Showing 1–43 of 43 results for author: Lombardi, J

Outcomes of Sub-Neptune Collisions

Authors: Tuhin Ghosh, Sourav Chatterjee, James C. Lombardi

Abstract: Observed high multiplicity planetary systems are often tightly packed. Numerical studies indicate that such systems are susceptible to dynamical instabilities. Dynamical instabilities in close-in tightly packed systems, similar to those found in abundance by Kepler, often lead to planet-planet collisions. For sub-Neptunes, the dominant type of observed exoplanets, the planetary mass is concentrate… ▽ More Observed high multiplicity planetary systems are often tightly packed. Numerical studies indicate that such systems are susceptible to dynamical instabilities. Dynamical instabilities in close-in tightly packed systems, similar to those found in abundance by Kepler, often lead to planet-planet collisions. For sub-Neptunes, the dominant type of observed exoplanets, the planetary mass is concentrated in a rocky core, but the volume is dominated by a low-density gaseous envelope. For these, using the traditional `sticky-sphere' assumption is questionable. Using both N-body integration and smoothed-particle hydrodynamics, we have simulated sub-Neptune collisions for a wide range in realistic kinematic properties such as impact parameters ($b^{\prime}$) and impact velocities ($v_{\rm{im}}$) to study the possible outcomes in detail. We find that the majority ($\sim 76\%$) of the collisions with kinematic properties similar to what is expected from dynamical instabilities in multiplanet systems may not lead to mergers of sub-Neptunes. Instead, the sub-Neptunes separate from each other, often with significant atmosphere loss. When mergers do occur, they can involve significant mass loss even from the core and can sometimes lead to complete disruption of one or both planets. Sub-Neptunes merge or disrupt if $b^{\prime}<b_{\rm{c}}$, a critical value dependent on $ v_{\rm{im}}/v_{\rm{esc}}$, where $v_{\rm{esc}}$ is the escape velocity from the surface of the hypothetical merged planet assuming sticky-sphere. For $ v_{\rm{im}}/v_{\rm{esc}}\lesssim2.5$, $b_{\rm{c}}\propto (v_{\rm{im}}/v_{\rm{esc}})^{-2}$, and collisions with $b^{\prime}<b_{\rm{c}}$ typically leads to mergers. On the other hand, for $ v_{\rm{im}}/v_{\rm{esc}}\gtrsim2.5$, $b_{\rm{c}}\propto v_{\rm{im}}/v_{\rm{esc}}$, and the collisions with $b^{\prime}<b_{\rm{c}}$ can result in complete destruction of one or both sub-Neptunes. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 15 pages, 7 figures, submitted to the Astrophysical Journal; comments welcome

Wind-Reprocessed Transients from Stellar-mass Black Hole Tidal Disruption Events

Authors: Kyle Kremer, Brenna Mockler, Anthony L. Piro, James C. Lombardi Jr

Abstract: Tidal disruptions of stars by stellar-mass black holes are expected to occur frequently in dense star clusters. Building upon previous studies that performed hydrodynamic simulations of these encounters, we explore the formation and long-term evolution of the thick, super-Eddington accretion disks formed. We build a disk model that includes fallback of material from the tidal disruption, accretion… ▽ More Tidal disruptions of stars by stellar-mass black holes are expected to occur frequently in dense star clusters. Building upon previous studies that performed hydrodynamic simulations of these encounters, we explore the formation and long-term evolution of the thick, super-Eddington accretion disks formed. We build a disk model that includes fallback of material from the tidal disruption, accretion onto the black hole, and disk mass losses through winds launched in association with the super-Eddington flow. We demonstrate that bright transients are expected when radiation from the central engine powered by accretion onto the black hole is reprocessed at large radii by the optically-thick disk wind. By combining hydrodynamic simulations of these disruption events with our disk+wind model, we compute light curves of these wind-reprocessed transients for a wide range of stellar masses and encounter penetration depths. We find typical peak bolometric luminosities of roughly $10^{41}-10^{44}\,$erg/s (depending mostly on accretion physics parameters) and temperatures of roughly $10^5-10^6\,$K, suggesting peak emission in the ultraviolet/blue bands. We predict all-sky surveys such as the Vera Rubin Observatory and ULTRASAT will detect up to thousands of these events per year in dense star clusters out to distances of several Gpc. △ Less

Submitted 21 July, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 16 Pages, 13 figures, 2 tables. Accepted for publication in MNRAS

Partial Tidal Disruption of Main-Sequence Stars by Intermediate-Mass Black Holes

Authors: Fulya Kıroğlu, James C. Lombardi Jr., Kyle Kremer, Giacomo Fragione, Shane Fogarty, Frederic A. Rasio

Abstract: We study close encounters of a $1\,M_{\odot}$ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, $M_{\rm BH}= 100-10^4\,M_{\odot}$. Possible outcomes vary from a small tidal perturbation for weak enc… ▽ More We study close encounters of a $1\,M_{\odot}$ middle-age main-sequence star (modeled using MESA) with massive black holes through hydrodynamic simulations, and explore in particular the dependence of the outcomes on the black hole mass. We consider here black holes in the intermediate-mass range, $M_{\rm BH}= 100-10^4\,M_{\odot}$. Possible outcomes vary from a small tidal perturbation for weak encounters all the way to partial or full disruption for stronger encounters. We find that stronger encounters lead to increased mass loss at the first pericenter passage, in many cases ejecting the partially disrupted star on an unbound orbit. For encounters that initially produce a bound system, with only partial stripping of the star, the fraction of mass stripped from the star increases with each subsequent pericenter passage and a stellar remnant of finite mass is ultimately ejected in all cases. The critical penetration depth that separates bound and unbound remnants has a dependence on the black hole mass when $M_{\rm BH} \lesssim 10^3\,M_{\odot}$. We also find that the number of successive close passages before ejection decreases as we go from the stellar-mass black hole to the intermediate-mass black hole regime. For instance, after an initial encounter right at the classical tidal disruption limit, a $1\,M_{\odot}$ star undergoes 16 (5) pericenter passages before ejection from a $10\,M_{\odot}$ ($100\,M_{\odot}$) black hole. Observations of periodic flares from these repeated close passages could in principle indicate signatures of a partial tidal disruption event. △ Less

Submitted 12 May, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 19 pages, 8 figures, 2 tables. Submitted to ApJ

Formation of Low-mass Black Holes and Single Millisecond Pulsars in Globular Clusters

Authors: Kyle Kremer, Claire S. Ye, Fulya Kıroğlu, James C. Lombardi Jr., Scott M. Ransom, Frederic A. Rasio

Abstract: Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in tidal disruption of the star. Using $N$-body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on host cluster properties. We find that tidal disruption events occur most frequently in… ▽ More Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in tidal disruption of the star. Using $N$-body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on host cluster properties. We find that tidal disruption events occur most frequently in core-collapsed globular clusters and that roughly $25\%$ of the disrupted stars are merger products (i.e., blue straggler stars). Using hydrodynamic simulations, we model the tidal disruptions themselves (over timescales of days) to determine the mass bound to the neutron star and the properties of the accretion disks formed. In general, we find that roughly $80-90\%$ of the initial stellar mass becomes bound to the neutron star following disruption. Additionally, we find that neutron stars receive impulsive kicks of up to about $20\,$km/s as a result of the asymmetry of unbound ejecta; these kicks place these neutron stars on elongated orbits within their host cluster, with apocenter distances well outside the cluster core. Finally, we model the evolution of the (hypercritical) accretion disks on longer timescales (days to years after disruption) to estimate the accretion rate onto the neutron stars and accompanying spin-up. As long as $\gtrsim1\%$ of the bound mass accretes onto the neutron star, millisecond spin periods can be attained. We argue the growing numbers of isolated millisecond pulsars observed in globular clusters may have formed, at least in part, through this mechanism. In the case of significant mass growth, some of these neutron stars may collapse to form low-mass ($\lesssim3\,M_{\odot}$) black holes. △ Less

Submitted 27 June, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 14 pages, 4 figures, 2 tables. Accepted for publication in ApJL

Temperature Distributions in the Protosatellite Disks of Uranus, Jupiter, Neptune and Saturn and Common Structure in Their Satellite Systems

Authors: James C. Lombardi Sr

Abstract: We observe similar structures in the orbital radii of satellite and ring systems of Uranus, Jupiter and Neptune. This stimulates an investigation that the evolution of these systems and Saturn's system is initiated by a common mechanism that involves the interaction of radiation with their subnebulae (protosatellite disks). A model is presented that is characterized by resonance created through st… ▽ More We observe similar structures in the orbital radii of satellite and ring systems of Uranus, Jupiter and Neptune. This stimulates an investigation that the evolution of these systems and Saturn's system is initiated by a common mechanism that involves the interaction of radiation with their subnebulae (protosatellite disks). A model is presented that is characterized by resonance created through stimulated radiative molecular association (SRMA) reactions. In this model thermal energy is extracted from a protoplanetary disk at specific distances from the protoplanet wherever there is a match between the local thermal energy of the disk and the energy of photons impinging on the disk. Radiation is supplied by a portion of the hydrogen molecule's spectrum for the Uranian, Jovian and Neptunian disks and a portion of the hydrogen atom's spectrum for the Saturnian disk. Findings shed light on the early evolution and structure of satellite systems including the complicated ring system of Saturn. They also link satellite orbital radii size distributions to shapes of temperature distributions (TD's) in protosatellite disks. Theoretically determined TD's in the protosatellite disks of Uranus (Mousis 2004) and Saturn (Mousis Gautier and Bockelee-Moran, 2002) are essential to the present investigation. △ Less

Submitted 5 March, 2022; originally announced March 2022.

Comments: 43 pages, 22 figures

Hydrodynamics of Collisions and Close Encounters between Stellar Black Holes and Main-sequence Stars

Authors: Kyle Kremer, James C. Lombardi Jr., Wenbin Lu, Anthony L. Piro, Frederic A. Rasio

Abstract: Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a… ▽ More Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a wide range of hydrodynamic simulations of close encounters between black holes and main-sequence stars that collectively cover the parameter space of interest, and we identify and classify the various possible outcomes. In the case of nearly head-on collisions, the star is completely disrupted with roughly half of the stellar material becoming bound to the black hole. For more distant encounters near the classical tidal-disruption radius, the star is only partially disrupted on the first pericenter passage. Depending upon the interaction details, the partially disrupted stellar remnant may be tidally captured by the black hole or become unbound (in some cases, receiving a sufficiently large impulsive kick from asymmetric mass loss to be ejected from its host cluster). In the former case, the star will undergo additional pericenter passages before ultimately being disrupted fully. Based on the properties of the material bound to the black hole at the end of our simulations (in particular, the total bound mass and angular momentum), we comment upon the expected accretion process and associated electromagnetic signatures that are likely to result. △ Less

Submitted 24 June, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 23 pages, 11 figures, 3 tables. Accepted for publication in ApJ

Simulating a stellar contact binary merger -- I. Stellar models

Authors: Roger W. M. Hatfull, Natalia Ivanova, James C. Lombardi Jr

Abstract: We study the initial conditions of a common envelope (CE) event resulting in a stellar merger. A merger's dynamics could be understood through its light curve, but no synthetic light curve has yet been created for the full evolution. Using the Smoothed Particle Hydrodynamics (SPH) code StarSmasher, we have created three-dimensional (3D) models of a 1.52 $M_\odot$ star that is a plausible donor in… ▽ More We study the initial conditions of a common envelope (CE) event resulting in a stellar merger. A merger's dynamics could be understood through its light curve, but no synthetic light curve has yet been created for the full evolution. Using the Smoothed Particle Hydrodynamics (SPH) code StarSmasher, we have created three-dimensional (3D) models of a 1.52 $M_\odot$ star that is a plausible donor in the V1309 Sco progenitor. The integrated total energy profiles of our 3D models match their initial one-dimensional (1D) models to within a 0.1 per cent difference in the top 0.1 $M_\odot$ of their envelopes. We have introduced a new method for obtaining radiative flux by linking intrinsically optically thick SPH particles to a single stellar envelope solution from a set of unique solutions. For the first time, we calculated our 3D models' effective temperatures to within a few per cent of the initial 1D models, and found a corresponding improvement in luminosity by a factor of $\gtrsim10^6$ compared to ray tracing. We let our highest resolution 3D model undergo Roche-lobe overflow with a 0.16 $M_\odot$ point-mass accretor ($P\simeq1.6$ days) and found a bolometric magnitude variability amplitude of $\sim0.3$ -- comparable to that of the V1309 Sco progenitor. Our 3D models are, in the top 0.1 $M_\odot$ of the envelope and in terms of total energy, the most accurate models so far of the V1309 Sco donor star. A dynamical simulation that uses the initial conditions we presented in this paper can be used to create the first ever synthetic CE evolution light curve. △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: 13 pages, 7 figures, 3 tables. Accepted for publication in MNRAS

Outcomes of Grazing Impacts Between Sub-Neptunes in Kepler Multis

Authors: Jason Hwang, Sourav Chatterjee, James Lombardi Jr., Jason Steffen, Frederic Rasio

Abstract: Studies of high-multiplicity, tightly-packed planetary systems suggest that dynamical instabilities are common and affect both the orbits and planet structures, where the compact orbits and typically low densities make physical collisions likely outcomes. Since the structure of many of these planets is such that the mass is dominated by a rocky core, but the volume is dominated by a tenuous gas en… ▽ More Studies of high-multiplicity, tightly-packed planetary systems suggest that dynamical instabilities are common and affect both the orbits and planet structures, where the compact orbits and typically low densities make physical collisions likely outcomes. Since the structure of many of these planets is such that the mass is dominated by a rocky core, but the volume is dominated by a tenuous gas envelope, the sticky-sphere approximation, used in dynamical integrators, may be a poor model for these collisions. We perform five sets of collision calculations, including detailed hydrodynamics, sampling mass ratios and core mass fractions typical in Kepler Multis. In our primary set of calculations, we use Kepler-36 as a nominal remnant system, as the two planets have a small dynamical separation and an extreme density ratio. We use an N-body code, Mercury 6.2, to integrate initially unstable systems and study the resultant collisions in detail. We use these collisions, focusing on grazing collisions, in combination with realistic planet models created using gas profiles from Modules for Experiments in Stellar Astrophysics and core profiles using equations of state from Seager et al. (2007), to perform hydrodynamic calculations, finding scatterings, mergers, and even a potential planet-planet binary. We dynamically integrate the remnant systems, examine the stability, and estimate the final densities, finding the remnant densities are sensitive to the core masses, and collisions result in generally more stable systems. We provide prescriptions for predicting the outcomes and modeling the changes in mass and orbits following collisions for general use in dynamical integrators. △ Less

Submitted 29 November, 2017; v1 submitted 6 July, 2017; originally announced July 2017.

Comments: 29 pages, 13 figures, Accepted for publication in ApJ

Dynamics and Collisional Evolution of Closely Packed Planetary Systems

Authors: Jason A. Hwang, Jason H. Steffen, James C. Lombardi Jr., Frederic A. Rasio

Abstract: High-multiplicity Kepler systems (referred to as Kepler multis) are often tightly packed and may be on the verge of instability. Many systems of this type could have experienced past instabilities, where the compact orbits and often low densities make physical collisions likely outcomes. We use numerical simulations to study the dynamical instabilities and planet-planet interactions in a synthetic… ▽ More High-multiplicity Kepler systems (referred to as Kepler multis) are often tightly packed and may be on the verge of instability. Many systems of this type could have experienced past instabilities, where the compact orbits and often low densities make physical collisions likely outcomes. We use numerical simulations to study the dynamical instabilities and planet-planet interactions in a synthetically generated sample of closely-packed, high-multiplicity systems. We focus specifically on systems resembling Kepler-11, a Kepler multi with six planets, and run a suite of dynamical integrations, sampling the initial orbital parameters around the nominal values reported in Lissauer et al. (2011a), finding that most of the realizations are unstable, resulting in orbit crossings and, eventually, collisions and mergers. We study in detail the dependence of stability on the orbital parameters of the planets and planet-pair characteristics to identify possible precursors to instability, compare the systems that emerge from dynamical instabilities to the observed Kepler sample (after applying observational corrections), and propose possible observable signatures of these instabilities. We examine the characteristics of each planet-planet collision, categorizing collisions by the degree of contact and collision energy, and find that grazing collisions are more common than direct impacts. Since the structure of many planets found in Kepler multis is such that the mass is dominated by a rocky core, but the volume is dominated by a low-density gaseous envelope, the sticky-sphere approximation may not be valid, and we present hydrodynamic calculations of planet-planet collisions clearly deviating from this approximation. Finally, we rerun a subset of our dynamical calculations using instead a modified prescription to handle collisions, finding, in general, higher multiplicity remnant systems. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 40 pages, 16 figures, accepted for publication in MNRAS

Micro - tidal disruption events by stellar compact objects and the production of ultra-long GRBs

Authors: Hagai B. Perets, Zhuo Li, James C. Lombardi Jr., Stephen R. Milcarek Jr

Abstract: We explore full/partial tidal disruption events (TDEs) of stars/planets by stellar compact objects (Black holes; BHs; or neutron stars; NSs), which we term micro-TDEs. Disruption of a star/planet with mass $M_{\star}$ may lead to the formation of a debris disk around the BH/NS. Efficient accretion of a fraction $(f_{acc}=0.1$ of the debris may then give rise to bright energetic long (… ▽ More We explore full/partial tidal disruption events (TDEs) of stars/planets by stellar compact objects (Black holes; BHs; or neutron stars; NSs), which we term micro-TDEs. Disruption of a star/planet with mass $M_{\star}$ may lead to the formation of a debris disk around the BH/NS. Efficient accretion of a fraction $(f_{acc}=0.1$ of the debris may then give rise to bright energetic long ($10^{3}-10^{4}\,s$), X-ray/Gamma-ray flares, with total energies of up to $(f_{acc}/0.1)\times10^{52}\,(M_{\star}/{0.6\,M_{\odot}})\,$ergs, possibly resembling ultra-long GRBs/XRFs. The energy of such flares depends on the poorly constrained accretion processes. Significantly fainter flares might be produced if most of the disk mass is blown away through strong outflows. We suggest three dynamical origins for such disruptions. In the first, a star/planet is tidally disrupted following a close random encounter with a BH/NS in a dense cluster. We estimate the BH (NS) micro-TDE rates from this scenario to be few$\times10^{-6}$ (few$\times10^{-7}$) ${\rm yr}^{-1}$ per Milky-Way galaxy. Another scenario involves the interaction of wide companions due to perturbations by stars in the field, likely producing comparable but lower rates. Finally, a third scenario involves a BH/NS which gain a natal velocity kick at birth, leading to a close encounter with a binary companion and the tidal disruption of that companion. Such events could be associated with a supernova, or even with a preceding GRB/XRF event, and would likely occur hours to days after the prompt explosion; the rates of such events could be larger than those obtained from the other scenarios, depending on the preceding complex binary stellar evolution. △ Less

Submitted 1 April, 2016; v1 submitted 24 February, 2016; originally announced February 2016.

Comments: ApJ, in press

The Influence of Neutrinos on r-Process Nucleosynthesis in the Ejecta of Black Hole-Neutron Star Mergers

Authors: Luke F. Roberts, Jonas Lippuner, Matthew D. Duez, Joshua A. Faber, Francois Foucart, James C. Lombardi Jr., Sandra Ning, Christian D. Ott, Marcelo Ponce

Abstract: During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and f… ▽ More During the merger of a black hole and a neutron star, baryonic mass can become unbound from the system. Because the ejected material is extremely neutron-rich, the r-process rapidly synthesizes heavy nuclides as the material expands and cools. In this work, we map general relativistic models of black hole-neutron star (BHNS) mergers into a Newtonian smoothed particle hydrodynamics (SPH) code and follow the evolution of the thermodynamics and morphology of the ejecta until the outflows become homologous. We investigate how the subsequent evolution depends on our mapping procedure and find that the results are robust. Using thermodynamic histories from the SPH particles, we then calculate the expected nucleosynthesis in these outflows while varying the level of neutrino irradiation coming from the postmerger accretion disk. We find that the ejected material robustly produces r-process nucleosynthesis even for unrealistically high neutrino luminosities, due to the rapid velocities of the outflow. Nonetheless, we find that neutrinos can have an impact on the detailed pat- tern of the r-process nucleosynthesis. Electron neutrinos are captured by neutrons to produce protons while neutron capture is occurring. The produced protons rapidly form low mass seed nuclei for the r-process. These low mass seeds are eventually incorporated into the first r-process peak at A~78, producing mainly Ge and Se. We consider the mechanism of this process in detail and discuss if it can impact galactic chemical evolution of the first peak r-process nuclei. △ Less

Submitted 6 April, 2017; v1 submitted 28 January, 2016; originally announced January 2016.

Comments: 14 pages, 12 figures, accepted for publication in MNRAS

Stimulated Radiative Molecular Association in the Early Solar System. II. Orbital Radii of the Planets and Other Satellites of the Sun

Authors: James C. Lombardi Sr

Abstract: In a previous investigation, the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn are shown to be directly related to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model was developed involving stimulated radiative molecular association (SRMA) reactions among photons and atoms in the protosatellite disks of the planets.… ▽ More In a previous investigation, the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn are shown to be directly related to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model was developed involving stimulated radiative molecular association (SRMA) reactions among photons and atoms in the protosatellite disks of the planets. In the present investigation, the previously developed model is applied to the planets and important satellites of the Sun. A key component of the model involves resonance associated with SRMA. Through this resonance, thermal energy is extracted from the protosun's protoplanetary disk at specific distances from the protosun wherever there is a match between the local thermal energy of the disk and the energy of photons impinging on the disk. Orbital radii of the planets and satellites are related to photon energies ($E_P$ values) in the spectrum of atomic hydrogen. An expression determined previously is used to relate $E_P$ values to temperatures in the disk. Results indicate the surface temperature of the protosun at the time when the evolution of the planets begins is higher than the surface temperature of a typical T Tauri star. The present investigation offers an explanation for the existence of the asteroid and classical Kuiper belts and predicts that a primordial belt once existed in the vicinity of Neptune. It also indicates that Uranus is formed from two protoplanets and is thus consistent with the theory that the large tilt of Uranus's axis of rotation was created by the collision of two bodies. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 16 pages, 6 figures, 2 tables

Stimulated Radiative Molecular Association in the Early Solar System: Orbital Radii of Satellites of Uranus, Jupiter, Neptune, and Saturn

Authors: James C. Lombardi Sr

Abstract: The present investigation relates the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model is developed involving stimulated radiative molecular association (SRMA) reactions among the photons and atoms in the protosatellite disks of the planets. In this model thermal ene… ▽ More The present investigation relates the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model is developed involving stimulated radiative molecular association (SRMA) reactions among the photons and atoms in the protosatellite disks of the planets. In this model thermal energy is extracted from each disk due to a resonance at radii where there is a match between the temperature in the disk and a photon energy. Matter accumulates at these radii, and satellites and rings are ultimately formed. Orbital radii of satellites of Uranus, Jupiter, and Neptune are related to photon energies ($E_{PM}$ values) in the spectrum of molecular hydrogen. Orbital radii of satellites of Saturn are related to photon energies ($E_{PA}$ values) in the spectrum of atomic hydrogen. The first hint that such relationships exist is found in the linearity of the graphs of orbital radii of uranian satellites vs. orbital radii of jovian satellites, as well as in the graphs of orbital radii of uranian satellites vs. orbital radii of neptunian satellites. An expression is determined which gives the temperature in protosatellite disks where the evolution of each satellite begins. This expression is used to find temperature distributions in the disks, which are found to be similar to distributions calculated by other investigators. △ Less

Submitted 31 May, 2015; originally announced June 2015.

Comments: 30 pages, 11 figures, submitted to Planetary and Space Science

Stability and Coalescence of Massive Twin Binaries

Authors: Jason A. Hwang, James C. Lombardi Jr., Frederic A. Rasio, Vassiliki Kalogera

Abstract: Massive stars are usually found in binaries, and binaries with periods less than 10 days may have a preference for near equal component masses. In this paper we investigate the evolution of these binaries all the way to contact and the possibility that these systems can be progenitors of double neutron star binaries. The small orbital separations of observed double neutron star binaries suggest th… ▽ More Massive stars are usually found in binaries, and binaries with periods less than 10 days may have a preference for near equal component masses. In this paper we investigate the evolution of these binaries all the way to contact and the possibility that these systems can be progenitors of double neutron star binaries. The small orbital separations of observed double neutron star binaries suggest that the progenitor systems underwent a common envelope phase at least once during their evolution. Bethe & Brown (1998) proposed that massive binary twins will undergo a common envelope evolution while both components are ascending the red giant branch or asymptotic giant branch simultaneously, also known as double-core evolution. Using models generated from the stellar evolution code Evolve Zero Age Main Sequence, we determine the range of mass ratios resulting in both components simultaneously ascending the RGB or AGB as a function of the difference in birth times, t. We find that, even for a generous t=5 Myr, the minimum mass ratio qmin=0.933 for an 8 Solar Mass primary and increases for larger primaries. We use a hydrodynamics code, StarSmasher, to study specifically the evolution of q=1 common envelope systems as a function of initial component mass, age, and orbital separation. We find the dynamical stability limit, the largest orbital separation where the binary becomes dynamically unstable, as a function of the component mass and age. Finally, we calculate the efficiency of ejecting matter during the inspiral phase to extrapolate the properties of the remnant binary from our numerical results, assuming the common envelope is completely ejected. We find that for the nominal core masses, there is a minimum orbital separation for a given component mass such that the helium cores survive common envelope evolution in a tightly bound binary and are viable progenitors for double neutron stars. △ Less

Submitted 21 May, 2015; originally announced May 2015.

Comments: 51 pages, 26 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 806, Issue 1, article id. 135, 23 pp. (2015)

Recombination energy in double white dwarf formation

Authors: Jose L. A. Nandez, Natalia Ivanova, James C. Lombardi Jr

Abstract: In this Letter we investigate the role of recombination energy during a common envelope event. We confirm that taking this energy into account helps to avoid the formation of the circumbinary envelope commonly found in previous studies. For the first time, we can model a complete common envelope event, with a clean compact double white dwarf binary system formed at the end. The resulting binary or… ▽ More In this Letter we investigate the role of recombination energy during a common envelope event. We confirm that taking this energy into account helps to avoid the formation of the circumbinary envelope commonly found in previous studies. For the first time, we can model a complete common envelope event, with a clean compact double white dwarf binary system formed at the end. The resulting binary orbit is almost perfectly circular. In addition to considering recombination energy, we also show that between 1/4 and 1/2 of the released orbital energy is taken away by the ejected material. We apply this new method to the case of the double-white dwarf system WD 1101+364, and we find that the progenitor system at the start of the common envelope event consisted of a $\sim1.5M_\odot$ red giant star in a $\sim 30$ day orbit with a white dwarf companion. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: 5 pages, 2 figures, 3 tables, accepted for publication in MNRAS Letters (acceptance date:8 March 2015)

Testing common envelopes on double white dwarf binaries

Authors: Jose L. A. Nandez, Natalia Ivanova, James C. Lombardi Jr

Abstract: The formation of a double white dwarf binary likely involves a common envelope (CE) event between a red giant and a white dwarf (WD) during the most recent episode of Roche lobe overflow mass transfer. We study the role of recombination energy with hydrodynamic simulations of such stellar interactions. We find that the recombination energy helps to expel the common envelope entirely, while if reco… ▽ More The formation of a double white dwarf binary likely involves a common envelope (CE) event between a red giant and a white dwarf (WD) during the most recent episode of Roche lobe overflow mass transfer. We study the role of recombination energy with hydrodynamic simulations of such stellar interactions. We find that the recombination energy helps to expel the common envelope entirely, while if recombination energy is not taken into account, a significant fraction of the common envelope remains bound. We apply our numerical methods to constrain the progenitor system for WD 1101+364 -- a double WD binary that has well-measured mass ratio of $q=0.87\pm0.03$ and an orbital period of 0.145 days. Our best-fit progenitor for the pre-common envelope donor is a 1.5 $M_\odot$ red giant. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: 5 pages, 3 figures, Contribution to the proceedings of "19th European White Dwarf Workshop", Montreal, Canada, August 11-15th, 2014, ASP Conference Series

An Efficient Radiative Cooling Approximation for Use in Hydrodynamic Simulations

Authors: James C. Lombardi Jr., William G. McInally, Joshua A. Faber

Abstract: To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive. Stamatellos and collaborators introduced a scheme for smoothed particle hydrodynamics (SPH) simulations based on the notion of polytropic pseudo-clouds that uses only… ▽ More To make relevant predictions about observable emission, hydrodynamical simulation codes must employ schemes that account for radiative losses, but the large dimensionality of accurate radiative transfer schemes is often prohibitive. Stamatellos and collaborators introduced a scheme for smoothed particle hydrodynamics (SPH) simulations based on the notion of polytropic pseudo-clouds that uses only local quantities to estimate cooling rates. The computational approach is extremely efficient and works well in cases close to spherical symmetry, such as in star formation problems. Unfortunately, the method, which takes the local gravitational potential as an input, can be inaccurate when applied to non-spherical configurations, limiting its usefulness when studying disks or stellar collisions, among other situations of interest. Here, we introduce the "pressure scale height method," which incorporates the fluid pressure scale height into the determination of column densities and cooling rates, and show that it produces more accurate results across a wide range of physical scenarios while retaining the computational efficiency of the original method. The tested models include spherical polytropes as well as disks with specified density and temperature profiles. We focus on applying our techniques within an SPH code, although our method can be implemented within any particle-based Lagrangian or grid-based Eulerian hydrodynamic scheme. Our new method may be applied in a broad range of situations, including within the realm of stellar interactions, collisions, and mergers. △ Less

Submitted 13 November, 2014; originally announced November 2014.

Comments: 12 pages, 7 figures, submitted to MNRAS

Implications of the delayed 2013 outburst of the ESO 243-49 HLX-1

Authors: O. Godet, J. Lombardi, F. Antonini, N. A. Webb, D. Barret, J. Vingless, M. Thomas

Abstract: After showing four outbursts spaced by $\sim 1$ year from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al. (2011), we proposed that the X-ray lightcurve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the… ▽ More After showing four outbursts spaced by $\sim 1$ year from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al. (2011), we proposed that the X-ray lightcurve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the IMBH in a highly eccentric orbit. In this scenario, the delay can be explained only if the orbital parameters can change suddenly from orbit to orbit. To investigate this, we ran Newtonian smooth particle hydrodynamical simulations starting with an incoming donor approaching an IMBH on a parabolic orbit. We survey a large parameter space by varying the star-to-BH mass ratio ($10^{-5}-10^{-3}$) and the periapsis separation $r_p$ from 2.2 to $2.7~r_t$ with $r_t$, the tidal radius. To model the donor, we choose several polytropes ($Γ= 5/2,~n=3/2$; $Γ=3/2,~n=2$; $Γ=5/3,~n=2$ & $n=3$). Once the system is formed, the orbital period decreases until reaching a minimum. Then, the period tends to increase over several periapsis passages due to tidal effects and increasing mass transfer, leading ultimately to the ejection of the donor. The development of stochastic fluctuations inside the donor could lead to sudden changes in the orbital period from orbit to orbit with the appropriate order of magnitude of what has been observed for HLX-1. Given the constraints on the BH mass ($M_{\rm BH} > 10^4~M_\odot$) and assuming that HLX-1 is currently near a minimum in period of $\sim 1$ yr, the donor has to be a white dwarf or a stripped giant core. We predict that if HLX-1 is indeed emerging from a minimum in orbital period, then the period would generally increase with each passage, although substantial stochastic fluctuations can be superposed on this trend. △ Less

Submitted 8 August, 2014; originally announced August 2014.

Comments: 22 pages, 10 figures, accepted for publication in ApJ (acceptance date: 2014-08-01)

V1309 SCO - Understanding a merger

Authors: Jose L. A. Nandez, Natalia Ivanova, James Lombardi Jr

Abstract: One of the two outcomes of a common envelope event is a merger of the two stars. To date, the best known case of a binary merger is the V1309 Sco outburst, where the orbital period was known and observed to decay up to the outburst. Using the hydrodynamical code StarSmasher, we study in detail which characteristics of the progenitor binary affect the outburst and produce the best match with observ… ▽ More One of the two outcomes of a common envelope event is a merger of the two stars. To date, the best known case of a binary merger is the V1309 Sco outburst, where the orbital period was known and observed to decay up to the outburst. Using the hydrodynamical code StarSmasher, we study in detail which characteristics of the progenitor binary affect the outburst and produce the best match with observations. We have developed a set of tools in order to quantify any common envelope event. We consider binaries consisting of a 1.52$M_\odot$ giant and a 0.16$M_\odot$ companion with $P_{\rm orb}\sim1.4$ days, varying the nature of the companion and its synchronization. We show that all considered progenitor binaries evolve towards the merger primarily due to the Darwin instability. The merger is accompanied by mass ejection that proceeds in several separate mass outbursts and takes away a few percent of the donor mass. This very small mass, nonetheless, is important as it is not only sufficient to explain the observed light-curve, but it also carries away up to 1/3 of the both initial total angular momentum and initial orbital energy. We find that all synchronized systems experience $L_2$ mass loss that operates during just a few days prior to merger and produces ring-shaped ejecta. The formed star is always a strongly heated radiative star that differentially rotates. We conclude that the case of a synchronized binary with a main-sequence companion gives the best match with observations of V1309 Sco. △ Less

Submitted 16 April, 2014; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: Published in ApJ, 24 pages, 23 figures, 2 figures for animations (online only)

Identification of the Long-Sought Common-Envelope Events

Authors: N. Ivanova, S. Justham, J. L. Avendano Nandez, J. C. Lombardi Jr

Abstract: Common-envelope events (CEEs), during which two stars temporarily orbit within a shared envelope, are believed to be vital for the formation of a wide range of close binaries. For decades, the only evidence that CEEs actually occur has been indirect, based on the existence of systems that could not be otherwise explained. Here we propose a direct observational signature of CEE arising from a physi… ▽ More Common-envelope events (CEEs), during which two stars temporarily orbit within a shared envelope, are believed to be vital for the formation of a wide range of close binaries. For decades, the only evidence that CEEs actually occur has been indirect, based on the existence of systems that could not be otherwise explained. Here we propose a direct observational signature of CEE arising from a physical model where emission from matter ejected in a CEE is controlled by a recombination front as the matter cools. The natural range of timescales and energies from this model, as well the expected colors, light-curve shapes, ejection velocities and event rate, match those of a recently-recognized class of red transient outbursts. △ Less

Submitted 24 January, 2013; originally announced January 2013.

Comments: 6 main and 22 supplemental pages, 5 total figures, one table and 2 movies. This is the authors version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science, vol 339, 2013, http://www.sciencemag.org/content/339/6118/433.abstract

Journal ref: Science, 339, 433 (2013)

On the origin of the metallicity dependence in dynamically formed extragalactic low-mass X-ray binaries

Authors: N. Ivanova, T. Fragos, D. -W. Kim, G. Fabbiano, J. L. Avendano Nandez, J. C. Lombardi, G. R. Sivakoff, R. Voss, A. Jordan

Abstract: Globular clusters (GCs) effectively produce dynamically-formed low-mass X-ray binaries (LMXBs). Observers detect ~100 times more LMXBs per stellar mass in GCs compared to stars in the fields of galaxies. It has also been observationally established that metal-rich GCs are about 3 times more likely to contain an X-ray source than their metal-poor counterparts. Recent observations have shown that th… ▽ More Globular clusters (GCs) effectively produce dynamically-formed low-mass X-ray binaries (LMXBs). Observers detect ~100 times more LMXBs per stellar mass in GCs compared to stars in the fields of galaxies. It has also been observationally established that metal-rich GCs are about 3 times more likely to contain an X-ray source than their metal-poor counterparts. Recent observations have shown that this ratio holds in extragalactic GCs for all bright X-ray sources with Lx between 2X10^{37} and 5X10^{38} erg/s. In this Letter, we propose that the observed metallicity dependence of LMXBs in extragalactic GCs can be explained by the differences in the number densities and average masses of red giants in populations of different metallicities. Red giants serve as seeds for the dynamical production of bright LMXBs via two channels - binary exchanges and physical collisions - and the increase of the number densities and masses of red giants boost LMXB production, leading to the observed difference. We also discuss a possible effect of the age difference in stellar populations of different metallicities. △ Less

Submitted 29 August, 2012; originally announced August 2012.

Comments: 5 pages, 2 figures, submitted to ApJL

Accretion disks around kicked black holes: Post-kick Dynamics

Authors: Marcelo Ponce, Joshua A. Faber, James C. Lombardi Jr

Abstract: Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion di… ▽ More Numerical calculations of merging black hole binaries indicate that asymmetric emission of gravitational radiation can kick the merged black hole at up to thousands of km/s, and a number of systems have been observed recently whose properties are consistent with an active galactic nucleus containing a supermassive black hole moving with substantial velocity with respect to its broader accretion disk. We study here the effect of an impulsive kick delivered to a black hole on the dynamical evolution of its accretion disk using a smoothed particle hydrodynamics code, focusing attention on the role played by the kick angle with respect to the orbital angular momentum vector of the pre-kicked disk. We find that for more vertical kicks, for which the angle between the kick and the normal vector to the disk $θ\lesssim 30^\circ$, a gap remains present in the inner disk, in accordance with the prediction from an analytic collisionless Keplerian disk model, while for more oblique kicks with $θ\gtrsim 45^\circ$, matter rapidly accretes toward the black hole. There is a systematic trend for higher potential luminosities for more oblique kick angles for a given black hole mass, disk mass and kick velocity, and we find large amplitude oscillations in time in the case of a kick oriented $60^\circ$ from the vertical. △ Less

Submitted 19 December, 2011; v1 submitted 8 July, 2011; originally announced July 2011.

Comments: Movies available as 'ancillary files'. Accepted for publication in ApJ

Journal ref: ApJ 745 (2012) 71

Twin Binaries: Studies of Stability, Mass Transfer, and Coalescence

Authors: James C. Lombardi Jr., William Holtzman, Katherine L. Dooley, Kyle Gearity, Vassiliki Kalogera, Frederic A. Rasio

Abstract: Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core ma… ▽ More Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. In addition to the formation of binary neutron stars, we also discuss the implications of our results for the production of planetary nebulae with double degenerate central binaries. △ Less

Submitted 14 June, 2011; v1 submitted 7 September, 2010; originally announced September 2010.

Comments: 17 pages, accepted to ApJ, final version includes discussion of planetary nebulae with central binaries and a new figure about shock heating, visualizations at http://webpub.allegheny.edu/employee/j/jalombar/movies/

Tidal breakup of binary stars at the Galactic Center. II. Hydrodynamic simulations

Authors: Fabio Antonini, James C. Lombardi, David Merritt

Abstract: In Paper I, we followed the evolution of binary stars as they orbited near the supermassive black hole (SMBH) at the Galactic center, noting the cases in which the two stars would come close enough together to collide. In this paper we replace the point-mass stars by fluid realizations, and use a smoothed-particle hydrodynamics (SPH) code to follow the close interactions. We model the binary compo… ▽ More In Paper I, we followed the evolution of binary stars as they orbited near the supermassive black hole (SMBH) at the Galactic center, noting the cases in which the two stars would come close enough together to collide. In this paper we replace the point-mass stars by fluid realizations, and use a smoothed-particle hydrodynamics (SPH) code to follow the close interactions. We model the binary components as main-sequence stars with initial masses of 1, 3 and 6 Solar masses, and with chemical composition profiles taken from stellar evolution codes. Outcomes of the close interactions include mergers, collisions that leave both stars intact, and ejection of one star at high velocity accompanied by capture of the other star into a tight orbit around the SMBH. For the first time, we follow the evolution of the collision products for many ($\gtrsim 100$) orbits around the SMBH. Stars that are initially too small to be tidally disrupted by the SMBH can be puffed up by close encounters or collisions, with the result that tidal stripping occurs in subsequent periapse passages. In these cases, mass loss occurs episodically, sometimes for hundreds of orbits before the star is completely disrupted. Repeated tidal flares, of either increasing or decreasing intensity, are a predicted consequence. In collisions involving a low-mass and a high-mass star, the merger product acquires a high core hydrogen abundance from the smaller star, effectively resetting the nuclear evolution "clock" to a younger age. Elements like Li, Be and B that can exist only in the outermost envelope of a star are severely depleted due to envelope ejection during collisions and due to tidal forces from the SMBH. In the absence of collisions, tidal spin-up of stars is only important in a narrow range of periapse distances, $r_t/2\lesssim r_per \lesssim r_t$ with $r_t$ the tidal disruption radius. △ Less

Submitted 27 February, 2011; v1 submitted 31 August, 2010; originally announced August 2010.

Comments: ApJ accepted, 22 pages, 19 figures. Version with high-resolution figures, and additional animations, available at this url: http://astrophysics.rit.edu/fantonini/tbbs2/

Formation of black-hole X-ray binaries in globular clusters

Authors: N. Ivanova, S. Chaichenets, J. Fregeau, C. O. Heinke, J. C. Lombardi Jr., T. Woods

Abstract: Inspired by the recent identification of the first candidate BH-WD X-ray binaries, where the compact accretors may be stellar-mass black hole candidates in extragalactic globular clusters, we explore how such binaries could be formed in a dynamical environment. We provide analyses of the formation rates via well known formation channels like binary exchange and physical collisions and propose th… ▽ More Inspired by the recent identification of the first candidate BH-WD X-ray binaries, where the compact accretors may be stellar-mass black hole candidates in extragalactic globular clusters, we explore how such binaries could be formed in a dynamical environment. We provide analyses of the formation rates via well known formation channels like binary exchange and physical collisions and propose that the only possibility to form BH-WD binaries is via coupling these usual formation channels with subsequent hardening and/or triple formation. Indeed, we find that the most important mechanism to make a BH-WD X-ray binary from an initially dynamically formed BH-WD binary is triple induced mass transfer via the Kozai mechanism. Even using the most optimistic estimates for the formation rates, we cannot match the observationally inferred production rates if black holes undergo significant evaporation from the cluster or form a completely detached subcluster of black holes. We estimate that at least 1% of all formed black holes, or presumably 10% of the black holes present in the core now, must be involved in interactions with the rest of the core stellar population. △ Less

Submitted 11 January, 2010; originally announced January 2010.

Comments: 10 pages, 2 figures, submitted to ApJ

On the onset of runaway stellar collisions in dense star clusters - II. Hydrodynamics of three-body interactions

Authors: Evghenii Gaburov, James Lombardi, Simon Portegies Zwart

Abstract: The onset of runaway stellar collisions in young star clusters is more likely to initiate with an encounter between a binary and a third star than between two single stars. Using the initial conditions of such three-star encounters from direct $N$-body simulations, we model the resulting interaction by means of Smoothed Particle Hydrodynamics (SPH). We find that, in the majority of the cases con… ▽ More The onset of runaway stellar collisions in young star clusters is more likely to initiate with an encounter between a binary and a third star than between two single stars. Using the initial conditions of such three-star encounters from direct $N$-body simulations, we model the resulting interaction by means of Smoothed Particle Hydrodynamics (SPH). We find that, in the majority of the cases considered, all three stars merge together, and in such three star mergers, the hydrodynamic simulations reveal that: (1) mass lost as ejecta can be a considerable fraction of the total mass in the system (up to $\sim25$%); (2) due to asymmetric mass loss, the collision product can sometimes receive a kick velocity that exceeds 10 km/s, large enough to allow the collision product to escape the core of the cluster; and (3) the energy of the ejected matter can be large enough (up to $\sim 3\times 10^{50}$ erg) to remove or disturb the inter cluster gas appreciably. △ Less

Submitted 6 April, 2009; originally announced April 2009.

Comments: 23 pages, 24 figures, Submitted to MNRAS

Mixing in massive stellar mergers

Authors: E. Gaburov, J. C. Lombardi, S. Portegies Zwart

Abstract: The early evolution of dense star clusters is possibly dominated by close interactions between stars, and physical collisions between stars may occur quite frequently. Simulating a stellar collision event can be an intensive numerical task, as detailed calculations of this process require hydrodynamic simulations in three dimensions. We present a computationally inexpensive method in which we ap… ▽ More The early evolution of dense star clusters is possibly dominated by close interactions between stars, and physical collisions between stars may occur quite frequently. Simulating a stellar collision event can be an intensive numerical task, as detailed calculations of this process require hydrodynamic simulations in three dimensions. We present a computationally inexpensive method in which we approximate the merger process, including shock heating, hydrodynamic mixing and mass loss, with a simple algorithm based on conservation laws and a basic qualitative understanding of the hydrodynamics of stellar mergers. The algorithm relies on Archimedes' principle to dictate the distribution of the fluid in the stable equilibrium situation. We calibrate and apply the method to mergers of massive stars, as these are expected to occur in young and dense star clusters. We find that without the effects of microscopic mixing, the temperature and chemical composition profiles in a collision product can become double-valued functions of enclosed mass. Such an unphysical situation is mended by simulating microscopic mixing as a post-collision effect. In this way we find that head-on collisions between stars of the same spectral type result in substantial mixing, while mergers between stars of different spectral type, such as type B and O stars ($\sim$10 and $\sim$40\msun respectively), are subject to relatively little hydrodynamic mixing. △ Less

Submitted 2 October, 2007; v1 submitted 20 July, 2007; originally announced July 2007.

Comments: Accepted by MNRAS

Stellar Collisions and Ultracompact X-ray Binary Formation

Authors: J. C. Lombardi Jr., Z. F. Proulx, K. L. Dooley, E. M. Theriault, N. Ivanova, F. A. Rasio

Abstract: (abridged) We report the results of SPH calculations of parabolic collisions between a subgiant or slightly evolved red-giant star and a neutron star (NS). Such collisions are likely to form ultracompact X-ray binaries (UCXBs) observed today in old globular clusters. In particular, we compute collisions of a 1.4 Msun NS with realistically modelled parent stars of initial masses 0.8 and 0.9 Msun,… ▽ More (abridged) We report the results of SPH calculations of parabolic collisions between a subgiant or slightly evolved red-giant star and a neutron star (NS). Such collisions are likely to form ultracompact X-ray binaries (UCXBs) observed today in old globular clusters. In particular, we compute collisions of a 1.4 Msun NS with realistically modelled parent stars of initial masses 0.8 and 0.9 Msun, each at three different evolutionary stages (corresponding to three different radii R). The distance of closest approach for the initial orbit varies from 0.04 R (nearly head-on) to 1.3 R (grazing). These collisions lead to the formation of a tight binary, composed of the NS and the subgiant or red-giant core, embedded in an extremely diffuse common envelope (CE) typically of mass ~0.1 to 0.3 Msun. Our calculations follow the binary for many hundreds of orbits, ensuring that the orbital parameters we determine at the end of the calculations are close to final. Some of the fluid initially in the envelope of the (sub)giant, from 0.003 to 0.023 Msun in the cases we considered, is left bound to the NS. The eccentricities of the resulting binaries range from about 0.2 for our most grazing collision to about 0.9 for the nearly head-on cases. In almost all the cases we consider, gravitational radiation alone will cause sufficiently fast orbital decay to form a UCXB within a Hubble time, and often on a much shorter timescale. Our hydrodynamics code implements the recent SPH equations of motion derived with a variational approach by Springel & Hernquist and by Monaghan. Numerical noise is reduced by enforcing an analytic constraint equation that relates the smoothing lengths and densities of SPH particles. We present tests of these new methods to help demonstrate their improved accuracy. △ Less

Submitted 23 November, 2005; v1 submitted 17 September, 2005; originally announced September 2005.

Comments: 41 pages, 17 figures, accepted by ApJ

Journal ref: Astrophys.J.640:441-458,2006

Formation of Ultracompact X-ray Binaries in Dense Star Clusters

Authors: N. Ivanova, F. A. Rasio, J. C. Lombardi, Jr., K. L. Dooley, Z. F. Proulx

Abstract: Bright, ultracompact X-ray binaries observed in dense star clusters, such as Galactic globular clusters, must have formed relatively recently, since their lifetimes as persistent bright sources are short (e.g., ~10^8 yr above 10^36 erg/s for a 1.4 Msun neutron star accreting from a degenerate helium companion with an initial mass of ~0.2 Msun). Therefore, we can use the present conditions in a c… ▽ More Bright, ultracompact X-ray binaries observed in dense star clusters, such as Galactic globular clusters, must have formed relatively recently, since their lifetimes as persistent bright sources are short (e.g., ~10^8 yr above 10^36 erg/s for a 1.4 Msun neutron star accreting from a degenerate helium companion with an initial mass of ~0.2 Msun). Therefore, we can use the present conditions in a cluster core to study possible dynamical formation processes for these sources. Here we show that direct physical collisions between neutron stars and red giants can provide a sufficient formation rate to explain the observed numbers of bright sources. These collisions produce tight, eccentric neutron star -- white dwarf binaries that decay to contact by gravitational radiation on timescales ~10^6-10^10 yr, usually shorter and often much shorter than the cluster age. △ Less

Submitted 10 February, 2005; v1 submitted 27 January, 2005; originally announced January 2005.

Comments: 4 pages, 1 figure, in emulateapj style. ApJ Letters, in press, includes corrections in proofs

Journal ref: Astrophys.J. 621 (2005) L109-L112

Hydrogen Spectra, Molecular Association and Orbital Radii in the Solar System

Authors: James C. Lombardi Sr.

Abstract: A relationship between the average orbital radii of the planets and their satellites in the solar system and the spectra of atomic and molecular hydrogen is identified and investigated. In this model, stimulated radiative association resonances develop early on in the disk of the protosun that cause the disk to cool at only certain radii, with each radius depending on a specific photon energy in… ▽ More A relationship between the average orbital radii of the planets and their satellites in the solar system and the spectra of atomic and molecular hydrogen is identified and investigated. In this model, stimulated radiative association resonances develop early on in the disk of the protosun that cause the disk to cool at only certain radii, with each radius depending on a specific photon energy in the atomic hydrogen spectrum. The planets then evolve from the relatively cool rings that are formed. Similar activity occurs in the formation of the satellite systems of the giant planets. The present investigation deals with the mechanism that generates rings from which the planets are formed. It does not deal with the evolution of the rings into planets. Many characteristics of the solar system are explained including the sizes of the orbital radii of the planets and their satellites, the tilt of Uranus's axis, the positions of the asteroid and Kuiper belts, the source of the scattered Kuiper belt objects, the positions of Saturn's main rings and the rings of Uranus, Jupiter, and Neptune. It also shows that a commonality exists in the structures of the solar system and the planetary systems that can be attributed to the common process that initiated their evolution. △ Less

Submitted 24 July, 2003; originally announced July 2003.

Comments: 44 pages including 8 tables and 8 figures

Modelling Collision Products of Triple-Star Mergers

Authors: James C. Lombardi, Ashley P. Thrall, Julia S. Deneva, Scott W. Fleming, Paul E. Grabowski

Abstract: In dense stellar clusters, binary-single and binary-binary encounters can ultimately lead to collisions involving two or more stars. A comprehensive survey of multi-star collisions would need to explore an enormous amount of parameter space, but here we focus on a number of representative cases involving low-mass main-sequence stars. Using both Smoothed Particle Hydrodynamics (SPH) calculations… ▽ More In dense stellar clusters, binary-single and binary-binary encounters can ultimately lead to collisions involving two or more stars. A comprehensive survey of multi-star collisions would need to explore an enormous amount of parameter space, but here we focus on a number of representative cases involving low-mass main-sequence stars. Using both Smoothed Particle Hydrodynamics (SPH) calculations and a much faster fluid sorting software package (MMAS), we study scenarios in which a newly formed product from an initial collision collides with a third parent star. By varying the order in which the parent stars collide, as well as the orbital parameters of the collision trajectories, we investigate how factors such as shock heating affect the chemical composition and structure profiles of the collision product. Our simulations and models indicate that the distribution of most chemical elements within the final product is not significantly affected by the order in which the stars collide, the direction of approach of the third parent star, or the periastron separations of the collisions. We find that the sizes of the products, and hence their collisional cross sections for subsequent encounters, are sensitive to the order and geometry of the collisions. For the cases that we consider, the radius of the product formed in the first (single-single star) collision ranges anywhere from roughly 2 to 30 times the sum of the radii of its parent stars. The final product formed in our triple-star collisions can easily be as large or larger than a typical red giant. We therefore expect the collisional cross section of a newly formed product to be greatly enhanced over that of a thermally relaxed star of the same mass. △ Less

Submitted 2 July, 2003; originally announced July 2003.

Comments: 20 pages, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.345:762,2003

Modest-2: A Summary

Authors: A. Sills, S. Deiters, P. Eggleton, M. Freitag, M. Giersz, D. Heggie, J. Hurley, P. Hut, N. Ivanova, R. Klessen, P. Kroupa, J. Lombardi, S. McMillan, S. Portegies Zwart, H. Zinnecker

Abstract: This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the… ▽ More This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified. △ Less

Submitted 11 March, 2003; v1 submitted 23 January, 2003; originally announced January 2003.

Comments: Accepted by New Astronomy

Journal ref: New Astron. 8 (2003) 605-628

MODEST-1: Integrating Stellar Evolution and Stellar Dynamics

Authors: Piet Hut, Michael M. Shara, Sverre J. Aarseth, Ralf S. Klessen, James C. Lombardi Jr., Junichiro Makino, Steve McMillan, Onno R. Pols, Peter J. Teuben, Ronald F. Webbink

Abstract: We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evol… ▽ More We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics. △ Less

Submitted 31 October, 2002; v1 submitted 15 July, 2002; originally announced July 2002.

Comments: workshop summary, to appear in New Astronomy, 36 pages, 3 figures; this slightly revised version mainly contains more references

Journal ref: New Astron. 8 (2003) 337-370

Stellar Collisions and the Interior Structure of Blue Stragglers

Authors: James C. Lombardi, Jessica Sawyer Warren, Frederic A. Rasio, Alison Sills, Aaron R. Warren

Abstract: Collisions of main sequence stars occur frequently in dense star clusters. In open and globular clusters, these collisions produce merger remnants that may be observed as blue stragglers. Detailed theoretical models of this process require lengthy hydrodynamic computations in three dimensions. However, a less computationally expensive approach, which we present here, is to approximate the merger… ▽ More Collisions of main sequence stars occur frequently in dense star clusters. In open and globular clusters, these collisions produce merger remnants that may be observed as blue stragglers. Detailed theoretical models of this process require lengthy hydrodynamic computations in three dimensions. However, a less computationally expensive approach, which we present here, is to approximate the merger process (including shock heating, hydrodynamic mixing, mass ejection, and angular momentum transfer) with simple algorithms based on conservation laws and a basic qualitative understanding of the hydrodynamics. These algorithms have been fine tuned through comparisons with the results of our previous hydrodynamic simulations. We find that the thermodynamic and chemical composition profiles of our simple models agree very well with those from recent SPH (smoothed particle hydrodynamics) calculations of stellar collisions, and the subsequent stellar evolution of our simple models also matches closely that of the more accurate hydrodynamic models. Our algorithms have been implemented in an easy to use software package, which we are making publicly available (see http://vassun.vassar.edu/~lombardi/mmas/). This software could be used in combination with realistic dynamical simulations of star clusters that must take into account stellar collisions. △ Less

Submitted 2 November, 2001; v1 submitted 19 July, 2001; originally announced July 2001.

Comments: This revised version has 37 pages, 13 figures, 4 tables; submitted to ApJ; for associated software package, see http://vassun.vassar.edu/~lombardi/mmas/ This revised version presents additional comparisons with SPH results and slightly improved merger recipes

Stellar Collisions and Blue Straggler Formation

Authors: James C. Lombardi, Frederic A. Rasio

Abstract: We review recent 3D hydrodynamic calculations of stellar collisions using the Smoothed Particle Hydrodynamics (SPH) method, and we discuss the implications of the results for the formation and evolution of blue stragglers in globular clusters. We also discuss the construction of simple analytic models for merger remnants, approximating the mass loss, shock heating, mixing, and angular momentum t… ▽ More We review recent 3D hydrodynamic calculations of stellar collisions using the Smoothed Particle Hydrodynamics (SPH) method, and we discuss the implications of the results for the formation and evolution of blue stragglers in globular clusters. We also discuss the construction of simple analytic models for merger remnants, approximating the mass loss, shock heating, mixing, and angular momentum transport during a collision with simple algorithms that can be calibrated using our 3D hydrodynamic results. The thermodynamic and chemical composition profiles predicted by these simple models are compared with those from our SPH simulations, demonstrating that our new models provide accurate representations of true collisional merger remnants. △ Less

Submitted 5 January, 2001; originally announced January 2001.

Comments: 29 pages, 7 figures, to appear in "Stellar Collisions, Mergers and their Consequences," ed. M. Shara (ASP Conference Series), for associated mpeg files, see http://research.amnh.org/~summers/scviz/lombardi.html

Evolution of Stellar Collision Products in Globular Clusters - II. Off-axis Collision

Authors: Alison Sills, Joshua A. Faber, James C. Lombardi, Jr., Frederic A. Rasio, Aaron Warren

Abstract: We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitte… ▽ More We continue our exploration of collisionally merged stars in the blue straggler region of the color-magnitude diagram. We report the results of new SPH calculations of parabolic collisions between two main-sequence stars, with the initial structure and composition profiles of the parent stars having been determined from stellar evolution calculations. Parallelization of the SPH code has permitted much higher numerical resolution of the hydrodynamics. We also present evolutionary tracks for the resulting collision products, which emerge as rapidly rotating blue stragglers. The rotating collision products are brighter, bluer and remain on the main sequence longer than their non-rotating counterparts. In addition, they retain their rapid rotation rates throughout their main sequence lifetime. Rotationally-induced mixing strongly affects the evolution of the collision products, although it is not sufficient to mix the entire star. We discuss the implications of these results for studies of blue straggler populations in clusters. This work shows that off-axis collision products cannot become blue stragglers unless they lose a large fraction of their initial angular momentum. The mechanism for this loss is not apparent, although some possibilities are discussed. △ Less

Submitted 16 August, 2000; originally announced August 2000.

Comments: 25 pages incl. 9 figures (one in colour). Submitted to ApJ

Tests of Spurious Transport in Smoothed Particle Hydrodynamics

Authors: James C. Lombardi, Alison Sills, Frederic A. Rasio, Stuart L. Shapiro

Abstract: We have performed a series of systematic tests to evaluate the effects of spurious transport in three-dimensional smoothed particle hydrodynamics (SPH) calculations. Our tests investigate (i) particle diffusion, (ii) shock heating, (iii) numerical viscosity, and (iv) angular momentum transport. The results are useful for quantifying the accuracy of the SPH scheme, especially for problems where s… ▽ More We have performed a series of systematic tests to evaluate the effects of spurious transport in three-dimensional smoothed particle hydrodynamics (SPH) calculations. Our tests investigate (i) particle diffusion, (ii) shock heating, (iii) numerical viscosity, and (iv) angular momentum transport. The results are useful for quantifying the accuracy of the SPH scheme, especially for problems where shear flows or shocks are present, as well as for problems where true hydrodynamic mixing is relevant. We examine the different forms of artificial viscosity (AV) which have been proposed by Monaghan, by Hernquist & Katz, and by Balsara. For each form, our tests suggest a single set of values for the AV parameters $α$ and $β$ (coefficients of the linear and quadratic terms) which are appropriate in a large number of situations. We also discuss how these parameters should be adjusted depending on the goals of the particular application. We find that both the Hernquist & Katz and Balsara forms introduce relatively small amounts of numerical viscosity. Furthermore, both Monaghan's and Balsara's AV do well at treating shocks and at limiting the amount of spurious mixing. For these reasons, we endorse the Balsara AV for use in a broad range of applications. △ Less

Submitted 28 July, 1998; originally announced July 1998.

Comments: 49 pages, 26 figures as 30 postscript files, submitted to The Journal of Computational Physics

Smoothed Particle Hydrodynamics Calculations of Stellar Interactions

Authors: Frederic A. Rasio, James C. Lombardi

Abstract: Smoothed Particle Hydrodynamics is a multidimensional Lagrangian method of numerical hydrodynamics that has been used to tackle a wide variety of problems in astrophysics. Here we develop the basic equations of the SPH scheme, and we discuss some of its numerical properties and limitations. As an illustration of typical astrophysical applications, we discuss recent calculations of stellar intera… ▽ More Smoothed Particle Hydrodynamics is a multidimensional Lagrangian method of numerical hydrodynamics that has been used to tackle a wide variety of problems in astrophysics. Here we develop the basic equations of the SPH scheme, and we discuss some of its numerical properties and limitations. As an illustration of typical astrophysical applications, we discuss recent calculations of stellar interactions, including collisions between main sequence stars and the coalescence of compact binaries. △ Less

Submitted 12 January, 2000; v1 submitted 7 May, 1998; originally announced May 1998.

Comments: 21 pages, invited review paper to appear in Journal of Computational and Applied Mathematics (JCAM). Fig. 2 is available from the authors in hardcopy form only. Revised, slightly abbreviated version

Journal ref: J. Comp. App. Math., 109 (1999) 213-230

Post-Newtonian Models of Binary Neutron Stars

Authors: James C. Lombardi, Frederic A. Rasio, Stuart L. Shapiro

Abstract: Using an energy variational method, we calculate quasi-equilibrium configurations of binary neutron stars modeled as compressible triaxial ellipsoids obeying a polytropic equation of state. Our energy functional includes terms both for the internal hydrodynamics of the stars and for the external orbital motion. We add the leading post-Newtonian (PN) corrections to the internal and gravitational… ▽ More Using an energy variational method, we calculate quasi-equilibrium configurations of binary neutron stars modeled as compressible triaxial ellipsoids obeying a polytropic equation of state. Our energy functional includes terms both for the internal hydrodynamics of the stars and for the external orbital motion. We add the leading post-Newtonian (PN) corrections to the internal and gravitational energies of the stars, and adopt hybrid orbital terms which are fully relativistic in the test-mass limit and always accurate to PN order. The total energy functional is varied to find quasi-equilibrium sequences for both corotating and irrotational binaries in circular orbits. We examine how the orbital frequency at the innermost stable circular orbit depends on the polytropic index n and the compactness parameter GM/Rc^2. We find that, for a given GM/Rc^2, the innermost stable circular orbit along an irrotational sequence is about 17% larger than the innermost secularly stable circular orbit along the corotating sequence when n=0.5, and 20% larger when n=1. We also examine the dependence of the maximum neutron star mass on the orbital frequency and find that, if PN tidal effects can be neglected, the maximum equilibrium mass increases as the orbital separation decreases. △ Less

Submitted 20 June, 1997; v1 submitted 27 May, 1997; originally announced May 1997.

Comments: 53 pages, LaTex, 9 figures as 10 postscript files, accepted by Phys. Rev. D, replaced version contains updated references

Journal ref: Phys.Rev. D56 (1997) 3416-3438

The Importance of Realistic Starting Models for Hydrodynamic Simulations of Stellar Collisions

Authors: Alison Sills, James C. Lombardi

Abstract: We demonstrate the necessity of using realistic stellar models taken from stellar evolution codes, as opposed to polytropes, for starting models in smoothed particle hydrodynamics calculations of collisions between main sequence stars. Evolved stars have mean molecular weight gradients, which affect their entropy profiles and therefore affect how they react during a collision. The structure of s… ▽ More We demonstrate the necessity of using realistic stellar models taken from stellar evolution codes, as opposed to polytropes, for starting models in smoothed particle hydrodynamics calculations of collisions between main sequence stars. Evolved stars have mean molecular weight gradients, which affect their entropy profiles and therefore affect how they react during a collision. The structure of stellar collision products of polytrope parent stars is significantly different from that of collision products of realistic parent models. These differences strongly affect the future evolution of the collision products, particularly products of collisions between unequal mass stars which have undergone significant chemical evolution. The use of polytropes as parent star models is likely to result in qualitatively mistaken results for the structure of the collision product. △ Less

Submitted 5 May, 1997; originally announced May 1997.

Comments: 13 pages, 4 encapsulated postscript figures, accepted by Astrophysical Journal (Letters)

Evolution of Stellar Collision Products in Globular Clusters -- I. Head-on Collisions

Authors: Alison Sills, James C. Lombardi, Charles D. Bailyn, Pierre Demarque, Frederic A. Rasio, Stuart L. Shapiro

Abstract: We explore the evolution of collisionally merged stars in the blue straggler region of the HR diagram. The starting models for our stellar evolution calculations are the results of the smoothed particle hydrodynamics (SPH) simulations of parabolic collisions between main sequence stars performed by Lombardi, Rasio and Shapiro (1996). Since SPH and stellar evolution codes employ different and oft… ▽ More We explore the evolution of collisionally merged stars in the blue straggler region of the HR diagram. The starting models for our stellar evolution calculations are the results of the smoothed particle hydrodynamics (SPH) simulations of parabolic collisions between main sequence stars performed by Lombardi, Rasio and Shapiro (1996). Since SPH and stellar evolution codes employ different and often contradictory approximations, it is necessary to treat the evolution of these products carefully. The mixture and disparity of the relevant timescales (hydrodynamic, thermal relaxation and nuclear burning) and of the important physical assumptions between the codes makes the combined analysis of the problem challenging, especially during the initial thermal relaxation of the star. In particular, the treatment of convection is important, and semiconvection must be modeled in some detail. The products of seven head-on collisions are evolved through their initial thermal relaxation, and then through the main sequence phase to the base of the giant branch. Their evolutionary tracks are presented. In contrast to the assumptions in previous work, these collision products do not develop substantial convective regions during their thermal relaxation, and therefore are not mixed significantly after the collision. △ Less

Submitted 5 May, 1997; originally announced May 1997.

Comments: 45 pages, 8 figures (encapsulated postscript), uses AASTeX, to be published in Astrophysical Journal volume 487

Collisions of Main-Sequence Stars and the Formation of Blue Stragglers in Globular Clusters

Authors: J. C. Lombardi, F. A. Rasio, S. L. Shapiro

Abstract: We report the results of new SPH calculations of parabolic collisions between two main-sequence stars in a globular cluster. Such collisions are directly relevant to the formation of blue stragglers. In particular, we consider parent stars of mass $M/M_{TO}=0.2,0.5,0.75$, and $1$, where $M_{TO}$ is the cluster turnoff mass (typically about $0.8\,M_\odot$). Low-mass stars (with $M=0.2 M_{TO}$ or… ▽ More We report the results of new SPH calculations of parabolic collisions between two main-sequence stars in a globular cluster. Such collisions are directly relevant to the formation of blue stragglers. In particular, we consider parent stars of mass $M/M_{TO}=0.2,0.5,0.75$, and $1$, where $M_{TO}$ is the cluster turnoff mass (typically about $0.8\,M_\odot$). Low-mass stars (with $M=0.2 M_{TO}$ or $0.5 M_{TO}$) are assumed to be fully convective and are therefore modeled as $n=1.5$ polytropes. Stars at the turnoff (with $M= M_{TO}$) are assumed to be mostly radiative and are modeled as $n=3$ polytropes. Intermediate-mass stars (with $M=0.75 M_{TO}$) are modeled as composite polytropes consisting of a radiative core with polytropic index $n=3$ and a convective envelope with $n=1.5$. We focus our study on the question of hydrodynamic mixing of helium and hydrogen, which plays a crucial role in determining the observable characteristics of blue stragglers. In all cases we find that there is negligible hydrodynamic mixing of helium into the outer envelope of the merger remnant. The amount of hydrogen mixed into the core of the merger depends strongly on the entropy profiles of the two colliding stars. In contrast to what has been done in previous studies, blue stragglers formed by direct stellar collisions should not necessarily be assumed to have initially homogeneous composition profiles. △ Less

Submitted 16 November, 1995; originally announced November 1995.

Comments: uses AAS latex macros v3, 66 pages with tables and figures, submitted to ApJ, complete postscript also available at http://ensor.mit.edu/~rasio/papers

ON BLUE STRAGGLER FORMATION BY DIRECT COLLISIONS OF MAIN-SEQUENCE STARS

Authors: James C. Lombardi, Jr., Frederic A. Rasio, Stuart L. Shapiro

Abstract: We report the results of new SPH calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close to the MS turn-off point in a globular cluster and are therefore modeled as $n=3$, $Γ=5/3$ polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic m… ▽ More We report the results of new SPH calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close to the MS turn-off point in a globular cluster and are therefore modeled as $n=3$, $Γ=5/3$ polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic mixing occurs between the dense, helium-rich inner cores and the outer envelopes. As a result, and in contrast to what has been assumed in previous studies, blue stragglers formed by direct stellar collisions are not necessarily expected to have anomalously high helium abundances in their envelopes or to have their cores replenished with fresh hydrogen fuel. △ Less

Submitted 28 February, 1995; originally announced February 1995.

Comments: uuencoded compressed postscript with figures, 10 pages (Figure 1 available upon request from best@sns.ias.edu)

Report number: IAS-SNS-AST 95/17

Journal ref: Astrophys.J. 445 (1995) L117