We analyze binary population models of double-neutron stars and compare results to the accurately measured orbital periods and eccentricities of the eight known such systems in our Galaxy. In contrast to past similar studies, we especially focus on the dominant evolutionary channels (we identify three); for the first time, we use a detailed understanding of the evolutionary history of three double neutron stars as actual constraints on the population models. We find that the evolutionary constraints derived from the double pulsar are particularly tight, and less than half of the examined models survive the full set of constraints. The top-likelihood surviving models yield constraints on the key binary evolution parameters, but most interestingly reveal (i) the need for electron-capture supernovae from relatively low-mass degenerate, progenitor cores, and (ii) the most likely evolutionary paths for the rest of the known double neutron stars. In particular, we find that J1913+16 likely went through a phase of Case BB mass transfer, and J1906+0746 and J1756-2251 are consistent with having been formed in electron-capture supernovae.

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.

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.

We model a circular mass-transferring binary system to calculate the exchange of angular momentum between stellar spins and the orbit due to direct impact of the mass transfer stream onto the surface of the accretor. We simulate mass transfer by calculating the ballistic motion of a point mass ejected from the $L_1$ point of the donor star, conserving the total linear and angular momentum of the system, and treating the stars as uniform density spheres with main sequence radii determined by their masses. We show that, contrary to previous assumptions in the literature, direct impact does not always act as a sink of orbital angular momentum and may in fact increase it by facilitating the transfer of angular momentum from the spin of the donor to the orbit. Here, we show an example of the exchange of angular momentum, as well as a measure of the orbital angular momentum changes for a variety of binary star systems with main sequence components.

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.

The rapid circularization and synchronization of the stellar components in an eccentric binary system at the onset of Roche lobe overflow (RLO) is a fundamental assumption common to all binary stellar evolution and population synthesis codes, even though the validity of this assumption is questionable both theoretically and observationally. Here we calculate the evolution of the orbital elements of an eccentric binary through the direct three-body integration of a massive particle ejected through the inner Lagrangian point of the donor star at periastron. The trajectory of this particle leads to three possible outcomes: direct accretion (DA) onto the companion star within a single orbit, self-accretion (SA) back onto the donor star within a single orbit, or a quasi-periodic orbit around the companion star. We calculate the secular evolution of the binary orbit in the first two cases and conclude that DA can increase or decrease the orbital semi-major axis and eccentricity, while SA always decreases the orbital both orbital elements. In cases where mass overflow contributes to circularizing the orbit, circularization can set in on timescales as short as a few per cent of the mass transfer timescale. In cases where mass overflow increases the eccentricity, the orbital evolution is governed by competition between mass overflow and tidal torques. In the absence of tidal torques, mass overflow resulting in DI can lead to substantially subsynchronously rotating donor stars. Contrary to common assumptions, DI furthermore does not always provide a strong sink of orbital angular momentum in close mass-transferring binaries; in fact we instead find that a significant part can be returned to the orbit during the particle orbit. The formulation presented here can be combined with stellar and binary evolution codes to generate a better picture of the evolution of eccentric, RLO binary star systems.

Since the discovery of the first double neutron star (DNS) system in 1975 by Hulse and Taylor, there are currently 8 confirmed DNS in our galaxy. For every system, the masses of both neutron stars, the orbital semi- major axis and eccentricity are measured, and proper motion is known for half of the systems. Using the orbital parameters and kinematic information, if available, as constraints for all system, we investigate the immediate progenitor mass of the second-born neutron star and the magnitude of the supernova kick it received at birth, with the primary goal to understand the core collapse mechanism leading to neutron star formation. Compared to earlier studies, we use a novel method to address the uncertainty related to the unknown radial velocity of the observed systems. For PSR B1534+12 and PSR B1913+16, the kick magnitudes are 150 - 270 km/s and 190 - 450 km/s (with 95% confidence) respectively, and the progenitor masses of the 2nd born neutron stars are 1.3 - 3.4 Msun and 1.4 - 5.0 Msun (95%), respectively. These suggest that the 2nd born neutron star was formed by an iron core collapse supernova in both systems. For PSR J0737-3039, on the other hand, the kick magnitude is only 5 - 120 km/s (95%), and the progenitor mass of the 2nd born neutron star is 1.3 - 1.9 Msun (95%). Because of the relatively low progenitor mass and kick magnitude, the formation of the 2nd born neutron star in PSR J0737-3039 is potentially connected to an electron capture supernova of a massive O - Ne - Mg white dwarf. For the remaining 5 Galactic DNS, the kick magnitude ranges from several tens to several hundreds of km/s, and the progenitor mass of the 2nd formed neutron star can be as low as ~1.5 Msun, or as high as ~8 Msun. Therefore in these systems, it is not clear which type of supernova is more likely to form the 2nd neutron star.

