Star-planet Systems Research Articles

ATMOSPHERE EXPANSION AND MASS LOSS OF CLOSE-ORBIT GIANT EXOPLANETS HEATED BY STELLAR XUV. II. EFFECTS OF PLANETARY MAGNETIC FIELD; STRUCTURING OF INNER MAGNETOSPHERE

This is the second paper in a series where we build a self-consistent model to simulate the mass-loss process of a close-orbit magnetized giant exoplanet, so-called hot Jupiter (HJ). In this paper we generalize the hydrodynamic (HD) model of an HJ expanding hydrogen atmosphere, proposed in the first paper, to include the effects of intrinsic planetary magnetic field. The proposed self-consistent axisymmetric 2D magnetohydrodynamics model incorporates radiative heating and ionization of the atmospheric gas, basic hydrogen chemistry for the appropriate account of major species composing HJ's upper atmosphere and related radiative energy deposition, and H3+ and Ly{\alpha} cooling processes. The model also takes into account a realistic solar-type X-ray/EUV spectrum for calculation of intensity and column density distribution of the radiative energy input, as well as gravitational and rotational forces acting in a tidally locked planet-star system. An interaction between the expanding atmospheric plasma and an intrinsic planetary magnetic dipole field leads to the formation of a current-carrying magnetodisk that plays an important role for topology and scaling of the planetary magnetosphere. A cyclic character of the magnetodisk behavior, composed of consequent phases of the disk formation followed by the magnetic reconnection with the ejection of a ring-type plasmoid, has been discovered and investigated. We found that the mass-loss rate of an HD 209458b analog planet is weakly affected by the equatorial surface field <0.3 G, but is suppressed by an order of magnitude at the field of 1 G.

Magnetic energy fluxes in close-in star-planet systems

AbstractMagnetic interactions between a close-in planet and its host star have been postulated to be a source of enhanced chromospheric emissions. We develop three dimensional global models of star-planet systems under the ideal magnetohydrodynamic (MHD) approximation to explore the impact of magnetic topology on the energy fluxes induced by the magnetic interaction. We conduct twin numerical experiments in which only the magnetic topology of the interaction is altered. We find that the Poynting flux varies by more than an order of magnitude when varying the magnetic topology from an aligned case to an anti-aligned case. This provides a simple and robust physical explanation for on/off enhanced chromospheric emissions induced by a close-in planet on time-scales of the order of days to years.

Planet migration and magnetic torques

AbstractThe possibility that magnetic torques may participate in close-in planet migration has recently been postulated. We develop three dimensional global models of magnetic star-planet interaction under the ideal magnetohydrodynamic (MHD) approximation to explore the impact of magnetic topology on the development of magnetic torques. We conduct twin numerical experiments in which only the magnetic topology of the interaction is altered. We find that magnetic torques can vary by roughly an order of magnitude when varying the magnetic topology from an aligned case to an anti-aligned case. Provided that the stellar magnetic field is strong enough, we find that magnetic migration time scales can be as fast as ~100 Myr. Hence, our model supports the idea that magnetic torques may participate in planet migration for some close-in star-planet systems.

Corot 310266512: A Light Curve With Primary, Secondary And Tertiary Eclipses

We present the photometric study of an interesting target in the CoRoT exoplanet database: CoRoT 310266512. Its light curve shows primary, secondary and tertiary eclipses that suggests the presence of at least three celestial bodies. The primary and secondary eclipses have the same orbital period, 7.42 days, and the tertiary eclipse has an orbital period of 3.27 days. Two of the tertiary eclipses fall within a primary eclipse and a secondary eclipse. The properties of the light curve indicate the presence of two physically separated systems. The primary and secondary eclipses corresponds to a binary system (System I). The tertiary eclipses correspond to a star-planet system or a star-dwarf system (System II). Some parameters of these two systems are obtained from JKTEBOP [1] program.

Tidal dissipation in stars and fluid planetary layers and its impact on the evolution of star-planet systems

