Papers for Monday, Apr 13 2026

We consider a scenario in which the Milky Way (MW) and M31 have had a previous pericentric passage, and investigate its compatibility with self-interacting dark matter (SIDM). Using initial conditions sampled from Local Group (LG) analogues in the IllustrisTNG simulation, we perform controlled re-simulations of the MW-M31 orbit, evolving the system under both standard cold dark matter (CDM) and various SIDM cross-sections. We find that the deep baryonic potential of the MW preconditions the halo's thermal structure, establishing an initial negative temperature gradient. This drives SIDM halos to bypass the standard core-formation phase and enter immediate core-collapse, resulting in monotonically increasing central densities. In full orbital simulations, the compact stellar component (disk/bulge) of the MW analog remains robust against tidal disruption for pericenter distances as close as $r_{\rm peri}\lesssim20$ kpc during an encounter at cosmic time $\sim8$ Gyr. The diffuse stellar halo is comparatively more susceptible, facing disruption for $r_{\rm peri}\lesssim100$ kpc. Our results demonstrate a dichotomy in structural evolution: the compact disk/bulge is sensitive to intrinsic SIDM thermodynamics but dynamically robust against the pericenter encounter, whereas the diffuse stellar halo is largely independent of the specific SIDM model but more vulnerable to orbital tidal disruptions.

We motivate the use of differentiable probabilistic programming techniques in order to account for the large model-space inherent to astrophysical $\gamma$-ray analyses. Targeting the longstanding Galactic Center $\gamma$-ray Excess (GCE) puzzle, we construct differentiable forward model and likelihood that make liberal use of GPU acceleration and vectorization in order to simultaneously account for a continuum of possible spatial morphologies consistent with the GCE emission in a fully probabilistic manner. Our setup allows for efficient inference over the large model space using variational methods. Beyond application to $\gamma$-ray data, a goal of this work is to showcase how differentiable probabilistic programming can be used as a tool to enable flexible analyses of astrophysical datasets.

We propose a novel cosmological scenario in which baryonic neutron stars could plausibly form in the early universe. If baryogenesis initially produces an excessively-large baryon asymmetry, $Y_B \gg 10^{-10},$ the baryonic mass inside the horizon can exceed the minimum neutron star mass before big bang nucleosynthesis (BBN). While this large asymmetry is present, non-relativistic baryons can dominate the universe and enhanced density perturbations on small scales can gravitationally collapse Hubble patches shortly after horizon re-entry. For some initial perturbations, just below the threshold for black hole formation, this collapse will be arrested only by nuclear pressure, possibly resulting in neutron star formation. Afterwards, there must be a large entropy injection to restore the observed baryon asymmetry, $Y_B \sim 10^{-10}$, and preserve the successful predictions of standard BBN. Unlike neutron stars that form from stellar collapse, primordial neutron stars can, in principle, be as light as $\sim 0.1 M_\odot$, limited only by the nuclear equation of state.

The Milky Way is surrounded by a hot diffuse circumgalactic medium (CGM) with temperatures of millions of degrees. Recent X-ray observations with the eROSITA satellite discovered a significant temperature asymmetry of this hot CGM, with the southern hemisphere being on average hotter than the northern one by a relative difference of ${\Delta} T/T \approx 12\%$, where $T$ is averaged over the entire CGM. In this Letter, we investigate whether the passage of the Magellanic Clouds can be responsible for this asymmetry by means of a hydrodynamical/N-body simulation. In the simulation, the Magellanic Clouds induce a relative motion of the Milky Way's disc of up to 40 km/s. This motion leads to compression of the CGM gas in the southern hemisphere, resulting in an overall temperature increase in that region. We estimate a south-north temperature difference of ${\Delta} T/T \approx 13-20\%$, consistent with the observations. We find that this temperature asymmetry is a recent phenomenon that began ~100 Myr ago.

Coronagraph designs which use photonic integrated circuits have the highest theoretical throughput for off-axis signals, and therefore the highest potential exoplanet yield for future high-contrast direct imaging campaigns. Using the rejected starlight, the photonic integrated circuit may also provide simultaneous wavefront sensing, allowing for the correction of non-common-path aberrations. This work considers how a photonic circuit should be configured to maximize its sensitivity to phase aberrations. Two cases are considered: in the first, the photonic circuit is coupled directly to an electric field in a piecewise manner, while in the second, the circuit is coupled to the field via an optical mode sorter. In either case, this work constructs a unitary matrix which can be applied by a photonic circuit to produce maximum sensitivity.

Whether the "Hubble sequence" of galaxy morphologies exists up to z~4 is still disputed, and one of the challenges is characterizing galaxy structure consistently across a wide range of redshifts. To enable a fair comparison across cosmic time, we constructed "absolute" images of galaxies spanning 0.15<z<4.5 and 8<log $M_{\star}$<11 from HST CANDELS and JWST CEERS surveys, by matching the effective resolution and surface brightness limit of galaxies, accounting for cosmological dimming and evolution in size and mass-to-light ratio. We measured the structural parameters of 2825 galaxies and used the UMAP technique to study the evolution of the morphological phase space. We find a continuous sequence spanning late-type to early-type galaxies, with no redshift gradient - indicating that a Hubble-like sequence is established by z~4. We show that our approach recovers a cleaner separation between early- and late-type galaxies than visual classifications. By tracing progenitors using empirical mass assembly histories, we find that progenitors of low-mass galaxies are predominantly star-forming disks at all epochs. Progenitors of massive galaxies follow two distinct paths: a stable star-forming disk population with little structural evolution, and an early-type population that builds up rapidly from irregular progenitors and quenches within a few Gyr, consistent with a compaction-driven quenching scenario.

The Galactic core-collapse supernova (SN) rate is estimated at $\approx 1-3$ per century; however, no optically visible SN has been discovered in the past 400 years. Although records of the last optically detected SN (Cassiopeia A) are debated, it is revealed today via its bright, variable mid-infrared (MIR) dust echoes -- offering the possibility of identifying dust-obscured, missed events via their dust echoes. We present the first all-sky, untargeted search for thermal dust echoes of luminous Galactic transients using difference imaging on 12 years of time-resolved NEOWISE co-adds (spanning $2009-2022$) followed by statistical detection of variable extended sources. We use echo features around Cas A, together with archival catalogs to train a convolutional neural network to classify transient candidates as dust echoes, point sources, artifacts, and high proper motion stars. Our model achieves $\approx 94$% accuracy in distinguishing echoes from other variable sources. Applying the classifier to $\approx 11$ million transient candidates, we search for spatial over-densities of echoes across the Galactic plane. We find that Cas A is the only region exhibiting echoes at the WISE sensitivity threshold of $W2$ surface brightness of $\approx 20$ Vega mag arcsec$^{-2}$ -- reflecting its unique combination of young age and luminous shock breakout. We present the largest catalog of time-resolved echo positions of Cas A (20477 within 10$^\circ$) that are being used for studies of the surrounding interstellar medium with the James Webb Space Telescope. Our results lay the groundwork for the imminent Roman space telescope surveys -- which will achieve $\approx 100\times$ higher sensitivity and $\approx 30\times$ better spatial resolution at wavelengths of $\lesssim 2.5\,\mu$m.

We report the discovery of GLIMPSED-329380, a z=6.64 galaxy behind Abell S1063, which shows signs of an extreme ionised outflow driven by an active galactic nucleus (AGN). The deep JWST/NIRSpec medium grating observations show spatially resolved structures of a host galaxy containing the very fast outflow and an AGN, which we analyse separately. The outflow, mainly traced by broad [O III]{\lambda}5008 and H{\alpha} emissions in the host, reaches a full-width half-maximum velocity of ~5500km/s, velocities only observed in AGN-dominated systems. From the Balmer decrement, we observe that while the narrow emission lines show no dust attenuation, the outflowing gas is dusty. We use emission lines diagnostics to infer gas abundances within the host galaxy. The oxygen abundance is 12+log(O/H) ~ 7.95 (~18% solar) and the host is slightly nitrogen-enriched with log(N/O) ~ -0.75. Despite its extreme velocity, the mass loading factor (<0.1) and the kinematic energy of the outflow (~10^43 erg/s) suggest limited impact on star formation. The AGN component shows many similarities with little red dots (LRDs): a characteristic "V-shape", exponential profile in hydrogen lines, numerous detection of forbidden [Fe II] lines, a Balmer break, and a broad absorption feature at ~4550 Å. This detection of a fast outflow in an LRD, rare in surveys dominated by low-resolution (e.g. PRISM) spectra, provides direct evidence of AGN activity in these systems.

We present the 1.6$-$28 $\mu$m spectra of the young protostar EC 53, obtained with JWST NIRSpec IFU and MIRI MRS during the quiescent and burst phases of its periodic brightness variations. To isolate ice absorption features, we modeled and removed the mid-infrared silicate dust absorption using a dedicated continuum-fitting procedure. In the optical depth spectrum, we first fit the broad H$_2$O ice features and then decomposed the major ice components, including NH$_3$, CO$_2$, CH$_3$OH, CO, and CH$_4$, by matching laboratory profiles for both pure and H$_2$O-mixed ices. The 4.62 $\mu$m and 6.85 $\mu$m bands are attributed to OCN$^-$ and NH$_4^+$ ions, respectively. Minor or tentative contributions from complex species (HCOOH, H$_2$CO, CH$_3$COOH, CH$_3$CHO, CH$_3$CH$_2$OH, and NH$_2$CHO) are also considered to our global ice analysis. The silicate-corrected spectra reveal no measurable change in any ice absorption band between the two phases, indicating that moderate and short-period accretion bursts in EC~53 do not significantly alter the physical or chemical state of the ices within its envelope. The derived abundances of these major species relative to H$_2$O significantly exceed the values typically observed toward other embedded protostars. Finally, we place the ice inventory of EC~53 in the context of other protostellar systems observed with JWST, highlighting that its chemically rich, thermally quiescent ice reservoir provides a benchmark for studying ice evolution under episodic accretion.

Most scientific research begins in the context of the then-contemporary state of knowledge of the field and moves toward a deeper understanding of the subject. This colloquium presents the Arecibo telescope$^{'}$s contribution to radio astronomy from the point of view of its contemporary impact and relevance and how that evolved over the 57 years of its long life. Sometimes, serendipitous discoveries at Arecibo and elsewhere brought revolutionary changes to the field. Further, significant upgrades to the reflector, optics, receiving, and data-taking systems helped clear a pathway for a leap ahead. Charting these movements through time and without any claim to completeness, this Colloquium presents a progression through Arecibo telescope$^{'}$s role in the history of radio astronomy.

