Papers for Wednesday, Oct 29 2025

Most stars form in dense stellar environments, where frequent close encounters can strongly perturb and reshape the early architecture of planetary systems. The solar system, with its rich population of distant comets, provides a natural laboratory to study these processes. We perform detailed numerical simulations using the LonelyPlanets framework that combines NBODY6++GPU and REBOUND, to explore the evolution of debris disks around solar system analogues embedded in stellar clusters. Two initial configurations are considered, an $Extended$ and a $Compact$ model, each containing four giant planets and either an extended or compact debris disk. We find that compact disks primarily form Kuiper belt and scattered disk-like populations through planet-disk interactions, while extended disks are more strongly shaped by stellar encounters, producing Oort cloud-like structures and interstellar comets with ejection velocities of 1-3 km/s. Stellar perturbations are most effective for encounter inclinations between $0^{\circ}$ and $30^{\circ}$, giving rise to distinct dynamical populations, like Sednoids, and inner Oort cloud analogues, and a characteristic tail in semi-major axis-eccentricity space. In coplanar encounters, the disk remains largely flattened, whereas polar flybys redistribute angular momentum vertically, producing nearly isotropic outer populations that resemble an emerging Oort cloud. Our results suggest that cometary reservoirs and interstellar objects are natural byproducts of planet-disk interactions and stellar flybys in dense clusters, linking the architecture of outer planetary systems to their birth environments.

The evolution of gaseous planets is a complex process influenced by various physical parameters and processes. In this study, we present critical modifications to the Modules for Experiments in Stellar Astrophysics (MESA) code to enhance its applicability to giant planet modelling. We introduce an equation of state specifically tailored for materials at planetary conditions. The equation of state considers the thermodynamic properties of hydrogen-helium mixtures and heavy elements, improving the accuracy of internal structure calculations. We also present modifications to the radiative opacity to allow the modelling of grains, clouds and opacity windows. Furthermore, we refine the treatment of convective mixing processes in MESA to better replicate convective mixing with the presence of composition gradients. Finally, we add a treatment for helium rain and settling. These modifications aim to enhance the predictive capabilities of MESA for giant planet evolution and are publicly available. We hope that these improvements will lead to a deeper understanding of giant planet evolution in the solar system and beyond.

Regolith simulants are essential for space research and technology development. Yet, their physical properties often differ from those of true planetary soil, particularly when compared to regolith properties in-situ, that experience notably reduced gravity. We focus on lunar regolith simulants and explore various techniques to modify existing simulants to replicate the mechanical/rheological behavior of Moon regolith in true lunar conditions. Our results are validated against data from in-situ tests conducted during the Luna and Apollo missions, enabling a direct comparison of physical properties of our enhanced simulants to true Lunar regolith, on the Moon. Analyzed in the Mohr-Coulomb model framework, the friction angle of most simulants is generally close to that of true regolith, but the measured cohesion is often higher on the Moon, notably due to the reduced gravity environment. We propose a method to increase the cohesion of an existing simulant and assess the mechanical behavior of our rheological regolith simulant using a standardized geotechnical, shear test. The experimental results are then directly compared to in-situ data, providing a quantitative basis for evaluating the fidelity of the enhanced simulants.

The nominal duration of Euclid's main mission is six years, but current best estimates indicate that the observatory has sufficient propellant to operate for up to ~14 years in total. In this work, we advocate dedicating six of these ~8 additional years to repeating the main survey, covering approximately one-third of the sky. This repetition would not only improve the sampling, signal-to-noise, quality, and depth of the survey, but -- most importantly -- would provide a six-year time baseline between two epochs if executed in the same sequence. The availability of multiple epochs would enable the derivation of proper motions for stars as faint as V~27, i.e., more than five magnitudes fainter than those measured by the Gaia mission. Although it may seem early to propose such a mission extension, in this work we quantitatively illustrate its immense scientific potential. We therefore intend to initiate the technical and scientific discussions early to ensure optimal planning. The here proposed extension would employ only the VIS channel -- owing to its superior astrometric capability and depth -- while simultaneously using NISP in slitless-spectroscopy mode to enhance the signal-to-noise ratio of first-epoch spectra that would also benefit of proper motions to identify and reject objects within the local Universe.

Merging galaxy clusters exhibit strong observational evidence for efficient particle acceleration in the intracluster medium (ICM), particularly in the form of synchrotron-emitting radio relics and halos. Cosmic ray (CR) electrons are likely accelerated (or re-accelerated) at merger and accretion shocks via diffusive shock acceleration (DSA). However, in the presence of the large diffusion coefficients one would naively expect in the rarefied, relatively unmagnetized ICM, this acceleration--in particular, the maximum proton energy ($E_{\rm max}$)--is limited by long acceleration times. On the other hand, recent work on CR transport suggests that the diffusion coefficient can be suppressed in ICM-like environments. In this picture, deviations from local thermodynamic equilibrium can trigger the mirror instability, creating plasma-scale magnetic structures, or "micromirrors," that efficiently scatter CRs. In this paper, we investigate the implications of micromirror confinement for shock acceleration in the ICM. We demonstrate that micromirrors enforce a minimum value of $E_{\rm max} \gtrsim 100$ GeV that does not rely on CR-driven magnetic field amplification. We also discuss micromirror confinement in the context of cosmological simulations and $\gamma$-ray observations, and present a simulation of a Coma-like merging cluster that self-consistently includes CR acceleration at shocks, with an effective diffusion coefficient set by micromirrors. We show that the introduction of micromirrors yields simulated galaxy clusters that remain consistent with $\gamma$-ray observations.

Distance-based methods involve the computation of distance values between features and are a well-established paradigm in machine learning. In anomaly detection, anomalies are identified by their large distance from normal data points. However, the performance of these methods often hinges on a single, user-selected distance metric (e.g., Euclidean), which may not be optimal for the complex, high-dimensional feature spaces common in astronomy. Here, we introduce a novel anomaly detection method, Distance Multi-Metric Anomaly Detection (DiMMAD), which uses an ensemble of distance metrics to find novelties. Using multiple distance metrics is effectively equivalent to using different geometries in the feature space. By using a robust ensemble of diverse distance metrics, we overcome the metric-selection problem, creating an anomaly score that is not reliant on any single definition of distance. We demonstrate this multi-metric approach as a tool for simple, interpretable scientific discovery on astronomical time series -- (1) with simulated data for the upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time, and (2) real data from the Zwicky Transient Facility. We find that DiMMAD excels at out-of-distribution anomaly detection -- anomalies in the data that might be new classes -- and beats other state-of-the-art methods in the goal of maximizing the diversity of new classes discovered. For rare in-distribution anomaly detection, DiMMAD performs similarly to other methods, but may allow for improved interpretability. All our code is open source: DiMMAD is implemented within DistClassiPy: this https URL, while all code to reproduce the results of this paper is available here: this https URL.

Self-gravitating galactic halos composed of self-interacting dark matter exhibit the formation of a highly dense core at the galactic center--a gravothermal collapse. Analytic models to describe this evolution have been developed and calibrated to numerical simulations initialized with isotropic particle velocity distributions, an assumption not necessarily warranted by the theory of halo formation. Here we study the dependence of the timescale for gravothermal collapse on the velocity distribution. To do so, we consider self-consistent initial conditions for halos with the same density distribution but with different velocity distributions. We consider models with constant anisotropy and with an anisotropy that increases with radius. The velocity distributions that we explore have collapse times that differ from that assuming isotropic distributions by more than a factor of two. We argue that these variations may depend on the global changes in velocity-dispersion profiles in these anisotropic models and not just on the degree of anisotropy.

We present a catalogue of compact sources detected in the SARAO MeerKAT 1.3 GHz Galactic Plane Survey (SMGPS). We extract 510599 compact sources, with areas less than five 8" beams, from the survey maps covering the regions $252°< l < 358^°$ and $2°< l < 61°$ at $|b| \leq 1.5^°$, which have an angular resolution of 8" and a sensitivity of $\sim$ 10-30 $\mu$Jy beam$^{-1}$. In this paper, we describe the source identification and characterisation methods, present the quality assurance of the catalogue, explore the nature of the catalogue sources, and present initial science highlights. We limit our catalogue to sources with a signal-to-noise ratio $\geq 5$, as the catalogue is $\sim$90 per cent complete, and has a false positive rate of less than 1 per cent at this threshold. The bulk of the catalogue sources are previously unknown to the literature, with the majority of unknown sources at sub-mJy levels. Initial science highlights from the catalogue include the detection of 213 radio quiet WISE HII region candidates, previously undetected in radio continuum studies. We show images that compare the SMGPS compact sources to CORNISH ultracompact HII regions, thus highlighting the sensitivity and unprecedented uv-coverage of the SMGPS, and the potential synergy of the SMGPS with other surveys.

Young stars host only a small fraction of the known exoplanet population because their photometric variability, magnetic activity, and frequent placement in dense, poorly-resolved regions hamper exoplanet detections. Yet, measuring planets at these ages is crucial since these phases are when dynamical processes that drive planetary migration are most active. We assess the expected yield of a hypothetical Nancy Grace Roman Space Telescope transit survey of the Rosette Nebula, a $\sim10\,\mathrm{Myr}$ star-forming region with a dense and diverse stellar population. Using the Roman Exposure Time Calculator to quantify sensitivity to Rosette members, we establish detection thresholds for companions and evaluate yields via Monte Carlo injection-recovery simulations. We predict the detection of $33 \pm 9$ young transiting exoplanets orbiting stellar hosts in a month-long survey, and $29 \pm 8$ in a two-week survey. The extended baseline primarily improves sensitivity to longer-period planets orbiting FGK stars, while most M dwarf detections are well-sampled within two weeks. Irrespective of the temporal baseline, transit detections are dominated by of 1-2\,$R_\oplus$ super-Earths and sub-Neptunes with $P\lesssim8$\,days. Such a sample would substantially expand the census of only three detected planets younger than $20\,\mathrm{Myr}$, probing an age regime in which planetary radii remain inflated, the stability of close-in orbits is uncertain, and planetary migration may still be ongoing. This survey offers a path to constrain early planetary evolution and establish prime follow-up targets for the James Webb Space Telescope, Vera Rubin Observatory, and the Habitable Worlds Observatory.

On August 18th, 2025, the LIGO--Virgo--KAGRA collaboration reported a sub-threshold gravitational wave candidate detection consistent with a sub-solar-mass neutron star merger, denoted S250818k. An optical transient, AT2025ulz, was discovered within the localization region. AT2025ulz initially appeared to meet the expected behavior of kilonova (KN) emission, the telltale signature of a binary neutron star merger. The transient subsequently rebrightened after $\sim$\,$5$ days and developed spectral features characteristic of a Type IIb supernova. In this work, we analyze the observations of the host galaxy of AT2025ulz obtained by the Dark Energy Spectroscopic Instrument (DESI). From the DESI spectrum, we obtain a secure redshift of $z = 0.084840 \pm 0.000006$, which places the transient within $2\sigma$ of the gravitational wave distance and results in an integral overlap between the gravitational wave alert and the transient location of $\log_{10}\mathcal{I} \approx 3.9-4.2$. Our analysis of the host galaxy's spectral energy distribution reveals a star-forming, dusty galaxy with stellar mass ${\sim} 10^{10}~M_\odot$, broadly consistent with the population of both short gamma-ray bursts and core-collapse supernova host galaxies. We also present our follow-up of DESI-selected candidate host galaxies using the Fraunhofer Telescope at the Wendelstein Observatory, and show the promise of DESI for associating or rejecting candidate electromagnetic counterparts to gravitational wave alerts. These results emphasize the value of DESI's extensive spectroscopic dataset in rapidly characterizing host galaxies, enabling spectroscopic host subtraction, and guiding targeted follow-up.

White dwarfs (WDs) are frequently observed to have strong magnetic fields up to $10^9$ G and expected to have a possible internal field as high as $\sim 10^{14}$ G. High internal fields can significantly deform a WD's equilibrium structure, generating a quadrupole moment. If the rotation axis is misaligned with the magnetic axis, the deformation can lead to the emission of continuous gravitational waves (CGWs). We examine the potential for detecting CGWs from magnetized WDs with future space-based detectors such as LISA, ALIA, DECIGO, Deci-Hz, BBO and TianQin. We model the field-induced deformation and compute the resulting GW strain, incorporating amplitude decay due to angular momentum loss from electromagnetic and gravitational radiation. This sets a timescale for detection -`active timescale' of $10^{5-6}$ yr, requiring observation while the object remains sufficiently young. Our results suggest that LISA could detect a few dozens of highly magnetized WDs across the Galaxy during its mission. As a specific case, we investigate ZTF J1901+1458- a compact, massive, fast-rotating, and strongly magnetized WD with spin period $\sim416$ s and inferred surface field $\sim10^{9}$ G. We find that this object would be detectable by LISA with four years of continuous data. This highlights the potential of CGW observations to probe magnetic field structure in WDs and their role in type Ia supernova progenitors.

