The steepness of the central density profiles of dark matter (DM) in low-mass galaxy halos (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhalos using stellar streams. We show that the present day morphological and dynamical properties of accreted globular cluster (GC) streams - those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) - are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy CDM subhalos produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhalos. A first comparison of MW streams "GD-1" and "Jhelum" (likely of accreted GC origin) with our simulations indicates a preference for cored subhalos. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhalos naturally possessed cored density profiles implying DM models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g., parallel structures and "spurs"). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming "gaps" in streams when the remnant of the parent subhalo later passes through the stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhalos previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

We present Forstand, a new code for constructing dynamical models of galaxies with the Schwarzschild orbit-superposition method. These models are constrained by line-of-sight kinematic observations and applicable to galaxies of all morphological types, including disks and triaxial rotating bars. Our implementation has several novel and improved features, is computationally efficient, and made publicly available. Using mock datasets taken from N-body simulations, we demonstrate that the pattern speed of a bar can be recovered with an accuracy of 10-20%, regardless of orientation, if the 3D shape of the galaxy is known or inferred correctly.

Direct detection experiments have set increasingly stringent limits on the cross section for spin-independent dark matter-nucleon interactions. In obtaining such limits, experiments primarily assume the standard halo model (SHM) as the distribution of dark matter in our Milky Way. Three astrophysical parameters are required to define the SHM: the local dark matter escape velocity, the local dark matter density and the circular velocity of the sun around the center of the galaxy. This paper studies the effect of the uncertainties in these three astrophysical parameters on the XENON1T exclusion limits using the publicly available DDCalc code. We compare limits obtained using the widely assumed escape velocity from the RAVE survey and the newly calculated escape velocity by Monari $et$ $al.$ using Gaia data. Our study finds that the astrophysical uncertainties are dominated by the uncertainty in the escape velocity (independent of the best fit value) at dark matter masses below 6 GeV and can lead to a variation of nearly 6 orders of magnitude in the exclusion limits at 4 GeV. Above a WIMP mass of 6 GeV, the uncertainty becomes dominated by the local dark matter density, leading to uncertainties of factors of $\sim$10 (3) at 6 (15) GeV WIMP mass in the exclusion limits. Additionally, this work finds that the updated best fit value for the escape velocity based on Gaia data leads to only very minor changes to the effects of the astrophysical uncertainties on the XENON1T exclusion limits.

The rotation curve of the Milky Way is commonly used to estimate the local dark matter density $\rho_{{\rm DM},\odot}$. However, the estimates are subject to the choice of the distribution of baryons needed in this type of studies. In this work we explore several Galactic mass models that differ in the distribution of baryons and dark matter, in order to determine $\rho_{{\rm DM},\odot}$. For this purpose we analyze the precise circular velocity curve measurement of the Milky Way up to $\sim 25$ kpc from the Galactic centre obtained from Gaia DR2 [1]. We find that the estimated value of $\rho_{{\rm DM},\odot}$ stays robust to reasonable changes in the spherical dark matter halo. However, we show that $\rho_{{\rm DM},\odot}$ is affected by the choice of the model for the underlying baryonic components. In particular, we find that $\rho_{{\rm DM},\odot}$ is mostly sensitive to uncertainties in the disk components of the Galaxy. We also show that, when choosing one particular baryonic model, the estimate of $\rho_{{\rm DM},\odot}$ has an uncertainty of only about $10\%$ of its best-fit value, but this uncertainty gets much bigger when we also consider the variation of the baryonic model. In particular, the rotation curve method does not allow to exclude the presence of an additional very thin component, that can increase $\rho_{{\rm DM},\odot}$ by more than a factor of 8 (the thin disk could even be made of dark matter). Therefore, we conclude that exclusively using the rotation curve of the Galaxy is not enough to provide a robust estimate of $\rho_{{\rm DM},\odot}$. For all the models that we study without the presence of an additional thin component, our resulting estimates of the local dark matter density take values in the range $\rho_{{\rm DM},\odot} \simeq \text{0.3--0.4}\,\mathrm{GeV/cm^3}$, consistent with many of the estimates in the literature.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

We present a Cepheid-based distance to the nearby Seyfert galaxy NGC\,6814 from {\it Hubble Space Telescope} observations. We obtained F555W and F814W imaging over the course of 12 visits with logarithmic time spacing in 2013 August$-$October. We detected and made photometric measurements for 16,469 unique sources across all images in both filters, from which we identify 90 excellent Cepheid candidates spanning a range of periods of $13-84$ days. We find evidence for incompleteness in the detection of candidates at periods <21 days. Based on the analysis of Cepheid candidates above the incompleteness limit, we determine a distance modulus for NGC\,6814 relative to the LMC of $\mu_{\rm rel\,LMC}=13.200^{+0.031}_{-0.031}$ mag. Adopting the recent constraint of the distance modulus to the LMC determined by Pietrzynski et al., we find $m-M=31.677^{+0.041}_{-0.041}$ which gives a distance of $21.65 \pm 0.41$ Mpc to NGC 6814.

A recently discovered young, high-velocity giant star J01020100-7122208 is a good candidate of hypervelocity star ejected from the Galactic center, although it has a bound orbit. If we assume that this star was ejected from the Galactic center, it can be used to constrain the Galactic potential, because the deviation of its orbit from a purely radial orbit informs us of the torque that this star has received after its ejection. Based on this assumption, we estimate the flattening of the dark matter halo of the Milky Way by using the Gaia DR2 data and the circular velocity data from Eilers et al. (2019). Our Bayesian analysis shows that the orbit of J01020100-7122208 favors a prolate dark matter halo within $\sim$ 10 kpc from the Galactic center. The posterior distribution of the density flattening $q$ shows a broad distribution at $q \gtrsim1$ and peaks at $q \simeq 1.5$. Also, 98.5\% of the posterior distribution is located at $q>1$, highly disfavoring an oblate halo.

We review the Schwarzschild orbit-superposition approach and present a new implementation of this method, which can deal with a large class of systems, including rotating barred disk galaxies. We discuss two conceptuals problems in this field: the intrinsic degeneracy of determining the potential from line-of-sight kinematics, and the non-uniqueness of deprojection and related biases in potential inference, especially acute for triaxial bars. When applied to mock datasets with known 3d shape, our method correctly recovers the pattern speed and other potential parameters. However, more work is needed to systematically address these two problems for real observational datasets.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

Tidally disrupted globular cluster streams are usually observed, and therefore perceived, as narrow, linear and one-dimensional structures in the 6D phase-space. Here we show that the GD-1 stellar stream ($\approx$ 30 pc wide), which is the tidal debris of a disrupted globular cluster, possesses a secondary diffuse and extended stellar component ($\approx$ 100 pc wide) around it, detected at >5$\sigma$ confidence level. Similar morphological properties are seen in synthetic streams that are produced from star clusters that are formed within dark matter sub-halos and then accrete onto a massive host galaxy. This lends credence to the idea that the progenitor of the highly retrograde GD-1 stream was originally formed outside of the Milky Way in a now defunct dark satellite galaxy. We deem that in future studies, this newly found $cocoon$ component may serve as a structural hallmark to distinguish between the in-situ and ex-situ (accreted) formed globular cluster streams.

Supermassive black holes are located at the center of most, if not all, massive galaxies. They follow close correlations with global properties of their host galaxies (scaling relations), and are thought to play a crucial role in galaxy evolution. Yet, we lack a complete understanding of fundamental aspects of their growth across cosmic time. In particular, we still do not understand: (1) whether black holes or their host galaxies grow faster and (2) what is the maximum mass that black holes can reach. The high angular resolution capability and sensitivity of 30-m class telescopes will revolutionize our understanding of the extreme end of the black hole and galaxy mass scale. With such facilities, we will be able to dynamically measure masses of the largest black holes and characterize galaxy properties out to redshift $z \sim 1.5$. Together with the evolution of black hole-galaxy scaling relations since $z \sim 1.5$, the maximum mass black hole will shed light on the main channels of black hole growth.

