Source author record

Gijs Nelemans

Gijs Nelemans appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2026arXiv

The first phase of mass transfer in low-mass binaries: neither stable nor a common envelope

The masses of the white dwarfs in a binary carry information about previous mass-transfer phases. The core mass -- radius relation of low-mass giants gives the size of the progenitor of a helium white dwarf at the moment it last filled its Roche lobe. Previously, we used this information for a few observed systems to propose a new mass-transfer type, based on an angular momentum balance. Our aim is to investigate if stable mass transfer instead of the angular momentum prescription is consistent with the observed double helium white dwarf masses. We reconstruct the progenitor evolution of observed double helium white dwarfs using the core mass -- radius relation and evaluate if the periods at the start of the second phases of mass transfer are consistent with the outcome of stable mass transfer. More generally, we calculate the mass distribution of double helium white dwarfs for three different progenitors scenarios: double common envelope (with parameter $αλ$), angular momentum prescription (with parameter $γ$) and stable mass transfer. We find that the observed systems are generally not consistent with stable mass transfer. Stable mass transfer leads to a tight correlation between the two white dwarf masses in a binary that is not consistent with the observed mass distribution. Double common envelope evolution is a particularly poor fit to the observations. The angular momentum prescription can populate the observed mass distribution, but not perfectly. We conclude that the first phase of mass transfer initiated on the red giant branch in low-mass systems does not generally proceed as stable mass transfer nor as common envelope, and thus is poorly understood. This may be related to the fact that for many observed binaries that have finished the first phase of mass transfer the orbit is eccentric, which is an unexpected outcome of mass transfer.

preprint2022arXiv

AutoSourceID-Light. Fast Optical Source Localization via U-Net and Laplacian of Gaussian

$\textbf{Aims}$. With the ever-increasing survey speed of optical wide-field telescopes and the importance of discovering transients when they are still young, rapid and reliable source localization is paramount. We present AutoSourceID-Light (ASID-L), an innovative framework that uses computer vision techniques that can naturally deal with large amounts of data and rapidly localize sources in optical images. $\textbf{Methods}$. We show that the AutoSourceID-Light algorithm based on U-shaped networks and enhanced with a Laplacian of Gaussian filter (Chen et al. 1987) enables outstanding performances in the localization of sources. A U-Net (Ronneberger et al. 2015) network discerns the sources in the images from many different artifacts and passes the result to a Laplacian of Gaussian filter that then estimates the exact location. $\textbf{Results}$. Application on optical images of the MeerLICHT telescope demonstrates the great speed and localization power of the method. We compare the results with the widely used SExtractor (Bertin & Arnouts 1996) and show the out-performances of our method. AutoSourceID-Light rapidly detects more sources not only in low and mid crowded fields, but particularly in areas with more than 150 sources per square arcminute.

preprint2022arXiv

Partial-envelope stripping and nuclear-timescale mass transfer from evolved supergiants at low metallicity

Stable mass transfer from a massive post-main sequence (post-MS) donor is thought to be a short-lived event of thermal-timescale mass transfer which strips the donor star of nearly its entire H-rich envelope, producing a hot helium star. This long-standing picture is based on stellar models with Hertzprung gap (HG) donors. Motivated by a finding that in low-metallicity binaries, post-MS mass transfer may instead be initiated by core-helium burning (CHeB) donors, we use the MESA code to compute grids of detailed massive binary models at three metallicities: Solar and that of the Large and the Small Magellanic Cloud (LMC, SMC). We find that metallicity strongly influences the course and outcome of mass transfer. We identify two novel types of post-MS mass transfer: (a) mass exchange on the long nuclear timescale that continues until the end of the CHeB phase, and (b) rapid mass transfer leading to detached binaries with mass-losers that are only partially stripped of their envelopes. In neither (a) or (b) does the donor become a fully stripped helium star by the end of CHeB. Boundaries between the different types of post-MS mass transfer are associated with the degree of rapid post-MS expansion and, for a given metallicity, are sensitive to the assumptions about internal mixing. At low metallicity, we predict fewer hot fully stripped stars formed through binary interactions and higher compactness of pre-supernova cores. Nuclear-timescale post-MS mass transfer suggests a strong preference for metal-poor host galaxies of ultra-luminous X-ray sources with black-hole (BH) accretors and massive donors, some of which might be the immediate progenitors of binary BH mergers. It also implies a population of interacting binaries with blue and yellow supergiant donors. Partially-stripped stars could potentially explain the puzzling nitrogen-enriched slowly-rotating (super)giants in the LMC.

preprint2022arXiv

Stellar response after stripping as a model for common-envelope outcomes

Binary neutron stars have been observed as millisecond pulsars, gravitational-wave sources, and as the progenitors of short gamma-ray bursts and kilonovae. Massive stellar binaries that evolve into merging double neutron stars are believed to experience a common-envelope episode. During this episode, the envelope of a giant star engulfs the whole binary. The energy transferred from the orbit to the envelope by drag forces or from other energy sources can eject the envelope from the binary system, leading to a stripped short-period binary. In this paper, we use one-dimensional single stellar evolution to explore the final stages of the common-envelope phase in progenitors of neutron star binaries. We consider an instantaneously stripped donor star as a proxy for the common-envelope phase and study the star's subsequent radial evolution. We determine a range of stripping boundaries which allow the star to avoid significant rapid re-expansion and which thus represent plausible boundaries for the termination of the common-envelope episode. We find that these boundaries lie above the maximum compression point, a commonly used location of the core/envelope boundary. We conclude that stars may retain fractions of a solar mass of hydrogen-rich material even after the common-envelope episode. If we consider orbital energy as the only energy source available, all of our models would overfill their Roche lobe after ejecting the envelope, whose binding energy includes gravitational, thermal, radiation, and recombination energy terms.

preprint2022arXiv

The Gravitational Wave Universe Toolbox: A software package to simulate observation of the Gravitational Wave Universe with different detectors