The first joint electromagentic (EM) and gravitational wave (GW) detection, known as GW170817, marked a critical juncture in our collective understanding of compact object mergers. However, it has now been 8 years since this discovery, and the search for a second EM-GW detection has yielded no robust discoveries. Recently, on August 18, 2025, the LIGO-Virgo-KAGRA collaboration reported a low-significance (high false alarm rate) binary neutron star merger candidate S250818k. Rapid optical follow-up revealed a single optical candidate AT2025ulz ($z=0.08484$) that initially appeared consistent with kilonova emission. We quickly initiated a set of observations with \textit{Swift}, \textit{XMM-Newton}, \textit{Chandra}, and the Very Large Array to search for non-thermal afterglow emission. Our deep X-ray and radio search rules out that the optical rebrightening of AT2025ulz is related to the afterglow onset, reinforcing its classification as a stripped-envelope supernova (SN 2025ulz). We derive constraints on the afterglow parameters for a hypothetical binary neutron star merger at the distance of AT2025ulz ($\approx 400$ Mpc) based on our X-ray and radio limits. We conclude that our observational campaign could exclude a GW170817-like afterglow out to viewing angles of $\theta_\textrm{v}\approx 12.5$ degrees. We briefly discuss the prospects for the future discovery of off-axis afterglows.

We present hydrodynamic simulations with the moving-mesh code AREPO of Tidal Disruption Encores (TDEEs) in nuclear star clusters (NSCs). TDEEs arise when a stellar-mass black hole (sBH) disrupts a star within the NSC, producing debris that is unbound from the sBH but remains gravitationally bound to the central massive black hole (MBH), leading to a delayed secondary flare. We find that the morphology and thermodynamics of the fallback material depend sensitively on the disruption geometry, MBH mass, and sBH-MBH separation. We identify two distinct morphological outcomes: ring encores, where debris circularize into a torus, and direct encores, where streams plunge toward the MBH, with encore luminosities peaking at times corresponding to the freefall timescale and one orbital period, respectively. Across all simulated cases, we find these events exhibit luminosities of $10^{40}-10^{42}$ erg/s with lightcurves characteristic of their morphology. Our work greatly improves the predictions of TDEE lightcurves and empowers observations to probe into NSC dynamics and sBH population while providing possible explanations for anomalous TDE-like flares.

On August 18, 2025, the LIGO-Virgo-KAGRA collaboration reported gravitational waves from a sub-threshold binary neutron star merger. If astrophysical, this event would have a surprisingly low chirp mass, suggesting that at least one neutron star was below a solar mass. The Zwicky Transient Facility mapped the coarse localization and discovered a transient, ZTF25abjmnps (AT2025ulz), that was spatially and temporally coincident with the gravitational wave trigger. The first week of follow-up suggested properties reminiscent of a GW170817-like kilonova. Subsequent follow-up suggests properties most similar to a young, stripped-envelope, Type IIb supernova. Although we cannot statistically rule out chance coincidence, we undertake due diligence analysis to explore the possible association between ZTF25abjmnps and S250818k. Theoretical models have been proposed wherein sub-solar neutron star(s) may form (and subsequently merge) via accretion disk fragmentation or core fission inside a core-collapse supernova i.e. a ``superkilonova". Here, we qualitatively discuss our multi-wavelength dataset in the context of the superkilonova picture. Future higher significance gravitational wave detections of sub-solar neutron star mergers with extensive electromagnetic follow-up would conclusively resolve this tantalizing multi-messenger association.

Magrathea is an open-source C++ code for modeling the internal structure of differentiated planets. The initial release, Huang et al. (2022), introduced the base solver, a modular framework for defining equations of state (EOSs) used within phase diagrams for each differentiated layer, and outlined a series of planned extensions. Many of those updates are now implemented. Magrathea v2 is a more versatile platform that supports a winder range of compositions, adds new tools for composition retrieval, and makes it easier for users to adapt the code to their own models.

The tidal torque theory (TTT) predicts that galaxy spins are correlated with the surrounding tidal field, reflecting how angular momentum is acquired during structure formation. We present a new observational test of this prediction using the final data release of the SDSS-IV MaNGA integral field spectroscopy survey, which enables direct spin measurements from stellar and ionized gas kinematics for nearly 10,000 galaxies. Unlike previous studies that relied on filament orientations, we directly utilize the three-dimensional tidal field reconstructed from the galaxy distribution, providing a physically defined reference frame for the analysis. We find that massive galaxies tend to align their spins parallel to the intermediate axis of the tidal field, consistent with the prediction of the TTT. They also tend to orient perpendicular to the major axis, hinting at the possible influence of baryonic processes. In contrast, low-mass galaxies exhibit the opposite trend, with a transition mass of $M_*\sim10^{10}-10^{10.5}\,M_\odot$. No significant alignment is detected with respect to the minor axis across all stellar masses. We further examine the dependence on morphology and environment, finding that S0 and early-type spiral galaxies exhibit stronger alignment signals than late-type spirals. The alignment trend becomes particularly pronounced in regions of high tidal anisotropy and high overdensity. A mutual information analysis identifies these environmental factors as the dominant drivers of the observed trends. Our results provide the most direct empirical evidence to date for the connection between galaxy spins and the cosmic tidal field.

We present high-resolution ($\sim$50$-$100 pc) Atacama Large Millimeter Array (ALMA) observations of $^{12}$CO(2-1) or $^{12}$CO(1-0) emission in seven local ($z$ $\lesssim$ 0.05) major mergers -- five of which are dual active galactic nuclei (AGN) systems, and two of which are single AGN systems. We model the molecular gas kinematics through rotating disk profiles using a Bayesian Markov chain Monte Carlo approach. The residuals were then used to isolate non-rotating components of the molecular gas-- the most likely contributor to future SMBH growth. We find that more massive SMBHs have higher surface densities of non-rotating molecular gas within their sphere of influence. This potential molecular gas supply, however, does not correlate with the current accretion efficiency of the SMBHs, suggesting that only a fraction of the observed non-rotating gas is currently reaching the SMBH. Finally, we tentatively find no significant differences in the nuclear molecular gas masses of single AGN and dual AGN hosts, both within the SMBH sphere of influence and within the central kiloparsec. Our results indicate that the probability of occurrence of the dual AGN phenomenon is likely dependent on AGN variability and/or obscuration rather than the availability of molecular gas in the nuclear regions.

We present the first detection of GeV $\gamma$-ray emission potentially associated with the pulsar wind nebula (PWN) hosted by the young core-collapse supernova remnant G292.0+1.8, based on a detailed time-resolved analysis of \textit{Fermi}-LAT data. By isolating the unpulsed component from the dominant magnetospheric radiation of PSR~J1124$-$5916, we successfully disentangle a candidate nebular emission in the GeV range, characterise its morphology and extract its spectrum. This identification places G292.0+1.8 among the few systems in which the pulsar and PWN contributions have been spectrally resolved at high energies, offering new insight into their respective emission mechanisms. We characterise the $\gamma$-ray spectrum of the pulsar and model the broadband spectral energy distribution (SED) of the PWN using radio, X-ray, and GeV data. The emission is well described by a single electron population with two spectral breaks: one intrinsic to the injection spectrum and another produced by synchrotron cooling in a magnetic field of $\sim$15~$\mu$G. Notably, the inferred magnetic field and the low TeV flux of the nebula resemble those of 3C~58, suggesting that similar low-field environments can arise in young PWNe. The high-energy portion of the SED is now tightly constrained by our GeV detection and existing TeV upper limits. Compared to our model, earlier predictions tend to underpredict the $\gamma$-ray flux, while others that succeed in reproducing the GeV component often overpredict the TeV emission. This mismatch underscores the challenges in modelling particle acceleration and radiation processes in young PWNe and establishes G292.0+1.8 as a valuable benchmark for testing and refining such models.

Sparse Autoencoders (SAEs) can efficiently identify candidate monosemantic features from pretrained neural networks for galaxy morphology. We demonstrate this on Euclid Q1 images using both supervised (Zoobot) and new self-supervised (MAE) models. Our publicly released MAE achieves superhuman image reconstruction performance. While a Principal Component Analysis (PCA) on the supervised model primarily identifies features already aligned with the Galaxy Zoo decision tree, SAEs can identify interpretable features outside of this framework. SAE features also show stronger alignment than PCA with Galaxy Zoo labels. Although challenges in interpretability remain, SAEs provide a powerful engine for discovering astrophysical phenomena beyond the confines of human-defined classification.

The Equation of State (EOS) of matter within neutron stars is a central topic in nuclear physics and astrophysics. A precise understanding of the composition and phase behavior of matter under such extreme conditions is crucial for uncovering the fundamental laws of the strong interaction. This study investigates hadron-quark hybrid stars using a two-flavor Nambu-Jona-Lasinio (NJL) model. As an effective theory, this model can describe the generation of dynamical quark masses and chiral symmetry restoration characteristic of dense quark matter. We construct the hybrid EOS by joining the BSR6 relativistic mean-field model for hadronic matter with the NJL model for quark matter. A quintic polynomial interpolation ensures a smooth ($C^2$ continuity) and thermodynamically consistent crossover between the phases. Based on this hybrid EOS, we solve the Tolman-Oppenheimer-Volkoff (TOV) equations to calculate macroscopic properties of neutron stars, such as the mass-radius ($M-R$) relationship and the tidal deformability parameter ($\Lambda$). By exploring key model parameters, we identify a region satisfying a wide range of multi-messenger constraints. Our resulting EOS supports a maximum mass consistent with PSR J0740+6620, while simultaneously predicting radii and tidal deformabilities for a $1.4M_{\odot}$ star that agree with NICER observations and limits from GW170817. This work thus presents a self-consistent model that resolves the tension between high-mass pulsars and small tidal deformabilities, deepening our understanding of the hadron-quark crossover.

We present the first systematic study of multi-domain map-to-map translation in galaxy formation simulations, leveraging deep generative models to predict diverse galactic properties. Using high-resolution magneto-hydrodynamical simulation data, we compare conditional generative adversarial networks and diffusion models under unified preprocessing and evaluation, optimizing architectures and attention mechanisms for physical fidelity on galactic scales. Our approach jointly addresses seven astrophysical domains - including dark matter, gas, neutral hydrogen, stellar mass, temperature, and magnetic field strength - while introducing physics-aware evaluation metrics that quantify structural realism beyond standard computer vision measures. We demonstrate that translation difficulty correlates with physical coupling, achieving near-perfect fidelity for mappings from gas to dark matter and mappings involving astro-chemical components such as total gas to HI content, while identifying fundamental challenges in weakly constrained tasks such as gas to stellar mass mappings. Our results establish GAN-based models as competitive counterparts to state-of-the-art diffusion approaches at a fraction of the computational cost (in training and inference), paving the way for scalable, physics-aware generative frameworks for forward modelling and observational reconstruction in the SKA era.

The cosmic dipole anomaly -- the observation of a significant mismatch between the dipole observed in the matter distribution and that expected given the kinematic interpretation of the cosmic microwave background dipole -- poses a serious challenge to the Cosmological Principle upon which the standard model of cosmology rests. Measurements of the dipole ($\ell=1$) in a given sample crucially depend on having control over other large-scale power ($\ell > 1$) so as to avoid biases, in particular those potentially caused by correlations among multipoles during fitting, and those by local source clustering. Currently, the most powerful catalogue that exhibits the cosmic dipole anomaly is the sample of 1.6~million mid-infrared quasars derived from CatWISE2020. We therefore analyse clustering properties of this catalogue by performing an inference analysis of large-scale multipoles in real space, and by computing its angular power spectrum on small scales to test for convergence with $\Lambda$CDM. After accounting for the known trend of the quasar number counts with ecliptic latitude, we find that any other large-scale power is consistent with noise, find no evidence for the presence of an octupole ($\ell=3$) in the data, and quantify the clustering dipole's proportion to be marginal. Our results therefore reaffirm the anomalously high dipole in the distribution of quasars.

The unprecedented volume and quality of data from space- and ground-based telescopes present an opportunity for machine learning to identify new classes of variable stars and peculiar systems that may have been overlooked by traditional methods. Extending prior methodological work, this study investigates the potential of an unsupervised learning approach to scale effectively to larger stellar populations, including objects in crowded fields, and without the need for pre-selected catalogues, specifically focusing on 13 405 sources selected from Gaia DR3 and lying in the selected region of the CMD. Our methodology incorporates unsupervised clustering techniques based primarily on statistical features extracted from Gaia DR3 epoch photometry. We used the t-distributed stochastic neighbour embedding (t-SNE) algorithm to identify variability classes, their subtypes, and spurious variability induced by instrumental effects. The clustering results revealed distinct groups, including hot subdwarfs, cataclysmic variables (CVs), eclipsing binaries, and objects in crowded fields, such as those in the Andromeda (M31) field. Several potential stellar subtypes also emerged within these clusters. Notably, objects previously labelled as RR Lyrae were found in an unexpected region of the CMD, potentially due to either unreliable astrometric measurements (e.g., due to binarity) or alternative evolutionary pathways. This study emphasises the robustness of the proposed method in finding variable objects in a large region of the Gaia CMD, including variable hot subdwarfs and CVs, while demonstrating its efficiency in detecting variability in extended stellar populations. The proposed unsupervised learning framework demonstrates scalability to large datasets and yields promising results in identifying stellar subclasses.

