The origin of chemically peculiar stars and non-zero eccentricity in evolved close binaries have been long-standing problems in binary stellar evolution. Answers to these questions may trace back to an intense mass transfer during AGB binary phase. In this work, we use AstroBEAR to solve the 3D radiation-hydrodynamic equations and calculate the mass transfer rate in asymptotic-giant-branch (AGB) binaries that undergo the wind-Roche-lobe-overflow or Bondi-Hoyle-Lyttleton (BHL) accretion. MESA produces the density and temperature of the boundary condition of the AGB star. To improve the resolution of the dynamics of a circumbinary disk, we implement an azimuthal angle-dependent 3D radiation transfer. We consider optically thin cooling and obtain the number density of the coolants by solving Saha equations. One of the goals of this work is to illustrate the transition from the wind-Roche-lobe-overflow to BHL accretion. Both circumbinary disks and spiral structure outflows can appear in the simulations. Circumbinary disks may form when the optical thickness in the equatorial region increases. The increase of the optical thickness is due to the deflected wind. The resulting mass transfer efficiency in our models is up to a factor of 8 times higher than what the standard BHL accretion scenario predicts, and the outflow gains up to 91% of its initial angular momentum when it reaches 1.3 binary separations. Consequently, some AGB binaries may undergo orbit shrinkage, and some will expand. The high mass transfer efficiency is closely related to the presence of the circumbinary disks.

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.

A class of optical transients known as Luminous Red Novae (LRNe) have recently been associated with mass ejections from binary stars undergoing common-envelope evolution. We use the population synthesis code COMPAS to explore the impact of a range of assumptions about the physics of common-envelope evolution on the properties of LRNe. In particular, we investigate the influence of various models for the energetics of LRNe on the expected event rate and light curve characteristics, and compare with the existing sample. We find that the Galactic rate of LRNe is $\sim 0.2$ yr$^{-1}$, in agreement with the observed rate. In our models, the luminosity function of Galactic LRNe covers multiple decades in luminosity and is dominated by signals from stellar mergers, consistent with observational constraints from iPTF and the Galactic sample of LRNe. We discuss how observations of the brightest LRNe may provide indirect evidence for the existence of massive ($> 40$ M$_\odot$) red supergiants. Such LRNe could be markers along the evolutionary pathway leading to the formation of double compact objects. We make predictions for the population of LRNe observable in future transient surveys with the Large Synoptic Survey Telescope and the Zwicky Transient Facility. In all plausible circumstances, we predict a selection-limited observable population dominated by bright, long-duration events caused by common envelope ejections. We show that the Large Synoptic Survey Telescope will observe $20$--$750$ LRNe per year, quickly constraining the luminosity function of LRNe and probing the physics of common-envelope events.

