Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics Earth and Planetary Astrophysics

astro-ph.EP Watch Earth and Planetary Astrophysics (sign in to use)

Wednesday 23 September

Earth and Planetary Astrophysics

Europa smooth plains formation: a SPH cryovolcanic model

Cryovolcanism is a fundamental geodynamical process that has been observed on several icy moons. Europa stands out as a primary target for investigating potential cryovolcanic activity and its role in shaping icy moon surfaces. This work characterizes the physical state of ascending cryomagma and evaluates the relevance of terrestrial analogues to explain Europa's surface features. After investigating cryomagma depressurization, we developed CryoPy, a numerical model based on Smoothed-Particle Hydrodynamics (SPH) for computational fluid dynamics and heat transfer. We simulated flow propagation and surface cooling using this model. Our simulations show that snow is a relevant physical analogue for cryolava, as the resulting flows closely match the smooth plains observed by Galileo in both extent and morphology. Additionally, we established the temporal thermal signatures of these events. Such active cryovolcanic events are within the detection limits of, primarily, the E-THEMIS infrared imager, providing a predictive framework for the Europa Clipper mission.
Read the full abstract

Earth and Planetary Astrophysics

Kelvin-Helmholtz vortices in the Martian +E hemisphere downstream of crustal magnetic anomalies

At unmagnetized planetary bodies, the Kelvin-Helmholtz (KH) instability at the solar-wind boundary is a key driver of atmospheric escape, fundamentally shaping long-term climate evolution. This shear-driven plasma mixing is theoretically predicted to be asymmetric, strongly favoring the hemisphere where the solar-wind electric field points toward the planet (-E hemisphere). Here we report the observation of a coherent train of fully developed KH vortices in the traditionally unfavorable +E hemisphere of Mars, using in-situ MAVEN measurements. Contemporaneous upstream constraints from Tianwen-1 confirm a steady interplanetary magnetic field, indicating sustained shear-driven growth rather than a transient boundary response. We find that these vortices developed within an expanded, weak-field boundary layer downstream of strong southern crustal magnetic anomalies. Our analysis suggests that regional boundary structuring can override the large-scale hemispheric preference to facilitate localized KH growth. These findings demonstrate that shear-driven escape pathways can operate in traditionally unfavorable regions, altering our understanding of atmospheric volatile loss at weakly magnetized planets.
Read the full abstract

All authors Junxi Niu, Tong Dang, Jiuhou Lei, Binzheng Zhang, Wenya Li, Binbin Tang, Junjie Chen, Sudong Xiao, Tielong Zhang

Earth and Planetary Astrophysics

Dynamical Behavior and Surface Features of Large Super-Fast Rotator Near-Earth Asteroid (436724) 2011 UW158

The sub-kilometer Super-Fast Rotator and Potentially Hazardous Asteroid (436724) 2011 UW158 exhibits an extreme dynamical environment driven by its rotational period of 0.61072 h. Using a 3-D polyhedral shape model, we investigate the surface dynamic stability under both cohesionless and cohesion-modified conditions, while also addressing families of periodic orbits, station-keeping maneuvers, and energy states in this rapidly rotating, irregular-body regime. Our results show that, with the inclusion of cohesion, the effective slopes decrease to ~45{\deg}, permitting the retention of boulders with maximum sizes of ~50 m near the equatorial region. The dominance of the centrifugal term in the effective potential shifts the exterior equilibrium points inward, creating highly unstable geopotential regions where loose regolith is susceptible to ejection, preventing long-term surface retention outside the polar regions. To quantify this instability, we introduce the equivalent instability speed ($v_{\mathrm{ei}}$), extending conventional escape speed analyses to the SFR regime. These results are consistent with significant internal cohesion in 2011 UW158, which enables the body to maintain structural integrity under rotational stresses exceeding the classical spin barrier. Additionally, thermophysical modeling reveals that the super-fast rotation produces a nearly uniform surface temperature distribution. Finally, we apply a spherical harmonics method to compute families of periodic orbits around the asteroid to assess the fuel cost of station-keeping maneuvers at reduced computational cost. Overall, this study characterizes the dynamical environment of 2011 UW158, contributing to the understanding of Super-Fast Rotator asteroids recently identified in LSST survey images.
Read the full abstract

All authors Caio Gomes, Andre Amarante, Filipe Monteiro, Alessandra Ferreira, Leonardo Braga, Allan Kardec de Almeida Jr., Duncan Lyster, Nicoli Rocha, Gabriel Damilano, Leonardo Barbosa

Earth and Planetary Astrophysics

The Dynamical Environment of Kepler-223 Planetary System

The Kepler space telescope was a successful mission that used the transit method. Multi-planetary systems are interesting because they are dynamically rich due to the planets' interactions. Kepler's mission was the first to find a system with four planets all in resonance, that is Kepler-223. The Kepler-223 system is of unique dynamical interest as the first exoplanetary system confirmed to host four sub-Neptune planets in a precise 3:4:6:8 resonant chain. Previous dynamical studies have struggled to fully recover the long-term phase protection mechanisms, often finding that higher-order resonant angles circulate rather than librate. In this work, we re-examine the dynamical state of Kepler-223 by performing high-resolution N-body integrations and constructing detailed stability maps in the (a,e) phase space, we demonstrate that the system is locked in a deeper resonant state than previously thought. Crucially, we find that all resonant arguments, exhibit robust libration over secular timescales. Our stability maps reveal that the planets reside in clearly defined stability islands that allow for higher eccentricities than previous estimates. We further compare these numerical islands with analytical resonance widths and find excellent agreement for the inner three planets, while the outermost planet (Kepler-223e) exhibits a stability region wider than first-order analytical theories predict, suggesting complex multi-body stabilization effects at the edge of the chain.
Read the full abstract

Tuesday 22 September

Earth and Planetary Astrophysics · Cross-listed from astro-ph.IM

South Pole-Aitken basin sample Return and eXploration (SPARX) Science Definition Team Report

