Search

Stellar active regions like spots and faculae can distort the shapes of spectral lines, inducing variations in the radial velocities that are often orders of magnitude larger than the signals from Earth-like planets. Efforts to mitigate these activity signals have hitherto focused on either the time or the velocity (wavelength) domains. We present a physics-driven Gaussian process (GP) framework to model activity signals directly in time series of line profiles or cross-correlation functions (CCFs). Unlike existing methods that correct activity signals in line profile time series, our approach exploits the time correlation between velocity (wavelength) bins in the line profile variations, and is based on a simplified but physically motivated model for the origin of these variations. When tested on both synthetic and real data sets with signal-to-noise ratios down to ∼100, our method was able to separate the planetary signal from the activity signal, even when their periods were identical. We also conducted injection/recovery tests using two years of realistically sampled HARPS-N solar data, demonstrating the ability of the method to accurately recover a signal induced by a 1.5-Earth mass planet with a semi-amplitude of 0.3 m s−1 and a period of 33 d during high solar activity. [ABSTRACT FROM AUTHOR]

Stellar magnetic activity induces both distortions and Doppler-shifts in the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the Sun's rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian processes (GPs). Planet signatures are still best retrieved with multidimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca  ii H & K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6–0.8 m s−1, likely to be dominated by signals induced by supergranulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign. [ABSTRACT FROM AUTHOR]

As the hunt for an Earth-like exoplanets has intensified in recent years, so has the effort to characterize and model the stellar signals that can hide or mimic small planetary signals. Stellar variability arises from a number of sources, including granulation, supergranulation, oscillations, and activity, all of which result in quasi-periodic or stochastic behaviour in photometric and/or radial velocity observations. Traditionally, the characterization of these signals has mostly been done in the frequency domain. However, the recent development of scalable Gaussian process regression methods makes direct time-domain modelling of stochastic processes a feasible and arguably preferable alternative, obviating the need to estimate the power spectral density of the data before modelling it. In this paper, we compare the two approaches using a series of experiments on simulated data. We show that frequency-domain modelling can lead to inaccurate results, especially when the time-sampling is irregular. By contrast, Gaussian process regression results are often more precise, and systematically more accurate, in both the regular and irregular time-sampling regimes. While this work was motivated by the analysis of radial velocity and photometry observations of main-sequence stars in the context of planet searches, we note that our results may also have applications for the study of other types of astrophysical variability such as quasi-periodic oscillations in X-ray binaries and active galactic nuclei variability. [ABSTRACT FROM AUTHOR]

