Engineering PapersSearch

Engineering topics

Hernquist, Lars

Publications and source records attributed to Hernquist, Lars.

36 records · Page 2

Heating of galactic disks by mergers

The paper considers dynamical aspects of mergers between disks and small satellite companions, with special attention given to: the self-consistent response of disks to satellite mergers, which generalizes an earlier investigation by Quinn & Goodman (1986); the evolution of satellite orbits and debris; and the long-term structure and kinematics of thickened disks as compared with observed properties of disk galaxies. For satellites with masses in the range 0.04-0.2 M sub disk on initially circular inclined prograde orbits merging with disks similar to the thin disk of the Galaxy, it is found that the dynamical coupling between disks and satellites leads to orbital decay. Disks respond to mergers by spreading both radially and vertically. The radial heating is driven largely by the spiral response of the dynamically cool disks to satellite perturbations.

Quinn, P. J.

Discreteness noise versus force errors in N-body simulations

A low accuracy in the force calculation per time step of a few percent for each particle pair is sufficient for collisionless N-body simulations. Higher accuracy is made meaningless by the dominant discreteness noise in the form of two-body relaxation, which can be reduced only by increasing the number of particles. Since an N-body simulation is a Monte Carlo procedure in which each particle-particle force is essentially random, i.e., carries an error of about 1000 percent, the only requirement is a systematic averaging-out of these intrinsic errors. We illustrate these assertions with two specific examples in which individual pairwise forces are deliberately allowed to carry significant errors: tree-codes on supercomputers and algorithms on special-purpose machines with low-precision hardware.

Hernquist, Lars

Gas distribution and starbursts in shell galaxies

Detailed maps of most elliptical galaxies reveal that, whereas the greatest part of their luminous mass originates from a smooth distribution with a surface brightness approximated by a de Vaucouleurs law, a small percentage of their light is contributed by low surface brightness distortions termed 'fine structures'. The sharp-edged features called 'shells' are successfully reproduced by merger and infall models involving accretion from less massive companions. In this context, dwarf spheroidal and compact disk galaxies are likely progenitors of these stellar phenomena. However, it is probable that the sources of shell-forming material also contain significant amounts of gas. This component may play an important role in constraining the formation and evolution of shell galaxies. To investigate the effects of the gaseous component, numerical simulations were performed to study the tidal disruption of dwarf galaxies containing both gas and stars by more massive primaries, and the evolution of the ensuing debris. The calculations were performed with a hybrid N-body/hydrodynamics code. Collisionless matter is evolved using a conventional N-body technique and gas is treated using smoothed particle hydrodynamics in which self-gravitating fluid elements are represented as particles evolving according to Lagrangian hydrodynamic equations. An isothermal equation of state is employed so the gas remains at a temperature 104 K. Owing to the large mass ratio between the primary and companion, the primary is modeled as a rigid potential and the self-gravity of both galaxies is neglected.

Weil, Melinda L.

Spokes in ring galaxies

We examine the response of self-gravitating primary galaxies consisting of dark matter halos and disks containing both stars and gas to collisions with less massive companions. The primaries were constructed using a technique which makes it possible to realize multi-component systems that are stable and virtually in precise equilibrium. A total of 65,536 particles were employed to represent the primary and 4096 to represent the companion. Half the particles in the primary comprise its halo and the other half its disk. Gas makes up 10 percent of the disk mass and is represented by 8192 of the disk particles. A system of units is used where the gravitational constant, total disk mass, and disk exponential scale length are unity. The primary motivation of the present study is to determine whether effects associated with dissipation and self-gravity can account for the unusual morphology of the Cartwheel galaxy.

Hernquist, Lars

Formation of dwarf galaxies in tidal tails

The results are reported of numerical simulations of encounters between disk galaxies, each modeled with a central bulge, an exponential disk, and a spheroidal dark-matter halo. It is found that dwarf systems form in material drawn out during the encounter; these objects can capture large amounts of moderately enriched gas but retain little dark matter from their parents' haloes. They should therefore have lower mass-to-light ratios than galaxies formed directly by the collapse of primordial matter.

Barnes, Joshua E.

A novel mechanism for creating double pulsars

Simulations of encounters between pairs of hard binaries, each containing a neutron star and a main-sequence star, reveal a new formation mechanism for double pulsars in dense cores of globular clusters. In many cases, the two normal stars are disrupted to form a common envelope around the pair of neutron stars, both of which will be spun up to become millisecond pulsars. We predict that a new class of pulsars, double millisecond pulsars, will be discovered in the cores of dense globular clusters. The genesis proceeds through a short-lived double-core common envelope phase, with the envelope ejected in a fast wind. It is possible that the progenitor may also undergo a double X-ray binary phase. Any circular, short-period double pulsar found in the galaxy would necessarily come from disrupted disk clusters, unlike Hulse-Taylor class pulsars or low-mass X-ray binaries which may be ejected from clusters or formed in the galaxy.

