We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative transfer and fluid dynamics. We use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate the atmospheres of the well-studied hot Jupiters HD~209458b and HD~189733b which both have dayside--nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH$_4$), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. 3D mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D - and even 2D - chemical mixing. We find that CH$_4$, carbon dioxide (CO$_2$) and ammonia (NH$_3$) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, detectable with near-future instruments.

We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of $0.001-0.075\,\mathrm{M_{\odot}}$ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H-He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by $\sim1-2\%$ at the stellar-substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the $3.5-5.5\,\mathrm{\mu m}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.

Low-mass stars are known to have magnetic fields that are believed to be of dynamo origin. Two complementary techniques are principally used to characterise them. Zeeman-Doppler imaging (ZDI) can determine the geometry of the large-scale magnetic field while Zeeman broadening can assess the total unsigned flux including that associated with small-scale structures such as spots. In this work, we study a sample of stars that have been previously mapped with ZDI. We show that the average unsigned magnetic flux follows an activity-rotation relation separating into saturated and unsaturated regimes. We also compare the average photospheric magnetic flux recovered by ZDI, $\langle B_V\rangle$, with that recovered by Zeeman broadening studies, $\langle B_I\rangle$. In line with previous studies, $\langle B_V\rangle$ ranges from a few % to $\sim$20% of $\langle B_I\rangle$. We show that a power law relationship between $\langle B_V\rangle$ and $\langle B_I\rangle$ exists and that ZDI recovers a larger fraction of the magnetic flux in more active stars. Using this relation, we improve on previous attempts to estimate filling factors, i.e. the fraction of the stellar surface covered with magnetic field, for stars mapped only with ZDI. Our estimated filling factors follow the well-known activity-rotation relation which is in agreement with filling factors obtained directly from Zeeman broadening studies. We discuss the possible implications of these results for flux tube expansion above the stellar surface and stellar wind models.

In the original manuscript arXiv:1710.10269, we presented a grid of forward model transmission spectra for hot Jupiter exoplanets. However, we recently identified an error in the treatment of rainout in our 1D atmosphere model ATMO. We explain the error, the correction, its validation and the changed conclusions in this erratum. The correction of this error led to changes in the equilibrium chemical abundances using rainout condensation and thereby the transmission spectra. We note that this error only affects the online library that includes rainout condensation, the library with local condensation (without rainout) is unaffected. The online library was updated with the correction in July 2018.

We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.

Future NASA concept missions that are currently under study, like Habitable Exoplanet Imaging Mission (HabEx) & Large Ultra-Violet Optical Infra Red (LUVOIR) Surveyor, would discover a large diversity of exoplanets. We propose here a classification scheme that distinguishes exoplanets into different categories based on their size and incident stellar flux, for the purpose of providing the expected number of exoplanets observed (yield) with direct imaging missions. The boundaries of this classification can be computed using the known chemical behavior of gases and condensates at different pressures and temperatures in a planetary atmosphere. In this study, we initially focus on condensation curves for sphalerite ZnS, H2O, CO2 and CH4. The order in which these species condense in a planetary atmosphere define the boundaries between different classes of planets. Broadly, the planets are divided into rocky (0.5 - 1.0RE), super-Earths (1.0- 1.75RE), sub-Neptunes (1.75-3.5RE), sub-Jovians (3.5 - 6.0RE) and Jovians (6-14.3RE) based on their planet sizes, and 'hot', 'warm' and 'cold' based on the incident stellar flux. We then calculate planet occurrence rates within these boundaries for different kinds of exoplanets, \eta_{planet}, using the community co-ordinated results of NASA's Exoplanet Program Analysis Group's Science Analysis Group-13 (SAG-13). These occurrence rate estimates are in turn used to estimate the expected exoplanet yields for direct imaging missions of different telescope diameter.

Based on optical high-resolution spectra obtained with CFHT/ESPaDOnS, we present new measurements of activity and magnetic field proxies of 442 low-mass K5-M7 dwarfs. The objects were analysed as potential targets to search for planetary-mass companions with the new spectropolarimeter and high-precision velocimeter, SPIRou. We have analysed their high-resolution spectra in an homogeneous way: circular polarisation, chromospheric features, and Zeeman broadening of the FeH infrared line. The complex relationship between these activity indicators is analysed: while no strong connection is found between the large-scale and small-scale magnetic fields, the latter relates with the non-thermal flux originating in the chromosphere. We then examine the relationship between various activity diagnostics and the optical radial-velocity jitter available in the literature, especially for planet host stars. We use this to derive for all stars an activity merit function (higher for quieter stars) with the goal of identifying the most favorable stars where the radial-velocity jitter is low enough for planet searches. We find that the main contributors to the RV jitter are the large-scale magnetic field and the chromospheric non-thermal emission. In addition, three stars (GJ 1289, GJ 793, and GJ 251) have been followed along their rotation using the spectropolarimetric mode, and we derive their magnetic topology. These very slow rotators are good representatives of future SPIRou targets. They are compared to other stars where the magnetic topology is also known. The poloidal component of the magnetic field is predominent in all three stars.