Sample return from South Pole-Aitken Basin (SPA) on the Moon has been recognized as the highest priority lunar science mission for well over two decades, appearing as a recommended New Frontiers class mission in the two previous Planetary Science Decadal Survey documents (2003-2012 New Frontiers, 2013-2022 Visions and Voyages), and as a prioritized mission for the Lunar Discovery and Exploration Program (LDEP) in the 2023-2032 Origins, Worlds, and Life Decadal Survey. The South Pole-Aitken basin presents the ideal target destination to test nearly 60 years of lunar science hypotheses owing to a unique combination of its size, antiquity, and location on the lunar farside. Convened by NASA in response to OWL Decadal Survey recommendations, the SPARX report provides analysis on prioritized science objectives, measurement requirements, implementation architectures, and a design reference mission for a South Pole-Aitken basin sample return mission. The mission profiles presented in this report leverage key advancements in remote surface exploration technology, the U.S. commercial space industry, and the programmatic focus on sustained human presence and exploration at the lunar south polar region. The science objectives and measurements presented in this report outline a mission that would transform lunar science, would address multiple high-level Decadal science questions at the only place in the Solar System where they can be credibly addressed, would pioneer a new paradigm for planetary surface exploration in the 21st century, would leverage the burgeoning space economy to pursue ground-breaking science on a more cost-efficient basis, and would unite robotic exploration with human exploration in a model for a future sustained lunar enterprise.
Read the full abstract

All authors Lauren M. Jozwiak, James T. Keane, Robert Citron, Natalie Curran, Steve Elardo, Caleb Fassett, Abigail Fraeman, Saira Hamid, Lena Heffern, Peter James, Daniel Moriarty, Kevin Righter, Jennifer Whitten, Hunter Vannier, Zoe Wilbur, Tomokatsu Morota, Hiroshi Nagaoka, Caroline-Emmanuelle Morisset, Patrick Hill, James Carpenter, Francesca MacDonald

Earth and Planetary Astrophysics · Cross-listed from astro-ph.IM

Increasing Sensitivity to Trailed Solar System Objects in Archival JWST NIRCam Imaging with Group Differencing

The JWST archive contains thousands of NIRCam exposures distributed across the sky. Small solar system bodies (SSSBs) are present in many of these fields, but are often trailed and more difficult to detect. The standard pipeline produces a single image by fitting a slope to a sequence of $N$ cumulative “up-the-ramp" groups of the detector. Thus, each JWST image comprises a time series of cumulative reads. We demonstrate that transforming the ramp into a series of group difference images reduces the per-image trail length of moving sources, producing an $(N-1)$-group photometrically calibrated time series from an $N$-group integration. This enables higher-cadence trajectory-based searches using as few as a single JWST exposure. A generalized least-squares fit of the group difference images along the trajectory of a moving source recovers the source flux at higher SNR when compared to the ramp fit, with the improvement scaling with the reduction in trailing. We validate our procedure on 13 moving objects detected in archival JWST NIRCam images. To estimate our sensitivity increase, we perform an injection/recovery simulation on the JWST imagery. Trailed moving objects recovered via group differencing show higher SNRs, improving completeness to moving sources compared to standard calibrations. For typical integration lengths of $N=4$ or $6$ for archival JWST/NIRCam exposures in the ecliptic, group differencing yields sensitivity gains between $\sim 1.1$–$1.4$ mag, based on simulation results. This method significantly increases the sensitivity of archival JWST surveys to trailed SSSBs.
Read the full abstract

Earth and Planetary Astrophysics · Cross-listed from astro-ph.SR

A Multiband Catalog of Infrared-Excess Sources in the Vela-Carina Region: VVVX, DECaPS, and GLIMPSE

Multiwavelength surveys of the Galactic midplane enable the identification of infrared-excess sources, signatures of circumstellar material, though high extinction and crowding hinder detection. We construct a homogeneous catalog of infrared-excess candidates in a $\sim$7.0 deg$^2$ field of the Vela-Carina region, characterizing their excess emission via spectral energy distribution (SED) analysis. We combine optical photometry from DECaPS, near-infrared photometry from VVVX, and mid-infrared photometry from Spitzer/GLIMPSE using hierarchical positional cross-matching. Candidates are selected via photometric quality and MIR color diagnostics, and their SEDs are analyzed with VOSA using theoretical atmosphere model grids. Photogeometric distances are adopted, while Gaia DR3 and NEOWISE variability provide complementary diagnostics. From the matched 249,719-source VVVX-GLIMPSE-DECaPS catalog, we identify 1,512 MIR-excess objects, of which 191 have complete 12-band photometry with S/N > 10, enabling detailed SED analysis. Photogeometric distances are statistically consistent with geometric distances for nearly the entire sample. SED analysis independently supports the MIR-excess nature of most sources, indicating emission beyond stellar photospheres. About 15% of sources are flagged variable by Gaia, while NEOWISE light curves reveal MIR variability of a few tenths of a magnitude. Comparisons with published YSO and variability-selected catalogs show limited overlap, particularly with variability-selected samples, indicating that MIR-excess selection complements other approaches to identifying candidate YSOs in the Vela-Carina region. The resulting catalog provides a well-defined MIR-excess sample with homogeneous multiband photometry, SED characterization, distances, and variability diagnostics, forming a foundation for future spectroscopic and time-domain studies.
Read the full abstract

All authors R. M. Torres, R. K. Saito, V. Fermiano, P. Esteves, B. W. Borges, D. Quispe, D. Minniti

Earth and Planetary Astrophysics · Cross-listed from astro-ph.IM

Laboratory Evidence that Methanol-Rich Ice Mantles Lower Methyl Formate Binding Energies

Complex organic molecules are widely detected in star-forming regions and are important precursors to prebiotic species. Because these molecules are thought to form primarily within icy grain mantles, their desorption is assumed to be regulated by the sublimation of the surrounding H$_2$O ice matrix. However, observations increasingly reveal gas-phase abundances of methyl formate (MF) at temperatures well below the sublimation temperature of water ice (T ~ 100K), challenging this view. To investigate how ice composition and phase influence MF desorption, we present laboratory temperature-programmed desorption measurements of MF from four astrophysically relevant substrates: amorphous and crystalline H$_2$O and CH$_3$OH. Desorption kinetics were analyzed using pre-exponential factors derived from transition-state theory with leading-edge analysis and non-negative least-squares inversion to recover binding-energy distributions. We find that MF binds more weakly to CH$_3$OH ice than to H$_2$O-rich substrates, with representative binding energies of ~6247K for amorphous H$_2$O, ~6213K for crystalline H$_2$O, and ~5469K for amorphous CH$_3$OH. MF deposited on crystalline CH$_3$OH exhibits coverage-independent leading edges consistent with multilayer or island desorption, yielding an effective binding energy of 5333K. Enhanced surface mobility and reduced trapping efficiency on methanol ice, combined with weaker MF binding, provide a natural explanation for the early appearance of gas-phase MF observed in star-forming environments, with characteristic desorption temperatures of ~85K on CH$_3$OH and ~95-100K on H$_2$O. Together, these effects support a two-step thermal desorption pathway in which MF associated with CH$_3$OH-rich ice can be released prior to bulk H$_2$O sublimation.
Read the full abstract

