Effects of Thermal Emission on the Transmission Spectra of Hot Jupiters

The atmosphere on the dayside of a highly irradiated close-in gas giant (also known as a hot Jupiter) absorbs a significant part of the incident stellar radiation which again gets re-emitted in the infrared wavelengths both from the day and the night sides of the planet. The re-emitted thermal radiation from the night side facing the observers during the transit event of such a planet contributes to the transmitted stellar radiation. We demonstrate that the transit spectra at the infrared region get altered significantly when such re-emitted thermal radiation of the planet is included. We assess the effects of the thermal emission of the hot Jupiters on the transit spectra by simulating observational spectroscopic data with corresponding errors from the different channels of the upcoming James Webb Space Telescope. We find that the effect is statistically significant with respect to the noise levels of those simulated data. Hence, we convey the important message that the planetary thermal re-emission must be taken into consideration in the retrieval models of transit spectra for hot Jupiters for a more accurate interpretation of the observed transit spectra.


INTRODUCTION
Transit spectroscopy is an essential tool for probing into the upper atmospheres of the close-in exoplanets (Charbonneau et al. 2002;Sing et al. 2016;Sengupta, Chakrabarty & Tinetti 2020). As an exoplanet transits across its parent star, a fraction of the starlight is transmitted through the planet's upper atmosphere. Light at specific wavelengths is preferentially absorbed, depending on the atmospheric chemical composition and physical properties. An accurate interpretation of the observed transit spectra, however, requires a self-consistent theoretical model that must incorporate the physical and chemical properties of the planetary atmospheres in sufficient detail, including scattering albedo (e.g. de Kok & Stam 2012;Sengupta, Chakrabarty & Tinetti 2020) and thermal re-emission from the night side (this work). Theoretical models for the planetary transit spectra presented by various groups are described in articles by Brown (2001) All of these models have only considered the effect of the total extinction coefficient and ignored the effect of scattering albedo in calculating the transmission depth at different wavelengths. Recently Sengupta, Chakrabarty & Tinetti (2020) demonstrated that the diffused reflection and transmission due to scattering of the transmitted starlight by the atoms and molecules affect the transmitted flux and hence, the transit depth significantly in the optical. Therefore, modeling the transit spectra correctly and consistently requires one to solve the multiple scattering radiative transfer equations. However, Sengupta, Chakrabarty & Tinetti (2020) show that effect of scattering albedo on transmission spectra is less than 10 ppm in the infrared region.
On the other hand, because of the extreme proximity to their parent stars, many gas giants are extremely hot. As they are tidally locked with their host stars, the dayside of such a planet is so hot due to the intense irradiation that its equilibrium temperature for zero albedo, T eq0 can reach as high as ∼4000 K (Gaudi et al. 2017). The heat is redistributed to the night side by the advection process and the night side of the planet facing the observer during transit also becomes hot depending on the heat re-circulation efficiency, ǫ of the planetary atmosphere. The average temperature T n of the night side of the planet can be estimated by the relation (Cowan & Agol 2011;Keating & Cowan 2017) T (1) where A B is the Bond albedo of the atmosphere. We define T n0 as the night-side average temperature for zero Bond albedo. Even with a small value of ǫ, the average night side temperature T n can be quite high if T eq0 is high enough. Some previous studies provide an estimation of the night-side temperature of a few exoplanets. For example, Keating & Cowan (2017) report T n = 1080 ± 11 K for WASP-43 b; Demory et al. (2016) report a night-side brightness temperature of 1380 ± 400 K for 55 Cancri e based on their observation in the 4.5-µm channel of the Spitzer Space Telescope Infrared Array Camera (IRAC); Arcangeli et al. (2019) report T n ≤ 1430 K for WASP-18 b at 3-σ level; etc. Such a hot region facing the observer would emit radiation in the infrared wavelengths which should be added up with the transmitted stellar radiation. Thus the transmitted flux would be affected by the re-emission of hot planets in the longer wavelengths.
In this paper, we demonstrate the effect of thermal re-emission on the transit spectra of the hot Jupiters of different size and surface gravity and with equilibrium temperature ranging from 1200K to 2400K. In Section 2 we provide the formalisms for calculating the transit depth with and without thermal re-emission. Section 3 outlines the detailed procedure followed to calculate the transmitted and the re-emitted flux from the hot-Jupiters. In Section 4 we discuss the 1D pressure-temperature grids calculated and the databases adopted for the calculation of abundance and absorption and scattering opacity for the modeling of the atmospheres of the hot-Jupiters. In Section 5 we discuss the results from all the case studies as well as the results from the testing of detectability of the effect of thermal emission on the transit spectra of hot Jupiters by simulating observational transit spectroscopic data to be observed from the upcoming James Webb Space Telescope (JWST) using the open-source Pandexo code (Batalha et al. 2017) available in the public domain 1 . In the last section, we conclude the key points.

