Many exoplanets have now been detected in orbits with ultra-short periods, very close to the Roche limit. Building upon our previous work, we study the possibility that mass loss through Roche lobe overflow (RLO) may affect the evolution of these planets, and could possibly transform a hot Jupiter into a lower-mass planet (hot Neptune or super-Earth). We focus here on systems in which the mass loss occurs slowly ("stable mass transfer" in the language of binary star evolution) and we compute their evolution in detail with the binary evolution code MESA. We include the effects of tides, RLO, irradiation and photo-evaporation of the planet, as well as the stellar wind and magnetic braking. Our calculations all start with a hot Jupiter close to its Roche limit, in orbit around a sun-like star. The initial orbital decay and onset of RLO are driven by tidal dissipation in the star. We confirm that such a system can indeed evolve to produce lower-mass planets in orbits of a few days. The RLO phase eventually ends and, depending on the details of the mass transfer and on the planetary core mass, the orbital period can remain around a few days for several Gyr. The remnant planets have a rocky core and some amount of envelope material, which is slowly removed via photo-evaporation at nearly constant orbital period; these have properties resembling many of the observed super-Earths and sub-Neptunes. For these remnant planets we also predict an anti-correlation between mass and orbital period; very low-mass planets ($M_{\rm pl}\,\lesssim\,5\,M_{\oplus}$) in ultra-short periods ($P_{\rm orb}$<1d) cannot be produced through this type of evolution.

Although not nearly as numerous as binaries with two white dwarfs, eccentric neutron star-white dwarf (NS-WD) binaries are important gravitational-wave (GW) sources for the next generation of space-based detectors sensitive to low frequency waves. Here we investigate periastron precession in these sources as a result of general relativistic, tidal, and rotational effects; such precession is expected to be detectable for at least some of the detected binaries of this type. Currently, two eccentric NS-WD binaries are known in the galactic field, PSR J1141-6545 and PSR B2303+46, both of which have orbits too wide to be relevant in their current state to GW observations. However, population synthesis studies predict the existence of a significant Galactic population of such systems. Though small in most of these systems, we find that tidally induced periastron precession becomes important when tides contribute to more than 3% of the total precession rate. For these systems, accounting for tides when analyzing periastron precession rate measurements can improve estimates of the WD component mass inferred and, in some cases, will prevent us from misclassifying the object. However, such systems are rare due to rapid orbital decay. To aid the inclusion of tidal effects when using periastron precession as a mass measurement tool, we derive a function that relates the WD radius and periastron precession constant to the WD mass.

Kepler-56 is a multi-planet system containing two coplanar inner planets that are in orbits misaligned with respect to the spin axis of the host star, and an outer planet. Various mechanisms have been proposed to explain the broad distribution of spin-orbit angles among exoplanets, and these theories fall under two broad categories. The first is based on dynamical interactions in a multi-body system, while the other assumes that disk migration is the driving mechanism in planetary configuration and that the star (or disk) is titled with respect to the planetary plane. Here we show that the large observed obliquity of Kepler-56 system is consistent with a dynamical origin. In addition, we use observations by Huber et al. (2013) to derive the obliquity's probability distribution function, thus improving the constrained lower limit. The outer planet may be the cause of the inner planets' large obliquities, and we give the probability distribution function of its inclination, which depends on the initial orbital configuration of the planetary system. We show that even in the presence of precise measurement of the true obliquity, one cannot distinguish the initial configurations. Finally we consider the fate of the system as the star continues to evolve beyond the main sequence, and we find that the obliquity of the system will not undergo major variations as the star climbs the red giant branch. We follow the evolution of the system and find that the innermost planet will be engulfed in ~129 Myr. Furthermore we put an upper limit of ~155 Myr for the engulfment of the second planet. This corresponds to ~ 3% of the current age of the star.

