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 and Gaia 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.4 $\pm$ 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 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.

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.

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.

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_\phi$. 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.

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.

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 $\gamma$ of the DG's DM density distribution, ranging from steep cusps ($\gamma=1$) to shallower density profiles ($\gamma<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 ($\gamma=1$) are favourable to forming a hard IMBH binary, whereas when $\gamma<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 $\gamma$. This suggests caution in employing cosmological simulations of galaxy formation to study BH dynamics in DGs.

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.

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.

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 \mu 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.

We present a new method for embedding a stellar disc in a cosmological dark matter halo and provide a worked example from a {\Lambda}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.

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.

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.

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.

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.

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.

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.

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.

We develop a framework for modelling the Milky Way using stellar streams and a wide range of photometric and kinematic observations. Through the use of mock data we demonstrate that a standard suite of Galactic observations leads to degeneracies in the inferred halo parameters. We then incorporate a GD-1-like stream into this suite using the orbit-fitting technique and show that the streams reduces the uncertainties in these parameters provided all observations are fit simultaneously. We also explore how the assumption of a disk-halo alignment can lead unphysical models. Our results may explain why some studies based on the Sagittarius stream find that the halo's intermediate axis is parallel to the disk spin axis even though such a configuration is highly unstable. Finally we show that both longer streams and multiple streams lead to improvements in our ability to infer the shape of our dark halo.

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

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.

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.

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.

Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33.

We present a mass model for the Milky Way, which is fit to observations of the Sagittarius stream together with constraints derived from a wide range of photometric and kinematic data. The model comprises a S\'{e}rsic bulge, an exponential disk, and an Einasto halo. Our Bayesian analysis is accomplished using an affine-invariant Markov chain Monte Carlo algorithm. We find that the best-fit dark matter halo is triaxial with axis ratios of $3.3\pm 0.7$ and $2.7\pm 0.4$ and with the short axis approximately aligned with the Sun-Galactic centre line. Our results are consistent with those presented in Law and Majewski (2010). Such a strongly aspherical halo is disfavoured by the standard cold dark matter scenario for structure formation.

We present results from a spectroscopic survey of the dwarf spheroidal And XXII and the two extended clusters EC1 and EC2. These three objects are candidate satellites of the Triangulum galaxy, M33, which itself is likely a satellite of M31. We use the DEep Imaging Multi-Object Spectrograph mounted on the Keck-II telescope to derive radial velocities for candidate member stars of these objects and thereby identify the stars that are most likely actual members. Eleven most probable stellar members (of 13 candidates) are found for AndXXII. We obtain an upper limit of sigma_v < 6.0 km s-1 for the velocity dispersion of AndXXII, [Fe/H] ~ -1.6 for its metallicity, and 255pc for the Plummer radius of its projected density profile. We construct a colour magnitude diagram for AndXXII and identify both the red giant branch and the horizontal branch. The position of the latter is used to derive a heliocentric distance to And XXII of 853 pm 26 kpc. The combination of the radial velocity, distance, and angular position of AndXXII indicates that it is a strong candidate for being the first known satellite of M33 and one of the very few examples of a galactic satellite of a satellite. N-body simulations imply that this conclusion is unchanged even if M31 and M33 had a strong encounter in the past few Gyr. We test the hypothesis that the extended clusters highlight tidally stripped galaxies by searching for an excess cloud of halo-like stars in their vicinity. We find such a cloud for the case of EC1 but not EC2. The three objects imply a dynamical mass for M33 that is consistent with previous estimates.

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.

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.

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.

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.

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

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 $\mu_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.

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.

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.

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.

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.

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.

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.

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo can not be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

We explore the secular dynamical evolution of an N-body model of M31 in the presence of a population of 100 dark matter satellites over 10 Gyr. The satellite population has structural and kinematic characteristics modelled to follow the predictions of Lambda-CDM cosmological simulations. Vertical disk heating is a small effect despite many interactions with the satellite population with only a 20% increase in vertical velocity dispersion sigma_z and the disk scale height z_d at the equivalent solar radius R = 2.5R_d . However, the stellar disk is noticeably flared after 10 Gyr with z_d nearly doubling at the disk edge. Azimuthal disk heating is much larger with sigma_R and sigma_z both increasing by 1.7x. However, in a control experiment without satellites dispersion increases by 1.5x suggesting that most of the effect is due to heating through scattering off of spiral structure excited by swing-amplified noise. Surprisingly, direct impacts of satellites on the disk can excite spiral structure with a significant amplitude and in some cases impacts close to the disk center also induce the bar instability. The large number of dark matter satellite impacts expected over a galaxy's lifetime may be a significant source of external perturbations for driving disk secular evolution.

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.

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.

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.

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.

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.

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.

