Showing 1–50 of 72 results for author: Widrow, L M

Linear operator theory of phase mixing

Authors: Keir Darling, Lawrence M. Widrow

Abstract: We study solutions of the collisionless Boltzmann equation (CBE) in a functional Koopman representation. This facilitates the use of linear spectral techniques characteristic of the analysis of Schrodinger-type equations. For illustrative purposes, we consider the classical phase mixing of a non-interacting distribution function in a quartic potential. Solutions are determined perturbatively relat… ▽ More We study solutions of the collisionless Boltzmann equation (CBE) in a functional Koopman representation. This facilitates the use of linear spectral techniques characteristic of the analysis of Schrodinger-type equations. For illustrative purposes, we consider the classical phase mixing of a non-interacting distribution function in a quartic potential. Solutions are determined perturbatively relative to a harmonic oscillator. We impose a form of coarse-graining by choosing a finite dimensional basis to represent the distribution function and time evolution operators, which sets a minimum length scale on phase space structure. We observe a relationship between the dimension of the representation and the multiplicity of the harmonic oscillator eigenvalues. System dynamics are understood in terms of degenerate subspaces of the linear operator spectra. Each subspace is associated with a mode of the harmonic oscillator, the first two being bending and breathing structures. The quartic potential splits the degenerate eigenvalues within each subspace. This facilitates the formation of spiral structure as deformations from the harmonic oscillator modes. We ultimately argue that this construction provides a promising avenue for study of self-interacting systems experiencing phase mixing, which is an outstanding problem in the context of the Gaia DR2 vertical phase space spirals. △ Less

Submitted 21 June, 2024; v1 submitted 25 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures. Changes: expanded discussion and conclusion, replaced Fig. 1 with new Fig. 1 and 2, moved definitions from Section 2 to Appendix, clarified wording regarding multiplicity and length scale, added steps to Section 5.1, added second order correction curves to Fig. 3, added equation 57, added Fig. 5 to show subspace evolution, changed dashed curves to solid in contour plots

Galactic Structure From Binary Pulsar Accelerations: Beyond Smooth Models

Authors: Thomas Donlon II, Sukanya Chakrabarti, Lawrence M. Widrow, Michael T. Lam, Philip Chang, Alice C. Quillen