Context. As the importance of Gravitational Wave (GW) Astrophysics increases rapidly, astronomers in different fields and with different backgrounds can have the need to get a quick idea of which GW source populations can be detected by which detectors and with what measurement uncertainties. Aims. The GW-Toolbox is an easy-to-use, flexible tool to simulate observations on the GW universe with different detectors, including ground-based interferometers (advanced LIGO, advanced VIRGO, KAGRA, Einstein Telescope, and also customised designs), space-borne interferometers (LISA and a customised design), pulsar timing arrays mimicking the current working ones (EPTA, PPTA, NANOGrav, IPTA) and future ones. We include a broad range of sources such as mergers of stellar mass compact objects, namely black holes, neutron stars and black hole-neutron stars; and supermassive black hole binaries mergers and inspirals, Galactic double white dwarfs in ultra-compact orbit, extreme mass ratio inspirals and Stochastic GW backgrounds. Methods. We collect methods to simulate source populations and determine their detectability with the various detectors. The paper aims at giving a comprehensive description on the algorithm and functionality of the GW-Toolbox. Results. The GW-Toolbox produces results that are consistent with more detailed calculations of the different source classes and can be accessed with a website interface (gw-universe.org) or as a python package (https://bitbucket.org/radboudradiolab/gwtoolbox). In the future, it will be upgraded with more functionality.

preprint2022arXiv

The GW-Universe Toolbox II: constraining the binary black hole population with second and third generation detectors

We employ the method used by the GW-Universe Toolbox to generate a synthetic catalogue of detection of stellar mass binary black hole (BBH) mergers. We study advanced LIGO (aLIGO) and Einstein Telescope (ET) as two representatives for the 2nd and 3rd generation GW observatories, and study how GW observations of BBHs can be used to constrain the merger rate as function of redshift and masses. We also simulate the observations from a detector that is half as sensitive as the ET at design which represents an early phase of ET. Two methods are used to obtain the constraints on the source population properties from the catalogues: 1. parametric differential merger rate model and applies a Bayesian inference on the parameters; and 2. non-parametric and uses weighted Kernel density estimators. The results show the overwhelming advantages of the 3rd generation detector over the 2nd generation for the study of BBH population properties, especially at a redshifts higher than ~2, where the merger rate is believed to peak. With the simulated aLIGO catalogue, the parametric Bayesian method can still give some constraints on the merger rate density and mass function beyond its detecting horizon, while the non-parametric method lose the constraining ability completely there. We also find that, despite the numbers of detection of the half-ET can be easily compatible with full ET after a longer observation duration, the catalogue from the full ET can still give much better constraints on the population properties, due to its smaller uncertainties on the physical parameters of the GW events.

preprint2020arXiv

Massive donors in interacting binaries: effect of metallicity

Metallicity is known to significantly affect the radial expansion of a massive star: the lower the metallicity, the more compact the star, especially during its post-MS evolution. We study this effect in the context of binary evolution. Using the stellar-evolution code MESA, we computed evolutionary tracks of stars at different metallicities, exploring variations of factors known to affect the radial expansion (e.g. semiconvection, overshooting, rotation). We find observational support for an evolution in which already at metallicity $0.2Z_{\odot}$ massive stars remain relatively compact during the Hertzprung-Gap (HG) phase and most of their expansion occurs during core-helium burning (CHeB). Consequently, we show that metallicity has a strong influence on the type of mass transfer evolution in binary systems. At solar metallicity, a case-B mass transfer is initiated shortly after the end of MS, and a giant donor is almost always a rapidly expanding HG star. At lower metallicity, the parameter space for mass transfer from a more evolved CHeB star increases dramatically. This means that envelope stripping and formation of helium stars in low-metallicity environments occurs later in the evolution of the donor, implying a much shorter duration of the Wolf-Rayet phase (even by an order of magnitude) and higher final core masses. This metallicity effect is independent of the impact of metallicity-dependent stellar winds. At very low metallicities, a significant fraction of massive stars in binaries engages in the first episode of mass transfer very late into their evolution, when they already have a well-developed CO core. The remaining lifetime ($< 10^4$ yr) is unlikely to be enough to strip the entire H-rich envelope. We also briefly discuss the extremely small parameter space for mass transfer from massive convective-envelope donors in the context of binary black hole merger formation.

preprint2020arXiv

The effect of the environment-dependent IMF on the formation and metallicities of stars over the cosmic history

Recent observational and theoretical studies indicate that the stellar initial mass function (IMF) varies systematically with the environment (star formation rate - SFR, metallicity). Although the exact dependence of the IMF on those properties is likely to change with improving observational constraints, the reported trend in the shape of the IMF appears robust. We present the first study aiming to evaluate the effect of the IMF variations on the measured cosmic SFR density (SFRD) as a function of metallicity and redshift, $f_{\rm SFR}$(Z,z). We also study the expected number and metallicity of white dwarf, neutron star and black hole progenitors under different IMF assumptions. Applying the empirically driven IMF variations described by the integrated galactic IMF (IGIMF) theory, we correct $f_{\rm SFR}$(Z,z) obtained by Chruslinska & Nelemans (2019) and find lower SFRD at high redshifts as well as a higher fraction of metal-poor stars being formed. In the local Universe, our calculation applying the IGIMF theory suggests more white dwarf and neutron star progenitors in comparison with the universal IMF scenario, while the number of black hole progenitors remains unaffected.

preprint2018arXiv

High mass X-ray binaries as progenitors of gravitational wave sources

X-ray binaries with black hole (BH) accretors and massive star donors at short orbital periods of a few days can evolve into close binary BH systems (BBH) that merge within the Hubble time through stable mass transfer evolution. From observational point of view, upon the Roche-lobe overflow, such systems will most likely appear as ultra-luminous X-ray sources (ULXs). To study this connection, we compute the mass transfer phase in systems with BH accretors and massive star donors ($M > 15 \,M_\odot$) at various orbital separations and metallicities using the MESA stellar evolution code. In the case of core-hydrogen and core-helium burning donors (cases A and C of mass transfer) we find the typical duration of super-Eddington mass transfer of up to $10^6$ and $10^5 \, \rm yr$ , with rates of $10^{-6}$ and $10^{-5} \,M_\odot \, \rm yr^{-1}$ , respectively. Given that roughly $0.5$ ULXs are found per unit of star formation rate ($\,M_\odot \, \rm yr^{-1}$), and assuming that $10\%$ of all the observed ULXs form merging BBH, we estimate the rate of BBH mergers from stable mass transfer evolution to be at most $10 {\rm ~Gpc}^{-3} {\rm ~yr}^{-1}$.

preprint2016arXiv

Building Blocks of the Milky Way's Accreted Spheroid

In the $Λ$CDM model of structure formation, a stellar spheroid grows by the assembly of smaller galaxies, the so-called building blocks. Combining the Munich-Groningen semi-analytical model of galaxy formation with the high resolution Aquarius simulations of dark matter haloes, we study the assembly history of the stellar spheroids of six Milky Way-mass galaxies, focussing on building block properties such as mass, age and metallicity. These properties are compared to those of the surviving satellites in the same models. We find that the building blocks have higher star formation rates on average, and this is especially the case for the more massive objects. At high redshift these dominate in star formation over the satellites, whose star formation timescales are longer on average. These differences ought to result in a larger $α$-element enhancement from Type II supernovae in the building blocks (compared to the satellites) by the time Type Ia supernovae would start to enrich them in iron, explaining the observational trends. Interestingly, there are some variations in the star formation timescales of the building blocks amongst the simulated haloes, indicating that [$α$/Fe] abundances in spheroids of other galaxies might differ from those in our own Milky Way.

preprint2016arXiv

The OmegaWhite Survey for short period variable stars II: An overview of results from the first four years

OmegaWhite is a wide-field, high cadence, synoptic survey targeting fields in the southern Galactic plane, with the aim of discovering short period variable stars. Our strategy is to take a series of 39 s exposures in the g band of a 1 square degree of sky lasting 2 h using the OmegaCAM wide field imager on the VLT Survey Telescope (VST). We give an overview of the initial 4 years of data which covers 134 square degrees and includes 12.3 million light curves. As the fields overlap with the VLT Survey Telescope Halpha Photometric Survey of the Galactic plane and Bulge (VPHAS+), we currently have $ugriHα$ photometry for ~1/3 of our fields. We find that a significant fraction of the light curves have been affected by the diffraction spikes of bright stars sweeping across stars within a few dozen of pixels over the two hour observing time interval due to the alt-az nature of the VST. We select candidate variable stars using a variety of variability statistics, followed by a manual verification stage. We present samples of several classes of short period variables, including: an ultra compact binary, a DQ white dwarf, a compact object with evidence of a 100 min rotation period, three CVs, one eclipsing binary with an 85 min period, a symbiotic binary which shows evidence of a 31 min photometric period, and a large sample of candidate delta Sct type stars including one with a 9.3 min period. Our overall goal is to cover 400 square degrees, and this study indicates we will find many more interesting short period variable stars as a result.

preprint2016arXiv

Tiling strategies for optical follow-up of gravitational wave triggers by wide field of view telescopes

Binary neutron stars are among the most promising candidates for joint gravitational-wave and electromagnetic astronomy. The goal of this work is to investigate the strategy of using gravitational wave sky-localizations for binary neutron star systems, to search for electromagnetic counterparts using wide field of view optical telescopes. We examine various strategies of scanning the gravitational wave sky-localizations on the mock 2015-16 gravitational-wave events. We propose an optimal tiling-strategy that would ensure the most economical coverage of the gravitational wave sky-localization, while keeping in mind the realistic constrains of transient optical astronomy. Our analysis reveals that the proposed tiling strategy improves the sky-localization coverage over naive contour-covering method. The improvement is more significant for observations conducted using larger field of view telescopes, or for observations conducted over smaller confidence interval of gravitational wave sky-localization probability distribution. Next, we investigate the performance of the tiling strategy for telescope arrays and compare their performance against monolithic giant field of view telescopes. We observed that distributing the field of view of the telescopes into arrays of multiple telescopes significantly improves the coverage efficiency by as much as 50% over a single large FOV telescope in 2016 localizations while scanning around 100 sq. degrees. Finally, we studied the ability of optical counterpart detection by various types of telescopes. In Our analysis for a range of wide field-of-view telescopes we found improvement in detection upon sacrificing coverage of localization in order to achieve greater observation depth for very large field-of-view - small aperture telescopes, especially if the intrinsic brightness of the optical counterparts are weak.

preprint2015arXiv

Constraining the formation of black-holes in short-period Black-Hole Low-Mass X-ray Binaries

The formation of stellar mass black holes is still very uncertain. Two main uncertainties are the amount of mass ejected in the supernova event (if any) and the magnitude of the natal kick the black hole receives at birth (if any). Repetto et al. (2012), studying the position of Galactic X-ray binaries containing black holes, found evidence for black holes receiving high natal kicks at birth. In this Paper we extend that study, taking into account the previous binary evolution of the sources as well. The seven short-period black-hole X-ray binaries that we use, are compact binaries consisting of a low-mass star orbiting a black hole in a period less than $1$ day. We trace their binary evolution backwards in time, from the current observed state of mass-transfer, to the moment the black hole was formed, and we add the extra information on the kinematics of the binaries. We find that several systems could be explained by no natal kick, just mass ejection, while for two systems (and possibly more) a high kick is required. So unless the latter have an alternative formation, such as within a globular cluster, we conclude that at least some black holes get high kicks. This challenges the standard picture that black hole kicks would be scaled down from neutron star kicks. Furthermore, we find that five systems could have formed with a non-zero natal kick but zero mass ejected (i.e. no supernova) at formation, as predicted by neutrino-driven natal kicks.

preprint2015arXiv

Constraining the Physics of AM Canum Venaticorum Systems with the Accretion Disk Instability Model

Recent work by Levitan et al has expanded the long-term photometric database for AM CVn stars. In particular, their outburst properties are well-correlated with orbital period, and allow constraints to be placed on the secular mass transfer rate between secondary and primary if one adopts the disk instability model for the outbursts. We use the observed range of outbursting behavior for AM CVn systems as a function of orbital period to place a constraint on mass transfer rate versus orbital period P. We infer a rate ~5 x 10^{-9} Msun/yr (P/1000 s)^{-5.2}. We show the functional form so obtained is consistent with the recurrence time-orbital period relation found by Levitan et al using a simple theory for the recurrence time. Also, we predict their steep dependence of outburst duration on orbital period will flatten considerably once the longer orbital period systems have more complete observations.

preprint2015arXiv

Stellar archaeology with Gaia: the Galactic white dwarf population

Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.

preprint2014arXiv

Binary white dwarfs in the halo of the Milky Way

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

preprint2014arXiv

Constraining parameters of white-dwarf binaries using gravitational-wave and electromagnetic observations

The space-based gravitational wave (GW) detector, \emph{evolved Laser Interferometer Space Antenna} (eLISA) is expected to observe millions of compact Galactic binaries that populate our Milky Way. GW measurements obtained from the eLISA detector are in many cases complimentary to possible electro-magnetic (EM) data. In our previous papers, we have shown that the EM data can significantly enhance our knowledge of the astrophysically relevant GW parameters of the Galactic binaries, such as the amplitude and inclination. This is possible due to the presence of some strong correlations between GW parameters that are measurable by both EM and GW observations, for example the inclination and sky position. In this paper, we quantify the constraints in the physical parameters of the white-dwarf binaries, i.e. the individual masses, chirp mass and the distance to the source that can be obtained by combining the full set of EM measurements such as the inclination, radial velocities, distances and/or individual masses with the GW measurements. We find the following $2-σ$ fractional uncertainties in the parameters of interest. The EM observations of distance constrains the the chirp mass to $\sim 15-25 %$, whereas EM data of a single-lined spectroscopic binary constrains the secondary mass and the distance with factors of 2 to $\sim 40 %$. The single-line spectroscopic data complemented with distance constrains the secondary mass to $\sim 25-30%$. Finally EM data on double-lined spectroscopic binary constrains the distance to $\sim 30%$. All of these constraints depend on the inclination and the signal strength of the binary systems. We also find that the EM information on distance and/or the radial velocity are the most useful in improving the estimate of the secondary mass,inclination and/or distance.

preprint2014arXiv

Localizing gravitational wave sources with optical telescopes and combining electromagnetic and gravitational wave data

Neutron star binaries, which are among the most promising sources for the direct detection of gravitational waves (GW) by ground based detectors, are also potential electromagnetic (EM) emitters. Gravitational waves will provide a new window to observe these events and hopefully give us glimpses of new astrophysics. In this paper, we discuss how EM information of these events can considerably improve GW parameter estimation both in terms of accuracy and computational power requirement. And then in return how GW sky localization can help EM astronomers in follow-up studies of sources which did not yield any prompt emission. We discuss how both EM source information and GW source localization can be used in a framework of multi-messenger astronomy. We illustrate how the large error regions in GW sky localizations can be handled in conducting optical astronomy in the advance detector era. We show some preliminary results in the context of an array of optical telescopes called BlackGEM, dedicated for optical follow-up of GW triggers, that is being constructed in La Silla, Chile and is expected to operate concurrent to the advanced GW detectors.

preprint2014arXiv

Measuring tides and binary parameters from gravitational wave data and eclipsing timings of detached white dwarf binaries

The discovery of the most compact detached white dwarf (WD) binary SDSS J065133.33+284423.3 has been discussed in terms of probing the tidal effects in white dwarfs. This system is also a verification source for the space-based gravitational wave (GW) detector, evolved Laser Interferometer Space Antenna (eLISA) which will observe short-period compact Galactic binaries with $P_{orb}\lesssim 5$ hrs. We address the prospects of doing tidal studies using eLISA binaries by showing the fractional uncertainties in the orbital decay rate and the rate of that decay, $\dot{f}, \ddot{f}$ expected from both the GW and EM data for some of the high-$f$ binaries. We find that $\dot{f}$ and $\ddot{f}$ can be measured using GW data only for the most massive WD binaries observed at high-frequencies. Form timing the eclipses for $\sim 10$ years, we find that $\dot{f}$ can be known to $\sim 0.1% $ for J0651. We find that from GW data alone, measuring the effects of tides in binaries is (almost) impossible. We also investigate the improvement in the knowledge of the binary parameters by combining GW amplitude and inclination with EM data with and without $\dot{f}$. In our previous work we found that EM data on distance constrained 2-$σ$ uncertainty in chirp mass to $15-25%$ whereas adding $\dot{f}$ reduces it to $0.11%$. EM data on $\dot{f}$ also constrains 2-$σ$ uncertainty in distance to $35%-19%$. EM data on primary mass constrains the secondary mass $m_2$ to factors of 2 to $\sim40%$ whereas adding $\dot{f}$ reduces this to $25%$. And finally using single-line spectroscopic constrains 2-$σ$ uncertainties in both the $m_2, d$ to factors of 2 to $\sim 40%$. Adding EM data on $\dot{f}$ reduces these 2-$σ$ uncertainties to $\leq 25%$ and $6%-19%$ respectively. Thus we find that EM measurements of $\dot{f}$ and radial velocity will be valuable in constraining binary parameters.

preprint2014arXiv

Observational clues to the progenitors of Type-Ia supernovae

Type-Ia supernovae (SNe Ia) are important distance indicators, element factories, cosmic-ray accelerators, kinetic-energy sources in galaxy evolution, and endpoints of stellar binary evolution. It has long been clear that a SN Ia must be the runaway thermonuclear explosion of a degenerate carbon-oxygen stellar core, most likely a white dwarf (WD). However, the specific progenitor systems of SNe Ia, and the processes that lead to their ignition, have not been identified. Two broad classes of progenitor binary systems have long been considered: single-degenerate (SD), in which a WD gains mass from a non-degenerate star; and double-degenerate (DD), involving the merger of two WDs. New theoretical work has enriched these possibilities with some interesting updates and variants. We review the significant recent observational progress in addressing the progenitor problem. We consider clues that have emerged from the observed properties of the various proposed progenitor populations, from studies of their sites, pre- and post-explosion, from analysis of the explosions themselves, and from the measurement of event rates. The recent nearby and well-studied event, SN 2011fe, has been particularly revealing. The observational results are not yet conclusive, and sometimes prone to competing theoretical interpretations. Nevertheless, it appears that DD progenitors, long considered the underdog option, could be behind some, if not all, SNe Ia. We point to some directions that may lead to future progress.

preprint2014arXiv

The coupled effect of tides and stellar winds on the evolution of compact binaries

We follow the evolution of compact binaries under the coupled effect of tides and stellar winds until the onset of Roche-lobe overflow. These binaries contain a compact object (either a black-hole, a neutron-star, or a planet) and a stellar component. We integrate the full set of tidal equations, which are based on Hut's model for tidal evolution, and we couple them with the angular momentum loss in a stellar wind. Our aim is twofold. Firstly, we wish to highlight some interesting evolutionary outcomes of the coupling. When tides are coupled with a non-massive stellar wind, one interesting outcome is that in certain types of binaries, the stellar spin tends to reach a quasi-equilibrium state, where the effect of tides and wind are counteracting each other. When tides are coupled with a massive wind, we parametrize the evolution in terms of the decoupling radius, at which the wind decouples from the star. Even for small decoupling radii this \emph{wind braking} can drive systems on the main sequence to Roche-lobe overflow that otherwise would fail to do so. Our second aim is to inspect whether simple timescale considerations are a good description of the evolution of the systems. We find that simple timescale considerations, which rely on neglecting the coupling between tides and stellar winds, do not accurately represent the true evolution of compact binaries. The outcome of the coupled evolution of the rotational and orbital elements can strongly differ from simple timescale considerations, as already pointed out by Barker & Ogilvie 2009 in the case of short-period planetary systems.

preprint2013arXiv

A determination of the space density and birth rate of hydrogen-line (DA) white dwarfs in the Galactic Plane, based on the UVEX survey

We present a determination of the average space density and birth rate of hydrogen-line (DA) white dwarfs within a radius of 1 kpc around the Sun, based on an observational sample of 360 candidate white dwarfs with g<19.5 and (g-r)<0.4, selected from the UV-excess Survey of the Northern Galactic Plane (UVEX), in combination with a theoretical white dwarf population that has been constructed to simulate the observations, including the effects of reddening and observational selection effects. The main uncertainty in the derivation of the white dwarf space density and current birth rate lies in the absolute photometric calibration and the photometric scatter of the observational data, which influences the classification method on colours, the completeness and the pollution. Corrections for these effects are applied. We derive an average space density of hydrogen-line (DA) white dwarfs with T_eff > 10,000K (M_V<12.2) of (3.8 +/- 1.1) x 1e-4 pc^-3, and an average DA white dwarf birth rate over the last 7e7 years of (5.4 + 1.5) x 1e-13 pc^-3 yr^-1. Additionally, we show that many estimates of the white dwarf space density from different studies are consistent with each other, and with our determination here.

preprint2013arXiv

On double-degenerate type Ia supernova progenitors as supersoft X-ray sources - A population synthesis analysis using SeBa

Context: The nature of the progenitors of type Ia supernova (SN) progenitors remains unclear. While it is usually agreed that single-degenerate progenitor systems would be luminous supersoft X-ray sources (SSSs), it was recently suggested that double-degenerate progenitors might also go through a SSS phase. Aims: We aim to examine the possibility of double-degenerate progenitor systems being SSSs, and place stringent upper limits on the maximally possible durations of any SSS phases and expected number of these systems in a galactic population. Method: We employ the binary population synthesis code SeBa to examine the mass-transfer characteristics of a possible SSS phase of double-degenerate type Ia SN progenitor systems for 1) the standard SeBa assumptions, and 2) an optimistic best-case scenario. The latter case establishes firm upper limits on the possible population of supersoft source double-degenerate type Ia SN progenitor systems. Results: Our results indicate that unlike what is expected for single-degenerate progenitor systems, the vast majority of the material accreted by either pure wind mass transfer or a combination of wind and RLOF mass transfer is helium rather than hydrogen. Even with extremely optimistic assumptions concerning the mass-transfer and retention efficiencies, the average mass accreted by systems that eventually become double-degenerate type Ia SNe is small. Consequently, the lengths of time that these systems may be SSSs are short, even under optimal conditions, and the expected number of such systems in a galactic population is negligible. Conclusions: The population of double-degenerate type Ia SN progenitors that are SSSs is at least an order of magnitude smaller than the population of single-degenerate progenitors expected to be SSSs, and the supersoft X-ray behaviour of double-degenerate systems typically ceases long before the SN explosions.