All authors Rachel E. Gross, Jeroen Terwisscha van Scheltinga, L. Ilsedore Cleeves, Catherine A. Dukes, Robin T. Garrod, Alexia Simon, Shane Sawyer, Adam Woodson

Earth and Planetary Astrophysics

Revisiting planetesimal formation from mm-sized grains in pressure bumps

We present evidence that the streaming instability (SI) can form planetesimals from millimeter grains inside axisymmetric pressure bumps, provided the bump is not too strongly reinforced. We conducted three-dimensional local shearing-box simulations of millimeter grains in a Gaussian pressure bump of amplitude $A = 0.6$, at a resolution of $160/H$. We varied only the Newtonian reinforcement timescale $t_{\rm reinf}$, the timescale on which the imposed bump is restored against particle back-reaction. Our reference run T1 uses $t_{\rm reinf} = 1\,\Omega^{-1}$. Particles pile up in a narrow axisymmetric band at the minimum-headwind point, cross the Roche density, and collapse through direct gravitational instability (GI). In our weakly reinforced run T100 with $t_{\rm reinf} = 100\,\Omega^{-1}$, the outcome is different. Particles drift through the bump roughly three times more slowly, and dense SI-like filaments develop across a wide radial region well before the Roche density is reached. These results are consistent with a residence-time criterion, which states that for the SI to grow, $t_{\rm cross} > t_{\rm grow}$, where $t_{\rm cross}$ is the time a particle spends inside the region where the ratio $Z/\Pi$ of the solid abundance $Z = \Sigma_p/\Sigma_g$ to the headwind parameter $\Pi = \Delta v/c_{\rm s}$ is high enough for strong clumping, and $t_{\rm grow}$ is the SI growth time. T1 violates this criterion, and T100 satisfies it. Our results suggest that the reinforcement timescale matters for millimeter grains because it helps set the crossing time, and therefore whether the SI has time to grow. We caution, however, that our reinforcement scheme is an idealised numerical construction rather than a physical bump-forming mechanism such as a planet, and a more definitive test of this picture will require simulations with a more physically motivated bump.
Read the full abstract

All authors Tabassum S Tanvir, Jacob B Simon, Daniel Carrera, Jeonghoon Lim

Earth and Planetary Astrophysics

M Dwarf Singles and Multis Have Dissimilar Metallicities Indicating Different Origins

It is unclear if host star metallicities differ between systems with one detected transiting planet ("singles") or multiple detected transiting planets ("multis") around FGK dwarfs. Here, we find strong evidence that they do for M dwarfs. We use a homogeneous sample of FGK and M dwarf chemical abundances to demonstrate that M dwarf multis are significantly more metal-poor than singles ($p$ = 0.00079), even when we only consider small (<4 $R_{\oplus}$) planets ($p$ = 0.0013). This result emerges at the late-K to M dwarf transition. We observe that M dwarf multis become more metal-poor as a function of increasing planet multiplicity. Additionally, planets in singles have significantly higher eccentricities and shorter orbital periods compared to planets in multis. We interpret these findings as evidence that M dwarf singles lead relatively active dynamical lives, and the different planetary architectures of singles and multis arise from distinct dynamics influenced by metallicity. This may not occur in FGK systems because their more massive protoplanetary disks may supply enough planet-forming material to instigate dynamics that lead to only one transiting planet, even when the host star metallicity is low. Another possibility is that multis are less compact around FGK dwarfs because they form from larger disks, which may often cause them to be observed as singles in transit data.
Read the full abstract

All authors Aida Behmard, Sheila Sagear, Yuxi Lu, Romy Rodríguez Martínez, Adrian M. Price-Whelan, Soichiro Hattori

Earth and Planetary Astrophysics

Simulated Effects of Radioactive Elements on Magnetic Evolution and Atmospheric Mass Loss of Earth-like Planets

Thousands of exoplanets have been discovered in recent decades, yet Earth remains the only known world with life. To assess the potential habitability of other planets, we must identify the importance of various physical conditions that may contribute to sustaining life. One such factor is the presence of a planetary magnetic field. A planet's magnetic field can protect its atmosphere from photoevaporation, and its long-term strength depends on interior heat flow driven by radiogenic heating from isotopes such as 40K, 232Th, 235U, and 238U. In this work, we use the VPLanet code to simulate the coupled thermal, magnetic, and atmospheric evolution of Earth-like planets with varying radiogenic heat budgets orbiting G-, K-, and M-dwarf stars. The range of isotope abundances tested is based on observed stellar elemental distributions for planet-hosting stars in the Milky Way. We find that variations in radiogenic heating slightly extend dynamo lifetimes and reduces atmosphere and surface water loss for K- and G-dwarf systems, while its effects are more negligible for M dwarfs. Our results show that, while stellar irradiation and planetary location primarily control photoevaporation rates, radiogenic heating contributes to planetary habitability in a small subset of cases and may help refine target selection for future surveys.
Read the full abstract

Earth and Planetary Astrophysics

JWST and HST Confirm Transit Timing Variations for the Ultra-Hot Jupiter TOI-2109 b

Ultra-hot Jupiters are a rare yet valuable subset of exoplanets that can help elucidate both the internal properties of their host stars as well as the orbital evolution of close-in planets. An extreme among extremes, TOI-2109 b is a super-Jupiter ($1.347 \pm 0.047 R_{\rm Jup}$, $5.02 \pm 0.75 M_{\rm Jup}$) with the shortest orbital period (0.67247479 $\pm$ 0.00000028 days) and the fastest predicted orbital decay rate of any known giant planet. Previous observations hinted at a quadratic long-term trend in the transit timing variations (TTVs) consistent with orbital decay, as well as short-term sinusoidal TTVs consistent with a near-resonant planetary companion that would indicate a dynamically quiescent migration pathway. Here, we combine archival ground-based, TESS, and CHEOPS photometry with new transit and eclipse data from HST and JWST to update the TTV analysis. We find very strong evidence for sinusoidal TTVs over a linear ephemeris ($\Delta$BIC of 33.6), consistent with a nearby planetary companion as suggested by J.-V. Harre et al. 2024, and rule out several non-planetary scenarios that could explain the sinusoidal TTVs. We also find strong evidence of a quadratic TTV trend ($\Delta$BIC of 6.75) consistent with an orbital decay rate of -3.46 $\pm$ 1.17 ms/yr and a stellar tidal quality factor of 2.13$^{+1.18}_{-0.61} \times10^7$, though we cannot rule out line-of-sight acceleration as the cause of the quadratic trend. If our results are confirmed, TOI-2109 b would be the only known exoplanet besides WASP-12 b with a decaying orbit and the most extreme of an ultra-rare class of close-in giant planets with nearby companions.
Read the full abstract