Using a population synthesis technique, we have calculated detailed models of the present-day field population of objects that have resulted from the merger of a giant primary and a main-sequence or brown dwarf secondary during common-envelope evolution. We used a grid of 116 stellar and 32 low-mass/brown dwarf models, a crude model of the merger process, and followed the angular momentum evolution of the binary orbit and the primary's rotation prior to merger, as well as the merged object's rotation after the merger. We find that present-day merged objects that are observable as giant stars or core-helium burning stars in our model population constitute between 0.24% and 0.33% of the initial population of ZAMS binaries, depending upon the input parameters chosen. The median projected rotational velocity of these merged objects is ~16 km/sec, an order of magnitude higher than the median projected rotational velocity in a model population of normal single stars calculated using the same stellar models and initial mass function. The masses of the merged objects are typically less than ~2 solar masses, with a median mass of 1.28 solar masses, which is slightly more than, but not significantly different from, their normal single star counterparts. The luminosities in our merged object population range from ~10-100 solar luminosities, with a strong peak in the luminosity distribution at ~60 solar luminosities, since the majority of the merged objects (57%) lie on the horizontal branch at the present epoch. The results of our population synthesis study are discussed in terms of possible observational counterparts either directly involving the high rotational velocity of the merger product or indirectly, via the effect of rotation on envelope abundances and on the amount and distribution of circumstellar matter.

We apply population synthesis techniques to calculate the present day population of post-common envelope binaries (PCEBs) for a range of theoretical models describing the common envelope (CE) phase. Adopting the canonical energy budget approach we consider models where the ejection efficiency, $\alpha_{\rmn{CE}}$ is either a constant, or a function of the secondary mass. We obtain the envelope binding energy from detailed stellar models of the progenitor primary, with and without the thermal and ionization energy, but we also test a commonly used analytical scaling. We also employ the alternative angular momentum budget approach, known as the $\gamma$-algorithm. We find that a constant, global value of $\alpha_{\rmn{CE}} \ga 0.1$ can adequately account for the observed population of PCEBs with late spectral-type secondaries. However, this prescription fails to reproduce IK Pegasi, which has a secondary with spectral type A8. We can account for IK Pegasi if we include thermal and ionization energy of the giant's envelope, or if we use the $\gamma$-algorithm. However, the $\gamma$-algorithm predicts local space densities that are 1 to 2 orders of magnitude greater than estimates from observations. In contrast, the canonical energy budget prescription with an initial mass ratio distribution that favours unequal initial mass ratios gives a local space density which is in good agreement with observations, and best reproduces the observed distribution of PCEBs. Finally, all models fail to reproduce the sharp decline for orbital periods, $P_{\rmn{orb}} \ga 1$ d in the orbital period distribution of observed PCEBs, even if we take into account selection effects against systems with long orbital periods and early spectral-type secondaries.

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, allowing for both mass and angular momentum loss from the system. Adopting a delta function mass transfer rate at the periastron of the binary orbit, we find that, depending on the initial binary properties at the onset of mass transfer, the orbital semi-major axis and eccentricity can either increase or decrease at a rate linearly proportional to the magnitude of the mass transfer rate at periastron. The range of initial binary mass ratios and eccentricities that leads to increasing orbital semi-major axes and eccentricities broadens with increasing degrees of mass loss from the system and narrows with increasing orbital angular momentum loss from the binary. Comparison with tidal evolution timescales shows that the usual assumption of rapid circularization at the onset of mass transfer in eccentric binaries is not justified, irrespective of the degree of systemic mass and angular momentum loss. This work extends our previous results for conservative mass transfer in eccentric binaries and can be incorporated into binary evolution and population synthesis codes to model non-conservative mass transfer in eccentric binaries.

We study tidal interactions in white dwarf binaries in the limiting case of quasi-static tides. The formalism is valid for arbitrary orbital eccentricities and therefore applicable to white dwarf binaries in the Galactic disk as well as globular clusters. In the quasi-static limit, the total perturbation of the gravitational potential shows a phase shift with respect to the position of the companion, the magnitude of which is determined primarily by the efficiency of energy dissipation through convective damping. We determine rates of secular evolution of the orbital elements and white dwarf rotational angular velocity for a 0.3 solar mass helium white dwarf in binaries with orbital frequencies in the LISA gravitational wave frequency band and companion masses ranging from 0.3 to 10^5 solar masses. The resulting tidal evolution time scales for the orbital semi-major axis are longer than a Hubble time, so that convective damping of quasi-static tides need not be considered in the construction of gravitational wave templates of white dwarf binaries in the LISA band. Spin-up of the white dwarf, on the other hand, can occur on time scales of less than 10Myr, provided that the white dwarf is initially rotating with a frequency much smaller than the orbital frequency. For semi-detached white dwarf binaries spin-up can occur on time scales of less than 1Myr. Nevertheless, the time scales remain longer than the orbital inspiral time scales due to gravitational radiation, so that the degree of asynchronism in these binaries increases. As a consequence, tidal forcing eventually occurs at forcing frequencies beyond the quasi-static tide approximation. For the shortest period binaries, energy dissipation is therefore expected to take place through dynamic tides and resonantly excited g-modes.

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.

