Paper

New observations and transit solutions of the inflated exoplanets HAT-P-40b and HAT-P-51b

, and

© 2021 National Astronomical Observatories, CAS and IOP Publishing Ltd.
, , Citation Diana P. Kjurkchieva et al 2021 Res. Astron. Astrophys. 21 044 DOI 10.1088/1674-4527/21/2/44

1674-4527/21/2/044

Abstract

We present high-precision photometric observations of the transiting exoplanets HAT-P-40b and HAT-P-51b by the Rozhen 2-m telescope. The newly-observed transit of HAT-P-40b is the first one with a complete curve. The orbital periods of the two targets were improved. We modeled the observed transits and found bigger stellar radii than those derived from the stellar models. The planet radii of HAT-P-40b and HAT-P-51b obtained from our transit solutions are bigger than the values calculated by the empirical relations for Jupiter-mass and Saturn-mass planets respectively. Their values reveal the highly-inflated nature of the two targets, especially that of HAT-P-51b. We established that the best transit solutions correspond to a quadratic limb-darkening law. The fitted limb-darkening coefficients of HAT-P-40 are close to the theoretical ones while those of HAT-P-51 are a little different. The precise astrometric Gaia distances of the two targets are smaller by 6%-7% than the calculated values from the stellar models. We propose the Gaia distances to be used for improvement of the stellar models as well as for more reliable calculation of the parameters of the known exoplanets.

Export citation and abstract BibTeX RIS

1. Introduction

Transiting extrasolar planets (TEPs) provide a possibility to determine their orbits, physical properties and internal structure.

Statistics reveal that the majority of TEPs have close orbits (a < 0.1 AU). Their circularity was attributed to tidal interaction (Goldreich & Soter 1966; Mazeh 2008; Pont et al. 2011).

Many of the known TEPs are giant planets, larger in size than the theoretical predictions (Mandushev et al. 2007; Hebb et al. 2009; Latham et al. 2010; Enoch et al. 2011; Smalley et al. 2012). It was found that the degree of planet inflation depends on their temperatures and orbital axes (Guillot 2005; Enoch et al. 2012).

In this paper we present the results from study of two inflated exoplanets.

2. The targets

The HATNet survey identified around 2000 exoplanet candidates and 67 were confirmed by spectroscopy. We chose to observe two of them, HAT-P-40b and HAT-P-51b, which have published data only from their discovery papers (Hartman et al. 2012, 2015).

The main goals of our observations were to improve the parameters and ephemerides of the two targets.

2.1. HAT-P-40b

HAT-P-40b was identified as a candidate TEP based on observations made with the HATNet wide-field photometric instruments during 2004–2005 (Hartman et al. 2012). It transits the F9V (V = 11.7 mag) star GSC 3607–01028 with an orbital period of ∼4.572 d. The follow-up photometric observations in Sloan i filter were obtained with the FLWO, FTN and BOS telescopes in 2010–2011. The spectroscopic observations revealed a radial velocity curve (RVC) with semi-amplitude of 58.1 m s−1.

Hartman et al. (2012) fitted the transit curves and RVC to derive period and initial epoch, orbital inclination, and relative stellar and relative planet radii (in units of planet orbital axis). The authors apply a Mandel & Agol (2002) model with fixed quadratic limb-darkening coefficients. The physical parameters of the star HAT-P-40, mass, radius, luminosity and distance, have been determined from a theoretical stellar model on the basis of its effective temperature and metallicity. Then the mass ratio and in turn the planet orbital axis are determined from the mass function. After that, the planet mass and radius are calculated.

Table 1 (column 2) reveals the obtained parameters by Hartman et al. (2012): period P (d); epoch T0 of transit center (HJD 2450000+); orbital semi-major axis a (AU); orbital inclination i (deg); stellar temperature Tst (K); stellar mass Mst (M); stellar radius Rst (R); distance d (pc); limb-darkening law (LDL); limb-darkening coefficients (LDC); planet mass Mp (MJ ); planet radius Rp (RJ ); planet density ρp (g cm−3); planet temperature Tp (K); and Safronov number Θ.

Table 1. Parameters of HAT-P-40

ParameterHartman et al. (2012)Our solution
P 4.457243(10)4.457223
T0 5813.17584(54)5813.17584*
a 0.0608(11)0.0608*
i 88.3(9)86.6(1)
Tst 6080(100)6080*
Mst 1.512(60)1.512*
Rst 2.206(61)2.248(4)
d (pc)501(16)464.5(4)*
LDLquadrquadr
LDC0.219*, 0.365* 0.31, 0.363
Mp 0.615(38)0.615*
Rp 1.730(62)1.6761(7)
ρp 0.15(1)0.165(2)
Tp 1770(33)* 1000–2600
Θ0.0286(2)0.0295(2)

*means fixed parameter.