Through tidal dissipation in a slowly spinning host star the orbits of many hot Jupiters may decay down to the Roche limit. We expect that in most cases the ensuing mass transfer will be stable. Using detailed numerical calculations we find that this evolution is quite rapid, potentially leading to complete removal of the gaseous envelope in a few Gyr, and leaving behind an exposed rocky core ("hot super-Earth"). Final orbital periods are quite sensitive to the details of the planet's mass-radius relation, and to the effects of irradiation and photo-evaporation, but could be as short as a few hours, or as long as several days. Our scenario predicts the existence of planets with intermediate masses ("hot Neptunes") that should be found precisely at their Roche limit and in the process of losing mass through Roche lobe overflow. The observed excess of small single-planet candidate systems observed by Kepler may also be the result of this process. If so, the properties of their host stars should track those of the hot Jupiters. Moreover, the number of systems that produced hot Jupiters could be 2-3 times larger than one would infer from contemporary observations.

Hot Jupiters formed through circularization of high-eccentricity orbits should be found at orbital separations $a$ exceeding $twice$ that of their Roche limit $a_{\rm R}$. Nevertheless, about a dozen giant planets have now been found well within this limit ($a_{\rm R}<a<2a_{\rm R}$), with one coming as close as 1.2$a_{\rm R}$. In this Letter, we show that orbital decay (starting beyond 2$a_{\rm R}$) driven by tidal dissipation in the star can naturally explain these objects. For a few systems (WASP-4 and 19), this explanation requires the linear reduction in convective tidal dissipation proposed originally by Zahn (1966, 1989) and verified by recent numerical simulations Penev et al. (2007), but rules out the quadratic prescription proposed by Goldreich & Nicholson (1977). Additionally, we find that WASP-19-like systems could potentially provide direct empirical constraints on tidal dissipation, as we could soon be able to measure their orbital decay through high precision transit timing measurements.

The extragalactic X-ray binary IC 10 X-1 has attracted attention as it is possibly the host of the most massive stellar-mass black-hole (BH) known to date. Here we consider all available observational constraints and construct its evolutionary history up to the instant just before the formation of the BH. Our analysis accounts for the simplest possible history that includes three evolutionary phases: binary orbital dynamics at core collapse, common envelope (CE) evolution, and evolution of the BH--helium star binary progenitor of the observed system. We derive the complete set of constraints on the progenitor system at various evolutionary stages. Specifically: right before the core collapse event, we find the mass of the BH immediate progenitor to be > 31 Msun (at 95% of confidence, same hereafter). The magnitude of the natal kick imparted to the BH is constrained to be < 130 km/s. Furthermore, we find that the "enthalpy" formalism recently suggested by Ivanova & Chaichenets is able to explain the existence of IC 10 X-1 without the need of invoking unreasonably high CE efficiencies. With this physically motivated formalism, we find that the CE efficiency required to explain the system is in the range of 0.6--1.

Two formation scenarios have been proposed to explain the tight orbits of hot Jupiters. They could be formed in orbits with a small inclination (with respect to the stellar spin) via disk migration, or in more highly inclined orbits via high-eccentricity migration, where gravitational interactions with a companion and tidal dissipation are at play. Here we target hot Jupiter systems where the misalignment $\lambda$ has been inferred observationally and we investigate whether their properties are consistent with high-eccentricity migration. Specifically, we study whether stellar tides can be responsible for the observed distribution of $\lambda$ and orbital separations. Improving on previous studies, we use detailed models for each star, thus accounting for how convection (and tidal dissipation) depends on stellar properties. In line with observations suggesting that hotter stars have higher $\lambda$, we find that $\lambda$ increases as the amount of stellar surface convection decreases. This trend supports the hypothesis that tides are the mechanism shaping the observed distribution of $\lambda$. Furthermore, we study the past orbital evolution of five representative systems, chosen to cover a variety of temperatures and misalignments. We consider various initial orbital configurations and integrate the equations describing the coupled evolution of the orbital separation, stellar spin, and misalignment. We account for stellar tides and wind mass loss, stellar evolution, and magnetic braking. We show that the current properties of these five representative systems can be explained naturally, given our current understanding of tidal dissipation and with physically motivated assumptions for the effects driving the orbital evolution.

