We conduct a rigorous examination of the nearby red supergiant $\alpha$ Orionis, or Betelgeuse, by drawing on the synthesis of new observational data and three different modeling techniques. Our observational results include the release of new, processed photometric measurements collected with the space-based SMEI instrument prior to Betelgeuse's recent, unprecedented dimming event. Our theoretical predictions include self-consistent results from multi-timescale evolutionary, oscillatory, and hydrodynamic simulations conducted with the Modules for Experiments in Stellar Astrophysics (MESA) software suite. Significant outcomes of our modeling efforts include a precise prediction for the star's radius: $750^{+62}_{-30}\,R_{\odot}$. In concert with additional constraints, this allows us to derive a new, independent distance estimate of $165^ {+16}_{-8}$\,pc and a parallax of $\pi=6.06^{+0.31}_{-0.52}$\,mas, in good agreement with \textit{Hipparcos} but less so with recent radio measurements. Seismic results from both perturbed hydrostatic and evolving hydrodynamic simulations constrain the period and driving mechanisms of Betelgeuse's dominant periodicities in new ways. Our analyses converge to the conclusion that Betelgeuse's $388$ day period is the result of pulsation in the fundamental mode, driven by the $\kappa$-mechanism. Grid-based hydrodynamic modeling reveals that the behavior of the oscillating envelope is mass-dependent, and likewise suggests that the non-linear pulsation excitation time could serve as a mass constraint. Our results corroborate recent conclusions that Betelgeuse is the outcome of a past merger. We place it definitively in the core helium-burning phase near the base of the red supergiant branch, and we report a present-day mass of $16.5$--$19 ~M_{\odot}$---slightly lower than typical literature values.

The last decade has seen a rapid development in asteroseismology thanks to the CoRoT and Kepler missions. With more detailed asteroseismic observations available, it is becoming possible to infer exactly how oscillations are driven and dissipated in solar-type stars. We have carried out three-dimensional (3D) stellar atmosphere simulations together with one-dimensional (1D) stellar structural models of key benchmark turn-off and subgiant stars to study this problem from a theoretical perspective. Mode excitation and damping rates are extracted from 3D and 1D stellar models based on analytical expressions. Mode velocity amplitudes are determined by the balance between stochastic excitation and linear damping, which then allows the estimation of the frequency of maximum oscillation power, $\nu_{\max}$, for the first time based on ab initio and parameter-free modelling. We have made detailed comparisons between our numerical results and observational data and achieved very encouraging agreement for all of our target stars. This opens the exciting prospect of using such realistic 3D hydrodynamical stellar models to predict solar-like oscillations across the HR-diagram, thereby enabling accurate estimates of stellar properties such as mass, radius and age.

We use self-similarity in N-body simulations of scale-free models to test for resolution dependence in the mass function and two-point correlation functions of dark matter halos. We use 1024$^3$ particle simulations performed with the ABACUS N-body code and define halos with two different algorithms, "friends of friends" (FOF) and ROCKSTAR. The FOF mass functions show a systematic deviation from self-similarity which is explained by a resolution dependence of the FOF mass assignment previously reported in the literature. Evidence for convergence is observed only starting from halos of several thousand particles, and mass functions are overestimated by as much as 20-25 percent for halos of 50 particles. The mass function of ROCKSTAR halos, on the other hand, shows good convergence from of order 50 to 100 particles per halo, with no detectable evidence at the one percent level of any systematic dependence for larger particle number. Tests show that the mass unbinding procedure in ROCKSTAR is the key factor in obtaining this much improved resolution. Applying the same analysis to the halo-halo two point correlation function, we find again strong evidence for convergence only for ROCKSTAR halos, at separations sufficiently large so that halos do not overlap. At these separations we can exclude dependence on resolution at the few percent level once halos have of order 50 to 100 particles. At smaller separations results are not converged even at significantly larger particle number, and bigger simulations would be required to establish the resolution required for convergence.

We demonstrate that testing for self-similarity in scale-free simulations provides an excellent tool to quantify the resolution at small scales of cosmological N-body simulations. Analysing two-point correlation functions measured in simulations using ABACUS, we show how observed deviations from self-similarity reveal the range of time and distance scales in which convergence is obtained. While the well-converged scales show accuracy below 1 percent, our results show that, with a small force softening length, the spatial resolution is essentially determined by the mass resolution. At later times the lower cut-off scale on convergence evolves in comoving units as $a^{-1/2}$ ($a$ being the scale factor), consistent with a hypothesis that it is set by two-body collisionality. A corollary of our results is that N-body simulations, particularly at high red-shift, contain a significant spatial range in which clustering appears converged with respect to the time-stepping and force softening but has not actually converged to the physical continuum result. The method developed can be applied to determine the resolution of any clustering statistic and extended to infer resolution limits for non-scale-free simulations.

We present 1D-MESA2HYDRO-3D, an open source, Python-based software tool that provides an accessible means of generating physically motivated initial conditions (ICs) for hydrodynamical simulations from 1-D stellar structure models. We test 1D-MESA2HYDRO-3D on five stellar models generated with the MESA stellar evolution code and verify its capacity as an IC generator with the Phantom smoothed-particle hydrodynamics code \citep{MESAIV, Phantom}. Consistency between the input density profiles, the 1D-MESA2HYDRO-3D-rendered particle distributions, and the state of the distributions after evolution over $10$ dynamical timescales is found for model stars ranging in structure and density from a radially extended supergiant to a white dwarf.

Most aspects of stellar evolution proceed far too slowly to be directly observable in a single star on human timescales. The thermally pulsing asymptotic giant branch is one exception. The combination of state-of-the-art modelling techniques with data assimilated from observations collected by amateur astronomers over many decades provide, for the first time, the opportunity to identify a star occupying precisely this evolutionary stage. In this study, we show that the rapid pulsation period change and associated reduction in radius in the bright, northern variable star T Ursae Minoris are caused by the recent onset of a thermal pulse. We demonstrate that T UMi transitioned into a double-mode pulsation state, and we exploit its asteroseismic features to constrain its fundamental stellar parameters. We use evolutionary models from MESA and linear pulsation models from GYRE to track simultaneously the structural and oscillatory evolution of models with varying mass. We apply a sophisticated iterative sampling scheme to achieve time resolution $\le10$ years at the onset of the relevant thermal pulses. We report initial mass of $2.0\pm0.15\, \mathrm{M}_\odot$ and an age of $1.17 \pm 0.21$ Gyr for T UMi. This is the most precise mass and age determination for a single asymptotic giant branch star ever obtained. The ultimate test of our models will be the continued observation of its evolution in real time: we predict that the pulsation periods in T UMi will continue shortening for a few decades before they rebound and begin to lengthen again, as the star expands in radius.

The relation between halo mass, M, and concentration, c, is a critical component in our understanding of the structure of dark matter halos. While numerous models for this relation have been proposed, almost none of them attempt to derive the evolution of the relation analytically. We build on previous efforts to model the c-M relation as a function of physical parameters such as the peak height, $\nu$, and the effective power spectrum slope, $n_{\rm eff}$, which capture the dependence of $c$ on halo mass, redshift, and cosmology. We present three major improvements over previous models. First, we derive an analytical expression for the c-M relation that is valid under the assumption of pseudo-evolution, i.e., assuming that the density profiles of halos are static in physical coordinates while the definition of their boundary evolves. We find that this ansatz is highly successful in describing the evolution of the low-mass end of the c-M relation. Second, we employ a new physical variable, the effective exponent of linear growth, $\alpha_{\rm eff}$, to parameterize deviations from an Einstein-de Sitter expansion history. Third, we combine an updated definition of $n_{\rm eff}$ with the additional dependence on $\alpha_{\rm eff}$ and propose a phenomenological extension of our analytical framework to include all halo masses. This semianalytical model matches simulated concentrations in both scale-free models and LambdaCDM to 5% accuracy with very few exceptions and differs significantly from all previously proposed models. We present a publicly available code to compute the predictions of our model in the python toolkit Colossus, including updated parameters for the model of Diemer and Kravtsov.

The evolution of self-gravitating systems, and long-range interacting systems more generally, from initial configurations far from dynamical equilibrium is often described as a simple two phase process: a first phase of violent relaxation bringing it to a quasi-stationary state in a few dynamical times, followed by a slow adiabatic evolution driven by collisional processes. In this context the complex spatial structure evident, e.g., in spiral galaxies is understood either in terms of instabilities of quasi-stationary states, or a result of dissipative non-gravitational interactions. We illustrate here, using numerical simulations, that purely self-gravitating systems evolving from quite simple initial configurations can in fact give rise easily to structures of this kind of which the lifetime can be large compared to the dynamical characteristic time, but short compared to the collisional relaxation time scale. More specifically, for a broad range of non-spherical and non-uniform rotating initial conditions, gravitational relaxation gives rise quite generically to long-lived non-stationary structures of a rich variety, characterized by spiral-like arms, bars and even ring-like structures in special cases. These structures are a feature of the intrinsically out-of-equilibrium nature of the system's collapse, associated with a part of the system's mass while the bulk is well virialized. They are characterized by predominantly radial motions in their outermost parts, but also incorporate an extended flattened region which rotates coherently about a well virialized core of triaxial shape with an approximately isotropic velocity dispersion. We discuss the possible relevance of these simple toy models to the observed structure of real galaxies emphasizing the difference between dissipative and dissipationless disc formation.

