Jovian-mass Planets Research Articles

Connections between Planetary Populations and Chemical Characteristics of Their Host Stars

Chemical anomalies in planet-hosting stars (PHSs) are studied in order to assess how the planetary nature and multiplicity affect the atmospheric chemical abundances of their host stars. We employ APOGEE DR17 to select thin-disk stars of the Milky Way, and crossmatch them with the Kepler Input Catalog to identify confirmed PHSs, which results in 227 PHSs with available chemical abundance ratios for six refractory elements. We also examine an ensemble of stars without planet signals, which are equivalent to the selected PHSs in terms of evolutionary stage and stellar parameters, to correct for Galactic chemical evolution effects, and derive the abundance gradient of refractory elements over the condensation temperature for the PHSs. Using the Galactic chemical evolution corrected abundances, we find that our PHSs do not show a significant difference in abundance slope from the stars without planets. However, when we examine the trends of the refractory elements of PHSs, based on the total number of their planets and their planet types, we find that the PHSs with giant planets are more depleted in refractory elements than those with rocky planets. Among the PHSs with rocky planets, the refractory depletion trends are potentially correlated with the terrestrial planets’ radii and multiplicity. In the cases of PHSs with giant planets, sub-Jovian PHSs demonstrate more depleted refractory trends than stars hosting Jovian-mass planets, raising questions on different planetary formation processes for Neptune-like and Jupiter-like planets.

OGLE-2017-BLG-0448Lb: A Low Mass–Ratio Wide-orbit Microlensing Planet?

The gravitational microlensing technique is most sensitive to planets in a Jupiter-like orbit and has detected more than 200 planets. However, only a few wide-orbit (s > 2) microlensing planets have been discovered, where s is the planet-to-host separation normalized to the angular Einstein ring radius, θ E. Here, we present the discovery and analysis of a strong candidate wide-orbit microlensing planet in the event OGLE-2017-BLG-0448. The whole light curve exhibits long-term residuals to the static binary-lens single-source model, so we investigate the residuals by adding the microlensing parallax, microlensing xallarap, an additional lens, or an additional source. For the first time, we observe a complex degeneracy between all four effects. The wide-orbit models with s ∼ 2.5 and a planet-to-host mass ratio of q ∼ 10−4 are significantly preferred, but we cannot rule out the close models with s ∼ 0.35 and q ∼ 10−3. A Bayesian analysis based on a Galactic model indicates that, despite the complicated degeneracy, the surviving wide-orbit models all contain a super-Earth-mass to Neptune-mass planet at a projected planet-host separation of ∼6 au and the surviving close-orbit models all consist of a Jovian-mass planet at ∼1 au. The host star is probably an M or K dwarf. We discuss the implications of this dimension-degeneracy disaster on microlensing light-curve analysis and its potential impact on statistical studies.

TOI-4201: An Early M Dwarf Hosting a Massive Transiting Jupiter Stretching Theories of Core Accretion* *This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory, Chile.

We confirm TOI-4201 b as a transiting Jovian-mass planet orbiting an early M dwarf discovered by the Transiting Exoplanet Survey Satellite. Using ground-based photometry and precise radial velocities from NEID and the Planet Finder Spectrograph, we measure a planet mass of 2.59−0.06+0.07 M J, making this one of the most massive planets transiting an M dwarf. The planet is ∼0.4% of the mass of its 0.63 M ⊙ host and may have a heavy-element mass comparable to the total dust mass contained in a typical class II disk. TOI-4201 b stretches our understanding of core accretion during the protoplanetary phase and the disk mass budget, necessitating giant planet formation to take place either much earlier in the disk lifetime or perhaps through alternative mechanisms like gravitational instability.

Searching for Planets Orbiting Fomalhaut with JWST/NIRCam

