The observable impact of runaway OB stars on protoplanetary discs

By: Gavin A. L. Coleman, Jinyoung Serena Kim, Thomas J. Haworth, Peter A. Hartman, Taylor C. Kalish

UV radiation from OB stars can drive ``external'' photoevaporative winds from discs in clusters, that have been shown to be important for disc evolution and planet formation. However, cluster dynamics can complicate the interpretation of this process. A significant fraction of OB stars are runaways, propagating at high velocity which might dominate over the wider cluster dynamics in setting the time variation of the UV field in part of the cl... more
UV radiation from OB stars can drive ``external'' photoevaporative winds from discs in clusters, that have been shown to be important for disc evolution and planet formation. However, cluster dynamics can complicate the interpretation of this process. A significant fraction of OB stars are runaways, propagating at high velocity which might dominate over the wider cluster dynamics in setting the time variation of the UV field in part of the cluster. We explore the impact of a runaway OB star on discs and the observational impact that may have. We find that discs exposed to even short periods of strong irradiation are significantly truncated, and only rebound slightly following the ``flyby'' of the UV source. This is predicted to leave an observable imprint on a disc population, with those downstream of the OB star vector being more massive and extended than those upstream. Because external photoevaporation acts quickly, this imprint is less susceptible to being washed out by cluster dynamics for faster runaway OB stars. The Gaia proper motion vector of the B star 42 Ori in NGC 1977 is transverse to the low mass stellar population and so may make a good region to search for this signature in resolved disc observations. less
Conditions for accretion favoring an unmelted Callisto and a
  differentiated Ganymede

By: Yannis Bennacer, Olivier Mousis, Marc Monnereau, Vincent Hue, Antoine Schneeberger

Analysis of Callisto's moments of inertia, derived from Galileo's gravity data, suggests that its structure is not fully differentiated. This possibly undifferentiated state contrasts sharply with the globally molten state inferred in its counterpart, Ganymede, and poses unique challenges to theories of the formation and evolution of the Galilean moons. During their formation, both moons experienced multiple heating mechanisms, including tida... more
Analysis of Callisto's moments of inertia, derived from Galileo's gravity data, suggests that its structure is not fully differentiated. This possibly undifferentiated state contrasts sharply with the globally molten state inferred in its counterpart, Ganymede, and poses unique challenges to theories of the formation and evolution of the Galilean moons. During their formation, both moons experienced multiple heating mechanisms, including tidal heating, radiogenic heating from short-lived radionuclides, accretional heating from impacts, and heat from the surrounding circumplanetary disk. Our study investigates the optimal conditions required to account for Callisto's partially differentiated state in contrast to Ganymede's complete differentiation. We investigate crucial accretion parameters, such as the timing of accretion onset, the duration of accretion, and the impactor size distribution. We find that the observed dichotomy between Ganymede and Callisto can be attributed to similar formation conditions, assuming an identical impactor size distribution and composition in the Jovian circumplanetary disk. The key differences in the formation of Ganymede and Callisto are the disk temperature at their respective formation locations and their final radii. Our results indicate that both moons accreted gradually over more than 2 Myr, concluding at least 5.5 Myr after the formation of calcium-aluminum-rich inclusions in the protosolar nebula. Our model demonstrates that Callisto can remain undifferentiated despite accreting a substantial influx of kilometer-sized impactors, potentially contributing up to 30% of the total mass inflow, while still allowing for the complete differentiation of Ganymede. less
KMT-2022-BLG-1818Lb,c: A Cold Super-Jupiter with a Saturn Sibling