A new digital optical module (DOM) has been developed for the proposed expansion to the IceCube detector at the South Pole, IceCube-Gen2. The "Gen2-DOM" has 4 times the integrated photon sensitivity of the current IceCube DOMs and has built off the design features of the IceCube Upgrade modules. The Gen2-DOM has up to 18 4" photomultiplier tubes (PMTs) in a borosilicate glass pressure vessel, arranged in a uniform 4$\pi$ angular distribution. The mechanical design has been optimized to fit into a reduced borehole diameter which, in turn, will reduce drilling costs during installation. Each PMT has a dedicated readout board, designed to increase sensitivity to high-energy events aligned with the science goals of IceCube-Gen2. Internal storage enables multi-level triggering schemes with reduced overall flow of data on the long cables. Twelve prototypes of the Gen2-DOM will be deployed in the IceCube Upgrade in the 2025-2026 austral summer. This article will focus on the current status of design development and initial performance testing results.

We present the results of experiments probing the retention of CO2 in crystalline water ice, frozen sodium chloride (NaCl) brines, and flash-frozen carbonated water using diffuse reflectance infrared spectroscopy. Characteristic absorptions alluding to the formation of clathrate hydrates in crystalline ices and frozen brines are observed. NaCl in frozen brines does not affect qualitatively affect the formation of clathrate hydrates. Generation and stability of clathrates in crystalline ice transiently subjected to pressure-temperature (P-T) conditions in the stability region is observed, despite conditions being unviable at the onset of freezing. Retention of CO2 in flash-frozen carbonated water is observed to be dependent on the temperature of the substrate during freezing. The state of CO2 retained in the resulting ices differs from clathrate hydrates, as inferred from the respective infrared spectra. Both mechanisms of CO2 retention are stable up to 140 K and under evacuated conditions. In the context of Europa, the P-T states traversed by the samples plausibly represent the typical conditions around endogenous CO2 if it is indeed transported from the subsurface ocean to the surface while being retained in ice/frozen brines and/or liquid emerging on the surface. However, the absorptions of CO2 in the laboratory infrared spectra do not match those detected on the leading side of Europa by the NIRSpec instrument on board JWST. Therefore, it is unlikely that the endogenous CO2 observed at the surface of Europa is sourced directly from the ocean, unless additional processes affect the observed bands of CO2 on Europa.

The method of weighted least squares is widely used in parameter estimation problems such as asteroid orbit determination. A key difficulty is the treatment of observational uncertainties, especially when combining heterogeneous datasets with differing precision. We propose a simple reweighting scheme that adjusts the contribution of each measurement group to ensure a statistically consistent least-squares solution. It consists of three steps: (i) estimating error standard deviations for each observational subset, (ii) rescaling their weights by the corresponding variances, and (iii) a weighted least-squares fit with the adjusted weights. We apply this to heliocentric orbit fitting of asteroids using ground-based astrometry and high-precision Gaia measurements. We validated the method by fitting each orbit to a restricted set and comparing with the complete set of measurements. For 7 objects, the reweighted solutions provide significantly improved agreement with older data. The most dramatic case is asteroid (21) Lutetia, where increasing the effective uncertainty of Gaia observations by a factor of 17 yields a substantially better fit, indicating the importance of accounting for systematic biases in high-precision datasets. We further apply the scheme to near-Earth asteroid 2024 YR4, grouping observations by visual magnitude. The reweighted orbit produces smaller uncertainty regions and a more stable solution, reducing predicted impact probabilities by roughly an order of magnitude. All computed probabilities remain below 0.5%, under the 1% International Asteroid Warning Network (IAWN) alert threshold. This reweighting procedure provides a practical way to combine heterogeneous measurements, improving the reliability of orbit determination and impact-risk assessment. The method is general and can be readily applied to other parameter estimation problems involving mixed-precision data.

The Hayabusa2 extended mission, nicknamed Hayabusa2# (# is pronounced SHARP, which stands for the Small Hazardous Asteroid Reconnaissance Probe), is JAXA's small body explorer to conduct science and engineering investigations in space. After the successful return to the Earth with the samples from the carbonaceous asteroid (162173) Ryugu on December 6, 2020, Hayabusa2 diverted away from Earth to start its decade-long extended mission. The major scope includes engineering demonstration of long-term maintenance strategies for spacecraft and operation systems and scientific investigations during various mission phases. Major scientific investigations include spacecraft-based telescopic observations of exoplanets and zodiacal dust observations during the cruise phase, flyby observations of the near-Earth asteroid (98943) Torifune in July 2026, and rendezvous observations of near-Earth asteroid 1998 KY26 in 2031. This study overviews Hayabusa2#'s flyby and the physical properties of Torifune. Although the flyby operation planning is still ongoing, the mission will attempt to fly by the target at a distance (from the asteroid's center) of ~1-10 km. The flyby speed is planned to be 5.25 km/s, while the encounter location is 0.81 au from the sun. The mission plans to fix the spacecraft's orientation during the flyby, only allowing for a very limited pointing change to attain higher resolution imaging. The mission will attempt to obtain science and engineering returns during the flyby. The planned investigations will offer stronger insights into material transport mechanisms in the inner solar system and a demonstration of planetary defense technologies.

We present a validation of our recently proposed non-conventional method, Constant Acceleration Accounted Perspective (CAAP), for estimating the instantaneous expansion speed of coronal mass ejection (CMEs), even when only single-point in situ observations are available. This validation is enabled by the radial alignment of SolO and Wind spacecraft (0.13 AU radial and 2.3 deg angular separation), providing simultaneous observations of the center (at Wind) and trailing edge (at SolO) of a CME associated magnetic cloud (MC) during 3-5 November 2021, allowing a direct measurement of its instantaneous expansion speed. These measurements are compared with CAAP-derived instantaneous expansion speed estimates at both spacecraft. The favorable spacecraft configuration also enables tracking the temporal evolution of CME substructures, including the shock, sheath, and MC. A discrepancy is noted between the low-inclination MC axis estimated from minimum variance analysis (MVA) and the highly inclined ENW-type MC axis suggested by visual inspection of in situ measurements. We also observe an apparent increase in the magnetic flux within the MC from SolO to Wind, indicating a noticeable deviation from magnetic flux conservation. During the CME's propagation from SolO to Wind, the shock becomes unexpectedly stronger at Wind, while the sheath thickness remains nearly the same, likely due to MC acceleration from back compression by a high-speed stream and ambient solar wind variability. Our results demonstrate the applicability of the CAAP method and the importance of accounting for temporal evolution in CME substructures for space weather studies.

cTreeBalls (cBalls for short) is a Python/C package useful to measure (2,3)-point clustering statistics. cBalls can efficiently calculate 3-point correlations of more than 200 million HEALPix pixels ( a full sky simulation with Nside = 4096) in less than 10 minutes on a single high-performance computing node, enabling a feasible analysis for the upcoming LSST data. It builds upon octree (Barnes & Hut, 1986) and kd-tree algorithms (Bentley, 1975), and supplies a user-friendly interface with flexible input/output (I/O) of catalogue data and measurement results, with the built program configurable through external parameter files and tracked through enhanced logging and warning/exception handling. For completeness and complementarity, methods for measuring two-point clustering statistics for periodic boxes are also included in the package. cTreeBalls was developed for its use in the Dark Energy Science Collaboration (DESC) of the Rubin Observatory Legacy Survey of Space and Time (LSST).

We present a theoretical model for the power spectrum and bispectrum of galaxy clustering that exploits the complementarity between small-scale power spectrum information and large-scale bispectrum measurements. We extend the FOLPS code by combining its one-loop EFT galaxy power spectrum with a tree-level galaxy bispectrum projected onto the tripolar spherical harmonics (Sugiyama) basis. To access additional small-scale information, we also consider a line-of-sight damping factor in both statistics, mirroring approaches commonly used in studies of redshift-space distortions. We test the model using DESI DR2 galaxy mocks. Even without damping, the joint analysis of the EFT power spectrum and bispectrum significantly improves constraints and reduces parameter degeneracies relative to power spectrum analyses alone. For LRG-like samples, including the damping further extends the range beyond $k\sim 0.3 \,h \text{Mpc}^{-1}$ in the power spectrum and $k \sim 0.24 \,h \text{Mpc}^{-1}$ in the bispectrum without introducing statistically significant parameter biases. This leads to up to $\sim 30\%$ tighter constraints on $A_s$ and $\omega_{cdm}$. For low signal-to-noise tracers such as QSOs, however, the damping parameters are weakly constrained and can absorb noise fluctuations, leading to shifts in inferred parameters. Similar limitations may arise in models where cosmological information is encoded in power-spectrum shape features degenerate with the damping, such as scenarios with massive neutrinos. In contrast, for $w_0w_a$CDM we obtain $15\%$ and $21\%$ tighter constraints on $w_0$ and $w_a$, respectively, yielding a deviation from constant dark energy at slightly more than the $1\sigma$ level using full-shape information alone. The code is publicly available at this https URL

Previous studies using p-H$_2$CO $J=3$--$2$ transitions at 218 GHz suggested widespread high-temperature gas exceeding 60 K and even 100 K in the CMZ, with heating mechanisms possibly related to cosmic rays or turbulent dissipation. However, at temperatures above 100 K, p-H$_2$CO $J=3$--$2$ line emission may lead to significant overestimates of kinetic temperature. This study combines o-H$_2$CO $J=5$--$4$ data from JCMT with p-H$_2$CO $J=3$--$2$ data from APEX to analyze three molecular clouds (The Brick, Sgr A1, and Sgr A2) with high temperatures. We used the non-LTE radiative transfer code RADEX to model spectral lines and constrain physical parameters with multiple line ratios, obtaining more reliable kinetic temperatures. Our results show that the previously reported extreme temperatures ($>100$ K) based on p-H$_2$CO $J=3$--$2$ line ratios are revised downward, with the average kinetic temperatures now constrained to 84--95 K using o-H$_2$CO $J=5$--$4$ line ratios, indicating systematic overestimation in the earlier studies. Further analysis reveals that the relationship between temperature and gas line width aligns more closely with predictions from models incorporating both high cosmic ray ionization rate and turbulent heating, suggesting that these molecular clouds are likely heated by a combination of cosmic-ray and turbulent dissipation mechanisms.