The bright, nearby binary $\alpha$ Centauri provides an excellent laboratory for testing stellar evolution models, as it is one of the few stellar systems for which we have high-precision classical (mass, radius, luminosity) and asteroseismic ($p$-mode) observations. Stellar models are created and fit to the classical and seismic observations of both stars by allowing for the free variation of convective mixing length parameter $\alpha_{\text{MLT}}$. This system is modeled using five different sets of assumptions about the physics governing the stellar models. There are 31 pairs of tracks (out of ${\sim} 150,000$ generated) which fit the classical, binary, and seismic observational constraints of the system within $3\,\sigma$. Models with each tested choice of input physics are found to be viable, but the optimal mixing lengths for Cen A and Cen B remain the same regardless of the physical prescription. The optimal mixing lengths are $\alpha_{\text{MLT,A}} /\alpha_{\odot}= 0.932$ and $\alpha_{\text{MLT,B}}/\alpha_{\odot} = 1.095$. That Cen A and Cen B require sub- and super-solar mixing lengths, respectively, to fit the observations is a trend consistent with recent findings, such as in Kervella et al. (2017), Joyce and Chaboyer (2018), and Viani et al. (2018). The optimal models find an age for $\alpha$ Centauri of $5.3 \pm 0.3$ Gyr.

Theoretical stellar evolution models are constructed and tailored to the best known, observationally derived characteristics of metal-poor ([Fe/H]$\sim-2.3$) stars representing a range of evolutionary phases: subgiant HD140283, globular cluster M92, and four single, main sequence stars with well-determined parallaxes: HIP46120, HIP54639, HIP106924, and WOLF1137. It is found that the use of a solar-calibrated value of the mixing length parameter $\alpha_{\text{MLT}}$ in models of these objects is ineffective at reproducing their observed properties. Empirically calibrated values of $\alpha_{\text{MLT}}$ are presented for each object, accounting for uncertainties in the input physics employed in the models. It is advocated that the implementation of an adaptive mixing length is necessary in order for stellar evolution models to maintain fidelity in the era of asteroseismic observations.

We present the results of a multisite photometric observing campaign on the rapidly oscillating Ap (roAp) star 2MASS 16400299-0737293 (J1640; $V=12.7$). We analyse photometric $B$ data to show the star pulsates at a frequency of $151.93$ d$^{-1}$ ($1758.45 \mu$Hz; $P=9.5$ min) with a peak-to-peak amplitude of 20.68 mmag, making it one of the highest amplitude roAp stars. No further pulsation modes are detected. The stellar rotation period is measured at $3.6747\pm0.0005$ d, and we show that rotational modulation due to spots is in anti-phase between broadband and $B$ observations. Analysis and modelling of the pulsation reveals this star to be pulsating in a distorted quadrupole mode, but with a strong spherically symmetric component. The pulsational phase variation in this star is suppressed, leading to the conclusion that the contribution of $\ell>2$ components dictate the shape of phase variations in roAp stars that pulsate in quadrupole modes. This is only the fourth time such a strong pulsation phase suppression has been observed, leading us to question the mechanisms at work in these stars. We classify J1640 as an A7 Vp SrEu(Cr) star through analysis of classification resolution spectra.

We describe how a simple class of out of equilibrium, rotating and asymmetrical mass distributions evolve under their self-gravity to produce a quasi-planar spiral structure surrounding a virialized core, qualitatively resembling a spiral galaxy. The spiral structure is transient, but can survive tens of dynamical times, and further reproduces qualitatively noted features of spiral galaxies as the predominance of trailing two-armed spirals and large pitch angles. As our models are highly idealized, a detailed comparison with observations is not appropriate, but generic features of the velocity distributions can be identified to be potential observational signatures of such a mechanism. Indeed, the mechanism leads generically to a characteristic transition from predominantly rotational motion, in a region outside the core, to radial ballistic motion in the outermost parts. Such radial motions are excluded in our Galaxy up to 15 kpc, but could be detected at larger scales in the future by GAIA. We explore the apparent motions seen by external observers of the velocity distributions of our toy galaxies, and find that it is difficult to distinguish them from those of a rotating disc with sub-dominant radial motions at levels typically inferred from observations. These simple models illustrate the possibility that the observed apparent motions of spiral galaxies might be explained by non-trivial non-stationary mass and velocity distributions without invoking a dark matter halo or modification of Newtonian gravity. In this scenario the observed phenomenological relation between the centripetal and gravitational acceleration of the visible baryonic mass could have a simple explanation.

The determination of the resolution of cosmological N-body simulations, i.e., the range of scales in which quantities measured in them represent accurately the continuum limit, is an important open question. We address it here using scale-free models, for which self-similarity provides a powerful tool to control resolution. Such models also provide a robust testing ground for the so-called stable clustering approximation, which gives simple predictions for them. Studying large N-body simulations of such models with different force smoothing, we find that these two issues are in fact very closely related: our conclusion is that the accuracy of two point statistics in the non-linear regime starts to degrade strongly around the scale at which their behaviour deviates from that predicted by the stable clustering hypothesis. Physically the association of the two scales is in fact simple to understand: stable clustering fails to be a good approximation when there are strong interactions of structures (in particular merging) and it is precisely such non-linear processes which are sensitive to fluctuations at the smaller scales affected by discretisation. Resolution may be further degraded if the short distance gravitational smoothing scale is larger than the scale to which stable clustering can propagate. We examine in detail the very different conclusions of studies by Smith et al. (2003) and Widrow et al. (2009} and find that the strong deviations from stable clustering reported by these works are the results of over-optimistic assumptions about scales resolved accurately by the measured power spectra, and the reliance on Fourier space analysis. We emphasise the much poorer resolution obtained with the power spectrum compared to the two point correlation function..

Simulations of purely self-gravitating N-body systems are often used in astrophysics and cosmology to study the collisionless limit of such systems. Their results for macroscopic quantities should then converge well for sufficiently large N. Using a study of the evolution from a simple space of spherical initial conditions - including a region characterised by so-called "radial orbit instability" - we illustrate that the values of N at which such convergence is obtained can vary enormously. In the family of initial conditions we study, good convergence can be obtained up to a few dynamical times with N $ \sim 10^3$ - just large enough to suppress two body relaxation - for certain initial conditions, while in other cases such convergence is not attained at this time even in our largest simulations with N $\sim 10^5$. The qualitative difference is due to the stability properties of fluctuations introduced by the N-body discretisation, of which the initial amplitude depends on N. We discuss briefly why the crucial role which such fluctuations can potentially play in the evolution of the N-body system could, in particular, constitute a serious problem in cosmological simulations of dark matter.

We explore the formation and relaxation of so-called quasi-stationary states (QSS) for particle distributions in three dimensions interacting via an attractive radial pair potential $V(r \rightarrow \infty) \sim 1/r^\gamma$ with $\gamma > 0$, and either a soft-core or hard-core regularization at small $r$. In the first part of the paper we generalize, for any spatial dimension $d \geq 2$, Chandrasekhar's approach for the case of gravity to obtain analytic estimates of the rate of collisional relaxation due to two body collisions. The resultant relaxation rates indicate an essential qualitative difference depending on the integrability of the pair force at large distances: for $\gamma >d-1$ the rate diverges in the large particle number $N$ (mean field) limit, unless a sufficiently large soft core is present; for $\gamma < d-1$, on the other hand, the rate vanishes in the same limit even in the absence of any regularization. In the second part of the paper we compare our analytical predictions with the results of extensive parallel numerical simulations in $d=3$, for a range of different exponents $\gamma$ and soft cores leading to the formation of QSS. We find, just as for the previously well studied case of gravity (which we also revisit), excellent agreement between the parametric dependence of the observed relaxation times and our analytic predictions. Further, as in the case of gravity, we find that the results indicate that, when large impact factors dominate, the appropriate cut-off is the size of the system (rather than, for example, the mean inter-particle distance). Our results provide strong evidence that the existence of QSS is robust only for long-range interactions with a large distance behavior $\gamma < d-1$; for $\gamma \geq d-1$ the existence of such states will be conditioned strongly on the short range properties of the interaction.

Isolated long-range interacting particle systems appear generically to relax to non-equilibrium states ("quasi-stationary states" or QSS) which are stationary in the thermodynamic limit. A fundamental open question concerns the "robustness" of these states when the system is not isolated. In this paper we explore, using both analytical and numerical approaches to a paradigmatic one dimensional model, the effect of a simple class of perturbations. We call them "internal local perturbations" in that the particle energies are perturbed at collisions in a way which depends only on the local properties. Our central finding is that the effect of the perturbations is to drive all the very different QSS we consider towards a unique QSS. The latter is thus independent of the initial conditions of the system, but determined instead by both the long-range forces and the details of the perturbations applied. Thus in the presence of such a perturbation the long-range system evolves to a unique non-equilibrium stationary state, completely different to its state in absence of the perturbation, and it remains in this state when the perturbation is removed. We argue that this result may be generic for long-range interacting systems subject to perturbations which are dependent on the local properties (e.g. spatial density or velocity distribution) of the system itself.

