BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael Joyce

Michael Joyce appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COcond-mat.stat-mechastro-phhep-phmath.AGmath.COastro-ph.GAcond-matgr-qchep-thmath.GRmath.RTphysics.ed-ph

Catalog footprint

What is connected

34works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

Accuracy of power spectra in dissipationless cosmological simulations

We exploit a suite of large \emph{N}-body simulations (up to N=$4096^3$) performed with \Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At the 1\% level it starts at early times from a fraction of $k_Λ$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, $\sim 2-3 k_Λ$ at the very latest times we evolve to. At the $5\%$ level, accuracy extends up to wavenumbers of order $5k_Λ$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple $n$-dependence, allowing a extrapolation to place conservative bounds on the accuracy of \emph{N}-body simulations of non-scale free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_Λ$ may need revision in light of our convergence analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantifying resolution in cosmological N-body simulations using self-similarity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Schubert polynomial analogues for degenerate involutions

We survey the recent study of involution Schubert polynomials and a modest generalization that we call degenerate involution Schubert polynomials. We cite several conditions when (degenerate) involution Schubert polynomials have simple factorization formulae. Such polynomials can be computed by traversing through chains in certain weak order posets, and we provide explicit descriptions of such chains in weak order for involutions and degenerate involutions. As an application, we give several examples of how certain multiplicity-free sums of Schubert polynomials factor completely into very simple linear factors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Attractor non-equilibrium stationary states in perturbed long-range interacting systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite $N$ corrections to Vlasov dynamics and the range of pair interactions

We explore the conditions on a pair interaction for the validity of the Vlasov equation to describe the dynamics of an interacting $N$ particle system in the large $N$ limit. Using a coarse-graining in phase space of the exact Klimontovich equation for the $N$ particle system, we evaluate, neglecting correlations of density fluctuations, the scalings with $N$ of the terms describing the corrections to the Vlasov equation for the coarse-grained one particle phase space density. Considering a generic interaction with radial pair force $F(r)$, with $F(r) \sim 1/r^γ$ at large scales, and regulated to a bounded behaviour below a "softening" scale $\varepsilon$, we find that there is an essential qualitative difference between the cases $γ< d$ and $γ> d$, i.e., depending on the integrability at large distances of the pair force. In the former case the corrections to the Vlasov dynamics for a given coarse-grained scale are essentially insensitive to the softening parameter $\varepsilon$, while for $γ> d$ the amplitude of these terms is directly regulated by $\varepsilon$, and thus by the small scale properties of the interaction. This corresponds to a simple physical criterion for a basic distinction between long-range ($γ\leq d $) and short range ($γ> d$) interactions, different to the canonical one ($γ\leq d +1$ or $γ> d +1$ ) based on thermodynamic analysis. This alternative classification, based on purely dynamical considerations, is relevant notably to understanding the conditions for the existence of so-called quasi-stationary states in long-range interacting systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the generation of triaxiality in the collapse of cold spherical self-gravitating systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chains in Weak Order Posets Associated to Involutions

The W-set of an element of a weak order poset is useful in the cohomological study of the closures of spherical subgroups in generalized flag varieties. We explicitly describe in a purely combinatorial manner the W-sets of the weak order posets of three different sets of involutions in the symmetric group, namely, the set of all involutions, the set of all fixed point free involutions, and the set of all involutions with signed fixed points (or "clans"). These distinguished sets of involutions parameterize Borel orbits in the classical symmetric spaces associated to the general linear group. In particular, we give a complete characterization of the maximal chains of an arbitrary lower order ideal in any of these three posets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Particle ejection during mergers of dark matter halos

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling quasi-stationary states in long range systems with dissipation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-similarity and stable clustering in a family of scale-free cosmologies

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 $κ$ fixing the ratio between the mass driving the expansion and the mass which clusters. As in the "usual" EdS model, corresponding to $κ=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 $γ_{sc} (n, κ)$ 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 $κ=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 $γ_{sc} (n, κ)$. 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 $γ_{sc} (n, κ) \approx 1$, corresponding to the case $n=-2$ for $κ=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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exponents of non-linear clustering in scale-free one dimensional cosmological simulations

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 γdepends on the initial conditions, characterized by the exponent n of the power spectrum of initial fluctuations, and on a single parameter κcontrolling the rate of expansion. The space of initial conditions/cosmology divides very clearly into two parts: (1) a region in which γdepends strongly on both n and κ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 γ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, κ) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Introduction of interactive learning into French university physics classrooms