A convolutional neural network (CNN) is used to construct a new catalog for solar flares based on high resolution (1-s cadence) Geostationary Operational Environmental Satellites (GOES) soft X-ray data. The CNN is trained to identify flare rise episodes. From 1 January 2018 to 22 August 2025, the algorithm detects 111,580 flare candidates, compared with 14,612 events in the corresponding GOES catalog. For each candidate, the probability of being a true positive is quantified by Bayesian inference based on the peak flux, rise time, and temporal coincidence with cataloged events where available. The flare size and waiting-time distributions are studied and compared with the GOES catalog. The CNN catalog shows a steeper power-law index for raw peak fluxes (-2.59 -\+ 0.02) than GOES (-2.25 -\+ 0.04), indicating the CNN's higher sensitivity to small events. After background correction, the indices are -1.97 -\+ 0.02 (CNN) and -2.05 -\+ 0.04 (GOES). The CNN catalog extends the power-law distribution of flare peak fluxes by one order of magnitude at the small-flux end compared with the GOES background-subtracted catalog. A Bayesian blocks analysis of the waiting-time distributions from the GOES and CNN catalogs indicates broad consistency with a piecewise Poisson process. We find that previously reported correlations between flare sizes and waiting times are significantly influenced by obscuration, that is, under-counting weaker or overlapping flares during periods of elevated flux. The new CNN catalog provides a foundation for complete and consistent studies of solar flare statistics.

In recent decades, the relevance of polarimetry in planetary sciences and astronomy has increased rapidly. Polarization is a fundamental property of light and can be modified by any scattering event. As such, polarization yields additional information that cannot be obtained by only assessing light's scalar properties. For instance, the polarization state of starlight scattered by planetary surfaces can provide useful insights on the composition, size, morphology, and porosity of regolith particles and might even indicate the presence of life. Beside being useful for characterization, polarimetry can also greatly enhance the detection of exoplanets. Here, polarization can be harnessed to enhance the contrast between the bright light of a star, which can be considered to be fully unpolarized, and the very dim but polarized light reflected by an exoplanet. In this paper, we discuss and review the current developments and advances in optical polarimetry and polarimetric instrumentation in Switzerland within the framework of the National Centre of Competence in Research PlanetS. We focus on their implications for the vast range of science cases that polarimetry can address within the research fields of planetary science and astronomy.

This work describes the context and approach for the detection of spectroscopic signatures from planets in the habitable zone of nearby stars. By understanding the limitations of current observatories, future telescopes can be understood, and their ability to characterize the atmospheres of exoplanets estimated. An example calculation is given for the signal-to-noise analysis for a planet like the current Earth of oxygen as a biosignature, and (an enhanced abundance) of hydrogen iodine as a technosignature. In the optimistic estimate, Earth is easily detected, O2 characterized in 20 hours, but signals from enhance HI are only visible after hundreds of hours, indicating the signals are too weak to realistically constrain.

Impacts play a fundamental role in shaping the physical and chemical properties of the objects in our Solar System. Given the challenges in replicating such collisions through laboratory experiments, computer simulations are an important tool to investigate their outcomes. Accurately modelling material properties such as shear strength, porosity, and the formation of cracks is crucial for understanding impacts on small bodies like asteroids and comets. Very large and massive objects are dominated by self-gravity and can be approximated as a fluid. In this regime the equation of state used to model the behaviour of the constituent materials plays a key role. However, for bodies of several hundred kilometres, which are already spheroidal due to self-gravity, shear strength must still be considered. This impact regime is most challenging to model and therefore often overlooked in publications. In this review we present different impact regimes and the relevant physics that must be included. We then discuss their application to a variety of Solar System objects and assess how recent observations and numerical simulations, focussing on the Smoothed Particle Hydrodynamics method, can be used to inform our understanding of impact processes and solar system formation.

We study the effects of using the optical properties of irregular hexahedral grains in photoionization models of circumstellar nebulae around evolved stars. Dust opacities for the irregular grains were obtained from the scattering properties available in the TAMUdust2020 database and these were implemented in the spectral synthesis code cloudy. A sample of photoionization models that use opacities from both spherical and irregular hexahedral grains across a standard MRN size distribution (0.005 to 0.25 um) was produced. We consider the optical properties of graphite, amorphous carbon and silicate dust grains and find that differences between the model nebula continua calculated using spherical and irregular dust grains increase with the grain size, especially for graphite. In particular, we find that the luminosities at the infrared peak for the hexahedral grain models can be up to 60% higher than those from the equivalent spherical grain models for the largest grains. This result suggests that traditional spherical grain assumptions may lead to an overestimate of the dust mass in photoionized nebulae.

We review the progresses made in global theoretical models of planetary system formation in the last decade using the example of the planetary system formation framework known as the Bern Model that has been continuously developed since before the beginning of the NCCR PlanetS. We highlight major developments and applications that have since been implemented, reflecting important recent advancements of planet formation theory overall, such as MHD wind-driven disk evolution, planetesimal evolution including fragmentation, dust evolution and pebble accretion, formation of planets in structured disks, interior structure models allowing for compositional gradients, as well as the analysis of the emerging planetary system architectures and the identification of different classes of architectures. We discuss how these new models impact the formation and evolution process and translate into different populations of planets and planetary systems. We also discuss the major strengths of the Bern Model, including successful predictions of the break in the planetary mass function at 30 MEarth, the prevalence of low-mass planets, the radius pile-up around 1 RJupiter, and the evaporation valley, with the recent New Generation Planetary Population Synthesis models with 100 seeds per disk providing quantitive matches to many RV-survey and Kepler diagnostics. This includes key characteristics of planetary system architectures. We also highlight the limitations of this model, some of them were addressed during the course of the NCCR PlanetS: the inclusion of the early phases of planet formation from dust to planetesimals, the hybrid pebble-planetesimals accretion of solids, simplified interior structure models, reliance on simplified parametrizations that may not encapsulate the full complexity of physical processes, and computational constraints.

High-precision high-fidelity spectrographs are the most powerful instruments for exoplanets detection and characterization. The sub-m/s radial-velocity precision, required to detect Earth-mass exoplanets, necessitates tackling all the sources of instrumental and stellar instabilities. We present the new high-precision high-fidelity spectrographs ESPRESSO, NIRPS, ANDES and RISTRETTO designed, developed, and operated with support of PlanetS.

The disk instability (DI) model for giant planet formation remains an attractive alternative in explaining the formation of giant planets at early times, giant planets at large radial distances, and giant planets orbiting M-stars. In this review, we present recent developments in the disk instability model including hydrodynamical as well as magneto-hydrodynamical (MHD) disk simulations, populations synthesis models, and simulations of clump-clump collisions. We also discuss advances in observations that can be used to constrain and test this formation scenario.

The origin of low surface brightness (LSB) galaxies remains a key open question in galaxy formation, reflecting the balance internal mechanisms and environmental influence. Using MaNGA integral-field spectroscopy, we investigate whether LSB and high surface brightness (HSB) galaxies of comparable stellar mass ($9 < \log M_\ast < 10$) occupy distinct environments or differ primarily through internal evolution. Our late-type sample comprises 113 central and 29 satellite LSB galaxies, and 374 central and 142 satellite HSB galaxies. We characterize environments on scales from 100 kpc to 10 Mpc, analyzing radial profiles of stellar mass surface density ($\Sigma_\ast$), star formation activity, and gas-phase metallicity. Central LSB and HSB galaxies inhabit similarly low-density large-scale ($>$200 kpc) environments, but LSB galaxies are more isolated on small scales ($\sim$100 kpc). Even after matching in stellar mass and environment, LSB galaxies show systematically lower $\Sigma_\ast$, $\Sigma_{SFR}$, and metallicities, often hosting diffuse, weakly star-forming bulges embedded in extended disks. These results indicate that LSB structure and star formation are not primarily governed by large-scale environment or halo mass. While secondary halo properties such as spin, concentration, or gas accretion history are often invoked, their environmental dependence appears weak. Instead, LSB-HSB differences for centrals likely reflect divergent assembly or interaction histories and internal processes -- such as angular momentum-driven disk evolution or inefficient gas conversion -- largely decoupled from large-scale environment. Nonetheless, environment still influences the observed star formation and chemical differences between central and satellite LSB galaxies.

We continued the search for single pulses (SPs) in the northern part of the all-sky High Time Resolution Universe survey, whose aim is to detect pulsars and other radio transients. This search is now about 21% complete and has yielded the first discovery of a fast radio burst (FRB) with the 100 m Effelsberg Radio Telescope. FRB20110220A was detected with an S/N-optimised dispersion measure of 501.0 pc/cm$^{3}$ and a width of 11.9 $\pm$ 3.5 ms, for a fluence of 0.6 $\pm$ 0.1 Jy ms. We obtained the first L-band detection of the rotating radio transient (RRAT) J2028+28, from which we obtained upper limits on the source's period and burst rate, as well as an improved position. We also discovered a new RRAT, J0404+53, which had previously been reported as an isolated SP candidate. Eight new SP trains and 272 faint isolated SP candidates were detected too. We used these candidates to demonstrate that their all-sky detection rates depend on Galactic latitude and longitude. This direction dependence suggests the existence of a faint Galactic SP population.

We extracted a list of 662 nearby (within $\sim16$ Mpc) HI-detection sources from the Five-hundred-meter Aperture Spherical radio Telescope (FAST) All Sky HI Survey (FASHI) and made a visual identification of them with optical counterpart. This inspection led to the discovery of 71 new dwarf galaxies. All of them are dwarf irregular galaxies with ongoing star formation. They are characterized by the following median parameters: visual magnitude of $g=17.8$ mag and color $(g-r)=0.29$ mag, HI-flux $S_\mathrm{HI}=718$ mJy km/s, HI-mass $M_\mathrm{HI}=3.7\times10^7$ $M_\odot$, as well as the HI line-width of $W_{50}=37$ km/s.