We report observations with the JWST/NIRCam coronagraph of the Fomalhaut (α PsA) system. This nearby A star hosts a complex debris disk system discovered by the IRAS satellite. Observations in F444W and F356W filters using the round 430R mask achieve a contrast ratio of ∼4 × 10−7 at 1″ and ∼4 × 10−8 outside of 3″. These observations reach a sensitivity limit of <1 M Jup across most of the disk region. Consistent with the hypothesis that Fomalhaut b is not a massive planet but is a dust cloud from a planetesimal collision, we do not detect it in either F356W or F444W (the latter band where a Jovian-sized planet should be bright). We have reliably detected 10 sources in and around Fomalhaut and its debris disk, all but one of which are coincident with Keck or Hubble Space Telescope sources seen in earlier coronagraphic imaging; we show them to be background objects, including the “Great Dust Cloud” identified in Mid-Infrared Instrument (MIRI) data. However, one of the objects, located at the edge of the inner dust disk seen in the MIRI images, has no obvious counterpart in imaging at earlier epochs and has a relatively red [F356W]–[F444W] > 0.7 mag (Vega) color. Whether this object is a background galaxy, brown dwarf, or a Jovian-mass planet in the Fomalhaut system will be determined by an approved Cycle 2 follow-up program. Finally, we set upper limits to any scattered light from the outer ring, placing a weak limit on the dust albedo at F356W and F444W.

Gas accretion onto Jupiter mass planets in discs with laminar accretion flows

Context. Numerous studies have shown that a gap-forming Jovian mass planet embedded in a protoplanetary disc, in which a turbulent viscosity operates, can accrete gas efficiently through the gap, and for typical parameters it doubles its mass in ~0.1 Myr. The planet also migrates inwards on a timescale that is closely related to the local viscous evolution timescale, which is also typically 0.1 Myr. These timescales are short compared to protoplanetary disc lifetimes, and raise questions about the origins of the cold gas giant exoplanets that have been discovered in abundance. It is understood that protoplanetary discs are unlikely to be globally turbulent, and instead they may launch magnetised winds such that accretion towards the star occurs in laminar accretion flows located in narrow layers near the surfaces of the disc. Aims. The aim of this study is to examine the rate at which gas accretes onto Jovian mass planets that are embedded in layered protoplanetary discs, and to compare the results with those obtained for viscous models. Methods. We use 3D hydrodynamical simulations of planets embedded in protoplanetary discs, in which a constant radial mass flux towards the star of ṁ = 10−8 M⊙ yr−1 is sustained. We consider a classical viscous α disc model, and also models in which an external torque is applied in narrow surface layers to mimic the effects of a magnetised wind. The accreting layers have a variety of depths, as parameterised by their column densities ΣA, and we consider values of ΣA in the range 0.1−10 g cm−2. Results. The viscous disc model gives results in agreement with previous studies. In accord with our recent work that examines the migration of Jovian mass planets in layered models, we find the accretion rate onto the planet in the layered models crucially depends on the ability of the planet to block the wind-induced mass flow towards the star. For ΣA = 10 g cm−2, the planet torque can block the mass flow in disc, accretion onto the planet is slow, and a mass doubling time of 10 Myr is obtained. For ΣA = 0.1 g cm−2, the flow is not blocked, accretion is fast, and the mass doubling time is 0.2 Myr. Conclusions. Our results show that although the radial mass flow through the layered disc models is always 10−8 M⊙ yr−1, adopting different values of ΣA leads to very different gas accretion rates onto embedded gas giant planets.

Evaporation before disruption: Comparing time-scales for Jovian planets in star-forming regions

ABSTRACT Simulations show that the orbits of planets are readily disrupted in dense star-forming regions; planets can also be exchanged between stars, or become free-floating and then be captured by other stars. However, dense star-forming regions also tend to be populous, containing massive stars that emit photoionizing radiation, which can evaporate the gas in protoplanetary discs. We analyse N-body simulations of star-forming regions containing Jovian-mass planets and determine the times when their orbits are altered, when they become free-floating, and when they are stolen or captured. Simultaneously, we perform calculations of the evolution of protoplanetary discs when exposed to FUV radiation fields from massive stars in the same star-forming regions. In almost half (44 per cent) of the planetary systems that are disrupted – either altered, captured, stolen or become free-floating, we find that the radius of the protoplanetary disc evolves inwards, or the gas in the disc is completely evaporated, before the planets’ orbits are disrupted. This implies that planets that are disrupted in dense, populous star-forming regions are more likely to be super-Earths or mini-Neptunes, as Jovian mass planets would not be able to form due to mass-loss from photoevaporation. Furthermore, the recent discoveries of distant Jovian mass planets around tightly packed terrestrial planets argue against their formation in populous star-forming regions, as photoevaporation would preclude gas giant planet formation at distances of more than a few au.

A New Third Planet and the Dynamical Architecture of the HD 33142 Planetary System* *Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programs 60.A-9700, 60.A-9036, 097.C-0090, 0100.C-0414, 0101.C-0232, 0102.C-0338, and MPG programs 088.C-0892, 099.A-9009, 0100.A-9006.