We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837 b, with an orbital period of ∼8.32 d. Utilizing data from the Transiting Exoplanet Survey Satellite and ground-based observations, we determine a planetary radius of |$0.818_{-0.024}^{+0.034}$| R J for TOI-837 b. Through detailed High Accuracy Radial Velocity Planet Searcher spectroscopic time series analysis, we derive a Doppler semi-amplitude of |$34.7_{-5.6}^{+5.3}$| |${\rm m\, s^{-1}}$|⁠ , corresponding to a planetary mass of |$0.379_{-0.061}^{+0.058}$| M J. The derived planetary properties suggest a substantial core of approximately 70 M⊕, constituting about 60 per cent of the planet's total mass. This finding poses a significant challenge to existing theoretical models of core formation. We propose that future atmospheric observations with JWST could provide insights into resolving ambiguities of TOI-837 b, offering new perspectives on its composition, formation, and evolution. [ABSTRACT FROM AUTHOR]

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-metre per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m s–1 in Sun-like stars. To overcome this, a multidimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S + LEAF 2 , enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to (a) measure rotation periods for 49 slow rotators in our sample, (b) explore the impact of these results on the spin-down of middle-aged late F, G, and K stars, and (c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations programme, and help fine-tune current stellar structure and evolution models. [ABSTRACT FROM AUTHOR]

The past two decades have seen a major expansion in the availability, size, and precision of time-domain data sets in astronomy. Owing to their unique combination of flexibility, mathematical simplicity, and comparative robustness, Gaussian processes (GPs) have emerged recently as the solution of choice to model stochastic signals in such data sets. In this review, we provide a brief introduction to the emergence of GPs in astronomy, present the underlying mathematical theory, and give practical advice considering the key modeling choices involved in GP regression. We then review applications of GPs to time-domain data sets in the astrophysical literature so far, from exoplanets to active galactic nuclei, showcasing the power and flexibility of the method. We provide worked examples using simulated data, with links to the source code; discuss the problem of computational cost and scalability; and give a snapshot of the current ecosystem of open-source GP software packages. In summary: GP regression is a conceptually simple but statistically principled and powerful tool for the analysis of astronomical time series. It is already widely used in some subfields, such as exoplanets, and gaining traction in many others, such as optical transients. Driven by further algorithmic and conceptual advances, we expect that GPs will continue to be an important tool for robust and interpretable time-domain astronomy for many years to come. [ABSTRACT FROM AUTHOR]

Detecting planetary signatures in radial velocity time-series of young stars is challenging due to their inherently strong stellar activity. However, it is possible to learn information about the properties of the stellar signal by using activity indicators measured from the same stellar spectra used to extract radial velocities. In this manuscript, we present a reanalysis of spectroscopic High Accuracy Radial Velocity Planet Searcher data of the young star K2-233, which hosts three transiting planets. We perform a multidimensional Gaussian process regression on the radial velocity and the activity indicators to characterize the planetary Doppler signals. We demonstrate, for the first time on a real data set, that the use of a multidimensional Gaussian process can boost the precision with which we measure the planetary signals compared to a one-dimensional Gaussian process applied to the radial velocities alone. We measure the semi-amplitudes of K2-233 b, c, and d as |$1.31_{-0.74}^{+0.81}$| , |$1.81_{-0.67}^{+0.71}$| ,  and |$2.72_{-0.70}^{+0.66}$| |${\rm m\, s^{-1}}$| , which translate into planetary masses of |$2.4_{-1.3}^{+1.5}$| , |$4.6_{-1.7}^{+1.8}$| , and |$10.3_{-2.6}^{+2.4}$| M⊕, respectively. These new mass measurements make K2-233 d a valuable target for transmission spectroscopy observations with JWST. K2-233 is the only young system with two detected inner planets below the radius valley and a third outer planet above it. This makes it an excellent target to perform comparative studies, to inform our theories of planet evolution, formation, migration, and atmospheric evolution. [ABSTRACT FROM AUTHOR]

We present a spectroscopic analysis of a 1-yr intensive monitoring campaign of the 22-Myr old planet-hosting M dwarf AU Mic using the HARPS spectrograph. In a companion paper, we reported detections of the planet radial velocity (RV) signatures of the two close-in transiting planets of the system, with respective semi-amplitudes of 5.8 ± 2.5 and 8.5 ± 2.5 m s−1 for AU Mic b and AU Mic c. Here, we perform an independent measurement of the RV semi-amplitude of AU Mic c using Doppler imaging to simultaneously model the activity-induced distortions and the planet-induced shifts in the line profiles. The resulting semi-amplitude of 13.3 ± 4.1 m s−1 for AU Mic c reinforces the idea that the planet features a surprisingly large inner density, in tension with current standard models of core accretion. Our brightness maps feature significantly higher spot coverage and lower level of differential rotation than the brightness maps obtained in late 2019 with the SPIRou spectropolarimeter, suggesting that the stellar magnetic activity has evolved dramatically over a ∼1-yr time span. Additionally, we report a 3σ detection of a modulation at 8.33 ± 0.04 d of the He  i D3 (5875.62 Å) emission flux, close to the 8.46-d orbital period of AU Mic b. The power of this emission (a few 1017 W) is consistent with 3D magnetohydrodynamical simulations of the interaction between stellar wind and the close-in planet if the latter hosts a magnetic field of ∼10 G. Spectropolarimetric observations of the star are needed to firmly elucidate the origin of the observed chromospheric variability. [ABSTRACT FROM AUTHOR]

Observations of exoplanet atmospheres in high resolution have the potential to resolve individual planetary absorption lines, despite the issues associated with ground-based observations. The removal of contaminating stellar and telluric absorption features is one of the most sensitive steps required to reveal the planetary spectrum and, while many different detrending methods exist, it remains difficult to directly compare the performance and efficacy of these methods. Additionally, though the standard cross-correlation method enables robust detection of specific atmospheric species, it only probes for features that are expected a priori. Here, we present a novel methodology using Gaussian process (GP) regression to directly model the components of high-resolution spectra, which partially addresses these issues. We use two archival CRyogenic Infra-Red Echelle Spectrograph (CRIRES)/Very Large Telescope (VLT) data sets as test cases, observations of the hot Jupiters HD 189733 b and 51 Pegasi b, recovering injected signals with average line contrast ratios of ∼4.37 × 10−3 and ∼1.39 × 10−3, and planet radial velocities Δ K p = 1.45 ± 1.53  |$\mathrm{km\, s^{-1}}$| and Δ K p = 0.12 ± 0.12  |$\mathrm{km\, s^{-1}}$| from the injection velocities, respectively. In addition, we demonstrate an application of the GP method to assess the impact of the detrending process on the planetary spectrum, by implementing injection-recovery tests. We show that standard detrending methods used in the literature negatively affect the amplitudes of absorption features in particular, which has the potential to render retrieval analyses inaccurate. Finally, we discuss possible limiting factors for the non-detections using this method, likely to be remedied by higher signal-to-noise data. [ABSTRACT FROM AUTHOR]

We present ground- and space-based photometric observations of TOI-270 (L231-32), a system of three transiting planets consisting of one super-Earth and two sub-Neptunes discovered by TESS around a bright (K-mag = 8.25) M3V dwarf. The planets orbit near low-order mean-motion resonances (5:3 and 2:1) and are thus expected to exhibit large transit timing variations (TTVs). Following an extensive observing campaign using eight different observatories between 2018 and 2020, we now report a clear detection of TTVs for planets c and d, with amplitudes of ∼10 min and a super-period of ∼3 yr, as well as significantly refined estimates of the radii and mean orbital periods of all three planets. Dynamical modelling of the TTVs alone puts strong constraints on the mass ratio of planets c and d and on their eccentricities. When incorporating recently published constraints from radial velocity observations, we obtain masses of |$M_{\mathrm{b}}=1.48\pm 0.18\, M_\oplus$|⁠ , |$M_{\mathrm{c}}=6.20\pm 0.31\, M_\oplus$|⁠ , and |$M_{\mathrm{d}}=4.20\pm 0.16\, M_\oplus$| for planets b, c, and d, respectively. We also detect small but significant eccentricities for all three planets : e b = 0.0167 ± 0.0084, e c = 0.0044 ± 0.0006, and e d = 0.0066 ± 0.0020. Our findings imply an Earth-like rocky composition for the inner planet, and Earth-like cores with an additional He/H2O atmosphere for the outer two. TOI-270 is now one of the best constrained systems of small transiting planets, and it remains an excellent target for atmospheric characterization. [ABSTRACT FROM AUTHOR]

The two most successful methods for exoplanet detection rely on the detection of planetary signals in photometric and radial velocity time-series. This depends on numerical techniques that exploit the synergy between data and theory to estimate planetary, orbital, and/or stellar parameters. In this work, we present a new version of the exoplanet modelling code pyaneti. This new release has a special emphasis on the modelling of stellar signals in radial velocity time-series. The code has a built-in multidimensional Gaussian process approach to modelling radial velocity and activity indicator time-series with different underlying covariance functions. This new version of the code also allows multiband and single transit modelling; it runs on Python 3 , and features overall improvements in performance. We describe the new implementation and provide tests to validate the new routines that have direct application to exoplanet detection and characterization. We have made the code public and freely available at https://github.com/oscaribv/pyaneti. We also present the codes citlalicue and citlalatonac that allow one to create synthetic photometric and spectroscopic time-series, respectively, with planetary and stellar-like signals. [ABSTRACT FROM AUTHOR]

Observational studies of exoplanets are suggestive of a ubiquitous presence of clouds. The current modelling techniques used in emission to account for the clouds tend to require prior knowledge of the cloud condensing species and often do not consider the scattering effects of the cloud. We explore the effects that thermal scattering has on the emission spectra by modelling a suite of hot Jupiter atmospheres with varying cloud single-scattering albedos (SSAs) and temperature profiles. We examine cases ranging from simple isothermal conditions to more complex structures and physically driven cloud modelling. We show that scattering from nightside clouds would lead to brightness temperatures that are cooler than the real atmospheric temperature if scattering is unaccounted for. We show that scattering can produce spectral signatures in the emission spectrum even for isothermal atmospheres. We identify the retrieval degeneracies and biases that arise in the context of simulated JWST spectra when the scattering from the clouds dominates the spectral shape. Finally, we propose a novel method of fitting the SSA spectrum of the cloud in emission retrievals, using a technique that does not require any prior knowledge of the cloud chemical or physical properties. [ABSTRACT FROM AUTHOR]

Binary supermassive black hole (SMBH) systems result from galaxy mergers, and will eventually coalesce due to gravitational wave (GW) emission if the binary separation can be reduced to ≲0.1 pc by other mechanisms. Here, we explore a gravitational self-lensing binary SMBH model for the sharp (duration ∼1 h), quasi-regular X-ray flares – dubbed quasi-periodic eruptions – recently observed from two low-mass active galactic nuclei: GSN 069 and RX J1301.9+2747. In our model, the binary is observed ∼edge-on, such that each SMBH gravitationally lenses light from the accretion disc surrounding the other SMBH twice per orbital period. The model can reproduce the flare spacings if the current eccentricity of RX J1301.9+2747 is ϵ0 ≳ 0.16, implying a merger within ∼1000 yr. However, we cannot reproduce the observed flare profiles with our current calculations. Model flares with the correct amplitude are ∼2/5 the observed duration, and model flares with the correct duration are ∼2/5 the observed amplitude. Our modelling yields three distinct behaviours of self-lensing binary systems that can be searched for in current and future X-ray and optical time-domain surveys: (i) periodic lensing flares, (ii) partial eclipses (caused by occultation of the background mini-disc by the foreground mini-disc), and (iii) partial eclipses with a very sharp in-eclipse lensing flare. Discovery of such features would constitute very strong evidence for the presence of a supermassive binary, and monitoring of the flare spacings will provide a measurement of periastron precession. [ABSTRACT FROM AUTHOR]

We present Mon-735, a detached double-lined eclipsing binary (EB) member of the ∼3 Myr old NGC 2264 star-forming region, detected by Spitzer. We simultaneously model the Spitzer light curves, follow-up Keck/HIRES radial velocities, and the system's spectral energy distribution to determine self-consistent masses, radii, and effective temperatures for both stars. We find that Mon-735 comprises two pre-main-sequence M dwarfs with component masses of M = 0.2918 ± 0.0099 and 0.2661 ± 0.0095 M⊙, radii of R = 0.762 ± 0.022 and 0.748 ± 0.023 R⊙, and effective temperatures of T eff = 3260 ± 73 and 3213 ± 73 K. The two stars travel on circular orbits around their common centre of mass in P = 1.9751388 ± 0.0000050 d. We compare our results for Mon-735, along with another EB in NGC 2264 (CoRoT 223992193), to the predictions of five stellar evolution models. These suggest that the lower mass EB system Mon-735 is older than CoRoT 223992193 in the mass–radius diagram (MRD) and, to a lesser extent, in the Hertzsprung–Russell diagram (HRD). The MRD ages of Mon-735 and CoRoT 223992193 are ∼7–9 and 4–6 Myr, respectively, with the two components in each EB system possessing consistent ages. [ABSTRACT FROM AUTHOR]

Exoplanet emission spectra are often modelled assuming that the hemisphere observed is well represented by a horizontally homogenized atmosphere. However, this approximation will likely fail for planets with a large temperature contrast in the James Webb Space Telescope (JWST) era, potentially leading to erroneous interpretations of spectra. We first develop an analytic formulation to quantify the signal-to-noise ratio and wavelength coverage necessary to disentangle temperature inhomogeneities from a hemispherically averaged spectrum. We find that for a given signal-to-noise ratio, observations at shorter wavelengths are better at detecting the presence of inhomogeneities. We then determine why the presence of an inhomogeneous thermal structure can lead to spurious molecular detections when assuming a fully homogenized planet in the retrieval process. Finally, we quantify more precisely the potential biases by modelling a suite of hot Jupiter spectra, varying the spatial contributions of a hot and a cold region, as would be observed by the different instruments of JWST /NIRSpec. We then retrieve the abundances and temperature profiles from the synthetic observations. We find that in most cases, assuming a homogeneous thermal structure when retrieving the atmospheric chemistry leads to biased results, and spurious molecular detection. Explicitly modelling the data using two profiles avoids these biases, and is statistically supported provided the wavelength coverage is wide enough, and crucially also spanning shorter wavelengths. For the high contrast used here, a single profile with a dilution factor performs as well as the two-profile case, with only one additional parameter compared to the 1D approach. [ABSTRACT FROM AUTHOR]

We present a novel retrieval technique that attempts to model phase curve observations of exoplanets more realistically and reliably, which we call the 2.5-dimensional (2.5D) approach. In our 2.5D approach we retrieve the vertical temperature profile and mean gaseous abundance of a planet at all longitudes and latitudes simultaneously, assuming that the temperature or composition, x , at a particular longitude and latitude (Λ, Φ) is given by |$x(\Lambda ,\Phi) = \bar{x} + (x(\Lambda ,0) - \bar{x})\cos ^n\Phi$|⁠ , where |$\bar{x}$| is the mean of the morning and evening terminator values of x (Λ, 0), and n is an assumed coefficient. We compare our new 2.5D scheme with the more traditional 1D approach, which assumes the same temperature profile and gaseous abundances at all points on the visible disc of a planet for each individual phase observation, using a set of synthetic phase curves generated from a GCM-based simulation. We find that our 2.5D model fits these data more realistically than the 1D approach, confining the hotter regions of the planet more closely to the dayside. We then apply both models to WASP-43b phase curve observations of HST /WFC3 and Spitzer /IRAC. We find that the dayside of WASP-43b is apparently much hotter than the nightside and show that this could be explained by the presence of a thick cloud on the nightside with a cloud top at pressure <0.2 bar. We further show that while the mole fraction of water vapour is reasonably well constrained to (1–10) × 10−4, the abundance of CO is very difficult to constrain with these data since it is degenerate with temperature and prone to possible systematic radiometric differences between the HST /WFC3 and Spitzer /IRAC observations. Hence, it is difficult to reliably constrain C/O. [ABSTRACT FROM AUTHOR]

Together with the stellar rotation, the spotted surfaces of low-mass magnetically active stars produce modulations in their brightness. These modulations can be resolved by photometric variability surveys, allowing direct measurements of stellar spin rates. In this proceedings, we present results of a multisite photometric survey dedicated to the measurement of spin rates in the 30 Myr cluster NGC 3766. Inside the framework of the Monitor Project, the cluster was monitored during 2014 in the i-band by the Wide Field Imager at the MPG/ESO 2.2-m telescope. Data from Gaia-DR2 and grizY photometry from DECam/CTIO were used to identify cluster members. We present spin rates measured for ⁓200 cluster members. [ABSTRACT FROM AUTHOR]

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasiperiodic -- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalizing over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine- fitting periodogram and an autocorrelation function method.We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these and many more, altogether 1102 Kepler objects of interest, and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterizing star--planet interactions. The code used to implement this method is available online. [ABSTRACT FROM AUTHOR]

KMOS (K-Band Multi Object Spectrograph) is a novel integral field spectrograph installed in the Very Large Telescope's (VLT's) ANTU unit. The instrument offers an ability to observe 24 2.8 arcsec x 2.8 arcsec subfields positionable within a 7.2 arcmin patrol field, each subfield producing a spectrum with a 14 x 14-pixel spatial resolution. The main science drivers for KMOS are the study of galaxies, star formation, and molecular clouds, but its ability to simultaneously measure spectra of multiple stars makes KMOS an interesting instrument for exoplanet atmosphere characterization via transmission spectroscopy. We set to test whether transmission spectroscopy is practical with KMOS, and what are the conditions required to achieve the photometric precision needed, based on observations of a partial transit of WASP- 19b, and full transits of GJ 1214b and HD 209458b. Our analysis uses the simultaneously observed comparison stars to reduce the effects from instrumental and atmospheric sources, and Gaussian processes to model the residual systematics. We show that KMOS can, in theory, deliver the photometric precision required for transmission spectroscopy. However, this is shown to require (a) pre-imaging to ensure accurate centring and (b) a very stable night with optimal observing conditions (seeing ~0.8 arcsec). Combining these two factors with the need to observe several transits, each with a sufficient out-of-transit baseline (and with the fact that similar or better precision can be reached with telescopes and instruments with smaller pressure), we conclude that transmission spectroscopy is not the optimal science case to take advantage of the abilities offered by KMOS and VLT. [ABSTRACT FROM AUTHOR]