FRB 20220912A is a highly active repeating fast radio burst (FRB) source, discovered by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using its real-time FRB detection system (CHIME/FRB). Here, we present results from a radio monitoring campaign of FRB 20220912A using CHIME, including ~200 hours of data collected by CHIME/Pulsar, spanning 1.5 years following the source's discovery. We present an analysis of a sample of 828 CHIME-detected bursts from FRB 20220912A, in the 400-800 MHz radio frequency band. The source remains highly active for ~10 weeks and has a bimodal wait-time distribution with peaks at $160^{+120}_{-70}$ ms and $306^{+14}_{-13}$ s. Assuming a radio efficiency factor of $10^{-4}$ and a beaming angle of 0.1, we estimate the total emitted energy from the source over the entire observing campaign to be $2 \times 10^{43}$ ergs. We report a 2.3$\sigma$ detection of a linear increase in the DM of $1.4 \pm 0.6$ pc cm$^{-3}$ yr$^{-1}$, with no significant trend in rotation measure (with a 3$\sigma$ upper limit of 13.4 rad m$^{-2}$ yr$^{-1}$). We contrast our findings with other active repeaters, which exhibit different DM and RM evolution to indicate that FRB 20220912A may reside in a unique local environment.

The nucleosynthetic heterogeneity between different asteroids and planets is well established. These isotopic variations manifest themselves at the part per millions level or larger, in isotopes that were synthesised in various stellar environments. To escape homogenisation, some of these isotopic signatures must have been preserved in dust, which ended up being heterogeneously distributed in the solar protoplanetary disc. The origin of the nucleosynthetic heterogeneity is still poorly constrained, potentially reflecting inherited isotope variations from the Sun's parental molecular cloud and/or processing and redistribution during the subsequent protoplanetary disc phase with thermal processing and size sorting as major processes. This chapter aims to provide a broad review of the dynamical, collisional, and thermal processes in protoplanetary discs -- from initial infall to gas dispersal -- that may have influenced the distribution and survival of the anomalous carrier phases, which finally accreted into asteroids and planets. While several of these mechanisms have been considered in past studies, they are often examined in isolation, which impedes the assessment of how their effects may be altered or amplified by additional disc processes. Size sorting in particular has received little attention, and here we highlight that this process likely occurred in the disc and can induce nucleosynthetic heterogeneity. By placing previous studies within the context of a comprehensive overview, we aim to clarify the broader physical framework in which anomalous carrier transport occurs and identify previously underexplored mechanisms that may have contributed to the final isotopic structure of the Solar System we see today.

The extragalactic background Light (EBL) from ultraviolet to infrared comprises the emission from all stars, galaxies, and actively accreting black holes in the observable Universe. A precise measurement of the EBL is critically important to probe models of star formation and galaxy evolution. The EBL can be measured via the absorption imprint left on the spectra of gamma-ray blazars. In this work, we rely on 15 years of {\it Fermi}-LAT data and 1576 blazars to measure the EBL optical depth in the $0<z<4.3$ range. We detect the EBL attenuation with $\sim23\sigma$ significance and measure the optical depth in 19 redshift bins, extending the coverage and improving on our previous results. This allows us to reconstruct the EBL evolution and find general consistency with recent EBL models. These results represent the most precise determination of the EBL with GeV $\gamma$ rays to date.

We study cosmological solutions of (pseudo)scalar theories with periodic potentials, in the presence of arbitrary cosmological fluids -- including a cosmological constant of either sign. Independently of the initial misalignment angle and field velocity, we derive an analytic bound that the axion mass parameter and decay constant fulfill as the universe decreases its acceleration rate, finding a natural application in models of thawing quintessence. As a first application, we illustrate the analytic handle our bound provides in bounding axion dark energy, after observational inputs from DESI and various supernovae data sets are taken into account. As a second application, we argue that our analytic bounds in combination with proposed quantum gravity constraints on axions exclude vast regions of parameter space. The combined constraints push the axion masses to be much larger than the Hubble scale, in tension with basic models of axion quintessence.

We apply discrete Morse theory, global topology, and persistent homology to characterize the impact of massive neutrinos on the multiscale cosmic web, focusing on filaments. The topology of the cosmic web is sensitive to neutrino imprints, and persistence diagrams provide more information than commonly used summary statistics by quantifying the longevity of topological features across densities. This scale-adaptive, parameter-free formalism is powerful, as massive neutrinos affect halos, walls, filaments, and voids in distinct ways. Within this framework, we simultaneously assess their impact on tracers and skeleton structures and capture their multiscale signals across cosmic time. Discrete Morse theory is also well suited for particle-based neutrino implementations, often affected by Poisson shot noise, as it preserves the salient features of the underlying smooth field. Using two independent sets of N-body simulations, we present filament statistics and persistence diagrams in massive-neutrino cosmologies. Our results show that neutrinos leave distinct imprints on filaments and skeleton connectivity, producing mass-dependent signatures most pronounced at high redshift (z~2) and detectable at the few-percent level for masses as small as $M_\nu \sim 0.1$ eV. Filaments thus provide an ideal environment for isolating neutrino effects. We also compare two implementations of massive neutrinos to assess systematics. Our study establishes a promising avenue for leveraging cosmic web topology, persistent homology, and environment-based statistics to constrain or directly detect neutrino mass and infer the mass hierarchy - a long-standing challenge in particle physics and a major objective of ongoing and upcoming galaxy redshift surveys (e.g., DES, DESI, Euclid, Rubin-LSST).

The exploration of planetary bodies in our Solar system and beyond relies on the processing and interpretation of large, spatio-temporally inconsistent, and heterogeneous datasets. Recent advances in machine learning (ML) provide unprecedented opportunities to address many fundamental challenges posed by these heterogeneous and hyper-dimensional datasets. This review chapter highlights innovative ML methodologies that were developed and used by NCCR PlanetS members to address three overarching challenges in (exo)planetary science. The first challenge is sequence modelling, which encompasses the intricate analysis of one-dimensional data such as time series of radial velocities and light curves, among other examples. Secondly, there is pattern recognition that involves studying correlations, leveraging convolutional neural networks for feature extraction, mapping and cross correlation among other examples., anomaly detection through variational autoencoders, and unsupervised clustering of mass spectrometric data. Lastly, there are generative models and emulation-based Bayesian analysis, which encompass the development of predictive models for planetary interior structure, employing Deep Neural Networks to understand planet formation mechanisms. These innovative ML methodologies herald a paradigm shift in the processing of data and numerical models that represent inherent challenges in planetary and exoplanetary science, paving the way for revolutionary discoveries and ideas in this field.

We present an extended optical monitoring of the quadruply-imaged gravitationally lensed quasar QSO 2237+0305, the Einstein Cross, including observations from different observatories in both hemispheres and using a new photometric technique. This technique uses a region far enough from the lens system to determine accurately the sky background level, and minimises contamination from the lensing galaxy by combining analytical and numerical modeling of its structure. The resulting light curves of the four quasar images describe variations across practically the entire optical spectrum and span about 9000 days in the $VRI$ bands. The multi-band microlensing variability is captured with an unprecedented level of detail, and a preliminary microlensing analysis reveals an almost linear scaling of source radius with wavelength, providing direct evidence for the wavelength-dependent structure of the region contributing to optical passband fluxes. Specifically, assuming a mean microlens mass $\langle M \rangle$ = 0.3 $\rm{M_{\odot}}$ and concentric Gaussian sources that move according to the velocity distribution peaks (speed and direction) reported in a previous microlensing analysis, we find that the half-light radius of the $g$-band source is 9.6 $\pm$ 2.7 lt-day and the size of the sources grows with wavelength with a power-law index of $\alpha$ = 0.94 $\pm$ 0.05. We conclude that these long-term light curves set stringent empirical constraints on models of quasar emission and microlensing physics.