By: Hongyu Li, Jiyuan Zhang, Cheongho Han, Weicheng Zang, Youn Kil Jung, Andrzej Udalski, Takahiro Sumi, Hongjing Yang, Renkun Kuang, Shude Mao, Michael D. Albrow, Sun-Ju Chung, Andrew Gould, Kyu-Ha Hwang, Yoon-Hyun Ryu, In-Gu Shin, Yossi Shvartzvald, Jennifer C. Yee, Sang-Mok Cha, Dong-Jin Kim, Seung-Lee Kim, Chung-Uk Lee, Dong-Joo Lee, Yongseok Lee, Byeong-Gon Park, Richard W. Pogge, Berto Monard, Yunyi Tang, Subo Dong, Zhuokai Liu, Grant Christie, Jennie McCormick, Tim Natusch, Qiyue Qian, Dan Maoz, Wei Zhu, Przemek Mróz, Michał K. Szymański, Jan Skowron, Radoslaw Poleski, Igor Soszyński, Paweł Pietrukowicz, Szymon Kozłowski, Krzysztof A. Rybicki, Patryk Iwanek, Krzysztof Ulaczyk, Marcin Wrona, Mariusz Gromadzki, Mateusz J. Mróz, Michał Jaroszyński, Fumio Abe, Ken Bando, David P. Bennett, Aparna Bhattacharya, Ian A. Bond, Akihiko Fukui, Ryusei Hamada, Shunya Hamada, Naoto Hamasak, Yuki Hirao, Stela Ishitani Silva, Naoki Koshimoto, Yutaka Matsubara, Shota Miyazaki, Yasushi Muraki, Tutumi Nagai, Kansuke Nunota, Greg Olmschenk, Clément Ranc, Nicholas J. Rattenbury, Yuki Satoh, Daisuke Suzuki, Sean Terry, Paul J. Tristram, Aikaterini Vandorou, Hibiki Yama

We present the discovery and analysis of the sixth microlensing two-planet system, KMT-2022-BLG-1818Lb,c, detected by a follow-up program targeting high-magnification events. Both planets are subject to the well-known ''Close/Wide'' degeneracy, although for the first planet, which has a super-Jovian mass ratio of $q_2 \simeq 5\times 10^{-3}$ in both solutions, the Close topology, with a normalized separation of $s\simeq 0.70$, is clearly pref... more
We present the discovery and analysis of the sixth microlensing two-planet system, KMT-2022-BLG-1818Lb,c, detected by a follow-up program targeting high-magnification events. Both planets are subject to the well-known ''Close/Wide'' degeneracy, although for the first planet, which has a super-Jovian mass ratio of $q_2 \simeq 5\times 10^{-3}$ in both solutions, the Close topology, with a normalized separation of $s\simeq 0.70$, is clearly preferred by $\Delta\chi^2=26$. However, contrary to all previous two-planet microlensing systems, the mass ratio for the second planet, $q_3$, is substantially (factor of $\sim 10$) different for the Close and Wide topologies of the first planet. While this degeneracy is resolved in the present case due to high-cadence follow-up observations, the appearance of this new degeneracy indicates the need for caution in the analysis of future two-planet systems. A Bayesian analysis suggests that the host is likely a K-dwarf star in the Galactic disk. The first planet is probably a super-Jupiter on a Jupiter-like orbit, while the second planet is a Saturn-class planet on either a Mercury-like or Saturn-like orbit. less
No Planet around the K Giant Star 42 Draconis

By: Artie P. Hatzes, Volker Perdelwitz, Marie Karjalainen, Jana Köhler, Michael Hartmann

Published radial velocity (RV) measurements of the K giant star 42 Dra reveal variations consistent with a 3.9 M_Jup mass companion in a 479-d orbit. This exoplanet can be confirmed if these variations are long-lived and coherent. Continued monitoring may also reveal other companions. We have acquired additional RV measurements of 42 Dra spanning fifteen years. Periodogram analyses were used to investigate the stability of the planet RV signa... more
Published radial velocity (RV) measurements of the K giant star 42 Dra reveal variations consistent with a 3.9 M_Jup mass companion in a 479-d orbit. This exoplanet can be confirmed if these variations are long-lived and coherent. Continued monitoring may also reveal other companions. We have acquired additional RV measurements of 42 Dra spanning fifteen years. Periodogram analyses were used to investigate the stability of the planet RV signal. We also investigated variations in the spectral line shapes using the bisector velocity span as well as infrared photometry from the COBE mission. The new RV measurements do not follow the published planet orbit. An orbital solution using the 2004 - 2011 data yields a period and eccentricity consistent with the published values, but the RV amplitude has decreased by a factor of four from the earlier measurements. Including some additional RV measurements taken between 2014 and 2018 reveal the presence of a second period at 530 d. The beating of this period with the one at 479-d may account for the observed amplitude variations. The planet hypothesis is conclusively ruled out by COBE/DIRBE 1.25 micron photometry that shows variations with the planet orbital period as well as an additional 170 d period. The amplitude variations in the RV as well the COBE/DIRBE photometry firmly establish that there is no giant planet around 42 Dra. The presence of multi-periodic variations suggests that these may be stellar oscillations, most likely oscillatory convection modes. These oscillations may account for some of the long period RV variations attributed to planets around K giant stars. This may skew the statistics of planet occurrence around intermediate mass stars. Long-term monitoring with excellent sampling is required to exclude amplitude variations in the long-periods found in radial velocity of K giant stars. less
Dust enrichment and growth in their earliest stages of protoplanetary
  disk formation

