Source author record

S. P. Littlefair

S. P. Littlefair 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

81works
7topics
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

81 published item(s)

preprint2022arXiv

A 62-minute orbital period black widow binary in a wide hierarchical triple

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution.

preprint2022arXiv

A tentative 114-minute orbital period challenges the ultra-compact nature of the X-ray binary 4U 1812-12

We present a detailed time-resolved photometric study of the ultra-compact X-ray binary candidate 4U 1812-12. The multicolor light curves obtained with HiPERCAM on the 10.4-m Gran Telescopio Canarias show an aprox 114 min modulation similar to a superhump. Under this interpretation, this period should lie very close to the orbital period of the system. Contrary to what its other observational properties suggest (namely, persistent dim luminosity, low optical-to-X-ray flux ratio and lack of hydrogen features in the optical spectrum), this implies that 4U1812-12 is most likely not an ultra-compact X-ray binary, which are usually defined as systems with orbital periods lower than 80 min. We discuss the nature of the system, showing that a scenario in which 4U 1812-12 is the progenitor of an ultra-compact X-ray binary may reconcile all the observables.

preprint2022arXiv

First detection of the outer edge of an AGN accretion disc: Very fast multiband optical variability of NGC 4395 with GTC/HiPERCAM and LT/IO:O

We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with disc reprocessing but not with reprocessing in the Broad Line Region (BLR). There is very little increase in lag with wavelength at long wavelengths, indicating an outer edge (Rout) to the reprocessor. We have compared truncated disc reprocessing models to the combined HiPERCAM and previous X-ray/UV lags. For the normally accepted mass of 3.6E5 solar, we obtain reasonable agreement with zero spin, Rout ~1700 Rg, and the DONE physically-motivated temperature-dependent disc colour correction factor (fcol). A smaller mass of 4E4 solar can only be accomodated if fcol=2.4, which is probably unrealistically high. Disc self gravity is probably unimportant in this low mass AGN but an obscuring wind may provide an edge. For the small mass the dust sublimation radius is similar to Rout, so the wind could be dusty. However for the more likely large mass the sublimation radius is further out so the optically-thick base of a line-driven gaseous wind is more likely. The inner edge of the BLR is close to Rout in both cases. These observations provide the first good evidence for a truncated AGN disc and caution that truncation should be included in reverberation lag modelling.

preprint2022arXiv

Multi-colour optical light curves of the companion star to the millisecond pulsar PSR J2051-0827

We present simultaneous, multi-colour optical light curves of the companion star to the black-widow pulsar PSR J2051-0827, obtained approximately 10 years apart using ULTRACAM and HiPERCAM, respectively. The ULTRACAM light curves confirm the previously reported asymmetry in which the leading hemisphere of the companion star appears to be brighter than the trailing hemisphere. The HiPERCAM light curves, however, do not show this asymmetry, demonstrating that whatever mechanism is responsible for it varies on timescales of a decade or less. We fit the symmetrical HiPERCAM light curves with a direct-heating model to derive the system parameters, finding an orbital inclination of $55.9^{+4.8}_{-4.1}$ degrees, in good agreement with radio-eclipse constraints. We find that approximately half of the pulsar's spin-down energy is converted to optical luminosity, resulting in temperatures ranging from approximately $5150^{+190}_{-190}$ K on the day side to $2750^{+130}_{-150}$ K on the night side of the companion star. The companion star is close to filling its Roche lobe ($f_{\rm RL} =0.88^{+0.02}_{-0.02}$) and has a mass of $0.039^{+0.010}_{-0.011}$ M$_{\odot}$, giving a mean density of $20.24^{+0.59}_{-0.44}$ g cm$^{-3}$ and an apsidal motion constant in the range $0.0036 < k_2 < 0.0047$. The companion mass and mean density values are consistent with those of brown dwarfs, but the apsidal motion constant implies a significantly more centrally-condensed internal structure than is typical for such objects.

preprint2022arXiv

Probing for the host galaxies of the fast X-ray transients XRT 000519 and XRT 110103

Over the past few years, $\sim$30 extragalactic fast X-ray transients (FXRTs) have been discovered, mainly in Chandra and XMM-Newton data. Their nature remains unclear, with proposed origins including a double neutron star merger, a tidal disruption event involving an intermediate-mass black hole and a white dwarf, or a supernova shock breakout. A decisive differentiation between these three promising mechanisms for their origin requires an understanding of the FXRT energetics, environments, and/or host properties. We present optical observations obtained with the Very Large Telescope for the FXRTs XRT 000519 and XRT 110103 and Gran Telescopio Canarias observations for XRT 000519 designed to search for host galaxies of these FXRTs. In the $g_s$, $r_s$ and $R$-band images, we detect an extended source on the North-West side of the $\sim$ $1^{\prime\prime}$ (68% confidence) error circle of the X-ray position of XRT 000519 with a Kron magnitude of $g_s=$26.29$\pm$0.09 (AB magnitude). We discuss the XRT 000519 association with the probable host candidate for various possible distances, and we conclude that if XRT 000519 is associated with the host candidate a supernova shock breakout scenario is likely excluded. No host galaxy is found near XRT 110103 down to a limiting magnitude of $R>25.8$.

preprint2021arXiv

Magnetic white dwarfs in post-common-envelope binaries

Magnitude-limited samples have shown that 20-25 per cent of cataclysmic variables contain white dwarfs with magnetic fields of Mega Gauss strength, in stark contrast to the approximately 5 per cent of single white dwarfs with similar magnetic field strengths. Moreover, the lack of identifiable progenitor systems for magnetic cataclysmic variables leads to considerable challenges when trying to understand how these systems form and evolve. Here we present a sample of six magnetic white dwarfs in detached binaries with low-mass stellar companions where we have constrained the stellar and binary parameters including, for the first time, reliable mass estimates for these magnetic white dwarfs. We find that they are systematically more massive than non-magnetic white dwarfs in detached binaries. These magnetic white dwarfs generally have cooling ages of more than 1 Gyr and reside in systems that are very close to Roche-lobe filling. Our findings are more consistent with these systems being temporarily detached cataclysmic variables, rather than pre-cataclysmic binaries, but we cannot rule out the latter possibility. We find that these systems can display unusual asymmetric light curves that may offer a way to identify them in larger numbers in future. Seven new candidate magnetic white dwarf systems are also presented, three of which have asymmetric light curves. Finally, we note that several newly identified magnetic systems have archival spectra where there is no clear evidence of magnetism, meaning that these binaries have been previously missed. Nevertheless, there remains a clear lack of younger detached magnetic white dwarf systems.

preprint2020arXiv

A new class of Roche lobe-filling hot subdwarf binaries

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD.

preprint2020arXiv

Spectroscopic and Photometric Periods of Six Ultracompact Accreting Binaries

Ultracompact accreting binary systems each consist of a stellar remnant accreting helium-enriched material from a compact donor star. Such binaries include two related sub-classes, AM CVn-type binaries and helium cataclysmic variables, in both of which the central star is a white dwarf. We present a spectroscopic and photometric study of six accreting binaries with orbital periods in the range of 40--70 min, including phase-resolved VLT spectroscopy and high-speed ULTRACAM photometry. Four of these are AM CVn systems and two are helium cataclysmic variables. For four of these binaries we are able to identify orbital periods (of which three are spectroscopic). SDSS J1505+0659 has an orbital period of 67.8 min, significantly longer than previously believed, and longer than any other known AM CVn binary. We identify a WISE infrared excess in SDSS J1505+0659 that we believe to be the first direct detection of an AM CVn donor star in a non-direct impacting binary. The mass ratio of SDSS J1505+0659 is consistent with a white dwarf donor. CRTS J1028-0819 has an orbital period of 52.1 min, the shortest period of any helium cataclysmic variable. MOA 2010-BLG-087 is co-aligned with a K-class star that dominates its spectrum. ASASSN-14ei and ASASSN-14mv both show a remarkable number of echo outbursts following superoutbursts (13 and 10 echo outbursts respectively). ASASSN-14ei shows an increased outburst rate over the years following its superoutburst, perhaps resulting from an increased accretion rate.

preprint2020arXiv

The first ultracompact Roche lobe-filling hot subdwarf binary

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star.

preprint2020arXiv

WD1032+011, an inflated brown dwarf in an old eclipsing binary with a white dwarf

We present the discovery of only the third brown dwarf known to eclipse a non-accreting white dwarf. Gaia parallax information and multi-colour photometry confirm that the white dwarf is cool (9950$\pm$150K) and has a low mass (0.45$\pm$0.05~MSun), and spectra and lightcurves suggest the brown dwarf has a mass of 0.067 $\pm$0.006 MSun (70 MJup) and a spectral type of L5 $\pm$1. The kinematics of the system show that the binary is likely to be a member of the thick disk and therefore at least 5 Gyr old. The high cadence lightcurves show that the brown dwarf is inflated, making it the first brown dwarf in an eclipsing white dwarf-brown dwarf binary to be so.

preprint2019arXiv

Characterising the i-band variability of YSOs over six orders of magnitude in timescale

We present an $i$-band photometric study of over 800 young stellar objects in the OB association Cep OB3b, which samples timescales from 1 minute to ten years. Using structure functions we show that on all timescales ($τ$) there is a monotonic decrease in variability from Class I to Class II through the transition disc (TD) systems to Class III, i.e. the more evolved systems are less variable. The Class Is show an approximately power-law increase ($τ^{0.8}$) in variability from timescales of a few minutes to ten years. The Class II, TDs and Class III systems show a qualitatively different behaviour with most showing a power-law increase in variability up to a timescale corresponding to the rotational period of the star, with little additional variability beyond that timescale. However, about a third of the Class IIs show lower overall variability, but their variability is still increasing at 10 years. This behaviour can be explained if all Class IIs have two primary components to their variability. The first is an underlying roughly power-law variability spectrum, which evidence from the infrared suggests is driven by accretion rate changes. The second component is approximately sinusoidal and results from the rotation of the star. We suggest that the systems with dominant longer-timescale variability have a smaller rotational modulation either because they are seen at low inclinations or have more complex magnetic field geometries. We derive a new way of calculating structure functions for large simulated datasets (the "fast structure function"), based on fast Fourier transforms.

preprint2019arXiv

Full orbital solution for the binary system in the northern Galactic disc microlensing event Gaia16aye

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

preprint2016arXiv

A LOFAR mini-survey for low-frequency radio emission from the nearest brown dwarfs

We have conducted a mini-survey for low-frequency radio emission from some of the closest brown dwarfs to the Sun with rapid rotation rates: SIMP J013656.5+093347, WISEPC J150649.97+702736.0, and WISEPA J174124.26+255319.5. We have placed robust 3-sigma upper limits on the flux density in the 111 - 169 MHz frequency range for these targets: WISE 1506: < 0.72 mJy; WISE 1741: < 0.87 mJy; SIMP 0136: < 0.66 mJy. At 8 hours of integration per target to achieve these limits, we find that systematic and detailed study of this class of object at LOFAR frequencies will require a substantial dedication of resources.

preprint2016arXiv

A radio pulsing white dwarf binary star