The formation of low-mass X-ray binaries (LMXBs) is an ongoing challenge in stellar evolution. The important subset of LMXBs are the binary systems with a neutron star (NS) accretor. In NS LMXBs with non-degenerate donors, the mass transfer is mainly driven by magnetic braking. The discrepancies between the observed mass transfer (MT) rates and the theoretical models were known for a while. Theory predictions of the MT rates are too weak and differ by an order of magnitude or more. Recently, we showed that with the standard magnetic braking, it is not possible to find progenitor binary systems such that they could reproduce -- at any time of their evolution -- most of the observed persistent NS LMXBs. In this ${\it Letter}$ we present a modified magnetic braking prescription, CARB (Convection And Rotation Boosted). CARB magnetic braking combines two recent improvements in understanding stellar magnetic fields and magnetized winds -- the dependence of the magnetic field strength on the outer convective zone and the dependence of the Alfv\`en radius on the donor's rotation. Using this new magnetic braking prescription, we can reproduce the observed mass transfer rates at the detected mass ratio and orbital period for all well-observed to-the-date Galactic persistent NS LMXBs. For the systems where the effective temperature of the donor stars is known, theory agrees with observations as well.

We present a population study of low- and intermediate-mass X-ray binaries (LMXBs) with neutron star accretors, performed using the detailed 1D stellar evolution code MESA. We identify all plausible Roche-lobe overflowing binaries at the start of mass transfer, and compare our theoretical mass transfer tracks to the population of well-studied Milky Way LMXBs. The mass transferring evolution depends on the accepted magnetic braking (MB) law for angular momentum loss. The most common MB prescription ("Skumanich MB") originated from observations of the time-dependence of rotational braking of Sun-type stars, where the angular momentum loss rate depends on the donor mass $M_d$, donor radius $R_d$, and rotation rate $\Omega$, $\dot{J} \propto M_d R_d^{\gamma} \Omega^3$. The functional form of the Skumanich MB can be also obtained theoretically assuming a radial magnetic field, isotropic isothermal winds, and boosting of the magnetic field by rotation. Here we show that this simple form of the Skumanich MB law gives mass transfer rates an order of magnitude too weak to explain most observed persistent LMXBs. This failure suggests that the standard Skumanich MB law should not be employed to interpret Galactic, or extragalactic, LMXB populations, with either detailed stellar codes or rapid binary population synthesis codes. We investigate modifications for the MB law, and find that including a scaling of the magnetic field strength with the convective turnover time, and a scaling of MB with the wind mass loss rate, can reproduce persistent LMXBs, and does a better job at reproducing transient LMXBs.

In this Letter we discuss what happens to hydrogen recombination energy that is released during regular common envelope (CE) events, as opposed to self-regulated CE events. We show that the amount of recombination energy that can be transferred away by either convection or radiation from the regions where recombination takes place is negligible. Instead, recombination energy is destined to be used either to help CE expansion, as a work term, or to accelerate local fluid elements. The exceptions are donors that initially have very high entropy material, $S/(k_B N_A)>37$ mol g$^{-1}$. The analysis and conclusions are independent of specific stellar models or evolutionary codes, and rely on fundamental properties of stellar matter such as the equation of state, Saha equation and opacities, as well as on stellar structure equations and the mixing length theory of convection.

We present 22 new (+3 confirmed) cataclysmic variables (CVs) in the non core-collapsed globular cluster 47 Tucanae (47 Tuc). The total number of CVs in the cluster is now 43, the largest sample in any globular cluster so far. For the identifications we used near-ultraviolet (NUV) and optical images from the Hubble Space Telescope, in combination with X-ray results from the Chandra X-ray Observatory. This allowed us to build the deepest NUV CV luminosity function of the cluster to date. We found that the CVs in 47 Tuc are more concentrated towards the cluster center than the main sequence turnoff stars. We compared our results to the CV populations of the core-collapsed globular clusters NGC 6397 and NGC 6752. We found that 47 Tuc has fewer bright CVs per unit mass than those two other clusters. That suggests that dynamical interactions in core-collapsed clusters play a major role creating new CVs. In 47 Tuc, the CV population is probably dominated by primordial and old dynamically formed systems. We estimated that the CVs in 47 Tuc have total masses of approx. 1.4 M_sun. We also found that the X-ray luminosity function of the CVs in the three clusters is bimodal. Additionally, we discuss a possible double degenerate system and an intriguing/unclassified object. Finally, we present four systems that could be millisecond pulsar companions given their X-ray and NUV/optical colors. For one of them we present very strong evidence for being an ablated companion. The other three could be CO- or He-WDs.

Dynamical interactions that take place between objects in dense stellar systems lead to frequent formation of exotic stellar objects, unusual binaries, and systems of higher multiplicity. They are most important for the formation of binaries with neutron stars and black holes, which are usually observationally revealed in mass-transferring binaries. Here we review the current understanding of compact object's retention, of the metallicity dependence on the formation of low-mass X-ray binaries with neutron stars, and how mass-transferring binaries with a black hole and a white dwarf can be formed. We discuss as well one old unsolved puzzle and two new puzzles posed by recent observations: what descendants do ultra-compact X-ray binaries produce, how are very compact triples formed, and how can black hole low-mass X-ray binaries acquire non-degenerate companions?

In this Letter we propose a formation channel for low-mass X-ray binaries with black hole accretors and non-degenerate donors via grazing tidal encounters with subgiants. We estimate that in a typically dense globular cluster with the core density of $10^5$ stars pc$^{-3}$, the formation rates are about one binary per Gyr per 50-100 retained black holes. The donors -- stripped subgiants -- will be strongly underluminous when compared to subgiant or giant branch stars of the same colors. The products of tidal stripping are underluminous by at least one magnitude for several hundred million years, when compared to normal stars of the same color, and differ from underluminous red stars that could be produced by non-catastrophic mass transfer in an ordinary binary. The dynamically formed binaries become quiescent low-mass X-ray binaries, with lifetimes of about a Gyr. The expected number of X-ray binaries is one per 50-200 retained black holes, while the expected number of strongly underluminous subsubgiant is about half this. The presence of strongly underluminous stars in a globular cluster may be indicative of black holes presence.

We review the fundamentals and the recent developments in understanding of common envelope physics. We report specifically on the progress that was made by the consideration of the recombination energy. This energy is found to be responsible for the complete envelope ejection in the case of a prompt binary formation, for the delayed dynamical ejections in the case of a self-regulated spiral-in, and for the steady recombination outflows during the transition between the plunge-in and the self-regulated spiral-in. Due to different ways how the recombination affects the common envelope during fast and slow spiral-ins, the apparent efficiency of the orbital energy use can be different between the two types of spiral-ins by a factor of ten. We also discuss the observational signatures of the common envelope events, their link a new class of astronomical transients, Luminous Red Novae, and to a plausible class of very luminous irregular variables.

After the initial fast spiral-in phase experienced by a common-envelope binary, the system may enter a slow, self-regulated phase, possibly lasting 100s of years, in which all the energy released by orbital decay can be efficiently transported to the surface, where it is radiated away. If the remaining envelope is to be removed during this phase, this removal must occur through some as-yet-undetermined mechanism. We carried out 1-d hydrodynamic simulations of a low-mass red giant undergoing a synthetic common-envelope event in such a slow spiral-in phase, using the stellar evolutionary code MESA. We simulated the heating of the envelope due to frictional dissipation from a binary companion's orbit in multiple configurations and investigated the response of the giant's envelope. We find that our model envelopes become dynamically unstable and develop large-amplitude pulsations, with periods in the range 3-20 years and very short growth time-scales of similar order. The shocks and associated rebounds that emerge as these pulsations grow are in some cases strong enough to dynamically eject shells of matter of up to 0.1 $\mathrm{M}_{\odot}$, $\sim 10$ % of the mass of the envelope, from the stellar surface at above escape velocity. These ejections are seen to repeat within a few decades, leading to a time-averaged mass-loss rate of order $10^{-3}$ $\mathrm{M}_{\odot} \: \mathrm{yr}^{-1}$ which is sufficiently high to represent a candidate mechanism for removing the entire envelope over the duration of the slow spiral-in phase.

We present a three-dimensional (3D) study of common envelope events (CEEs) to provide a foundation for future one-dimensional (1D) methods to model the self-regulated phase of a CEE. The considered CEEs with a low-mass red giant end with one of three different outcomes -- merger, slow spiral-in, or prompt formation of a binary. To understand which physical processes determine different outcomes, and to evaluate how well 1D simulations model the self-regulated phase of a CEE, we introduce tools that map our 3D models to 1D profiles. We discuss the differences in the angular momentum and energy redistribution in 1D and 3D codes. We identified four types of ejection processes: the pre-plunge-in ejection, the outflow during the plunge-in, the outflow driven by recombination, and the ejection triggered by a contraction of the circumbinary envelope. Significant mass is lost in all cases, including the mergers. Therefore a self-regulated spiral-in can start only with a strongly reduced envelope mass. We derive the condition to start a recombination outflow, which can proceed either as a runaway or a stationary outflow. We show that the way the energy of the inspiraling companion is added to the envelope in 1D studies intensifies the envelope's entropy increase, alters the start of the recombination outflow, and leads to different outcomes in 1D and 3D studies. The steady recombination outflow may dispel most of the envelope in all slow spiral-in cases, making the existence of a long-term self-regulated phase debatable, at least for low-mass giant donors.

The mass transfer in binaries with massive donors and compact companions, when the donors rapidly evolve after their main sequence, is one of the dominant formation channels of merging double stellar-mass black hole binaries. This mass transfer was previously postulated to be unstable and was expected to lead to a common envelope event. The common envelope event then would end with either double black hole formation, or with the merger of the two stars. We re-visit the stability of this mass transfer, and find that for a large range of the binary orbital separations this mass transfer is stable. This newly found stability allows us to reconcile the theoretical rate for double black hole binary mergers predicted by population synthesis studies, and the empirical rate obtained by LIGO. Futhermore, the stability of the mass transfer leads to the formation of ultra-luminous X-ray sources. The theoretically predicted formation rates of ultra-luminous X-ray sources powered by a stellar-mass BH, as well as the range of produced X-ray luminosity, can explain the observed bright ultra-luminous X-ray sources.

Common envelope events are important interactions between two binary stars that lead to the formation of close binary systems. We present here a systematic three-dimensional study in which we model common envelope events with low-mass giant donors. The results allow us to revise the energy formalism that is usually used to determine common envelope event outcomes. We show that the energy budget for this type of system should include the recombination energy, and that it also must take into account that a significant fraction of the released orbital energy is taken away by the ejecta. We provide three ways in which our results can be used by binary population synthesis studies: a relation that links the observed post-common envelope binary with the initial binary parameters, a fitting formula for the $\alpha_{\rm ce}\lambda$ parameter of the standard energy formalism, and a revised energy formalism that takes into account both the recombination energy and the energy that is taken away by the ejecta.

Sco X-1 is a low-mass X-ray binary (LMXB) that has one of the most precisely determined set of binary parameters such as the mass accretion rate, companions mass ratio and the orbital period. For this system, as well as for a large fraction of other well-studied LMXBs, the observationally-inferred mass accretion rate is known to strongly exceed the theoretically expected mass transfer rate. We suggest that this discrepancy can be solved by applying a modified magnetic braking prescription, which accounts for increased wind mass loss in evolved stars compared to main sequence stars. Using our mass transfer framework based on {\tt MESA}, we explore a large range of binaries at the onset of the mass transfer. We identify the subset of binaries for which the mass transfer tracks cross the Sco X-1 values for the mass ratio and the orbital period. We confirm that no solution can be found for which the standard magnetic braking can provide the observed accretion rates, while wind-boosted magnetic braking can provide the observed accretion rates for many progenitor binaries that evolve to the observed orbital period and mass ratio.

We report the discovery of the likely white dwarf companions to radio millisecond pulsars 47 Tuc Q and 47 Tuc S in the globular cluster 47 Tucanae. These blue stars were found in near-ultraviolet images from the Hubble Space Telescope for which we derived accurate absolute astrometry, and are located at positions consistent with the radio coordinates to within 0.016 arcsec (0.2sigma). We present near-ultraviolet and optical colours for the previously identified companion to millisecond pulsar 47 Tuc U, and we unambiguously confirm the tentative prior identifications of the optical counterparts to 47 Tuc T and 47 Tuc Y. For the latter, we present its radio-timing solution for the first time. We find that all five near-ultraviolet counterparts have U300-B390 colours that are consistent with He white dwarf cooling models for masses ~0.16-0.3 Msun and cooling ages within ~0.1-6 Gyr. The Ha-R625 colours of 47 Tuc U and 47 Tuc T indicate the presence of a strong Ha absorption line, as expected for white dwarfs with an H envelope.

The stability of mass transfer in binaries with convective giant donors remains an open question in modern astrophysics. There is a significant discrepancy between what the existing methods predict for a response to mass loss of the giant itself, as well as for the mass transfer rate during the Roche lobe overflow. Here we show that the recombination energy in the superadiabatic layer plays an important and hitherto unaccounted-for role in he donor's response to mass loss, in particular on its luminosity and effective temperature. Our improved optically thick nozzle method to calculate the mass transfer rate via $L_1$ allows us to evolve binary systems for a substantial Roche lobe overflow. We propose a new, strengthened criterion for the mass transfer instability, basing it on whether the donor experiences overflow through its outer Lagrangian point. We find that with the new criterion, if the donor has a well-developed outer convective envelope, the critical initial mass ratio for which a binary would evolve stably through the conservative mass transfer varies from $1.5$ to $2.2$, which is about twice as large as previously believed. In underdeveloped giants with shallow convective envelopes this critical ratio may be even larger. When the convective envelope is still growing, and in particular for most cases of massive donors, the critical mass ratio gradually decreases to this value, from that of radiative donors.

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.

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.

The energy budget in common-envelope events (CEEs) is not well understood, with substantial uncertainty even over to what extent the recombination energy stored in ionised hydrogen and helium might be used to help envelope ejection. We investigate the reaction of a red-giant envelope to heating which mimics limiting cases of energy input provided by the orbital decay of a binary during a CEE, specifically during the post-plunge-in phase during which the spiral-in has been argued to occur on a time-scale longer than dynamical. We show that the outcome of such a CEE depends less on the total amount of energy by which the envelope is heated than on how rapidly the energy was transferred to the envelope and on where the envelope was heated. The envelope always becomes dynamically unstable before receiving net heat energy equal to the envelope's initial binding energy. We find two types of outcome, both of which likely lead to at least partial envelope ejection: "runaway" solutions in which the expansion of the radius becomes undeniably dynamical, and superficially "self-regulated" solutions, in which the expansion of the stellar radius stops but a significant fraction of the envelope becomes formally dynamically unstable. Almost the entire reservoir of initial helium recombination energy is used for envelope expansion. Hydrogen recombination is less energetically useful, but is nonetheless important for the development of the dynamical instabilities. However, this result requires the companion to have already plunged deep into the envelope; therefore this release of recombination energy does not help to explain wide post-common-envelope orbits.

One of the principal mechanisms that is responsible for the origin of blue stragglers is mass transfer that takes place while one of the binary companions overfills its Roche lobe. In this Chapter, we overview the theoretical understanding of mass transfer via Roche lobe overflow: classification, how both the donor and of the accretor respond to the mass transfer on different timescales (adiabatic response, equilibrium response, superadiabatic response, time-dependent response) for different types of their envelopes (convective and radiative). These responses, as well as the assumption on how liberal the process is, are discussed in terms of the stability of the ensuing mass transfer. The predictions of the theory of mass transfer via Roche lobe overflow are then briefly compared with the observed mass-transferring systems with both degenerate and non-degenerate donors. We conclude with the discussion which cases of mass transfer and which primordial binaries could be responsible for blue stragglers formation via Roche lobe overflow, as well as how this can be enhanced for blue stragglers formed in globular clusters.

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.

Ultracompact X-ray binaries (UCXBs) are thought to have relatively simple binary evolution post-contact, leading to clear predictions of their luminosity function. We test these predictions by studying the long-term behavior of known UCXBs in our Galaxy, principally using data from the MAXI All-Sky Survey and the Galactic Bulge Scans with RXTE's PCA instrument. Strong luminosity variations are common (and well-documented) among persistent UCXBs, which requires an explanation other than the disk instability mechanism. We measure the luminosity function of known UCXBs in the Milky Way, which extends to lower luminosities than some proposed theoretical luminosity functions of UCXBs. The difference between field and globular cluster X-ray luminosity functions in other galaxies cannot be explained by an increased fraction of UCXBs in globular clusters. Instead, our measured luminosity function suggests that UCXBs only make up a small fraction of the X-ray binaries above a few times 10^36 ergs/s in both old field populations and globular clusters.

We study the mass transfer rates and disk stability conditions of ultracompact X-ray binaries (UCXBs) using empirical time-averaged X-ray luminosities from Paper I (Cartwright et al. 2013) and compiled information from the literature. The majority of UCXBs are consistent with evolutionary tracks for white dwarf donors. Three UCXBs with orbital periods longer than 40 minutes have mass transfer rates above 10^{-10} Msun/year, inconsistent with white dwarf donor tracks. We show that if helium star donors can retain their initial high entropy, they can explain the observed mass transfer rates of these UCXBs. Several UCXBs show persistent luminosities apparently below the disk instability limit for irradiated He accretion disks. We point out that a predominantly C and/or O disk (as observed in the optical spectra of several) lowers the disk instability limit, explaining this disagreement. The orbital period and low time-averaged mass transfer rate of 2S 0918-549 provide evidence that the donor star is a low-entropy C/O white dwarf, consistent with optical spectra. We combine existing information to constrain the masses of the donors in 4U 1916-053 (0.064+-0.010 Msun) and 4U 1626-67 (<0.036 Msun for a 1.4 Msun neutron star). We show that 4U 1626-67 is indeed persistent, and not undergoing a transient outburst, leaving He star models as the best explanation for the donor.

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.

In dense stellar systems, dynamical interactions between objects lead to frequent formation of exotic stellar objects, unusual binaries, and systems of higher multiplicity. They are especially important for the formation of low mass X-ray binaries (LMXBs), which are not only formed 100 times more efficiently than in the field, but also have a puzzling dependence on metallicity. In this contribution we review how compact objects are formed and retained, the mechanisms of dynamical formation and the specifics of the evolution of mass-transferring binaries with neutron stars and black holes in globular clusters - those two kinds of compact objects have different favored paths to become luminous in X-ray. We describe how stellar evolution affects ostensibly purely dynamical formation, producing the observed metallicity dependence for LMXBs. We also discuss the next puzzle to be solved on our journey to understand the link between LMXBs and millisecond pulsars formation.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

We present comprehensive observational results of the metallicity effect on the fraction of globular clusters (GC) that contain low-mass X-ray binaries (LMXB), by utilizing all available data obtained with Chandra for LMXBs and HST ACS for GCs. Our primary sample consists of old elliptical galaxies selected from the ACS Virgo and Fornax surveys. To improve statistics at both the lowest and highest X-ray luminosity, we also use previously reported results from other galaxies. It is well known that the LMXB fraction is considerably higher in red, metal-rich, than in blue, metal-poor GCs. In this paper, we test whether this metallicity effect is X-ray luminosity-dependent, and find that the effect holds uniformly in a wide luminosity range. This result is statistically significant (at >= 3 sigma) in LMXBs with luminosities in the range LX = 2 x 10^37 - 5 x 10^38 erg s-1, where the ratio of LMXB fractions in metal-rich to metal-poor GCs is R = 3.4 +- 0.5. A similar ratio is also found at lower (down to 10^36 erg s-1) and higher luminosities (up to the ULX regime), but with less significance (~2 sigma confidence). Because different types of LMXBs dominate in different luminosities, our finding requires a new explanation for the metallicity effect in dynamically formed LMXBs. We confirm that the metallicity effect is not affected by other factors such as stellar age, GC mass, stellar encounter rate, and galacto-centric distance.

Ultra-compact X-ray binaries (UCXBs) are accreting systems with periods less than 1 hour, which qualifies them to contain a degenerate donor-companion. One would expect such systems to have the easiest theoretical explanation, compared to other kinds of X-ray binaries. Nonetheless, current theory fails to explain high mass transfer (MT) rates in three recently well observed long-period UCXBs. We find that this range of MT rates can be maintained if the donor is a remnant of an out-of-thermal-equilibrium naked core of a giant which was revealed in a very recent episode of a common envelope (CE) event.

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.

In interacting binaries, comparison of a donor star's radial response to mass loss with the response of its Roche radius determines whether mass loss persists and, if so, determines the timescale and stability of the ensuing evolutionary phase. For giants with deep convective envelopes, the canonical description holds that once mass transfer begins it typically proceeds catastrophically on the dynamical timescale, as the star cannot lose sufficient heat in order to avoid expansion. However, we demonstrate that the local thermal timescale of the envelope's superadiabatic outer surface layer remains comparable to that of mass loss in most cases of "dynamical" mass loss. We argue therefore that if mass loss proceeds on a timescale longer than this, then even a deep convective envelope will not dramatically expand, as the surface layer will have time to relax thermally and reconstitute itself. We demonstrate that in general the polytropic approximation gives much too strict a criterion for stability, and discuss the dependence of the donor's response on its radius in addition to its core mass. In general, we find that the effective response of the donor on rapid timescales cannot be determined accurately without detailed evolutionary calculations.

The common envelope event is one of the most important and uncertain evolutionary stages that lead to formation of compact binaries. While the problem is almost 30 years old, its theoretical foundation did not progress much from the first proposed consideration. For many years, the simple estimate provided by \alpha \lambda-formalism has been intensively used by population synthesis studies and, not surprisingly, frequently contradicted observations. In recent years, the advancements in our studies of stellar structure, progress of the numerical techniques for hydrodynamical simulations as well as increase of the computer power and new observations started to bring improvements to our understanding of the common envelope phase. We review main physical processes taking place during the common envelope phase from the theoretical point of view and how they affect the values of classical formal parameters. In particular, we discuss the energy budget problem -- what are the energy sources, sinks and what is the condition for the envelope to disperse, as well as the importance of choosing the definition of the remnant core to the common envelope outcome.

One class of compact binaries of special interest is that of double white dwarfs (DWDs). For many of these systems, the exact nature of the evolutionary channels by which they form remains uncertain. The canonical explanation calls for the progenitor binary system to undergo two subsequent mass-transfer events, both of which are unstable and lead to a common envelope (CE) phase. However, it has been shown that if both CE events obey the standard alpha prescription, it is not possible to reproduce all of the observed systems. As an alternative prescription, the gamma-formalism was proposed, which parametrizes the fraction of angular momentum carried away with mass loss, in contrast to the alpha prescription, which parametrizes energy loss. We demonstrate that the gamma-prescription is also inadequate in describing the evolution of an arbitrary DWD binary; clearly we require a deeper understanding of the physical mechanisms underlying their formation. We then present a detailed model for the evolution of Red Giant -- Main Sequence binaries during the first episode of mass transfer, and demonstrate that their evolution into DWDs need not arise through two phases of dynamical mass loss. Instead, the first episode of dramatic mass loss may be stable, non-conservative mass transfer. The second phase is then well described by the alpha prescription. We find that the considered progenitors can reproduce the properties of the observed helium DWDs in which the younger component is the more massive.

In this Letter we discuss a modification to the criterion for the common envelope (CE) event to result in envelope dispersion. We emphasize that the current energy criterion for the CE phase is not sufficient for an instability of the CE, nor for an ejection. However, in some cases, stellar envelopes undergo stationary mass outflows, which are likely to occur during the slow spiral-in stage of the CE event. We propose the condition for such outflows, in a manner similar to the currently standard $\alpha_{\rm CE}\lambda$-prescription but with an addition of $P/\rho$ term in the energy balance equation, accounting therefore for the enthalpy of the envelope rather than merely the gas internal energy. This produces a significant correction, which might help to dispense with an unphysically high value of energy efficiency parameter during CE phase, currently required in the binary population synthesis studies to make the production of low-mass X-ray binaries (LMXBs) with a black hole companion to match the observations.

Although several dozen double white dwarfs (DWDs) have been observed, for many the exact nature of the evolutionary channel(s) by which they form remains uncertain. The canonical explanation calls for the progenitor binary system to undergo two subsequent mass-transfer events, both of which are unstable and lead to a common envelope (CE). However, it has been shown that if both CE events obey the standard alpha-prescription (parametrizing energy loss), it is not possible to reproduce all of the observed systems. The gamma-prescription was proposed as an alternative to this description, instead parametrizing the fraction of angular momentum carried away in dynamical-timescale mass loss. In this paper, we analyze simultaneous energy and angular-momentum conservation, and show that the gamma-prescription cannot adequately describe a CE event for an arbitrary binary, nor can the first phase of mass loss always be understood in general as a dynamical-timescale event. We consider in detail the first episode of mass transfer in binary systems with initially low companion masses, with a primary mass in the range 1.0--1.3 solar masses and an initial mass ratio between the secondary and primary stars of 0.83--0.92. In these systems, the first episode of dramatic mass loss may be stable, non-conservative mass transfer. This strips the donor's envelope and dramatically raises the mass ratio; the considered progenitor binary systems can then evolve into DWDs after passing through a single CE during the second episode of mass loss. We find that such a mechanism reproduces the properties of the observed DWD systems which have an older component with less than approximately 0.46 solar masses and mass ratios between the younger and older WDs greater than 1.

In contrast to the field, the binaries in dense stellar systems are frequently not primordial, and could be either dynamically formed or significantly altered from their primordial states. Destruction and formation of binaries occur in parallel all the time. The destruction, which constantly removes soft binaries from a binary pool, works as an energy sink and could be a reason for cluster entering the binary-burning phase. The true binary fraction is greater than observed, as a result, the observable binary fraction evolves differently from the predictions. Combined measurements of binary fractions in globular clusters suggest that most of the clusters are still core-contracting. The formation, on other hand, affects most the more evolutionary advanced stars, which significantly enhances the population of X-ray sources in globular clusters. The formation of binaries with a compact objects proceeds mainly through physical collisions, binary-binary and single-binary encounters; however, it is the dynamical formation of triples and multiple encounters that principally determine whether the formed binary will become an X-ray source.

One of the most important and uncertain stages in the binary evolution is the common envelope (CE) event. Significant attention has been devoted in the literature so far to the energy balance during the CE event, expected to determine the outcome. However this question is intrinsically coupled with the problem of what is left from the donor star after the CE and its immediate evolution. In this paper we argue that an important stage has been overlooked: post-CE remnant thermal readjustment phase. We propose a methodology for unambiguously defining the post-CE remnant mass after it has been thermally readjusted, namely by calling the core boundary the radius in the hydrogen shell corresponding to the local maximum of the sonic velocity. We argue that the important consequences of the thermal readjustment phase are: (i) a change in the energy budget requirement for the CE binaries and (ii) a companion spin-up and chemical enrichment, as a result of the mass transfer that occurs during the remnant thermal readjustment (TR). More CE binaries are expected to merge. If the companion is a neutron star, it will be mildly recycled during the TR phase. The mass transfer during the TR phase is much stronger than the accretion rate during the common envelope, and therefore satisfies the condition for a hypercritical accretion better. We also argue that the TR phase is responsible for a production of mildly recycled pulsars in double neutron stars.

Although many double white dwarfs (DWDs) have been observed, the evolutionary channel by which they are formed from low-mass/long-period red-giant-main-sequence (RG-MS) binaries remains uncertain. The canonical explanations involve some variant of double common-envelope (CE) evolution, however it has been found that such a mechanism cannot produce the observed distribution. We present a model for the initial episode of mass transfer (MT) in RG-MS binaries, and demonstrate that their evolution into double white dwarfs need not arise through a double-CE process, as long as the initial primary's core mass (Md,c) does not exceed 0.46M$_{\odot}$. Instead, the first episode of dramatic mass loss may be stable, non-conservative MT. We find a lower bound on the fraction of transferred mass that must be lost from the system in order to provide for MT, and demonstrate the feasibility of this channel in producing observed low-mass (with M$_{d,c}$ < 0.46M$_{\odot}$) DWD systems.

We present a new mechanism for the ejection of a common envelope in a massive binary, where the energy source is nuclear energy rather than orbital energy. This can occur during the slow merger of a massive primary with a secondary of 1-3 Msun when the primary has already completed helium core burning. We show that, in the final merging phase, hydrogen-rich material from the secondary can be injected into the helium-burning shell of the primary. This leads to a nuclear runaway and the explosive ejection of both the hydrogen and the helium layer, producing a close binary containing a CO star and a low-mass companion. We argue that this presents a viable scenario to produce short-period black-hole binaries and long-duration gamma-ray bursts (LGRBs). We estimate a LGRB rate of about 1.e-6 per year at solar metallicity, which implies that this may account for a significant fraction of all LGRBs, and that this rate should be higher at lower metallicity.

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.

Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to increase with time over a range of initial cluster central densities for initial binary fractions <~ 90%. The dominant processes driving the evolution of the core binary fraction are mass segregation of binaries into the cluster core and preferential destruction of binaries there. On a global scale, these effects and the preferential tidal stripping of single stars tend to roughly balance, leading to overall cluster binary fractions that are roughly constant with time. Our findings suggest that the current hard binary fraction near the half-mass radius is a good indicator of the hard primordial binary fraction. However, the relationship between the true binary fraction and the fraction of main-sequence stars in binaries (which is typically what observers measure) is non-linear and rather complicated. We also consider the importance of soft binaries, which not only modify the evolution of the binary fraction, but can drastically change the evolution of the cluster as a whole. Finally, we describe in some detail the recent addition of single and binary stellar evolution to our cluster evolution code.

In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor effective temperature, these kinematical constraints harbor a wealth of information on the system's past evolution. Here, we consider all this available information to reconstruct the full evolutionary history of the black hole X-ray binary XTE J1118+480, assuming that the system originated in the Galactic disk and the donor has solar metallicity. This analysis accounts for four evolutionary phases: mass transfer through the ongoing X-ray phase, tidal evolution before the onset of Roche-lobe overflow, motion through the Galactic potential after the formation of the black hole, and binary orbital dynamics due to explosive mass loss and possibly a black hole natal kick at the time of core collapse. We find that right after black hole formation, the system consists of a ~6.0-10.0 solar masses black hole and a ~1.0-1.6 solar masses main-sequence star. We also find that that an asymmetric natal kick is not only plausible but required for the formation of this system, and derive a lower and upper limit on the black hole natal kick velocity magnitude of 80 km/s and 310 km/s, respectively.

In this paper, the second of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing neutron stars (NSs) in dense globular clusters (GCs). For this study, 70 dense clusters were simulated independently, with a total stellar mass ~2x10^7Msun, exceeding the total mass of all dense GCs in our Galaxy. We find that, in order to reproduce the empirically derived formation rate of low-mass X-ray binaries (LMXBs), we must assume that NSs can be formed via electron-capture supernovae (ECS) with typical natal kicks smaller than in core-collapse supernovae. Our results explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different GCs. We predict that the number of quiescent LMXBs in different GCs should not have a strong metallicity dependence. In our cluster model the following mass-gaining events create populations of MSPs that do not match the observations: (i) accretion during a common envelope event with a NS formed through ECS, and (ii) mass transfer (MT) from a WD donor. Some processes lead only to a mild recycling. In addition, for MSPs, we distinguish low-magnetic-field (long-lived) and high-magnetic-field (short-lived) populations. With this distinction and by considering only those mass-gaining events that appear to lead to NS recycling, we obtain good agreement of our models with the numbers and characteristics of observed MSPs in 47 Tuc and Terzan 5, as well as with the cumulative statistics for MSPs detected in GCs of different dynamical properties. We find that significant production of merging double NSs potentially detectable as short gamma-ray bursts occurs only in very dense, most likely core-collapsed GCs. (abridged)

We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentricity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could provide unique dynamical identifications of these systems in globular clusters.

Using a population synthesis approach, we compute the total merger rate in the local Universe for double neutron stars, double black holes, and black hole -- neutron star binaries. These compact binaries are the prime source candidates for gravitational-wave detection by LIGO and VIRGO. We account for mergers originating from field populations and from dense stellar clusters. For both populations we use the same treatment of stellar evolution. Our results indicate that the merger rates of double neutron stars and black hole -- neutron star binaries are strongly dominated by field populations, while merging black hole binaries are formed much more effectively in dense stellar clusters. The overall merger rate of double compact objects depends sensitively on the (largely unknown) initial mass fraction contained in dense clusters f_cl. For f_cl < 0.0001, the Advanced LIGO detection rate will be dominated by field populations of double neutron star mergers, with a small but significant number of detections ~20 yr^-1. However for a higher mass fraction in clusters, f_cl > 0.001, the detection rate will be dominated by numerous mergers of double black holes originating from dense clusters, and it will be considerably higher, ~25 - 300 yr^-1. In addition, we show that, once mergers of double black holes are detected, it is easy to differentiate between systems formed in the field and in dense clusters, since the chirp mass distributions are strikingly different. Finally, we point out that there may exist a population of merging black hole binaries in intergalactic space.

Dynamical interactions that occur between objects in dense stellar systems are particularly important for the question of formation of X-ray binaries. We present results of numerical simulations of 70 globular clusters with different dynamical properties and a total stellar mass of 2*10^7 Msun. We find that in order to retain enough neutron stars to match observations we must assume that NSs can be formed via electron-capture supernovae. Our simulations explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different globular clusters. For millisecond pulsars, we obtain good agreement between our models and the numbers and characteristics of observed pulsars in the clusters Terzan 5 and 47 Tuc

In this contribution we discuss how neutron stars are produced and retained in globular clusters, outlining the most important dynamical channels and evolutionary events that affect thepopulation of mass-transferring binaries with neutron stars and result in the formation of recycled pulsars. We confirm the importance of electron-capture supernovae in globular clusters as the major supplier of retained neutron stars.By comparing the observed millisecond pulsar population and the results obtained from simulations, we discuss several constraints on the evolution of mass-transferring systems.In particular, we find that in our cluster model the following mass-gaining events create populations of MSPs that do not match the observations (with respect to binary periods and companion masses or the number of produced systems) and therefore likely do not lead to NSs spun up to millisecond periods: (i) accretion during a common envelope event with a NS formed through accretion-induced collapse, and (ii) mass transfer from a WD donor. By restricting ourselves to the evolutionary and dynamical paths that most likely lead to neutron star recycling, we obtain good agreement between our models and the numbers and characteristics of observed millisecond pulsars in the clusters Terzan 5 and 47 Tuc.

We present a comprehensive description of the population synthesis code StarTrack. The original code has been significantly modified and updated. Special emphasis is placed here on processes leading to the formation and further evolution of compact objects (white dwarfs, neutron stars, and black holes). Both single and binary star populations are considered. The code now incorporates detailed calculations of all mass-transfer phases, a full implementation of orbital evolution due to tides, as well as the most recent estimates of magnetic braking. This updated version of StarTrack can be used for a wide variety of problems, with relevance to many current and planned observatories, e.g., studies of X-ray binaries (Chandra, XMM-Newton), gravitational radiation sources (LIGO, LISA), and gamma-ray burst progenitors (HETE-II, Swift). The code has already been used in studies of Galactic and extra-galactic X-ray binary populations, black holes in young star clusters, Type Ia supernova progenitors, and double compact object populations. Here we describe in detail the input physics, we present the code calibration and tests, and we outline our current studies in the context of X-ray binary populations.

In recent years, an increasing number of proper motions have been measured for Galactic X-ray binaries. When supplemented with accurate determinations of the component masses, orbital period, and donor luminosity and effective temperature, these kinematical constraints harbor a wealth of information on the systems' past evolution. The constraints on compact object progenitors and kicks derived from this are of immense value for understanding compact object formation and exposing common threads and fundamental differences between black hole and neutron star formation. Here, we present the results of such an analysis for the black hole X-ray binary XTE J1118+480. We present results from modeling the mass transfer phase, following the motion in the Galaxy back to the birth site of the black hole, and examining the dynamics of symmetric and asymmetric core-collapses of the black hole progenitor.

This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes.

In this Letter we examine the idea that a subset of short-period black-hole low-mass X-ray binaries could be powered by the mass transfer from pre-main sequence donors. As the star contracts towards the main sequence, the strong magnetic fields operate the magnetic braking which dissipates the orbital angular momentum, driving the binary to contact. We show that the periods and apparent donor spectral classes of the X-ray binaries with a pre-main sequence donor agree better with the available observations of black hole X-ray binaries than those of binaries with a main-sequence donor. This mechanism also explains, without the need for additional hypotheses, the roughly primordial abundance of Li detected in donor companions of black hole X-ray candidates in our Galaxy.

We present a theoretical study of double compact objects as potential short/hard gamma-ray burst (GRB) progenitors. An updated population synthesis code StarTrack is used to calculate properties of double neutron stars and black-hole neutron star binaries. We obtain their formation rates, estimate merger times and finally predict their most likely merger locations and afterglow properties for different types of host galaxies. Our results serve for a direct comparison with the recent HETE-II and SWIFT observations of several short bursts, for which afterglows and host galaxies were detected. We also discuss the possible constraints these observations put on the evolutionary models of double compact object formation. We emphasize that our double compact object models can successfully reproduce at the same time short GRBs within both young, star-forming galaxies (e.g., GRB 050709 and GRB 051221A), as well as within old, elliptical hosts (e.g., GRB 050724 and probably GRB 050509B)

In this paper, the first of a series, we study the stellar dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We examine the processes leading to the creation of X-ray binaries such as cataclysmic variables and AM CVn systems. Using numerical simulations, we identify the dominant formation channels and we predict the expected numbers and characteristics of detectable systems, emphasizing how the cluster sources differ from the field population. We explore the dependence of formation rates on cluster properties and we explain in particular why the distribution of cataclysmic variables has only a weak dependence on cluster density. We also discuss the frequency of dwarf nova outbursts in globular clusters and their connection with moderately strong white dwarf magnetic fields. We examine the rate of Type Ia supernovae via both single and double degenerate channels in clusters and we argue that those rates may contribute to the total SN Ia rate in elliptical galaxies. Considering coalescing white dwarf binaries we discuss possible constraints on the common envelope evolution of their progenitors and we derive theoretical expectations for gravitational wave detection by LISA.

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

We model the dynamical evolution of primordial black holes (BHs) in dense star clusters using a simplified treatment of stellar dynamics in which the BHs are assumed to remain concentrated in an inner core, completely decoupled from the background stars. Dynamical interactions involving BH binaries are computed exactly and are generated according to a Monte Carlo prescription. Recoil and ejections lead to complete evaporation of the BH core on a timescale ~10^9 yr for typical globular cluster parameters. Orbital decay driven by gravitational radiation can make binaries merge and, in some cases, successive mergers can lead to significant BH growth. Our highly simplified treatment of the cluster dynamics allows us to study a large number of models and to compute statistical distributions of outcomes, such as the probability of massive BH growth and retention in a cluster. We find that, in most models, there is a significant probability (~20-80%) of BH growth with final masses > 100 M_{\sun}. In one case, a final BH formed with mass ~ 620 M_{\sun}. However, if the typical merger recoil speed (due to asymmetric emission of gravitational radiation) significantly exceeds the cluster escape speed, no growth ever occurs. Independent of the recoil speed, we find that BH-BH mergers enhanced by dynamical interactions in cluster cores present an important source of gravitational waves for ground--based laser interferometers. Under optimistic conditions, the total rate of detections by Advanced LIGO, could be as high as a few tens of events per year from inspiraling BHs from clusters.

We propose that the observed difference in the formation rates of bright low-mass X-ray binaries in metal-rich and metal-poor globular clusters can be explained by taking into account the difference in the stellar structure of main sequence donors with masses between ~0.85 Msun and ~1.25 Msun at different metallicities. This difference is caused by the absence of an outer convective zone in metal-poor main sequence stars in this mass range. In the result, magnetic braking, a powerful mechanism of orbital shrinkage, does not operate and dynamically formed main sequence - neutron star binaries fail to start mass transfer or appear as bright low-mass X-ray binaries.

We propose that the shape of the upper-end X-ray luminosity function observed in elliptical galaxies for point sources is a footprint of the black-hole (BH) mass spectrum among old X-ray transients formed in the galaxies. We show that this underlying BH mass spectrum is modified by a weighting factor that is related to the transient duty cycle and it generally depends on the BH mass and XRB donor type (main-sequence, red-giant, or white-dwarf donors). A duty cycle that depends on the binary mass-transfer rate relative to the critical rate for transient behavior is most probably favored. We also find that the derived BH mass spectrum slope depends on the strength of angular momentum loss due to magnetic braking for main-sequence donors. More specifically, we find that, for ``standard'' magnetic braking, BH XRBs with red-giant donors dominate the upper-end XLF; for weaker magnetic braking prescriptions main-sequence donors are found to be dominant. In both cases the BH mass spectrum has a differential slope of ~2.5 and an upper BH mass cut-off at ~20 Msun is needed to understand the very brightest of the BH XRBs in elliptical galaxies. We expect that our result will help to constrain binary population synthesis models and the adopted relations between black holes and the masses of their progenitors.

Observations of galactic and extra-galactic globular clusters have shown that on average metal-rich clusters are ~3 times as likely to contain a bright X-ray source than their metal-poor counterparts. We propose that this can be explained by taking into account the difference in the stellar structure of main sequence donors with masses between ~0.85 Msun and ~1.25 Msun at different metallicities. Metal-poor main sequence stars in this mass range do not have an outer convective zone while metal-rich stars do. The absence of this zone turns off magnetic braking, a powerful mechanism of orbital shrinkage, leading to the failure of dynamically formed main sequence - neutron star binaries to start mass transfer or appear as bright low-mass X-ray binaries.

We analyze the upper-end X-ray luminosity function (XLF) observed in elliptical galaxies for point sources. We propose that the observed XLF is dominated by transient BH systems in outburst and the XLF shape reflects the black hole (BH) mass spectrum among old X-ray transients. The BH mass spectrum -- XLF connection depends on a weighting factor that is related to the transient duty cycle and depends on the host-galaxy age, the BH mass and the donor type (main sequence, red giant, or white dwarf). We argue that the assumption of a constant duty cycle for all systems leads to results inconsistent with current observations. The type of dominant donors in the upper-end XLF depends on what type of magnetic braking operates: in the case of ``standard'' magnetic braking, BH X-ray binaries with red-giant donors dominate, and in the case of weaker magnetic braking prescriptions main sequence donors dominate.

In dense stellar systems, frequent dynamical interactions between binaries lead to the formation of multiple systems. In this contribution we discuss the dynamical formation of hierarchically stable triples: the formation rate, main characteristics of dynamically formed population of triples and the impact of the triples formation on the population of close binaries. In particular, we estimate how much the population of blue stragglers and compact binaries could be affected.

While many observed ultra-luminous X-ray sources (ULXs, Lx > 10^39 erg s^-1) could be extragalactic X-ray binaries (XRBs) emitting close to the Eddington limit, the highest-luminosity ULXs (Lx > 3x10^39 erg s^-1) exceed the isotropic Eddington luminosity for even high-stellar-mass accreting black hole XRBs. It has been suggested that these highest-luminosity ULXs may contain accreting intermediate-mass black hole (IMBH) binaries. We consider this hypothesis for dense, young (about 100 Myr) stellar clusters where we assume that a 50-500 solar mass central IMBH has formed through runaway growth of a massive star. Using numerical simulations of the dynamics and evolution of the central black hole's captured companions, we obtain estimates of the incidence of mass transfer phases and possible ULX activity throughout the IMBH's evolutionary history. We find that, although it is common for the central black hole to acquire binary companions, there is a very low probability that these interacting binaries will become observable ULX sources.

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.

In recent years proper motion measurements have been added to the set of observational constraints on the current properties of Galactic X-ray binaries. We develop an analysis that allows us to consider all this available information and reconstruct the full evolutionary history of X-ray binaries back to the time of compact object formation. This analysis accounts for mass transfer through the ongoing X-ray phase, tidal circularization before the onset of Roche-lobe overflow, motion through the Galactic potential after the formation of the compact object, and binary orbital dynamics and hydrodynamic modeling of the core collapse. We apply the analysis to the soft X-ray transient GRO J1655-40 and, for the first time, use its full 3D peculiar velocity constraints right after core collapse instead of lower limits on the current space velocity given by the present-day radial velocity. We find that the system has remained within 200pc from the Galactic plane throughout its entire life time and that the mass loss and a kick possibly associated with the black hole (BH) formation imparted a kick velocity of 45-115 km/s to the binary's center of mass. Right after BH formation, the system consists of a 3.5-6.3 Msun BH and a 2.3-4 Msun main-sequence star. At the onset of the X-ray phase the donor is still on the main sequence. We find that a symmetric BH formation event cannot be formally excluded, but that the associated system parameters are only marginally consistent with the currently observed binary properties. BH formation mechanisms involving an asymmetric supernova explosion with associated BH kick velocities of a few tens of km/s, on the other hand, satisfy the constraints much more comfortably. We also derive an upper limit on the BH kick magnitude of 210 km/s. (abridged)

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code (StarTrack), and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing 3-body and 4-body interactions (Fewbody). We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5-10%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.

The population of compact binaries in dense stellar systems is affected strongly by frequent dynamical interactions between stars and their interplay with the stellar evolution. In this contribution, we consider these effects on binaries with a white dwarf accretor, in particular cataclysmic variables and AM CVns. We examine which processes can successfully lead to the creation of such X-ray binaries. Using numerical simulations, we identify predominant formation channels and predict the expected numbers of detectable systems. We discuss also why the distribution of cataclysmic variables has a weaker dependence upon the cluster density than the distribution of quiescent low-mass X-ray binaries and why dwarf nova outbursts may not occur among globular cluster cataclysmic variables.

We investigate the dynamical formation and evolution of binaries containing neutron stars in dense globular clusters. Our numerical simulations combine a simple Monte Carlo prescription for stellar dynamics, a sophisticated binary population synthesis code, and a small-N-body integrator for computing 3-body and 4-body interactions. Our results suggest that there is no ``retention problem,'' i.e., that, under standard assumptions, globular clusters can retain enough neutron stars to produce the observed numbers of millisecond pulsars. We also identify the dominant evolutionary and dynamical channels through which globular clusters produce their two main types of binary millisecond pulsars

We consider observational tests for the nature of Ultraluminous X-ray sources (ULXs). These must distinguish between thermal-timescale mass transfer on to stellar-mass black holes leading to anisotropic X-ray emission, and accretion on to intermediate-mass black holes. We suggest that long-term transient behavior via the thermal-viscous disk instability could discriminate between these two possibilities for ULXs in regions of young stellar populations. Thermal-timescale mass transfer generally produces stable disks and persistent X-ray emission. In contrast, mass transfer from massive stars to black holes produces unstable disks and thus transient behavior, provided that the black hole mass exceeds some minimum value. This minimum mass depends primarily on the donor mass and evolutionary state. We show that it exceeds 50 solar masses for a large fraction (greater than 90%) of the mass-transfer lifetime for the most likely donors in young clusters. Thus if long-term monitoring reveals a large transient fraction among ULXs in a young stellar population, these systems would be good candidates for intermediate-mass black holes in a statistical sense; information about the donor star is needed to make this identification secure in any individual case. A transient ULX population would imply a much larger population of quiescent systems of the same type.

In dense stellar systems the frequent dynamical interactions between stars play a crucial role in the formation and evolution of compact binaries. We study these processes using a novel approach combining a state-of-the-art binary population synthesis code with a simple treatment of dynamical interactions in dense star cluster cores. Here we focus on the dynamical and evolutionary processes leading to the formation of compact binaries containing white dwarfs in dense globular clusters. We demonstrate that dynamics can increase by factors ~ 2 - 100 the production rates of interesting binaries such as cataclysmic variables, ``nonflickerers'' (He white dwarfs with a heavier dark companion), merging white dwarf binaries with total masses above the Chandrasekhar limit, and white dwarf binaries emitting gravitational waves in the LISA band.

We study the evolution of binary stars in globular clusters using a novel approach combining a state-of-the-art population synthesis code with a simple treatment of dynamical interactions in the dense cluster core. We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.

The evolution of binaries consisting of evolved main sequence stars (1 < M_d/Msun < 3.5) with white dwarf companions (0.7 < M_wd/Msun < 1.2) is investigated through the thermal mass transfer phase. Taking into account the stabilizing effect of a strong, optically thick wind from the accreting white dwarf surface, we have explored the formation of several evolutionary groups of systems for progenitors with initial orbital periods of 1 and 2 days. The numerical results show that CO white dwarfs can accrete sufficient mass to evolve to a Type Ia supernova and ONeMg white dwarfs can be built up to undergo accretion induced collapse for donors more massive than about 2 Msun. For donors less massive than ~2 Msun the system can evolve to form a He and CO or ONeMg white dwarf pair. In addition, sufficient helium can be accumulated (~0.1 Msun) in systems characterized by 1.6 < M_d/Msun < 1.9 and 0.8 < M_wd/Msun < 1 such that sub Chandrasekhar mass models for Type Ia supernovae, involving off center helium ignition, are possible for progenitor systems evolving via the Case A mass transfer phase. For systems characterized by mass ratios > 3 the system likely merges as a result of the occurrence of a delayed dynamical mass transfer instability. A semi-analytical model is developed to delineate these phases which can be easily incorporated in population synthesis studies of these systems.

We present a description for the angular momentum loss rate due to magnetic braking for late type stars taking into account recent observational data on the relationship between stellar activity and rotation. The analysis is based on an idealized two component coronal model subject to constraints imposed on the variation of the coronal gas density with rotation period inferred from the observed variation of X-ray luminosity, L_x, with rotation rate, Omega, (L_x \propto Omega^2) for single rotating dwarfs. An application of the model to high rotation rates leads to a gradual turnover of the X-ray luminosity which is similar to the saturation recently observed in rapidly rotating dwarfs. The resulting angular momentum loss rate, \dot J, depends on Omega in the form \dot J \propto Omega^beta where beta ~ 3 for slow rotators and ~ 1.3 for fast rotators. The relation at high rotation rates significantly differs from the power law exponent for slowly rotating stars, depressing the angular momentum loss rate without necessarily requiring the saturation of the magnetic field. The application of this model to the evolution of cataclysmic variable binary systems leads to mass transfer rates that are in approximate accord with those observed in comparison to rates based on either a Skumanich law or an empirical law based on beta=1.

We have calculated the evolution of 60 model binary systems consisting of helium stars in the mass range of M_He= 2.5-6Msun with a 1.4Msun neutron star companion to investigate the formation of double neutron star systems.Orbital periods ranging from 0.09 to 2 days are considered, corresponding to Roche lobe overflow starting from the helium main sequence to after the ignition of carbon burning in the core. We have also examined the evolution into a common envelope phase via secular instability, delayed dynamical instability, and the consequence of matter filling the neutron star's Roche lobe. The survival of some close He-star neutron-star binaries through the last mass transfer episode (either dynamically stable or unstable mass transfer phase) leads to the formation of extremely short-period double neutron star systems (with P<~0.1days). In addition, we find that systems throughout the entire calculated mass range can evolve into a common envelope phase, depending on the orbital period at the onset of mass transfer. The critical orbital period below which common envelope evolution occurs generally increases with M_He. In addition, a common envelope phase may occur during a short time for systems characterized by orbital periods of 0.1-0.5 days at low He-star masses (~ 2.6-3.3Msun). The existence of a short-period population of double neutron stars increases the predicted detection rate of inspiral events by ground-based gravitational-wave detectors and impacts their merger location in host galaxies and their possible role as gamma-ray burst progenitors. We use a set of population synthesis calculations and investigate the implications of the mass-transfer results for the orbital properties of DNS populations.

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 study the complete merger of two massive stars inside a common envelope and the subsequent evolution of the merger product, a rapidly rotating massive supergiant. Three qualitatively different types of mergers have been identified and investigated in detail, and the post-merger evolution has been followed to the immediate presupernova stage. The ``quiet merger'' case does not lead to significant changes in composition, and the star remains a red supergiant. In the case of a ``moderate merger'', the star may become a blue supergiant and end its evolution as a blue supergiant, depending on the core to total mass ratio (as may be appropriate for the progenitor of SN 1987A). In the case of the most effective ``explosive merger'', the merger product stays a red giant. In last two cases, the He abundance in the envelope is increased drastically, but significant s-processing is mainly expected in the ``explosive merger'' case.

Subdwarf B (sdB) stars (and related sdO/sdOB stars) are believed to be helium core-burning objects with very thin hydrogen-rich envelopes. In recent years it has become increasingly clear from observational surveys that a large fraction of these objects are members of binary systems. To better understand their formation, we here present the results of a detailed investigation of the three main binary evolution channels that can lead to the formation of sdB stars: the common envelope (CE) ejection channel, the stable Roche lobe overflow (RLOF) channel and the double helium white dwarfs (WDs) merger channel. We obtained the conditions for the formation of sdB stars from each of these channels using detailed stellar and binary evolution calculations where we modelled the detailed evolution of sdB stars and carried out simplified binary population synthesis simulations. The observed period distribution of sdB stars in compact binaries strongly constrains the CE ejection parameters. We also present the distribution of sdB stars in the $T_{\rm eff}$ - $\log g$ diagram, the Hertzsprung-Russell diagram and the distribution of mass functions.

We study the nucleosynthesis and the induced mixing during the merging of massive stars inside a common envelope. The systems of interest are close binaries, initially consisting of a massive red supergiant and a main-sequence companion of a few solar masses. We apply parameterized results based on hydrodynamical simulations to model the stream-core interaction and the response of the star in a standard stellar-evolution code. Preliminary results are presented illustrating the possibility of unusual nucleosynthesis and post-merging dredge-up which can cause composition anomalies in the supergiant's envelope.

We present detailed simulations of the interaction of a stream emanating from a mass-losing secondary with the core of a massive supergiant in the slow merger of the two stars inside a common envelope. The dynamics of the stream can be divided into a ballistic phase, starting at the L_1 point, and a hydrodynamical phase where the stream interacts strongly with the core. Considering the merger of a 1 and 5Msun star with a 20Msun evolved supergiant, we present two-dimensional hydrodynamical simulations using the PROMETHEUS code to demonstrate how the penetration depth and post-impact conditions depend on the initial properties of stream material (e.g. entropy, angular momentum, stream width) and the properties of the core (e.g. density structure and rotation rate). Using these results, we present a fitting formula for the entropy generated in the stream--core interaction and a recipe for the determination of the penetration depth based on a modified Bernoulli integral.

We study the merging of massive stars inside a common envelope for binary systems consisting of a red supergiant with a mass of 15-20 Msun and a main-sequence companion of 1-5 Msun. We are particularly interested in the stage when the secondary, having overfilled its Roche lobe inside the common envelope, starts to transfer mass to the core of the primary at a very high mass-transfer rate and the subsequent nucleo-synthesis in the core-impact region. Using a parametrized model for the structure of the envelope at this stage, we perform 2-dimensional hydrodynamical calculations with the Munich Prometheus code to calculate the dynamics of the stream emanating from the secondary and its impact on the core of the primary. We find that, for the lower end of the estimated mass-transfer rate, low-entropy, hydrogen-rich material can penetrate deep into the primary core where nucleosynthesis through the hot CNO cycle can take place and that the associated neutron exposure may be sufficiently high for significant s-processing. For mass-transfer rates at the high end of our estimated range and higher densities in the stream, the stream impact can lead to the dredge-up of helium, but the neutron production is too low for significant s-processing.