Based on recently-taken and archival HARPS, FEROS, and HIRES radial velocities (RVs), we present evidence for a new planet orbiting the first ascent red giant star HD 33142 (with an improved mass estimate of M1.52 ± 0.03 M ⊙), already known to host two planets. We confirm the Jovian-mass planets HD 33142b and c, with periods of P b = days and P c = days and minimum dynamical masses of = M Jup and = MJup, respectively. Furthermore, our periodogram analysis of the precise RVs shows strong evidence for a short-period Doppler signal in the residuals of a two-planet Keplerian fit, which we interpret as a third, Saturn-mass planet with = M Jup in a close-in orbit with an orbital period of P d = days. We study the dynamical behavior of the three-planet system configuration with an N-body integration scheme, finding it to be long-term stable with the planets alternating between low and moderate eccentricity episodes. We also perform N-body simulations, including stellar evolution and second-order dynamical effects such as planet–stellar tides and stellar mass loss on the way to the white dwarf phase. We find that planets HD 33142b, c, and d are likely to be engulfed near the tip of the red giant branch phase due to tidal migration. These results make the HD 33142 system an essential benchmark for planet population statistics of the multiple-planet systems found around evolved stars.

SOLES I: The Spin–Orbit Alignment of K2-140 b

Abstract Obliquity measurements for stars hosting relatively long-period giant planets with weak star-planet tidal interactions may play a key role in distinguishing between formation theories for shorter-period hot Jupiters. Few such obliquity measurements have been made to date due to the relatively small sample of known wide-orbiting, transiting Jovian-mass planets and the challenging nature of these targets, which tend to have long transit durations and orbit faint stars. We report a measurement of the Rossiter–McLaughlin effect across the transit of K2-140 b, a Jupiter-mass planet with period P = 6.57 days orbiting a V = 12.6 star. We find that K2-140 is an aligned system with projected spin–orbit angle λ = 0.5° ± 9.7°, suggesting a dynamically cool formation history. This observation builds toward a population of tidally detached giant planet spin–orbit angles that will enable a direct comparison with the distribution of close-orbiting hot-Jupiter orbital configurations, elucidating the prevalent formation mechanisms of each group.

Observational Signature of Tightly Wound Spirals Driven by Buoyancy Resonances in Protoplanetary Disks

Abstract Besides the spirals induced by the Lindblad resonances, planets can generate a family of tightly wound spirals through buoyancy resonances. The excitation of buoyancy resonances depends on the thermal relaxation timescale of the gas. By computing timescales of various processes associated with thermal relaxation, namely, radiation, diffusion, and gas–dust collision, we show that the thermal relaxation in protoplanetary disks’ surface layers (Z/R ≳ 0.1) and outer disks (R ≳ 100 au) is limited by infrequent gas–dust collisions. The use of the isothermal equation of state or rapid cooling, common in protoplanetary disk simulations, is therefore not justified. Using three-dimensional hydrodynamic simulations, we show that the collision-limited slow thermal relaxation provides favorable conditions for buoyancy resonances to develop. Buoyancy resonances produce predominantly vertical motions, whose magnitude at the 12CO emission surface is of the order of 100 m s−1 for Jovian-mass planets, sufficiently large to detect using molecular line observations with ALMA. We generate synthetic observations and describe characteristic features of buoyancy resonances in Keplerian-subtracted moment maps and velocity channel maps. Based on the morphology and magnitude of the perturbation, we propose that the tightly wound spirals observed in TW Hya could be driven by a (sub-)Jovian-mass planet at 90 au. We discuss how non-Keplerian motions driven by buoyancy resonances can be distinguished from those driven by other origins. We argue that observations of multiple lines tracing different heights, with sufficiently high spatial/spectral resolution and sensitivity to separate the emission arising from the near and far sides of the disk, will help constrain the origin of non-Keplerian motions.

KMT-2019-BLG-1339L: An M Dwarf with a Giant Planet or a Companion near the Planet/Brown Dwarf Boundary