White dwarfs are compact stars, similar in size to Earth but ~200,000 times more massive. Isolated white dwarfs emit most of their power from ultraviolet to near-infrared wavelengths, but when in close orbits with less dense stars, white dwarfs can strip material from their companions, and the resulting mass transfer can generate atomic line and X-ray emission, as well as near- and mid-infrared radiation if the white dwarf is magnetic. However, even in binaries, white dwarfs are rarely detected at far-infrared or radio frequencies. Here we report the discovery of a white dwarf / cool star binary that emits from X-ray to radio wavelengths. The star, AR Scorpii (henceforth AR Sco), was classified in the early 1970s as a delta-Scuti star, a common variety of periodic variable star. Our observations reveal instead a 3.56 hr period close binary, pulsing in brightness on a period of 1.97 min. The pulses are so intense that AR Sco's optical flux can increase by a factor of four within 30 s, and they are detectable at radio frequencies, the first such detection for any white dwarf system. They reflect the spin of a magnetic white dwarf which we find to be slowing down on a 10^7 yr timescale. The spin-down power is an order of magnitude larger than that seen in electromagnetic radiation, which, together with an absence of obvious signs of accretion, suggests that AR Sco is primarily spin-powered. Although the pulsations are driven by the white dwarf's spin, they originate in large part from the cool star. AR Sco's broad-band spectrum is characteristic of synchrotron radiation, requiring relativistic electrons. These must either originate from near the white dwarf or be generated in situ at the M star through direct interaction with the white dwarf's magnetosphere.

preprint2016arXiv

Flares, wind and nebulae: the 2015 December mini-outburst of V404 Cygni

After more than 26 years in quiescence, the black hole transient V404 Cyg went into a luminous outburst in June 2015, and additional activity was detected in late December of the same year. Here, we present an optical spectroscopic follow-up of the December mini-outburst, together with X-ray, optical and radio monitoring that spanned more than a month. Strong flares with gradually increasing intensity are detected in the three spectral ranges during the ~10 days following the Swift trigger. Our optical spectra reveal the presence of a fast outflowing wind, as implied by the detection of a P-Cyg profile (He I - 5876 A) with a terminal velocity of ~2500 km/s. Nebular-like spectra -- with an H_alpha equivalent width of ~500 A -- are also observed. All these features are similar to those seen during the main June 2015 outburst. Thus, the fast optical wind simultaneous with the radio jet is most likely present in every V404 Cyg outburst. Finally, we report on the detection of a strong radio flare in late January 2016, when X-ray and optical monitoring had stopped due to Sun constraints.

preprint2016arXiv

Furiously Fast and Red: Sub-second Optical Flaring in V404 Cyg during the 2015 Outburst Peak

We present observations of rapid (sub-second) optical flux variability in V404 Cyg during its 2015 June outburst. Simultaneous three-band observations with the ULTRACAM fast imager on four nights show steep power spectra dominated by slow variations on ~100-1000s timescales. Near the peak of the outburst on June 26, a dramatic change occurs and additional, persistent sub-second optical flaring appears close in time to giant radio and X-ray flaring. The flares reach peak optical luminosities of ~few x 10^{36} erg/s. Some are unresolved down to a time resolution of 24 milliseconds. Whereas the fast flares are stronger in the red, the slow variations are bluer when brighter. The redder slopes, emitted power, and characteristic timescales of the fast flares can be explained as optically-thin synchrotron emission from a compact jet arising on size scales ~140-500 Gravitational radii (with a possible additional contribution by a thermal particle distribution). The origin of the slower variations is unclear. The optical continuum spectral slopes are strongly affected by dereddening uncertainties and contamination by strong H{alpha} emission, but the variations of these slopes follow relatively stable loci as a function of flux. Cross-correlating the slow variations between the different bands shows asymmetries on all nights consistent with a small red skew (i.e., red lag). X-ray reprocessing and non-thermal emission could both contribute to these. These data reveal a complex mix of components over five decades in timescale during the outburst.

preprint2016arXiv

High-speed photometry of the disintegrating planetesimals at WD1145+017: evidence for rapid dynamical evolution

We obtained high-speed photometry of the disintegrating planetesimals orbiting the white dwarf WD1145+017, spanning a period of four weeks. The light curves show a dramatic evolution of the system since the first observations obtained about seven months ago. Multiple transit events are detected in every light curve, which have varying durations(~3-12min) and depths (~10-60%). The time-averaged extinction is ~11%, much higher than at the time of the Kepler observations. The shortest-duration transits require that the occulting cloud of debris has a few times the size of the white dwarf, longer events are often resolved into the superposition of several individual transits. The transits evolve on time scales of days, both in shape and in depth, with most of them gradually appearing and disappearing over the course of the observing campaign. Several transits can be tracked across multiple nights, all of them recur on periods of ~4.49h, indicating multiple planetary debris fragments on nearly identical orbits. Identifying the specific origin of these bodies within this planetary system, and the evolution leading to their current orbits remains a challenging problem.

preprint2016arXiv

HiPERCAM: A high-speed, quintuple-beam CCD camera for the study of rapid variability in the Universe

HiPERCAM is a high-speed camera for the study of rapid variability in the Universe. The project is funded by a 3.5MEuro European Research Council Advanced Grant. HiPERCAM builds on the success of our previous instrument, ULTRACAM, with very significant improvements in performance thanks to the use of the latest technologies. HiPERCAM will use 4 dichroic beamsplitters to image simultaneously in 5 optical channels covering the u'g'r'i'z' bands. Frame rates of over 1000 per second will be achievable using an ESO CCD controller (NGC), with every frame GPS timestamped. The detectors are custom-made, frame-transfer CCDs from e2v, with 4 low-noise (2.5e-) outputs, mounted in small thermoelectrically-cooled heads operated at 180 K, resulting in virtually no dark current. The two reddest CCDs will be deep-depletion devices with anti-etaloning, providing high quantum efficiencies across the red part of the spectrum with no fringing. The instrument will also incorporate scintillation noise correction via the conjugate-plane photometry technique. The opto-mechanical chassis will make use of additive manufacturing techniques in metal to make a light-weight, rigid and temperature-invariant structure. First light is expected on the 4.2m William Herschel Telescope on La Palma in 2017 (on which the field of view will be 10' with a 0.3"/pixel scale), with subsequent use planned on the 10.4m Gran Telescopio Canarias on La Palma (on which the field of view will be 4' with a 0.11"/pixel scale) and the 3.5m New Technology Telescope in Chile.

preprint2016arXiv

Long-term eclipse timing of white dwarf binaries: an observational hint of a magnetic mechanism at work

We present a long-term programme for timing the eclipses of white dwarfs in close binaries to measure apparent and/or real variations in their orbital periods. Our programme includes 67 close binaries, both detached and semi-detached and with M-dwarfs, K-dwarfs, brown dwarfs or white dwarfs secondaries. In total, we have observed more than 650 white dwarf eclipses. We use this sample to search for orbital period variations and aim to identify the underlying cause of these variations. We find that the probability of observing orbital period variations increases significantly with the observational baseline. In particular, all binaries with baselines exceeding 10 yrs, with secondaries of spectral type K2 -- M5.5, show variations in the eclipse arrival times that in most cases amount to several minutes. In addition, among those with baselines shorter than 10 yrs, binaries with late spectral type (>M6), brown dwarf or white dwarf secondaries appear to show no orbital period variations. This is in agreement with the so-called Applegate mechanism, which proposes that magnetic cycles in the secondary stars can drive variability in the binary orbits. We also present new eclipse times of NN Ser, which are still compatible with the previously published circumbinary planetary system model, although only with the addition of a quadratic term to the ephemeris. Finally, we conclude that we are limited by the relatively short observational baseline for many of the binaries in the eclipse timing programme, and therefore cannot yet draw robust conclusions about the cause of orbital period variations in evolved, white dwarf binaries.

preprint2016arXiv

Long-Term Photometry of IC 348 with the YETI Network

We present long-term photometric observations of the young open cluster IC 348 with a baseline time-scale of 2.4 yr. Our study was conducted with several telescopes from the Young Exoplanet Transit Initiative (YETI) network in the Bessel $R$ band to find periodic variability of young stars. We identified 87 stars in IC 348 to be periodically variable; 33 of them were unreported before. Additionally, we detected 61 periodic non-members of which 41 are new discoveries. Our wide field of view was the key to those numerous newly found variable stars. The distribution of rotation periods in IC 348 has always been of special interest. We investigate it further with our newly detected periods but we cannot find a statistically significant bimodality. We also report the detection of a close eclipsing binary in IC 348 composed of a low-mass stellar component ($M \gtrsim 0.09\,\mathrm{M}_{\odot}$) and a K0 pre-main sequence star ($M \approx 2.7\,\mathrm{M}_{\odot}$). Furthermore, we discovered three detached binaries among the background stars in our field of view and confirmed the period of a fourth one.

preprint2016arXiv

Search for transiting exoplanets and variable stars in the open cluster NGC 7243

We report results of the first five observing campaigns for the open stellar cluster NGC 7243 in the frame of project Young Exoplanet Transit Initiative (YETI). The project focuses on the monitoring of young and nearby stellar clusters, with the aim to detect young transiting exoplanets, and to study other variability phenomena on time-scales from minutes to years. After five observing campaigns and additional observations during 2013 and 2014, a clear and repeating transit-like signal was detected in the light curve of J221550.6+495611. Furthermore, we detected and analysed 37 new eclipsing binary stars in the studied region. The best fit parameters and light curves of all systems are given. Finally, we detected and analysed 26 new, presumably pulsating variable stars in the studied region. The follow-up investigation of these objects, including spectroscopic measurements of the exoplanet candidate, is currently planned.

preprint2016arXiv

The crowded magnetosphere of the post common envelope binary QS Virginis

We present high speed photometry and high resolution spectroscopy of the eclipsing post common envelope binary QS Virginis (QS Vir). Our UVES spectra span multiple orbits over more than a year and reveal the presence of several large prominences passing in front of both the M star and its white dwarf companion, allowing us to triangulate their positions. Despite showing small variations on a timescale of days, they persist for more than a year and may last decades. One large prominence extends almost three stellar radii from the M star. Roche tomography reveals that the M star is heavily spotted and that these spots are long-lived and in relatively fixed locations, preferentially found on the hemisphere facing the white dwarf. We also determine precise binary and physical parameters for the system. We find that the 14,220 +/- 350K white dwarf is relatively massive, 0.782 +/- 0.013Ms, and has a radius of 0.01068 +/- 0.00007Rs, consistent with evolutionary models. The tidally distorted M star has a mass of 0.382 +/- 0.006Ms and a radius of 0.381 +/- 0.003Rs, also consistent with evolutionary models. We find that the magnesium absorption line from the white dwarf is broader than expected. This could be due to rotation (implying a spin period of only ~700 seconds), or due to a weak (~100kG) magnetic field, we favour the latter interpretation. Since the M star's radius is still within its Roche lobe and there is no evidence that its over-inflated we conclude that QS Vir is most likely a pre-cataclysmic binary just about to become semi-detached.

preprint2016arXiv

Transmission spectroscopy of the inflated exoplanet WASP-52b, and evidence for a bright region on the stellar surface