During the violent relaxation of a self-gravitating system a significant fraction of its mass may be ejected. If the time varying gravitational field also breaks spherical symmetry this mass can potentially carry angular momentum. Thus starting initial configurations with zero angular momentum can in principle lead to a bound virialized system with non-zero angular momentum. We explore here, using numerical simulations, how much angular momentum can be generated in a virialized structure in this way, starting from configurations of cold particles which are very close to spherically symmetric. For initial configurations in which spherical symmetry is broken only by the Poissonian fluctuations associated with the finite particle number $N$, with $N$ in range $10^3$ to $10^5$, we find that the relaxed structures have standard "spin" parameters $\lambda \sim 10^{-3}$, and decreasing slowly with $N$. For slightly ellipsoidal initial conditions, in which the finite-$N$ fluctuations break the residual reflection symmetries, we observe values $\lambda \sim 10^{-2}$, of the same order of magnitude as those reported for elliptical galaxies. The net angular momentum vector is typically aligned close to normal to the major semi-axis of the triaxial relaxed structure, and also with that of the ejected mass. This simple mechanism may provide an alternative, or complement, to "tidal torque theory" for understanding the origin of angular momentum in astrophysical structures.

Synthetic RGBB magnitudes are generated with the most recent theoretical stellar evolution models computed with the Dartmouth Stellar Evolution Program (DSEP) code. They are compared to the observational work of Nataf et al., who present RGBB magnitudes for 72 globular clusters. A DSEP model using a chemical composition with enhanced $\alpha$ capture [$\alpha$/Fe] $ =+0.4$ and an age of 13 Gyr shows agreement with observations over metallicities ranging from [Fe/H] = $0$ to [Fe/H] $\approx-1.5$, with discrepancy emerging at lower metallicities.

Initially cold and spherically symmetric self-gravitating systems may give rise to a virial equilibrium state which is far from spherically symmetric, and typically triaxial. We focus here on how the degree of symmetry breaking in the final state depends on the initial density profile. We note that the most asymmetric structures result when, during the collapse phase, there is a strong injection of energy preferentially into the particles which are localized initially in the outer shells. These particles are still collapsing when the others, initially located in the inner part, are already re-expanding; the motion of particles in a time varying potential allow them to gain kinetic energy --- in some cases enough to be ejected from the system. We show that this mechanism of energy gain amplifies the initial small deviations from perfect spherical symmetry due to finite $N$ fluctuations. This amplification is more efficient when the initial density profile depends on radius, because particles have a greater spread of fall times compared to a uniform density profile, for which very close to symmetric final states are obtained}. These effects lead to a distinctive correlation of the orientation of the final structure with the distribution of ejected mass, and also with the initial (very small) angular fluctuations.

We study non-linear gravitational clustering from cold gaussian power-law initial conditions in a family of scale-free EdS models, characterized by a free parameter $\kappa$ fixing the ratio between the mass driving the expansion and the mass which clusters. As in the "usual" EdS model, corresponding to $\kappa=1$, self-similarity provides a powerful instrument to delimit the physically relevant clustering resolved by a simulation. Likewise, if stable clustering applies, it implies scale-free non-linear clustering. We derive the corresponding exponent $\gamma_{sc} (n, \kappa)$ of the two point correlation function. We then report the results of extensive N-body simulations, of comparable size to those previously reported in the literature for the case $\kappa=1$, and performed with an appropriate modification of the GADGET2 code. We observe in all cases self-similarity in the two point correlations, down to a lower cut-off which decreases monotonically in time in comoving coordinates. The self-similar part of the non-linear correlation function is fitted well in all cases by a single power-law with an exponent in good agreement with $\gamma_{sc} (n, \kappa)$. Our results thus indicate that stable clustering provides an excellent approximation to the non-linear correlation function over the resolved self-similar scales, at least down to $\gamma_{sc} (n, \kappa) \approx 1$, corresponding to the case $n=-2$ for $\kappa=1$. We conclude, in contrast notably with the results of Smith et al. (2003), that a clear identification of the breakdown of stable clustering in self-similar models - and the possible existence of a "universal" region in which non-linear clustering becomes independent of initial conditions - remains an important open problem, which should be addressed further in significantly larger simulations.

Dark matter halos are built from accretion and merging. During merging some of the dark matter particles may be ejected with velocities higher than the escape velocity. We use both N-body simulations and single-particle smooth-field simulations to demonstrate that rapid changes to the mean field potential are responsible for such ejection, and in particular that dynamical friction plays no significant role in it. Studying a range of minor mergers, we find that typically between 5-15% of the particles from the smaller of the two merging structures are ejected. We also find that the ejected particles originate essentially from the small halo, and more specifically are particles in the small halo which pass later through the region in which the merging occurs.

Hamiltonian systems with long-range interactions give rise to long lived out of equilibrium macroscopic states, so-called quasi-stationary states. We show here that, in a suitably generalized form, this result remains valid for many such systems in the presence of dissipation. Using an appropriate mean-field kinetic description, we show that models with dissipation due to a viscous damping or due to inelastic collisions admit "scaling quasi-stationary states", i.e., states which are quasi-stationary in rescaled variables. A numerical study of one dimensional self-gravitating systems confirms both the relevance of these solutions, and gives indications of their regime of validity in line with theoretical predictions. We underline that the velocity distributions never show any tendency to evolve towards a Maxwell-Boltzmann form.

One dimensional versions of cosmological N-body simulations have been shown to share many qualitative behaviours of the three dimensional problem. They can resolve a large range of time and length scales, and admit exact numerical integration. We use such models to study how non-linear clustering depends on initial conditions and cosmology. More specifically, we consider a family of models which, like the 3D EdS model, lead for power-law initial conditions to self-similar clustering characterized in the strongly non-linear regime by power-law behaviour of the two point correlation function. We study how the corresponding exponent \gamma depends on the initial conditions, characterized by the exponent n of the power spectrum of initial fluctuations, and on a single parameter \kappa controlling the rate of expansion. The space of initial conditions/cosmology divides very clearly into two parts: (1) a region in which \gamma depends strongly on both n and \kappa and where it agrees very well with a simple generalisation of the so-called stable clustering hypothesis in three dimensions, and (2) a region in which \gamma is more or less independent of both the spectrum and the expansion of the universe. We explain the observed location of the boundary in (n, \kappa) space dividing the "stable clustering" region from the "universal" region. We compare and contrast our findings to results in three dimensions, and discuss in particular the light they may throw on the question of "universality" of non-linear clustering in this context.

Beyond convergence studies and comparison of different codes, there are essentially no controls on the accuracy in the non-linear regime of cosmological N body simulations, even in the dissipationless limit. We propose and explore here a simple test which has not been previously employed: when cosmological codes are used to simulate an isolated overdensity, they should reproduce, in physical coordinates, those obtained in open boundary conditions without expansion. In particular, the desired collisionless nature of the simulations can be probed by testing for stability in physical coordinates of virialized equilibria. We investigate and illustrate the test using a suite of simulations in an Einstein de Sitter cosmology from initial conditions which rapidly settle to virial equilibrium. We find that the criterion of stable clustering allows one to determine, for given particle number N in the "halo" and force smoothing, a maximum red-shift range over which the collisionless limit may be represented with desired accuracy. We also compare our results to the so-called Layzer Irvine test, showing that it provides a weaker, but very useful, tool to constrain the choice of numerical parameters. Finally we outline in some detail how these methods could be employed to test the choice of the numerical parameters used in a cosmological simulation.

We study quasi-stationary states (QSS) resulting from violent relaxation in the one-dimensional self-gravitating "sheet model", revisiting in particular the question of the adequacy of the theory of Lynden-Bell (LB) to describe them. For "waterbag" initial conditions characterized by a single phase space density, the prediction of this theory is, in this model, a function of only one parameter, which can conveniently be chosen to be the ratio of the energy to that in the degenerate limit. Studying a class of such initial conditions in which the shape of the initial waterbag is varied, we find that the LB predictions are reasonably good always in the low energy region, while at higher energies (i.e. in the non-degenerate limit) they are generally not even qualitatively correct, although certain initial conditions can still be found where they are as good as at low energy. We find notably that, in line with what has been observed by Levin et al. in some other models, when LB theory does not work the QSS are always characterized by the presence of a {\it degenerate} core, which these authors explain as the result of dynamical resonances. In short LB theory appears to be a good approximation only when violent relaxation is sufficiently "gentle", and otherwise a degenerate core-halo structure results.