THE TRANSIT DEPTH WITH PLANETARY THERMAL RE-EMISSION
The transit spectra of an exoplanet are expressed in term of transit depth D(λ) which is the difference in stellar flux during out of transit and during the transit of the planet and normalized to the unblocked stellar flux. When the planetary radiation does not contribute to the stellar flux, it can be written as (e.g. Kempton et al. 2017;Sengupta, Chakrabarty & Tinetti 2020) where, D N E is the transit depth with no thermal radiation from the transitting planet, R P A is the sum of the base radius R P and the atmospheric height of the planet, F in and F out are the in-transit and out-of-transit stellar flux respectively. In case of pure transmission when thermal emission of the planet is ignored, F out is the stellar flux F * . F P is the portion of the stellar flux which is transmitted through the upper atmosphere of the planet and undergoes absorption and scattering through the medium. The above equation can also be written as: For a transiting planet with a hot night side facing the observer, the re-emitted radiation flux F T h is added to the observed fluxes. Hence, the transit spectra including the effect of the re-emission from the planet can be expressed as: A transit spectrum is usually produced from observation by calculating the transit depth at different wavelengths normalized with respect to the baseline flux, i.e., the flux observed right before ingress or immediately after egress.
One of the methods is to calculate the wavelength-dependent transit depth from the light curves at different wavelength bins extracted from the time-series spectra of the host stars observed during a transit event using spacebased instruments like HST+STIS, HST+WFC3, etc. (Charbonneau et al. 2002;Gibson et al. 2012;Berta et al. 2012;Deming et al. 2013;Sing et al. 2016, etc.) or ground-based instruments like VLT+FORS, VLT+FORS2, Magel-lan+MMIRS, GEMINI-N+GMOS, GTC+OSIRIS etc (Bean et al. 2011;Stevenson et al. 2014;Nikolov et al. 2016;Pallé et al. 2016;Huitson et al. 2017, etc.). The wavelength-dependent transit depth can also be calculated directly from the photometric light curves at different wavelength bands observed using ground-based instruments or space-based instruments like Spitzer Space Telescope Infrared Array Camera (IRAC) etc. during a transit event (Tinetti et al. 2007;Sing et al. 2016). In either case, the normalized wavelength-dependent transit depth (equivalently, the wavelength-dependent radius) can be more accurately modeled with the incorporation of the planetary thermal emission F T h . However, the impact on the accuracy by omitting the effect of F T h might be, in some cases (e.g., at wavelength shorter than 2 µm), small compared to the uncertainties in the data. We discuss the significance of the effect of F T h with respect to the uncertainties in the observational data elaborately in Section 5.
Finally, Equation 2 and Equation 4 gives As evident from Equation 5 the contribution from the thermal re-emission of the planet reduces the transit depth.

