Source author record

Lawrence M. Widrow

Lawrence M. Widrow 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

25works
4topics
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

25 published item(s)

preprint2022arXiv

Estimate of the Mass and Radial Profile of the Orphan-Chenab Stream's Dwarf Galaxy Progenitor Using MilkyWay@home

We fit the mass and radial profile of the Orphan-Chenab Stream's (OCS) dwarf galaxy progenitor by using turnoff stars in the Sloan Digital Sky Survey (SDSS) and the Dark Energy Camera (DEC) to constrain N-body simulations of the OCS progenitor falling into the Milky Way on the 1.5 PetaFLOPS MilkyWay@home distributed supercomputer. We infer the internal structure of the OCS's progenitor under the assumption that it was a spherically symmetric dwarf galaxy comprised of a stellar system embedded in an extended dark matter halo. We optimize the evolution time, the baryonic and dark matter scale radii, and the baryonic and dark matter masses of the progenitor using a differential evolution algorithm. The likelihood score for each set of parameters is determined by comparing the simulated tidal stream to the angular distribution of OCS stars observed in the sky. We fit the total mass of the OCS's progenitor to ($2.0\pm0.3$) $\times 10^7 M_\odot$ with a mass-to-light ratio of $γ=73.5\pm10.6$ and ($1.1\pm0.2$)$\times10^6M_{\odot}$ within 300 pc of its center. Within the progenitor's half-light radius, we estimate total a mass of ($4.0\pm1.0$)$\times10^5M_{\odot}$. We also fit the current sky position of the progenitor's remnant to be $(α,δ)=((166.0\pm0.9)^\circ,(-11.1\pm2.5)^\circ)$ and show that it is gravitationally unbound at the present time. The measured progenitor mass is on the low end of previous measurements, and if confirmed lowers the mass range of ultrafaint dwarf galaxies. Our optimization assumes a fixed Milky Way potential, OCS orbit, and radial profile for the progenitor, ignoring the impact of the Large Magellanic Cloud.

preprint2016arXiv

Galactoseismology and the Local Density of Dark Matter

We model vertical breathing mode perturbations in the Milky Way's stellar disc and study their effects on estimates of the local dark matter density, surface density, and vertical force. Evidence for these perturbations, which involve compression and expansion of the Galactic disc perpendicular to its midplane, come from the SEGUE, RAVE, and LAMOST surveys. We show that their existence may lead to systematic errors of $10\%$ or greater in the vertical force $K_z(z)$ at $|z|=1.1\,{\rm kpc}$. These errors translate to $\gtrsim 25\%$ errors in estimates of the local dark matter density. Using different mono-abundant subpopulations as tracers offers a way out: if the inferences from all tracers in the Gaia era agree, then the dark matter determination will be robust. Disagreement in the inferences from different tracers will signal the breakdown of the unperturbed model and perhaps provide the means for determining the nature of the perturbation.

preprint2016arXiv

Simulating a slow bar in the low surface brightness galaxy UGC 628

We present a disc-halo N-body model of the low surface brightness galaxy UGC 628, one of the few systems that harbours a "slow" bar with a ratio of corotation radius to bar length of $\mathcal{R} \equiv R_c/a_b \sim 2$. We select our initial conditions using SDSS DR10 photometry, a physically motivated radially variable mass-to-light ratio profile, and rotation curve data from the literature. A global bar instability grows in our submaximal disc model, and the disc morphology and dynamics agree broadly with the photometry and kinematics of UGC 628 at times between peak bar strength and the onset of buckling. Prior to bar formation, the disc and halo contribute roughly equally to the potential in the galaxy's inner region, giving the disc enough self gravity for bar modes to grow. After bar formation there is significant mass redistribution, creating a baryon dominated inner and dark matter dominated outer disc. This implies that, unlike most other low surface brightness galaxies, UGC 628 is not dark matter dominated everywhere. Our model nonetheless implies that UGC 628 falls on same the relationship between dark matter fraction and rotation velocity found for high surface brightness galaxies, and lends credence to the argument that the disc mass fraction measured at the location where its contribution to the potential peaks is not a reliable indicator of its dynamical importance at all radii.

preprint2015arXiv

Estimating the Galactic Mass Profile in the Presence of Incomplete Data