Studies of a class of infinite one dimensional self-gravitating systems have highlighted that, on the one hand, the spatial clustering which develops may have scale invariant (fractal) properties, and, on the other, that they display "self-similar" properties in their temporal evolution. The relevance of these results to three dimensional cosmological simulations has remained unclear. We show here that the measured exponents characterizing the scale-invariant non-linear clustering are in excellent agreement with those derived from an appropriately generalized "stable-clustering" hypothesis. Further an analysis in terms of "halos" selected with a friend-of-friend algorithm reveals that such structures are, statistically, virialized across the range of scales corresponding to scale-invariance. Thus the strongly non-linear clustering in these models is accurately described as a virialized fractal structure, very much in line with the "clustering hierarchy" which Peebles originally envisaged qualitatively as associated with stable clustering. If transposed to three dimensions these results would imply, notably, that cold dark matter halos (or even subhalos) are 1) not well modeled as smooth objects, and 2) that the supposed "universality" of their profiles is, like apparent smoothness, an artefact of poor numerical resolution.

Borexino, a liquid scintillator detector at LNGS, is designed for the detection of neutrinos and antineutrinos from the Sun, supernovae, nuclear reactors, and the Earth. The feeble nature of these signals requires a strong suppression of backgrounds below a few MeV. Very low intrinsic radiogenic contamination of all detector components needs to be accompanied by the efficient identification of muons and of muon-induced backgrounds. Muons produce unstable nuclei by spallation processes along their trajectory through the detector whose decays can mimic the expected signals; for isotopes with half-lives longer than a few seconds, the dead time induced by a muon-related veto becomes unacceptably long, unless its application can be restricted to a sub-volume along the muon track. Consequently, not only the identification of muons with very high efficiency but also a precise reconstruction of their tracks is of primary importance for the physics program of the experiment. The Borexino inner detector is surrounded by an outer water-Cherenkov detector that plays a fundamental role in accomplishing this task. The detector design principles and their implementation are described. The strategies adopted to identify muons are reviewed and their efficiency is evaluated. The overall muon veto efficiency is found to be 99.992% or better. Ad-hoc track reconstruction algorithms developed are presented. Their performance is tested against muon events of known direction such as those from the CNGS neutrino beam, test tracks available from a dedicated External Muon Tracker and cosmic muons whose angular distribution reflects the local overburden profile. The achieved angular resolution is 3-5 deg and the lateral resolution is 35-50 cm, depending on the impact parameter of the crossing muon. The methods implemented to efficiently tag cosmogenic neutrons are also presented.

We report results of a study of the Newtonian dynamics of N self-gravitating particles which start in a quasi-uniform spherical configuration, without initial velocities. These initial conditions would lead to a density singularity at the origin at a finite time when N \rightarrow \infty, but this singularity is regulated at any finite N (by the associated density fluctuations). While previous studies have focussed on the behaviour as a function of N of the minimal size reached during the contracting phase, we examine in particular the size and energy of the virialized halo which results. We find the unexpected result that the structure decreases in size as N increases, scaling in proportion to N^{-1/3}, a behaviour which is associated with an ejection of kinetic energy during violent relaxation which grows in proportion to N^{1/3}. This latter scaling may be qualitatively understood, and if it represents the asymptotic behaviour in N implies that this ejected energy is unbounded above. We discuss also tests we have performed which indicate that this ejection is a mean-field phenomenon (i.e. a result of collisionless dynamics).

We formalize a classification of pair interactions based on the convergence properties of the {\it forces} acting on particles as a function of system size. We do so by considering the behavior of the probability distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that the size of the system is taken to infinity at constant particle density, i.e., in the "usual" thermodynamic limit. For a pair interaction potential V(r) with V(r) \rightarrow \infty) \sim 1/r^a defining a {\it bounded} pair force, we show that P(F) converges continuously to a well-defined and rapidly decreasing PDF if and only if the {\it pair force} is absolutely integrable, i.e., for a > d-1, where d is the spatial dimension. We refer to this case as {\it dynamically short-range}, because the dominant contribution to the force on a typical particle in this limit arises from particles in a finite neighborhood around it. For the {\it dynamically long-range} case, i.e., a \leq d-1, on the other hand, the dominant contribution to the force comes from the mean field due to the bulk, which becomes undefined in this limit. We discuss also how, for a \leq d-1 (and notably, for the case of gravity, a=d-2) P(F) may, in some cases, be defined in a weaker sense. This involves a regularization of the force summation which is generalization of the procedure employed to define gravitational forces in an infinite static homogeneous universe. We explain that the relevant classification in this context is, however, that which divides pair forces with a > d-2 (or a < d-2), for which the PDF of the {\it difference in forces} is defined (or not defined) in the infinite system limit, without any regularization. In the former case dynamics can, as for the (marginal) case of gravity, be defined consistently in an infinite uniform system.

We revisit the issue of relaxation to thermal equilibrium in the so-called "sheet model", i.e., particles in one dimension interacting by attractive forces independent of their separation. We show that this relaxation may be very clearly detected and characterized by following the evolution of order parameters defined by appropriately normalized moments of the phase space distribution which probe its entanglement in space and velocity coordinates. For a class of quasi-stationary states which result from the violent relaxation of rectangular waterbag initial conditions, characterized by their virial ratio R_0, we show that relaxation occurs on a time scale which (i) scales approximately linearly in the particle number N, and (ii) shows also a strong dependence on R_0, with quasi-stationary states from colder initial conditions relaxing much more rapidly. The temporal evolution of the order parameter may be well described by a stretched exponential function. We study finally the correlation of the relaxation times with the amplitude of fluctuations in the relaxing quasi-stationary states, as well as the relation between temporal and ensemble averages.

"Quasi-stationary" states are approximately time-independent out of equilibrium states which have been observed in a variety of systems of particles interacting by long-range interactions. We investigate here the conditions of their occurrence for a generic pair interaction V(r \rightarrow \infty) \sim 1/r^a with a > 0, in d>1 dimensions. We generalize analytic calculations known for gravity in d=3 to determine the scaling parametric dependences of their relaxation rates due to two body collisions, and report extensive numerical simulations testing their validity. Our results lead to the conclusion that, for a < d-1, the existence of quasi-stationary states is ensured by the large distance behavior of the interaction alone, while for a > d-1 it is conditioned on the short distance properties of the interaction, requiring the presence of a sufficiently large soft-core in the interaction potential.

The dynamics of infinite, asymptotically uniform, distributions of self-gravitating particles in one spatial dimension provides a simple toy model for the analogous three dimensional problem. We focus here on a limitation of such models as treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (i.e. the definition requires explicitly the breaking of statistical translational invariance). The problem arises because naive background subtraction (due to expansion, or by "Jeans' swindle" for the static case), applied as in three dimensions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a discussion by Kiessling, we show that the problem may be resolved by defining the force in infinite point distributions as the limit of an exponentially screened pair interaction. We show that this prescription gives a well defined (finite) force acting on particles in a class of perturbed infinite lattices, which are the point processes relevant to cosmological N-body simulations. For identical particles the dynamics of the simplest toy model is equivalent to that of an infinite set of points with inverted harmonic oscillator potentials which bounce elastically when they collide. We discuss previous results in the literature, and present new results for the specific case of this simplest (static) model starting from "shuffled lattice" initial conditions. These show qualitative properties (notably its "self-similarity") of the evolution very similar to those in the analogous simulations in three dimensions, which in turn resemble those in the expanding universe.

When an open system of classical point particles interacting by Newtonian gravity collapses and relaxes violently, an arbitrary amount of energy may in principle be carried away by particles which escape to infinity. We investigate here, using numerical simulations, how this released energy and other related quantities (notably the binding energy and size of the virialized structure) depends on the initial conditions, for the one parameter family of starting configurations given by randomly distributing N cold particles in a spherical volume. Previous studies have established that the minimal size reached by the system scales approximately as N^{-1/3}, a behaviour which follows trivially when the growth of perturbations (which regularize the singularity of the cold collapse in the infinite N limit) are assumed to be unaffected by the boundaries. Our study shows that the energy ejected grows approximately in proportion to N^{1/3}, while the fraction of the initial mass ejected grows only very slowly with N, approximately logarithmically, in the range of N simulated. We examine in detail the mechanism of this mass and energy ejection, showing explicitly that it arises from the interplay of the growth of perturbations with the finite size of the system. A net lag of particles compared to their uniform spherical collapse trajectories develops first at the boundaries and then propagates into the volume during the collapse. Particles in the outer shells are then ejected as they scatter through the time dependent potential of an already re-expanding central core. Using modified initial configurations we explore the importance of fluctuations at different scales, and discreteness (i.e. non-Vlasov) effects in the dynamics.

