Impact of a non-Gaussian density field on Sunyaev–Zeldovich observables

The main statistical properties of the Sunyaev–Zeldovich (S–Z) effect – the power spectrum, cluster number counts and angular correlation function – are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts. We do so for the usual Press & Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a  χ²  distributed density field. These three S–Z observables are found to be very significantly dependent on the choice of the mass function. The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S–Z observables in terms of cluster mass and redshift. The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function. A semiquantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.     

An analytic model for the epoch of halo creation

In this paper we describe the Bayesian link between the cosmological mass function and the distribution of times at which isolated haloes of a given mass exist. By assuming that clumps of dark matter undergo monotonic growth on the time-scales of interest, this distribution of times is also the distribution of 'creation' times of the haloes. This monotonic growth is an inevitable aspect of gravitational instability. The spherical top-hat collapse model is used to estimate the rate at which clumps of dark matter collapse. This gives the prior for the creation time given no information about halo mass. Applying Bayes' theorem then allows any mass function to be converted into a distribution of times at which haloes of a given mass are created. This general result covers both Gaussian and non-Gaussian models. We also demonstrate how the mass function and the creation time distribution can be combined to give a joint density function, and discuss the relation between the time distribution of major merger events and the formula calculated. Finally, we determine the creation time of haloes within three N -body simulations, and compare the link between the mass function and creation rate with the analytic theory.     

Abundance distribution of heavy nuclides up to a subnuclear density and the equation of state for gravitationally collapsing stars

With the use of the empirical equations for the binding and internal energies of heavy nuclides with density dependence, the abundance distribution of neutron-excess nuclides which appear up to some subnuclear density along the collapsing trajectory of a stellar core is calculated with an assumption of the -equilibrium condition. Free nucleons are simply assumed to be non-interacting, degenerate fermions. Using the abundance distribution thus derived, the quantities affecting the hydrodynamics of the core, such as the average mass number and the mass fraction of heavy nuclides, the mass fraction of free neutrons, the number of electrons per baryon, the average excitation energy per nucleus, and the entropy, the pressure and the adiabatic index of the system are then calculated. It is found that we can obtain fairly reliable values of these physical quantities by introducing only the magic and near-magic number nuclides in the calculation, and that the results are hardly affected by the difference of the semi-empirical nuclear mass formula we use to derive the binding energies of unknown nuclides and the formula of the grand partition function of these nuclides.     

The jeans problem for a thin galaxy in steady state

Given the potential, the equation of continuity and Poisson’s equation are solved for the variation perpendicular to the Galactic plane for a thin galaxy in a steady state. Simple expressions are obtained for the joint density function for the mass density and velocity, and for the distribution function for the velocity and its moments. These results are applied using a potential suggested by Woolley and Stewart (1967) and Whitley (1977), which is Camm’s potential due to an isothermal gas stratified in parallel layers, plus the potential due to the Galactic centre. The resulting velocity distribution is quite close to the normal distribution usually assumed and cannot be distinguished from it by the data. The mass density distribution fits the log (relative density) curves better than the Camm solution, especially at large distances from the Galactic plane. A formu1a, valid under conditions more genera1 than usual, is found for the total mass density in the neighbourhood of the Sun.     

Gas phase processes affecting galactic evolution

Gas processes affecting star formation are reviewed with an emphasis on gravitational and magnetic instabilities as a source of turbulence. Gravitational instabilities are pervasive in a multi-phase medium, even for sub-threshold column densities, suggesting that only an ISM with a pure-warm phase can stop star formation. The instabilities generate turbulence, and this turbulence influences the structure and timing of star formation through its effect on the gas distribution and density. The final trigger for star formation is usually direct compression by another star or cluster. The star formation rate is apparently independent of the detailed mechanisms for star formation, and determined primarily by the total mass of gas in a dense form. If the density distribution function is a log-normal, as suggested by turbulence simulations, then this dense gas mass can be calculated and the star formation rate determined from first principles. The results suggest that only 10^-4 of the ISM mass actively participates in the star formation process and that this fraction does so because its density is larger than 10⁵ cm^-3, at which point several key processes affecting dynamical equilibrium begin to break down. This revised version was published online in August 2006 with corrections to the Cover Date.     