All authors Michael Greklek-McKeon, Julie Inglis, Nicole L. Wallack, Shreyas Vissapragada, Ian Wong, John R. T. Allen, Yayaati Chachan, Fei Dai, Guangwei Fu, Tiffany Kataria, Thaddeus D. Komacek, Johanna K. Teske

Earth and Planetary Astrophysics

Revealing the Shiny Nature of Lava Worlds: Small Planet Phase Curves with TESS

Recent thermal emission observations suggest that `lava worlds'–Earth-sized exoplanets hot enough to melt their silicate surfaces–have atmospheres on a population-level. This evidence comes from measured dayside brightness temperatures much lower than expected for a dark, maximally hot bare rock. However, these results are degenerate between a thick atmosphere advecting heat to the planet's nightside or a high albedo reflecting incoming radiation. This degeneracy can be resolved through visible light phase curve observations, which can constrain how much light is reflected by the planet. Here, we present a population-level analysis of the visible-light TESS phase curves of $R\lesssim5\,R_{\oplus}$ exoplanets. We select the population by defining a `Reflection Spectroscopy Metric' meant to quantify the relative signal-to-noise of reflected light observations. We report new $>3\,\sigma$ detections of the TESS phase curve signals of LTT 9779 b and the lava worlds TOI-2431 b, TOI-561 b, 55 Cancri e, K2-141 b, and TOI-1444 b. We use these results and previous results from Kepler and CHEOPS to homogeneously derive geometric albedos for these planets, revealing that lava worlds generally show high albedos compared to the low geometric albedos measured for hot Jupiters at similar temperatures. These results support an emerging picture of reflective silicate clouds on these worlds cooling their daysides. The upcoming PLATO mission is expected to tightly constrain the reflected light phase curves of lava worlds.
Read the full abstract

All authors Brandon Park Coy, Michael Zhang, Fei Dai, Qiao Xue

Earth and Planetary Astrophysics

Characterization of Planetary Systems with Astrometry

Astrometry is the oldest method used to search for extrasolar planets. This technique provides key information about a planetary system, like the true mass of planets and its three-dimensional orbital architecture. Nowadays, the space-based mission Gaia, on one hand, and ground-based interferometry at infrared and radio wavelengths, on the other hand, can achieve astrometric accuracies of tens of micro-arseconds. Thanks to this unprecedented level of accuracy, the number of astrometrically discovered exoplanets has started to grow. In this chapter, I review the basic principles for orbital characterization of unseen planetary companions using astrometric data. I also review the properties of the exoplanets discovered so far with astrometric methods and comment on the contributions that this technique can make to the exoplanet field.
Read the full abstract

Earth and Planetary Astrophysics

High-resolution optical spectroscopy reveals the distinctive volatile composition of the interstellar comet 3I/ATLAS

Comets preserve volatile material from the early stages of planetary-system formation. The discovery of the interstellar comet 3I/ATLAS enables a direct comparison between a comet formed around another star and comets in the Solar System. We characterise the coma composition of 3I/ATLAS across perihelion using homogeneous high-resolution near-UV and optical spectroscopy with UVES on the ESO Very Large Telescope, from August 2025 (r$_h$=3.14 au) pre-perihelion to February 2026 (r$_h$=4.34 au) post-perihelion. From flux-calibrated spectra, we measured OH, NH, CN, C$_3$, CH and C$_2$ emissions, derived production rates, examined abundance ratios and heliocentric-distance trends, constrained NH$_2$, and compared OH-based water production with estimates from forbidden oxygen lines. Gas production rates show broadly symmetric pre- and post-perihelion behaviour, with steeper heliocentric-distance dependences than typically observed in Solar System comets. NH and NH$_2$ are strongly depleted, indicating a severe deficiency of ammonia-related volatiles. The C$_2$/CN ratio indicates carbon-chain depletion and evolves with heliocentric distance, consistent with coma-driven effects such as hyperactivity rather than intrinsic compositional changes. Water production from forbidden oxygen lines exceeds OH-based estimates except near perihelion, suggesting varying contributions from H$_2$O and CO$_2$ to coma chemistry. Relative abundance correlations place 3I/ATLAS outside the canonical Solar System comet cluster and closest to strongly carbon-chain depleted comets, including G-Z-type comets, while showing strong depletion of NH- and C$_3$-bearing species relative to CN. 3I/ATLAS therefore extends the known diversity of cometary compositions, demonstrating that volatile reservoirs in other planetary systems can differ substantially from those sampled by Solar System comets.
Read the full abstract

All authors K. Aravind, E. Jehin, C. Opitom, J. Manfroid, D. Hutsemékers, D. Bodewits, M. Bannister, R. C. Dorsey, L. Ferellec, H. Kawakita, M. Lippi, C. Snodgrass

Earth and Planetary Astrophysics

NGTS clusters survey - VII. An Enigmatic Short-Period Circumsecondary Disk Candidate in Orion

Young stellar systems provide key insights into the processes governing star/planet formation, often exhibiting complex and evolving photometric variability. In this paper, we present an analysis of the object around NOI 105872, observed by the Next Generation Transit Survey (NGTS) to eclipse a young ($4.1\pm0.4$ Myr) M-type star in the Orion Nebula Cluster. The system displays deep, asymmetric eclipses with a period of $\sim0.69$ days, alongside variability in morphology. We combine photometric observations spanning approximately a decade from multiple ground- and space-based facilities to investigate the nature of this system. We explore a range of physical scenarios to reproduce the eclipse profiles and their evolution. After analysis, our preferred explanation is a circumsecondary disk surrounding a companion in orbit around the M-type primary. We adopt a six-parameter disk model and fit the eclipse profiles independently to each dataset. The best-fitting solution yields a disk radius of $\sim 0.005 \mathrm{au}$ or $0.98 R_{*}$, with broadly consistent orientations across epochs, but with variations that may be attributed to disk precession or intrinsic stellar variability. Constraints from the Hill radius imply a companion mass consistent with a very low-mass stellar companion. These results establish NOI 105872 as a candidate for a circumsecondary disk system in a young stellar environment. Continued photometric monitoring and spectroscopic follow-up will be essential to confirm this interpretation. If confirmed, systems such as NOI 105872 provide valuable benchmarks for understanding disk dynamics, companion formation, and the origin of complex eclipse phenomena in young stars.
Read the full abstract

