Mars finished forming while the solar nebula was still present, and acquired its primordial atmosphere from this reservoir. The absence of a detectable cometary xenon signature in the present-day Martian atmosphere suggests that the capture of solar nebular gas was significant enough to dilute later cometary contributions. By quantifying the mass of cometary material efficiently retained on Mars, we place a lower bound on the mass of the primordial Martian atmosphere. To test the robustness of our conclusions, we use cometary bombardment data from two independent studies conducted within a solar system evolutionary model consistent with its current structure. Our calculations show that, even under the most conservative scenario, the minimal mass of the primordial martian atmospheres would yield a surface pressure of no less than 2.9 bar. Such a massive nebular envelope is consistent with recent models in which atmospheric capture is strongly enhanced by the presence of heavier species on Mars - due to outgassing or redox buffering with a magma ocean.

We present new VLTI/GRAVITY astrometry and updated orbit fits for the directly imaged companions YSES 1 b and HR 2562 B, substellar objects straddling the planet-brown dwarf boundary. Using high-precision astrometry, radial velocity (RV) data, and proper motions, we derive revised orbital parameters with orbitize! arXiv:1910.01756. For YSES 1 b, the inclusion of GRAVITY astrometry and a relative radial velocity measurement from arXiv:2409.16660 overcomes the traditional challenge of constraining eccentricities for distant companions, enabling the first orbit fit and yielding a constrained eccentricity of 0.44 (0.20). This represents the first full orbit fit for the system. Additionally, we calculate a median line-of-sight stellar obliquity of 12 (+11, -8) degrees, providing further insight into the system's dynamical architecture. For HR 2562 B, our analysis agrees with arXiv:2302.04893, confirming a low-eccentricity orbit (0.34 (0.20)) and an inclination of 87 (1) degrees. We find HR 2562 B's orbit to be nearly coplanar with the debris disk, with a mutual inclination of 3.7 (0.3) degrees. For both YSES 1 b and HR 2562 B the lower eccentricities favor an in situ formation scenario over extreme scattering or cloud fragmentation.

The stellar Rossby number, a dimensionless parameter quantifying the influence of Coriolis forces on convective motions, plays a pivotal role in understanding magnetic stellar evolution. In this work, we explore the connection between the Rossby number and potential dynamo mechanisms in Sun-like stars, as well as its dependence on fundamental stellar properties. We present a novel, detailed asteroseismic calibration of the convective turnover time, incorporating for the first time Gaia photometry alongside surface gravity, effective temperature, and stellar metallicity. Our analysis employs an expanded sample of more than 150 stars, including targets from the Kepler LEGACY and KOI surveys, as well as more evolved stars observed by TESS and K2. This sample spans evolutionary stages from the main sequence to the early red giant branch (RGB), enabling a comprehensive investigation of Rossby number trends across stellar evolution.

