CyberSec Research

Understanding the diurnal behavior of lee wave clouds on Mars provides critical insight into the planet's mesoscale atmospheric dynamics and their interaction with surface topography. Lee wave clouds exhibit distinct spatial and temporal patterns that vary over the Martian day. In this study, we investigate the diurnal distribution and frequency of lee wave cloud activity during Martian Year (MY) 36 using observations from the EXI instrument aboard the Emirates Mars Mission (EMM) "Hope" spacecraft. A total of 50 lee wave events were identified, with a pronounced peak in activity during afternoon hours between solar longitudes (Ls) 270 deg and 360 deg. Our analysis reveals a seasonal and local-time dependence for these clouds, providing a comparative framework with previous mission datasets. These findings not only enhance the current understanding of Martian weather processes but also support future efforts to model and predict terrain-induced cloud dynamics across key locations on Mars.

We recently released 10 years of HARPS-N solar telescope and the goal of this manuscript is to present the different optimisations made to the data reduction, to describe data curation, and to perform some analyses that demonstrate the extreme RV precision of those data. By analysing all the HARPS-N wavelength solutions over 13 years, we bring to light instrumental systematics at the 1 m/s level. After correction, we demonstrate a peak-to-peak precision on the HARPS-N wavelength solution better than 0.75 m/s over 13 years. We then carefully curate the decade of HARPS-N re-reduced solar observations by rejecting 30% of the data affected either by clouds, bad atmospheric conditions or well-understood instrumental systematics. Finally, we correct the curated data for spurious sub-m/s RV effects caused by erroneous instrumental drift measurements and by changes in the spectral blaze function over time. After curation and correction, a total of 109,466 HARPS-N solar spectra and respective RVs over a decade are available. The median photon-noise precision of the RV data is 0.28 m/s and, on daily timescales, the median RV rms is 0.49 m/s, similar to the level imposed by stellar granulation signals. On 10-year timescales, the large RV rms of 2.95 m/s results from the RV signature of the Sun's magnetic cycle. When modelling this long-term effect using the Magnesium II activity index, we demonstrate a long-term RV precision of 0.41 m/s. We also analysed contemporaneous HARPS-N and NEID solar RVs and found the data from both instruments to be of similar quality and precision, with an overall RV differece rms of 0.79 m/s. This decade of high-cadence HARPS-N solar observations with short- and long-term precision below 1 m/s represents a crucial dataset to further understand stellar activity signals in solar-type stars , and to advance other science cases requiring such an extreme precision.

V883 Ori is an FU-Orionis-type outburst system characterized by a shoulder at 50-70 au in its ALMA band 6 and 7 intensity profiles. Previously, this feature was attributed to dust pile-up from pebble disintegration at the water snowline. However, recent multi-wavelength observations show continuity in the spectral index across the expected snowline region, disfavoring abrupt changes in grain properties. Moreover, extended water emission is detected beyond 80 au, pointing to a snowline further out. This Letter aims to explain both features with a model in which the snowline is receding. We construct a 2D disk model that solves the cooling and subsequent vapor recondensation during the post-outburst dimming phase. Our results show that both the intensity shoulder and the extended water emission are natural relics of a retreating snowline: the shoulder arises from excess surface density generated by vapor recondensation at the moving condensation front, while the outer water vapor reservoir persists due to the long recondensation timescales of $10^{2}-10^{3}$ yr at the disk atmosphere. As V883 Ori continues to fade, we predict that the intensity shoulder will migrate inward by an observationally significant amount of 10 au over about 25 years.

In binary systems with a strongly misaligned disk, the central binary stars can travel a significant vertical distance above and below the disk's orbital plane. This can cause large changes in illumination of the disk over the course of the binary orbital period. We use both analytic and radiative transfer models to examine the effect of changes in stellar illumination on the appearance of the disk, particularly in the case of the polar disk HD 98800B. We find that the observed flux from the disk can vary significantly over the binary orbital period, producing a periodically varying lightcurve which peaks twice each binary orbit. The amount of flux variation is strongly influenced by the disk geometry. We suggest that these flux variations produce several observable signatures, and that these observables may provide constraints on different properties of the disk such as its vertical structure, geometry, and cooling rate.