All authors Niamh Mallaghan, Ernst J. W. de Mooij, Christopher A. Watson, Edward Gillen, Matthew Kenworthy, Louise D. Nielsen, David R. Anderson, Matthew P. Battley, Edward M. Bryant, Matthew R. Burleigh, Benjamin M. J. Cadell, Sarah L. Casewell, Alexander Chaushev, V. S. Dhillon, Elsa Ducrot, Jorge Fernandez Fernandez, Samuel Gill, Michael Gillon, Michael R. Goad, Faith Hawthorn, Katlyn L. Hobbs, James McCormac, Maximiliano Moyano, Catriona A. Murray, Toby Rodel, Suman Saha, Ramotholo R. Sefako, John Southworth, Mathilde Timmermans, Amaury H. M. J. Triaud, Jose Vines, Richard G. West, Peter J. Wheatley, Tafadzwa Zivave

Earth and Planetary Astrophysics

Molecules in the shadows of the protoplanetary disk HD 143006: Evidence for C/O $\gtrsim 1$ from simultaneously modeling fourteen ALMA lines

The inner disk of the T Tauri star HD 143006 is inclined and casts a broad shadow on a large portion of the outer disk. We use this extraordinary geometry to study the impact of stellar ultraviolet (UV) irradiation on disk chemistry, temperature, and line formation. We present new ALMA observations of this disk in CO, $^{13}$CO, C$^{18}$O, CN, HCN, HCO$^+$, CS and C$_3$H$_2$ and combine these with previously-published data to get a rich dataset, which is modeled with ProDiMo models for the illuminated and shadowed side by including an optional inclined inner disk to partially block stellar illumination of the outer disk. The observed line flux contrasts between the illuminated and the shadowed disk sides are small. We explain this by introducing “equilibrated models” which consider a slow heating/cooling relaxation of gas parcels orbiting in and out of the shadow. We fit the radial line intensity profiles of all 14 ALMA lines for both halves of the disk, the SED beyond 20$\,\mu$m, and two ALMA continuum images, with an overall reduced $\chi$ of about 2.3, in one combined model. Line intensities only increase slightly on the illuminated side because the gas warms up slowly, although fast UV photochemistry significantly reduces the observable molecular concentrations. For the lines of N$_2$H$^+$ and CS, which probe deeper layers where the gas and dust temperatures are coupled, slow UV-desorption of ices such as NH$_3$ and H$_2$S plays a crucial role. This can result in a brightening of these lines on the shadowed side. The carbon-to-oxygen ratio in the outer disk of HD 143006 is slightly greater than one ($\sim 1.05$). The disk mass is fitted to be as low as $0.0012 M_{\odot}$, but almost equally well-fitting models with a fixed disk mass of $0.021 M_\odot$ can be obtained after re-adjusting the disk shape, dust, and gas heating/cooling parameters. (abridged)
Read the full abstract

All authors S. E. van Terwisga, P. Woitke, A. M. Weiss, A. D. Bosman, L. Trapman

Earth and Planetary Astrophysics

Transport and Thermochemical Kinetics of Alkali Species in Hot Jupiter Atmospheres: Implications for Magnetic Models

Interactions between the planetary magnetic field of a hot Jupiter and the winds within its atmosphere are moderated by the electron and ion abundances. For a solar metallicity atmosphere at pressures between 10 and $10^{-4}$ bar, the primary sources for these charged species are expected to be the potassium and sodium atoms, and models of magnetic effects in hot Jupiter atmospheres have typically assumed that charged species exist within the atmosphere in their equilibrium abundances. In this work, we investigate the thermochemical kinetics of alkali species within a hot Jupiter atmosphere, thus allowing for the thermal ionisation fraction to depart from equilibrium, and demonstrate that quenching of ions occurs between 1 and 10 mbar. This increase in ionisation fraction on the nightside and at mid-latitudes results in stronger magnetic interaction, here modelled using the magnetic drag approximation, further slowing the jet and altering the global circulation. We additionally show that while quenching occurs, large magnetic resistivities resulting from low electron abundances continue to exist on the nightside between 10 and 100 mbar, in tension with magnetohydrodynamic models that assume horizontally uniform magnetic resistivities to simplify the explicit integration of the induction equation. At lower pressures ($\sim$ 0.1 mbar) where the electrons have become horizontally homogenised, the magnetic Reynolds number approaches or exceeds unity, signalling the breakdown of the magnetic drag approximation. We thus stress that chemical kinetics represents a dynamically important component of modelling hot Jupiter atmospheres while at the same time further complicating the treatment of magnetic effects.
Read the full abstract

All authors D. A. Christie, T. M. Evans-Soma, N. J. Mayne, E. Hébrard, M. Zamyatina, J. E. Owen, K. Kohary

Earth and Planetary Astrophysics

Molecular mapping of an exoplanet with JWST: NH3 detection in the temperate super-Jupiter Epsilon Indi Ab

Epsilon Indi Ab is currently the coldest, and one of the closest, directly imaged exoplanets (T_eff = 275 K, d=3.6 pc), offering a unique opportunity to test our understanding of the physical and chemical properties of gas giants across the full temperature range spanned by hotter directly imaged exoplanets and the cooler giants of the Solar System. We aim to characterize the atmosphere of Epsilon Indi Ab using infrared spectroscopic observations. We present JWST spectroscopic observations of Epsilon Indi Ab obtained with MIRI in the Medium Resolution Spectrograph mode (4.9-27.9 microns, October 2025) and NIRSpec in the Integral Field Unit mode (2.87-5.27 microns, May 2026). The data are highly contaminated by the stellar light from the bright host star. We therefore adopt a high-pass filtering and cross-correlation approach based on petitRADTRANS atmospheric models to detect the planetary signal and specific molecules. The cross-correlation analysis enables a robust detection of the planet with both instruments (up to S/N = 15.2). Using molecular mapping, we securely identify the presence of NH3 (maximum S/N = 15.6), H2O (maximum S/N = 10.8), CH4 (S/N = 3.8) and marginally detect CO2 (S/N = 2.5). These first results confirm the potential of infrared spectroscopy combined to cross-correlation techniques to characterize the atmosphere of directly imaged exoplanets with JWST, even the in case where the contamination by the host star is critical.
Read the full abstract

