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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 ($\Gamma = 5/2,~n=3/2$; $\Gamma=3/2,~n=2$; $\Gamma=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.

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.

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.

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.

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 $\theta\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 $\theta\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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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 $\alpha$ and $\beta$ (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.

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.

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.

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.

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.