We interpret simulations of secularly-evolving disc galaxies through orbit morphology. Using a new algorithm that measures the volume of orbits in real space using a tessellation, we rapidly isolate commensurate (resonant) orbits. We identify phase-space regions occupied by different orbital families. Compared to spectral methods, the tessellation algorithm can identify resonant orbits within a few dynamical periods, crucial for understanding an evolving galaxy model. The flexible methodology accepts arbitrary potentials, enabling detailed descriptions of the orbital families. We apply the machinery to four different potential models, including two barred models, and fully characterise the orbital membership. We identify key differences in the content of orbit families, emphasizing the presence of orbit families indicative of the bar evolutionary state and the shape of the dark matter halo. We use the characterisation of orbits to investigate the shortcomings of analytic and self-consistent studies, comparing our findings to the evolutionary epochs in self-consistent barred galaxy simulations. Using insight from our orbit analysis, we present a new observational metric that uses spatial and kinematic information from integral field spectrometers that may reveal signatures of commensurabilities and allow for a differentiation between models.

We track the angular momentum transfer in n-body simulations of barred galaxies by measuring torques to understand the dynamical mechanisms responsible for the evolution of the bar-disc-dark matter halo system. We find evidence for three distinct phases of barred galaxy evolution: assembly, secular growth, and steady-state equilibrium. Using a decomposition of the disc into orbital families, we track bar mass and angular momentum through time and correlate the quantities with the phases of evolution. We follow the angular momentum transfer between particles and identify the dominant torque channels. We find that the halo model mediates the assembly and growth of the bar for a high central density halo, and the outer disc mediates the assembly and growth of the bar in a low central density halo model. Both galaxies exhibit a steady-state equilibrium phase where the bar is neither lengthening nor slowing. The steady-state equilibrium results from the balance of torque between particles that are gaining and losing angular momentum. We propose observational metrics for barred galaxies that can be used to help determine the evolutionary phase of a barred galaxy, and discuss the implications of the phases for galaxy evolution as a whole.

We study the mechanisms and evolutionary phases of bar formation in n-body simulations of a stellar disc and dark matter halo system using harmonic basis function expansion analysis to characterize the dynamical mechanisms in bar evolution. We correlate orbit families with phases of bar evolution by using empirical orthogonal functions that act as a spatial filter and form the gravitational potential basis. In both models we find evidence for three phases in evolution with unique harmonic signatures. We recover known analytic results, such as bar slowdown owing to angular momentum transfer. We also find new dynamical mechanisms for bar evolution: a steady-state equilibrium configuration and harmonic interaction resulting in harmonic mode locking, both of which may be observable. Additionally, we find that ellipse fitting may severely overestimate measurements of bar length by a factor of two relative to the measurements based on orbits that comprise the true backbone supporting the bar feature. The bias will lead to overestimates of both bar mass and bar pattern speed, affecting inferences about the evolution of bars in the real universe, such as the fraction of bars with fast pattern speeds. We propose a direct observational technique to compute the radial extent of trapped orbits and determine a dynamical length for the bar.

We explore the application of Bayesian image analysis to infer the properties of an SDSS early-type galaxy sample including AGN. We use GALPHAT (Yoon et al. 2010) with a Bayes-factor model comparison to photometrically infer an AGN population and verify this using spectroscopic signatures. Our combined posterior sample for the SDSS sample reveals distinct low and high concentration modes after the point-source flux is modeled. This suggests that ETG parameters are intrinsically bimodal. The bimodal signature was weak when analyzed by GALFIT (Peng et al. 2002, 2010). This led us to create several ensembles of synthetic images to investigate the bias of inferred structural parameters and compare with GALFIT. GALPHAT inferences are less biased, especially for high-concentration profiles: GALPHAT S\'ersic index $n$, $r_{e}$ and MAG deviate from the true values by $6\%$, $7.6\%$ and $-0.03 \,\mathrm{mag}$, respectively, while GALFIT deviates by $15\%$, $22\%$ and $-0.09$\, mag, respectively. In addition, we explore the reliability for the photometric detection of AGN using Bayes factors. For our SDSS sample with $r_{e}\ge 7.92\,$arcsec, we correctly identify central point sources with $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 5$ for $n\le6$ and $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 3$ for $n>6$. The magnitude range increases and classification error decreases with increasing resolution, suggesting that this approach will excel for upcoming high-resolution surveys. Future work will extend this to models that test hypotheses of galaxy evolution through the cosmic time.

We investigate the complex interactions between the stellar disc and the dark-matter halo during bar formation and evolution using N-body simulations with fine temporal resolution and optimally chosen spatial resolution. We find that the forming stellar bar traps dark matter in the vicinity of the stellar bar into bar-supporting orbits. We call this feature the shadow bar. The shadow bar modifies both the location and magnitude of the angular momentum transfer between the disc and dark matter halo and adds 10 per cent to the mass of the stellar bar over 4 Gyr. The shadow bar is potentially observable by its density and velocity signature in spheroid stars and by direct dark matter detection experiments. Numerical tests demonstrate that the shadow bar can diminish the rate of angular momentum transport from the bar to the dark matter halo by more than a factor of three over the rate predicted by dynamical friction with an untrapped dark halo, and thus provides a possible physical explanation for the observed prevalence of fast bars in nature.

Over a handful of rotation periods, dynamical processes in barred galaxies induce non-axisymmetric structure in dark matter halos. Using n-body simulations of a Milky Way-like barred galaxy, we identify both a trapped dark-matter component, a shadow bar, and a strong response wake in the dark-matter distribution that affects the predicted dark-matter detection rates for current experiments. The presence of a baryonic disk together with well-known dynamical processes (e.g. spiral structure and bar instabilities) increase the dark matter density in the disk plane. We find that the magnitude of the combined stellar and shadow bar evolution, when isolated from the effect of the axisymmetric gravitational potential of the disk, accounts for >30% of this overall increase in disk-plane density. This is significantly larger that of previously claimed deviations from the standard halo model. The dark-matter density and kinematic wakes driven by the Milky Way bar increase the detectability of dark matter overall, especially for the experiments with higher $v_{min}$. These astrophysical features increase the detection rate by more than a factor of two when compared to the standard halo model and by a factor of ten for experiments with high minimum recoil energy thresholds. These same features increase (decrease) the annual modulation for low (high) minimum recoil energy experiments. We present physical arguments for why these dynamics are generic for barred galaxies such as the Milky Way rather than contingent on a specific galaxy model.

Our understanding of the mechanisms governing the structure and secular evolution galaxies assume nearly integrable Hamiltonians with regular orbits; our perturbation theories are founded on the averaging theorem for isolated resonances. On the other hand, it is well-known that dynamical systems with many degrees of freedom are irregular in all but special cases. The best developed framework for studying the breakdown of regularity and the onset is the Kolmogorov-Arnold-Moser (KAM) theory. Here, we use a numerical version of the KAM procedure to construct regular orbits (tori) and locate irregular orbits (broken tori). Irregular orbits are most often classified in astronomical dynamics by their exponential divergence using Lyapunov exponents. Although their computation is numerically challenging, the procedure is straightforward and they are often used to estimate the measure of regularity. The numerical KAM approach has several advantages: 1) it provides the morphology of perturbed orbits; 2) its constructive nature allows the tori to be used as basis for studying secular evolution; 3) for broken tori, clues to the cause of the irregularity may be found by studying the largest, diverging Fourier terms; and 4) it is more likely to detect weak chaos and orbits close to bifurcation. Conversely, it is not a general technique and works most cleanly for small perturbations. We develop a perturbation theory that includes chaos by retaining an arbitrary number of interacting terms rather than eliminating all but one using the averaging theorem. The companion papers show that models with significant stochasticity seem to be the rule, not the exception.

Owing to the pioneering work of Contopoulos, a strongly barred galaxy is known to have irregular orbits in the vicinity of the bar. By definition, irregular orbits can not be represented by action-angle tori everywhere in phase space. This thwarts perturbation theory and complicates our understanding of their role in galaxy structure and evolution. This paper provides a qualitative introduction to a new method based on KAM theory for investigating the morphology of regular and irregular orbits based on direct computation of tori described in Paper 1 and applies it to a galaxy disc bar. Using this method, we find that much of the phase space inside of the bar radius becomes chaotic for strong bars, excepting a small region in phase space between the ILR and corotation resonances for orbits of moderate ellipticity. This helps explain the preponderance of moderately eccentric bar-supporting orbits as the bar strength increases. This also suggests that bar strength may be limited by chaos! The chaos results from stochastic layers that form around primary resonances owing to separatrix splitting. Most investigations of orbit regularity are performed using numerical computation of Lyapunov exponents or related indices. We show that Lyapunov exponents poorly diagnose the degree of stochasticity in this problem; the island structure in the stochastic sheaths allow orbit to change morphology while presenting anomalously small Lyapunov exponent values (i.e. weak chaos). For example, a weakly chaotic orbit may appear to change its morphology spontaneously, while appearing regular except during the change itself. The numerical KAM approach sensitively detects these dynamics and provides a model Hamiltonian for further investigation. It may underpredict the number of broken tori for strong perturbations.

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