Tidal dissipation in stars and planets is one of the key physical mechanisms that drive the evolution of planetary systems. It intrinsically depends on the nature of the tidal response of celestial bodies, which is directly linked to their internal structure and friction. Indeed, it is highly resonant in the case of fluids. In this work, we present a local analytical modeling of tidal gravito-inertial waves, which can be excited in stars and fluid planetary layers. This model allows us to understand the properties of their resonant dissipation as a function of the excitation frequencies, the rotation, the stratification, and the viscous and thermal properties of the studied fluid regions. Next, we introduce such a complex tidal dissipation frequency-spectra in a celestial mechanics numerical code to give a qualitative overview of its impact on the evolution of planetary systems. We consider the example of a two-body coplanar system with a punctual perturber orbiting a central fluid body. We demonstrate how the viscous dissipation of tidal waves can lead to a strongly erratic orbital evolution. Finally, we characterize such a non-regular dynamics as a function of the properties of resonances, which have been determined thanks to our local fluid model.

TRANSIENT CHAOS AND FRACTAL STRUCTURES IN PLANETARY FEEDING ZONES

The circular restricted three body problem is investigated in the context of accretion and scattering processes. In our model a large number of identical non-interacting mass-less planetesimals are considered in planar case orbiting a star-planet system. This description allows us to investigate in dynamical systems approach the gravitational scattering and possible captures of the particles by the forming planetary embryo. Although the problem serves a large variety of complex motion, the results can be easily interpreted because of the low dimensionality of the phase space. We show that initial conditions define isolated regions of the disk, where accretion or escape of the planetesimals occur, these have, in fact, a fractal structure. The fractal geometry of these "basins" implies that the dynamics is very complex. Based on the calculated escape rates and escape times, it is also demonstrated that the planetary accretion rate is exponential for short times and follows a power-law for longer integration. A new numerical calculation of the maximum mass that a planet can reach (described by the expression of the isolation mass) is also derived.

The pulsating sdB+M eclipsing system NY Virginis and its circumbinary planets

We searched for circumbinary planets orbiting NY Vir in historical eclipse times including our long-term CCD data. Sixty-eight times of minimum light with accuracies better than 10 s were used for the ephemeris computations. The best fit to those timings indicated that the orbital period of NY Vir has varied due to a combination of two sinusoids with periods of $P_3$=8.2 yr and $P_4$=27.0 yr and semi-amplitudes of $K_3$=6.9 s and $K_4$=27.3 s, respectively. The periodic variations most likely arise from a pair of light-time effects due to the presence of third and fourth bodies that are gravitationally bound to the eclipsing pair. We have derived the orbital parameters and the minimum masses, $M_3 \sin i_3$ = 2.8 M$\rm_{Jup}$ and $M_4 \sin i_4$ = 4.5 M$\rm_{Jup}$, of both objects. A dynamical analysis suggests that the outer companion is less likely to orbit the binary on a circular orbit. Instead we show that future timing data might push its eccentricity to moderate values for which the system exhibits long-term stability. The results demonstrate that NY Vir is probably a star-planet system, which consists of a very close binary star and two giant planets. The period ratio $P_{3}/P_{4}$ suggests that a long-term gravitational interaction between them would result in capture into a nearly 3:10 mean motion resonance. When the presence of the circumbinary planets is verified and understood more comprehensively, the formation and evolution of this planetary system should be advanced greatly.

Asteroseismic inference on the spin-orbit misalignment and stellar parameters of HAT-P-7

The measurement of obliquities in star-planet systems is of great importance for the understanding of planet system formation and evolution. The bright and well studied HAT-P-7 system is intriguing as several Rossiter-McLaughlin (RM) measurements found a large projected obliquity in this system, but it was so far not possible to determine if the orbit is polar and/or retrograde. The goal of this study is to measure the stellar inclination and hereby the full 3D obliquity of the HAT-P-7 system instead of only the 2D projection as measured by the RM effect. In addition we provide an updated set of stellar parameters for the star. We use the full set of available observations from Kepler spanning Q0-Q17 to produce the power spectrum of HAT-P-7. We extract oscillation mode frequencies via an MCMC peak-bagging routine, and use the results from this to estimate the stellar inclination angle. Combining this with the projected obliquity from RM and the inclination of the orbital plane allows us to determine the stellar obliquity. We use asteroseismology to model the star from the extracted frequencies using two different approaches to the modelling where either the MESA or the GARSTEC stellar evolution codes are adopted. Using our updated asteroseismic modelling we find, i.a., the following stellar parameters for HAT-P-7: M=1.51{+0.04}{-0.05}Msun, $R=2.00{+0.01}{-0.02}Rsun, and age = 2.07{+0.28}{-0.23} Gyr. Our asteroseismic modelling offers a high precision on the stellar parameters, for instance is the uncertainty on age of the order ~11%. For the stellar inclination we estimate i_star<36.5 deg., which translates to an obliquity between 83 and 111 deg. We find that the planet HAT-P-7b is likely retrograde in its orbit, and that the orbit is close to being polar. The new parameters for the star gives an updated planetary density of 0.65+-0.03 g cm^{-3}, which is lower than previous estimates.