By: E. I. Vorobyov, V. G. Elbakyan, A. Skliarevskii, V. Akimkin, I. Kulikov

Aims. Dust enrichment and growth during the initial stages of protoplanetary disk formation were numerically investigated. A particular objective was to determine the effects of various growth barriers, which were mimicked by setting a series of upper permissible limits on maximum dust sizes. Methods. We used the ngFEOSAD code to simulate the three-dimensional dynamics of gas and dust in the polytropic approximation starting from the gravitat... more
Aims. Dust enrichment and growth during the initial stages of protoplanetary disk formation were numerically investigated. A particular objective was to determine the effects of various growth barriers, which were mimicked by setting a series of upper permissible limits on maximum dust sizes. Methods. We used the ngFEOSAD code to simulate the three-dimensional dynamics of gas and dust in the polytropic approximation starting from the gravitational collapse of a slowly rotating Bonnor-Ebert sphere to $\approx 12$ kyr after the first hydrostatic core and disk formation. Results. We found that dust growth starts in the contracting cloud in the evolution stage that precedes disk formation and the disk begins to form in an environment that is already enriched in grown dust. The efficiency of dust growth in the disk is limited by dust growth barriers. For dust grains with maximum size < 100 $\mu$m these are likely electrostatic or bouncing barriers, and for larger grains the fragmentation and drift barriers play the major role. The disk midplane becomes quickly enriched with dust, while the vertically integrated distribution of dust shows notable local variations around the canonical 1:100 dust-to-gas mass ratio. These positive and negative deviations are likely caused by local hydrodynamic flows, since the globally integrated dust-to-gas ratio deviates negligibly from the initial 1:100 value. We note that care should be taken when considering models with a fixed dust size, as it may attain a profound negative radial gradient already in the very early stages of disk formation. Models with a constant Stokes number may be preferable in this context. Conclusions. The early dust enrichment and growth may facilitate planet formation as suggested by observations of protoplanetary disk substructures. less
The Influence of Passing Field Stars on the Solar System's Dynamical
  Future

By: Nathan A. Kaib, Sean N. Raymond

The long-term dynamical future of the Sun's planets has been simulated and statistically analyzed in great detail, but most prior work considers the solar system as completely isolated, neglecting the potential influence of field star passages. To understand the dynamical significance of field star encounters, we simulate several thousand realizations of the modern solar system in the presence of passing field stars for 5 Gyrs. We find that t... more
The long-term dynamical future of the Sun's planets has been simulated and statistically analyzed in great detail, but most prior work considers the solar system as completely isolated, neglecting the potential influence of field star passages. To understand the dynamical significance of field star encounters, we simulate several thousand realizations of the modern solar system in the presence of passing field stars for 5 Gyrs. We find that the impulse gradient of the strongest stellar encounter largely determines the net dynamical effect of field stars. Because the expected strength of such an encounter is uncertain by multiple orders of magnitude, the possible significance of field stars can be large. Our simulations indicate that isolated models of the solar system can underestimate the degree of our giant planets' future secular orbital changes by over an order of magnitude. In addition, our planets and Pluto are significantly less stable than previously thought. Field stars transform Pluto from a completely stable object over 5 Gyrs to one with a ~5% instability probability. Furthermore, field stars increase the odds of Mercury's instability by ~50-80%. We also find a ~0.3% chance that Mars will be lost through collision or ejection and a ~0.2% probability that Earth will be involved in a planetary collision or ejected. Compared to previously studied instabilities in isolated solar systems models, those induced by field stars are much more likely to involve the loss of multiple planets. In addition, they typically happen sooner in our solar system's future, making field star passages the most likely cause of instability for the next 4-4.5 Gyrs. less
4 SciCasts by .
Double Hot Jupiters Through ZLK Migration