In many areas of astronomy, spectra of different objects are co-added or stacked to improve signal-to-noise and reveal population-level characteristics. As the number of exoplanets with measured transmission spectra grows, it becomes important to understand when stacking spectra from different exoplanets is appropriate and what stacked spectra represent physically. Stacking will be particularly valuable for long-period planets, where repeated observations of the same planet are time-consuming. Here, we show that stacked exoplanet transmission spectra are approximately mathematically equivalent to spectra generated from the geometric mean of each planet's abundance ratios. We test this by comparing stacked and geometric mean spectra across grids of forward models over JWST's NIRSpec/G395H wavelength range (2.8-5.2$\mu$m). For two dominant species (e.g., H$_2$O and CO$_2$), the geometric mean accurately reflects the stacked spectrum if abundance ratios are self-similar across planets. Introducing a third species (e.g., CH$_4$) makes temperature a critical factor, with stacking becoming inappropriate across the CO/CH$_4$ boundary. Surface gravity exerts only a minor influence when stacking within comparable planetary regimes. We further assess the number of stacked, distinct sub-Neptunes with high-metallicity atmospheres and low-pressure cloud decks required to rule out a flat spectrum at $>5\sigma$, as a function of both cloud deck pressure and per-planet spectral precision. These results provide guidance on when stacking is useful and on how to interpret stacked exoplanet spectra in the era of population studies of exoplanets.

Undifferentiated asteroids, particularly the parent bodies of carbon-rich chondrite groups, might be promising candidates for future space resource utilization due to their primitive composition and potential to host valuable metals and rare earth elements. However, our understanding of their bulk elemental composition remains limited, as most data are derived from reflectance spectra with low mineralogical resolution. Sample return missions have started to change that, as returned materials are already available to study. Still the available meteorites provide a valuable source of information about the diversity of undifferentiated asteroids in the interplanetary space. To improve compositional insights, we conducted Inductively Coupled Plasma Mass Spectrometry (ICP-MS) and ICP-AES (Inductively coupled Plasma Atomic Emission Spectroscopy) analyses on a representative suite of carbonaceous chondrites. These meteorites, considered analogs of undifferentiated asteroids, preserve materials from the early solar system and provide a geochemical record of their parent bodies. Our results highlight the abundance and distribution of transition metals, siderophile elements, and rare earth elements across several chondrite groups. These findings support the view that C-type asteroids may serve as viable sources of critical materials, while also informing future mission planning, extraction strategies, and the development of new technologies for low-gravity resource operations.

The present investigation is directed at exploring southern polar ionospheric responses to intense space weather events and their correlations with plasma convection and auroral precipitation. The main phases of six geomagnetic storms occurring in the year 2023 (ascending phase of the present solar cycle) are considered for this study. The ionospheric Total Electron Content (TEC) measurements derived from GPS receivers covering the Antarctic region are used for probing the electron density perturbations during these events. Auroral precipitation maps are shown to illustrate the locations of the GPS stations relative to particle precipitation. SuperDARN maps are shown to understand the effects of plasma convection over these locations. Correlation between the enhanced TEC observations with the auroral precipitation (R $\sim$ 0.31) and the plasma convection (R $\sim$ 0.88) reveals that the latter is more responsible for causing significant enhancements in the diurnal maximum values of TEC over the Antarctic region in comparison to the former. Therefore, this work shows correlation studies between two physical processes and ionospheric density enhancements over the under-explored south polar region under strong levels of geomagnetic activity during 2023.

Thermal emission spectra provide key insights into the atmospheric composition and especially the temperature structure of an exoplanet. With broader wavelength coverage, sensitivity and higher resolution, JWST has enabled robust constraints on these properties, including detections of photochemical products. This advances the need for retrieval frameworks capable of navigating complex parameter spaces for accurate data interpretation. In this work, we introduce the emission retrieval module of NEXOTRANS, which employs both one- and two-stream radiative transfer approximations and leverages Bayesian and machine learning techniques for retrievals. It also incorporates approximate disequilibrium chemistry models to infer photochemical species like SO2. We applied NEXOTRANS to the JWST NIRCam and MIRI emission observations of WASP-69b, covering the 2-12 microns range. The retrievals place robust constraints on the volume mixing ratios (VMR) of H2O, CO2, CO, CH4, and potential SO2. The best-fit model, i.e, free chemistry combined with non-uniform aerosol coverage, yields a log(VMR) = -3.78 (+0.15/-0.17) for H2O and -5.77 (+0.09/-0.10) for CO2 which has a sharp absorption at 4.3 micron. The second best-fit model, the hybrid equilibrium chemistry (utilizing equilibrium chemistry-grids) combined with non-uniform aerosol yields a C/O of 0.42 (+0.17/-0.13) and a metallicity of log[M/H] = 1.24 (+0.17/-0.14), corresponding to approximately 17.38 times the solar value. This hybrid chemistry retrieval also constrain SO2 with a log(VMR) = -4.85 (+0.28/-0.29), indicating possible absorption features in the 7-8 microns range. These results highlight NEXOTRANS's capability to significantly advance JWST emission spectra interpretation, offering broader insights into exoplanetary atmospheres.