Hartman et al. (2012) concluded that HAT-P-40b falls in the population of large radius (R > 1.5 RJ ), sub-Jupiter-mass planets, which are referred to as highly inflated planets.

2.2. HAT-P-51b

HAT-P-51b was identified as a candidate TEP based on HATNet observations during 2008–2011. It transits the G0V (V = 13.4 mag) star GSC 2296–00637 with an orbital period of ∼4.218 d. The follow-up photometric observations in Sloan i filter were obtained with the FLWO telescope in 2012. The spectroscopic observations revealed RVC with semi-amplitude of 40 m s−1.

Following the same procedure as for HAT-P-40, Hartman et al. (2015) obtained the parameter values given in Table 2 (second column). As a result, they classified HAT-P-51b as a hot Saturn-like planet that falls near the upper envelope of the distribution of points in the mass-radius diagram for transiting planets with 0.1 MJ < M < 10 MJ .

Table 2. Parameters of HAT-P-51

ParameterHartman et al. (2015)Our Solution
P 4.2180278(59)4.218022(1)
T0 6194.12204(40)6194.12204*
a 0.05069(49)0.05069*
i 88.48(0.57)87.4(1)
Tst 5450(50)5450*
Mst 0.976(28)0.976*
Rst 1.041(32)1.090(1)
d (pc)470(16)440(4)*
LDLquadrquadr
LDC0.335*, 0.2989*0.24, 0.28
Mp 0.309(18)0.309*
Rp 1.293(54)1.414(1)
ρp 0.178(24)0.138(1)
Tp 1192(21)*200–1600
Θ0.0247(18)0.0227(5)

3. Observations

Our photometric observations in R filter were carried out with the Rozhen 2-m telescope with mean photometric precision of 0.002 mag. HAT-P-40 was observed on 2019 Oct 10 with 15 s exposures. HAT-P-51 was targeted on 2018 Nov 2 with 35 s exposures.

The standard procedure was utilized for calibration of the images (de-biasing and flat-fielding) by the software MaxImDL. We tested several sets of reference stars (similar in brightness to or brighter than the targets) in order to obtain the most precise photometry. The R magnitudes (Table 3) of the chosen reference stars were obtained using their Sloan g and i magnitudes (APASS DR9, Henden et al. 2015) and the transformations of Jester et al. (2005). The photometry was performed by MaxImDL and the obtained data were cleaned of trends (Figs. 12). The shape of the transit curve of HATP-40b seems slightly serrated. This is due to the atmospheric fluctuations during that night (notable for the small transit depth).

Fig. 1

Fig. 1  Top: the Rozhen transit of HAT-P-40b and the synthetic curve corresponding to the best solution; Bottom: the residuals of the fit.

Standard image
Fig. 2

Fig. 2 The same as Fig. 1 for HAT-P-51b.

Standard image

Table 3. Coordinates and R Magnitudes of the Comparison Stars

TargetComparison starsRA (2000)DEC (2000) R
HAT-P-402MASS J01273318+385537501 27 33.2+38 55 37.614.34
 2MASS J01275452+385412701 27 54.5+38 54 12.713.98
HAT-P-51TYC 2296–996–101 24 10.6+32 47 23.912.34
 SDSS J012418.72+324925.201 24 18.8+32 49 25.415.29

4. Model of the observed transits

Our photometric data were modeled via the code TAC maker 1.1.1 (Kjurkchieva et al. 2013, 2014). It does not use any simplifications of the configuration (dark planet, linear trajectory, etc.) and may fit data by linear, quadratic, square-root and logarithmic stellar LDL.

We fixed the values of stellar temperature Tst, orbital radius a and initial epoch T0 to the values from the previous solutions (Tables 12). The fitted parameters were the period P, stellar radius Rst, planet radius Rp , orbital inclination i, planet temperature Tp and LDC. We varied these parameters around their values from the previous solutions (Tables 12) to search for the minimum of χ2. It should be noted that we fitted the planet temperature Tp freely, instead of fixing it to Teq values from the previous solutions.

We searched for transit solutions for different LDLs: linear, quadratic, square-root and logarithmic.

The results of our best transit solutions are expressed in the last columns of Tables 12. The synthetic curves are displayed in Figures 12 as continuous lines.

Considering the known values of masses Mp and Mst (Tables 12), we calculated the planet densities ρp and Safronov numbers Θ.

5. Analysis of the results

5.1. HAT-P-40

We think our transit solution of HAT-P-40b is more accurate than that of Hartman et al. (2012) because it is based on more precise data as well as it is a result of reproducing a complete transit curve as opposed to the partial transits fitted by the discoverers. The new observations led to improving the orbital period of HAT-P-40 (Table 1).

The derived stellar radius from the transit solution is bigger by 2% than that determined by Hartman et al. (2012) from stellar models (Table 1).