We report that LAMOST-HVS1 is a massive hyper-runaway subgiant star with mass of 8.3 Msun and super-Solar metallicity, ejected from the inner stellar disk of the Milky Way $\sim$ 33 Myr ago with the intrinsic ejection velocity of $568^{+19}_{-17}$ km/s (corrected for the streaming motion of the disk), based on the proper motion data from Gaia Data Release 2 (DR2) and high-resolution spectroscopy. The extremely large ejection velocity indicates that this star was not ejected by the supernova explosion of the binary companion. Rather, it was probably ejected by a 3- or 4-body dynamical interaction with more massive objects in a high-density environment. Such a high-density environment may be attained at the core region of a young massive cluster with mass of $\gtrsim 10^4$ Msun. The ejection agent that took part in the ejection of LAMOST-HVS1 may well be an intermediate mass black hole ($\gtrsim$ 100 Msun), a very massive star ($\gtrsim$ 100 Msun), or multiple ordinary massive stars ($\gtrsim$ 30 Msun). Based on the flight time and the ejection location of LAMOST-HVS1, we argue that its ejection agent or its natal star cluster is currently located near the Norma spiral arm. The natal star cluster of LAMOST-HVS1 may be an undiscovered young massive cluster near the Norma spiral arm.

The time irreversibility and fast relaxation of collapsing $N$-body gravitating systems (as opposed to the time reversibility of the equations of motion for individual stars or particles) are traditionally attributed to information loss due to coarse-graining in the observation. We show that this subjective element is not necessary once one takes into consideration the fundamental fact that these systems are discrete, i.e. composed of a finite number $N$ of stars or particles. We show that a connection can be made between entropy estimates for discrete systems and the Nyquist-Shannon sampling criterion. Specifically, given a sample with $N$ points in a space of $d$ dimensions, the Nyquist-Shannon criterion constrains the size of the smallest structures defined by a function in the continuum that can be uniquely associated with the discrete sample. When applied to an $N$-body system, this theorem sets a lower limit to the size of phase-space structures (in the continuum) that can be resolved in the discrete data. As a consequence, the finite $N$ system tends to a uniform distribution after a relaxation time that typically scales as $N^{1/d}$. This provides an explanation for the fast achievement of a stationary state in collapsing $N$-body gravitating systems such as galaxies and star clusters, without the need to advocate for the subjective effect of coarse-graining.

We investigate the old problem of the fast relaxation of collisionless $N$-body systems which are collapsing or perturbed, emphasizing the importance of (non-collisional) discreteness effects. We integrate orbit ensembles in fixed external potentials, estimating the entropy of the ensemble to analyze the time evolution of the distribution function. We show that these estimates capture the correct physical behavior expected from the 2nd Law of Thermodynamics, without any spurious entropy production. For self-consistent (i.e. stationary) samples, the entropy is conserved, while for non-self-consistent samples, it increases within a few dynamical times up to a maximum where it stabilizes (even in integrable potentials). Our results shed light on the main ingredients for this fast collisionless relaxation. The fundamental ingredient is the discreteness (finite $N$) of gravitational systems in any potential. Additionally, in non-integrable potentials, the presence of chaotic orbits accelerates the entropy production. Contrary to the traditional violent relaxation scenario, our results indicate that a time-dependent potential is not necessary for this fast relaxation. For the first time, in connection with the Nyquist-Shannon theorem we determine the $N$-dependence of this discreteness-driven relaxation, deriving a typical relaxation time $T/\tau_{cr}\approx 0.1 N^{1/6}$, with slightly weaker $N$-dependencies for non-integrable potentials with substantial fractions of chaotic orbits. This timescale is much smaller than the collisional timescale even for small-$N$ systems such as open clusters and represents an upper limit for real collisionless $N$-body systems. Additionally, our results reinforce the conclusion draw in Beraldo e Silva et al. (2017) that the Vlasov equation does not provide an adequate kinetic description of the fast relaxation of collapsing collisionless $N$-body systems.

We present a full dynamical model of the Fornax dwarf spheroidal galaxy obtained with the spherically symmetric Schwarzschild orbit superposition method applied to the largest kinematic data set presently available. We modelled the total mass content of the dwarf with the mass-to-light ratio $\Upsilon$ varying with radius and found that Fornax is with high probability embedded in a massive and extended dark matter halo. We estimated the total mass contained within 1$\,$kpc to be $M(<1\,$kpc$)=1.25^{+0.06}_{-0.13} \times 10^8\,$M$_{\odot}$. The data are consistent with the constant mass-to-light ratio, i.e. the mass-follows-light model, only at 3$\sigma$ level, but still require a high value of $\Upsilon \approx 11.2\,$M$_{\odot}/$L$_{\odot}$. Our results are in general agreement with previous estimates of the dynamical mass profile of Fornax. As the Schwarzschild method does not require any assumptions on the orbital anisotropy of the stars, we obtained a profile of the anisotropy parameter $\beta$ as an output of our modelling. The derived anisotropy is close to zero in the centre of the galaxy and decreases with radius, but remains consistent with isotropic orbits at all radii at 1$\sigma$ confidence level.

Gravitational 3-body interaction among binary stars and the supermassive black hole (SMBH) at the center of the Milky Way occasionally ejects a hypervelocity star (HVS) with a velocity of ~1000 km/s. Due to the ejection location, such a HVS initially has negligible azimuthal angular momentum Lz ~ 0 kpc km/s. Even if the halo is mildly triaxial, Lz of a recently ejected nearby HVS remains negligible, since its flight time from the Galactic Center is too short to accumulate noticeable torque. However, if we make a wrong assumption about the Solar position and velocity, such a HVS would apparently have noticeable non-zero azimuthal angular momentum, due to the wrong reflex motion of the Sun. Conversely, with precise astrometric data for a nearby HVS, we can measure the Solar position and velocity by assuming that the HVS has zero azimuthal angular momentum. Based on this idea, here we propose a method to estimate the Galactocentric distance of the Sun R0 and the Galactocentric Solar azimuthal velocity Vsun by using a HVS. We demonstrate with mock data for a nearby HVS candidate that the Gaia astrometric data, along with the currently available constraint on Vsun/R0 from the proper motion measurement of Sgr A*, can constrain R0 and Vsun with uncertainties of ~0.27 kpc and ~7.8 km/s (or fractional uncertainties of 3 percent), respectively. Our method will be a promising tool to constrain (R0, Vsun), given that Gaia is expected to discover many nearby HVSs in the near future.

We report the discovery of 30 stars with extreme space velocities ($>$ 480 km/s) in the Gaia-DR2 archive. These stars are a subset of 1743 stars with high-precision parallax, large tangential velocity ($v_{tan}>$ 300 km/s), and measured line-of-sight velocity in DR2. By tracing the orbits of the stars back in time, we find at least one of them is consistent with having been ejected by the supermassive black hole at the Galactic Center. Another star has an orbit that passed near the Large Magellanic Cloud (LMC) about 200 Myr ago. Unlike previously discovered blue hypervelocity stars, our sample is metal-poor (-1.5 $<$ [Fe/H] $<$ -1.0) and quite old ($>$ 1 Gyr). We discuss possible mechanisms for accelerating old stars to such extreme velocities. The high observed space density of this population, relative to potential acceleration mechanisms, implies that these stars are probably bound to the Milky Way (MW). If they are bound, the discovery of this population would require a local escape speed of around $\sim$ 600 km/s and consequently imply a virial mass of $M_{200} \sim 1.4 \times 10^{12} M_\odot$ for the MW.

