Engineering topics

Hernquist, Lars

Publications and source records attributed to Hernquist, Lars.

At least 37 records · Page 2

NASA NTRS1994

Dense stellar cores in merger remnants

We use numerical models which include star formation to analyze the mass profiles of remnants formed by mergers of disk galaxies. During a merger, dissipation in gas and ensuing star formation leave behind a dense stellar core in the remnant. Rather than joining smoothly onto a de Vaucouleurs profile, the starburst population leads to a sharp break in the surface density profile at a few percent of the effective radius. While our results are preliminary, the lack of such signatures in most elliptical galaxies suggests that mergers of gas-rich disk galaxies may not have contributed greatly to the populations of present-day ellipticals.

Mihos, J. Christopher↗

NASA NTRS1994

Kinematic misalignments in remnants of multiple mergers

Observations indicate that the rotation and minor axes in elliptical galaxies are usually well aligned. However, computer models of dissipationless mergers of pairs of disk galaxies have failed to establish such a correlation, implying that most ellipticals probably did not originate in this manner. Here we consider remnants produced from repeated merging in small groups of a half-dozen or so disk galaxies and find that the internal kinematics of such objects are more representative of actual elliptical galaxies. Although our coverage of parameter space is not complete, the remnants produced in our simulations invariably possess small kinematic misalignments. These results suggest that the majority of elliptical galaxies have had evolutionary histories more complex than those envisaged by simple merger hypotheses.

Weil, Melinda L.↗

NASA NTRS1994

Ultraluminous starbursts in major mergers

We use numerical simulation to investigate the triggering of starbursts in merging disk galaxies. The properties of the merger-driven starbursts are sensitive to the structure of the progenitor galaxies; specifically, to the amount of material in a dense central bulge. Galaxies without bulges develop bars shortly after their first close passage, driving significant gas inflow and subsequent starbursts in the centers of the galaxies. These starbursts significantly deplete the star-forming gas, so that only relatively weak starbursts arise during the final merger. By contrast, models of galaxies with central bulges show that a bulge acts to stabilize the galaxies against inflow and starbursts until the galaxies actually merge. At this time, strong dissipation leads to the formation of a massive central gas mass and an ensuing star-formation rate two orders of magnitude greater than that in our isolated disk models. These starbursts are very short in duration, typically approximately 50 Myr, suggesting that the rarity of ultraluminous infrared galaxies is a result of their being in a very short evolutionary phase, rather than special and rare formation conditions. The fact that these mergers display many of the properties of ultraluminous infrared galaxies--tidal features, double nuclei, massive compact gas concentrations, and extreme levels of starburst activity--suggets that merger-driven starbursts can explain the emission from many ultraluminous infrared galaxies without an active nucleus.

Mihos, J. Christopher↗

NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗

NASA NTRS1994

Induced population gradients in galaxy merger remnants

We use numerical models to investigate the nature of induced population gradients in merger remnants. As noted by White, we find that the population mixing in stellar-dynamical mergers is rather moderate, leading to metallicity gradients in the remnant that are closely tied to gradients in the progenitor galaxies. Furthermore, the resultant metallicity gradients are poorly fitted by power laws, showing significant nonlinearities within an effective radius. If constant metallicity bulges are added to the progenitor disks, the strong relation between final and initial metallicity gradients is diluted, but the final gradients are still not well fitted by power laws. Detailed studies of the shape of the metallicity gradients in elliptical galaxies should help determine the degree to which stellar-dynamical mergers have contributed to the population of present-day ellipticals. We also consider the effects of centrally concentrated starbursts on metallicity gradients. Using simple models for metallicity enhancement, we find that addition of a metal-rich starburst population can reproduce the magnitude and shape of metallicity gradients observed in elliptical galaxies. However, even in such cases the metallicity gradients in the merger remnant should steepen significantly beyond an effective radius. Modelling such merger-induced starbursts with population synthesis techniques, we find that if the starburst is characterized by solar metallicity, it should be detectable as a bluing of the nuclear regions for several Gyr. However, if the starburst is metal-rich, the reddening effects of higher metallicity makes the broad-band signature of the starburst much more difficult to detect.

