Space–time foam and cosmic-ray interactions

It has been proposed that propagation of cosmic-rays at extreme-energy may be sensitive to Lorentz-violating metric fluctuations (“foam”). We investigate the changes in interaction thresholds for cosmic-rays and gamma-rays interacting on the CMB and IR backgrounds, for a class of stochastic models of space–time foam. The strength of the foam is characterized by the factor (E/M_P)^a, where a is a phenomenological suppression parameter. We find that there exists a critical value of a (dependent on the particular reaction: a_crit3 for cosmic-rays, 1 for gamma-rays), below which the threshold energy can only be lowered, and above which the threshold energy may be raised, but at most by a factor of two. Thus, it does not appear possible in this class of models to extend cosmic-ray spectra significantly beyond their classical absorption energies. However, the lower thresholds resulting from foam may have signatures in the cosmic-ray spectrum. In the context of this foam model, we find that cosmic-ray energies cannot exceed the fundamental Planck scale, and so set a lower bound of 10⁸ TeV for the scale of gravity. We also find that suppression of p→pπ⁰ and γ→e⁻e⁺ “decays” favors values aa_crit. Finally, we comment on the apparent non-conservation of particle energy–momentum, and speculate on its re-emergence as dark energy in the foamy vacuum.     

Basic N-body modelling of the evolution of globular clusters – I. Time scaling

We consider the use of N -body simulations for studying the evolution of rich star clusters (i.e. globular clusters).The dynamical processes included in this study are restricted to gravitational (point-mass) interactions, the steady tidal field of a galaxy, and instantaneous mass loss resulting from stellar evolution. With evolution driven by these mechanisms, it is known that clusters fall roughly into two broad classes: those that dissipate promptly in the tidal field, as a result of mass loss; and those that survive long enough for their evolution to become dominated by two-body relaxation.
The time-scales of the processes we consider scale in different ways with the number of stars in the simulation, and the main aim of the paper is to suggest how the scaling of a simulation should be done so that the results are representative of the evolution of a 'real' cluster. We investigate three different ways of scaling time. One of these is appropriate to the first type of cluster, i.e. those that dissipate rapidly; similarly, a second scaling is appropriate only to the second (relaxation-dominated) type. We also develop a hybrid scaling, which is a satisfactory compromise for both types of cluster. Finally we present evidence that the widely used Fokker–Planck method produces models that are in good agreement with N -body models of those clusters that are relaxation-dominated, at least for N -body models with several thousand particles, but that the Fokker–Planck models evolve too fast for clusters that dissipate promptly.     

Cosmological parameter estimation and the spectral index from inflation

Accurate estimation of cosmological parameters from microwave background anisotropies requires high-accuracy understanding of the cosmological model. Normally, a power-law spectrum of density perturbations is assumed, in which case the spectral index n can be measured to around ± 0.004 using microwave anisotropy satellites such as MAP Planck . However, inflationary models generically predict that the spectral index n of the density perturbation spectrum will be scale-dependent. We carry out a detailed investigation of the measurability of this scale dependence by Planck , including the influence of polarization on the parameter estimation. We also estimate the increase in the uncertainty in all other parameters if the scale dependence has to be included. This increase applies even if the scale dependence is too small to be measured, unless it is assumed absent. We study the implications for inflation models, beginning with a brief examination of the generic slow-roll inflation situation, and then move to a detailed examination of a recently devised hybrid inflation model for which the scale dependence of n may be observable.     

Extragalactic source counts and contributions to the anisotropies of the cosmic microwave background: predictions for the Planck Surveyor mission

We present predictions for the counts of extragalactic sources, the contributions to fluctuations and their angular power spectrum in each channel foreseen for the Planck Surveyor (formerly COBRAS/SAMBA ) mission. The contribution to fluctuations owing to clustering of both radio and far-IR sources is found to be generally small in comparison with the Poisson term; however the relative importance of the clustering contribution increases and may eventually become dominant if sources are identified and subtracted down to faint flux limits. The central Planck frequency bands are expected to be 'clean': at high galactic latitude (| b | > 20°), where the reduced galactic noise does not prevent the detection of the extragalactic signal, only a tiny fraction of pixels is found to be contaminated by discrete extragalactic sources. Moreover, the 'flat' angular power spectrum of fluctuations resulting from extragalactic sources substantially differs from that of primordial fluctuations; therefore, the removal of contaminating signals is eased even at frequencies where point sources give a sizeable contribution to the foreground noise.     