We have measured the transmission spectrum of the extremely inflated hot Jupiter WASP-52b using simultaneous photometric observations in SDSS u', g' and a filter centred on the sodium doublet (NaI) with the ULTRACAM instrument mounted on the 4.2m William Herschel Telescope. We find that Rayleigh scattering is not the dominant source of opacity within the planetary atmosphere and find a transmission spectrum more consistent with wavelength-independent opacity such as from clouds. We detect an in-transit anomaly that we attribute to the presence of stellar activity and find that this feature can be more simply modelled as a bright region on the stellar surface akin to Solar faculae rather than spots. A spot model requires a significantly larger planet/star radius ratio than that found in previous studies. Our results highlight the precision that can be achieved by ground-based photometry with errors in the scaled planetary radii of less than one atmospheric scale height, comparable to HST observations.

preprint2016arXiv

Using Gaussian processes to model light curves in the presence of flickering: the eclipsing cataclysmic variable ASASSN-14ag

The majority of cataclysmic variable (CV) stars contain a stochastic noise component in their light curves, commonly referred to as flickering. This can significantly affect the morphology of CV eclipses and increases the difficulty in obtaining accurate system parameters with reliable errors through eclipse modelling. Here we introduce a new approach to eclipse modelling, which models CV flickering with the help of Gaussian processes (GPs). A parameterised eclipse model - with an additional GP component - is simultaneously fit to 8 eclipses of the dwarf nova ASASSN-14ag and system parameters determined. We obtain a mass ratio $q$ = 0.149 $\pm$ 0.016 and inclination $i$ = 83.4 $^{+0.9}_{-0.6}$ $^{\circ}$. The white dwarf and donor masses were found to be $M_{w}$ = 0.63 $\pm$ 0.04 $M_{\odot}$ and $M_{d}$ = 0.093 $^{+0.015}_{-0.012}$ $M_{\odot}$, respectively. A white dwarf temperature $T_{w}$ = 14000 $^{+2200}_{-2000}$ K and distance $d$ = 146 $^{+24}_{-20}$ pc were determined through multicolour photometry. We find GPs to be an effective way of modelling flickering in CV light curves and plan to use this new eclipse modelling approach going forward.

preprint2015arXiv

A large, long-lived structure near the trojan L5 point in the post common-envelope binary SDSS J1021+1744

SDSS J1021+1744 is a detached, eclipsing white dwarf / M dwarf binary discovered in the Sloan Digital Sky Survey. Outside the primary eclipse, the light curves of such systems are usually smooth and characterised by low-level variations caused by tidal distortion and heating of the M star component. Early data on SDSS J1021+1744 obtained in June 2012 was unusual in showing a dip in flux of uncertain origin shortly after the white dwarf's eclipse. Here we present high-time resolution, multi-wavelength observations of 35 more eclipses over 1.3 years, showing that the dip has a lifetime extending over many orbits. Moreover the "dip" is in fact a series of dips that vary in depth, number and position, although they are always placed in the phase interval 1.06 to 1.26 after the white dwarf's eclipse, near the L5 point in this system. Since SDSS J1021+1744 is a detached binary, it follows that the dips are caused by the transit of the white dwarf by material around the Lagrangian L5 point. A possible interpretation is that they are the signatures of prominences, a phenomenon already known from H-alpha observations of rapidly rotating single stars as well as binaries. What makes SDSS J1021+1744 peculiar is that the material is dense enough to block continuum light. The dips appear to have finally faded out around 2015 May after the first detection by Parsons et al. in 2012, suggesting a lifetime of years.

preprint2015arXiv

Fourteen new eclipsing white dwarf plus main-sequence binaries from the SDSS and Catalina surveys

We report on the search for new eclipsing white dwarf plus main-sequence (WDMS) binaries in the light curves of the Catalina surveys. We use a colour selected list of almost 2000 candidate WDMS systems from the Sloan Digital Sky Survey, specifically designed to identify WDMS systems with cool white dwarfs and/or early M type main-sequence stars. We identify a total of 17 eclipsing systems, 14 of which are new discoveries. We also find 3 candidate eclipsing systems, 2 main-sequence eclipsing binaries and 22 non-eclipsing close binaries. Our newly discovered systems generally have optical fluxes dominated by the main-sequence components, which have earlier spectral types than the majority of previously discovered eclipsing systems. We find a large number of ellipsoidally variable binaries with similar periods, near 4 hours, and spectral types M2--3, which are very close to Roche-lobe filling. We also find that the fraction of eclipsing systems is lower than found in previous studies and likely reflects a lower close binary fraction among WDMS binaries with early M-type main-sequence stars due to their enhanced angular momentum loss compared to fully convective late M type stars, hence causing them to become cataclysmic variables quicker and disappear from the WDMS sample. Our systems bring the total number of known detached, eclipsing WDMS binaries to 71.

preprint2015arXiv

Magnetospherically driven optical and radio aurorae at the end of the stellar main sequence

Aurorae are detected from all the magnetized planets in our Solar System, including Earth. They are powered by magnetospheric current systems that lead to the precipitation of energetic electrons into the high-latitude regions of the upper atmosphere. In the case of the gas-giant planets, these aurorae include highly polarized radio emission at kilohertz and megahertz frequencies produced by the precipitating electrons, as well as continuum and line emission in the infrared, optical, ultraviolet and X-ray parts of the spectrum, associated with the collisional excitation and heating of the hydrogen-dominated atmosphere. Here we report simultaneous radio and optical spectroscopic observations of an object at the end of the stellar main sequence, located right at the boundary between stars and brown dwarfs, from which we have detected radio and optical auroral emissions both powered by magnetospheric currents. Whereas the magnetic activity of stars like our Sun is powered by processes that occur in their lower atmospheres, these aurorae are powered by processes originating much further out in the magnetosphere of the dwarf star that couple energy into the lower atmosphere. The dissipated power is at least four orders of magnitude larger than what is produced in the Jovian magnetosphere, revealing aurorae to be a potentially ubiquitous signature of large-scale magnetospheres that can scale to luminosities far greater than those observed in our Solar System. These magnetospheric current systems may also play a part in powering some of the weather phenomena reported on brown dwarfs.

preprint2015arXiv

PHL 1445: An eclipsing cataclysmic variable with a substellar donor near the period minimum

We present high-speed, three-colour photometry of the eclipsing dwarf nova PHL 1445, which, with an orbital period of 76.3 min, lies just below the period minimum of ~82 min for cataclysmic variable stars. Averaging four eclipses reveals resolved eclipses of the white dwarf and bright spot. We determined the system parameters by fitting a parameterised eclipse model to the averaged lightcurve. We obtain a mass ratio of q = 0.087 +- 0.006 and inclination i = 85.2 +- 0.9 degrees. The primary and donor masses were found to be Mw = 0.73 +- 0.03 Msun and Md = 0.064 +- 0.005 Msun, respectively. Through multicolour photometry a temperature of the white dwarf of Tw = 13200 +- 700 K and a distance of 220 +- 50 pc were determined. The evolutionary state of PHL 1445 is uncertain. We are able to rule out a significantly evolved donor, but not one that is slightly evolved. Formation with a brown dwarf donor is plausible; though the brown dwarf would need to be no older than 600 Myrs at the start of mass transfer, requiring an extremely low mass ratio (q = 0.025) progenitor system. PHL 1445 joins SDSS 1433 as a sub-period minimum CV with a substellar donor. These existence of two such systems raises an alternative possibility; that current estimates for the intrinsic scatter and/or position of the period minimum may be in error.

preprint2015arXiv

pt5m - a 0.5m robotic telescope on La Palma

pt5m is a 0.5m robotic telescope located on the roof of the 4.2m William Herschel Telescope (WHT) building, at the Roque de los Muchachos Observatory, La Palma. Using a 5-position filter wheel and CCD detector, and bespoke control software, pt5m provides a high quality robotic observing facility. The telescope first began robotic observing in 2012, and is now contributing to transient follow-up and time-resolved astronomical studies. In this paper we present the scientific motivation behind pt5m, as well as the specifications and unique features of the facility. We also present an example of the science we have performed with pt5m, where we measure the radius of the transiting exoplanet WASP-33b. We find a planetary radius of 1.603 +/- 0.014 R(J).

preprint2015arXiv

Rotation periods and astrometric motions of the Luhman 16AB brown dwarfs by high-resolution lucky-imaging monitoring

Context. Photometric monitoring of the variability of brown dwarfs can provide useful information about the structure of clouds in their cold atmospheres. The brown-dwarf binary system Luhman 16AB is an interesting target for such a study, as its components stand at the L/T transition and show high levels of variability. Luhman 16AB is also the third closest system to the Solar system, allowing precise astrometric investigations with ground-based facilities. Aims. The aim of the work is to estimate the rotation period and study the astrometric motion of both components. Methods. We have monitored Luhman 16AB over a period of two years with the lucky-imaging camera mounted on the Danish 1.54m telescope at La Silla, through a special i+z long-pass filter, which allowed us to clearly resolve the two brown dwarfs into single objects. An intense monitoring of the target was also performed over 16 nights, in which we observed a peak-to-peak variability of 0.20 \pm 0.02 mag and 0.34 \pm 0.02 mag for Luhman 16A and 16B, respectively. Results. We used the 16-night time-series data to estimate the rotation period of the two components. We found that Luhman 16B rotates with a period of 5.1 \pm 0.1 hr, in very good agreement with previous measurements. For Luhman 16A, we report that it rotates slower than its companion and, even though we were not able to get a robust determination, our data indicate a rotation period of roughly 8 hr. This implies that the rotation axes of the two components are well aligned and suggests a scenario in which the two objects underwent the same accretion process. The 2-year complete dataset was used to study the astrometric motion of Luhman 16AB. We predict a motion of the system that is not consistent with a previous estimate based on two months of monitoring, but cannot confirm or refute the presence of additional planetary-mass bodies in the system.

preprint2015arXiv

The First Millimeter Detection of a Non-Accreting Ultracool Dwarf

The well-studied M9 dwarf TVLM 513-46546 is a rapid rotator (P_rot ~ 2 hr) hosting a stable, dipolar magnetic field of ~3 kG surface strength. Here we report its detection with ALMA at 95 GHz at a mean flux density of $56 \pm 12$ uJy, making it the first ultracool dwarf detected in the millimeter band, excluding young, disk-bearing objects. We also report flux density measurements from unpublished archival VLA data and new optical monitoring data from the Liverpool Telescope. The ALMA data are consistent with a power-law radio spectrum that extends continuously between centimeter and millimeter wavelengths. We argue that the emission is due to the synchrotron process, excluding thermal, free-free, and electron cyclotron maser emission as possible sources. During the interval of the ALMA observation that phases with the maximum of the object's optical variability, the flux density is higher at a ~1.8 sigma significance level. These early results show how ALMA opens a new window for studying the magnetic activity of ultracool dwarfs, particularly shedding light on the particle acceleration mechanism operating in their immediate surroundings.

preprint2015arXiv

Total eclipse of the heart: The AM CVn Gaia14aae / ASSASN-14cn

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

preprint2014arXiv

A parameter study of the eclipsing CV in the Kepler field, KIS J192748.53+444724.5