In this paper we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an empirical recipe for the evolving satellite mass and the corresponding mass loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direct orbital integration, or by using the action-angle formalism. The model naturally produces streaky features within the stream. These are formed due to the radial oscillation of the progenitor and the bursts of stars emitted near pericenter, rather than clumping at particular oscillation phases as sometimes suggested. When detectable, these streaky features are a reliable diagnostic for the stream's direction of motion and encode other information on the progenitor and its orbit. We show several tests of the recipe in alternate potentials, including a case with a chaotic progenitor orbit which displays a marked effect on the width of the stream. Although the specific ejection recipe may need adjusting when elements such as the orbit or satellite density profile are changed significantly, our examples suggest that model tidal streams can be quickly and accurately generated by models of this general type for use in Bayesian sampling.

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

We describe a hybrid Direct Simulation Monte Carlo (DSMC) code for simultaneously solving the collisional Boltzmann equation for gas and the collisionless Boltzmann equation for stars and dark matter for problems important to galaxy evolution. This project is motivated by the need to understand the controlling dynamics at interfaces between gases of widely differing densities and temperature, i.e. multiphase media. While more expensive than hydrodynamics, the kinetic approach does not suffer from discontinuities and it applies when the continuum limit does not, such as in the collapse of galaxy clusters and at the interface between coronal halo gas and a thin neutral gas layer. Finally, the momentum flux is carried, self-consistently, by particles and this approach explicitly resolves and thereby captures shocks. The DSMC method splits the solution into two pieces: 1) the evolution of the phase-space flow without collisions; and 2) the evolution governed the collision term alone without phase-space flow. This splitting approach makes DSMC an ideal match to existing particle-based n-body codes. If the mean free path becomes very small compared to any scale of interest, the method abandons simulated particle collisions and simply adopts the relaxed solution in each interaction cell consistent with the overall energy and momentum fluxes. This is functionally equivalent to solving the Navier-Stokes equations on a mesh. Our implementation is tested using the Sod shock tube problem and the non-linear development of an Kelvin-Helmholtz unstable shear layer.

We use the Direct Simulation Monte Carlo (DSMC) method combined with an n-body code to study the dynamics of the interaction between a gas-rich spiral galaxy and intracluster or intragroup medium, often known as the ram pressure scenario. The advantage of this gas kinetic approach over traditional hydrodynamics is explicit treatment of the interface between the hot and cold, dense and rarefied media typical of astrophysical flows and the explicit conservation of energy and momentum and the interface. This approach yields some new physical insight. Owing to the shock and backward wave that forms at the point ICM--ISM contact, ICM gas is compressed, heated and slowed. The shock morphology is Mach-disk-like. In the outer galaxy, the hot turbulent post-shock gas flows around the galaxy disk, while heating and ablating the initially cool disk gas. The outer gas and angular momentum are lost to the flow. In the inner galaxy, the hot gas pressurizes the neutral ISM gas causing a strong two-phase instability. As a result, the momentum of the wind is no longer impulsively communicated to the cold gas as assumed in the Gunn-Gott (1972) formula, but oozes through the porous disk, transferring its linear momentum to the disk en masse. The escaping gas mixture has a net positive angular momentum and forms a slowly rotating sheath. The shear flow caused by the post-shock ICM flowing through the porous multiphase ISM creates a strong Kelvin-Helmholtz instability in the disk that results in Cartwheel-like ring and spoke morphology.

Recent theoretical work has implicated fast bar formation modes and subsequent evolution as the creation mechanism for the observed peanut-shaped bulges in some edge-on disk galaxies. We demonstrate an N-body simulation of a disk undergoing a contrasting slow mode of bar growth, unsubjected to a buckling instability, which nonetheless grows the 4:1 orbit family responsible for a peanut-shaped bulge. We also present a simulation with fast mode bar growth, which exhibits thickening similar to other work. A novel orbit classification method that finds dynamically distinct families is presented, allowing for a detailed analysis of angular momentum transfer channels within the disk.

This paper introduces the Bayesian Inference Engine (BIE), a general parallel, optimised software package for parameter inference and model selection. This package is motivated by the analysis needs of modern astronomical surveys and the need to organise and reuse expensive derived data. The BIE is the first platform for computational statistics designed explicitly to enable Bayesian update and model comparison for astronomical problems. Bayesian update is based on the representation of high-dimensional posterior distributions using metric-ball-tree based kernel density estimation. Among its algorithmic offerings, the BIE emphasises hybrid tempered MCMC schemes that robustly sample multimodal posterior distributions in high-dimensional parameter spaces. Moreover, the BIE is implements a full persistence or serialisation system that stores the full byte-level image of the running inference and previously characterised posterior distributions for later use. Two new algorithms to compute the marginal likelihood from the posterior distribution, developed for and implemented in the BIE, enable model comparison for complex models and data sets. Finally, the BIE was designed to be a collaborative platform for applying Bayesian methodology to astronomy. It includes an extensible object-oriented and easily extended framework that implements every aspect of the Bayesian inference. By providing a variety of statistical algorithms for all phases of the inference problem, a scientist may explore a variety of approaches with a single model and data implementation. Additional technical details and download details are available from http://www.astro.umass.edu/bie. The BIE is distributed under the GNU GPL.

M31 has a giant stream of stars extending far to the south and a great deal of other tidal debris in its halo, much of which is thought to be directly associated with the southern stream. We model this structure by means of Bayesian sampling of parameter space, where each sample uses an N-body simulation of a satellite disrupting in M31's potential. We combine constraints on stellar surface densities from the Isaac Newton Telescope survey of M31 with kinematic data and photometric distances. This combination of data tightly constrains the model, indicating a stellar mass at last pericentric passage of log(M_s / Msun) = 9.5+-0.1, comparable to the LMC. Any existing remnant of the satellite is expected to lie in the NE Shelf region beside M31's disk, at velocities more negative than M31's disk in this region. This rules out the prominent satellites M32 or NGC 205 as the progenitor, but an overdensity recently discovered in M31's NE disk sits at the edge of the progenitor locations found in the model. M31's virial mass is constrained in this model to be log(M200) = 12.3+-0.1, alleviating the previous tension between observational virial mass estimates and expectations from the general galactic population and the timing argument. The techniques used in this paper, which should be more generally applicable, are a powerful method of extracting physical inferences from observational data on tidal debris structures.

Weinberg (2012) described a constructive algorithm for computing the marginal likelihood, Z, from a Markov chain simulation of the posterior distribution. Its key point is: the choice of an integration subdomain that eliminates subvolumes with poor sampling owing to low tail-values of posterior probability. Conversely, this same idea may be used to choose the subdomain that optimizes the accuracy of Z. Here, we explore using the simulated distribution to define a small region of high posterior probability, followed by a numerical integration of the sample in the selected region using the volume tessellation algorithm described in Weinberg (2012). Even more promising is the resampling of this small region followed by a naive Monte Carlo integration. The new enhanced algorithm is computationally trivial and leads to a dramatic improvement in accuracy. For example, this application of the new algorithm to a four-component mixture with random locations in 16 dimensions yields accurate evaluation of Z with 5% errors. This enables Bayes-factor model selection for real-world problems that have been infeasible with previous methods.

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $\Lambda$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

We introduce a new galaxy image decomposition tool, GALPHAT (GALaxy PHotometric ATtributes), to provide full posterior probability distributions and reliable confidence intervals for all model parameters. GALPHAT is designed to yield a high speed and accurate likelihood computation, using grid interpolation and Fourier rotation. We benchmark this approach using an ensemble of simulated Sersic model galaxies over a wide range of observational conditions: the signal-to-noise ratio S/N, the ratio of galaxy size to the PSF and the image size, and errors in the assumed PSF; and a range of structural parameters: the half-light radius $r_e$ and the Sersic index $n$. We characterise the strength of parameter covariance in Sersic model, which increases with S/N and $n$, and the results strongly motivate the need for the full posterior probability distribution in galaxy morphology analyses and later inferences. The test results for simulated galaxies successfully demonstrate that, with a careful choice of Markov chain Monte Carlo algorithms and fast model image generation, GALPHAT is a powerful analysis tool for reliably inferring morphological parameters from a large ensemble of galaxies over a wide range of different observational conditions. (abridged)