Nanophase metallic iron ( $\mathrm{npFe}^0$ ) is a key indicator of space weathering on the lunar surface, primarily attributed to solar wind irradiation and micrometeoroid impacts. Recent discoveries of hematite ( $\mathrm{Fe}_2 \mathrm{O}_3$ ), a highly oxidized form of iron, in the lunar polar regions challenge the prevailing understanding of the Moon's reducing environment. This study, using ReaxFF molecular dynamics simulations of micrometeoroid impacts on fayalite ( $\mathrm{Fe}_2 \mathrm{SiO}_4$ ), investigates the atomistic mechanisms leading to both reduced and oxidized iron species. Our simulations reveals that the high-temperature and pressure conditions at the impact crater surface produces a reduced iron environment while providing a transient oxygen-rich environment in the expanding plume. Our findings bridge previously disparate observations-linking impact-driven $\mathrm{npFe}^0$ formation to the puzzling presence of oxidized iron phases on the Moon, completing the observed strong dichotomous distribution of hematite between the nearside and farside of the Moon. These findings highlight that micrometeoroid impacts, by simultaneously generating spatially distinct redox environments, provide a formation mechanism that reconciles the ubiquitous identification of nanophase metallic iron ( $\mathrm{npFe}^0$ ) in returned lunar samples with $\mathrm{Fe}^{3+}$ signatures detected by remote sensing. This underscores the dynamic nature of space weathering processes. For a more nuanced understanding of regolith evolution, we should also consider the presence of different generations or types of $\mathrm{npFe}{ }^0$, such as those formed from solar wind reduction versus impact disproportionation.

Radiative transfer calculations are essential for modeling planetary atmospheres. However, standard methods are computationally demanding and impose accuracy-speed trade-offs. High computational costs force numerical simplifications in large models (e.g., General Circulation Models) that degrade the accuracy of the simulation. Radiative transfer calculations are an ideal candidate for machine learning emulation: fundamentally, it is a well-defined physical mapping from a static atmospheric profile to the resulting fluxes, and high-fidelity training data can be created from first principles calculations. We developed a radiative transfer emulator using an encoder-only transformer neural network architecture, trained on 1D profiles representative of solar-composition hot Jupiter atmospheres. Our emulator reproduced bolometric two-stream layer fluxes with mean test set errors of ~1% compared to the traditional method and achieved speedups of 100x. Emulating radiative transfer with machine learning opens up the possibility for faster and more accurate routines within planetary atmospheric models such as GCMs.

The existence of giant extrasolar planets on short-period orbits ("hot Jupiters") represents a challenge to theories of planet formation. A leading explanation invokes perturbations from distant companions, i.e., the Eccentric Kozai-Lidov (EKL) mechanism, which can excite the eccentricities of initially wide-orbiting planets to values of order unity. The resulting tidal dissipation at periastron shrinks and circularizes the orbits to their observed configurations. While observations of orbital misalignment and distant companions support this scenario, theoretical models have struggled to reproduce the observed hot Jupiter occurrence rate. Population synthesis studies often predict that many source "cold Jupiters" are destroyed by tidal disruption during highly eccentric passages. We revisit this question with improved treatments of the mass loss and angular momentum return experienced by tidally perturbed planets. Numerical studies are performed by combining secular dynamical evolution with planetary structural evolution using Modules for Experiments in Stellar Astrophysics (MESA). We also use an analytical approach to estimate rates of tidal disruption and hot Jupiter survival. Our new population synthesis studies of giant planets in stellar binaries show that improved treatment of tidal mass loss enhances hot Jupiter survival by a factor of $\sim2-3$, yielding occurrence rates ($\gtrsim 0.5\%$ around FGK stars) consistent with observations. Angular momentum return from mass accreted onto the star may also produce a pileup of hot Jupiters near three-day orbital periods that is in statistical agreement with observations. These results suggest that EKL-driven high-eccentricity migration, when combined with realistic planetary mass loss, may be a dominant channel for hot Jupiter formation.