We propose a physically motivated and self-consistent prescription for the modeling of transient neutron star (NS) low-mass X-ray binary (LMXB) properties, such as duty cycle (DC), outburst duration and recurrence time. We apply this prescription to the population synthesis (PS) models of field LMXBs presented by Fragos et al. (2008), and compare the transient LMXB population to the Chandra X-ray survey of the two elliptical galaxies NGC 3379 and NGC 4278, which revealed several transient sources (Brassington et al., 2008, 2009). We are able to exclude models with a constant DC for all transient systems, while models with a variable DC based on the properties of each system are consistent with the observed transient populations. We predict that the majority of the observed transient sources in these two galaxies are LMXBs with red giant donors. Our comparison suggests that LMXBs formed through evolution of primordial field binaries are dominant in globular cluster (GC) poor elliptical galaxies, while they still have a significant contribution in GC rich ones.

We quantify an evolutionary channel for single sdB stars based on mergers of binaries containing a red giant star and a lower mass main sequence or brown dwarf companion in our Galaxy. Population synthesis calculations that follow mergers during the common envelope phase of evolution of such systems reveal a population of rapidly rotating horizontal branch stars with a distribution of core masses between 0.32 Mo - 0.7 Mo that is strongly peaked between 0.47 Mo - 0.54 Mo. The high rotation rates in these stars are a natural consequence of the orbital angular momentum deposition during the merger and the subsequent stellar contraction of the merged object from the tip of the red giant branch. We suggest that centrifugally enhanced mass loss facilitated by the rapid rotation of these stars may lead to the formation of single sdB stars for some of these objects.

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.

We apply population synthesis techniques to calculate the present day number of two types of white dwarf-main sequence star (WDMS) binaries within the cataclysmic variable period gap. The first are post-common envelope binaries with secondary stars that have masses between 0.17 and 0.36 Msun (gPCEBs), such that they will commence mass transfer within the period gap. The second type are systems that were CVs at some point in their past, but detached once they evolved down in orbital period to ~3 h as a consequence of disrupted magnetic braking, and are crossing the period gap via gravitational radiation (dCVs). Full population synthesis calculations are performed where we assume either constant, global values of the common envelope ejection efficiency, or consider the ejection efficiency as a function of secondary mass. Several forms of magnetic braking are also considered. We predict an excess of dCVs over gPCEBs within the period gap of ~4 to ~13 assuming an ejection efficiency between 0.1 and 0.6, and a flat initial mass ratio distribution. This excess is revealed as a prominent peak at the location of the period gap in the orbital period distribution of the combined gPCEB and dCV population. We suggest that if such a feature is observed in the orbital period distribution of an observed sample of short orbital period WDMS binaries, this would strongly corroborate the disruption of magnetic braking.

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.

In globular clusters, dynamical interactions give rise to a population of eccentric double white dwarfs detectable by the Laser Interferometer Space Antenna (LISA) up to the Large Magellanic Cloud. In this Letter, we explore the detectability of periastron precession in these systems with LISA. Unlike previous investigations, we consider contributions due to tidal and rotational distortions of the binary components in addition to general relativistic contributions to the periastron precession. At orbital frequencies above a few mHz, we find that tides and stellar rotation dominate, opening up a possibly unique window to the study of the interior and structure of white dwarfs.

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.

We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter when moving through the envelope of a giant companion.

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass transfer in eccentric binaries, assuming conservation of total system mass and orbital angular momentum. Assuming a delta function mass transfer rate centered at periastron, we find rates of secular change of the orbital semi-major axis and eccentricity which are linearly proportional to the magnitude of the mass transfer rate at periastron. The rates can be positive as well as negative, so that the semi-major axis and eccentricity can increase as well as decrease in time. Adopting a delta-function mass-transfer rate of $10^{-9} M_\sun {\rm yr}^{-1}$ at periastron yields orbital evolution timescales ranging from a few Myr to a Hubble time or more, depending on the binary mass ratio and orbital eccentricity. Comparison with orbital evolution timescales due to dissipative tides furthermore shows that tides cannot, in all cases, circularize the orbit rapidly enough to justify the often adopted assumption of instantaneous circularization at the onset of mass transfer. The formalism presented can be incorporated in binary evolution and population synthesis codes to create a self-consistent treatment of mass transfer in eccentric binaries.

