Source author record

Martin D. Weinberg

Martin D. Weinberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

A Comprehensive Perturbative Formalism for Phase-Mixing in Perturbed Disks. I. Phase spirals in an Infinite, Isothermal Slab

Galactic disks are highly responsive systems that often undergo external perturbations and subsequent collisionless equilibration, predominantly via phase-mixing. We use linear perturbation theory to study the response of infinite isothermal slab analogues of disks to perturbations with diverse spatio-temporal characteristics. Without self-gravity of the response, the dominant Fourier modes that get excited in a disk are the bending and breathing modes, which, due to vertical phase-mixing, trigger local phase-space spirals that are one- and two-armed, respectively. We demonstrate how the lateral streaming motion of slab stars causes phase spirals to damp out over time. The ratio of the perturbation timescale ($τ_{\mathrm{P}}$) to the local, vertical oscillation time ($τ_z$) ultimately decides which of the two modes is excited. Faster, more impulsive ($τ_{\mathrm{P}} < τ_z$) and slower, more adiabatic ($τ_{\mathrm{P}} > τ_z$) perturbations excite stronger breathing and bending modes, respectively, although the response to very slow perturbations is exponentially suppressed. For encounters with satellite galaxies, this translates to more distant and more perpendicular encounters triggering stronger bending modes. We compute the direct response of the Milky Way disk to several of its satellite galaxies, and find that recent encounters with all of them excite bending modes in the Solar neighborhood. The encounter with Sagittarius triggers a response that is at least $1-2$ orders of magnitude larger than that due to any other satellite, including the Large Magellanic Cloud. We briefly discuss how ignoring the presence of a dark matter halo and the self-gravity of the response might impact our conclusions.

preprint2021arXiv

EXP: N-body integration using basis function expansions

We present the N-body simulation techniques in EXP. EXP uses empirically-chosen basis functions to expand the potential field of an ensemble of particles. Unlike other basis function expansions, the derived basis functions are adapted to an input mass distribution, enabling accurate expansion of highly non-spherical objects, such as galactic discs. We measure the force accuracy in three models, one based on a spherical or aspherical halo, one based on an exponential disc, and one based on a bar-based disc model. We find that EXP is as accurate as a direct-summation or tree-based calculation, and in some ways is better, while being considerably less computationally intensive. We discuss optimising the computation of the basis function representation. We also detail numerical improvements for performing orbit integrations, including timesteps.

preprint2020arXiv

Separatrix Divergence of Stellar Streams in Galactic Potentials

Flattened axisymmetric galactic potentials are known to host minor orbit families surrounding orbits with commensurable frequencies. The behavior of orbits that belong to these orbit families is fundamentally different than that of typical orbits with non-commensurable frequencies. We investigate the evolution of stellar streams on orbits near the boundaries between orbit families (separatrices) in a flattened axisymmetric potential. We demonstrate that the separatrix divides these streams into two groups of stars that belong to two different orbit families, and that as a result, these streams diffuse more rapidly than streams that evolve elsewhere in the potential. We utilize Hamiltonian perturbation theory to estimate both the timescale of this effect and the likelihood of a stream evolving close enough to a separatrix to be affected by it. We analyze two prior reports of stream-fanning in simulations with triaxial potentials, and conclude that at least one of them is caused by separatrix divergence. These results lay the foundation for a method of mapping the orbit families of galactic potentials using the morphology of stellar streams. Comparing these predictions with the currently known distribution of streams in the Milky Way presents a new way of constraining the shape of our Galaxy's potential and distribution of dark matter.

preprint2016arXiv

Dark Matter Trapping by Stellar Bars: The Shadow Bar

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.

preprint2015arXiv

Practical application of KAM theory to galactic dynamics: I. Motivation and methodology

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.

preprint2015arXiv

Practical application of KAM theory to galactic dynamics: II. Application to weakly chaotic orbits in barred galaxies

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.

preprint2015arXiv

Star Formation and Stellar Mass Assembly in Dark Matter Halos: From Giants to Dwarfs

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.

preprint2014arXiv

An Empirical Model for the Star Formation History in Dark Matter Halos

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.

preprint2014arXiv

Bayesian inferences of galaxy formation from the K-band luminosity and HI mass functions of galaxies: constraining star formation and feedback

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.

preprint2014arXiv

Generation of mock tidal streams

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.

preprint2013arXiv

A remarkably simple and accurate method for computing the Bayes Factor from a Markov chain Monte Carlo Simulation of the Posterior Distribution in high dimension

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.

preprint2013arXiv

Computational statistics using the Bayesian Inference Engine

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.

preprint2013arXiv

Creation of Peanut-Shaped Bulges via the Slow Mode of Bar Growth

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.

preprint2013arXiv

Direct Simulation Monte Carlo for astrophysical flows: I. Motivation and methodology

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.

preprint2013arXiv

Direct Simulation Monte Carlo for astrophysical flows: II. Ram pressure dynamics

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.

preprint2013arXiv

Inferring the Andromeda Galaxy's mass from its giant southern stream with Bayesian simulation sampling

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.

preprint2012arXiv

Bayesian inference of galaxy formation from the K-band luminosity function of galaxies: tensions between theory and observation

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.

preprint2011arXiv

A Bayesian approach to the semi-analytic model of galaxy formation: methodology

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 $Λ$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.

preprint2011arXiv

On the algorithms of radiative cooling in semi-analytic models

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.

preprint2010arXiv

Computing the Bayesian Factor from a Markov chain Monte Carlo Simulation of the Posterior Distribution

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]

preprint2009arXiv

The dynamics of satellite disruption in cold dark matter haloes

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.

preprint2003arXiv

The Optical and Near-Infrared Properties of Galaxies: I. Luminosity and Stellar Mass Functions

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]

preprint2000arXiv

Structure of the Large Magellanic Cloud from 2MASS

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 $σ_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.

preprint1995arXiv

Minimum relaxation n-body simulations using orthogonal series force computation

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.