PyTranSpot\texttt{PyTranSpot} - A tool for multiband light curve modeling of planetary transits and stellar spots

Several studies have shown that stellar activity features, such as occulted and non-occulted starspots, can affect the measurement of transit parameters biasing studies of transit timing variations and transmission spectra. We present $\texttt{PyTranSpot}$, which we designed to model multiband transit light curves showing starspot anomalies, inferring both transit and spot parameters. The code follows a pixellation approach to model the star with its corresponding limb darkening, spots, and transiting planet on a two dimensional Cartesian coordinate grid. We combine $\texttt{PyTranSpot}$ with an MCMC framework to study and derive exoplanet transmission spectra, which provides statistically robust values for the physical properties and uncertainties of a transiting star-planet system. We validate $\texttt{PyTranSpot}$'s performance by analyzing eleven synthetic light curves of four different star-planet systems and 20 transit light curves of the well-studied WASP-41b system. We also investigate the impact of starspots on transit parameters and derive wavelength dependent transit depth values for WASP-41b covering a range of 6200-9200 $\AA$, indicating a flat transmission spectrum.Comment: 17 pages, 22 figures; accepted for publication in Astronomy & Astrophysic

Vehicular traffic flow at an intersection with the possibility of turning

We have developed a Nagel-Schreckenberg cellular automata model for describing of vehicular traffic flow at a single intersection. A set of traffic lights operating in fixed-time scheme controls the traffic flow. Open boundary condition is applied to the streets each of which conduct a uni-directional flow. Streets are single-lane and cars can turn upon reaching to the intersection with prescribed probabilities. Extensive Monte Carlo simulations are carried out to find the model flow characteristics. In particular, we investigate the flows dependence on the signalisation parameters, turning probabilities and input rates. It is shown that for each set of parameters, there exist a plateau region inside which the total outflow from the intersection remains almost constant. We also compute total waiting time of vehicles per cycle behind red lights for various control parameters.Comment: 8 pages, 17 eps figures, Late

Opening a new window to other worlds with spectropolarimetry

A high level of diversity has already been observed among the planets of our own Solar System. As such, one expects extrasolar planets to present a wide range of distinctive features, therefore the characterisation of Earth- and super Earth-like planets is becoming of key importance in scientific research. The SEARCH (Spectropolarimetric Exoplanet AtmospheRe CHaracerisation) mission proposal of this paper represents one possible approach to realising these objectives. The mission goals of SEARCH include the detailed characterisation of a wide variety of exoplanets, ranging from terrestrial planets to gas giants. More specifically, SEARCH will determine atmospheric properties such as cloud coverage, surface pressure and atmospheric composition, and may also be capable of identifying basic surface features. To resolve a planet with a semi major axis of down to 1.4AU and 30pc distant SEARCH will have a mirror system consisting of two segments, with elliptical rim, cut out of a parabolic mirror. This will yield an effective diameter of 9 meters along one axis. A phase mask coronagraph along with an integral spectrograph will be used to overcome the contrast ratio of star to planet light. Such a mission would provide invaluable data on the diversity present in extrasolar planetary systems and much more could be learned from the similarities and differences compared to our own Solar System. This would allow our theories of planetary formation, atmospheric accretion and evolution to be tested, and our understanding of regions such as the outer limit of the Habitable Zone to be further improved.Comment: 23 pages, accepted for publication in Experimental Astronom

A scenario of planet erosion by coronal radiation