Computation of the marginal likelihood from a simulated posterior distribution is central to Bayesian model selection but is computationally difficult. I argue that the marginal likelihood can be reliably computed from a posterior sample by careful attention to the numerics of the probability integral. Posing the expression for the marginal likelihood as a Lebesgue integral, we may convert the harmonic mean approximation from a sample statistic to a quadrature rule. As a quadrature, the harmonic mean approximation suffers from enormous truncation error as consequence . In addition, I demonstrate that the integral expression for the harmonic-mean approximation converges slowly at best for high-dimensional problems with uninformative prior distributions. These observations lead to two computationally-modest families of quadrature algorithms that use the full generality sample posterior but without the instability. The first algorithm automatically eliminates the part of the sample that contributes large truncation error. The second algorithm uses the posterior sample to assign probability to a partition of the sample space and performs the marginal likelihood integral directly. This eliminates convergence issues. The first algorithm is analogous to standard quadrature but can only be applied for convergent problems. The second is a hybrid of cubature: it uses the posterior to discover and tessellate the subset of that sample space was explored and uses quantiles to compute a representive field value. Neither algorithm makes strong assumptions about the shape of the posterior distribution and neither is sensitive outliers. [abridged]

We investigate the physical mechanisms of tidal heating and satellite disruption in cold dark matter host haloes using N-body simulations based on cosmological initial conditions. We show the importance of resonant shocks and resonant torques with the host halo to satellite heating. A resonant shock (torque) couples the radial (tangential) motion of a satellite in its orbit to its phase space. For a satellite on a circular orbit, an ILR-like resonance dominates the heating and this heating results in continuous satellite mass loss. We estimate the requirements for simulations to achieve these dynamics using perturbation theory. Both resonant shocks and resonant torques affect satellites on eccentric orbits. We demonstrate that satellite mass loss is an outside-in process in energy space; a satellite's stars and gas are thus protected by their own halo against tidal stripping. We simulate the evolution of a halo similar to the Large Magellanic Cloud (LMC) in our Galactic dark matter halo and conclude that the LMC stars have not yet been stripped. Finally, we present a simple algorithm for estimating the evolution of satellite mass that includes both shock heating and resonant torques.

We introduce a novel image decomposition package, GALPHAT, that provides robust estimates of galaxy surface brightness profiles using Bayesian Markov Chain Monte Carlo. The GALPHAT-determined posterior distribution of parameters enables us to assign rigorous statistical confidence intervals to maximum a posteriori estimates and to test complex galaxy formation and evolution hypotheses. We describe the GALPHAT algorithm, assess its performance using test image data, and demonstrate that it has sufficient speed for production analysis of a large galaxy sample. Finally we briefly introduce our ongoing science program to study the distribution of galaxy structural properties in the local universe using GALPHAT.

Using a representative sample of 911 central galaxies (CENs) from the SDSS DR4 group catalogue, we study how the structure of the most massive members in groups and clusters depend on (1) galaxy stellar mass (Mstar), (2) dark matter halo mass of the host group (Mhalo), and (3) their halo-centric position. We establish and thoroughly test a GALFIT-based pipeline to fit 2D Sersic models to SDSS data. We find that the fitting results are most sensitive to the background sky level determination and strongly recommend using the SDSS global value. We find that uncertainties in the background translate into a strong covariance between the total magnitude, half-light size (r50), and Sersic index (n), especially for bright/massive galaxies. We find that n depends strongly on Mstar for CENs, but only weakly or not at all on Mhalo. Less (more) massive CENs tend to be disk (spheroid)-like over the full Mhalo range. Likewise, there is a clear r50-Mstar relation for CENs, with separate slopes for disks and spheroids. When comparing CENs with satellite galaxies (SATs), we find that low mass (<10e10.75 Msun/h^2) SATs have larger median n than CENs of similar Mstar. Low mass, late-type SATs have moderately smaller r50 than late-type CENs of the same Mstar. However, we find no size differences between spheroid-like CENs and SATs, and no structural differences between CENs and SATs matched in both mass and colour. The similarity of massive SATs and CENs shows that this distinction has no significant impact on the structure of spheroids. We conclude that Mstar is the most fundamental property determining the basic structure of a galaxy. The lack of a clear n-Mhalo relation rules out a distinct group mass for producing spheroids, and the responsible morphological transformation processes must occur at the centres of groups spanning a wide range of masses. (abridged)

(Abbreviated) We investigate the dynamical mechanisms responsible for producing tidal tails from dwarf satellites using N-body simulations. We identify two important dynamical co-conspirators: 1) the points where the attractive force of the host halo and satellite are balanced do not occur at equal distances from the satellite centre or at the same equipotential value for massive satellites, breaking the morphological symmetry of the leading and trailing tails; and 2) the escaped ejecta in the leading (trailing) tail continues to be decelerated (accelerated) by the satellite's gravity leading to large offsets of the ejecta orbits from the satellite orbit. The effect of the satellite's self gravity decreases only weakly with a decreasing ratio of satellite mass to host halo mass, demonstrating the importance of these effects over a wide range of subhalo masses. Not only will the morphology of the leading and trailing tails for massive satellites be different, but the observed radial velocities of the tails will be displaced from that of the satellite orbit; both the displacement and the peak radial velocity is proportional to satellite mass. If the tails are assumed to follow the progenitor satellite orbits, the tails from satellites with masses greater than 0.0001 of the host halo virial mass in a spherical halo will appear to indicate a flattened halo. Therefore, a constraint on the Milky Way halo shape using tidal streams requires mass-dependent modelling. Similarly, we compute the the distribution of tail orbits both in E_{r}-r^{-2} space and in E-L_{z} space, advocated for identifying satellite stream relics. The acceleration of ejecta by a massive satellite during escape spreads the velocity distribution and obscures the signature of a well-defined ``moving group'' in phase space.

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter halos in the $\Lambda$CDM paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though halos and disks may assemble simultaneously. We further demonstrate the validity of the simple approximation for $\Lambda$CDM halos with isotropic velocity distributions using three-dimensional N-body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low-angular momentum orbits in the vicinity of the disk and the growth time scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modeling the response of dark matter halos to the growth of a disk.

We examine the outer Galactic HI disk for deviations from the b=0 plane by constructing maps of disk surface density, mean height, and thickness. We find that the Galactic warp is well described by a vertical offset plus two Fourier modes of frequency 1 and 2, all of which grow with Galactocentric radius. A perturbation theory calculation demonstrates that the tidal influence of the Magellanic Clouds creates distortions in the dark matter halo, which combine to produce disk warp amplitudes comparable to the observations. We also construct a detailed map of the perturbed surface density of HI in the outer disk demonstrating that the Galaxy is a non-axisymmetric multi-armed spiral. Overdensities in the surface density are coincident with regions of reduced gas thickness.

We have used the HST/ACS images to identify 4700 Lyman break galaxies (LBGs) in GOODS. We present the results from a parametric analysis of the 2-D surface brightness profiles, for 1333 LBGs at z > 2.5 with rest-frame UV(1600 Angstrom) AB magnitude < -20.5. Based on the Sersic index, n, which measures the profile shape, we find that about 40% of LBGs at z=3 have light profiles close to exponential, and only 30% have the high concentrations seen for spheroids. About 30% of LBGs appear to have multiple cores or disturbed morphologies suggestive of close pairs or on-going mergers. The fraction of spheroid-like (n > 2.5) LBGs decreases by about 15% from z = 5 to 3. A comparison of LBGs with the starburst galaxies at z = 1.2, shows that the fraction of spheroid-like profiles is about 20% higher among LBGs. The ellipticity distribution for LBGs exhibits a pronounced skew towards high ellipticities (> 0.5), which cannot be explained by morphologies similar to the local disks and spheroids viewed at random orientations. The peak of the distribution evolves toward lower ellipticities, from 0.7 at z = 4 to 0.5 at z = 3. At z = 1.2 the distribution is relatively flat as seen among the present-day galaxies. The dominance of elongated morphologies among LBGs suggests that in a significant fraction of them we may be witnessing star-formation in clumps along gas-rich filaments, or the earliest gas-rich bars that encompass essentially the entire visible galaxy. Similar features are found to be ubiquitous in hydrodynamical simulations in which galaxy formation at high redshifts occurs in filamentary inflows of dynamically cold gas within the dark matter halos, and involves gas- rich mergers.

N-body simulations predict that CDM halo-assembly occurs in two phases: 1) a fast accretion phase with a rapidly deepening potential well; and 2) a slow accretion phase characterised by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM halos of the NFW form and provides physical insight into the properties of the mass accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropised by fluctuations in the gravitational potential during the fast accretion phase, we show that gravitational collapse in this phase leads to an inner profile rho(r) ~ r^{-1}. Slow accretion onto an established potential well leads to an outer profile with rho(r) ~ r^{-3}. The concentration of a halo is determined by the fraction of mass that is accreted during the fast accretion phase. Using an ensemble of realistic mass accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time, and hence the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N-body simulations. Using a simple analytic model that captures much of the important physics we show that the inner r^{-1} profile of CDM halos is a natural result of hierarchical mass assembly with a initial phase of rapid accretion.