Planck and re-ionization history: a model selection view

We use Bayesian model selection tools to forecast the Planck satellite's ability to distinguish between different models for the re-ionization history of the Universe, using the large angular scale signal in the cosmic microwave background polarization spectrum. We find that Planck is not expected to be able to distinguish between an instantaneous re-ionization model and a two-parameter smooth re-ionization model, except for extreme values of the additional re-ionization parameter. If it cannot, then it will be unable to distinguish between different two-parameter models either. However, Bayesian model averaging will be needed to obtain unbiased estimates of the optical depth to re-ionization. We also generalize our results to a hypothetical future cosmic variance limited microwave anisotropy survey, where the outlook is more optimistic.     

Thermodynamic Modeling of Developed Structural Turbulence Taking into Account Fluctuations of Energy Dissipation

A thermodynamic approach to the construction of a phenomenological macroscopic model of developed turbulence in a compressible fluid is considered with regard for the formation of space–time dissipative structures. A set of random variables were introduced into the model as internal parameters of the turbulent–chaos subsystem. This allowed us to obtain, by methods of nonequilibrium thermodynamics, the kinetic Fokker–Planck equation in the configuration space. This equation serves to determine the temporary evolution of the conditional probability distribution function of structural parameters pertaining to the cascade process of fragmentation of large-scale eddies and temperature inhomogeneities and to analyze Markovian stochastic processes of transition from one nonequilibrium stationary turbulent-motion state to another as a result of successive loss of stability caused by a change in the governing parameters. An alternative method for investigating the mechanisms of such transitions, based on the stochastic Langevin-type equation intimately related to the derived kinetic equation, is also considered. Some postulates and physical and mathematical assumptions used in the thermodynamic model of structurized turbulence are discussed in detail. In particular, we considered, using the deterministic transport equation for conditional means, the cardinal problem of the developed approach—the possibility of the existence of asymptotically stable stationary states of the turbulent-chaos subsystem. Also proposed is the nonequilibrium thermodynamic potential for internal coordinates, which extends the well-known Boltzmann–Planck relationship for equilibrium states to the nonequilibrium stationary states of the representing ensemble. This potential is shown to be the Lyapunov function for such states. The relation is also explored between the internal intermittence in the inertial interval of scales and the fluctuations of the energy of dissipation. This study is aimed at constructing representative models of natural space environments. It develops a synergetic approach to modeling the structurized turbulence of astrophysical and geophysical systems, which was proposed by the author in previous papers (Kolesnichenko, 2002, 2003).     

GAUGE: the GrAnd Unification and Gravity Explorer

The GAUGE (GrAnd Unification and Gravity Explorer) mission proposes to use a drag-free spacecraft platform onto which a number of experiments are attached. They are designed to address a number of key issues at the interface between gravity and unification with the other forces of nature. The equivalence principle is to be probed with both a high-precision test using classical macroscopic test bodies, and, to lower precision, using microscopic test bodies via cold-atom interferometry. These two equivalence principle tests will explore string-dilaton theories and the effect of space–time fluctuations respectively. The macroscopic test bodies will also be used for intermediate-range inverse-square law and an axion-like spin-coupling search. The microscopic test bodies offer the prospect of extending the range of tests to also include short-range inverse-square law and spin-coupling measurements as well as looking for evidence of quantum decoherence due to space–time fluctuations at the Planck scale.     

Evolution of multimass globular clusters in the Galactic tidal field with the effects of velocity anisotropy