The effects of discreteness arising from the use of the N-body method on the accuracy of simulations of cosmological structure formation are not currently well understood. After a discussion of how the relevant discretisation parameters introduced should be extrapolated to recover the Vlasov-Poisson limit, we study numerically, and with analytical methods we have developed recently, the central issue of how finite particle density affects the precision of results. In particular we focus on the power spectrum at wavenumbers around and above the Nyquist wavenumber, in simulations in which the force resolution is taken smaller than the initial interparticle spacing. Using simulations of identical theoretical initial conditions sampled on four different "pre-initial" configurations (three different Bravais lattices, and a glass) we obtain a {\it lower bound} on the real discreteness error. With the guidance of our analytical results, we establish with confidence that the measured dispersion is not contaminated either by finite box size effects or by subtle numerical effects. Our results show notably that, at wavenumbers {\it below} the Nyquist wavenumber, the dispersion increases monotonically in time throughout the simulation, while the same is true above the Nyquist wavenumber once non-linearity sets in. For normalizations typical of cosmological simulations, we find lower bounds on errors at the Nyquist wavenumber of order of a percent, and larger above this scale. The only way this error may be reduced below these levels at these scales, and indeed convergence to the physical limit firmly established, is by extrapolation, at fixed values of the other relevant parameters, to the regime in which the mean comoving interparticle distance becomes less than the force smoothing scale.

In cosmology numerical simulations of structure formation are now of central importance, as they are the sole instrument for providing detailed predictions of current cosmological models for a whole class of important constraining observations. These simulations are essentially molecular dynamics simulations of N (>> 1), now up to of order several billion) particles interacting through their self-gravity. While their aim is to produce the Vlasov limit, which describes the underlying (``cold dark matter'') models, the degree to which they actually do produce this limit is currently understood, at best, only very qualitatively, and there is an acknowledged need for ``a theory of discreteness errors''. In this talk I will describe, for non-cosmologists, both the simulations and the underlying theoretical models, and will then focus on the issue of discreteness, describing some recent progress in addressing this question quantitatively.

In two recent articles a detailed study has been presented of the out of equilibrium dynamics of an infinite system of self-gravitating points initially located on a randomly perturbed lattice. In this article we extend the treatment of the early time phase during which strong non-linear correlations first develop, prior to the onset of ``self-similar'' scaling in the two point correlation function. We establish more directly, using appropriate modifications of the numerical integration, that the development of these correlations can be well described by an approximation of the evolution in two phases: a first perturbative phase in which particles' displacements are small compared to the lattice spacing, and a subsequent phase in which particles interact only with their nearest neighbor. For the range of initial amplitudes considered we show that the first phase can be well approximated as a transformation of the perturbed lattice configuration into a Poisson distribution at the relevant scales. This appears to explain the ``universality'' of the spatial dependence of the asymptotic non-linear clustering observed from both shuffled lattice and Poisson initial conditions.

The usual thermodynamic limit for systems of classical self-gravitating point particles becomes well defined, as a {\it dynamical} problem, using a simple physical prescription for the calculation of the force, equivalent to the so-called ``Jeans' swindle''. The relation of the resulting intrinsically out of equilibrium problem, of particles evolving from prescribed uniform initial conditions in an infinite space, to the one studied in current cosmological models (in an expanding universe) is explained. We then describe results of a numerical study of the dynamical evolution of such a system, starting from a simple class of infinite ``shuffled lattice'' initial conditions. The clustering, which develops in time starting from scales around the grid scale, is qualitatively very similar to that seen in cosmological simulations, which begin from lattices with applied correlated displacements and incorporate an expanding spatial background. From very soon after the formation of the first non-linear structures, a spatio-temporal scaling relation describes well the evolution of the two-point correlations. At larger times the dynamics of these correlations converges to what is termed ``self-similar'' evolution in cosmology, in which the time dependence in the scaling relation is specified entirely by that of the linearized fluid theory. We show how this statistical mechanical ``toy model'' can be useful in addressing various questions about these systems which are relevant in cosmology. Some of these questions are closely analagous to those currently studied in the literature on long range interactions, notably the relation of the evolution of the particle system to that in the Vlasov limit and the nature of approximately quasi-stationary states.

We consider the construction of point processes from tilings, with equal volume tiles, of d-dimensional Euclidean space. We show that one can generate, with simple algorithms ascribing one or more points to each tile, point processes which are "superhomogeneous'' (or "hyperuniform''), i.e., for which the structure factor S(k) vanishes when the wavenumber k tends to zero. The exponent of the leading small-k behavior depends in a simple manner on the nature of the correlation properties of the specific tiling and on the conservation of the mass moments of the tiles. Assigning one point to the center of mass of each tile gives the exponent \gamma=4 for any tiling in which the shapes and orientations of the tiles are short-range correlated. Smaller exponents, in the range 4-d<\gamma<4 (and thus always superhomogeneous for d\leq 4), may be obtained in the case that the latter quantities have long-range correlations. Assigning more than one point to each tile in an appropriate way, we show that one can obtain arbitrarily higher exponents in both cases. We illustrate our results with explicit constructions using known deterministic tilings, as well as some simple stochastic tilings for which we can calculate S(k) exactly. Our results provide, we believe, the first explicit analytical construction of point processes with \gamma > 4. Applications to condensed matter physics, and also to cosmology, are briefly discussed.

We study how the two-point density correlation properties of a point particle distribution are modified when each particle is divided, by a stochastic process, into an equal number of identical "daughter" particles. We consider generically that there may be non-trivial correlations in the displacement fields describing the positions of the different daughters of the same "mother" particle, and then treat separately the cases in which there are, or are not, correlations also between the displacements of daughters belonging to different mothers. For both cases exact formulae are derived relating the structure factor (power spectrum) of the daughter distribution to that of the mother. These results can be considered as a generalization of the analogous equations obtained in ref. [1] (cond-mat/0409594) for the case of stochastic displacement fields applied to particle distributions. An application of the present results is that they give explicit algorithms for generating, starting from regular lattice arrays, stochastic particle distributions with an arbitrarily high degree of large-scale uniformity.

We apply a recently developed perturbative formalism which describes the evolution under their self-gravity of particles displaced from a perfect lattice to quantify precisely, up to shell crossing, the effects of discreteness in dissipationless cosmological N-body simulations. We give simple expressions, explicitly dependent on the particle density, for the evolution of power in each mode as a function of red-shift. For typical starting red-shifts the effect of finite particle number is to {\it slow down} slightly the growth of power compared to that in the fluid limit (e.g. by about ten percent at half the Nyquist frequency), and to induce also dispersion in the growth as a function of direction at a comparable level. In the limit that the initial red-shift tends to infinity, at fixed particle density, the evolution in fact diverges from that in the fluid limit (described by the Zeldovich approximation). Contrary to widely held belief, this means that a simulation started at a red-shift much higher than the red-shift of shell crossing actually gives a worse, rather than a better, result. We also study how these effects are modified when there is a small-scale regularization of the gravitational force. We show that such a smoothing may reduce the anisotropy of the discreteness effects, but it then {\it increases} their average effect. This behaviour illustrates the fact that the discreteness effects described here are distinct from those usually considered in this context, due to two-body collisions. Indeed the characteristic time for divergence from the collisionless limit is proportional to $N^{2/3}$, rather than $N/ \log N$ in the latter case.

The relation between the results of cosmological N-body simulations, and the continuum theoretical models they simulate, is currently not understood in a way which allows a quantification of N dependent effects. In this first of a series of papers on this issue, we consider the quantification of such effects in the initial conditions of such simulations. A general formalism developed in [1] allows us to write down an exact expression for the power spectrum of the point distributions generated by the standard algorithm for generating such initial conditions. Expanded perturbatively in the amplitude of the input (i.e. theoretical, continuum) power spectrum, we obtain at linear order the input power spectrum, plus two terms which arise from discreteness and contribute at large wavenumbers. For cosmological type power spectra, one obtains as expected, the input spectrum for wavenumbers k smaller than that characteristic of the discreteness. The comparison of real space correlation properties is more subtle because the discreteness corrections are not as strongly localised in real space. For cosmological type spectra the theoretical mass variance in spheres and two point correlation function are well approximated {\it above} a finite distance. For typical initial amplitudes this distance is a few times the inter-particle distance, but it diverges as this amplitude (or, equivalently, the initial red-shift of the cosmological simulation) goes to zero, at fixed particle density. We discuss briefly the physical significance of these discreteness terms in the initial conditions, in particular with respect to the definition of the continuum limit of N-body simulations.

We study, using numerical simulations, the dynamical evolution of self-gravitating point particles in static euclidean space, starting from a simple class of infinite ``shuffled lattice'' initial conditions. These are obtained by applying independently to each particle on an infinite perfect lattice a small random displacement, and are characterized by a power spectrum (structure factor) of density fluctuations which is quadratic in the wave number k, at small k. For a specified form of the probability distribution function of the ``shuffling'' applied to each particle, and zero initial velocities, these initial configurations are characterized by a single relevant parameter: the variance $\delta^2$ of the ``shuffling'' normalized in units of the lattice spacing $\ell$. The clustering, which develops in time starting from scales around $\ell$, is qualitatively very similar to that seen in cosmological simulations, which begin from lattices with applied correlated displacements and incorporate an expanding spatial background. From very soon after the formation of the first non-linear structures, a spatio-temporal scaling relation describes well the evolution of the two-point correlations. At larger times the dynamics of these correlations converges to what is termed ``self-similar'' evolution in cosmology, in which the time dependence in the scaling relation is specified entirely by that of the linearized fluid theory. Comparing simulations with different $\delta$, different resolution, but identical large scale fluctuations, we are able to identify and study features of the dynamics of the system in the transient phase leading to this behavior. In this phase, the discrete nature of the system explicitly plays an essential role.