THE EFFECTS OF REFRACTION ON TRANSIT TRANSMISSION SPECTROSCOPY: APPLICATION TO EARTH-LIKE EXOPLANETS

We quantify the effects of refraction in transit transmission spectroscopy on spectral absorption features and on temporal variations that could be used to obtain altitude-dependent spectra for planets orbiting stars of different stellar types. We validate our model against altitude-dependent transmission spectra of the Earth from ATMOS and against lunar eclipse spectra from Palle et al. (2009). We perform detectability studies to show the potential effects of refraction on hypothetical observations of Earth analogs with the James Webb Space Telescope (JWST) Near-Infrared Spectrograph (NIRSPEC). Due to refraction, there will be a maximum tangent pressure level that can be probed during transit for each given planet-star system. We show that because of refraction, for an Earth-analog planet orbiting in the habitable zone of a Sun-like star only the top 0.3 bars of the atmosphere can be probed, leading to a decrease in the signal to noise ratio (SNR) of absorption features by 60%, while for an Earth-analog planet orbiting in the habitable zone of an M5V star it is possible to probe almost the entire atmosphere with minimal decreases in SNR. We also show that refraction can result in temporal variations in the transit transmission spectrum which may provide a way to obtain altitude-dependent spectra of exoplanet atmospheres. Additionally, the variations prior to ingress and subsequent to egress provide a way to probe pressures greater than the maximum tangent pressure that can be probed during transit. Therefore, probing the maximum range of atmospheric altitudes, and in particular the near-surface environment of an Earth-analog exoplanet, will require looking at out-of-transit refracted light in addition to the in-transit spectrum.

Analysis and interpretation of 15 quarters ofKeplerdata of the disintegrating planet KIC 12557548 b

The Kepler object KIC 12557548 shows irregular eclipsing behaviour with a constant 15.685 hr period, but strongly varying transit depth. In this paper we fit individual eclipses, in addition to fitting binned light curves, to learn more about the process underlying the eclipse depth variation. Additionally, we put forward observational constraints that any model of this planet-star system will have to match. We find two quiescent spells of ~30 orbital periods each where the transit depth is <0.1%, followed by relatively deep transits. Additionally, we find periods of on-off behaviour where >0.5% deep transits are followed by apparently no transit at all. Apart from these isolated events we find neither significant correlation between consecutive transit depths nor a correlation between transit depth and stellar intensity. We find a three-sigma upper limit for the secondary eclipse of 4.9*10^-5, consistent with a planet candidate with a radius of less than 4600 km. Using the short cadence data we find that a 1-D exponential dust tail model is insufficient to explain the data. We improved our model to a 2-D, two-component dust model with an opaque core and an exponential tail. Using this model we fit individual eclipses observed in short cadence mode. We find an improved fit of the data, quantifying earlier suggestions by Budaj (2013) of the necessity of at least two components. We find that deep transits have most absorption in the tail, and not in a disk-shaped, opaque coma, but the transit depth and the total absorption show no correlation with the tail length.

Planet packing in circumbinary systems