Here we present CAFein, a new computational tool for investigating radiative dissipation of dynamic tides in close binaries and of non-adiabatic, non-radial stellar oscillations in isolated stars in the linear regime. For the latter, CAFein computes the non-adiabatic eigenfrequencies and eigenfunctions of detailed stellar models. The code is based on the so-called Riccati method, a numerical algorithm that has been successfully applied to a variety of stellar pulsators, and which doesn't suffer of the major drawbacks of commonly used shooting and relaxation schemes. Here we present an extension of the Riccati method to investigate dynamic tides in close binaries. We demonstrate CAFein's capabilities as a stellar pulsation code both in the adiabatic and non-adiabatic regime, by reproducing previously published eigenfrequencies of a polytrope, and by successfully identifying the unstable modes of a stellar model in the $\beta$ Cephei/SPB region of the Hertzsprung-Russell diagram. Finally, we verify CAFein's behavior in the dynamic tides regime by investigating the effects of dynamic tides on the eigenfunctions and orbital and spin evolution of massive Main Sequence stars in eccentric binaries, and of hot Jupiter host stars. The plethora of asteroseismic data provided by the NASA's Kepler satellite, some of which include the direct detection of tidally excited stellar oscillations, make CAFein quite timely. Furthermore, the increasing number of observed short-period detached double white dwarfs (WD) and the observed orbital decay in the tightest of such binaries open up a new possibility of investigating WD interiors through the effects of tides on their orbital evolution

In recent years, the Extremely Low Mass White Dwarf (ELM WD) survey has quintupled the number of known close, detached double WD binaries (DWD). The tightest such DWD, SDSS J065133.33+284423.3 (J0651), harbors a He WD eclipsing a C/O WD every $\simeq\,12\,$min. The orbital decay of this source was recently measured to be consistent with general relativistic (GR) radiation. Here we investigate the role of dynamic tides in a J0651-Like binary and we uncover the potentially new phenomenon of "antiresonance" locking. In the most probable scenario of an asynchronous binary at birth, we find that dynamic tides play a key role in explaining the measured GR-driven orbital decay, as they lock the system at stable antiresonances with the star's eigenfrequencies. We show how such locking is naturally achieved and how, while locked at an antiresonance, GR drives the evolution of the orbital separation, while dynamic tides act to synchronize the spin of the He WD with the companion's orbital motion, but \emph{only on the GR timescale}. Given the relevant orbital and spin evolution timescales, the system is clearly on its way to synchronism, if not already synchronized.

The paucity of known High-Mass X-Ray Binaries (HMXB) with naked He donor stars (hereafter He star) in the Galaxy has been noted over the years as a surprising fact, given the significant number of Galactic HMXBs containing H-rich donors, which are expected to be their progenitors. This contrast has further sharpened in light of recent observations uncovering a preponderance of HMXBs hosting loosely bound Be donors orbiting neutron stars (NS), which would be expected to naturally evolve into He-HMXBs through dynamical mass transfer onto the NS and a common-envelope (CE) phase. Hence, reconciling the large population of Be-HMXBs with the observation of only one He-HMXB can help constrain the dynamics of CE physics. Here, we use detailed stellar structure and evolution models and show that binary mergers of HMXBs during CE events must be common in order to resolve the tension between these observed populations. We find that, quantitatively, this scenario remains consistent with the typically adopted energy parameterization of CE evolution, yielding expected populations which are not at odds with current observations. However, future observations which better constrain the underlying population of loosely bound O/B-NS binaries are likely to place significant constraints on the efficiency of CE ejection.