CALCULATIONS OF TRANSMISSION AND EMISSION FLUX
To calculate the transmitted stellar flux F P that passes through the atmosphere of a hot Jupiter, we first, calculate the reduced stellar intensity that suffers absorption and scattering in the planetary atmosphere and then integrate over the angle subtended by the annular region of the atmosphere. Sengupta, Chakrabarty & Tinetti (2020) have shown that an accurate approach to calculate the reduced intensity is to solve the multi-scattering radiative transfer equations that incorporate the diffused reflection and transmission of radiation due to scattering. The radiative transfer equations including diffused reflection and transmission for a plane-parallel geometry can be expressed as (e.g. Chandrasekhar 1960; Sengupta, Chakrabarty & Tinetti 2020) where I(τ LOS , µ, λ) is the specific intensity of the diffused radiation field along the direction µ = cos θ, θ being the angle between the axis of symmetry and the ray path, F * is the incident stellar flux in the direction −µ 0 , ω is the albedo for single scattering i.e. the ratio of scattering co-efficient to the extinction coefficient, p(µ, µ ′ ) is the scattering phase function that describes the angular distribution of the photon before and after scattering and τ LOS is the optical depth along the line of sight to the observer. The detail method for calculating τ LOS is given in Sengupta, Chakrabarty & Tinetti (2020). We adopt the Rayleigh phase function (e.g. Chandrasekhar 1960) for cloudfree atmospheres. For cloudy atmospheres, we use the Mie phase function. This is elaborated in Section 5.5. The planetary thermal emission is important only at longer wavelengths where the effect of scattering albedo is negligible, less than 10 ppm (Sengupta, Chakrabarty & Tinetti 2020), even in the presence of atmospheric clouds. Therefore, in the present work, we have ignored the effect of scattering albedo in the calculation of F P .
On the other hand, to calculate the flux of the planetary thermal radiation F T h , we treat the planet as a self-luminous object and solve the radiative transfer equations in the following form: where, I(τ, µ, λ) denotes the specific intensity of the thermal radiation field along the direction µ, τ is the optical depth in the radial direction, B(τ, λ) is the Planck function corresponding to the temperature of the atmospheric layer with optical depth τ at a particular wavelength. The above radiative transfer equations are solved by using the Discrete Space Theory method (Peraiah & Grant 1973). The numerical method is described in Sengupta, Chakrabarty & Tinetti (2020) and in Sengupta & Marley (2009). T eq0 =1200 K, g=30 m/s 2 T eq0 =1600 K, g=15 m/s 2 T eq0 =1600 K, g=30 m/s 2 T eq0 =1600 K, g=60 m/s 2 T eq0 =2000 K, g=15 m/s 2 T eq0 =2000 K, g=30 m/s 2 T eq0 =2000 K, g=60 m/s 2 T eq0 =2000 K, g=100 m/s 2 T eq0 =2400 K, g=30 m/s 2 Temperature (K) Pressure (bar) Figure 1. Pressure-Temperature profiles for different Teq0 and g adopted in this work. For all the cases, the effect of TiO and VO is included.