The recent discovery of planets orbiting main sequence binaries will provide crucial constraints for theories of binary and planet formation. The formation pathway for these planets is complicated by uncertainties in the formation mechanism of the host stars. In this paper, we compare the dynamical states of single and binary star planetary systems. Specifically, we pose two questions: (1) What does it mean for a circumbinary system to be dynamically packed? (2) How many systems are required to differentiate between a population of packed or sparse planets? We determine when circumbinary systems become dynamically unstable as a function of the separation between the host-stars and the inner planet, and the first and second planets. We show that these represent unique stability constraints compared to single-star systems. We find that although the existing Kepler data is insufficient to distinguish between a population of packed or sparse circumbinary systems, a more thorough study of circumbinary TTVs combined with an order of magnitude increase in the number of systems may prove conclusive. Future space missions such as TESS provide the best opportunity for increasing the sample size.

Eccentricity Dependence on Iron Abundance

AbstractThe occurrence and eccentricity distribution of planets as a function of period is significantly different for iron-rich and iron-poor planet systems. We find that iron-poor stars with planets having periods between 525 and 600 days have higher eccentricity than such systems outside this range. If whole planet pollution causes the correlation of giant planet eccentricity with stellar iron abundance, then this cluster could be due to a paucity of pollution in this period range. Newly reported patterns of planet occurrence must result from planet system architectural features such as the snow line, followed by subsequent migration. Different results favor pollution or higher initial iron abundance causing the higher occurrence fraction of giant planets hosted by iron-rich stars, but the two explanations could be complementary. Relations between planet and stellar parameters are a major product of planet-finding, which promise further insights into star-planet system formation and evolution. Collaborators are sought to study these patterns. We expect a spirited debate over the relative contributions of initial abundances, disk accretion, and whole planet accretion.

Photospheric activity, rotation, and magnetic interaction in LHS 6343 A

Context. The Kepler mission has recently discovered a brown dwarf companion transiting one member of the M4V+M5V visual binary system LHS 6343 AB with an orbital period of 12.71 days. Aims. The particular interest of this transiting system lies in the synchronicity between the transits of the brown dwarf C component and the main modulation observed in the light curve, which is assumed to be caused by rotating starspots on the A component. We model the activity of this star by deriving maps of the active regions that allow us to study stellar rotation and the possible interaction with the brown dwarf companion. Methods. An average transit profile was derived, and the photometric perturbations due to spots occulted during transits are removed to derive more precise transit parameters. We applied a maximum entropy spot model to fit the out-of-transit optical modulation as observed by Kepler during an uninterrupted interval of 500 days. It assumes that stellar active regions consist of cool spots and bright faculae whose visibility is modulated by stellar rotation. Results. Thanks to the extended photometric time series, we refine the determination of the transit parameters and find evidence of spots that are occulted by the brown dwarf during its transits. The modelling of the out-of-transit light curve of LHS 6343 A reveals several starspots rotating with a slightly longer period than the orbital period of the brown dwarf, i.e., 13.13 +- 0.02 days. No signature attributable to differential rotation is observed. We find evidence of a persistent active longitude on the M dwarf preceding the sub- companion point by 100 deg and lasting for at least 500 days. This can be relevant for understanding how magnetic interaction works in low-mass binary and star-planet systems.

GJ 1214 reviewed

GJ 1214 is orbited by a transiting super-Earth-mass planet. It is a primary target for ongoing efforts to understand the emerging population of super-Earth-mass planets around M dwarfs. We present new precision astrometric measurements, a re-analysis of HARPS radial velocity measurements, and medium-resolution infrared spectroscopy of GJ 1214. We combine these measurements with recent transit follow-up observations and new catalog photometry to provide a comprehensive update of the star-planet properties. The distance is obtained with 0.6% relative uncertainty using CAPScam astrometry. The new value increases the nominal distance to the star by ~10% and is significantly more precise than previous measurements. Updated Doppler measurements combined with published transit observations significantly refine the constraints on the orbital solution. The analysis of the infrared spectrum and photometry confirm that the star is enriched in metals compared to the Sun. Using all this information, combined with empirical mass-luminosity relations for low mass stars, we derive updated values for the bulk properties of the star-planet system. We also use infrared absolute fluxes to estimate the stellar radius and to re-derive the star-planet properties. Both approaches provide very consistent values for the system. Our analysis shows indicates that the favoured mean density of GJ 1214b is 1.6 +/-0.6 g cm^{-3}. We illustrate how fundamental properties of M dwarfs are of paramount importance in the proper characterization of the low mass planetary candidates orbiting them. Given that the distance is now known to better than 1%, interferometric measurements of the stellar radius, additional high precision Doppler observations, and/or or detection of the secondary transit (occultation), are necessary to further improve the constraints on the GJ 1214 star-planet properties.

