Source author record

Edmund Bertschinger

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2016arXiv

Women in physics in the United States: recruitment and retention

Initiatives to increase the number, persistence, and success of women in physics in the US reach pre-teen girls through senior women. Programs exist at both the local and national levels. In addition, researchers have investigated issues related to gender equity in physics and physics education. Anecdotal evidence suggests increased media coverage of the underrepresentation of women in science. All of these efforts are both motivated and made more effective by the collection and presentation of data on the presence, persistence, and promise of women in physics.

preprint2015arXiv

An Extension of the Faddeev-Jackiw Technique to Fields in Curved Spacetimes

The Legendre transformation on singular Lagrangians, e.g. Lagrangians representing gauge theories, fails due to the presence of constraints. The Faddeev-Jackiw technique, which offers an alternative to that of Dirac, is a symplectic approach to calculating a Hamiltonian paired with a well-defined initial value problem when working with a singular Lagrangian. This phase space coordinate reduction was generalized by Barcelos-Neto and Wotzasek to simplify its application. We present an extension of the Faddeev-Jackiw technique for constraint reduction in gauge field theories and non-gauge field theories that are coupled to a curved spacetime that is described by General Relativity. A major difference from previous formulations is that we do not explicitly construct the symplectic matrix, as that is not necessary. We find that the technique is a useful tool that avoids some of the subtle complications of the Dirac approach to constraints. We apply this formulation to the Ginzburg-Landau action and provide a calculation of its Hamiltonian and Poisson brackets in a curved spacetime.

preprint2012arXiv

Scale-Dependent Growth from a Transition in Dark Energy Dynamics

We investigate the observational consequences of the quintessence field rolling to and oscillating near a minimum in its potential, "if" it happens close to the present epoch (z<0.2). We show that in a class of models, the oscillations lead to a rapid growth of the field fluctuations and the gravitational potential on subhorizon scales. The growth in the gravitational potential occurs on timescales << H^(1). This effect is present even when the quintessence parameters are chosen to reproduce an expansion history consistent with observations. For linearized fluctuations, we find that although the gravitational potential power spectrum is enhanced in a scale-dependent manner, the shape of the dark matter/galaxy power spectrum is not significantly affected. We find that the best constraints on such a transition in the quintessence field is provided via the integrated Sachs-Wolfe (ISW) effect in the CMB temperature power spectrum. Going beyond the linearized regime, the quintessence field can fragment into large, localized, long lived excitations (oscillons) with sizes comparable to galaxy clusters; this fragmentation could provide additional observational constraints. Two quoted "signatures" of modified gravity are a scale-dependent growth of the gravitational potential and a difference between the matter power spectrum inferred from measurements of lensing and galaxy clustering. Here, both effects are achieved by a minimally coupled scalar field in general relativity with a canonical kinetic term.

preprint2011arXiv

One Gravitational Potential or Two? Forecasts and Tests

The metric of a perturbed Robertson-Walker spacetime is characterized by three functions: a scale-factor giving the expansion history and two potentials which generalize the single potential of Newtonian gravity. The Newtonian potential induces peculiar velocities and, from these, the growth of matter fluctuations. Massless particles respond equally to the Newtonian potential and to a curvature potential. The difference of the two potentials, called the gravitational slip, is predicted to be very small in general relativity but can be substantial in modified gravity theories. The two potentials can be measured, and gravity tested on cosmological scales, by combining weak gravitational lensing or the Integrated Sachs-Wolfe effect with galaxy peculiar velocities or clustering.

preprint2010arXiv

Seen and unseen tidal caustics in the Andromeda galaxy

Indirect detection of high-energy particles from dark matter interactions is a promising avenue for learning more about dark matter, but is hampered by the frequent coincidence of high-energy astrophysical sources of such particles with putative high-density regions of dark matter. We calculate the boost factor and gamma-ray flux from dark matter associated with two shell-like caustics of luminous tidal debris recently discovered around the Andromeda galaxy, under the assumption that dark matter is its own supersymmetric antiparticle. These shell features could be a good candidate for indirect detection of dark matter via gamma rays because they are located far from the primary confusion sources at the galaxy's center, and because the shapes of the shells indicate that most of the mass has piled up near apocenter. Using a numerical estimator specifically calibrated to estimate densities in N-body representations with sharp features and a previously determined N-body model of the shells, we find that the largest boost factors do occur in the shells but are only a few percent. We also find that the gamma-ray flux is an order of magnitude too low to be detected with Fermi for likely dark matter parameters, and about 2 orders of magnitude less than the signal that would have come from the dwarf galaxy that produces the shells in the N-body model. We further show that the radial density profiles and relative radial spacing of the shells, in either dark or luminous matter, is relatively insensitive to the details of the potential of the host galaxy but depends in a predictable way on the velocity dispersion of the progenitor galaxy.

preprint2010arXiv

Self-Similar Spherical Collapse with Tidal Torque