We study the evolution under their self-gravity of infinite ``shuffled lattice'' particle distributions, focussing specifically on the comparison of this evolution with that of ``daughter'' particle distributions, defined by a simple coarse-graining procedure. We consider both the case that such coarse-grainings are performed (i) on the initial conditions, and (ii) at a finite time with a specific additional prescription. In numerical simulations we observe that, to a first approximation, these coarse-grainings represent well the evolution of the two-point correlation properties over a significant range of scales. We note, in particular, that the form of the two-point correlation function in the original system, when it is evolving in the asymptotic ``self-similar'' regime, may be reproduced well in a daughter coarse-grained system in which the dynamics are still dominated by two-body (nearest neighbor) interactions. Using analytical results on the early time evolution of these systems, however, we show that small observed differences between the evolved system and its coarse-grainings at the initial time will in fact diverge as the ratio of the coarse-graining scale to the original inter-particle distance increases. The second coarse-graining studied, performed at a finite time in a specified manner, circumvents this problem. It also makes more physically transparent why gravitational dynamics from these initial conditions tends toward a ``self-similar'' evolution. We finally discuss the precise definition of a limit in which a continuum (specifically Vlasov-like) description of the observed linear and non-linear evolution should be applicable.

N-body simulations sample their initial conditions on an initial particle distribution, which for cosmological simulations is usually a glass or grid, whilst a Poisson distribution is used for galaxy models, spherical collapse etc. These pre-initial conditions have inherent correlations, noise due to discreteness and preferential alignments, whilst the glass distribution is poorly defined and computationally expensive to construct. We present a novel particle distribution which can be useful as a pre-initial condition for N-body simulations, using a simple construction based on a ``quaquaversal'' tiling of space. This distribution has little preferred orientation (i.e. is statistically isotropic), has a rapidly vanishing large scale power-spectrum (P(k) ~ k^4), and is trivial to create. It should be particularly useful for warm dark matter and cold collapse simulations.

We study the statistics of the force felt by a particle in the class of spatially correlated distribution of identical point-like particles, interacting via a $1/r^2$ pair force (i.e. gravitational or Coulomb), and obtained by randomly perturbing an infinite perfect lattice. In the first part we specify the conditions under which the force on a particle is a well defined stochastic quantity. We then study the small displacements approximation, giving both the limitations of its validity, and, when it is valid, an expression for the force variance. In the second part of the paper we extend to this class of particle distributions the method introduced by Chandrasekhar to study the force probability density function in the homogeneous Poisson particle distribution. In this way we can derive an approximate expression for the probability distribution of the force over the full range of perturbations of the lattice, i.e., from very small (compared to the lattice spacing) to very large where the Poisson limit is recovered. We show in particular the qualitative change in the large-force tail of the force distribution between these two limits. Excellent accuracy of our analytic results is found on detailed comparison with results from numerical simulations. These results provide basic statistical information about the fluctuations of the interactions (i) of the masses in self-gravitating systems like those encountered in the context of cosmological N-body simulations, and (ii) of the charges in the ordered phase of the One Component Plasma.

We present a perturbative treatment of the evolution under their mutual self-gravity of particles displaced off an infinite perfect lattice, both for a static space and for a homogeneously expanding space as in cosmological N-body simulations. The treatment, analogous to that of perturbations to a crystal in solid state physics, can be seen as a discrete (i.e. particle) generalization of the perturbative solution in the Lagrangian formalism of a self-gravitating fluid. Working to linear order, we show explicitly that this fluid evolution is recovered in the limit that the initial perturbations are restricted to modes of wavelength much larger than the lattice spacing. The full spectrum of eigenvalues of the simple cubic lattice contains both oscillatory modes and unstable modes which grow slightly faster than in the fluid limit. A detailed comparison of our perturbative treatment, at linear order, with full numerical simulations is presented, for two very different classes of initial perturbation spectra. We find that the range of validity is similar to that of the perturbative fluid approximation (i.e. up to close to ``shell-crossing''), but that the accuracy in tracing the evolution is superior. The formalism provides a powerful tool to systematically calculate discreteness effects at early times in cosmological N-body simulations.

We discuss some of the basic implications of recent results on galaxy correlations published by the SDSS collaboration. In particular we focus on the evidence which has been recently presented for the scale and nature of the transition to homogeneity in the galaxy distribution, and results which describe the dependence of clustering on luminosity. The two questions are in fact strictly entangled, as the stability of the measure of the amplitude of the correlation function depends on the scale at which the mean density becomes well defined. We note that the recent results which indicate the convergence to well defined homogeneity in a volume equivalent to that of a sphere of radius 70 Mpc/h, place in doubt previous detections of ``luminosity bias'' from measures of the amplitude of the correlation function. We emphasize that the way to resolve these issues is to first use, in volume limited samples corresponding to different ranges of luminosity, the unnormalized two point statistics to establish the scale (and value) at which the mean density becomes well defined. We note also that the recent SDSS results for these statistics are in good agreement with those obtained by us through analyses of many previous samples, confirming in particular that the galaxy distribution is well described by a fractal dimension D ~ 2 up to a scale of at least 20 Mpc/h. We discuss critically the agreement of this new data with current theoretical models.

We investigate the possibility of generating initial conditions for cosmological N-body simulations by simulating a system whose correlations at thermal equilibrium approximate well those of cosmological density perturbations. The system is an appropriately modified version of the standard ``one component plasma'' (OCP). We show first how a well-known semi-analytic method can be used to determine the potential required to produce the desired correlations, and then verify our results for some cosmological type spectra with simulations of the full molecular dynamics. The advantage of the method, compared to the standard one, is that it gives by construction an accurate representation of both the real and reciprocal space correlation properties of the theoretical model. Furthermore the distributions are also statistically homogeneous and isotropic. We discuss briefly the modifications needed to implement the method to produce configurations appropriate for large N-body simulations in cosmology, and also the generation of initial velocities in this context.

We apply a simple linearization, well known in solid state physics, to approximate the evolution at early times of cosmological N-body simulations of gravity. In the limit that the initial perturbations, applied to an infinite perfect lattice, are at wavelengths much greater than the lattice spacing $l$ the evolution is exactly that of a pressureless self-gravitating fluid treated in the analagous (Lagrangian) linearization, with the Zeldovich approximation as a sub-class of asymptotic solutions. Our less restricted approximation allows one to trace the evolution of the discrete distribution until the time when particles approach one another (i.e. ``shell crossing''). We calculate modifications of the fluid evolution, explicitly dependent on $l$ i.e. discreteness effects in the N body simulations. We note that these effects become increasingly important as the initial red-shift is increased at fixed $l$. The possible advantages of using a body centred cubic, rather than simple cubic, lattice are pointed out.

We present an analysis of different sets of gravitational N-body simulations, all describing the dynamics of discrete particles with a small initial velocity dispersion. They encompass very different initial particle configurations, different numerical algorithms for the computation of the force, with or without the space expansion of cosmological models. Despite these differences we find in all cases that the non-linear clustering which results is essentially the same, with a well-defined simple power-law behaviour in the two-point correlations in the range from a few times the lower cut-off in the gravitational force to the scale at which fluctuations are of order one. We argue, presenting quantitative evidence, that this apparently universal behaviour can be understood by the domination of the small scale contribution to the gravitational force, coming initially from nearest neighbor particles.

A well known argument in cosmology gives that the power spectrum (or structure function) $P(k)$ of mass density fluctuations produced from a uniform initial state by physics which is causal (i.e. moves matter and momentum only up to a finite scale) has the behaviour $P(k) \propto k^4$ at small $k$. Noting the assumption of analyticity at $k=0$ of $P(k)$ in the standard derivation of this result, we introduce a class of solvable one dimensional models which allows us to study the relation between the behaviour of $P(k)$ at small $k$ and the properties of the probability distribution $f(l)$ for the spatial extent $l$ of mass and momentum conserving fluctuations. We find that the $k^4$ behaviour is obtained in the case that the first {\it six} moments of $f(l)$ are finite. Interestingly the condition that the fluctuations be localised - taken to correspond to the convergence of the first two moments of $f(l)$ - imposes only the weaker constraint $P(k) \propto k^n$ with $n$ anywhere in the range $0< n \leq 4$. We interpret this result to suggest that the causality bound will be loosened in this way if quantum fluctuations are permitted.