We study the evolution of globular clusters with mass spectra under the influence of the steady Galactic tidal field, including the effects of velocity anisotropy. Similarly to single-mass models, velocity anisotropy develops as the cluster evolves, but the degree of anisotropy is much smaller than in isolated clusters. Except for very early epochs of the cluster evolution, the velocity distributions of nearly all mass components become tangentially anisotropic at the outer parts. We examine how the mass function (MF) changes in time. Specifically, we find that the power-law index of the MF decreases monotonically with the total mass of the cluster, in agreement with previous findings based on isotropic models or N -body studies. This is also consistent with the behaviour of the observed slopes of MFs for a limited number of clusters. We attempt to compare our results with multimass King models, although it is almost impossible to fit the entire density profiles for all mass components. When the MF is fixed, the central densities of individual components show significant differences between Fokker–Planck and King models. We obtain 'best-fitting' multimass King models, for which the central density of individual components as well as the total density distribution agrees with the Fokker–Planck models by adjusting the MF. The MFs obtained in this way closely resemble the MF within the half-mass radius of the Fokker–Planck result. Also, we find that the local MFs predicted by Fokker–Planck calculations vary more rapidly with radius than best-fitting multimass King models. The projected velocity profiles for anisotropic models show significant flattening toward the tidal radius compared with the isotropic model. This is caused by the fact that the tangential velocity dispersion becomes dominant at the outer parts. Such a behaviour of velocity profile appears to be consistent with the observed profiles of the collapsed cluster M15.     

Fast parameter estimation from the cosmic microwave background power spectrum

The statistical properties of a map of the primary fluctuations in the cosmic microwave background (CMB) may be specified to high accuracy by a few thousand power spectra measurements, provided the fluctuations are Gaussian, yet the number of parameters relevant for the CMB is probably no more than ∼10–20. Consequently, there is a large degree of redundancy in the power spectrum data. In this paper, we show that the moped data compression technique can reduce the CMB power spectrum measurements to ∼10–20 numbers (one for each parameter), from which the cosmological parameters can be estimated virtually as accurately as from the complete power spectrum. Combined with recent advances in the speed of generation of theoretical power spectra, this offers opportunities for very fast parameter estimation from real and simulated CMB skies. The evaluation of the likelihood itself, at Planck resolution, is speeded up by factors up to ∼10⁸, ensuring that this step will not be the dominant part of the data analysis pipeline.     

Solar variability and stochastic effects on climate

The effect of fluctuations on a simple energy balance model is examined. A new, long characteristic time scale referring to the passage between different stable climatic states is identified. It is shown that a weak external forcing whose period is comparable to this scale enables the system to switch between different states with a high probability. The connection with glaciation cycles is pointed out.Proceedings of the 14th ESLAB Symposium on Physics of Solar Variations, 16–19 September 1980, Scheveningen, The Netherlands.     

Global fits for the spectral index of the cosmological curvature perturbation

Best-fitting values of the spectral index of the curvature perturbation are presented, assuming the ΛCDM cosmology. Apart from the spectral index, the parameters are the Hubble parameter, the total matter density and the baryon density. The data points are intended to represent all measurements that are likely to affect the result significantly. The cosmic microwave anisotropy is represented by the COBE normalization, and heights of the first and second peaks are given by the latest Boomerang and Maxima data. The slope of the galaxy correlation function and the matter density contrast on the 8  h ⁻¹ Mpc scale are each represented by a data point, as are the expected values of the Hubble parameter and matter density. The 'low-deuterium' nucleosynthesis value of the baryon density provides a final data point, the fit giving a value higher by about one standard deviation. The reionization epoch is calculated from the model by assuming that it corresponds to the collapse of a fraction f ≳10⁻⁴ of matter. We consider the case of a scale-independent spectral index, and also the scale-dependent spectral index predicted by running mass models of inflation. In the former case, the result is compared with the prediction of models of inflation based on effective field theory, in which the field value is small on the Planck scale. A detailed comparison is made with other fits, and other approaches to the comparison with theory.     

Processing and compression of noise-dominated data: application to the cosmic microwave background data on board the Planck satellite