N-body simulations have revealed a wealth of information about dark matter halos however their results are largely empirical. Using analytic means, we attempt to shed light on simulation results by generalizing the self-similar secondary infall model to include tidal torque. In this first of two papers, we describe our halo formation model and compare our results to empirical mass profiles inspired by N-body simulations. Each halo is determined by four parameters. One parameter sets the mass scale and the other three define how particles within a mass shell are torqued throughout evolution. We choose torque parameters motivated by tidal torque theory and N-body simulations and analytically calculate the structure of the halo in different radial regimes. We find that angular momentum plays an important role in determining the density profile at small radii. For cosmological initial conditions, the density profile on small scales is set by the time rate of change of the angular momentum of particles as well as the halo mass. On intermediate scales, however, $ρ\propto r^{-2}$, while $ρ\propto r^{-3}$ close to the virial radius.

preprint2010arXiv

Velocity Structure of Self-Similar Spherically Collapsed Halos

Using a generalized self-similar secondary infall model, which accounts for tidal torques acting on the halo, we analyze the velocity profiles of halos in order to gain intuition for N-body simulation results. We analytically calculate the asymptotic behavior of the internal radial and tangential kinetic energy profiles in different radial regimes. We then numerically compute the velocity anisotropy and pseudo-phase-space density profiles and compare them to recent N-body simulations. For cosmological initial conditions, we find both numerically and analytically that the anisotropy profile asymptotes at small radii to a constant set by model parameters. It rises on intermediate scales as the velocity dispersion becomes more radially dominated and then drops off at radii larger than the virial radius where the radial velocity dispersion vanishes in our model. The pseudo-phase-space density is universal on intermediate and large scales. However, its asymptotic slope on small scales depends on the halo mass and on how mass shells are torqued after turnaround. The results largely confirm N-body simulations but show some differences that are likely due to our assumption of a one-dimensional phase space manifold.

preprint2009arXiv

Hierarchy in the Phase Space and Dark Matter Astronomy

We develop a theoretical framework for describing the hierarchical structure of the phase space of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo model, the stable clustering hypothesis is applied for the first time here to the small-scale phase space structure. As an application, we show that the particle dark matter annihilation signal could be up to two orders of magnitude larger than that of the smooth halo within the Galactic virial radius. The local boost is inversely proportional to the smooth halo density, and thus is O(1) within the solar radius, which could translate into interesting signatures for dark matter direct detection experiments: The temporal correlation of dark matter detection can change by a factor of 2 in the span of 10 years, while there will be significant correlations in the velocity space of dark matter particles. This can introduce O(1) uncertainty in the direction of local dark matter wind, which was believed to be a benchmark of directional dark matter searches or the annual modulation signal.

preprint2008arXiv

Hierarchical Phase Space Structure of Dark Matter Haloes: Tidal debris, Caustics, and Dark Matter annihilation

Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. The density of debris is not constant, but rather can grow due to formation of caustics at the apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. Understanding this phase structure can be of significant importance for the interpretation of many astrophysical observations and in particular dark matter detection experiments. With this purpose in mind, we develop a general theoretical framework to describe the hierarchical structure of the phase space of cold dark matter haloes. We do not make any assumption of spherical symmetry and/or smooth and continuous accretion. Instead, working with correlation functions in the action-angle space, we can fully account for the hierarchical structure (predicting a two-point correlation function ~ (ΔJ)^{-1.6} in the action space), as well as the primordial discreteness of the phase space. As an application, we estimate the boost to the dark matter annihilation signal due to the structure of the phase space within virial radius: the boost due to the hierarchical tidal debris is of order unity, whereas the primordial discreteness of the phase structure can boost the total annihilation signal by up to an order of magnitude. The latter is dominated by the regions beyond 20% of the virial radius, and is largest for the recently formed haloes with the least degree of phase mixing.

preprint2007arXiv

Variational Integrators for the Gravitational N-Body Problem

This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the algorithm to handle individual time steps; this introduces fifth-order errors in angular momentum conservation and symplecticity. We show that using adaptive block power of two timesteps does not increase the error in symplecticity. In contrast to other high-order, symplectic, individual timestep, momentum-preserving algorithms, the algorithm takes only forward timesteps. We compare a code integrating an N-body system using the algorithm with a direct-summation force calculation to standard stellar cluster simulation codes. We find that our algorithm has about 1.5 orders of magnitude better symplecticity and momentum conservation errors than standard algorithms for equivalent numbers of force evaluations and equivalent energy conservation errors.

preprint2004arXiv

A tensor formalism for transfer and Compton scattering of polarized light

A novel covariant formalism for the treatment of the transfer and Compton scattering of partially polarized light is presented. This was initially developed to aid in the computation of relativistic corrections to the polarization generated by the Sunyaev-Zeldovich effect (demonstrated in a companion paper), but it is of more general utility. In this approach, the polarization state of a light beam is described by a tensor constructed from the time average of quadratic products of the electric field components in a local observer frame. This leads naturally to a covariant description which is ideal for calculations involving the boosting of polarized light beams between Lorentz frames, and is more flexible than the traditional Stokes parameter approach in which a separate set of polarization basis vectors is required for each photon. The covariant kinetic equation for Compton scattering of partially polarized light by relativistic electrons is obtained in the tensor formalism by a heuristic semi-classical line of reasoning. The kinetic equation is derived first in the electron rest frame in the Thomson limit, and then is generalized to account for electron recoil and allow for scattering from an arbitrary distribution of electrons.