We investigate the existence and properties of equipotential surfaces and Lagrangian points in non-synchronous, eccentric binary star and planetary systems under the assumption of quasi-static equilibrium. We adopt a binary potential that accounts for non-synchronous rotation and eccentric orbits, and calculate the positions of the Lagrangian points as functions of the mass ratio, the degree of asynchronism, the orbital eccentricity, and the position of the stars or planets in their relative orbit. We find that the geometry of the equipotential surfaces may facilitate non-conservative mass transfer in non-synchronous, eccentric binary star and planetary systems, especially if the component stars or planets are rotating super-synchronously at the periastron of their relative orbit. We also calculate the volume-equivalent radius of the Roche lobe as a function of the four parameters mentioned above. Contrary to common practice, we find that replacing the radius of a circular orbit in the fitting formula of Eggleton (1983) with the instantaneous distance between the components of eccentric binary or planetary systems does not always lead to a good approximation to the volume-equivalent radius of the Roche-lobe. We therefore provide generalized analytic fitting formulae for the volume-equivalent Roche lobe radius appropriate for non-synchronous, eccentric binary star and planetary systems. These formulae are accurate to better than 1% throughout the relevant 2-dimensional parameter space that covers a dynamic range of 16 and 6 orders of magnitude in the two dimensions.

Current observations of double neutron stars provide us with a wealth of information that we can use to investigate their evolutionary history and the physical conditions of neutron star formation. Understanding this history and formation conditions further allow us to make theoretical predictions for the formation of other double compact objects with one or two black hole components and assess the detectability of such systems by ground-based gravitational-wave interferometers. In this paper we summarize our group's body of work in the past few years and we place our conclusions and current understanding in the framework of other work in this area of astrophysical research.

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.

The population of non magnetic cataclysmic variables evolving under the influence of a circumbinary disk is investigated for systems above the upper edge of the period gap at orbital periods greater than 2.75hr. For a fractional mass input rate into the disk, corresponding to 3e-4 of the mass transfer rate, the model systems exhibit a bounce at orbital periods greater than 2.75hr. The simulations reveal that (1) some systems can exist as dwarf nova type systems throughout their lifetime, (2) dwarf nova type systems can evolve into nova-like systems as their mass transfer rate increases with increasing circumbinary disk mass, and (3) nova-like systems can evolve back into dwarf nova systems during their postbounce evolution to longer orbital periods. Among these subclasses, nova-like cataclysmic variables would be the best candidates to search for circumbinary disks at wavelengths greater than 10 micron. The theoretical orbital period distribution of our population synthesis model is in reasonable accord with the combined population of dwarf novae and nova-like systems above the period gap, suggesting the possibility that systems with unevolved donors need not detach and evolve below the period gap as in the disrupted magnetic braking model. The resulting population furthermore reveals the possible presence of systems with small mass ratios (a property of systems exhibiting superhump phenomena at long orbital periods) and a preference of O/Ne/Mg white dwarfs in dwarf nova systems in comparison to nova-like systems. The importance of observational bias in accounting for the differing populations is examined, and it is shown that an understanding of these effects is necessary in order to confront the theoretical distributions with the observed ones in a meaningful manner. (abridged)

We revisit the formation of PSR J0737-3039, taking into account the most recent observational constraints. We show that the most likely kick velocity and progenitor parameters depend strongly on the consideration of the full five-dimensional PDF for the magnitude and direction of the kick velocity imparted to pulsar B at birth, the mass of pulsar B's pre-supernova helium star progenitor, and the pre-supernova orbital separation, and on the adopted prior assumptions. The priors consist of the transverse systemic velocity, the age of the system, and the treatment of the unknown radial velocity. Since the latter cannot be determined from observation, we adopt a statistical approach and use theoretical radial-velocity distributions obtained from population synthesis calculations for coalescing double neutron stars. We find that the prior assumptions about the pre-supernova helium star mass affect the derived most likely parameters significantly: when the minimum helium star mass required for neutron star formation is assumed to be 2.1Msun, the most likely kick velocity ranges from 70-180km/s; when masses lower than 2.1Msun are assumed to allow neutron star formation, the most likely kick velocity can be as low as a few km/s, although the majority of the considered models still yield most likely kick velocities of 50-170km/s. We also show that the proximity of the double pulsar to the Galactic plane and the small proper motion do not pose stringent constraints on the kick velocity and progenitor mass of pulsar B. Instead, the constraints imposed by the orbital dynamics of asymmetric supernova explosions turn out to be much more restrictive. We conclude that the currently available observational constraints cannot be used to favor a specific core-collapse and neutron star formation mechanism. (abridged)

Using a comprehensive binary population synthesis scheme, we investigate the statistical properties of a sample of eclipsing binaries that is detectable by an idealised extrasolar planet transit survey with specifications broadly similar to those of the SuperWASP (Wide Angle Search for Planets) project. In this idealised survey the total number of detectable single stars in the Galactic disc is of the order of 10^6-10^7, while, for a flat initial mass ratio distribution, the total number of detectable eclipsing binaries is of the order of 10^4-10^5. The majority of the population of detectable single stars is made up of main-sequence stars (60%), horizontal-branch stars (20%), and giant-branch stars (10%). The largest contributions to the population of detectable eclipsing binaries stem from detached double main-sequence star binaries (60%), detached giant-branch main-sequence star binaries (20%), and detached horizontal-branch main-sequence star binaries (10%). The ratio of the number of eclipsing binaries to the number of single stars detectable by the idealised SuperWASP survey varies by less than a factor of 2.5 across the sky, and decreases with increasing Galactic latitude. It is found to be largest in the direction of the Galactic longitude l=-7.5deg and the Galactic latitude b=-22.5deg. We also show that the fractions of systems in different subgroups of eclipsing binaries are sensitive to the adopted initial mass ratio distribution, which is one of the poorest constrained input parameters in present-day binary population synthesis calculations. This suggests that once statistically meaningful results from transit surveys are available, they will be able to significantly improve the predictive power of population synthesis studies of interacting binaries and related objects. (abridged)