We present high-speed, three-colour photometry of the eclipsing dwarf nova KIS J192748.53+444724.5 which is located in the Kepler field. Our data reveal sharp features corresponding to the eclipses of the accreting white dwarf followed by the bright spot where the gas stream joins the accretion disc. We determine the system parameters via a parameterized model of the eclipse fitted to the observed lightcurve. We obtain a mass ratio of q = 0.570 +/- 0.011 and an orbital inclination of 84.6 +/- 0.3 degrees. The primary mass is M_w = 0.69 +/- 0.07 Msun. The donor star's mass and radius are found to be M_d = 0.39 +/- 0.04 Msun and R_d = 0.43 +/- 0.01 Rsun, respectively. From the fluxes of the white dwarf eclipse we find a white dwarf temperature of T_w = 23000 +/- 3000 K, and a photometric distance to the system of 1600 +/- 200 pc, neglecting the effects of interstellar reddening. The white dwarf temperature in KISJ1927 implies the white dwarf is accreting at an average rate of Mdot = (1.4 +/- 0.8) x10e-9 Msun/yr, in agreement with estimates of the secular mass loss rate from the donor.

preprint2014arXiv

Defocussed Transmission Spectroscopy: A potential detection of sodium in the atmosphere of WASP-12b

We report on a pilot study of a novel observing technique, defocussed transmission spectroscopy, and its application to the study of exoplanet atmospheres using ground-based platforms. Similar to defocussed photometry, defocussed transmission spectroscopy has an added advantage over normal spectroscopy in that it reduces systematic errors due to flat-fielding, PSF variations, slit-jaw imperfections and other effects associated with ground-based observations. For one of the planetary systems studied, WASP-12b, we report a tentative detection of additional Na absorption of 0.12+/-0.03[+0.03]% during transit using a 2A wavelength mask. After consideration of a systematic that occurs mid-transit, it is likely that the true depth is actually closer to 0.15%. This is a similar level of absorption reported in the atmosphere of HD209458b (0.135+/-0.017%, Snellen et al. 2008). Finally, we outline methods that will improve the technique during future observations, based on our findings from this pilot study.

preprint2014arXiv

Discovery of ZZ Cetis in detached white dwarf plus main-sequence binaries

We present the first results of a dedicated search for pulsating white dwarfs (WDs) in detached white dwarf plus main-sequence binaries. Candidate systems were selected from a catalogue of WD+MS binaries, based on the surface gravities and effective temperatures of the WDs. We observed a total of 26 systems using ULTRACAM mounted on ESO's 3.5m New Technology Telescope (NTT) at La Silla. Our photometric observations reveal pulsations in seven WDs of our sample, including the first pulsating white dwarf with a main-sequence companion in a post common envelope binary, SDSSJ1136+0409. Asteroseismology of these new pulsating systems will provide crucial insight into how binary interactions, particularly the common envelope phase, affect the internal structure and evolution of WDs. In addition, our observations have revealed the partially eclipsing nature of one of our targets, SDSSJ1223-0056.

preprint2014arXiv

Multiwaveband photometry of the irradiated brown dwarf WD0137-349B

WD0137-349 is a white dwarf-brown dwarf binary system in a 116 minute orbit. We present radial velocity observations and multiwaveband photometry from V, R and I in the optical, to J, H and Ks in the near-IR and [3.6], [4.5], [5.8] and [8.0] microns in the mid-IR. The photometry and lightcurves show variability in all wavebands, with the amplitude peaking at [4.5] microns, where the system is also brightest. Fluxes and brightness temperatures were computed for the heated and unheated atmosphere of the brown dwarf (WD0137-349B) using synthetic spectra of the white dwarf using model atmosphere simulations. We show that the flux from the brown dwarf dayside is brighter than expected in the Ks and [4.5] micron bands when compared to models of irradiated brown dwarfs with full energy circulation and suggest this over-luminosity may be attributed to H2 fluorescence or H3+ being generated in the atmosphere by the UV irradiation.

preprint2014arXiv

Precise parameters for both white dwarfs in the eclipsing binary CSS 41177

We present ULTRACAM photometry and X-Shooter spectroscopy of the eclipsing double white dwarf binary CSS 41177, the only such system that is also a double-lined spectroscopic binary. Combined modelling of the light curves and radial velocities yield masses and radii for both white dwarfs without the need to assume mass-radius relations. We find that the primary white dwarf has a mass of M1 = 0.38(2) Msun and a radius of R1 = 0.0222(4) Rsun. The secondary white dwarf's mass and radius are M2 = 0.32(1) Msun and R2 = 0.0207(4) Rsun, and its temperature and surface gravity (T2 = 11678(313) K, log(g2) = 7.32(2)) put it close to the white dwarf instability strip. However, we find no evidence for pulsations to roughly 0.5% relative amplitude. Both masses and radii are consistent with helium white dwarf models with thin hydrogen envelopes of 0.0001 Mstar. The two stars will merge in 1.14 Gyr due to angular momentum loss via gravitational wave emission.

preprint2014arXiv

Spectroscopically confirmed brown dwarf members of Coma Berenices and the Hyades

We have obtained low and medium resolution spectra of 9 brown dwarf candidate members of Coma Berenices and the Hyades using SpEX on the NASA InfaRed Telescope Facility and LIRIS on the William Herschel Telescope. We conclude that 7 of these objects are indeed late M or early L dwarfs, and that two are likely members of Coma Berenices, and four of the Hyades. Two objects, cbd40 and Hy3 are suggested to be a field L dwarfs, although there is also a possibility that Hy3 is an unresolved binary belonging to the cluster. These objects have masses between 71 and 53 M$_{\rm Jup}$, close to the hydrogen burning boundary for these clusters, however only an optical detection of Lithium can confirm if they are truly substellar.

preprint2014arXiv

The substellar companion in the eclipsing white dwarf binary SDSS J141126.20+200911.1

We present high time resolution SDSS-$g'$ and SDSS-$z'$ light curves of the primary eclipse in SDSS J141126.20+200911.1, together with time-resolved X-Shooter spectroscopy and near-infrared $JHK_{s}$ photometry. Our observations confirm the substellar nature of the companion, making SDSS J141126.20+200911.1 the first eclipsing white dwarf/brown dwarf binary known. We measure a (white dwarf model dependent) mass and radius for the brown dwarf companion of $M_{2} = 0.050 \pm 0.002$ $M_{\odot}$ and $R_{2} = 0.072 \pm 0.004$ $M_{\odot}$, respectively. The lack of a robust detection of the companion light in the $z'$-band eclipse constrains the spectral type of the companion to be later than L5. Comparing the NIR photometry to the expected white dwarf flux reveals a clear $K_s$-band excess, suggesting a spectral type in the range L7-T1. The radius measurement is consistent with the predictions of evolutionary models, and suggests a system age in excess of three Gyr. The low companion mass is inconsistent with the inferred spectral type of L7-T1, instead predicting a spectral type nearer T5. This indicates that irradiation of the companion in SDSS J1411 could be causing a significant temperature increase, at least on one hemisphere.

preprint2014arXiv

ULTRASPEC: a high-speed imaging photometer on the 2.4-m Thai National Telescope

ULTRASPEC is a high-speed imaging photometer mounted permanently at one of the Nasmyth focii of the 2.4-m Thai National Telescope (TNT) on Doi Inthanon, Thailand's highest mountain. ULTRASPEC employs a 1024x1024 pixel frame-transfer, electron-multiplying CCD (EMCCD) in conjunction with re-imaging optics to image a field of 7.7'x7.7' at (windowed) frame rates of up to ~200 Hz. The EMCCD has two outputs - a normal output that provides a readout noise of 2.3 e- and an avalanche output that can provide essentially zero readout noise. A six-position filter wheel enables narrow-band and broad-band imaging over the wavelength range 330-1000 nm. The instrument saw first light on the TNT in November 2013 and will be used to study rapid variability in the Universe. In this paper we describe the scientific motivation behind ULTRASPEC, present an outline of its design and report on its measured performance on the TNT.

preprint2013arXiv

A J-band detection of the sub-stellar mass donor in SDSS J1433+1011

We present time-resolved J-band spectroscopy of the short period cataclysmic variable SDSS J143317.78+101123.3. We detect absorption lines from the sub-stellar donor star in this system, which contributes 38 +/- 5% to the J-band light. From the relative strengths of the absorption lines in the J-band, we estimate the spectral type of the donor star to be L2 +/- 1. These data are the first spectroscopic detection of a confirmed sub-stellar donor in a cataclysmic variable, and the spectral type is consistent with that expected from semi-empirical evolutionary models. Using skew mapping, we have been able to derive an estimate for the radial velocity of the donor of Kd = 520 +/- 60 km/s. This value is consistent with, though much less precise than, predictions from mass determinations found via photometric fitting of the eclipse light curves.

preprint2013arXiv

A lithium depletion boundary age of 22 Myr for NGC 1960

We present a deep Cousins RI photometric survey of the open cluster NGC 1960, complete to R_C \simeq 22, I_C \simeq 21, that is used to select a sample of very low-mass cluster candidates. Gemini spectroscopy of a subset of these is used to confirm membership and locate the age-dependent "lithium depletion boundary" (LDB) --the luminosity at which lithium remains unburned in its low-mass stars. The LDB implies a cluster age of 22 +/-4 Myr and is quite insensitive to choice of evolutionary model. NGC 1960 is the youngest cluster for which a LDB age has been estimated and possesses a well populated upper main sequence and a rich low-mass pre-main sequence. The LDB age determined here agrees well with precise age estimates made for the same cluster based on isochrone fits to its high- and low-mass populations. The concordance between these three age estimation techniques, that rely on different facets of stellar astrophysics at very different masses, is an important step towards calibrating the absolute ages of young open clusters and lends confidence to ages determined using any one of them.

preprint2013arXiv

A magnetic white dwarf in a detached eclipsing binary

SDSS J030308.35+005444.1 is a close, detached, eclipsing white dwarf plus M dwarf binary which shows a large infrared excess which has been interpreted in terms of a circumbinary dust disk. In this paper we present optical and near-infrared photometric and spectroscopic data for this system. At optical wavelengths we observe heated pole caps from the white dwarf caused by accretion of wind material from the main-sequence star on to the white dwarf. At near-infrared wavelengths we see the eclipse of two poles on the surface of the white dwarf by the main-sequence star, indicating that the white dwarf is magnetic. Our spectroscopic observations reveal Zeeman split emission lines in the hydrogen Balmer series, which we use to measure the magnetic field strength as 8MG. This measurement indicates that the cyclotron lines are located in the infrared, naturally explaining the infrared excess without the need for a circumbinary dust disk. We also detect magnetically-confined material located roughly midway between the two stars. Using measurements of the radial velocity amplitude and rotational broadening of the M star we constrain the physical parameters of the system, a first for a magnetic white dwarf, and the location of the poles on the surface of the white dwarf. SDSS J030308.35+005444.1 is a pre-cataclysmic variable that will likely evolve into an intermediate polar in ~1Gyr.

preprint2013arXiv

Finding compact hot subdwarf binaries in the Galactic disc

We started a new project which aims to find compact hot subdwarf binaries at low Galactic latitudes. Targets are selected from several photometric surveys and a spectroscopic follow-up campaign to find radial velocity variations on timescales as short as tens of minutes has been started. Once radial variations are detected phase-resolved spectroscopy is obtained to measure the radial velocity curve and the mass function of the system. The observing strategy is described and the discovery of two short period hot subdwarf binaries is presented. UVEXJ032855.25+503529.8 contains a hot subdwarf B star (sdB) orbited by a cool M-dwarf in a P=0.11017 days orbit. The lightcurve shows a strong reflection effect but no eclipses are visible. HS 1741+2133 is a short period (P=0.20 days) sdB most likely with a white dwarf (WD) companion.

preprint2013arXiv