preprint2013arXiv

Population synthesis of triple systems in the context of mergers of carbon-oxygen white dwarfs

Hierarchical triple systems are common among field stars yet their long-term evolution is poorly understood theoretically. In such systems Kozai cycles can be induced in the inner binary system during which the inner orbit eccentricity and the inclination between both binary orbits vary periodically. These cycles, combined with tidal friction and gravitational wave emission, can significantly affect the inner binary evolution. To investigate these effects quantitatively we perform a population synthesis study of triple systems and focus on evolutionary paths that lead to mergers of carbon-oxygen (CO) white dwarfs (WDs), which constitute an important candidate progenitor channel for type Ia supernovae (SNe Ia). We approach this problem by Monte Carlo sampling from observation-based distributions of systems within the primary mass range 1.0 - 6.5 M_Sun and inner orbit semi-major axes a_1 and eccentricities e_1 satisfying a_1 (1-e_1^2) > 12 AU, i.e. non-interacting in the absence of a tertiary component. We evolve these systems by means of a newly developed algorithm that couples secular triple dynamics with an existing binary population synthesis code. We find that the tertiary significantly alters the inner binary evolution in about 24% of all sampled systems. In particular, we find several channels leading to CO WD mergers. Amongst these is a novel channel in which a collision occurs in wide inner binary systems as a result of extremely high eccentricities induced by Kozai cycles. With assumptions on which CO WD mergers lead to a SN Ia explosion we estimate the SNe Ia delay time distribution resulting from triples and compare to a binary population synthesis study and to observations. Although we find that our triple rate is low, we have determined a lower limit of the triple-induced SNe Ia rate and further study is needed that includes triples with initially tighter inner orbits.