A powerful method to measure the mass profile of a galaxy is through the velocities of tracer particles distributed through its halo. Transforming this kind of data accurately to a mass profile M(r), however, is not a trivial problem. In particular, limited or incomplete data may substantially affect the analysis. In this paper we develop a Bayesian method to deal with incomplete data effectively; we have a hybrid-Gibbs sampler that treats the unknown velocity components of tracers as parameters in the model. We explore the effectiveness of our model using simulated data, and then apply our method to the Milky Way using velocity and position data from globular clusters and dwarf galaxies. We find that in general, missing velocity components have little effect on the total mass estimate. However, the results are quite sensitive to the outer globular cluster Pal 3. Using a basic Hernquist model with an isotropic velocity dispersion, we obtain credible regions for the cumulative mass profile M(r) of the Milky Way, and provide estimates for the model parameters with 95 percent Bayesian credible intervals. The mass contained within 260 kpc is 1.37x10^12 solar masses, with a 95 percent credible interval of (1.27,1.51)x10^12 solar masses. The Hernquist parameters for the total mass and scale radius are 1.55 (+0.18/-0.13)x10^12 solar masses and 16.9 (+4.8/-4.1) kpc, where the uncertainties span the 95 percent credible intervals. The code we developed for this work, Galactic Mass Estimator (GME), will be available as an open source package in the R Project for Statistical Computing.

preprint2015arXiv

Vertical Oscillations of Fluid and Stellar Disks

A satellite galaxy or dark matter subhalo that passes through a stellar disk may excite coherent oscillations in the disk perpendicular to its plane. We determine the properties of these modes for various self-gravitating plane symmetric systems (Spitzer sheets) using the matrix method of Kalnajs. In particular, we find an infinite series of modes for the case of a barotropic fluid. In general, for a collisionless system, there is a double series of modes, which include normal modes and/or Landau-damped oscillations depending on the phase space distribution function of the stars. Even Landau-damped oscillations may decay slowly enough to persist for several hundred Myr. We discuss the implications of these results for the recently discovered vertical perturbations in the kinematics of solar neighborhood stars and for broader questions surrounding secular phenomena such as spiral structure in disk galaxies.

preprint2014arXiv

Bending and Breathing Modes of the Galactic Disk

We explore the hypothesis that a passing satellite or dark matter subhalo has excited coherent oscillations of the Milky Way's stellar disk in the direction perpendicular to the Galactic midplane. This work is motivated by recent observations of spatially dependent bulk vertical motions within ~ kpc of the Sun. A satellite can transfer a fraction of its orbital energy to the disk stars as it plunges through the Galactic midplane thereby heating and thickening the disk. Bulk motions arise during the early stages of such an event when the disk is still in an unrelaxed state. We present simple toy-model calculations and simulations of disk-satellite interactions, which show that the response of the disk depends on the relative velocity of the satellite. When the component of the satellite's velocity perpendicular to the disk is small compared with that of the stars, the perturbation is predominantly a bending mode. Conversely, breathing and higher order modes are excited when the vertical velocity of the satellite is larger than that of the stars. We argue that the compression and rarefaction motions seen in three different surveys are in fact breathing mode perturbations of the Galactic disk.

preprint2013arXiv

A New Approach to Detailed Structural Decomposition from the SPLASH and PHAT Surveys: Kicked-up Disk Stars in the Andromeda Galaxy?

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 4 < R_proj < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ~20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < m_814W < 22) stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sersic, exponential, and cored power-law, respectively, and the LFs are modeled as broken power-laws. We present probability distributions for each of 32 parameters describing the SB profiles and LFs of the three subcomponents. We find that the number of stars with a disk-like LF is ~5% larger than the the number with disk-like (dynamically cold) kinematics, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a "kicked-up disk" halo population in M31. The fraction of kicked-up disk stars is consistent with that found in simulations. We also find evidence for a radially varying disk LF, consistent with a negative metallicity gradient in the stellar disk.

preprint2013arXiv

A New Approach to Detailed Structural Decomposition: Kicked-up Disk Stars in Andromeda's Halo?

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 0.4 < R < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ~20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < I < 22) stars from the Panchromatic Hubble Andromeda Treasury (PHAT) survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sersic, exponential, and cored power-law, respectively, and the LFs are modeled as broken power-laws. We find that the number of stars with a disk-like LF is ~5% larger than the number in the dynamically cold component, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a "kicked-up disk" population in the stellar halo of M31. The fraction of kicked-up disk stars is consistent with that found in simulations. See Dorman et al. (2013, submitted) for more information.