We report on a project to introduce interactive learning strategies (ILS) to physics classes at the Université Pierre et Marie Curie (UPMC), one of the leading science universities in France. In Spring 2012, instructors in two large introductory classes, first-year, second-semester mechanics, and second-year introductory E&M, enrolling approximately 500 and 250 students respectively, introduced ILS into some sections of each class. The specific ILS utilized were Think-Pair-Share questions and Peer Instruction in the main lecture classrooms, and UW Tutorials for Introductory Physics in recitation sections. Pre- and post-instruction assessments (FCI and CSEM respectively) were given, along with a series of demographics questions. We were able to compare the results of the FCI and CSEM between interactive and non-interactive classes taught simultaneously with the same curriculum. We also analyzed final exam results, as well as the results of student and instructor attitude surveys between classes. In our analysis, we argue that Multiple Linear Regression modeling is superior to other common analysis tools, including normalized gain. Our results show that ILS are effective at improving student learning by all measures used: research-validated concept inventories and final exam scores, on both conceptual and traditional problem-solving questions. Multiple Linear Regression analysis reveals that interactivity in the classroom is a significant predictor of student learning, showing a similar or stronger relationship with student learning than such ascribed characteristics as parents' education, and achieved characteristics such as GPA and hours studied per week. Analysis of student and instructors attitudes shows that both groups believe that ILS improve student learning in the physics classroom, and increases student engagement and motivation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of isolated overdensities as a control on cosmological N body simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ordered Bell numbers, Hermite polynomials, Skew Young Tableaux, and Borel orbits

We give three interpretations of the number $b$ of orbits of the Borel subgroup of upper triangular matrices on the variety $\ms{X}$ of complete quadrics. First, we show that $b$ is equal to the number of standard Young tableaux on skew-diagrams. Then, we relate $b$ to certain values of a modified Hermite polynomial. Third, we relate $b$ to a certain cell decomposition on $\ms{X}$ previously studied by De Concini, Springer, and Strickland. Using these, we give asymptotic estimates for $b$ as the dimension of the quadrics increases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-linear gravitational clustering of cold matter in an expanding universe: indications from 1D toy models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quasi-stationary states in the self-gravitating sheet model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A dynamical classification of the range of pair interactions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational force in an infinite one-dimensional Poisson distribution

We consider the statistical properties of the gravitational field F in an infinite one-dimensional homogeneous Poisson distribution of particles, using an exponential cut-off of the pair interaction to control and study the divergences which arise. Deriving an exact analytic expression for the probability density function (PDF) P(F), we show that it is badly defined in the limit in which the well known Holtzmark distribution is obtained in the analogous three-dimensional case. A well defined P(F) may, however, be obtained in the infinite range limit by an appropriate renormalization of the coupling strength, giving a Gaussian form. Calculating the spatial correlation properties we show that this latter procedure has a trivial physical meaning. Finally we calculate the PDF and correlation properties of differences of forces (at separate spatial points), which are well defined without any renormalization. We explain that the convergence of these quantities is in fact sufficient to allow a physically meaningful infinite system limit to be defined for the clustering dynamics from Poissonian initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quasi-stationary states and the range of pair interactions

"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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relaxation to thermal equilibrium in the self-gravitating sheet model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

1-d gravity in infinite point distributions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Universality of power law correlations in gravitational clustering

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.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The Quantum Boltzmann Equation in a Nontrivial Background

This talk is a status report on our study of quantum transport equations relevant for baryogenesis computations. Our main finding is that, as a consequence of localization in space, the quasiparticle picture of the plasma dynamics breaks down at first non-trivial order in gradient expansion. While in this talk we focus on bosons, we expect that a similar picture holds for fermions. We then argue that the quasiparticle picture is recovered in the adiabatic limit of frequent scattering.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Cosmology with a Primordial Scaling Field

A weakly coupled scalar field $Φ$ with a simple exponential potential $V=M_P^4\exp(-λΦ/M_P)$ where $M_P$ is the reduced Planck mass, and $λ> 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 $Ω_ϕ$ (determined by $λ$) to the energy density. Such fields arise generically in particle physics theories involving compactified dimensions, with values of $λ$ which give a cosmologically relevant $Ω_ϕ$. 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 $Ω=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 $Ω_ϕ\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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Turning Around the Sphaleron Bound: Electroweak Baryogenesis in an Alternative Post-inflationary Cosmology

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.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Electroweak Baryogenesis and the Expansion Rate of the Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Electroweak Baryogenesis from a Classical Force

We describe a new effect which produces baryons at a first order electroweak phase transition. It operates when there is a CP-violating field present on propagating bubble walls. The novel aspect is that it involves a purely classical force, which alters the motion of particles across the wall and through diffusion creates a chiral asymmetry in front of the wall. We develop a technique for computing the baryon asymmetry using the Boltzmann equation, and a fluid approximation which allows us to model strong scattering effects. The final formula for the baryon asymmetry has a remarkably simple form.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Efficient Electroweak Baryogenesis from Lepton Transport

One mechanism for generating a baryon asymmetry at the electroweak phase transition involves propagation of particle asymmetries generated by reflection from the bubble walls into the unbroken phase. Hitherto attention has focussed on top quarks because of their large mass and thus effective scattering from bubble walls. In this paper we point out that leptons may be more efficient mediators of this type of electroweak baryogenesis, particularly for thicker bubble walls favored by perturbative calculations. We carry out an analytic calculation of each stage of the mechanism and conclude it produces a baryon asymmetry compatible with observations for a wide range of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Inflation and squeezed quantum states

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.

0 citations0 reviewsNo code linkedOpen access