All authors R. Kravtchenko, O. Berné, I. Schroetter, P. Amiot, F. Debras

Earth and Planetary Astrophysics · Cross-listed from astro-ph.IM

Fitting Moving Objects in Up-The-Ramp Data with Applications to the Roman Space Telescope and JWST

A moving object breaks the fundamental property of constant per-pixel count rates in an astronomical image read out up-the-ramp. In this paper, we show how to fit a moving object's path across a detector as that detector is read out nondestructively. We write the full likelihood function for every pixel subject to a constant count rate plus a time-dependent count rate due to a moving source. Assuming the moving source to be point-like and assuming the effective point-spread function to be known, we are left with four parameters that enter the likelihood nonlinearly: two for position and two for velocity. All remaining parameters can be optimized using closed-form expressions. Our approach extracts maximal information on a moving source's position and speed and enables the source to be accurately removed from the image. We investigate the dependence of flux, position, and velocity precision on the target's speed and the readout pattern. We also find a small, positive bias on the recovered flux due to the need to fit for an uncertain position and speed. Our approach can be used for space-based images with minor Solar system bodies in the foreground, e.g.~from Roman and JWST, or for ground-based observations with satellites in the foreground. We demonstrate the promise of our method with a fit to an asteroid track observed serendipitously by the NIRISS instrument on JWST, comparing it to the performance of the JWST pipeline. Python code implementing our approach is available at https://github.com/t-brandt/moving_source. The total computational cost to fit the track of a moving object is $\sim$1 second on a 2023 Macbook Pro.
Read the full abstract

Earth and Planetary Astrophysics · Cross-listed from math.DS

Identifiable Uniqueness of Impulsively Forced, Inert, and Stabilized Body Trajectories

An unusual object in the sky raises two related questions: what kind of object it is and what physical process drives its motion. Departures from a gravity-only trajectory can provide evidence, but the same apparent acceleration may be generated by distinct dynamical mechanisms. We formulate this as an identification problem and establish conditions for identifiable uniqueness among three reference classes: inert dynamics, impulsive forcing, and active stabilization such as a controlled spacecraft. The distinction rests on how each class responds to disturbances. Although different mechanisms can produce similar average trajectories, the specified reference models imply different stochastic responses around those trajectories. We show when these differences uniquely determine the dynamical class from the observed trajectory and when several classes remain observationally compatible.
Read the full abstract

Earth and Planetary Astrophysics

Possible Observational Survey Strategies to Maximize Exoplanet Yields: A Report from the Survey Strategies Task Group in the Exoplanet Science Yields Working Group

We present a variety of possible observational sequences from planet detection to characterization (i.e., survey strategies) for the initial exoEarth detection and characterization survey phase of the Habitable Worlds Observatory. An array of survey strategies were collated by the HWO Exoplanet Survey Strategies Task Group; by investigating and comparing expected resultant exoEarth yields for multiple survey strategies, we can better understand the most efficient strategy to characterize multiple planets in unconventional ways. We present eight varied survey strategies, presented as steps of Detection, Photometry, Orbit Measurement, and/or Characterization, including discussion of possible shortcomings and benefits of each strategy. We present exoEarth yield calculations from the Altruistic Yield Optimizer (AYO) and the Exoplanet Open-Source Imaging Mission Simulator (EXOSIMS) for select strategies, with initial considerations of thermal emission and detector noise. Herein, a yield is the expectation value given the astrophysical assumptions of the number of planets that meet the observational criteria during the survey. We find that AYO and EXOSIMS find the highest yields for H2O at 0.9 microns, although the calculated yields vary. When allowing for wavelength optimization within AYO, we find that H2O is most optimal to observe at 0.9 microns when given the full range of VIS and NIR values, with an expected exoEarth yield of 29.1. H2O and O2 dual characterization is possible at the expense of yield planet loss, with a calculated possible yield of 20.9 in the VIS. The characterization of CH4 and CO2 in the NIR result in the lowest exoEarth yields at 19.9 and 6.4 possible expected planets, respectively. Future work updating AYO and EXOSIMS is necessary to explore further survey strategies that follow different observation procedures, some of which is currently underway.
Read the full abstract

All authors Natasha Latouf, Corey Spohn, Kevin Fogarty, Chris Stark, Dmitry Savransky n=and Eleonora Alei, Anthony Boccaletti, Gagandeep Kaur, Avi Mandell, Rhonda Morgan, Aki Roberge, Jean-Baptiste Ruffio, Jean Schneider, Edward Schwieterman, Armen Tokadijan, Daniel Viudez-Moreiras

Earth and Planetary Astrophysics

Large-scale chaos in restricted hierarchical triples driven by short-range forces

Context: The eccentric von Zeipel-Lidov-Kozai effect, which is widely applied to diverse astrophysical settings, can drive the inner binary to extremely high eccentricities, where short-range forces such as general relativity (GR) become prominent. Aims: Poincare surfaces of section show that GR effects reshape phase-space structures, giving rise to a widespread and large-scale chaotic sea. This work aims to uncover the underlying mechanism of GR-enabled chaos. Methods: The dynamical structures are studied within an adiabatic framework, where an adiabatic invariant is constructed. Phase portraits, defined by the level curves of this invariant, establish analytical boundaries for the chaotic domains. Numerical simulations are subsequently performed to map orbital flips, maximum eccentricities, and Fast Lyapunov Indicators (FLI) across the initial eccentricity-inclination parameter space. Results: Phase portrait analysis demonstrates that GR precession introduces an uncertainty zone near polar inclinations, where trajectories crossing this zone inherently become chaotic. Leveraging this framework, the chaotic boundaries across the full parameter space are analytically derived, yielding excellent agreement with numerical maps of flipping orbits, maximum eccentricity, and FLI. Conclusions: In the presence of short-range forces, flipping orbits accompanied by extreme eccentricity excitation are fundamentally chaotic, and the underlying mechanism of large-scale chaos originates from the periodic passage of trajectories through the uncertainty zone in phase space over secular timescales.
Read the full abstract