We present an analysis of Spitzer/Infrared Array Camera primary transit and secondary eclipse light curves measured for HD 209458b, using Gaussian process models to marginalize over the intrapixel sensitivity variations in the 3.6 and 4.5 µm channels and the ramp effect in the 5.8 and 8.0 µm channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not confirm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better fit. We find that a solar-abundance clear-atmosphere model without a thermal inversion underpredicts the measured emission in the 4.5 µm channel, which may suggest the atmosphere is depleted in carbon monoxide. An acceptable fit to the emission data can be achieved by assuming that the planet radiates as an isothermal blackbody with a temperature of 1484 ± 18 K. [ABSTRACT FROM AUTHOR]

Among the available methods for dating stars, gyrochronology is a powerful one because it requires knowledge of only the star's mass and rotation period. Gyrochronology relations have previously been calibrated using young clusters, with the Sun providing the only age dependence, and are therefore poorly calibrated at late ages. We used rotation period measurements of 310 Kepler stars with asteroseismic ages, 50 stars from the Hyades and Coma Berenices clusters and 6 field stars (including the Sun) with precise age measurements to calibrate the gyrochronology relation, whilst fully accounting for measurement uncertainties in all observable quantities. We calibrated a relation of the form P = An x (B -- V -- c)b, where P is rotation period in days, A is age in Myr, B and V are magnitudes and a, b and n are the free parameters of our model. We found a = 0.40-0.05+0.3, b = -0.02+0.05 and n = -0.09+0.02 Markov Chain Monte Carlo methods were used to explore the posterior probability distribution functions of the gyrochronology parameters and we carefully checked the effects of leaving out parts of our sample, leading us to find that no single relation between rotation period, colour and age can adequately describe all the subsets of our data. The Kepler asteroseismic stars, cluster stars and local field stars cannot all be described by the same gyrochronology relation. The Kepler asteroseismic stars may be subject to observational biases; however, the clusters show unexpected deviations from the predicted behaviour, providing concerns for the overall reliability of gyrochronology as a dating method. [ABSTRACT FROM AUTHOR]