We present and analyze follow-up, higher resolution ($R$ $\sim$ 70) $H$ and $K$ band integral field spectroscopy of the superjovian exoplanet HIP 99770 b with SCExAO/CHARIS. Our new data recover the companion at a high signal-to-noise ratio in both bandpasses and more than double the astrometric baseline for its orbital motion. Jointly modeling HIP 99770 b's position and the star's astrometry from \textit{Hipparcos} and \textit{Gaia} yields orbital parameters consistent with those from the discovery paper, albeit with smaller errors, and a slight preference for a smaller semimajor axis ($\sim$15.7--15.8 au)and a larger eccentricity ($\sim$0.28--0.29), disfavoring a circular orbit. We revise its dynamical mass slightly downwards to 15.0$_{-4.4}^{+4.5}$ $M_{\rm Jup}$ for a flat prior and 13.1$_{-5.2}^{+4.8}$ $M_{\rm Jup}$ for a more standard log-uniform mass prior, where the inclusion of its relative radial-velocity measurement is primarily responsible for these changes. \textcolor{red}{We find consistent results for HIP 99770 b's dynamical mass including recent VLTI/GRAVITY astrometry, albeit with a slightly smaller, better constrained eccentricity of $e$ $\sim$ 0.22$^{+0.10}_{-0.13}$}. HIP 99770 b is a $\sim$ 1300 K object at the L/T transition with a gravity intermediate between that of the HR 8799 planets and older, more massive field brown dwarfs with similar temperatures but with hints of equilibrium chemistry. HIP 99770 b is particularly well suited for spectroscopic follow up with Roman CGI during the technology demonstration phase at 730 nm to further constrain its metallicity and chemistry; JWST thermal infrared observations could likewise explore the planet's carbon chemistry, metallicity, and clouds.

The infall of the Large Magellanic Cloud (LMC) into the Milky Way's halo impacts the distribution of stars and dark matter in our Galaxy. Mapping the observational consequences of this encounter can inform us about the properties of both galaxies, details of their interaction, and possibly distinguish between different dark matter models. N-body simulations predict a localized overdensity trailing the LMC's orbit both in baryonic and dark matter, known as the wake. We collected wide-field, deep near-infrared, and optical photometry using VIRCAM and DECam across four fields along the expected wake, covering the sky region expected to span most of its predicted density contrast. We identify over 400 stars comprising two different tracers - near main sequence turn-off stars and red giants - that map the halo between 60-100 kpc, deriving stellar halo densities as a function of sky position and Galactocentric radius. We detect (1) a break in the halo radial density profile at 70 kpc not seen in Northern halo studies, and (2) a clear halo overdensity starting also at 70 kpc, with density contrast increasing steadily toward the expected current location of the wake. If this overdensity is the LMC wake, its peak density contrast is as pronounced as the most massive LMC model considered. Contamination from unidentified substructures may bias our wake detections, so wider-area surveys with similar depth are needed for confirmation.

Since the first detection of gravitational waves (GWs) in 2015, the International Gravitational-wave Network has made substantial strides in improving the sensitivity of ground-based detectors. Despite these advancements, many GW signals remain below the detection threshold due to environmental noise that limits sensitivity. In recent years, algorithms such as DeepClean have been developed to estimate and remove contamination from various noise sources, addressing linear, non-linear, and non-stationary coupling mechanisms. In this paper, we present noise reduction in the Virgo detector using DeepClean, serving as a preliminary step toward integrating Virgo into the online noise reduction pipeline for the O5 observing run. Our results demonstrate the applicability of DeepClean in Virgo O3b data, where noise was reconstructed from a total of 225 auxiliary witness channels. These channels were divided into 13 subsets, each corresponding to a specific frequency band, with training and subtraction performed layer-wise in a sequential manner. We observe that the subtraction improves the binary neutron star inspiral range by up to 1.3 Mpc, representing an approximately 2.5% increase. To ensure robust validation, we conduct an injection study with binary black hole waveforms. Matched-filter analyses of the injections showed an average improvement of 1.7% in the recovered signal-to-noise ratio, while parameter estimation confirmed that DeepClean introduces no bias in the recovered parameters. The successful demonstration provides a pathway for online non-linear noise subtraction in Virgo in the future observing runs.

The Laser Interferometer Space Antenna (LISA) is expected to detect thousands of individually resolved gravitational wave sources, overlapping in time and frequency, on top of unresolved astrophysical and/or primordial backgrounds. Disentangling resolved sources from backgrounds and extracting their parameters in a computationally intensive "global fit" is normally regarded as a necessary step toward reconstructing the properties of the underlying astrophysical populations. Here, we show that it is possible to infer the properties of the most numerous population of LISA sources - Galactic double white dwarfs - directly from the frequency (or, equivalently, time) strain series, by using a simulation-based approach that bypasses the global fit entirely. By training a normalizing flow on a custom-designed compression of simulated LISA frequency series from the Galactic double white dwarf population, we demonstrate how to infer the posterior distribution of population parameters (e.g., mass function, frequency, and spatial distributions). This allows for extracting information on the population parameters from both resolved and unresolved sources simultaneously and in a computationally efficient manner. Our approach to target population properties directly can be readily extended to other source classes (e.g., massive and stellar-mass black holes, extreme mass ratio inspirals), provided fast simulations are available, and to scenarios involving non-Gaussian or non-stationary noise (e.g., data gaps).

Extreme mass ratio inspirals (EMRIs) are among the most interesting gravitational wave (GW) sources for space-borne GW detectors. However, successful GW data analysis remains challenging due to many issues, ranging from the difficulty of modeling accurate waveforms, to the impractically large template bank required by the traditional matched filtering search method. In this work, we introduce a proof-of-principle approach for EMRI detection based on convolutional neural networks (CNNs). We demonstrate the performance with simulated EMRI signals buried in Gaussian noise. We show that over a wide range of physical parameters, the network is effective for EMRI systems with a signal-to-noise ratio larger than 50, and the performance is most strongly related to the signal-to-noise ratio. The method also shows good generalization ability towards different waveform models. Our study reveals the potential applicability of machine learning technology like CNNs towards more realistic EMRI data analysis.

The LIGO-Virgo-KAGRA Collaboration recently reported an exceptional gravitational-wave event, GW231123. This gravitational-wave signal was assumed to be generated from the merger of a binary black hole system, with source frame masses of $137^{+22}_{-17}~\textup{M}_\odot$ and $103^{+20}_{-52}~ \textup{M}_\odot$ (90\% credible intervals). As seen by the two LIGO detectors, the signal has only $\sim 5$ cycles, between 30 and 80 Hz, over $\sim 10$ ms. It is of critical importance to confirm the origin of this signal. Here we present the results of a Bayesian model comparison to test whether the gravitational-wave signal was actually generated by a binary black hole merger, or emitted from cusps or kinks on a cosmic string. We find significant evidence for a binary black hole merger origin of the signal.

Giant low-surface-brightness disk galaxies (gLSBGs) are rare objects with disk radii up-to 160 kpc and dynamical masses of an order of up to 10$^{12}$ $M_{\odot}$. Their very existence challenges currently accepted theories of galaxy formation and evolution, as it is difficult to build such large, dynamically cold disks through mergers without destroying them. We present deep MUSE mosaic observations of two nearby gLSBGs with compact elliptical satellites: UGC 1382, which hosts a globally counter-rotating gaseous disk, and AGC 192040, which does not. We analyze properties of ionized gas and present spatially resolved kinematics and metallicity maps; as well as stellar population analysis for the central regions of the galaxies. The radial gradients of gas-phase metallicities are flat for both galaxies. Our estimates of the effective oxygen yield suggest 'passive' gas in the outskirts of both stellar systems that is not involved in star formation. Our observational data indicate that both galaxies experienced mergers several Gyrs ago. However, the scenarios of formation of giant disks appear to be slightly different for these two systems. For AGC 192040 we propose the gas accretion from the filament followed by the intermediate-mass ratio merger with the companion on a prograde orbit. For UGC 1382 multiple gas-rich mergers with companions on retrograde orbits are preferred by the data.

The HyperScout-H (HS-H) instrument is one of the payloads aboard ESA's Hera spacecraft. Hera is a planetary defence mission that aims to provide a detailed characterization of the near-Earth binary asteroid (65803) Didymos-Dimorphos after the NASA/DART mission impact. HS-H is a versatile dual-use payload, functioning as a hyperspectral imager that captures both images and spectral data within the 0.65--0.95 $\mu$m wavelength range. The observations from this instrument will offer key insights regarding the composition of the two bodies Didymos and Dimorphos, space weathering effects, and the potential presence of exogenous material on these asteroids. Thanks to its wide field of view ($\approx 15.5^\circ \times 8.3^\circ$ in paraxial approximation), HS-H will be able to monitor the system's orbital dynamic and dust environment. At the same time, both components of this binary asteroid remain in the field of view for most of the asteroid phase of the mission. These results also complement the data obtained from other instruments in characterizing the geomorphological units. The data that will be obtained by HS-H will enable the creation of maps highlighting key spectral features, such as taxonomic classification, spectral slope, and band parameters. This article presents the pre-flight calibration of the instrument, outlines the science objectives, and discusses the expected investigations. The instrument's capabilities are demonstrated through laboratory observations of two meteorite samples and a dedicated software toolbox was developed specifically for processing the instrument's data.

Scaling up interferometry to 8m collectors should smooth-out the optical piston perturbations and allow a slow fringe tracker to obtain high precision correction on faint targets. In practice, the GRAVITY fringe tracker still observes high frequency OPD components that limit the exposure time, its precision and limiting magnitude. Perturbations seem to come from mechanical vibrations in the train of mirrors. As part of the GRAVITY+ efforts, accelerometers were added to all the mirrors of the coudé train to compensate in real-time the optical path using the main delay lines. We show their effectiveness on vibrations peaks between 40 and 200Hz and outline prospects for the upgrade of the deformable mirrors and the beam-compressor differential delay lines.

Galaxy mergers are believed to play an important role in triggering rapid supermassive black hole (SMBH) growth. As merging nuclei approach each other, the physical properties of the participating galaxies and the associated SMBH growth are expected to evolve significantly. This study measures and characterizes these physical properties throughout the merger sequence. We constructed multiwavelength Spectral Energy Distributions (SEDs) from hard X-rays to the far-infrared (FIR) for a sample of 72 nearby Active Galactic Nuclei (AGN) host galaxies. The sample comprises 64 interacting systems, including single AGNs in mergers and dual AGNs, with nuclear separations $\leq$30 kpc, as well as eight isolated active galaxies with merging features. We carefully adapted available photometric measurements at each wavelength to account for their complex morphologies and varying spatial resolutions, to perform SED fitting using CIGALE, aimed to derive critical physical properties. Our results reveal that merging galaxies hosting AGN(s) show deviations from the star-forming main sequence, and a wide range of star formation rates (SFRs). Both AGN activity and star formation are significantly influenced by the merger process, but these effects are more prominent in major, mass ratios <4:1, interactions. We find that the projected nuclear separation is not a good tracer of the merger stage. Instead, morphological classification accurately assesses the merger progression. Based on this morphological analysis, late-stage mergers exhibit elevated SFRs (5.1$\times$), AGN luminosities (2.4$\times$), and nuclear obscuration (2.8$\times$) compared to earlier stages, supporting previous findings and reinforcing the link between merger-driven galaxy evolution and SMBH growth.

Flexible and efficient noise characterization is crucial for the precise estimation of gravitational wave parameters. We introduce a fast and accurate Bayesian method for estimating the power spectral density (PSD) of long, stationary time series tailored specifically for LISA data analysis. Our approach models the PSD as a geometric mean of a parametric and a nonparametric component, combining the computational efficiency of parametric models with the flexibility to capture deviations from theoretical expectations. The nonparametric component is expressed by a mixture of penalized B-splines. Adaptive, data-driven knot placement performed once during initialization eliminates computationally expensive reversible-jump Markov Chain Monte Carlo, while hierarchical roughness penalty priors prevent overfitting. This design yields stable, flexible PSD estimates with runtimes of minutes instead of hours. Validation on simulated autoregressive AR(4) data demonstrates estimator consistency. It shows that well-matched parametric components reduce the integrated absolute error compared to an uninformative baseline, requiring fewer spline knots to achieve comparable accuracy. Applied to a year of simulated LISA $X$-channel noise, our method achieves relative integrated absolute errors of $\mathcal{O}(10^{-2})$ with computation times less than three minutes, which makes it suitable for iterative analysis pipelines and multi-year mission datasets.

To date, galaxy image simulations for weak lensing surveys usually approximate the light profiles of all galaxies as a single or double Sérsic profile, neglecting the influence of galaxy substructures and morphologies deviating from such a simplified parametric characterization. While this approximation may be sufficient for previous data sets, the stringent cosmic shear calibration requirements and the high quality of the data in the upcoming Euclid survey demand a consideration of the effects that realistic galaxy substructures have on shear measurement biases. Here we present a novel deep learning-based method to create such simulated galaxies directly from HST data. We first build and validate a convolutional neural network based on the wavelet scattering transform to learn noise-free representations independent of the point-spread function of HST galaxy images that can be injected into simulations of images from Euclid's optical instrument VIS without introducing noise correlations during PSF convolution or shearing. Then, we demonstrate the generation of new galaxy images by sampling from the model randomly and conditionally. Next, we quantify the cosmic shear bias from complex galaxy shapes in Euclid-like simulations by comparing the shear measurement biases between a sample of model objects and their best-fit double-Sérsic counterparts. Using the KSB shape measurement algorithm, we find a multiplicative bias difference between these branches with realistic morphologies and parametric profiles on the order of $6.9\times 10^{-3}$ for a realistic magnitude-Sérsic index distribution. Moreover, we find clear detection bias differences between full image scenes simulated with parametric and realistic galaxies, leading to a bias difference of $4.0\times 10^{-3}$ independent of the shape measurement method. This makes it relevant for stage IV weak lensing surveys such as Euclid.

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

Direct numerical simulations are used to study the interaction of a stream of small heavy inertial particles with the laminar and turbulent wakes of an immobile sphere facing an incompressible uniform inflow. Particles that do not collide with the obstacle but move past it, are found to form preferential concentrations both in the sphere boundary layer and in its wake. In the laminar case, the upstream diverging flow pattern is responsible for particle clustering on a cylinder that extends far downstream the sphere. The interior of this surface contains no particles and can be seen as a shadow of the large obstacle. Such concentration profiles are also present in the case of turbulent wakes but show a finite extension. The sphere shadow is followed by a region around the axis of symmetry where the concentration is higher than the average. It originates from a resonant centrifugal expulsion of particles from shed vortices. The consequence of this concentration mechanism on monodisperse inter-particle collisions is also briefly discussed. They are enhanced by both the increased concentration and the presence of large velocity differences between particles in the wake.

ASTEP (Antarctica Search for Transiting ExoPlanets) is a pilot project that aims at searching and characterizing transiting exoplanets from Dome C in Antarctica and to qualify this site for photometry in the visible. Two instruments were installed at Dome C and ran for six winters in total. The analysis of the collected data is nearly complete. We present the operation of the instruments, and the technical challenges, limitations, and possible solutions in light of the data quality. The instruments performed continuous observations during the winters. Human interventions are required mainly for regular inspection and ice dust removal. A defrosting system is efficient at preventing and removing ice on the mirrors. The PSF FWHM is 4.5 arcsec on average which is 2.5 times larger than the specification, and is highly variable; the causes are the poor ground-level seeing, the turbulent plumes generated by the heating system, and to a lower extent the imperfect optical alignment and focusing, and some astigmatism. We propose solutions for each of these aspects that would largely increase the PSF stability. The astrometric and guiding precisions are satisfactory and would deserve only minor improvements. Major issues are encountered with the camera shutter which did not close properly after two winters; we minimized this issue by heating the shutter and by developing specific image calibration algorithms. Finally, we summarize the site testing and science results obtained with ASTEP. Overall, the ASTEP experiment will serve as a basis to design and operate future optical and near-infrared telescopes in Antarctica.