Finding Exoplanets Orbiting Young Active Stars. I. Technique

Stellar activity, such as starspots, can induce radial velocity (RV) variations that can mask or even mimic the RV signature of orbiting exoplanets. For this reason RV exoplanet surveys have been unsuccessful when searching for planets around young, active stars and are therefore failing to explore an important regime which can help to reveal how planets form and migrate. This paper describes a new technique to remove spot signatures from the stellar line-profiles of moderately rotating, active stars (vsini ranging from 10 to 50 km/s). By doing so it allows planetary RV signals to be uncovered. We used simulated models of a G5V type star with differing dark spots on its surface along with archive data of the known active star HD49933 to validate our method. The results showed that starspots could be effectively cleaned from the line-profiles so that the stellar RV jitter was reduced by more than 80%. Applying this procedure to the same models and HD49933 data, but with fake planets injected, enabled the effective removal of starspots so that Jupiter mass planets on short orbital periods were successfully recovered. These results show that this approach can be useful in the search for hot-Jupiter planets that orbit around young, active stars with a vsini of 10 - 50 km/s.

preprint2013arXiv

KIC 4544587: an Eccentric, Short Period Binary System with delta Sct Pulsations and Tidally Excited Modes

We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the delta Scuti and gamma Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software PHOEBE to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98+/-0.07 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do in fact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes.

preprint2013arXiv

Optical transmission photometry of the highly inflated exoplanet WASP-17b

We present ground-based high-precision observations of the transit of WASP-17b using the multi-band photometer ULTRACAM on ESO's NTT in the context of performing transmission spectrophotometry of this highly inflated exoplanet. Our choice of filters (SDSS u', g' and r' bands) is designed to probe for the presence of opacity sources in the upper atmosphere. We find evidence for a wavelength dependence in the planet radius in the form of enhanced absorption in the SDSS r' band, consistent with a previously detected broad sodium feature. We present a new independent measurement of the planetary radius at Rpl = 1.97 +/- 0.06 Rjup, which confirms this planet as the most inflated exoplanet known to date. Our measurements are most consistent with an atmospheric profile devoid of enhanced TiO opacity, previously predicted to be present for this planet.

preprint2013arXiv

Pre-main-sequence isochrones -- II. Revising star and planet formation timescales

We have derived ages for 13 young (<30 Myr) star-forming regions and find they are up to a factor two older than the ages typically adopted in the literature. This result has wide-ranging implications, including that circumstellar discs survive longer (~10-12 Myr) and that the average Class I lifetime is greater (~1 Myr) than currently believed. For each star-forming region we derived two ages from colour-magnitude diagrams. First we fitted models of the evolution between the zero-age main-sequence and terminal-age main-sequence to derive a homogeneous set of main-sequence ages, distances and reddenings with statistically meaningful uncertainties. Our second age for each star-forming region was derived by fitting pre-main-sequence stars to new semi-empirical model isochrones. For the first time (for a set of clusters younger than 50 Myr) we find broad agreement between these two ages, and since these are derived from two distinct mass regimes that rely on different aspects of stellar physics, it gives us confidence in the new age scale. This agreement is largely due to our adoption of empirical colour-Teff relations and bolometric corrections for pre-main-sequence stars cooler than 4000 K. The revised ages for the star-forming regions in our sample are: ~2 Myr for NGC 6611 (Eagle Nebula; M 16), IC 5146 (Cocoon Nebula), NGC 6530 (Lagoon Nebula; M 8), and NGC 2244 (Rosette Nebula); ~6 Myr for σ Ori, Cep OB3b, and IC 348; ~10 Myr for λ Ori (Collinder 69); ~11 Myr for NGC 2169; ~12 Myr for NGC 2362; ~13 Myr for NGC 7160; ~14 Myr for χ Per (NGC 884); and ~20 Myr for NGC 1960 (M 36).

preprint2013arXiv

Seventy six T dwarfs from the UKIDSS LAS: benchmarks, kinematics and an updated space density

We report the discovery of 76 new T dwarfs from the UKIDSS Large Area Survey (LAS). Near-infrared broad and narrow-band photometry and spectroscopy are presented for the new objects, along with WISE and warm-Spitzer photometry. Proper motions for 128 UKIDSS T dwarfs are presented from a new two epoch LAS proper motion catalogue. We use these motions to identify two new benchmark systems: LHS 6176AB, a T8+M4 pair and HD118865AB, a T5.5+F8 pair. Using age constraints from the primaries and evolutionary models to constrain the radii we have estimated their physical properties from their bolometric luminosity. We compare the colours and properties of known benchmark T dwarfs to the latest model atmospheres and draw two principal conclusions. Firstly, it appears that the H - [4.5] and J - W2 colours are more sensitive to metallicity than has previously been recognised, such that differences in metallicity may dominate over differences in Teff when considering relative properties of cool objects using these colours. Secondly, the previously noted apparent dominance of young objects in the late-T dwarf sample is no longer apparent when using the new model grids and the expanded sample of late-T dwarfs and benchmarks. This is supported by the apparently similar distribution of late-T dwarfs and earlier-type T dwarfs on reduced proper motion diagrams that we present. Finally, we present updated space densities for the late-T dwarfs, and compare our values to simulation predictions and those from WISE.

preprint2013arXiv

The planets around NN Ser: still there

We present 25 new eclipse times of the white dwarf binary NN Ser taken with the high-speed camera ULTRACAM on the WHT and NTT, the RISE camera on the Liverpool Telescope, and HAWK-I on the VLT to test the two-planet model proposed to explain variations in its eclipse times measured over the last 25 years. The planetary model survives the test with flying colours, correctly predicting a progressive lag in eclipse times of 36 seconds that has set in since 2010 compared to the previous 8 years of precise times. Allowing both orbits to be eccentric, we find orbital periods of 7.9 +/- 0.5 yr and 15.3 +/- 0.3 yr, and masses of 2.3 +/- 0.5 Mjup and 7.3 +/- 0.3 Mjup. We also find dynamically long-lived orbits consistent with the data, associated with 2:1 and 5:2 period ratios. The data scatter by 0.07 seconds relative to the best-fit model, by some margin the most precise of any of the proposed eclipsing compact object planet hosts. Despite the high precision, degeneracy in the orbit fits prevents a significant measurement of a period change of the binary and of N-body effects. Finally, we point out a major flaw with a previous dynamical stability analysis of NN Ser, and by extension, with a number of analyses of similar systems.

preprint2013arXiv

Timing variations in the secondary eclipse of NN Ser

The eclipsing white dwarf plus main-sequence binary NN Serpentis provides one of the most convincing cases for the existence of circumbinary planets around evolved binaries. The exquisite timing precision provided by the deep eclipse of the white dwarf has revealed complex variations in the eclipse arrival times over the last few decades. These variations have been interpreted as the influence of two planets in orbit around the binary. Recent studies have proved that such a system is dynamically stable over the current lifetime of the binary. However, the existence of such planets is by no means proven and several alternative mechanisms have been proposed that could drive similar variations. One of these is apsidal precession, which causes the eclipse times of eccentric binaries to vary sinusoidally on many year timescales. In this paper we present timing data for the secondary eclipse of NN Ser and show that they follow the same trend seen in the primary eclipse times, ruling out apsidal precession as a possible cause for the variations. This result leaves no alternatives to the planetary interpretation for the observed period variations, although we still do not consider their existence as proven. Our data limits the eccentricity of NN Ser to e<0.001. We also detect a 3.3+/-1.0 second delay in the arrival times of the secondary eclipses relative to the best planetary model. This delay is consistent with the expected 2.84+/-0.04 second Romer delay of the binary, and is the first time this effect has been detected in a white dwarf plus M dwarf system.

preprint2013arXiv

Transmission photometry of WASP-12b: simultaneous measurement of the planetary radius in three bands

Transmission spectroscopy has been successfully used from both the ground and in space to characterise the atmospheres of transiting exoplanets. This technique is challenging from the ground because ground-based spectrographs tend not to be designed to be photometrically stable, and effects such as variable slit losses cause significant systematic uncertainties. An alternative approach is to use simultaneous photometric observations in multiple wavebands to determine wavelength dependent transit depth differences. We report an application of this technique to one of the hottest known exoplanets, WASP-12b, using the triple-beam camera ULTRACAM. We obtained simultaneous light curves in Sloan u', and two narrow band filters centered on 4169 and 6010 angstroms, with FWHMs 52 and 118 angstroms respectively. We fit these light curves with a photometric model and determine the planetary radius in the three different bands. Our data show no evidence for a difference in planetary radius over the wavelength range we study, and are consistent with an atmosphere that is dominated by Rayleigh scattering from a high altitude haze, as well as more complicated atmosphere models which include the effects of molecules such as TiO. Our planetary radius measurements have an average precision of 2.6 per cent, compared to the ~1.4 - 2.4 per cent radius differences predicted by the models over this wavelength range. We also find a consistent time of ingress and egress across our three wavebands, in contrast to the early ingress which has been reported for this system at shorter wavelengths.

preprint2013arXiv

ULTRACAM z'-band Detection of the Secondary Eclipse of WASP-12b

We present z'-band secondary eclipse photometry of the highly irradiated hot Jupiter WASP-12b using ULTRACAM on the 4.2m William Herschel Telescope. We measure a decrease in flux of δ= 0.130+/-0.013% during the passage of the planet behind the star, which is significantly deeper than the previous measurement at this wavelength (0.082+/-0.015%, López-Morales et al. 2010). Our secondary eclipse is best fit with a mid-eclipse phase, ϕ, that is compatible with a circular orbit ϕ= 0.501+/-0.002, in agreement with previous results (Croll et al. 2011). In combination with existing data, our eclipse depth measurement allows us to constrain the characteristics of the planet's atmosphere, which is consistent with a carbon-rich model, with no evidence for a strong thermal inversion. If the difference in eclipse depth reported here compared to that of López-Morales et al. (2010) is of physical origin, as opposed to due to systematics, it may be caused by temporal variability in the flux, due to atmospheric dynamics.

preprint2012arXiv

A Radial Velocity Study of CTCV J1300-3052

We present time-resolved spectroscopy of the eclipsing, short period cataclysmic variable CTCV J1300-3052. Using absorption features from the secondary star, we determine the radial velocity semi-amplitude of the secondary star to be K2 = 378 \pm 6 km/s, and its projected rotational velocity to be v sin i = 125 \pm 7 km/s. Using these parameters and Monte Carlo techniques, we obtain masses of M1 = 0.79 \pm 0.05 MSun for the white dwarf primary and M2 = 0.198 \pm 0.029 MSun for the M-type secondary star. These parameters are found to be in excellent agreement with previous mass determinations found via photometric fitting techniques, supporting the accuracy and validity of photometric mass determinations in short period CVs.

preprint2012arXiv

An accurate mass and radius measurement for an ultracool white dwarf