preprint2013arXiv

A Vast Thin Plane of Co-rotating Dwarf Galaxies Orbiting the Andromeda Galaxy

Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way's disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology.

preprint2013arXiv

The effect of bars on the M-sigma relation: offset, scatter and residuals correlations

We analyse a set of collisionless disc galaxy simulations to study the consequences of bar formation and evolution on the M-sigma relation of supermassive black holes. The redistribution of angular momentum driven by bars leads to a mass increase within the central region, raising the velocity dispersion of the bulge, sigma, on average by ~12% and as much as ~20%. If a disc galaxy with a SMBH satisfying the M-sigma relation forms a bar, and the SMBH does not grow in the process, then the increase in sigma moves the galaxy off the M-sigma relation. We explore various effects that can affect this result including contamination from the disc and anisotropy. The displacement from the M-sigma relation for individual model barred galaxies correlates with both M(B)/M(B+D) and beta_phi(B+D) measured within the effective radius of the bulge. Overall, this process leads to an M-sigma for barred galaxies offset from that of unbarred galaxies, as well as an increase in its scatter. We assemble samples of unbarred and barred galaxies with classical bulges and find tentative hints of an offset between the two consistent with the predicted. Including all barred galaxies, rather than just those with a classical bulge, leads to a significantly larger offset.

preprint2012arXiv

Galactoseismology: Discovery of Vertical Waves in the Galactic Disk

We present evidence for a Galactic North-South asymmetry in the number density and bulk velocity of solar neighborhood stars. The number density profile, which is derived from main-sequence stars in the Sloan Digital Sky Survey, shows a (North - South)/(North + South) deficit at |z| ~ 400 pc and an excess at |z| ~ 800 pc. The bulk velocity profile, which is derived from the Sloan Extension for Galactic Understanding and Exploration, shows a gradual trend across the Galactic midplane as well as smaller-scale features. We speculate that the North-South asymmetry, which has the appearance of a wavelike perturbation, is intrinsic to the disk. We explore the physics of this phenomenon through an analysis of the linearized Boltzmann and Poisson equations and through one-dimensional simulations. The perturbation may be excited by the passage of a satellite galaxy or dark matter subhalo through the Galactic disk, in which case we are witnessing a recent disk-heating event.

preprint2011arXiv

Magnetic fields in the Large-Scale Structure of the Universe

Magnetic fields appear to be ubiquitous in astrophysical environments. Their existence in the intracluster medium is established through observations of synchrotron emission and Faraday rotation. On the other hand, the nature of magnetic fields outside of clusters, where observations are scarce and controversial, remains largely unknown. In this chapter, we review recent developments in our understanding of the nature and origin of intergalactic magnetic fields, and in particular, intercluster fields. A plausible scenario for the origin of galactic and intergalactic magnetic fields is for seed fields, created in the early universe, to be amplified by turbulent flows induced during the formation of the large scale structure. We present several mechanisms for the generation of seed fields both before and after recombination. We then discuss the evolution and role of magnetic fields during the formation of the first starts. We describe the turbulent amplification of seed fields during the formation of large scale structure and the nature of the magnetic fields that arise. Finally, we discuss implications of intergalactic magnetic fields.

preprint2011arXiv

Peaks above the Maxwellian Sea: A New Approach to Finding Substructure in N-Body Haloes

We describe a new algorithm for finding substructures within dark matter haloes from N-body simulations. The algorithm relies upon the fact that dynamically distinct substructures in a halo will have a {\em local} velocity distribution that differs significantly from the mean, i.e. smooth background halo. We characterize the large-scale mean field using a coarsely grained cell-based approach, while a kernel smoothing process is used to determined the local velocity distribution. Comparing the ratio of these two estimates allows us to identify particles which are strongly cluster in velocity space relative to the background and thus resident in substructure. From this population of outliers, groups are identified using a Friends-of-Friends-like approach. False positives are rejected using Poisson noise arguments. This approach does not require a search of the full phase-space structure of a halo, a non-trivial task, and is thus computationally advantageous. We apply our algorithm to several test cases and show that it identifies not only subhaloes, bound overdensities in phase-space, but can recover tidal streams with a high purity. Our method can even find streams which do not appear significantly overdense in either physical or phase-space.

preprint2011arXiv