In recent years, accurate observational constraints become available for an increasing number of Galactic X-ray binaries. Together with proper motion measurements, we could reconstruct the full evolutionary history of X-ray binaries back to the time of compact object formation. In this paper, we present the first study of the persistent X-ray source Cygnus\;X-1 that takes into account of all available observational constraints. Our analysis accounts for three evolutionary phases: orbital evolution and motion through the Galactic potential after the formation of black hole (BH), and binary orbital dynamics at the time of core collapse. We find that the mass of the BH immediate progenitor is $15.0 - 20.0$ M$_\sun$, and at the time of core collapse, the BH has potentially received a small kick velocity of $\le 77$ km s$^{-1}$ at 95% confidence. If the BH progenitor mass is less than $\sim 17$ M$_\sun$, a non zero natal kick velocity is required to explain the currently observed properties of Cygnus\;X-1. Since the BH has only accreted mass from its companion's stellar wind, the negligible amount of accreted mass is impossible to explain the observationally inferred BH spin of $a_* > 0.95$, and the origin of this extreme BH spin must be connected to the BH formation itself. Right after the BH formation, we find that the BH companion is a $19.8 - 22.6$ M$_\sun$ main sequence star, orbiting the BH at a period of $4.7 - 5.2$ days. Furthermore, recent observations show that the BH companion is currently super-synchronized. This super-synchronism indicates that the strength of tides exerted on the BH companion should be weaker by a factor of at least two compared to the usually adopted strength.

Galactic interacting double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the GW detector LISA, with more than 10^4 binaries expected to be detected over the mission's lifetime. Part of this population is expected to be eccentric, and here we investigate the potential for constraining the white dwarf (WD) properties through apsidal precession in these binaries. We analyze the tidal, rotational, and general relativistic contributions to apsidal precession by using detailed He WD models, where the evolution of the star's interior is followed throughout the cooling phase. In agreement with previous studies of zero-temperature WDs, we find that apsidal precession in eccentric DWDs can lead to a detectable shift in the emitted GW signal when binaries with cool (old) components are considered. This shift increases significantly for hot (young) WDs. We find that apsidal motion in hot (cool) DWDs is dominated by tides at orbital frequencies above ~10^{-4}Hz (10^{- 3}$Hz). The analysis of apsidal precession in these sources while ignoring the tidal component would lead to an extreme bias in the mass determination, and could lead us to misidentify WDs as neutron stars or black holes. We use the detailed WD models to show that for older, cold WDs, there is a unique relationship that ties the radius and apsidal precession constant to the WD masses, therefore allowing tides to be used as a tool to constrain the source masses.

Galactic short period double white dwarfs (DWD) are guaranteed gravitational wave (GW) sources for the next generation of space-based interferometers sensitive to low-frequency GWs (10^{-4}- 1 Hz). Here we investigate the possibility of constraining the white dwarf (WD) properties through measurements of apsidal precession in eccentric binaries. We analyze the general relativistic (GR), tidal, and rotational contributions to apsidal precession by using detailed He WD models. We find that apsidal precession can lead to a detectable shift in the emitted GW signal, the effect being stronger (weaker) for binaries hosting hot (cool) WDs. We find that in hot (cool) DWDs tides dominate the precession at orbital frequencies above ~0.01 mHz (~1 mHz). Analyzing the apsidal precession of these sources only accounting for GR would potentially lead to an extreme overestimate of the component masses. Finally, we derive a relation that ties the radius and apsidal precession constant of cool WD components to their masses, therefore allowing tides to be used as an additional mass measurement tool.