We present kinematics of 35 highly r-process-enhanced ([Eu/Fe] >= +0.7) metal-poor (-3.8 < [Fe/H] < -1.4) field stars. We calculate six-dimensional positions and velocities, evaluate energies and integrals of motion, and compute orbits for each of these stars using parallaxes and proper motions from the second Gaia data release and published radial velocities. All of these stars have halo kinematics. Most stars (66%) remain in the inner regions of the halo (< 13 kpc), and many (51%) have orbits that pass within 2.6 kpc of the Galactic center. Several stars (20%) have orbits that extend beyond 20 kpc, including one with an orbital apocenter larger than the Milky Way virial radius. We apply three clustering methods to search for structure in phase space, and we identify eight groups. No abundances are considered in the clustering process, but the [Fe/H] dispersions of the groups are smaller than would be expected by random chance. The orbital properties, clustering in phase space and metallicity, and lack of highly r-process-enhanced stars on disk-like orbits indicate that such stars likely were accreted from disrupted satellites. Comparison with the galaxy luminosity-metallicity relation suggests M_V > -9 for most of the progenitor satellites, characteristic of ultra-faint or low-luminosity classical dwarf spheroidal galaxies. Environments with low rates of star formation and Fe production, rather than the nature of the r-process site, may be key to obtaining the [Eu/Fe] ratios found in highly r-process-enhanced stars.

In our previous work we confirmed the reliability of the spherically symmetric Schwarzschild orbit-superposition method to recover the mass and velocity anisotropy profiles of spherical dwarf galaxies. Here we investigate the effect of its application to intrinsically non-spherical objects. For this purpose we use a model of a dwarf spheroidal galaxy formed in a numerical simulation of a major merger of two disky dwarfs. The shape of the stellar component of the merger remnant is axisymmetric and prolate which allows us to identify and measure the bias caused by observing the spheroidal galaxy along different directions, especially the longest and shortest principal axis. The modelling is based on mock data generated from the remnant that are observationally available for dwarfs: projected positions and line-of-sight velocities of the stars. In order to obtain a reliable tool while keeping the number of parameters low we parametrize the total mass distribution as a radius-dependent mass-to-light ratio with just two free parameters we aim to constrain. Our study shows that if the total density profile is known, the true, radially increasing anisotropy profile can be well recovered for the observations along the longest axis whereas the data along the shortest axis lead to the inference an incorrect, isotropic model. On the other hand, if the density profile is derived from the method as well, the anisotropy is always underestimated but the total mass profile is well recovered for the data along the shortest axis whereas for the longest axis the mass content is overestimated.

The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We present a new study concerning the application of the Schwarzschild orbit superposition method to model spherical galaxies. The method aims to recover the mass and the orbit anisotropy parameter profiles of the objects using measurements of positions and line-of-sight velocities usually available for resolved stellar populations of dwarf galaxies in the Local Group. To test the reliability of the method, we used different sets of mock data extracted from four numerical realizations of dark matter haloes. The models shared the same density profile but differed in anisotropy profiles, covering a wide range of possibilities, from constant to increasing and decreasing with radius. The tests were done in two steps, first assuming that the mass profile of the dwarf is known and employing the method to retrieve the anisotropy only, and then varying also the mass distribution. We used two kinds of data samples: unrealistically large ones based on over 270 000 particles from the numerical realizations and small ones matching the amount of data available for the Fornax dwarf. For the large data samples we recover both the mass and the anisotropy profiles with very high accuracy. For the realistically small ones we also find a reasonably good agreement between the fitted and the input anisotropies, however the total density profiles can be significantly biased as a result of their oversensitivity to the available data. Our results therefore provide convincing evidence in favour of the applicability of the Schwarzschild method to break the mass-anisotropy degeneracy in dwarf galaxies.