Internal waves and tides in star-planet systems

Internal waves and tidal interactions constitute powerful mechanisms for angular momentum exchanges in star-planet systems. In this lecture, we review the state of the art of their modelling and we discuss their impact on the evolution of stellar systems. In this context, we first show how internal waves are able to deeply modify the secular transport of angular momentum in stellar interiors in the whole Hertzsprung-Russelll diagram. Then, we show how tidal waves, as well as the related exchanges of angular momentum, strongly depend on the internal structure of stars and planets.

WTS: A near-infrared transit survey

The WFCAM Transit Survey is a transiting planet survey running on the United Kingdom Infrared Telescope targeting M dwarf stars in the near infrared. The survey has been operating since 2007 and gathering photometric time series of about 15000M dwarf stars brighter than J = 17mag. We identified and followed-up planet candidates from the most complete field and found two hot Jupiters around non-M dwarf hosts (WTS-1b & WTS-2b) but found no planets around the M dwarfs. 1. OBSERVATIONS AND DATA QUALITY The WFCAM Transit Survey (WTS) is an ongoing campaign survey on the United Kingdom Infrared Telescope (UKIRT) operating in the near-infrared, in J-band since 2007. The WTS targets 4 superfields, two in the winter and two in the summer sky and each covering 1.6 square degrees in size. A field is made up of 8 pawprints with slightly overlapping regions at the edges. The WTS operates as a backup program of the large UKIRT surveys (such as the UKIDSS surveys and UHS). The WTS runs in queue mode using 1 or 0.5 hour long observing blocks resulting 4 or 2 epochs per block respectively. The exposure time is 90s. Table 1 summarizes main properties of the fields and the current number of epochs. We expect to complete 1 2-3 fields before the closure of UKIRT in 2013. The main science goal to identify planets orbiting around M dwarfs or earlier stars. Other science topics includes but not limited to characterizing of low mass eclipsing binaries; studying M dwarf variability in the J band; identifying high proper motion sources e.g. cool white dwarfs, brown dwarfs, TNOs. Figure 1 summarizes our light curve noise properties in the 19hr field. The saturation appears at J = 13mag. We reach 3-4mmag precision at the bright end and at J = 17mag where the sky noise dominates we have a precision of 2%. Mainly because of follow-up possibilities, we limit our main science goal of finding planets to stellar sources brighter than J = 17magnitude, but for other science topics we used fainter sources, too. For reference we show the signal depth for a few star-planet systems. In Figure 2 we present a basic sensitivity diagram, showing the effect of the irregular observing pattern and number of epochs in the different fields. We used a simple simulation of a Neptune sized planet transiting around an M0 star. The simulation provides us an unbiased sample of different transit lengths and timings derived from random geometry of inclination (we allow grazing systems, too).

WTS-2 b: Too close for comfort?

We report the discovery of WTS-2 b, a typical hot Jupiter in an unusually close 1.02-day orbit to a K-dwarf star. This is the second planet to be discovered in the infrared light curves of the WFCAM Transit Survey (WTS) and is only one-and-a-half times the separation from its host star at which is would be destroyed by Roche lobe overflow. The predicted remaining lifetime of the planet is just 38 Myrs, assuming a tidal dissipation quality factor of Q'* = 106. The magnitude of Q'* is largely unconstrained by observations, thus WTS-2 b provides a useful calibration point for theories describing how frictional processes within a host star affect the tidal orbital evolution of its companion giant planets. It is expected that stars with large convective envelopes are more efficient at dissipating the orbital energy of the planet, and WTS-2 b provides an observational constraint in the sparsely populated K-dwarf regime. In addition, despite its relatively faint magnitude, the favourable size ratio of the WTS-2 star-planet system and the predicted hot equilibrium temperature of the planet will make it possible to characterise the planet's atmosphere via secondary eclipse measurements using existing ground-based instrumentation.