We show that a Magellanic Cloud origin for the warp of the Milky Way can explain most quantitative features of the outer HI layer recently identified by Levine, Blitz & Heiles (2005). We construct a model similar to that of Weinberg (1998) that produces distortions in the dark matter halo, and we calculate the combined effect of these dark-halo distortions and the direct tidal forcing by the Magellanic Clouds on the disk warp in the linear regime. The interaction of the dark matter halo with the disk and resonances between the orbit of the Clouds and the disk account for the large amplitudes observed for the vertical m=0,1,2 harmonics. The observations lead to six constraints on warp forcing mechanisms and our model reasonably approximates all six. The disk is shown to be very dynamic, constantly changing its shape as the Clouds proceed along their orbit. We discuss the challenges to MOND placed by the observations.

This paper explores resonance-driven secular evolution between a bar and dark-matter halo using N-body simulations. We make direct comparisons to our analytic theory (Weinberg & Katz 2005) to demonstrate the great difficulty that an N-body simulation has representing these dynamics for realistic astronomical interactions. In a dark-matter halo, the bar's angular momentum is coupled to the central density cusp (if present) by the Inner Lindblad Resonance. Owing to this angular momentum transfer and self-consistent re-equilibration, strong realistic bars WILL modify the cusp profile, lowering the central densities within about 30% of the bar radius in a few bar orbits. Past results to the contrary (Sellwood 2006, McMillan & Dehnen 2005) may be the result of weak bars or numerical artifacts. The magnitude depends on many factors and we illustrate the sensitivity of the response to the dark-matter profile, the bar shape and mass, and the galaxy's evolutionary history. For example, if the bar length is comparable to the size of a central dark-matter core, the bar may exchange angular momentum without changing its pattern speed significantly. We emphasise that this apparently simple example of secular evolution is remarkably subtle in detail and conclude that an N-body exploration of any astronomical scenario requires a deep investigation into the underlying dynamical mechanisms for that particular problem to set the necessary requirements for the simulation parameters and method (e.g. particle number and Poisson solver). Simply put, N-body simulations do not divinely reveal truth and hence their results are not infallible. They are unlikely to provide useful insight on their own, particularly for the study of even more complex secular processes such as the production of pseudo-bulges and disk heating.

We cross correlate the well-defined and very complete spectroscopic Main Galaxy Sample (MGS) of 156,000 bright (r<17.5 mag) galaxies from the SDSS with 2MASS sources to explore the nature and completeness of the 2MASS K-band selection of nearby galaxies. 2MASS detects 90% of the MGS brighter than r=17 mag. For r<16, 93.1% of the MGS is found in the 2MASS Extended Source Catalog (XSC). These detections span the representative range of optical and near-IR galaxy properties, but with a surface brightness-dependent bias to preferentially miss the most blue and low-concentration sources, consistent with the most morphologically late-type galaxy population. An XSC completeness of 97.5% is achievable at bright magnitudes, with blue LSBs being the only major source of incompleteness, if one follows our careful matching criteria and weeds out spurious SDSS sources. We conclude that the rapid drop in XSC completeness at r>16 reflects the sharp surface-brightness limit of the extended source detection algorithm in 2MASS. As a result, the r>16 galaxies found in the XSC are over-representative in red early types and under-representative in blue latetypes. At r>16 the XSC suffers an additional selection effect from the 2-3" spatial resolution limit of 2MASS. Therefore, 2MASS continues to detect 90% of of the MGS at 16<r<17, but with a growing fraction found in the Point Source Catalog (PSC) only. Overall, one third of the MGS is detected in the 2MASS PSC but not the XSC. A combined K<13.57 and r<16 selection provides the most representative inventory of galaxies in the local cosmos with near-IR and optical measurements, and 90.8% completeness. Using data from SDSS-DR2, this sample contains 19,156 galaxies with a median redshift of 0.052. (abridged)

We have used infrared and radio observations to search for a dwarf galaxy associated with the high-velocity cloud (HVC) known as Complex H. Complex H is a large (> 400 deg^2) and probably nearby (d = 27 kpc) HVC whose location in the Galactic plane has hampered previous investigations of its stellar content. The HI mass of the cloud is 2.0 x 10^7 (d/27 kpc)^2 Msun, making Complex H one of the most massive HVCs if its distance is > ~20 kpc. Virtually all similar HI clouds in other galaxy groups are associated with LSB dwarf galaxies. We selected mid-IR MSX sources in the direction of Complex H that appeared likely to be star-forming regions and observed them at the wavelength of the CO 1-0 transition. 59 of the 60 observed sources show emission at Milky Way velocities, and we detected no emission at velocities consistent with that of Complex H. We use these observations to set an upper limit on the ongoing star formation rate in the HVC of < 5 x 10^-4 Msun yr^-1. We also searched the 2MASS database for evidence of any dwarf-galaxy-like stellar population in the direction of the HVC and found no trace of a distant red giant population, with an upper limit on the stellar mass of ~10^6 Msun. Given the lack of evidence for either current star formation or an evolved population, we conclude that Complex H cannot be a dwarf galaxy with properties similar to those of known dwarfs. Complex H is therefore one of the most massive known HI clouds without stars. If Complex H is self-gravitating, then it is one of the few known dark galaxy candidates. These findings may offer observational support for the idea that the Cold Dark Matter substructure problem is related to the difficulty of forming stars in low-mass dark matter halos; alternatively, Complex H could be an example of a cold accretion flow onto the Milky Way. [slightly abridged]

Only through resonances can non-axisymmetric features such as spiral arms and bars exert torques over large scales and change the overall structure of a near-equilibrium galaxy. We describe the resonant interaction mechanism in detail and derive explicit criteria for the particle number required to simulate these dynamical processes accurately using N-body simulations and illustrate them with numerical experiments. To do this, we perform direct numerical solution of perturbation theory and make detailed comparisons with N-body simulations. The criteria include: sufficient particle coverage in phase space near the resonance and enough particles to minimize gravitational potential fluctuations that will change the dynamics of the resonant encounter. Some of our more surprising findings are as follows. First, the Inner-Lindblad-like resonance (ILR), responsible for coupling the bar to the central halo cusp, requires almost 10^9 equal mass particles within the virial radius for a Milky-Way-like bar in an NFW profile. Second, orbits that linger near the resonance receive more angular momentum than orbits that move through the resonance quickly. Small-scale fluctuations present in state-of-the-art particle-particle simulations can knock orbits out of resonance, preventing them from lingering and, thereby, decrease the torque. The required particle numbers are sufficiently high for scenarios of interest that apparent convergence in particle number is misleading: the convergence is in the noise-dominated regime. State-of-the-art simulations are not adequate to follow all aspects of secular evolution driven by the bar-halo interaction. We present a procedure to test the requirements for individual N-body codes for the actual problem of interest. [abridged]

We calculate the angular correlation function of galaxies in the Two Micron All Sky Survey. We minimize the possible contamination by stars, dust, seeing and sky brightness by studying their cross correlation with galaxy density, and limiting the galaxy sample accordingly. We measure the correlation function at scales between 1-18 arcdegs using a half million galaxies. We find a best fit power law to the correlation function has a slope of 0.76 and an amplitude of 0.11. However, there are statistically significant oscillations around this power law. The largest oscillation occurs at about 0.8 degrees, corresponding to 2.8 h^{-1} Mpc at the median redshift of our survey, as expected in halo occupation distribution descriptions of galaxy clustering. We invert the angular correlation function using Singular Value Decomposition to measure the three-dimensional power spectrum and find that it too is in good agreement with previous measurements. A dip seen in the power spectrum at small wavenumber k is statistically consistent with CDM-type power spectra. A fit of CDM-type power spectra to k < 0.2 h Mpc^{-1} give constraints of \Gamma_{eff}=0.116 and \sigma_8=0.96. This suggest a K_s-band linear bias of 1.1+/-0.2. This \Gamma_{eff} is different from the WMAP CMB derived value. On small scales the power-law shape of our power spectrum is shallower than that derived for the SDSS. These facts together imply a biasing of these different galaxies that might be nonlinear, that might be either waveband or luminosity dependent, and that might have a nonlocal origin.

Lynden-Bell & Kalnajs (1972) presented a useful formula for computing the long-range torque between spiral arms and the disk at large. The derivation uses second-order perturbation theory and assumes that the perturbation slowly grows over a very long time: the time-asymptotic limit. This formula has been widely used to predict the angular momentum transport between spiral arms and stellar bars between disks and dark-matter halos. However, this paper shows that the LBK time-asymptotic limit is not appropriate because the characteristic evolution time for galaxies is too close to the relevant dynamical times. We demonstrate that transients, not present in the time-asymptotic formula, can play a major role in the evolution for realistic astronomical time scales. A generalisation for arbitrary time dependence is presented and illustrated by the bar--halo and satellite--halo interaction. The natural time dependence in bar-driven halo evolution causes quantitative differences in the overall torque and qualitative differences in the physical- and phase-space location of angular momentum transfer. The time-dependent theory predicts that four principal resonances dominate the torque at different times and accurately predicts the results of an N-body simulation. In addition, we show that the Inner Lindblad Resonance (ILR) is responsible for the peak angular momentum exchange but, due to the time dependence, the changes occur over a broad range of energies, radii and frequencies. We describe the implication of these findings for the satellite--halo interaction using a simple model and end with a discussion of possible impact on other aspects secular galaxy evolution.

