Engineering topics

Hernquist, Lars

Publications and source records attributed to Hernquist, Lars.

At least 19 records

NASA NTRS1995

Double core evolution. 7: The infall of a neutron star through the envelope of its massive star companion

Binary systems with properties similar to those of high-mass X-ray binaries are evolved through the common envelope phase. Three-dimensional simulations show that the timescale of the infall phase of the neutron star depends upon the evolutionary state of its massive companion. We find that tidal torques more effectively accelerate common envelope evolution for companions in their late core helium-burning stage and that the infall phase is rapid (approximately several initial orbital periods). For less evolved companions the decay of the orbit is longer; however, once the neutron star is deeply embedded within the companion's envelope the timescale for orbital decay decreases rapidly. As the neutron star encounters the high-density region surrounding the helium core of its massive companion, the rate of energy loss from the orbit increases dramatically leading to either partial or nearly total envelope ejection. The outcome of the common envelope phase depends upon the structure of the evolved companion. In particular, it is found that the entire common envelope can be ejected by the interaction of the neutron star with a red supergiant companion in binaries with orbital periods similar to those of long-period Be X-ray binaries. For orbital periods greater than or approximately equal to 0.8-2 yr (for companions of mass 12-24 solar mass) it is likely that a binary will survive the common envelope phase. For these systems, the structure of the progenitor star is characterized by a steep density gradient above the helium core, and the common envelope phase ends with a spin up of the envelope to within 50%-60% of corotation and with a slow mass outflow. The efficiency of mass ejection is found to be approximately 30%-40%. For less evolved companions, there is insufficient energy in the orbit to unbind the common envelope and only a fraction of it is ejected. Since the timescale for orbital decay is always shorter than the mass-loss timescale from the common envelope, the two cores will likely merge to form a Thorne-Zytkow object. Implications for the origin of Cyg X-3, an X-ray source consisting of a Wolf-Rayet star and a compact companion, and for the fate of the remnant binary consisting of a helium star and a neutron star are briefly discussed.

Terman, James L.↗

NASA NTRS1995

Using X-rays to determine which compact groups are illusory

If the large-scale galaxy distribution is filamentary, as suggested by some observations and recent hydrodynamical simulations, then lengthwise views of filaments will apparently produce compact groups (CGs) that are in reality stretched out along the line of sight. This possibility has been advocated recently by Hernquist, Katz, & Weinberg (1995). Here we propose a test for this hypothesis using X-ray emission from CGs. The observable quantity Q identical with L(sub x)a(sup 3)(sup p)/(L(sup 2)(sup g))T(sup 1/2)(sub x) should be proportional to the axis ratio of the group, a/c, where a and c are the short and long axes of a prolate distribution, a(sub p) is the radius of the group projected onto the sky, L(sub x) is the bolometric X-ray luminosity, L(sub g) is the group blue luminosity, and T(sub x) is the gas temperature. We find that the distribution of Q is consistent with the notion that many CGs with unusually small values of a/c are frauds, i.e., that the values of Q are anomalously small. Absent other information, it is equally possible that CGs are very gas-poor relative to rich clusters; however, this can be tested using the Sunyaev-Zeldovich effect. If the groups have a close to normal ratio of gas to total mass, but are simply stretched out along the line of sight, a Sunyaev-Zeldovich signal should be detectable.

Ostriker, Jeremiah P.↗

NASA NTRS1995

Physically detached 'compact groups'

A small fraction of galaxies appear to reside in dense compact groups, whose inferred crossing times are much shorter than a Hubble time. These short crossing times have led to considerable disagreement among researchers attempting to deduce the dynamical state of these systems. In this paper, we suggest that many of the observed groups are not physically bound but are chance projections of galaxies well separated along the line of sight. Unlike earlier similar proposals, ours does not require that the galaxies in the compact group be members of a more diffuse, but physically bound entity. The probability of physically separated galaxies projecting into an apparent compact group is nonnegligible if most galaxies are distributed in thin filaments. We illustrate this general point with a specific example: a simulation of a cold dark matter universe, in which hydrodynamic effects are included to identify galaxies. The simulated galaxy distribution is filamentary and end-on views of these filaments produce apparent galaxy associations that have sizes and velocity dispersions similar to those of observed compact groups. The frequency of such projections is sufficient, in principle, to explain the observed space density of groups in the Hickson catalog. We discuss the implications of our proposal for the formation and evolution of groups and elliptical galaxies. The proposal can be tested by using redshift-independent distance estimators to measure the line-of-sight spatial extent of nearby compact groups.

Hernquist, Lars↗

NASA NTRS1995

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

NASA NTRS1994

Star-forming galaxy models: Blending star formation into TREESPH

We have incorporated star-formation algorithms into a hybrid N-body/smoothed particle hydrodynamics code (TREESPH) in order to describe the star forming properties of disk galaxies over timescales of a few billion years. The models employ a Schmidt law of index n approximately 1.5 to calculate star-formation rates, and explicitly include the energy and metallicity feedback into the Interstellar Medium (ISM). Modeling the newly formed stellar population is achieved through the use of hybrid SPH/young star particles which gradually convert from gaseous to collisionless particles, avoiding the computational difficulties involved in creating new particles. The models are shown to reproduce well the star-forming properties of disk galaxies, such as the morphology, rate of star formation, and evolution of the global star-formation rate and disk gas content. As an example of the technique, we model an encounter between a disk galaxy and a small companion which gives rise to a ring galaxy reminiscent of the Cartwheel (AM 0035-35). The primary galaxy in this encounter experiences two phases of star forming activity: an initial period during the expansion of the ring, and a delayed phase as shocked material in the ring falls back into the central regions.

Mihos, J. Christopher↗

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↗