Abstract We analyze KMT-2019-BLG-1339, a microlensing event with an obvious but incompletely resolved brief anomaly feature around the peak of the light curve. Although the origin of the anomaly is identified to be a companion to the lens with a low mass ratio q, the interpretation is subject to two different degeneracy types. The first type is the ambiguity in ρ, representing the angular source radius scaled to the angular radius of the Einstein ring, θ E, and the other is the s ↔ s −1 degeneracy. The former type, “finite-source degeneracy,” causes ambiguities in both s and q, while the latter induces an ambiguity only in s. Here, s denotes the separation (in units of θ E) in projection between the lens components. We estimate that the lens components have masses and according to the two solutions subject to the finite-source degeneracy, indicating that the lens comprises an M dwarf and a companion with a mass around the planet/brown dwarf boundary or a Jovian-mass planet. It is possible to lift the finite-source degeneracy by conducting future observations utilizing a high-resolution instrument because the relative lens–source proper motion predicted by the solutions is widely different.

Constraints on a Putative Planet Sculpting the V4046 Sagittarii Circumbinary Disk

Abstract We analyze the highest-resolution millimeter continuum and near-infrared (NIR) scattered-light images presented to date of the circumbinary disk orbiting V4046 Sgr, a ∼20-Myr-old actively accreting, close binary T Tauri star system located a mere 72.4 pc from Earth. We observed the disk with the Atacama Large Millimeter/submillimeter Array (ALMA) at 870 μm during Cycle 4, and we analyze these data in conjunction with archival NIR (H band) polarimetric images obtained with SPHERE/IRDIS on the ESO Very Large Telescope. At 0.″3 (20 au) resolution, the 870 μm image reveals a marginally resolved ring that peaks at ∼32 au and has an extension of ∼90 au. We infer a lower limit on a dust mass of ∼60.0 M ⊕ within the 870 μm ring, and confirm that the ring is well aligned with the larger -scale gaseous disk. A second, inner dust ring is also tentatively detected in the ALMA observations; its position appears coincident with the inner (∼14 au radius) ring detected in scattered light. Using synthetic 870 μm and H-band images obtained from disk–planet interaction simulations, we attempt to constrain the mass of the putative planet orbiting at 20 au. Our trials suggest that a circumbinary Jovian-mass planet may be responsible for generating the dust ring and gap structures detected within the disk. We discuss the longevity of the gas-rich disk orbiting V4046 Sgr in the context of the binary nature of the system.

Exploring the formation by core accretion and the luminosity evolution of directly imaged planets

Context. A low-mass companion to the two-solar mass star HIP 65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. Aims. We want to derive statistical constraints on the mass and initial entropy of HIP 65426 b and to explore possible formation pathways of directly imaged objects within the core-accretion paradigm, focusing on HIP 65426 b. Methods. Constraints on the planet’s mass and post-formation entropy are derived from its age and luminosity combined with cooling models. For the first time, the results of population synthesis are also used to inform the results. Then a formation model that includes N-body dynamics with several embryos per disc is used to study possible formation histories and the properties of possible additional companions. Finally, the outcomes of two- and three-planet scattering in the post-disc phase are analysed, taking tides into account for small-pericentre orbits. Results. The mass of HIP 65426 b is found to be mp = 9.9−1.8+1.1 MJ using the hot population and mp = 10.9−2.0+1.4 MJ with the cold-nominal population. We find that core formation at small separations from the star followed by outward scattering and runaway accretion at a few hundred astronomical units succeeds in reproducing the mass and separation of HIP 65426 b. Alternatively, systems having two or more giant planets close enough to be on an unstable orbit at disc dispersal are likely to end up with one planet on a wide HIP 65426 b-like orbit with a relatively high eccentricity (≳ 0.5). Conclusions. If this scattering scenario explains its formation, HIP 65426 b is predicted to have a high eccentricity and to be accompanied by one or several roughly Jovian-mass planets at smaller semi-major axes, which also could have a high eccentricity. This could be tested by further direct-imaging as well as radial-velocity observations.

Shock-generating Planetesimals Perturbed by a Giant Planet in a Gas Disk

Abstract We examined the excitations of planetesimals caused by the resonances of a giant planet in a protoplanetary gas disk. The highly excited planetesimals generate bow shocks, the mechanism of which results in chondrule formation, crystallization of silicate dust, and evaporation of icy planetesimals. The planetesimals beyond 2:1 resonance migrate owing to the gas drag and obtain the maximum eccentricity around 3:1 resonance, which is located at approximately half the planetary distance. The eccentricity depends on the parameters of the planetesimals and the Jovian planet, such as size and location, and the gas density of the disk. The maximum relative velocity of a 100 km sized planetesimal with respect to the gas disk reaches up to ∼12 km s−1 in the case of Jupiter owing to secular resonance, which occurs because of the disk’s gravity. We find that if a Jovian-mass planet is located within 10 au, the planetesimals larger than 100 km gain sufficient velocity to cause the melting of chondrule precursors and crystallization of the silicate. The maximum velocity is higher for large planetesimals and eccentric planets. Planetesimals are trapped temporarily in the resonances and continue to have high speed over ≳1 Myr after the formation of a Jovian planet. This duration fits into the timescale of chondrule formation suggested by the isotopic data. The evaporation of icy planetesimals occurs when a Jovian planet is located within 15 au. This mechanism can be a new indicator of planet formation in exosystems if some molecules ejected from icy planetesimals are detected.