CoRoT-7 is an magnitude K-star whose light curve shows transits with a depth of 0.3 mmag and a period of 0.854 d, superimposed on variability at the 1 per cent level, due to the modulation of evolving active regions with the star's 23-d rotation period. In this paper, we revisit the published HARPS radial-velocity (RV) measurements of the object, which were previously used to estimate the companion mass, but have been the subject of ongoing debate. We build a realistic model of the star's activity during the HARPS observations, by fitting simultaneously the linewidth (as measured by the width of the cross-correlation function) and the line bisector, and use it to evaluate the contribution of activity to the RV variations. The data show clear evidence of errors above the level of the formal uncertainties, which are accounted for neither by activity nor by any plausible planet model and which increase rapidly with a decreasing signal-to-noise ratio (S/N) of the spectra. We cite evidence of similar systematics in mid-S/N spectra of other targets obtained with HARPS and other high-precision RV spectrographs, and discuss possible sources. Allowing for these, we re-evaluate the semi-amplitude of the CoRoT-7b signal, finding m s, a tentative detection with a much reduced significance (1.2σ) compared to previous estimates. We also argue that the combined presence of activity and additional errors precludes a meaningful search for additional low-mass companions, despite previous claims to the contrary. Taken at face value, our analysis points to a lower density for CoRoT-7b, the 1σ mass range spanning and allowing for a wide range of bulk compositions. In particular, an ice-rich composition is compatible with the RV constraints. More generally, this study highlights the importance of a realistic treatment of both activity and uncertainties, particularly in the medium S/N regime, which applies to most small planet candidates from CoRoT and Kepler. [ABSTRACT FROM AUTHOR]