The First Magnetic Fields

We review current ideas on the origin of galactic and extragalactic magnetic fields. We begin by summarizing observations of magnetic fields at cosmological redshifts and on cosmological scales. These observations translate into constraints on the strength and scale magnetic fields must have during the early stages of galaxy formation in order to seed the galactic dynamo. We examine mechanisms for the generation of magnetic fields that operate prior during inflation and during subsequent phase transitions such as electroweak symmetry breaking and the quark-hadron phase transition. The implications of strong primordial magnetic fields for the reionization epoch as well as the first generation of stars is discussed in detail. The exotic, early-Universe mechanisms are contrasted with astrophysical processes that generate fields after recombination. For example, a Biermann-type battery can operate in a proto-galaxy during the early stages of structure formation. Moreover, magnetic fields in either an early generation of stars or active galactic nuclei can be dispersed into the intergalactic medium.

preprint2011arXiv

The Luminosity Profile and Structural Parameters of the Andromeda Galaxy

We have constructed an extended composite luminosity profile for the Andromeda galaxy, M31, and have decomposed it into three basic luminous structural components: a bulge, a disk and a halo. The dust-free Spitzer/IRAC imaging and extended spatial coverage of ground-based optical imaging and deep star counts allow us to map M31's structure from its center to 22 kpc along the major axis. We apply different decomposition methods for the 1D luminosity profiles and 2D images. These include non-linear least-squares and Bayesian Monte-Carlo Markov-chain analyses. The basic photometric model for M31 has a Sersic bulge with shape index n = 2.2 +/- 0.3 and effective radius R_e = 1.0 +/- 0.2 kpc, a dust-free exponential disk of scale length R_d = 5.3 +/- 0.5 kpc; the parameter errors reflect the range between various decomposition methods. The bulge parameter, n, is rather insensitive to bandpass effects and its value suggests a first rapid formation via mergers followed by secular growth from the disk. The M31 halo has a 2D power-law index ~ -2.5 +/- 0.2. We find that the M31 bulge light is mostly dominant over the range R_minor < 1.2 kpc. The disk takes over in the range 1.2 kpc < R_minor < 9 kpc, whereas the halo dominates at R_minor > 9 kpc. The stellar nucleus, bulge, disk, and halo components each contribute roughly 0.05%, 23%, 73% and 4% of the total light of M31 out to 200 kpc along the minor axis. Nominal errors for the structural parameters of the M31 bulge, disk and halo amount to 20%. If M31 and the Milky Way are at all typical, faint stellar halos should be routinely detected in galaxy surveys reaching below mu_i ~ 27 mag/arcsec^2. We stress that our results rely on this photometric analysis alone. Structural parameters may change when other fundamental constraints, such as those provided by abundance gradients and stellar kinematics, are considered simultaneously (Abridged).

preprint2010arXiv

Dynamical Models for NGC 6503 using a Markov Chain Monte Carlo Technique

We use Bayesian statistics and Markov chain Monte Carlo (MCMC) techniques to construct dynamical models for the spiral galaxy NGC 6503. The constraints include surface brightness profiles which display a Freeman Type II structure; HI and ionized gas rotation curves; the stellar rotation, which is nearly coincident with the ionized gas curve; and the line of sight stellar dispersion, with a sigma-drop at the centre. The galaxy models consist of a Sersic bulge, an exponential disc with an optional inner truncation and a cosmologically motivated dark halo. The Bayesian/MCMC technique yields the joint posterior probability distribution function for the input parameters. We examine several interpretations of the data: the Type II surface brightness profile may be due to dust extinction, to an inner truncated disc or to a ring of bright stars; and we test separate fits to the gas and stellar rotation curves to determine if the gas traces the gravitational potential. We test each of these scenarios for bar stability, ruling out dust extinction. We also find that the gas likely does not trace the gravitational potential, since the predicted stellar rotation curve is then inconsistent with the observed stellar rotation curve. The disc is well fit by an inner-truncated profile, but the possibility of ring formation by a bar to reproduce the Type II profile is also a realistic model. We further find that the halo must have a cuspy profile with gamma > 1; the bulge has a lower M/L than the disc, suggesting a star forming component in the centre of the galaxy; and the bulge, as expected for this late type galaxy, has a low Sersic index with n_b~1-2, suggesting a formation history dominated by secular evolution.

preprint2010arXiv