Dynamical Constraints on Nontransiting Planets Orbiting TRAPPIST-1

Abstract We derive lower bounds on the orbital distance and inclination of a putative planet beyond the transiting seven planets of TRAPPIST-1, for a range of masses ranging from 0.08 M Jup to 3.5 M Jup. While the outer architecture of this system will ultimately be constrained by radial velocity measurements over time, we present dynamical constraints from the remarkably coplanar configuration of the seven transiting planets, which is sensitive to modestly inclined perturbers. We find that the observed configuration is unlikely if a Jovian-mass planet inclined by ≥3° to the transiting planet exists within 0.53 au, exceeding any constraints from transit timing variations (TTV) induced in the known planets from an undetected perturber. Our results will inform RV programs targeting TRAPPIST-1, and for near coplanar outer planets, tighter constraints are anticipated for radial velocity (RV) precisions of ≲140 m s−1. At higher inclinations, putative planets are ruled out to greater orbital distances with orbital periods up to a few years.

Outer Architecture of Kepler-11: Constraints from Coplanarity

Abstract The outer architectures of Kepler’s compact systems of multiple transiting planets remain poorly constrained, and few of these systems have lower bounds on the orbital distance of any massive outer planets. We infer a minimum orbital distance and upper limits on the inclination of a hypothetical Jovian-mass planet orbiting exterior to the six transiting planets at Kepler-11. Our constraints are derived from dynamical models together with observations provided by the Kepler mission. First, the lack of transit timing variations (TTV) in the outermost transiting planet Kepler-11 g imply that the system does not contain a Jovian-mass perturber within 2 au from the star. Second, we test under what initial conditions a Jovian-mass planet moderately inclined from the transiting planets would make their co-transiting configuration unlikely. The transiting planets are secularly coupled and exhibit small mutual inclinations over long timescales, although the outermost transiting planet, Kepler-11 g, is weakly coupled to the inner five. We rule out a Jovian-mass planet on a 3° inclination within 3.0 au, and higher inclinations out to farther orbital distances, unless an undetected planet exists orbiting in the dynamical gap between Kepler-11 f and Kepler-11 g. Our constraints depend little on whether we assume the six transiting planets of Kepler-11 were initially perfectly coplanar or whether a minimum initial mutual inclination between the transiting planets is adopted based on the measured impact parameters of the transiting planets.

THREE TEMPERATE NEPTUNES ORBITING NEARBY STARS*

ABSTRACT We present the discovery of three modestly irradiated, roughly Neptune-mass planets orbiting three nearby Solar-type stars. HD 42618 b has a minimum mass of 15.4 ± 2.4 , a semimajor axis of 0.55 au, an equilibrium temperature of 337 K, and is the first planet discovered to orbit the solar analogue host star, HD 42618. We also discover new planets orbiting the known exoplanet host stars HD 164922 and HD 143761 (ρ CrB). The new planet orbiting HD 164922 has a minimum mass of 12.9 ± 1.6 and orbits interior to the previously known Jovian mass planet orbiting at 2.1 au. HD 164922 c has a semimajor axis of 0.34 au and an equilibrium temperature of 418 K. HD 143761 c orbits with a semimajor axis of 0.44 au, has a minimum mass of 25 ± 2 , and is the warmest of the three new planets with an equilibrium temperature of 445 K. It orbits exterior to the previously known warm Jupiter in the system. A transit search using space-based CoRoT data and ground-based photometry from the Automated Photometric Telescopes (APTs) at Fairborn Observatory failed to detect any transits, but the precise, high-cadence APT photometry helped to disentangle planetary-reflex motion from stellar activity. These planets were discovered as part of an ongoing radial velocity survey of bright, nearby, chromospherically inactive stars using the Automated Planet Finder (APF) telescope at Lick Observatory. The high-cadence APF data combined with nearly two decades of radial velocity data from Keck Observatory and gives unprecedented sensitivity to both short-period low-mass, and long-period intermediate-mass planets.