Deep optical imaging is used to study time-domain properties of young brown dwarfs in sigma Orionis over typical rotational time-scales and to search for new substellar and planetary-mass cluster members. We used VIMOS at the VLT to monitor a 24 arcmin x 16 arcmin field in the I-band. Using the individual images from this run we investigated the photometric time series of nine substellar cluster members with masses from 10 to 60 Mj. The deep stacked image shows cluster members down to $\approx$ 5 Mj. We search for new planetary mass objects by combining our deep I-band photometry with public J-band magnitudes and by examining the nearby environment of known very low mass members for possible companions. We find two brown dwarfs, with significantly variable, aperiodic light curves, both with masses around 50 Mj, one of which was previously unknown to be variable. The physical mechanism responsible for the observed variability is likely different for the two objects. The variability of the first object, a single-lined spectroscopic binary, is most likely linked to its accretion disc, the second may be caused by variable extinction by large grains. We find five new candidate members from the colour-magnitude diagram and three from a search for companions within 2000 au. We rule all eight sources out as potential members based on non-stellar shape and/or infrared colours. The I-band photometry is made available as a public dataset. We present two variable brown dwarfs. One is consistent with ongoing accretion, the other, exhibiting apparent transient variability without presence of an accretion disc. Our analysis confirms the existing census of substellar cluster members down to $\approx$ 7 Mj. The zero result from our companion search agrees with the low occurrence rate of wide companions to brown dwarfs found in other works.

Magnetic activity and rotation are known to be intimately linked for low-mass stars. Understanding rotation evolution over the stellar lifetime is therefore an important goal within stellar astrophysics. In recent years, there has been increased focus on how the complexity of the stellar magnetic field affects the rate of angular momentum-loss from a star. This is a topic that Zeeman-Doppler imaging (ZDI), a technique that is capable of reconstructing the large-scale magnetic field topology of a star, can uniquely address. Using a potential field source surface model, we estimate the open flux, mass loss-rate and angular momentum-loss rates for a sample of 66 stars that have been mapped with ZDI. We show that the open flux of a star is predominantly determined by the dipolar component of its magnetic field for our choice of source surface radius. We also show that, on the main sequence, the open flux, mass- and angular momentum-loss rates increase with decreasing Rossby number. The exception to this rule is stars less massive than $0.3M_{\odot}$. Previous work suggests that low mass M dwarfs may possess either strong, ordered and dipolar fields or weak and complex fields. This range of field strengths results in a large spread of angular momentum-loss rates for these stars and has important consequences for their spin down behaviour. Additionally, our models do not predict a transition in the mass-loss rates at the so called wind dividing line noted from Ly$\alpha$ studies.

Hot Jupiters are giant Jupiter-like exoplanets that orbit 100x closer to their host stars than Jupiter does to the Sun. These planets presumably form in the outer part of the primordial disc from which both the central star and surrounding planets are born, then migrate inwards and yet avoid falling into their host star. It is however unclear whether this occurs early in the lives of hot Jupiters, when still embedded within protoplanetary discs, or later, once multiple planets are formed and interact. Although numerous hot Jupiters were detected around mature Sun-like stars, their existence has not yet been firmly demonstrated for young stars, whose magnetic activity is so intense that it overshadows the radial velocity signal that close-in giant planets can induce. Here we show that hot Jupiters around young stars can be revealed from extended sets of high-resolution spectra. Once filtered-out from the activity, radial velocities of V830 Tau derived from new data collected in late 2015 exhibit a sine wave of period 4.93 d and semi-amplitude 75 m/ s, detected with a false alarm probability <0.03%. We find that this signal is fully unrelated to the 2.741-d rotation period of V830 Tau and we attribute it to the presence of a 0.77 Jupiter mass planet orbiting at a distance of 0.057 au from the host star. Our result demonstrates that hot Jupiters can migrate inwards in <2 Myr, most likely as a result of planet-disc interactions, and thus yields strong support to the theory of giant planet migration in gaseous protoplanetary discs.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