Abstract: We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρ_0=0.062\pm0.017$ \m… ▽ More We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($ρ_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($ρ_{0,\textrm{DM}}=-0.010\pm0.018$ \msun/pc$^3$); these values are similar to, but have smaller uncertainties than previous pulsar timing measurements of these quantities. Here, we provide for the first time, values for the Oort constants and the slope of the rotation curve from direct acceleration measurements. We find that $A=15.4\pm2.6$ km/s/kpc and $B=-13.1\pm2.6$ km/s/kpc (consistent with results from \textit{Gaia}), and the slope of the rotation curve near the Sun is $-2\pm5$ km/s/kpc. We show that the Galactic acceleration field is clearly asymmetric, but due to data limitations it is not yet clear which physical processes drive this asymmetry. We provide updated models of the Galactic potential that account for various sources of disequilibrium; these models are incompatible with commonly used kinematic potentials. This indicates that use of kinematically derived Galactic potentials in precision tests (e.g., in tests of general relativity with pulsar timing) may be subject to larger uncertainties than reported. The acceleration data indicates that the mass of the Galaxy within the Solar circle is $2.3 \times 10^{11}$ M$_\odot$, roughly twice as large as currently accepted models. Additionally, the residuals of the acceleration data compared to existing Galactic models have a dependence on radial position; this trend can be explained if the Sun has an additional acceleration away from the Galactic center. △ Less

Submitted 20 May, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

The Formation of Bars and Warps in Rotating Halos

Authors: Robin Joshi, Lawrence M. Widrow

Abstract: We investigate the effects of halo kinematics on the dynamics of stellar discs by simulating the evolution of isolated disc-halo systems from equilibrium initial conditions. Our main results come from four simulations where the initial disc is identical and the halo is either treated as a rigid potential or is live with isotropic orbits or orbits that preferentially rotate with or counter to the d… ▽ More We investigate the effects of halo kinematics on the dynamics of stellar discs by simulating the evolution of isolated disc-halo systems from equilibrium initial conditions. Our main results come from four simulations where the initial disc is identical and the halo is either treated as a rigid potential or is live with isotropic orbits or orbits that preferentially rotate with or counter to the disc. We confirm previous results that bar formation is more vigorous in models with a live halo than a rigid one and is further enhanced when halo orbits preferentially rotate with the disc. We discuss two types of buckling events with different symmetries about the mid plane, one that occurs just as the bar is forming and the other well after the bar has been established. We also show that warps are most easily excited and maintained when the halo is counter-rotating with the disc, in agreement with theoretical predictions. Our most novel result is the discovery of a rotating halo instability, which causes the disc and halo cusp to spiral outward from the centre of mass of the system whether the halo rotates with the disc or counter to it and also occurs in a disc bulge halo system that does not form a bar. We provide a heuristic linear model that captures the essential dynamics of the instability. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Comments: 14 pages, 13 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Fuzzy Dark Matter Dynamics and the Quasiparticle Hypothesis

Authors: Boris Zupancic, Lawrence M. Widrow

Abstract: Dark matter may be composed of ultra-light bosons whose de Broglie wavelength in galaxies is of order 1 kpc. The standard model for this fuzzy dark matter (FDM) is a complex scalar field that obeys the Schrödinger-Poisson equations. The wavelike nature of FDM leads to fluctuations in the gravitational field that can pump energy into the stellar components of a galaxy. Heuristic arguments and theor… ▽ More Dark matter may be composed of ultra-light bosons whose de Broglie wavelength in galaxies is of order 1 kpc. The standard model for this fuzzy dark matter (FDM) is a complex scalar field that obeys the Schrödinger-Poisson equations. The wavelike nature of FDM leads to fluctuations in the gravitational field that can pump energy into the stellar components of a galaxy. Heuristic arguments and theoretical analyses suggest that these fluctuations can be modelled by replacing FDM with a system of quasiparticles (QPs). We test this hypothesis by comparing self-consistent simulations of a Schrödinger field with those using a system of QPs in one spatial dimension. Simulations of pure FDM systems allow us to derive a phenomenological relation between the number of QPs that is required to model FDM with a given de Broglie wavelength. We also simulate systems of FDM and stars and find that the FDM pumps energy into the stars whether it is described by QPs or a Schrödinger field with the FDM adiabatically contracting and the stellar system adiabatically expanding. However, we find that QPs overestimate dynamical heating. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 9 pages, 10 figures, accepted for publication MNRAS

Swing Amplification and the Gaia Phase Spirals

Authors: Lawrence M. Widrow

Abstract: We explore the interplay between in-plane and vertical dynamics in stellar discs within the framework of the shearing box approximation. Julian and Toomre used the shearing sheet to show that leading density waves are amplified as they swing into a trailing ones. We extend their formalism into the dimension perpendicular to the disc and obtain explicit solutions for the response of a disc to an im… ▽ More We explore the interplay between in-plane and vertical dynamics in stellar discs within the framework of the shearing box approximation. Julian and Toomre used the shearing sheet to show that leading density waves are amplified as they swing into a trailing ones. We extend their formalism into the dimension perpendicular to the disc and obtain explicit solutions for the response of a disc to an impulsive, external excitation. An excitation that is is symmetric about the mid plane produces a density/breathing wave as well as two-armed phase spirals in the vertical phase space plane. On the other hand, an excitation that is antisymmetric about the mid plane leads to a bending wave and single-armed phase spirals. In either case, self-gravity plays a crucial role in driving the evolution of the disturbance and determining the amplitude and pitch angle of the ensuing spirals. We also show that when the disc is excited by a co-rotating cloud, it develops stationary phase spirals in the wake of the cloud. The results call into question simple kinematic arguments that have been used to determine the age of the phase spirals seen in the Gaia survey △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 10 pages, 11 figures, submitted to MNRAS

Mapping Milky Way disk perturbations in stellar number density and vertical velocity using Gaia DR3

Authors: Axel Widmark, Lawrence M. Widrow, Aneesh Naik

Abstract: We have mapped the number density and mean vertical velocity of the Milky Way's stellar disk out to roughly two kiloparsecs from the Sun using Gaia Data Release 3 (DR3) and complementary photo-astrometric distance information from StarHorse. For the number counts, we carefully masked spatial regions that are compromised by open clusters, great distances, or dust extinction and used Gaussian proces… ▽ More We have mapped the number density and mean vertical velocity of the Milky Way's stellar disk out to roughly two kiloparsecs from the Sun using Gaia Data Release 3 (DR3) and complementary photo-astrometric distance information from StarHorse. For the number counts, we carefully masked spatial regions that are compromised by open clusters, great distances, or dust extinction and used Gaussian processes to arrive at a smooth, non-parametric estimate for the underlying number density field. We find that the number density and velocity fields depart significantly from an axisymmetric and mirror-symmetric model. These departures, which include projections of the Gaia phase-space spiral, signal the presence of local disturbances in the disk. We identify two features that are present in both stellar number density and mean vertical velocity. One of these features appears to be associated with the Local Spiral Arm. It is most prominent at small heights and is largely symmetric across the mid-plane of the disk. The density and velocity field perturbations are phase-shifted by roughly a quarter wavelength, suggesting a breathing mode that is propagating in the direction of Galactic longitude $l\sim 270$ deg. The second feature is a gradient in the stellar number density and mean vertical velocity with respect to Galactocentric radius. This feature, which extends across the entire region of our analysis, may be associated with the extension of the Galactic warp into the Solar neighbourhood in combination with more localised bending waves. △ Less

Submitted 31 October, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 13 pages, 11 figures, appendices add 12 pages, 11 figures; accepted for publication in A&A

Journal ref: A&A 668, A95 (2022)

Gaussian Process Model for the Local Stellar Velocity Field from Gaia Data Release 2

Authors: Patrick Nelson, Lawrence M. Widrow

Abstract: We model the local stellar velocity field using position and velocity measurements for 4M stars from the second data release of Gaia. We determine the components of the mean or bulk velocity in ~27k spatially-defined bins. Our assumption is that these quantities constitute a Gaussian process where the correlation between the bulk velocity at different locations is described by a simple covariance… ▽ More We model the local stellar velocity field using position and velocity measurements for 4M stars from the second data release of Gaia. We determine the components of the mean or bulk velocity in ~27k spatially-defined bins. Our assumption is that these quantities constitute a Gaussian process where the correlation between the bulk velocity at different locations is described by a simple covariance function or kernel. We use a sparse Gaussian process algorithm based on inducing points to construct a non-parametric, smooth, and differentiable model for the underlying velocity field. We estimate the Oort constants A, B, C, and K and find values in excellent agreement with previous results. Maps of the velocity field within 2 kpc of the Sun reveal complicated substructures, which provide clear evidence that the local disk is in a state of disequilibrium. We present the first 3D map of the divergence of the stellar velocity field and identify regions of the disk that may be undergoing compression and rarefaction. △ Less

Submitted 14 June, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 10 pages, 10 figures; typos corrected, figure format

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

Authors: Eric J. Mendelsohn, Heidi Jo Newberg, Siddhartha Shelton, Lawrence M. Widrow, Jeffery M. Thompson, Carl J. Grillmair

Abstract: 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 a… ▽ More 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. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 33 pages, 53 figures, to be published by ApJ

Inferring Time-Dependent Distribution Functions from Kinematic Snapshots

Authors: Keir Darling, Lawrence M. Widrow

Abstract: We propose a method for constructing the time-dependent phase space distribution function (DF) of a collisionless system from an isolated kinematic snapshot. In general, the problem of mapping a single snapshot to a time-dependent function is intractable. Here we assume a finite series representation of the DF, constructed from the spectrum of the system's Koopman operator. This reduces the origin… ▽ More We propose a method for constructing the time-dependent phase space distribution function (DF) of a collisionless system from an isolated kinematic snapshot. In general, the problem of mapping a single snapshot to a time-dependent function is intractable. Here we assume a finite series representation of the DF, constructed from the spectrum of the system's Koopman operator. This reduces the original problem to one of mapping a kinematic snapshot to a discrete spectrum rather than to a time-dependent function. We implement this mapping with a convolutional neural network (CNN). The method is demonstrated on two example models: the quantum simple harmonic oscillator, and a self-gravitating isothermal plane. The latter system exhibits phase space spiral structure similar to that observed in Gaia Data Release 2. △ Less

Submitted 1 July, 2021; v1 submitted 6 February, 2021; originally announced February 2021.

Comments: 14 pages, 11 figures

The Stellar Distribution Function and Local Vertical Potential from Gaia DR2

Authors: Haochuan Li, Lawrence M. Widrow

Abstract: We develop a novel method to simultaneously determine the vertical potential, force and stellar $z-v_z$ phase space distribution function (DF) in our local patch of the Galaxy. We assume that the Solar Neighborhood can be treated as a one-dimensional system in dynamical equilibrium and directly fit the number density in the $z-v_z$ plane to what we call the Rational Linear DF (RLDF) model. This mo… ▽ More We develop a novel method to simultaneously determine the vertical potential, force and stellar $z-v_z$ phase space distribution function (DF) in our local patch of the Galaxy. We assume that the Solar Neighborhood can be treated as a one-dimensional system in dynamical equilibrium and directly fit the number density in the $z-v_z$ plane to what we call the Rational Linear DF (RLDF) model. This model can be regarded as a continuous sum of isothermal DFs though it has only one more parameter than the isothermal model. We apply our method to a sample of giant stars from Gaia Data Release 2 and show that the RLDF provides an excellent fit to the data. The well-known phase space spiral emerges in the residual map of the $z-v_z$ plane. We use the best-fit potential to plot the residuals in terms of the frequency and angle of vertical oscillations and show that the spiral maps into a straight line. From its slope, we estimate that the phase spirals were generated by a perturbation $\sim540$ Myr years ago. We also determine the differential surface density as a function of vertical velocity dispersion, a.k.a. the vertical temperature distribution. The result is qualitatively similar to what was previously found for SDSS/SEGUE G dwarfs. Finally, we address parameter degeneracies and the validity of the 1D approximation. Particularly, the mid-plane density derived from a cold subsample, where the 1D approximation is more secure, is closer to literature values than that derived from the sample as a whole. △ Less

Submitted 1 March, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: 15 pages, 16 figures

The Milky Way's Shell Structure Reveals the Time of a Radial Collision

Authors: Thomas Donlon II, Heidi Jo Newberg, Robyn Sanderson, Lawrence M. Widrow

Abstract: We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenit… ▽ More We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenitor dwarf galaxy passed through the Galactic center 2.7 +/- 0.2 Gyr ago. Based on the time of collision, it is also possible that the VRM is related to the phenomenon that created phase-space spirals in the vertical motion of the disk and/or the Splash, and could have caused a burst of star formation in the inner disk. We analyze phase mixing in a collection of radial merger N-body simulations, and find that shell structure similar to that observed in Milky Way data disappears by 5 Gyr after collision with the Galactic center. The method used to calculate the merger time of the VRM was able to reliably recover the correct merger times for these simulations. Previous work supports the idea that the VRM and the Gaia Sausage/Gaia-Enceladus Merger are the same. However, the Gaia Sausage is widely believed to be 8--11 Gyr old. The disparate ages could be reconciled if the larger age is associated with an infall time when the progenitor crossed the virial radius; we do not constrain the time at which the progenitor became bound to the Milky Way. Alternatively, the Gaia Sausage could be younger than previously thought. △ Less

Submitted 9 September, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

GalactICS with Gas

Authors: N. Deg, L. M. Widrow, T. Randriamampandry, C. Carignan

Abstract: We present a new version of the GalactICS code that can generate self-consistent equilibrium galaxy models with a two-component stellar disc and a gas disc as well as a centrally-concentrated bulge and extended dark halo. The models can serve as initial conditions for simulations of isolated galaxies that include both hydrodynamics and collisionless dynamics. We test the code by evolving a pair of… ▽ More We present a new version of the GalactICS code that can generate self-consistent equilibrium galaxy models with a two-component stellar disc and a gas disc as well as a centrally-concentrated bulge and extended dark halo. The models can serve as initial conditions for simulations of isolated galaxies that include both hydrodynamics and collisionless dynamics. We test the code by evolving a pair of simple gas disc-halo models, which differ only in the initial temperature of the gas component. The models are similar to the ones considered in the Wang et al. except that here, the halo is live whereas they included the halo as a fixed potential. We find that the basic structural properties of the models, such as the rotation curve and surface density profiles, are well-preserved over 1.5 Gyr. We also construct a Milky Way model that includes thin and thick stellar disc components, a gas disc, a bulge, and a dark halo. Bar formation occurs in all disc-like components at about 1 Gyr. The bar is strongest in the thin disc while the gas disc contains the most prominent spiral features. The length of the bar in our model is comparable to what has been inferred for the Galactic bar. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: 10 pages, 13 figures, accepted to MNRAS

Eigenfunctions of Galactic Phase Space Spirals from Dynamic Mode Decomposition

Authors: Keir Darling, Lawrence M. Widrow

Abstract: We investigate the spatiotemporal structure of simulations of the homogeneous slab and isothermal plane models for the vertical motion in the Galactic disc. We use Dynamic Mode Decomposition (DMD) to compute eigenfunctions of the simulated distribution functions for both models, referred to as DMD modes. In the case of the homogeneous slab, we compare the DMD modes to the analytic normal modes of… ▽ More We investigate the spatiotemporal structure of simulations of the homogeneous slab and isothermal plane models for the vertical motion in the Galactic disc. We use Dynamic Mode Decomposition (DMD) to compute eigenfunctions of the simulated distribution functions for both models, referred to as DMD modes. In the case of the homogeneous slab, we compare the DMD modes to the analytic normal modes of the system to evaluate the feasibility of DMD in collisionless self gravitating systems. This is followed by the isothermal plane model, where we focus on the effect of self gravity on phase mixing. We compute DMD modes of the system for varying relative dominance of mutual interaction and external potential, so as to study the corresponding variance in mode structure and lifetime. We find that there is a regime of relative dominance, at approximately $ 4:1 $ external potential to mutual interaction where the DMD modes are spirals in the $ (z,v_z) $ plane, and are nearly un-damped. This leads to the proposition that a system undergoing phase mixing in the presence of weak to moderate self gravity can have persisting spiral structure in the form of such modes. We then conclude with the conjecture that such a mechanism may be at work in the phase space spirals observed in Gaia Data Release 2, and that studying more complex simulations with DMD may aid in understanding both the timing and form of the perturbation that lead to the observed spirals. △ Less

Submitted 5 September, 2019; v1 submitted 18 April, 2019; originally announced April 2019.

Comments: 8 pages, 6 figures

The large-scale structure of the halo of the Andromeda galaxy II. Hierarchical structure in the Pan-Andromeda Archaeological Survey

Authors: Alan W. McConnachie, Rodrigo Ibata, Nicolas Martin, Annette M. N. Ferguson, Michelle Collins, Stephen Gwyn, Mike Irwin, Geraint F. Lewis, A. Dougal Mackey, Tim Davidge, Veronica Arias, Anthony Conn, Patrick Cote, Denija Crnojevic, Avon Huxor, Jorge Penarrubia, Chelsea Spengler, Nial Tanvir, David Valls-Gabaud, Arif Babul, Pauline Barmby, Nicholas F. Bate, Edouard Bernard, Scott Chapman, Aaron Dotter , et al. (7 additional authors not shown)

Abstract: The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections betwee… ▽ More The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products. △ Less

Submitted 25 October, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: Accepted for publication in the Astrophysical Journal. 51 pages, 24 figures, 5 tables. Some figures have degraded resolution. All PAndAS data products are available via the CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/pandas/query.html where you can also find a version of the paper with full resolution figures

Can Stellar Discs in a Cosmological Setting Avoid Forming Strong Bars?

Authors: Jacob S. Bauer, Lawrence M. Widrow

Abstract: We investigate the connection between the vertical structure of stellar discs and the formation of bars using high-resolution simulations of galaxies in isolation and in the cosmological context. In particular, we simulate a suite of isolated galaxy models that have the same Toomre Q parameter and swing amplification parameter but that differ in the vertical scale height and velocity dispersion. W… ▽ More We investigate the connection between the vertical structure of stellar discs and the formation of bars using high-resolution simulations of galaxies in isolation and in the cosmological context. In particular, we simulate a suite of isolated galaxy models that have the same Toomre Q parameter and swing amplification parameter but that differ in the vertical scale height and velocity dispersion. We find that the onset of bar formation occurs more slowly in models with thicker discs. Moreover, thicker discs and also discs evolved in simulations with larger force softening also appear to be more resilient to buckling, which acts to regulate the length and strength of bars. We also simulate disc-halo systems in the cosmological environment using a disc-insertion technique developed in a previous paper. In this case, bar formation is driven by the stochastic effects of a triaxial halo and subhalo-disc interactions and the initial growth of bars appears to be relatively insensitive to the thickness of the disc. On the other hand, thin discs in cosmological halos do appear to be more susceptible to buckling than thick ones and therefore bar strength correlates with disc thickness as in the isolated case. More to the point, one can form discs in cosmological simulations with relatively weak bars or no bars at all provided the discs as thin as the discs we observe and the softening length is smaller than the disc scale height. △ Less

Submitted 31 August, 2018; originally announced September 2018.

Comments: 16 pages, 13 figures. Submitted to MNRAS

Emergence of the Gaia Phase Space Spirals from Bending Waves

Authors: Keir Darling, Lawrence M. Widrow

Abstract: We discuss the physical mechanism by which pure vertical bending waves in a stellar disc evolve to form phase space spirals similar to those discovered by Antoja et al. ( arXiv:1804.10196) in Gaia Data Release 2. These spirals were found by projecting Solar Neighbourhood stars onto the $z-v_z$ plane. Faint spirals appear in the number density of stars projected onto the $z-v_z$ plane, which can be… ▽ More We discuss the physical mechanism by which pure vertical bending waves in a stellar disc evolve to form phase space spirals similar to those discovered by Antoja et al. ( arXiv:1804.10196) in Gaia Data Release 2. These spirals were found by projecting Solar Neighbourhood stars onto the $z-v_z$ plane. Faint spirals appear in the number density of stars projected onto the $z-v_z$ plane, which can be explained by a simple model for phase wrapping. More prominent spirals are seen when bins across the $z-v_z$ plane are coloured by median $v_R$ or $v_φ$. We use both toy model and fully self-consistent simulations to show that the spirals develop naturally from vertical bending oscillations of a stellar disc. The underlying physics follows from the observation that the vertical energy of a star (essentially, its "radius" in the $z-v_z$ plane) correlates with its angular momentum or, alternatively, guiding radius. Moreover, at fixed physical radius, the guiding radius determines the azimuthal velocity. Together, these properties imply the link between in-plane and vertical motion that lead directly to the Gaia spirals. We show that the cubic $R-z$ coupling term in the effective potential is crucial for understanding the morphology of the spirals. This suggests that phase space spirals might be a powerful probe of the Galactic potential. In addition, we argue that self-gravity is necessary to properly model the evolution of the bending waves and their attendant phase space spirals. △ Less

Submitted 14 November, 2018; v1 submitted 30 July, 2018; originally announced July 2018.

Comments: 7 Pages, 9 Figures, submitted to MNRAS on July 30, 2018

Simulating non-axisymmetric flows in disc galaxies

Authors: T. H. Randriamampandry, N. Deg, C. Carignan, L. M. Widrow

Abstract: We present a two-step method to simulate and study non-circular motions in strongly barred galaxies. The first step is to constrain the initial parameters using a Bayesian analysis of each galaxy's azimuthally averaged rotation curve, the 3.6 $\rm μm$ surface brightness, and the gas surface density. The second step is to generate equilibrium models using the GalactICS code and evolve them via GADG… ▽ More We present a two-step method to simulate and study non-circular motions in strongly barred galaxies. The first step is to constrain the initial parameters using a Bayesian analysis of each galaxy's azimuthally averaged rotation curve, the 3.6 $\rm μm$ surface brightness, and the gas surface density. The second step is to generate equilibrium models using the GalactICS code and evolve them via GADGET-2. The bar strengths and mock velocity maps of the resulting snapshots are compared to observations in order to determine the best representation of the galaxy. We test our method on the unbarred galaxy NGC 3621 and the barred galaxies NGC 1300 and NGC 1530. NGC 3621 provides a validation of our method of generating initial conditions. NGC 1530 has an intermediate bar orientation that allows for a comparison to DiskFit. Finally NGC 1300 has a bar oriented parallel to the galaxy's major axis, where other algorithms tend to fail. Our models for NGC 3621 and NGC 1530 are comparable to those obtained using commonly available algorithms. Moreover, we have produced one of the first mass distribution models for NGC 1300. △ Less

Submitted 15 July, 2018; originally announced July 2018.

Comments: 19 pages, 15 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 618, A106 (2018)

Formation of LISA Black Hole Binaries in Merging Dwarf Galaxies: the Imprint of Dark Matter

Authors: Tomas Tamfal, Pedro R. Capelo, Stelios Kazantzidis, Lucio Mayer, Douglas Potter, Joachim Stadel, Lawrence M. Widrow

Abstract: Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (… ▽ More Theoretical models for the expected merger rates of intermediate-mass black holes (IMBHs) are vital for planned gravitational-wave detection experiments such as the Laser Interferometer Space Antenna (LISA). Using collisionless $N$-body simulations of dwarf galaxy (DG) mergers, we examine how the orbital decay of IMBHs and the efficiency of IMBH binary formation depend on the central dark matter (DM) density profile of the merging DGs. Specifically, we explore various asymptotic inner slopes $γ$ of the DG's DM density distribution, ranging from steep cusps ($γ=1$) to shallower density profiles ($γ<1$), motivated by well-known baryonic-feedback effects as well as by DM models that differ from cold DM at the scales of DGs. We find that the inner DM slope is crucial for the formation (or lack thereof) of an IMBH binary; only mergers between DGs with cuspy DM profiles ($γ=1$) are favourable to forming a hard IMBH binary, whereas when $γ<1$ the IMBHs stall at a separation of 50-100 pc. Consequently, the rate of LISA signals from IMBH coalescence will be determined by the fraction of DGs with a cuspy DM profile. Conversely, the LISA event rates at IMBH mass scales offer in principle a novel way to place constraints on the inner structure of DM halos in DGs and address the core-cusp controversy. We also show that, with spatial resolutions of $\sim$0.1 kpc, as often adopted in cosmological simulations, all IMBHs stall, independent of $γ$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs. △ Less

Submitted 27 September, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 9 pages, 4 figures, 1 table, accepted for publication in ApJL

Bending Waves in the Milky Way's disc from halo substructure

Authors: Matthew H. Chequers, Lawrence M. Widrow, Keir Darling

Abstract: We use $N$-body simulations to investigate the excitation of bending waves in a Milky Way-like disc-bulge-halo system. The dark matter halo consists of a smooth component and a population of subhaloes while the disc is composed of thin and thick components. Also considered is a control simulation where all of the halo mass is smoothly distributed. We find that bending waves are more vigorously exc… ▽ More We use $N$-body simulations to investigate the excitation of bending waves in a Milky Way-like disc-bulge-halo system. The dark matter halo consists of a smooth component and a population of subhaloes while the disc is composed of thin and thick components. Also considered is a control simulation where all of the halo mass is smoothly distributed. We find that bending waves are more vigorously excited in the thin disc than the thick one and that they are strongest in the outer regions of the disc, especially at late times. By way of a Fourier decomposition, we find that the complicated pattern of bending across the disc can be described as a superposition of waves, which concentrate along two branches in the radius-rotational frequency plane. These branches correspond to vertical resonance curves as predicted by a WKB analysis. Bending waves in the simulation with substructure have a higher amplitude than those in the smooth-halo simulation, though the frequency-radius characteristics of the waves in the two simulations are very similar. A cross correlation analysis of vertical displacement and bulk vertical velocity suggests that the waves oscillate largely as simple plane waves. We suggest that the wave-like features in astrometric surveys such as the Second Data Release from \textit{Gaia} may be due to long-lived waves of a dynamically active disc rather than, or in addition to, perturbations from a recent satellite-disc encounter. △ Less

Submitted 3 August, 2018; v1 submitted 31 May, 2018; originally announced May 2018.

Comments: 16 pages, 14 figures, 1 table, accepted for publication in MNRAS on Aug, 2, 2018

Disc-Halo Interactions in ΛCDM

Authors: Jacob S. Bauer, Lawrence M. Widrow, Denis Erkal

Abstract: We present a new method for embedding a stellar disc in a cosmological dark matter halo and provide a worked example from a ΛCDM zoom-in simulation. The disc is inserted into the halo at a redshift z = 3 as a zero-mass rigid body. Its mass and size are then increased adiabatically while its position, velocity, and orientation are determined from rigid-body dynamics. At z = 1, the rigid disc is rep… ▽ More We present a new method for embedding a stellar disc in a cosmological dark matter halo and provide a worked example from a ΛCDM zoom-in simulation. The disc is inserted into the halo at a redshift z = 3 as a zero-mass rigid body. Its mass and size are then increased adiabatically while its position, velocity, and orientation are determined from rigid-body dynamics. At z = 1, the rigid disc is replaced by an N-body disc whose particles sample a three-integral distribution function (DF). The simulation then proceeds to z = 0 with live disc and halo particles. By comparison, other methods assume one or more of the following: the centre of the rigid disc during the growth phase is pinned to the minimum of the halo potential, the orientation of the rigid disc is fixed, or the live N-body disc is constructed from a two rather than three-integral DF. In general, the presence of a disc makes the halo rounder, more centrally concentrated, and smoother, especially in the innermost regions. We find that methods in which the disc is pinned to the minimum of the halo potential tend to overestimate the amount of adiabatic contraction. Additionally, the effect of the disc on the subhalo distribution appears to be rather insensitive to the disc insertion method. The live disc in our simulation develops a bar that is consistent with the bars seen in late-type spiral galaxies. In addition, particles from the disc are launched or "kicked up" to high galactic latitudes. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: Accepted to MNRAS 2018 January 08. Received 2018 January 04; in original form 2017 August 17

Spontaneous generation of bending waves in isolated Milky Way-like discs

Authors: Matthew H. Chequers, Lawrence M. Widrow

Abstract: We study the spontaneous generation and evolution of bending waves in $N$-body simulations of two isolated Milky Way-like galaxy models. The models differ by their disc-to-halo mass ratios, and hence by their susceptibility to the formation of a bar and spiral structure. Seeded from shot noise in the particle distribution, bending waves rapidly form in both models and persist for many billions of… ▽ More We study the spontaneous generation and evolution of bending waves in $N$-body simulations of two isolated Milky Way-like galaxy models. The models differ by their disc-to-halo mass ratios, and hence by their susceptibility to the formation of a bar and spiral structure. Seeded from shot noise in the particle distribution, bending waves rapidly form in both models and persist for many billions of years. Waves at intermediate radii manifest as corrugated structures in vertical position and velocity that are tightly wound, morphologically leading, and dominated by the $m=1$ azimuthal Fourier component. A spectral analysis of the waves suggests they are a superposition of modes from two continuous branches in the Galactocentric radius-rotational frequency plane. The lower-frequency branch is dominant and is responsible for the corrugated, leading, and warped structure. Over time, power in this branch migrates outward, lending credence to an inside-out formation scenario for the warp. Our power spectra qualitatively agree with results from linear perturbation theory and a WKB analysis, both of which include self-gravity. Thus, we conclude that the waves in our simulations are self-gravitating and not purely kinematic. These waves are reminiscent of the wave-like pattern recently found in Galactic star counts from the Sloan Digital Sky Survey and smoothly transition to a warp near the disc's edge. Velocity measurements from \textit{Gaia} data will be instrumental in testing the true wave nature of the corrugations. We also compile a list of "minimum requirements" needed to observe bending waves in external galaxies. △ Less

Submitted 22 August, 2017; originally announced August 2017.

Comments: Accepted for publication in MNRAS, 15 pages, 10 figures, 1 table

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

Authors: Matthew H. Chequers, Kristine Spekkens, Lawrence M. Widrow, Colleen Gilhuly

Abstract: 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 gl… ▽ More 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. △ Less

Submitted 29 August, 2016; originally announced August 2016.

Comments: 9 pages, 6 figures, 1 table, accepted for publication in MNRAS

Galactoseismology and the Local Density of Dark Matter

Authors: Nilanjan Banik, Lawrence M. Widrow, Scott Dodelson

Abstract: 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… ▽ More 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. △ Less

Submitted 6 October, 2016; v1 submitted 10 August, 2016; originally announced August 2016.

Comments: 9 pages, 11 figures, submitted to MNRAS, minor typos corrected, references added and conclusion improved

Report number: FERMILAB-PUB-16-315-A-AE-PPD-T

Journal ref: MNRAS (February 01, 2017) 464 (4): 3775-3783

Estimating the Galactic Mass Profile in the Presence of Incomplete Data

Authors: Gwendolyn M. Eadie, William E. Harris, Lawrence M. Widrow

Abstract: 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;… ▽ More 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. △ Less

Submitted 24 March, 2015; originally announced March 2015.

Comments: Accepted to ApJ; 14 pages, 7 figures, 2 tables

Journal ref: ApJ 806:54 (2015) 14pp

Vertical Oscillations of Fluid and Stellar Disks

Authors: Lawrence M. Widrow, Gage Bonner

Abstract: 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, fo… ▽ More 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. △ Less

Submitted 19 March, 2015; originally announced March 2015.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

Bending and Breathing Modes of the Galactic Disk

Authors: Lawrence M. Widrow, Jarrett Barber, Matthew H. Chequers, Edward Cheng

Abstract: 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 plung… ▽ More 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. △ Less

Submitted 15 April, 2014; originally announced April 2014.

Comments: 12 pages, 12 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2013, 440, 1971

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

Authors: Claire E. Dorman, Lawrence M. Widrow, Puragra Guhathakurta, the PHAT Collaboration

Abstract: 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 bri… ▽ More 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. △ Less

Submitted 16 October, 2013; originally announced October 2013.

Comments: 4 pages, 3 figures. To appear in the proceedings of the conference "Structure and Dynamics of Disk Galaxies." Full-length paper accepted to ApJ

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

Authors: Claire E. Dorman, Lawrence M. Widrow, Puragra Guhathakurta, Anil C. Seth, Daniel Foreman-Mackey, Eric F. Bell, Julianne J. Dalcanton, Karoline M. Gilbert, Evan D. Skillman, Benjamin F. Williams

Abstract: 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… ▽ More 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. △ Less

Submitted 15 October, 2013; originally announced October 2013.

Comments: 25 pages, 20 figures, accepted for publication in ApJ

Galaxy Masses

Authors: S. Courteau, M. Cappellari, R. S. de Jong, A. A. Dutton, E. Emsellem, H. Hoekstra, L. V. E. Koopmans, G. A. Mamon, Claudia Maraston, T. Treu, L. M. Widrow

Abstract: Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing me… ▽ More Galaxy masses play a fundamental role in our understanding of structure formation models. This review addresses the variety and reliability of mass estimators that pertain to stars, gas, and dark matter. The different sections on masses from stellar populations, dynamical masses of gas-rich and gas-poor galaxies, with some attention paid to our Milky Way, and masses from weak and strong lensing methods, all provide review material on galaxy masses in a self-consistent manner. △ Less

Submitted 28 April, 2014; v1 submitted 12 September, 2013; originally announced September 2013.

Comments: 145 pages, 28 figures, to appear in Reviews of Modern Physics. Figure 22 is missing here, and Figs. 15, 26-28 are at low resolution. This version has a slightly different title and some typos fixed in Chapter 5. For the full review with figures, please consult: http://www.astro.queensu.ca/~courteau/GalaxyMasses_28apr2014.pdf

Journal ref: Rev Mod Phys 86 (2014), 47-119

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

Authors: Markus Hartmann, Victor P. Debattista, David R. Cole, Monica Valluri, Lawrence M. Widrow, Juntai Shen

Abstract: 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 sati… ▽ More 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. △ Less

Submitted 10 September, 2013; originally announced September 2013.

Comments: submitted to MNRAS, 19 pages, 13 figures, 5 tables, comments are very welcome!

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

Authors: Rodrigo A. Ibata, Geraint F. Lewis, Anthony R. Conn, Michael J. Irwin, Alan W. McConnachie, Scott C. Chapman, Michelle L. Collins, Mark Fardal, Annette M. N. Ferguson, Neil G. Ibata, A. Dougal Mackey, Nicolas F. Martin, Julio Navarro, R. Michael Rich, David Valls-Gabaud, Lawrence M. Widrow

Abstract: 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, a… ▽ More 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. △ Less

Submitted 3 January, 2013; originally announced January 2013.

Comments: Published in the 3rd Jan 2013 issue of Nature. 19 pages, 4 figures, 1 three-dimensional interactive figure. To view and manipulate the 3-D figure, an Adobe Reader browser plug-in is required; alternatively save to disk and view with Adobe Reader

Journal ref: Nature 493, 62-65 (2013)

Galactoseismology: Discovery of Vertical Waves in the Galactic Disk

Authors: Lawrence M. Widrow, Susan Gardner, Brian Yanny, Scott Dodelson, Hsin-Yu Chen

Abstract: 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 Gal… ▽ More 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. △ Less

Submitted 30 March, 2012; originally announced March 2012.

Comments: 10 pages, 5 figures, accepted for publication in the Astrophysical Journal Letters

Magnetic fields in the Large-Scale Structure of the Universe

Authors: Dongsu Ryu, Dominik R. G. Schleicher, Rudolf A. Treumann, Christos G. Tsagas, Lawrence M. Widrow

Abstract: 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 unders… ▽ More 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. △ Less

Submitted 19 September, 2011; originally announced September 2011.

Comments: Accepted for publication in Space Science Reviews. Pdf with full resolution figures can be downloaded from http://canopus.cnu.ac.kr/ryu/cosmic-mag2.pdf

The First Magnetic Fields

Authors: Lawrence M. Widrow, Dongsu Ryu, Dominik Schleicher, Kandaswamy Subramanian, Christos G. Tsagas, Rudolf A. Treumann

Abstract: 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… ▽ More 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. △ Less

Submitted 19 September, 2011; originally announced September 2011.

Comments: Accepted for publication in Space Science Reviews. Pdf can be also downloaded from http://canopus.cnu.ac.kr/ryu/cosmic-mag1.pdf

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

Authors: Pascal J. Elahi, Robert J. Thacker, Lawrence M. Widrow

Abstract: 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, whil… ▽ More 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. △ Less

Submitted 21 July, 2011; originally announced July 2011.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

The Luminosity Profile and Structural Parameters of the Andromeda Galaxy

Authors: Stephane Courteau, Lawrence M. Widrow, Michael McDonald, Puragra Guhathakurta, Karoline M. Gilbert, Yucong Zhu, Rachael Lynn Beaton, Steven R. Majewski

Abstract: 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 d… ▽ More 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). △ Less

Submitted 17 June, 2011; originally announced June 2011.

Comments: 43 pages, 10 figures. Accepted for publication at ApJ

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

Authors: Alan W. McConnachie, Annette M. N. Ferguson, Michael J. Irwin, John Dubinski, Lawrence M. Widrow, Aaron Dotter, Rodrigo Ibata, Geraint F. Lewis

Abstract: 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 ste… ▽ More 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. △ Less

Submitted 14 September, 2010; originally announced September 2010.

Comments: Accepted for publication in ApJ. 17 figures. ApJ preprint format

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

Authors: David Puglielli, Lawrence M. Widrow, Stéphane Courteau

Abstract: 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… ▽ More 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. △ Less

Submitted 8 April, 2010; v1 submitted 3 February, 2010; originally announced February 2010.

Comments: 57 pages, 22 figures, ApJ accepted. Minor changes made to match accepted article.

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

Authors: Nicolas F. Martin, Alan W. McConnachie, Mike Irwin, Lawrence M. Widrow, Annette M. N. Ferguson, Rodrigo A. Ibata, John Dubinski, Arif Babul, Scott Chapman, Mark Fardal, Geraint F. Lewis, Julio Navarro, R. Michael Rich

Abstract: 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-H… ▽ More 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. △ Less

Submitted 17 September, 2009; v1 submitted 2 September, 2009; originally announced September 2009.

Comments: 10 pages, 6 figures, accepted for publication in ApJ; v2: minor typographical corrections

Journal ref: Astrophys.J.705:758-765,2009

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

Authors: Alan W. McConnachie, Michael J. Irwin, Rodrigo A. Ibata, John Dubinski, Lawrence M. Widrow, Nicolas F. Martin, Patrick Cote, Aaron L. Dotter, Julio F. Navarro, Annette M. N. Ferguson, Thomas H. Puzia, Geraint F. Lewis, Arif Babul, Pauline Barmby, Olivier Bienayme, Scott C. Chapman, Robert Cockcroft, Michelle L. M. Collins, Mark A. Fardal, William E. Harris, Avon Huxor, A. Dougal Mackey, Jorge Penarrubia, R. Michael Rich, Harvey B. Richer , et al. (4 additional authors not shown)

Abstract: 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 format… ▽ More 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. △ Less

Submitted 2 September, 2009; originally announced September 2009.

Comments: Published in Nature. Supplementary movie available at https://www.astrosci.ca/users/alan/PANDAS/Latest%20news%3A%20movie%20of%20orbit.html

Journal ref: Nature 461:66-69,2009

Can substructure in the Galactic Halo explain the ATIC and PAMELA results?

Authors: P. J. Elahi, L. M. Widrow, R. J. Thacker

Abstract: Recently, ATIC and PAMELA measured an anomalously large flux of leptonic cosmic rays which may arise from dark matter self-annihilation. While the annihilation signal predicted for a smooth halo is 10^2-10^3 times smaller than the measured excess, the signal can be boosted by the presence of subhalos. We investigate the feasibility of large boost factors using a new Monte Carlo calculation techn… ▽ More Recently, ATIC and PAMELA measured an anomalously large flux of leptonic cosmic rays which may arise from dark matter self-annihilation. While the annihilation signal predicted for a smooth halo is 10^2-10^3 times smaller than the measured excess, the signal can be boosted by the presence of subhalos. We investigate the feasibility of large boost factors using a new Monte Carlo calculation technique that is constrained by previous simulation work on halo substructure. The model accounts for the observed decrease in the amount of substructure with decreasing halo mass and the scatter in halo structural parameters such as the density concentration parameter. Our results suggest that boost factors of the order of 10^2 are ruled out at more than 14 sigma. We conclude that substructure alone, at least with commonly assumed annihilation cross-sections, cannot explain the anomalous flux measured by ATIC and PAMELA. △ Less

Submitted 23 October, 2009; v1 submitted 23 June, 2009; originally announced June 2009.

Comments: 6 pages, 2 figures, submitted to Phys Rev D

Journal ref: Phys.Rev.D80:123513,2009

Power spectrum for the small-scale Universe

Authors: Lawrence M. Widrow, Pascal J. Elahi, Robert J. Thacker, Mark Richardson, Evan Scannapieco

Abstract: 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 prohibi… ▽ More 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. △ Less

Submitted 7 July, 2009; v1 submitted 28 January, 2009; originally announced January 2009.

Comments: 30 pages including 10 figures; accepted for publication in MNRAS

Subhaloes in Scale-Free Cosmologies

Authors: P. J. Elahi, R. J. Thacker, L. M. Widrow, E. Scannapieco

Abstract: We explore the dependence of the subhalo mass function on the spectral index n of the linear matter power spectrum using scale-free Einstein-de Sitter simulations with n=-1 and n=-2.5. We carefully consider finite volume effects that may call into question previous simulations of n<-2 power spectra. Subhaloes are found using a 6D friends-of-friends algorithm in all haloes originating from high-s… ▽ More We explore the dependence of the subhalo mass function on the spectral index n of the linear matter power spectrum using scale-free Einstein-de Sitter simulations with n=-1 and n=-2.5. We carefully consider finite volume effects that may call into question previous simulations of n<-2 power spectra. Subhaloes are found using a 6D friends-of-friends algorithm in all haloes originating from high-sigma peaks. For n=-1, we find that the cumulative subhalo mass function is independent of the parameters used in the subhalo finding algorithm and is consistent with the subhalo mass function found in LCDM simulations. In particular, the subhalo mass function is well fit by a power-law with an index of alpha=-0.9, that is the mass function has roughly equal mass in subhaloes per logarithmic interval in subhalo mass. Conversely, for n=-2.5, the algorithm parameters affect the subhalo mass function since subhaloes are more triaxial with less well defined boundaries. We find that the index alpha is generally larger with alpha>=-0.75. We infer that although the subhalo mass function appears to be independent of n so long as n>=-2, it begins to flatten as n->-3. Thus, the common practice of using alpha=-1.0 may greatly overestimate the number of subhaloes at the smallest scales in the CDM hierarchy. △ Less

Submitted 9 March, 2009; v1 submitted 2 November, 2008; originally announced November 2008.

Comments: 13 pages, 8 figures, to be published in MNRAS

Dynamical Models for Disk Galaxies with Triaxial Halos

Authors: Lawrence M. Widrow

Abstract: We construct self-consistent dynamical models for disk galaxies with triaxial, cuspy halos. We begin with an equilibrium, axisymmetric, disk-bulge-halo system and apply an artificial acceleration to the halo particles. By design, this acceleration conserves energy and thereby preserving the system's differential energy distribution even as its phase space distribution function is altered. The ha… ▽ More We construct self-consistent dynamical models for disk galaxies with triaxial, cuspy halos. We begin with an equilibrium, axisymmetric, disk-bulge-halo system and apply an artificial acceleration to the halo particles. By design, this acceleration conserves energy and thereby preserving the system's differential energy distribution even as its phase space distribution function is altered. The halo becomes triaxial but its spherically-averaged density profile remains largely unchanged. The final system is in equilibrium, to a very good approximation, so long as the halo's shape changes adiabatically. The disk and bulge are ``live'' while the halo is being deformed; they respond to the changing gravitational potential but also influence the deformation of the halo. We test the hypothesis that halo triaxiality can explain the rotation curves of low surface brightness galaxies by modelling the galaxy F568-3. △ Less

Submitted 31 January, 2008; originally announced January 2008.

Comments: 23 pages including 9 figures

Dynamical Blueprints for Galaxies

Authors: Lawrence M. Widrow, Brent Pym, John Dubinski

Abstract: 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 co… ▽ More 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. △ Less

Submitted 22 January, 2008; originally announced January 2008.

Comments: 54 pages, 20 figures

Journal ref: Astrophys.J.679:1239-1259,2008

The Vcirc-sigma0 Relation of Galaxies

Authors: Stephane Courteau, Michael McDonald, Lawrence M. Widrow

Abstract: Courteau et al. (2007a) reported on the dependence of the ratio of a galaxy's maximum circular velocity, Vcirc, to its central velocity dispersion, sigma0, on morphology, or equivalently total light concentration. This Vcirc-sigma0 concentration relation, which involves details about the local and global galaxy physics, poses a fundamental challenge for galaxy structure models. Furthermore, not… ▽ More Courteau et al. (2007a) reported on the dependence of the ratio of a galaxy's maximum circular velocity, Vcirc, to its central velocity dispersion, sigma0, on morphology, or equivalently total light concentration. This Vcirc-sigma0 concentration relation, which involves details about the local and global galaxy physics, poses a fundamental challenge for galaxy structure models. Furthermore, not only must these models reproduce the Vcirc-sigma0 relation and its various dependences, they must simultaneously match other fundamental scaling relations such as the velocity-size-luminosity and color-luminosity relations of galaxies. We focus here on the interpretation of parameters that enter the Vcirc-sigma0 relation to enable proper data-model comparisons and follow-up studies by galaxy modelers and observers. △ Less

Submitted 23 September, 2007; originally announced September 2007.

Comments: 5 pages (including questions), 1 figure. To appear in the Proceedings of the IAU Symposium 245, "Formation and Evolution of Galaxy Bulges", held at Oxford, U.K., July 2007, Eds. M. Bureau, E. Athanassoula, B. Barbuy

The Bulge-Halo Connection in Galaxies: A Physical Interpretation of the Vcirc-sigma_0 Relation

Authors: Stephane Courteau, Michael McDonald, Lawrence M. Widrow, Jon Holtzman

Abstract: We explore the dependence of the ratio of a galaxy's circular velocity, Vcirc, to its central velocity dispersion, sigma_0, on morphology, or equivalently total light concentration. Such a dependence is expected if light traces the mass. Over the full range of galaxy types, masses and brightnesses, and assuming that the gas velocity traces the circular velocity, we find that galaxies obey the re… ▽ More We explore the dependence of the ratio of a galaxy's circular velocity, Vcirc, to its central velocity dispersion, sigma_0, on morphology, or equivalently total light concentration. Such a dependence is expected if light traces the mass. Over the full range of galaxy types, masses and brightnesses, and assuming that the gas velocity traces the circular velocity, we find that galaxies obey the relation log(Vcirc/sigma_0)= 0.63-0.11*C28 where C28=5log(r80/r20) and the radii are measured at 80 percent and 20 percent of the total light. Massive galaxies scatter about the Vcirc = sqrt(2)*sigma_0 line for isothermal stellar systems. Disk galaxies follow the simple relation Vcirc/sigma_0=2(1-B/T), where B/T is the bulge-to-total light ratio. For pure disks, C28~2.8, B/T -> 0, and Vcirc~=2*sigma_0. Self-consistent equilibrium galaxy models from Widrow & Dubinski (2005) constrained to match the size-luminosity and velocity-luminosity relations of disk galaxies fail to match the observed Vcirc/sigma_0 distribution. Furthermore, the matching of dynamical models for Vcirc(r)/sigma(r) with observations of dwarf and elliptical galaxies suffers from limited radial coverage and relatively large error bars; for dwarf systems, however, kinematical measurements at the galaxy center and optical edge suggest Vcirc(Rmax) > 2*sigma_0 (in contrast with past assumptions that Vcirc = sqrt(2)*sigma_0 for dwarfs.) The Vcirc-sigma_0-C28 relation has direct implications for galaxy formation and dynamical models, galaxy scaling relations, the mass function of galaxies, and the links between respective formation and evolution processes for a galaxy's central massive object, bulge, and dark matter halo. △ Less

Submitted 1 November, 2006; v1 submitted 24 October, 2006; originally announced October 2006.

Comments: Accepted for publication in ApJL. Current version matches ApJL page requirement

Journal ref: Astrophys.J.655:L21-L24,2007

Substructure around M31 : Evolution and Effects

Authors: Jean-Rene Gauthier, John Dubinski, Lawrence M. Widrow

Abstract: We investigate the evolution of a population of 100 dark matter satellites orbiting in the gravitational potential of a realistic model of M31. We find that after 10 Gyr, seven subhalos are completely disrupted by the tidal field of the host galaxy. The remaining satellites suffer heavy mass loss and overall, 75% of the mass initially in the subhalo system is tidally stripped. Not surprisingly,… ▽ More We investigate the evolution of a population of 100 dark matter satellites orbiting in the gravitational potential of a realistic model of M31. We find that after 10 Gyr, seven subhalos are completely disrupted by the tidal field of the host galaxy. The remaining satellites suffer heavy mass loss and overall, 75% of the mass initially in the subhalo system is tidally stripped. Not surprisingly, satellites with pericentric radius less than 30 kpc suffer the greatest stripping and leave a complex structure of tails and streams of debris around the host galaxy. Assuming that the most bound particles in each subhalo are kinematic tracers of stars, we find that the halo stellar population resulting from the tidal debris follows an r^{-3.5} density profile at large radii. We construct B-band photometric maps of stars coming from disrupted satellites and find conspicuous features similar both in morphology and brightness to the observed Giant Stream around Andromeda. An assumed star formation efficiency of 5-10% in the simulated satellite galaxies results in good agreement with the number of M31 satellites, the V-band surface brightness distribution, and the brightness of the Giant Stream. During the first 5 Gyr, the bombardment of the satellites heats and thickens the disk by a small amount. At about 5 Gyr, satellite interations induce the formation of a strong bar which, in turn, leads to a significant increase in the velocity dispersion of the disk. △ Less

Submitted 12 June, 2006; v1 submitted 31 May, 2006; originally announced June 2006.

Comments: 45 pages, 18 figures. To be submitted to the Astrophysical Journal, version 2.0 : scale height value corrected, references added, and some figures have been modified

Journal ref: Astrophys.J.653:1180-1193,2006

On Universal Halos and the Radial Orbit Instability

Authors: Joseph D. MacMillan, Lawrence M. Widrow, Richard N. Henriksen

Abstract: The radial orbit instability drives dark matter halos toward a universal structure. This conclusion, first noted by Huss, Jain, and Steinmetz, is explored in detail through a series of numerical experiments involving the collapse of an isolated halo into the non-linear regime. The role played by the radial orbit instability in generating the density profile, shape, and orbit structure is careful… ▽ More The radial orbit instability drives dark matter halos toward a universal structure. This conclusion, first noted by Huss, Jain, and Steinmetz, is explored in detail through a series of numerical experiments involving the collapse of an isolated halo into the non-linear regime. The role played by the radial orbit instability in generating the density profile, shape, and orbit structure is carefully analyzed and, in all cases, the instability leads to universality independent of initial conditions. New insights into the underlying physics of the radial orbit instability are presented. △ Less

Submitted 19 April, 2006; originally announced April 2006.

Comments: 31 pages, 11 figures, submitted to the Astrophysical Journal

Journal ref: Astrophys.J.653:43-52,2006

HST Imaging of MEGA Microlensing Candidates in M31

Authors: Patrick Cseresnjes, Arlin P. S. Crotts, Jelte T. A. de Jong, Alex Bergier, Edward A. Baltz, Geza Gyuk, Konrad Kuijken, Lawrence M. Widrow

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2005; v1 submitted 12 July, 2005; originally announced July 2005.

Comments: Accepted for publication in ApJL. Higher resolution version available at http://www.astro.columbia.edu/~patrick/hst/hst_ml.pdf

Journal ref: Astrophys.J.633:L105-L108,2005