The present-day bulk elemental composition of an exoplanet can provide insight into a planet's formation and evolutionary history. Such information is now being measured for dozens of planets with state-of-the-art facilities using Bayesian atmosphere retrievals. We collect measurements of exoplanet composition of gas giants into a Library of Exoplanet Atmospheric Composition Measurements for comparison on a population level. We develop an open-source toolkit, ExoComp, to standardize between solar abundance, metallicity, and C/O ratio definitions. We find a systematic enhancement in the metallicity of exoplanets compared to T-dwarf and stellar populations, a strict bound in C/O between 0 and 1, and statistically significant differences between measurements from direct, eclipse, and transmission spectroscopy. In particular, the transit spectroscopy population exhibits a systematically lower C/O ratio compared to planets observed with eclipse and direct spectroscopy. While such differences may be astrophysical signals, we discuss many of the challenges and subtleties of such a comparison. We characterize the mass-metallicity trend, finding a slope consistent between planets measured in transit versus eclipse, but offset in metallicity. Compared to the Solar System and constraints from interior modeling, gas giant atmospheres appear to exhibit a steeper mass-metallicity trend. We hope that the tools available in ExoComp and the data in the Library of Exoplanet Atmospheric Composition Measurements can enhance the science return of the wide-array of space- and ground-based exoplanet science being undertaken by the community.

The properties of stars and planets are shaped by the initial conditions of their natal clouds. However, the spatial scales over which the initial conditions can exert a significant influence are not well constrained. We report the first evidence for parsec-scale spatial correlations of stellar magnetospheric inclinations ($i_{\rm mag}$), observed in the Lupus low-mass star forming region. Applying consensus clustering with a hierarchical density-based clustering algorithm, we demonstrate that the detected spatial dependencies are stable against perturbations by measurement uncertainties. The $i_{\rm mag}$ correlation scales are on the order of ~3 pc, which aligns with the typical scales of the Lupus molecular cloud filaments. Our results reveal a connection between large-scale forces -- in the form of expanding shells from the Upper Scorpius and Upper-Centaurus-Lupus regions -- and sub-au scale system configurations. We find that Lupus has a non-uniform $i_{\rm mag}$ distribution and suggest that this results from the preferential elongation of protostellar cores along filamentary axes. Non-uniformity would have significant implications for exoplanet occurrence rate calculations, so future work should explore the longevity of these biases driven by the star-cloud connection.

We present a fully integrated model of comet evolution that couples thermal and compositional processes with dynamical processes continuously, from formation to present-day activity. The combined code takes into account changes in orbital parameters that define the heliocentric distance as a function of time, which is fed into the thermal/compositional evolution code. The latter includes a set of volatile species, gas flow through the porous interior, crystallization of amorphous ice, sublimation and refreezing of volatiles in the pores. We follow the evolution of three models, with radii of 2, 10 and 50 km for 4.6 Gyr, through different dynamical epochs, starting in the vicinity of Neptune, moving to the Oort Cloud and after a long sojourn there, back inward to the planetary region. The initial composition includes a mixture of CO, CO2 ices, amorphous water ice with trapped CO and CO2, and dust.We find that the CO ice is completely depleted in the small object, but preserved in the larger ones from a depth of 500 m to the center, while the CO2 and the amorphous ice are entirely preserved. Of crucial importance is the change in CO abundance profiles during the cooling phase, as the objects migrate to the OC. Upon return from the Oort Cloud, the activity is driven by CO sublimation at large heliocentric distances (up to 50 au), by CO2 inward of 13 au and by gas released from crystallizing amorphous ice at about 7 au. We test the effect of radioactive heating by long-lived isotopes and find that it is negligible. Considering sub-solar temperatures and limited active areas, we show that CO2 production rates can exceed the detection limit as far out as 25 au.

The search for life beyond our Solar system has been a long and difficult endeavour. The majority of current efforts are focused on the potential detection of biosignatures. However, their detection and interpretation are extremely challenging. Technosignatures appear as an attractive alternative, given their expected univocal interpretation. In recent years, the number of publications discussing them have skyrocketted, both in their more rigurous and speculative sides. In this article, we explore the 28.8 years of archival radial velocity data of $\zeta^2$ Ret with the aim of detecting the proposed giant planet Calpamos, suspected source of a signal of technological origin. We performed a global model fitting the radial velocity data along with activity indicators and modelled the stellar magnetic cycle and rotation. The analysis rules out the presence of the proposed planet, as well as of any other planets more massive than 2-20 $\mathrm{M}_\oplus$ $m_{p}$ sin $i$, depending on orbital period. We show that the previously identified long-period RV signal is definitively caused by the magnetic cycle of the star.