We present results of a microlensing survey toward the Andromeda Galaxy (M31) carried out during four observing seasons at the Isaac Newton Telescope (INT). This survey is part of the larger microlensing survey toward M31 performed by the Microlensing Exploration of the Galaxy and Andromeda (MEGA) collaboration. Using a fully automated search algorithm, we indentify 14 candidate microlensing events, three of which are reported here for the first time. Observations obtained at the Mayall telescope are combined with the INT data to produce composite lightcurves for these candidates. The results from the survey are compared with theoretical predictions for the number and distribution of events. These predictions are based on a Monte Carlo calculation of the detection efficiency and disk-bulge-halo models for M31. The models provide the full phase-space distribution functions for the lens and source populations and are motivated by dynamical and observational considerations. They include differential extinction and span a wide range of parameter space characterised primarily by the mass-to-light ratios for the disk and bulge. For most models, the observed event rate is consistent with the rate predicted for self-lensing -- a MACHO halo fraction of 30% or higher can be ruled at the 95% confidence level. The event distribution does show a large near-far asymmetry hinting at a halo contribution to the microlensing signal. Two candidate events are located at particularly large projected radii on the far side of the disk. These events are difficult to explain by self lensing and only somewhat easier to explain by MACHO lensing. A possibility is that one of these is due to a lens in a giant stellar stream.

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 $\delta$-function mass distribution), brown dwarfs for disk lenses, and stellar masses for bulge lenses.

We describe a new set of self-consistent, equilibrium disk galaxy models that incorporate an exponential disk, a Hernquist model bulge, an NFW halo and a central supermassive black hole. The models are derived from explicit distribution functions for each component and the large number of parameters permit detailed modeling of actual galaxies. We present techniques that use structural and kinematic data such as radial surface brightness profiles, rotation curves and bulge velocity dispersion profiles to find the best-fit models for the Milky Way and M31. Through N-body realizations of these models we explore their stability against the formation of bars. The models permit the study of a wide range of dynamical phenomenon with a high degree of realism.

We present the first M31 candidate microlensing events from the Microlensing Exploration of the Galaxy and Andromeda (MEGA) survey. MEGA uses several telescopes to detect microlensing towards the nearby Andromeda galaxy, M31, in order to establish whether massive compact objects are a significant contribution to the mass budget of the dark halo of M31. The results presented here are based on observations with the Isaac Newton Telescope on La Palma, during the 1999/00 and 2000/01 observing seasons. In this data set, 14 variable sources consistent with microlensing have been detected, 12 of which are new and 2 have been reported previously by the POINT-AGAPE group. A preliminary analysis of the spatial and timescale distributions of the candidate events support their microlensing nature. We compare the spatial distributions of the candidate events and of long-period variable stars, assuming the chances of finding a long-period variable and a microlensing event are comparable. The spatial distribution of our candidate microlensing events is more far/near side asymmetric than expected from the detected long-period variable distribution. The current analysis is preliminary and the asymmetry not highly significant, but the spatial distribution of candidate microlenses is suggestive of the presence of a microlensing halo.

Terrestrial dark matter detection experiments probe the velocity-space distribution of dark matter particles in the vicinity of the Earth. We present a novel method, to be used in conjunction with standard cosmological simulations of hierarchical clustering, that allows one to extract a truly local velocity-space distribution in exquisite detail. Preliminary results suggest a new picture for this distribution which is decidedly non-Maxwellian but instead is characterized by randomly positioned peaks in velocity space. We discuss the implications of these results for both WIMP and axion detection experiments.

We present a suite of semi-analytic disk-bulge-halo models for the Andromeda galaxy (M31) which satisfy three fundamental conditions: (1) internal self-consistency; (2) consistency with observational data; and (3) stability of the disk against the formation of a central bar. The models are chosen from a set first constructed by Kuijken and Dubinski. We develop an algorithm to search the parameter space for this set in order to best match observations of the M31 rotation curve, inner velocity dispersion profile, and surface brightness profile. Models are obtained for a large range of bulge and disk masses; we find that the disk mass must be of order 8 * 10^10 M_sun and that the preferred value for the bulge mass is 2.5 * 10^10 M_sun. N-body simulations are carried out to test the stability of our models against the formation of a bar within the disk. We also calculate the baryon fraction and halo concentration parameter for a subset of our models and show that the results are consistent with the predictions from cosmological theories of structure formation. In addition, we describe how gravitational microlensing surveys and dynamical studies of globular clusters and satellites can further constrain the models.