The basic theory of dynamic tides in close binaries is reviewed. Particular attention is paid to resonances between dynamic tides and free oscillation modes and to the role of the apsidal-motion rate in probing the internal structure of binary components. The discussed effects are generally applicable to stars across the entire Hertzsprung-Russell diagram, including the binary OB-stars discussed at this meeting.

The population synthesis of cataclysmic variables below the period is investigated. A grid of detailed binary evolutionary sequences has been calculated and included in the simulations to take account of additional angular momentum losses beyond that associated with gravitational radiation and mass loss, due to nova outbursts, from the system. As a specific example, we consider the effect of a circumbinary disk to gain insight into the ingredients necessary to reproduce the observed orbital period distribution. The resulting distributions show that the period minimum lies at about 80 minutes with the number of systems monotonically increasing with increasing orbital period to a maximum near 90 minutes. There is no evidence for an accumulation of systems at the period minimum which is a common feature of simulations in which only gravitational radiation losses are considered. The period distribution is found to be fairly flat for orbital periods ranging from about 85 to 120 minutes. The steepness of the lower edge of the period gap can be reproduced, for example, by an input of systems at periods near 2.25 hrs due to a flow of cataclysmic variable binary systems from orbital periods longer than 2.75 hrs. The good agreement with the cumulated distribution function of observed systems within the framework of our model indicates that the angular momentum loss by a circumbinary disk or a mechanism which mimics its features coupled with a weighting factor to account for selection effects in the discovery of such systems and a flow of systems from above the period gap to below the period gap are important ingredients for understanding the overall period distribution of cataclysmic variable binary systems.

We examine the nature of RXJ0806.3+1527 and show that it is possible to reconcile the observed period decrease and X-ray luminosity with the transfer of mass between two white dwarfs provided that: either the system is (i) still in the early and short-lived (less than ~100yr) stages of mass transfer due to atmospheric Roche-lobe overflow, or (ii) in a standard, long-term, quasi-stationary mass-transfer phase that is significantly (~90%) non-conservative and the conversion of accretion energy to X-rays is quite inefficient. In either of the two cases and for a wide range of physical parameters, we find that orbital angular momentum is lost from the system at a rate that is a factor of a few (less than ~4) higher than the rate associated with the emission of gravitational waves. Although the physical origin of this extra angular momentum loss is not clear at present, it should be taken into account in the consideration of RXJ0806.3+1527 as a verification Galactic source for LISA.

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)

In view of the recent measurement of the scintillation velocity for PSR J0737-3039, we examine the complete set of constraints imposed on the pulsar B natal kicks (magnitude and orientation) and predict the most favorable pulsar kick velocity and spin tilt for both isotropic and polar kicks. Our analysis takes into account both currently unknown parameters: the orientation of the orbital plane on the sky (Omega) and the radial component of the systemic velocity (V_r). Assuming that the system's peculiar velocity is entirely due to the second supernova explosion, we find that the system may have crossed the Galactic plane multiple times since the birth of the second neutron star and that the post-supernova peculiar velocity could have been as high as 1200km/s. We also confirm the absolute lower and upper limits on the physical parameters derived in our earlier study. For specific combinations of the two unknown parameters Omega and V_r, however, we find much tighter constraints on the pre-supernova binary configuration and natal kicks imparted to pulsar B, as well as on the age of system. Once Omega is measured in the coming year, it will be straightforward to use the results presented here to further constrain the natal kicks and the spin-tilt predictions. We complete our comprehensive study and derive similar constraints and spin-tilt predictions for PSR B1534+12, where the only free parameter is V_r. Lastly, for PSR B1913+16, we update the progenitor and kick constraints using the measured pulsar spin tilt and allowing for Roche-lobe overflow from the progenitor of the pulsar companion.

We calculate the evolution of planets undergoing a strong tidal encounter using smoothed particle hydrodynamics (SPH), for a range of periastron separations. We find that outside the Roche limit, the evolution of the planet is well-described by the standard model of linear, non-radial, adiabatic oscillations. If the planet passes within the Roche limit at periastron, however, mass can be stripped from it, but in no case do we find enough energy transferred to the planet to lead to complete disruption. In light of the three new extrasolar planets discovered with periods shorter than two days, we argue that the shortest-period cases observed in the period-mass relation may be explained by a model whereby planets undergo strong tidal encounters with stars, after either being scattered by dynamical interactions into highly eccentric orbits, or tidally captured from nearly parabolic orbits. Although this scenario does provide a natural explanation for the edge found for planets at twice the Roche limit, it does not explain how such planets will survive the inevitable expansion that results from energy injection during tidal circularization.