We present moderate resolution (~6 km/s) spectroscopy of 284 M giant candidates selected from the Two Micron All Sky Survey photometry. Radial velocities (RVs) are presented for stars mainly in the south, with a number having positions consistent with association to the trailing tidal tail of the Sagittarius (Sgr) dwarf galaxy. The latter show a clear RV trend with orbital longitude, as expected from models of the orbit and destruction of Sgr. A minimum 8 kpc width of the trailing stream about the Sgr orbital midplane is implied by verified RV members. The coldness of this stream (dispersion ~10 km/s) provides upper limits on the combined contributions of stream heating by a lumpy Galactic halo and the intrinsic dispersion of released stars, which is a function of the Sgr core mass. The Sgr trailing arm is consistent with a Galactic halo containing one dominant, LMC-like lump, however some lumpier halos are not ruled out. An upper limit to the total M/L of the Sgr core is 21 in solar units. A second structure that roughly mimics expectations for wrapped, leading Sgr arm debris crosses the trailing arm in the Southern Hemisphere; however, this may also be an unrelated tidal feature. Among the <13 kpc M giants toward the South Galactic Pole are some with large RVs that identify them as halo stars, perhaps part of the Sgr leading arm near the Sun. The positions and RVs of Southern Hemisphere M giants are compared with those of southern globular clusters potentially stripped from the Sgr system and support for association of Pal 2 and Pal 12 with Sgr debris is found. Our discussion includes description of a masked-filtered cross-correlation methodology that achieves better than 1/20 of a resolution element RVs in moderate resolution spectra.

2MASS has provided a three-dimensional map of the >360 degree, wrapped tidal tails of the Sagittarius (Sgr) dwarf spheroidal galaxy, as traced by M giant stars. With the inclusion of radial velocity data for stars along these tails, strong constraints exist for dynamical models of the Milky Way-Sgr interaction. N-body simulations of Sgr disruption with model parameters spanning a range of initial conditions (e.g., Sgr mass and orbit, Galactic rotation curve, halo flattening) are used to find parameterizations that match almost every extant observational constraint of the Sgr system. We discuss the implications of the Sgr data and models for the orbit, mass and M/L of the Sgr bound core as well as the strength, flattening, and lumpiness of the Milky Way potential.

We present the first all-sky view of the Sagittarius (Sgr) dwarf galaxy mapped by M giant star tracers detected in the complete Two Micron All-Sky Survey (2MASS). The main body is fit with a King profile of 30 deg limiting radius, but with a break in the density profile from stars in tidal tails. We argue that much of the observed structure beyond the 224' core radius may be unbound as the satellite undergoes catastrophic disruption. A striking, >150 deg trailing tidal tail extends from the Sgr center and arcs across the South Galactic Hemisphere. A prominent leading debris arm extends from the Sgr center northward of the Galactic plane to an ~40 kpc apoGalacticon, loops towards the North Galactic Cap (NGC) and descends back towards the Galactic plane, foreshortened and covering the NGC. The Sgr tails lie along a well-defined orbital plane that shows little precession, which supports the notion of a nearly spherical Galactic potential. The Sun lies near the path of leading Sgr debris; thus, former Sgr stars may be near or in the solar neighborhood. The number of M giants in the Sgr tails is >15% that within the King limiting radius of the Sgr center. That several gigayear old M giants are so widespread along the Sgr tidal arms not only places limits on the dynamical age of these arms but poses a timing problem that bears on the recent binding energy of the Sgr core and that is naturally explained by recent and catastrophic mass loss. Sgr appears to contribute >75% of the high latitude, halo M giants; no evidence for M giant tidal debris from the Magellanic Clouds is found. Generally good correspondence is found between the M giant, all-sky map of the Sgr system and all previously published detections of potential Sgr debris with the exception of Sgr carbon stars -- which must be subluminous to resolve the discrepancy.

We use a large sample of galaxies from the Two Micron All Sky Survey (2MASS) and the Sloan Digital Sky Survey (SDSS) to calculate galaxy luminosity and stellar mass functions in the local Universe. We estimate k-corrections, evolution corrections, and stellar mass-to-light ratios (M/Ls) by fitting the galaxy colors with simple models. Our optical and near-infrared luminosity functions agree with most recent literature optical and near-infrared determinations within the uncertainties. We argue that 2MASS is biased against low surface brightness galaxies, and use SDSS plus our knowledge of stellar populations to estimate that the true K-band luminosity function has a steeper faint end slope than the direct estimate. Assuming a universally-applicable stellar initial mass function (IMF), we derive the stellar mass function of galaxies. The faint end slope slope for the stellar mass function is steeper than -1.1, reflecting the low stellar M/Ls characteristic of low-mass galaxies. We estimate an upper limit to the stellar mass density in the local Universe Omega* h = 0.0020+/-0.0006 by assuming an IMF as rich in low-mass stars as allowed by observations of galaxy dynamics. Finally, we find that the characteristic luminosity or mass of early-type galaxies is larger than for later types, and the faint end slope is steeper for later types than for earlier types. Accounting for typing uncertainties, we estimate that at least half, and perhaps as much as 3/4, of the stellar mass in the Universe is in early-type galaxies. We present also SDSS/2MASS color-M/L correlations, an updated discussion of near-infrared stellar M/L estimates, and the volume-corrected distribution of g and K-band stellar M/Ls as a function of stellar mass. [Abridged]

The evolution of a stellar bar transforms not only the galactic disk, but also the host dark matter halo. We present high resolution, fully self-consistent N-body simulations that clearly demonstrate that dark matter halo central density cusps flatten as the bar torques the halo. This effect is independent of the bar formation mode and occurs even for rather short bars. The halo and bar evolution is mediated by resonant interactions between orbits in the halo and the bar pattern speed, as predicted by linear Hamiltonian perturbation theory. The bar lengthens and slows as it loses angular momentum, a process that occurs even in rather warm disks. We demonstrate that the bar and halo response can be critically underestimated for experiments that are unable to resolve the relevant resonant dynamics; this occurs when the phase space in the resonant region is under sampled or plagued by noise.

The unprecedented sky coverage and photometric uniformity of the Two Micron All Sky Survey (2MASS) provides a rich resource for investigating the galaxies populating the local Universe. A full characterization of the large-scale clustering distribution is important for theoretical studies of structure formation. 2MASS offers an all-sky view of the local galaxy population at 2.15 micron, unbiased by young stellar light and minimally affected by dust. We use 2MASS to map the local distribution of galaxies, identifying the largest structures in the nearby universe. The inhomogeneity of these structures causes an acceleration on the Local Group of galaxies, which can be seen in the dipole of the Cosmic Microwave Background (CMB). We find that the direction of the 2MASS clustering dipole is 11 degrees from the CMB dipole, confirming that the local galaxy distribution accelerates the Local Group. From the magnitude of the dipole we find a value of the linear bias parameter b=1.37 +/- 0.3 in the K_s-band. The 2MASS clustering dipole is 19 degrees from the latest measurement of the dipole using galaxies detected by the Infrared Astronomical Satellite (IRAS) suggesting that bias may be non-linear in some wavebands.

We estimate the baryonic (stellar+cold gas) mass function of galaxies in the local Universe by assigning a complete sample of Two Micron All Sky Survey and Sloan Digital Sky Survey galaxies a gas fraction based on a statistical sample of the entire population, under the assumption of a universally-applicable stellar initial mass function. The baryonic mass function is similar to the stellar mass function at the high mass end, and has a reasonably steep faint-end slope owing to the typically high cold gas fractions and low stellar mass-to-light ratios characteristic of low-mass galaxies. The Schechter Function fit parameters are phi* = 0.0108(6) per cubic Mpc per dex of mass, M* = 5.3(3)x10^10 solar masses, and alpha = -1.21(5), with formal error estimates given in parentheses (for a Hubble constant of 100 km/s per Mpc). We show that the neutral and molecular hydrogen mass functions derived using this indirect route are in agreement with direct estimates, validating our indirect method. Integrating under the baryonic mass function and incorporating all sources of uncertainty, we find that the baryonic (stellar+cold gas) mass density implied by this estimate is Omega_cold baryon = 0.0024+0.0007-0.0014, or 8% +4% -5% of the Big Bang nucleosynthesis expectation.