preprint2013arXiv

The effect of common-envelope evolution on the visible population of post-common-envelope binaries

Context. An important ingredient in binary evolution is the common-envelope (CE) phase. Although this phase is believed to be responsible for the formation of many close binaries, the process is not well understood. Aims. We investigate the characteristics of the population of post-common-envelope binaries (PCEB). As the evolution of these binaries and their stellar components are relatively simple, this population can be directly used to constraint CE evolution. Methods. We use the binary population synthesis code SeBa to simulate the current-day population of PCEBs in the Galaxy. We incorporate the selection effects in our model that are inherent to the general PCEB population and that are specific to the SDSS survey, which enables a direct comparison for the first time between the synthetic and observed population of visible PCEBs. Results. We find that selection effects do not play a significant role on the period distribution of visible PCEBs. To explain the observed dearth of long-period systems, the α-CE efficiency of the main evolutionary channel must be low. In the main channel, the CE is initiated by a red giant as it fills its Roche lobe in a dynamically unstable way. Other evolutionary paths cannot be constrained more. Additionally our model reproduces well the observed space density, the fraction of visible PCEBs amongst white dwarf (WD)- main sequence (MS) binaries, and the WD mass versus MS mass distribution, but overestimates the fraction of PCEBs with helium WD companions.

preprint2013arXiv