In this talk, we review some recent advances in the theory of dynamic tides in close binaries. We particularly focus on the effects of resonances of dynamic tides with free oscillation modes and on the role of dynamic tides in the comparison of theoretically predicted and observationally inferred apsidal-motion rates.

We present an updated set of constraints for the progenitor of PSR J0737-3039 and for the natal kicks imparted to pulsar B taking into account both the evolutionary and kinematic history of the double neutron star. For this purpose, we use recently reported scintillation velocity measurements to trace the motion of the system in the Galaxy backwards in time as a function of the unknown orientation Omega of the systemic velocity projected on plane of the sky as well as the unknown radial velocity Vr. The absolute limits on the orbital separation and the mass of pulsar B's helium star progenitor just before its supernova explosion are 1.2Rsun < A0 < 1.7Rsun and 2.1Msun < M0 < 4.7Msun. The kick velocity is constrained to be between 60km/s and 1660km/s and to be misaligned from the pre-SN orbital angular axis (which could be associated with pulsar B's spin axis) by least 25 degrees. We also derive probability distribution functions for the kick velocity imparted to pulsar B and for the misalignment angle between pulsar A's spin and the post-supernova orbital angular momentum for both isotropic and polar kicks. The most probable values of both quantities depend sensitively on the unknown radial velocity. In particular, tilt angles lower than 30-50 degrees tend to be favored for current radial velocities of less than 500km/s in absolute value, while tilt angles higher than 120 degrees tend to be associated with radial velocities in excess of 1000km/s in absolute value.

We present results of a binary population synthesis study on the orbital period distribution of wide binary millisecond pulsars forming through four evolutionary channels. In three of the channels, the progenitor of the millisecond pulsar undergoes a common envelope phase prior to the supernova explosion which gives birth to the neutron star. In the fourth channel, the primary avoids the common-envelope phase and forms a neutron star when it ascends the asymptotic giant branch. The four formation channels yield an orbital period distribution which typically shows a short-period peak below 10 days, a long-period peak around 100 days, and a cut-off near 200 days. The agreement with the orbital period distribution of observed binary millisecond pulsars in the Galactic disk is best when the common-envelope ejection is efficient, the mass-transfer phase responsible for spinning up the pulsar is highly non-conservative, and no or moderate supernova kicks are imparted to neutron stars at birth.

We considered the formation of detached white dwarf main-sequence star (WDMS) binaries through seven evolutionary channels subdivided according to the evolutionary process that gives rise to the formation of the white dwarf or its helium-star progenitor: dynamically stable Roche-lobe overflow (Algol-type evolution), dynamically unstable Roche-lobe overflow (common-envelope evolution), or stellar winds (single star evolution). We examine the sensitivity of the population to changes in the amount of mass lost from the system during stable Roche-lobe overflow, the common-envelope ejection efficiency, and the initial mass ratio or initial secondary mass distribution. In the case of a flat initial mass ratio distribution, the local space density of WDMS binaries is of the order of 10^{-3}/pc^3. This number decreases to 10^{-4}/pc^3 when the initial mass ratio distribution is approximately proportional to the inverse of the initial mass ratio. More than 75% of the WDMS binary population stems from wide systems in which both components evolve as if they were single stars. The remaining part of the population is dominated by systems in which the white dwarf is formed in a common-envelope phase. The birthrate of WDMS binaries forming through a common-envelope phase is about 10 times larger than the birthrate of WDMS binaries forming through a stable Roche-lobe overflow phase. The ratio of the number of helium white dwarf systems to the number of carbon/oxygen or oxygen/neon/magnesium white dwarf systems derived from large samples of observed WDMS binaries by, e.g., future planet-search missions such as SuperWASP, COROT, and Kepler may furthermore constrain the common-envelope ejection efficiency.

A strongly relativistic binary pulsar has been recently discovered with the 64m Parkes telescope (Burgay et al. 2003). Here we use the measured properties of this binary (masses and orbital characteristics as well as age estimates), and we derive the complete set of constraints imposed on the physical properties of the binary pulsar progenitor right before the second supernova explosion. We find that: (i) according to our current understanding of neutron-star formation, the helium-rich progenitor of the second neutron star is most likely overflowing its Roche lobe; (ii) the neutron-star kick magnitude is constrained in the range 60-1560 km/s, with the most probable value being equal to 150 km/s. While the first conclusion is in agreement with Dewi & van den Heuvel (2003), our upper limit on the kick magnitude is significantly larger than that derived by these authors. We find that the difference arises because Dewi & van den Heuvel (2003) inadvertently neglected to consider kicks directed out of the pre-supernova orbital plane.