The Milky Way bulge shows a box/peanut or X-shaped bulge (hereafter BP/X) when viewed in infrared or microwave bands. We examine orbits in an N-body model of a barred disk galaxy that is scaled to match the kinematics of the Milky Way (MW) bulge. We generate maps of projected stellar surface density, unsharp masked images, 3D excess-mass distributions (showing mass outside ellipsoids), line-of-sight number count distributions, and 2D line-of-sight kinematics for the simulation as well as co-added orbit families, in order to identify the orbits primarily responsible for the BP/X shape. We estimate that between 19-23\% of the mass of the bar is associated with the BP/X shape and that most bar orbits contribute to this shape which is clearly seen in projected surface density maps and 3D excess mass for non-resonant box orbits, "banana" orbits, "fish/pretzel" orbits and "brezel" orbits. {We find that nearly all bar orbit families contribute some mass to the 3D BP/X-shape. All co-added orbit families show a bifurcation in stellar number count distribution with heliocentric distance that resembles the bifurcation observed in red clump stars in the MW. However, only the box orbit family shows an increasing separation of peaks with increasing galactic latitude $|b|$, similar to that observed.} Our analysis shows that no single orbit family fully explains all the observed features associated with the MW's BP/X shaped bulge, but collectively the non-resonant boxes and various resonant boxlet orbits contribute at different distances from the center to produce this feature. We propose that since box orbits have three incommensurable orbital fundamental frequencies, their 3-dimensional shapes are highly flexible and, like Lissajous figures, this family of orbits is most easily able to adapt to evolution in the shape of the underlying potential.

Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(\sigma_r, \sigma_\theta, \sigma_\phi)$, or by the velocity anisotropy parameter $\beta=1-(\sigma_\theta^2+\sigma_\phi^2)/(2\sigma_r^2)$. Direct measurements of $(\sigma_r, \sigma_\theta, \sigma_\phi)$ consistently suggest $\beta =0.5$-$0.7$ for nearby halo stars. In contrast, the value of $\beta$ at large Galactocentric radius $r$ is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate $\beta$ for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply $\beta<0$ at $r \gtrsim 20 \;\rm{kpc}$, which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalogue analyses to show that the estimates of $\beta$ based on local fitting methods are reliable only at $r \leq 15 \;\rm{kpc}$ with the current sample size ($\sim10^3$ stars at a given radius). As $r$ increases, the line-of-sight velocity (corrected for the Solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so that it becomes increasingly more difficult to estimate tangential velocity dispersion $(\sigma_\theta, \sigma_\phi)$ from line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at $r \geq 20\;\rm{kpc}$.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.

An application of active contours without edges is presented as an efficient and effective means of extracting and characterizing coronal holes. Coronal holes are regions of low-density plasma on the Sun with open magnetic field lines. As the source of the fast solar wind, the detection and characterization of these regions is important for both testing theories of their formation and evolution and from a space weather perspective. Coronal holes are detected in full disk extreme ultraviolet (EUV) images of the corona obtained with the Solar Dynamics Observatory Atmospheric Imaging Assembly (SDO/AIA). The proposed method detects coronal boundaries without determining any fixed intensity value in the data. Instead, the active contour segmentation employs an energy-minimization in which coronal holes are assumed to have more homogeneous intensities than surrounding active regions and quiet Sun. The segmented coronal holes tend to correspond to unipolar magnetic regions, are consistent with concurrent solar wind observations, and qualitatively match the coronal holes segmented by other methods. The means to identify a coronal hole without specification of a final intensity threshold may allow this algorithm to be more robust across multiple datasets, regardless of data type, resolution, and quality.

Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. Streams even in weakly chaotic regions look very different from those in regular regions. We find that streams can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be $>$1000 orbital periods ($>$1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.

Dark Stars are stellar objects made (almost entirely) of hydrogen and helium, but powered by the heat from Dark Matter annihilation, rather than by fusion. They are in hydrostatic and thermal equilibrium, but with an unusual power source. Weakly Interacting Massive Particles (WIMPs), among the best candidates for dark matter, can be their own antimatter and can annihilate inside the star, thereby providing a heat source. Although dark matter constitutes only $\lesssim 0.1\%$ of the stellar mass, this amount is sufficient to power the star for millions to billions of years. Thus, the first phase of stellar evolution in the history of the Universe may have been dark stars. We review how dark stars come into existence, how they grow as long as dark matter fuel persists, and their stellar structure and evolution. The studies were done in two different ways, first assuming polytropic interiors and more recently using the MESA stellar evolution code; the basic results are the same. Dark stars are giant, puffy ($\sim$ 10 AU) and cool (surface temperatures $\sim$10,000 K) objects. We follow the evolution of dark stars from their inception at $\sim 1 M_\odot$ as they accrete mass from their surroundings to become supermassive stars, some even reaching masses $> 10^6 M_\odot$ and luminosities $>10^{10} L_\odot$, making them detectable with the upcoming James Webb Space Telescope. Once the dark matter runs out and the dark star dies, it may collapse to a black hole; thus dark stars may provide seeds for the supermassive black holes observed throughout the Universe and at early times. Other sites for dark star formation may exist in the Universe today in regions of high dark matter density such as the centers of galaxies. The current review briefly discusses dark stars existing today, but focuses on the early generation of dark stars.

Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration - metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (> 50%) of the stars in the outer disk; they populate the high metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[alpha/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2+/-1.2 kpc) and time (8.8+/-0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our Milky Way.

The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of r-process material in Ret 2 (mean [Eu/Fe]=+1.69+/-0.05) found in three of these stars (mean [Fe/H]=-2.88+/-0.10). The abundances closely match the r-process pattern found in the well-studied metal-poor halo star CS22892-052. Such r-process-enhanced stars have not been found in any other ultra-faint dwarf galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 ([Fe/H]=-3.42+/-0.20) contains only trace amounts of Sr ([Sr/Fe]=-1.73+/-0.43) and no detectable heavier elements. One r-process enhanced star is also enhanced in C (natal [C/Fe]=+1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and r-process enhancements are decoupled. The r-process-deficient star is enhanced in Mg ([Mg/Fe]=+0.81+/-0.14), and the other three stars show normal levels of alpha-enhancement (mean [Mg/Fe]=+0.34+/-0.03). The abundances of other alpha and Fe-group elements closely resemble those in ultra-faint dwarf galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large r-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae.

We apply the Schwarzschild orbit superposition method to mock data in order to investigate the accuracy of recovering the profile of the orbit anisotropy. The mock data come from four numerical realizations of dark matter haloes with well defined anisotropy profiles. We show that when assuming a correct mass distribution we are able to determine the anisotropy with high precision and clearly distinguish between the models.

We examine a large random sample of orbits in self-consistent simulations of N-body bars. Orbits in the bars are classified both visually and with a new automated orbit classification method based on frequency analysis. The well known prograde x1 orbit family originates from the same parent orbit as the box orbits in stationary and rotating triaxial ellipsoids. However only a small fraction of bar orbits ~4% have predominately prograde motion like their periodic parent orbit. Most bar orbits arising from the x1 orbit have little net angular momentum in the bar frame making them equivalent to box orbits in rotating triaxial potentials. A small fraction of bar orbits (~7%) are long axis tubes that behave exactly like those in triaxial ellipsoids:they are tipped about the intermediate-axis due to the Coriolis force, with the sense of tipping determined by the sign of their angular momentum about the long axis. No orbits parented by prograde periodic x2 orbits are found in the pure bar model, but a tiny population (~2%) of short axis tube orbits parented by retrograde x4 orbits are found. When a central point mass representing a supermassive black hole (SMBH) is grown adiabatically at the center of the bar, those orbits that lie in the immediate vicinity of the SMBH are transformed into precessing Keplerian orbits (PKOs) which belong the same major families (short axis tubes, long axis tubes and boxes) occupying the bar at larger radii. During the growth of a SMBH the inflow of mass and outward transport of angular momentum transforms some x1 and long axis tube orbits into prograde short axis tubes. This study has important implications for future attempts to constrain the masses of SMBHs in barred galaxies using orbit based methods like the Schwarzschild orbit superposition scheme and for understanding the observed features in barred galaxies.

We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This sub- structure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified IMF and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, on the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a 'sawtooth' shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically-increasing AMR.

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

Using N-body simulations, we study the properties of a bar induced in a dwarf galaxy as a result of tidal interaction with the Milky Way. The dwarf is initially composed of a disk embedded in a dark matter halo and we follow its evolution on a typical orbit for 10 Gyr. It undergoes an evolution typical of tidally stirred dwarf galaxies: it loses mass, the stellar component transforms from a disk to a spheroid and the rotation of the stars is partially replaced by random motions. A tidally induced bar forms at the first pericentre passage and survives until the end of the evolution. Fourier decomposition of the distribution of stars reveals that only even modes are significant and preserve a hierarchy so that the bar mode is always the strongest. They show a characteristic profile with a maximum, similar to simulated bars forming in isolated galaxies and observed bars in real galaxies. We adopt the maximum of the bar mode as a measure of the bar strength and we estimate the bar length by comparing the density profiles along the bar and perpendicular to it. Both the bar strength and the bar length decrease with time, mainly at pericentres, as a result of tidal torques acting at those times and not to secular evolution. The pattern speed of the bar varies significantly on a time scale of 1 Gyr and seems to be fully controlled by the orientation of the tidal torque from the Milky Way. The bar is never tidally locked but we discover a hint of a 5/2 orbital resonance between the third and fourth pericentre passage. The speed of the bar expressed in terms of the ratio of the bar length to the corotation radius decreases in the long run so that the bar changes from initially rather fast to slow in the later stages. The boxy/peanut shape is clearly visible in the edge-on view between the second and third pericentre and its occurrence is preceded by a short period of buckling instability.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

We present a revised measurement of the mass of the central black hole (Mbh) in the Seyfert 1 galaxy NGC 4151. The new stellar dynamical mass measurement is derived by applying an axisymmetric orbit-superposition code to near-infrared integral field data obtained using adaptive optics with the Gemini NIFS spectrograph. When our models attempt to fit both the NIFS kinematics and additional low spatial resolution kinematics, our results depend sensitively on how chi-squared is computed--probably a consequence of complex bar kinematics that manifest immediately outside the nuclear region. The most robust results are obtained when only the high spatial resolution kinematic constraints in the nuclear region are included in the fit. Our best estimates for the BH mass and H-band mass-to-light ratio are Mbh~(3.76+/-1.15)E7 Msun (1-sigma error) and M/L(H-band)~0.34+/-0.03 Msun/Lsun (3-sigma error), respectively (the quoted errors reflect the model uncertainties). Our BH mass measurement is consistent with estimates from both reverberation mapping (3.57[+0.45/-0.37]E7 Msun) and gas kinematics (3.0[+0.75/-2.2]E7 Msun; 1-sigma errors), and our best-fit mass-to-light ratio is consistent with the photometric estimate of M/L(H-band)=0.4+/-0.2 Msun/Lsun. The NIFS kinematics give a central bulge velocity dispersion sigma_c=116+/-3 km/s, bringing this object slightly closer to the M-sigma relation for quiescent galaxies. Although NGC 4151 is one of only a few Seyfert 1 galaxies in which it is possible to obtain a direct dynamical BH mass measurement--and thus, an independent calibration of the reverberation mapping mass scale--the complex bar kinematics makes it less than ideally suited for this purpose.

The optical light profiles of nearby early type galaxies are known to exhibit a smooth transition from nuclear light deficits to nuclear light excesses with decreasing galaxy mass, with as much as 80 per cent of the galaxies with stellar masses below 10^10 Msun hosting a massive nuclear star cluster. At the same time, while all massive galaxies are thought to harbor nuclear super-massive black holes (SMBHs), observational evidence for SMBHs is slim at the low end of the mass function. Here, we explore the environmental dependence of the nucleation fraction by comparing two homogeneous samples of nearby field vs. cluster early type galaxies with uniform Hubble Space Telescope (HST) coverage. Existing Chandra X-ray Telescope data for both samples yield complementary information on low-level accretion onto nuclear SMBHs. Specifically, we report on dual-band (F475W & F850LP) Advanced Camera for Surveys (ACS) imaging data for 28 out of the 103 field early type galaxies that compose the AMUSE-Field Chandra survey, and compare our results against the companion HST and Chandra surveys for a sample of 100 Virgo cluster early types (ACS Virgo Cluster and AMUSE-Virgo surveys, respectively). We model the two-dimensional light profiles of the field targets to identify and characterize NSCs, and find a field nucleation fraction of 26% (+17%, -11%; at the 1-sigma level), consistent with the measured Virgo nucleation fraction of 30% (+17%, -12%), across a comparable mass distribution. Coupled with the Chandra result that SMBH activity is higher for the field, our findings indicate that, since the last epoch of star formation, the funneling of gas to the nuclear regions has been inhibited more effectively for Virgo galaxies, arguably via ram pressure stripping.

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

We use collisionless $N$-body simulations to determine how the growth of a supermassive black hole (SMBH) influences the nuclear kinematics in both barred and unbarred galaxies. In the presence of a bar, the increase in the velocity dispersion $\sigma$ (within the effective radius) due to the growth of an SMBH is on average $\lesssim 10%$, whereas the increase is only $\lesssim 4%$ in an unbarred galaxy. In a barred galaxy, the increase results from a combination of three separate factors (a) orientation and inclination effects; (b) angular momentum transport by the bar that results in an increase in the central mass density; (c) an increase in the vertical and radial velocity anisotropy of stars in the vicinity of the SMBH. In contrast the growth of the SMBH in an unbarred galaxy causes the velocity distribution in the inner part of the nucleus to become less radially anisotropic. The increase in $\sigma$ following the growth of the SMBH is insensitive to a variation of a factor of 10 in the final mass of the SMBH, showing that it is the growth process rather than the actual SMBH mass that alters bar evolution in a way that increases $\sigma$. We argue that using an axisymmetric stellar dynamical modeling code to measure SMBH masses in barred galaxies could result in a slight overestimate of the derived $M_{BH}$, especially if a constant M/L ratio is assumed. We conclude that the growth of a black hole in the presence of a bar could result in an increase in $\sigma$ which is roughly of $4-8%$ larger than the increase that occurs in an axisymmetric system. While the increase in $\sigma$ due to SMBH growth in a barred galaxy might partially account for the claimed offset of barred galaxies and pseudo bulges from the $M_{BH}-\sigma$ relation obtained for elliptical galaxies and classical bulges in unbarred galaxies, it is inadequate to account for all of the offset.

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.

The relations observed between supermassive black holes and their host galaxies suggest a fundamental link in the processes that cause these two objects to evolve. A more comprehensive understanding of these relations could be gained by increasing the number of supermassive black hole mass (M) measurements. This can be achieved, in part, by continuing to model the stellar dynamics at the centers of galactic bulges using data of the highest possible spatial resolution. Consequently, we present here an atlas of galaxies in the Space Telescope Imaging Spectrograph (STIS) data archive that may have spectra suitable for new M estimates. Archived STIS G750M data for all non-barred galactic bulges are co-aligned and combined, where appropriate, and the radial signal-to-noise ratios calculated. The line-of-sight velocity distributions from the CaII triplet are then determined using a maximum penalized likelihood method. We find 19 out of 42 galaxies may provide useful new M estimates since they are found to have data that is comparable in quality with data that has been used in the past to estimate M. However, we find no relation between the signal-to-noise ratio in the previously analyzed spectra and the uncertainties of the black hole masses derived from the spectra. We also find that there is a very limited number of appropriately observed stellar templates in the archive from which to estimate the effects of template mismatching.

Models of the Sagittarius Stream have consistently found that the Milky Way disc is oriented such that its short axis is along the intermediate axis of the triaxial dark matter halo. We attempt to build models of disc galaxies in such an `intermediate-axis orientation'. We do this with three models. In the first two cases we simply rigidly grow a disc in a triaxial halo such that the disc ends up perpendicular to the global intermediate axis. We also attempt to coax a disc to form in an intermediate-axis orientation by producing a gas+dark matter triaxial system with gas angular momentum about the intermediate axis. In all cases we fail to produce systems which remain with stellar angular momentum aligned with the halo's intermediate axis, even when the disc's potential flattens the inner halo such that the disc is everywhere perpendicular to the halo's local minor axis. For one of these unstable simulations we show that the potential is even rounder than the models of the Milky Way potential in the region probed by the Sagittarius Stream. We conclude that the Milky Way's disc is very unlikely to be in an intermediate axis orientation. However we find that a disc can persist off one of the principal planes of the potential. We propose that the disc of the Milky Way must be tilted relative to the principal axes of the dark matter halo. Direct confirmation of this prediction would constitute a critical test of Modified Newtonian Dynamics.

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

We present velocity-resolved reverberation results for five active galactic nuclei. We recovered velocity-delay maps using the maximum-entropy method for four objects: Mrk 335, Mrk 1501, 3C120, and PG2130+099. For the fifth, Mrk 6, we were only able to measure mean time delays in different velocity bins of the Hbeta emission line. The four velocity-delay maps show unique dynamical signatures for each object. For 3C120, the Balmer lines show kinematic signatures consistent with both an inclined disk and infalling gas, but the HeII 4686 emission line is suggestive only of inflow. The Balmer lines in Mrk 335, Mrk 1501, and PG 2130+099 show signs of infalling gas, but the HeII emission in Mrk 335 is consistent with an inclined disk. We also see tentative evidence of combined virial motion and infalling gas from the velocity-binned analysis of Mrk 6. The maps for 3C120 and Mrk 335 are two of the most clearly defined velocity-delay maps to date. These maps constitute a large increase in the number of objects for which we have resolved velocity-delay maps and provide evidence supporting the reliability of reverberation-based black hole mass measurements.

We present the results from a detailed analysis of photometric and spectrophotometric data on five Seyfert 1 galaxies observed as a part of a recent reverberation mapping program. The data were collected at several observatories over a 140-day span beginning in 2010 August and ending in 2011 January. We obtained high sampling-rate light curves for Mrk 335, Mrk 1501, 3C120, Mrk 6, and PG2130+099, from which we have measured the time lag between variations in the 5100 Angstrom continuum and the H-beta broad emission line. We then used these measurements to calculate the mass of the supermassive black hole at the center of each of these galaxies. Our new measurements substantially improve previous measurements of MBH and the size of the broad line-emitting region for four sources and add a measurement for one new object. Our new measurements are consistent with photoionization physics regulating the location of the broad line region in active galactic nuclei.

We present the first results from a detailed analysis of photometric and spectrophotometric data on the narrow-line Seyfert 1 galaxy Mrk 335, collected over a 120-day span in the fall of 2010. From these data we measure the lag in the He II 4686 broad emission line relative to the optical continuum to be 2.7 \pm 0.6 days and the lag in the H\beta 4861 broad emission line to be 13.9 \pm 0.9 days. Combined with the line width, the He II lag yields a black hole mass, MBH = (2.6 \pm 0.8)\times 10^7 Msun. This measurement is consistent with measurements made using the H\beta 4861 line, suggesting that the He II emission originates in the same region as H\beta, but at a much smaller radius. This constitutes the first robust lag measurement for a high-ionization line in a narrow-line Seyfert 1 galaxy.

We study the capability of the James Webb Space Telescope (JWST) to detect Supermassive Dark Stars (SMDS). If the first stars are powered by dark matter heating in triaxial dark matter haloes, they may grow to be very large and very bright, visible in deep imaging with JWST and even Hubble Space Telescope (HST). We use HST surveys to place bounds on the numbers of SMDSs that may be detected in future JWST imaging surveys. We showed that SMDS in the mass range $10^6-10^7 M_\odot$ are bright enough to be detected in all the wavelength bands of the NIRCam on JWST . If SMDSs exist at z ~10, 12, and 14, they will be detectable as J-band, H-band, or K-band dropouts, respectively. With a total survey area of 150 arcmin^2 (assuming a multi-year deep parallel survey with JWST), we find that typically the number of $10^6 M_\odot$ SMDSs found as H or K-band dropouts is ~10^5\fsmds, where the fraction of early DM haloes hosting DS is likely to be small, \fsmds<<1. If the SDMS survive down to z=10 where HST bounds apply, then the observable number of SMDSs as H or K-band dropouts with JWST is ~1-30. While individual SMDS are bright enough to be detected by JWST, standard PopIII stars are not, and would only be detected in first galaxies with total stellar masses of ~$10^6-10^8 M_\odot$. Differentiating first galaxies at z>10 from SMDSs would be possible with spectroscopy: the SMDS (which are too cool produce significant nebular emission) will have only absorption lines while the galaxies are likely to produce emission lines as well. Of particular interest would be the 1640 HeII emission line as well as H{\alpha} lines which would be signatures of early galaxies rather than SMDSs. The detection of SMDSs would not only provide alternative evidence for WIMPs but would also provide possible seeds for the formation of supermassive black holes that power QSOs at z~6.

Resolved surveys of the Milky Way's stellar halo can obtain all 6 phase space coordinates of tens of thousands of individual stars, making it possible to compute their 3-dimensional orbits. Spectral analysis of large numbers of halo orbits can be used to construct frequency maps which are a compact, yet informative representation of their phase space distribution function (DF). Such maps can be used to infer the major types of orbit families that constitute the DF of stellar halo and their relative abundances. The structure of the frequency maps, especially the resonant orbits, reflects the formation history and shape of the dark matter potential and its orientation relative to the disk. The application of frequency analysis to cosmological hydrodynamic simulations of disk galaxies shows that the orbital families occupied by halo stars and dark matter particles are very similar, implying that stellar halo orbits can be used to constrain the DF of the dark matter halo, possibly impacting future direct dark matter detection experiments. An application of these methods to a sample of \sim 16,000 Milky Way halo and thick disk stars from the SDSS-SEGUE survey yields a frequency map with strong evidence for resonant trapping of halo stars by the Milky Way disk, in a manner predicted by controlled simulations in which the disk grows adiabatically. The application of frequency analysis methods to current and future phase space data for Milky Way halo stars will provide new insights into the formation history of the dierent components of the Galaxy and the DF of the halo.

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function (ABRIDGED).

We survey the properties of all orbit families in the rotating frame of a family of realistic triaxial potentials with central supermassive black holes (SMBHs). In such galaxies, most regular box orbits (vital for maintaining triaxiality) are associated with resonances which occupy two-dimensional surfaces in configuration space. For slow figure rotation all orbit families are largely stable. At intermediate pattern speeds a significant fraction of the resonant box orbits as well as inner long-axis tubes are destabilized by the "envelope doubling" that arises from the Coriolis forces and are driven into the destabilizing center. Thus, for pattern rotation periods .2 Gyr < Tp < 5 Gyr, the two orbit families that are most important for maintaining triaxiality are highly chaotic. As pattern speed increases there is also a sharp decrease in the overall fraction of prograde short-axis tubes and a corresponding increase in the retrograde variety. At the highest pattern speeds (close to that of triaxial bars), box-like orbits undergo a sudden transition to a new family of stable retrograde loop-like orbits, which resemble orbits in three-dimensional bars, and circulate about the short axis. Our analysis implies that triaxial systems (with central cusps and SMBHs) can either have high pattern speeds like fast bars or low patten speeds like triaxial elliptical galaxies or dark matter halos found in N-body simulations. Intermediate pattern speeds produce a high level of stochasticity in both the box and inner long-axis tube orbit families implying that stable triaxial systems are unlikely to have such pattern speeds.

The first stars to form in the history of the universe may have been powered by dark matter annihilation rather than by fusion. This new phase of stellar evolution may have lasted millions to billions of years. These dark stars can grow to be very large, > 10^5 solar masses, and are relatively cool (~10^4 K). They are also very bright, being potentially detectable in the upcoming James Webb Space Telescope or even the Hubble Space Telescope. Once the dark matter runs out, the dark stars have a short fusion phase, before collapsing into black holes (BH). The resulting BH could serve as seeds for the (unexplained) supermassive black holes at high redshift and at the centers of galaxies.

The first phase of stellar evolution in the history of the Universe may be Dark Stars, powered by dark matter heating rather than by nuclear fusion. Weakly Interacting Massive Particles, which may be their own antipartners, collect inside the first stars and annihilate to produce a heat source that can power the stars for millions to billions of years. In this paper we show that these objects can grow to be supermassive dark stars (SMDS) with masses $\gtrsim (10^5-10^7) \msun$. The growth continues as long as dark matter heating persists, since dark stars are large and cool (surface temperature $\lesssim 5\times 10^4$K) and do not emit enough ionizing photons to prevent further accretion of baryons onto the star. The dark matter may be provided by two mechanisms: (1) gravitational attraction of dark matter particles on a variety of orbits not previously considered, and (2) capture of WIMPs due to elastic scattering. Once the dark matter fuel is exhausted, the SMDS becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that may provide seeds for supermassive black holes in the Universe. SMDS are very bright, with luminosities exceeding $(10^9-10^{11}) L_\odot$. We demonstrate that for several reasonable parameters, these objects will be detectable with JWST. Such an observational discovery would confirm the existence of a new phase of stellar evolution powered by dark matter.

Using spectral methods, we analyse the orbital structure of prolate/triaxial dark matter (DM) halos in N-body simulations to understand the processes that drive the evolution of shapes of DM halos and elliptical galaxies in which central masses are grown. A longstanding issue is whether the change in the shapes of DM halos is the result of chaotic scattering of box orbits, or whether they change shape adiabatically in response to the evolving galactic potential. We use orbital frequencies to classify orbits, to quantify orbital shapes, and to identify resonant orbits and chaotic orbits. The frequency-based method overcomes the limitations of Lyapunov exponents which are sensitive to numerical discreteness effects. Regardless of the distribution of the baryonic component, the shape of a DM halo changes primarily due to changes in the shapes of individual orbits within a given family. Orbits with small pericentric radii are more likely to change both their orbital type and shape than orbits with large pericentric radii. Whether the evolution is regular (and reversible) or chaotic (and irreversible), depends primarily on the radial distribution of the baryonic component. The growth of an extended baryonic component of any shape results in a regular rather than chaotic change in orbital populations. In contrast the growth of a massive and compact central component results in chaotic scattering of a significant fraction of both box and long-axis tube orbits. The growth of a disk causes a significant fraction of halo particles to become trapped by major global orbital resonances. Despite the fact that shape of a DM halo is always quite oblate following the growth of a central baryonic component, a significant fraction of its orbit population has characteristics of its triaxial or prolate progenitor (ABRIDGED).

We summarize recent developments in the use of spectral methods for analyzing large numbers of orbits in N-body simulations to obtain insights into the global phase space structure of dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand the physical mechanism by which the shapes of dark matter halos evolve in response to the growth of central baryonic components. Halos change shape primarily because individual orbits change their shapes adiabatically in response to the growth of a baryonic component, with those at small radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for centrophilic box orbits.

We study the evolution of phase-space density during the hierarchical structure formation of LCDM halos. We compute both a spherically-averaged surrogate for phase-space density (Q) and the coarse-grained distribution function f(x,v) for dark matter particles that lie within~2 virial radii of four Milky-Way-sized dark matter halos. The estimated f(x,v) spans over four decades at any radius. Dark matter particles that end up within two virial radii of a Milky-Way-sized DM halo at $z=0$ have an approximately Gaussian distribution in log(f) at early redshifts, but the distribution becomes increasingly skewed at lower redshifts. The value corresponding to the peak of the Gaussian decreases as the evolution progresses and is well described by a power-law in (1+z). The highest values of f are found at the centers of dark matter halos and subhalos, where f can be an order of magnitude higher than in the center of the main halo. The power-law Q(r) profile likely reflects the distribution of entropy (K = sigma^2/rho^{2/3} \propto r^{1.2}), which dark matter acquires as it is accreted onto a growing halo. The estimated f(x, v), on the other hand, exhibits a more complicated behavior. Although the median coarse-grained phase-space density profile F(r) can be approximated by a power-law in the inner regions of halos and at larger radii the profile flattens significantly. This is because phase-space density averaged on small scales is sensitive to the high-f material associated with surviving subhalos, as well as relatively unmixed material (probably in streams) resulting from disrupted subhalos, which contribute a sizable fraction of matter at large radii. (ABRIDGED)

We use dissipationless N-body simulations to investigate the evolution of the true coarse-grained phase-space density distribution f(x,v) in equal-mass mergers between dark matter (DM) halos. The halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and we employ the phase-space density estimator ``Enbid'' developed by Sharma & Steinmetz to compute f(x,v). The adopted force resolution allows robust phase-space density profile estimates in the inner ~1% of the virial radii of the simulated systems. We confirm that mergers result in a decrease of the coarse-grained phase-space density in accordance with expectations from Mixing Theorems for collisionless systems. We demonstrate that binary mergers between identical DM halos produce remnants that retain excellent memories of the inner slopes and overall shapes of the phase-space density distribution of their progenitors. The robustness of the phase-space density profiles holds for a range of orbital energies, and a variety of encounter configurations including sequences of several consecutive merger events, designed to mimic hierarchical merging, and collisions occurring at different cosmological epochs. If the progenitor halos are constructed with appreciably different asymptotic power-law indices, we find that the inner slope and overall shape of the phase-space density distribution of the remnant are substantially closer to that of the initial system with the steepest central density cusp. These results explicitly demonstrate that mixing is incomplete in equal-mass mergers between DM halos, as it does not erase memory of the progenitor properties. Our results also confirm the recent analytical predictions of Dehnen (2005) regarding the preservation of merging self-gravitating central density cusps.

In the past decade, surveys of the stellar component of the Galaxy have revealed a number of streams from tidally disrupted dwarf galaxies and globular clusters. Simulations of hierarchical structure formation in LCDM cosmologies predict that the dark matter halo of a galaxy like the Milky Way contains hundreds of subhalos with masses of ~10^8 solar masses and greater, and it has been suggested that the existence of coherent tidal streams is incompatible with the expected abundance of substructure. We investigate the effects of dark matter substructure on tidal streams by simulating the disruption of a self-gravitating satellite on a wide range of orbits in different host models both with and without substructure. We find that the halo shape and the specific orbital path more strongly determine the overall degree of disruption of the satellite than does the presence or absence of substructure, i.e., the changes in the large-scale properties of the tidal debris due to substructure are small compared to variations in the debris from different orbits in a smooth potential. Substructure typically leads to an increase in the degree of clumpiness of the tidal debris in sky projection, and in some cases a more compact distribution in line-of-sight velocity. Substructure also leads to differences in the location of sections of debris compared to the results of the smooth halo model, which may have important implications for the interpretation of observed tidal streams. A unique signature of the presence of substructure in the halo which may be detectable by upcoming surveys is identified. We conclude, however, that predicted levels of substructure are consistent with a detection of a coherent tidal stream from a dwarf galaxy.

We present a stellar dynamical estimate of the black hole (BH) mass in the Seyfert 1 galaxy, NGC 4151. We analyze ground-based spectroscopy as well as imaging data from the ground and space, and we construct 3-integral axisymmetric models in order to constrain the BH mass and mass-to-light ratio. The dynamical models depend on the assumed inclination of the kinematic symmetry axis of the stellar bulge. In the case where the bulge is assumed to be viewed edge-on, the kinematical data give only an upper limit to the mass of the BH of ~4e7 M_sun (1 sigma). If the bulge kinematic axis is assumed to have the same inclination as the symmetry axis of the large-scale galaxy disk (i.e., 23 degrees relative to the line of sight), a best-fit dynamical mass between 4-5e7 M_sun is obtained. However, because of the poor quality of the fit when the bulge is assumed to be inclined (as determined by the noisiness of the chi^2 surface and its minimum value), and because we lack spectroscopic data that clearly resolves the BH sphere of influence, we consider our measurements to be tentative estimates of the dynamical BH mass. With this preliminary result, NGC 4151 is now among the small sample of galaxies in which the BH mass has been constrained from two independent techniques, and the mass values we find for both bulge inclinations are in reasonable agreement with the recent estimate from reverberation mapping (4.57[+0.57/-0.47]e7 M_sun) published by Bentz et al.

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.

Hubble Space Telescope (HST) images and spectra of the nucleated dwarf elliptical galaxy NGC 205 are combined with 3-integral axisymmetric dynamical models to constrain the mass (M_BH) of a putative nuclear black hole. This is only the second attempt, after M33, to use resolved stellar kinematics to search for a nuclear black hole with mass below 10^6 solar masses. We are unable to identify a best-fit value of M_BH in NGC 205; however, the data impose a upper limit of 2.2x10^4 M_sun (1sigma confidence) and and upper limit of 3.8x10^4 M_sun (3sigma confidence). This upper limit is consistent with the extrapolation of the M_BH-sigma relation to the M_BH < 10^6 M_sunregime. If we assume that NGC 205 and M33 both contain nuclear black holes, the upper limits on M_BH in the two galaxies imply a slope of ~5.5 or greater for the M_BH-sigma relation. We use our 3-integral models to evaluate the relaxation time (T_r) and stellar collision time (T_coll) in NGC 205; T_r~10^8 yr or less in the nucleus and T_coll~10^11 yr. The low value of T_r is consistent with core collapse having already occurred, but we are unable to draw conclusions from nuclear morphology about the presence or absence of a massive black hole.

We investigate the ability of three-integral, axisymmetric, orbit-based modeling algorithms to recover the parameters defining the gravitational potential (M/L ratio and black hole mass Mh) in spheroidal stellar systems using stellar kinematical data. We show that the potential estimation problem is generically under-determined when applied to long-slit kinematical data of the kind used in most black hole mass determinations to date. A range of parameters (M/L, Mh) can provide equally good fits to the data, making it impossible to assign best-fit values. We illustrate the indeterminacy using a variety of data sets derived from realistic models as well as published observations of the galaxy M32. In the case of M32, our reanalysis demonstrates that data published prior to 2000 are equally consistent with Mh in the range 1.5x10^6-5x10^6 solar masses, with no preferred value in that range. While the HST/STIS data for this galaxy may overcome the degeneracy in Mh, HST data for most galaxies do not resolve the black hole's sphere of influence and in these galaxies the degree of degeneracy allowed by the data may be substantial. We investigate the effect on the degeneracy of enforcing smoothness (regularization) constraints. However we find no indication that the true potential can be recovered simply by enforcing smoothness. For a given smoothing level, all solutions in the minimum-chisquare valley exhibit similar levels of noise. These experiments affirm that the indeterminacy is real and not an artifact associated with non-smooth solutions. (Abridged)

We have obtained optical long-slit spectroscopy of the nucleus of M32 using the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. The stellar rotation velocity and velocity dispersion, as well as the full line-of-sight velocity distribution (LOSVD), were determined as a function of position along the slit using two independent spectral deconvolution algorithms. We see three clear kinematical signatures of the nuclear black hole: a sudden upturn, at ~0.3 arc seconds from the center, in the stellar velocity dispersions; a flat or rising rotation curve into the center; and strong, non-Gaussian wings on the central LOSVD. The central velocity dispersion is ~130 km/s (Gaussian fit) or ~175 km/s (corrected for the wings). Both the velocity dispersion spike and the shape of the central LOSVD are consistent with the presence of a supermassive compact object in M32 with a mass in the range 2-5 x 10^6 solar masses. These data are a significant improvement on previous stellar kinematical data, making M32 the first galaxy for which the imprint of the black hole's gravitation on the stellar velocities has been observed with a resolution comparable to that of gas-dynamical studies.

The orbits of stars in galaxies are generically chaotic: the chaotic behavior arises in part from the intrinsically grainy nature of a potential that is composed of point masses. Even if the potential is assumed to be smooth, however, orbits in non-axisymmetric galaxies can be chaotic due to the presence of central density cusps or black holes. The chaotic nature of orbits implies that perturbations will grow exponentially and this in turn is expected to result in a diffusion in phase space. We show that the degree of orbital evolution is not well predicted by the growth rate of infinitesimal perturbations, i.e. by the Liapunov exponent. A more useful criterion is whether perturbations continue to grow exponentially until their scale is of order the size of the system. We illustrate these ideas in a potential consisting of N fixed point masses. Liapunov exponents are large for all values of N, but orbits become increasingly regular in their behavior as N increases; the reason is that the exponential divergence saturates at smaller and smaller distances as N is increased. The objects which impede diffusion are the invariant tori. When there are few stable tori, an ensemble of chaotic orbits evolves rapidly toward a nearly stationary state. This mixing process occurs on timescales of a few crossing times in triaxial potentials containing massive central singularities, consistent with the rapid evolution observed in N-body simulations of galaxies with central black holes.

The maximally compact representation of a regular orbit is in terms of its action-angle variables. Computing the map between a trajectory's Cartesian coordinates and its action-angle variables is called torus construction. This article reviews various approaches to torus construction and their application to galactic dynamics.

Box orbits in triaxial potentials are generically thin, that is, they lie close in phase space to a resonant orbit satisfying a relation of the form l\omega_1 +m\omega_2+n\omega_3=0 between the three fundamental frequencies. Resonant orbits are confined to a membrane; they play roughly the same role, in three dimensions, that closed orbits play in two. Stable resonant orbits avoid the center of the potential; orbits that are thick enough to pass near the center are typically stochastic. Very near the center, where the gravitational potential is dominated by the black hole, resonant orbits continue to exist, including at least one family whose elongation is parallel to the long axes of the triaxial figure.

Figure rotation substantially increases the fraction of stochastic orbits in triaxial systems. This increase is most dramatic in systems with shallow cusps showing that it is not a direct consequence of scattering by a central density cusp or black hole. In a recent study of stationary triaxial potentials (Valluri & Merritt 1998) it was found that the most important elements that define the structure of phase space are the two dimensional resonant tori. The increase in the fraction of stochastic orbits in models with figure rotation is a direct consequence of the destabilization of these resonant tori. The presence of a large fraction of stochastic orbits in a triaxial bulge will result in the evolution of its shape from triaxial to axisymmetric. The timescales for evolution can be as short as a few crossing times in the bulges of galaxies and evolution is accelerated by figure rotation. This suggests that low luminosity ellipticals and the bulges of early type spirals are likely to be predominantly axisymmetric.

We use Laskar's frequency mapping technique to study the dynamics of triaxial galaxies with central density cusps and nuclear black holes. For ensembles of 10^4 orbits, we numerically compute the three fundamental frequencies of the motion, allowing us to map out the Arnold web. We also compute diffusion rates of stochastic orbits in frequency space. The objects of fundamental importance in structuring phase space are found to be the 3-dimensional resonant tori; even when stable, such tori are not necessarily associated with periodic orbits as in systems with only two degrees of freedom. Boxlike orbits are generically stochastic, but some tube orbits are stochastic as well. The spectrum of diffusion rates for box-like orbits at a given energy is well approximated as a power law over at least six decades. Models with high central concentrations -- steep central cusps or massive black holes -- exhibit the most stochasticity. A black hole with a mass of 0.3% the mass of the galaxy is as effective as the steepest central density cusp at inducing stochastic diffusion. There is a transition to global stochasticity in box-like phase space when the mass of a central black hole exceeds 2% the galaxy mass. We predict a greater average degree of dynamical evolution in faint ellipticals, due to their high central densities and short crossing times. The evolution time is estimated to be shorter than a galaxy lifetime for absolute magnitudes fainter than about -19 or -20, consistent with the observed change in many elliptical galaxy properties at this luminosity.

The motion of stars in the gravitational potential of a triaxial galaxy is generically chaotic. However, the timescale over which the chaos manifests itself in the orbital motion is a strong function of the degree of central concentration of the galaxy. Here, chaotic diffusion rates are presented for orbits in triaxial models with a range of central density slopes and nuclear black-hole masses. Typical diffusion times are found to be less than a galaxy lifetime in triaxial models where the density increases more rapidly than 1/r at the center, or which contain black holes with masses that exceed roughly 0.1% of the galaxy mass. When the mass of a central black hole exceeds roughly 0.02 times the mass of the galaxy, there is a transition to global stochasticity and the galaxy evolves to an axisymmetric shape in little more than a crossing time. This rapid evolution may provide a negative feedback mechanism that limits the mass of nuclear black holes to a few percent of the stellar mass of a galaxy.

We have examined the stability of a sequence of oblate elliptical galaxy models having the Stackel form suggested by Kuz'min & Kutuzov. We have employed the 2-integral DFs given by Dejonghe & de Zeeuw for which flattened non-rotating models are characterized by counter-streaming motion and are radially cool; we introduce net rotation in some models by changing the sign of the z-component of angular momentum for a fraction of the particles. We have found that all non-rotating and slowly rotating members of this sequence rounder than E7 are stable, and that even maximally rotating models rounder than E4 are stable. In the absence of strong rotation, the most disruptive instability, and the last to be stabilized by increasing thickness, is a lopsided (m=1) mode. This instability appears to be driven by counter-rotation in radially cool models. Its vigour is lessened as rotation is increased, but it remains strong even in models with net angular momentum 90% of that of a maximally rotating model before finally disappearing in maximally rotating models. Strongly rotating models are more unstable to bar-forming modes which afflict maximally rotating models with c/a >~ 0.5, but this mode is quickly stabilized by moderate fractions of counter-rotating particles. Bending instabilities appear not to be very important; they are detectable in the inner parts of the flatter models, but are less vigorous and more easily stabilized than the lop-sided or bar modes in every case. We briefly discuss the possible relevance of the lop-sided instability to the existence of many lop-sided disk galaxies.

We report results of a search for optical emission line gas in the Hickson compact group 62. From narrow band (H-alpha) imaging observations we estimate 10^4 Msun of ionized hydrogen gas within the central regions of three of the galaxies in the group. There has been considerable interest in this group since ROSAT detected an extended halo of diffuse X-ray emission, possibly with a cooling flow at the center. The optical emission line luminosities from the three galaxies are comparable to each other (~10^{39} ergs/s) and resemble those seen in normal X-ray luminous early type galaxies. There is no evidence for any excess emission in the galaxy which lies at the center of the X-ray cooling flow. We find that the most likely source of excitation for the emission line nebulosities is photoionization by young stars.

We investigate the timescales for stochasticity and chaotic mixing in a family of triaxial potentials that mimic the distribution of light in elliptical galaxies. Some of the models include central point masses designed to represent nuclear black holes. Most of the boxlike orbits are found to be stochastic, with mean Liapunov times that are 3-6 times the period of the long-axis orbit. In models with large cores or small black holes, the stochastic orbits mimic regular box orbits for hundreds of oscillations at least. However a small core radius or significant black hole mass causes most of the stochastic orbits to diffuse through phase space on the same timescale, visiting a significant fraction of the volume beneath the equipotential surface. We estimate timescales for chaotic mixing in the more strongly stochastic models by evolving ensembles of 10^4 points until their distribution reaches a nearly steady state. Mixing initially takes place rapidly, with characteristic times of 10-30 dynamical times, as the phase points fill a region similar in shape to that of a box orbit. Subsequent mixing is slower, with characteristic times of hundreds of orbital times. Mixing rates were found to be enhanced by the addition of modest force perturbations. The consequences for the structure and evolution of elliptical galaxies are discussed.