We present a systematic comparison of three independent machine learning (ML)-based searches for strong gravitational lenses applied to the Dark Energy Survey (Jacobs et al. 2019a,b; Rojas et al. 2022; Gonzalez et al. 2025). Each search employs a distinct ML architecture and training strategy, allowing us to evaluate their relative performance, completeness, and complementarity. Using a visually inspected sample of 1651 systems previously reported as lens candidates, we assess how each model scores these systems and quantify their agreement with expert classifications. The three models show progressive improvement in performance, with F1-scores of 0.31, 0.35, and 0.54 for Jacobs, Rojas, and Gonzalez, respectively. Their completeness for moderate- to high-confidence lens candidates follows a similar trend (31%, 52%, and 70%). When combined, the models recover 82% of all such systems, highlighting their strong complementarity. Additionally, we explore ensemble strategies: average, median, linear regression, decision trees, random forests, and an Independent Bayesian method. We find that all but averaging achieve higher maximum F1 scores than the best individual model, with some ensemble methods improving precision by up to a factor of six. These results demonstrate that combining multiple, diverse ML classifiers can substantially improve the completeness of lens samples while drastically reducing false positives, offering practical guidance for optimizing future ML-based strong lens searches in wide-field surveys.

Exoplanetary systems exhibit a broad range of architectures which, in turn, enable a variety of dynamical environments. Many of the known planetary systems do not transit the host star, and so we measure the minimum masses of their planets, making it difficult to fully assess the dynamical environment within the system. Astrometry can resolve the mass ambiguity and thus allow a more complete dynamical analysis of systems to be conducted. Gaia-4 and Gaia-5 are two such systems, whose study with radial velocities and data from the Gaia mission revealed that each star harbors a massive planet on a highly eccentric orbit. In this work, we provide the results of a dynamical analysis of each system, including calculations of the Habitable Zone (HZ), from which we show that the presence of the known companions largely exclude the presence of planets within the HZ. We discuss the diagnostics of potential past planet-planet scattering events, and the occurrence of similar systems whereby a giant planet on an eccentric orbit can substantially disrupt orbital integrity of terrestrial planets. These "wrecking ball" systems have an impact on the target selection for planned direct imaging missions that seek to identify potentially habitable environments.

In this two-paper series, we present a straightforward mathematical model for synthesizing quasar absorption line profiles from sight lines through idealized, spatial-kinematic models of the circumgalactic medium (CGM) and their host galaxies. In Paper I, we developed the spatial components of the galaxy/CGM structures (disk, halo, wind, accretion) and their 3D velocity fields. We derived the formalism for arbitrary observed orientation of the galaxy/CGM model and quasar line of sight positioning. In this paper, following a brief review of Paper I, we present the formalism for populating the galaxy/CGM structures with multiphase photoionized and collisionally ionized gas and for generating HI and metal-line absorption profiles. Example absorption line systems through a fiducial galaxy/CGM model are presented. These flexible spatial-kinematic absorption models (SKAMs) can be directly applied to and/or easily modified/expanded for studying individual or ensembles of observed absorption line systems, for exploring various competing theoretical scenarios of the baryon cycle as studied through quasar absorption line systems, and/or serving as pedagogical tools for developing physical intuition. We briefly describe a SKAM GUI that is in early stages of development.

Across the Galactic disk lies a diverse population of X-ray sources, with the fainter end remaining poorly understood due to past survey sensitivity limits. We aim to classify and characterize faint X-ray sources detected in the eROSITA All-Sky Survey (eRASS1) towards the inner Galactic disk ($350^\circ < l < 360^\circ$, $-1^\circ < b < 1^\circ$) using deeper XMM-Newton observations (typical exposure of $\sim 20\,\text{ks}$). We analyzed 189 eRASS1 sources, combining X-ray spectral fitting ($0.2$--$10\,\text{keV}$) with Gaia astrometric and photometric data for robust classification. Our results show that the eRASS1 catalog towards the Galactic disk is overwhelmingly dominated by coronal sources ($\sim 74\%$), primarily active stars and binaries, with $\sim 8\%$ being wind-powered massive stars and $\sim 18\%$ being accreting compact objects. We propose an empirical hardness-ratio cut ($\text{HR} > -0.2$) to efficiently isolate these non-coronal sources. By stacking the classified population and comparing with the Galactic Ridge X-ray Emission (GRXE), we estimate that $\sim 6\%$ of the GRXE flux in the $0.5$--$2.0\,\text{keV}$ band is resolved into point sources above the eRASS1 flux limit ($\sim 5\times 10^{-14}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$). This resolved soft-band emission is dominated by active stars, while hard-band flux originates primarily from X-ray binaries. We conclude that the eRASS1 catalog retains a non-negligible population of compact objects that can be effectively distinguished using X-ray color selection.