Mihos, J. Christopher↗

NASA NTRS1994

Structure of merger remnants. 4: Isophotal shapes

This paper examines the shapes of isophotes of galaxy merger remnants. More specifically, we perform a series of numerical experiments to study galaxy mergers. The simulations explore a variety of encounter geometries, types of progenitor galaxies, and particle numbers. We 'observe' each of the remnants from 64 viewpoints to estimate how the isophotal shapes vary with the orientation of the remnant. Also, by comparing the results from the various simulations, we learn how encounter geometry and the structure of the progenitors can affect the shape of the remnants. The encounter geometry, the structure of the progenitors, and the orientation of the remnant play deciding roles in the shape of the 'observed' isophotes, so much so that these simulated merger remnants have both 'boxy' and 'disky' isophotes and ellipticities ranging from E0 to E7. Thus, the shapes seen in our merger simulation remnants span a similar range to observed isophotal shapes of real galaxies. Furthermore, to estimate the errors in this analysis, we introduce the statistical bootstrap. Bootstrapping is used rather widely in observational astronomy when dealing with small samples; however, its use is practically unheard of for tackling theoretical problems. We believe that bootstrapping can be extremely useful for dealing with the small samples found in numerical simulations; therefore, we present the rudiments and basis of the technique with emphasis on its use in N-body calculations. Also, by comparing the errors for the smaller simulations with those of the larger ones, we speculate on the number of particles required to accurately explore isophotal shapes in simulations.

Heyl, Jeremy S.↗

NASA NTRS1994

Triggering of starbursts in galaxies by minor mergers

Using numerical simulation, we explore the triggering of starburst activity in disk galaxies which accrete low-mass dwarf companions. In response to the tidal perturbation of an infalling satellite, a disk galaxy develops a strong two-armed spiral pattern, which in turn drives large quantities of disk gas into its central regions. The global star formation rate stays constant during the early stages of an accretion, before rising rapidly by an order of magnitude when the central gas density becomes very large. The associated central starburst is quite compact. Models which include a bulge component in the disk galaxy show that the presence of a bulge can suppress the radial gas flow and limit the strength of the associated starburst, depending on the overall mass profile. The fact that such relatively common 'minor' mergers may trigger strong starburst activity suggests that many disk galaxies may have experienced starbursts at some point in their lifetime. Implications for galaxy evolution and formation are discussed.

Mihos, J. Christopher↗

NASA NTRS1994

Double-core evolution. 5: Three-dimensional effects in the merger of a red giant with a dwarf companion

The evolution of the common envelope phase of a binary system consisting of a 4.67 solar mass red giant and a 0.94 solar mass dwarf is studied using smoothed particle hydrodynamics. We demonstrate that the three-dimensional effects associated with the gravitational tidal torques lead to a rapid decay of the orbit on timescales approximately less than 1 yr. The relative orbit of the two cores in the common envelope is initally eccentric and tends to circularize as the orbital separation of the two cores decreases. The angular momentum lost from the orbital motion is distributed throughout the common envelope, and the double core does not evolve to a state of co-rotation for the evolutionary time followed. The energy dissipated from the relative orbit and deposited in the common envelope results in the ejection of approximately 13% of the mass of the envelope. The mass is ejected in all directions, but there is a preference for mass ejection in the orbital plane of the binary system. For example, approximately 80% of the ejected mass lies within 30 deg of the binary orbital plane. Because gravitational forces are long range, most of the energy and angular momentum is imparted to a small fraction of the common envelope resulting in an efficiency of the mass ejection process of approximately 15%. The core of the red giant executes significant displacement with respect to the center of mass of the system and contributes nearly equally to the total energy dissipation rate during the latter phases of the evolution. The degree of departure from synchronism of the initial binary system can be an important property of the system which can affect the outcome of the common envelope phase.

Terman, James L.↗

NASA NTRS1993

Primordial black holes in globular clusters