The photometric properties of a vast stellar substructure in the outskirts of M33

We have surveyed $\sim40$sq.degrees surrounding M33 with CFHT MegaCam in the g and i filters, as part of the Pan-Andromeda Archaeological Survey. Our observations are deep enough to resolve the top 4mags of the red giant branch population in this galaxy. We have previously shown that the disk of M33 is surrounded by a large, irregular, low-surface brightness substructure. Here, we quantify the stellar populations and structure of this feature using the PAndAS data. We show that the stellar populations of this feature are consistent with an old population with $<[Fe/H]>\sim-1.6$dex and an interquartile range in metallicity of $\sim0.5$dex. We construct a surface brightness map of M33 that traces this feature to $μ_V\simeq33$mags\,arcsec$^{-2}$. At these low surface brightness levels, the structure extends to projected radii of $\sim40$kpc from the center of M33 in both the north-west and south-east quadrants of the galaxy. Overall, the structure has an "S-shaped" appearance that broadly aligns with the orientation of the HI disk warp. We calculate a lower limit to the integrated luminosity of the structure of $-12.7\pm0.5$mags, comparable to a bright dwarf galaxy such as Fornax or AndII and slightly less than $1\$ of the total luminosity of M33. Further, we show that there is tentative evidence for a distortion in the distribution of young stars near the edge of the HI disk that occurs at similar azimuth to the warp in HI. The data also hint at a low-level, extended stellar component at larger radius that may be a M33 halo component. We revisit studies of M33 and its stellar populations in light of these new results, and we discuss possible formation scenarios for the vast stellar structure. Our favored model is that of the tidal disruption of M33 in its orbit around M31.

preprint2009arXiv

PAndAS' cubs: discovery of two new dwarf galaxies in the surroundings of the Andromeda and Triangulum galaxies

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

preprint2009arXiv

Power spectrum for the small-scale Universe

The first objects to arise in a cold dark matter universe present a daunting challenge for models of structure formation. In the ultra small-scale limit, CDM structures form nearly simultaneously across a wide range of scales. Hierarchical clustering no longer provides a guiding principle for theoretical analyses and the computation time required to carry out credible simulations becomes prohibitively high. To gain insight into this problem, we perform high-resolution (N=720^3 - 1584^3) simulations of an Einstein-de Sitter cosmology where the initial power spectrum is P(k) propto k^n, with -2.5 < n < -1. Self-similar scaling is established for n=-1 and n=-2 more convincingly than in previous, lower-resolution simulations and for the first time, self-similar scaling is established for an n=-2.25 simulation. However, finite box-size effects induce departures from self-similar scaling in our n=-2.5 simulation. We compare our results with the predictions for the power spectrum from (one-loop) perturbation theory and demonstrate that the renormalization group approach suggested by McDonald improves perturbation theory's ability to predict the power spectrum in the quasilinear regime. In the nonlinear regime, our power spectra differ significantly from the widely used fitting formulae of Peacock & Dodds and Smith et al. and a new fitting formula is presented. Implications of our results for the stable clustering hypothesis vs. halo model debate are discussed. Our power spectra are inconsistent with predictions of the stable clustering hypothesis in the high-k limit and lend credence to the halo model. Nevertheless, the fitting formula advocated in this paper is purely empirical and not derived from a specific formulation of the halo model.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

preprint2008arXiv

Dynamical Blueprints for Galaxies

We present an axisymmetric, equilibrium model for late-type galaxies which consists of an exponential disk, a Sersic bulge, and a cuspy dark halo. The model is specified by a phase space distribution function which, in turn, depends on the integrals of motion. Bayesian statistics and the Markov Chain Monte Carlo method are used to tailor the model to satisfy observational data and theoretical constraints. By way of example, we construct a chain of 10^5 models for the Milky Way designed to fit a wide range of photometric and kinematic observations. From this chain, we calculate the probability distribution function of important Galactic parameters such as the Sersic index of the bulge, the disk scale length, and the disk, bulge, and halo masses. We also calculate the probability distribution function of the local dark matter velocity dispersion and density, two quantities of paramount significance for terrestrial dark matter detection experiments. Though the Milky Way models in our chain all satisfy the prescribed observational constraints, they vary considerably in key structural parameters and therefore respond differently to non-axisymmetric perturbations. We simulate the evolution of twenty-five models which have different Toomre Q and Goldreich-Tremaine X parameters. Virtually all of these models form a bar, though some, more quickly than others. The bar pattern speeds are ~ 40 - 50 km/s/kpc at the time when they form and then decrease, presumably due to coupling of the bar with the halo. Since the Galactic bar has a pattern speed ~50 km/s/kpc we conclude that it must have formed recently.