Sigurdsson, Steinn

Structure of merger remnants. I - Bulgeless progenitors

The study examines mergers of identical galaxies consisting of self-gravitating disks and halos in the context of the suggestion that such events may form elliptical galaxies. It is shown that the luminous remnants of such mergers do indeed share many common properties with observed ellipticals. Specifically, the end states of the simulations considered rotate slowly in regions of relatively high surface density, having typical values of less than about 0.2 there. Morphologically, the remnants display a variety of structures, including shells and loops comprising loosely bound material and boxy and disky isophotes. The luminous matter is well-fitted by ellipsoidal generalizations of Hernquists's (1990, 1992) model for elliptical galaxies, implying that the surface brightness profiles are essentially de Vaucouleurs-like over a large radial interval. It is proposed that mergers of pure stellar disks do not represent an attractive mechanism for the production of massive elliptical galaxies.

Hernquist, Lars

Gravity and count probabilities in an expanding universe

The time evolution of nonlinear clustering on large scales in cold dark matter, hot dark matter, and white noise models of the universe is investigated using N-body simulations performed with a tree code. Count probabilities in cubic cells are determined as functions of the cell size and the clustering state (redshift), and comparisons are made with various theoretical models. We isolate the features that appear to be the result of gravitational instability, those that depend on the initial conditions, and those that are likely a consequence of numerical limitations. More specifically, we study the development of skewness, kurtosis, and the fifth moment in relation to variance, the dependence of the void probability on time as well as on sparseness of sampling, and the overall shape of the count probability distribution. Implications of our results for theoretical and observational studies are discussed.

Bouchet, Francois R.

Bar-spheroid interaction in galaxies

N-body simulation and linear analysis is employed to investigate the secular evolution of barred galaxies, with emphasis on the interaction between bars and spheroidal components of galaxies. This interaction is argued to drive secular transfer of angular momentum from bars to spheroids, primarily through resonant coupling. A moderately strong bar, having mass within corotation about 0.3 times the enclosed spheroid mass, is predicted to shed all its angular momentum typically in less than about 10 exp 9 yr. Even shorter depletion time scales are found for relatively more massive bars. It is suggested either that spheroids around barred galaxies are structured so as to inhibit strong coupling with bars, or that bars can form by unknown processes long after disks are established. The present models reinforce the notion that bars can drive secular evolution in galaxies.

Hernquist, Lars

Galaxies and gas in a cold dark matter universe

We use a combined gravity/hydrodynamics code to simulate the formation of structure in a random 22 Mpc cube of a cold dark matter universe. Adiabatic compression and shocks heat much of the gas to temperatures of 10 exp 6 - 10 exp 7 K, but a fraction of the gas cools radiatively to about 10 exp 4 K and condenses into discrete, highly overdense lumps. We identify these lumps with galaxies. The high-mass end of their baryonic mass function fits the form of the observed galaxy luminosity function. They retain independent identities after their dark halos merge, so gravitational clustering produces groups of galaxies embedded in relatively smooth envelopes of hot gas and dark matter. The galaxy correlation function is approximately an r exp -2.1 power law from separations of 35 kpc to 7 Mpc. Galaxy fluctuations are biased relative to dark matter fluctuations by a factor b about 1.5. We find no significant 'velocity bias' between galaxies and dark matter particles. However, virial analysis of the simulation's richest group leads to an estimated Omega of about 0.3, even though the simulation adopts Omega = 1.

Katz, Neal

Formation of shells in major mergers

Numerical simulations are used to study the fate of disks in mergers between equal-mass galaxies. Contrary to popular belief, mergers between equal-mass galaxies can form shells, loops, and ripples. Material that was originally in the outer disk of the premerger spirals falls into the remnant late in the merger event, long after the inner region of the remnant has relaxed. Thus, its evolution is similar to that of an accreted dwarf satellite that forms shells through 'phase wrapping'. However, the mechanism described in this letter avoids a number of difficulties with the accretion model; specifically, it can explain the observed correlation between properties of the host galaxies and shell alignments and luminosities. In view of this, we argue that many shell systems may have originated through 'major' mergers of comparable-mass galaxies rather than via 'minor' mergers or accretions. In particular, the shell elliptical NGC 3923 appears to be a good candidate for shell formation by a major merger.

Hernquist, Lars

Statistical mechanics of violent relaxation