The direct, empirical determination of the local value of the Hubble constant (H0) has markedly advanced thanks to improved instrumentation, measurement techniques, and distance estimators. However, combining determinations from different estimators is non-trivial, due to correlated calibrations and different analysis methodologies. Using covariance weighting and leveraging the broad and comprehensive community of experts, we constructed a rigorous and transparent Distance Network (DN) to find a consensus value and uncertainty for the local H0. All critically reviewed the available data sets, spanning parallaxes, detached eclipsing binaries, masers, Cepheids, the TRGB, Miras, JAGB stars, SN Ia, Surface Brightness Fluctuations, SN II, the Fundamental Plane, and Tully-Fisher relations and voted for indicators to define a `baseline' DN and others to assess robustness and sensitivity of the results. We provide open-source software and data products to support full transparency and future extensions of this effort. Our conclusions: 1) Local H0 is robustly determined, with first-rank indicators internally consistent within their uncertainties; 2) A covariance-weighted combination yields an uncertainty of 1.1% (baseline) or 0.9% (all estimators); 3) The contribution from SNe Ia is consistent across four current compilations of optical magnitudes or using NIR-only magnitudes; 4) Removing either Cepheids or TRGB has minimal effect; 5) Replacing SNe Ia with galaxy-based indicators changes H0 by less than 0.1 km/s/Mpc, while doubling its uncertainty; 6) The baseline result is H0=73.50+/-0.81 km/s/Mpc. Compared to early Universe results, our result differs by 7.1sigma from flat {\Lambda}CDM with Planck+SPT+ACT and 5.0 sigma with BBN+BAO (DESI2). A networked approach is invaluable for enabling further progress in accuracy and precision without overreliance on any single method, sample or group.

This review is focused on recent {\sl Self-Organized Criticality (SOC)} literature of astrophysical phenomena, covering the last decade of (2015-2025), while previous SOC literature (1987-2014) is reviewed elsewhere. The selection of literature is mostly based on searches with the NASA-supported {\sl Astrophysics Data System (ADS)}. The discussed astrophysical SOC phenomena are subdivided into solar flares, solar atmosphere (photosphere, chromosphere, corona), heliospheric systems (coronal mass ejections, solar wind, solar energetic particles), planetary systems (asteroids and small bodies, lunar cratering, Saturnian ring systems, magnetospheric systems), stellar flares, and galactic systems (pulsar glitches, gamma ray bursts, soft gamma-ray repeaters, supergiant fast X-ray transients, fast transient radio bursts, magnetars, blazars, black holes).

Throughout their lives, short period exoplanets (<100 days) experience X-ray and extreme-UV (XUV) stellar irradiation that can heat and photoionize planets' upper atmospheres, driving transonic outflows. This photoevaporative mass loss plays a role in both evolution and observed demographics; however, mass loss rates are not currently directly observable and can only be inferred from models. To that end, we present an open-source fast 1D, XUV multi-frequency, multispecies, steady-state, hydrodynamic Parker Wind photoevaporation relaxation model based on Murray-Clay et al. (2009,arXiv:0811.0006). The model can move smoothly between high and low flux regimes and accepts custom multi-frequency stellar spectra. While the inclusion of high-energy X-rays increases mass loss rates ($\dot{M}$), metals decrease $\dot{M}$, and the net result for a typical hot Jupiter is a similar $\dot{M}$, but a hotter, faster, and more gradually ionized wind. We find that mulitfrequency photons (e.g., 13.6-2000eV) are absorbed over a broader range of heights in the atmosphere resulting in a wind-launch radius, $R_{XUV}$, that is of order 10 nanobars for all but the highest surface gravity planets. Grids of H/He solar metallicity atmospheres reveal that, for typical hot Jupiters like HD 209458b, $R_{XUV}$~1.1-1.8$R_P$ for low-fluxes, meaning that the energy-limited mass loss rate, $\dot{M}_{Elim}(R)$, computed at $R=R_P$ is a good approximation. However, for planets with low escape velocities, like many sub-Neptunes and super-Earths, $R_{XUV}$ can be >>$R_P$, making it necessary to use $\dot{M}_{Elim}(R=R_{XUV})$ to avoid significantly underestimating mass loss rates. For both high escape velocities and large incident fluxes, radiative cooling is significant and energy-limited mass loss overestimates $\dot{M}$.

We present a sample of 212 early-type galaxies (ETGs) at redshifts $0.25 < z < 0.75$. We combine deep integral-field spectroscopy from the MUSE-DEEP survey with high-resolution HST imaging to study the structure, kinematics, and stellar populations of these galaxies. We measure spatially resolved stellar kinematics and use the specific angular momentum proxy, $\lambda_R$, to classify galaxies into fast and slow rotators. We find a slow rotator fraction consistent with local Universe samples, suggesting little evolution in the massive ETG population since $z \sim 1$. The kinematic and photometric axes of fast rotators are generally well-aligned, similar to their local counterparts. We find that global stellar population properties, such as age, metallicity, and mass-to-light ratio ($M_*/L$), correlate strongly with the central velocity dispersion ($\sigma_\mathrm{e}$), following trends established for local ETGs. Slow rotators are typically more massive, have higher $\sigma_\mathrm{e}$, and are more metal-rich than fast rotators. Our findings indicate that the fundamental structural, kinematic, and stellar population scaling relations of massive ETGs were already in place by $z \sim 0.75$, suggesting their evolutionary pathways have remained stable over the last $\sim 7$ Gyr.

Understanding how the internal kinematics of massive galaxies evolve is key to constraining the physical processes that drive their assembly. We investigate the evolution of rotational support in massive ($\log M_{\ast}/M_{\odot} \geq 10.6$) early-type galaxies (ETGs) over the past $\sim$7 Gyr. We use MUSE integral-field spectroscopic (IFS) data for 212 ETGs at intermediate redshift ($0.25 < z < 0.75$) from the MAGNUS sample. We compare their kinematics to a carefully matched local sample of 787 ETGs ($z \leq 0.05$) from the MaNGA survey. Using the specific stellar angular momentum proxy, $\lambda_R$, we quantify the balance between ordered rotation and random motions. We derive intrinsic $\lambda_R$ values by applying a uniform correction for seeing and point-spread function (PSF) effects to both samples. We find a significant evolutionary trend: the intermediate-redshift ETGs are systematically more rotationally supported than their local counterparts. The median PSF-corrected $\lambda_R$ for the MAGNUS sample is $0.48 \pm 0.05$, substantially higher than the median of $0.34 \pm 0.03$ for the matched MaNGA sample. This corresponds to a positive slope in the $\lambda_R-z$ relation of $\mathrm{d} \lambda_R / \mathrm{d} z = 0.3 \pm 0.04$ for the combined sample. The decline in rotational support is most pronounced for the most massive galaxies ($\log M_{\ast}/M_{\odot} > 11.3$). Our results provide robust evidence that massive ETGs have undergone significant kinematic evolution, losing angular momentum as they evolve towards the present day, consistent with theoretical models where processes such as dry mergers play a crucial role in shaping the dynamical state of galaxies.

A new chapter is opening in the theory of core-collapse supernovae and neutron star mergers as simulations of these events begin to incorporate fast flavor conversion (FFC) and other forms of neutrino flavor mixing. Using numerical experiments, we show that the approximations of FFC that have been implemented so far are limited by at least two of three factors: (1) approximating continuous evolution as a discrete sequence of instabilities, (2) using spatially homogeneous asymptotic states, and (3) assuming that FFC must be accompanied by instability. The factors we identify in this work will be important considerations as the research area progresses from initial exploratory studies to more quantitatively precise assessments.

In this study, we analyze microlensing events from the 2023 and 2024 observing seasons to identify cases likely caused by binary systems composed of BDs. By applying criteria that the binary-lens events exhibit well-resolved caustics, short time scales ($t_{\rm E} \lesssim 9$ days), and have small angular Einstein radii ($\theta_{\rm E} \lesssim 0.17$~mas), we identify six candidate binary BD events: MOA-2023-BLG-331, KMT-2023-BLG-2019, KMT-2024-BLG-1005, KMT-2024-BLG-1518, MOA-2024-BLG-181, and KMT-2024-BLG-2486. Analysis of these events leads to models that provide precise estimates for both lensing observables, $t_{\rm E}$ and $\theta_{\rm E}$. We estimate the masses of the binary components through Bayesian analysis, utilizing the constraints from $t_{\rm E}$ and $\theta_{\rm E}$. The results show that for the events KMT-2024-BLG-1005, KMT-2024-BLG-1518, MOA-2024-BLG-181, and KMT-2024-BLG-2486, the probability that both binary components lie within the BD mass range exceeds 50\%, indicating a high likelihood that the lenses of these events are binary BDs. In contrast, for MOA-2023-BLG-331L and KMT-2023-BLG-2019L, the probabilities that the lower-mass components of the binary lenses lie within the BD mass range exceed 50\%, while the probabilities for the heavier components are below 50\%, suggesting that these systems are more likely to consist of a low-mass M dwarf and a BD. The brown-dwarf nature of the binary candidates can ultimately be confirmed by combining the measured lens-source relative proper motions with high-resolution imaging taken at a later time.

Gamma Ray Bursts (GRBs) are classified as long (LGRBs) and short (SGRBs) with collapsars and compact object mergers (neutron star (NS)-NS or NS-Black hole) as progenitors, respectively. The former are expected to follow the cosmic star formation rate (SFR), while the latter follows a delayed version of the SFR. However, this division has come under question in several ways, the most prominent being the observational evidence of a significant excess of LGRBs at low redhifts by several independent investigations, summarized in arXiv:2305.15081. This could indicate that the progenitors of low-redshift LGRBs, whose formation rates are delayed, (similar to that of SGRBs) are compact mergers rather than collapsars. Two recent observations of low-redshift LGRBs show associations with kilonovae, a clear feature of compact mergers. Most results showing this separation are based on analyses of small (less than 200) samples of LGRBs with measured redshifts. The aim of this paper is to use a larger sample of LGRBs. The number of LGRBs with measured redshifts has increased by more than a factor of 2 over the last decade. To this data set we add a sample of LGRBs whose redshifts are estimated using a machine learning (ML) method (arXiv:2410.13985). This, in addition to increasing the sample size, reduces the observational selection bias arising from the process of redshift measurement. To account for this bias, we use the non-parametric, non-binning Efron-Petrosian method to establish the degree of correlation between luminosity and redshift, the luminosity evolution, which then allows us to use the Lynden-Bell $C^-$ method to obtain a complete description of the luminosity function. We find similar low redshift excess for the larger sample with measured redshifts. Adding the sources with ML-estimated redshifts, which tend to have more sources in mid-range redshifts, the excess is reduced.

We propose to use the unprecedented polarization sensitivity of PRIMA's PRIMAger Polarization Imager and its high resolution in Band 1 (92 $\mu$m) to map magnetic fields across two contrasting molecular cloud environments: the well-studied Perseus cloud and the isolated Musca filament. This comparative study will leverage the existing VLA radio observations that provide line-of-sight magnetic field component of the Perseus cloud, along with upcoming POSSUM survey results for Musca, to construct the first detailed 3D magnetic field vector maps at sub-parsec resolution. Perseus, with its known formation history through interstellar structure (e.g., bubble) interactions, will reveal how magnetic fields evolve during active star formation phases, while Musca, an isolated filament with lower star-formation activity, will show magnetic field morphology in early evolutionary stages. With PRIMA's resolution of 0.01 pc for Perseus and $<0.01$ pc for Musca, we will resolve magnetic field structures at scales necessary for understanding cloud fragmentation and star formation efficiency, and the roles that magnetic fields play in these processes. Our survey will require approximately 1438 hours to cover 42 deg$^2$ of Perseus and 12 deg$^2$ of Musca. With this, we aim to provide the first comprehensive view of how environment and evolutionary state may influence magnetic field evolution in molecular clouds and how magnetic fields influence cloud formation, fragmentation, and star formation.

The Large Binocular Telescope (LBT) is a world-leading astronomical observatory, where the italian partnership has played an important role in increasing the telescope's productivity, both through an optimized observing strategy and through peer-reviewed publications that are well recognized by the international astronomical community. This manuscript provides an updated overview of the active and past instruments at LBT, together with the key usage statistics. In particular, we analyze the operational performance recorded in the LBT Italia night logs during INAF's observing time and assess the scientific impact of each instrument. Between 2014 and 2025 LBT Italia produced an average of 14 referred publication per year, based on an annual average of 311 hours of on-sky time. This corresponds to approximately 2.2 nights of telescope time per publication. The results of this analysis are placed in an international context to evaluate the competitiveness of LBT, and we outline future perspective for scientific exploitation.

We investigate how the external medium surrounding prestellar cores affects the star formation process by conducting three-dimensional resistive magnetohydrodynamic simulations. The initial cores follow critical Bonnor-Ebert profiles and are embedded in environments with different ambient densities. The simulations follow the evolution at least until the envelope mass within a radius equal to twice the critical Bonnor-Ebert radius drops to 35% of the initial cloud mass. We reveal that in environments with higher external density, enhanced mass inflow from the envelope leads to Bondi-like accretion as the protostellar mass increases. The continued inflow substantially increases the final stellar mass, resulting in star formation efficiencies that appear to exceed unity in dense environments. The external medium also influences the evolution of circumstellar disks and protostellar outflows: with the high-density external medium, disks grow rapidly but their mass becomes smaller relative to the protostellar mass, and the outflow is sustained over a long duration. However, the ratio of angular momentum removed by outflows and magnetic braking to that introduced by inflowing gas decreases with increasing external density. These results suggest that the density of the external medium regulates not only protostellar mass growth but also the inflow-outflow balance and angular momentum transport in magnetized, rotating star-forming cores.

A tidal disruption event (TDE) occurs when a star passes within the tidal radius of a supermassive black hole (SMBH). In TDEs it is expected that the orbital angular momentum of the disrupted star is generally misaligned with the SMBH spin axis, which should result in a misaligned super-Eddington disk precessing around the SMBH spin axis due to the Lense-Thirring effect. In this paper, we investigate the distinct observational signatures produced from such TDE disks, by performing radiative transfer calculations upon previous super-Eddington disk simulations. We demonstrate that the precession of the disk and wind drive time-dependent obscuration and reprocessing of X-ray radiation. Depending on the orientation of the viewing angle of the observer and the tilt angle of the disk, four main types of variability are induced: 1) The smooth-TDEs: The emissions from these TDEs show no fluctuations; 2) The dimmer: The main emission type (X-ray or optical) stays the same, with small to moderate modulations of brightness; 3) The blinker: X-ray and optical emissions take turns to dominate in one cycle of precession, with dramatic changes in the X-ray fluxes. 4) The siren: X-ray and optical emissions take over each other twice per cycle, possibly with two different peak X-ray fluxes within one cycle. In all three scenarios, we observe an inverse correlation between X-ray and optical emissions. Our model provides a unified physical framework for interpreting TDE multi-wavelength variability through disk precession dynamics and gives an alternative interpretation to the interesting case of J045650.3-203750 which was suggested to be a repeated partial TDE previously.

X-ray bursts from the magnetar can lead to the formation of fireballs trapped by the magnetic field and co-rotating with the star. The fireball emission could occasionally be eclipsed by the magnetar, especially when the burst duration is comparable to the magnetar's spin period. In this work, we discover a peculiar type of burst whose light curve has a plateau-like feature among the long bursts of the magnetar SGR J1935+2154. Based on these bursts, we identified four burst candidates with eclipse-like characteristics. By fitting their light curves with the eclipse fireball model, the viewing angle of the magnetar relative to its spin axis is estimated to be $17^\circ \pm 10^\circ$. The distances from the fireballs to the magnetar are found to be more than five times the magnetar's radius, indicating that the fireballs are suspended in the magnetosphere rather than adhering to the magnetar surface. We also find this configuration is well consistent with the implication of the cyclotron resonance scattering feature in their spectra. Our results suggest that some intermediate X-ray bursts of SGR 1935+2154 may originate from magnetic reconnection within the magnetosphere rather than the starquake.

The design of an indigenous Lunar Navigation Satellite System (LNSS) is receiving growing attention due to the surge in planned lunar missions and the limited accessibility of Earth-based Global Navigation Satellite Systems (GNSS) in the cislunar environment. Several studies have explored LNSS architecture using geometric analysis in both near and distant lunar orbits. The existing LNSS optimization efforts have primarily focused on global lunar coverage using analytical station-keeping models with low accuracy. Furthermore, current south pole-focused research is restricted to Elliptical Lunar Frozen Orbits (ELFOs) and lacks comprehensive optimization approach. Additionally, integration with Earth GNSS systems for ephemeris computation and time synchronization has not been adequately addressed in prior studies. In this work, we present a comprehensive LNSS mission design framework based on evolutionary multi-objective optimization integrated with a high-fidelity numerical lunar orbit propagation model. The optimization simultaneously considers navigation performance in the lunar south pole region, semi-analytical continuous station-keeping maneuver model for realistic dV estimate, and GPS-LNSS integration analysis parameters. The resulting Pareto front offers a diverse set of LNSS configurations that balance coverage, accuracy, and dV requirements. Notably, the optimization identifies diverse non-frozen elliptical orbit solutions that achieve over 90% south pole coverage with acceptable navigation accuracy using as few as six satellites and dV of less than 0.4 km/s per satellite per year. This represents a significant reduction in constellation size compared to previous studies, offering a cost-effective yet operationally efficient solution for future LNSS missions.

Addressing the spatial uncertainty and spectral blending challenges in CSST slitless spectroscopy, we present a deep learning-driven, end-to-end framework based on the You Only Look Once (YOLO) models. This approach directly detects, classifies, and analyzes spectral traces from raw 2D images, bypassing traditional, error-accumulating pipelines. YOLOv5 effectively detects both compact zero-order and extended first-order traces even in highly crowded fields. Building on this, YOLO11 integrates source classification (star/galaxy) and discrete astrophysical parameter estimation (e.g., redshift bins), showcasing complete spectral trace analysis without other manual preprocessing. Our framework processes large images rapidly, learning spectral-spatial features holistically to minimize errors. We achieve high trace detection precision (YOLOv5) and demonstrate successful quasar identification and binned redshift estimation (YOLO11). This study establishes machine learning as a paradigm shift in slitless spectroscopy, unifying detection, classification, and preliminary parameter estimation in a scalable system. Future research will concentrate on direct, continuous prediction of astrophysical parameters from raw spectral traces.

We apply an extended Alcock-Paczyński (AP) test to the Sloan Digital Sky Survey data to constrain the dark energy models with the Chevallier-Polarski-Linder (CPL) parametrization of the dark energy equation of state. The extended AP test method uses the full shape of redshift-space two-point correlation funcion(CF) as the standard shape in order to measure the expansion history of the universe. We calibrate the standard shape by using the cosmology-dependent nonlinear evolution of the CF shape in the Multiverse simulations. Further validation of the method and calibration of possible systematics are performed based on mock samples from the Horizon Run 4 simulation. Using the AP test alone, we constrain the flat CDM plus CPL-type dark energy model (flat $w^{\rm CPL}$CDM) to have $\Omega_m=0.289_{-0.029}^{+0.031}$, $w_0=-0.798_{-0.102}^{+0.192}$ and $w_a=-0.165_{-0.945}^{+0.610}$. The result does not show evidence for a dynamically evolving dark energy model. When combined with other results from the low-redshift universe, such as the PantheonPlus supernova compilation and DESI BAO data, the constraint on $w_a$ becomes $w_a=-0.124_{-0.368}^{+0.334}$, which is still consistent with zero.

We present the decadal brightening of non-thermal emission flux in the northeastern filament of the young supernova remnant Cassiopeia A (Cas A), which highlights dynamic processes in the forward shock. This filament, characterized by the highest particle acceleration rate among Cas A's outer shells, offers an exceptional opportunity to investigate underlying astrophysical mechanisms. Since 2000, the non-thermal flux has increased by several tens of percent before plateauing, while the spectral shape has remained largely unchanged. Over the past two decades, the filament's morphology has evolved significantly, splitting into two distinct sections. Detailed analysis reveals contrasting behaviors: one section shows a flux increase followed by saturation, while the other maintains a steady flux. These differences likely arise from nonlinear effects, including magnetohydrodynamic interactions influenced by magnetic field orientation, interactions with surrounding material, and complex fluid dynamics associated with young supernova remnants. The localized evolution of this filament, captured with high spatial resolution, provides critical insights into the temporal dynamics of non-thermal particles and the generation of cosmic rays from asymmetric supernova explosions.

High rates of stable mass transfer (MT) occur for some binary star systems, resulting in luminous transients and circumbinary outflows. We perform hydrodynamical simulations of a $10 \ M_\odot$ donor star and a $5\ M_\odot$ point mass secondary, incorporating approximate effects of radiative cooling. By varying the orbital separation of the system, we probe MT rates between $10^{-5}$ and $10^{-1} M_\odot$/yr. Mass flows from the donor into an accretion disk, with significant equatorially-concentrated outflows through the outer Lagrange point L2 occurring for MT rates $\gtrsim 10^{-3} M_\odot$/yr, while the MT remaining mostly conservative for lower MT rates. In all cases, any outflowing gas approximately carries the specific angular momentum of L2. The gas cooling luminosity $L$ and temperature increases with MT rate, with $L \sim 10^{5} L_\odot$ and $T \sim 10^4 \, {\rm K}$ for simulations featuring the strongest outflows, with contributions from both the accretion disk and circumbinary outflow. The most luminous transients associated with mass outflows will be rare due to the high MT rate requirement, but generate significant optical emission from both the accretor's disk and the circumbinary outflow.

Stellar mass is a fundamental quantity that determines the properties and evolution of stars. However, estimating stellar masses in star-forming regions is challenging because young stars are obscured by dense gas and the regions are highly inhomogeneous, making spherical dynamical estimates unreliable. Supervised machine learning could link such complex structures to stellar mass, but it requires large, high-quality labeled datasets from high-resolution magneto-hydrodynamical (MHD) simulations, which are computationally expensive. We address this by pretraining a vision transformer on one million synthetic fractal images using the self-supervised framework DINOv2, and then applying the frozen model to limited high-resolution MHD simulations. Our results demonstrate that synthetic pretraining improves frozen-feature regression stellar mass predictions, with the pretrained model performing slightly better than a supervised model trained on the same limited simulations. Principal component analysis of the extracted features further reveals semantically meaningful structures, suggesting that the model enables unsupervised segmentation of star-forming regions without the need for labeled data or fine-tuning.

Neutron stars provide unique laboratories for probing physics of dense nuclear matter under extreme conditions. Their thermal and luminosity evolution reflects key internal properties such as the equation of state (EoS), nucleon superfluidity and superconductivity, envelope composition, and magnetic field, and so on. Recent observations [\textit{e.g.}, V. Abramkin \textit{et al.,} ApJ \textbf{924}, 128 (2022)] have revealed unexpectedly warm old neutron stars, which cannot be explained by standard neutrino-photon cooling models. The failure of the standard cooling models implies the presence of additional internal heating mechanism. Building on the previous study [M. Fujiwara \textit{et al}., JCAP \textbf{03}, 051 (2024)], which proposed vortex creep heating (VCH) from the frictional motion of superfluid vortices as a viable mechanism, we extend the cooling framework to include both VCH and direct Urca (DUrca) processes. These are implemented in our code to explore their combined impact, particularly for massive neutron stars where DUrca operates. By varying rotational parameters ($P$, $\dot{P}$, $P_0$), EoS models (APR, BSk24), pairing gaps, and envelope compositions, we examine how heating-cooling interplay shapes the temperature evolution. Our results show that VCH can substantially mitigate the rapid cooling driven by DUrca, offering new evolutionary pathways for massive neutron stars.

Hycean planets are hypothetical exoplanets characterized by $H_2O$ oceans and $H_2$-rich atmospheres. These planets are high-priority targets for biosignature searches, as they combine abundant surface liquid water with easy-to-characterize $H_2$-rich atmospheres. Perhaps their most unusual climate feature is convective inhibition, which can dramatically alter a planet's temperature structure. However, so far hycean planets have mostly been investigated using 1D models that do not account for convective inhibition, and its effects are still poorly understood. This work develops pen-and-paper theory to analyze the effects of moist convective inhibition on hycean planets. The theory is tested and verified against a 1D radiative-convective model. We show that hycean planets near the onset of convective inhibition can exhibit either bistability or oscillations, due to the inhibition layer's trapping of heat and moisture. Meanwhile, hot hycean planets exhibit multistability, in which the inhibition layer and surface climate show multiple stable equilibria due to the lack of constraints on the water cycle inside the inhibition layer. The water cycle inside the inhibition layer is influenced by numerous processes that are challenging to resolve in 1D, including turbulent diffusion, convective overshoot and large-scale circulations. Our results demonstrate that hycean planets have unexpectedly rich climate dynamics. Meanwhile, previous claims about hycean planets should be treated with caution until confirmed with more self-consistent 1D and 3D models; this includes the claim that K2-18b might be habitable, and the proposal to infer $H_2O$ oceans on sub-Neptunes from JWST measurements of chemical species in their upper atmospheres.

Growing observational evidence suggests that Class II protoplanetary disks may undergo substantial interactions with their environment in the form of late infall. This mass inflow predominantly manifests itself in the form of so-called streamers: filaments and arcs of gas connecting large-scale, extended gas structures to disk scales. Prevalent late infall has far-reaching consequences for planet formation theory, challenging the long-standing treatment of evolved disks in isolation. In this work, we investigate the emergence of late-infall streamers in different formation scenarios, their morphology and multiplicity, as well as their dependence on environmental conditions. We conduct this investigation by performing 3D hydrodynamical simulation using the grid-based code FARGO3D, which we post-process to obtain synthetic observations using the Monte Carlo radiative transfer code RADMC3D. We find that, while a late infall event in the form of a single encounter with a "cloudlet" of gas can produce a streamer via an interplay between the fall-back of bound material and shocks, such features dissipate quickly, on a timescale of ~10 kyr. Furthermore, we find that streamers can also form naturally in a turbulent, dense environment without the need for such encounters, which could act to reconcile short-lived streamers with ubiquitous detection of these structures. Here, we find multiple co-existing streamers for a disk velocity relative to the interstellar medium of $v_\mathrm{sys}=0.5~\mathrm{km}~\mathrm{s}^{-1}$ and a turbulent velocity dispersion of $\sigma_\mathrm{turb}=0.5~\mathrm{km}~\mathrm{s}^{-1}$. We find considerable dependence of the streamer morphology on the environment, which may act as a utility to constrain the physical conditions of the gas surrounding planet-forming disk, and therefore the conditions under which planets form.

In-situ observations of the fast solar wind in the inner-heliosphere show that minor ions and ion sub-populations often exhibit distinct drift velocities. Both alpha particles and proton beams stream at speeds that rarely exceed the local Alfvén speed relative to the core protons, suggesting the presence of instabilities that constrain their maximum drift. We aim to propose a mechanism that generates an alpha-particle beam through non-linear Landau damping, primarily driven by the relative super-Alfvénic drift between protons and alpha particles. To investigate this process, we perform one-dimensional, fully kinetic particle-in-cell simulations of a non-equilibrium multi-species plasma, complemented by its linear theory to validate the model during the linear phase. Our results provide clear evidence that the system evolves by producing an alpha-particle beam, thereby suggesting a local mechanism for alpha-beam generation via non-linear Landau damping.

[abbreviated] Long Secondary Periods (LSPs) are observed in about one third of pulsating red giants yet remain unexplained. Four key observational constraints anchor the discussion: (i) a roughly 30 percent occurrence rate in semi-regular variable AGB stars (SRVs), with a much lower rate or absence in regularly pulsating Mira-type AGB stars (Miras), (ii) about 50 percent of LSP stars show a secondary mid-infrared minimum, (iii) Keplerian fits to radial-velocity curves favor argument of periastron greater than 180 degrees, and (iv) the radial-velocity to light-curve phase lag clusters around minus 90 degrees. We test whether a close-in, eccentric, low-mass companion that spends only part of its orbit within the giant's dust-formation zone can match all four empirical constraints. Guided by observed radial-velocity amplitudes and periods of about 500 to 1500 days, we adopt companion masses of 0.08 to 0.25 solar masses, orbital separations of 1.5 to 3 AU, eccentricities up to 0.6, and dust condensation radii of about 2.5 to 3 AU for SRVs (larger for Miras via scaling with luminosity). We compute the fraction of the orbit spent outside the condensation radius and apply line-of-sight criteria for detectability. We test the first three empirical constraints analytically, then model the phase offset using three-dimensional hydrodynamical simulations. Our scenario explains the observed excess of systems with argument of periastron greater than 180 degrees and yields LSP detectability of about 27 percent for SRVs and about 2.5 percent for Miras, with a conditional secondary mid-infrared eclipse fraction of about 44 percent. The hydrodynamical models place the optical-depth maximum downstream of the companion near apastron and then shift it forward by roughly 90 to 225 degrees later in the orbit, consistent with the observed radial-velocity to light-curve offsets.

Cosmological hydrodynamical simulations are essential tools for studying the formation and evolution of galaxies and their central supermassive black holes. While they reproduce many key observed properties of galaxies, their limited volumes have hindered comprehensive studies of the AGN and quasar populations. In this work, we leverage the FLAMINGO simulation suite, focusing on its large $(2.8$ $\mathrm{Gpc})^3$ volume, to investigate two key observables of quasar activity: the quasar luminosity function (QLF) and quasar clustering. FLAMINGO reproduces the observed QLF at low redshift ($z \lesssim 1$) and for faint quasars ($L_\mathrm{bol} \lesssim 10^{45}$ $\mathrm{erg s^{-1}}$), but significantly underpredicts the abundance of bright quasars at $z \approx 1$-$3$. Introducing a 0.75 dex log-normal luminosity scatter to represent unresolved small-scale variability boosts the number of bright quasars by upscattering lower-luminosity systems, thereby improving agreement with observations at the bright end. A decomposition of the QLF by black hole mass reveals that this boost is primarily driven by low-mass black holes radiating above the Eddington limit. Nevertheless, limitations remain in fully reproducing the rise and decline of the bright quasar population over cosmic time and in matching the black hole masses inferred from quasar spectra. Thanks to FLAMINGO's large volume, we can robustly sample rare, luminous quasars and measure their spatial clustering for $\log_{10} L_\mathrm{bol}/\mathrm{erg s^{-1}} \gtrsim 45.5$. The simulation reproduces the observed clustering across $0 \lesssim z \lesssim 3$, and the reduced luminosity dependence introduced by scatter aligns with observational trends. However, it underpredicts the clustering strength at $z \approx 4$, consistent with other models and possibly reflecting high-redshift observational uncertainties.

Quasar absorption systems not only affect the way quasars are selected, but also serve as key probes of galaxies, providing insight into their chemical evolution and interstellar medium (ISM). Recently, a method based on Gaia astrometric measurements has aided the selection of quasars reddened by dust hitherto overlooked. We conducted a spectroscopic study using VLT/X-Shooter on one such dust-reddened quasar, Q2310-3358. This quasar, at $z = 2.3908\pm0.0003$, is associated with a Damped Lyman-alpha absorber (DLA) at nearly the same redshift $2.4007\pm0.0003$, with a neutral hydrogen column density of $\log N(\mathrm{H\,I}) = 21.214 \pm 0.003$. The DLA is very metal-rich (close to the Solar metallicity after correction for depletion on dust grains). Its properties align with the metal-to-dust ratio and the mass-metallicity relation established in previous large samples of DLAs. Surprisingly, given its proximity to the quasar in redshift, the absorber has strong cold gas characteristics, including CI and H$_2$. Based on the derived kinetic temperature of $71^{+28}_{-15}$~K, we infer the presence of a strong UV radiation field, which in turn suggests that the quasar and the DLA are in close proximity, i.e. part of the same galaxy and not just different objects in the same overdensity of galaxies. We use the line ratios of the CI fine-structure lines to constrain the density of the cold gas, yielding $n_{\rm H} \sim 10^{3}~\mathrm{cm}^{-3}$. Our analysis extends the understanding of $z_{abs} \approx z_{em}$ absorption line systems and provides valuable constraints on the interplay between dust, metals, and neutral gas in the ISM of early galaxies.

We report on the jet structure in SS 433 based on X-ray high resolution spectroscopy with the XRISM/Resolve. The source was observed over 5 days covering both inside and outside an eclipse of the compact object by the companion star. Doppler-shifted, ionized Fe and Ni K emission lines were resolved, as well as lower-energy lines including Si and S K lines. Time-resolved spectral analysis showed that Fe and Ni K line widths were $1020 \pm 40$ km s$^{-1}$ (corresponding the 1$\sigma$ width) in the eclipse phase, gradually increased during the egress, and reached $1740 \pm 30$ km s$^{-1}$ outside the eclipse. A time-averaged spectrum outside the eclipse confirmed that the Fe and Ni K lines in 5.5-9 keV are significantly broader than the Si and S K$\alpha$ emission lines in 2--4 keV. Specifically, the width in 5.5-9 keV was measured to be $1900 \pm 80$ km s$^{-1}$, whereas the width in 2-4 keV is $1300^{+300}_{-400}$ km s$^{-1}$ for the approaching (blueshifted) jet component. These results indicate that radial velocity dispersion of the jet plasma in SS 433 decreases as it moves outward. We interpret this variation as progressive jet collimation along its axis, as suggested by Namiki et al. (2003), or a decrease in turbulence in the jet plasma flow within the X-ray emitting region. We also detected a clear difference in velocity dispersion between the approaching and receding (redshifted) jet components in the 5.5-9 keV band outside eclipse. The receding jet exhibited a smaller velocity dispersion ($1400 \pm 200$ km s$^{-1}$) than the approaching jet. Since the observation was conducted when the approaching jet was tilted toward the observer, this may suggest that the receding jet was more extensively occulted by the accretion disk.

We present a detailed chemical analysis study of 16 candidate metal-poor stars, previously identified with 2dF+AAOmega, using X-Shooter spectra and the Korg 1D LTE spectral synthesis code. We generally confirm the earlier metallicity estimates and reveal four EMP ([Fe/H] $< -3$) stars in the current sample. Two of these stars, including the most metal-poor at [Fe/H] = $-3.83 \pm 0.07$, are kinematically associated with the GSE accretion event, increasing the number of known GSE stars with [Fe/H] $< -3.5$ to eight. From the X-Shooter spectra we also determine abundances for 16 elements, with the element-to-iron abundance ratios generally consistent with high-resolution studies of Galactic halo stars. Within the sample, we identify three peculiar stars: the first is a GSE nitrogen enhanced metal poor (NEMP; [N/Fe] $= 1.62 \pm 0.10$ and [C/Fe] $= 0.27 \pm 0.08$) star with unusually high Na ([Na/Fe] $= 2.28 \pm 0.07$) and Li (A(Li)$_{\rm 3DNLTE} = 1.90 \pm 0.08$) abundances, but which lacks any enhancement in [Al/Fe] or [Mg/Fe]. The second is a halo r-II star significantly enhanced in Sr ([Sr/Ba] $= 0.37 \pm 0.08$), suggesting mixture of r-process and s-process enrichment, uncommon for r-II stars. Whilst the third is a halo star very depleted in N ([N/Fe] $< -1.13$), with low C ([C/Fe] $= -0.36 \pm 0.08$) and otherwise 'normal' [X/Fe] abundances, suggesting enrichment with Type II supernova that proceeds enrichment from massive AGB stars. This study reveals the substantial degree of chemical diversity in the stellar populations which assembled the early Milky Way.

Multiple populations (MPs) is a intra-star cluster phenomenon consisting in star-to- star variation of the abundance of some light chemical elements. They have been observed in many star clusters, most of them old globular clusters, populating the Milky Way and other satellite galaxies. Since the study of MPs became more system- atic, different astrophysical parameters have been claimed to be the main responsible for its occurrence. However, at the present time, no attempt would seem to have solved this conundrum. This work deals with a potential trigger of the MPs phenomenon, based on the gathered observational evidence of the existence of MPs in some star clusters and the absence of its in others. We found that star clusters with MPs mostly formed during time intervals of intense star formation activity in a galaxy, for instance during the galaxy formation epoch, a close galaxy encounter, etc. At those time inter- vals where relative peaks in the galaxy star formation rate occur, star clusters with masses above a lower mass limit harbour MPs. This lower star cluster mass limit would marginally depend on the star cluster age.

In this two-paper series, we present a straightforward mathematical model for synthesizing quasar absorption line profiles from sight lines through idealized, spatial-kinematic models of the circumgalactic medium (CGM) and their host galaxies. Here, in Paper I, we develop the spatial geometries of multiple galaxy/CGM structures and populate these structures with 3D velocity fields. For arbitrary viewing angles and galaxy-quasar impact parameters, we derive observer coordinate-based expressions for the perceived azimuthal angle and galaxy inclination and a generalized scalar expression for the line-of-sight velocity as a function of position along the line of sight. We motivate and develop four idealized galaxy/CGM spatial-kinematic structures based on empirical data and theoretical predictions: (1) a rotating galactic disk/extra-planar gas, (2) a static or dynamic spherical halo, (3) an outflowing bi-polar galactic wind, and (4) an inward spiraling flared planar accretion. Using a small set of free parameters, the spatial geometries and velocity fields can be adjusted and explored, including velocity gradients, wind stalling, and accretion trajectories. These spatial-kinematic models are designed to be flexible and easily modified and can be tailored for studying individual galaxy-absorber pairs or galaxy group environments; they can be applied to real-world observations or hydrodynamic simulations of the baryon cycle as studied through quasar absorption line systems. These models also serve as tools for developing physical intuition. In Paper II, we will present the formalism for populating the galaxy/CGM structures with multiphase photoionized and collisionally ionized gas and for generating absorption profiles for ions of interest.

We present half-light radius measurements for the dwarf satellites of Andromeda, based on multi-epoch imaging from the Isaac Newton Telescope (INT) Monitoring Survey of Local Group dwarf galaxies. This analysis is conducted within a larger study to identify long-period variable (LPV) stars in these galaxies. The survey was performed with the Wide Field Camera on the 2.5-m INT and covers multiple epochs obtained between 2015 and 2018 in the i (Sloan) and V (Harris) bands. To determine the half-light radii, we derived surface brightness and number density profiles for each system and fitted them with Exponential, Plummer, and Sérsic models. The resulting half-light radii are in good agreement with literature values but reveal subtle variations linked to differences in stellar distribution and morphology. Distances were independently estimated using the Tip of the Red Giant Branch (TRGB) method, yielding values consistent with previous determinations. The complete photometric and variability catalogs will be made publicly available through CDS/VizieR, providing a valuable resource and foundation for future studies of the structure, stellar populations, and evolution of Andromeda's dwarf companions.

We selected Bulge stars from APOGEE DR17 cross-matched with astrometric data from \textit{Gaia} DR3. Bulge stars were divided into sub-samples with line-of-sight velocity dispersion analyzed and the peaks of MDF were detected by both Gaussian Mixture Models (GMM) and \texttt{this http URL\_peaks}. GMM is also conducted to kinematically distinguish the metal-poor and metal-rich populations. Analyses were put on the Bulge stars (including retrograde stars), their elemental abundances, and the [Mg/Mn]-[Al/Fe] plane to investigate potential accreted components. Finally, the shapes (X-shaped/boxy) of Bulge stars with different metallicities were analyzed through least-squares fitting based on the analytical Bulge models. By studying the kinematic, elemental and structural dependences on metallicity for Bulge stars, our findings are concluded as follows: 1. Six peaks are detected in the Bulge MDF, encompassing values reported in previous studies, suggesting a complex composition of Bulge populations. 2. An inversion relationship is well-observed in metal-rich sub-samples, while absent in metal-poor sub-samples. 3. Metal-poor populations exhibit larger dispersions than metal-rich stars (which is also revealed by GMM decomposition), suggesting that metal-rich stars are kinematically coherent. 4. Retrograde stars are confined to $\sim1$ kpc of the Galactic center, with their relative fraction decreasing at higher [Fe/H] -- a trend potentially linked to the ``spin-up'' process of Galactic disks. 5. Metal-rich Bulge stars with [Al/Fe] $<-0.15$ are likely associated with from disk accreted substructure, while all elemental planes exhibit bimodality but Na abundances rise monotonically with metallicity. 6. In general, stars with all metallicities support a boxy profile.

High-energy neutrino astronomy has advanced rapidly in recent years, with IceCube, KM3NeT, and Baikal-GVD establishing a diffuse astrophysical flux and pointing to promising source candidates. These achievements mark the transition from first detections to detailed source studies, motivating next-generation detectors with larger volumes, improved angular resolution, and full neutrino-flavour sensitivity. We present a performance study of large underwater neutrino telescopes, taking the proposed TRIDENT array in the South China Sea as a case study, with a focus on comparing the performance of various detector configurations against the TRIDENT baseline design. Both track-like events primarily from muon neutrinos, which provide precise directional information, and cascade events from all flavours, which offer superior energy resolution, diffuse-source sensitivity, and all-sky flavour coverage, are included to achieve a balanced performance across source types. The time to discover potential astrophysical sources with both track- and cascade-like events is used as the figure of merit to compare a variety of detector design choices. Our results show that, for a fixed number of optical modules, simply enlarging the instrumented volume does not inherently lead to improved performance, while taller strings can provide modest gains across all detector channels, within engineering constraints. Distributing dense clusters of strings over a large volume is found to generally worsen discovery potential compared to the baseline layout. Finally, the optical properties of the sea-water emerge as the key factor dictating the optimisation of detector layout, highlighting the need for in-situ measurements and early deployment of optical modules to guide the final array configuration.

Asteroseismic modelling of high-amplitude $\delta$ Scuti (HADS) variables critically depends on their accurate classification, which provides robust constraints on stellar physical parameters. As a foundational step in this direction, we present a detailed analysis of the pulsational behavior of 50 HADS star candidates using high-precision photometric data from the Transiting Exoplanet Survey Satellite (TESS). We confirm 46 as genuine HADS variables, with 40 stars having their dominant frequencies identified for the first time. Among them, 7 pulsate solely in the fundamental mode; 21 exhibit the fundamental mode alongside at least one low-amplitude nonradial mode; 5 are pure double-mode pulsators (the fundamental and first-overtone modes), 13 show double-mode pulsations accompanied by additional low-amplitude nonradial modes. The remaining four stars are classified as other types of variables: two (TIC 69546708 and TIC 110937533) are confirmed hot subdwarfs, one (TIC 8765832) is a cataclysmic variable, and one (TIC 32302937) is a likely hot subluminous star, but it requires further spectroscopic confirmation. We investigate the period-luminosity (PL) relation, also known as the Leavitt law, for these 46 confirmed HADS stars, deriving a revised relation: $M_{V}= (-3.31 \pm0.39)~\mathrm{log}~P+(-1.68 \pm 0.39)$. This result is consistent with previous studies. Their distribution in the Hertzsprung-Russell diagram indicates that HADS stars are not strictly confined to be within a narrow instability strip previously found, but can extend beyond it, with a distribution toward lower temperatures. The refined classifications presented here establish a high-quality sample for precise asteroseismic modelling and enhance the potential for future machine-learning-assisted searches and classifications of HADS stars in large-scale photometric surveys.

Planetary nebulae (PNe) are valuable spatial and kinematic tracers of nearby galaxies. In this mini-review, I focus on their role in tracing the halo and intra-cluster/intra-group light assembly in groups and clusters of galaxies within 100~Mpc and, in particular, the link between characteristic PN metrics such as the $\alpha$-parameter and the PN luminosity function and changes from the underlying in-situ to ex-situ stellar populations. These results from nearby groups and clusters are placed into context with integral-field spectroscopic surveys of galaxies, which allow the co-spatial measurement of PN and stellar population properties. I provide an outlook on upcoming instrumentation that will provide new opportunities for the study of extragalactic PN populations. I address the challenges of reconciling observations of extragalactic PN populations with predictions from stellar evolution models and how revised late-stellar-evolution models have alleviated some of the tensions between observations and theory.

The origin of ultra-high-energy cosmic rays (UHECRs) remains an open questions in astrophysics. We explore two primary scenarios for the distribution of UHECR sources, assuming that their production rate follows either the cosmic star-formation-rate or stellar-mass density. By jointly fitting the UHECR energy spectrum and mass composition measured by the Pierre Auger Observatory above the ankle (10^{18.7} eV), we derive constraints on the acceleration mechanisms, source energetics, and elemental abundances at escape. Using these constraints, we generate sky maps above 40 EeV based on a catalog of over 400,000 galaxies out to 350 Mpc, providing a near-infrared flux-limited sample that maps the two stellar-activity tracers across the full sky. A crucial factor in understanding UHECR propagation is the influence of large-scale cosmic structures, particularly galaxy clusters, the largest gravitationally bound systems in the Universe, which are filled with magnetized diffuse plasma. Intermittent sources hosted in galaxies within such structures, coupled with cosmic magnetic fields, shape the observed UHECR arrival directions and provide insights into the burst rate of the sources. We show that these environments can significantly impact UHECR transport, making them particularly opaque to heavy nuclei. Additionally, we compute the expected secondary neutrino and photon fluxes from UHECR interactions in these environments and compare them with current experimental limits, constraining the maximum energy that particles can achieve. Finally, we assess the compatibility of these constraints with astrophysical candidates, identifying long gamma-ray bursts as the most promising sources.

Modern simulations and observations in Astronomy & Cosmology (A&C) produce massively large data volumes, posing significant challenges for storage, access and data analysis. A long-standing bottleneck in high-performance computing, especially now in the exascale era, has been the requirement to write these large datasets to disks, which limits the performance. A promising solution to this challenge is in-situ processing, where analysis and visualization are performed concurrently with the simulation itself, bypassing the storage of the simulation data. In this work, we present new results from an approach for in-situ processing based on Hecuba, a framework that provides a highly distributed database for streaming A&C simulation data directly into the visualization pipeline to make possible on-line visualization. By integrating Hecuba with the high-performance cosmological simulator ChaNGa, we enable real-time, in-situ visualization of N-body simulation results using tools such as ParaView and VisIVO.

The Astronomy & Cosmology (A&C) community is presently witnessing an unprecedented growth in the quality and quantity of data coming from simulations and observations. Writing results of numerical simulations to disk files has long been a bottleneck in high-performance computing. To access effectively and extract the scientific content of such large-scale data sets appropriate tools and techniques are needed. This is especially true for visualization tools, where petascale data size problems cannot be visualized without some data filtering, which reduces either the resolution or the amount of data volume managed by the visualization tool. A solution to this problem is to run the analysis and visualization concurrently (in-situ) with the simulation and bypass the storage of the full results. In particular we use Hecuba, a framework offering a highly distributed database to stream A\&C simulation data for on-line visualization. We will demonstrate the Hecuba platform integration with the Changa high performant cosmological simulator and the in-situ visualization of its N-body results with the ParaView and VisIVO tools.

COMCUBE-S (Compton Telescope CubeSat Swarm) is a proposed mission aimed at understanding the radiation mechanisms of ultra-relativistic jets from Gamma-Ray Bursts (GRBs). It consists of a swarm of 16U CubeSats carrying a state-of-the-art Compton polarimeter and a BGO spectrometer to perform timing, spectroscopic and polarimetric measurements of the prompt emission from GRBs. The mission is currently in a feasibility study phase (Phase A) with the European Space Agency to prepare an in-orbit demonstration. Here, we present the simulation work used to optimise the design and operational concept of the microsatellite constellation, as well as estimate the mission performance in terms of GRB detection rate and polarimetry. We used the MEGAlib software to simulate the response function of the gamma-ray instruments, together with a detailed model for the background particle and radiation fluxes in low-Earth orbit. We also developed a synthetic GRB population model to best estimate the detection rate. These simulations show that COMCUBE-S will detect about 2 GRBs per day, which is significantly higher than that of all past and current GRB missions. Furthermore, simulated performance for linear polarisation measurements shows that COMCUBE-S will be able to uniquely distinguish between competing models of the GRB prompt emission, thereby shedding new light on some of the most fundamental aspects of GRB physics.

We stack 3.75 Megaseconds of early XRISM Resolve observations of ten galaxy clusters to search for unidentified spectral lines in the $E=$ 2.5-15 keV band (rest frame), including the $E=3.5$ keV line reported in earlier, low spectral resolution studies of cluster samples. Such an emission line may originate from the decay of the sterile neutrino, a warm dark matter (DM) candidate. No unidentified lines are detected in our stacked cluster spectrum, with the $3\sigma$ upper limit on the $m_{\rm s}\sim$ 7.1 keV DM particle decay rate (which corresponds to a $E=3.55$ keV emission line) of $\Gamma \sim 1.0 \times 10^{-27}$ s$^{-1}$. This upper limit is 3-4 times lower than the one derived by Hitomi Collaboration et al. (2017) from the Perseus observation, but still 5 times higher than the XMM-Newton detection reported by Bulbul et al. (2014) in the stacked cluster sample. XRISM Resolve, with its high spectral resolution but a small field of view, may reach the sensitivity needed to test the XMM-Newton cluster sample detection by combining several years worth of future cluster observations.

Carbonaceous nano-grains are a significant component of interstellar dust and dominate the mid-infrared emission of photodissociation regions (PDRs). We study the evolution of nano-grains across the illuminated edge of the Horsehead PDR, especially their abundance and size properties. This work is part of the Physics and Chemistry of PDR Fronts program studying dust and gas in PDRs with JWST. We use NIRCam+MIRI photometric bands and NIRSpec+MRS spectroscopy to map the illuminated edge. We model dust emission using the THEMIS dust model with the SOC radiative transfer code. Detailed modeling of high angular resolution JWST data allows us to obtain constraints on nano-grain properties. We find that diffuse ISM dust cannot account for the observed data, requiring evolved grains. A sharp density increase is observed at the illuminated edge, consistent with ALMA observations revealing a sharp transition between molecular and ionized gas. Although the PDR length could not be directly determined, we estimate an upper limit of approximately 0.015 pc. This implies a lower limit on small grain abundance (greater than 0.003), showing small grains are not depleted at the Horsehead edge, unlike in the Orion Bar. Our findings indicate a high-density environment and less steep size distribution for nano-grains at the illuminated edge versus the diffuse ISM. This implies nano-grain destruction mechanisms might be less efficient in the Horsehead's moderate-UV field than in more intense PDRs. These results support a model where nano-grain population recovery is slower in moderate-UV environments, leading to a unique dust size distribution at the edge of the Horsehead Nebula.

Reconstructions of solar spectral irradiance - especially in the ultraviolet (UV) range - are crucial for understanding Earth's climate system. Although total solar irradiance (TSI) has been thoroughly investigated, the spectral composition of solar radiation offers a deeper insight into its interactions with the atmosphere, biosphere, and climate. UV radiation, in particular, plays a key role in stratospheric chemistry and the dynamics of stratospheric ozone. Reconstructing solar irradiance over the past centuries requires accounting for both the cyclic modulation of active-region coverage associated with the 11-year solar cycle and the longer-term secular trends, including their centennial variability. This study utilizes an empirical framework, based on a 1000-year record of Open Solar Flux, to characterize the various temporal components of solar irradiance variability. We then combine these components to reconstruct Solar UV irradiance variations in spectral bands crucial for Earth's atmospheric studies.

Observations of $\gamma$-rays from diffuse gas provide the opportunity to study the distribution of high energy particles in different astrophysical environments. In the circumgalactic medium (CGM) and the intracluster medium (ICM), it is expected that relativistic cosmic rays collide with thermal particles and produce $\gamma$-rays through pion decay. The $\gamma$-ray luminosity of a plasma depends on where cosmic rays are: if they are in denser gas, they produce more $\gamma$-rays. In this work, we study how different cosmic-ray transport mechanisms impact the $\gamma$-ray luminosity of a turbulent, multiphase medium formed from an initially diffuse medium. Two quantities set the luminosity: the average cosmic-ray energy density and the correlation of cosmic-ray energy and gas density. Overall, cosmic rays must escape cold dense regions in order to produce less $\gamma$-ray emission and be consistent with observations. Our simulations with fast transport mechanisms (either diffusion or streaming) are degenerate: they each produce a lower $\gamma$-ray luminosity than slow transport simulations by two orders of magnitude. This result means that fast transport (particularly in dense clumps) is necessary for simulations to agree with the dearth of observations of $\gamma$-ray emission from diffuse gas like the CGM and ICM. We also show the significant difference in luminosity is the result of cosmic-ray reacceleration. This reacceleration is different from the turbulent reacceleration described by Ptuskin (1988). Instead, condensing, cold clouds drive a significant increase in the average cosmic-ray energy and, as a result, the $\gamma$-ray luminosity.

The presence of rings and moons around exoplanets is likely to be one of the next great discoveries in exoplanet research. Using theories developed for the Solar System, we explore the possibility of coupled ring-moon cycles around exoplanets and what these processes mean for the observability of these features. Around Neptune- and Earth-like planets, we find that ring-moon cycles are capable of producing long-lived rings of comparable and greater relative transit depths than Saturn's. In multi-planet systems, secular spin-orbit resonances can provide the necessary planetary obliquity for these rings to contribute noticeably to transit lightcurves. We model the geometry of a ring's cross-section at various angles in comparison to the cross-section of a transiting planet to determine whether the ring may be detectable during the planet's transit. Ringed planets have also been proposed as an alternative to inflated planetary radii seen in transit, leading to abnormally low observed densities. Ring-moon cycles can produce late-forming and sometimes long-lived rings that can have the potential of explaining at least some of these observations. We also discuss some inconsistencies in the calculation of exoplanet oblateness due to rotation that we have come across in the course of this work.

We present the largest ground-based catalogue of M-dwarf flares to date, comprising 1,229 time-resolved events identified in Zwicky Transient Facility Data Release 17. Using high-cadence ZTF observations collected between April 2018 and September 2020, we analyzed over 93 million variable light curves containing 4.1 billion photometric measurements. Flare candidates were identified through a machine-learning pipeline trained on simulated light curves generated by injecting TESS-based flare templates into ZTF data and then refined through an extensive post-filtering stage combining an additional classifier, metadata checks, and human inspection. For 655 flares with reliable Gaia-based distances and well-sampled light curves, we derived bolometric energies ranging from 10^31 to 10^35 erg. A clear correlation is observed between flare frequency and spectral subtype, with a sharp increase toward later M dwarfs, particularly near M4-M5, coinciding with the transition to full convection. Using 680 flaring stars with known vertical distances from the Galactic plane, we find that the fraction of flaring stars decreases with increasing Galactic height. The resulting catalogue provides the most comprehensive ground-based sample of M-dwarf flares and establishes a framework for flare detection and classification in upcoming wide-field surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time.

In the era of Planck cosmology, the inflationary paradigm is best fitted towards the cosmological attractor scenarios, including the induced inflation, universal attractors, conformal attractors, and special attractors that are cataloged as $\xi$-models and $\alpha$-models. The recent hint from the ACT results pushes the scalar spectral index closer to the scale-invariant Harrison-Zeldovich spectrum, calling for a theoretical paradigm shift towards a Harrison-Zeldovich attractor, which is difficult to realize in the standard single-field slow-roll inflationary scenario. In this Letter, we achieve the Harrison-Zeldovich attractor scenario via nonminimal derivative coupling, attracting the monomial inflation, hilltop inflation, and $\alpha$-attractor E-model towards the Harrison-Zeldovich spectrum.

Late-time light curve plateaus in tidal disruption events (TDEs) are often approximated as flat and time-independent. This simplification is motivated by theoretical modeling of spreading late time TDE disks, which often predicts slow light curve evolution. However, if time evolution can be detected, late-time light curves will yield more information than has been previously accessible. In this work, we re-examine late-time TDE data to test how well the flat plateau assumption holds. We use Markov Chain Monte Carlo to estimate the maximum likelihood for a family of theory-agnostic models and apply the Akaike information criterion to find that that roughly one third of our sample favors evolving plateaus, one third favors truly flat plateaus, and one third shows no statistically significant evidence for any plateau. Next, we refit the TDEs that exhibit statistically significant plateaus using a magnetically elevated $\alpha$-disk model, motivated by the lack of clear thermal instability in late time TDE light curves. From these model-dependent fits, we obtain estimates for the supermassive black hole (SMBH) mass, the mass of the disrupted star, and the $\alpha$ parameter itself. Fitted $\alpha$ values range from $10^{-3}$ to 0.4 (the mean fitted $\alpha=10^{-1.8}$, with scatter of 0.6 dex), broadly consistent with results from magnetohydrodynamic simulations. Finally, we estimate the timescales of disk precession in magnetically elevated TDE models. Theoretically, we find that disk precession times may be orders of magnitude shorter than in unmagnetized Shakura-Sunyaev disks, and grow in time as $T_{\rm prec}\propto t^{35/36}$; empirically, by using fitted $\alpha$ parameters, we estimate that late time disks may experience $\sim$few-10 precession cycles.

Highly-irradiated giant exoplanets known as "ultra-hot Jupiters" are anticipated to exhibit large variations of atmospheric temperature and chemistry as a function of longitude, latitude, and altitude. Previous observations have hinted at these variations, but the existing data have been fundamentally restricted to probing hemisphere-integrated spectra, thereby providing only coarse information on atmospheric gradients. Here we present a spectroscopic eclipse map of an extrasolar planet, resolving the atmosphere in multiple dimensions simultaneously. We analyze a secondary eclipse of the ultra-hot Jupiter WASP-18b observed with the NIRISS instrument on JWST. The mapping reveals weaker longitudinal temperature gradients than were predicted by theoretical models, indicating the importance of hydrogen dissociation and/or nightside clouds in shaping global thermal emission. Additionally, we identify two thermally distinct regions of the planet's atmosphere: a "hotspot" surrounding the substellar point and a "ring" near the dayside limbs. The hotspot region shows a strongly inverted thermal structure due to the presence of optical absorbers and a water abundance marginally lower than the hemispheric average, in accordance with theoretical predictions. The ring region shows colder temperatures and poorly constrained chemical abundances. Similar future analyses will reveal three-dimensional thermal, chemical, and dynamical properties of a broad range of exoplanet atmospheres.

Despite the well-established presence of amorphous compounds on planetary bodies such as the Moon and Mercury due to space weathering, the direct effect of atomic arrangement on the surface binding energies (SBEs) of elements on these bodies remains largely unexplored. Accurate SBE values are essential for reliably predicting sputtering yields and the energy distribution of ejecta. Here, we use molecular dynamics simulations to quantify SBEs for the different elements sputtered from amorphous atomic arrangements of the plagioclase feldspar end members, albite and anorthite, and compare to their crystalline counterparts. Results show that while the mean elemental SBEs from amorphous surfaces are not significantly different from their crystalline counterparts, the random orientation in amorphous structures gives rise to a spectrum of bonding configurations, resulting in a distribution of SBEs with a wider range. This contrasts with the clearly discretized set of SBE values associated with the ordered atomic structure of crystalline surfaces. We then consider sputtering by H, He, and a solar wind combination of 96% H and 4% He. For each of these cases, we demonstrate that there is minimal difference (<10% for albite and <20% for anorthite) between the sputtering yields of amorphous and crystalline surfaces. We attribute these results to the presence of the same elemental bonds across different atomic arrangements, which leads to similar mean SBEs and, consequently, comparable sputtering yields.

We study Big Bang Nucleosynthesis (BBN) constraints on heavy QCD axions. BBN offers a powerful probe of new physics that modifies the neutron-to-proton ratio during the process, thanks to the precisely measured primordial Helium-4 abundance. A heavy QCD axion provides an attractive target for this probe, because not only is it a well-motivated hypothetical particle by the strong CP problem, but also it dominantly decays to hadrons if kinematically allowed. A range of its lifetime is thus excluded where the hadronic decays would significantly alter the neutron-to-proton ratio. We compute axion-induced modification of the neutron-to-proton ratio, and obtain robust upper bounds on the axion lifetimes, as low as 0.017 s for the axion mass higher than 300 MeV. Remarkably, this is stronger than projected future CMB bounds via $N_{\rm eff}$. Our bounds are largely insensitive to uncertainties in hadronic cross sections and the axion's branching fractions into various hadrons, as well as to the precise value of the initial axion abundance. We also incorporate, for the first time, several key improvements, such as scattering processes by energetic $K_L$ and secondary hadrons, that can also be important for studying general hadronic injections during BBN, not limited to those from axion decays.

The growing evidence for nano-hertz gravitational waves, from NANOGrav and other observations, may be pointing to a cosmological first-order phase transition at temperatures of $\mathcal{O}(10-100)\;\mathrm{MeV}$. Such an interpretation requires dynamics beyond the Standard Model in this energy range. If so, it may well be the case that core-collapse supernova explosions would recreate the first-order phase transition leaving a unique imprint on the spectrum of neutrinos emitted in the initial few seconds. This scenario is also suggestive of a low-mass seesaw mechanism to explain neutrino masses. We outline the prospects for future observations of Galactic supernovae to uncover the signals of this scenario, which could get further confirmation with additional pulsar timing array data establishing the primordial origin of the observed nano-hertz gravitational waves.

We study the scalar potentials that arise from higher curvature corrections in general $f(R)$ theories of gravity and their connection to a dynamical species scale. Starting from general considerations in arbitrary dimensions, we show that at large field values, the scalar potential generated by an infinite series of curvature terms and the field dependent species scale arising from circle compactification both decay exponentially, in complementary ways. We identify conditions under which these two effects precisely balance out, giving rise to exponentially flat, plateau-like potentials. We additionally find a precise embedding of Starobinsky inflation consistent with the Swampland program, and we discuss possible implications the mechanism proposed could have for M and string theory.

The QCD axion can successfully solve the strong CP problem under the condition that the Peccei-Quinn symmetry is respected to extremely high standard. We explore a class of gauge theories that accommodate a high-quality axion, known as the Barr-Seckel models, paying special attention to the cosmology of topological defects. In models with domain wall number equal to unity, we show that axion strings with winding number larger than one can always be decomposed into a number of unit-winding axion strings and pure gauge strings. This mechanism enables the axion string-wall network to be destroyed timely to render a viable cosmology. The subsequent decay of gauge strings into gravitational waves generically produces a double-plateau in frequency space allowing the mechanism to be tested by upcoming experiments.

The unstably bound, critical null geodesics of the Kerr spacetime form a distinguished class of orbits whose properties govern observables such as the photon ring and the high-frequency component of black-hole ringdown. This set of orbits defines a codimension-two submanifold of the null-geodesic phase space known as the photon shell. In this work we investigate the photon shell's intrinsic symplectic geometry. Using the induced symplectic form, we construct the canonical volume form on the shell and compute the differential phase-space volume it encloses as a function of radius -- equivalently, the radial density of states. In the near-extremal limit the photon shell bifurcates into near-horizon and far-region components; we find that approximately $3\%$ of the shell's phase-space volume resides in the near-horizon component. We also analyze a thickening of the photon shell that includes near-critical orbits, and compute its differential phase-space volume. Beyond their intrinsic theoretical interest, these results may inform the interpretation of high-resolution observations of spinning black holes.

We present a comprehensive study of secluded dark matter (DM) $\chi$, where the relic abundance is set by annihilations into lighter dark mediators $\phi$ that couple only feebly to the Standard Model (SM). In contrast to canonical WIMPs, which are now strongly constrained by direct and indirect searches, secluded models still achieve the observed relic abundance via thermal freeze-out into hidden-sector mediators, while predicting highly suppressed present-day signals. We analyze three minimal models: (i) a $U(1)_X$ gauge boson ($A'$) with kinetic mixing; (ii) a scalar DM candidate $S$ with a scalar mediator $K$ that has a trilinear vertex; and (iii) a Dirac fermion $\chi$ whose mass arises from a Higgs-mixed singlet $H_p$. For each model we derive annihilation and scattering rates in both WIMP-like and secluded regimes, and solve the Boltzmann equations: a single-species equation for the WIMP case and a coupled $\chi$-$\phi$ system for the secluded case to account for possible early departure of the mediator from thermal equilibrium with the SM bath. In this regard, we provide explicit lower limits on the portal coupling $\epsilon$ required to keep the mediator in thermal equilibrium with the SM bath and to ensure mediator decay before BBN. We show that for portal couplings $\epsilon \ll 10^{-3}$ the relic density is dominantly controlled by DM annihilation into mediator pairs, while spin-independent scattering lies well below current limits and remains viable even for future experiments approaching the irreducible neutrino background floor. Indirect constraints are typically weak due to $p$-wave suppression, off-resonance $s$-channels, and cascade spectra controlled by $\epsilon^2$. Finally, we highlight the most promising collider tests, which remain sensitive despite tiny portal couplings.

We present a theory-agnostic framework to search for extra-fundamental modes in the ringdown phase of black hole mergers. These modes, distinct from standard General Relativity (GR) deviations, originate from modifications of the theory of gravity. Applying our method to the events GW150914, GW190521, and GW200129, we find no significant evidence for extra modes beyond those expected in GR.

A gravitational-wave background (GWB) arising from the superposition of numerous unresolved gravitational-wave signals has yet to be detected. Potential contributing sources to such a background include compact binary coalescences (CBCs) and core-collapse supernovae (CCSNe). In this work, we place upper limits on the gravitational-wave energy emitted by CCSNe using cross-correlation measurements made with Advanced LIGO and Advanced Virgo detectors during their third observing run (O3). Specifically, we obtain a $95\%$ credibility upper limit of $0.01~ {M_\odot c^2}$ while accounting for the contribution from CBC sources to a GWB. This result improves on previous constraint obtained from initial LIGO data by approximately two orders of magnitude. We also explore the detection prospects of third-generation ground-based detectors such as the Einstein Telescope and Cosmic Explorer for both individual CCSNe events and the GWB. Our results show that single events are likely to be detected prior to the GWB.

We introduce a higher-derivative action for dark matter whose energy-momentum tensor describes an imperfect fluid with nonzero pressure, energy flux, and anisotropic stress. In the limit where the higher-derivative couplings are switched off, the energy-momentum tensor reduces to pressureless dust. A systematic derivation follows from extending the singular conformal transformation used in the mimetic dark matter scenario to include higher-derivative terms while the resulting action is general and does not rely on the mimetic framework. On a homogeneous cosmological background, the dynamics coincides with that of pressureless dust, while in the presence of inhomogeneities the higher-derivative terms generate nonzero acceleration and vorticity, making it possible to avoid the formation of caustic singularities even if the strong energy condition satisfies. In particular, within the mimetic realization these terms resolve the usual caustic pathology of mimetic dark matter.

We present the invariant characterization of the region where vacuum breakdown into electron-positron ($e^+e^-$) pairs occurs due to an overcritical electric field, the dyadoregion, in the case of a Kerr black hole (BH) in the presence of an external, asymptotically uniform test magnetic field aligned with the BH rotation axis, using the Wald solution. We calculate the dyadoregion morphology, the electromagnetic energy available for the pairs, the pair-creation rate, the number density of pairs, the average energy per pair, and their energy density and pressure. These results provide initial conditions for simulating the subsequent dynamics of the pair-produced plasma and astrophysical applications in the context of high-energy transients involving BHs in strong electromagnetic fields.

The symmetry energy expansion was developed to connect isospin symmetric matter probed in nuclear experiments to asymmetric matter found in neutron stars. Using the isospin asymmetry derived from the Gell-Mann-Nishijima formula, we derive the symmetry energy expansion for quark matter that has unique properties compared to hadronic matter. To test our methods, we use perturbative Quantum Chromodynamics (pQCD) calculations at next-to-leading-order, where realistic quark masses can be included. We find that pQCD at electroweak equilibrium is not isospin symmetric but rather obtains a small skewness term in the symmetry energy expansion. We predict that if equations of state for nuclear matter must match pQCD results, then a non-monotonic dip in the symmetry energy would appear.

We present a new model of the dark sector involving Dirac fermion dark matter, with axial coupling to a dark photon which provides a portal to Standard Model particles. In the non-relativistic limit, this implies that the dominant effective operator relevant to direct detection is ${\cal O}_8$. The resulting event rate for direct detection is suppressed by either the dark matter velocity or the momentum transfer. In this scenario there are much wider regions of the dark parameter space that are consistent with all of the existing constraints associated with thermal relic density, direct detection and collider searches.

Modified gravity theories with a nonminimal coupling between curvature and matter offer a compelling alternative to dark energy and dark matter by introducing an explicit interaction between matter and curvature invariants. Two of the main consequences of such an interaction are the emergence of an additional force and the non-conservation of the energy--momentum tensor, which can be interpreted as an energy exchange between matter and geometry. By adopting this interpretation, one can then take advantage of many different approaches in order to investigate the phenomenon of gravitationally induced particle creation. One of these approaches relies on the so-called irreversible thermodynamics of open systems formalism. By considering the scalar--tensor formulation of one of these theories, we derive the corresponding particle creation rate, creation pressure, and entropy production, demonstrating that irreversible particle creation can drive a late-time de Sitter acceleration through a negative creation pressure, providing a natural alternative to the cosmological constant. Furthermore, we demonstrate that the generalized second law of thermodynamics holds: the total entropy, from both the apparent horizon and enclosed matter, increases monotonically and saturates in the de Sitter phase, imposing constraints on the allowed particle production dynamics. Furthermore, we present brief reviews of other theoretical descriptions of matter creation processes. Specifically, we consider approaches based on the Boltzmann equation and quantum-based aspects and discuss the generalization of the Klein--Gordon equation, as well as the problem of its quantization in time-varying gravitational fields.

A recent study shows that incorporating a new term into the thermodynamic potential density, as required by the thermodynamic consistency criterion, can effectively resolve the thermodynamic inconsistency problems of the conventional perturbative QCD model. This additional term plays a crucial role in resolving inconsistencies at relatively low densities and becomes negligible at extremely high densities. Within this revised perturbative QCD model, we find that if we require only that the energy per baryon of up-down ($ud$) quark matter exceeds 930 MeV so as not to contradict the standard nuclear physics, the maximum mass of an $ud$ quark star allowed by the revised perturbative QCD model can reach up to 2.17 $M_{\odot}$. From this perspective, the observed 2.14 $M_{\odot}$ pulsar PSR J0740+6620 may be an $ud$ quark star. However, if we further impose the constraint that the tidal deformability of a 1.4 $M_{\odot}$ $ud$ quark star must be consistent with the GW170817 event, the maximum mass allowed by the revised perturbative QCD model would decrease to no more than 2.08 $M_{\odot}$. Consequently, our results suggest that the compact object with a mass of 2.50-2.67 $M_{\odot}$, as observed in the GW190814 event, cannot be an $ud$ quark star, according to the revised perturbative QCD model.

By critically evaluating higher-order nonlinear effects to the isotope shifts (ISs) in the low-lying transition frequencies of the singly charged calcium ion, stringent constraint on the electron-neutron coupling due to a hypothetical boson describing physics beyond the Standard Model is inferred. It shows an order magnitude difference compared to the previously reported limit demonstrating importance of higher-order effects in the analysis of nonlinearity in the King's plot. The first-order IS parameters and enhancement factor ($D$) were evaluated using two complementary approaches in the relativistic coupled-cluster theory framework: namely finite-field (FF) and analytical response (AR) approaches. Extraction of the second-order IS parameters in the FF approach show numerical instabilities, so they are determined in the AR approach. Comparison of these factors with previous calculation shows substantial differences in the magnitudes. However, $D$ values from both the FF and AR approaches display excellent agreement. We also show explicitly roles of electron correlation effects in the evaluation of $D$ values accurately.

Accurately observing the rarefied media of the upper atmosphere, exosphere, and planetary and solar system environments and beyond requires highly sensitive metrological techniques. We present the operating concept and architecture of an in-situ sensing solution based on the dynamics of a levitated nanoparticle (levitodynamics). It can detect and measure impacts of individual particles in rarefied media. Dubbed `LEVITAS', our sensor consists of a dispenser of dielectric nanoparticles and optical trapping of a single nanoparticle in the focus of a laser beam. The trapped nanoparticle constitutes a harmonic oscillator at frequencies in the kilohertz range whose position can be tracked at the standard quantum limit by interferometric detection of the laser photons it scatters. Here, we simulate microcanonical impacts on the nanoparticle and show that the density, velocity, temperature, and composition of the surrounding medium can be estimated accurately. We illustrate the performance of LEVITAS in circumstances ranging from low Earth orbit out to exospheric distances, across which individual impacts can be detected at favourable rates. Furthermore, LEVITAS may be employed to accurately measure highly rarefied neutral distributions within vastly different areas of momentum space. This we demonstrate by simulating the measurement of high-velocity neutral gas particles from the interstellar medium penetrating the heliosphere and flowing through our solar system.