The unprecedented sky coverage and photometric uniformity of the Two Micron All Sky Survey (2MASS) provides a rich resource for obtaining a detailed understanding of the galaxies populating our local (z<0.1) Universe. A full characterization of the physical structure of nearby galaxies is essential for theoretical and observational studies of galaxy evolution and structure formation. We have begun a quantified description of the internal structure and morphology of 10,000 bright (10<Ks<11) 2MASS galaxies through multi-component model fits to the 2D surface brightness profiles using GIM2D. From our initial Monte Carlo tests on 77 galaxies drawn at random from the RC3, we find that the model derived structural parameter errors due to sky uncertainies are typically less than 10%.

Using data from the Two Micron All-Sky Survey (2MASS), we identify a population of infrared carbon (IR C) stars with J-K >= 2 in the Milky Way. These stars are shown to trace the stellar bar previously identified in IR and optical surveys. The properties of C stars strongly suggest that they are of intermediate age. We conclude that the bar is likely to have formed more recently than 3 Gyr ago, and must be younger than 6 Gyr. Implications and further tests of this conclusion are briefly discussed.

I review the arguments for the importance of halo structure in driving galaxy evolution and coupling a galaxy to its environment. We begin with a general discussion of the key dynamics and examples of structure dominated by modes. We find that simulations with large numbers of particles (N > 1e6) are required to resolve the dynamics. Finally, I will describe some new results which demonstrates that a disk bar can produce cores in a cuspy CDM dark-matter profile within a gigayear. An inner Lindblad-like resonance couples the rotating bar to halo orbits at all radii through the cusp, rapidly flattening it. This resonance disappears for profiles with cores and is responsible for a qualitative difference in bar-driven halo evolution with and without a cusp. Although the bar gives up the angular momentum in its pattern to make the core, the formation epoch is rich in accretion events to recreate or trigger a classic stellar bar. The evolution of the cuspy inner halo by the first-generation bar paves the way for a long-lived subsequent bar with low torque and a stable pattern speed.

Simulations predict that the dark matter halos of galaxies should have central cusps, while those inferred from observed galaxies do not have cusps. We demonstrate, using both linear perturbation theory and n-body simulations, that a disk bar, which should be ubiquitous in forming galaxies, can produce cores in cuspy CDM dark matter profiles within five bar orbital times. Simulations of forming galaxies suggest that one of Milky Way size could have a 10 kpc primordial bar; this bar will remove the cusp out to approximately 5 kpc in approximately 1.5 gigayears, while the disk only loses approximately 8% of its original angular momentum. An inner Lindblad-like resonance couples the rotating bar to orbits at all radii through the cusp, transferring the bar pattern angular momentum to the dark matter cusp, rapidly flattening it. This resonance disappears for profiles with cores and is responsible for a qualitative difference in bar driven halo evolution with and without a cusp. This bar induced evolution will have a profound effect on the structure and evolution of almost all galaxies. Hence, both to understand galaxy formation and evolution and to make predictions from theory it is necessary to resolve these dynamical processes. Unfortunately, correctly resolving these important dynamical processes in ab initio calculations of galaxy formation is a daunting task, requiring at least 4,000,000 halo particles using our SCF code, and probably requiring many times more particles when using noisier tree, direct summation, or grid based techniques, the usual methods employed in such calculations.

These lectures will cover methods for studying the evolution of galaxies since their formation. Because the properties of a galaxy depend on its history, an understanding of galaxy evolution requires that we understand the dynamical interplay between all components. The first part will emphasize n-body simulation methods which minimize sampling noise. These techniques are based on harmonic expansions and scale linearly with the number of bodies, similar to Fourier transform solutions used in cosmological simulations. Although fast, until recently they were only efficiently used for small number of geometries and background profiles. These same techniques may be used to study the modes and response of a galaxy to an arbitrary perturbation. In particular, I will describe the modal spectra of stellar systems and role of damped modes which are generic to stellar systems in interactions and appear to play a significant role in determining the common structures that we see. The general development leads indirectly to guidelines for the number of particles necessary to adequately represent the gravitational field such that the modal spectrum is resolvable. I will then apply these same excitation to understanding the importance of noise to galaxy evolution.

Much of our effort in understanding the long-term evolution and morphology of the Milky Way and other galaxies has focused on the equilibrium of its luminous disk. However, the interplay between all components, seen and unseen, is a major cause of observed features and drives evolution. I will review the key underlying dynamics, and in a number of examples, show how this leads to lopsidedness and offset nuclei, may trigger bars and cause warps. Indeed, the Milky Way like most spiral galaxies show exhibit many of these features. In addition, the mechanisms suggest that observed morphology depends on the properties of the galaxy and only weakly on any particular disturbance. Because of this convergence, understanding a galaxy's history will be subtle and require the level of detail that study of the Milky Way can provide.