By: Yurou Liu, Tiger Lu, Malena Rice

The von Zeipel-Lidov-Kozai (ZLK) mechanism with tidal friction has been demonstrated as a promising avenue to generate hot Jupiters in stellar binary systems. Previous population studies of hot Jupiter formation have largely examined this mechanism in systems comprised of three bodies: two stars and one planet. However, because stars in a binary system form in similar environments with comparable metallicities, the formation of a single hot J... more
The von Zeipel-Lidov-Kozai (ZLK) mechanism with tidal friction has been demonstrated as a promising avenue to generate hot Jupiters in stellar binary systems. Previous population studies of hot Jupiter formation have largely examined this mechanism in systems comprised of three bodies: two stars and one planet. However, because stars in a binary system form in similar environments with comparable metallicities, the formation of a single hot Jupiter in such a system may imply that the conditions are more likely met for the companion star, as well. We investigate the ZLK mechanism with tidal friction as a potential mechanism to produce double hot Jupiter systems in stellar binaries. Using N-body simulations, we characterize the evolution of two cold Jupiters, each orbiting one star in a binary system, undergoing mirrored ZLK migration. We then examine the robustness of this mechanism to asymmetries in stellar masses, planet masses, and planet orbital inclinations relative to the binary plane. We predict that, under the assumptions that (1) most hot Jupiters in binary star systems form through ZLK migration of primordially formed cold Jupiters and (2) if one star in a binary system forms a cold Jupiter, the second does as well, a comprehensive search could identify double hot Jupiters in up to ~9% of the close- to moderate- separation $a<2000$ AU) binary systems that already host a known hot Jupiter. We also argue that a blind search for ZLK-migrated double hot Jupiters should prioritize twin stellar binaries with pericenter approaches of a few hundred AU. less
Statistical geochemical constraints on present-day water outgassing as a
  source of secondary atmospheres on the TRAPPIST-1 exoplanets

By: Trent B. Thomas, Victoria S. Meadows, Joshua Krissansen-Totton, Megan T. Gialluca, Nicholas F. Wogan, David C. Catling

The TRAPPIST-1 planetary system is observationally favorable for studying if planets orbiting M stars can retain atmospheres and host habitable conditions. Recent JWST secondary eclipse observations of TRAPPIST-1 c rule out a thick \ch{CO2} atmosphere but do not rule out atmospheric water vapor or its photochemical product, oxygen. Given the high expected escape rate, maintenance of atmospheric water vapor would require a present-day water so... more
The TRAPPIST-1 planetary system is observationally favorable for studying if planets orbiting M stars can retain atmospheres and host habitable conditions. Recent JWST secondary eclipse observations of TRAPPIST-1 c rule out a thick \ch{CO2} atmosphere but do not rule out atmospheric water vapor or its photochemical product, oxygen. Given the high expected escape rate, maintenance of atmospheric water vapor would require a present-day water source, such as volcanic outgassing. Here, we simulate water outgassing on the TRAPPIST-1 planets over a broad phase space based on solar system terrestrial bodies. We then apply two filters based on observation and geochemistry that narrow this phase space and constrain the plausible outgassing scenarios. For all seven TRAPPIST-1 planets, we find that the water outgassing rate is most likely $\sim$0.03x Earth's but has upper limits of $\sim$8x Earth's. The allowed range also implies low, Mars-like magma emplacement rates and relatively dry, Earth-like mantles, although mantle water mass fractions up to 1 wt\% are possible. We also present scenarios with magma emplacement rates similar to Mars, Earth, and Io, resulting in different preferred mantle water content and outgassing rates. We find that water outgassing rates are potentially high enough to balance water escape rates, providing a theoretical pathway for the TRAPPIST-1 planets to maintain surface water or water-vapor-containing atmospheres over long timescales. The bounds on outgassing rates and interior properties can be used in atmospheric chemistry and escape models to contextualize future observations of the TRAPPIST-1 planets, and may be applicable to other terrestrial exoplanets. less
4 SciCasts by .
A potential exomoon from the predicted planet obliquity of $β$
  Pictoris b

