Papers for Thursday, Feb 12 2026

Artificial intelligence offers much promise, but its use in scientific research should be restrained so that the primary aim of academia -- advancing knowledge for humans -- is safeguarded.

Solar vortices are fundamental components of solar atmospheric dynamics, serving as natural laboratories for magnetic field twisting, energy concentration and transport, wave guidance, and plasma coupling across atmospheric layers. Numerical and observational studies show that solar vortices are intimately connected to key physical processes including magnetic reconnection, atmospheric heating, turbulence, and wave generation. This white paper, prepared for the UK Space Frontiers 2035 call, outline five high-priority scientific questions addressing vortex generation mechanisms, cross-layer coupling, magnetic restructuring, collective wave-guidance structures, and their role in triggering explosive events and modulating the solar wind. Key observations and capabilities required to make significant advancements over the coming decade are identified. The UK solar physics community has established world-leading expertise in vortex dynamics, combining strengths in high-resolution observations, MHD turbulence theory, numerical modelling, and space instrumentation. UK researchers have made foundational contributions to Solar Orbiter, delivered critical systems for DKIST, and maintain active involvement in MUSE and SOLAR-C EUVST missions. Our technical approach centres on developing next-generation instrumentation: a multi-band, space-qualified system employing four tunable Fabry-Pérot Interferometers providing diffraction-limited, high-cadence spectropolarimetric coverage from the deep photosphere to the low corona. This capability will be validated through a staged mission architecture beginning with balloon-borne demonstrators. Continuing this effort over the coming decade is vital to maintain UK leadership in this field and achieve the goals of roadmap for solar system research.

Reconstructing the early Universe from the evolved present-day Universe is a challenging and computationally demanding problem in modern astrophysics. We devise a novel generative framework, Cosmo3DFlow, designed to address dimensionality and sparsity, the critical bottlenecks inherent in current state-of-the-art methods for cosmological inference. By integrating 3D Discrete Wavelet Transform (DWT) with flow matching, we effectively represent high-dimensional cosmological structures. The Wavelet Transform addresses the ``void problem'' by translating spatial emptiness into spectral sparsity. It decouples high-frequency details from low-frequency structures through spatial compression, and wavelet-space velocity fields facilitate stable ordinary differential equation (ODE) solvers with large step sizes. Using large-scale cosmological $N$-body simulations, at $128^3$ resolution, we achieve up to $50\times$ faster sampling than diffusion models, combining a $10\times$ reduction in integration steps with lower per-step computational cost from wavelet compression. Our results enable initial conditions to be sampled in seconds, compared to minutes for previous methods.

Recently, Jones et al. [arXiv:2310.12859] claimed strong evidence for the statistical anisotropy of the universe. The claim is based on a joint analysis of four different anomaly tests of the cosmic microwave background data, each of which is known to be anomalous, with a lower level of significance. They reported a combined $p$-value of about $3\times 10^{-8}$, which is more than a $5\sigma$ level of significance. We observe that statistical anisotropy is not even relevant for two of the four considered tests, which seems sufficient to invalidate the authors' claim. Furthermore, even if one reinterprets the claim as evidence against $\Lambda$CDM rather than statistical anisotropy, we argue that this result significantly suffers from the look-elsewhere effect. Assuming a set of independent (i.e., uncorrelated) tests, we show that if the four tests with the smallest $p$-values are cherry-picked from 10 independent tests, the $p$-value reported by Jones et al. corresponds to only $3\sigma$ significance. If there are 27 independent tests, the significance falls to $2\sigma$. These numbers, however, overstate our argument, since the four tests used by Jones et al. are slightly correlated. Determining the correlation of Jones et al.'s tests by comparing their joint $p$-value with the product of the four separate $p$-values, we find that about 16 or 50 tests are sufficient to reduce the significance of Jones et al.'s results to 3$\sigma$ or 2$\sigma$ significance, respectively. We also provide a list of anomaly tests discussed in the literature (and propose a few generalizations), suggesting that very plausibly 16 (or even 50) independent tests have been published, and possibly many more have been considered but not published. We conclude that the current data is consistent with the $\Lambda$CDM model and, in particular, with statistical isotropy.

Tidal disruption events (TDEs) have traditionally been discovered in optical sky surveys through targeted searches of nuclear transients. However, it is expected that some TDEs will occur outside the galaxy nucleus, arising from wandering black holes originating in galaxy mergers. Here we present observations of TDE 2025abcr, the first optical TDE discovered in the outskirts of a host galaxy. The TDE was identified by a custom 'off-nuclear' implementation of the ML classifier $\texttt{tdescore}$, which classifies new ZTF transients based on their lightcurves. Follow-up observations confirm that TDE 2025abcr is a TDE-H+He, occurring 9.5$"$ (10.3 kpc projected distance) from the nucleus of a massive galaxy ($\mathrm{M}_{\star}$ = $10^{11.18 \pm 0.03}\mathrm{M}_{\odot}$) with a central black hole mass of $10^{8.82 \pm 0.65}\mathrm{M}_{\odot}$. TDE 2025abcr itself was likely disrupted by a much lighter black hole ($10^{6.09\pm0.53}\mathrm{M}_{\odot}$, as estimated with peak luminosity scaling relations). The black hole was either dynamically ejected from the nucleus or lies at the center of a very faint tidally-stripped dwarf galaxy undergoing a minor merger. Late-time observations of TDE 2025abcr could confirm the origin of this apparent 'orphan' black hole. The rate of highly offset ($\gtrsim$3 kpc) TDEs can be constrained to $<$10% of the nuclear TDE rate, but our discovery implies that many dozens of similar sources will be detected by the Vera C. Rubin each year with resolvable offsets.

At time of writing, large language models (LLMs) are beginning to obtain the ability to design, execute, write up, and referee scientific projects on the data-science side of astrophysics. What implications does this have for our profession? In this white paper, I list - and argue for - a set of facts or "points of agreement" about what astrophysics is, or should be; these include considerations of novelty, people-centrism, trust, and (the lack of) clinical value. I then list and discuss every possible benefit that astrophysics can be seen as bringing to us, and to science, and to universities, and to the world; these include considerations of love, weaponry, and personal (and personnel) development. I conclude with a discussion of two possible (extreme and bad) policy recommendations related to the use of LLMs in astrophysics, dubbed "let-them-cook" and "ban-and-punish." I argue strongly against both of these; it is not going to be easy to develop or adopt good moderate policies.

The initial mass function (IMF) describes the distribution of stellar masses in a population of newly born stars and is amongst the most fundamental concepts in astrophysics. It is not only the direct result of the star formation process but it also explains the evolution of galaxies' luminosities, metal yields, star-formation efficiencies, and supernova production rates. Because most stars exist in binary systems, however, a full statistical account of stellar mass requires not the IMF but rather the joint distribution of a binary population's primary- and secondary-star masses. This joint distribution must respect the IMF of the stars from which the population has been assembled as well as the distribution of mass ratios that results from the assembly mechanism. Despite its importance, this joint distribution is known only in the case of random pairing. Here we present a method for constructing it in the general case. We also illustrate the use of our method by recovering the known result for random pairing and by finding the previously unknown result for uniform pairing.

In this work, we present an updated prescription of contemporary tidal dissipation theory adapted for rapid binary population synthesis. Our simplified expressions encode the dependence of tidal dissipation on stellar structure, stratification, and tidal forcing frequency, while remaining computationally efficient. We implement these prescriptions in the rapid population synthesis code COMPAS, and demonstrate the self-consistent coupling of tides with stellar evolution and binary properties such as orbital periods, spins, and eccentricities for several representative binary systems. When compared with commonly used tidal prescriptions, our equilibrium tidal dissipation efficiencies can be stronger by 1-2 orders of magnitude for low mass main sequence and giant type stars, and dynamical tides can be stronger by 1-7 orders of magnitude due to the explicit dependence on internal stellar structure and the presence of inertial wave dissipation. Despite our simplistic approach, our models agree with detailed stellar simulations to within an order of magnitude across tidal dissipation mechanisms.

Stellar collisions can occur frequently in dense cluster environments, and play a crucial role in producing exotic phenomena from blue stragglers in globular clusters to high-energy transients in galactic nuclei. Successive collisions and mergers of massive stars could also lead to the formation of massive black holes, serving as seeds for supermassive black hole in the early universe. While analytic fitting formulae exist for predicting collision outcomes, they do not generalize across different energy scales or stellar evolutionary phases. Smoothed particle hydrodynamics (SPH) simulations are often used to compute the outcomes of stellar collisions, but, even at low resolution, their computational cost makes running on-the-fly calculations during an $N$-body simulation quite challenging. Here we present a new grid of $27,720$ SPH calculations of main-sequence star collisions, spanning a wide range of masses, ages, relative velocities, and impact parameters. Using this grid, we train machine learning models to predict both collision outcomes (merger vs disruption, or flyby) and final remnant masses. We compare the performance of nearest neighbors, support vector machines, and neural networks, achieving classification balanced accuracy of $98.4\%$, and regression relative errors as low as $0.11\%$ and $0.15\%$ for the final stars $1$ and $2$, respectively. We make our trained models publicly available as part of the package collAIder, enabling rapid predictions of stellar collision outcomes in $N$-body models of dense star cluster dynamics.

We investigate how different galaxy formation models impact the stellar properties of brightest group galaxies (BGGs) in four cosmological simulations: ROMULUS, SIMBA, SIMBA-C, and OBSIDIAN. The stellar masses, specific star formation rates, and mass-weighted stellar ages of the simulated BGGs are analysed alongside those of observed BGGs from X-ray-selected galaxy groups in the COSMOS field. We find that the global properties and underlying evolutionary pathways of simulated BGG populations are strongly impacted by the strength and mechanism of their respective active galactic nucleus (AGN) feedback models, which play a critical role in regulating the growth of massive galaxies. OBSIDIAN's sophisticated three-regime AGN feedback model achieves the highest overall agreement with COSMOS observations, matching stellar property distributions, quenched fractions, and the evolution of star formation in increasingly massive systems. We find evidence suggesting that BGG populations of OBSIDIAN and COSMOS undergo a gradual decline in star formation with stellar mass, in contrast to SIMBA and SIMBA-C, which display rapid quenching linked to the onset of powerful AGN jet feedback. By comparison, ROMULUS produces highly star-forming, under-quenched BGGs due to the inefficiency of its thermal AGN feedback in preventing cooling flows from fuelling BGG growth. The success of the OBSIDIAN simulation demonstrates the importance of physically motivated subgrid prescriptions for realistically capturing the processes that shape BGGs and their dynamic group environments.