Studies of cool white dwarfs in the solar neighbourhood have placed a limit on the age of the Galactic disk of 8-9 billion years. However, determining their cooling ages requires the knowledge of their effective temperatures, masses, radii, and atmospheric composition. So far, these parameters could only be inferred for a small number of ultracool white dwarfs for which an accurate distance is known, by fitting their spectral energy distributions (SEDs) in conjunction with a theoretical mass-radius relation. However, the mass-radius relation remains largely untested, and the derived cooling ages are hence model-dependent. Here we report direct measurements of the mass and radius of an ultracool white dwarf in the double-lined eclipsing binary SDSS J013851.54-001621.6. We find M(WD)=0.529+/-0.010Msol and R(WD)=0.0131+/-0.0003Rsol. Our measurements are consistent with the mass-radius relation and we determine a robust cooling age of 9.5 billion years for the 3570K white dwarf. We find that the mass and radius of the low mass companion star, M(sec)=0.132+/-0.003Msol and R(sec)=0.165+/-0.001Rsol, are in agreement with evolutionary models. We also find evidence that this >9.5 Gyr old M5 star is still active, far beyond the activity lifetime for a star of its spectral type. This is likely caused by the high tidally-enforced rotation rate of the star. The companion star is close to filling its Roche lobe and the system will evolve into a cataclysmic variable in only 70 Myr. Our direct measurements demonstrate that this system can be used to calibrate ultracool white dwarf atmospheric models.

preprint2012arXiv

Eclipsing Post Common Envelope Binaries from the Catalina Surveys

We analyse the Catalina Real-time Transient Survey light curves of 835 spectroscopically confirmed white dwarf plus main-sequence binaries from the Sloan Digital Sky Survey (SDSS) with g<19, in search of new eclipsing systems. We identify 29 eclipsing systems, 12 of which were previously unknown. This brings the total number of eclipsing white dwarf plus main-sequence binaries to 49. Our set of new eclipsing systems contains two with periods of 1.9 and 2.3 days making them the longest period eclipsing white dwarf binaries known. We also identify one system which shows very large ellipsoidal modulation (almost 0.3 magnitudes), implying that the system is both very close to Roche-lobe overflow and at high inclination. However, our follow up photometry failed to firmly detect an eclipse meaning that this system either contains a cool white dwarf and hence the eclipse is very shallow and undetectable in our red-sensitive photometry or that it is non-eclipsing. Radial velocity measurements for the main-sequence stars in three of our newly identified eclipsing systems imply that their white dwarf masses are lower than those inferred from modelling their SDSS spectra. 13 non-eclipsing post common envelope binaries were also identified, from either reflection or ellipsoidal modulation effects. The white dwarfs in our newly discovered eclipsing systems span a wide range of parameters, including; low mass (~0.3Msol), very hot (80,000K) and a DC white dwarf. The spectral types of the main-sequence stars range from M2 to M6. This makes our sample ideal for testing white dwarf and low-mass star mass-radius relationships as well as close binary evolution.

preprint2012arXiv

Pre-main-sequence isochrones -- I. The Pleiades benchmark

We present a critical assessment of commonly used pre-main-sequence isochrones by comparing their predictions to a set of well-calibrated colour-magnitude diagrams of the Pleiades in the wavelength range 0.4 to 2.5 microns. Our analysis shows that for temperatures less than 4000 K the models systematically overestimate the flux by a factor two at 0.5 microns, though this decreases with wavelength, becoming negligible at 2.2 microns. In optical colours this will result in the ages for stars younger than 10 Myr being underestimated by factors between two and three. We show that using observations of standard stars to transform the data into a standard system can introduce significant errors in the positioning of pre-main-sequences in colour-magnitude diagrams. Therefore we have compared the models to the data in the natural photometric system in which the observations were taken. Thus we have constructed and tested a model of the system responses for the Wide-Field Camera on the Isaac Newton Telescope. As a benchmark test for the development of pre-main-sequence models we provide both our system responses and the Pleiades sequence.

preprint2012arXiv

z'-band Ground-Based Detection of the Secondary Eclipse of WASP-19b

We present the ground-based detection of the secondary eclipse of the transiting exoplanet WASP-19b. The observations were made in the Sloan z'-band using the ULTRACAM triple-beam CCD camera mounted on the NTT. The measurement shows a 0.088\pm0.019% eclipse depth, matching previous predictions based on H- and K-band measurements. We discuss in detail our approach to the removal of errors arising due to systematics in the data set, in addition to fitting a model transit to our data. This fit returns an eclipse centre, T0, of 2455578.7676 HJD, consistent with a circular orbit. Our measurement of the secondary eclipse depth is also compared to model atmospheres of WASP-19b, and is found to be consistent with previous measurements at longer wavelengths for the model atmospheres we investigated.

preprint2011arXiv

A J-band detection of the donor star in the dwarf nova OY Carinae, and an optical detection of its `iron curtain'

Purely photometric models can be used to determine the binary parameters of eclipsing cataclysmic variables with a high degree of precision. However, the photometric method relies on a number of assumptions, and to date there have been very few independent checks of this method in the literature. We present time-resolved spectroscopy of the P=90.9 min eclipsing cataclysmic variable OY Carinae obtained with X-shooter on the VLT, in which we detect the donor star from K I lines in the J-band. We measure the radial velocity amplitude of the donor star K2 = 470.0 +/- 2.7 km/s, consistent with predictions based upon the photometric method (470 +/- 7 km/s). Additionally, the spectra obtained in the UVB arm of X-shooter show a series of Fe I and Fe II lines with a phase and velocity consistent with an origin in the accretion disc. This is the first unambiguous detection at optical wavelengths of the `iron curtain' of disc material which has been previously reported to veil the white dwarf in this system. The velocities of these lines do not track the white dwarf, reflecting a distortion of the outer disc that we see also in Doppler images. This is evidence for considerable radial motion in the outer disk, at up to 90 km/s towards and away from the white dwarf.

preprint2011arXiv

A precision study of two eclipsing white dwarf plus M dwarf binaries

We use a combination of X-shooter spectroscopy, ULTRACAM high-speed photometry and SOFI near-infrared photometry to measure the masses and radii of both components of the eclipsing post common envelope binaries SDSS J1212-0123 and GK Vir. For both systems we measure the gravitational redshift of the white dwarf and combine it with light curve model fits to determine the inclinations, masses and radii. For SDSS J1212-0123 we find a white dwarf mass and radius of 0.439 +/- 0.002 Msun and 0.0168 +/- 0.0003 Rsun, and a secondary star mass and radius of 0.273 +/- 0.002 Msun and 0.306 +/- 0.007 Rsun. For GK Vir we find a white dwarf mass and radius of 0.564 +/- 0.014 Msun and 0.0170 +/- 0.0004 Rsun, and a secondary star mass and radius of 0.116 +/- 0.003 Msun and 0.155 +/- 0.003 Rsun. The mass and radius of the white dwarf in GK Vir are consistent with evolutionary models for a 50,000K carbon-oxygen core white dwarf. Although the mass and radius of the white dwarf in SDSS J1212-0123 are consistent with carbon-oxygen core models, evolutionary models imply that a white dwarf with such a low mass and in a short period binary must have a helium core. The mass and radius measurements are consistent with helium core models but only if the white dwarf has a very thin hydrogen envelope, which has not been predicted by evolutionary models. The mass and radius of the secondary star in GK Vir are consistent with evolutionary models after correcting for the effects of irradiation by the white dwarf. The secondary star in SDSS J1212-0123 has a radius ~9 per cent larger than predicted.

preprint2011arXiv

A search for optical bursts from RRAT J1819-1458: II. Simultaneous ULTRACAM-Lovell Telescope observations

The Rotating RAdio Transient (RRAT) J1819-1458 exhibits ~3 ms bursts in the radio every ~3 min, implying that it is visible for only ~1 per day. Assuming that the optical light behaves in a similar manner, long exposures of the field would be relatively insensitive due to the accumulation of sky photons. A much better way of detecting optical emission from J1819-1458 would then be to observe with a high-speed optical camera simultaneously with radio observations, and co-add only those optical frames coincident with the dispersion-corrected radio bursts. We present the results of such a search, using simultaneous ULTRACAM and Lovell Telescope observations. We find no evidence for optical bursts in J1819-1458 at magnitudes brighter than i'=19.3 (5-sigma limit). This is nearly 3 magnitudes fainter than the previous burst limit, which had no simultaneous radio observations.

preprint2011arXiv

Accretion-induced luminosity spreads in young clusters: evidence from stellar rotation

We present an analysis of the rotation of young stars in the associations Cepheus OB3b, NGC 2264, NGC 2362 and the Orion Nebula Cluster (ONC). We discover a correlation between rotation rate and position in a colour-magnitude diagram (CMD) such that stars which lie above an empirically determined median pre-main sequence rotate more rapidly than stars which lie below this sequence. The same correlation is seen, with a high degree of statistical significance, in each association studied here. If position within the CMD is interpreted as being due to genuine age spreads within a cluster, then the stars above the median pre-main sequence would be the youngest stars. This would in turn imply that the most rapidly rotating stars in an association are the youngest, and hence those with the largest moments of inertia and highest likelihood of ongoing accretion. Such a result does not fit naturally into the existing picture of angular momentum evolution in young stars, where the stars are braked effectively by their accretion discs until the disc disperses. Instead, we argue that, for a given association of young stars, position within the CMD is not primarily a function of age, but of accretion history. We show that this hypothesis could explain the correlation we observe between rotation rate and position within the CMD.

preprint2011arXiv

Cataclysmic Variables below the Period Gap: Mass Determinations of 14 Eclipsing Systems

We present high-speed, three-colour photometry of the eclipsing cataclysmic variables CTCV 1300, CTCV 2354 and SDSS 1152. All three systems are below the observed "period gap" for cataclysmic variables. For each system we determine the system parameters by fitting a parameterised model to the observed eclipse light curve by chi-squared minimisation. We also present an updated analysis of all other eclipsing systems previously analysed by our group. New donor masses are generally between 1 and 2 sigma of those originally published, with the exception of SDSS 1502 and DV UMa. We note that the donor mass of SDSS 1501 has been revised upwards by 0.024Msun. This system was previously identified as having evolved passed the minimum orbital period for cataclysmic variables, but the new mass determination suggests otherwise. Our new analysis confirms that SDSS 1035 and SDSS 1433 have evolved past the period minimum for cataclysmic variables, corroborating our earlier studies. We find that the radii of donor stars are oversized when compared to theoretical models, by approximately 10 percent. We show that this can be explained by invoking either enhanced angular momentum loss, or by taking into account the effects of star spots. We are unable to favour one cause over the other, as we lack enough precise mass determinations for systems with orbital periods between 100 and 130 minutes, where evolutionary tracks begin to diverge significantly. We also find a strong tendency towards high white dwarf masses within our sample, and no evidence for any He-core white dwarfs. The dominance of high mass white dwarfs implies that erosion of the white dwarf during the nova outburst must be negligible, or that not all of the mass accreted is ejected during nova cycles, resulting in the white dwarf growing in mass. (Abridged)

preprint2011arXiv

No wide spread of stellar ages in the Orion Nebula Cluster

