2026-07-14 | | Total: 20
Using a 1D mass and heat transport model, we numerically integrate heat flow and gas transport in a porous body exposed to a transient heating event while embedded in a protoplanetary disk. When small icy grains are heated, volatiles sublimate, enriching the disk with volatile gases. When a porous body enters this heated, volatile-rich environment, volatile gases diffuse throughout the cool, porous body and deposit ice where the partial pressure of a volatile exceeds its vapor pressure. We simulate sublimation and deposition fronts of water, carbon dioxide, and carbon monoxide. Our simulations show that an ice deposition front forms and moves deeper into the porous body as the body warms. The amount of nebular gas deposited in an initially dry body is usually extremely low; however, in an initially icy body, an ice deposition front contains locally sublimated volatiles. In this case, the front can increase the ice volume fraction (by a factor of 2) in a thin layer below the surface. We find that the propagation speed, propagation strength, and final depth of an ice deposition front primarily depend on pore size. We propose that nebular heating events can alter the subsurface morphology and physical properties of porous icy objects embedded in a protoplanetary disk.
Stellar activity is a major complication in the detection and characterization of exoplanets by both radial velocities and transits, and the upcoming Habitable Worlds Observatory (HWO) invites us to also consider its effect on direct imaging. Spectra of directly-imaged planets can vary with the activity of their host stars because the face of the star we see is not the same as the face reflected by the planet. This discrepancy could potentially result in inaccurate measurements of the planet's radius and unexpected, externally-caused variability in its contrast spectrum with the star. To assess the scientific capabilities and requirements of HWO, it is important that the magnitude of these effects be quantified. We present results of model retrievals of Earth-like exoplanets observed with an HWO-style survey, as they would appear when affected by starspots, using the ExoVista code for spectrum generation and the BARBIE code for spectroscopic retrieval. Both Solar-type stellar activity and highly active rapid rotators are considered and compared with an idealized quiescent host star. In the quiescent case, ${\rm SNR}\approx5$ is needed to detect atmospheric water vapor at 0.9 microns and ${\rm SNR}\approx13$ at 0.74 microns. We find that for Solar-type activity, the effect on retrievals will be negligible, but it could present problems for certain highly-active stars at limiting geometries. For an extreme case with a single large spot with 10\% coverage not visible to the observer, the SNR required for water detection increases to ${\rm SNR}\approx8$ at 0.9 microns and ${\rm SNR}\approx18$ at 0.74 microns. It also decreases the accuracy of the retrieved albedo, resulting in a value $\sim2/3$ of the true value. In light of these results, we estimate the impact that stellar variability and starspots may have on an HWO-style survey.
Observations of tens of secondary CO gas discs associated with cold exo-Kuiper belts together with other arguments have led Kral et al (2024) to propose that water ice could also sublimate in exo-asteroid belts, suggesting a new pathway for the delivery of water to terrestrial planets, including Earth. We aim to model such water vapour discs and to characterise their physical properties across a range of extrasolar systems with different host stars. We further investigate the implications for the accretion of this water by potential planets located in the inner regions of these systems. We adapt and extend the model of Kral et al (2024) to follow the outgassing, photodissociation, and viscous evolution of water vapour discs. We perform a suite of simulations exploring the parameter space, focusing on the stellar mass, the mass of the parent belt, and its orbital location. We additionally include an inner planet to estimate the mass of water accreted as a function of disc properties and system architecture. We find that systems hosting Sun's mass (and higher) stars produce water vapour very efficiently, sublimating nearly all of the ice initially present in the belt. In most cases, the bulk of the gas mass is generated early, when the stellar luminosity is highest. The amount of water accreted by inner planets can approach the initial ice mass of the belt, leading to planets with water inventories comparable to or exceeding those of Earth, potentially creating ocean planets. We find that water outgassing occurs early after the protoplanetary disc dissipates in systems containing exo-asteroid belts. ALMA, JWST and ELT are capable of detecting this water vapour for several tens of millions of years, even in relatively low-mass water discs. Hence, if such water gas discs are present, they should be detectable with current facilities.
Brown dwarfs are massive analogues of extrasolar giant planets. Compared to exoplanets whose observations are generally limited by the presence of their bright host star, brown dwarfs are ideal targets for studying substellar atmospheric physics, chemistry and dynamics. Previous observations and simulations of their atmospheres suggest preferential cloud formation around the equator, associated with an equator-pole thermal gradient. Here we show that this atmospheric structure should induce latitudinal chemical variations detectable by the Doppler effect. We introduce a new method - Differential Molecular Rotational Broadening - which consists in comparing the apparent rotational broadening of individual molecules from high-resolution spectra. Application of this approach to VLT-CRIRES observations for different molecules (CO, H2O, CH4 and NH3) in the atmosphere of the brown dwarf DENIS J0255-4700 confirms the existence of latitudinal chemical variations. Our data suggest a depletion of CH4 and NH3 at low latitudes, consistent with an equatorial cloud belt. Our method could be applied to multiple brown dwarfs and exoplanets to map their atmospheres and to study various atmospheric processes.
HD 98800 is a young hierarchical quadruple system composed of the tight binaries AaAb and BaBb on a wide, highly inclined outer orbit. The B subsystem hosts a circumbinary disc in a polar configuration, and the geometry of the system offers a rare opportunity to observe the passage of the disc around BaBb in front of AaAb. We aim to update the overall orbital solution to offer precise time windows of the main occultation features by combining the orbit with the state-of-the-art knowledge of the disc's structure. We combine new and published radial-velocity measurements and multi-wavelength astrometry for the outer AB orbit and both inner subsystems, AaAb and BaBb, in a joint orbital fit. The revised solution is consistent with previous dynamical masses and improved the outer orbit, reducing the uncertainties in the period and periastron epoch by about a factor of two. It also narrows the predicted crossing-phase windows to 5--15 days at the 1σ level, improving the timing of ingress, egress, and the first cavity-crossing predictions. These results provide a more accurate timing framework for future observations of the occultation, although the predicted epochs remain model-dependent because of uncertainties in the disc structure.
The atmospheric composition of rocky exoplanets offers an important tool for constraining the properties of the interior of this type of planet, beyond what is possible from measurements of their mass and radius alone. However, the interpretation of these observations requires an understanding of the complex interplay of a larger number of coupled planetary and atmospheric processes. This review provides an overview of the current state of knowledge regarding rocky exoplanet atmospheres, beginning with their formation and escape mechanisms. We specifically highlight the importance of long-term interaction between the atmosphere, the surface, and the interior on rocky planets. Furthermore, this review addresses the influence of biological activity and photochemical reactions on the atmospheric compositions. Consequently, establishing how these different processes contribute to shaping the atmospheres of rocky exoplanets during their evolution is fundamental for the characterization of these planets with future space missions and ground-based surveys.
The majority of the Solar System planets are sources of bright radio emission, driven by energetic electrons trapped within each planet's magnetic field. Detection of this emission from exoplanets provides a unique opportunity to characterise their magnetic fields, which is key to determining the atmospheric evolution of exoplanets. However, a conclusive detection of radio emission from an exoplanet remains at large, primarily due to a lack of sensitivity at low radio frequencies. On the other hand, planet-like radio signatures have been detected on objects called ultracool dwarfs (UCDs) for over two decades. UCDs are of comparable sizes to Jupiter, but are more massive. They also possess similar interior structures to Jupiter, the region where magnetic fields are generated. Therefore, UCDs are ideal targets to study to advance our understanding of how magnetic fields manifest at planetary scales. In this Chapter, we outline the revolutionary role that the Square Kilometre Array will play in the study of exoplanets and UCDs. We anticipate that it will facilitate the first detection of radio emission from giant exoplanets with strong magnetic fields, and will deliver thousands of detections of UCDs within a few hundred parsecs. Combined with very long baseline interferometry, we also expect that astrometric monitoring will enable the detection of planets of a few Earth masses orbiting nearby radio-emitting UCDs. These findings will open a new window into how planets form and evolve in extrasolar systems.
The enrichment of heavy isotopes of volatile elements in the Martian atmosphere indicates that Mars lost a large portion of its atmosphere through escape to space. Recent atmospheric measurements by ExoMars Trace Gas Orbiter (TGO) have suggested that the vertical profiles of oxygen isotopic compositions are influenced by chemical reactions involving isotopic fractionation. However, their quantitative impacts have not yet been fully evaluated. In this study, we develop a 1D photochemical model that incorporates oxygen isotopic fractionation associated with CO$_2$ photolysis and O$_3$ formation to investigate the vertical profiles of oxygen isotopic compositions. Our calculations show that CO is depleted in heavy oxygen isotopes relative to CO$_2$, reaching $δ^{18}$O $\sim -25$ per mil and $δ^{17}$O $\sim -15$ per mil, primarily due to isotopic fractionation during CO$_2$ photolysis. The vertical profiles of oxygen and carbon isotopic compositions are in good agreement between our model and the TGO measurements. O$_3$ is strongly enriched in $^{18}$O and $^{17}$O, reaching $δ^{18}$O $\sim 100$ per mil and $δ^{17}$O $\sim 50$ per mil as a consequence of the isotopic fractionation during its formation, whereas atomic oxygen is highly depleted in the heavy oxygen isotopes with $δ^{18}$O $\lesssim -100$ per mil and $δ^{17}$O $\lesssim -50$ per mil so as to compensate for their enrichment in O$_3$. These chemical fractionation processes can deplete the heavy oxygen isotopes in species that escape from the upper atmosphere, and thereby enhance the isotopic fractionation associated with oxygen escape to space. Such fractionated isotopic compositions of escaping oxygen may be detectable by the Martian Moons eXploration (MMX) mission.
Growing interest in space interferometry for detecting bio-signatures of exoplanets has led to the development of several low-cost mission concepts involving Earth-orbiting formation flying techniques to perform nulling interferometry of incoming light from exoplanetary systems. Pursuing these developments, this work proposes a multi-objective dynamic programming scheme to optimize the balance between fuel expenditure and scientific outcomes of low-cost formation flying space interferometry missions by leveraging several insights from mission design. This scheme accounts for nominal concept of operations constraints, fuel expenditure, and observability conditions to produce a globally optimal Pareto front policy for exoplanetary system observation, and represents one of the first applications of multi-objective optimization to distributed space telescope astronomy.
Phobos evolves in a highly dynamical environment where surface-material motion is controlled by the combined effects of self-gravity, time-dependent Martian tides, and inertial forces. In such a low-gravity regime, the displacement of loose material, cannot be inferred from topographic slope alone, making a dynamical approach essential for interpreting Phobos' surface morphology and for supporting the Martian Moons eXploration (MMX) mission led by JAXA. Here, using our RAVEL code, we apply a dynamical model that combines the surface acceleration field with friction on a digital terrain model of Phobos to compute surface regolith trajectories. The model does not aim to predict the triggering of slope failure. Instead, it addresses where material would preferentially move once motion is initiated. This reveals large scale coherent dynamical regions and a sparse network of preferred regolith transport routes, termed here Regolith Migration Pathways (RMPs). The final positions of the RMPs correlate with smooth, low-relief terrains and spectrally neutral units, consistent with depositional mantles formed by long-term regolith infill, whereas rough, high-standing areas with abundant small craters and blue spectral slopes tend to correspond to dynamically active or denuded source regions. In contrast, spectrally red terrains are generally associated with dynamically quiet, morphologically rough surfaces where our model predicts negligible regolith motion, suggesting older, less frequently reworked units. Taken together, these patterns indicate that much of Phobos' surface morphology and spectral heterogeneity can be explained by long-term regolith redistribution driven by the surface acceleration field along RMPs. We provide a 3D morphodynamic atlas of RMPs across Phobos' surface, which will be useful for constraining the geographical provenance of samples to be collected by the MMX spacecraft.
We propose a novel muon Moonshot concept for lunar subsurface tomography based on upward-going muons originated from the lunar regolith. Unlike the Earth, the Moon lacks an atmosphere, leaving a dense regolith below and a near-vacuum environment above. Consequently, while most downward-going hadrons are absorbed before decaying, upward-going hadrons escaping the regolith can decay in flight, producing a significant source of lunar muons. These muons are detectable by instruments on the lunar surface or in near-lunar orbit. We perform Monte Carlo simulations to investigate their energy spectra, angular distributions, and integrated fluxes under various theoretical and detector configurations. The results indicate that the lunar muon flux is sensitive to detector altitude under a flat-terrain assumption, demonstrating its potential as a novel non-invasive probe of shallow subsurface voids. We also present case studies on detecting underground cavities and water resources, with cavity-induced flux variations observable in less than a minute and weaker water signals distinguishable after about ten minutes of data collection, and discuss potential implementations in future lunar missions.
Exoplanet atmosphere characterization heavily relies on molecular spectroscopic data. Despite efforts to obtain comprehensive spectral libraries for the chemical characterization of exoplanet atmospheres, large gaps remain, particularly for larger molecules and higher frequencies at high spectral resolution. One key example is the methane (CH4) optical spectrum. CH4, the simplest hydrocarbon, is a crucial species for exoplanet atmosphere characterization and a possible biosignature. However, until now, high-resolution linelists at optical wavelengths for CH4 have been very challenging to obtain either experimentally or computationally, leaving the high resolution spectrum of CH4 uncharacterised across most of the visible spectrum. This restricts exploration of CH4 absorption in the optical regime, as upcoming instruments such as ELT-ANDES and VLT-RISTRETTO will start probing the atmospheres of ever smaller exoplanets in optical wavelengths. To address this spectroscopic data limitation, we observed Titan's optical spectrum, dominated by CH4 absorption, at the highest spectral resolution to date with VLT-ESPRESSO. From it, we produced an empirical, low-temperature high-resolution (R ~ 190000) linelist of CH4 in optical wavelengths which we present here, with thousands of previously unidentified lines. We employ this CH4 linelist (RRS-2026) to build a template suitable for high resolution cross-correlation spectroscopy (HRCCS) studies, a first for CH4 in optical wavelengths. With this new linelist, we performed the first HRCCS detection of CH4 in the atmospheres of Titan and Jupiter using optical high resolution spectra. This work sets the stage for the search for CH4 in exoplanet atmospheres through HRCCS with current and future ground-based high-resolution optical spectrographs, showcasing how Solar System observations provide useful products for exoplanet research.
Binary stars are ubiquitous, yet it remains unclear how wide-orbit stellar companions influence the formation of hot Jupiters, particularly around M dwarfs. Here, we first report the discovery of TOI-5628Ab, a giant planet transiting a mid-type M dwarf ($M_\ast=0.36\pm0.02\ M_\odot$) every 4.34 days, accompanied by an associated white dwarf TOI-5628B ($M_{\rm WD}=0.59\pm0.16\ M_\odot$) at a projected distance of about 2,500 AU. Using TESS, ground-based photometry and SPIRou RVs, we constrain the planet radius to $0.74\pm0.04\ R_J$ and mass to $0.09\pm0.04\ M_J$, with a $3σ$ upper limit of $0.22\ M_J$. Building on this system, we further conduct a homogeneous systematic search for co-moving stellar companions with projected semi-major axis between 100 and 10,000 AU around all M dwarfs with confirmed giant planets with periods smaller than 10 days and radii larger than 0.7 $R_J$, as well as a group of field M stars with stellar properties similar to the planet sample, based on the stellar kinematics from Gaia DR3. We measure a stellar multiplicity of $34.2\pm9.5\%$ for M dwarfs hosting short-period giant planets, which is substantially higher than the fraction of $5.3\pm3.7\%$ for the field M stars by approximately a factor of 6. Our results suggest that wide-orbit stellar companions tend to promote the formation of short-period giant planets around M stars with masses $0.21 \leq M_\ast\leq 0.64\ M_\odot$, and high-eccentricity migration may play an important role in producing such systems.
We study transit light curves of exoplanets, incorporating gravitational lensing effects, and explore the possibility of estimating planetary masses using only the transit light curve. We analyze existing data for planets with relatively large orbital distances ($>1$ AU) from Kepler and TESS to evaluate how well their masses can be constrained. For six exoplanets, we derive 2$σ$ upper limits of 200-2000 Jupiter masses. We also construct mock data to assess the feasibility of determining masses for exoplanets with larger orbital radii than currently known. Our results show that a precision of $3\times 10^{-6}$ is required to recover the true masses of planets with 10 (5) Jupiter masses at orbital radii of $a=20$ (100) AU. Furthermore, we demonstrate that neglecting gravitational effects can lead to underestimation of exoplanet radii, particularly for planets with large orbital separations. Using mock data with a precision of $2\times 10^{-4}$, assuming a 30 minutes cadence with Kepler for a $K_p=14$ mag star, we find that the estimated radii can be underestimated by 1-20% for 5 and 10 Jupiter-mass planets with $a=$ 10-100 AU. These results highlight the importance of including microlensing effects when modeling transit light curves of wide-orbit planets.
The Moon Mineralogy Mapper (M3) observed a widespread 3-micron absorption feature on the lunar surface, commonly attributed to surficial OH and/or H2O. The feature varies with latitude, local time, and lunar phase, suggesting that at least part of the optically active reservoir is dynamic. We compare phase- and local-time-resolved M3 abundance estimates with a global solar wind proton implantation, hydrogen diffusion, and H2 exosphere model. The model tracks implanted H in the upper grain rim, thermally activated retention and loss, recombinative H2 release, and phase-driven source variations associated with the solar wind, magnetosheath, and magnetotail. We compute the surface abundance using a weighted activation-energy-bin representation, in which the activation-energy distribution is discretized into fixed probability-weighted bins rather than sampled stochastically by Monte Carlo particles. This implementation reduces particle sampling noise in low-flux regions. A refined grid of Gaussian effective activation-energy distributions shows that the low- and mid-latitude M3 trends are best matched by distributions centered near Ec = 0.520 eV with width Ew = 0.090 eV. The preferred case gives a combined low- and mid-latitude RMSE of approximately 27 ppm and a small mean bias of only a few ppm. These fitted parameters should be interpreted as an effective response of the M3-sensitive retained H/OH reservoir, not as a unique mineralogical activation energy. The model reproduces the dominant local-time and phase-dependent abundance structure, including the reduced retained abundance near local noon, while model-data differences near plasma-transition intervals and high-latitude behavior may reflect additional sampling, saturation, or surface-plasma processes.
Observations with the {\it James Webb Space Telescope} (JWST) can, in principle, detect moons and rotational oblateness of giant exoplanets through subtle distortions of transit light curves. The most favorable planets are expected to be on wide orbits ($\gtrsim$0.3~AU) where moons and rapid rotation are more likely to survive tidal evolution. No unambiguous detections have yet been reported. Here, we forecast the number of systems with sufficiently favorable properties to allow for secure detections, using JWST noise models, analytic detectability scalings, giant-planet occurrence rates, and the Gaia star catalog. For planets orbiting 0.9--1.6$\,M_\odot$ stars and a noise model based on demonstrated JWST performance, single-transit observations should be capable of detecting Jupiter-like rotational oblateness in several known systems and of order 10 systems yet to be discovered, if obliquities are typically $\gtrsim$10$^\circ$. A similar number of systems are favorable for Ganymede-sized moons, if such moons are common. The yields can increase to tens or hundreds of systems if lower-mass host stars are included or if JWST can achieve photon-limited performance. Time-correlated noise on 1--10 hr timescales can strongly suppress these yields; a noise floor of a few tens of parts per million is enough to hide oblateness or moons in many otherwise favorable systems. Successful searches will therefore require both a more complete census of long-period transiting giant planets and low levels of instrumental systematics and stellar variability.
Obtaining accurate fluxes for faint sources from the ground at near- and mid-infrared wavelengths is challenging because of rapidly varying atmospheric absorption, and is often limited by the lack of a nearby spectro-photometric calibrator. We present the STellar Absolute Reference Spectroscopic Flux Library (STARSFLUX), an all-sky catalogue of 64,484 calibrator spectra spanning 0.3--30 $μ$m, based on the Mid-infrared stellar Diameters and Fluxes compilation Catalogue (MDFC). STARSFLUX combines Gaia DR3 stellar parameters, multi-band photometry from space- and ground-based surveys, and synthetic NewEra PHOENIX atmosphere models. For each star, we fit the observed SED with an interpolated model spectrum, angular diameter, and extinction to produce an absolutely flux-calibrated spectrum. The derived angular diameters can also be used to calibrate interferometric observations. We validate the catalogue using three independent tests. First, STARSFLUX angular radii agree closely with Gaia DR3 radii, with a median $|ΔR|/R_{\rm Gaia,DR3} \simeq 4.8$ per cent. Second, integrated L-band (2.8--4.2 $μ$m) fluxes for 12 stars in common with the Cohen infrared standards agree with Cohen values with a mean absolute percentage difference of $4.7 \pm 2.6$ per cent. Third, comparison with Gaia DR3 BP/RP spectra shows that the near-UV/visible/near-IR absolute spectro-photometric fluxes agree to approximately 3 per cent. The FITS spectra are available at https://home.strw.leidenuniv.nl/~gamez/ and will be submitted to VizieR.
This geometric study evaluates a compact seven-petal, 6:5-resonant tulip-shaped orbit constellation for lunar south-pole positioning, navigation, and timing (PNT) and direct-to-Earth relay services. The tulip-shaped orbits are compared against elliptical lunar frozen orbit (ELFO) constellations over the NASA LunaNet Service Volume 2 (SV2), covering lunar latitudes south of -75 degrees. We compare a six-satellite tulip baseline with a minimum-cost five-satellite variant; both use the same shared three-body orbit and differ only in satellite count and along-track phasing. Performance is scored against three Initial Operating Capability C (IOC-C) metrics: line-of-sight (LOS) link availability, Lunar Augmented Navigation System (LANS) geometric dilution of precision (GDOP) below 6, and daily extravehicular activity (EVA) usable-PNT windows. Both tulip constellations satisfy all three IOC-C metrics across SV2. The six-satellite configuration meets requirements with wide margin: 75% worst-point daily GDOP availability and 18 hours of daily EVA support. The five-satellite variant also passes, but with thinner margin: 44% availability and 10 hours of EVA support. Unlike the ELFO configurations, each spacecraft in the tulip-shaped constellation maintains continuous Earth line of sight, providing persistent geometric opportunity for direct single-hop Earth relay. Because all spacecraft share a single three-body orbit, initial phasing and post-failure reconstitution reduce to along-track drift maneuvers between neighboring orbits, with screening-level estimates indicating low maneuver cost.
Solar sails present an attractive but challenging propulsion method for large-scale, long-duration trajectory design problems. In 2025, the 13th Global Trajectory Optimization Competition (GTOC13) presented a trajectory design problem involving an exoplanetary solar sailing spacecraft in the fictional Altaira system, where the goal is to collect scientific return from flybys of planets, comets, and asteroids. High-scoring solutions combine combinatorial gravity assist tour design with continuous solar sail trajectory optimization. This paper presents the solution approach developed by the team `TheAntipodes' during GTOC13. The approach combines several search and optimization stages: (1) trade studies to identify competitive entry opportunities, (2) large-scale beam search over ballistic gravity assist tours to identify beneficial planetary structures, (3) resonant targeting strategies for Vulcan flyby sequences, and (4) multi-leg solar sail trajectory refinement using sequential convex programming (SCP). A key component of the refinement process is the use of a lossless control-convex solar sail formulation, which allows for large portions of the trajectory, including all gravity assist geometry and flyby timing, to be optimized simultaneously to maximize score. The resulting trajectory placed third, with a score of 337.878 from 133 scoring flybys, and exhibited a structure broadly similar to those of the other high-scoring solutions. This demonstrates the scalability of methods such as SCP for very large trajectory design problems.
We present results from a spectroscopic survey of 32 FGK stars hosting brown dwarfs, using high-resolution optical spectra (R = 130,000 and 50,000) obtained with the PEPSI spectrograph on the Large Binocular Telescope. The primary goal of this survey is to determine precise stellar parameters and abundances for 11 elements (C, O, Mg, Si, Ca, Al, Ti, Fe, Y, S, and N) in these systems. We employ spectral synthesis within the BACCHUS framework to derive precise stellar properties and elemental abundance ratios. For our average S/N $>$ 200 data, we achieve a typical error of 42 K in T$_\mathrm{eff}$ and $\sim$0.03 dex for [Fe/H]. We observe a significant dispersion from a solar C/O ratio among the sample of brown dwarf host stars that host primarily wide-orbit brown dwarfs. Using established theoretical chemical frameworks, we discuss the implications of the observed Mg/Si and Ca/Al ratios for cloud properties in the brown dwarf companions. Finally, we evaluate the applicability of the [Y/Mg] stellar clock for our sample and discuss the broader implications of our results. This work provides a timely and uniform abundance analysis of host stars, supporting extended wavelength brown dwarf observations in the era of JWST.