The stellar content of galaxies is tightly connected to the mass and growth of their host dark matter halos. Observational constraints on this relation remain limited, particularly for low-mass groups, leaving uncertainties in how galaxies assemble their stars across halo mass scales. Accurately measuring the brightest central galaxy (BCG) stellar-to-halo mass relation (SHMR) over a wide mass range is therefore crucial for understanding galaxy formation and the role of feedback processes. Here we present the SHMR spanning $M_{\rm halo} \sim 10^{12}$-$10^{15}\,M_\odot$, using halo masses derived from eROSITA eRASS1 X-ray data and BCG stellar masses based on SDSS photometry. By stacking X-ray spectra of optically selected groups, we recover robust average halo gas temperatures for each bin, which are then converted to halo masses via the $M$-$T_X$ relation. We find that the SHMR peaks near $M_{\rm halo} \sim 10^{12}\,M_\odot$, with a declining stellar fraction at higher masses. This trend reflects a combination of processes that reduce the efficiency of stellar mass growth in massive halos, such as AGN feedback, reduced cooling efficiency, and the increasing dominance of ex-situ assembly, while halos continue to grow through mergers and accretion. Our measurements are consistent over the full mass range with previous observational studies, including weak lensing, X-ray analyses of individual clusters, and kinematical and dynamical methods. Comparisons with hydrodynamical simulations show good agreement at low masses but reveal significant discrepancies in the normalization at cluster scales, highlighting the sensitivity of BCG stellar growth to feedback prescriptions and halo assembly history. These results provide the first X-ray-based observational SHMR covering three orders of magnitude in halo mass, establish a robust benchmark for testing galaxy formation models.

We present the second set of discoveries from Keck Observations in the INfrared of Taurus and $\rho$ Oph Exoplanets And Ultracool dwarfs (KOINTREAU), an adaptive optics survey of young stars in the Taurus and $\rho$ Oph star-forming regions using Keck/NIRC2 in conjunction with the Keck infrared pyramid wavefront sensor. We have discovered two faint comoving companions to young stars ISO-Oph 96 and 2MASS J16262785-2625152. The companion to ISO-Oph 96, KOINTREAU-3b, is at a projected separation of 340 au (2.49"). Using our NIRC2 photometry and evolutionary models, and assuming that the companion has the same extinction as its host star, we infer that KOINTREAU-3b has a mass of $3.4\pm0.7$ M$_{\rm Jup}$. The companion to 2MASS J16262785-2625152, KOINTREAU-4b, has a projected separation of 180 au (1.25") and could have a mass of either $11.5^{+1.2}_{-1.6}$ M$_{\rm Jup}$ or $15.3^{+0.7}_{-0.8}$ M$_{\rm Jup}$, depending on whether the host star is a member of $\rho$ Oph or Upper Sco.

We present the deepest and highest-resolution radio continuum imaging of the Galactic globular cluster Terzan 5, one of the most crowded locations in the radio sky. In these new 2$-$4 GHz Karl G. Jansky Very Large Array images, we detect 38 of the 49 confirmed pulsars, including extensive multi-frequency eclipse mapping of the luminous redback Ter5A. Nonetheless, there is still a large amount of diffuse residual flux from pulsars that are fainter than our 2.5 GHz continuum detection limit of $\sim 11\,\mu$Jy. Using a range of approaches including image-based simulations, we model the fluxes of the detected pulsars together with the residual flux. We find a minimum total population of $N\sim250$ detectable pulsars in Terzan 5 and perhaps substantially more, though the luminosity function remains very uncertain. Consideration of the $\gamma$-ray properties of the cluster, though also not unambiguous to interpret, leads to consistent conclusions. These pulsar population estimates are larger than inferred from previous work and highlight Terzan 5 as a keystone target for next-generation radio facilities.

We present a homogeneous spectroscopic dataset of 22,339 stars in 78 Milky Way dwarf galaxy satellites and globular clusters. All data were taken with the Keck II telescope and Deep Extragalactic Imaging Multiobject Spectrograph (DEIMOS) spectrograph using the 1200G grating (spectral resolution R~6000). Based on a uniform data reduction of 411 DEIMOS masks, we present a catalog of individual stellar radial velocities, equivalent width-based [Fe/H] metallicities, and membership estimates. The Milky Way satellites range from M_V = 2 to -14 (M* = 10^1.5 to 10^7.5 Msun); the majority of individual stars presented in these systems have magnitudes 17 > r > 22. The data were reduced to 1D spectra using PypeIt, which provides near Poisson statistics-level sky subtraction. Radial velocities were determined via dmost, a forward modeling method first presented here, which combines both synthetic telluric and stellar templates to determine stellar radial velocities. We assess the accuracy and precision our method via comparison to thousands of repeat measurements and literature values. We determine a velocity error floor of 1.1 km/s and a CaII triplet-based metallicity error floor of 0.1 dex. We calculate internal velocity dispersions and compare to literature values, demonstrating 20-50% improved precision over the literature in most cases. In a companion paper, we use our homogeneous catalogs to explore properties of these Milky Way satellites, including previously unpublished measurements in several systems including Bootes II and Draco II. We provide full access to the data catalogs to enable further studies.

Population-level studies of Milky Way satellites used to constrain dark matter or the threshold of galaxy formation often rely on velocity dispersions and metallicities derived from heterogeneous spectroscopic analyses. Systematic differences between data reduction pipelines and membership criteria can masquerade as astrophysical signals, or obscure real trends. Here, we present the largest self-consistent sample of spectroscopically-derived quantities for Milky Way satellite galaxies and globular clusters based on a homogeneous re-analysis of individual stars observed with the Keck/DEIMOS spectrograph. We determine enclosed dynamical masses, mean [Fe/H] metallicities, and metallicity dispersions for 67 systems with 10 or more member stars. At a given stellar mass, systems classified as satellite galaxies are well separated from globular clusters in their dynamical mass and mass-to-light ratios. The average enclosed mass densities of satellite galaxies agree with semi-analytic CDM model predictions. For satellite galaxies, we observe a break in the stellar mass-metallicity relation near log M*/Msun = 4 (M_V ~ -4.5). Above this stellar mass, satellite galaxies show the well-known tight trend (0.16 dex scatter in [Fe/H]) of decreasing metallicity with stellar mass; below log M*/Msun = 4, the mass-metallicity relation flattens and/or increases in scatter. Satellite galaxies have internal metallicity scatter between 0.3-0.4 dex across our stellar mass range. These uniform measurements will enable tighter constraints on the stellar mass-halo mass relation, improved J-factor estimates for dark matter searches, and lay a foundation for interpreting the flood of new Milky Way satellites expected in the LSST/Roman/Euclid era.