The wide luminosity dispersion seen for stars at a given effective temperature in the H-R diagrams of young clusters and star forming regions is often interpreted as due to significant (~10 Myr) spreads in stellar contraction age. In the scenario where most stars are born with circumstellar discs, and that disc signatures decay monotonically (on average) over timescales of only a few Myr, then any such age spread should lead to clear differences in the age distributions of stars with and without discs. We have investigated large samples of stars in the Orion Nebula Cluster (ONC) using three methods to diagnose disc presence from infrared measurements. We find no significant difference in the mean ages or age distributions of stars with and without discs, consistent with expectations for a coeval population. Using a simple quantitative model we show that any real age spread must be smaller than the median disc lifetime. For a log-normal age distribution, there is an upper limit of <0.14 dex (at 99% confidence) to any real age dispersion, compared to the ~=0.4 dex implied by the H-R diagram. If the mean age of the ONC is 2.5 Myr, this would mean at least 95% of its low-mass stellar population has ages between 1.3--4.8 Myr. We suggest that the observed luminosity dispersion is caused by a combination of observational uncertainties and physical mechanisms that disorder the conventional relationship between luminosity and age for pre main-sequence stars. This means that individual stellar ages from the H-R diagram are unreliable and cannot be used to directly infer a star formation history. Irrespective of what causes the wide luminosity dispersion, the finding that any real age dispersion is less than the median disc lifetime argues strongly against star formation scenarios for the ONC lasting longer than a few Myr.

preprint2011arXiv

On the alignment of debris disks and their host stars' rotation axis -implications for spin-orbit misalignment in exoplanetary systems

It has been widely thought that measuring the misalignment angle between the orbital plane of a transiting exoplanet and the spin of its host star was a good discriminator between different migration processes for hot-Jupiters. Specifically, well-aligned hot-Jupiter systems (as measured by the Rossiter-McLaughlin effect) were thought to have formed via migration through interaction with a viscous disk, while misaligned systems were thought to have undergone a more violent dynamical history. These conclusions were based on the assumption that the planet-forming disk was well-aligned with the host star. Recent work by a number of authors has challenged this assumption by proposing mechanisms that act to drive the star-disk interaction out of alignment during the pre-main sequence phase. We have estimated the stellar rotation axis of a sample of stars which host spatially resolved debris disks. Comparison of our derived stellar rotation axis inclination angles with the geometrically measured debris-disk inclinations shows no evidence for a misalignment between the two.

preprint2011arXiv

The first observation of optical pulsations from a soft gamma repeater: SGR 0501+4516

We present high-speed optical photometry of the soft gamma repeater SGR 0501+4516, obtained with ULTRACAM on two consecutive nights approximately 4 months after the source was discovered via its gamma-ray bursts. We detect SGR 0501+4516 at a magnitude of i' = 24.4+/-0.1. We present the first measurement of optical pulsations from an SGR, deriving a period of 5.7622+/-0.0003 s, in excellent agreement with the X-ray spin period of the neutron star. We compare the morphologies of the optical pulse profile with the X-ray and infrared pulse profiles; we find that the optical, infrared and harder X-rays share similar double-peaked morphologies, but the softer X-rays exhibit only a single-peaked morphology, indicative of a different origin. The optical pulsations appear to be in phase with the X-ray pulsations and exhibit a root-mean-square pulsed fraction of 52+/-7%, approximately a factor of two greater than in the X-rays. Our results find a natural explanation within the context of the magnetar model for SGRs.

preprint2011arXiv

The photometric period in ES Ceti

We present ULTRACAM photometry of ES Cet, an ultracompact binary with a 620s orbital period. The mass transfer in systems such as this one is thought to be driven by gravitational radiation, which causes the binary to evolve to longer periods since the semi-degenerate donor star expands in size as it loses mass. We supplement these ULTRACAM+WHT data with observations made with smaller telescopes around the world over a nine year baseline. All of the observations show variation on the orbital period, and by timing this variation we track the period evolution of this system. We do not detect any significant departure from a linear ephemeris, implying a donor star that is of small mass and close to a fully degenerate state. This finding favours the double white dwarf formation channel for this AM CVn star. An alternative explanation is that the system is in the relatively short-lived phase in which the mass transfer rate climbs towards its long-term value.

preprint2011arXiv

The Planets around the Post-Common Envelope Binary NN Serpentis

We have detected 2 circumbinary planets around the close binary system NN Serpentis using the orbital time delay effect measured via the sharp eclipses of the white dwarf primary. The present pre-cataclysmic binary was formed when the original - 2 M primary expanded into a red giant, causing the secondary star to drop from its original orbit at a separation of about 1.4 A.U. down to its current separation at 0.0043 A.U. A quasi-adiabatic evolution of the circumbinary planets' orbits during the common-envelope phase would have placed them in unstable configurations, suggesting that they may have suffered significant orbital drag effects and were originally in much larger orbits. Alternatively, they may have been created as 2nd-generation planets during the last million years from the substantial amount of material lost during the creation of the binary, making them the youngest planets known. Either solution shows how little we actually understand about planetary formation.

preprint2011arXiv

The shortest period detached white dwarf + main-sequence binary

We present high-speed ULTRACAM and SOFI photometry and X-shooter spectroscopy of the recently discovered 94 minute orbital period eclipsing white dwarf / main-sequence binary SDSS J085746.18+034255.3 (CSS 03170) and use these observations to measure the system parameters. We detect a shallow secondary eclipse and hence are able to determine an orbital inclination of 85.5 +/- 0.2 deg. The white dwarf has a mass of 0.51 +/- 0.05 Msun and a radius of 0.0247 +/- 0.0008 Rsun. With a temperature of 35,300 +/- 400K the white dwarf is highly over-inflated if it has a carbon-oxygen core, however if it has a helium core then its mass and radius are consistent with evolutionary models. Therefore, the white dwarf in SDSS J085746.18+034255.3 is most likely a helium core white dwarf with a mass close to the upper limit expected from evolution. The main-sequence star is an M8 dwarf with a mass of 0.09 +/- 0.01 Msun and a radius of 0.110 +/- 0.004 Rsun placing it close to the hydrogen burning limit. The system emerged from a common envelope ~20 million years ago and will reach a semi-detached configuration in ~400 million years, becoming a cataclysmic variable with a period of 66 minutes, below the period minimum.

preprint2010arXiv

Estimating the masses of extra-solar planets

All extra-solar planet masses that have been derived spectroscopically are lower limits since the inclination of the orbit to our line-of-sight is unknown except for transiting systems. It is, however, possible to determine the inclination angle, i, between the rotation axis of a star and an observer's line-of-sight from measurements of the projected equatorial velocity (v sin i), the stellar rotation period (P_rot) and the stellar radius (R_star). This allows the removal of the sin i dependency of spectroscopically derived extra-solar planet masses under the assumption that the planetary orbits lie perpendicular to the stellar rotation axis. We have carried out an extensive literature search and present a catalogue of v sin i, P_rot, and R_star estimates for exoplanet host stars. In addition, we have used Hipparcos parallaxes and the Barnes-Evans relationship to further supplement the R_star estimates obtained from the literature. Using this catalogue, we have obtained sin i estimates using a Markov-chain Monte Carlo analysis. This allows proper 1-sigma two-tailed confidence limits to be placed on the derived sin i's along with the transit probability for each planet to be determined. While a small proportion of systems yield sin i's significantly greater than 1, most likely due to poor P_rot estimations, the large majority are acceptable. We are further encouraged by the cases where we have data on transiting systems, as the technique indicates inclinations of ~90 degrees and high transit probabilities. In total, we estimate the true masses of 133 extra-solar planets. Of these, only 6 have revised masses that place them above the 13 Jupiter mass deuterium burning limit. Our work reveals a population of high-mass planets with low eccentricities and we speculate that these may represent the signature of different planetary formation mechanisms at work.

preprint2010arXiv

Orbital Period Variations in Eclipsing Post Common Envelope Binaries

We present high speed ULTRACAM photometry of the eclipsing post common envelope binaries DE CVn, GK Vir, NN Ser, QS Vir, RR Cae, RX J2130.6+4710, SDSS 0110+1326 and SDSS 0303+0054 and use these data to measure precise mid-eclipse times in order to detect any period variations. We detect a large (~ 250 sec) departure from linearity in the eclipse times of QS Vir which Applegate's mechanism fails to reproduce by an order of magnitude. The only mechanism able to drive this period change is a third body in a highly elliptical orbit. However, the planetary/sub-stellar companion previously suggested to exist in this system is ruled out by our data. Our eclipse times show that the period decrease detected in NN Ser is continuing, with magnetic braking or a third body the only mechanisms able to explain this change. The planetary/sub-stellar companion previously suggested to exist in NN Ser is also ruled out by our data. Our precise eclipse times also lead to improved ephemerides for DE CVn and GK Vir. The width of a primary eclipse is directly related to the size of the secondary star and variations in the size of this star could be an indication of Applegate's mechanism or Wilson (starspot) depressions which can cause jitter in the O-C curves. We measure the width of primary eclipses for the systems NN Ser and GK Vir over several years but find no definitive variations in the radii of the secondary stars. However, our data are precise enough (dRsec / Rsec < 10^-5) to show the effects of Applegate's mechanism in the future. We find no evidence of Wilson depressions in either system. We also find tentative indications that flaring rates of the secondary stars depend on their mass rather than rotation rates.

preprint2010arXiv

Probing GRB environments with time variability: ULTRASPEC fast imaging of GRB 080210

We present high time resolution (1.09 s) photometry of GRB 080210 obtained with ULTRASPEC mounted on the ESO/3.6-m telescope, starting 68.22 min after the burst and lasting for 26.45 min. The light curve is smooth on both short (down to 2.18 s) and long time scales, confirmed by a featureless power spectrum. On top of the fireball power-law decay, bumps and wiggles at different time scales can, in principle, be produced by density fluctuations in the circumburst medium, substructures in the jet or by refreshed shocks. Comparing our constraints with variability limits derived from kinematic arguments, we exclude under-density fluctuations producing flux dips larger than 1 per cent with time scales \Deltat > 9.2 min (2 per cent on \Deltat > 2.3 min for many fluctuating regions). In addition, we study the afterglow VLT/FORS2 spectrum, the optical-to-X-ray spectral energy distribution (SED) and the time decay. The SED is best fit with a broken power law with slopes βopt = 0.71 \pm 0.01 and βX = 1.59 \pm 0.07, in disagreement with the fireball model, suggesting a non-standard afterglow for GRB 080210. We find AV = 0.18 \pm 0.03 mag optical extinction due to SMC-like dust and an excess X-ray absorption of log (NH/cm-2) = 21.58 +0.18 -0.26 assuming Solar abundances. The spectral analysis reveals a damped Lyα absorber (log (NH I /cm-2) = 21.90 \pm 0.10) with a low metallicity ([X/H] = -1.21 \pm 0.16), likely associated with the interstellar medium of the GRB host galaxy (z = 2.641).

preprint2010arXiv

SDSS J0926+3624: The shortest period eclipsing binary star

With orbital periods of the order of tens of minutes or less, the AM Canum Venaticorum stars are ultracompact, hydrogen deficient binaries with the shortest periods of any binary subclass, and are expected to be among the strongest gravitational wave sources in the sky. To date, the only known eclipsing source of this type is the P = 28 min binary SDSS J0926+3624. We present multiband, high time resolution light curves of this system, collected with WHT/ULTRACAM in 2006 and 2009. We supplement these data with additional observations made with LT/RISE, XMM_Newton and the Catalina Real-Time Transient Survey. From light curve models we determine the mass ratio to be q = M2 / M1 = 0.041 +/- 0.002 and the inclination to be 82.6 +/- 0.3 deg. We calculate the mass of the primary white dwarf to be 0.85 +/- 0.04 solar masses and the donor to be 0.035 +/- 0.003 solar masses, implying a partially degenerate state for this component. We observe superhump variations that are characteristic of an elliptical, precessing accretion disc. Our determination of the superhump period excess is in agreement with the established relationship between this parameter and the mass ratio, and is the most precise calibration of this relationship at low q. We also observe a quasi-periodic oscillation in the 2006 data, and we examine the outbursting behaviour of the system over a 4.5 year period.