Surveys for exoplanetary transits are usually limited not by photon noise but rather by the amount of red noise in their data. In particular, although the CoRoT space-based survey data are being carefully scrutinized, significant new sources of systematic noises are still being discovered. Recently, a magnitude-dependant systematic effect was discovered in the CoRoT data by Mazeh et al. and a phenomenological correction was proposed. Here we tie the observed effect to a particular type of effect, and in the process generalize the popular Sysrem algorithm to include external parameters in a simultaneous solution with the unknown effects. We show that a post-processing scheme based on this algorithm performs well and indeed allows for the detection of new transit-like signals that were not previously detected. [ABSTRACT FROM AUTHOR]

We report on the results of a time-series photometric survey of M50 (NGC 2323), a open cluster, carried out using the Cerro Tololo Inter-American Observatory (CTIO) 4-m Blanco telescope and Mosaic-II detector as part of the Monitor project. Rotation periods were derived for 812 candidate cluster members over the mass range . The rotation period distributions show a clear mass-dependent morphology, statistically indistinguishable from those in NGC 2516 and M35 taken from the literature. Due to the availability of data from three observing runs separated by ∼10 and 1 month time-scales, we are able to demonstrate clear evidence for evolution of the photometric amplitudes, and hence spot patterns, over the 10 month gap. We are not able to constrain the time-scales for these effects in detail due to limitations imposed by the large gaps in our sampling, which also prevent the use of the phase information. [ABSTRACT FROM AUTHOR]