preprint1995arXiv

Constraining Peaks in Gaussian Primordial Density Fields: an Application of the Hoffman-Ribak Method

We develop an algorithm for setting up initial Gaussian random density and velocity fields containing one or more peaks or dips, in an arbitrary cosmological scenario. The intention is to generate appropriate initial conditions for cosmological N-body simulations that focus on the evolution of the progenitors of the present-day galaxies and clusters. The procedure is an application of the direct and accurate prescription of Hoffman & Ribak (1991) for generating constrained random fields. For each peak a total of 21 physical characteristics can be specified, including its scale, position, density Hessian, velocity, and velocity gradient. The velocity (or, equivalently, gravity) field constrants are based on a generalization of the formalism developed by Bardeen et al. (1986). The resulting density field is sculpted such that it induces the desired amount of net gravitational and tidal forces. We provide a detailed mathematical presentation of the formalism. Afterwards we provide analytical estimates of the likelihood of the imposed constraints. Amongst others, it is shown that the tidal field has a strong tendency to align itself along the principal axes of the mass tensor. The method is illustrated by means of some concrete examples. In addition to the illustration of constraint-field correlation functions and how they add up to the mean fields, followed by illustrations of the variance characteristics of field realizations, we concentrate in particular on the consequences of imposing gravitational field constraints (or, equivalent in the linear regime for growing mode fluctuations, peculiar velocity field constraints).

preprint1995arXiv

COSMICS: Cosmological Initial Conditions and Microwave Anisotropy Codes

COSMICS is a package of fortran programs useful for computing transfer functions and microwave background anisotropy for cosmological models, and for generating gaussian random initial conditions for nonlinear structure formation simulations of such models. Four programs are provided: {\bf linger\_con} and {\bf linger\_syn} integrate the linearized equations of general relativity, matter, and radiation in conformal Newtonian and synchronous gauge, respectively; {\bf deltat} integrates the photon transfer functions computed by the linger codes to produce photon anisotropy power spectra; and {\bf grafic} tabulates normalized matter power spectra and produces constrained or unconstrained samples of the matter density field. Version 1.0 of COSMICS is available at http://arcturus.mit.edu/cosmics/ . The current release gives fortran-77 programs that run on workstations and vectorized supercomputers. Unix makefiles are included that make it simple to build and test the package. A future release will include portable parallel versions of the linger codes using standard message-passing libraries.

preprint1995arXiv

Cosmological Perturbation Theory in the Synchronous and Conformal Newtonian Gauges

This paper presents a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. It differs from others in the literature in that we give, in both gauges, a complete discussion of all particle species that are relevant to any flat cold dark matter (CDM), hot dark matter (HDM), or CDM+HDM models (including a possible cosmological constant). The particles considered include CDM, baryons, photons, massless neutrinos, and massive neutrinos (an HDM candidate), where the CDM and baryons are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. Particular care is applied to the massive neutrino component, which has been either ignored or approximated crudely in previous works. Isentropic initial conditions on super-horizon scales are derived. The coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric and density perturbations are then solved numerically in both gauges for the standard CDM model and two CDM+HDM models with neutrino mass densities $\onu=0.2$ and 0.3, assuming a scale-invariant, adiabatic spectrum of primordial fluctuations. We also give the full details of the cosmic microwave background anisotropy, and present the first accurate calculations of the angular power spectra in the two CDM+HDM models including photon polarization, higher neutrino multipole moments, and helium recombination. The numerical programs for both gauges are available at http://arcturus.mit.edu/cosmics/ .

preprint1995arXiv

Local Approximations to the Gravitational Collapse of Cold Matter

We investigate three different local approximations for nonlinear gravitational instability in the framework of cosmological Lagrangian fluid dynamics of cold dust. They include the Zel'dovich approximation (ZA), the ``non-magnetic'' approximation of Bertschinger \& Jain (1994, NMA), and a new ``local tidal'' approximation (LTA). The LTA is exact for any perturbations whose gravitational and velocity equipotentials have the same constant shape with time, including spherical, cylindrical, and plane-parallel perturbations. We tested all three local approximations with the collapse of a homogeneous triaxial ellipsoid, for which an exact solution exists for an ellipsoid embedded in empty space and an excellent approximation is known in the cosmological context. We find that the LTA is significantly more accurate in general than the ZA and the NMA. Like the ZA, but unlike the NMA, the LTA generically leads to pancake collapse. For a randomly chosen mass element in an Einstein-de Sitter universe, assuming a Gaussian random field of initial density fluctuations, the LTA predicts that at least 78\% of initially underdense regions collapse owing to nonlinear effects of shear and tides.