PRESSURE-TEMPERATURE GRIDS AND THE OPACITY AND ABUDANCE DATA
The atmospheric pressure-temperature structure plays an important role not only in determining the optical depth of the medium but also in estimating the thermal radiation of the planet. In order to calculate the pressure-temperature profiles, we use the FORTRAN implementation of the analytical models of non-Grey irradiated planets provided by Parmentier & Guillot (2014); Parmentier et al. (2015). This code is available in public domain 2 . It uses the functional form for Rosseland opacity derived by Valencia et al. (2013). The Rosseland opacities of Freedman, Marley, & Lodders (2008) are adopted in the derivation. For T eq0 ≥ 1700K, the pressure-temperature (P-T) profiles show temperature inversion when TiO and VO are included. Previous studies, including, e.g., Sengupta, Chakrabarty & Tinetti (2020) find that the hot-Jupiters are almost opaque to the transmitted flux at pressures level below 1 bar. Therefore, in the present work, we have considered the P-T profiles up to 1 bar pressure level such that the base radii R P of the hot Jupiters considered in the calculations of the transmission spectra are located at 1 bar pressure level.
However, the thermal radiation of a hot planet emerges from a deeper layer of the atmosphere and so we have considered the P-T profiles of hot Jupiters with T eq0 ranging between 1200K to 2400K and the surface gravity over a range of 15-100 m/s 2 at 10 bar pressure level. Figure 1 shows the P-T profiles for all the case-studies with different T eq0 .
For all the calculations we have adopted solar metallicity and solar system abundance for the atoms and the molecules present in the atmospheres. We have considered 28 molecular and atomic species as mentioned in Sengupta, Chakrabarty & Tinetti (2020). The abundance for all these atoms and molecules has been calculated using the abundance database given in the open-source package "Exo-Transmit" (Kempton et al. 2017) available in the public domain 3 . In this package, the abundances for major atmospheric constituents as a function of temperature and pressure are calculated based on the solar system abundances of Lodders (2003). Also, to calculate the absorption and scattering coefficients we have used the opacity database from the same package. These opacities are based on the molecular databases of Freedman, Marley, & Lodders (2008) and Freedman et al. (2014). The line lists used to generate the molecular opacities are also tabulated by Lupu et al. (2014). Both the abundance and opacity databases are available over a broad grid of pressure and temperature from which we have interpolated to our particular P-T profiles. We have adopted the equation of state (EOS) for rain-out condensation. We have also calculated the opacity due to atmospheric clouds comprising mainly of amorphous Forsterite by using Mie theory (see Section 5.5).
Modeling an atmosphere with a high day-night temperature contrast is not straightforward because 1D pressuretemperature (P-T) profile may not be adequate in such a scenario and therefore one requires a 3D P-T mapping using Global Circulation Model (GCM) or a limb-averaged P-T profile (Kataria et al. 2015(Kataria et al. , 2016Evans et al. 2018). In order to avoid such complications, we have assumed a hot-Jupiter atmosphere with ǫ = 1 (T n0 = T eq0 ) i.e. a globally averaged P-T profile in all the cases as our main motive is to demonstrate the effect of the night-side temperature T n . This corresponds to f = 0.25, where f is the flux parameter as defined in Burrows, Sudarsky, & Hubbard (2003) and also, used in Guillot (2010) The night side of the planetary atmosphere may also get warmed up by the release of the internal energy characterized by the internal temperature (T int ) and may cause thermal radiation (self-emission) but this radiation should be insignificant as compared to the thermal re-emission from planet older than 100 Myr (Burrows et al. 1997) in age.

CASE STUDIES: SIMULATION AND TESTING OF DETECTABILITY
We checked the detectability of the effect of thermal emission on the transit spectra by using the simulated data. We also investigate the extent to which the transit spectra D E depend on the planetary properties such as T n , R P /R * , the atmospheric clouds and the surface gravity as well as the spectral types of the host stars. The following subsections describe the results from these case studies.

Detectability with JWST
To understand the significance of the effect of thermal emission from transiting hot Jupiters on the transit spectra, we need to compare the difference between D N E and D E with the noise levels of the actual observed data. For that purpose, we have focused on the observing capability of the upcoming James Webb Space Telescope (JWST) as this mission is going to be at the forefront of exoplanet characterization. Considering the contribution of thermal emission from hot Jupiters to the transit spectra, we use our model calculations for D E and, simulate some observational transit spectra with error-bars that can be observed using the IR instruments of JWST. The simulation is done by using the open-source code Pandexo (Batalha et al. 2017). A host star of spectral type G2V with J-mag=8 and a saturation level equal to 70 % of the full well potential are considered for this simulation. The spectra with noiselevels with a resolution of R ∼ 50 are calculated by combining the data over 4 observed transit events, each with a duration of 2 hours (T 14 ). We have considered the instrument modes viz. NIRSpec G140M and NIRSpec G235M for the wavelength region 1-3 µm and NIRSpec G395M and MIRI LRS (slitless) for the wavelength region 3-8 µm (see Table 1 of Batalha et al. (2017)) and the corresponding simulated transit depth and noise levels are denoted by D G140M ± σ G140M , D G235M ± σ G235M , D 395M ± σ G395M and D LRS ± σ LRS respectively. For each pair of instrument modes we present the comparison of the models (D N E vs D E ), constructed using different sets of planetary paramters, with the simulated data along with their residuals in Figure 2-5. In Figure 2 and Figure 4 the model parameters are (i) T n = 1600K, R P = 1 R J , R * = 1 R ⊙ , (ii) T n = 2000 K, R P = 1 R J , R * = 1 R ⊙ , (iii) T n = 1600 K, R P = 1 R J , R * = 0.8 R ⊙ , and (iv) T n = 2000 K, R P = 1R J , R * = 0.8R ⊙ . In Figure 3 and Figure  Transit epth (%) T n =1600 K, R P =1 R J , R * =0.8 R 1 T n =1600 K, R P =1 R J , R * =0.8 R 1 χ 2 NE = 2.02 χ 2 E = 1.03 T n =2000 K, R P =1 R J , R * =0.8 R 1 T n =2000 K, R P =1 R J , R * =0.   the models (top sub-panels) and the mean of the ratio of (D N E − D E ) to the 1-σ noise-levels of the above modes (viz. σ G140M , σ G235M , σ G395M , and σ LRS ). These figures show that the difference between D N E and D E for R P /R * = 1 is of no or extremely low significance (< 20 ppm and < 4-σ). Also, at wavelengths up to 2 µm, the difference between D N E and D E is significant (> 25 ppm and ∼ 22-σ) only for T n = 2000 K, R P = 1.4R J , R * = 0.8R ⊙ . For all other combinations of T n , R P , and R * , the difference is of no or low significance (< 25 ppm and < 4-σ). However, for wavelengths longer than 2 µm, we find that the difference between the models increases with increasing T n and R P /R * and reaches up to 500 ppm (330-σ) for T n = 2000 K, R P = 1.4R J , R * = 0.8R ⊙ . The chi-square values shown in these figures imply that for higher values of T n and R P /R * and wavelength 2 µm, the simulated data are fitted well with the model transit spectra only when  Figure 4. Simulated observational data of transit depth with emission combined over 4 observed transit events using instrument modes NIRSpec G395M and MIRI LRS (slitelss) of JWST with error-bars, viz. DG395M ± σG395M (red) and DLRS ± σLRS (blue) respectively, are shown in the top panels assuming a G2V host star with J-band magnitude = 8. Models without and with thermal emission (DNE in magenta and DE in green respectively) for different values of Tn and R * with RP = 1 RJ are shown with corresponding chi-square values, keeping g fixed at g=30 m/s 2 . The bottom panels show the difference between the model without emission and the simulated observational data (red and blue triangles with magenta error-bars) as compared to the difference between the the model with emission and the same simulated observational data (red and blue circles with green error-bars). Also, the mean of the ratio of the difference between the two models to the noise levels of the two modes, viz.  the instruments on-board JWST, the effect of thermal emission from hot Jupiters must be taken into consideration in the retrieval model for transit spectra.

Host Stars of Different Spectral Types
The transit depth with no planetary thermal emission, D N E , is absolutely independent of the stellar spectrum, as evident from Equation 6. The factor FP F * solely depends on the physical and chemical properties of the planetary atmospheres and provides only the reduction in the stellar flux due to absorption. However, Equation 7 suggests that, when planetary emission is included, the transit depth, D E , for planets with the same T n , becomes dependent on the flux of the host star. Consequently, it depends on the spectral types of the host stars. Figure 6 displays the difference between D N E and D E for planets with R P /R * = 1.4R J /R ⊙ (∼0.144), g=30m/s 2 and T n =1600K orbiting stars of spectral types F5V, G5V, K5V, and M5V. Also, the 1-σ noise levels of the JWST instrument modes NIRSpec G140M, NIRSpec G235M, NIRSpec G395M, and MIRI LRS (slitless) for number of observed transits equal to 2 and 4 and host stars with J-band magnitude of 8 and 10 are shown in this figure. The model spectra show that the transit depths D E for stars of different spectral types differ significantly at wavelengths longer than 3 µm. The cooler the host stars are, the more significant the difference is between the models with respect to the noise levels, as evident from Equation 5.  Figure 6. Difference between the models of transit depth without and with thermal emission from the hot-Jupiters with Tn = 1600 K, g=30 m/s 2 , RP /R * = 1.4RJ /R⊙ (∼0.144), orbiting around stars of different spectral types. Transit depth without emission is independent of the host star spectral type. The 1-σ noise-levels are shown in dashed lines from left to right for the JWST channels NIRSpec G140M, NIRSpec G235M, NIRSpec G395M, and MIRI LRS (slitless) respectively. The red and black dashed lines correspond to noise-levels for number of observed transits equal to 2 and 4 respectively. The thick and thin dashed lines correspond to noise-levels for host stars with J-band magnitude of 8 and 10 respectively.

Average Night-Side Temperature and Planetary Size
In order to investigate the effect of T n , we calculate the difference between D N E and D E at different values of T n0 e.g., 1200K, 1600K, 2000K and 2400K. Figure 13 of Parmentier & Guillot (2014) demonstrates that the Bond albedo of planets with high equilibrium temperature is extremely low (< 0.01) for solar composition. Consequently, from Equation 1 it follows that, T n ≈ T n0 . These values of T n correspond to atmospheric scale heights of 138 km, 186 km, 254 km and 387 km respectively. Figure 7 shows the difference between D N E and D E for these values of T n and for different values of R P and R * e.g., (i) R P = 1R J , R * = 1R ⊙ , (ii) R P = 1R J , R * = 0.8R ⊙ , (iii) R P = 1.4R J , R * = 1R ⊙ and (iv) R P = 1.4R J , R * = 0.8R ⊙ . It is clear from the figure that with the increase in T n , the difference between D N E and D E increases because the thermal re-emission from the planet increases. The factor R P /R * also strongly dictates the significance of the difference with respect the noise-levels. This happens due to the fact that, with increasing R P /R * , the ratio of the thermal luminosity of the planet to the luminosity of the host star increases. Obviously, for a fixed planetary radius, the difference in D N E and D E increases with the decrease in the size of the host star.
Also, the 1-σ noise levels of the JWST instruments NIRSpec G140M, NIRSpec G235M, NIRSpec G395M and MIRI LRS (slitless) for number of observed transit equal to 2 and 4 and host stars with J-band magnitude of 8 and 10 are shown in Figure 7. This helps us comprehend the significance of the difference between the models with respect to the noise levels for different numbers of observed transits and different host star J-band magnitudes. However, it can be safely ascertained that for higher values of R P /R * (see bottom 2 panels of Figure 7) and for T n > 1200 K, the deviation of D E from the standard model of transmission spectra (D N E ) at wavelength beyond 2 µm is so significant that observed transit spectra can be misinterpreted by the standard model by 10-300 σ (representing a change up to 0.5% in transit depth).

Surface Gravity
We have calculated D N E and D E for g=15, 30, 60 and 100 m/s 2 as shown in Figure 8 for a fixed values of T n =2000K, R P = 1.4R J and R * = 1R ⊙ . These values of g correspond to atmospheric scale heights, estimated by using T n =2000K, of 508 km, 254 km, 127 km and 76 km respectively. We find that with increasing g the transmission flux F P decreases. However, the value of g has almost no effect on the calculation of F T h and hence, the difference between D N E and D E is almost independent of g, the surface gravity of the planet.

Atmospheric Clouds
Clouds and hazes are a ubiquitous feature in the planetary atmospheres. For hot exoplanets, silicates may condensate in the upper atmosphere. Gas giant planets with comparatively cooler night sides can have thick atmospheric clouds that may affect the spectra in the optical and near infrared wavelength region. However, at higher day or night temperatures, clouds may either completely evaporate or may form a thin layer of haze in the uppermost atmosphere. As discussed in the previous subsection, the effect of thermal emission is significant only at wavelengths longer than 2µm and for a night-side temperature T n > 1200 K. Therefore, even in the presence of a thin layer of haze, we don't expect the transmission spectra of planets having T n > 1200 K to be affected in the infra-red region where thermal re-emission is important. Nevertheless, we have investigated the effects of the thin clouds on the transmission spectra of a planet with T n = 1600 K. For this purpose, we have considered a simple model for thin haze in the uppermost atmosphere. The formalism is adopted from Griffith, Yelle, & Marley (1998); Saumon et al. (2000). We consider grains of amorphous Forsterite (Mg 2 SO 4 ) of mean diameter 0.5 µm as the dominant constituent of the cloud located within a thin region of the atmosphere bound by a base and a deck. Within this region, the sizes of the particles follow a log-normal distribution and the vertical density distribution of the cloud particles follows the relation n(P ) = n 0 P P 0 where, n(P) is the number density of cloud particles at pressure level P, P 0 is the pressure at the base radius, and n 0 is a free parameter with the dimension of number density. The details of the model adopted can be found in Sengupta, Chakrabarty & Tinetti (2020). The deck and base of the haze are fixed at 0.1 Pa and 100 Pa pressure levels respectively. We use the Mie theory of scattering to calculate the wavelength-dependent scattering coefficients, extinction coefficients, and the phase functions at different pressure atmospheric depth (van de Hulst 1957; Hansen & Travis 1974;Fowler 1983;Bohren & Huffman 1983;Sengupta & Marley 2009;Sengupta, Chakrabarty & Tinetti 2020, etc.). Figure 9 shows D N E (in top panel) and the difference between D N E and D E (bottom panel) for cloud models with n 0 = 1000 cm −3 and n 0 = 5000 cm −3 . We compare the results with that of a cloud-free atmosphere. Although the transit depth without thermal emission, D N E , can alter depending on the cloud structure and opacity, the difference between D N E and D E does not change much as the emission flux (F T h ) is not affected significantly in the infra-red 0 20 40 60 80 100 120 140 160 R P = 1R J , R * = 1R ⊙ R P = 1R J , R * = 1R ⊙ R P = 1R J , R * = 1R ⊙ R P = 1R J , R * = 1R ⊙ T n =2400⊙K T n =2000⊙K T n =1600⊙K T n =1200⊙K R P = 1R J , R * = 0.8R ⊙ R P = 1R J , R * = 0.8R ⊙ R P = 1R J , R * = 0.8R ⊙ R P = 1R J , R * = 0.  Figure 7. Difference between the models of transit depth without and with planetary thermal emission for different values of Tn, RP and R * , keeoing g fixed at 30m/s 2 . The 1-σ noise-levels are shown in dashed lines from left to right for the JWST channels NIRSpec G140M, NIRSpec G235M, NIRSpec G395M and MIRI LRS (slitless) respectively. The red and black dashed lines correspond to noise-levels for number of observed transits equal to 2 and 4 respectively. The thick and thin dashed lines correspond to noise-levels for host stars with J-band magnitude of 8 and 10 respectively.
region by the presence of cloud. Figure 9 shows that the difference between D N E and D E does not change at all with clouds. Hence, the clouds do not play an important role in determining the transit depth at the infra-red wavelength region of ultra-hot Jupiters.

CONCLUSION
We demonstrate that the effect of thermal re-emission from the night side of hot Jupiters on the transit spectra can be significant at the infrared wavelength region if the equilibrium temperature of the planet is higher than about 1200K . Top -Models of transit depth without and with thermal emission for different values of g and for Tn=2000K, RP =1.4RJ and R * =1R⊙. Bottom -Difference between the above models for each value g which shows no dependance on g. The 1-σ noise-levels are shown in dashed lines from left to right for the JWST channels NIRSpec G140M, NIRSpec G235M, NIRSpec G395M and MIRI LRS (slitless) respectively. The red and black dashed lines correspond to noise-levels for number of observed transits equal to 2 and 4 respectively. The thick and thin dashed lines correspond to noise-levels for host stars with J-band magnitude of 8 and 10 respectively. and if the planet is large enough in size such that R P /R * > 0.1. The contribution of planetary thermal emission to the transit spectra can significantly exceed the total noise budget (photon noise plus readout noise) of the IR instruments on-board the upcoming JWST that will perform transit spectroscopy. Hence, a retrieval model that does not include planetary thermal emission would overestimate the transit depth and thus can lead to a wrong interpretation of the planetary properties of the hot Jupiters. Therefore, for a consistent and accurate interpretation of the observed transit spectra, it is essential to include the diffused reflection and transmission due to scattering in the optical and near infrared wavelength region and the thermal re-emission at the near and mid infrared region of hot gas giant planets. Both need the solutions of the multi-scattering radiative transfer equations.
We thank the reviewer for several valuable comments and suggestions. No cloud Clouds with n 0 =1000 cm −3 Clouds with n 0 =5000 cm −3 Figure 9. Top -Models of transit depth with thermal emission with different cloud abundance as well as without any cloud for Tn=1600K, g=30 m/s 2 , RP =1.4RJ and R * =1R⊙. Bottom -Difference between the models without and with thermal emission from the planets for the above cloud abundances as well as no cloud.