SPIRou is a near-infrared spectropolarimeter and high-precision radial-velocity instrument, to be mounted on the 3.6m Canada-France-Hawaii telescope ontop Maunakea and to be offered to the CFHT community from 2018. It focuses on two main scientific objectives : (i) the search and study of Earth-like planets around M dwarfs, especially in their habitable zone and (ii) the study of stellar and planetary formation in the presence of stellar magnetic field. The SPIRou characteristics (complete coverage of the near infrared wavelengths, high resolution, high stability and efficiency, polarimetry) also allow many other programs, e.g., magnetic fields and atmospheres of M dwarfs and brown dwarfs, star-planet interactions, formation and characterization of massive stars, dynamics and atmospheric chemistry of planets in the solar system.

Wood et al suggested that mass-loss rate is a function of X-ray flux ($\dot{M} \propto F_x^{1.34}$) for dwarf stars with $F_x \lesssim F_{x,6} \equiv 10^6$ erg cm$^{-2}$ s$^{-1}$. However, more active stars do not obey this relation. These authors suggested that the break at $F_{x,6}$ could be caused by significant changes in magnetic field topology that would inhibit stellar wind generation. Here, we investigate this hypothesis by analysing the stars in Wood et al's sample that had their surface magnetic fields reconstructed through Zeeman-Doppler Imaging (ZDI). Although the solar-like outliers in the $\dot{M}$ -- $F_x$ relation have higher fractional toroidal magnetic energy, we do not find evidence of a sharp transition in magnetic topology at $F_{x,6}$. To confirm this, further wind measurements and ZDI observations at both sides of the break are required. As active stars can jump between states with highly toroidal to highly poloidal fields, we expect significant scatter in magnetic field topology to exist for stars with $F_x \gtrsim F_{x,6}$. This strengthens the importance of multi-epoch ZDI observations. Finally, we show that there is a correlation between $F_x$ and magnetic energy, which implies that $\dot{M}$ -- magnetic energy relation has the same qualitative behaviour as the original $\dot{M}$ -- $F_x$ relation. No break is seen in any of the $F_x$ -- magnetic energy

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star LkCa4 within the MaTYSSE programme, involving ESPaDOnS at the Canada-France-Hawaii Telescope. Despite an age of only 2Myr and a similarity with prototypical classical T Tauri stars, LkCa4 shows no evidence for accretion and probes an interesting transition stage for star and planet formation. Large profile distortions and Zeeman signatures are detected in the unpolarized and circularly-polarized lines of LkCa4 using Least-Squares Deconvolution (LSD), indicating the presence of brightness inhomogeneities and magnetic fields at the surface of LkCa4. Using tomographic imaging, we reconstruct brightness and magnetic maps of LkCa4 from sets of unpolarized and circularly-polarized LSD profiles. The large-scale field is strong and mainly axisymmetric, featuring a ~2kG poloidal component and a ~1kG toroidal component encircling the star at equatorial latitudes - the latter making LkCa4 markedly different from classical TTauri stars of similar mass and age. The brightness map includes a dark spot overlapping the magnetic pole and a bright region at mid latitudes - providing a good match to the contemporaneous photometry. We also find that differential rotation at the surface of LkCa4 is small, typically ~5.5x weaker than that of the Sun, and compatible with solid-body rotation. Using our tomographic modelling, we are able to filter out the activity jitter in the RV curve of LkCa4 (of full amplitude 4.3km/s) down to a rms precision of 0.055km/s. Looking for hot Jupiters around young Sun-like stars thus appears feasible, even though we find no evidence for such planets around LkCa4.

We present Rossiter-McLaughlin observations of WASP-13b and WASP-32b and determine the sky-projected angle between the normal of the planetary orbit and the stellar rotation axis ($\lambda$). WASP-13b and WASP-32b both have prograde orbits and are consistent with alignment with measured sky-projected angles of $\lambda={8^{\circ}}^{+13}_{-12}$ and $\lambda={-2^{\circ}}^{+17}_{-19}$, respectively. Both WASP-13 and WASP-32 have $T_{\mathrm{eff}}<6250$K and therefore these systems support the general trend that aligned planetary systems are preferentially found orbiting cool host stars. A Lomb-Scargle periodogram analysis was carried out on archival SuperWASP data for both systems. A statistically significant stellar rotation period detection (above 99.9\% confidence) was identified for the WASP-32 system with $P_{\mathrm{rot}}=11.6 \pm 1.0 $ days. This rotation period is in agreement with the predicted stellar rotation period calculated from the stellar radius, $R_{\star}$, and $v \sin i$ if a stellar inclination of $i_{\star}=90^{\circ}$ is assumed. With the determined rotation period, the true 3D angle between the stellar rotation axis and the planetary orbit, $\psi$, was found to be $\psi=11^{\circ} \pm 14$. We conclude with a discussion on the alignment of systems around cool host stars with $T_{\mathrm{eff}}<6150$K by calculating the tidal dissipation timescale. We find that systems with short tidal dissipation timescales are preferentially aligned and systems with long tidal dissipation timescales have a broad range of obliquities.