We present the results of a systematic search for transiting planets in a ∼5 Myr open cluster, NGC 2362. We observed ∼1200 candidate cluster members, of which ∼475 are believed to be genuine cluster members, for a total of ∼100 h. We identify 15 light curves with reductions in flux that pass all our detection criteria, and six of the candidates have occultation depths compatible with a planetary companion. The variability in these six light curves would require very large planets to reproduce the observed transit depth. If we assume that none of our candidates are, in fact, planets then we can place upper limits on the fraction of stars with hot Jupiters (HJs) in NGC 2362. We obtain 99 per cent confidence upper limits of 0.22 and 0.70 on the fraction of stars with HJs ( fp) for 1–3 and 3–10 d orbits, respectively, assuming all HJs have a planetary radius of 1.5 RJup. These upper limits represent observational constraints on the number of stars with HJs at an age ≲10 Myr, when the vast majority of stars are thought to have lost their protoplanetary discs. Finally, we extend our results to the entire Monitor project, a survey searching young, open clusters for planetary transits, and find that the survey as currently designed should be capable of placing upper limits on fp near the observed values of fp in the solar neighbourhood. [ABSTRACT FROM AUTHOR]

With more than 80 transits observed in the CoRoT lightcurve with a cadence of 32 s, CoRoT-Exo-2b provides an excellent case to search for the secondary eclipse of the planet, with an expected signal of less than 10−4 in relative flux. The activity of the star causes a modulation on the flux that makes the detection of this signal challenging. We describe the technique used to seek for the secondary eclipse, that leads to a tentative 2.5σ detection of a 5.5×10−5 eclipse. If the effect of the spots are not taken into account, the times of transit centers will also be affected. They could lead to an erroneous detection of periodic transit timing variations of ~20 s and with a 7.45 d period. By measuring the transit central times at different depths of the transit (transit bisectors), we show that there are no such periodic variations in the CoRoT-Exo-2b O-C residuals larger than ~10 s. [ABSTRACT FROM PUBLISHER]

We report on the results of a time-series photometric survey of NGC 2362, carried out using the CTIO 4-m Blanco telescope and Mosaic-II detector as part of the Monitor project. Rotation periods were derived for 271 candidate cluster members over the mass range . The rotation period distributions show a clear mass-dependent morphology, qualitatively similar to that in NGC 2264, as would be expected from the age of this cluster. Using models of angular momentum evolution, we show that angular momentum losses over the ∼1–5 Myr age range appear to be needed in order to reproduce the evolution of the slowest rotators in the sample from the ONC to NGC 2362, as found by many previous studies. By incorporating Spitzer IRAC mid-infrared (mid-IR) measurements, we found that three to four objects showing mid-IR excesses indicative of the presence of circumstellar discs were all slow rotators, as would be expected in the disc regulation paradigm for early pre-main-sequence angular momentum evolution, but this result is not statistically significant at present, given the extremely limited sample size. [ABSTRACT FROM AUTHOR]