Threshold biasing of a Gaussian random field gives a linear amplification of the reduced two point correlation function at large distances. We show that for standard cosmological models this does not translate into a linear amplification of the power spectrum (PS) at small k. For standard CDM type models this means that the``turn-over'' at small k of the original PS disappears in the PS of the biased field for the physically relevant range of the threshold parameter. In real space this difference is manifest in the asymptotic behaviour of the normalised mass variance in spheres of radius R, which changes from the ``super-homogeneous'' behaviour to a Poisson-like behaviour.This qualitative change results from the intrinsic stochasticity of the threshold sampling. While our quantitative results are specific to the simplest threshold biasing model, we argue that our qualitative conclusions should be valid generically for any biasing mechanism involving a scale-dependent amplification of the correlation function. One implication is that the real-space correlation function will be a better instrument to probe for the underlying Harrison Zeldovich spectrum in the distribution of visible matter, as the characteristic asymptotic negative power-law \xi (r) \sim -r^{-4} tail is undistorted by biasing.

The initial conditions describing seed fluctuations for the formation of structure in standard cosmological models, i.e.the Harrison-Zeldovich distribution, have very characteristic ``super-homogeneous'' properties: they are statistically translation invariant, isotropic, and the variance of the mass fluctuations in a region of volume V grows slower than V. We discuss the geometrical construction of distributions of points in ${\bf R}^3$ with similar properties encountered in tiling and in statistical physics, e.g. the Gibbs distribution of a one-component system of charged particles in a uniform background (OCP). Modifications of the OCP can produce equilibrium correlations of the kind assumed in the cosmological context. We then describe how such systems can be used for the generation of initial conditions in gravitational $N$-body simulations.

We analyse with simple real-space statistics the Virgo consortium's cosmological N-body simulations. Significant clustering rapidly develops well below the initial mean interparticle separation \Lambda_i, where the gravitational force on a particle is dominated by that with its nearest neighbours. A power-law behaviour in the two point correlation function emerges, which in the subsequent evolution is continuously amplified and shifted to larger scales, in a roughly self-similar manner. We conclude that the density fluctuations at the smallest scales due to the particle-like nature of the distribution being evolved are thus essential in the development of these correlations, and not solely, as usually supposed, the very small continuous (fluid-like) fluctuations at scales larger than \Lambda_i >.

In this lecture we address three different but related aspects of the initial continuous fluctuation field in standard cosmological models. Firstly we discuss the properties of the so-called Harrison-Zeldovich like spectra. This power spectrum is a fundamental feature of all current standard cosmological models. In a simple classification of all stationary stochastic processes into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with $P(0)=0$, belong. In statistical physics language they are well described as glass-like. Secondly, the initial continuous density field with such small amplitude correlated Gaussian fluctuations must be discretised in order to set up the initial particle distribution used in gravitational N-body simulations. We discuss the main issues related to the effects of discretisation, particularly concerning the effect of particle induced fluctuations on the statistical properties of the initial conditions and on the dynamical evolution of gravitational clustering.

After reviewing some basic relevant properties of stationary stochastic processes (SSP), we discuss the properties of the so-called Harrison-Zeldovich like spectra of mass density perturbations. These correlations are a fundamental feature of all current standard cosmological models. Examining them in real space we note they imply a "sub-poissonian" normalised variance in spheres $\sigma_M^2(R) \sim R^{-4} \ln R$. In particular this latter behaviour is at the limit of the most rapid decay ($\sim R^{-4}$) of this quantity possible for any stochastic distribution (continuous or discrete). In a simple classification of all SSP into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with $P(0)=0$, belong. In statistical physics language they are well described as lattice or glass-like. We illustrate their properties through two simple examples: (i) the ``shuffled'' lattice and the One Component Plasma at thermal equilibrium.

After reviewing the basic relevant properties of stationary stochastic processes (SSP), defining basic terms and quantities, we discuss the properties of the so-called Harrison-Zeldovich like spectra. These correlations, usually characterized exclusively in k-space (i.e. in terms of power spectra P(k)), are a fundamental feature of all current standard cosmological models. Examining them in real space we note their characteristics to be a {\it negative} power law tail \xi(r) \sim - r^{-4} and a {\it sub-poissonian} normalised variance in spheres \sigma^2(R) \sim R^{-4} \ln R. We note in particular that this latter behaviour is at the limit of the most rapid decay (\sim R^{-4}) of this quantity possible for any stochastic distribution (continuous or discrete). This very particular characteristic is usually obscured in cosmology by the use of Gaussian spheres. In a simple classification of all SSP into three categories, we highlight with the name ``super-homogeneous'' the properties of the class to which models like this, with P(0)=0, belong. In statistical physics language they are well described as glass-like. They do not have either ``scale-invariant'' features, in the sense of critical phenomena, nor fractal properties. We illustrate their properties with some simple examples, in particular that of a ``shuffled'' lattice.

In most direct estimates of the mass density (visible or dark) of the Universe, a central input parameter is the luminosity density of the Universe. Here we consider the measurement of this luminosity density from red-shift surveys, as a function of the yet undetermined characteristic scale R_H at which the spatial distribution of visible matter tends to a well defined homogeneity. Making the canonical assumption that the cluster mass to luminosity ratio M/L is the universal one, we can estimate the total mass density as a function \Omega_m(R_H,M/L). Taking the highest estimated cluster value M/L ~300h and a conservative lower limit R_H > 20 Mpc/h, we obtain the upper bound \Omega_m < 0.1 . We note that for values of the homogeneity scale R_H in the range R_H ~ (90 +/- 45) hMpc, the value of \Omega_m may be compatible with the nucleosynthesis inferred density in baryons.

Scalar fields can play a dominant role in the dynamics of the Universe until shortly before nucleosynthesis. Examples are provided by domination by a kinetic mode of a scalar field, which may be both the inflaton and the late time `quintessence', and also by more conventional models of reheating. The resultant modification to the pre-nucleosynthesis expansion rate can allow solely an asymmetry in right handed electrons to produce a net baryon asymmetry when reprocessed by the anomalous B+L violating processes of the standard model. The production of such a source asymmetry - what we term `electrogenesis' - requires no additional B or L violation beyond that in the standard model. We consider a specific model for its generation, by a simple perturbative out of equilibrium decay of Higgs like scalar fields with CP-violating Yukawa couplings to the standard model leptons. We show that, because of the much enhanced expansion rate, such a mechanism can easily produce an adequate asymmetry from scalars with masses as low as 1 TeV. Kinetic mode domination is strongly favoured because it evades large entropy release which dilutes the asymmetry. We also discuss briefly the effect of the abelian hypercharge anomaly.

