Svensk Sammanfattning Vetenskap om exoplaneter är väldig kompetitiv. Teleskopobservationer är svåra och dyra. Korrekt selektion av objekt och observationsschema är viktiga för att detektera biomarkörer utan att missa dem. Att planera och schemalägga är där- för lika viktiga som observationerna. För att kunna detektera biomarkörer i en exoplanets atmosfär observeras planeten när den passerar framför stjärnan den cirkulerar runt om. Teleskopen och dess placering ger riktlinjer till observation- sprocessen, så väl som planetens och stjärnans egenskaper. I den första delen har vi diskuterat litteratur om exoplanetvetenskap, speciellt om hur man kan bestämma atmosfärens sammansättning och modellering av planeternas atmosfärer. I den an- dra delen presenterar vi en sådan programvara, utvecklad för CRIRES+. Vi visar några användningsexempel och presenterar förslag av planeter som kan observeras över året. Vi jämför resultaten med ett annat liknande verktyg och diskuterar skillnaderna.
1 Abstract
In recent decades, thousands of exoplanets have been discovered. The next step is to characterize the observed planets in terms of their radii, masses, density, physical conditions and composition of their atmospheres. Several space-based observatories such as TESS and CHEOPS have started determining the first three observables but characterization of exoplanetary atmospheres is waiting for observation campaigns with instruments like CRIRES+ at the VLT and NIRSpec on the JWST. To ensure the ef- ficiency of data acquisition, careful planning of observations is necessary. In this project we developed a software tool to select and rank candidates based on the feasibility of observations of atmospheric features during transits with CRIRES+. We also review different techniques to retrieve transmission spectra from transit observations and modeling of exoplanet atmospheres in order to clarify the re- quirements for the data. Our CRIRES-planning-tool is built on astronomical observation planning methods from astropy and astroplan and the exposure time calculator designed for CRIRES+ by ESO and UU. We conclude that observations for atmospheric characterization with CRIRES+ are fea- sible. However, we observed that for a robust candidate selection, careful iterative tuning of proposed constraints is required.
2 Contents
1 Introduction 4
2 Theory and literature review5 2.1 Exoplanet transit...... 5 2.1.1 Photometric lightcurve and transmission spectroscopy...... 10 2.2 Atmospheric retrieval methods for transmission spectroscopy...... 13 2.2.1 A novel planetary model-independent transmission spectrum recovery method.. 13 2.2.2 Cross correlation methods to retrieve the transmission spectrum of an exoplane- tary atmosphere...... 16 2.2.3 Optimal estimation and bayesian inference...... 17 2.3 Planetary atmospheric models...... 20 2.3.1 Parametric forward models...... 23 2.3.2 Scattering, albedos and the two-stream approximation...... 28 2.3.3 Self-consistent one-dimensional planetary models...... 36 2.3.4 General circulation models and three dimensional radiative hydrodynamics mod- els, detailed cloud and dust physics...... 44
3 Summary and conclusion of the literature review 51
4 Description and results of the CRIRES-planning-tool 54 4.1 Functionalities and methods of CRIRES-planning-tool...... 55 4.2 Constraints for candidates...... 59 4.3 Signal-to-noise ratio and the exposure time calculator...... 61 4.4 Input parameters to the ETC, restrictions and future improvements...... 62 4.5 Test and comparison of results...... 63 4.6 Results of detectable transits for one year...... 64
5 Discussion 69
6 Conclusion 70
3 1 Introduction
Today, the atmospheric composition of the planets in our Solar System has been roughly determined (de Pater and Lissauer, 2010). Besides, discoveries of exoplanets have revealed that the architecture of our Solar System is not typical at all (Wright et al., 2011). Therefore, we may expect very different planetary compositions, and atmospheric properties than can be found in our Solar System. On the other side, investigating their compositions will bring us closer to answer one of the most interesting questions about our own existence. Are we alone in the universe? However, we shall see, before we are able to answer this question we still have a long way to go. The next step will be to characterize planetary atmospheres and learn about the habitability of these remote and different worlds. To shed light onto the question of atmospheric composition of exoplanets, the CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES)(Kaeufl et al., 2004) was developed. An upgraded version of CRIRES, now called CRIRES+, is mounted at Paranal observatory on the very large telescope, (VLT), unit telescope 3 (UT3). To plan observations with CRIRES+ we developed a software tool called (CRIRES-planning-tool). Although it was developed to plan transit observations, it may easily be extended to other kinds of observations with CRIRES+. With the upgrade CRIRES+ has now enhanced observation efficiency, increased wavelength coverage, an upgrade of the focal plane detector array to increase the number of accessible diffraction orders, new wavelength calibration methods to reach wavelength precision of 5 m/s and a new spectropolarimetric unit (Follert et al., 2014). Previous scientific highlights with CRIRES were measuring the length of an exoplanet day for the first time (Snellen et al., 2014), creating the first weather map for the nearest brown dwarf to Earth (Crossfield et al., 2014) or finding the first superstorm on an exoplanet, HD209458 b and measuring its mass (Snellen et al., 2010a) to name but a few.
To plan transit observations of exoplanets we used the python libraries astropy,(Price-Whelan et al., 2018) and astroplan,(Morris et al., 2018), which provide useful methods to assess times and coordi- nates for observability at observatories all over the world. To calculate the signal-to-noise ratio (SNR) for an observation we use the Exposure Time Calculator (ETC), developed by the European Southern Observatory (ESO) in Garching, Munich, Germany together with Uppsala University (UU)1.
CRIRES-planning-tool is written in python and selects candidates from the NASA Exoplanet Archive (Akeson et al., 2013) to assess which candidates can be observed when, what signal-to-noise ratio can be reached for CRIRES+, and ranks the planets after criteria deduced from these properties. The output of the software tool is used for decision making, which targets should be observed and when, and provides an important basis to plan for best efficiency in time and data quality. This manuscript provides an overview of the structure of the software tool and describes and discusses some of the used methods. The software comes with an extensive documentation. The README file provides information about the different functionalities, accessibility, structure, installation, and further development possibilities of the CRIRES-planning-tool. The tool is available on GitHub2.
The thesis contains two major parts. The first part (chapter2) gives an overview of the present state of exoplanetary astrophysical techniques, both in observation and theory. In section 2.1 we introduce the reader to the technique of transit photometry, highlighting its strengths and limitations, and derive the basic equations for transmission spectroscopy. In section 2.2 we present two important methods for transmission spectra retrieval, and two techniques to compare retrieved spectra with synthesized spectra
1https://etctestpub.eso.org/observing/etc/crires 2https://github.com/jonaszubindu/CRIRES-planning-tool
4 from planetary models. In section 2.3 we provide an extensive overview about the different physics treated by a planetary model, such as radiative transfer, thermodynamics, convection, chemistry, clouds, and hydrodynamics and present examples to the different kinds of planetary models developed to study exoplanetary atmospheres. At the end of this part we will give a short summary and conclusion of the reviewed topics.
In the second part (chapter4), we shall provide the reader with a brief introduction to the CRIRES-planning-tool. In sections 4.1, and 4.2, we explain the functionalities and used methods to select candidates for obser- vations and assess their observability from the Paranal observatory, in section 4.3 we give an overview about the ETC and the computation of signal-to-noise ratio, and in section 4.4 explain each parameter used for the ETC. In section 4.5 we compare our tool with two other transit observation planning tools, highlighting differences between these tools and verifying our own findings. Finally in 4.6 we present preliminary results of possible transit observations over the course of one year. A discussion of the CRIRES-planning-tool can be found in chapter5, and a short conclusion of the findings is given in6.
2 Theory and literature review
One idea to determine the composition of an exoplanet’s atmosphere is to retrieve the transmission spectrum of the star light passing through the atmosphere and investigate the spectrum for known molecular transitions. Since planetary atmospheres consist of molecules, and in case of hot Jupiters, atoms, one aims at observing only molecular bands. Subtracting the stellar spectrum, and if the observatory is based on Earth subtracting the telluric spectrum, the exoplanet’s transmission spectra should reveal its atmospheric composition and allow the deduction of the prevailing conditions at the surface of the planet. However, one is faced with several challenges. We will present these challenges in the following chapter one by one, and solutions or methods to circumnavigate these issues and retrieve as much information as possible about the composition and structure of an exoplanet’s atmosphere. We will discuss these solutions and point out where more research or improvement of the mentioned techniques is necessary.
2.1 Exoplanet transit Transits of exoplanets in front of different stars have been observed thousands of times by now and are increasing by the month. The technique is based on the photometric observation of a star and investigation of its lightcurve for a periodically occurring dip in the received brightness or flux. The number of exoplanets that can be observed this way is limited by the inclination of their orbit with respect to the observer. Assuming a planet of mass Mp and radius Rp transiting in front of its host star with M? and R? for stellar mass and radius the orbit can be projected onto the plane of the observer. Any point of the orbit of the exoplanet has a projected distance rsky to the barycenter of the orbit. The distance r of a planet from the barycenter is given by
a(1 − e2) r = , (1) 1 + e cos f and can be projected onto the plane of the observer, such that the projected distance rsky takes the form q 2 2 rsky = r 1 − sin (ω + f) · sin i . (2)
5 a is the semi-major axis of the planet’s orbit, e the eccentricity, f the true anomaly, i the orbital inclination with respect to the projected plane of the observer, and ω is the argument of periapse. The term under the square root represents the projection of the orbit onto the√ observer’s plane. The 2 2 projected distance expressed in Cartesian coordinates takes the form rsky = X + Y , where the X axis is pointing in the direction of the descending node of the planet’s orbit. Following the conventions of Winn(2010), X can be inferred directly from orbital mechanics and the projection term in rsky:
X = −r · cos (ω + f) (3)
π Applying X = 0 we get the condition for transit mid-time ftra = + 2 − ω. Similarly the condition for π occultation of the planet by the star is given by focc = − 2 − ω. From a photometric measurement of a transit and occultation one can constrain the parameter space for ω and the eccentricity e. If possible, the transit observations can be combined with the radial velocity method to measure e directly and in this way also infer ω. With the exception of a and the inclination i, the parameter space of an orbit is covered. Winn(2010) also deduce the probability of an observable transiting planet from the penumbra and antumbra cones of a transiting planet:
R ± R 1 + e sin ω p = ? p . (4) tra a 1 − e2
The penumbra is represented here by the + sign and includes grazing orbits (not full planetary disk occulting the star) and − represents fully enclosed transits. The probability ptra to observe a transit is therefore determined by the orbital parameters a, e, and ω. Using population studies for e (Xie et al., 2016) and integrating out ω, we can estimate how many exoplanets we should be able to observe. For our purpose we do not want to have an estimate of how many exoplanets might be observable in the future (Brakensiek and Ragozzine, 2016), however, we would like to calculate the expectancy to find another Earth-twin around a Sun-like star, i. e. a G-dwarf. Plugging in the same values as for Earth in (4), we can compute the expected number of systems we must observe to find an Earth-sized planet at 1 AU around a G-dwarf. With the fraction η of systems with an Earth like planet we must observe N > 200 · η−1 systems to find a transiting Earth-twin. Using estimations for the relative abundance of G-dwarfs (Mamajek, 2016), the distribution of semi-major axis a around G-dwarfs (Han et al., 2014 (accessed July 20, 2020) (see Figure1), and the distribution of eccentricities e (Xie et al., 2016) of all exoplanets, we estimate η:
η = fa × fG × fe ≈ 0.00054, (5) 27 where fa = 1497 stands for the fraction of exoplanets found at 1 AU restricted to G-dwarf host stars, fG ≈ 0.06 the fraction of G-dwarf stars in our galactic neighbourhood and fe ≈ 0.5 the fraction of planets with eccentricity e ≈ 0. Plugging in the numbers from above we get N ≈ 3700000. This means on average out of 3700000 stars, 1 star is a G-dwarf with a transiting planet the size of Earth at 1 AU and a low eccentricity e below 0.1. Filtering the exoplanets archive (Han et al., 2014 (accessed July 20, 2020) for transiting planets around G-dwarfs with 0.8 < a < 1.2 and e < 0.1, we have not found any planets fulfilling the given constraints except 6 planets with undetermined eccentricity. Expecting to find an observable Earth twin in the near future might be ambitious to unrealistic. We will see later also that the required signal-to-noise ratio (S/N) to observe a Earth-twin around a G-dwarf is way beyond present detection limits. Concluding from this, we are forced to look for extraterrestrial life or hospitable environments for life elsewhere than what we are used to.
6 exoplanets.org | 7/20/2020
200
150
100 Distribution
50
0 0.01 0.1 1 10
Semi-Major Axis [Astronomical Units (AU)]
Figure 1: Distribution of the semi-major axis of exoplanets around G-dwarfs, pro- duced with (Han et al., 2014 (accessed July 20, 2020).
Measuring the lightcurve and observing the transit of an exoplanet we can infer two other important parameters of an eclipsing system, namely the period and the transit duration. Both are important to predict transit observability. We denote the period with P and the transit duration with Ttot. This should not be confused with the full transit duration Tfull, which stands for the time the complete planetary disk is occulting the star. The time difference from Ttot to the beginning of Tfull and vice versa is called the ingress and egress time during a transit (see Figure2). Using these properties we can predict future transits of an exoplanet in question:
tc(n) = tc(0) + P n, (6) where tc is the time of conjunction or for our purposes the transit mid-time. The uncertainty in tc(n) is mostly sensitive to the uncertainty in P. Using Gaussian error propagation we get
p 2 2 2 ∆tc(n) = ∆tc(0) + n ∆P ≈ n∆P. (7) Other errors that may influence the predictions and planning of transit observations are influences from other bodies in the observed system, which are hard to constrain. However, provide an excellent way of discovering the presence of other unobserved objects. If we also would like to observe occultations we can get another estimate for the eccentricity e and the observer’s celestial longitude ω from the time difference between transit and occultation ∆Tc P 4 ∆T ≈ 1 + e cos ω . (8) c 2 π
7 Additionally to observe exoplanets from an observatory on Earth, we must include the barycentric correction, which accounts for the light travel time from the observed system to the barycenter of the Solar System compared to the travel time to the observer. In the following expression UTC stands for universal time coordinated and if not mentioned otherwise, we always refer to UTC. BJD stands for the barycentric julian date and JD for the julian date. ~r is the position of the observer with respect to the barycenter of the Solar System, nˆ the unit vector from the observer to the observed system and c is the speed of light. Then the barycentric correction for objects outside of our Solar System can be written as (Eastman et al., 2010) ~r · nˆ BJD ' JD + . (9) UTC UTC c
The correction term is interchangeable and we can compute JDUTC , to compute our predicted transit times. This method was used to convert the transit observation times tc from the barycentric frame to the observers frame, the Paranal observatory.
If one measures the radial velocity of the host star and thereby knows the eccentricity e, the velocity semi-amplitude K? and the period P , one can estimate the mass of the planet Mp (Murray and Correia, 2010; Winn, 2010): √ M K 1 − e2 P 1/3 p = ? . (10) 2/3 (Mp + M?) sin i 2πG G is Newton’s gravitational constant. Using the transit method together with the radial velocity method the sin i degeneracy is approximately broken since in case of a transit sin i ≈ 1. However, normally 2/3 we are faced with the situation that Mp M? and we can only measure Mp/M? . To measure the mass of the planet the mass of the star must be determined by other means, such as spectral type, luminosity, and spectroscopically determined properties, such as effective temperature Teff , surface gravity g?, and metallicity [Fe/H]. The surface gravity gp at the surface of the planet can be derived by using the expression for the radial velocity semi amplitude of the star K? (Murray and Correia, 2010)
Mp 2πa/P sin i Mp 2πa/P sin i K? = √ ≈ √ (11) Mp + M? 1 − e2 M? 1 − e2 and Kepler’s third law. We obtain the surface gravity in terms of Rp/a, P, K?, and i: √ 2 Mp 2π 1 − e K? gp = G 2 = 2 . (12) Rp P (Rp/a) sin i
Measuring the photometric lightcurve during a transit we can obtain the depth of the lightcurve δtra, which provides us with further information about the observed planet star system, such as the ratio between the radius of the planet and the star k = Rp/R?: 2 Ip(ttra) δtra ≈ k 1 − . (13) I?
In case of negligible nightside emission Ip = 0 compared to the disk averaged stellar intensity I?, δtra 2 can be set equal to k , and the stellar radius R? can be computed from stellar models or approximately −1/3 by R? ∝ M? . Thereby one can also solve the degeneracy between Rp and a, since by using Kepler’s third law again we can compute a using M?. Notice that without making assumptions about the stellar properties of the host star, the presented treatment to obtain Rp is always degenerate with either R?, M?, and a, or sin i.
8 Winn(2010) present another way to approximate the stellar radius and the orbital inclination from the shape of the transit lightcurve. However, in a real observation of a lightcurve we are faced with limb darkening, which means the blocked fraction of the stellar flux decreases at the edges of the stellar disk, resulting in an increase of the received flux. To obtain information about the shape of the lightcurve such that one may determine transit, ingress and egress time, one needs to correct the lightcurves for limb darkening.
Knutson et al.(2007) propose a way to fit the lightcurves in different bandpasses simoultanously for the shape of a limb darkening corrected lightcurve, using the non-linear limb darkening law presented by Mandel and Agol(2002): 4 X n/2 Iˆ(r) = 1 − cn(1 − µ ) (14) n=1 Iˆ(r) is the specific intensity at position r on the stellar disk, where Iˆ(0) = 1 and 0 ≤ r ≤ 1, and 0 stands for the center of the√ stellar disk, while 1 stands for the outer edge of the stellar disc. r is related to µ through µ = cos θ = 1 − r2 , where θ represents the angle between the line of sight and and the line of emergence of the flux, standing perpendicular on the stellar surface. cn are the four coefficients to fit the limb darkening from stellar models or from photometric transit lightcurves. A more frequently encountered limb darkening law is the quadratic limb darkening with c1 = c3 = 0, c2 = γ1 + 2γ2 and c4 = −γ2, and the new fitting coefficients γ1, 2. The choice of using the quadratic over the full treatment is due to higher efficiency and less complexity in fitting lightcurves. Mandel and Agol(2002) present the full treatment of possible transit scenarios and the following lightcurve geometry with application of limb darkening. Knutson et al.(2007) apply the non-linear limb darkening for small planets k ≤ 0.1 which returns a set of equations to fit the limb darkening coefficients simultaneously with the decoupled parameters k and i. The calculations are lengthy and complicated and we forgo presenting them here and only show the expressions derived by Winn(2010). Another technique developed to fit the specific intensity dependent on limb darkening was proposed by Aronson and Piskunov(2018) and basically follows a similar idea as described by Mandel and Agol(2002).
The impact parameter b is depicted in Figure2, and is fundamentally related to the geometry of the orbit of the transiting planet: a cos i 1 − e2 btra = , (15) R? 1 + e sin ω where the + in the denominator must be replaced with − for occultations bocc. If we can retrieve the limb darkening corrected shape of the lightcurve, we can decouple R?/a from b using the time of ingress and egress (Ttot − Tfull) and computing the scaled stellar radius R?/a for non-grazing orbits Rp R? a and small eccentricities in terms of Ttot and Tfull:
√ 2 2 √ 2 2 (1 − δtra ) − (Tfull/Ttot) (1 + δtra ) b = 2 (16) 1 − (Tfull/Ttot) q 2 2 R π Ttot − Tfull 1 + e sin ω ? = √ (17) a 1/4 P 2 2δtra 1 − e
These expressions can be derived by computing Ttot and Tfull, also depicted in Figure2, in terms of the 2 shape of the lightcurve and the orbit and inverting for b and R?/a. Note that we have now independent expressions to compute i and R?/a. Still, to lift the degeneracy between Rp, R?, M? and a we need
9 Figure 2: Illustration of the different quantities introduced to describe the lightcurve of a transiting planet in front of its host star. T represents Tfull, Ttot would cover the times from tI until tIV , tI−II is the ingress time, and tIII−IV the egress time. b is called the impact parameter of the transit and is closely related to the orbital parameters of the transiting planet (Winn, 2010) . some input parameters for the star, either from fitted observational data, with stellar models, or from other theoretical assumptions for the star. Using photometric lightcurves together with radial velocity observations enables us to constrain the parameter space of the planet’s orbit, make good estimates for the radius of the planet Rp and the surface gravity gp and can be used as inputs for planetary climate and atmospheric models. However, the observed planetary radius Rp from a lightcurve is not well defined since cloudy or hazy atmospheres can alter the fraction of blocked stellar flux to total flux and defining the surface of a gas planet might be ambiguous. Nevertheless, the different values for Rp we may get by observing the lightcurve at different wavelengths is at the base of the idea of transmission spectroscopy of exoplanetary atmospheres, which is further elaborated in the next section.
2.1.1 Photometric lightcurve and transmission spectroscopy Using spectroscopy to probe transmission spectra of a planetary atmosphere needs some assumptions first. We define the radius of the planet Rp as the opaque part of the planet in all wavelengths, deduced from a photometric lightcurve as derived above. This means that the planetary radius Rp may be set at the top of an opaque cloud layer or at the optical depth τ ∼ 1, integrated over all wavelengths.
10 The photometric lightcurve, observed at frequency ν is defined by the fraction of the received flux to the stellar flux:
2 k αtra(t), transits Fν(t) Ip(t) f(t)ν ≡ = 1 + − 0, outside eclipse (18) Fν?(t) I? k2α (t) Iνp(t) , occultations occ Iν? where we have used the definition for the monochromatic flux (Pradhan and Nahar, 2011)
Z Z 1 Fν = Iν cos θdΩ = 2π Iνµdµ (19) ∂Σ −1
2 with Iν = Ip(t)+I? the received intensity in [ergs/(cm s str Hz)]. Here, the polar angle θ stands for the angle under which the light was emitted from the surface dA on the star. The monochromatic flux Fν 2 2 is given in [ergs/(cm s Hz)] or [W/(m s Hz)]. αtra and αocc describe the geometry of the area blocking the star light during a transit or the area adding to emitted light from the planet during occultations and out of transit. Mandel and Agol(2002) give the full analytic expressions for all possible cases that might be encountered doing transit observations. The different cases are derived from geometrical considerations.
Measuring the radius of the planet Rp at a certain wavelength might return a different radius than obtained from photometry. We rewrite the new radius as Rν = Rp +δRν depending on the wavelength, respectively the frequency ν. Then the part adding to the new occulting planetary disk can be described 2 by an annulus 2πRpδRν + πδRν with thickness δRν. The thickness of the annulus δRν can be related to the pressure scale height of the atmosphere observed at wavelength ν via the optical depth τν at ν. The optical depth is the integrated opacity along a path C weighted with the density ρ: Z τν = κν(s)ρ(s)ds. (20) C
The opacity can be related to the atomic absorption coefficient αν: κν = ανn, where n is the relative number density of a certain species. The expression for the absorption coefficient α is often seen to describe a single absorption line, whereas the opacity κ is used in general for radiative transfer with its continuous contributions from absorption and emission as well. Depending on the geometry and the physics we are investigating, κν or αν can be dependent on the integration path C. We will now show how the thickness of the annulus δRν can be related to the absorption at ν. For more details, see section 2.3.3 and (Gray, 2005).
Assume we have an atmopshere consisting of a single species. For a grazing light ray traveling through the annulus of the semi-transparent atmosphere at an impact parameter xν, we can compute the tangential optical depth τν according to Webb and Wormleaton(2001) and Brown(2001). The geometry is depcited in Figure3. In this derivation the geometry does not take into account any effects of refraction.
The tangential optical depth in this geometry is given by Z ∞ (r + Rp)dr τν(x) = 2 κνµn(r) . (21) p 2 2 x (r + Rp) − (x + Rp)
11 Figure 3: Geometry for a grazing light ray traveling through the semi-transparent atmosphere of a planet at impact parameter x, maximum atmospheric height H and planetary radius Rp (Webb and Wormleaton, 2001) . with the molecular mass µ of the species in question. In case of hydrostatic equilibrium we can use a barometric law for the number density of the species with height r: k T n(r) = n e−r/h, h = B , (22) 0 µg where the pressure scale height h is dependent on the surface gravity g, under the assumption of H Rp, such that g does not change much with altitude, the temperature T of the atmosphere, and kB is the Boltzmann factor. This may be generalized to non-isothermal atmospheres with different temperature dependencies in different layers and non-homogeneous processes such as strong irradiation effects on one side of the planet in case of tidally locked planets. Assuming the absorption coefficient to be relatively constant along a tangential line of sight, the integral above can be evaluated using the modified Bessel function K1((x + Rp)/h) (Abramowitz and Stegun, 1965): (x + R ) τ (x) = 2n κ (x + R )eRp/hK p (23) ν 0 ν p 1 h In the absence of absorption lines we can estimate the overall pressure scale height h to one order of magnitude accuracy (Benneke and Seager, 2012; Miller-Ricci et al., 2008) by measuring the impact parameter x at two different wavelengths under the assumption τ(x1,2) = τ(δR1,2) = 1. Equating both expressions and solving for h, we get x1 − x2 h = q , (24) ln α1 x1 α2 x2 where κ1 and κ2 are the opacity or scattering cross section at the two different wavelengths λ1 and λ2 respectively.
12 To investigate spectral features, we need to relate the measured impact parameter x to the decrease of the stellar flux by the semi-transparent atmosphere: δF?,ν. The specific photometric light curve fν (18) in terms of the blocked stellar flux δF?,ν can be expressed as:
δF?,ν fν = 1 − , (25) F?,ν where the second term is given by (Brown, 2001) Z ∞ δF?,ν 1 0 0 0 = 2 2π(Rp + x )[1 − exp (−τν(x ))]dx , (26) F?,ν πR? 0 and one integrates over all impact parameters. τ is the integral along every line of sight at each impact parameter. The factor 2π arises from integrating around the whole annulus visible by the observer, neglecting effects of rotation of the planet, the star or wind in the atmosphere of the planet. Following the treatment of Brown(2001) and the method to measure the transmission spectrum of an exoplanetary atmosphere by Aronson and Waldén(2015), we may equate the expression for δFν with ia · P , where the latter describes the spectrum of the semi-transparent atmosphere expressed in terms of the measured impact parameter xν: Z ∞ 2 2 0 0 0 ip&a(·2πRpxν + xν + Rp) = 2π(Rp + x )Iν(µ(φ))[1 − exp (−τν(x ))]dx (27) 0 Z ip&a(φ) = Iν(µ(φ))dµ(φ) (28) Ap&a(φ) where ip&a denotes the specific intensity from the star covered by the planet and its atmosphere. For I(µ(φ)) one can use a quadratic Limb darkening law or spectral models. The dependency of µ on the orbital parameter φ can be looked up in Figure4. It describes the relation between the position of the planet on its orbit around the star to the observer’s view of the position of the planet. What we have not treated in our derivation of δF/F is the influence of the refractive index and the night side emission of the planet. Brown(2001) point out that the path of a light ray through the atmosphere is not perfectly straight. However to compute the complete obscurred area by the atmosphere of the planet they explain that the minimum height of a tangential light ray not altered by refraction is relevant. They do not discuss night side emission, and Aronson and Waldén(2015) do not discuss the influence of refraction or night side emission and neglect it. Using the previous arguments we may now present their novel retrieval method.
2.2 Atmospheric retrieval methods for transmission spectroscopy 2.2.1 A novel planetary model-independent transmission spectrum recovery method Aronson and Waldén(2015) describe a method to recover transmission spectra of exoplanetary atmo- spheres in the near-infrared from Earth based observatories with high-resolution spectrometry. The method should be applied on data from CRIRES+ in the future. The main idea is that even in presence of strong telluric absorption CRIRES+ will be able to measure in between the telluric lines and observe features shifting with time due to a time dependent Doppler shift, while the planet moves across the surface of the star. The full signal received by a high-resolution spectrograph can be expressed as received flux per exposure normalized to the exposure time: 2 ˜ π(Rp + xν) Sn(ν, t) = F (ν, t, vs) − ip&a(φ, vs) 2 · f(vs, vp) · T (ν, t), (29) πR?
13 with F (ν, t) the stellar flux and T (ν, t) the telluric transmission and f(vp, vs) a Doppler correction factor for the movement of the planet with respect to the star, the rotation of the star, the movement of the Earth, and the rotation of the Earth. vs and vp are the velocities of the star and the planet relative to Earth depending on their present position during the transit. They are dependent also on the orbital parameters φ at t and the shape of the orbit (a, ω, e, i). The geometry is depicted in Figure4. Additional Doppler shifts can be introduced by planetary winds across the terminator (the only part of the atmosphere observable with transmission spectroscopy) and the rotation of the planet, however, a reasonable assumption is that the planets are tidally locked and so the influence of the planets rotation is small. On the other hand, horizontal winds arise from the strong stellar irradiation onto a tidally locked planet and dissipate some of the energy to the night side. Snellen et al.(2010b) used CRIRES, the predecessor of CRIRES+, to measure the high altitude winds of the gas giant HD 209458 b with transmission spectroscopy and a cross correlation technique that we will present in the next section 2.2.2. They showed that these winds can reach speeds of several 100 ms−1, close to the sound speed of the transported gases and can therefore alter the spectrum through line broadening and overall wavelength shifts significantly. To obtain the complete received signal one needs to convolve the spectrum with the instrument profile Γ. It can be approximated by a Gaussian with full-width at half maximum equal to the spectral resolution of the spectrograph and a noise function for each pixel. The 1 noise function follows a Poisson distribution with a width of λ = S/N , where S/N is the signal-to-noise ratio. The signal is then finally multiplied with a continuum normalization η. Applying these factors we get h i Sn,syn(ν, t) = S˜n(ν, t) ⊗ Γ · (1 − N(ν, φ)) · η. (30)
The continuum normalization is necessary due to the limited wavelength range of high-precision spectrographs and generally to correct for any in-
� strumental changes in the signal over time. Con- � � ventionally the continuum is normalized via fit- ting of the recovered spectrum or removing all broadband signal contributions and computing the relative line strengths to the total stellar flux (de Kok et al., 2013). Aronson and Waldén (2015) use sensitivity curves, derived from mea- suring the Solar spectrum before and after the science exposure to compare with a theoretical Solar spectrum and normalize the measured spec-
� &