The Galactic Bulge Survey: completion of the X-ray survey observations

We provide the Chandra source list for the last ~quarter of the area covered by the Galactic Bulge Survey (GBS). The GBS targets two strips of 6\degr x 1\degr (12 square degrees in total), one above (1\degr<b<2\degr) and one below (-2\degr<b<-1\degr) the Galactic plane in the direction of the Galactic Center at X-ray, optical and near-infrared wavelengths. For the X-ray part of the survey we use 2 ks per Chandra pointing. We find 424 X-ray sources in the 63 Chandra observations we report on here. These sources are in addition to the 1216 X-ray sources discovered in the first part of the GBS survey described before. We discuss the characteristics and the X-ray variability of the brightest of the sources as well as the radio properties from existing radio surveys. We point out an interesting asymmetry in the number of X-ray sources as a function of their Galactic l and b coordinates which is probably caused by differences in average extinction towards the different parts of the GBS survey area.

preprint2013arXiv

Ultra-compact X-ray binaries in the Galactic Bulge

We model the present-day number and characteristics of Ultra-compact X-ray binaries (UCXBs) in the Galactic Bulge. The main objective is to compare the results with the forthcoming data from the Galactic Bulge Survey (GBS) in order to learn about formation scenarios, notably core collapse supernova rate, common envelope evolution and dynamical interactions, as well as the binary initial mass function. We use a binary population synthesis code and detailed stellar evolutionary tracks, combined with observations of the 11 known UCXBs. We predict a few hundred faint UCXBs with periods of 20-40 minutes that would be detectable by the GBS.

preprint2013arXiv

Using electromagnetic observations to aid gravitational-wave parameter estimation of compact binaries observed with LISA II: The effect of knowing the sky position

In this follow-up paper, we continue our study of the effect of using knowledge from electromagnetic observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new \textit{Laser Interferometer Space Antenna} in the low-frequency range, $10^{-4} \:\mathrm{Hz}<f<1 \:\mathrm{Hz}$. In the first paper, we have shown that the strong correlation between amplitude and inclination can be used for mildly inclined binaries to improve the uncertainty in amplitude, and that this correlation depends on the inclination of the system. In this paper we investigate the overall effect of the other orientation parameters, namely the sky position and the polarisation angle. We find that after the inclination, the ecliptic latitude of the source has the strongest effect in determining the GW parameter uncertainties. We ascertain that the strong correlation we found previously, only depends on the inclination of the source and not on the other orientation parameters. We find that knowing the sky position of the source from electromagnetic data can reduce the GW parameter uncertainty up to a factor of $\sim 2$, depending on the inclination and the ecliptic latitude of the system. Knowing the sky position and inclination can reduce the uncertainty in amplitude by a factor larger than 40. We also find that unphysical errors in the inclinations, which we found when using the Fisher matrix, can affect the corresponding uncertainties in the amplitudes, which need to be corrected.

preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

preprint2012arXiv

Gravitational-wave emission from compact Galactic binaries

Compact Galactic binaries where at least one member is a white dwarf or neutron star constitute the majority of individually detectable sources for future low-frequency space-based gravitational-wave (GW) observatories; they also form an unresolved continuum, the dominant Galactic foreground at frequencies below a few mHz. Due to the paucity of electromagnetic observations, the majority of studies of Galactic-binary populations so far have been based on population-synthesis simulations. However, recent surveys have reported several new detections of white-dwarf binaries, providing new constraints for population estimates. In this article, we evaluate the impact of revised local densities of interacting white-dwarf binaries on future GW observations. Specifically: we consider five scenarios that explain these densities with different assumptions on the formation of interacting systems; we simulate corresponding populations of detached and interacting white-dwarf binaries; we estimate the number of individually detectable GW sources and the magnitude of the confusion-noise foreground, as observed by space-based detectors with 5- and 1-Mkm arms. We confirm earlier estimates of thousands of detached-binary detections, but project only few ten to few hundred detections of interacting systems. This reduction is partly due to our assessment of detection prospects, based on the iterative identification and subtraction of bright sources with respect to both instrument and confusion noise. We also confirm earlier estimates for the confusion-noise foreground, except in one scenario that explains smaller local densities of interacting systems with smaller numbers of progenitor detached systems.

preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

preprint2012arXiv

Obscuration of Supersoft X-ray Sources by Circumbinary Material - A Way to Hide Type Ia Supernova Progenitors?

The progenitors of supernovae (SNe) type Ia are usually assumed to be either a single white dwarf (WD) accreting from a non-degenerate companion (the SD channel) or the result of two merging WDs (DD channel). However, no consensus currently exists as to which progenitor scenario is the correct one, or whether the observed SN Ia rate is produced by a combination of both channels. Unlike a DD progenitor a SD progenitor is expected to emit supersoft X-rays for a prolonged period of time (~1 Myr) as a result of the burning of accreted matter on the surface of the WD. An argument against the SD channel as a significant producer of SNe type Ia has been the lack of observed supersoft X-ray sources (SSS) and the lower-than-expected integrated soft X-ray flux from elliptical galaxies. We wish to determine if it is possible to obscure the supersoft X-ray emission from a nuclear burning white dwarf in an accreting single degenerate binary system. In case of obscured systems we wish to determine their general observational characteristics. We examine the emergent X-ray emission from a canonical SSS system surrounded by a spherically symmetric configuration of material, assuming a black body spectrum with T_BB=50 eV and L=10^38 erg/s. The circumbinary material is assumed to be of solar chemical abundances, and we leave the mechanism behind the mass loss into the circumbinary region unspecified. If steadily accreting, nuclear burning WDs are canonical SSS our analysis suggests that they can be obscured by relatively modest circumbinary mass loss rates. This may explain the discrepancy of SSS compared to the SN Ia rate inferred from observations if the SD progenitor scenario contributes significantly to the SN Ia rate. Recycled emissions from obscured systems may be visible in other wavebands than X-rays. It may also explain the lack of observed SSS in symbiotic binary systems.

preprint2012arXiv

Prospects for observing ultra-compact binaries with space-based gravitational wave interferometers and optical telescopes

Space-based gravitational wave interferometers are sensitive to the galactic population of ultra-compact binaries. An important subset of the ultra-compact binary population are those stars that can be individually resolved by both gravitational wave interferometers and electromagnetic telescopes. The aim of this paper is to quantify the multi-messenger potential of space-based interferometers with arm-lengths between 1 and 5 Gm. The Fisher Information Matrix is used to estimate the number of binaries from a model of the Milky Way which are localized on the sky by the gravitational wave detector to within 1 and 10 square degrees and bright enough to be detected by a magnitude limited survey. We find, depending on the choice of GW detector characteristics, limiting magnitude, and observing strategy, that up to several hundred gravitational wave sources could be detected in electromagnetic follow-up observations.

preprint2012arXiv

Radio sources in the Chandra Galactic Bulge survey

We discuss radio sources in the Chandra Galactic Bulge survey region. By cross-matching the X-ray sources in this field with the NVSS archival data, we find 12 candidate matches. We present a classification scheme for radio/X-ray matches in surveys taken in or near the Galactic Plane, taking into account other multi-wavelength data. We show that none of the matches found here is likely to be due to coronal activity from normal stars because the radio to X-ray flux ratios are systematically too high. We show that one of the sources could be a radio pulsar, and that one could be a planetary nebula, but that the bulk of the sources are likely to be background active galactic nuclei (AGN), with many confirmed through a variety of approaches. Several of the AGN are bright enough in the near infrared (and presumably in the optical) to use as probes of the interstellar medium in the inner Galaxy.

preprint2012arXiv

Spectroscopic follow-up of UV-excess objects selected from the UVEX survey

We present the results of the first spectroscopic follow-up of 132 optically blue UV-excess sources selected from the UV-excess survey of the Northern Galactic Plane (UVEX). The UV-excess spectra are classified into different populations and grids of model spectra are fit to determine spectral types, temperatures, surface gravities and reddening. From this initial spectroscopic follow-up 95% of the UV-excess candidates turn out to be genuine UV-excess sources such as white dwarfs, white dwarf binaries, subdwarfs type O and B, emission line stars and QSOs. The remaining sources are classified as slightly reddened main-sequence stars with spectral types later than A0V. The fraction of DA white dwarfs is 47% with reddening smaller than E(B-V)<0.7 mag. Relations between the different populations and their UVEX photometry, Galac- tic latitude and reddening are shown. A larger fraction of UVEX white dwarfs is found at magnitudes fainter than g>17 and Galactic latitude smaller than |b|<4 compared to main-sequence stars, blue horizontal branch stars and subdwarfs.

preprint2012arXiv

Supernova Type Ia progenitors from merging double white dwarfs: Using a new population synthesis model

The study of Type Ia supernovae (SNIa) has lead to greatly improved insights into many fields in astrophysics, however a theoretical explanation of the origin of these events is still lacking. We investigate the potential contribution to the SNIa rate from the population of merging double carbon-oxygen white dwarfs. We aim to develope a model that fits the observed SNIa progenitors as well as the observed close double white dwarf population. We differentiate between two scenarios for the common envelope (CE) evolution; the alpha-formalism based on the energy equation and the gamma-formalism that is based on the angular momentum equation. In one model we apply the alpha-formalism always. In the second model the gamma-formalism is applied, unless the binary contains a compact object or the CE is triggered by a tidal instability for which the alpha-formalism is used. The binary population synthesis code SeBa was used to evolve binary systems from the zero-age main sequence to the formation of double white dwarfs and subsequent mergers. SeBa has been thoroughly updated since the last publication of the content of the code. The limited sample of observed double white dwarfs is better represented by the simulated population using the gamma-formalism than the alpha-formalism. For both CE formalisms, we find that although the morphology of the simulated delay time distribution matches that of the observations within the errors, the normalisation and time-integrated rate per stellar mass are a factor 7-12 lower than observed. Furthermore, the characteristics of the simulated populations of merging double carbon-oxygen white dwarfs are discussed and put in the context of alternative SNIa models for merging double white dwarfs.

preprint2012arXiv

Theoretical Delay Time Distributions

We briefly discuss the method of population synthesis to calculate theoretical delay time distributions of type Ia supernova progenitors. We also compare the results of the different research groups and conclude that although one of the main differences in the results for single degenerate progenitors is the retention efficiency with which accreted hydrogen is added to the white dwarf core, this cannot explain all the differences.

preprint2012arXiv

Upper limits on bolometric luminosities of ten type Ia supernova progenitors from Chandra observations

We present an analysis of Chandra observations of the position of ten nearby (< 25 Mpc) type Ia supernovae, taken before the explosions. No sources corresponding to progenitors were found in any of the observations. We calculated upper limits on the bolometric luminosities of the progenitors assuming black-body X-ray spectra with temperatures of 30-150 eV. This is inspired by the fact that luminous super-soft X-ray sources have been suggested as the direct progenitors of type Ia supernovae. The upper limits of two supernovae in our sample are comparable to the luminosities of the brightest observed super-soft sources, ruling out such sources as the progenitors of these supernovae. In contrast to Liu et al (2012) we find that for SN2011fe we can rule out Eddington luminosity systems for black body temperatures as low as 40 eV. Our findings are consistent with statistical studies comparing the observed type Ia supernova rate to the number of super-soft sources or the integrated X-ray luminosity in external galaxies. This suggest that either the progenitors of type Ia supernovae are not accreting, nuclear burning white dwarfs, or that they do not look like the classical super-soft sources, e.g. because they are obscured.