Determination of the area and mass distribution of orbital debris fragments

An important factor in modeling the orbital debris environment is the loss rate of debris due to atmospheric drag and luni/solar perturbations. An accurate knowledge of the area-to-mass ratio of debris fragments is required for the calculation of the effect of atmospheric drag. In general, this factor is unknown and assumed values are used. However, this ratio can be calculated for fragments for which changes in the orbital elements due to atmospheric drag as a function of time are known. This is the inverse of the technique used to determine the atmospheric density from the decay of satellites with accurately known area-to-mass ratios. These kinds of propagation programs are routinely used in predicting the decay of an orbiting vehicle. In this work the area-to-mass ratio of about 2600 fragments arising from the breakup of 24 artificial satellites have been determined. An analysis of the data on about 200 objects (rocket bodies, scientific satellites, etc.) with known mass, size, and shape has also been made. The value of the radar cross-section (RCS), as measured by the Eglin radar operating at 70 cm wavelength, has been correlated to the effective area of these objects. The measurements of the area-to-mass ratio of these objects then provide a calibration of the actual to the calculated mass. It has been shown that the debris mean mass, m, is related to the mean effective area, A, by a power law relation, m = k A ^1.86. However, for a given effective area the mass distribution is very broad. Moreover, the cumulative mass distribution, N(>m), can be expressed as N(>m) = D(m + b), where D, b, and c are constants. The asymptotic slope, c, of low intensity explosions is on the average lower than the slope for high intensity explosions, but there is considerable spread of this slope in each class. Part of the flattening, as indicated by the finite value of the parameter, b, can be understood as arising out of the spread in the RCS values due to the tumbling motion of the fragments and effects related to the detectability of the fragment by the Eglin radar. It has been established that the mass in a given breakup calculated using this technique is in good agreement with the expected mass value. These results can be used in modeling the breakups of other artificial earth satellites and safety analysis.     

Some static relativistic compact charged fluid spheres in general relativity

In this work a new family of relativistic models of electrically charged compact star has been obtained by solving Einstein–Maxwell field equations with preferred form of one of the metric potentials and a suitable form of electric charge distribution function. The resulting equation of state (EOS) has been calculated. The relativistic stellar structure for matter distribution obtained in this work may reasonably models an electrically charged compact star whose energy density associated with the electric fields is on the same order of magnitude as the energy density of fluid matter itself (e.g. electrically charged bare strange stars). Based on the analytic model developed in the present work, the values of the relevant physical quantities have been calculated by assuming the estimated masses and radii of some well known strange star candidates like X-ray pulsar Her X-1, millisecond X-ray pulsar SAX J 1808.4-3658, and 4U 1820-30.     

The cloud-in-cloud problem for non-Gaussian density fields

The cloud-in-cloud problem is studied in the context of the extension to non-Gaussian density fields of the PS approach for the calculation of the mass function. As an example of a non-Gaussian probability distribution function (PDF), we consider the chi-square distribution with various degrees of freedom. We generate density fields in cubic boxes with periodic boundary conditions, and then determine the number of points considered collapsed at each scale through a hierarchy of smoothing windows. We find that the mass function we obtain differs from that predicted using the extended PS formalism, particularly for low values of σ and for those PDFs that differ most from a Gaussian.     

The statistics of weak lensing at small angular scales: probability distribution function

Weak gravitational lensing surveys have the potential to probe mass density fluctuation in the Universe directly. Recent studies have shown that it is possible to model the statistics of the convergence field at small angular scales by modelling the statistics of the underlying density field in the highly non-linear regime. We propose a new method to model the complete probability distribution function of the convergence field as a function of smoothing angle and source redshift. The model relies on a hierarchical ansatz for the behaviour of higher order correlations of the density field. We compare our results with ray-tracing simulations and find very good agreement over a range of smoothing angles. Whereas the density probability distribution function is not sensitive to the cosmological model, the probability distribution function for the convergence can be used to constrain both the power spectrum and cosmological parameters.     