Context: According to theory, high-energy emission from the coronae of cool stars can severely erode the atmospheres of orbiting planets. No observational tests of the long term effects of erosion have yet been made. Aims: To analyze the current distribution of planetary mass with X-ray irradiation of the atmospheres in order to make an observational assessment of the effects of erosion by coronal radiation. Methods: We study a large sample of planet-hosting stars with XMM-Newton, Chandra and ROSAT; make a careful identification of X-ray counterparts; and fit their spectra to make accurately measurements of the stellar X-ray flux. Results: The distribution of the planetary masses with X-ray flux suggests that erosion has taken place: most surviving massive planets, (M_p sin i >1.5 M_J), have been exposed to lower accumulated irradiation. Heavy erosion during the initial stages of stellar evolution is followed by a phase of much weaker erosion. A line dividing these two phases could be present, showing a strong dependence on planet mass. Although a larger sample will be required to establish a well-defined erosion line, the distribution found is very suggestive. Conclusions: The distribution of planetary mass with X-ray flux is consistent with a scenario in which planet atmospheres have suffered the effects of erosion by coronal X-ray and EUV emission. The erosion line is an observational constraint to models of atmospheric erosion.Comment: A&A 511, L8 (2010). 4 pages, 3 figures, 1 online table (included). Language edited; corrected a wrong unit conversion (g/s -> M_J/Gyr); corrected values in column 12 of Table 1 (slightly underestimated in first version), and Figure 2 updated accordingl

Transit timing analysis of CoRoT-1b

CoRoT, the pioneer space-based transit search, steadily provides thousands of high-precision light curves with continuous time sampling over periods of up to 5 months. The transits of a planet perturbed by an additional object are not strictly periodic. By studying the transit timing variations (TTVs), additional objects can be detected in the system. A transit timing analysis of CoRoT-1b is carried out to constrain the existence of additional planets in the system. We used data obtained by an improved version of the CoRoT data pipeline (version 2.0). Individual transits were fitted to determine the mid-transit times, and we analyzed the derived $O-C$ diagram. N-body integrations were used to place limits on secondary planets. No periodic timing variations with a period shorter than the observational window (55 days) are found. The presence of an Earth-mass Trojan is not likely. A planet of mass greater than $\sim 1$ Earth mass can be ruled out by the present data if the object is in a 2:1 (exterior) mean motion resonance with CoRoT-1b. Considering initially circular orbits: (i) super-Earths (less than 10 Earth-masses) are excluded for periods less than about 3.5 days, (ii) Saturn-like planets can be ruled out for periods less than about 5 days, (iii) Jupiter-like planets should have a minimum orbital period of about 6.5 days.Comment: 6 pages, accepted at A&

Estimation of the XUV radiation onto close planets and their evaporation

Context: The current distribution of planet mass vs. incident stellar X-ray flux supports the idea that photoevaporation of the atmosphere may take place in close-in planets. Integrated effects have to be accounted for. A proper calculation of the mass loss rate due to photoevaporation requires to estimate the total irradiation from the whole XUV range. Aims: The purpose of this paper is to extend the analysis of the photoevaporation in planetary atmospheres from the accessible X-rays to the mostly unobserved EUV range by using the coronal models of stars to calculate the EUV contribution to the stellar spectra. The mass evolution of planets can be traced assuming that thermal losses dominate the mass loss of their atmospheres. Methods: We determine coronal models for 82 stars with exoplanets that have X-ray observations available. Then a synthetic spectrum is produced for the whole XUV range (~1-912 {\AA}). The determination of the EUV stellar flux, calibrated with real EUV data, allows us to calculate the accumulated effects of the XUV irradiation on the planet atmosphere with time, as well as the mass evolution for planets with known density. Results: We calibrate for the first time a relation of the EUV luminosity with stellar age valid for late-type stars. In a sample of 109 exoplanets, few planets with masses larger than ~1.5 Mj receive high XUV flux, suggesting that intense photoevaporation takes place in a short period of time, as previously found in X-rays. The scenario is also consistent with the observed distribution of planet masses with density. The accumulated effects of photoevaporation over time indicate that HD 209458b may have lost 0.2 Mj since an age of 20 Myr. Conclusions: Coronal radiation produces rapid photoevaporation of the atmospheres of planets close to young late-type stars. More complex models are needed to explain fully the observations.Comment: Accepted by A&A. 10 pages, 8 figures, 7 Tables (2 online). Additional online material includes 7 pages, 6 figures and 6 tables, all include

The Kepler-11 system: evolution of the stellar high-energy emission and {initial planetary} atmospheric mass fractions

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.Comment: 8 pages, 3 figure

Suppressed Far-UV stellar activity and low planetary mass-loss in the WASP-18 system