By: Michael Poon, Hanno Rein, Dang Pham

Planet obliquity is the alignment or misalignment of a planet spin axis relative to its orbit normal. In a multiplanet system, this obliquity is a valuable signature of planet formation and evolutionary history. The young $\beta$ Pictoris system hosts two coplanar super-Jupiters and upcoming JWST observations of this system will constrain the obliquity of the outer planet, $\beta$ Pictoris b. This will be the first planet obliquity measurem... more
Planet obliquity is the alignment or misalignment of a planet spin axis relative to its orbit normal. In a multiplanet system, this obliquity is a valuable signature of planet formation and evolutionary history. The young $\beta$ Pictoris system hosts two coplanar super-Jupiters and upcoming JWST observations of this system will constrain the obliquity of the outer planet, $\beta$ Pictoris b. This will be the first planet obliquity measurement in an extrasolar, multiplanet system. First, we show that this new planet obliquity is likely misaligned by using a wide range of simulated observations in combination with published measurements of the system. Motivated by current explanations for the tilted planet obliquities in the Solar System, we consider collisions and secular spin-orbit resonances. While collisions are unlikely to occur, secular spin-orbit resonance modified by the presence of an exomoon around the outer planet can excite a large obliquity. The largest induced obliquities ($\sim 60^\circ$) occur for moons with at least a Neptune-mass and a semimajor axis of $0.03-0.05~\mathrm{au}$ ($40-70$ planet radii). For certain orbital alignments, such a moon may observably transit the planet (transit depth of $3-7\%$, orbital period of $3-7$ weeks). Thus, a nonzero obliquity detection of $\beta$ Pictoris b implies that it may host a large exomoon. Although we focus on the $\beta$ Pictoris system, the idea that the presence of exomoons can excite high obliquities is very general and applicable to other exoplanetary systems. less
Interior and Gravity Field Models for Uranus Suggest Mixed-composition
  Interior: Implications for the Uranus Orbiter and Probe

By: Zifan Lin, Sara Seager, Benjamin P. Weiss

The interior composition and structure of Uranus are ambiguous. It is unclear whether Uranus is composed of fully differentiated layers dominated by an icy mantle or has smooth compositional gradients. The Uranus Orbiter and Probe (UOP), the next NASA Flagship mission prioritized by the Planetary Science and Astrobiology Survey 2023-2032, will constrain the planet's interior by measuring its gravity and magnetic fields. To characterize the ... more
The interior composition and structure of Uranus are ambiguous. It is unclear whether Uranus is composed of fully differentiated layers dominated by an icy mantle or has smooth compositional gradients. The Uranus Orbiter and Probe (UOP), the next NASA Flagship mission prioritized by the Planetary Science and Astrobiology Survey 2023-2032, will constrain the planet's interior by measuring its gravity and magnetic fields. To characterize the Uranian interior, here we present CORGI, a newly developed planetary interior and gravity model. We confirm that high degrees of mixing are required for Uranus interior models to be consistent with the $J_2$ and $J_4$ gravity harmonics measured by Voyager 2. Empirical models, which have smooth density profiles that require extensive mixing, can reproduce the Voyager 2 measurements. Distinct-layer models with mantles composed of H$_2$O-H/He or H$_2$O-CH$_4$-NH$_3$ mixtures are consistent with the Voyager 2 measurements if the heavy element mass fraction, $Z$, in the mantle $\lesssim85\%$, or if atmospheric $Z$ $\gtrsim25\%$. Our gravity harmonics model shows that UOP $J_2$ and $J_4$ measurements can distinguish between high ($Z\geq25\%$) and low ($Z=12.5\%$) atmospheric metallicity scenarios. The UOP can robustly constrain $J_6$ and potentially $J_8$ given polar orbits within rings. An ice-rich composition can naturally explain the source of Uranus' magnetic field. However, because the physical properties of rock-ice mixtures are poorly known, magnetic field generation by a rock-rich composition cannot be ruled out. Future experiments and simulations on realistic planetary building materials will be essential for refining Uranus interior models. less
No dilute core produced in simulations of giant impacts onto Jupiter