Multiresolution analysis is applied to the problem of halo identification in cosmological N-body simulations. The procedure makes use of a discrete wavelet transform known as the algorithme a trous and segmentation analysis. It has the ability to find subhalos in the dense regions of a parent halo and can discern the multiple levels of substructure expected in the hierarchical clustering scenario. As an illustration, a 500,000 particle dark matter halo is analyzed and over 600 subhalos are found. Statistical properties of the subhalo population are discussed.

A variety of observations suggest that magnetic fields are present in all galaxies and galaxy clusters. These fields are characterized by a modest strength (10^{-7}-10^{-5} G) and huge spatial scale (~Mpc). It is generally assumed that magnetic fields in spiral galaxies arise from the combined action of differential rotation and helical turbulence, a process known as the alpha-omega dynamo. However fundamental questions concerning the nature of the dynamo as well as the origin of the seed fields necessary to prime it remain unclear. Moreover, the standard alpha-omega dynamo does not explain the existence of magnetic fields in elliptical galaxies and clusters. The author summarizes what is known observationally about magnetic fields in galaxies, clusters, superclusters, and beyond. He then reviews the standard dynamo paradigm, the challenges that have been leveled against it, and several alternative scenarios. He concludes with a discussion of astrophysical and early Universe candidates for seed fields.