The parameter space favourable for the resonant excitation of free oscillation modes by dynamic tides in close binary components is explored using qualitative considerations to estimate the order of magnitude of the tidal force and the frequency range covered by the tidally induced oscillations. The investigation is valid for slowly rotating stars with masses in the interval between 2 and 20 solar masses, and an evolutionary stage ranging from the beginning to the end of the main sequence. Oscillation modes with eigenfrequencies of the order of five times the inverse of the star's dynamical time scale - i.e. the lowest-order p-modes, the f-mode, and the lowest-order g-modes - are found to be outside the favourable parameter space since their resonant excitation requires orbital eccentricities that are too high for the binary to stay detached when the components pass through the periastron of their relative orbit. Resonances between dynamic tides and g-modes with frequencies of the order of half of the inverse of the star's dynamical time scale on the other hand are found to be favourable for orbital periods up to 200 times the star's dynamical time scale, provided that the binary mass ratio q is larger than 1/3 and the orbital eccentricity e is larger than approximately 0.25. This favourable range comes down to orbital periods of up to 5-12 days in the case of 2-20 solar mass zero-age main-sequence binary components, and orbital periods of up to 21-70 days in the case of terminal main-sequence binary components.

We explore the population of candidate pre-low-mass X-ray binaries in which a neutron star accretes mass from the wind of a low-mass companion (mass < 2Msun) in the framework of a binary population synthesis study. The simulated accretion-luminosity distribution shows a primary peak close to 1e31 erg/s and a secondary peak near 1e28 erg/s. The relative contribution of the two peaks depends primarily on the magnitude of the kick velocity imparted to the neutron star at birth. The secondary peak is negligible for average kick velocities larger than 200 km/s, but becomes dominant for average kick velocities smaller than 50 km/s. Regardless of the relative contributions of the two peaks, our calculations suggest that pre-low-mass X-ray binaries may provide a non-negligible contribution to the population of discrete low-luminosity X-ray sources in the Galaxy.

We report new radial velocity observations of GP Vel/HD77581, the optical companion to the eclipsing X-ray pulsar Vela X-1. Using data spanning more than two complete orbits of the system, we detect evidence for tidally induced non-radial oscillations on the surface of GP Vel, apparent as peaks in the power spectrum of the residuals to the radial velocity curve fit. By removing the effect of these oscillations (to first order) and binning the radial velocities, we have determined the semi-amplitude of the radial velocity curve of GP Vel to be K_o=22.6+/-1.5 km/s. Given the accurately measured semi-amplitude of the pulsar's orbit, the mass ratio of the system is 0.081+/-0.005. We are able to set upper and lower limits on the masses of the component stars as follows. Assuming GP Vel fills its Roche lobe then the inclination angle of the system, i=70.1+/-2.6 deg. In this case we obtain the masses of the two stars as M_x=2.27 +/-0.17 M_sun for the neutron star and M_o=27.9+/-1.3 M_sun for GP Vel. Conversely, assuming the inclination angle is i=90 deg, the ratio of the radius of GP Vel to the radius of its Roche lobe is beta=0.89+/-0.03 and the masses of the two stars are M_x=1.88+/-0.13 M_sun and M_o=23.1+/-0.2 M_sun. A range of solutions between these two sets of limits is also possible, corresponding to other combinations of i and beta. In addition, we note that if the zero phase of the radial velocity curve is allowed as a free parameter, rather than constrained by the X-ray ephemeris, a significantly improved fit is obtained with an amplitude of 21.2+/-0.7 km/s and a phase shift of 0.033+/-0.007 in true anomaly. The apparent shift in the zero phase of the radial velocity curve may indicate the presence of an additional radial velocity component at the orbital period.

The paper is devoted to the analysis of the systematic deviations between the classical formula for the rate of secular apsidal motion in close binaries and the formula established within the framework of the theory of dynamic tides. The basic behaviour is found to be governed by the dependency of the forcing frequencies on the orbital period and the rotational angular velocity, irrespective of the stellar model. For binaries with short orbital periods a good agreement between the classical apsidal-motion formula and the formula established within the framework of the theory of dynamic tides is only possible for certain values of the rotational angular velocity.

This investigation is devoted to the effects of nonadiabatic resonant dynamic tides generated in a uniformly rotating stellar component of a close binary. The companion is considered to move in a fixed Keplerian orbit, and the effects of the centrifugal force and the Coriolis force are neglected. Semi-analytical solutions for the linear, nonadiabatic resonant dynamic tides are derived by means of a two-time variable expansion procedure. The solution at the lowest order of approximation consists of the resonantly excited oscillation mode and displays a phase shift with respect to the tide-generating potential. Expressions are established for the secular variations of the semi-major axis, the orbital eccentricity, and the star's angular velocity of rotation caused by the phase shift. The orders of magnitude of these secular variations are considerably larger than those derived earlier by Zahn (1977) for the limiting case of dynamic tides with small frequencies. For a 5 solar mass ZAMS star, an orbital eccentricity e = 0.5, and orbital periods in the range from 2 to 5 days, numerous resonances of dynamic tides with second-degree lower-order gravity modes are seen to induce secular variations of the semi-major axis, the orbital eccentricity, and the star's angular velocity of rotation with time scales shorter than the star's nuclear life time.