By: T. D. Sandnes, V. R. Eke, J. A. Kegerreis, R. J. Massey, L. F. A. Teodoro

A giant impact has been proposed as a possible formation mechanism for Jupiter's dilute core - the planet's inferred internal structure in which the transition between its core of heavy elements and its predominantly hydrogen-helium envelope is gradual rather than a discrete interface. A past simulation suggested that a head-on impact of a 10 $M_\oplus$ planet into an almost fully formed, differentiated Jupiter could lead to a post-impact p... more
A giant impact has been proposed as a possible formation mechanism for Jupiter's dilute core - the planet's inferred internal structure in which the transition between its core of heavy elements and its predominantly hydrogen-helium envelope is gradual rather than a discrete interface. A past simulation suggested that a head-on impact of a 10 $M_\oplus$ planet into an almost fully formed, differentiated Jupiter could lead to a post-impact planet with a smooth compositional gradient and a central heavy-element fraction as low as $Z\approx0.5$. Here, we present simulations of giant impacts onto Jupiter using improved numerical methods to reassess the feasibility of this scenario. We use the REMIX smoothed particle hydrodynamics (SPH) formulation, which has been newly developed to improve the treatment of mixing in SPH simulations, in particular between dissimilar materials. We perform a suite of giant impact simulations to probe the effects of impact speed, impact angle, pre-impact planet structure, and material equations of state on the evolution of heavy elements during a giant impact onto Jupiter. In all of our simulations, heavy elements re-settle over short timescales to form a differentiated core, even in cases where the core is initially disrupted into a transiently mixed state. A dilute core is not produced in any of our simulations. Our results, combined with recent observations that indicate that Saturn also has a dilute core, suggest that such structures are produced as part of the extended formation and evolution of giant planets, rather than through extreme, low-likelihood giant impacts. less
Metallicity Dependence of Giant Planets around M Dwarfs

By: Tianjun Gan, Christopher A. Theissen, Sharon X. Wang, Adam J. Burgasser, Shude Mao

We investigate the stellar metallicity ([Fe/H] and [M/H]) dependence of giant planets around M dwarfs by comparing the metallicity distribution of 746 field M dwarfs without known giant planets with a sample of 22 M dwarfs hosting confirmed giant planets. All metallicity measurements are homogeneously obtained through the same methodology based on the near-infrared spectra collected with a single instrument SpeX mounted on the NASA Infrared... more
We investigate the stellar metallicity ([Fe/H] and [M/H]) dependence of giant planets around M dwarfs by comparing the metallicity distribution of 746 field M dwarfs without known giant planets with a sample of 22 M dwarfs hosting confirmed giant planets. All metallicity measurements are homogeneously obtained through the same methodology based on the near-infrared spectra collected with a single instrument SpeX mounted on the NASA Infrared Telescope Facility. We find that 1) giant planets favor metal-rich M dwarfs at a 4-5$\sigma$ confidence level, depending on the band of spectra used to derive metallicity; 2) hot ($a/R_\ast\leq 20$) and warm ($a/R_\ast> 20$) Jupiters do not show a significant difference in the metallicity distribution. Our results suggest that giant planets around M and FGK stars, which are already known to prefer metal-rich hosts, probably have a similar formation channel. In particular, hot and warm Jupiters around M dwarfs may have the same origin as they have indistinguishable metallicity distributions. With the refined stellar and planetary parameters, we examine the stellar metallicities and the masses of giant planets where we find no significant correlation. M dwarfs with multiple giant planets or with a single giant planet have similar stellar metallicities. Mid-to-late type M stars hosting gas giants do not show an apparent preference to higher metallicities compared with those early-M dwarfs with gas giants and field M dwarfs. less
TOI-4504: Exceptionally large Transit Timing Variations induced by two
  resonant warm gas giants in a three planet system