We propose a functional that is extremized through violent relaxation. It is based on the Ansatz that the wave-particle scattering during violent dynamical processes can be approximated as a sequence of discrete scattering events that occur near a particle's perigalacticon. This functional has an extremum whose structure closely resembles that of spheroidal stellar systems such as elliptical galaxies. The results described here, therefore, provide a simple framework for understanding the physical nature of violent relaxation and support the view that galaxies are structured in accord with fundamental statistical principles.

Spergel, David N.

Starbursts in the nuclei of shell galaxies

Simulations of shell-forming galaxy mergers are reported in which gas and stars are tracked separately to study the implications of observations that indicate fine structures in elliptical galaxies. The hybrid N-body/hydrodynamics code by Hernquist and Katz (1989) is employed neglecting the self-gravity of galaxies in tidal disruptions by massive primaries of small dwarf galaxies. The evolution of the galactic debris is examined showing that the gas and stars are effectively segregated with the stars forming sharp-edged features and the gas forming disks in the nucleus of the primary. Star formation is theorized to occur in the nuclei of shell galaxies due to shocks in the gas, and the added presence of a black hole can lead to an active galactic nucleus. The simulation is of interest to the interpretation of such observations as those of NGC1275 in which dormant elliptical galaxies are activated by mergers between large galaxies and gas-rich satellites.

Hernquist, Lars

A self-consistent field method for galactic dynamics

The present study describes an algorithm for evolving collisionless stellar systems in order to investigate the evolution of systems with density profiles like the R exp 1/4 law, using only a few terms in the expansions. A good fit is obtained for a truncated isothermal distribution, which renders the method appropriate for galaxies with flat rotation curves. Calculations employing N of about 10 exp 6-7 are straightforward on existing supercomputers, making possible simulations having significantly smoother fields than with direct methods such as tree-codes. Orbits are found in a given static or time-dependent gravitational field; the potential, phi(r, t) is revised from the resultant density, rho(r, t). Possible scientific uses of this technique are discussed, including tidal perturbations of dwarf galaxies, the adiabatic growth of central masses in spheroidal galaxies, instabilities in realistic galaxy models, and secular processes in galactic evolution.

Hernquist, Lars

Dynamics of interacting galaxies

The dynamics of interacting galaxies as observed in the present epoch is reviewed. Topics discussed include numerical methods for modeling galactic interactions, the signatures of interacting galaxies due to tidal forces, and events that add mass to a galaxy. The review also covers major mergers between systems of comparable mass, forms of activity triggered or induced by galactic interactions, and galaxies' return to normality and related cosmological issues. Finally, some questions that yet have to be answered are examined.

Barnes, Joshua E.

A post-core-collapse model for the nucleus of M33

A model for the post-collapse core of M33 is proposed. By analogy with globular clusters, the formation of tight binary systems powers the reexpansion of the core, and it is expected that low-mass X-ray binaries form at a rate about equal to that of all galactic globular clusters combined. About a dozen such binaries should be present, and their combined emission may explain the large and enigmatic unresolved X-ray emission from the nucleus of M33. In addition, tidal formation of cataclysmic binaries may lead to a nova rate of more than about one per century. Numerous blue stragglers have probably formed through mergers caused by stellar collisions, but their density is probably too low to explain the blue color of the nucleus. Young stars are the likely cause for the color, and their presence may complicate the simple dynamical picture presented here.

Hernquist, Lars

Origin of kinematic subsystems in elliptical galaxies

The formation of a counterrotating central gas disk in a merger of two gas-rich disk galaxies of equal mass is demonstrated. Such a structure may well account for the unusual gas kinematics found in the merger remnant NGC7252. Continued star formation in such gaseous disks may produce central components with decoupled kinematics resembling the cores of some elliptical galaxies.

Hernquist, Lars

The fueling of active galaxies

Collisions of galaxies are often invoked to explain violent phenomena in the universe. The dynamics of interacting galaxies is intrinsically three-dimensional and involves both gas and stellar dynamics. In general, a computational approach is needed to model galactic collisions. Galaxy encounters are studied using a hybrid N-body/hydrodynamics code, capable of integrating systems of stars, gas, and dark matter in a fully self-consistent manner. These experiments demonstrate that gravitational coupling between gas and stars in galactic interactions can drive most of the gas throughout a galaxy into the nucleus of a merger remnant. The high densities in these gas concentrations are likely to result in strong bursts of star formation. Hence, this process may explain the nuclear starbursts in some systems of interacting galaxies. Further collapse of these gas concentrations can trigger even more intense activity if some gas is eventually accreted by a supermassive black hole. Such an evolutionary sequence may account for some quasars and active galactic nuclei.

Hernquist, Lars