On the orbital evolution of a pair of giant planets in mean motion resonance

Pairs of extrasolar giant planets in a mean motion commensurability are common with 2:1 resonance occurring most frequently. Disc-planet interaction provides a mechanism for their origin. However, the time scale on which this could operate in particular cases is unclear. We perform 2D and 3D numerical simulations of pairs of giant planets in a protoplanetary disc as they form and maintain a mean motion commensurability. We consider systems with current parameters similar to those of HD 155358, 24 Sextantis and HD 60532, and disc models of varying mass, decreasing mass corresponding to increasing age. For the lowest mass discs, systems with planets in the Jovian mass range migrate inwards maintaining a 2:1 commensurability. Systems with the inner planet currently at around 1 au from the central star could have originated at a few au and migrated inwards on a time scale comparable to protoplanetary disc lifetimes. Systems of larger mass planets such as HD 60532 attain 3:1 resonance as observed. For a given mass accretion rate, results are insensitive to the disc model for the range of viscosity prescriptions adopted, there being good agreement between 2D and 3D simulations. However, in a higher mass disc a pair of Jovian mass planets passes through 2:1 resonance before attaining a temporary phase lasting a few thousand orbits in an unstable 5:3 resonance prior to undergoing a scattering. Thus finding systems in this commensurability is unlikely.

GAS GAPS IN THE PROTOPLANETARY DISK AROUND THE YOUNG PROTOSTAR HL TAU

ABSTRACT We have analyzed the HCO+ (1–0) data of the Class I–II protostar, HL Tau, obtained from the Atacama Large Millimeter/submillimeter Array long baseline campaign. We generated the HCO+ image cube at an angular resolution of ∼0.″07 (∼10 au) and performed azimuthal averaging on the image cube to enhance the signal-to-noise ratio and measure the radial profile of the HCO+ integrated intensity. Two gaps at radii of ∼28 and ∼69 au and a central cavity are identified in the radial intensity profile. The inner HCO+ gap is coincident with the millimeter continuum gap at a radius of 32 au. The outer HCO+ gap is located at the millimeter continuum bright ring at a radius of 69 au and overlaps with the two millimeter continuum gaps at radii of 64 and 74 au. On the contrary, the presence of the central cavity is likely due to the high optical depth of the 3 mm continuum emission and not the depletion of the HCO+ gas. We derived the HCO+ column density profile from its intensity profile. From the column density profile, the FWHM widths of the inner and outer HCO+ gaps are both estimated to be ∼14 au, and their depths are estimated to be ∼2.4 and ∼5.0. These results are consistent with the expectation from the gaps opened by forming (sub-)Jovian mass planets, while placing tight constraints on the theoretical models solely incorporating the variation of dust properties and grain sizes.

3.6 AND 4.5μm PHASE CURVES OF THE HIGHLY IRRADIATED ECCENTRIC HOT JUPITER WASP-14b

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

Multiple mean motion resonances in the HR 8799 planetary system

HR 8799 is a nearby star hosting at least four ~10 Jovian mass planets in wide orbits up to ~70au, detected through the direct, high-contrast infrared imaging. Large companions and debris disks reported interior to ~10au, and exterior to ~100au indicate massive protoplanetary disc in the past. The dynamical state of the HR 8799 system is not yet fully resolved, due to limited astrometric data covering tiny orbital arcs. We construct a new, orbital model of the HR 8799 system, assuming rapid migration of the planets after their formation in wider orbits. We found that the HR 8799 planets are likely involved in double Laplace resonance, 1e:2d:4c:8b MMR. Quasi-circular planetary orbits are coplanar with the stellar equator and inclined by ~25 degrees to the sky plane. This best-fit orbital configuration matches astrometry, debris disk models, and mass estimates from cooling models. The multiple MMR is stable for the age of the star ~160Myr, for at least 1 Gyr unless significant perturbations to the N-body dynamics are present. We predict four configurations with the fifth hypothetical innermost planet HR 8799f in ~9.7au, or ~7.5au orbit, extending the MMR chain to triple Laplace resonance 1f:2e:4d:8c:16b MMR or to the 1f:3e:6d:12c:24b MMR, respectively. Our findings may establish strong boundary conditions for the system formation and its early history.