preprint2010arXiv

Two planets orbiting the recently formed post-common envelope binary NN Serpentis

Planets orbiting post-common envelope binaries provide fundamental information on planet formation and evolution. We searched for such planets in NN Ser ab, an eclipsing short-period binary that shows long-term eclipse time variations. Using published, reanalysed, and new mid-eclipse times of NN Ser ab obtained between 1988 and 2010, we find excellent agreement with the light-travel-time effect by two additional bodies superposed on the linear ephemeris of the binary. Our multi-parameter fits accompanied by N-body simulations yield a best fit for the objects NN Ser (ab)c and d locked in a 2:1 mean motion resonance, with orbital periods P_c=15.5 yrs and P_d=7.7 yrs, masses M_c sin i_c = 6.9 M_Jup and M_d sin i_d = 2.2 M_Jup, and eccentricities e_c=0 and e_d=0.20. A secondary chi**2 minimum corresponds to an alternative solution with a period ratio of 5:2. We estimate that the progenitor binary consisted of an A star with ~2 M_Sun and the present M dwarf secondary at an orbital separation of ~1.5 AU. The survival of two planets through the common-envelope phase that created the present white dwarf requires fine tuning between the gravitational force and the drag force experienced by them in the expanding envelope. The alternative is a second-generation origin in a circumbinary disk created at the end of this phase. In that case, the planets would be extremely young with ages not exceeding the cooling age of the white dwarf of 10**6 yrs.

preprint2009arXiv

Physical properties of IP Pegasi: an eclipsing dwarf nova with an unusually cool white dwarf

We present high speed photometric observations of the eclipsing dwarf nova IP Peg taken with the triple-beam camera ULTRACAM mounted on the William Herschel Telescope. The primary eclipse in this system was observed twice in 2004, and then a further sixteen times over a three week period in 2005. Our observations were simultaneous in the Sloan u', g' and r' bands. By phase-folding and averaging our data we make the first significant detection of the white dwarf ingress in this system and find the phase width of the white dwarf eclipse to be 0.0935 +/- 0.0003, significantly higher than the previous best value of between 0.0863 and 0.0918. The mass ratio is found to be q = M2 /M1 = 0.48 +/- 0.01, consistent with previous measurements, but we find the inclination to be 83.8 +/- 0.5 deg, significantly higher than previously reported. We find the radius of the white dwarf to be 0.0063 +/- 0.0003 solar radii, implying a white dwarf mass of 1.16 +/- 0.02 solar masses. The donor mass is 0.55 +/- 0.02 solar masses. The white dwarf temperature is more difficult to determine, since the white dwarf is seen to vary significantly in flux, even between consecutive eclipses. This is seen particularly in the u'-band, and is probably the result of absorption by disc material. Our best estimate of the temperature is 10,000 - 15,000K, which is much lower than would be expected for a CV with this period, and implies a mean accretion rate of less than 5 times 10^-11 solar masses per year, more than 40 times lower than the expected rate.

preprint2009arXiv

Pre-main-sequence variability across the radiative-convective gap

We use I band imaging to perform a variability survey of the 13 Myr-old cluster h Per. We find a significant fraction of the cluster members to be variable. Most importantly, we find that variable members lie almost entirely on the convective side of the gap in the cluster sequence between fully convective stars and those which have a radiative core. This result is consistent with a scenario in which the magnetic field changes topology when the star changes from being fully convective, to one containing a radiative core. When the star is convective the magnetic field appears dominated by large-scale structures, resulting in global-size spots that drive the observed variability. For those stars with radiative cores we observe a marked absence of variability due to spots, which suggests a switch to a magnetic field dominated by smaller-scale structures, resulting in many smaller spots and thus less apparent variability. This implies that wide field variability surveys may only be sensitive to fully convective stars. On the one hand this reduces the chances of picking out young groups (since the convective stars are the lower mass and therefore fainter objects), but conversely the absolute magnitude of the head of the convective sequence provides a straightforward measure of age for those groups which are discovered.

preprint2009arXiv

Precise mass and radius values for the white dwarf and low mass M dwarf in the pre-cataclysmic binary NN Serpentis

We derive precise system parameters for the pre-cataclysmic binary, NN Ser. From light curve fitting we find an orbital inclination of i = 89.6 +/- 0.2 deg. From the HeII absorption line we find K_{WD}= 62.3 +/- 1.9 km/s. The irradiation-induced emission lines from the surface of the secondary star give a range of observed radial velocities. The corrected values give a radial velocity of K_{sec}= 301 +/- 3 km/s, with an error dominated by the systematic effects of the model. This leads to a binary separation of a = 0.934 +/- 0.009 R_{sun}, radii of R_{WD} = 0.0211 +/- 0.0002 R_{sun} and R_{sec} = 0.149 +/- 0.002 R_{sun} and masses of M_{WD} = 0.535 +/- 0.012 M_{sun} and M_{sec} = 0.111 +/- 0.004 M_{sun}. The masses and radii of both components of NN Ser were measured independently of any mass-radius relation. For the white dwarf, the measured mass, radius and temperature show excellent agreement with a `thick' hydrogen layer of fractional mass M_{H}/{M}_{WD} = 10^{-4}. The measured radius of the secondary star is 10% larger than predicted by models, however, correcting for irradiation accounts for most of this inconsistency, hence the secondary star in NN Ser is one of the first precisely measured very low mass objects to show good agreement with models. ULTRACAM r', i' and z' photometry taken during the primary eclipse determines the colours of the secondary star as (r'-i')_{sec}= 1.4 +/- 0.1 and (i'-z')_{sec} = 0.8 +/- 0.1 which corresponds to a spectral type of M4 +/- 0.5. This is consistent with the derived mass, demonstrating that there is no detectable heating of the unirradiated face, despite intercepting radiative energy from the white dwarf which exceeds its own luminosity by over a factor of 20.

preprint2009arXiv

Rotation of young stars in Cepheus OB3b

We present a photometric study of I-band variability in the young association Cepheus OB3b. The study is sensitive to periodic variability on timescales of less than a day, to more than 20 days. After rejection of contaminating objects using V, I, R and narrowband H-alpha photometry, we find 475 objects with measured rotation periods, which are very likely pre-main-sequence members of the Cep OB3b star forming region. We revise the distance and age to Cep OB3b, putting it on the self-consistent age and distance ladder of Mayne & Naylor (2008). This yields a distance modulus of 8.8 +/- 0.2 mags, corresponding to a distance of 580 +/- 60 pc, and an age of 4-5Myrs. The rotation period distribution confirms the general picture of rotational evolution in young stars, exhibiting both the correlation between accretion (determined in this case through narrowband H-alpha photometry) and rotation expected from disc locking, and the dependence of rotation upon mass that is seen in other star forming regions. However, this mass dependence is much weaker in our data than found in other studies. Comparison to the similarly aged NGC 2362 shows that the low-mass stars in Cep OB3b are rotating much more slowly. This points to a possible link between star forming environment and rotation properties. Such a link would call into question models of stellar angular momentum evolution, which assume that the rotational period distributions of young clusters and associations can be assembled into an evolutionary sequence, thus ignoring environmental effects.

preprint2009arXiv

SDSS unveils a population of intrinsically faint cataclysmic variables at the minimum orbital period

We discuss the properties of 137 cataclysmic variables (CVs) which are included in the Sloan Digital Sky Survey (SDSS) spectroscopic data base, and for which accurate orbital periods have been measured. 92 of these systems are new discoveries from SDSS and were followed-up in more detail over the past few years. 45 systems were previously identified as CVs because of the detection of optical outbursts and/or X-ray emission, and subsequently re-identified from the SDSS spectroscopy. The period distribution of the SDSS CVs differs dramatically from that of all the previously known CVs, in particular it contains a significant accumulation of systems in the orbital period range 80--86 min. We identify this feature as the elusive "period minimum spike" predicted by CV population models, which resolves a long-standing discrepancy between compact binary evolution theory and observations. We show that this spike is almost entirely due to the large number of CVs with very low accretion activity identified by SDSS. The optical spectra of these systems are dominated by emission from the white dwarf photosphere, and display little or no spectroscopic signature from the donor stars, suggesting very low-mass companion stars. We determine the average absolute magnitude of these low-luminosity CVs at the period minimum to be <M_g>=11.6+-0.7. Comparison of the SDSS CV sample to the CVs found in the Hamburg Quasar Survey and the Palomar Green Survey suggests that the depth of SDSS is the key ingredient resulting in the discovery of a large number of intrinsically faint short-period systems.

preprint2004arXiv

Contamination and exclusion in the sigma Orionis young group

We present radial velocities for 38 low-mass candidate members of the sigma Orionis young group. We have measured their radial velocities by cross-correlation of high resolution (R~6000) AF2/WYFFOS spectra of the gravity sensitive NaI doublet at 8183, 8195Angstroms. The total sample contained 117 objects of which 54 have sufficient signal-to-noise to detect NaI at an equivalent width of 3Angstroms, however we only detect NaI in 38 of these. This implies that very low-mass members of this young group display weaker NaI absorption than similarly aged objects in the Upper Scorpius OB association. We develop a technique to assess membership using radial velocities with a range of uncertainties that does not bias the selection when large uncertainties are present. The resulting membership probabilities are used to assess the issue of exclusion in photometric selections, and we find that very few members are likely to be excluded by such techniques. We also assess the level of contamination in the expected pre-main sequence region of colour-magnitude space brighter than I = 17. We find that contamination by non-members in the expected PMS region of the colour-magnitude diagram is small. We conclude that although radial velocity alone is insufficient to confirm membership, high signal-to-noise observations of the NaI doublet provide the opportunity to use the strength of NaI absorption in concert with radial velocities to asses membership down to the lowest masses, where Lithium absorption no longer distinguishes youth.

preprint2003arXiv

Evidence for high accretion-rates in Weak-Line T Tauri stars?

We have discovered T Tauri stars which show startling spectral variability between observations seperated by 20 years. In spectra published by Bouvier & Appenzeller (1992) these objects showed very weak H-alpha emission, broad CaII absorption and so called ``composite spectra'', where the spectral type inferred from the blue region is earlier than that inferred from the red. We present here new spectroscopy which shows that all four stars now exhibit strong H-alpha emission, narrow CaII emission and a spectral type which is consistent at all wavelengths. We propose a scheme to understand these changes whereby the composite spectra of these stars can be explained by a period of active accretion onto the central, young star. In this scheme the composite spectrum consists of a contribution from the stellar photosphere and a contribution from a hot, optically thick, accretion component. The optically thick nature of the accretion flow explains the weakness of the H-alpha emission during this phase. Within this scheme, the change to a single spectral type at all wavelengths and emergence of strong H-alpha emission are consistent with the accretion columns becoming optically thin, as the accretion rate drops. There is a strong analogy here with the dwarf novae class of interacting binaries, which show similar behaviour during the decline from outbursts of high mass-transfer rate. The most important consequence of this interpretation is that these objects bring into question the association of Weak-Line T Tauri stars (WTTs) with non-accreting or discless objects. In light of this result we consider the justification for this paradigm.