The obtained planet radius of 1.676 RJ is bigger than the value of 1.62 RJ calculated by the empirical relation of Enoch et al. (2012) for a Jupiter-mass planet. Thus, the solution of the newly-observed transit confirms the conclusion of Hartman et al. (2012) that HAT-P-40b belongs to the family of highly inflated Jupiters with 0.4 MJ < M < 1.5 MJ and R > 1.5 RJ .

The varying of Tp in the range 1000–2600 K led to changing of χ2 by only 0.0007%. This means that the quality of our data does not allow determination of the planet temperature within this range.

The χ2 values of the obtained best solutions for several LDLs differed by around 2%. The absolute minimum of χ2 corresponded to the quadratic LDL (Table 1).

We found essential sensibility of the solution quality in terms of the LDC. It confirms the conclusion of Espinoza & Jordán (2015) that the fitting of the LDC in the transit-solution procedure leads to better results than fixing them to theoretical values. In the case of HAT-P-40, the fitted LDC for quadratic LDL u1 = 0.310 and u2 = 0.363 (Table 1) turned out close to the theoretical values u1 = 0.329 and u2 = 0.294 for filter R (Claret et al. 1995).

The precise astrometric Gaia distance of HAT-P-40 of 464.5 pc (Bailer-Jones et al. 2018) is smaller by around 7% than the calculated value by Hartman et al. (2012) based on the star luminosity obtained from the stellar models (Table 1). This means that probably the luminosity of HAT-P-40 estimated by the stellar models is bigger than the true value by the corresponding amount.

5.2. HAT-P-51

Our transit solution should be more accurate than that of Hartman et al. (2015) because it is based on more precise photometric data. The new observations led to improving the orbital period of HAT-P-51 (Table 1).

The stellar radius derived from our transit solution is bigger by 5% than that determined by Hartman et al. (2015) from stellar models (Table 2).

The planet radius of 1.414 RJ obtained from our transit solution is considerably bigger than the value of 0.923 RJ calculated by the empirical relation of Enoch et al. (2012) for a Saturn-mass planet.

The varying of the planet temperature Tp in the range 200–1600 K led to changing of χ2 by only 0.01%. This means that the quality of our data does not allow precise determination of the planet temperature within this range.

The χ2 values of the obtained best solutions for different LDLs differed by around 2%. The absolute minimum of χ2 corresponded to the quadratic LDL of HAT-P-51. The fitted LDC for a quadratic LDL u1 = 0.24 and u2 = 0.28 (Table 2) turned out a little different from the theoretical values u1 = 0.423 and u2 = 0.242 for filter R (Claret et al. 1995).

The Gaia distance of HAT-P-51 of 440 pc (Bailer-Jones et al. 2018) is smaller by around 6% than the calculated value by Hartman et al. (2015) based on the star luminosity obtained from the stellar models (Table 2). This means that the luminosity of HAT-P-51 estimated by the stellar models is bigger than the true value.

6. Conclusions

We presented new high-quality transit observations of the exoplanets HAT-P-40b and HAT-P-51b. Their modeling led to the following results.

(1) The orbital periods of the two targets were improved.

(2) The stellar radii derived from our transit solutions are bigger than those calculated from the stellar models by Hartman et al. (2012, 2015).

(3) The planet radii of HAT-P-40b and HAT-P-51b obtained from our transit solutions are bigger than the values calculated by the empirical relations of Enoch et al. (2012) for Jupiter-mass and Saturn-mass planets respectively. This confirms the highly-inflated nature of the two targets, especially for HAT-P-51b.

(4) The best transit solutions correspond to a quadratic LDL. We found that the fitted LDC of HAT-P-40 are close to the theoretical ones whereas those of HAT-P-51 are a little different. This means that the models of stellar atmospheres applied for estimation of LDC give quite reliable results taking into account the high sensitivity of the transit solution quality for the LDC.

(5) The Gaia distances of the two targets are smaller by 6–7% than the calculated values from the stellar models. This implies that the true stellar luminosities (and radii) of HAT-P-40 and HAT-P-51 are smaller than those determined from the stellar models. Hence, the stellar models and corresponding global stellar parameters would be improved on the basis of the Gaia distances. Moreover, this would allow more reliable determination of the parameters of the known exoplanets.

Acknowledgements

The research was supported partly by project DN08/20 of Scientific Foundation of the Bulgarian Ministry of Education and Science, project D01-277/16.12.2019 of the Bulgarian Ministry of Education and Science and project RD RD-08-116/2020 of Shumen University. The authors are very grateful to the anonymous referee for the valuable comments.

This work has made use of data from the European Space Agency (ESA) mission Gaia, processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research also has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and NASA Astrophysics Data System Abstract Service.

Please wait… references are loading.
10.1088/1674-4527/21/2/44