preprint2005arXiv

HST Imaging of MEGA Microlensing Candidates in M31

We investigate $HST$/ACS and WFPC2 images at the positions of five candidate microlensing events from a large survey of variability in M31 (MEGA). Three closely match unresolved sources, and two produce only flux upper limits. All are confined to regions of the color-magnitude diagram where stellar variability is unlikely to be easily confused with microlensing. Red variable stars cannot explain these events (although background supernova are possible for two). If these lenses arise in M31's halo, they are due to masses $0.15 < m / M_\odot < 0.49$ (95% certainty, for a $δ$-function mass distribution), brown dwarfs for disk lenses, and stellar masses for bulge lenses.

preprint1995arXiv

Large Scale Structure Tests of Warm Dark Matter

Warm dark matter (WDM) is an intriguing model of structure formation from the point of view of both cosmology and particle physics. We consider a one-parameter family of WDM models. The linear power spectra for these models is calculated and compared with the corresponding spectra for cold dark matter (CDM), hot dark matter (HDM) and mixed dark matter (MDM) as well as the power spectrum derived from observations. Our linear analyses suggest that a model universe dominated by a particle whose mass to temperature ratio $m_x/T_x$ is increased by a factor of two as compared with the standard HDM neutrino gives a reasonable fit to the data on large $(>8h^{-1} {\rm ~Mpc})$ scales. $N$-body simulations for this particular WDM model show features of both HDM and CDM. As in HDM, the first objects to collapse are large pancake-like structures. The final matter distribution is rather smooth and structures as small as galaxy halos are excluded. However, there appear to be virialized rich clusters evident in the CDM but not the HDM simulations. Unfortunately, a simple comparison of the matter distribution and its statistical properties with observations indicates that WDM, like CDM, has too much power at small scales. This is particularly evident in the small-scale pairwize velocity dispersion. The cluster multiplicity function has the wrong shape with too many rich clusters being produced, though this conclusion is based on the simple assumption that light traces mass in groups of galaxies.

preprint1995arXiv

Measuring Hubble's Constant in our Inhomogeneous Universe

Recent observations of Cepheids in the Virgo cluster have bolstered the evidence that supports a Hubble constant in 70-90 km/s/Mpc range. This evidence, by and large, probes the expansion of the Universe within 100 Mpc. We investigate the possibility that the expansion rate within this region is systematically higher than the true expansion rate due to the presence of a local, large underdense region or void. We begin by calculating the expected deviations between the locally measured Hubble constant and the true Hubble constant for a variety of models. We also discuss the expected correlations between these deviations and mass fluctuation for the sample volume. We find that the fluctuations are small for the standard cold dark matter as well as mixed dark matter models but can be substantial in a number of interesting and viable nonstandard scenarios. However, deviations in the Hubble flow for a region of radius 200 Mpc are small for virtually all reasonable models. Therefore, methods based on supernovae or the Sunyaev-Zel'dovich effect, which can probe 200 Mpc scales, will be essential in determining the true Hubble constant. We discuss, in detail, the fluctuations induced in the cosmic background radiation by voids at the last scattering surface. In addition, we discuss the dipole and quadrupole fluctuations one would expect if the void enclosing us is aspherical or if we lie off-center.

preprint1994arXiv

On Stationary, Self-Similar Distributions of a Collisionless, Self-Gravitating, Gas

We study systematically stationary solutions to the coupled Vlasov and Poisson equations which have `self-similar' or scaling symmetry in phase space. In particular, we find analytically {\it all} spherically symmetric distribution functions where the mass density and gravitational potential are strict power laws in $r$, the distance from the symmetry point. We treat as special cases, systems built from purely radial orbits and systems that are isotropic in velocity space. We then discuss systems with arbitrary velocity space anisotropy finding a new and very general class of distribution functions. These distributions may prove useful in modelling galaxies. Distribution functions in cylindrical and planar geometries are also discussed. Finally, we study spatially spheroidal systems that again exhibit strict power-law behaviour for the density and potential and find results in agreement with results published recently.