We investigate the dynamical instability of a self-gravitating thermal system in the quantum regime, where Fermi degeneracy pressure becomes significant. Using a truncated Fermi-Dirac distribution and solving the Tolman-Oppenheimer-Volkoff equation, we identify marginally stable configurations following Chandrasekhar's criterion. While Fermi pressure stabilizes a system against gravitational collapse in Newtonian gravity, in general relativity it can instead drive the instability, enabling collapse even at low temperatures. We discuss implications for the formation of massive black holes in the early Universe through the gravothermal collapse of dark matter.

Frontier AI agents show increasing promise as scientific research assistants, and may eventually be useful for extended, open-ended research workflows. However, in order to use agents for novel research, we must first assess the underlying faithfulness and correctness of their work. To evaluate agents as research assistants, we introduce ReplicationBench, an evaluation framework that tests whether agents can replicate entire research papers drawn from the astrophysics literature. Astrophysics, where research relies heavily on archival data and computational study while requiring little real-world experimentation, is a particularly useful testbed for AI agents in scientific research. We split each paper into tasks which require agents to replicate the paper's core contributions, including the experimental setup, derivations, data analysis, and codebase. Each task is co-developed with the original paper authors and targets a key scientific result, enabling objective evaluation of both faithfulness (adherence to original methods) and correctness (technical accuracy of results). ReplicationBench is extremely challenging for current frontier language models: even the best-performing language models score under 20%. We analyze ReplicationBench trajectories in collaboration with domain experts and find a rich, diverse set of failure modes for agents in scientific research. ReplicationBench establishes the first benchmark of paper-scale, expert-validated astrophysics research tasks, reveals insights about agent performance generalizable to other domains of data-driven science, and provides a scalable framework for measuring AI agents' reliability in scientific research.

We study cosmological gravitational particle production (CGPP) in Higgs inflation, wherein the inflaton is a scalar field with quartic self-coupling $\lambda$ and a nonminimal coupling to gravity $\xi$, and which may, but need not be, the Higgs boson of the Standard Model (SM). We find an explosive particle production peaked on a characteristic comoving wavenumber $k\sim \xi^{2/3} a H$ with a peak occupation number that scales with $\xi$. This new peak in production can easily dominate over the conventional (minimally coupled inflation) CGPP even for modest values of $\xi$. The results apply for a wide range of $\xi$, e.g., as low as $\xi=10$, which can be realized for the Standard Model Higgs given suitable RG flow of the quartic coupling. We discuss implications for late time relics such as dark matter.