The clustering of galaxies is well characterized by fractal properties, with the presence of an eventual cross-over to homogeneity still a matter of considerable debate. In this letter we discuss the cosmological implications of a fractal distribution of matter, with a possible cross-over to homogeneity at an undetermined scale R_{homo}. Contrary to what is generally assumed, we show that, even when R_{homo} -> \infty, this possibility can be treated consistently within the framework of the expanding universe solutions of Friedmann. The fractal is a perturbation to an open cosmology in which the leading homogeneous component is the cosmic background radiation (CBR). This cosmology, inspired by the observed galaxy distributions, provides a simple explanation for the recent data which indicate the absence of deceleration in the expansion (q_o \approx 0). Correspondingly the `age problem' is also resolved. Further we show that the model can be extended back from the curvature dominated arbitrarily deep into the radiation dominated era, and we discuss qualitatively the modifications to the physics of the anisotropy of the CBR, nucleosynthesis and structure formation.

We study the fermion propagator in a spatially varying classical background field, and show that, contrary to common wisdom, it may get nontrivial gradient corrections already at the first order in derivative expansion. This occurs whenever the fermion self-energy acquires a spatially (or temporally) varying pseudoscalar term, a simple example of which is given by a complex mass term m(x)= m_R + i gamma_5 m_I. Such effective mass terms arise for example in extensions of the Standard Model during the electroweak transition, and they are crucial in providing the CP-violation necessary for electroweak baryogenesis.

In a recent analysis of number counts in the ESP survey Scaramella et al. (1998) claim to find evidence for a cross-over to homogeneity at large scales, and against a fractal behaviour with dimension $D \approx 2$. In this comment we note firstly that, if such a cross-over exists as described by the authors, the scale characterizing it is ~ 100 - 300 Mpc/h. This invalidates the ``standard'' analysis of the same catalogue given elsewhere by the authors which results in a ``correlation length'' of only r_0 = 4 Mpc/h. Furthermore we show that the evidences for a cross-over to homogeneity rely on the choice of cosmological model, and most crucially on the so called K corrections. We show that the D ~ 3 behaviour seen in the K-corrected data of Scaramella et al. is in fact unstable, increasing systematically towards D=4 as a function of the absolute magnitude limit. This behaviour can be quantitatively explained as the effect of an unphysical K-correction in the relevant range of red-shift (z ~ 0.1- 0.3). A more consistent interpretation of the number counts is that D is in the range 2 - 2.5, depending on the cosmological model, consistent with the continuation of the fractal D ~ 2 behaviour observed at scales up to ~100 Mpc/h. This implies a smaller K-correction. Given, however, the uncertainty in the effect of intrinsic fluctuations on the number counts statistic, and its sensitivity on these large scales to the uncertain K corrections, we conclude that it is premature to put a definitive constraint on the galaxy distribution using the ESP data alone.

We report our analysis of the properties of galaxy clustering for a new redshift sample of galaxies, the CfA2-South catalog, using statistical methods which do not rely on the assumption of homogeneity. We find that, up to ~ 20 Mpc/h, which is the largest scale to which correlation properties can be reliably inferred, the galaxy clustering is scale-invariant and characterized by a fractal dimension D=1.9 \pm 0.1. Further there is no statistical evidence for homogeneity at any of the larger scales (up to ~150 Mpc/h) probed more weakly by the catalog. These results means that characteristic ``correlation lengths'' for the clustering of galaxies derived using standards methods of analysis are not meaningful. Further the results are very consistent with those obtained from many other catalogs using the methods adopted here, which show the D =2 fractal continuing to beyond 100 Mpc/h. The incompleteness of the relevant data conjectured by various authors to give rise to such behaviour is therefore proved to have no significant effect (up to 20 Mpc/h) on the measured correlations.

A weakly coupled scalar field $\Phi$ with a simple exponential potential $V=M_P^4\exp(-\lambda\Phi/M_P)$ where $M_P$ is the reduced Planck mass, and $\lambda > 2$, has an attractor solution in a radiation or matter dominated universe in which it mimics the scaling of the dominant component, contributing a fixed fraction $\Omega_\phi$ (determined by $\lambda$) to the energy density. Such fields arise generically in particle physics theories involving compactified dimensions, with values of $\lambda$ which give a cosmologically relevant $\Omega_\phi$. For natural initial conditions on the scalar field in the early universe the attractor solution is established long before the epoch of structure formation, and in contrast to the solutions used in other scalar field cosmologies, it is one which does not involve an energy scale for the scalar field characteristic of late times . We study in some detail the evolution of matter and radiation perturbations in a standard inflation-motivated $\Omega=1$ dark-matter dominated cosmology with this extra field. Using a full Einstein-Boltzmann calculation we compare observable quantities with current data. We find that, for $\Omega_\phi\simeq 0.08-0.12$, these models are consistent with large angle cosmic microwave background anisotropies as detected by COBE, the linear mass variance as compiled from galaxy surveys, big bang nucleosynthesis, the abundance of rich clusters and constraints from the Lyman-$\alpha$ systems at high redshift. Given the simplicity of the model, its theoretical motivation and its success in matching observations, we argue that it should be taken on a par with other currently viable models of structure formation.

We calculate the baryon asymmetry generated at the electroweak phase transition in the minimal supersymmetric standard model, treating the particles in a WKB approximation in the bubble wall background. A set of diffusion equations for the particle species relevant to baryon generation, including source terms arising from the CP violation associated with the complex phase $\delta$ of the $\mu$ parameter, are derived from Boltzmann equations, and solved. The conclusion is that $\delta$ must be $\gsim 0.1$ to generate a baryon asymmetry consistent with nucleosynthesis. We compare our results to several other recent computations of the effect, arguing that some are overestimates.

The usual sphaleron bound and the statement of the impossibility of baryon production at a second order phase transition or analytic cross-over are reformulated in the first part of the paper as requirements of the expansion rate of the Universe at the electroweak scale. With an (exact or effective) additional contribution to the energy density scaling as 1/a^6, which dominates until just before nucleosynthesis, the observed baryon asymmetry may be produced at the electroweak scale in simple extensions of the Minimal Standard Model, even in the case that the phase transition is not first order. We focus our attention on one such cosmology, in which the Universe goes through a period termed `kination' in which its energy is dominated by the kinetic energy of a scalar field. The required kinetic energy dominated modes can occur either as a field rolls down an exponential (or steeper) potential, or in the oscillation of a field about the minimum of a steep power-law potential. We implement in detail the former case with a single exponential field first driving inflation, and then rolling into a kinetic energy dominated mode. Reheating is achieved using an alternative to the usual mechanism due to Spokoiny, in which the Universe is `reheated' by particle creation in the expanding background. Density perturbations of the magnitude required for structure formation may also be generated. We show that the analogous model for the power-law potential cannot be consistently implemented. In models with inflation driven by a second field and the usual mechanism of reheating (by decay of the inflaton) the required kinetic energy dominated cosmology is viable in both types of potential.

Radiation domination at the electroweak epoch is a simplifying assumption, but one for which there is no observational basis. Treating the expansion rate as a variable, I re-examine electroweak baryogenesis in various scenarios. At a first order phase transition the main effect is on the sphaleron bound, which becomes a lower bound on the expansion rate in any given theory. At a second-order or cross-over phase transition, the created baryon asymmetry is directly proportional to the expansion rate. I sketch an alternative post-inflationary cosmology, in which the kinetic energy of a scalar field dominates the Universe until shortly before nucleosynthesis, and argue that the observed baryon asymmetry could be produced in this case even at an analytic cross-over.

A scalar field with an exponential potential has the particular property that it is attracted into a solution in which its energy scales as the dominant component (radiation or matter) of the Universe, contributing a fixed fraction of the total energy density. We study the growth of perturbations in a CDM dominated $\Omega=1$ universe with this extra field, with an initial flat spectrum of adiabatic fluctuations. The observational constraints from structure formation are satisfied as well, or better, than in other models, with a contribution to the energy density from the scalar field $\Omega_\phi \sim 0.1$ which is small enough to be consistent with entry into the attractor prior to nucleosynthesis.

In the standard model there are charges with abelian anomaly only (e.g. right-handed electron number) which are effectively conserved in the early universe until some time shortly before the electroweak scale. A state at finite chemical potential of such a charge, possibly arising due to asymmetries produced at the GUT scale, is unstable to the generation of hypercharge magnetic field. Quite large magnetic fields ($\sim 10^{22}$ gauss at $T\sim 100$ GeV with typical inhomogeneity scale $ \sim \frac{ 10^6}{T}$) can be generated. These fields may be of cosmological interest, potentially acting as seeds for amplification to larger scale magnetic fields through non-linear mechanisms. Previously derived bounds on exotic $B-L$ violating operators may also be evaded.

The standard requirement for the production of baryons at the electroweak phase transition, that the phase transition be first order and the sphaleron bound be satisfied, is predicated on the assumption of a radiation dominated universe at that epoch. One simple alternative - domination by the energy in a kinetic mode of a scalar field which scales as $1/a^6$ - gives a significantly weakened sphaleron bound for the preservation of a baryon asymmetry produced at a first-order phase transition, and allows the possibility that the observed baryon asymmetry be produced when the phase transition is second-order or cross-over. Such a phase of `kination' at the electroweak scale can occur in various ways as a scalar field evolves in an exponential potential after inflation.

We investigate baryogenesis at a first order electroweak phase transition in the presence of a CP violating condensate on the bubble walls, in the regime in which the bubble walls are `thick', in the sense that fermions interact with the plasma many times as the bubble wall passes. Such a condensate is present in multi-Higgs extensions of the standard model and may be formed via an instability in the minimal standard model. We concentrate on particles with typical thermal energies in the plasma, whose interactions with the wall are accurately described by the WKB approximation, in which a classical chiral force is evident. The deviations from thermal equilibrium produced by motion of the wall are then treated using a classical Boltzmann equation which we solve in a fluid approximation. From the resulting equations we find two effects important for baryogenesis: (i) a classical chiral force term due to the $CP$ violating background, and (ii) a term arising from hypercharge violating interactions which are pushed out of equilibrium by the background{field. Provided the wall propagates slower than the speed of sound, both terms lead to the diffusion of a chiral asymmetry in front of the wall. This can produce a baryon asymmetry of the observed magnitude for typical wall velocities and thicknesses.

We investigate `non-local' schemes for baryogenesis at a first order electroweak phase transition, in which the effects of a $CP$ violating condensate on the bubble wall propagate into the unbroken phase where the sphaleron rate is unsupressed. Such a condensate exists in multi-Higgs extensions of the standard model, and may exist due to an instability in the minimal standard model. In this paper we first discuss the general problem of determining the $CP$ violating perturbations, distinguishing two regimes (quantum and classical). We then give an analytic treatment of quantum mechanical reflection in the thin wall regime, in which interactions with the plasma can be neglected as a particle propagates across the wall. We focus on leptons because of their much weaker coupling to the plasma. We argue that they are likely to be accurately described by this calculation, but quarks are not. The relative magnitude of the baryon asymmetry produced for different fermions depends on their relative Yukawa couplings ({\it not} their zero temperature masses), their transport properties and their interactions. We calculate the baryon asymmetry for various parameter ranges and conclude that asymmetries comparable with observations can be generated.

The inflationary cosmology is analyzed from the point of view of squeezed quantum states. As noted by Grishchuk and Sidorov, the amplification of quantum fluctuations into macroscopic perturbations which occurs during cosmic inflation is a process of quantum squeezing. We carefully develop the squeezed state formalism and derive the equations that govern the evolution of a gaussian initial state. We derive the power spectrum of density perturbations for a simple inflationary model and discuss its features. We conclude that the squeezed state formalism provides an interesting framework within which to study the amplification process, but,in disagreement with the claims of Grishchuk and Sidorov, that it does {\em not} provide us with any new physical results.