preprint2012arXiv

Using electromagnetic observations to aid gravitational-wave parameter estimation of compact binaries observed with LISA

We present a first-stage study of the effect of using knowledge from electromagnetic (EM) observations in the gravitational wave (GW) data analysis of Galactic binaries that are predicted to be observed by the new Laser Interferometer Space Antenna in the low-frequency range, $10^{-4} \mathrm{Hz}<f<1 \mathrm{Hz}$. In particular, we examine the extent to which the accuracy of GW parameter estimation improves if we use available information from EM data. We do this by investigating whether correlations exist between the GW parameters that describe these binaries and whether some of these parameters are also available from EM observations. We used verification binaries, which are known as the guaranteed sources for eLISA and will test the functioning of the instrument. We find that of the seven parameters that characterise such a binary, only a few are correlated. The most useful result is the strong correlation between amplitude and inclination, which can be used to constrain the parameter uncertainty in amplitude by making use of the constraint of inclination from EM measurements. The improvement can be up to a factor of $\sim6.5$, but depends on the signal-to-noise ratio of the source data. Moreover, we find that this strong correlation depends on the inclination. For mildly face-on binaries ($ι\lesssim 45^{\circ}$), EM data on inclination can improve the estimate of the GW amplitude by a significant factor. However, for edge-on binaries ($ι\sim 90^{\circ}$), the inclination can be determined accurately from GW data alone, thus GW data can be used to select systems that will likely be eclipsing binaries for EM follow-up.

preprint2011arXiv

Stellar variability on time-scales of minutes: results from the first 5 years of the Rapid Temporal Survey (RATS)

The Rapid Temporal Survey (RATS) explores the faint, variable sky. Our observations search a parameter space which, until now, has never been exploited from the ground. Our strategy involves observing the sky close to the Galactic plane with wide-field CCD cameras. An exposure is obtained approximately every minute with the total observation of each field lasting around 2 hours. In this paper we present the first 6 epochs of observations which were taken over 5 years from 2003--2008 and cover over 31 square degrees of which 16.2 is within 10° of the Galactic plane. The number of stars contained in these data is over 3.0 x10^6. We have developed a method of combining the output of two variability tests in order to detect variability on time-scales ranging from a few minutes to a few hours. Using this technique we find 1.2 x 10^5 variables -- equal to 4.1 per cent of stars in our data. Follow-up spectroscopic observations have allowed us to identify the nature of a fraction of these sources. These include a pulsating white dwarf which appears to have a hot companion, a number of stars with A-type spectra that vary on a period in the range 20--35 min. Our primary goal is the discovery of new AM CVn systems: we find 66 sources which appear to show periodic modulation on a time-scales less than 40 min and a colour consistent with the known AM CVn systems. Of those sources for which we have spectra of, none appears to be an AM CVn system, although we have 12 candidate AM CVn systems with periods less than 25 min for which spectra are still required. Although our numbers are not strongly constraining, they are consistent with the predictions of Nelemans et al.

preprint2011arXiv

The formation of the eccentric-orbit millisecond pulsar J1903+0327 and the origin of single millisecond pulsars

The millisecond pulsar J1903+0327 is accompanied by an ordinary G-dwarf star in an unusually wide ($P_{\rm orb} \simeq 95.2$\,days) and eccentric ($e \simeq 0.44$) orbit. The standard model for producing MSPs fails to explain the orbital characteristics of this extraordinary binary, and alternative binary models are unable to explain the observables. We present a triple-star model for producing MSPs in relatively wide eccentric binaries with a normal (main-sequence) stellar companion. We start from a stable triple system consisting of a Low-Mass X-ray Binary (LMXB) with an orbital period of at least 1 day, accompanied by a G-dwarf in a wide and possibly eccentric orbit. Variations in the initial conditions naturally provide a satisfactory explanation for the unexplained triple component in the eclipsing soft X-ray transient 4U~2129+47 or the cataclysmic variable EC 19314-5915. The best explanation for J1903, however, results from the expansion of the orbit of the LMXB, driven by the mass transfer from the evolving donor star to its neutron star companion, which causes the triple eventually to becomes dynamically unstable. Using numerical computations we show that, depending on the precise system configuration at the moment the triple becomes dynamically unstable, the ejection of each of the three components is possible. If the donor star of the LMXB is ejected, a system resembling J1903, will result. If the neutron star is ejected, a single MSP results. This model therefore also provides a straightforward mechanism for forming single MSP in the Galactic disk. We conclude that the Galaxy contains some 30--300 binaries with characteristics similar to J1903, and about an order of magnitude fewer single millisecond pulsars produced with the proposed triple scenario.

preprint2009arXiv

A census of AM CVn stars: three new candidates and one confirmed 48.3-minute binary