Effect of hot jupiters on the variability of stellar radiation

The effect of hot jupiters on the variability of light from 14 bright stars with planets that are listed in the Pulkovo Spectrophotometric Catalog (PSC) and the Hipparcos photometric catalog is studied. The spectral range 320-1080 nm and the Vt and Bt photometric magnitudes are examined. The photometric observations reveal microvariability of the star τBoo in the Vt and Bt bands with amplitudes of ~0m.018 and ~0m.020 and periods of ~3.4 d and ~3.2 d, respectively. The coincidence of the periods of the variability with the orbital period of τBoo b around the star indicates that the microvariability is associated with the orbit of the planet. Microvariability of the star 51 Peg is also observed in the Bt band with an amplitude of 0m.023, essentially the same as the magnitude of the decrease in the brightness as the planet crosses the star’s disk, but with a period of 11.6 d, a factor of 3 longer than the planet’s orbital period of 4.2 d. Apparently the planet 51 Peg b crosses the star’s disk once in three orbital cycles. This could indicate the presence of an additional body (another planet?) in the star-planet system that perturbs the planet’s orbit. According to spectrophotometric observations, no effect caused by even the very hottest jupiters was found.

A stringent upper limit to 18 cm radio emission from the extrasolar planet systemτBoötis

Context: It has been speculated for many years that some extrasolar planets may emit strong cyclotron emission at low radio frequencies in the range 10-100 MHz. Despite several attempts no such emission has yet been seen. Aims: The hot Jupiter system tau Bootis is one of the nearest (d=15 pc) exoplanets known to date. The gravitational influence of this massive hot Jupiter (M=6 M_jup) has locked the star-planet system, making the star rotate in P~3.3 days, similar to the orbital period of the planet. From the well established correlation between stellar rotation and radio luminosity, it is conceivable that the tau Bootis system emits strong radio emission at significantly higher frequencies than currently probed, which we aimed to investigate with this work. Methods: We observed tau Bootis with the Westerbork Synthesis Radio Telescope (WSRT) at a frequency of 1.7 GHz. for 12 hours in spectral line mode, reaching a noise level of 42 microJy/beam at the position of the target. Results: No 18cm radio emission is detected from tau Bootis, resulting in a 3 sigma upper limit of 0.13 mJy, corresponding to a 18cm radio luminosity of <3.7e13 erg/s/Hz. We observe tau Bootis to be two orders of magnitude fainter than expected from the stellar relation between radio luminosity and rotation velocity. Conclusions: This implies that either the tau Bootis system is underluminous in the radio compared to similar fast-rotating stars, or that we happened to observe the target during a low state of radio emission.

GLANCING VIEWS OF THE EARTH: FROM A LUNAR ECLIPSE TO AN EXOPLANETARY TRANSIT

It has been posited that lunar eclipse observations may help predict the in-transit signature of Earth-like extrasolar planets. However, a comparative analysis of the two phenomena addressing in detail the transport of stellar light through the planet's atmosphere has not yet been presented. Here, we proceed with the investigation of both phenomena by making use of a common formulation. Our starting point is a set of previously unpublished near-infrared spectra collected at various phases during the August 2008 lunar eclipse. We then take the formulation to the limit of an infinitely distant observer in order to investigate the in-transit signature of the Earth-Sun system as being observed from outside our Solar System. The refraction-bending of sunlight rays that pass through the Earth's atmosphere is a critical factor in the illumination of the eclipsed Moon. Likewise, refraction will have an impact on the in-transit transmission spectrum for specific planet-star systems depending on the refractive properties of the planet's atmosphere, the stellar size and the planet's orbital distance. For the Earth-Sun system, at mid-transit, refraction prevents the remote observer's access to the lower ~12-14 km of the atmosphere and, thus, also to the bulk of the spectroscopically-active atmospheric gases. We demonstrate that the effective optical radius of the Earth in transit is modulated by refraction and varies by ~12 km from mid-transit to 2nd contact. The refractive nature of atmospheres, a property which is rarely accounted for in published investigations, will pose additional challenges to the characterization of Earth-like extrasolar planets. Refraction may have a lesser impact for Earth-like extrasolar planets within the habitable zone of some M-type stars.