We report on the results of an I-band time-series photometric survey of NGC 2547 using the MPG/ESO 2.2-m telescope with Wide Field Imager, achieving better than 1 per cent photometric precision per data point over . Candidate cluster members were selected from a V versus V− I colour–magnitude diagram over (covering masses from down to below the brown dwarf limit), finding 800 candidates, of which we expect ∼330 to be real cluster members, taking into account contamination from the field (which is most severe at the extremes of our mass range). Searching for periodic variations in these gave 176 detections over the mass range . The rotation period distributions were found to show a clear mass-dependent morphology, qualitatively intermediate between the distributions obtained from similar surveys in NGC 2362 and 2516, as would be expected from the age of this cluster. Models of the rotational evolution were investigated, finding that the evolution from NGC 2362 to 2547 was qualitatively reproduced (given the uncertainty in the age of NGC 2547) by solid body and core-envelope decoupled models from our earlier NGC 2516 study without need for significant modification. [ABSTRACT FROM AUTHOR]

Intrinsic stellar variability can hinder the detection of shallow transits, particularly in space-based data. Therefore, this variability has to be filtered out before running the transit search. Unfortunately, filtering out the low frequency signal of the stellar variability also modifies the transit shape. This results in errors in the measured transit depth and duration used to derive the planet radius, and orbital inclination. We present an evaluation of the magnitude of this effect based on 20 simulated light curves from the CoRoT blind exercise 2 (BT2). We then present an iterative filter which uses the strictly periodic nature of the transits to separate them from other forms of variability, so as to recover the original transit shape before deriving the planet parameters. On average with this filter, we improve the estimation of the transit depth and duration by 15% and 10% respectively. [ABSTRACT FROM PUBLISHER]

We report the discovery of a low-mass pre-main-sequence (PMS) eclipsing binary (EB) with a 5.3 d orbital period. JW 380 was detected as part of a high-cadence time-resolved photometric survey (the Monitor project) using the 2.5-m Isaac Newton Telescope and Wide Field Camera for a survey of a single field in the Orion nebula cluster (ONC) region in V and i bands. The star is assigned a 99 per cent membership probability from proper motion measurements, and radial velocity observations indicate a systemic velocity within 1σ of that of the ONC. Modelling of the combined light and radial velocity curves of the system gave stellar radii of and for the primary and the secondary, with a significant third light contribution which is also visible as a third peak in the cross-correlation functions used to derive radial velocities. The masses and radii appear to be consistent with stellar models for 2–3 Myr age from several authors, within the present observational errors. These observations probe an important region of mass–radius parameter space, where there are currently only a handful of known PMS EB systems with precise measurements available in the literature. [ABSTRACT FROM AUTHOR]

We present an extension of the formalism recently proposed by Pepper and Gaudi to evaluate the yield of transit surveys in homogeneous stellar systems, incorporating the impact of correlated noise on transit time-scales on the detectability of transits, and simultaneously incorporating the magnitude limits imposed by the need for radial velocity (RV) follow-up of transit candidates. New expressions are derived for the different contributions to the noise budget on transit time-scales and the least-squares detection statistic for box-shaped transits, and their behaviour as a function of stellar mass is re-examined. Correlated noise that is constant with apparent stellar magnitude implies a steep decrease in detection probability at the high-mass end which, when considered jointly with the RV requirements, can severely limit the potential of otherwise promising surveys in star clusters. However, we find that small-aperture, wide-field surveys may detect hot Neptunes whose RV signal can be measured with present-day instrumentation in very nearby (<100 pc) clusters. [ABSTRACT FROM AUTHOR]

We report on the results of an i-band time-series photometric survey of NGC 2516 using the Cerro Tololo Inter-American Observatory (CTIO) 4-m Blanco telescope and 8k Mosaic-II detector, achieving better than 1 per cent photometric precision per data point over . Candidate cluster members were selected from a V versus colour–magnitude diagram over (covering masses from down to below the brown dwarf limit), finding 1685 candidates, of which we expect ∼1000 to be real cluster members, taking into account contamination from the field (which is most severe at the extremes of our mass range). Searching for periodic variations in these gave 362 detections over the mass range . The rotation period distributions were found to show a remarkable morphology as a function of mass, with the fastest rotators bounded by , and the slowest rotators for bounded by a line of , with those for following a flatter relation closer to . Models of the rotational evolution were investigated, finding that the evolution of the fastest rotators was well reproduced by a conventional solid body model with a mass-dependent saturation velocity, whereas core–envelope decoupling was needed to reproduce the evolution of the slowest rotators. None of our models were able to simultaneously reproduce the behaviour of both populations. [ABSTRACT FROM AUTHOR]