We present results from the application of a global photometric calibration (GPC) procedure to calibration data from the first 2 years of The Two Micron All Sky Survey (2MASS). The GPC algorithm uses photometry of both primary standards and moderately bright `tracer' stars in 35 2MASS calibration fields. During the first two years of the Survey, each standard was observed on approximately 50 nights, with about 900 individual measurements. Based on the photometry of primary standard stars and secondary tracer stars and under the assumption that the nightly zeropoint drift is linear, GPC ties together all calibration fields and all survey nights simultaneously, producing a globally optimized solution. Calibration solutions for the Northern and Southern hemisphere observatories are found separately, and are tested for global consistency based on common fields near the celestial equator. Several results from the GPC are presented, including establishing candidate secondary standards, monitoring of near-infrared atmospheric extinction coefficients, and verification of global validity of the standards. The solution gives long-term averages of the atmospheric extinction coefficients, A_J=0.096, A_H=0.026, A_{K_s}=0.066 (North) and A_J=0.092, A_H=0.031, A_{K_s}=0.065 (South), with formal error of 0.001. The residuals show small seasonal variations, most likely due to changing atmospheric content of water vapor. Extension of the GPC to approximately 100 field stars in each of the 35 calibration fields yields a catalog of more than two thousand photometric standards ranging from 10th to 14th magnitude, with photometry that is globally consistent to $\sim 1%$.

We present a theory for describing the evolution of a galaxy caused by stochastic events such as weak mergers, transient spiral structure, orbiting blobs, etc. This noise excites large-scale patterns that drives the evolution of the galactic density profile. In dark-matter haloes, the repeated stochastic perturbations preferentially ring the lowest-order modes of the halo with only a very weak dependence on the details of their source. Shaped by these modes, the profile quickly takes on a nearly self-similar form. We show that this form has the features of the ``universal profile'' reported by Navarro, Frenk, & White independent of initial conditions in a companion paper. In this sense, this noise-driven process is a near-equilibrium form of violent relaxation.

Using the theory describing the evolution of a galaxy halo due to stochastic fluctuations developed in the companion paper, we show that a halo quickly evolves toward the same self-similar profile, independent of its initial profile and concentration. The self-similar part of profile takes the form of a double power law with inner and outer exponents taking the values near -1.5 and -3 respectively. The precise value of the inner exponent depends on the magnitude and duration of the noisy epoch and most likely on form of the inner profile to start. The outer exponent is the result of evolution dominated by the external l=1 multipole resulting from the inner halo's response to noise. Three different noise processes are studied: (1) a bombardment by blobs of mass small compared to the halo mass (`shrapnel'); (2) orbital evolution of substructure by dynamical friction (`satellites'); and (3) noise caused by the orbit of blobs in the halo (`black holes'). The power spectra in the shrapnel and satellite cases is continuous and results in the double power law form, independent of initial conditions. The power spectrum for black holes is discrete and has a different form with a much slower rate of evolution. A generic prediction of this study is that noise from transient processes will drive evolution toward the same double power law with only weak constraints on the noise source and initial conditions.

We derive structural parameters and evidence for extended tidal debris from star count and preliminary standard candle analyses of the Large Magellanic Cloud based on Two Micron All Sky Survey (2MASS) data. The full-sky coverage and low extinction in K_s presents an ideal sample for structural analysis of the LMC. The star count surface densities and deprojected inclination for both young and older populations are consistent with previous work. We use the full areal coverage and large LMC diameter to Galactrocentric distance ratio to infer the same value for the disk inclination based on perspective. A standard candle analysis based on a sample of carbon long-period variables (LPV) in a narrow color range, 1.6<J-K_s<1.7 allows us to probe the three-dimensional structure of the LMC along the line of sight. The intrinsic brightness distribution of carbon LPVs in selected fields implies that $\sigma_M\simlt 0.2^m$ for this color cut. The sample provides a {\it direct} determination of the LMC disk inclination: $42.3^\circ\pm 7.2^\circ$. Distinct features in the photometric distribution suggest several distinct populations. We interpret this as the presence of an extended stellar component of the LMC, which may be as thick as 14 kpc, and intervening tidal debris at roughly 15 kpc from the LMC.

We present a morphological analysis of the feature-rich 2MASS LMC color-magnitude diagram, identifying Galactic and LMC populations and estimating the density of LMC populations alone. We also present the projected spatial distributions of various stellar populations. Major populations are identified based on matching morphological features of the CMD with expected positions of known populations, isochrone fits, and analysis of the projected spatial distributions. The LMC populations along the first-ascent RGB and AGB are quantified. We find the RGB tip at $K_s=12.3\pm0.1$. Preliminary isochrone analysis is done for giant populations in the bar and the outer regions of the Cloud. We find no significant differences in metallicities and ages between the fields. The observed LMC giant branch is well-fit by published tracks in the CIT/CTIO system with a distance modulus of $\mu=18.5\pm0.1$, reddening $E_{B-V}=0.15-0.20$, metallicity $Z=0.004^{+0.002}_{-0.001}$ and age 3-13 Gyr. Analysis of deep 2MASS engineering data with six times the standard exposure produces similar estimates.

We study the internal response of a galaxy to an unbound encounter and present a survey of orbital parameters covering typical encounters in different galactic environments. Overall, we conclude that relatively weak encounters by low-mass interloping galaxies can cause observable distortions in the primaries. The resulting asymmetries may persist long after the interloper is evident. We focus our attention on the production of structure in dark halos and in cluster ellipticals. Any distortion produced in a dark halo can distort the embedded stellar disk, possibly leading to the formation of lopsided and warped disks. Features excited in an elliptical are directly observable and we predict that asymmetries in the morphologies of these systems can be produced by relatively small perturbers. We use these detailed results to predict the distribution of the $A$ parameter defined by Abraham et al. and the shift between the center of mass of the primary system and the position of the peak of density for a range of environments. We find that high-density, low-velocity dispersion environments are more likely to host galaxies with significant asymmetries. Our distribution for the $A$ parameter is in good agreement with the range spanned by the observed values for local galaxy clusters and for distant galaxies in the MDS and in the HDF.Finally, we propose a generalized asymmetry parameter $A(r)$ which provides detailed information on the radial structure of the asymmetry produced by the mechanism explored in our work. (abridged)

A combination of analytic models and n-body simulations implies that the structural evolution of the Large Magellanic Cloud (LMC) is dominated by its dynamical interaction with the Milky Way. Although expected at some level, the scope of the involvement has significant observational consequences. First, LMC disk orbits are torqued out of the disk plane, thickening the disk and populating a spheroid. The torque results from direct forcing by the Milky Way tide and, indirectly, from the drag between the LMC disk and its halo resulting from the induced precession of the LMC disk. The latter is a newly reported mechanism that can affect all satellite interations. However, the overall torque can not isotropize the stellar orbits and their kinematics remains disk-like. Such a kinematic signature is observed for nearly all LMC populations. The extended disk distribution is predicted to increase the microlensing toward the LMC. Second, the disk's binding energy slowly decreases during this process, puffing up and priming the outer regions for subsequent tidal stripping. Because the tidally stripped debris will be spatially extended, the distribution of stripped stars is much more extended than the HI Magellanic Stream. This is consistent with upper limits to stellar densities in the gas stream and suggests a different strategy for detecting the stripped stars. And, finally, the mass loss over several LMC orbits is predicted by n-body simulation and the debris extends to tens of kiloparsecs from the tidal boundary. Although the overall space density of the stripped stars is low, possible existence of such intervening populations have been recently reported and may be detectable using 2MASS.

The satellite-disk interaction provides limits on halo properties in two ways: (1) physical arguments motivate the excitation of observable Galactic disk structure in the presence of a massive halo, although precise limits on halo parameters are scenario-dependent; (2) conversely, the Milky Way as a whole has significant dynamical effect on LMC structure and this interaction also leads to halo limits. Together, these scenarios give strong corroboration of our current gravitational mass estimates and suggests a rapidly evolving LMC.

An expansion of a density field or particle distribution in basis functions which solve the Poisson equation both provides an easily parallelized n-body force algorithm and simplifies perturbation theories. The expansion converges quickly and provides the highest computational advantage if the lowest-order potential-density pair in the basis looks like the unperturbed galaxy or stellar system. Unfortunately, there are only a handful of such basis in the literature which limits this advantage. This paper presents an algorithm for deriving these bases to match a wide variety of galaxy models. The method is based on efficient numerical solution of the Sturm-Liouville equation and can be used for any geometry with a separable Laplacian. Two cases are described in detail. First for the spherical case, the lowest order basis function pair may be chosen to be exactly that of the underlying model. The profile may be cuspy or have a core and truncated or of infinite extent. Secondly, the method yields a three-dimensional cylindrical basis appropriate for studying galaxian disks. In this case, the vertical and radial bases are coupled; the lowest order radial part of the basis function can be chosen to match the underlying profile only in the disk plane. Practically, this basis is still a very good match to the overall disk profile and converges in a small number of terms.

Gravitational amplification of Poisson noise in stellar systems increases the dipole noise power by roughly a factor of six and the quadrupole noise by 50% for a King model profile. The dipole noise is amplified by a factor of fifteen for the core-free Hernquist model. The predictions are computed using the dressed-particle formalism of Rostoker & Rosenbluth (1960) and are demonstrated by n-body simulation. This result implies that a collisionless n-body simulation is impossible; the fluctuation noise which causes relaxation is an intrinic part of self gravity. Applied to dark matter halos of disk galaxies, particle numbers of at least 10^6 will be necessary to suppress this noise at a level that does not dominate or significantly affect the disk response. Conversely, halos are most likely far from phase-mixed equilibrium and the resulting noise spectrum may seed or excite observed structure such as warps, spiral arms and bars. Fluctuation noise, similar to that due to a population of 10^6 solar-mass black holes can produce observable warping and possibly excite or seed other disk structure.

This paper explores the dynamical interaction between extended halo and spheroid components of a disk galaxy and an environmental disturbance. One finds that resonant interaction between a galaxy and passing interlopers or satellite companions can carry the disturbance inward, deep inside the halo, where it can perturb the disk. Applied to the Milky Way for example, the LMC and SMC appear to be sufficient to cause the observed Galactic warp and possibly seed other asymmetries. A polar satellite orbit produces the largest warp and therefore the inferred LMC orbit is nearly optimal for maximum warp production. Both the magnitude and morphology of the response depend on the details of the disk and halo models. The results and prognosis for N-body simulations are discussed. Discreteness noise in the halo, similar to that due to a population of $10^6M_\odot$ black holes, can produce observable warping.

Previous work reported a bar signature in color-selected IRAS variable stars. Here, we estimate the source density of these variables while consistently accounting for spatial incompleteness in data using a likelihood approach. The existence of the bar is confirmed with shoulder of 4 kpc and position angle of 19 +/- 1 degrees. The ratio of non-axisymmetric to axisymmetric components gives similar estimate for the bar size 3.3 +/- 0.1 kpc and position angle 24 +/- 2 degrees. We estimate a scale length 4.00 +/- 0.55 kpc for the IRAS variable population, suggesting that these stars represent the old disk population. We use this density reconstruction to estimate the optical depth to microlensing for the large-scale bar in the Galactic disk. We find an enhancement over an equivalent axisymmetric disk by 30% but still too small to account for the MACHO result. In addition, we predict a significant asymmetry at positive and negative longitudes along lines of sight through the end of the bar with optical depths comparable to that in Baade's window. An infrared microlensing survey may be a sensitive tool for detecting or constraining structural asymmetries. More generally, this is a pilot study for Bayesian star count analyses. Bayesian approach allows the assessment of prior probabilities to the unknown parameters of the model; the resulting likelihood function is straightforwardly modified to incorporate all available data.

We study the dynamical evolution of disk and halo globular clusters in the Milky Way using a series of Fokker-Planck calculations combined with parametric statistical models. Our sample of 113 clusters with velocity data is predicted to descend from an initial population of 250 clusters, implying more than a factor of two decrease in population size due to evolution. Approximately 200 of these clusters are in a halo component and 50 in a disk component. The estimated initial halo population follows a coreless R^{-3.38} density profile in good agreement with current estimates for the distribution of halo field stars. The observed core in the present-day distribution of halo clusters results from the rapid evaporation of clusters in the inner regions of the Galaxy. The initial halo population is also predicted to have a radially biased orbit distribution in rough agreement with the observed kinematics of halo field stars. The isotropy of the present-day halo cluster distribution results from the evaporation of clusters on elongated orbits. Similarly, the initial disk component has a nearly isotropic initial distribution that becomes more tangentially biased with time. However, the inferred initial characteristics of the disk component do not match the kinematics of the rapidly rotating thin or thick disk stellar populations. These characteristics may be more indicative of the flattened halo component discussed by Zinn (1993). Detailed examination of cluster evolution confirms the importance of disk heating. Clusters on low-inclination orbits experience the strongest disk heating because of optimal matches in resonant frequencies. Disk heating on high-inclination orbits is weaker but still dominates over spheroidal heating.

Evolution and disruption of galaxies orbiting in the gravitational field of a larger cluster galaxy are driven by three coupled mechanisms: 1) the heating due to its time dependent motion in the primary; 2) mass loss due to the tidal strain field; and 3) orbital decay. Previous work demonstrated that tidal heating is effective well inside the impulse approximation limit. Not only does the overall energy increase over previous predictions, but the work is done deep inside the secondary galaxy, e.g. at or inside the half mass radius in most cases. Here, these ideas applied to cannibalization of elliptical galaxies with fundamental-plane parameters. In summary, satellites which can fall to the center of a cluster giant by dynamical friction are evaporated by internal heating by the time they reach the center. This suggests that true merger-produced multiple nuclei giants should be rare. Specifically, secondaries with mass ratios as small as 1\% on any initial orbit evaporate and those on eccentric orbits with mass ratios as small as 0.1\% evolve significantly and nearly evaporate in a galaxian age. Captured satellites with mass ratios smaller than roughly 1\% have insufficient time to decay to the center. After many accretion events, the model predicts that the merged system has a profile similar to that of the original primary with a weak increase in concentration.

We study the combined effects of relaxation, tidal heating and binary heating on globular cluster evolution, exploring the physical consequences of external effects and examining evolutionary trends in the Milky Way population. Our analysis demonstrates that heating on circular and low-eccentricity orbits can dominate cluster evolution. The results also predict rapid evolution on eccentric orbits either due to strong relaxation caused by the high densities needed for tidal limitation or due to efficient bulge shocking of low density clusters. The combination of effects leads to strong evolution of the population as a whole. For example, within the solar circle, tidally-limited 10^5 M_sun clusters lose at least 40% of their mass in 10 Gyr. At high eccentricity most of these clusters evaporate completely. Bulge shocking disrupts clusters within 40 kpc which have less than 80% of their mass within their pericentric inner Lagrange point. Our results are consistent with suggestions that the shape of the cluster luminosity function results from evaporation and disruption of low mass clusters; they further predict that the net velocity dispersion of the cluster system in the inner Galaxy has decreased with time. Preliminary constraints on formation models are also discussed. We conclude that the observed cluster system has largely been shaped by dynamical selection.

The globular cluster population in M87 has decreased measurably through dynamical evolution caused by relaxation, binary heating and time-dependent tidal perturbation. For fundamental plane ellipticals in general, cluster populations evolve more rapidly in smaller galaxies because of the higher mass density. A simple evolutionary model reproduces the observed trend in specific frequency with luminosity for an initially constant relationship. Fits of theoretically evolved populations to M87 cluster data from McLaughlin et al. (1994) show the following: 1) dynamical effects drive evolution in the initial mass and space distributions and can account for the large core in the spatial profile as well as producing radial-dependence in the mass spectrum; 2) evolution reduces S_N by 50% within 16 kpc and 35% within 50 kpc, implying that S_N was initially 26 in this region. We estimate that 15% of the `missing' clusters lie below the detection threshold with mass less than 10^5 M_sun.

I selectively review the various dynamical scenarios that have been explored to date, especially those that illustrate the conundrums. In short, although the existence of asymmetries are convincing enough, the interpretation remains ambiguous. A coherent picture for the Milky Way asymmetries is an obvious lack; each mechanism is considered independently of all others. With this motivation, I will present a nonstandard picture which has the potential to explain some of the signatures and a few of the conundrums. Moreover, it illustrates global features are dynamically connected and I believe that it is time to revisit the Milky Way in this context. Finally, I will end with a wish list for future work---both observational and theoretical---designed to help us pin down the Milky Way.

This report describes a modification of the orthogonal function Poisson solver for n-body simulations that minimizes relaxation caused by small particle number fluctuations. With the standard algorithm, the noise leading to relaxation can be reduced by making the expansion basis similar to the particle distribution and by carefully choosing the maximum order in the expansion. The proposed algorithm accomplishes both tasks simultaneously while the simulation is running. This procedure is asymptotically equivalent to expanding in an orthogonal series which is matched to the distribution to start and truncating at low order. Because the modified algorithm adapts to a time-evolving distribution, it has advantage over a fixed basis. The required changes to the standard algorithm are minor and do not affect its overall structure or scalability. Tests show that the overhead in CPU time is small in practical applications. The decrease in relaxation rate is demonstrated for both axisymmetric and non-axisymmetric systems and the robustness of the algorithm is demonstrated by following the evolution of unstable generalized polytropes. Finally, the empirically based moment analysis which leads to the uncorrelated basis is an ideal tool for investigating structure and modes in n-body simulations and an example is provided.

Previous attempts at disturbing the galactic disk by the Magellanic Clouds relied on direct tidal forcing. However, by allowing the halo to actively respond rather than remain a rigid contributor to the rotation curve, the Clouds may produce a wake in the halo which then distorts the disk. Recent work reported here suggests that the Magellanic Clouds use this mechanism to produce disk distortions sufficient to account for both the radial location, position angle and sign of the HI warp and observed anomalies in stellar kinematics towards the galactic anticenter and LSR motion.

The adiabatic criterion, widely used in astronomical dynamics, is based on the harmonic oscillator. It asserts that the change in action under a slowly varying perturbation is exponentially small. Recent mathematical results precisely define the conditions for invariance show that this model does not apply in general. In particular, a slowly varying perturbation may cause significant evolution stellar dynamical systems even if its time scale is longer than any internal orbital time scale. This additional `heating' may have serious implications for the evolution of star clusters and dwarf galaxies which are subject to long-term environmental forces. The mathematical developments leading to these results are reviewed, and the conditions for applicability to and further implications for stellar systems are discussed. Companion papers present a computational method for a general time-dependent disturbance and detailed example.

A new theory of gravitational shocking based on time-dependent perturbation theory shows that the changes in energy and angular momentum due to a slowly varying disturbance are not exponentially small for stellar dynamical systems in general. It predicts significant shock heating by slowly varying perturbations previously thought to be negligible according to the adiabatic criterion. The theory extends the scenarios traditionally computed only with the impulse approximation and is applicable to a wide class of disturbances. The approach is applied specifically to the problem of disk shocking of star clusters.

The previous two companion papers demonstrate that slowly varying perturbations do not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with $R_g\lta8\kpc$: a concentration $c=1.5$ cluster with $R_g=8\kpc$ loses up to $95\%$ of its initial mass in $15\Gyr$. Only those with concentration $c\lta1.3$ survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will increase this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, $x$, which may be considerably larger than their initial values.

A perturber may excite a coherent mode in a star cluster or galaxy. If the stellar system is stable, it is commonly assumed that such a mode will be strongly damped and therefore of little practical consequence other than redistributing momentum and energy deposited by the perturber. This paper demonstrates that this assumption is false; weakly damped modes exist and may persist long enough to have observable consequences. To do this, a method for investigating the dispersion relation for spherical stellar systems and for locating weakly damped modes in particular is developed and applied to King models of varying concentration. This leads to the following remarkable result: King models exhibit {\it very} weakly damped $m=1$ modes over a wide range of concentration ($0.67\le c\le1.5$ have been examined). The predicted damping time is tens to hundreds of crossing times. This mode causes the peak density to shift from and slowly revolve about the initial center. The existence of the mode is supported by n-body simulation. Higher order modes and possible astronomical consequences are discussed. Weakly damped modes, for example, may provide a natural explanation for observed discrepancies between density and kinematic centers in galaxies, the location of velocity cusps due to massive black holes, and $m=1$ disturbances of disks embedded in massive halos. Gravitational shocking may excite the $m=1$ mode in globular clusters, which could modify their subsequent evolution and displace the positions of exotic remnants.

There have been several recent suggestions that the Milky Way has rotating bar-like features based on HI and star count data. In this paper, I show that such features cause distinctive stellar kinematic signatures near OLR and ILR. The effects of these resonances may be observable far from the peak density of the pattern and relatively nearby the solar position. The details of the kinematic signatures depend on the evolutionary history of the `bar' and therefore velocity data, both systemic and velocity dispersion, may be used to probe the evolutionary history as well as the present state of the Galaxy. Kinematic models for a variety of simple scenarios are presented. Models with evolving pattern speeds show significantly stronger dispersion signatures than those with static pattern speeds, suggesting that useful observational constraints are possible. The models are applied to the proposed rotating spheroid and bar models; we find: 1) none of these models chosen to represent the proposed large-scale rotating spheroid are consistent with the stellar kinematics; and 2) a Galactic bar with semimajor axis of $3\kpc$ will cause a large increase in velocity dispersion in the vicinity of OLR ($\sim5\kpc$) with little change in the net radial motion and such a signature is suggested by K-giant velocity data. Potential future observations and analyses are discussed.