We address the problem of encoding and compressing data dominated by noise. Information is decomposed into 'reference' sequences plus arrays containing noisy differences susceptible to being described by a known probability distribution. One can then give reliable estimates of the optimal compression rates by estimating the corresponding Shannon entropy. As a working example, this idea is applied to an idealized model of the cosmic microwave background (CMB) data on board the Planck satellite. Data reduction is a critical issue in space missions because the total information that can be downloaded to Earth is sometimes limited by telemetry allocation. Similar limitations might arise in remotely operated ground based telescopes. This download-rate limitation could reduce the amount of diagnostics sent on the stability of the instruments and, as a consequence, curb the final sensitivity of the scientific signal. Our proposal for Planck consists of taking differences of consecutive circles at a given sky pointing. To a good approximation, these differences could be made independent of the external signal, so that they are dominated by thermal (white) instrumental noise, which is simpler to model than the sky signal. Similar approaches can be found in other individual applications. Generic simulations and analytical predictions show that high compression rates,     can be obtained with minor or zero loss of sensitivity. Possible effects of digital distortion are also analysed. The proposed scheme is flexible and reliable enough to be optimized in relation to other critical aspects of the corresponding application. For Planck , this study constitutes an important step towards a more realistic modelling of the final sensitivity of the CMB temperature anisotropy maps.     

On the Reliability of Merger-Trees and the Mass-Growth Histories of Dark Matter Haloes

We have used merger-trees realizations to study the formation of dark matter haloes. The construction of merger-trees is based on three different pictures about the formation of structures in the Universe. These pictures include the spherical collapse (SC), the ellipsoidal collapse (EC) and the non-radial collapse (NR). The reliability of merger-trees has been examined comparing their predictions related to the distribution of the number of progenitors, as well as the distribution of formation times, with the predictions of analytical relations. The comparison yields a very satisfactory agreement. Subsequently, the mass-growth histories (MGH) of haloes have been studied and their formation scale factors have been derived. This derivation has been based on two different definitions that are (a) the scale factor when the halo reaches half its present day mass and (b) the scale factor when the mass-growth rate falls below some specific value. Formation scale factors follow approximately power laws of mass. It has also been shown that MGHs are in good agreement with models proposed in the literature that are based on the results of N-body simulations. The agreement is found to be excellent for small haloes but, at the early epochs of the formation of large haloes, MGHs seem to be steeper than those predicted by the models based on N-body simulations. This rapid growth of mass of heavy haloes is likely to be related to a steeper central density profile indicated by the results of some N-body simulations.     

Was the Planck era of the universe a ‘cold’ or a ‘hot’ stage of the Big Bang?

Two models are presented of the Planck era in the Big Bang.The first of them starts from the state of an essentially cold and non-expanding universe (assumed to be closed). It is shown that the conditions in the Planck era rendered the total annihilation of the matter with the Planck density absolutely necessary. The annihilation was, further, a necessary condition for the expansion of the post-Planck universe, which, by itself, also was a necessity. The hot stage of the Big Bang in the post-Planck era was, in fact, a necessary consequence resulting from the scheme.The second model is motivated by the weak points of the first one. It is concluded that, owing to the extreme physical conditions, the cosmic medium within the Planck era could only consist of ultrahighly energetic photons.Finally, it is proposed that the second model is to be preferred.     

综合脉冲星时的小波分析算法

由单脉冲星定义的脉冲星时受几种噪声源的影响,为削弱这些影响得到一种更稳定的时间尺度,采取将多颗脉冲星定义的单脉冲星时进行综合的分析方法得到综合脉冲星时．用两种方法:经典加权算法和小波分解算法,对PSRB1855 09和PSR B1937 21两颗脉冲星进行综合计算并做出比较,经典加权算法无法兼顾脉冲星不同频率上的稳定度,而用小波分析的方法对PSR B1855 09和PSR B1937 21两颗脉冲星进行综合,能够有效抑制不同频率上噪声的影响,达到更好的效果．     

Detection/estimation of the modulus of a vector. Application to point-source detection in polarization data