Photometric redshifts (photo-$z$'s) will be crucial for studies of galaxy evolution, large-scale structure, and transients with the Nancy Grace Roman Space Telescope. Deep learning methods leverage pixel-level information from ground-based images to achieve the best photo-$z$'s for low-redshift galaxies, but their efficacy at higher redshifts with deep, space-based imaging remains largely untested. We used Hubble Space Telescope CANDELS optical and near-infrared imaging to evaluate fully-supervised, self-supervised, and semi-supervised deep learning photo-$z$ algorithms out to $z\sim3$. Compared to template-based and classical machine learning photometry methods, the fully-supervised and semi-supervised models achieved better performance. Our new semi-supervised model, PITA (Photo-$z$ Inference with a Triple-loss Algorithm), outperformed all others by learning from unlabeled and labeled data through a three-part loss function that incorporates images and colors for all objects as well as redshifts when available. PITA produces a latent space that varies smoothly in magnitude, color, and redshift, resulting in the best photo-$z$ performance even when the redshift training set was significantly reduced. In contrast, the self-supervised approach produced a latent space with significant color and redshift fluctuations that hindered photo-$z$ inference. Looking forward to Roman, we recommend using semi supervised deep learning to take full advantage of the information contained in the hundreds of millions of high-resolution images and color measurements, together with the limited redshift measurements available, to achieve the most accurate photo-$z$ estimates for both faint and bright sources.

IceCube recently reported the observation of TeV neutrinos from the nearby Seyfert galaxy NGC~1068, and the corresponding neutrino flux is significantly higher than the upper limit implied by observations of GeV-TeV gamma rays. This suggests that neutrinos are produced near the supermassive black hole, where the radiation density is high enough to obscure gamma rays. We use a set of muon neutrinos with interaction vertices inside the detector, which have good sensitivity to sources in the Southern sky, from IceCube data recorded between 2011 and 2021. We then search for individual and collective neutrino signals from 14 Seyfert galaxies in the Southern Sky selected from the Swift Burst Alert Telescope (BAT) AGN Spectroscopic Survey. Using the correlations between keV X-rays and TeV neutrinos predicted by disk-corona models, and assuming production characteristics similar to NGC~1068, a collective neutrino signal search reveals an excess of $6.7_{-3.2}^{+4.0}$ events, which is inconsistent with background expectations at the 3$\sigma$ level of significance. In this paper, we present new independent evidence that Seyfert galaxies contribute to the extragalactic flux of high-energy neutrinos.

A common-envelope (CE) phase occurs when a star engulfs its companion and is widely considered the primary channel for producing Luminous Red Novae (LRNe). In this study, we combine binary-population synthesis with stellar-evolution calculations to systematically estimate the mass, velocity, and launching radius of ejecta produced during coalescence across a range of binary configurations. Our aim is to quantify how unstable mass-transfer dynamics in binaries at various evolutionary stages shape CE outcomes, enabling a predictive framework for modeling the LRN luminosity function. We find a bimodal distribution of plateau luminosities with significant implications for binary mass stability criteria that can be tested with forthcoming LSST observations. This bimodality emerges from differing mass-ejection outcomes during common-envelope interactions, which can lead either to stellar mergers, often accompanied by tidal disruption of the companion, or to successful envelope ejection. Although our predicted plateau luminosities and timescales broadly match existing observations, the models underpredict the number of LRNe with long-duration plateaus ($t_p \gtrsim 100\, \text{d}$) by about a third. We propose that these long-duration events arise from highly extended progenitors whose envelopes are ejected over multiple orbits (i.e., non-impulsively), producing relatively faint, long-lived transients. By constraining ejecta properties and incorporating pre-outburst progenitor imaging, we show how our models can clarify the physical processes that drive unstable mass transfer in these events. Finally, we argue that common-envelope interactions involving white-dwarf accretors can yield exotic outcomes, including red giants containing embedded white dwarfs that resemble Thorne-Żytków objects (TŻOs), along with calcium-rich supernovae that preserve hydrogen envelopes.

'Spider' binary systems - black widow and redback compact binaries differentiated by their companion's mass and nature - are an important type of pulsar system exhibiting a rich empirical phenomenology, including radio eclipses, optical light curves from a heated companion, as well as non-thermal X-ray and GeV orbital light curves and spectra. Multi-wavelength observations have now resulted in the detection of >~50 of these systems in which a millisecond pulsar heats and ablates its low-mass companion via its intense pulsar wind. Broadband observations have established the presence of relativistic leptons that have been accelerated in the pulsar magnetosphere and near the intrabinary shock, as well as a hot companion, presenting an ideal environment for the creation of orbitally-modulated inverse Compton fluxes that should be within reach of current and future Cherenkov telescopes. We have included an updated synchrotron kernel, different parametric injection spectral shapes, and several intrabinary shock geometries in our emission code to improve our predictions of the expected TeV signatures from spider binaries. Our updated phase-dependent spectral and energy-dependent light curve outputs may aid in constraining particle energetics, wind properties, shock geometry, and system inclination of several spider binaries.

The magnetic field of the Sun drives a wide range of eruptive phenomena, from small-scale nanoflares to large flares and coronal mass ejections (CMEs). While direct observations of solar activity cover only the past few decades, indirect evidence indicates that the Sun can occasionally produce events orders of magnitude stronger than any recorded ones in the modern era. Two complementary lines of evidence exist. First, extreme solar particle events (ESPEs) have been inferred from prominent spikes in cosmogenic isotope concentrations preserved in precisely dated natural archives such as tree rings and ice cores over the past 15 millennia. Second, high-precision space-borne photometry has revealed superflares on thousands of stars similar to the Sun. Whether these solar and stellar extremes are physically related remains an open question. We summarise the present state of understanding and discuss physical mechanisms that may link them. Although superflares and ESPEs are both extremely energetic manifestations of magnetic energy storage and release, their relationship does not appear to be one-to-one. Their occurrence and energetics likely depend on how magnetic flux and topology govern the partitioning of released energy between radiation, mass ejection, and particle acceleration.

Observational extragalactic catalogs over wide sky areas are essential for uncovering the large-scale structure of the Universe. They allow, among others, cosmological studies and density analyses that impose strong constraints on models of galaxy formation and evolution. By taking advantage of the wide field images and the 12 optical bands of the Southern Photometric Local Universe Survey (S-PLUS), we aim at providing a catalog of galaxies located, in projection, towards the Fornax galaxy cluster, within $\sim$ 5 virial radii in right ascension (R.A.) and $\sim$ 3 virial radius in declination (Dec) around NGC,1399, the dominant galaxy of the cluster. Such a catalog will allow unprecedented large-scale structure studies in that sky region. Supervised deep learning algorithms have been developed, utilizing neural networks complemented with dimensionality reduction techniques, to classify and separate spurious objects, stars and galaxies in a photometric catalog previously built for the S-PLUS Fornax Project (S+FP). That catalog was built using a combination of SExtractor configurations optimized for galaxy detection and characterization. A catalog of 119,580 galaxies was obtained in the direction of the Fornax cluster containing photometric information in the 12 optical bands of S-PLUS complemented with GALEX (UV), VHS-VISTA (NIR) and AllWISE (MIR) data. We estimate photometric redshifts ({\sigma}_ NMAD $\sim$ 0.0219) with a lower limit of z_ lim $\sim$ 0.03. Stellar masses, star formation rates (SFRs) and D4000_N index estimates were obtained through a machine learning approach, by matching S-PLUS photometric data to SDSS spectroscopic data. The completeness of the catalog (72%) was calculated by comparing with mock catalogs ...

As a byproduct of sub-Neptune formation, planetary embryos with high eccentricity can remain in outer orbits, near 1 au from the star. In this work, we investigate the long-term evolution of systems consisting of close-in sub-Neptunes (SNs) and outer high-eccentricity embryos. Our analysis focuses on collisions between SNs and embryos, particularly their atmospheric mass loss. We performed N-body simulations for various initial eccentricities and numbers of embryos. We analyzed the impact-induced atmospheric loss using post-processing methods, finding that the embryos and SNs collide at high speeds on timescales of several million years, leading to the loss of the SNs' atmospheres. Depending on the embryos' eccentricity and the orbital radius of the SNs, the impact velocity can be quite high, ranging from 2 to 5 times the escape velocity. On average, about 15%-30% of the atmosphere is dissipated per collision, so after 3-6 collisions, the atmospheric mass of an SN is reduced to about 1/3 of its initial value. Collisions between SNs and embryos can thus explain the presence of planets within the radius gap. Depending upon the initial eccentricity and the number of remaining embryos, additional collisions can occur, potentially accounting for the formation of the radius gap. This study also indicates that collisions between remaining embryos and SNs may help to explain the observed rarity of SNs with atmospheric mass fractions greater than 10%, commonly termed the "radius cliff."

Discovered almost 10 years ago, the giant planet 51 Eridani b is one of the least separated (0.2 arcsec) and faintest (J = 19.74 mag) directly imaged exoplanets known to date. Its atmospheric properties have been thoroughly investigated through low- and medium-resolution spectroscopic observations, enabling robust characterization of the planet's bulk parameters. However, the planet's intrinsically high contrast renders high-resolution spectroscopic observations difficult, despite their potential to yield key measurements essential for a more comprehensive characterization. This study seeks to constrain the planet's radial velocity, enabling a full 3D orbital solution when integrated with previous measurements. We have obtained 4 high-contrast high-resolution (R = 140000) spectroscopic datasets of the planet, collected over a two-year interval with the HiRISE visitor instrument at the VLT to derive the planet's radial velocity. Using self-consistent models of atmosphere, we were able to derive the radial velocity of the planet at each of the 4 epochs. These radial velocity measurements were then used in combination with all existing relative astrometry in order to constrain the orbit of the planet. Our radial velocity measurements allow us to break the degeneracy along the line of sight, making it now possible the unambiguous interpretation of the phase curve of the companion. We further constrain the orbital parameters, particularly the eccentricity, for which we derive e = 0.55 (-0.07, +0.03). The relatively high eccentricity indicates that the system has experienced dynamical interactions induced by an external perturber. We place constraints on the mass and semi-major axis of a hypothetical, unseen outer planet capable of producing the observed high eccentricities.

A recent paper on SOL2012-06-03 reported the detection for the first time of two distinct phases of $>$100 MeV $\gamma$-radiation indicating separate acceleration processes. But such two-phase emission has been seen before and was first observed in SOL1982-06-03. The second phase is known as Late Phase Gamma-Ray Emission (LPGRE) and was cataloged for $>$40 solar eruptions, including SOL2012-06-03. Here we provide evidence that the second SOL2012-06-03 $\pi$-decay peak is the onset of LPGRE that lasted for $>$8 min. Its delay from the impulsive X-ray peak is consistent with the time it would take flare-produced suprathermal protons to overtake the expanding CME and be accelerated by its shock. The high accelerated ion-to-electron ratio in SOL2012-06-03 and other LPGRE events is consistent with the ratio observed in gradual SEP events produced by shocks and is inconsistent with ratios typically found in impulsive flares and solar energetic particle events produced by reconnection.

Carbon is an essential element for a habitable world. Inner (r < 3 au) disk planetary carbon compositions are strongly influenced by supply and survival of carbonaceous solids. Here we trace the journey of carbon from the interstellar medium to the processes leading to planet formation. The review highlights the following central aspects: -Organics forming in evolved star envelopes are supplemented by aromatic molecules forming in the dense ISM to represent the seeds of (hydro)carbon supply through pervasive pebble drift to rocky planets and sub-Neptune cores. -Within the protoplanetary disk the sharp gradient in the C/Si content of Solar System bodies and mineral geochemistry outlines a tale of carbon loss from pebbles to within planetesimals and planets, and from planetary atmospheres. -Within two planet formation paradigms (pebble and planetesimal accretion) a range of planetary carbon content is possible that is strongly influenced by early (< 0.5 Myr) formation of a pressure bump that titrates drift. Overall, it is unlikely that the carbon architecture of our Solar System applies to all systems. In the absence of giant planets, carbon-rich rocky worlds and sub-Neptunes may be common. We outline observations that support their presence and discuss habitability of terrestrial worlds.

Context: The Canis Major (CMa) star-forming region, a remote molecular cloud complex within the recently discovered Radcliffe Wave, remains under-explored in the literature. Aims: We revisit the stellar census in the CMa region, characterizing its stellar population, kinematics, and age using recent astrometric and photometric data from the third data release of the Gaia space mission (Gaia DR3). Methods: We conducted a membership analysis of Gaia DR3 sources across a 16 deg$^2$ field encompassing the youngest subgroups in CMa. This new stellar census, combined with spectroscopic observations, allowed us to investigate the structure, kinematics, and age of this region. Results: We identified 1531 objects as members of the CMa region, confirming 401 previously known members and introducing 1130 new candidate members. These objects have magnitudes ranging from 10 to 18 mag in the G band from Gaia DR3. We identified two subgroups of CMa stars in our sample labelled as Cluster A and Cluster B. They are located at roughly the same distance ($d_{A} = 1150^{+79}_{-88}$ pc and $d_{B} = 1183^{+103}_{-108}$ pc) and exhibit similar space motions that can be derived thanks to the precise radial velocities obtained in this study. The subgroups have a mean isochronal age of about 2-3 Myr. However, based on infrared photometry we show that Cluster A has a higher fraction of disc-bearing stars suggesting that it could be somewhat younger than Cluster B. Conclusions: Our analysis provides new insights into the stellar population of the Canis Major region, by identifying new members, characterizing their kinematics, and assessing their evolutionary stages. Future studies incorporating additional data from upcoming Gaia data releases, multi-wavelength and high-resolution spectroscopic observations will be essential to further advance our understanding of the history of star formation in this region.

Recent large-scale multi-band photometric surveys now enable elemental-abundance estimates for millions of stars with accuracies approaching those of low- to medium-resolution spectroscopy. Using [Fe/H] and [C/Fe] estimates derived from the Javalambre Photometric Local Universe Survey (J-PLUS) DR3 and the Southern Photometric Local Universe Survey (S-PLUS) DR4, which together cover $\sim$6,200 deg$^2$ of the sky, we identify large numbers of carbon-enhanced metal-poor (CEMP) stars in the Milky Way. After applying data-quality cuts and evolutionary corrections to the carbon-abundance estimates, we construct a combined J/S-PLUS sample of $\sim$6.40 million stars and identify $\sim$104,900 CEMP candidates, roughly twice the number of CEMP candidates identified from Gaia XP spectra by Lucey et al. We photometrically confirm that the absolute carbon abundance $A$(C) separates CEMP stars into two primary groups, CEMP-no and CEMP-$s$ stars, consistent with previous spectroscopic studies. We also recover CEMP morphological Groups I-III in the Yoon-Beers diagram, as well as the recently proposed Group IV, and show that it is statistically distinct even in photometric data. A cumulative frequency analysis confirms that the CEMP fraction increases toward lower metallicity and that CEMP-no stars dominate in the most metal-poor regime. By comparing frequencies with and without Group IV stars, we assess their relation to CEMP-no and CEMP-$s$ stars, and examine CEMP distributions across different Galactic components. The resulting catalog provides a substantial sample for future spectroscopic follow-up, in particular to constrain the likely origin(s) of the Group IV stars.

Fast radio bursts (FRBs) emit broad band radio wave radiation that may, in rare cases, encode atomic hydrogen (HI) absorption signals produced as they traverse the interstellar medium of their host galaxies. Combining such signals with high resolution HI emission maps offers a unique opportunity to probe the dynamics of neutral gas at cosmological distances through constraints of the HI excitation temperature $T_{spin}$, which characterises the balance of neutral gas phases and the underlying thermal processes within these galactic environments. While no absorption signal has been recorded in an FRB to date, we demonstrate a proof of concept with the bright (F = 35 Jy ms) and narrow (0.2 ms) FRB 20211127I detected by ASKAP. We find a 3${\sigma}$ upper limit on the integrated optical depth in the pulse-averaged spectrum of 33 km s$^{-1}$, and, based on the HI emission observed in a 3 hr MeerKAT L-band observation, subsequently find a lower limit on $T_{spin}$ of 26 K. While this test case provides little constraining power, we find that narrow, non-repeating FRBs with fluences greater than 20/70/150 Jy ms observed with all dishes with the current MeerKAT/ASKAP/DSA telescopes can probe integrated optical depths below 5 km s$^{-1}$. Furthermore, we highlight that utilising FAST's incredible sensitivity to stack thousands of bursts from hyperactive repeaters also provides a plausible avenue through which HI absorption, and hence $T_{spin}$, can be measured. Finally, we discuss how HI absorption can address several modern challenges in FRB science, providing a physical anchor for locating bursts within their host galaxies and helping to disentangle the host contribution to dispersion and scattering.

Context: JWST has enabled transmission spectroscopy at unprecedented precision, but stellar heterogeneities (spots and faculae) remain a dominant contamination source that can bias atmospheric retrievals if uncorrected. Aims: We present a fast, unsupervised methodology to reduce stellar contamination and instrument-specific noise in exoplanet transmission spectra using denoising autoencoders, improving the reliability of retrieved atmospheric parameters. Methods: We design and train denoising autoencoder architectures on large synthetic datasets of terrestrial (TRAPPIST-1e analogues) and sub-Neptune (K2-18b analogues) planets. Reconstruction quality is evaluated with the $\chi^2$ statistic over a wide range of signal-to-noise ratios, and atmospheric retrieval experiments on contaminated spectra are used to compare against standard correction approaches in accuracy and computational cost. Results: The autoencoders reconstruct uncontaminated spectra while preserving key molecular features, even at low S/N. In retrieval tests, pre-processing with denoising autoencoders reduces bias in inferred abundances relative to uncorrected baselines and matches the accuracy of simultaneous stellar-contamination fitting while reducing computational time by a factor of three to six. Conclusions: Denoising autoencoders provide an efficient alternative to conventional correction strategies and are promising components of future atmospheric characterization pipelines for both rocky and gaseous exoplanets.

With the improvements in radio interferometry sensitivity, the quest for coherent radio emission from exoplanets and ultra-cool dwarfs, which is indicative of their magnetic fields, has gained significant momentum in recent years. We investigated the relatively unexplored possibility of radio emission from wide-orbit brown dwarf companions, which may radiate through rapid rotation, as in isolated ultra-cool dwarfs, or via interactions between their extended magnetospheres and the host star's wind. We analysed $\sim 60$ hours of Upgraded Giant Metrewave Radio Telescope and Karl G. Jansky Very Large Array data for a set of well-characterized systems previously unobserved at 0.3-2 GHz. The targets include companions orbiting the G-type stars HD 26161 and BD-004475, the K-type HD 153557A and $\nu$ Oph, and the M dwarfs GJ 3626 and 2MJ01225093-2439505. No detections were obtained with 3$\sigma$ upper limits down to $\sim 25\,\mu$Jy/beam in Stokes V in the best cases. The light-curve analysis also revealed no evidence of short ($\gtrsim$ minutes), intense ($\gtrsim$ mJy) radio bursts. The upper limits provide tentative constraints on model parameters. However, the effects of model uncertainties, limited observational coverage, and intrinsic variability or beaming of the emission must be considered. The improvement in sensitivity of the next-generation radio interferometers will likely allow to go below the expected flux range over a much larger range of free parameters.

We investigate the total mass density slope evolution in massive early-type galaxies (ETGs) over the last 6.5 billion years ($0 < z < 0.75$). We perform a detailed dynamical analysis of approximately 200 ETGs spanning the redshift range $0.24 < z < 0.75$, utilizing spatially resolved stellar kinematics derived from high signal-to-noise ratio (S/N) MUSE-DEEP spectroscopy and surface brightness models from high-resolution HST imaging. We constrain mass distributions using the Jeans Anisotropic Modeling (JAM) technique coupled with Multi-Gaussian Expansion (MGE) method. To rigorously constrain evolutionary trends, we combine this intermediate-redshift dataset with a local ETG sample ($z \sim 0.05$) from the MaNGA survey. We adopt dynamical constraints for the local sample derived using an identical homogeneous methodology, ensuring a strictly consistent comparison. We found that the total density profiles of the intermediate-redshift ETG sample are approximately isothermal and exhibit a median mass-weighted total density slope, $<\gamma_{\rm T}>=2.19 \pm 0.01$ at $<z>=0.44$, which is shallower than the local baseline of $<\gamma_{\rm T}> = 2.26 \pm 0.01$ at $<z>=0.04$. This structural shift corresponds to a redshift gradient of $\mathrm{d} \gamma_{\rm T}/\mathrm{d} z \approx -0.20 \pm 0.03$, detected at $\sim$5-$\sigma$ significance. We demonstrate that this trend is robust against model assumptions and persists even when restricting the analysis to high-velocity dispersion systems ($\sigma_e > 150$ km/s). Our findings are consistent with previous lensing-based studies and in tension with cosmological simulations. The observed steepening suggests that dissipative processes, such as gas-rich accretion and mergers, must play a non-negligible role in the late-stage assembly of massive ETGs.

Earth's climate is influenced by over a dozen feedbacks, but only three dominate its long-term climate behavior. Models of the exoplanet habitable zone (HZ) assume that this is similar for other Earth-like planets. We used dynamical simulations to study Earth-like planets with a fourth, (potentially strong) generalized climate feedback. Across over 20,000 climate simulations, we find that the addition of the fourth feedback produces novel behaviors, including runaway and chaotic climate trajectories, that are more diverse than one would expect based on Earth's climate configuration. Non-negligible fourth feedbacks -- if negative -- would not lessen the probability of planets with temperate climates. However, positive fourth feedbacks decrease the fraction of exo-Earth candidates that are long-term habitable. Therefore, strong fourth feedbacks will alter (and mostly shrink) the boundaries of the classical habitable zone. When combined with occurrence rates of Earth-sized planets around sun-like stars, our results imply that the fraction of stars hosting rocky planets with temperate climates may be substantially lower than classical estimates under Earth-like climate assumptions. Our results are subject to the validity of the model assumptions and not intended to represent conclusive predictions about exoplanet populations but rather to demonstrate the potential climate diversity that emerges from non-Earth-like model configurations. Our conclusions provide context on sample sizes and science questions for next-generation exoplanet surveys.

Coronal mass ejections (CMEs) are the primary drivers of adverse space-weather events, yet their initiation and onset prediction remain insufficiently understood due to the complexity of the magnetic topology and physical processes in real solar source regions. Here, based on fully observational-data-driven magnetohydrodynamic simulation, we successfully reproduce the initiation of a CME originating from the super active region AR 13663, with only a one-minute time lag between the flare peak in observations and the velocity peak of the rising flux rope in the simulation. Moreover, the eruptive structure exhibits a multi-stage kinematic evolution: an initial slow acceleration, a plateau at a nearly stationary height, and a subsequent impulsive acceleration. These stages correspond to torus instability, the downward tension force exerted by the overlying toroidal field, and fast magnetic reconnection, respectively. Our results highlight the inherently multistage nature of CME initiation in real events. In configurations with strong overlying toroidal fields, the downward toroidal-field-induced tension force can suppress the rise of the flux rope and produce a plateau phase at a nearly stable height, even when torus instability occurs. In contrast, the subsequent fast magnetic reconnection beneath the flux rope can drive the impulsive eruption more effectively. The close agreement between the observed and simulated peak times over one minute demonstrates the strong potential of our data-driven model for predicting CME onset.

Symbiotic stars are interacting binaries consisting of a red giant and typically a white dwarf, important as potential Type Ia supernova progenitors. Despite theoretical predictions that white dwarfs in symbiotic systems should often be magnetic, direct evidence has been elusive. We report the discovery of a $682.5 \pm 7$ s periodicity in the XMM-Newton X-ray light curve that we interpret as the spin period of the WD in CH Cygni. This detection provides strong evidence for a magnetic white dwarf in CH Cygni, making it only the second WD symbiotic star with confirmed X-ray pulsations after R Aquarii. Our discovery supports the magnetic propeller model previously proposed for CH Cygni's jet activity and state transitions.

We report the first detection of $\gamma$-ray emission from the young radio galaxy JVAS~J1311+1658, classified as a compact symmetric object (CSO). This detection is characterized by a recent GeV $\gamma$-ray flare identified in Fermi-LAT data during MJD~60032.6--60132.6, with a $\gamma$-ray source detected at a significance level of $\sim6.2\sigma$. The average 0.1--300~GeV flux is measured to be $(1.6 \pm 0.6)\times10^{-8}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$, with a photon spectral index of $\Gamma = 2.15 \pm 0.185$. We find that a radiative model of the radio lobes significantly underestimates the observed $\gamma$-ray emission. The strong flux and short-term variability over $\sim$100 days suggest that the emission likely originates from newly launched sub-kiloparsec-scale jets at the core. This detection provides a unique window into the extreme environments and early-stage jet activity of young radio galaxies, offering insights into their initial evolution and the formation of relativistic jets in the earliest phases of galaxy growth.

Star-forming edge-on galaxies often exhibit extended halo radiation in multiple bands, providing ideal laboratories for studying the transfer of matter from the disk to the halo. We investigate the transport of cosmic-ray electrons (CREs) and the associated galactic wind, and assess their impact on the surrounding medium in NGC 3044. We obtained the NGC 3044 total intensity image at 943 MHz from the Australian SKA Pathfinder (ASKAP) observations with a resolution of 16 arcsec and an rms noise of 20 $\mu$Jy beam$^{-1}$. The sensitivity is higher than the previous observations at similar frequencies. We find that the ASKAP intensity profiles perpendicular to the disk can be fit with two exponential components. The scale heights of the thin and thick disks are $0.43 \pm 0.13$ kpc and $1.91 \pm 0.26$ kpc, respectively. By jointly fitting total intensity and spectral index profiles with one-dimensional advection and diffusion models, we find that CREs are advected outward from the disk with the velocity increasing with height in a power law. Beyond $\sim3$ kpc, the velocity exceeds the escape speed of $\sim400$ km s$^{-1}$, indicating a strong wind. We further identify a possible superbubble of radius $\sim3$ kpc filled with soft X-ray emitting hot gas and surrounded by an HI shell and a bright H$\alpha$ rim. These results demonstrate that radio continuum observations provide a powerful probe of cosmic-ray-driven winds in normal star-forming spiral galaxies.

Based on the classical unification, blazars, namely BL Lacertae objects (BL Lacs) and flat-spectrum radio quasars (FSRQs), are believed to correspond with radio galaxies when observed at small jet viewing angles. In this paper, we aim to compile a sample of Fermi blazars and radio galaxies to provide new insights towards a unified accretion and ejection scenario between aligned and misaligned radio-loud active galactic nuclei (AGNs), by considering their optical emission-line classifications (low- and high-excitation radio galaxies, LERGs, HERGs), which are more representative of their accretion states. We adopted statistical analyses of accretion properties and high-energy beaming patterns for both Fermi blazars and radio galaxies to investigate a unified accretion-ejection scenario. In the gamma-ray luminosity-photon index plane, HERGs populate the region of higher luminosities and softer photon indices, akin to FSRQs, whereas LERGs fill lower luminosities with harder photon indices, analogous to BL Lacs. This parallel segregation indicates that LERGs and HERGs represent the misaligned counterparts of BL Lacs and FSRQs, respectively. The unified picture is further supported by the Compton dominance-photon index diagram, where FSRQs and HERGs dominated by external Compton (EC) emissions are distinctly separated from BL Lacs and LERGs governed by synchrotron self-Compton (SSC) emissions. Similarly, the diagram of accretion rate versus gamma-ray photon index reveals two distinct accretion-ejection states: a low-accretion-rate branch (BL Lacs and LERGs) is associated with the SSC model, and a high-accretion-rate branch (FSRQs and HERGs) is linked to the EC model. These results strongly strengthen the idea of a unified accretion and ejection paradigm between blazars and radio galaxies into two distinct states.

The cosmic distance scale is built on multiple different techniques for estimating distances in space that are often connected and dependent on multiple measurements and assumptions. Double white dwarf binaries (DWDs) are common objects and are expected to produce gravitational wave (GW) signals that can be observed with space-based detectors such as LISA. By analyzing these signals we should be able to estimate the distance and sky location of the source. Previous studies have done this for circular binaries which, while they are abundant, have, in general, weaker signals than eccentric binaries and it is not possible to differentiate whether a circular binary is in the field or in a dense environment such as a globular cluster (GC). In this paper we used eccentric binaries from MOCCA GC simulations, simulated the GW signal from each binary at locations related to GCs in the Milky Way and estimated the precision on the distance and the sky location of the source. We find that distances can be estimated with higher precision than current day methods even with low eccentricity binaries and higher eccentricity further increases this precision. Although the probability of finding a tight and eccentric DWD is far lower than a circular one, we can expect to find at least a few in the dense environments of the Milky Way, such as GCs. These estimations would be independent measurements with high precision to objects inside dense environments, such as GCs inside the Milky Way and the Magellanic Clouds.

We investigate the effects of different dust and gas opacity descriptions on the structure and evolution of the inner regions of protoplanetary discs. The influence on the episodic instability of the inner rim is hereby of central interest. 2D axisymmetric radiation hydrodynamic models are employed to simulate the evolution of the inner disc over several thousand years. Our simulations greatly expand on previous models by implementing detailed opacity descriptions in terms of their mean and frequency-dependent values, allowing us to also consider binned frequency-dependent irradiation. The adaptive opacity description significantly affects the structure of the inner disc rim, with gas opacities exerting the greatest influence. The resulting effects include shifts in the position of both the dust sublimation front and the dead zone inner edge, a significantly altered temperature in the dust-free region and the manifestation of an equilibrium temperature degeneracy as a sharp temperature transition. The episodic instability due to MRI activation in the dead zone still occurs, but at lower inner disc densities. While the gas opacities set the initial conditions for the instability, the evolution of the outburst itself is mainly governed by the dust opacities. The analysis of criteria for non-axisymmetric instabilities reveals possible breaking of the density peaks produced by the burst. However, due to the periodicity of the instability, the inner edge itself may remain stable throughout quiescent phases according to linear criteria. Although the thermal structure of the inner disc is crucially affected by different opacity descriptions, the mechanism of the periodic instability of the DZIE remains active and is only marginally influenced by gas opacities. The observational consequences of the severely altered temperatures may be significant and require further investigation.

The linear structure reported by van Dokkum et al. (2023) has been proposed as either a massive stellar wake produced by a runaway supermassive black hole (SMBH) or a bulgeless edge-on galaxy. New JWST/NIRSpec IFU observations target the tip of the structure, where a SMBH would produce a bow shock, whereas a normal galaxy would host an HII region. Using standard BPT diagrams ([OIII]5007/Hb vs [NII]6583/Ha and [OIII]5007/Hb vs [OII]6716,6731/Ha), we find that the line ratios at the tip fall on the locus of low-metallicity low-extinction HII regions. This region does not overlap with loci typical of shocks in merging galaxies. Thus, these results are consistent with the interpretation that the linear structure is a star-forming galaxy, with the bright knot representing one of its HII regions.

High expansion velocity carbon stars (HVCs) are a rare class of evolved stars whose circumstellar envelopes (CSEs) combine C-rich chemistry with unusually high expansion velocities typical of O-rich massive evolved stars. AFGL2233 has been proposed as a high-mass evolved object that exhausted hot-bottom burning. Studying its chemistry is essential to understand the nature and evolution of these objects. We characterize the chemical composition and isotopic ratios of the CSE of AFGL2233 and investigate chemical peculiarities, including the presence of N- and O-bearing species in a C-rich environment. We carried out a complete line survey at 3 mm and 1 mm using the IRAM 30m telescope, complemented by Herschel/HIFI FIR observations and interferometric maps of SiO, C2H, and HCN. Molecular emission was analyzed using rotational diagrams and radiative transfer modeling under the LVG approximation. Column densities and fractional abundances were derived for more than 30 molecular species, including isotopologues, and compared with other evolved stars. The Gaia DR3 distance of 1.236 kpc implies a luminosity of ~2 Lsun, consistent with an initial mass of 4.5-9 Msun. The molecular inventory confirms C-rich chemistry but reveals unusually high abundances of NH3, H2O, and SiN. The isotopic ratios vary among species, with 12C/13C ranging from 7 to 55. The C2H/C4H ratio is abnormally high compared with C-rich AGB stars. The presence of SiN and high NH3 may indicate N-enrichment or the influence of a companion. AFGL2233 is likely a high-mass AGB or super-AGB star with a complex evolutionary history involving nucleosynthesis, shocks, and possible binary interaction.

In the search for metal-poor stars, large spectroscopic surveys are an invaluable tool. However, the spectra of metal-poor stars can be difficult to analyse because of the relative lack of available lines, which can also lead to misclassification. We aim to identify the stars observed by the APOGEE survey that are below the metallicity limit of APOGEE's analysis. For the highest confidence candidates, we classify the orbital properties of the stars to investigate whether their orbital distribution matches what we would expect for stars that are this metal poor and to select especially interesting targets for spectroscopic follow-up purposes. We examined the properties derived by APOGEE for metal-poor stars from the literature to find signatures of stars with a metallicity below the range of the grid used for spectral analysis. The calibrated APOGEE stellar parameters provide a clear signature of the most metal-poor stars in the survey, indicated by null values for their metallicities while having effective temperatures and surface gravities determined by the pipeline. From comparison with the literature, we find that, within a temperature range of 3700 - 6700 K and above a threshold of S/N > 30, the vast majority of APOGEE stars without calibrated metallicities are very metal poor. The radial velocities provided by APOGEE, Gaia DR3 positions and astrometry along with spectrophotometric distances derived in this work allowed for the computation of their orbits. In this work, we select 289 very metal-poor red giant stars (likely below = -2.5) from the APOGEE results. Our sample contains 16 very metal-poor member candidates of the Magellanic Clouds, 14 very metal-poor stars with orbits confined to the inner Galaxy, and 13 inner Galaxy halo interlopers.

H-band spectra contain very limited spectral information for stars at the most metal-poor tail ( Fe/H < -2.5) because the available Fe lines in FGK stars in this wavelength range are weak. The first paper in this series successfully identified a sample of 327 very metal-poor stars (with [Fe/H] < -2) from the APOGEE database, 289 of which are on the red giant branch. The spectra of these stars were not properly analysed by the APOGEE main pipeline because they are very metal poor. In this work, we measure metallicities for these stars using the abundances of the elements Mg and Si. We demonstrate that the absorption lines of the elements Mg and Si are of good quality despite the challenging combination of (low) metallicity, wavelength regime, spectral resolution, and signal-to-noise ratios available for these spectra. A specialised pipeline was designed to measure the abundance of Mg and Si in APOGEE spectra and yielded a robust estimate of the overall metallicity. In order to provide reliable measurements, we tested three different sets of assumptions for Mg and Si enhancement. We present Mg and Si abundances as well as overall metallicities for 327 stars, all of which had previously gotten null values from the main APOGEE pipeline for either the calibrated M/H or [Fe/H] . The typical uncertainties for our measurements are 0.2 dex. We found five stars in our sample with unusual [Si/Mg] abundances higher than 0.5, and we connect this signature to globular cluster stars, and this might be related to specific supernova events. Our data suggest a concentration of high [Si/Mg] stars in the Sextans dwarf galaxy. Other dwarf galaxies are found to agree well with results in the literature. Our derived metallicities range between -3.1 $\leq$ [M/H] $\leq$ -2.25, thereby pushing the metal-poor tail of APOGEE results down by 0.6 dex.

Beyond cosmic noon, galaxies appear as faint whispers amid noise, yet this epoch is key to understanding massive galaxy assembly. ALMA's sensitivity to cold dust and [C II] emission allows us to probe their interstellar medium, but faint signals make robust denoising essential. We evaluate and benchmark denoising strategies including Principal Component Analysis, Independent Component Analysis, sparse unsupervised representations: iterative soft thresholding with 2D-1D wavelets, and supervised deep learning with a 3D U-Net, to identify techniques that suppress noise while preserving flux and morphology across peak SNRs of 2.5-8, applied to (i) synthetic spectral cubes of rotating toy disk galaxies, (ii) synthetic [C II] IFU cubes from FIRE simulations, and (iii) ALMA [C II] observations of CRISTAL galaxies and W2246-0526. Performance is assessed via RMSE, conservation of flux and spectra, noise reduction, and SNR improvement of the central galaxy. For synthetic cubes: PCA and ICA provide marginal improvement; IST reduces noise effectively at moderate SNRs but can suppress emission at low SNRs; and the U-Net outperforms IST, though it can produce quantifiable hallucinations at lower-SNRs. For moderate-SNR observations (ALMA-CRISTAL), U-Net and IST achieve comparable performance, conserving >91% flux and increasing SNR by >6. However, for observations with complex morphologies absent in the training set (W2246), the U-Net underperforms relative to IST, recovering ~80% flux, while IST robustly conserves flux and improves SNR by ~3, highlighting generalisation challenges and the need for physically-motivated training priors. We conclude that IST is a robust unsupervised denoiser for moderate-SNR data, and a synthetically trained U-Net generalises effectively to real data, dependent on training priors. This framework offers a pathway for transferable denoising for ALMA, VLT/MUSE, and JWST.

Pollux is a high-resolution spectropolarimeter proposed by an European consortium for HWO. The current design of Pollux features four spectropolarimetric channels, three of which are in the UV range. For the near-UV (NUV) [236-472 nm] and mid-UV (MUV) [118-236 nm] channels, the polarimeters consist of waveplates and prisms made of MgF2, a birefringent material. However, no such birefringent material can be used for the far-UV (FUV) channel [100-123 nm]. Therefore, the polarimeter for this FUV channel is composed solely of mirrors in an innovative assembly. In this talk, we aim to detail the architecture of the test bench that will allow us to validate the performance of these different polarimeters, as part of the HWO GOMaP. Given that we are working in the vacuum ultraviolet (VUV) range, the test bench operates in a vacuum chamber in a clean room. We will discuss the adaptable architecture of the bench based on wavelength and the measurement methodology that we will implement to test if the polarimeters achieve the precision of $10^{-3}$ required for the Pollux instrument. With this test bench, we will successfully increase the Technology Readiness Level (TRL) of UV spectropolarimeters and, for the first time, develop a means to test FUV spectropolarimetry.

Thermal instability plays a crucial role in the dynamics of astrophysical plasmas. Building upon the foundational work of Field (1965) and the subsequent analysis by Waters and Proga (2019), this study revisits the characteristics of thermal instability in a uniform, non-magnetic medium. Our primary aim is to reevaluate and expand the understanding of instability domains, focusing on the classification and characteristics of thermal and acoustic modes in the presence of heating, radiative cooling, and thermal conduction. Except for Spitzers expression for parallel thermal conductivity, the heating and cooling processes are unspecified. Additionally, we investigate the existence of isobaric and isochoric thermal modes across the extreme limits of very short and very long wavelengths, as well as at intermediate wavelengths. We perform an in-depth analysis of the dispersion relation for an infinite, uniform hydrodynamic medium, as initially derived by Field (1965). This approach enables the generation of growth-rate and dispersion diagrams, providing insight into the behaviour of thermal instability across different wavelength ranges. With the inclusion of thermal conduction, our study refines the classification of the seven instability regions previously outlined by Waters and Proga (2019). Our findings confirm that their classification holds when the Field length is smaller than or comparable to the thermal wavelength. For larger Field lengths, a simplified classification becomes impractical. Furthermore, we discuss the potential implications of the catastrophic cooling instability (Waters and Stricklan 2025) in the context of cool coronal rain formation. Our work not only validates certain classifications, but also proposes new insights into complex instability behaviours, thereby enhancing the theoretical framework established by Field and others.

A robust stellar initial mass function (IMF) is crucial in any studies related to star formation. However, the direct measurement of the stellar IMF is confined to the local universe, limited by the resolving power of telescopes. Recently, a new method for accessing the stellar IMF beyond the local universe has been developed. The observed detection rate of transient lensed stars -- individual, massive, thus luminous stars in strongly lensed galaxies that are temporarily detectable upon stellar microlensing -- can serve as a probe to break the IMF-star formation history degeneracy in studies utilizing spectral energy distribution fitting, hence providing a window to look at the IMF at a subsample of gravitationally lensed galaxies. In this proceeding, I summarize the contributed talk given at IAUS402 entitled the same as this contribution and highlight some key results, which currently show no evidence for a top-heavy IMF in $z \approx 1$ galaxies.

We report the discovery of a new massive giant planet, HD 164604 c ($a_c = 5.556_{-0.10}^{+0.093}$ au, $e_c = 0.196_{-0.078}^{+0.078}$ and $m_c = 9.5_{-1.25}^{+1.2}$ or $7.6_{-1.0}^{+1.0}\,M_{\rm Jup}$), orbiting a K3.5 dwarf, The result is based on the combined analysis of high-precision radial-velocity data, Hipparcos, and Gaia DR2 and DR3 astrometry. We refine the orbital parameters of the inner planet HD 164604 b to $a_b = 1.362_{-0.012}^{+0.012}$ au, $e_b = 0.479_{-0.021}^{+0.027}$, and $m_b = 13.2_{-1.5}^{+1.8}\,M_{\rm Jup}$ (or $8.8_{-1.5}^{+1.9}\,M_{\rm Jup}$). Depending on the two possible orbital orientations of HD 164604 c, the true mutual inclination between the two planets is $\psi_{bc}=5.0^{+3.7}_{-2.2}$$^\circ$ (prograde) or $162.1^{+7.1}_{-4.7}$$^\circ$ (retrograde). Long-term N-body integrations show that most orbits with the retrograde configuration remain dynamically stable for at least 10 Myr, while orbits with the prograde motion might rapidly evolve into chaos or lead to ejection. The retrograde architecture points to a violent dynamical history, possibly involving von Zeipel-Lidov-Kozai cycles or three-body scattering, while the prograde scenario might be consistent with coplanar and mild disk migration. Future Gaia DR4 astrometry will break the inclination degeneracy and distinguish between prograde and retrograde orbits for HD 164604 c.

The dense Galactic Center environment is expected to host compact binary inspirals detectable by future space-borne gravitational wave (GW) observatories (e.g., LISA, TianQin, Taiji) in the millihertz band. Aided by information from these facilities, next-generation ground-based GW detectors (e.g., Cosmic Explorer, Einstein Telescope) can potentially capture gravitational radiation in the hectohertz band from rapidly spinning neutron star (NS) components in such binaries. These Galactic Center systems are thus anticipated to act as dual-line (i.e., low-frequency inspiral and high-frequency spin) GW sources. However, the formation channels of these systems remain largely unexplored. In this \textit{Letter}, we propose that the von Zeipel-Lidov-Kozai (ZLK) effect can enhance the formation of dual-line GW sources in hierarchical triples involving the Galactic supermassive black hole, Sgr A*. We show that ZLK-driven oscillations in the eccentricity and inclination of the inner binary can modulate the GW emission from both the binary inspiral and the individual NS spins. This effect boosts the expected dual-line source count by a factor of $\sim 3$, from rare to $\mathcal{O}(1)$ in 4 years, making dual-line observations substantially more probable. Our results demonstrate that the ZLK effect provides an important formation channel for Galactic dual-line GW sources.

Symbiotic stars (SySts) are interacting binaries composed of a red giant transferring material to a hot compact star, typically a white dwarf. Although only about 300 systems are confirmed, the Galactic population is estimated at 1.2 x 10^3 - 1.5 x 10^4, indicating that most remain undiscovered. We identify new SySts using a machine-learning approach that combines Gaia DR3, 2MASS, and WISE photometry, parallaxes, and the pseudo-equivalent width of H alpha. A Random Forest model was trained on 166 confirmed S-type SySts and 1600 non-symbiotic stars, applying SMOTE to mitigate class imbalance. The model achieved an F1-score of 89% for the symbiotic class. Applied to 2.5 x 10^6 color-selected sources, it identified 990 candidates with probabilities more than 70%. We further refined the sample using physically motivated cuts on effective temperature, surface gravity, metallicity, and SkyMapper photometry, yielding 12 high-confidence candidates. These objects show cool temperatures, low surface gravities, near-solar metallicity, H alpha emission, moderate-to-high luminosities, and UV excess consistent with S-type SySts. Validation on recently confirmed systems recovered 92.3%, demonstrating the robustness and generalizability of our method.

The X-ray output of low-metallicity starburst galaxies is a key component of stellar feedback, tracing the processes responsible for gas ionization and chemical enrichment. The integrated X-ray luminosity ($L_X$) from high-mass X-ray binaries in star-forming galaxies scales with both the star formation rate (SFR) and host-galaxy metallicity $Z$. Due to the inverse correlation between $L_X/\mathrm{SFR}$ and $Z$, the contribution of X-ray binaries to the ionizing photon budget is enhanced in metal-poor systems and may ionize He II in the surrounding interstellar medium, powering nebular He II $\lambda4686$ emission. The blue compact dwarf galaxy ESO 338-IG04 (ESO 338-4) provides a nearby laboratory for studying stellar feedback in a low-metallicity starburst, combining vigorous recent star formation, low metallicity ($12+\log(\mathrm{O/H})\approx7.9$), and a rich population of massive stellar clusters. We characterize the X-ray emission of ESO 338-4 and its halo using new deep Chandra and XMM-Newton observations. We analyze X-ray spectra, light curves, and images to constrain the nature of its X-ray sources. We identify five ultra-luminous X-ray sources (ULXs) and diffuse hot gas surrounding the galaxy. Two ULXs are spatially associated with stellar clusters. The total X-ray luminosity exceeds $10^{41}\,\mathrm{erg\,s^{-1}}$. The brightest source, ULX1, shows variability on day timescales and lacks a stellar-cluster counterpart. Photoionization modeling shows that X-ray sources significantly impact the ionizing photon budget; models with ULX1 as the ionizing source predict nebular He II $\lambda4686$ luminosities of $\sim10^{39}\,\mathrm{erg\,s^{-1}}$.

Recent cosmological tensions pose difficulties for $\Lambda$CDM. Forthcoming facilities will be able to check whether these tensions result from systematic effects or indeed with the $\Lambda$CDM model itself. However, these new data will primarily probe gravitational interactions and provide only limited information about non-gravitational interactions. Distinguishing between competing models that make similar predictions yet rely on fundamentally different principles, therefore requires suitably diverse physical tests. Observational constraints on spacetime variations of fundamental constants fill this need. The fine-structure constant, $\alpha = e^2/\hbar c$, can be measured using absorption systems towards bright quasars using the Many Multiplet method, and using atomic doublets from line emitting gas in galaxies. A spectroscopic facility such as the WST could produce more than 100,000 new measurements of $\alpha$ from quasars together with a million measurements from galaxies. When combined with other probes, such a large and homogeneous dataset of $\alpha$ measurements would provide unprecedented constraints on physics beyond $\Lambda$CDM.

The correlation between inner super-Earths (SEs) and outer cold Jupiters (CJs) provides an important constraint on the formation and dynamical evolution of planetary systems. Previous studies have suggested a positive connection between these two populations, particularly around metal-rich stars, and proposed that nearly all CJ-hosting stars may also harbor inner SEs. In this work, we use TESS transits to investigate the occurrence of hot SEs in systems with known CJs detected by radial velocity (RV). Out of a statistical sample of 132 CJ systems, we identify five transiting hot super-Earths ($1$-$4R_\oplus$, $P<10\mathrm{d}$) around four stars, including one new candidate (TOI-6965.01) around HD 50554. To enable statistical analysis, we first validate the two candidates around HD 50554 using TESS photometry, archival RV measurements, and Gaia astrometry. After accounting for detection sensitivity and geometric transit probability, we find that the presence of CJs enhances the occurrence rate of hot super-Earths by a factor of $8.1^{+4.3}_{-3.2}$ relative to field stars, with the case of no enhancement ruled out at the 99.9% confidence level. Taking into account the average multiplicity of hot SEs, we find that about 87% of CJ systems host at least one hot SE. Our results provide strong supporting evidence for a positive SE-CJ correlation. We also briefly explore the correlation around metal-poor hosts and for specific sub-populations (e.g., warm super-Earths or cold super-Jupiters).

The improved sensitivity of Gravitational-Wave detectors and the development of eccentric waveform models enable us to explore the growing catalog of gravitational-wave events with measurable eccentricity. This opens new opportunities to gain insight into the formation channels and evolutionary pathways of compact binary systems using eccentricity. However, most recent population analyses have been limited to quasi-circular binaries, primarily due to constraints in waveform modeling and sensitivity estimates. We are now entering an era where both of these limitations are being addressed, allowing for a more comprehensive investigation of eccentric binary populations. In this work, we perform the first population inference analysis that simultaneously fits the mass, spin, redshift, and eccentricity distribution. Specifically, we use source-parameter estimation provided by the Rapid Iterative FiTting (RIFT) framework using the SEOBNRv5EHM waveform model, and a default O4a population model extended to include eccentricity. We find population properties broadly consistent with conclusions obtained in previous analyses assuming quasi-circular binaries. Consistent with our conclusions about each event, we bound the branching ratio for eccentric events to be below $0.0239$ at $90\%$ confidence with our fiducial eccentricity mixture models. Using four different parametric population models for eccentricity, we point out that the rate of eccentric events is weakly constrained by observations and highly model-dependent.

The Habitable Worlds Observatory (HWO), NASA's next flagship science mission, follows in the tradition of the Nancy Grace Roman Space Telescope and other preceding great observatories. HWO will directly image and characterize Earth-like exoplanet and their atmospheres, with the capability to detect biosignatures and potentially answer the question of whether we are we alone. HWO will also serve as a powerful general astrophysics observatory, enabling breakthroughs in galaxy evolution, stellar astrophysics, and dark matter studies. Currently in pre-formulation, the project has established Exploratory Analytic Cases (EACs), a series of architectural concept designs used to assess the mission's demanding science objectives while exploring challenging engineering parameters. This paper describes the first three EACs, starting with observing strategies and error budget formulation and then progressing to design formulations, trade studies and lessons learned; this paper also discusses the integrated modeling pipeline, a key multidisciplinary system-level analysis capability, and analysis findings as applied to the first EAC. These activities set the stage for the follow on EACs 4 and 5, which will further explore the trade space and prepare for the baseline design that will support the Mission Concept Review (MCR).

The environments of deflectors in strong lensing systems affect our ability to test cosmological models and constrain evolutionary properties of galaxies. Here we measure the deflector scale (Einstein mass) and deflector environment (halo mass) of 89 spectroscopically confirmed strong lenses in the ASTRO3D Galaxy Evolution With Lenses (AGEL) survey. We classify deflector scale by measuring $\theta_{\rm{E}}$ to determine the mass enclosed by the Einstein radius, $M(<\theta_{\rm{E}})$. We quantify deflector environment by using photometric redshifts to determine the galaxy surface density to the fifth-nearest neighbor $\Sigma_5(z)$. We find that 47.2% of our deflectors are embedded in cluster environments, whereas only 9.0% have cluster-scale Einstein radii (masses). We measure a weak correlation ($r = 0.38$) between Einstein mass and $\Sigma_5(z)$, suggesting that the assumption of single galaxy-scale deflectors in lens modeling is overly-simplified. We hypothesize that the weak correlation results from galaxy-scale bias in the original AGEL selection and the observational challenge of detecting faint arcs with large Einstein radii. Comparing number densities, $N_{\rm{gal}}$, between AGEL and control fields, we find that AGEL deflectors are in systematically denser environments. Our study provides a method to identify strong lenses as a function of deflector environment and approximate the impact of large-scale environment in lens modeling. We provide the measured lensing parameters for our 89 AGEL systems as well as $z_{\rm{phot}}$ and $r$-mag (AB) maps of the line-of-sight.

In this work, we attempted to apply the dipole fitting (DF) method to galaxy clusters to search for cosmic anisotropic signals, and to construct a statistical isotropic analysis scheme for them. Compared to Type Ia supernova (SNe Ia), the galaxy clusters offer a significant advantage in terms of spatial distribution. This advantage makes the anisotropic signals obtained from them more reliable. From 313 galaxy clusters (Chandra + XMM-Newton), we find two preferred directions (l, b) = (${257.82^{\circ}}_{-52.88}^{+58.01}$, $-31.30{^{\circ}}_{-39.46}^{+35.92}$) and ($80.89{^{\circ}}_{-52.46}^{+60.97}$, $31.75{^{\circ}}_{-40.16}^{+35.19}$). The former to a direction where the universe is expanding at a faster rate than the surrounding area, while the latter to a slower rate of expansion. The corresponding magnitude of anisotropy is $|A|$ = 5.2 $\sim$ 5.4 $\times$ 10$^{-4}$. The results of statistical isotropy analyses give $\sim$1.0$\sigma$ confidence level. From the reanalyses based on the subsamples including Chandra, XMM-Newton, low reshift (LR, $z < 0.10$), high redshift (HR, $z > 0.10$) datasets, we find that the observation equipment and sample redshift can affect the preferred direction, anisotropic magnitude, and statistical significance of anisotropy. The XMM-Newton dataset gives a statistical significance of 2.26$\sigma$ (Mock) and 2.86$\sigma$ (Iso) which are much higher than that from Chandra and the total datasets. The magnitude of anisotropy $|A|$ from HR dataset is larger than that from LR dataset. Overall, our results indicate the presence of anisotropic signals in galaxy clusters, which must be taken seriously. Further test is still needed to better understand these signals.

Black hole jets represent one of the most extreme manifestations of astrophysical processes, linking accretion physics, relativistic magnetohydrodynamics, and large-scale feedback in galaxies and clusters. Despite decades of observational and theoretical work, the mechanisms governing jet launching, collimation, and energy dissipation remain open questions. In this article, we discuss how upcoming facilities such as the Event Horizon Telescope (EHT), the Cherenkov Telescope Array (CTA), the Vera C. Rubin Observatory (LSST), and the Whole Earth Blazar Telescope (WEBT) will provide unprecedented constraints on jet dynamics, variability, and multi-wavelength signatures. Furthermore, we highlight theoretical challenges, including the role of magnetically arrested disks (MADs), plasma microphysics, and general relativistic magnetohydrodynamic (GRMHD) simulations in shaping our understanding of jet formation. By combining high-resolution imaging, time-domain surveys, and advanced simulations, the next decade promises transformative progress in unveiling the physics of black hole jets.

We derive updated cosmological bounds on light axion-like particles (ALPs) coupled to leptons or photons, using a full phase-space treatment of their production from the primordial thermal plasma. The ALP phase-space distribution, obtained by solving the momentum-dependent Boltzmann equation for the relevant production processes, is consistently propagated into the computation of cosmological observables, allowing us to assess the impact of non-thermal spectral distortions on the effective number of relativistic species, $\Delta N_{\rm eff}$. Using state-of-the-art measurements of the cosmic microwave background from Planck, the Atacama Cosmology Telescope, and the South Pole Telescope, complemented with Big Bang Nucleosynthesis determinations of primordial deuterium and helium abundances, we obtain the following 95\% credible limits on the ALP decay constant: $f_a > 1.63 \times 10^6 \, {\rm GeV}$, $9.41 \times 10^6 \, {\rm GeV}$ and $8.06 \times 10^4 \, {\rm GeV}$ for ALPs coupled to electrons, muons and taus, respectively. For the ALP-photon coupling we find $g_{a\gamma} < 1.98 \times 10^{-8} \, {\rm GeV}^{-1}$. Including baryon acoustic oscillation data from the Dark Energy Spectroscopic Instrument mildly relaxes the constraints, in line with previous analyses of extra relativistic degrees of freedom. Finally, we present forecasts for the LiteBIRD$+$Simons Observatory and LiteBIRD$+$CMB-HD configurations, discussing the importance of an exact phase-space treatment for robust cosmological bounds on ALP interactions.

Recent hydrodynamical simulations have identified a disappearing thermal X-ray signature in massive black hole binaries (MBHBs) embedded in circumbinary disks, arising from the tidal truncation and depletion of minidiscs shortly before merger. This feature has been proposed as a promising electromagnetic counterpart to MBHB mergers detectable by LISA. In this work, we examine whether non-thermal X-ray emission powered by magnetic reconnection could obscure or modify this thermal X-ray drop. We construct semi-analytic models for both the thermal X-ray emission from minidiscs and the non-thermal synchrotron emission produced by reconnection in magnetically dominated black hole magnetospheres. Evaluating these models across the MBHB mass range relevant for LISA, we find that for physically motivated magnetic field strengths and accretion rates, the non-thermal X-ray luminosity remains several orders of magnitude below the thermal component throughout the inspiral. Even under optimistic assumptions that enhance the non-thermal emission, it remains significantly subdominant. We further incorporate the magnetospheric balding framework to model the decay of non-thermal emission near merger, finding that reconnection-powered X-ray emission fades on short, mass-scaled timescales once the external magnetic flux supply is disrupted. Taken together, our results indicate that non-thermal emission is unlikely to mask the disappearing thermal X-ray signature, reinforcing its robustness as an electromagnetic counterpart to MBHB mergers and its potential utility for multi-messenger studies with LISA.

LHAASO J2108+5157 is one of the few ultra-high energy gamma-ray sources in the LHAASO catalogue without secure counterpart at longer wavelengths. Several Galactic scenarios have been proposed, including an evolved supernova remnant and a pulsar wind nebula. Yet, no shocked gas, shell-like structure, or compact pulsar candidate has been identified. Follow-up observations with VERITAS and the LST-1 prototype have not firmly clarified its nature. A recent microquasar candidate from GMRT radio data remains uncertain. Here we present the first dedicated near-infrared study of the field, combining deep JHKs imaging with narrow band observations targeting the H2 v=1-0 S(1) line. Our observations were initially planned to encompass the full source region, but now only partially cover the latest updated position and size of LHAASO J2108+5157. We find no evidence of shocked emission, extended nebular structures, or an accreting compact object signature in the covered field. The GMRT radio source, despite its jet-like morphology, exhibits near-infrared properties incompatible with both a Galactic microquasar and a nearby radio galaxy, discouraging an association with the gamma-ray emission. Our analysis reveals no convincing counterpart consistent within the positional uncertainty, leaving LHAASO J2108+5157 as an enigmatic ultra-high energy emitter that requires deeper observations.

Minimal atomic dark matter with its distinctive cooling mechanisms offers an instructive framework for understanding the potential impact of dark matter on small-scale structure formation and early cosmology. The model consists of two fermions with opposite charges under a hidden Abelian gauge symmetry $U(1)_{D}$ and masses $m_{p_{D}}$ and $m_{e_{D}}$, respectively. Analogous to hydrogen in the Standard Model, these fermions interact via their own electromagnetic-like force, with a dark fine structure constant denoted by $\alpha_{D}$, and can bind into neutral atomic (and molecular) dark states. Previous work has largely focused on the benchmark scenario where the dark sector mirrors ordinary matter, with $m_{e_{D}}$ near the electron mass, $m_{p_{D}}$ near the proton mass, and $\alpha_{D}\sim 1/137$. We extend this analysis by investigating dark recombination and cooling physics across the full parameter space of masses and couplings. Combining Cosmic Microwave Background (CMB) measurements from Planck and ACT with BAO and Pantheon+ data, we place new constraints on the atomic dark matter parameter space, identifying regions where acoustic damping and recombination dynamics leave observable imprints on the CMB.

Gravitational production of decoupled scalars during inflationary and post-inflationary phases is efficient and can lead to over-production. We study this production with various reheating scenarios such as a generic power-law inflaton potential $V_{\rm inf}\propto \phi^k$ as well as a multi-stage reheating scenario. We derive constraints on the scalar self-interaction coupling $\lambda_s$, the mass $m_s$, and coefficients of quantum gravity-induced operators. We find that the constraints depend sensitively on the reheating dynamics. Our analysis demonstrates that universal gravity effects do not necessarily spoil the predictivity of non-thermal dark matter scenarios with $k < 4$ and low reheating temperatures, as an extended reheating phase dilutes gravitationally-produced relics. For $k > 4$, on the other hand, the relic abundance is enhanced during the reheating phase, leading to stringent constraints on the scalar. In multi-stage reheating, we show that the enhancement/dilution effect of subsequent reheating phases factorises.

We show that scalar field fluctuations alone can drive cosmic acceleration, provided the universe is spatially closed and the Compton wavelength of the field exceeds the radius of curvature. This mechanism may open new perspectives on inflation and dark energy, which could arise from a gas of sufficiently light bosons in a closed universe.

Motivated by recent accurate measurements of disk/jet coprecessions around some galactic supermassive black holes, the accelerations experienced by an uncharged, spinless object in the Kerr metric, written in harmonic coordinates, are analytically calculated up to the formal second post-Newtonian order. To such a level, some new accelerations make their appearance. They are proportional to even and odd powers of the hole's angular momentum. Their counterparts are not known where the primary is a material body. After expressing them in a coordinate-independent, vector form valid for any orientations of the hole's spin axis in space, their orbital effects are perturbatively worked out in terms of the particle's Keplerian orbital elements. The resulting expressions, averaged over one orbital revolution, are valid for generic shapes and inclinations of the orbit. The orbital plane's precession proportional to the first power of the hole's angular momentum and to the reciprocal of the fourth power of the speed of light amounts to about twenty per cent of the corresponding Lense-Thirring effect. The latter is believed to be the cause of the accurately measured disk/jet precessional phenomenology, currently measured to a few per cent accuracy. Although at a lesser extent, also the precession proportional to the second power of the hole's spin and to the reciprocal of the fourth power of the speed of light is measurable. Allowed domains in the parameter space of the jet precession around M87$^\ast$ are displayed.

Voyager spacecraft (V1 and V2) provide unique in situ measurements of perturbations propagating beyond the heliopause through the very local interstellar medium (VLISM), including the shocks and pressure fronts whose origin is debated. In particular, a jump in magnetic field strength, observed by V1 in 2020.4 at 149.3 au from the Sun, was followed by a distinct "hump" and persistently strong magnetic field, both requiring theoretical explanation. This paper offers an interpretation of those observations using a self-consistent, MHD model of the solar wind - LISM interaction driven by the OMNI and interplanetary scintillation data combined with a turbulence analysis of Voyager data. Our simulations convincingly demonstrate that global, solar-cycle-driven compressions, on hitting the heliopause, can reproduce those puzzling V1 observations. They appear to be associated with solar cycle 24, whereas similar interstellar magnetic field structures can occur once per cycle. The turbulence analysis reveals time-dependent magnetic compressibility that persists up to 165 au at scales below 10 days. Turbulence intermittency at scales below 1 hour is mostly confined to specific intervals, possibly associated with a broad foreshock region. The apparent disappearance of intermittency since 2022 reflects the turbulence weakening rather than a fundamental change in VLISM properties. We predict that V1 will record relatively strong magnetic field strengths until $\sim$2030, followed by weaker, infrequent perturbations. At V2, we expect multiple solar-driven compressions before 2026, followed by a major event induced by solar cycle 25 around 2030. New Horizons is expected to cross the termination shock at 80$\pm$ au in 2031.

The origin of cosmic magnetic fields is widely attributed to the amplification of weak seed fields by turbulent dynamos. However, a critical understanding gap remains between the microscopic generation of these seeds and the macroscopic onset of the dynamo. Current kinetic models, often constrained to electron scales, predict premature saturation via magnetic trapping, leaving the generated fields potentially too weak and small-scale to effectively prime magnetohydrodynamic (MHD) processes. Here, using high-resolution kinetic simulations with a realistic mass ratio, we reveal the physics of this unexplored ion-kinetic regime. Under generalized continuous shear driving, used to simulate ubiquitous macroscopic flows, we demonstrate that the saturation of electron instabilities is not the endpoint but a precursor to a distinct, ion-dominated evolution. Massive ions, sustaining the velocity shear, trigger a subsequent filamentation instability that accesses the vast ion kinetic energy reservoir. This mechanism amplifies the magnetic energy by orders of magnitude beyond the electron-saturation limit, expanding the field coherence to ion scales. Our results establish ion kinetics as the essential ''missing link'' that bridges the divide between microscopic plasma instabilities and macroscopic cosmic dynamos.

Unidentified Falling Objects (UFOs) refer to sporadic beam losses observed during LHC operation. The prevailing explanation is that micrometer-sized dust particles released from the beam screen produce beam losses through interactions with the protons. However, the release mechanism of these particles remains unknown. We propose that roughly $(1-10)$% of UFOs may be caused by axion quark nuggets (AQNs), macroscopic dark matter (DM) candidates with masses of order $(5-1000)\,$g. The AQN model naturally relates the dark- and visible-matter abundances ($\Omega_\mathrm{DM}\sim\Omega_\mathrm{visible}$) and provides a mechanism for generating the baryon-antibaryon asymmetry, with DM composed of both matter and antimatter AQNs. When passing underground within approximately 100km of the LHC, an antimatter AQN generates acoustic waves strong enough to trigger multiple UFO events within $2\,$s. If three correlated UFOs (placed at different locations along the LHC ring) are detected, the signal-to-noise ratio can exceed 5 across the entire allowed AQN mass range for a measurement time of about 360 hours. Practically, the LHC can serve as a large broadband acoustic detector for AQNs.

We propose a system to optimize parametric designs subject to radiation pressure, \ie the effect of light on the motion of objects. This is most relevant in the design of spacecraft, where radiation pressure presents the dominant non-conservative forcing mechanism, which is the case beyond approximately 800 km altitude. Despite its importance, the high computational cost of high-fidelity radiation pressure modeling has limited its use in large-scale spacecraft design, optimization, and space situational awareness applications. We enable this by offering three innovations in the simulation, in representation and in optimization: First, a practical computer graphics-inspired Monte-Carlo (MC) simulation of radiation pressure. The simulation is highly parallel, uses importance sampling and next-event estimation to reduce variance and allows simulating an entire family of designs instead of a single spacecraft as in previous work. Second, we introduce neural networks as a representation of forces from design parameters. This neural proxy model, learned from simulations, is inherently differentiable and can query forces orders of magnitude faster than a full MC simulation. Third, and finally, we demonstrate optimizing inverse radiation pressure designs, such as finding geometry, material or operation parameters that minimizes travel time, maximizes proximity given a desired end-point, minimize thruster fuel, trains mission control policies or allocated compute budget in extraterrestrial compute.

Binding energies (BEs) of adsorbates on interstellar dust grains critically control adsorption, desorption, diffusion, and surface reactivity, and therefore strongly influence astrochemical models of star- and planet-forming regions. While recent computational studies increasingly report full distributions of BEs rather than single representative values, these distributions are typically derived for either bare grain surfaces or thick water-ice mantles. In this work, we bridge these regimes by systematically investigating the BE distributions of water on partially and fully ice-covered dust grain surfaces. We employ machine-learning interatomic potentials (MLIPs) based on graph neural networks to model water adsorption on graphene and on the Mg-terminated (010) surface of forsterite, representing carbonaceous and silicate grains, respectively. The models enable extensive sampling of adsorption sites on water clusters, monolayers, and bilayers generated under both crystalline (thermally processed) and amorphous (low-temperature) growth conditions. At submonolayer coverage, the chemical nature of the underlying grain strongly affects both ice morphology and binding energies, with Mg-O interactions on silicate surfaces producing particularly deep binding sites. From monolayer coverage onward, adsorption on both substrates is dominated by hydrogen bonding within the ice, reducing the influence of the grain material. Across all coverages, amorphous ice structures systematically shift the BE distributions toward stronger binding compared to crystalline ice, introducing highly stable defect and pocket sites. These results demonstrate that BE distributions in the submonolayer to few-layer ice regime are broad and highly surface dependent, and they provide physically motivated input for next-generation astrochemical models incorporating surface heterogeneity.

Anisotropic, non-thermal, and multi-temperature distributed particle momenta are commonly observed in collisionless space plasmas, such as the solar wind. Using Liouville's theorem, we argue that anisotropic compression or expansion of the plasma, followed by a relaxation of the resulting anisotropic stress must lead to non-equilibrium states that are either anisotropic, non-thermal distribution functions, different electron and ion temperatures, or a combination of these effects. We present arguments showing that a plasma in thermal equilibrium undergoing anisotropic compression or expansion cannot return to thermal equilibrium in the absence of particle collisions. Since most astrophysical plasmas are practically collisionless and experience significant anisotropic compression or expansion, we expect anisotropic, non-thermal, and multi-temperature particle distributions to be ubiquitous, in agreement with solar wind measurements.

The early universe grants access to energy scales far beyond those achievable in terrestrial experiments and allows unstable Standard Model particles to play an active dynamical role. In this contribution, we focus on recent studies aimed at quantifying the potential of the early universe to probe the properties and interactions of axions. The discussion is organized around four classes of axion scenarios, ordered from long to short lifetimes: (i) stable or long-lived axions contributing to dark radiation; (ii) stable or long-lived axions produced out-of-equilibrium and constituting dark matter; (iii) metastable axions whose decays inject energy into the primordial plasma and leave observable signatures in the global 21 cm signal; and (iv) very short-lived axions that act only as portals to additional degrees of freedom. Together, these scenarios highlight the interplay between axion phenomenology and early universe cosmology and demonstrate the potential of cosmological data to probe axions over a broad range of masses and lifetimes.