It has recently been recognized that significant numbers of medium-mass back holes (of order 10 solar masses) should form in globular clusters during the early stages of their evolution. Here we explore the dynamical and observational consequences of the presence of such a primordial black-hole population in a globular cluster. The holes initially segregate to the cluster cores, where they form binary and multiple black-hole systems. The subsequent dynamical evolution of the black-hole population ejects most of the holes on a relatively short timescale: a typical cluster will retain between zero and four black holes in its core, and possibly a few black holes in its halo. The presence of binary, triple, and quadruple black-hole systems in cluster cores will disrupt main-sequence and giant stellar binaries; this may account for the observed anomalies in the distribution of binaries in globular clusters. Furthermore, tidal interactions between a multiple black-hole system and a red giant star can remove much of the red giant's stellar envelope, which may explain the puzzling absence of larger red giants in the cores of some very dense clusters.

Sigurdsson, Steinn↗

NASA NTRS1993

Structure of merger remnants. II - Progenitors with rotating bulges

Mergers of identical galaxies consisting of self-gravitating disks, bulges, and halos are examined in the context of the suggestion that such events may form elliptical galaxies. Unlike earlier studies, the simulations reported here include effects arising from intrinsic spin of bulges. It is found that the disks and bulges are able to redistribute their angular momentum so that the luminous remnants rotate slowly near their centers. In addition, if the bulges are sufficiently concentrated, the core radii of the remnants are significantly reduced relative to those of end-states formed in mergers between pure stellar disks to the extent that the remnants share structural properties with observed elliptical galaxies. Nevertheless, it does appear that stellar-dynamical mergers between spiral progenitors will represent a viable mechanism for the production of massive elliptical galaxies only if sufficient mass resides already in dense, spheroidal components. These results suggest that any ellipticals formed in this manner having featureless light profiles were victims of a 'disk-bulge conspiracy' analogous to the disk-halo conspiracy thought to give rise to smooth rotation curves in spiral galaxies. Possible observational signatures of mergers are discussed, along with implications of the findings for our understanding of galaxy formation and evolution.

Hernquist, Lars↗

NASA NTRS1993

N-body realizations of compound galaxies

A prescription for constructing N-body models of galaxies consisting of more than one component is described. Spatial density profiles are realized exactly, but the phase space distribution is approximated using moments of the collisionless Boltzmann equation. While this approach is not fully rigorous, empirical tests suggest that it is adequate for studies of, e.g., interacting galaxies and the forced response of galaxies to imposed perturbations such as bars.

Hernquist, Lars↗

NASA NTRS1993

Segregation of gas and stars in shell galaxies

Using a code which is capable of evolving composite systems of collisionless matter and gas, we explore mergers like those thought responsible for the shells seen around many elliptical galaxies. If a small companion containing both gas and stars is accreted by a more massive primary, the stellar and gaseous debris are rapidly segregated: while the stars are free to oscillate back and forth in the primary's potential, thereby forming shells, the oppositely directed flows near the center of the primary effectively dissipate the orbital kinetic energy of the gas. Consequently, the gas settles into compact disks or rings in the nucleus of the primary, depending on orbital parameters. We note implications of these findings for the production of polar rings and the onset of nuclear activity in galaxies.

Weil, Melinda L.↗

NASA NTRS1993

Some cautionary remarks about smoothed particle hydrodynamics

Potential difficulties with smoothed particle hydrodynamics are discussed. In particular, empirical tests are used to demonstrate that the errors resulting from the use of variable smoothing can be much larger than commonly believed. Fortunately, however, these errors, which are normally small, do not appear to promote instability on small scales, such as fragmentation in self-gravitating fluids. Still, while SPH remains a useful tool for many problems of astrophysical interest, a rigorous formulation of it, which is adaptive but still satisfies conservation properties, is clearly wanting.

Hernquist, Lars↗

NASA NTRS1993

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

NASA NTRS1993

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↗

NASA NTRS1993

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

NASA NTRS1993

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↗

NASA NTRS1992

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