Monday 21 September

Earth and Planetary Astrophysics

1I/'Oumuamua-like objects's SFDs with older, current and future surveys

Interstellar Objects (ISOs) may provide direct evidence of planetesimals ejected from early-stage solar systems. Studying these objects and predicting their detection through various surveys is vital for enhancing our understanding of planetary formation. We show how ISO size distribution (SFD), f_D~D^-q, controls expected detections in current and future surveys. We focus on TESS and the Vera Rubin Observatory LSST while analyzing Pan-STARRS and ATLAS. TESS is expected to detect ~1 ISO from shallow SFDs (q<2.5). On the other hand, LSST is likely to detect ~0.2-100 ISO yr^-1 over its lifetime. ATLAS is likely to detect one non-cometary object for populations with q<3. We further estimate ~0-8 detections of cometary like-3I/ATLAS objects by LSST in 10 years and ~1 cometary object detected by TESS over its mission time, independent of the nucleus size distribution. By considering both 1I/'Oumuamua and 3I/ATLAS from the same SFD, our results suggest that more asteroid-like objects are produced rather than comets during the early stages of planetary formation. Imminent detections of ISOs by LSST will help constrain the slopes of the SFDs for asteroid-like and cometary ISOs better, and, if applicable, differentiate between them.
Read the full abstract

Earth and Planetary Astrophysics

A Very Fast Symplectic Integrator for Planetary Systems: WHFast512 in x86 Assembly

Long-term integrations of planetary systems are used in a variety of astrophysical applications, for example to determine a system's stability. In this paper, we describe significant speed and accuracy improvements to the WHFast512 integrator. By completely rewriting WHFast512 in x86 assembly, we gain fine-grained control over CPU instructions and are able to keep the entire simulation state in CPU registers, only writing to memory when an output is required. We also reduce the number of branches in the code, improve the Kepler solver's speed and accuracy with a convergence check, switch the Hamiltonian splitting to Jacobi coordinates, implement symplectic correctors, and add support for ejection and collision detection. As a benchmark we run simulations of the Solar System with all eight planets and general relativistic corrections for 5 Gyr. With our improvements, we can integrate this system in 9 hours, more than 8x faster than the standard version of WHFast. This makes WHFast512 by far the fastest N-body integrator for this kind of simulation. We run extensive tests to verify the integrator's accuracy, with a particular focus on removing any long-term bias, and show that the accuracy is now five orders of magnitude better for typical setups compared to the earlier version of WHFast512. The new WHFast512 integrator is freely available in the REBOUND integrator package.
Read the full abstract

Earth and Planetary Astrophysics

Inside the Eirenesphere: The Interplay of Porosity, Heat Flux and Mineralogy in Exoplanetary Aquable and Habitable Layers

The classical circumstellar habitable zone restricts the search for life to planetary surfaces where stellar irradiation sustains liquid water, overlooking vast subsurface environments. Here, we present a geophysical model integrating internal radial structure, mineralogy, radiogenic heat, and pressure-dependent porosity to quantify the three-dimensional habitable volume-the eirenesphere-of rocky exoplanets. We distinguish between aquability (thermodynamic water stability) and habitability, which requires temperatures and pressures within extremophile biological limits, plus sufficient porosity for fluid circulation. To compare planetary capacities, we introduce the Eirenesphere Volumetric Index (EVI), measuring average subsurface biosphere volume. Applying this framework, we find that Earth's current state (EVI$\approx 0.06$ terrestrial oceans) represents only a moderate regime. Instead, mature super-Earths with high geothermal activity provide the most extensive environments for deep biospheres. Crustal mineralogy exerts a first-order control: thermally insulating felsic crusts sustain significantly larger eirenespheres than primary mafic lithologies. Tracking secular cooling reveals internal habitability is an evolutionary property; young planets host confined biospheres due to steep thermal gradients, whereas mature worlds maximize habitable volumes over billions of years. Finally, we demonstrate that around solar-like stars, subsurface habitability persists out to $5-7$ au, effectively decoupling life's potential from surface radiative balance. Extrapolating to the Milky Way reveals a staggering galactic capacity for subsurface life, on the order of billions of terrestrial oceans. Ultimately, this framework shifts the astrobiological paradigm from a surface-dependent phenomenon to an intrinsic planetary property, offering a metric to prioritize targets.
Read the full abstract

Earth and Planetary Astrophysics

Probing Exoplanetary Rings with Asterodensity Profiling: A PhotoRing Analysis of Kepler-51

The Kepler-51 system hosts transiting super-puff planets whose extremely low bulk densities challenge structural models, possibly because unresolved rings inflate their inferred radii. We test this hypothesis by combining asterodensity profiling with TTV-aware transit fits of Kepler long-cadence photometry. For Kepler-51b we measure a density ratio $\Psi\equiv\rho_{\star,\mathrm{obs}}/\rho_{\star,\mathrm{true}} = 1.18\pm0.06$, a $2.9\sigma$ offset from $\Psi=1$; for Kepler-51d we find $\Psi=1.04\pm0.08$ ($0.6\sigma$), consistent with no anomaly. Unmodelled eccentricity is disfavoured for planet b by TTV-derived dynamical constraints ($p=0.007$) but remains viable for planet d ($p=0.33$). The offset of Kepler-51b is also compatible with the PhotoRing (PR) effect of unmodeled rings. Using a geometric ringed-transit forward model and Bayesian retrieval, we identify ring configurations that reproduce the transit observables and the measured density ratios. The preferred solutions are consistent with compact, moderately inclined rings that raise the true bulk densities by factors of $2-5$. Conversely, the absence of a large PR offset excludes extended, optically thick, high-obliquity rings at $99\%$ credibility for b and $70\%$ for d from archival photometry alone. These objects may thus be "puff ringed planets" although recent JWST observations of Kepler-51d show no direct signature of an opaque ring. While the required impact parameters and stellar densities show tension with independent constraints -attributable to Kepler's 30-minute cadence- this study presents the first complete application of the PR effect, offered as a methodological demonstration and exclusion experiment rather than as evidence that rings are present. All code is publicly available to support applications to other systems.
Read the full abstract

All authors Jorge I. Zuluaga, Sebastian Numpaque, David M. Kipping, Jaime A. Alvarado-Montes