Given a set of images, whose pixel values can be considered as the components of a vector, it is interesting to estimate the modulus of such a vector in some localized areas corresponding to a compact signal. For instance, the detection/estimation of a polarized signal in compact sources immersed in a background is relevant in some fields like astrophysics. We develop two different techniques, one based on the Neyman–Pearson lemma, the Neyman–Pearson filter (NPF), and another based on pre-filtering before fusion, the filtered fusion (FF), to deal with the problem of detection of the source and estimation of the polarization given two or three images corresponding to the different components of polarization (two for linear polarization, three including circular polarization). For the case of linear polarization, we have performed numerical simulations on two-dimensional patches to test these filters following two different approaches (a blind and a non-blind detection), considering extragalactic point sources immersed in cosmic microwave background (CMB) and non-stationary noise with the conditions of the 70 GHz Planck channel. The FF outperforms the NPF, especially for low fluxes. We can detect with the FF extragalactic sources in a high noise zone with fluxes      Jy for (blind/non-blind) detection and in a low noise zone with fluxes      Jy for (blind/non-blind) detection with low errors in the estimated flux and position.     

The build-up of haloes within Press–Schechter theory

Modelling the build-up of haloes is important for linking the formation of galaxies with cosmological models. A simple model of halo growth is provided by Press–Schechter (PS) theory, where the initial field of density fluctuations is smoothed using spherically symmetric filters centred on a given position to obtain information about the likelihood of later collapse on varying scales. In this paper the predicted halo mass growth is compared for three filter shapes: Gaussian, top-hat and sharp k -space. Preliminary work is also presented analysing the build-up of haloes within numerical simulations using a friends-of-friends group finder. The best-fit to the simulation mass function was obtained using PS theory with a top-hat filter. By comparing both the backwards conditional mass function, which gives the distribution of halo progenitors, and the distribution of halo mergers in time, the build-up of haloes in the simulations is shown to be better fitted by PS theory with a sharp k -space filter. This strengthens previous work, which also found the build-up of haloes in simulations to be well matched to PS theory with a sharp k -space filter by providing a direct comparison of different filters and by extending the statistical tools used to analyse halo mass growth. The usefulness of this work is illustrated by showing that the cosmological evolution in the proportion of haloes that have undergone recent merger is predicted to be independent of mass and power spectrum and to only depend upon cosmology. Recent results from observations of field galaxies are shown to match the evolution expected, but are not sufficiently accurate to distinguish usefully between cosmological parameters.     

Dynamical evolution of rotating stellar systems – II. Post-collapse, equal-mass system

We present the first post-core-collapse models of initially rotating star clusters, using the numerical solution of an orbit-averaged 2D Fokker–Planck equation. Based on the code developed by Einsel & Spurzem, we have improved the speed and the stability and included the steady three-body binary heating source. We have confirmed that rotating clusters, whether they are in a tidal field or not, evolve significantly faster than non-rotating ones. Consequences for the observed shapes, density distribution and kinematic properties of young and old star clusters are discussed. The results are compared with gaseous and 1D Fokker–Planck models in the non-rotating case.     

Evolution of galaxy clustering

We study the evolution of the correlation function of dark matter haloes in the CDM class of models. We show that the halo correlation function does not evolve in proportion with the correlation function of the underlying mass distribution. The earliest haloes to collapse, which correspond to rare peaks in the density field, cluster very strongly. The amplitude of the halo correlation function decreases from its initial, large, value. This decrease continues until the average peaks have collapsed, after which the amplitude grows slowly. This behaviour is shown to be generic and the epoch of minimum amplitude depends only on the rms  fluctuations in mass at the relevant scale and, to a much smaller extent, on the slope of the power spectrum at that scale. We discuss the relevance of this result for interpretation of observations of galaxy and quasar clustering.     

The correlation of peaks in the microwave background

We present accurate small-angle predictions of the correlation function of hotspots in the microwave background radiation for Gaussian theories such as those predicted in most inflation models. The correlation function of peaks above a certain threshold depends only on the threshold and the power spectrum of temperature fluctuations. Since these are both potentially observable quantities in a microwave background map, there are no adjustable parameters in the predictions. These correlations should therefore provide a powerful test of the Gaussian hypothesis, and provide a useful discriminant between inflation and topological defect models such as the cosmic string model. The correlations have a number of oscillatory features, which should be detectable at high signal-to-noise ratio with future satellite experiments such as MAP and Planck .