We present three new candidate AM CVn binaries, plus one confirmed new system, from a spectroscopic survey of color-selected objects from the Sloan Digital Sky Survey. All four systems were found from their helium emission lines in low-resolution spectra taken on the Hale telescope at Palomar, and the Nordic Optical Telescope and the William Herschel Telescope on La Palma. The ultra-compact binary nature of SDSS J090221.35+381941.9 was confirmed using phase-resolved spectroscopy at the Keck-I telescope. From the characteristic radial velocity `S-wave' observed in the helium emission lines we measure an orbital period of 48.31 +/- 0.08 min. The continuum emission can be described with a blackbody or a helium white dwarf atmosphere of T_eff ~ 15,000K, in agreement with theoretical cooling models for relatively massive accretors and/or donors. The absence in the spectrum of broad helium absorption lines from the accreting white dwarf suggests that the accreting white dwarf cannot be much hotter than 15,000K, or that an additional component such as the accretion disk contributes substantially to the optical flux. Two of the candidate systems, SDSS J152509.57+360054.5 and SDSS J172102.48+273301.2, do show helium absorption in the blue part of their spectra in addition to the characteristic helium emission lines. This, in combination with the high effective temperatures of ~18,000K and ~16,000K suggests both two be at orbital periods below ~40min. The third candidate, SDSS J164228.06+193410.0, exhibits remarkably strong helium emission on top of a relatively cool (T_eff~12,000K) continuum, indicating an orbital period above ~50min.

preprint2009arXiv

Impact of helium diffusion and helium-flash-induced carbon production on gravity-mode pulsations in subdwarf B stars

Realistic stellar models are essential to the forward modelling approach in asteroseismology. For practicality however, certain model assumptions are also required. For example, in the case of subdwarf B stars, one usually starts with zero-age horizontal branch structures without following the progenitor evolution. We analyse the effects of common assumptions in subdwarf B models on the g-mode pulsational properties. We investigate if and how the pulsation periods are affected by the H-profile in the core-envelope transition zone. Furthermore, the effects of C-production and convective mixing during the core helium flash are evaluated. Finally, we reanalyse the effects of stellar opacities on the mode excitation in subdwarf B stars. We find that helium settling causes a shift in the theoretical blue edge of the g-mode instability domain to higher effective temperatures. This results in a closer match to the observed instability strip of long-period sdB pulsators, particularly for l<=3 modes. We show further that the g-mode spectrum is extremely sensitive to the H-profile in the core-envelope transition zone. If atomic diffusion is efficient, details of the initial shape of the profile become less important in the course of evolution. Diffusion broadens the chemical gradients, and results in less effective mode trapping and different pulsation periods. Furthermore, we report on the possible consequences of the He-flash for the g-modes. The outer edge of a flash-induced convective region introduces an additional chemical transition in the stellar models, and the corresponding spike in the Brunt-Vaisala frequency produces a complicated mode trapping signature in the period spacings.

preprint2009arXiv

RAT J1953+1859: a dwarf nova discovered through high amplitude QPOs in quiescence

We report the discovery of an accreting binary, RAT J1953+1859, made during the RApid Temporal Survey (RATS) on the Isaac Newton Telescope. It showed high amplitude (0.3 mag) quasi-periodic oscillations on a timescale of ~20 mins. Further observations made using the Nordic Optical Telescope showed it to be ~4 mag brighter than in the discovery images. These photometric observations, together with radial velocity data taken using the William Herschel Telescope, point to an orbital period of ~90 mins. These data suggest that RAT J1953+1859 is a dwarf novae of the SU UMa type. What makes RAT J1953+1859 unusual is that it is the first such system to be discovered as a result of high amplitude QPOs during quiescence. This suggests that high-cadence wide-field surveys could be another means to discover cataclysmic variables as a result of their short period variability.

preprint2009arXiv

The influence of short term variations in AM CVn systems on LISA measurements

We study the effect of short term variations of the evolution of AM CVn systems on their gravitational wave emissions and in particular LISA observations. We model the systems according to their equilibrium mass-transfer evolution as driven by gravitational wave emission and tidal interaction, and determine their reaction to a sudden perturbation of the system. This is inspired by the suggestion to explain the orbital period evolution of the ultra-compact binary systems V407 Vul and RX-J0806+1527 by non-equilibrium mass transfer. The characteristics of the emitted gravitational wave signal are deduced from a Taylor expansion of a Newtonian quadrupolar emission model, and the changes in signal structure as visible to the LISA mission are determined. We show that short term variations can significantly change the higher order terms in the expansion, and thus lead to spurious (non) detection of frequency derivatives. This may hamper the estimation of the parameters of the system, in particular their masses and distances. However, we find that overall detection is still secured as signals still can be described by general templates. We conclude that a better modelling of the effects of short term variations is needed to prepare the community for astrophysical evaluations of real gravitational wave data of AM CVn systems.

preprint2005arXiv

Simulation of the White Dwarf -- White Dwarf galactic background in the LISA data

LISA (Laser Interferometer Space Antenna) is a proposed space mission, which will use coherent laser beams exchanged between three remote spacecraft to detect and study low-frequency cosmic gravitational radiation. In the low-part of its frequency band, the LISA strain sensitivity will be dominated by the incoherent superposition of hundreds of millions of gravitational wave signals radiated by inspiraling white-dwarf binaries present in our own galaxy. In order to estimate the magnitude of the LISA response to this background, we have simulated a synthesized population that recently appeared in the literature. We find the amplitude of the galactic white-dwarf binary background in the LISA data to be modulated in time, reaching a minimum equal to about twice that of the LISA noise for a period of about two months around the time when the Sun-LISA direction is roughly oriented towards the Autumn equinox. Since the galactic white-dwarfs background will be observed by LISA not as a stationary but rather as a cyclostationary random process with a period of one year, we summarize the theory of cyclostationary random processes and present the corresponding generalized spectral method needed to characterize such process. We find that, by measuring the generalized spectral components of the white-dwarf background, LISA will be able to infer properties of the distribution of the white-dwarfs binary systems present in our Galaxy.

preprint2005arXiv

The White Dwarf -- White Dwarf galactic background in the LISA data

LISA (Laser Interferometer Space Antenna) is a proposed space mission, which will use coherent laser beams exchanged between three remote spacecraft to detect and study low-frequency cosmic gravitational radiation. In the low-part of its frequency band, the LISA strain sensitivity will be dominated by the incoherent superposition of hundreds of millions of gravitational wave signals radiated by inspiraling white-dwarf binaries present in our own galaxy. In order to estimate the magnitude of the LISA response to this background, we have simulated a synthesized population that recently appeared in the literature. We find the amplitude of the galactic white-dwarf binary background in the LISA data to be modulated in time, reaching a minimum equal to about twice that of the LISA noise for a period of about two months around the time when the Sun-LISA direction is roughly oriented towards the Autumn equinox. Since the galactic white-dwarfs background will be observed by LISA not as a stationary but rather as a cyclostationary random process with a period of one year, we summarize the theory of cyclostationary random processes, present the corresponding generalized spectral method needed to characterize such process, and make a comparison between our analytic results and those obtained by applying our method to the simulated data. We find that, by measuring the generalized spectral components of the white-dwarf background, LISA will be able to infer properties of the distribution of the white-dwarfs binary systems present in our Galaxy.

preprint2000arXiv

Population synthesis for double white dwarfs I.Close detached systems

We model the population of double white dwarfs in the Galaxy and find a better agreement with observations compared to earlier studies, due to two modifications. The first is the treatment of the first phase of unstable mass transfer and the second the modelling of the cooling of the white dwarfs. A satisfactory agreement with observations of the local sample of white dwarfs is achieved if we assume that the initial binary fraction is ~ 50% and that the lowest mass white dwarfs (M < 0.3 Msun) cool faster than the most recently published cooling models predict. With this model we find a Galactic birth rate of close double white dwarfs of 0.05 yr^{-1}, a birth rate of AM CVn systems of 0.005 yr^{-1}, a merger rate of pairs with a combined mass exceeding the Chandrasekhar limit (which may be progenitors of SNe Ia) of 0.003 yr^{-1} and a formation rate of planetary nebulae of 1 yr^{-1}. We estimate the total number of double white dwarfs in the Galaxy as 2.5 10^8. In an observable sample with a limiting magnitude V_lim = 15 we predict the presence of ~855 white dwarfs of which ~220 are close pairs. Of these 10 are double CO white dwarfs of which one has a combined mass exceeding the Chandrasekhar limit and will merge within a Hubble time.

preprint2000arXiv

The population of close double white dwarfs in the Galaxy

We present a new model for the Galactic population of close double white dwarfs. The model accounts for the suggestion of the avoidance of a substantial spiral-in during mass transfer between a giant and a main-sequence star of comparable mass and for detailed cooling models. It agrees well with the observations of the local sample of white dwarfs if the initial binary fraction is close to 50% and an ad hoc assumption is made that white dwarfs with mass less than about 0.3 solar mass cool faster than the models suggest. About 1000 white dwarfs brighter than V=15 have to be surveyed for detection of a pair which has total mass greater than the Chandrasekhar mass and will merge within 10 Gyr.