Four evolutionary channels leading to the formation of wide binary millisecond pulsars are investigated. The majority of binary millisecond pulsars are found to descend from systems in which the most massive component undergoes a common-envelope phase prior to the supernova explosion leading to the birth of the neutron star. The orbital period distribution of simulated samples of wide binary millisecond pulsars is compared with the observed distribution of Galactic binary millisecond pulsars for a variety of parameters describing the formation and evolution of binaries. The distribution functions typically show a short-period peak below 10 days and a long-period peak around 100 days. The observed distribution is best reproduced by models with highly non-conservative mass transfer, common-envelope efficiencies equal to or larger than unity, a critical mass ratio for the delayed dynamical instability larger than 3, and no or moderate supernova kicks at the birth of the neutron star. Few systems are found with orbital periods longer than 200 days, irrespective of the accretion efficiency of neutron stars. This occurs as a result of the upper limit on the initial orbital periods beyond which the binary avoids the common-envelope phase prior to the supernova explosion of the primary.

We have discovered both intermediate-order gravity mode and low-order pressure mode pulsation in the same star, HD 209295. It is therefore both a Gamma Doradus and a Delta Scuti star, which makes it the first pulsating star to be a member of two classes. The star is a single-lined spectroscopic binary with an orbital period of 3.10575 d and an eccentricity of 0.352. Weak pulsational signals are found in both the radial velocity and line-profile variations, allowing us to show that the two highest-amplitude Gamma Doradus pulsation modes are consistent with l=1 and |m|=1. In our 280 h of BVI multi-site photometry we detected ten frequencies in the light variations, one in the Delta Scuti regime and nine in the Gamma Doradus domain. Five of the Gamma Doradus frequencies are exact integer multiples of the orbital frequency. This observation leads us to suspect they are tidally excited. Results of theoretical modeling (stability analysis, tidal excitation) were consistent with the observations. We could not detect the secondary component of the system in infrared photometry, suggesting that it may not be a main-sequence star. Archival data of HD 209295 show a strong ultraviolet excess, the origin of which is not known. The orbit of the primary is consistent with a secondary mass of M > 1.04 Msun indicative of a neutron star or a white dwarf companion.

A theoretical framework for the determination of tidally induced radial-velocity variations in a component of a close binary is presented. Both the free and the forced oscillations of the component are treated as linear, isentropic perturbations of a spherically symmetric star. Resonances between dynamic tides and free oscillation modes are taken into account by means of the formalism developed by Smeyers et al. (1998). The amplitude of the tidally induced radial-velocity variations seen by the observer depends on the orbital eccentricity and on the orbital inclination. The amplitude increases with increasing orbital eccentricity and is most sensitive to the value of the orbital inclination when 20 < i < 70. In the case of a 5 solar mass ZAMS star with a 1.4 solar mass compact companion, it is shown that resonant dynamic tides can lead to radial-velocity variations with amplitudes large enough to be detected in observations. The shape of the tidally induced radial-velocity curves varies from very irregular for orbital periods away from any resonances with free oscillation modes to sinusoidal for orbital periods close to a resonance with a free oscillation mode. Our investigation is concluded with an application to the slowly pulsating B star HD 177863 showing the possibility of resonant excitation of a high-order second-degree g-mode in this star.

The validity of the classical formula for the rate of secular apsidal motion in close binaries is investigated for a sequence of models of a 5 solar mass star ranging from the last stages of the C12 -> N14 reaction to the phase where hydrogen is exhausted in the core. For binaries with short orbital periods, the apsidal-motion rates predicted by the classical formula deviate from the rates determined within the framework of the theory of dynamic tides due to the effects of the compressibility of the stellar fluid and due to resonances of dynamic tides with free oscillation modes of the component stars (Smeyers & Willems 2001). As the star evolves on the main sequence, the deviations caused by the compressibility of the stellar fluid increase with increasing radius of the star. The additional deviations caused by the resonances are largest near the end of the core-hydrogen burning phase. Both of these deviations increase with increasing values of the orbital eccentricity.

We report the discovery of both intermediate-order gravity mode and low-order pressure mode pulsation in the same star, HD 209295. It is therefore both a gamma Doradus and a delta Scuti star, which makes it the first confirmed member of two classes of pulsating star. This object is located in a close binary system with an unknown, but likely degenerate companion in an eccentric orbit, and some of the gamma Doradus pulsation frequencies are exact integer multiples of the orbital frequency. We suggest that these pulsations are tidally excited. HD 209295 may be the progenitor of an intermediate-mass X-Ray binary.