Earth and Planetary Astrophysics · Cross-listed from astro-ph.SR

Elemental Abundances of Cool Stars: A Spectroscopic Framework Applied to Five Planet-Host Stars

Stellar elemental abundances provide important constraints on the chemical environments from which planets form. M and late-K dwarfs are prime targets for the characterization of small exoplanets, yet their complex spectra make accurate elemental abundance measurements challenging. In this work, we present a spectral synthesis-based framework to derive detailed abundances for cool stars from high-resolution, near-infrared spectra. We apply our framework to SPIRou YJHK spectra of five planet-host stars–GJ 9827, K2-18, TOI-1201, TOI-1685, and Barnard's Star–and derive abundances of C, O, Na, Mg, Al, Si, K, Ca, Ti, Cr, Mn, and Fe for all five stars, plus Sc, V, and Ni for a subset, with typical uncertainties of $\lesssim 0.1$ dex. Our recovered metallicities and elemental abundances are broadly consistent with expected trends and previous literature analyses of these stars. We derive stellar ratios of C/O, Mg/Si, Fe/Mg, and Fe/O and find that all five systems are broadly consistent with solar ratios, with the exception of the metal-poor Barnard's Star ([M/H] $= -0.33 \pm 0.06$ dex), which is substantially iron-poor ([Fe/H] $= -0.722 \pm 0.085$ dex). Our results establish stellar chemical baselines for interior structure and atmospheric studies of these planetary systems with the Hubble Space Telescope (HST) and James Webb Space Telescope (JWST), while highlighting the remaining model-dependent limitations of precise abundance analyses for cool stars.
Read the full abstract

All authors Nicole Gromek, Drew Weisserman, Ryan Cloutier, Alexandrine L'Heureux, Charles Cadieux

Earth and Planetary Astrophysics

Feeding the Circumplanetary Disk: 3D Simulations of Dust Filtration and Accretion in PDS 70 c

How circumplanetary disks (CPDs) capture and retain solids is central to constraining the timescale for formation of rocky satellites and interpreting submillimeter continuum observations around giant planets. Here, we investigate the transport of gas and dust from the circumstellar disk into the CPD, using PDS 70 as our fiducial example. Based on observation-driven parameters, we perform high-resolution 3D adaptive mesh refinement (AMR) hydrodynamic simulations including a multifluid dust component to study dust accretion onto planets with masses of 1 $M_\mathrm{J}$ and 2.5 $M_\mathrm{J}$. We find that the pressure maximum at the gap edge imposes strong, size-dependent dust filtration, drastically lowering the solid content of the accreting flow. The net dust-to-gas mass ratio of the material accreting onto the planet is reduced by roughly two orders of magnitude relative to the outer disk. Only small grains ($\lesssim 61 \mu$m for the 1 $M_\mathrm{J}$ case and $\lesssim 10 \mu$m for the 2.5 $M_\mathrm{J}$ case) are able to accrete efficiently onto the CPD. Despite this filtering, we show that a continuous inflow of small grains can still deliver sufficient mass to build the observed PDS 70 c CPD or a Galilean-like satellite system within a few million years. Because more massive planets more effectively prevent the accretion of large grains, the dust that reaches the CPD is dominated by small particles with low millimeter-wave opacities. Consequently, in the absence of grain growth, interpreting millimeter continuum measurements of CPDs around massive giant planets may require invoking larger total dust masses than typically assumed.
Read the full abstract

All authors Charles H. Gardner, Andrea Isella, Hui Li, Shengtai Li, Gennaro D'Angelo, Adam M. Dempsey

Earth and Planetary Astrophysics

Upside down: GJ3090 b the first retrograde exoplanet around an M dwarf detected with NIRPS

The angle between stellar spin axis and planetary orbital plane can provide key insights into the formation and dynamical evolution of planetary systems. In particular, this measurement in multi-planet systems can be used to further discriminate between different competing migration scenarios. We present six transit observations of the sub-Neptune GJ3090 b obtained with NIRPS and HARPS. GJ3090 b is the inner planet of a confirmed multi-planet system orbiting an M dwarf (M2). Using high spectral resolution and high temporal cadence spectroscopic observations, we analyzed the Rossiter-McLaughlin (RM) effect induced by GJ3090 b to determine its orbital obliquity. Through the RM revolutions technique, we find that the planet is on a retrograde orbit with a derived 3D obliquity of $\psi = 136^{+24}_{-18}\,\mathrm{deg}$. We find no evidence of massive outer planetary or wide stellar binary companions, which disfavors scenarios involving gravitational perturbations from a massive body and instead points toward a primordial misalignment of the protoplanetary disk. Our results establish GJ3090 b as the first planet on a retrograde orbit discovered around an M dwarf and the first highly misaligned confirmed multi-planet system without a known massive companion. We further propose late secondary disk accretion around GJ3090, in which the disk is not expected to be aligned with the stellar spin axis, followed by disk-driven migration as the most likely mechanisms to explain the observed architecture. This work also illustrates the capability of the RM revolutions technique when applied to near-infrared data to probe orbital architectures of the smallest planets around M dwarfs that have remained mostly inaccessible.
Read the full abstract

All authors Y. Carteret, V. Bourrier, L. Parc, A. Winter, R. Allart, I. de Castro Leão, F. Bouchy, V. M. Passegger, C. Cadieux, P. Lamontagne, É. Artigau, F. Baron, S. C. C. Barros, B. Benneke, X. Bonfils, M. Bryan, B. L. Canto Martins, R. Cloutier, E. Cristo, J. Gagné, N. J. Cook, N. B. Cowan, J. R. De Medeiros, X. Delfosse, R. Doyon, X. Dumusque, D. Ehrenreich, J. I. González Hernández, D. Lafrenière, C. Lovis, L. Malo, C. Melo, L. Mignon, C. Mordasini, F. Pepe, R. Rebolo, J. Rowe, N. C. Santos, D. Ségransan, A. Suárez Mascareño, S. Udry, D. Valencia, G. Wade, K. Al Moulla, R. G. de Souza Barros de Amorim, L. Dang, E. Deibert, D. O. Fontinele, T. Forveille, Y. G. C. Frensch, R. de Lima Gomes, D. Mounzer, S. Pelletier, R. Rosener, B. N. Skinner, A. Srivastava, A. K. Stefanov, V. Vaulato, J. P. Wardenier, D. Weisserman