基于Shapelets基函数的引力透镜质量重构方法

引力透镜效应是探测星系团物质分布的有效方法之一.目前,利用引力透镜数据重构星系团质量分布的主流方法可以分为两大类,即参数法和非参数法.在实际研究工作中,受限于质量模型假设和计算分辨率等方面的影响,现有的重构算法仍有诸多亟需解决的问题.基于Shapelets基函数的引力透镜质量重构方法通过基函数来实现引力透镜质量重构,使用Shapelets基函数分解引力透镜势,以引力透镜中多重像的位置和背景星系椭率畸变为限制条件来迭代求解基函数系数从而得到透镜体的质量分布.通过拟合一个模拟的NFW (Navarro,Frenk and White)透镜系统测试了新方法的可行性,结果表明新方法可以在整体上重构出透镜体的质量分布,并拟合出接近真实的源位置,能够为星系团质量测量提供一套灵活且高效的重构算法.     

Suprathermal particle distributions in space physics: Kappa distributions and entropy

The energization of a charged test-particle of mass m in contact with a large ensemble of charged particles of mass M at equilibrium is studied with the Fokker-Planck equation for Coulomb collisions and a quasi-linear diffusion operator for wave-particle interactions. The features of the nonequilibrium steady state velocity distribution of the test-particle system is studied as a function of the mass ratio m/M, and the relative strengths of the wave-particle interactions and Coulomb collisions. It is shown that the steady distribution function is not necessarily a Kappa distribution. The temperature of heavy minor ions given by the model is shown to vary linearly with the mass ratio as observed in the solar wind. The time evolution of the distribution function with and without the energization by wave-particle interactions is calculated and it is demonstrated that the Kullback relative entropy rather than the Tsallis nonextensive entropy rationalizes the results obtained.     

A simple model for the evolution of supermassive black holes and the quasar population

An empirically motivated model is presented for accretion-dominated growth of supermassive black holes (SMBH) in galaxies, and the implications are studied for the evolution of the quasar population in the Universe. We investigate the core aspects of the quasar population, including space density evolution, evolution of the characteristic luminosity, plausible minimum masses of quasars, the mass function of SMBH and their formation epoch distribution. Our model suggests that the characteristic luminosity in the quasar luminosity function arises primarily as a consequence of a characteristic mass scale above which there is a systematic separation between the black hole and the halo merging rates. At lower mass scales, black hole merging closely tracks the merging of dark haloes. When combined with a declining efficiency of black hole formation with redshift, the model can reproduce the quasar luminosity function over a wide range of redshifts. The observed space density evolution of quasars is well described by formation rates of SMBH above  ∼10⁸  M_⊙  . The inferred mass density of SMBH agrees with that found independently from estimates of the SMBH mass function derived empirically from the quasar luminosity function.     

Mass Estimation in the Globular Cluster M4

Using the kinematical data we try to solve the Jeans equation to provide mass estimation of the globular cluster M4. The proper motion data provides two independent velocity dispersion profiles which we need in this estimation. Moreover, we have calculated the density distribution function and checked the anisotropy of the velocity dispersion. The comparison of the dynamical mass with the estimation stellar mass for M4 gives the best estimation of the mass inside a spherical shell centered on the cluster, with a radius corresponding to 840 arcsec on the sky, as 4-5× 10⁴ M . This revised version was published online in July 2006 with corrections to the Cover Date.     

X-ray group and cluster mass profiles in MOND: unexplained mass on the group scale