We have begun a large-scale photometric survey of nearby open clusters and star-forming regions, the Monitor project, aiming to measure time-series photometry for >10 000 cluster members over >10 deg2 of sky, to find low-mass eclipsing binary and planet systems. We describe the software pipeline we have developed for this project, showing that we can achieve peak rms accuracy over the entire data set of better than ∼2 mmag using aperture photometry, with rms <1 per cent over ∼4 mag, in data from 2- and 4-m class telescopes with wide-field mosaic cameras. We investigate the noise properties of our data, finding correlated ‘red’ noise at the ∼1–1.5 mmag level in bright stars, over transit-like time-scales of 2.5 h. An important source of correlated noise in aperture photometry is image blending, which produces variations correlated with the seeing. We present a simple blend index based on fitting polynomials to these variations, and find that subtracting the fit from the data provides a method to reduce their amplitude, in lieu of using techniques, such as point spread function fitting photometry, which tackle their cause. Finally, we use thesysrem algorithm to search for any further systematic effects. [ABSTRACT FROM AUTHOR]

We report on the results of a V- and i-band time-series photometric survey of M34 (NGC 1039) using the Wide Field Camera (WFC) on the Isaac Newton Telescope (INT), achieving better than 1 per cent precision per data point for . Candidate cluster members were selected from a V versus V− I colour–magnitude diagram over , finding 714 candidates, of which we expect ∼400 to be real cluster members (taking into account contamination from the field). The mass function was computed, and found to be consistent with a lognormal distribution in d N/d log M. Searching for periodic variable objects in the candidate members gave 105 detections over the mass range . The distribution of rotation periods for was found to peak at ∼7 d, with a tail of fast rotators down to periods of . For we found a peak at short periods, with a lack of slow rotators (e.g. ), consistent with the work of other authors at very low masses. Our results are interpreted in the context of previous work, finding that we reproduce the same general features in the rotational period distributions. A number of rapid rotators were found with velocities ∼ a factor of 2 lower than in the Pleiades, consistent with models of angular momentum evolution assuming solid body rotation without needing to invoke core–envelope decoupling. [ABSTRACT FROM AUTHOR]

Spectroscopic observations of transiting exoplanets have provided the first indications of their atmospheric structure and composition. Optimal estimation retrievals have been successfully applied to solar system planets to determine the temperature, composition and aerosol properties of their atmospheres, and have recently been applied to exoplanets. We show the effectiveness of the technique when combined with simulated observations from the proposed space telescope EChO, and also discuss the difficulty of constraining a complex system with sparse data and large uncertainties, using the super-Earth GJ 1214b as an example. [ABSTRACT FROM PUBLISHER]

We investigate the variability properties of main-sequence stars in the first month of Kepler data, using a new astrophysically robust systematics correction. We find that 36% appear more variable than the Sun, and confirm the trend of increasing variability with decreasing effective temperature. We define low- and high-variability samples, with a cut at twice the level of the active Sun, and compare properties of the stars belonging to each sample. We find tentative evidence that the more active stars have lower proper motions. The frequency content of the variability shows clear evidence for periodic or quasi-periodic behaviour in 16% of stars, and highlights significant differences in the nature of variability between spectral types. Most A and F stars have short periods (< 2 days) and highly sinusoidal variability, suggestive of pulsations, whilst G, K and M stars tend to have longer periods (> 5 days, with a trend towards longer periods at later spectral types) and show a mixture of periodic and stochastic variability, indicative of activity. Finally, we use autoregressive models to characterise the stochastic component of the variability, and show that its typical amplitude and time-scale increase towards later spectral types, which we interpret as an increase in the characteristic size and life-time of active regions. Full details will be published shortly. [ABSTRACT FROM PUBLISHER]

We present a re-analysis of archival HST/NICMOS transmission spectroscopy of the exoplanet system, HD 189733, from which detections of several molecules have been claimed. As expected, we can replicate the transmission spectrum previously published when we use an identical model for the systematic effects, although the uncertainties are larger as we use a residual permutation algorithm in an effort to account for instrumental systematics. We also find that the transmission spectrum is considerably altered when slightly changing the instrument model, and conclude that the NICMOS transmission spectrum is too dependent on the method used to remove systematics to be considered a robust detection of molecular species, given that there is no physical reason to believe that the baseline flux should be modelled as a linear function of any chosen set of parameters. [ABSTRACT FROM PUBLISHER]

Virtually all exoplanet detection and characterisation methods are based on time-domain data. This invited talk gave an overview of some recent results in the field, highlighting some of the time-series-specific challenges encountered along the way. In particular it focussed on planetary transits: how to detect shallow, rare transits in noisy data, and how to model them with extreme accuracy to extract information about the transiting planet's atmosphere. Space-based transit surveys also constitute an extraordinary goldmine of information on stellar variability, and the talk touched briefly upon some recent statistical work in that field. [ABSTRACT FROM PUBLISHER]