We present the discovery of four new transiting hot jupiters, detected mainly from SuperWASP-North and SOPHIE observations. These new planets, WASP-52b, WASP-58b, WASP-59b, and WASP-60b, have orbital periods ranging from 1.7 to 7.9 days, masses between 0.46 and 0.94 M_Jup, and radii between 0.73 and 1.49 R_Jup. Their G1 to K5 dwarf host stars have V magnitudes in the range 11.7-13.0. The depths of the transits are between 0.6 and 2.7%, depending on the target. With their large radii, WASP-52b and 58b are new cases of low-density, inflated planets, whereas WASP-59b is likely to have a large, dense core. WASP-60 shows shallow transits. In the case of WASP-52 we also detected the Rossiter-McLaughlin anomaly via time-resolved spectroscopy of a transit. We measured the sky-projected obliquity lambda = 24 (+17/-9) degrees, indicating that WASP-52b orbits in the same direction as its host star is rotating and that this prograde orbit is slightly misaligned with the stellar equator. These four new planetary systems increase our statistics on hot jupiters, and provide new targets for follow-up studies.

We investigate the effects of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b. We use a simplified dynamical model and a robust chemical network, as opposed to previous studies which have used a three dimensional circulation model coupled to a simple chemical kinetics scheme. The temperature structure and distribution of the main atmospheric constituents are calculated in the limit of an atmosphere that rotates as a solid body with an equatorial rotation rate of 1 km/s. Such motion mimics a uniform zonal wind which resembles the equatorial superrotation structure found by three dimensional circulation models. The uneven heating of this tidally locked planet causes, even in the presence of such a strong zonal wind, large temperature contrasts between the dayside and nightside, of up to 800 K. This would result in important longitudinal variations of some molecular abundances if the atmosphere were at chemical equilibrium. The zonal wind, however, acts as a powerful disequilibrium process. We identify the existence of a pressure level of transition between two regimes, which may be located between 100 and 0.1 mbar depending on the molecule. Below this transition layer, chemical equilibrium holds, while above it, the zonal wind tends to homogenize the chemical composition of the atmosphere, bringing molecular abundances in the limb and nightside regions close to chemical equilibrium values characteristic of the dayside, i.e. producing an horizontal quenching effect in the abundances. Reasoning based on timescales arguments indicates that horizontal and vertical mixing are likely to compete in HD 209458b's atmosphere, producing a complex distribution where molecular abundances are quenched horizontally to dayside values and vertically to chemical equilibrium values characteristic of deep layers.

The early evolution of Earth's atmosphere and the origin of life took place at a time when physical conditions at the Earth where radically different from its present state. The radiative input from the Sun was much enhanced in the high-energy spectral domain, and in order to model early planetary atmospheres in detail, a knowledge of the solar radiative input is needed. We present an investigation of the atmospheric parameters, state of evolution and high-energy fluxes of the nearby star kap^1 Cet, previously thought to have properties resembling those of the early Sun. Atmospheric parameters were derived from the excitation/ionization equilibrium of Fe I and Fe II, profile fitting of Halpha and the spectral energy distribution. The UV irradiance was derived from FUSE and HST data, and the absolute chromospheric flux from the Halpha line core. From careful spectral analysis and the comparison of different methods we propose for kap^1 Cet the following atmospheric parameters: Teff = 5665+/-30 K (Halpha profile and energy distribution), log g = 4.49+/-0.05 dex (evolutionary and spectroscopic) and [Fe/H] = +0.10+/-0.05 dex (Fe II lines). The UV radiative properties of kap^1 Cet indicate that its flux is some 35% lower than the current Sun's between 210 and 300 nm, it matches the Sun's at 170 nm and increases to at least 2-7 times higher than the Sun's between 110 and 140 nm. The use of several indicators ascribes an age to kap^1 Cet in the interval ~0.4-0.8 Gyr and the analysis of the theoretical HR diagram suggests a mass ~1.04 Msun. This star is thus a very close analog of the Sun when life arose on Earth and Mars is thought to have lost its surface bodies of liquid water. Photochemical models indicate that the enhanced UV emission leads to a significant increase in photodissociation rates compared with those commonly assumed of the early Earth. Our results show that reliable calculations of the chemical composition of early planetary atmospheres need to account for the stronger solar photodissociating UV irradiation.