Spectral observations of 3I/ATLAS (C/2025 N1) with JWST/NIRSpec and SPHEREx reveal an extreme CO2 enrichment (CO2/H2O = 7.6+-0.3) that is 4.5 sigma above solar system comet trends and among the highest ever recorded. This unprecedented composition, combined with substantial absolute CO levels (CO/H2O = 1.65+-0.09) and red spectral slopes, provides direct evidence for galactic cosmic ray (GCR) processing of the outer layers of the interstellar comet nucleus. Laboratory experiments demonstrate that GCR irradiation efficiently converts CO to CO2 while synthesizing organic-rich crusts, suggesting that the outer layers of 3I/ATLAS consist of irradiated material which properties are consistent with the observed composition of 3I/ATLAS coma and with its observed spectral reddening. Estimates of the erosion rate of 3I/ATLAS indicate that current outgassing samples the GCR-processed zone only (depth ~15-20 m), never reaching pristine interior material. Outgassing of pristine material after perihelion remains possible, though it is considered unlikely. This represents a paradigm shift: long-residence interstellar objects primarily reveal GCR-processed material rather than pristine material representative of their primordial formation environments. With 3I/ATLAS approaching perihelion in October 2025, immediate follow-up observations are critical to confirm this interpretation and establish GCR processing as a fundamental evolutionary pathway for interstellar objects.

The discovery of hot Jupiters has challenged the classical planet formation theory. Although various formation mechanisms have been proposed, the dominant channel and relative contributions remain unclear. Furthermore, hot Jupiters offer a unique opportunity to test tidal theory and measure the fundamental tidal quality factor, which is yet to be well-constrained. In this work, based on a hot Jupiter sample around single Sun-like stars with kinematic properties, {we find that the declining trend of their frequency is broken with a ridge at about 2 Gyr, providing direct evidence that hot Jupiters are formed with multiple origins of different timescales. By fitting with the theoretical expectations, we provide a constraint of tidal factor for Sun-like stars, which aligns well with the detected number of hot Jupiters with orbital decay. Moreover, we simultaneously constrain the relative importance of different channels: although the majority of hot Jupiters are formed early, within several tenths of Gyr via 'Early' models (e.g., in-situ formation, disk migration, planet-planet scattering and Kozai-Lidov interaction), a significant portion (about 40%) should be formed late on a relatively long timescale extending up to several Gyr mainly via the secular chaos mechanism, further supported by the obliquity distribution of 'late-arrived' hot Jupiters. Our findings provide a unified framework that reconciles hot Jupiter demographics and long-term evolution with multichannel formation.

A poleward-thinning ice shell can drive circulation in the subsurface oceans of icy moons by imposing a meridional temperature gradient--colder at the equator than the pole--through the freezing point suppression due to pressure. This temperature gradient sets a buoyancy gradient, whose sign depends on the thermal expansion coefficient determined by ocean salinity. Together with vertical mixing, this buoyancy forcing shapes key oceanic features, including zonal currents in thermal wind balance, baroclinic instability of those currents, meridional heat transport by eddies, and vertical stratification. We use high-resolution numerical simulations to explore how variations in ice shell thickness affect these processes. Our simulations span a wide range of topographic slopes, pole-to-equator temperature differences, and vertical mixing strengths, for both fresh and salty oceans. We find that baroclinic eddies dominate large-scale circulation and meridional heat transport, consistent with studies assuming a flat ice-ocean interface. However, sloped topography introduces new effects: when lighter water overlies denser water along the slope, circulation weakens as a stratified layer thickens beneath the poles. Conversely, when denser water lies beneath the poles, circulation strengthens as topography increases the available potential energy. We develop a scaling framework that predicts heat transport and stratification across all simulations. Applying this framework to Enceladus, Europa, and Titan, we infer ocean heat fluxes, stratification, and tidal energy dissipation and showing large-scale circulation constrains tidal heating and links future observations of ice thickness and rotation to subsurface ocean dynamics.