JWST has revealed a large population of massive black holes (BHs) in the early Universe with unusual properties which mark them as distinct from low-redshift active galactic nuclei. Such findings have prompted the development of new models of BH formation and growth, and of their co-evolution with host galaxies. Linking the gas-phase metallicity of BH environments to seed masses is key to understanding which evolutionary pathways could explain the population of JWST-discovered BHs. We present new high-resolution JWST NIRSpec/IFU observations covering the rest-frame optical emission lines of a Little Red Dot (LRD) at $z=3.55$, known as The Cliff, from the `Red Unknowns: Bright Infrared Extragalactic Survey' (RUBIES). We find evidence for low metallicity ($Z=0.017\pm0.004 \ Z_\odot$) based on the low narrow-line [OIII]$\lambda5007$/H$\beta$ ratio, supported by the non-detection of low-ionisation emission lines such as [OII]$\lambda\lambda3727,3729$ and [NII]$\lambda\lambda6548,6583$. We find that the observed properties of The Cliff, including its overmassive BH, can be reproduced by some simulations of black hole growth and evolution down to $z\sim3.5$. However, these simulation runs require high seed masses ($10^4 - 10^5\ M_\odot$) and appear as rarely in the simulation volume as in the RUBIES survey volume over redshifts $3<z<4$, highlighting the unusual nature of The Cliff. Future simulations and numerical models will help to uncover how such a metal poor system managed to develop a massive black hole and persist to such low redshift.

Stacked prism lenses (SPLs) are a type of refractive X-ray optics currently under development with the potential to greatly improve on current X-ray telescope optics in terms of focal length, angular resolution, efficiency and scalability. For this work, SPLs are manufactured using two-photon polymerization (2PP), with production being significantly faster and with higher geometric fidelity than previous methods. Preliminary laboratory tests show improved efficiency compared to previous manufacturing methods and promising optical capabilities. Two-photon polymerization is shown to be a reliable method for producing SPLs, and when challenges around printing time and assembly are addressed, the path towards an SPL X-ray telescope lies open.

In this chapter, we summarize the underlying numerical methods needed for efficient $N$-body integration of planetary systems. We discuss how symplectic integrators have been developed to tackle the complementary problems of long-term orbital integration and short-term collisional interactions. The public code GENGA, a parallel GPU/CPU planet formation and orbital dynamics simulation code, was developed to unify these methods and take full advantage of the newest available computing hardware. We present state-of-the-art N-body simulations performed with GENGA in a comparative study regarding the basic properties that emerge during the late stages of the terrestrial planet formation process. We show that in modern N-body simulations the commonly used acceleration factor f, used to speed up the collisional growth of planets in simulations, should be avoided since it can lead to distorted chemical composition of the planets. We make a detailed comparison of low to high-resolution simulations, showing that the formation time scale depends on the size of the initial planetesimals. These simulations also show that terrestrial planets can form resonant chains without the need of orbital migration due to gas effects.

Oxygen and carbon-rich AGB stars - and objects directly polluted by them - are excellent laboratories to investigate the nucleosynthesis and mixing processes occurring during the later phases of the of low- and intermediate-mass star evolution. The determination of the abundances of several s-elements is a key tool for constraining theoretical AGB models. This contribution discusses the main results, recent advances, and current problems on this subject.

Utilising the optical imaging Fourier transform spectrograph SITELLE, the Star-formation, Ionized Gas and Nebular Abundances Legacy Survey (SIGNALS) is designed to study the connection between star-forming regions and their environments. Targeting $31$ local star-forming galaxies, its data products also lend themselves to planetary nebula (PN) surveys. We present here a new pipeline to find PNe using automated emission-line diagnostics and morphology tests, that is able to distinguish PNe from contaminants with an accuracy similar to that of past visual methods. We also perform thorough completeness tests using mock PNe inserted into the data cubes with full spectra. We apply these tools to a pilot sample of two dwarf irregular galaxies from the SIGNALS survey, NGC 4214 and NGC 4449, with other galaxies to follow. For these two galaxies, we identify $25$ PNe (including $6$ new discoveries) and $23$ PNe (including $13$ new discoveries), respectively, and calculate PN luminosity function distances of $3.09^{+0.25}_{-0.46}$ and $3.91^{+0.33}_{-0.52}$ Mpc, respectively, the latter consistent with previous estimates. We also calculate the bolometric PN specific frequency of our galaxies ($\alpha_\mathrm{bol}$), as well as a newly defined $V$-band PN specific frequency ($\alpha_\mathrm{V}$) based solely on the galaxies' total luminosities in that band.

We aim to provide a systematic and quantitative description of the hydrogenation network connecting HCN and HNC to methylamine on interstellar water ices, while identifying dominant pathways and bottlenecks. To this end, we performed a comprehensive quantum-chemical investigation of H-addition, H-abstraction, reactions with H2, and water-assisted H-transfer isomerization, covering intermediates linking HCN and HNC to CH3NH2. Calculations were carried out on amorphous solid water clusters of 14 molecules. Using benchmarked density functional theory, we derived activation barriers, elucidated mechanisms, and determined the binding energy distribution of H2CN and CNH2, also assessing deuterium substitution effects. H-addition reactions generally involve activation barriers, except for radical species. Considering both barrier heights and tunneling crossover temperatures, the most favorable sequence originates from HNC rather than HCN. The network evolves toward methanimine (H2CNH), the central species, or the singlet carbene HC:NH2, from which further hydrogenation leads to methylamine. Along these paths, several reactions are barrierless, while some H-abstraction processes compete with addition. Reactions involving H2 are uncommon, as most are endoergic. Deuterium substitution weakly affects classical barriers but significantly influences tunneling efficiencies. Our results support efficient formation of methanimine and methylamine from HNC on cold interstellar ices, with methanimine acting as a chemical sink, whereas HCN is less reactive and more likely to persist. These findings provide quantitative constraints for astrochemical models.

The 2025 UK National Astronomy Meeting (NAM) in Durham played host to a session titled "Unseen Astronomy", involving a variety of astronomy researchers in diverse fields. This unique meeting focussed on a number of novel projects exploring alternatives to purely visual means of display in Astronomy, encompassing spheres of education, communication and research, and straddling both accessible and general use applications. The successful inclusion of such a session at a major conference reflects the explosion of interest in multimodal astronomy in recent years, and hints at its transformative potential. Here, I aim to outline and motivate the topic of multi-modal science and consider its exciting potential. I will discuss this in the context of our own work in the area, the community building being undertaken to bring together researchers considering multi-modality, and efforts to impact astronomy at large.

Most detected transiting planets have orbits which would fit within the one of Mercury, exposing them to intense stellar irradiation and interactions that significantly alter their properties. In contrast, colder planets with longer orbital periods are less affected, offering crucial insights into their formation and migration histories. Characterizing transiting warm and temperate planets is a key missing piece in the exoplanet puzzle. Dedicated photometric and spectroscopic follow-up of transiting events detected in space-based photometric data opened the way to detecting long-period transiting exoplanets. The wealth of information available for these transiting planets makes them golden targets for in-depth characterization. For giant planets, combining precise masses, radii, and ages with state-of-the-art planetary evolution models allows the estimation of their planetary bulk compositions, a crucial element to explore their formation and evolution pathways. Furthermore, these planets are compelling candidates for hosting moons and circumplanetary rings-features that could illuminate dynamical histories, satellite formation processes, and even potential habitable environments.

Transit timing variations (TTV) are considered a tool for constraining the masses of transiting planets in the absence of radial-velocity data. Although theoretical studies have long revealed that TTV mass determinations intrinsically suffer from degeneracies, existing analyses of TTV data typically report a single-mode solution under a model with a specified number of planets. This is because fitting TTV curves in the high-dimensional solution space of TTV posterior is extremely challenging; even locating a single solution requires substantial computational resources. We developed an efficient mode-first searching algorithm that can locate multiple solutions in a single MCMC run. We applied this algorithm to Kepler-9 b and c, which have the highest-quality TTV data. We found that the observed TTV can be reproduced by many combinations of planetary masses spanning a broad range, rather than the previously assumed precise determination. The mass of Kepler-9 b can range from 31.6 to 47.1 $M_{\oplus}$, while that of Kepler-9 c can range from 21.8 to 32.3 $M_{\oplus}$, and even more broadly under looser constraints. These degenerate solutions follow a linear relationship under a tight mass ratio between the two planets, consistent with previous theoretical predictions. Furthermore, we demonstrate that achieving a globally converged posterior distribution for Kepler-9's TTV is impossible using a sampling algorithm that preserves the Markovian property. This underscores the need for caution when interpreting results from sampling algorithms that lack mathematical guarantees of global convergence.

Nitrogen is a key element for building habitable worlds, yet only a small fraction of the available N-budget of planet-forming disks has been detected. In particular, the lack of any IR NH$_3$ detection is striking, as this molecule is predicted to be rather abundant in the warm, inner regions of protoplanetary disks, and therefore potentially readily incorporated into (giant) planets' atmospheres. We present a combined modeling and observational study of N-bearing molecules in planet-forming disks, using detailed thermo-chemical disk models that investigate the sensitivity of N-containing molecules to the bulk elemental composition of the disk. Our models predict a strong increase in HCN flux with high C/H, and conversely a strong increase in flux from NO when O/H is high. The flux from NH$_3$ is not very sensitive to O/H, but does decrease at high C/H due to competition with HCN. However, the absolute NH$_3$ flux predicted by our model is not large enough to be detected with JWST-MIRI, even when N/H is enhanced by an order of magnitude. The flux from NO, on the other hand, is potentially detectable, and could therefore provide further insights into the N-budget of the inner disk. Using a cross-correlation technique, we search for NH$_3$ and NO detections in three disks, GW Lup, Sz 98, and V1094 Sco. We do not find any NH$_3$ detections, and only one tentative NO detection in V1094 Sco, though this needs further study to be confirmed. Additionally, we demonstrate that future facilities in the FIR may provide a better opportunity to detect NH$_3$ and thereby draw a comparison to the NH$_3$ budget known to be present in interstellar ices.

We present ALMA detections (or stringent upper limits) of the [CII] 158 $\mu m$ emission line and underlying dust continuum from five massive quenched galaxies (QGs) at 2<z<4.7. We find extreme variations in the molecular gas fractions ($\rm{f_g=M_{mol}/M_{\star}}$), spanning 0.1%-25%, if a standard $\rm{\alpha_{[CII]}}$ applies. We attempt a first empirical calibration of $\rm{\alpha_{[CII]}}$ with respect to dust continuum in a $z=2$ lensed QG and with respect to CO(3-2) in a $z=3.1$ QG, finding no evidence of strong deviations from the standard value. Dust continuum measurements, coupled with JWST/MIRI fluxes, suggest higher dust temperatures compared to expectations from $z<2$ QGs, reaching $T_{d}\sim40-50 \,K$ in two galaxies. Coupled with remarkably high total infrared luminosities (LIR) not explained by observed JWST colors not by energy balance based on literature dust extinction measurements, and with [CII] deficits down to $\rm{[CII]/LIR\sim 2\times10^{-4}}$ typical of (Ultra)Luminous Infrared Galaxies, our findings point to additional dust-heating mechanisms other than dust-absorbed stellar radiation. Surprisingly, JWST/NIRCam and ALMA imaging reveal widespread disturbed stellar morphologies and offsets/tails in dust and gas, indicative of ongoing interactions. While larger samples are needed to assess how common these features are in high-z QGs, these findings support a merger-driven origin for the phenomenology observed in these systems, with key similarities with respect to local post-starburst galaxies where low-velocity shocks and turbulence also inject energy into the residual ISM.

The relationship between the star formation rate surface density and the molecular gas surface density in galaxies is key to understanding galaxy evolution. To investigate the molecular Kennicutt-Schmidt (K-S) relation and its dependence on gas density, we analyze a uniform sample of 36 nearby galaxies from the AMISS survey, focusing on the CO(1-0), CO(2-1), and CO(3-2) transitions, which trace progressively denser and warmer molecular gas. Using statistical methods that combine binning with Markov Chain Monte Carlo (MCMC) fitting, we derive the slope, scatter, and intercept of the $\Sigma_{\mathrm{SFR}}$-$\Sigma_{\mathrm{CO}}$ relation for each transition. We find power-law slopes of 1.26, 1.14, and 1.07 for CO(1-0), CO(2-1), and CO(3-2), respectively, consistent with a trend toward increasingly linear star formation relations at higher-J transitions. This behavior supports the idea that denser gas is more directly linked to ongoing star formation and is consistent with previous findings of near-linear correlations between HCN or high-J CO luminosities and global SFR. The observed trend suggests an underlying relation between gas and SFR volume densities with a power-law index of $\sim$1.5, indicating enhanced star formation efficiency in denser environments. These findings underscore the critical role of dense gas in regulating star formation and highlight the importance of tracer selection and excitation conditions when interpreting the K-S relation across different environments.

Many exoplanetary systems are multiplanet configurations whose long-term dynamics are governed by N-body gravitational interactions. Consequently, their detection signatures cannot be adequately described by Keplerian orbits. Accurately interpreting the observational data of these systems -- including radial velocity (RV), astrometry, and transit timing variations (TTVs) -- requires N-body integration. Motivated by this need, we developed a Bayesian fitting framework that couples N-body integration with Markov chain Monte Carlo (MCMC) to retrieve the system parameters of multiplanet systems. The code, named \texttt{Nii-body}, integrates an adaptive Runge--Kutta--Fehlberg 7(8) (RKF78) solver with an automated parallel tempering MCMC algorithm. Using simplified synthetic astrometric observations, we evaluated the efficiency and robustness of \texttt{Nii-body}'s N-body orbit retrieval on an idealized two-planet model, demonstrating its potential for future application to real observational data. The N-body fitting workflow can be readily extended to RV, TTVs, or combined datasets, providing a versatile engine for high-precision orbital inference in multiplanet systems.

This chapter reviews the current state of observational and theoretical efforts in the characterization of exoplanet atmospheres, with a focus on developments enabled through the Swiss National Centre for Competence in Research (NCCR) PlanetS. It covers the essential physical and chemical processes that govern atmospheric dynamics, radiative transfer, chemistry, and cloud formation in exoplanets and brown dwarfs. The review discusses the modeling approaches used to simulate these processes, ranging from simplified 1D models to fully coupled 3D general circulation models. Atmospheric retrieval frameworks are presented as tools for inferring atmospheric properties from observational data, highlighting both classical Bayesian techniques and emerging machine learning methods. Observational strategies using instruments like HST, JWST, and ground-based high-resolution spectrographs are also examined. Special emphasis is placed on the interplay between theory and observation, and how developments in modeling, data analysis, and instrumentation collectively advance our understanding of planetary atmospheres beyond the Solar System.

The extragalactic background light (EBL), the cumulative radiation from all extragalactic sources, traces galaxy formation and cosmic evolution. High-energy $\gamma$ rays attenuated via pair production with EBL photons are a powerful probe of the EBL. In this work, we use very-high-energy (VHE; $E_\gamma > 100\,\mathrm{GeV}$) $\gamma$ rays to measure the local EBL intensity and test its consistency with galaxy counts and direct measurements. Our analysis employs a sample of 268 spectra from 45 sources observed with Imaging Atmospheric Cherenkov telescopes. A model-dependent study shows seven EBL templates require only $\le 10\%$ rescaling to fit the observed $\gamma$-ray attenuation. The galaxy-count-anchored model gives the closest match. We then derive template-marginalized TeV optical depths from a representative model subset. We combine them with \textit{Fermi}-LAT GeV measurements to reconstruct the EBL at $z = 0$ using empirical and physically motivated models. The two reconstructions agree and follow the integrated galaxy light to within $2$--$3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$ (typically $<25\%$) over $0.5$--$30\,\mu$m. Both are consistent with low-zodiacal-light observations, including outer solar system and dark cloud measurements. In contrast, the near-IR excess reported by IRTS and CIBER exceeds our reconstructed intensity by $3$--$5\sigma$, implying an additional $\gtrsim 5$--$10\,\mathrm{nW\,m^{-2}\,sr^{-1}}$ incompatible with the $\gamma$-ray optical depths. Combined with GeV constraints on EBL evolution to $z \simeq 4$, these TeV optical depths provide a VHE-anchored determination of the local EBL intensity. The agreement with galaxy counts and deep-space measurements indicates that known galaxy populations account for most of the optical and near-IR background, leaving limited room for an additional diffuse component.

''Little Red Dots'' (LRDs) are broad-line sources at high redshift, initially identified by their compact morphologies, red colours and prominent Balmer breaks. The origin of their optical-to-near-infrared continua is debated, with proposed explanations ranging from direct recombination emission to thermalised blackbodies from stellar-like atmospheres. Here we report evidence for Paschen jumps in a subset of LRDs, consistent with free-bound recombination to hydrogen $n=3$. The Paschen and Brackett continuum shapes across the sample are consistent with minimally reddened emission from low-temperature gas with $T_e\lesssim10\,000$ K, while the presence of Paschen jump signatures limits scenarios in which the emission is thermalised. Further, the extreme H$\alpha$ equivalent widths and the tight observed correlation between H$\alpha$ and the continuum follow naturally if both originate in recombination emission. This provides an observational upper limit on the contribution of any direct AGN accretion component and any stellar-atmosphere-like component, as well as on the fraction of line emission that can be thermalised as it traverses the cocoon. Ultimately, nebular radiative-transfer models provide a self-consistent explanation of the continuum, line strengths and line profiles without requiring multiple separately fitted components.

Gamma-ray astronomy from hundreds of GeV to PeV is confined to ground-based experiments that detect air showers induced by $\gamma$-rays entering Earth's atmosphere. While particle detector arrays feature huge detection areas, accurately reconstructing the primary particle properties is difficult due to the sparse sampling of the air shower and its intrinsic fluctuations. In this work, using simulations of a future water-Cherenkov array, we investigate two end-to-end deep learning approaches based on the transformer architecture with different computational complexities that utilize calibrated raw data. We benchmark both methods against well-established methods in the field in terms of $\gamma$-hadron separation, angular, core, and energy reconstruction. Our results show significant improvements across the whole energy range, particularly at low and intermediate energies. This work is the first to consistently demonstrate improved performance in both event reconstruction and $\gamma$-hadron separation using a single architecture.

The 1D flux power spectrum ($P_{\mathrm{1D}}$) of the Ly$\alpha$ forest provides an exceptionally high-resolution probe of structure formation down to small scales ($k\approx1-10~\text{$h~$Mpc$^{-1}$}$). These scales carry the imprints of massive neutrinos, warm dark matter, and the running of the primordial power spectrum spectral index. The effective field theory (EFT) is a promising perturbative approach to systematically and efficiently describe the Ly$\alpha$ forest, but it faces challenges in its application to $P_{\mathrm{1D}}$, as many EFT parameters become degenerate when projected along the line of sight. In addition, this projection generates new stochastic terms from the integration over small-scale modes. In this work, we address these issues by compressing the EFT model space using the Fisher matrix formalism and linearizing the resulting compression directions, enabling analytic template marginalization and significantly reducing the computational cost of likelihood evaluation. We use hydrodynamical simulations to obtain a baseline estimate of EFT parameters, and use the DESI DR1 $P_{\mathrm{1D}}$ measurements to derive compression directions. We then marginalize over deviations from the baseline using these compression directions and forecast the constraining power of our formalism. We find that even in conservative scenarios where each data redshift bin requires its own set of EFT parameters, the cosmological constraints saturate with the linear bias, two leading-order 1D stochastic terms, and three principal combinations of the remaining EFT templates. In this case, our forecasted precision of the amplitude ($\Delta^2_p$) and the logarithmic slope ($n_p$) of the linear matter power spectrum at the pivot scale ($k_p=0.7~\text{Mpc}^{-1}$) is $10\%$ and $2.0\%$, respectively, which is similar to emulator-based analyses that include observational data systematics.

This study evaluates the effect of proposed constellations -- ranging from current deployments to mega-constellations and very bright reflector concepts -- on direct trail losses, diffuse background, and scattered sky brightness. We use a numerical model for Mie and Rayleigh scattering in the V band, adapted from moonlight sky-brightness calculations and validated against observations of moonlight and stellar background light. This is combined with the SatConAnalytic package to quantify scattered light, diffuse light from undetected satellites, and direct losses from detected trails. Constellations comprising approximately 60,000 satellites that adhere to the V_550km > 7 recommendation exert a negligible effect on sky brightness, contributing only about 10^-4 of the natural dark sky. Conversely, mega-constellations with 10^6 satellites render trails pervasive. Bright satellites, such those from AST SpaceMobile, significantly impact saturating detectors even when their number is moderate. Extremely bright satellites pose a far more severe threat: a 5000-satellite Reflect Orbital-like constellation elevates the scattered sky background by 20%-30%, and a population of 50,000 increases it by 200%-300%. The constellations currently proposed for launch, over 1,700,000 objects and including satellites brighter than V_550km = 7, would substantially degrade observations. Maintaining satellite brightness below V_550km = 7 is important for all instruments, but critical for safeguarding saturating instruments, such as the VRO LSST camera and for limiting sky-background pollution. Even under this constraint, the total satellite population must remain below ~100,000 satellites to ensure that field-of-view losses do not exceed typical technical downtime.

Recent XRISM observations of active galactic nuclei such as PDS 456 have revealed ``forests'' of absorption lines best modeled by five distinct absorption zones with varying large blueshifts. We propose a model in which these relativistic blueshifts originate from the motion of the accretion disc itself, rather than from a clumpy super-Eddington outflow at hundreds of gravitational radii $r_g\equiv GM/c^2$. We demonstrate that thin rings of absorbing material lying just above the accretion disc at varying radii can produce the observed energy shifts and separations of the absorption zones. In this model, the PDS 456 transmission spectrum is well reproduced by rings with widths $\Delta r\lesssim1r_g$ at locations between the black hole's innermost stable circular orbit (ISCO) and $\approx15r_g$. This model suggests that the absorption forests seen in XRISM observations can probe the surface structure of the innermost ($\lesssim15r_g$) regions of quasar accretion discs.

The discovery of circumbinary planets (CBPs) has advanced our understanding of planet formation and dynamical evolution in complex environments. However, the population of such planets remains small, leading their underlying physical properties to be loosely constrained. In this work, we have developed a semi-automated framework to identify planetary transit events in light curves of eclipsing binaries observed by the Transiting Exoplanet Survey Satellite (TESS). Our search method, ${\tt mono-cbp}$, removes stellar eclipses and applies a custom detrending procedure, searching for individual transit events and applying automated vetting procedures to filter false positive signals. We searched a sample of binaries from the TESS Eclipsing Binary Catalogue, yielding one candidate transit event. ${\tt mono-cbp}$ was also tested on the known population of transiting CBPs, using the Kepler long-cadence photometry for the Kepler transiting CBPs and the TESS Full Frame Image photometry for the TESS CBPs. Excluding transits that are shallower than the intrinsic noise of the Kepler/TESS data, ${\tt mono-cbp}$ achieved a recovery rate of $\geq50$ per cent for each planet, reaching >75 per cent for 9 of the 14 planets. To test the limits of our framework, we injected simulated transit profiles with varying depth and duration into our sample of TESS light curves, finding that our recovery rate is a strong function of transit duration and the metrics used to filter false positive signals. This framework may be applied to large samples of TESS eclipsing binaries with little computational burden and to photometry from future space-based photometric surveys.

The Wide-field Spectroscopic Telescope (WST) is a proposed 12-meter segmented facility optimized for seeing limited observations in the visible and designed to operate both a high-multiplex multi-object spectrograph and a panoramic integral field spectrograph (IFS). The WST IFS concept builds on instruments such as MUSE at the VLT (Very Large Telescope), using field splitters and image slicers to reformat a large field into pseudo-slits feeding spectrographs with two optimized spectral channels. This paper presents the spectrograph architecture developed for the WST IFS, aiming to achieve high through put and image quality over a wide wavelength range in a cost-effective manner. We investigate the use of curved detectors as a means to simplify the spectrograph layout, reduce aberrations, and potentially improve efficiency. This study establishes a promising baseline for the IFS spectrographs and assesses the benefits of incorporating curved sensors that can guide the development of future large-scale integral field spectrographs.

Low-mass galaxies with intense starbursts exhibit spectra dominated by extreme nebular emission and faint stellar continua. These extreme emission-line galaxies (EELGs) are key laboratories to study star formation, feedback, and ionizing photon escape in low-metallicity environments. We exploit the DESI survey to assemble the k-Means of Extreme Nebulae from DEsi outLiers (k-MENDEL), a statistically robust sample of ~16,000 EELGs at 0.01 < z < 0.96 selected via automatic k-means classification. Using SED fitting and Te-based metallicities, we characterize EELGs including "blueberry" and "green pea" galaxies, spanning stellar masses of 10^6-10^10 Msun and SFRs of 0.1-100 Msun/yr. k-MENDEL extends previous SDSS samples toward higher redshifts and lower metallicities (12+log(O/H) ~ 7.0-8.5). EELGs lie systematically above the star-forming main sequence, with sSFRs up to ~100 Gyr^-1. They follow a shallower mass-metallicity relation offset by 0.3-0.5 dex from local relations, closely resembling young galaxies observed with JWST at z > 3-10. The large intrinsic metallicity scatter, even after projecting along the fundamental metallicity relation, indicates strong departures from simple "bathtub" models, suggesting massive inflows of metal-poor gas followed by strong feedback. While ~6% of the sample shows AGN-like signatures, the most extreme star-forming systems reach high ionization (O32 ~ 5-60) comparable to confirmed Lyman-continuum emitters. Our results support the interpretation of EELGs as short-lived, non-equilibrium phases in the evolution of low-mass galaxies and highlight their importance as nearby analogs of galaxies likely driving cosmic reionization (Abridged).

The vertical redistribution of materials in the lunar regolith - ranging from continuously produced space-weathering products to sporadic pulses of supernova- or kilonova-derived isotopes - remains a fundamental problem in planetary science. We present a unified stochastic model of regolith gardening induced by the impact flux. Treating gardening as a competition between impact-driven advection and diffusion predicts the maturity profiles of Apollo cores over more than two orders of magnitude in time ($1.4 \times 10^7$ to $4.5 \times 10^8$ years). This model describes well the depth profiles of live Fe60 in Apollo regolith samples, suggesting that supernova dust capture is independent of native iron abundance, and is consistent with a uniform influx at the latitudes of the Apollo landing sites. We extend our model to predict lunar signals for live r-process species that might originate from supernovae or kilonovae: Pu244 tied to terrestrial detections, and I129, Hf182, and Cm247 based on r-process calculations. The Pu244/Fe60 depth profile can probe the origin of Pu244, motivating searches in Artemis regolith samples down to depths O(100) cm.

Planet formation remains a fundamentally important yet poorly understood process. Protoplanetary disks, the birthplaces of planetary systems, exhibit a wide range of substructures that are increasingly interpreted as signatures of interactions with forming planets. However, the direct detection rate of protoplanets within these disks remains low, leaving critical gaps in our understanding of the physical mechanisms driving their formation and early evolution. In this chapter, we review recent efforts by the high-contrast imaging community to directly observe forming protoplanets and their immediate environments. These observations aim to provide key constraints on thermal and accretion processes, planetary growth, and the formation of circumplanetary disks and satellite systems. We also propose a path forward for deriving observational estimates of the planet mass-to-radius ratio ($M_p/R_p$), a crucial parameter for distinguishing between competing formation models and understanding the thermal evolution of young planets. Finally, we highlight how upcoming instruments on the Extremely Large Telescope (ELT), with their unprecedented combination of high spatial and spectral resolution, will transform our ability to probe planet formation at the smallest and most critical scales.

We combine the highest resolution N-body simulation of a $\sim 10^{12}\, M_\odot$ $\Lambda$CDM halo (Aquarius-A) with the {\sc GALFORM} galaxy formation semianalytic model to study the full satellite population expected in a MW-like galaxy. The model assumes that galaxies only form in subhalos whose peak circular velocity exceeds the H-cooling threshold, all of which are well resolved in the simulation. The number of luminous subhalos ever accreted into the main halo is thus well defined, and implies that the total number of MW satellites, down to arbitrarily low luminosity, should not exceed a few hundred. The model tracks satellites even after their halos cease to be resolved ("orphan" galaxies), and includes a novel treatment of dark matter and stellar tidal stripping which takes into account that all $\Lambda$CDM subhalos survive until the present because of their cuspy inner density profiles. After accounting for tides, our results match well the massive end of the observed MW satellite mass function and predict that a large number of ultra-faint dwarfs are missing from the current MW satellite census. The missing UFDs are predicted to avoid the innermost regions of the host, and to have properties that overlap with those of the many ultra-faint compact MW satellites (UFCSs) discovered recently, with properties intermediate between globular clusters and dwarf galaxies. Our results suggest that many UFCS systems are dark matter-dominated dwarfs with velocity dispersions between $1-3$km/s, which have survived disruption because they reside in the dense cusp of $\Lambda$CDM subhalos. UFCSs should have mean densities of order $10^{10}$-$10^{11}\,M_\odot/$kpc$^3$, higher than those of more extended ultra-faint systems. If confirmed, our results would provide support for the cuspy nature of $\Lambda$CDM dark matter halos and for the hydrogen-cooling threshold for galaxy formation.

The inverse energy cascade in turbulent Taylor-Couette flow is studied in line with the results of the large eddy simulation. The simulation results show that the inverse energy cascade first occurs within the core region of the flow channel of the Taylor-Couette flow at higher Reynolds number. It is uncovered that this phenomenon is induced by the pulsed zero shear stress resulting from the singularities of the Navier-Stokes equation. In the core area between the two cylinders, the shear stress is nearly zero at higher Reynolds number. The turbulence generated there has high turbulent energy due to discontinuity of the tangential velocity. Since the energy transfer between the fluid layers is inhibited due to the low shear stress, the turbulent energy cannot be transferred along the radial direction, and small-scale vortices with high turbulent energy are produced. These small-scale vortices are located with the large-scale vortices and cannot be dissipated owing to low shear stress. A peak in the energy spectrum at middle frequency (or wave number) is formed due to the concentration of the small-scale vortices. As the number of the singular points of the Navier-Stokes equation increases with the increasing Reynolds number, the region with zero shear stress expands along the radial direction, intensifying nonlinear instability and energy accumulation. This, in turn, leads to more prominent peaks in the energy spectrum, resulting in a more pronounced inverse energy cascade.

Net primary productivity (NPP) forms the basis of biological carbon pump, but its estimates in high-latitude regions remain highly uncertain despite its disproportional importance for the global carbon sink. Optical satellites are limited by cloud cover, low irradiance, and shallow light penetration, with uncertainties further exacerbated by the lack of in situ validations and regional model tuning for NPP measurements. This study compared two satellite-based models, a global (VGPM) and a regionally tuned (BIO) NPP model, with a time series of in situ NPP. Using a high-frequency, depth-resolved moored profiler in the subpolar Northwest Atlantic (56°N) in 2016, in situ NPP was estimated by daily bio-optical profiles and prior measurement of photosynthesis-irradiance (P-I) parameters. Our findings indicated that satellite-derived estimates of depth-integrated NPP were overestimated by a factor of 2.5 to 4. However, the reasons for the discrepancies varied between the VGPM and BIO model. VGPM used global photosynthetic parameters with a simplified depth assumption, leading to an unrealistic vertical structure for depth-integrated NPP, despite its surface values were lower than in situ estimates. A major phytoplankton bloom in June-July was missed by VGPM, likely due to the use of non-regionally calibrated OCI Chl-a, which led to an underestimation of biomass. In contrast, the BIO model used regionally tuned POLY4 Chl-a products, and the differences in the assignment of P-I parameters accounted for the remaining discrepancies. This study showed the possibility to reach good agreement between satellite and in situ NPPs if the challenge of P-I assignment can be overcome. We recommend further studies to investigate discrepancies of NPP estimates in high-latitude regions, focusing on data sources and model choices, as well as improving regional model calibration to enhance NPP accuracy.

The Love numbers of a gravitating body are response coefficients encoding its tidal deformability. In compact binary systems, they appear in the gravitational waveform during the inspiral phase and will be measurable by upcoming gravitational-wave observatories. This review provides a comprehensive and pedagogical account of the theoretical foundations of Love numbers and surveys the most recent advances in the study of tidal effects in compact objects, with particular emphasis on black holes and neutron stars. We begin with a gentle introduction to tidal effects in Newtonian gravity, leading into a discussion of how to robustly define tidal responses in general relativity using the effective field theory and post-Newtonian frameworks. After an overview of the perturbation theory of black holes and neutron stars, we review the computation of Love numbers and dissipative response coefficients in a wide range of settings, including the static, dynamical, and nonlinear tidal responses of Kerr black holes and neutron stars in four-dimensional general relativity. We further discuss the extension of these results to charged black holes and a number of "new physics" scenarios, including higher-dimensional black holes and other black objects, (anti-) de Sitter black holes, supergravity black holes, and theories beyond general relativity. Finally we provide an overview of the tantalizing zoo of hidden symmetries of general relativity that have been uncovered in the attempt to explain the famous vanishing of static black hole Love numbers.

We study the loss of hyperbolicity of perturbation equations for black hole solutions of scalar Gauss-Bonnet gravity. We consider a class of coupling functions allowing for static black hole solutions with arbitrary small masses. For masses below a minimum value, such solutions become unphysical, because the perturbation equations become elliptic; this arguably corresponds to the loss of validity of the effective field theory. We analyse the dependence of this minimum mass on the parameters of the theory, finding that with an appropriate choice of the coupling function, such mass can be chosen arbitrarily small. However, this does not correspond to larger deviations from general relativity, since observable quantities like the black hole scalar charge are bounded by above.

The tidal response of compact objects provides a powerful probe of their internal structure and of the surrounding gravitational field. We provide a comprehensive and unified overview of tidal effects in black holes, neutron stars, and exotic compact objects, with emphasis on both static and dynamical responses to external fields, encoded in Love numbers and dissipation numbers. We discuss the vanishing of static bosonic Love numbers for black holes in vacuum General Relativity, their modifications in alternative theories, in non-standard models of compact objects, and in the presence of matter, as well as their role in testing deviations from Einstein's theory and environmental effects. A fundamental distinction between bosonic and fermionic perturbations is highlighted, as the latter yield nonzero static Love numbers even for black holes in General Relativity. For neutron stars, we overview the dependence of tidal Love numbers on the equation of state, the emergence of quasi-universal relations, and the impact of rotation, nonlinearities, and dynamical effects. Exotic compact objects typically feature nonvanishing static tidal Love numbers -- a striking observational signature that differentiates them from black holes. We further review how tidal effects influence the gravitational-wave signals from binary inspirals, and explore their implications for gravitational-wave astronomy. In particular, we stress their significance for current and future detectors as tools to test General Relativity, constrain the nuclear equation of state, and probe the fundamental nature of compact objects and their environments.

Gravitational-wave detectors are affected by short-duration non-Gaussian noise transients, commonly referred to as glitches, which can obscure astrophysical signals and complicate downstream analyses. While recent work has demonstrated the effectiveness of deep learning models for glitch classification using image-based time-frequency representations, comparatively less attention has been given to systematic evaluations of machine-learning architectures operating directly on tabular glitch metadata. In this work, we present a comprehensive benchmark of classical and deep learning models for multiclass glitch classification using numerical features derived from the Gravity Spy dataset. We compare gradient-boosted decision trees with a diverse set of neural architectures, including multilayer perceptrons, attention-based models, and neural decision ensembles, and evaluate them in terms of classification performance, inference efficiency, parameter efficiency, data-scaling behavior, and cross-model interpretability alignment. We find that while tree-based methods remain strong baselines for tabular data, several deep learning models achieve competitive performance with substantially fewer parameters and exhibit distinct inductive biases and scaling behavior. A cross-model attribution analysis further reveals partially consistent feature-importance hierarchies across architectures, providing new insight into interpretability structure in tabular models. These results clarify trade-offs between performance, complexity, data efficiency, and interpretability in tabular gravitational-wave analyses and provide practical guidance for deploying machine-learning models in detector characterization pipelines.

Fast and reliable inference of gravitational-wave source parameters is crucial for analyzing large catalogs that are reaching the size of hundreds of detections, and for identifying short-lived electromagnetic counterparts. Neural posterior estimation has emerged as a powerful inference method, where the model is trained on simulated gravitational-wave data at considerable computational cost, but thereafter enables extremely fast and inexpensive inference at test time. Here, we extend this approach by incorporating domain-specific physical insights and methods in the model architecture. These include compressing the data by heterodyning against a reference waveform chosen via approximate likelihood maximization, removing parameter degeneracies through tailored coordinate systems, and eliminating known multimodalities by folding the parameter space. As a result, the network is approximately equivariant to changes in the source parameters, and achieves a reduced training cost and improved model interpretability. Our implementation, called labrador, can be trained end-to-end on a 1-day timescale on $\sim 10^2$ CPU cores and a V100 GPU, achieving a median importance-sampling efficiency of 1% on quadrupolar, aligned-spin signals in a broad mass range (chirp mass $\mathcal M \in 1\text{-}50\,\mathrm{M}_\odot$, mass ratio $q > 0.1$). labrador is the first neural inference code to achieve extensive coverage of long-duration signals with secondary masses $m_2 < 10\,\mathrm{M}_\odot$, rendered possible by its equivariance property. Among our novel contributions is a numerically stable procedure that enables neural posterior estimation when the simulation and inference priors differ.

In arXiv:2603.18175, the authors argue, based on numerical studies of particular cases, that the quantum damping of cosmological shear in a modified loop quantum cosmological model (mLQC-I) that was recently found in arXiv:2510.14021 is not generic and that the universe never becomes truly classical. In this brief Note, we revisit these claims by carefully examining the underlying assumptions and the class of initial conditions considered. We show that the examples analyzed in arXiv:2603.18175 correspond to configurations that do not represent physically admissible collapsing Bianchi I universes, as they involve mixed expanding-contracting directions and lead to effectively lower-dimensional post-bounce geometries. Restricting to physically relevant initial conditions corresponding to genuine three-dimensional contraction, we find that the quantum damping of cosmological shear is a robust dynamical feature. This conclusion is supported by both numerical and perturbative analyses, which demonstrate that the post-bounce evolution admits an isotropic attractor, with anisotropies decaying exponentially and independently of the matter content, provided that the weak energy condition is satisfied. We further outline a plausible post-bounce mechanism for the onset of classicalization.

The motion of the Solar System with respect to the cosmic rest frame induces a kinematic dipole in the stochastic gravitational-wave background (GWB). Detecting this signal with space-based interferometers would provide an independent measurement of our peculiar velocity and a GW probe of cosmic anisotropies. We present a fully analytic derivation of the response of the \emph{Laser Interferometer Space Antenna} (LISA) to a kinematic dipole, and construct an optimal estimator for its detection. We show that the dipolar response is governed by a single frequency-dependent function fixed by symmetry, and we compute its behaviour across the LISA band. Using Fisher forecasts, we find that for a scale-invariant background detectability requires $h^2\Omega_{\rm GW} \gtrsim 5\times 10^{-8}$ for \emph{fiducial} LISA, and $h^2\Omega_{\rm GW} \gtrsim 5\times 10^{-10}$ for a detector with characteristic instrumental-noise amplitudes improved by an order of magnitude. Prospects are more favorable for signals with richer frequency profile. We also explore the potential of the kinematic dipole to break degeneracies, particularly in the presence of strong galactic foregrounds or noise features that closely mimic the primordial signal.

With the Event Horizon Telescope and future Very Long Baseline Interferometry arrays poised to image supermassive black holes, there is an urgent need to understand dynamic aspects of small-scale structure near the supermassive black hole. In this study, we introduce the relative magnification factor to characterize point sources distributed on the surface of the accretion disk near a black hole. We investigate the influence of source motion on this factor, comparing static sources with those corotating with the disk. In contrast to the static case, which can be well-understood in the standard framework of gravitational lensing, corotating sources exhibit significant distortions in the distribution of the magnification factor on both the image and source planes, indicating that the caustic structure is substantially modulated by source motion. This magnification factor pattern encodes signatures of the kinematics of accretion flow when the time-delayed images are incorporated. This potentially offers a novel probe for investigating the interplay between spacetime geometry and properties of accretion flow.

We propose a covariant, gauge-independent construction of foliation-based scalar-tensor theories, yielding diffeomorphism-invariant operators involving only gradients on the hypersurfaces where the scalar field is constant, assumed to be spacelike. This defines a basis of independent invariants up to four derivatives of $\phi$, including the first nontrivial parity-odd pseudoscalar at this order, with a straightforward extension to higher derivatives. Our framework goes beyond degenerate higher-order scalar-tensor (DHOST) theories and provides a nonlinear extension of U-DHOST (where $\nabla_\mu\phi$ is supposed to be timelike) directly in covariant form, without using unitary gauge as a starting point or imposing degeneracy a priori. After minimal coupling to gravity, we analyze the theory through its Hamiltonian constraint structure and linear cosmological perturbations about an FLRW background, and show that it propagates three physical degrees of freedom.

We develop a self-consistent analytical two-fluid framework for plasma evolution in the short-time regime, elucidating the fundamental mechanism underlying the coupled generation of flow and magnetic fields. We show that consistency between ion momentum and mass conservation imposes a structural constraint on the system: the total pressure must satisfy the Laplace equation, $\nabla^2 P = 0$. This constraint enables a class of exact analytical solutions in which pressure gradients simultaneously drive plasma flow and generate magnetic fields through a Biermann-type mechanism. Using representative parameters, we obtain magnetic-field strengths and flow velocities consistent with both laser-produced plasmas and large-scale astrophysical systems. This framework provides a unified description of pressure-driven magnetogenesis and plasma flow in the short-time regime.

We provide the first investigation of the solar production of symmetrons, a well-motivated class of screened scalar fields with density dependent couplings to the Standard Model, and their subsequent absorption in underground direct detection experiments. We compute the flux of symmetrons produced through photon conversion in the magnetic field of the solar tachocline, and constrain the resulting luminosity to not exceed 3% of the observed solar output. Even under the conservative assumption that production occurs only in the tachocline, this criterion yields robust astrophysical bounds on previously uncharted regions of symmetron parameter space, and predicts a keV-scale symmetron spectrum at Earth. We then derive the corresponding absorption signal in liquid xenon detectors, where symmetrons can interact with electrons through both conformal and disformal couplings. Using binned data from XENONnT, we obtain new direct-detection limits that are complementary to the solar luminosity constraint, further tightening the viable symmetron parameter space. Our results demonstrate that the Sun provides a testable, previously unexploited, source of symmetrons, and highlight that the interplay of astrophysical and laboratory searches offers a powerful strategy for probing screened scalar theories.

Measurement of $N_{\rm eff}$ in the CMB (Cosmic Microwave Background) observations, like Planck 2018 and BBN (Big Bang Nucleosynthesis) has already set stringent constraints on the interaction strength of light particles beyond the Standard Model (BSM). Despite such negligible couplings of such BSM particles to the visible sector, they are inevitably produced in the early universe through gravity-mediated processes. If a sizable density of light particles survives around CMB formation, they may act as dark radiation (DR) contributing to $N_{\rm eff}$ at CMB epoch. In this work, we study the production of such light BSM particles through the gravity-mediated scatterings in an effective field theory (EFT) setup assuming that all non-gravitational couplings of the BSM particle are negligible. Since the production is sensitive to the spin of the produced particle, we perform a concrete analysis for two representative cases: scalar dark Higgs DR and vector dark photons this http URL the Planck 2018 observations, we find constraints on the reheating temperature ($T_{\rm RH}$) and background equation of state ($w_\Phi$) during reheating in such scenarios featuring dark Higgs and dark photon. A comparative discussion involving gravity-mediated production of Dirac right-handed neutrinos ($\nu_R$) and light axion-like particles (ALP) is also presented. Finally, for completeness, we also analyze the scenario where the production occurs through a generic spin-2 mediator characterized by an effective scale $\Lambda$ delineating the parameter space that is currently ruled out from Planck-2018 and can be probed by the future CMB experiments like LiteBird, Simon Observatory, CMB-S4, CMB-HD.