In the cold dark matter model of structure formation, galaxies are assembled hierarchically from mergers and the accretion of subclumps. This process is expected to leave residual substructure in the Galactic dark halo, including partially disrupted clumps and their associated tidal debris. We develop a model for such halo substructure and study its implications for dark matter (WIMP and axion) detection experiments. We combine the Press-Schechter model for the distribution of halo subclump masses with N-body simulations of the evolution and disruption of individual clumps as they orbit through the evolving Galaxy to derive the probability that the Earth is passing through a subclump or stream of a given density. Our results suggest that it is likely that the local complement of dark matter particles includes a 1-5% contribution from a single clump. The implications for dark matter detection experiments are significant, since the disrupted clump is composed of a `cold' flow of high-velocity particles. We describe the distinctive features due to halo clumps that would be seen in the energy and angular spectra of detection experiments. The annual modulation of these features would have a different signature and phase from that for a smooth halo and, in principle, would allow one to discern the direction of motion of the clump relative to the Galactic center.

January 2000 completes the first season of intensive, wide-field observations of microlensing and stellar variability in M31 by MEGA (Microlensing Exploration of the Galaxy and Andromeda) at the Isaac Newton 2.5m Telescope, the KPNO 4m, and the 1.3m and 2.5m telescopes of MDM Observatory. In preliminary analysis, we detect ~50000 variable objects, including some consistent with microlensing events. We present the level of sensitivity to be reached in our planned three-year program to test for the presence of a significant halo microlensing population in M31, as well as its spatial distribution and mass-function. We also discuss our application of image subtraction to these wide fields and HST WFPC2 Snapshot followup observations to confirm candidates identified from previous years' surveys. We present intermediate results from our smaller-field survey, on the MDM 1.3m and Vatican Advanced Technology 1.8m Telescope, from 1994-1998, wherein we have discovered 8 additional probable microlensing events, over about one-half the time base of the project, in addition to confirming three of our original 6 microlensing candidates from 1995.

Various analytic expressions have been proposed for the density profile of dark matter halos. We consider six of these expressions for which the density profile has a power-law fall-off rho proportional to r^{-3} at large radii and a power-law cusp rho proportional to r^{-gamma} (gamma=0, 1/2, 1, 3/2) at small radii. The phase-space distribution function for these models is calculated assuming spherical symmetry and either an isotropic velocity dispersion tensor or an anisotropic dispersion tensor of the type proposed by Osipkov and Merritt. The differential energy distribution for these models is also derived. Several applications are discussed including the analysis of dark matter search experiments and the study of halo formation in a cosmological setting. Analytic fitting formulae for some of the models are provided.

We demonstrate that the Biermann battery mechanism for the creation of large scale magnetic fields can arise in a simple model protogalaxy. Analytic calculations and numerical simulations follow explicitly the generation of vorticity (and hence magnetic field) at the outward-moving shock that develops as the protogalactic perturbation collapses. Shear angular momentum then distorts this field into a dipole-like configuration. The magnitude of the field created in the fully formed disk galaxy is estimated to be 10^(-17) Gauss, approximately what is needed as a seed for the galactic dynamo.

Using the VLA, we have performed the first observational test of dark matter in the form of cold, primordial fractal clouds, as envisioned by Pfenniger, Combes, & Martinet (1994) and Pfenniger & Combes (1994). We show that, after a Hubble Time, primordial fractal clouds will convert most of their HI to H2, but a small fraction of HI remains which is optically thick. This opens up a new window for detecting dark matter which may exist in this form. The detectability of such gas depends on its filling factor and temperature and therefore should be observable in absorption against a background source, with observations of sufficient sensitivity and resolution. The current VLA observations have made a first step towards this goal by taking advantage of a fortuitous alignment between the extension of the HI disk of the nearby galaxy, NGC 3079, and a background quasar Q 0957+561. Our observations probe 28 independent beams against the quasar and all of velocity space between the extension of a flat rotation curve and a Keplerian decline for the halo region of NGC 3079. We do not detect any absorption features and investigate, in detail, the implication of this result for the hypothesis that dark matter is in the form of fractal clouds. While our analysis does not rule out fractal clouds as dark matter, it does lay out the groundwork for future, more sensitive observations and we consider what form these might take to probe the range of possible cloud properties more deeply. It is interesting that the observations can rule out cold optically thin HI gas, if it exists, to a limit of 0.001% of the dark matter. In contrast, the existence of cold HI in a fractal hierarchy would be an efficient way of hiding dark matter.

Navarro, Frenk, and White have suggested that the density profiles of simulated dark matter halos have a ``universal'' shape so that a given halo can be characterized by a single free parameter which fixes its mass. In this paper, we revisit the spherical infall model in the hope of recognizing in detail the existence and origin of any such universality. A system of particles is followed from linear perturbation, through first shell crossing, then through an accretion or infall phase, and finally to virialization. During the accretion phase, the system relaxes through a combination of phase mixing, phase space instability, and moderate violent relation. It is driven quickly, by the flow of mass through its surface, toward self-similar evolution. The self-similar solution plays its usual role of intermediate attractor and can be recognized from a virial-type theorem in scaled variables and from our numerical simulations. The transition to final equilibrium state once infall has ceased is relatively gentle, an observation which leads to an approximate form for the distribution function of the final system. The infall phase fixes the density profile in intermediate regions of the halo to be close to r^{-2}. We make contact with the standard hierarchical clustering scenario and explain how modifications of the self-similar infall model might lead to density profiles in agreement with those found in numerical simulations.

We conduct gravitational microlensing experiments in a galaxy taken from a cosmological N-body simulation. Hypothetical observers measure the optical depth and event rate toward hypothetical LMCs and compare their results with model predictions. Since we control the accuracy and sophistication of the model, we can determine how good it has to be for statistical errors to dominate over systematic ones. Several thousand independent microlensing experiments are performed. When the ``best-fit'' triaxial model for the mass distribution of the halo is used, the agreement between the measured and predicted optical depths is quite good: by and large the discrepancies are consistent with statistical fluctuations. If, on the other hand, a spherical model is used, systematic errors dominate. Even with our ``best-fit'' model, there are a few rare experiments where the deviation between the measured and predicted optical depths cannot be understood in terms of statistical fluctuations. In these experiments there is typically a clump of particles crossing the line of sight to the hypothetical LMC. These clumps can be either gravitationally bound systems or transient phenomena in a galaxy that is still undergoing phase mixing. Substructure of this type, if present in the Galactic distribution of Machos, can lead to large systematic errors in the analysis of microlensing experiments. We also describe how hypothetical WIMP and axion detection experiments might be conducted in a simulated N-body galaxy.

We consider the evolution of an initially cold, spherically symmetric system of self-gravitating particles. At early times the system undergoes self-similar collapse of the type described by Fillmore & Goldreich and Bertschinger. This stage of phase mixing soon gives way to a period of violent relaxation driven by an instability in the similarity solution. The onset of violent relaxation is illustrated by numerical simulations and supported by non-linear analysis of the ``scaled'' collisionless Boltzmann and Poisson equations.

The Schr\"odinger Method is a novel approach for modeling numerically self-gravitating, collisionless systems that may have certain advantages over N-body and phase space methods. In particular, smoothing is part of the dynamics and not just the force calculation. This paper describes test-bed simulations which illustrate the viability of the Schr\"odinger Method. We develop the techniques necessary to handle ``hot'' systems as well as spherically symmetric systems. We demonstrate that the method can maintain a stable, equilibrium star cluster by constructing and then evolving a Plummer sphere. We also consider nonequilibrium initial conditions and follow the system as it attempts to reach virial equilibrium through phase mixing and violent relaxation. Finally we make a few remarks concerning the dynamics of axions and other bosonic dark matter candidates. The Schr\"odinger Method, in principle, provides an exact treatment of these fields. However such ``scalar field'' simulations are feasible and warrented only if the deBroglie wavelength of the particle is comparable to the size of the system of interest, a situation that is almost certainly not the case for axions in the Galaxy. Therefore the Schr\"odinger Method treats axions in the same way as other collisionless particles. We challenge recent claims in the literature that axions in the Galaxy form soliton stars.

We propose a new numerical technique for following the evolution of a self-gravitating collisionless system in general relativity. Matter is modeled as a scalar field obeying the coupled Klein-Gordon and Einstein equations. A phase space distribution function, constructed using covariant coherent states, obeys the relativistic Vlasov equation provided the de Broglie wavelength for the field is very much smaller than the scales of interest. We illustrate the method by solving for the evolution of a system of particles in a static, plane-symmetric, background spacetime.

It is widely believed that dark matter halos are flattened, that is closer to oblate than prolate. The evidence cited is based largely on observations of galaxies which do not look anything like our own and on numerical simulations which use ad hoc initial conditions. Given what we believe to be a ``reasonable doubt'' concerning the shape of dark Galactic halo we calculate the optical depth and event rate for microlensing of stars in the LMC assuming a wide range of models that include both prolate and oblate halos. We find, in agreement with previous analysis, that the optical depth for a spherical (E0) halo and for an oblate (E6) halo are roughly the same, essentially because two competing effects cancel approximately. However the optical depth for an E6 prolate halo is reduced by ~35%. This means that an all-Macho prolate halo with reasonable parameters for the Galaxy is consistent with the published microlensing event rate.

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

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

We investigate a new hybrid-model universe containing two types of dark matter, one ``warm'' and the other ``hot''. The hot component is an ordinary light neutrino with mass $\sim 25h^2$~eV while the warm component is a sterile neutrino with mass $\sim 700h^2$~eV. The two types of dark matter arise entirely within the neutrino sector and do not require separate physical origins. We calculate the linear transfer functions for a representative sample of warm-plus-hot models. The transfer functions, and results from several observational tests of structure formation, are compared with those for the cold-plus-hot models that have been studied extensively in the literature. On the basis of these tests, we conclude that warm-plus-hot dark matter is essentially indistinguishable from cold-plus-hot dark matter, and therefore provides a viable scenario for large scale structure. We demonstrate that a neutrino mass matrix can be constructed which provides the requisite dark matter constituents, while remaining consistent with all cosmological bounds.

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

We propose that the recently reported MACHO/EROS events correspond to gravitational amplification by dark clouds rather than compact objects. These clouds must be very dense with $M\sim 0.1 ~M_\odot$ and $R\la 10^{14}~{\rm cm}$. In all likelihood, the clouds will be members of a family of objects with different sizes and masses. We therefore expect events of longer duration than the ones reported by the MACHO and EROS groups but with light curves very different from the ones derived assuming point mass lenses. We suggest that one such event has already been observed in radio measurements of the quasar 1502+106. The abundances of free electrons, metals, complex molecules, and dust grains are constrained to be very small suggesting that the clouds are formed from a primordial mixture of hydrogen and helium. Cosmic rays and background UV radiation ionize a halo around the cloud. Radio waves from distant sources will be scattered by the electrons in this halo, an effect which may have already been observed in quasars such as 1502+106. We argue that dark clouds are a viable alternative to compact objects for baryonic dark matter in the halo.

The simplest model that can accomodate a viable nonbaryonic dark matter candidate is the standard electroweak theory with the addition of right-handed or sterile neutrinos. We reexamine this model and find that the sterile neutrinos can be either hot, warm, or cold dark matter. Since their only direct coupling is to left-handed or active neutrinos, the most efficient production mechanism is via neutrino oscillations. If the production rate is always less than the expansion rate, then these neutrinos will never be in thermal equilibrium. However, enough of them may be produced so that they provide the missing mass necessary for closure. We consider a single generation of neutrino fields $\left (\nu_L,\,\nu_R\right )$ with a Dirac mass, $\mu$, and a Majorana mass for the right-handed components only, $M$. For $M\gg \mu$ we show that the number density of sterile neutrinos is proportional to $\mu^2/M$ so that the energy density today is {\it independent of} $M$. However $M$ is crucial in determining the large scale structure of the Universe. In particular, $M\simeq 0.1-1.0 {\rm ~keV}$ leads to warm dark matter and a structure formation scenario that may have some advantages over both the standard hot and cold dark matter scenarios.