M33 X-7 is among the most massive X-Ray binary stellar systems known, hosting a rapidly spinning 15.65 Msun black hole orbiting an underluminous 70 Msun Main Sequence companion in a slightly eccentric 3.45 day orbit. Although post-main-sequence mass transfer explains the masses and tight orbit, it leaves unexplained the observed X-Ray luminosity, star's underluminosity, black hole's spin, and eccentricity. A common envelope phase, or rotational mixing, could explain the orbit, but the former would lead to a merger and the latter to an overluminous companion. A merger would also ensue if mass transfer to the black hole were invoked for its spin-up. Here we report that, if M33 X-7 started as a primary of 85-99 Msun and a secondary of 28-32 Msun, in a 2.8-3.1 day orbit, its observed properties can be consistently explained. In this model, the Main Sequence primary transferred part of its envelope to the secondary and lost the rest in a wind; it ended its life as a ~16 Msun He star with a Fe-Ni core which collapsed to a black hole (with or without an accompanying supernova). The release of binding energy and, possibly, collapse asymmetries "kicked" the nascent black hole into an eccentric orbit. Wind accretion explains the X-Ray luminosity, while the black hole spin can be natal.

Black hole (BH) X-ray binaries (XRBs) are X-ray luminous binary systems comprising a BH accreting matter from a companion star. Understanding their origins sheds light on the still not well understood physics of BH formation. M33 X-7 hosts one of the most massive stellar-mass BH among all XRBs known to date, a 15.65 Msun BH orbiting a 70 Msun companion star in a 3.45 day orbit. The high masses of the two components and the tight orbit relative to the large H-rich stellar component challenge our understanding of the typically invoked BH-XRBs formation channels. The measured underluminosity of the optical component further complicates the picture. A solution to the evolutionary history of this system that can account for all its observed properties has yet to be presented, and here we propose the first scenario that is consistent with the complete set of current observational constraints. In our model, M33 X-7 started its life hosting a 85-99 Msun primary and a 28-32 Msun companion in a Keplerian orbit of 2.8-3.1 days. In order to form a BH of 15.65 Msun, the initially most massive component transferred part of its envelope to the companion star and lost the rest in a strong stellar wind. During this dynamically stable mass transfer phase the companion accreted matter, to become the presently underluminous 70 Msun star.

We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.

We explore the formation and evolution of the black hole X-ray binary system M33 X-7. In particular, we examine whether accounting for systematic errors in the stellar parameters inherent to single star models, as well as the uncertainty in the distance to M33, can explain the discrepancy between the observed and expected luminosity of the ~70 solar masses companion star. Our analysis assumes no prior interactions between the companion star and the black hole progenitor. We use four different stellar evolution codes, modified to include a variety of current stellar wind prescriptions. For the models satisfying the observational constraints on the donor star's effective temperature and luminosity, we recalculate the black hole mass, the orbital separation, and the mean X-ray luminosity. Our best model, satisfying simultaneously all observational constraints except the observationally inferred companion mass, consists of a ~13 solar masses black hole and a ~54 solar masses companion star. We conclude that a star with the observed mass and luminosity can not be explained via single star evolution models, and that a prior interaction between the companion star and the black hole progenitor should be taken into account.

One of the primary goals when studying stellar systems with neutron stars has been to reveal the physical properties of progenitors and understand how neutron star spins and birth kicks are determined. Over the years a consensus understanding had been developed, but recently some of the basic elements of this understanding are being challenged by current observations of some binary systems and their theoretical interpretation. In what follows we review such recent developments and highlight how they are interconnected; we particularly emphasize some of the assumptions and caveats of theoretical interpretations and examine their validity (e.g., in connection to the unknown radial velocities of pulsars or the nuances of multi-dimensional statistical analysis). The emerging picture does not erase our earlier understanding; instead it broadens it as it reveals additional pathways for neutron star formation and evolution: neutron stars probably form at the end of both core collapse of Fe cores of massive stars and electron-capture supernovae of ONeMg cores of lower-mass stars; birth kicks are required to be high (well in excess of 100 km/s) for some neutron stars and low (< 100 km/s) for others depending on the formation process; and spin up may occur not just through Roche-lobe overflow but also through wind accretion or phases of hypercritical accretion during common envelope evolution.