Although very successful in explaining the observed conspiracy between the baryonic distribution and the gravitational field in spiral galaxies without resorting to dark matter (DM), the modified Newtonian dynamics (MOND) paradigm still requires DM in X-ray bright systems. Here, to get a handle on the distribution and importance of this DM, and thus on its possible form, we deconstruct the mass profiles of 26 X-ray emitting systems in MOND, with temperatures ranging from 0.5 to 9 keV. Initially, we compute the MOND dynamical mass as a function of radius, then subtract the known gas mass along with a component of galaxies which include the cD galaxy with   M / L _K = 1  . Next, we test the compatibility of the required DM with ordinary massive neutrinos at the experimental limit of detection  ( m _ν= 2 eV)  , with density given by the Tremaine–Gunn limit. Even by considering that the neutrino density stays constant and maximal within the central 100 or 150 kpc (which is the absolute upper limit of a possible neutrino contribution there), we show that these neutrinos can never account for the required DM within this region. The natural corollary of this finding is that, whereas clusters  ( T ≳ 3 keV)  might have most of their mass accounted for if ordinary neutrinos have a 2 eV mass, groups  ( T ≲ 2 keV)  cannot be explained by a 2 eV neutrino contribution. This means that, for instance, cluster baryonic dark matter (CBDM, Milgrom) or even sterile neutrinos would present a more satisfactory solution to the problem of missing mass in MOND X-ray emitting systems.     

Kinetic models for the exospheres of Jupiter and Saturn

The exospheric theory based on the Kappa velocity distribution function (VDF) is used to model the exosphere of the giant planets Jupiter and Saturn. Such Kappa velocity distribution functions with an enhanced population of suprathermal particles are indeed often observed in space plasmas and in the space environment of the planets. The suprathermal particles have significant effects on the escape flux, density and temperature profiles of the particles in the exosphere of the giant planets. The polar wind flux becomes several orders larger when suprathermal electrons are considered, so that the planetary ionosphere becomes then a significant source for their inner magnetosphere. Moreover, the number density of the particles decreases slower as a function of the altitude when a Kappa distribution is considered instead of a Maxwellian one. Two-dimensional maps of density are calculated for typical values of the temperatures. The exospheric formalism is also applied to study the escape flux from the exospheres of Io and Titan, respectively, moons of Jupiter and Saturn.     

Turbulent Fragmentation and Star Formation

We review the main results from recent numerical simulations of turbulent fragmentation and star formation. Specifically, we discuss the observed scaling relationships, the “quiescent” (subsonic) nature of many star-forming cores, their energy balance, their synthesized polarized dust emission, the ages of stars associated with the molecular gas from which they have formed, the mass spectra of clumps, and the density and column density probability distribution function of the gas. We then give a critical discussion on recent attempts to explain and/or predict the star formation efficiency and the stellar initial mass function from the statistical nature of turbulent fields. Finally, it appears that turbulent fragmentation alone cannot account for the final stages of fragmentation: although the turbulent velocity field is able to produce filaments, the spatial distribution of cores in such filaments is better explained in terms of gravitational fragmentation.     

Inertia tensor and tidal transfer of angular momentum for special ellipsoidal mass distributions

This paper deals with ellipsoidal mass distributions made by concentric, coaxial, and similar shells, whose density is specified for two special situations both including the whole range between the limiting cases of homogeneous systems, on one side, and Roche generalized systems, on the other side. The components of the related inertia tensor are calculated and the dependence of the spin growth on the density distribution is shown, concerning tidal transfer of angular momentum to both virialized proto-galaxies and encountering galaxies.     

On the motion of a mass point inside a rotating inhomogeneous ellipsoidal body

The plane motion of a mass point inside an inhomogeneous rotating ellipsoidal body with a homothetic density distribution is considered. The force function of the problem is expanded in terms of the ellipsoid's second eccentricities up to the fourth order, which are taken as small parameters. We derive an expression for the perturbing function and solve the equations of perturbed motion in orbital elements.     

Shaping of Planetary Nebulae by Dust Driven AGB Winds

We have investigated the interacting winds model (IWM) in which the shapes of elliptical Planetary Nebulae (PNe) can be explained by the asymmetric mass loss produced by a rotating AGB star. The mass loss mechanism is based on a dust driven wind calculated for stationary situations. Already for small rotation rates of the AGB star we get a significantly angle-dependent mass loss which is concentrated towards the equatorial plane. This pole to equator density variation in the space surrounding the star influences the shape of the later developed PN. We compare these theoretical shapes with observed PNe and for some objects with well known quantities our model can fit the observations quite well. This revised version was published online in July 2006 with corrections to the Cover Date.