WASP-18 hosts a massive, very close-in Jupiter-like planet. Despite its young age (R′HK activity parameter lies slightly below the basal level; there is no significant time-variability in the log R′HK value; there is no detection of the star in the X-rays. We present results of far-UV observations of WASP-18 obtained with COS on board of HST aimed at explaining this anomaly. From the star’s spectral energy distribution, we infer the extinction (E(B − V) ≈ 0.01mag) and then the ISM column density for a number of ions, concluding that ISM absorption is not the origin of the anomaly. We measure the flux of the four stellar emission features detected in the COS spectrum (C II, C III, C IV, Si IV). Comparing the C II/C IV flux ratio measured for WASP-18 with that derived from spectra of nearby stars with known age, we see that the far-UV spectrum of WASP-18 resembles that of old (>5Gyr), inactive stars, in stark contrast with its young age. We conclude that WASP-18 has an intrinsically low activity level, possibly caused by star-planet tidal interaction, as suggested by previous studies. Re-scaling the solar irradiance reference spectrum to match the flux of the Si IV line, yields an XUV integrated flux at the planet orbit of 10.2 erg s−1 cm−2. We employ the rescaled XUV solar fluxes to model of the planetary upper atmosphere, deriving an extremely low thermal mass-loss rate of 10−20MJ Gyr−1. For such high-mass planets, thermal escape is not energy limited, but driven by Jeans escape

Temporal variations in the evaporating atmosphere of the exoplanet HD 189733b

Atmospheric escape has been detected from the exoplanet HD 209458b through transit observations of the hydrogen Lyman-alpha line. Here we present spectrally resolved Lyman-alpha transit observations of the exoplanet HD 189733b at two different epochs. These HST/STIS observations show for the first time, that there are significant temporal variations in the physical conditions of an evaporating planetary atmosphere. While atmospheric hydrogen is not detected in the first epoch observations, it is observed at the second epoch, producing a transit absorption depth of 14.4+/-3.6% between velocities of -230 to -140 km/s. Contrary to HD 209458b, these high velocities cannot arise from radiation pressure alone and require an additional acceleration mechanism, such as interactions with stellar wind protons. The observed absorption can be explained by an atmospheric escape rate of neutral hydrogen atoms of about 10^9 g/s, a stellar wind with a velocity of 190 km/s and a temperature of ~10^5K. An X-ray flare from the active star seen with Swift/XRT 8 hours before the second-epoch observation supports the idea that the observed changes within the upper atmosphere of the planet can be caused by variations in the stellar wind properties, or by variations in the stellar energy input to the planetary escaping gas (or a mix of the two effects). These observations provide the first indication of interaction between the exoplanet's atmosphere and stellar variations.Comment: To be published in A&A Letters, June 28, 201

CoRoT-22 b: a validated 4.9 RE exoplanet in 10-day orbit

The CoRoT satellite has provided high-precision photometric light curves for more than 163,000 stars and found several hundreds of transiting systems compatible with a planetary scenario. If ground-based velocimetric observations are the best way to identify the actual planets among many possible configurations of eclipsing binary systems, recent transit surveys have shown that it is not always within reach of the radial-velocity detection limits. In this paper, we present a transiting exoplanet candidate discovered by CoRoT whose nature cannot be established from ground-based observations, and where extensive analyses are used to validate the planet scenario. They are based on observing constraints from radial-velocity spectroscopy, adaptive optics imaging and the CoRoT transit shape, as well as from priors on stellar populations, planet and multiple stellar systems frequency. We use the fully Bayesian approach developed in the PASTIS analysis software, and conclude that the planet scenario is at least 1400 times more probable than any other false positive scenario. The primary star is a metallic solar-like dwarf, with Ms = 1.099+-0.049 Msun and Rs = 1.136 (+0.038,-0.090) Rsun . The validated planet has a radius of Rp = 4.88 (+0.17,-0.39) RE and mass less than 49 ME. Its mean density is smaller than 2.56 g/cm^3 and orbital period is 9.7566+-0.0012 days. This object, called CoRoT-22 b, adds to a large number of validated Kepler planets. These planets do not have a proper measurement of the mass but allow statistical characterization of the exoplanet population