By: Michaela Vítková, Rafael Brahm, Trifon Trifonov, Petr Kabáth, Andrés Jordán, Thomas Henning, Melissa J. Hobson, Jan Eberhardt, Marcelo Tala Pinto, Felipe I. Rojas, Nestor Espinoza, Martin Schlecker, Matías I. Jones, Maximiliano Moyano, Susana Eyheramendy, Carl Ziegler, Jack J. Lissauer, Andrew Vanderburg, Karen A. Collins, Bill Wohler, David Watanabe, George R. Ricker, Roland Vanderspek, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Marek Skarka

We present a joint analysis of TTVs and Doppler data for the transiting exoplanet system TOI-4504. TOI-4504 c is a warm Jupiter-mass planet that exhibits the largest known transit timing variations (TTVs), with a peak-to-node amplitude of $\sim$ 2 days, the largest value ever observed, and a super-period of $\sim$ 930 d. TOI-4504 b and c were identified in public TESS data, while the TTVs observed in TOI-4504 c, together with radial velocit... more
We present a joint analysis of TTVs and Doppler data for the transiting exoplanet system TOI-4504. TOI-4504 c is a warm Jupiter-mass planet that exhibits the largest known transit timing variations (TTVs), with a peak-to-node amplitude of $\sim$ 2 days, the largest value ever observed, and a super-period of $\sim$ 930 d. TOI-4504 b and c were identified in public TESS data, while the TTVs observed in TOI-4504 c, together with radial velocity (RV) data collected with FEROS, allowed us to uncover a third, non-transiting planet in this system, TOI-4504 d. We were able to detect transits of TOI-4504 b in the TESS data with a period of 2.4261$\pm 0.0001$ days and derive a radius of 2.69$\pm 0.19$ R$_{\oplus}$. The RV scatter of TOI-4504 was too large to constrain the mass of TOI-4504 b, but the RV signals of TOI-4504 c \& d were sufficiently large to measure their masses. The TTV+RV dynamical model we apply confirms TOI-4504 c as a warm Jupiter planet with an osculating period of 82.54$\pm 0.02$ d, mass of 3.77$\pm 0.18$ M$_{\rm J}$ and a radius of 0.99$\pm 0.05$ R$_{\rm J}$, while the non-transiting planet TOI-4504 d, has an orbital period of 40.56$\pm 0.04$ days and mass of 1.42$_{-0.06}^{+0.07}$ M$_{\rm J}$. We present the discovery of a system with three exoplanets: a hot sub-Neptune and two warm Jupiter planets. The gas giant pair is stable and likely locked in a first-order 2:1 mean-motion resonance (MMR). The TOI-4504 system is an important addition to MMR pairs, whose increasing occurrence supports a smooth migration into a resonant configuration during the protoplanetary disk phase. less
Eta-Earth Revisited I: A Formula for Estimating the Maximum Number of
  Earth-like Habitats

By: Helmut Lammer, Manuel Scherf, Laurenz Sproß

In this hypothesis article, we discuss the basic requirements of planetary environments where aerobe organisms can grow and survive, including atmospheric limitations of millimeter-to-meter-sized biological animal life based on physical limits, and O$_2$, N$_2$, and CO$_2$ toxicity levels. By assuming that animal-like extraterrestrial organisms adhere to similar limits, we define Earth-like Habitats ($\eta_{\rm EH}$) as rocky exoplanets in ... more
In this hypothesis article, we discuss the basic requirements of planetary environments where aerobe organisms can grow and survive, including atmospheric limitations of millimeter-to-meter-sized biological animal life based on physical limits, and O$_2$, N$_2$, and CO$_2$ toxicity levels. By assuming that animal-like extraterrestrial organisms adhere to similar limits, we define Earth-like Habitats ($\eta_{\rm EH}$) as rocky exoplanets in the Habitable Zone of Complex Life that host N$_2$-O$_2$-dominated atmospheres with minor amounts of CO$_2$, at which advanced animal-like life can in principle evolve and exist. We then derive a new formula that can be used to estimate the maximum occurrence rate of such Earth-like Habitats in the Galaxy. This contains realistic probabilistic arguments that can be fine-tuned and constrained by atmospheric characterization with future space and ground-based telescopes. As an example, we briefly discuss two specific requirements feeding into our new formula that, although not quantifiable at present, will become scientifically quantifiable in the upcoming decades due to future observations of exoplanets and their atmospheres. less