The inference of astrophysical and cosmological properties from the Lyman-α forest conventionally relies on summary statistics of the transmission field that carry useful but limited information. We present a deep learning framework for inference from the Lyman-α forest at the field level. This framework consists of a 1D residual convolutional neural network (ResNet) that extracts spectral features and performs regression on thermal parameters of the intergalactic medium that characterize the power-law temperature-density relation. We trained this supervised machinery using a large set of mock absorption spectra from NYX hydrodynamic simulations at z = 2.2 with a range of thermal parameter combinations (labels). We employed Bayesian optimization to find an optimal set of hyperparameters for our network, and then employed a committee of 20 neural networks for increased statistical robustness of the network inference. In addition to the parameter point predictions, our machine also provides a self-consistent estimate of their covariance matrix with which we constructed a pipeline for inferring the posterior distribution of the parameters. We compared the results of our framework with the traditional summary based approach, namely the power spectrum and the probability density function (PDF) of transmission, in terms of the area of the 68% credibility regions as our figure of merit (FoM). In our study of the information content of perfect (noise- and systematics-free) Lyα forest spectral datasets, we find a significant tightening of the posterior constraints – factors of 10.92 and 3.30 in FoM over the power spectrum only and jointly with PDF, respectively – which is the consequence of recovering the relevant parts of information that are not carried by the classical summary statistics.

Auto-chemotaxis, the directed movement of cells along gradients in chemicals they secrete, is central to the formation of complex spatiotemporal patterns in biological systems. Since the introduction of the Keller-Segel model, numerous variants have been analyzed, revealing phenomena such as coarsening of aggregates, stable aggregate sizes, and spatiotemporally chaotic dynamics. Here, we consider general mass-conserving Keller-Segel models, that is, models without cell growth and death, and analyze the generic long-time dynamics of the chemotactic aggregates. Building on and extending our previous work, which demonstrated that chemotactic aggregation can be understood through a generalized Maxwell construction balancing density fluxes and reactive turnover, we use singular perturbation theory to derive the full rates of mass competition between well-separated aggregates. We analyze how this mass-competition process drives coarsening in both diffusion- and reaction-limited regimes, with the diffusion-limited rate aligning with our previous quasi-steady-state analyses. Our results generalize earlier mathematical findings, demonstrating that coarsening is driven by self-amplifying mass transport and aggregate coalescence. Additionally, we provide a linear stability analysis of the lateral instability, predicting it through a nullcline-slope criterion similar to the curvature criterion in spinodal decomposition. Overall, our findings suggest that chemotactic aggregates behave similarly to phase-separating droplets, providing a robust framework for understanding the coarse-grained dynamics of auto-chemotactic cell populations and providing a basis for the analysis of more complex multi-species chemotactic systems.

Context. Observational evidence has accumulated in recent years, showing that the Galactic bulge includes two populations, a metal-poor one and a metal-rich one, which in addition to having different metallicities show different alpha over iron abundances, spatial distribution, and kinematics. While the metal-rich, barred component has been fairly well characterized, the metal-poor, spheroidal component has been more elusive and harder to describe. RR Lyrae variables are clean tracers of the old bulge component, and they are, on average, more metal-poor than red clump stars. Aims. In the present paper, we provide a new catalog of 16488 ab-type RR Lyrae variables in the bulge region within |l|≲10° and |b|≲2.8°, extracted from multi-epoch Point Spread Function photometry performed on VISTA Variable in the Vía Láctea survey data. We used the catalog to constrain the shape of the old, metal-poor, bulge stellar population. Methods. The identification of ab-type RR Lyrae among a large sample of candidate variables of different types has been performed via a combination of a Random Forest classifier and visual inspection. We optimized this process in such a way to extract a clean catalog with high purity, although for this reason its completeness, close to the midplane, is lower compared to a few other near-infrared catalogs covering the same region of the sky. Results. We used the present catalog to derive the shape of their distribution around the Galactic center, resulting in an elongated spheroid with projected axis ratio of b/a~0.7 and an inclination angle of ϕ~20 degrees. We discuss how observational biases, such as errors on the distances and a nonuniform sampling in longitude, affect both the present measurements and previous ones, especially those based on red clump stars. Because the latter have not been taken into account before, we refrain from a quantitative comparison between these shape parameters and those derived for the main Galactic bar. Nonetheless, qualitatively, taking into account observational biases would lower the estimated ellipticity of the bar derived from red clump stars, and hence reduce the difference with the present results. Conclusions. We publish a high-purity RRab sample for future studies of the oldest Galactic bulge population, close to the midplane. We explore different choices for the period-luminosity-metallicity relation, highlighting how some of them introduce spurious trends of distance with either the period or the metallicity, or both. We provide evidence that they trace a structure that is less elongated than the main bar, though we also highlight some biases of these kind of studies not discussed before. ★Based on observations taken within the ESO VISTA Public Survey VVV, Program ID 179.B-2002.

Extragalactic and galactic cosmic rays scatter with the cosmic neutrino background during propagation to Earth, yielding a flux of relic neutrinos boosted to larger energies. If an overdensity of relic neutrinos is present in galaxies, and neutrinos are massive enough, this flux might be detectable by high-energy neutrino experiments. For a lightest neutrino of mass <inline-formula><mml:math display="inline"><mml:msub><mml:mi>m</mml:mi><mml:mi>ν</mml:mi></mml:msub><mml:mo>∼</mml:mo><mml:mn>0.1</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>eV</mml:mi></mml:math></inline-formula>, we find an upper limit on the local relic neutrino overdensity of <inline-formula><mml:math display="inline"><mml:mo>∼</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>13</mml:mn></mml:msup></mml:math></inline-formula> and an upper limit on the relic neutrino overdensity at TXS <inline-formula><mml:math display="inline"><mml:mrow><mml:mn>0506</mml:mn><mml:mo>+</mml:mo><mml:mn>056</mml:mn></mml:mrow></mml:math></inline-formula> of <inline-formula><mml:math display="inline"><mml:mo>∼</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>10</mml:mn></mml:msup></mml:math></inline-formula>. Future experiments like GRAND or IceCube-Gen2 could improve these bounds by orders of magnitude.

It is well established that maximizing the information extracted from upcoming and ongoing stage-IV weak-lensing surveys requires higher order summary statistics that complement the standard two-point statistics. In this work, we focus on weak-lensing peak statistics to test two popular modified gravity models, <inline-formula><tex-math id="TM0001" notation="LaTeX">$f(R)$</tex-math></inline-formula> and nDGP, using the FORGE and BRIDGE weak-lensing simulations, respectively. From these simulations, we measure the peak statistics as a function of both cosmological and modified gravity parameters simultaneously. Our findings indicate that the peak abundance is sensitive to the strength of modified gravity, while the peak two-point correlation function is sensitive to the nature of the screening mechanism in a modified gravity model. We combine these simulated statistics with a Gaussian Process Regression emulator and a Gaussian likelihood to generate stage-IV forecast posterior distributions for the modified gravity models. We demonstrate that, assuming small scales can be correctly modelled, peak statistics can be used to distinguish general relativity from <inline-formula><tex-math id="TM0002" notation="LaTeX">$f(R)$</tex-math></inline-formula> and nDGP models at the 2σ level with a stage-IV survey area of <inline-formula><tex-math id="TM0003" notation="LaTeX">$300$</tex-math></inline-formula> and <inline-formula><tex-math id="TM0004" notation="LaTeX">$1000 \, \rm {deg}^2$</tex-math></inline-formula>, respectively. Finally, we show that peak statistics can constrain <inline-formula><tex-math id="TM0005" notation="LaTeX">$\log _{10}\left(|f_{R0}|\right) = -6$</tex-math></inline-formula> per cent to 2 per cent precision, and <inline-formula><tex-math id="TM0006" notation="LaTeX">$\log _{10}(H_0 r_c) = 0.5$</tex-math></inline-formula> per cent to 25 per cent precision.

This research introduces an innovative application of physics-informed neural networks (PINNs) to tackle the intricate challenges of radiative transfer (RT) modelling in exoplanetary atmospheres, with a special focus on efficiently handling scattering phenomena. Traditional RT models often simplify scattering as absorption, leading to inaccuracies. Our approach utilizes PINNs, noted for their ability to incorporate the governing differential equations of RT directly into their loss function, thus offering a more precise yet potentially fast modelling technique. The core of our method involves the development of a parametrized PINN tailored for a modified RT equation, enhancing its adaptability to various atmospheric scenarios. We focus on RT in transiting exoplanet atmospheres using a simplified 1D isothermal model with pressure-dependent coefficients for absorption and Rayleigh scattering. In scenarios of pure absorption, the PINN demonstrates its effectiveness in predicting transmission spectra for diverse absorption profiles. For Rayleigh scattering, the network successfully computes the RT equation, addressing both direct and diffuse stellar light components. While our preliminary results with simplified models are promising, indicating the potential of PINNs in improving RT calculations, we acknowledge the errors stemming from our approximations as well as the challenges in applying this technique to more complex atmospheric conditions. Specifically, extending our approach to atmospheres with intricate temperature-pressure profiles and varying scattering properties, such as those introduced by clouds and hazes, remains a significant area for future development.

We study the theory of a scalar in the fundamental representation of the internal supergroup $SU(N|M)$. Remarkably, for $M=N+1$ its tree-level mass does not receive quantum corrections at one loop from either self-coupling or interactions with gauge bosons and fermions. This property comes at the price of introducing both degrees of freedom with wrong statistics and with wrong sign kinetic terms. We detail a method to break $SU(N|M)$ down to its bosonic subgroup through a Higgs-like mechanism, allowing for the partial decoupling of the dangerous modes, and study the associated vacuum structure up to one loop.

If primordial black holes (PBHs) of asteroidal mass make up the entire dark matter they could be detectable through their gravitational influence in the solar system. In this work, we study the perturbations that PBHs induce on the orbits of planets. Detailed numerical simulations of the solar system, embedded in a halo of PBHs, are performed. Using the Earth-Mars distance as an observational probe, we show that the perturbations are below the current detection limits and thus PBHs are not directly constrained by solar system ephemerides. We estimate that an improvement in the measurement accuracy by more than an order of magnitude or the extraction of signals well below the noise level are required to detect the gravitational influence of PBHs in the solar system in the foreseeable future.

We present constraints on the $f(R)$ gravity model using a sample of 1,005 galaxy clusters in the redshift range $0.25 - 1.78$ that have been selected through the thermal Sunyaev-Zel'dovich effect (tSZE) from South Pole Telescope (SPT) data and subjected to optical and near-infrared confirmation with the Multi-component Matched Filter (MCMF) algorithm. We employ weak gravitational lensing mass calibration from the Dark Energy Survey (DES) Year 3 data for 688 clusters at $z < 0.95$ and from the Hubble Space Telescope (HST) for 39 clusters with $0.6 < z < 1.7$. Our cluster sample is a powerful probe of $f(R)$ gravity, because this model predicts a scale-dependent enhancement in the growth of structure, which impacts the halo mass function (HMF) at cluster mass scales. To account for these modified gravity effects on the HMF, our analysis employs a semi-analytical approach calibrated with numerical simulations. Combining calibrated cluster counts with primary cosmic microwave background (CMB) temperature and polarization anisotropy measurements from the Planck2018 release, we derive robust constraints on the $f(R)$ parameter $f_{R0}$. Our results, $\log_{10} |f_{R0}| < -5.32$ at the 95 % credible level, are the tightest current constraints on $f(R)$ gravity from cosmological scales. This upper limit rules out $f(R)$-like deviations from general relativity that result in more than a $\sim$20 % enhancement of the cluster population on mass scales $M_\mathrm{200c}>3\times10^{14}M_\odot$.

Context. The forward modelling of galaxy surveys has recently gathered interest as one of the primary methods to achieve the required precision on the estimate of the redshift distributions for stage IV surveys, allowing them to perform cosmological tests with unprecedented accuracy. One of the key aspects of forward modelling a galaxy survey is the connection between the physical properties drawn from a galaxy population model and the intrinsic galaxy spectral energy distributions (SEDs), achieved through stellar population synthesis (SPS) codes (e.g. FSPS). However, SPS requires a large number of detailed assumptions on the constituents of galaxies, for which the model choice or parameter values are currently uncertain. Aims. In this work, we perform a sensitivity study of the impact that the variations of the SED modelling choices have on the mean and scatter of the tomographic galaxy redshift distributions. Methods. We assumed the PROSPECTOR-β model as the fiducial input galaxy population model and used its SPS parameters to build 9-bands ugriZYJHK_s observed-frame magnitudes of a fiducial sample of galaxies. We then built samples of galaxy magnitudes by varying one SED modelling choice at a time. We modelled the colour-redshift relation of these galaxy samples using the self-organising map (SOM) approach that optimally groups similar redshifts galaxies by their multidimensional colours. We placed galaxies in the SOM cells according to their simulated observed-frame colours and used their cell assignment to build colour-selected tomographic bins. Finally, we compared each variant's binned redshift distributions against the estimates obtained for the original PROSPECTOR-β model. Results. We find that the SED components related to the initial mass function, as well as the active galactic nuclei, the gas physics, and the attenuation law substantially bias the mean and the scatter of the tomographic redshift distributions with respect to those estimated with the fiducial model. Conclusions. For the uncertainty of these choices currently present in the literature and regardless of the applied stellar mass function based re-weighting strategy, the bias in the mean and the scatter of the tomographic redshift distributions are greater than the precision requirements set by next-generation Stage IV galaxy surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and Euclid.

Context. The ESO public survey VISTA Variables in the Vía Láctea (VVV) surveyed the inner Galactic bulge and the adjacent southern Galactic disk from 2009–2015. Upon its conclusion, the complementary VVV extended (VVVX) survey has expanded both the temporal as well as spatial coverage of the original VVV area, widening it from 562 to 1700 sq. deg., as well as providing additional epochs in JHK_s filters from 2016–2023. Aims. With the completion of VVVX observations during the first semester of 2023, we present here the observing strategy, a description of data quality and access, and the legacy of VVVX. Methods. VVVX took ~2000 h, covering about 4% of the sky in the bulge and southern disk. VVVX covered most of the gaps left between the VVV and the VISTA Hemisphere Survey (VHS) areas and extended the VVV time baseline in the obscured regions affected by high extinction and hence hidden from optical observations. Results. VVVX provides a deep JHK_s catalogue of ≳1.5 × 10⁹ point sources, as well as a K_s band catalogue of ~10⁷ variable sources. Within the existing VVV area, we produced a 5D map of the surveyed region by combining positions, distances, and proper motions of well-understood distance indicators such as red clump stars, RR Lyrae, and Cepheid variables. Conclusions. In March 2023 we successfully finished the VVVX survey observations that started in 2016, an accomplishment for ESO Paranal Observatory upon 4200 h of observations for VVV+VVVX. The VVV+VVVX catalogues complement those from the Gaia mission at low Galactic latitudes and provide spectroscopic targets for the forthcoming ESO high-multiplex spectrographs MOONS and 4MOST. ★Based on observations taken within the ESO VISTA Public Survey VVV and VVVX, Programmes ID 179.B-2002 and 198.B-2004, respectively.

We derive a cutting rule for equal-time in-in correlators including cosmological correlators based on Keldysh $r/a$ basis, which decomposes diagrams into fully retarded functions and cut-propagators consisting of Wightman functions. Our derivation relies only on basic assumptions such as unitarity, locality, and the causal structure of the in-in formalism, and therefore holds for theories with arbitrary particle contents and local interactions at any loop order. As an application, we show that non-local cosmological collider signals arise solely from cut-propagators under the assumption of microcausality. Since the cut-propagators do not contain (anti-)time-ordering theta functions, the conformal time integrals are factorized, simplifying practical calculations.

Hyperluminous infrared galaxies (HyLIRGs) are the rarest and most extreme starbursts and found only in the distant Universe (z ≳ 1). They have intrinsic infrared (IR) luminosities L_IR ≥ 10¹³ L_⊙ and are commonly found to be major mergers. Recently, the Planck All-Sky Survey to Analyze Gravitationally-lensed Extreme Starbursts project (PASSAGES) searched ~10⁴ deg² of the sky and found ~20 HyLIRGs. We describe a detailed study of PJ0116-24, the brightest (μL_IR ≈ 2.6 × 10¹⁴ L_⊙, magnified with μ ≈ 17) Einstein-ring HyLIRG in the southern sky, at z = 2.125, with observations from the near-IR integral-field spectrograph VLT/ERIS and the submillimetre interferometer ALMA. We detected Hα, Hβ, [N II] and [S II] lines and obtained an extreme Balmer decrement (Hα/Hβ ≈ 8.73 ± 1.14). We modelled the molecular-gas and ionized-gas kinematics with CO(3-2) and Hα data at ~100-300 pc and (sub)kiloparsec delensed scales, respectively, finding consistent regular rotation. We found PJ0116-24 to be highly rotationally supported (v_rot/σ_{0, mol. gas} ≈ 9.4) with a richer gaseous substructure than other known HyLIRGs. Our results imply that PJ0116-24 is an intrinsically massive (M_baryon ≈ 10^11.3 M_⊙) and rare starbursty disk (star-formation rate, SFR = 1,490 M_⊙ yr⁻¹) probably undergoing secular evolution. This indicates that the maximal SFR (≳1,000 M_⊙ yr⁻¹) predicted by simulations could occur during a galaxy's secular evolution, away from major mergers.

Over the past decade, advancement of observational capabilities, specifically the Atacama Large Millimeter/submillimeter Array (ALMA) and Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instruments, alongside theoretical innovations like pebble accretion, have reshaped our understanding of planet formation and the physics of protoplanetary disks. Despite this progress, mysteries persist along the winded path of micrometer-sized dust, from the interstellar medium, through transport and growth in the protoplanetary disk, to becoming gravitationally bound bodies. This review outlines our current knowledge of dust evolution in circumstellar disks, yielding the following insights: ▪ Theoretical and laboratory studies have accurately predicted the growth of dust particles to sizes that are susceptible to accumulation through transport processes like radial drift and settling. ▪ Critical uncertainties in that process remain the level of turbulence, the threshold collision velocities at which dust growth stalls, and the evolution of dust porosity. ▪ Symmetric and asymmetric substructures are widespread. Dust traps appear to be solving several long-standing issues in planet formation models, and they are observationally consistent with being sites of active planetesimal formation. ▪ In some instances, planets have been identified as the causes behind substructures. This underlines the need to study earlier stages of disks to understand how planets can form so rapidly. Theoretical and laboratory studies have accurately predicted the growth of dust particles to sizes that are susceptible to accumulation through transport processes like radial drift and settling. Critical uncertainties in that process remain the level of turbulence, the threshold collision velocities at which dust growth stalls, and the evolution of dust porosity. Symmetric and asymmetric substructures are widespread. Dust traps appear to be solving several long-standing issues in planet formation models, and they are observationally consistent with being sites of active planetesimal formation. In some instances, planets have been identified as the causes behind substructures. This underlines the need to study earlier stages of disks to understand how planets can form so rapidly. In the future, better probes of the physical conditions in optically thick regions, including densities, turbulence strength, kinematics, and particle properties, will be essential for unraveling the physical processes at play.

Axion quark nuggets (AQN) are hypothetical, macroscopically large objects with a mass greater than a few grams and sub-micrometer size, formed during the quark-hadron transition. Originating from the axion field, they offer a possible resolution of the similarity between visible and dark components of the Universe, i.e. Ω_DM ∼ Ω_visible and observed matter-antimatter asymmetry. These composite objects behave as cold dark matter, interacting with ordinary matter and resulting in pervasive electromagnetic radiation throughout the Universe. This work aims to predict the electromagnetic signature in large-scale structures from this AQN-baryon interaction, accounting for thermal and non-thermal radiations. We use Magneticum hydrodynamical simulations to describe the realistic distribution and dynamics of gas and dark matter at cosmological scales. We construct a light cone encompassing a 1.4 square degree area on the sky, extending up to redshift z = 5.4, and we calculate the electromagnetic signature across a wide range of frequencies from radio, starting at ν ∼ 1 GHz, up to a few keV X-ray energies. We find that the AQNs electromagnetic signature is characterized by global (monopole) and fluctuation signals. The amplitude of both signals strongly depends on the average nugget mass and the ionization level of the baryonic environment, allowing us to identify a most optimistic scenario and a minimal configuration. The signal of our most optimistic scenario is often near the sensitivity limit of existing instruments, such as FIRAS in the ν = [100-500] GHz range and the South Pole Telescope for high-resolution ℓ > 4000 at ν = 95 GHz. Fluctuations in the Extra-galactic Background Light caused by the axion quark nuggets in the most optimistic scenario can also be tested with space-based imagers Euclid and James Webb Space Telescope. In general, our minimal configuration is still out of reach of existing instruments, but future experiments might be able to pose some constraints. We conclude that the axion quark nuggets model represents a viable model for dark matter, which does not violate the canons of cosmology nor existing observations. A reanalysis of existing data sets could provide some evidence of axion quark nuggets if our most optimistic configuration is correct. The best chances for testing the model reside in 1) ultra-deep infrared and optical surveys, 2) future experiments to probe the frequency spectrum of the cosmic microwave background, and 3) low-frequency (1 GHz < ν < 100 GHz) and high-resolution (ℓ ≳ 10⁴) observations.

We develop parametrizations of eight of the lowest Born-Oppenheimer potentials for quarkonium hybrid mesons as functions of the separation <inline-formula><mml:math display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula> of the static quark and antiquark sources. The parameters are determined by fitting results calculated using pure <inline-formula><mml:math display="inline"><mml:mi>S</mml:mi><mml:mi>U</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>3</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> lattice gauge theory. The parametrizations have the correct limiting behavior at small <inline-formula><mml:math display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula>, where the potentials form multiplets associated with gluelumps. They have the correct limiting behavior at large <inline-formula><mml:math display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula>, where the potentials form multiplets associated with excitations of a relativistic string. There is a narrow avoided crossing in the small-<inline-formula><mml:math display="inline"><mml:mi>r</mml:mi></mml:math></inline-formula> region between two potentials with the same Born-Oppenheimer quantum numbers.

Context. Planets form in the disks surrounding young stars. The time at which the planet formation process begins is still an open question. Annular substructures such as rings and gaps in disks are intertwined with planet formation, and thus their presence or absence is commonly used to investigate the onset of this process. Aims. Current observations show that a limited number of disks surrounding protostars exhibit annular substructures, all of them in the Class I stage. The lack of observed features in most of these sources may indicate a late emergence of substructures, but it could also be an artifact of these disks being optically thick. To mitigate the problem of optical depth, we investigated substructures within a very young Class 0 disk characterized by low inclination using observations at longer wavelengths. Methods. We used 3 mm ALMA observations tracing dust emission at a resolution of 7 au to search for evidence of annular substructures in the disk around the deeply embedded Class 0 protostar Oph A SM1. Results. The observations reveal a nearly face-on disk (inclination ∼ 16°) extending up to 40 au. The radial intensity profile shows a clear deviation from a smooth profile near 30 au, which we interpret as the presence of either a gap at 28 au or a ring at 34 au with Gaussian widths of σ = 1.4_‑1.2^+2.3 au and σ = 3.9_‑1.9^+2.0 au, respectively. Crucially, the 3 mm emission at the location of the possible gap or ring is determined to be optically thin, precluding the possibility that this feature in the intensity profile is due to the emission being optically thick. Conclusions. Annular substructures resembling those in the more evolved Class I and II disks could indeed be present in the Class 0 stage, which is earlier than suggested by previous observations. Similar observations of embedded disks in which the high-optical-depth problem can be mitigated are clearly needed to better constrain the onset of substructures in the embedded stages.

Observed low-mass galaxies with nuclear star clusters (NSCs) can host accreting massive black holes (MBH). We present simulations of dwarf galaxies ($M_{\mathrm{baryon}} \sim 0.6 - 2.4 \times 10^8 \rm \, M_\odot$) at solar mass resolution ($0.5\rm \, M_\odot < m_{\mathrm{gas}} < 4 \rm \, M_\odot$) with a multi-phase interstellar medium (ISM) and investigate the impact of NSCs on MBH growth and nuclear star formation (SF). The Griffin simulation model includes non-equilibrium low temperature cooling, chemistry and the effect of HII regions and supernovae (SN) from massive stars. Individual stars are sampled down to 0.08 $\rm M_\odot$ and their non-softened gravitational interactions with MBHs are computed with the regularised Ketju integrator. MBHs with masses in the range of $10^2 - 10^5 \, \rm M_\odot$ are represented by accreting sink particles without feedback. We find that the presence of NSCs boost nuclear SF (i.e. NSC growth) and MBH accretion by funneling gas to the central few parsecs. Low-mass MBHs grow more rapidly on $\sim 600$ Myr timescales, exceeding their Eddington rates at peak accretion. MBH accretion and nuclear SF is episodic (i.e. leads to multiple stellar generations), coeval and regulated by SN explosions. On 40 - 60 Myr timescales the first SN of each episode terminates MBH accretion and nuclear SF. Without NSCs, low-mass MBHs do not grow and MBH accretion and reduced nuclear SF become irregular and uncorrelated. This study gives the first insights into the possible co-evolution of MBHs and NSCs in low-mass galaxies and highlights the importance of considering dense NSCs in galactic studies of MBH growth.

Strongly lensed supernovae (SNe) are a rare class of transient that can offer tight cosmological constraints that are complementary to methods from other astronomical events. We present a follow-up study of one recently-discovered strongly lensed SN, the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. We measure updated, template-subtracted photometry for SN Zwicky and derive improved time delays and magnifications. This is possible because SNe are transient, fading away after reaching their peak brightness. Specifically, we measure point spread function (PSF) photometry for all four images of SN Zwicky in three Hubble Space Telescope WFC3/UVIS passbands (F475W, F625W, F814W) and one WFC3/IR passband (F160W), with template images taken $\sim 11$ months after the epoch in which the SN images appear. We find consistency to within $2\sigma$ between lens model predicted time delays ($\lesssim1$ day), and measured time delays with HST colors ($\lesssim2$ days), including the uncertainty from chromatic microlensing that may arise from stars in the lensing galaxy. The standardizable nature of SNe Ia allows us to estimate absolute magnifications for the four images, with images A and C being elevated in magnification compared to lens model predictions by about $6\sigma$ and $3\sigma$ respectively, confirming previous work. We show that millilensing or differential dust extinction is unable to explain these discrepancies and find evidence for the existence of microlensing in images A, C, and potentially D, that may contribute to the anomalous magnification.

The 1-point matter density probability distribution function (PDF) captures some of the non-Gaussian information lost in standard 2-point statistics. The matter PDF can be well predicted at mildly non-linear scales using large deviations theory. This work extends those predictions to biased tracers like dark matter halos and the galaxies they host. We model the conditional PDF of tracer counts given matter density using a tracer bias and stochasticity model previously used for photometric data. We find accurate parametrisations for tracer bias with a smoothing scale-independent 2-parameter Gaussian Lagrangian bias model and a quadratic shot noise. We relate those bias and stochasticity parameters to the one for the power spectrum and tracer-matter covariances. We validate the model against the Quijote suite of N-body simulations and find excellent agreement for both halo and galaxy density PDFs and their cosmology dependence. We demonstrate the constraining power of the tracer PDFs and their complementarity to power spectra through a Fisher forecast. We focus on the cosmological parameters $\sigma_8$ and $\Omega_m$ as well as linear bias parameters, finding that the strength of the tracer PDF lies in disentangling tracer bias from cosmology. Our results show promise for applications to spectroscopic clustering data when augmented with a redshift space distortion model

Studying the orbital motion of stars around Sagittarius A* in the Galactic Center provides a unique opportunity to probe the gravitational potential near the supermassive black hole at the heart of our Galaxy. Interferometric data obtained with the GRAVITY instrument at the Very Large Telescope Interferometer (VLTI) since 2016 has allowed us to achieve unprecedented precision in tracking the orbits of these stars. GRAVITY data have been key to detecting the in-plane, prograde Schwarzschild precession of the orbit of the star S2, as predicted by General Relativity. By combining astrometric and spectroscopic data from multiple stars, including S2, S29, S38, and S55 - for which we have data around their time of pericenter passage with GRAVITY - we can now strengthen the significance of this detection to an approximately $10 \sigma$ confidence level. The prograde precession of S2's orbit provides valuable insights into the potential presence of an extended mass distribution surrounding Sagittarius A*, which could consist of a dynamically relaxed stellar cusp comprised of old stars and stellar remnants, along with a possible dark matter spike. Our analysis, based on two plausible density profiles - a power-law and a Plummer profile - constrains the enclosed mass within the orbit of S2 to be consistent with zero, establishing an upper limit of approximately $1200 \, M_\odot$ with a $1 \sigma$ confidence level. This significantly improves our constraints on the mass distribution in the Galactic Center. Our upper limit is very close to the expected value from numerical simulations for a stellar cusp in the Galactic Center, leaving little room for a significant enhancement of dark matter density near Sagittarius A*.

We perform canonical quantization of General Relativity, as an effective quantum field theory below the Planck scale, within the BRST-invariant framework. We show that the promotion of constraints to dynamical equations of motion for auxiliary fields leads to the healthy Hamiltonian flow. In particular, we show that the classical properties of Einstein's gravity, such as vanishing Hamiltonian modulo boundary contribution, is realized merely as an expectation value in appropriate physical states. Most importantly, the physicality is shown not to entail trivial time-evolution for correlation functions. In the present approach we quantize the theory once and for all around the Minkowaski vacuum and treat other would-be classical backgrounds as BRST-invariant coherent states. This is especially important for cosmological spacetimes as it uncovers features that are not visible in ordinary semi-classical treatment. The Poincaré invariance of the vacuum, essential for our quantization, provides strong motivation for spontaneously-broken supersymmetry.

We describe a software package, TomOpt, developed to optimise the geometrical layout and specifications of detectors designed for tomography by scattering of cosmic-ray muons. The software exploits differentiable programming for the modeling of muon interactions with detectors and scanned volumes, the inference of volume properties, and the optimisation cycle performing the loss minimisation. In doing so, we provide the first demonstration of end-to-end-differentiable and inference-aware optimisation of particle physics instruments. We study the performance of the software on a relevant benchmark scenario and discuss its potential applications. Our code is available on Github (Strong et al 2024 available at: github.com/GilesStrong/tomopt).

The shape of ²⁸⁶No is investigated with the relativistic density functional theory on a three-dimensional lattice space without any symmetry restriction in a microscopic and self-consistent way. It is found that the ground state of ²⁸⁶No has a pure non-axial octupole shape and coexists with a tetrahedral isomeric state. The energy difference between the two states is only 0.12 MeV, and they are separated by a potential barrier of about 0.5 MeV. The occurrence of the octupole correlations is analyzed with the evolution of the single-particle levels near the Fermi surface driven by the octupole deformations.

We present a machine learning search for high-redshift ($5.0 < z < 6.5$) quasars using the combined photometric data from the DESI Imaging Legacy Surveys and the WISE survey. We explore the imputation of missing values for high-redshift quasars, discuss the feature selections, compare different machine learning algorithms, and investigate the selections of class ensemble for the training sample, then we find that the random forest model is very effective in separating the high-redshift quasars from various contaminators. The 11-class random forest model can achieve a precision of $96.43\%$ and a recall of $91.53\%$ for high-redshift quasars for the test set. We demonstrate that the completeness of the high-redshift quasars can reach as high as $82.20\%$. The final catalog consists of 216,949 high-redshift quasar candidates with 476 high probable ones in the entire Legacy Surveys DR9 footprint, and we make the catalog publicly available. Using MUSE and DESI-EDR public spectra, we find that 14 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for MUSE, and 20 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for DESI-EDR. Additionally, we estimate photometric redshift for the high-redshift quasar candidates using random forest regression model with a high precision.

In this thesis, I present my work developing two instruments for studying cosmic and solar particles and the secondary radiation created by them. The RadMap Telescope is a compact radiation monitor for characterizing the environment inside the International Space Station. The Lunar Cosmic-Ray and Neutron Spectrometer is a versatile instrument designed to detect areas of increased sub-surface hydrogen abundance via neutron spectroscopy, which we use to search for water-ice deposits in the Moon's polar regions.

Small-scale winds driven from accretion discs surrounding active galactic nuclei (AGN) are expected to launch kpc-scale outflows into their host galaxies. However, the ways in which the structure of the interstellar medium (ISM) affects the multiphase content and impact of the outflow remain uncertain. We present a series of numerical experiments featuring a realistic small-scale AGN wind with velocity $5\times 10^3 \!-\! 10^4\rm {\ km\ s^{-1}}$ interacting with an isolated galaxy disc with a manually controlled clumpy ISM, followed at sub-pc resolution. Our simulations are performed with AREPO and probe a wide range of AGN luminosities ($L_{\rm {AGN}}{=} 10^{43-47}\rm {\ erg\ s^{-1}}$) and ISM substructures. In homogeneous discs, the AGN wind sweeps up an outflowing, cooling shell, where the emerging cold phase dominates the mass and kinetic energy budgets, reaching a momentum flux $\dot{p} \approx 7\ L/c$. However, when the ISM is clumpy, outflow properties are profoundly different. They contain small, long-lived ($\gtrsim 5\ \rm {Myr}$), cold ($T{\lesssim }10^{4.5}{\rm {\ K}}$) cloudlets entrained in the faster, hot outflow phase, which are only present in the outflow if radiative cooling is included in the simulation. While the cold phase dominates the mass of the outflow, most of the kinetic luminosity is now carried by a tenuous, hot phase with $T \gtrsim 10^7 \, \rm K$. While the hot phases reach momentum fluxes $\dot{p} \approx (1 - 5)\ L/c$, energy-driven bubbles couple to the cold phase inefficiently, producing modest momentum fluxes $\dot{p} \lesssim L/c$ in the fast-outflowing cold gas. These low momentum fluxes could lead to the outflows being misclassified as momentum-driven using common observational diagnostics. We also show predictions for scaling relations between outflow properties and AGN luminosity and discuss the challenges in constraining outflow driving mechanisms and kinetic coupling efficiencies using observed quantities.

We study the generation of gravitational waves (GWs) during a first-order cosmological phase transition (PT) using the recently introduced Higgsless approach to numerically evaluate the fluid motion induced by the PT. We present for the first time spectra from strong first-order PTs ($\alpha = 0.5$), alongside weak ($\alpha = 0.0046$) and intermediate ($\alpha = 0.05$) transitions previously considered in the literature. We test the regime of applicability of the stationary source assumption, characteristic of the sound-shell model, and show that it agrees with our numerical results when the kinetic energy, sourcing GWs, does not decay with time. However, we find in general that for intermediate and strong PTs, the kinetic energy in our simulations decays following a power law in time, and provide a theoretical framework that extends the stationary assumption to one that allows to include the time evolution of the source. This decay of the kinetic energy, potentially determined by non-linear dynamics and hence, related to the production of vorticity, modifies the usually assumed linear growth with the source duration to an integral over time of the kinetic energy fraction, effectively reducing the growth rate. We validate the novel theoretical model with the results of our simulations covering a broad range of wall velocities. We provide templates for the GW amplitude and spectral shape for a broad range of PT parameters.

The Euclid mission, designed to map the geometry of the dark Universe, presents an unprecedented opportunity for advancing our understanding of the cosmos through its photometric galaxy cluster survey. This paper focuses on enhancing the precision of halo bias (HB) predictions, which is crucial for deriving cosmological constraints from the clustering of galaxy clusters. Our study is based on the peak-background split (PBS) model linked to the halo mass function (HMF); it extends with a parametric correction to precisely align with results from an extended set of $N$-body simulations carried out with the OpenGADGET3 code. Employing simulations with fixed and paired initial conditions, we meticulously analyze the matter-halo cross-spectrum and model its covariance using a large number of mock catalogs generated with Lagrangian Perturbation Theory simulations with the PINOCCHIO code. This ensures a comprehensive understanding of the uncertainties in our HB calibration. Our findings indicate that the calibrated HB model is remarkably resilient against changes in cosmological parameters including those involving massive neutrinos. The robustness and adaptability of our calibrated HB model provide an important contribution to the cosmological exploitation of the cluster surveys to be provided by the Euclid mission. This study highlights the necessity of continuously refining the calibration of cosmological tools like the HB to match the advancing quality of observational data. As we project the impact of our model on cosmological constraints, we find that, given the sensitivity of the Euclid survey, a miscalibration of the HB could introduce biases in cluster cosmology analyses. Our work fills this critical gap, ensuring the HB calibration matches the expected precision of the Euclid survey. The implementation of our model is publicly available in github.com/TiagoBsCastro/CCToolkit.

We investigate the redshift evolution of the concentration-mass relationship of dark matter haloes in state-of-the-art cosmological hydrodynamic simulations and their dark-matter-only counterparts. By combining the IllustrisTNG suite and the novel MillenniumTNG simulation, our analysis encompasses a wide range of box size ($50 - 740 \: \rm cMpc$) and mass resolution ($8.5 \times 10^4 - 3.1 \times 10^7 \: \rm M_{\odot}$ per baryonic mass element). This enables us to study the impact of baryons on the concentration-mass relationship in the redshift interval $0<z<7$ over an unprecedented halo mass range, extending from dwarf galaxies to superclusters ($\sim 10^{9.5}-10^{15.5} \, \rm M_{\odot}$). We find that the presence of baryons increases the steepness of the concentration-mass relationship at higher redshift, and demonstrate that this is driven by adiabatic contraction of the profile, due to gas accretion at early times, which promotes star formation in the inner regions of haloes. At lower redshift, when the effects of feedback start to become important, baryons decrease the concentration of haloes below the mass scale $\sim 10^{11.5} \, \rm M_{\odot}$. Through a rigorous information criterion test, we show that broken power-law models accurately represent the redshift evolution of the concentration-mass relationship, and of the relative difference in the total mass of haloes induced by the presence of baryons. We provide the best-fit parameters of our empirical formulae, enabling their application to models that mimic baryonic effects in dark-matter-only simulations over six decades in halo mass in the redshift range $0<z<7$.

We introduce the parity-odd power (POP) spectra, a novel set of observables for probing parity violation in cosmological N-point statistics. POP spectra are derived from composite fields obtained by applying non-linear transformations, involving also gradients, curls, and filtering functions, to a scalar field. This compresses the parity-odd trispectrum into a power spectrum. These new statistics offer several advantages: they are computationally fast to construct, estimating their covariance is less demanding compared to estimating that of the full parity-odd trispectrum, and they are simple to model theoretically. We measure the POP spectra on simulations of a scalar field with a specific parity-odd trispectrum shape. We compare these measurements to semi-analytic theoretical calculations and find agreement. We also explore extensions and generalizations of these parity-odd observables.

Understanding star formation in galaxies requires resolving the physical scale on which star formation often occurs: the scale of star clusters. We present a multiwavelength, eight-parsec resolution study of star formation in the circumnuclear star cluster and molecular gas rings of the early-type spiral NGC 1386. The cluster ring formed simultaneously ~4 Myr ago. The clusters have similar properties in terms of mass and star formation rate, resembling those of H II regions in the Milky Way disc. The molecular CO gas resolves into long filaments, which define a secondary ring detached from the cluster ring. Most clusters are in CO voids. Their separation with respect to the CO filaments is reminiscent of that seen in galaxy spiral arms. By analogy, we propose that a density wave through the disc of this galaxy may have produced this gap in the central kpc. The CO filaments fragment into strings of dense, unresolved clouds with no evidence of a stellar counterpart. These clouds may be the sites of a future population of clusters in the ring. The free-fall time of these clouds, ~10 Myr, is close to the orbital time of the CO ring. This coincidence could lead to a synchronous bursting ring, as is the case for the current ring. The inward spiralling morphology of the CO filaments and co-spatiality with equivalent kpc-scale dust filaments are suggestive of their role as matter carriers from the galaxy outskirts to feed the molecular ring and a moderately active nucleus.

The accuracy of reaction theories used to extract properties of exotic nuclei from scattering experiments is often unknown or not quantified, but of utmost importance when, e.g., constraining the equation of state of asymmetric nuclear matter from observables as the neutron-skin thickness. In order to test the Glauber multiple-scattering model, the total interaction cross section of Image 1 on carbon targets was measured at initial beam energies of 400, 550, 650, 800, and 1000 MeV/nucleon. The measurements were performed during the first experiment of the newly constructed R³B (Reaction with Relativistic Radioactive Beams) experiment after the start of FAIR Phase-0 at the GSI/FAIR facility with beam energies of 400, 550, 650, 800, and 1000 MeV/nucleon. The combination of the large-acceptance dipole magnet GLAD and a newly designed and highly efficient Time-of-Flight detector enabled a precise transmission measurement with several target thicknesses for each initial beam energy with an experimental uncertainty of ±0.4%. A comparison with the Glauber model revealed a discrepancy of around 3.1% at higher beam energies, which will serve as a crucial baseline for the model-dependent uncertainty in future fragmentation experiments.

Multiply lensed images of a same source experience a relative time delay in the arrival of photons due to the path length difference and the different gravitational potentials the photons travel through. This effect can be used to measure absolute distances and the Hubble constant (H₀) and is known as time-delay cosmography. The method is independent of the local distance ladder and early-universe physics and provides a precise and competitive measurement of H₀. With upcoming observatories, time-delay cosmography can provide a 1% precision measurement of H₀ and can decisively shed light on the current reported `Hubble tension'. This manuscript details the general methodology developed over the past decades in time-delay cosmography, discusses recent advances and results, and, foremost, provides a foundation and outlook for the next decade in providing accurate and ever more precise measurements with increased sample size and improved observational techniques.

We identify axion miniclusters collapsing in the radiation-dominated era and follow them to redshift z=99 with N-body simulations. We find that the majority of the densest miniclusters end up in the center of larger minicluster halos at late times. Soon after their formation, the miniclusters exhibit Navarro-Frenk-White (NFW) profiles but they subsequently develop a steeper inner slope approaching ρ∼r-2 on small scales. Using the so far most highly resolved axion structure formation simulation with 20483 particles we examine the structure of previously studied minicluster halos. While the density profiles of their subhalos are NFW-like we confirm that a modified NFW profile with a steeper inner slope provides a better description for minicluster halos with masses above ∼10-12M⊙. We show that miniclusters with a higher central density might be in contrast to pure NFW halos dense enough to induce gravitational microlensing. Likewise, more compact minicluster halos will have immediate implications for direct and indirect axion detection.

Neutron stars (NSs) are observed with high space velocities and elliptical orbits in binaries. The magnitude of these effects points to natal kicks that originate from asymmetries during the supernova (SN) explosions. Using a growing set of long-time 3D SN simulations with the PROMETHEUS-VERTEX code, we explore the interplay of NS kicks that are induced by asymmetric neutrino emission and by asymmetric mass ejection. Anisotropic neutrino emission can arise from a large-amplitude dipolar convection asymmetry inside the proto-NS (PNS) termed LESA (Lepton-number Emission Self-sustained Asymmetry) and from aspherical accretion downflows around the PNS, which can lead to anisotropic neutrino emission (absorption/scattering) with a neutrino-induced NS kick roughly opposite to (aligned with) the kick by asymmetric mass ejection. In massive progenitors, hydrodynamic kicks can reach up to more than 1300 km s^‑1, whereas our calculated neutrino kicks reach (55–140) km s^‑1 (estimated upper bounds of (170–265) km s^‑1) and only ∼(10–50) km s^‑1, if LESA is the main cause of asymmetric neutrino emission. Therefore, hydrodynamic NS kicks dominate in explosions of high-mass progenitors, whereas LESA-induced neutrino kicks dominate for NSs born in low-energy SNe of the lowest-mass progenitors, when these explode nearly spherically. Our models suggest that the Crab pulsar with its velocity of ∼160 km s^‑1, if born in the low-energy explosion of a low-mass, single-star progenitor, should have received a hydrodynamic kick in a considerably asymmetric explosion. Black holes, if formed by the collapse of short-lived PNSs and solely kicked by anisotropic neutrino emission, obtain velocities of only some km s^‑1.

LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2 μK-arcmin and a resolution of 0.5° at 100 GHz. Its primary goal is to measure the tensor-toscalar ratio r with an uncertainty δr = 0.001, including systematic errors and margin. If r ≥ 0.01, LiteBIRD expects to achieve a > 5σ detection in the l = 2-10 and l = 11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects.

We utilize the Magneticum suite of state-of-the-art hydrodynamical, as well as dark-matter-only simulations to investigate the effects of baryonic physics on cosmic voids in the highest-resolution study of its kind. This includes the size, shape and inner density distributions of voids, as well as their radial density and velocity profiles traced by (sub-) halos, baryonic and cold dark matter particles. Our results reveal observationally insignificant effects that slightly increase with the inner densities of voids and are exclusively relevant on scales of only a few Mpc. Most notably, we identify deviations in the distributions of baryons and cold dark matter around halo-defined voids, relevant for weak lensing studies. In contrast, we find that voids identified in cold dark matter, as well as in halos of fixed tracer density exhibit nearly indistinguishable distributions and profiles between hydrodynamical and dark-matter-only simulations, consolidating the universality and robustness of the latter for comparisons of void statistics with observations in upcoming surveys. This corroborates that voids are the components of the cosmic web that are least affected by baryonic physics, further enhancing their use as cosmological probes.

Collective nuclear excitations, like giant resonances, are sensitive to nuclear deformation, as evidenced by alterations in their excitation energies and transition strength distributions. A common theoretical framework to study these collective modes, the random-phase approximation (RPA), has to deal with large dimensions spanned by all possible particle-hole configurations satisfying certain symmetries. It is the aim of this work to establish a new theoretical framework to study the impact of deformation on spin-isospin excitations, that is able to provide fast and reliable solutions of the RPA equations. The nuclear ground state is determined with the axially deformed relativistic Hartree-Bogoliubov (RHB) model based on relativistic point-coupling energy density functionals (EDFs). To study the excitations in the charge-exchange channel, an axially deformed proton-neutron relativistic quasiparticle RPA (pnRQRPA) is developed in the linear response approach. After benchmarking the axially deformed pnRQRPA in the spherical limit, a study of spin-isospin excitations including Fermi, Gamow-Teller (GT), and spin-dipole (SD) is performed for selected pf-shell nuclei. For GT transitions, it is demonstrated that deformation leads to a considerable fragmentation of the strength function. A mechanism inducing the fragmentation is studied by decomposing the total strength to different projections of total angular momentum K and constraining the nuclear shape to either spherical, prolate, or oblate. A similar fragmentation is also observed for SD transitions, although somewhat moderated by the complex structure of these transitions, while, as expected, the Fermi strength is almost shape independent. The axially deformed pnRQRPA introduced in this work open perspectives for the future studies of deformation effects on astrophysically relevant weak interaction processes, in particular beta decay and electron capture.

A comparison of theoretical and experimental values of the scalar spin-spin interaction ($J$-coupling) in tritium deuteride molecules yield constraints for nucleon-nucleon exotic interactions of the dimensionless coupling strengths $g_Vg_V$, $g_Ag_A$ and $g_pg_p$, corresponding to the exchange of an vector, axial-vector, and pseudoscalar (axionlike) boson. The couplings between proton ($p$) and nucleon ($N$), denoted by $g_V^p g_V^N$, $g_p^p g_p^N$ are constrained to be less than $1.4 \times 10^{-6}$ and $2.7\times 10^{-6}$, respectively, for boson masses around 5 keV. The coupling constant $g_A^p g_A^N$ is constrained to be less than $1.0 \times 10^{-18}$ for boson masses $\leq 100$ eV. It is noteworthy that this study represents the first instance in which constraints on $g_V g_V$ have been established through the analysis of the potential term $V_2 + V_3$ for both tritium deuteride and hydrogen deuteride molecules.

A recent evaluation of three-loop nonplanar Feynman integrals contributing to Higgs plus jet production has established their dependence on two novel symbol letters. We show that the resulting alphabet is described by a $G_2$ cluster algebra, enlarging the $C_2$ cluster algebra found to cover all previously known integrals relevant for this process. The cluster algebra connection we find reveals new adjacency relations, which significantly reduce the function space dimension of the non-planar triple ladder integral. These adjacencies may be understood in part by embedding $G_2$ inside higher-rank cluster algebras.

Novel interactions beyond the four known fundamental forces in nature (electromagnetic, gravitational, strong and weak interactions), may arise due to "new physics" beyond the standard model, manifesting as a "fifth force". This review is focused on spin-dependent fifth forces mediated by exotic bosons such as spin-0 axions and axionlike particles and spin-1 Z' bosons, dark photons, or paraphotons. Many of these exotic bosons are candidates to explain the nature of dark matter and dark energy, and their interactions may violate fundamental symmetries. Spin-dependent interactions between fermions mediated by the exchange of exotic bosons have been investigated in a variety of experiments, particularly at the low-energy frontier. Experimental methods and tools used to search for exotic spin-dependent interactions, such as atomic comagnetometers, torsion balances, nitrogen-vacancy spin sensors, and precision atomic and molecular spectroscopy, are described. A complete set of interaction potentials, derived based on quantum field theory with minimal assumptions and characterized in terms of reduced coupling constants, are presented. A comprehensive summary of existing experimental and observational constraints on exotic spin-dependent interactions is given, illustrating the current research landscape and promising directions of further research.

Large angular scale surveys in the absence of atmosphere are essential for measuring the primordial $B$-mode power spectrum of the Cosmic Microwave Background (CMB). Since this proposed measurement is about three to four orders of magnitude fainter than the temperature anisotropies of the CMB, in-flight calibration of the instruments and active suppression of systematic effects are crucial. We investigate the effect of changing the parameters of the scanning strategy on the in-flight calibration effectiveness, the suppression of the systematic effects themselves, and the ability to distinguish systematic effects by null-tests. Next-generation missions such as LiteBIRD, modulated by a Half-Wave Plate (HWP), will be able to observe polarisation using a single detector, eliminating the need to combine several detectors to measure polarisation, as done in many previous experiments and hence avoiding the consequent systematic effects. While the HWP is expected to suppress many systematic effects, some of them will remain. We use an analytical approach to comprehensively address the mitigation of these systematic effects and identify the characteristics of scanning strategies that are the most effective for implementing a variety of calibration strategies in the multi-dimensional space of common spacecraft scan parameters. We also present Falcons, a fast spacecraft scanning simulator that we developed to investigate this scanning parameter space.

We compute the electron self-energy in Quantum Electrodynamics to three loops in terms of iterated integrals over kernels of elliptic type. We make use of the differential equations method, augmented by an $\epsilon$-factorized basis, which allows us to gain full control over the differential forms appearing in the iterated integrals to all orders in the dimensional regulator. We obtain compact analytic expressions, for which we provide generalized series expansion representations that allow us to evaluate the result numerically for all values of the electron momentum squared. As a by product, we also obtain $\epsilon$-resummed results for the self-energy in the on-shell limit $p^2 = m^2$, which we use to recompute the known three-loop renormalization constants in the on-shell scheme.

Nested sampling (NS) is a stochastic method for computing the log-evidence of a Bayesian problem. It relies on stochastic estimates of prior volumes enclosed by likelihood contours, which limits the accuracy of the log-evidence calculation. We propose to transform the prior volume estimation into a Bayesian inference problem, which allows us to incorporate a smoothness assumption for likelihood-prior volume relations. As a result, we aim to increase the accuracy of the volume estimates and thus improve the overall log-evidence calculation using NS. The method presented works as a post-processing step for NS and provides posterior samples of the likelihood-prior-volume relation, from which the log-evidence can be calculated. We demonstrate an implementation of the algorithm and compare its results with plain NS on two synthetic datasets for which the underlying evidence is known. We find a significant improvement in accuracy for runs with less than one hundred active samples in NS, but are prone to numerical problems beyond this point.

We address how the interplay between the finite availability and carrying capacity of particles at different parts of a spatially extended system can control the steady state currents and density profiles in the one-dimensional current-carrying lanes connecting the different parts of the system. To study this, we set up a minimal model consisting of two particle reservoirs of the same finite carrying capacity connected by two equally sized anti-parallel asymmetric exclusion processes (TASEP). We focus on the steady-state currents and particle density profiles in the two TASEP lanes. The ensuing phases and the phase diagrams, which can be remarkably complex, are parametrized by the model parameters defining particle exchange between the TASEP lanes and the reservoirs and the filling fraction of the particles that determine the total resources available. These parameters may be tuned to make the densities of the two TASEP lanes globally uniform or piece-wise continuous in the form of a combination of a single localized domain wall and a spatially constant density or a pair of delocalized domain walls. Our model reveals that the two reservoirs can be preferentially populated or depopulated in the steady states.

Active matter systems evade the constraints of thermal equilibrium, leading to the emergence of intriguing collective behavior. A paradigmatic example is given by motor-filament mixtures, where the motion of motor proteins drives alignment and sliding interactions between filaments and their self-organization into macroscopic structures. After defining a microscopic model for these systems, we derive continuum equations, exhibiting the formation of active supramolecular assemblies such as micelles, bilayers, and foams. The transition between these structures is driven by a branching instability, which destabilizes the orientational order within the micelles, leading to the growth of bilayers at high microtubule densities. Additionally, we identify a fingering instability, modulating the shape of the micelle interface at high motor densities. We study the role of various mechanisms in these two instabilities, such as contractility, active splay, and anchoring, allowing for generalization beyond the system considered here.

How can Transformers model and learn enumerative geometry? What is a robust procedure for using Transformers in abductive knowledge discovery within a mathematician-machine collaboration? In this work, we introduce a new paradigm in computational enumerative geometry in analyzing the $\psi$-class intersection numbers on the moduli space of curves. By formulating the enumerative problem as a continuous optimization task, we develop a Transformer-based model for computing $\psi$-class intersection numbers based on the underlying quantum Airy structure. For a finite range of genera, our model is capable of regressing intersection numbers that span an extremely wide range of values, from $10^{-45}$ to $10^{45}$. To provide a proper inductive bias for capturing the recursive behavior of intersection numbers, we propose a new activation function, Dynamic Range Activator (DRA). Moreover, given the severe heteroscedasticity of $\psi$-class intersections and the required precision, we quantify the uncertainty of the predictions using Conformal Prediction with a dynamic sliding window that is aware of the number of marked points. Next, we go beyond merely computing intersection numbers and explore the enumerative "world-model" of the Transformers. Through a series of causal inference and correlational interpretability analyses, we demonstrate that Transformers are actually modeling Virasoro constraints in a purely data-driven manner. Additionally, we provide evidence for the comprehension of several values appearing in the large genus asymptotic of $\psi$-class intersection numbers through abductive hypothesis testing.

We consider the dimensional reduction of N=(2,0) conformal supergravity in six dimensions on a two-torus to N=4 conformal supergravity in four dimensions. At the level of kinematics, the six-dimensional Weyl multiplet is shown to reduce to a mixture of the N=4 Weyl and vector multiplets, which can be reinterpreted as a new off-shell multiplet of N=4 conformal supergravity. Similar multiplets have been constructed in other settings and are referred to as dilaton Weyl multiplets. We derive it here for the first time in a maximally supersymmetric context in four dimensions. Furthermore, we present the non-linear relations between all the six- and four-dimensional bosonic and fermionic fields, that are obtained by comparing the off-shell supersymmetry transformation rules.

Modern hydrodynamic simulations of core-collapse supernovae and neutron-star mergers require knowledge not only of the equilibrium properties of strongly interacting matter, but also of the system's response to perturbations, encoded in various transport coefficients. Using perturbative and holographic tools, we derive here an improved weak-coupling and a new strong-coupling result for the most important transport coefficient of unpaired quark matter, its bulk viscosity. These results are combined in a simple analytic pocket formula for the quantity that is rooted in perturbative quantum chromodynamics at high densities but takes into account nonperturbative holographic input at neutron-star densities, where the system is strongly coupled. This expression can be used in the modeling of unpaired quark matter at astrophysically relevant temperatures and densities.

In this work we present the study of $p\Lambda$ and $pp\Lambda$ scattering processes using femtoscopic correlation functions. This observable has been recently used to access the low-energy interaction of hadrons emitted in the final state of high-energy collisions, delivering unprecedented precision information of the interaction among strange hadrons. The formalism for particle pairs is well established and it relates the measured correlation functions with the scattering wave function and the emission source. In the present work we analyze the $NN\Lambda$ scattering in free space and relate the corresponding wave function to the $pp\Lambda$ correlation measurement performed by the ALICE collaboration. The three-body problem is solved using the hyperspherical adiabatic basis. Regarding the $p\Lambda$ and $pp\Lambda$ interactions, different models are used and their impact on the correlation function is studied. The three body force considered in this work is anchored to describe the binding energy of the hypertriton and to give a good description of the two four-body hypernuclei. As a main result we have observed a huge, low-energy peak in the $pp\Lambda$ correlation function, mainly produced by the $J^\pi=1/2^+$ three-body state. The study of this peak from an experimental as well as a theoretical point of view will provide important constraints to the two- and three-body interactions.

We explore the potential application of quantum computers to the examination of lattice holography, which extends to the strongly coupled bulk theory regime. With adiabatic evolution, we compute the ground state of a spin system on a (2 +1 )-dimensional hyperbolic lattice, and measure the spin-spin correlation function on the boundary. Notably, we observe that with achievable resources for coming quantum devices, the correlation function demonstrates an approximate scale-invariant behavior, aligning with the pivotal theoretical predictions of the anti-de Sitter/conformal field theory correspondence.

We argue that the Standard Model is accompanied by a new pseudo-scalar degree of freedom, $\eta_w$-meson, which cancels the topological susceptibility of the electroweak vacuum and gets its mass from this effect. The prediction is based on the analyticity properties of the Chern-Simons correlator combined with the basic features of gravity. Depending on the quality-level of the $U(1)_{B+L}$-symmetry, $\eta_w$ emerges as a $B+L$ pseudo-Goldstone boson or as a Stückelberg $2$-form of the electroweak gauge redundancy. An intriguing scenario of the first category is the emergence of $\eta_w$ in the form of the phase of a $U(1)_{B+L}$-violating fermion condensate triggered by the instantons, somewhat similarly to $\eta'$-meson in QCD. Regardless of its particular origin, the presence of $\eta_w$-meson in the theory appears to be a matter of consistency.

The projected sensitivity of the effective electron neutrino-mass measurement with the KATRIN experiment is below 0.3 eV (90 % CL) after five years of data acquisition. The sensitivity is affected by the increased rate of the background electrons from KATRIN's main spectrometer. A special shifted-analysing-plane (SAP) configuration was developed to reduce this background by a factor of two. The complex layout of electromagnetic fields in the SAP configuration requires a robust method of estimating these fields. We present in this paper a dedicated calibration measurement of the fields using conversion electrons of gaseous $^\mathrm{83m}$Kr, which enables the neutrino-mass measurements in the SAP configuration.

The amount of turbulent pressure in galaxy clusters is still debated, especially as for the impact of the dynamical state and the hydro-method used for simulations. We study the turbulent pressure fraction in the intra cluster medium of massive galaxy clusters. We aim to understand the impact of the hydrodynamical scheme, analysis method, and dynamical state on the final properties of galaxy clusters from cosmological simulations. We perform non-radiative simulations of a set of zoom-in regions of seven galaxy clusters with Meshless Finite Mass (MFM) and Smoothed Particle Hydrodynamics (SPH). We use three different analysis methods based on: $(i)$ the deviation from hydrostatic equilibrium, $(ii)$ the solenoidal velocity component obtained by a Helmholtz-Hodge decomposition, and $(iii)$ the small-scale velocity obtained through a multi-scale filtering approach. We split the sample of simulated clusters into active and relaxed clusters. Our simulations predict an increased turbulent pressure fraction for active compared to relaxed clusters. This is especially visible for the velocity-based methods. For these, we also find increased turbulence for the MFM simulations compared to SPH, consistent with findings from more idealized simulations. The predicted non-thermal pressure fraction varies between a few percent for relaxed clusters and $\approx13\%$ for active ones within the cluster center and increases towards the outskirts. No clear trend with redshift is visible. Our analysis quantitatively assesses the importance played by the hydrodynamical scheme and the analysis method to determine the non-thermal/turbulent pressure fraction. While our setup is relatively simple (non-radiative runs), our simulations show agreement with previous, more idealized simulations, and make a step further toward the understanding of turbulence.

Context. The Gaia mission has delivered hundreds of thousands of variable star light curves in multiple wavelengths. Recent work demonstrates that these light curves can be used to identify (non-)radial pulsations in OBAF-type stars, despite their irregular cadence and low light curve precision, of the order of a few millimagnitudes. With the considerably more precise TESS photometry, we revisited these candidate pulsators to conclusively ascertain the nature of their variability.
Aims: We seek to re-classify the Gaia light curves with the first two years of TESS photometry for a sample of 58 970 p- and g-mode pulsators, encompassing γ Dor, δ Scuti, slowly pulsating B, and β Cep variables. From the TESS data, we seek to assess the quality of Gaia's classification of non-radial pulsators, which is based on sparse, years-long light curves of millimagnitude precision. We also supply four new catalogues containing the confirmed pulsators, along with their dominant and secondary pulsation frequencies, the number of independent mode frequencies, and a ranking according to their usefulness for future asteroseismic ensemble analysis.
Methods: We first analysed the TESS light curves independent of their Gaia classification by pre-whitening all dominant pulsation modes down to a 1% false alarm probability. Using this, in combination with a feature-based random forest classifier, we identified different variability types across the sample.
Results: We find that the Gaia photometry is exceptionally accurate for detecting the dominant and secondary frequencies, reaching approximately 80% accuracy in frequency for p- and g-mode pulsators. The majority of Gaia classifications are consistent with the classifications from the TESS data, illustrating the power of the low-cadence Gaia photometry for pulsation studies. We find that the sample of g-mode pulsators forms a continuous group of variable stars along the main sequence across B, A, and F spectral types, implying that the mode excitation mechanisms for all these pulsators need to be updated with improved physics. Finally, we provide a rank-ordered table of pulsators according to their asteroseismic potential for follow-up studies, based on the number of sectors they have been observed in, their classification probability, and the number of independent modes found in the TESS light curves from the nominal mission.
Conclusions: Our catalogue offers a major increase in the number of confirmed g-mode pulsators with an identified dominant mode suitable for follow-up TESS ensemble asteroseismology of such stars.

Full Tables 1 and 3 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/688/A93.

The results of intermediate data products can be found at Zenodo https://doi.org/10.5281/zenodo.12578015.

Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap–ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring owing to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid growth and possibly planetesimal production via the streaming instability. We present evidence of dust trapping in the bright ring of the planet-hosting disk Elias 2-24, from the analysis of 1.3 and 3 mm Atacama Large Millimeter/submillimeter Array observations at high spatial resolution (0.″029, 4.0 au). We leverage the high spatial resolution to demonstrate that larger grains are more efficiently trapped and place constraints on the local turbulence (8 × 10^‑4 < α _turb < 0.03) and the gas-to-dust ratio (Σ _g /Σ _d < 30) in the ring. Using a scattering-included marginal probability analysis, we measure a total dust disk mass of <inline-formula> <mml:math overflow="scroll"><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mi>dust</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mrow><mml:mn>13.8</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>0.5</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>0.7</mml:mn></mml:mrow></mml:msubsup><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>4</mml:mn></mml:mrow></mml:msup><mml:mspace width="0.33em"></mml:mspace><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>&#8857;</mml:mo></mml:mrow></mml:msub></mml:math> </inline-formula>. We also show that at the orbital radius of the proposed perturber the gap is cleared of material down to a flux contrast of 10^‑3 of the peak flux in the disk.

The Orion Nebula Cluster (ONC) is the closest site of very young (∼1 Myr) massive star formation The ONC hosts more than 1600 young and X-ray bright stars with masses ranging from ∼0.1–35 M _⊙. The Chandra HETGS Orion Legacy Project observed the ONC with the Chandra High Energy Transmission Grating Spectrometer (HETGS) for 2.1 Ms. We describe the spectral extraction and cleaning processes necessary to separate overlapping spectra. We obtained 36 high-resolution spectra, which include a high-brilliance X-ray spectrum of θ ¹ Ori C with over 100 highly significant X-ray lines. The lines show Doppler broadening between 300 and 400 km s^‑1. Higher spectral diffraction orders allow us to resolve line components of high Z He-like triplets in θ ¹ Ori C with unprecedented spectral resolution. Long-term light curves spanning ∼20 yr show all stars to be highly variable, including the massive stars. Spectral fitting with thermal coronal emission line models reveals that most sources show column densities of up to a few times 10²² cm^‑2 and high coronal temperatures of 10–90 MK. We observe a bifurcation of the high-temperature component where some stars show a high component of 40 MK, while others show above 60 MK, indicating heavy flaring activity. Some lines are resolved with Doppler broadening above our threshold of ∼200 km s^‑1, up to 500 km s^‑1. This data set represents the largest collection of HETGS high-resolution X-ray spectra from young pre-main-sequence stars in a single star-forming region to date.

Context. Planet formation models are necessary to understand the origins of diverse planetary systems. Circumstellar disc substructures have been proposed as preferred locations of planet formation, but a complete formation scenario has not been covered by a single model so far.
Aims: We aim to study the formation of giant planets facilitated by disc substructure and starting with sub-micron-sized dust.
Methods: We connect dust coagulation and drift, planetesimal formation, N-body gravity, pebble accretion, planet migration, planetary gas accretion, and gap opening in one consistent modelling framework.
Results: We find rapid formation of multiple gas giants from the initial disc substructure. The migration trap near the substructure allows for the formation of cold gas giants. A new pressure maximum is created at the outer edge of the planetary gap, which triggers the next generation of planet formation resulting in a compact chain of giant planets. A high planet formation efficiency is achieved, as the first gas giants are effective at preventing dust from drifting further inwards, which preserves material for planet formation.
Conclusions: Sequential planet formation is a promising framework to explain the formation of chains of gas and ice giants.

Numerical simulations indicate that the clustering of dark matter halos is not only dependent on the halo masses but has a secondary dependence on other properties, such as the assembly history of the halo. This phenomenon, known as the halo assembly bias (HAB), has been found mostly on galaxy scales; observational evidence on larger scales is scarce. In this work, we propose a novel method for exploring HAB on cluster scales using large samples of superclusters. Leveraging the largest-ever X-ray galaxy cluster and supercluster samples obtained from the first SRG/eROSITA all-sky survey, we constructed two subsamples of galaxy clusters that consist of supercluster members and isolated clusters, respectively. After correcting for the selection effects on redshift, mass, and survey depth, we computed the excess in the concentration of the intracluster gas of isolated clusters with respect to supercluster members, defined as δc_gas ≡ c_{gas, ISO}/c_{gas, SC} − 1, to investigate the environmental effect on the concentration of clusters, a sign of HAB on cluster scales. We find that the average gas mass concentration of isolated clusters is a few percent higher than that of supercluster members, with a maximum significance of 2.8σ. The result for δc_gas varies with the overdensity ratio, f, in supercluster identification, cluster mass proxies, and mass and redshift ranges but remains positive in almost all the measurements. We measure slightly larger δc_gas when adopting a higher f for supercluster identification. The δc_gas is also higher for low-mass and low-redshift clusters. We performed weak lensing analyses to compare the total mass concentration of the two classes and find a similar trend in total mass concentration as obtained from the gas mass concentration. Our results are consistent with the prediction of HAB on cluster scales, where halos located in denser environments are less concentrated; this trend is stronger for halos with lower masses and at lower redshifts. These phenomena can be explained by the fact that clusters in denser environments, such as superclusters, have experienced more mergers than isolated clusters in their assembling history. This work paves the way to explore HAB with X-ray superclusters and demonstrates that large samples of superclusters with X-ray and weak-lensing data can advance our understanding of the evolution of the large-scale structure.

We present a high-resolution numerical study of the sinking and merging of massive black holes (MBHs) with masses in the range of $10^3 - 10^7 \, \mathrm{M}_\odot$ in multiple minor mergers of low-mass dark matter haloes without and with galaxies ($4\times 10^8 \, \mathrm{M}_\odot \lesssim {M}_{\mathrm{halo}} \lesssim 2\times 10^{10} \, \mathrm{M}_\odot)$. The KETJU simulation code, a combination of the GADGET tree solver with accurate regularized integration, uses unsoftened forces between the star/dark matter components and the MBHs for an accurate treatment of dynamical friction and scattering of dark matter/stars by MBH binaries or multiples. Post-Newtonian corrections up to order 3.5 for MBH interactions allow for coalescence by gravitational wave emission and gravitational recoil kicks. Low-mass MBHs ($\lesssim 10^5 \, \mathrm{M}_\odot$) hardly sink to the centre or merge. Sinking MBHs have various complex evolution paths - binaries, triplets, free-floating MBHs, and dynamically or recoil ejected MBHs. Collisional interactions with dark matter alone can drive MBHs to coalescence. The highest mass MBHs of $\gtrsim 10^6 \, \rm M_\odot$ mostly sink to the centre and trigger the scouring of dark matter and stellar cores. The scouring can transform a centrally baryon-dominated system into a dark-matter-dominated system. Our idealized high-resolution study highlights the difficulty to bring in and keep low-mass MBHs in the centres of low-mass haloes/galaxies - a remaining challenge for merger assisted MBH seed growth mechanisms.

The braking torque that dictates the timing properties of magnetars is closely tied to the large-scale dipolar magnetic field on their surface. The formation of this field has been a topic of ongoing debate. One proposed mechanism, based on macroscopic principles, involves an inverse cascade within the neutron star's crust. However, this phenomenon has not been observed in realistic simulations. In this study, we provide compelling evidence supporting the feasibility of the inverse cascading process in the presence of an initial helical magnetic field within realistic neutron star crusts and discuss its contribution to the amplification of the large-scale magnetic field. Our findings, derived from a systematic investigation that considers various coordinate systems, peak wavenumber positions, crustal thicknesses, magnetic boundary conditions, and magnetic Lundquist numbers, reveal that the specific geometry of the crustal domain - with its extreme aspect ratio - requires an initial peak wavenumber from small-scale structures for the inverse cascade to occur. However, this extreme aspect ratio limits the inverse cascade to magnetic field structures on scales comparable to the neutron star's crust, making the formation of a large-scale dipolar surface field unlikely. Despite this limitation, the inverse cascade can significantly impact the magnetic field evolution in the interior of the crust, potentially explaining the observed characteristics of highly magnetized objects with weak surface dipolar fields, such as low-field magnetars or central compact objects.

Context. Large millimeter surveys of star-forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. The ever-increasing number of these surveys has led to a flourishing of population study, a valuable tool and approach that is spreading in ever more fields. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles.
Aims: This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available.
Methods: We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport, varying the initial conditions of the disk and substructure to find the best match with the observed distributions. Disks both with and without substructure have been examined. We obtained the simulated population distribution for the disk sizes, millimeter fluxes, and spectral indices by post-processing the resulting disk profiles (surface density, maximum grain size, and disk temperature).
Results: We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure; namely, a perturbation that is strong enough to be able to trap particles, that is formed early in the evolution of the disk, and that is within 0.4 Myr. Agreement is reached by relatively high initial disk masses (10^−2.3 M_⋆ ⩽ M_disk ⩽ 10^−0.5 M_⋆) and moderate levels of turbulence (10^−3.5 ⩽ α ⩽ 10^−2.5). Other disk parameters play a weaker role. Only opacities with a high absorption efficiency can reproduce the observed spectral indices.
Conclusions: Disk population synthesis is a precious tool for investigating and constraining the parameters and initial conditions of planet-forming disks. The generally low observed spectral indices call for significant substructure, like that which planets in the mass range of Saturn to a few Jupiters would induce, to already be present before 0.4 Myr. Our results indicate that substructure, which so far has only been assessed in individual disks, is likely ubiquitous and extends to the whole population, and imply that most "smooth" disks hide unresolved substructure.

The collective properties of star clusters are investigated using a simulation of the collision between two dwarf galaxies. The characteristic power law of the cluster mass function, N(M), with a logarithmic slope dN/dM ~ -1, is present from cluster birth and remains throughout the simulation. The maximum mass of a young cluster scales with the star formation rate (SFR). The relative average minimum separation, R(M)= N(M)^{1/p}D_min(M)/D(M_low), for average minimum distance D_min(M) between clusters of mass M, and for lowest mass, M_low, measured in projection (p=2) or three dimensions (p=3), has a negative slope, d log R/d log M ~ -0.2, for all masses and ages. This agrees with observations of R(M) in low-mass galaxies studied previously. Like the slope of N(M), R(M) is apparently a property of cluster birth for dwarf galaxies that does not depend on SFR or time. The negative slope for R(M) implies that massive clusters are more concentrated relative to lower mass clusters throughout the entire mass range. Cluster growth through coalescence is also investigated. The ratio of the kinetic to potential energy of all near-neighbor clusters is generally large, but a tail of low values in the distribution of this ratio suggests that a fraction of the clusters merge, ~8% by number throughout the ~300 Myr of the simulation and up to 60% by mass for young clusters in their first 10 Myr, scaling with the SFR above a certain threshold.

Understanding the explosion mechanism and hydrodynamic evolution of core-collapse supernovae (SNe) is a long-standing quest in astronomy. The asymmetries caused by the explosion are encoded into the line profiles which appear in the nebular phase of the SN evolution - with particularly clean imprints in He star explosions. Here, we carry out nine different supernova simulations of He-core progenitors, exploding them in 3D with parametrically varied neutrino luminosities using the PROMETHEUS-HOTB code, hydrodynamically evolving the models to the homologous phase. We then compute nebular phase spectra with the 3D Non-Local Thermodynamic Equilibrium spectral synthesis code EXTRASS (EXplosive TRAnsient Spectral Simulator). We study how line widths and shifts depend on progenitor mass, explosion energy, and viewing angle. We compare the predicted line profile properties against a large set of Type Ib observations, and discuss the degree to which current neutrino-driven explosions can match observationally inferred asymmetries. With self-consistent 3D modelling - circumventing the difficulties of representing $^{56}$Ni mixing and clumping accurately in 1D models - we find that neither low-mass He cores exploding with high energies nor high-mass cores exploding with low energies contribute to the Type Ib SN population. Models which have line profile widths in agreement with this population give sufficiently large centroid shifts for calcium emission lines. Calcium is more strongly affected by explosion asymmetries connected to the neutron star kicks than oxygen and magnesium. Lastly, we turn to the near-infrared spectra from our models to investigate the potential of using this regime to look for the presence of He in the nebular phase.

We present ultraviolet, optical and near-infrared photometric and optical spectroscopic observations of the luminous, fast blue optical transient (LFBOT), CSS161010:045834-081803 (CSS161010). The transient was found in a low-redshift (z=0.033) dwarf galaxy. The light curves of CSS161010 are characterized by an extremely fast evolution and blue colours. The V-band light curve shows that CSS161010 reaches an absolute peak of M$_{V}^{max}=-20.66\pm0.06$ mag in 3.8 days from the start of the outburst. After maximum, CSS161010 follows a power-law decline $\propto t^{-2.8\pm0.1}$ at all optical bands. These photometric properties are comparable to those of well-observed LFBOTs such as AT 2018cow, AT 2020mrf and AT 2020xnd. However, unlike these objects, the spectra of CSS161010 show a remarkable transformation from a blue and featureless continuum to spectra dominated by very broad, entirely blueshifted hydrogen emission lines of velocities of up to 10% of the speed of light. The persistent blueshifted emission and the lack of any emission at the rest wavelength of CSS161010 are unique features not seen in any transient before CSS161010. The combined observational properties of CSS161010 and its dwarf galaxy host favour the tidal disruption of a star by an intermediate-mass black hole as its origin.

The Hubble tension is one of the most relevant unsolved problems in cosmology today. Strongly gravitationally lensed transient objects, such as strongly lensed supernovae, are an independent and competitive probe that can be used to determine the Hubble constant. In this context, the time delay between different images of lensed supernovae is a key ingredient. We present a method, to retrieve time delays and the amount of differential dust extinction between multiple images of lensed type IIP supernovae through their color curves, which display a kink in the time evolution. With multiple realistic mock color curves based on an observed unlensed supernova from the Carnegie Supernova Project, we demonstrate that we can retrieve the time delay with uncertainties of $\pm$1.0 days for light curves with 2-day cadence and 35% missing data due to weather losses. The differential dust extinction is more susceptible to uncertainties, because it depends on imposing the correct extinction law. Further we also investigate the kink structure in the color curves for different rest-frame wavelength bands, particularly rest-frame UV light curves from SWIFT, finding sufficiently strong kinks for our method to work for typical lensed SN redshifts that would redshift the kink feature to optical wavelengths. With the upcoming Rubin Observatory Legacy Survey of Space and Time, hundreds of strongly lensed supernovae will be detected and our new method for lensed SN IIP is readily applicable to provide delays.

We present the discovery of large radio shells around a massive pair of interacting galaxies and extended diffuse X-ray emission within the shells. The radio data were obtained with the Australian Square Kilometre Array Pathfinder (ASKAP) in two frequency bands centred at 944 MHz and 1.4 GHz, respectively, while the X-ray data are from the XMM-Newton observatory. The host galaxy pair, which consists of the early-type galaxies ESO 184-G042 and LEDA 418116, is part of a loose group at a distance of only 75 Mpc (redshift $z = 0.017$). The observed outer radio shells (diameter ${\sim}$145 kpc) and ridge-like central emission of the system, ASKAP J1914-5433 (Physalis), are likely associated with merger shocks during the formation of the central galaxy (ESO 184-G042) and resemble the new class of odd radio circles (ORCs). This is supported by the brightest X-ray emission found offset from the centre of the Physalis system, instead centred at the less massive galaxy, LEDA 418116. The host galaxy pair is embedded in an irregular envelope of diffuse light, highlighting ongoing interactions. We complement our combined radio and X-ray study with high-resolution simulations of the circumgalactic medium (CGM) around galaxy mergers from the Magneticum project to analyse the evolutionary state of the Physalis system. We argue that ORCs/radio shells could be produced by a combination of energy release from the central active galactic nucleus and subsequent lightening up in radio emission by merger shocks travelling through the CGM of these systems.

Despite recent progress in the spectroscopic characterization of individual exoplanets, the atmospheres of key ultra-hot Jupiters (UHJs) still lack comprehensive investigations. These include WASP-178b, one of the most irradiated UHJs known to date. We observed the dayside emission signal of this planet with CRIRES⁺ in the spectral K band. By applying the cross-correlation technique and a Bayesian retrieval framework to the high-resolution spectra, we identified the emission signature of ¹²CO (S/N = 8.9) and H₂O (S/N = 4.9), and a strong atmospheric thermal inversion. A joint retrieval with space-based secondary eclipse measurements from TESS and CHEOPS allowed us to refine our results on the thermal profile and thus to constrain the atmospheric chemistry, yielding a solar to super-solar metallicity (1.4 ± 1.6 dex) and a solar C/O ratio (0.6 ± 0.2). We infer a significant excess of spectral line broadening and identify a slight Doppler-shift between the ¹²CO and H₂O signals. These findings provide strong evidence for a super-rotating atmospheric flow pattern and suggest the possible existence of chemical inhomogeneities across the planetary dayside hemisphere. In addition, the inclusion of photometric data in our retrieval allows us to account for stellar light reflected by the planetary atmosphere, resulting in an upper limit on the geometric albedo (0.23). The successful characterization of WASP-178b's atmosphere through a joint analysis of CRIRES⁺, TESS, and CHEOPS observations highlights the potential of combined studies with space- and ground-based instruments and represents a promising avenue for advancing our understanding of exoplanet atmospheres.

Context: Photoevaporation is an important process for protoplanetary disc dispersal but there has so far been a lack of consensus from simulations over the mass-loss rates and the most important part of the high-energy spectrum for driving the wind. Aims: We aim to isolate the origins of these discrepancies through carefully-benchmarked hydrodynamic simulations of X-ray photoevaporation with time-dependent thermochemistry calculated on the fly. Methods: We conduct hydrodynamic simulations with pluto where the thermochemistry is calculated using prizmo. We explore the contribution of certain key microphysical processes and the impact of using different spectra used previously in literature studies. Results: We find that additional cooling results from the excitation of O by neutral H, which leads to dramatically reduced mass-loss across the disc compared to previous X-ray photoevaporation models, with an integrated rate of 10^-9 Msun/yr. Such rates would allow for longer-lived discs than previously expected from population synthesis. An alternative spectrum with less soft X-ray produces mass-loss rates around a factor of 2-3 times lower. The chemistry is significantly out of equilibrium, with the survival of H2 into the wind aided by advection. This leads to its role as the dominant coolant at 10s au - thus stabilising a larger radial temperature gradient across the wind - as well as providing a possible wind tracer.

Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC 1333 IRAS 4A, located in the Perseus molecular cloud. Using high-resolution (~ 100 au) ALMA observations, we investigate the [D₂CO]/[HDCO] ratio along the cavity walls of the outflows emanating from IRAS 4A1. Our analysis reveals a consistent decrease in the deuteration ratio (from ~ 60-20% to ~ 10%) with increasing distance from the protostar (from ~ 2000 au to ~ 4000 au). Given the large measured [D₂CO]/[HDCO], both HDCO and D₂CO are likely injected by the shocks along the cavity walls into the gas-phase from the dust mantles, formed in the previous prestellar phase. We propose that the observed [D₂CO]/[HDCO] decrease is due to the density profile of the prestellar core from which NGC 1333 IRAS 4A was born. When considering the chemical processes at the base of formaldehyde deuteration, the IRAS 4A's prestellar precursor had a predominantly flat density profile within 3000 au and a decrease of density beyond this radius.

The energy released by AGN has the potential to heat or remove the gas of the ISM, thus likely impacting the cold molecular gas reservoir of host galaxies at first, with star formation following on longer timescales. Previous works on high-z galaxies have yielded conflicting results, possibly due to selection biases and other systematics. To provide a reliable benchmark for galaxy evolution models at cosmic noon (z=1-3), two surveys were conceived: SUPER and KASHz, both targeting unbiased X-ray-selected AGN at z>1 that span a wide bolometric luminosity range. In this paper, we assess the effects of AGN feedback on the molecular gas content of host galaxies in a statistically robust, uniformly selected, coherently analyzed sample of AGN at z=1-2.6, drawn from the KASHz and SUPER surveys. By using ALMA data in combination with dedicated SED modeling, we retrieve CO and FIR luminosity as well as $M_*$ of SUPER and KASHz AGN. We selected non-active galaxies from PHIBBS, ASPECS and multiple ALMA/NOEMA surveys of sub-mm galaxies. By matching the samples in z, $M_*$ and $L_{FIR}$, we compared the properties of AGN and non-active galaxies within a Bayesian framework. We find that AGN hosts at given $L_{FIR}$ are on average CO depleted compared to non-active galaxies, confirming what was previously found in the SUPER survey. Moreover, the molecular gas fraction distributions of AGN and non-active galaxies are statistically different, with that of of AGN being skewed to lower values. Our results indicate that AGN can indeed reduce the total cold molecular gas reservoir of their host galaxies. Lastly, by comparing our results with predictions from three cosmological simulations (TNG, Eagle and Simba) filtered to match the observed properties, we confirm already known discrepancies and highlight new ones between observations and simulations.[Abridged]

Recent hydrodynamical simulations of the late stages of supernova remnant (SNR) evolution have revealed that as they merge with the ambient medium, SNRs implode, leading to the formation of dense clouds in their center. While being highly chemically enriched by their host SNR, these clouds appear to have similar properties as giant molecular clouds, which are believed to be the main site of star formation. Here, we develop a simple model in order to estimate the efficiency of the star formation that might be triggered by the implosion of SNRs. We separately consider two cases: cyclic star formation, maintained by the episodic driving of feedback from new generations of stars, and a single burst of star formation, triggered by a single explosion. We find that in the cyclic case, star formation is inefficient, with implosion-triggered star formation contributing a few percent of the observed star formation efficiency per freefall timescale. In the single-burst case, higher star formation efficiencies can be obtained. However, while the implosion-triggered process might not contribute much to the overall star formation, due to the high chemical enrichment of the birth clouds, it can explain the formation of a significant fraction of metal-rich stars.

We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early time (δ t < 2 days) spectra show transient, narrow emission lines from shock ionization of confined (r < 10¹⁵ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of H I, He I/II, C IV, and N III/IV/V from the CSM persist on a characteristic timescale (t _IIn) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early time IIn-like features in addition to 35 "comparison" SNe with no evidence of early time IIn-like features, all with ultraviolet observations. The total sample includes 50 unpublished objects with a total of 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both t _IIn and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through the matching of peak multiband absolute magnitudes, rise times, t _IIn, and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: <inline-formula> <mml:math overflow="scroll"><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>IIn</mml:mi></mml:mrow></mml:msub><mml:mo>≈</mml:mo><mml:mn>3.8</mml:mn><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>&#775;</mml:mo></mml:mrow></mml:mover><mml:mrow><mml:mo stretchy="true">/</mml:mo></mml:mrow></mml:math> </inline-formula>(0.01 M _&#8857; yr^‑1)] days.

We present a new method to analytically prove global stability in ghost-ridden dynamical systems. Our proposal encompasses all prior results and consequentially extends them. In particular, we show that stability can follow from a conserved quantity that is unbounded from below, contrary to expectation. Novel examples illustrate all our results. Our findings take root on a careful examination of the literature, here comprehensively reviewed for the first time. This work lays the mathematical basis for ulterior extensions to field theory and quantization, and it constitutes a gateway for inter-disciplinary research in dynamics and integrability.

This study investigates the clustering and bias of Luminous Red Galaxies (LRG) in the BOSS-LOWZ, -CMASS, -COMB, and eBOSS samples, using two types of simulated lightcones: (i) high-fidelity lightcones from UCHUUN-body simulation, employing SHAM technique to assign LRG to (sub)haloes, and (ii) 16 000 covariance lightcones from GLAM-UCHUUN-body simulations, including LRG using HOD data from UCHUU. Our results indicate that UCHUU and GLAM lightcones closely replicate BOSS/eBOSS data, reproducing correlation function and power spectrum across scales from redshifts 0.2 to 1.0, from 2 to $150 \,h^{-1}\,\mathrm{Mpc}$ in configuration space, from 0.005 to $0.7\, h\,\mathrm{Mpc}^{-1}$ in Fourier space, and across different LRG stellar masses. Furthermore, comparing with existing MD-PATCHY and EZMOCK BOSS/eBOSS lightcones based on approximate methods, our GLAM-UCHUU lightcones provide more precise clustering estimates. We identify significant deviations from observations within $20 \,h^{-1}\,\mathrm{Mpc}$ scales in MD-PATCHY and EZMOCK, with our covariance matrices indicating that these methods underestimate errors by between 10 per cent and 60 per cent. Lastly, we explore the impact of cosmology on galaxy clustering. Our findings suggest that, given the current level of uncertainties in BOSS/eBOSS data, distinguishing models with and without massive neutrino effects on large-scale structure (LSS) is challenging. This paper highlights the UCHUU and GLAM-UCHUU simulations' robustness in verifying the accuracy of Planck cosmological parameters, providing a strong foundation for enhancing lightcone construction in future LSS surveys. We also demonstrate that generating thousands of galaxy lightcones is feasible using N-body simulations with adequate mass and force resolution.

More than 30% of white dwarfs exhibit atmospheric metals, which are understood to be from recent or ongoing accretion of circumstellar debris. In cool white dwarfs, surface motions should rapidly homogenise photospheric abundances, and the accreted heavy elements should diffuse inward on a timescale much longer than that for surface mixing. The recent discovery of a metal scar on WD 0816-310 implies its B ≈ 140 kG magnetic field has impeded surface mixing of metals near the visible magnetic pole. Here, we report the discovery of a second magnetic, metal-polluted white dwarf, WD 2138-332, which exhibits periodic variability in longitudinal field, metal line strength, and broadband photometry. All three variable quantities have the same period, and show remarkable correlations: the published light curves have a brightness minimum exactly when the longitudinal field and line strength have a maximum, and a maximum when the longitudinal field and line strength have a minimum. The simplest interpretation of the line strength variability is that there is an enhanced metal concentration around one pole of the magnetic field; however, the variable line-blanketing cannot account for the observed multi-band light curves. More theoretical work is required to understand the efficiency of horizontal mixing of the accreted metal atoms, and the origin of photometric variability. Because both magnetic, metal-polluted white dwarfs that have been monitored to date show that metal line strengths vary in phase with the longitudinal field, we suggest that metal scars around magnetic poles may be a common feature of metal-polluted white dwarfs.

Frequentist parameter inference using profile likelihoods has received increased attention in the cosmology literature recently since it can give important complementary information to Bayesian credible intervals. Here, we give a pedagogical review to frequentist parameter inference in cosmology with particular focus on when the graphical profile likelihood construction gives meaningful constraints, i.e.\ confidence intervals with correct coverage. This construction rests on the assumption of the asymptotic limit of a large data set such as in \textit{Wilks' theorem}. We assess the validity of this assumption in the context of three cosmological models with \textit{Planck} 2018 \texttt{Plik\_lite} data: While our tests for the $\Lambda$CDM model indicate that the profile likelihood method gives correct coverage, $\Lambda$CDM with the sum of neutrino masses as a free parameter appears consistent with a Gaussian near a boundary motivating the use of the boundary-corrected or Feldman-Cousins graphical method; for $w_0$CDM with the equation of state of dark energy, $w_0$, as a free parameter, we find indication of a violation of the assumptions. Finally, we compare frequentist and Bayesian constraints of these models. Our results motivate care when using the graphical profile likelihood method in cosmology.

The $\delta N$ formalism has been the major computational tool to study the superhorizon evolution of the scalar type perturbation sourced by scalar fields. Recently, this formalism was generalized to compute an arbitrary scalar, vector, and tensor type perturbations, including the gravitational waves (GWs), sourced by an arbitrary bosonic fields. In this paper, we explain how to use the generalized $\delta N$ formalism (the g$\delta N$ formalism), considering a model with U(1) gauge fields as a concrete example. Several new findings on this model and prospects on future gravitational wave experiments are also discussed, including the condition for the two linear polarizations of GWs to have different amplitudes. This paper provides a detailed explanation of our previous paper published in Physical Review Letters. We also discuss the Weinberg's adiabatic mode for an anisotropic background, showing a qualitative difference from the one for the FLRW background.

Galaxy chemical enrichment mechanisms have primarily been constrained by [$\alpha$/Fe] and [Fe/H] measurements of individual stars and integrated light from stellar populations. However such measurements are limited at higher redshifts (z>1). Recently, we proposed an analogous diagram of the oxygen-to-argon abundance ratio, log(O/Ar), vs Ar abundance, 12+log(Ar/H), as a new diagnostic window for emission nebulae. In this Letter, using robust line flux measurements including temperature sensitive auroral lines, we present direct determination of O and Ar abundances in nine SFGs from JWST/NIRSPEC spectra at z$\sim$1.3-7.7, and two more with Keck/MOSFIRE spectra at z$\sim$2.2. Utilising their positions on the log(O/Ar) vs 12+log(Ar/H) plane, we present the first inference of galaxy chemical enrichment mechanisms from an ensemble of galaxies. The SFGs at z$\sim$1.3-3.4 are consistent with the solar neighbourhood galactic chemical enrichment models of the Milky Way Galaxy that are driven by core-collapse and Type Ia supernovae. Such enrichment mechanisms thus occur at least out to z$\sim$3.4. However, the highest-redshift SFGs (z$\sim$3.6-7.7) have very low log(O/Ar) values, revealing a different enrichment process at z>3.6. Such low log(O/Ar) values may be caused by a rapid but intermittent star-formation and/or additional sources. The new diagnostic window for SFGs enables us to reveal the unique fingerprints of galaxy chemical enrichment out to cosmic dawn.

Modern cosmological research in large-scale structure has witnessed an increasing number of machine-learning applications. Among them, convolutional neural networks (CNNs) have received substantial attention due to their outstanding performance in image classification, cosmological parameter inference, and various other tasks. However, many models based on CNNs are criticized as "black boxes" due to the difficulties in relating their outputs intuitively and quantitatively to the cosmological fields under investigation. To overcome this challenge, we present the Cosmological Correlator Convolutional Neural Network (C3NN)—a fusion of CNN architecture and cosmological N-point correlation functions (NPCFs). We demonstrate that its output can be expressed explicitly in terms of the analytically tractable NPCFs. Together with other auxiliary algorithms, we can open the "black box" by quantitatively ranking different orders of the interpretable outputs based on their contribution to classification tasks. As a proof of concept, we demonstrate this by applying our framework to a series of binary classification tasks using Gaussian and log-normal random fields and relating its outputs to the NPCFs describing the two fields. Furthermore, we exhibit the model's ability to distinguish different dark energy scenarios (w ₀ = ‑0.95 and ‑1.05) using N-body simulated weak-lensing convergence maps and discuss the physical implications coming from their interpretability. With these tests, we show that C3NN combines advanced aspects of machine learning architectures with the framework of cosmological NPCFs, thereby making it an exciting tool to extract physical insights in a robust and explainable way from observational data.

For decades, the boundary of cosmic filaments has been a subject of debate. In this work, we determine the physically motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at $z=0$ in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from $z=4$ to 1, and then slightly increases. Moreover, we find that the radius of the filament depends on the length of the filament, the distance from the connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experienced rapid contraction before $z=1$, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves.

Context. Recent observations suggest that planet formation starts early, in protostellar disks of ≤10⁵ yr, which are characterized by strong interactions with the environment, such as through accretion streamers and molecular outflows.
Aims: To investigate the impact of such phenomena on the physical and chemical properties of a disk, it is key to understand what chemistry planets inherit from their natal environment.
Methods: In the context of the ALMA large program Fifty AU Study of the chemistry in the disk/envelope system of solar-like protostars (FAUST), we present observations on scales from ∼1500 au to ∼60 au of H₂CO, HDCO, and D₂CO toward the young planet-forming disk IRS 63.
Results: The H₂CO probes the gas in the disk as well as in a large scale streamer (∼1500 au) impacting onto the southeast disk side. We detected for the first time deuterated formaldehyde, HDCO and D₂CO, in a planet-forming disk and HDCO in the streamer that is feeding it. These detections allowed us to estimate the deuterium fractionation of H₂CO in the disk: [HDCO]/[H₂CO] ∼ 0.1 − 0.3 and [D₂CO]/[H₂CO] ∼ 0.1. Interestingly, while HDCO follows the H₂CO distribution in the disk and in the streamer, the distribution of D₂CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the southeast disk side, where the streamer crashes onto the disk. In addition, D₂CO was detected in two spots along the blue- and redshifted outflow. This suggests that (i) in the disk, HDCO formation is dominated by gas-phase reactions in a manner similar to H₂CO, while (ii) D₂CO is mainly formed on the grain mantles during the prestellar phase and/or in the disk itself and is at present released in the gas phase in the shocks driven by the streamer and the outflow.
Conclusions: These findings testify to the key role of streamers in the buildup of the disk concerning both the final mass available for planet formation and its chemical composition.

In the present study we employ three distinct physically motivated speed of sound bounds in order to construct hybrid models where the high density phase is described by the maximally stiff equation of state. In particular, we consider the bounds related to special relativity, relativistic kinetic theory and conformality. The low density hadronic phase is described by a state-of-the-art microscopic relativistic Brueckner-Hartree-Fock theory. This work aims to access the effect of the different speed of sound constraints on the relevant parameter space of the key parameters of first-order phase transitions by utilising recent astronomical data. This involves a systematic analysis that also includes two distinct schemes for the construction of hybrid models, namely the Maxwell and Gibbs methods. Finally, a relevant discussion is conducted on the possible occurrence of a thermodynamic inconsistency that is related to the stability of the high density phase over hadronic matter at large densities.

The mass ranges allowed for primordial black holes (PBHs) to constitute all of dark matter (DM) are broadly constrained. However, these constraints rely on the standard semiclassical approximation which assumes that the evaporation process is self-similar. Quantum effects such as memory burden take the evaporation process out of the semiclassical regime latest by the time the black hole loses half of its mass. What happens beyond this time is currently not known. However, theoretical evidence based on prototype models indicates that the evaporation slows down, thereby extending the lifetime of a black hole. This modifies the mass ranges constrained, in particular, by big bang nucleosynthesis (BBN) and cosmic microwave background spectral distortions. We show that previous constraints are largely relaxed when the PBH lifetime is extended, making it possible for PBHs to constitute all of DM in previously excluded mass ranges. In particular, this is the case for PBHs lighter than 10⁹ g that enter the memory burden stage before BBN and are still present today as DM.

In recent years the amount of publicly available astronomical data has increased exponentially, with a remarkable example being large scale multiepoch photometric surveys. This wealth of data poses challenges to the classical methodologies commonly employed in the study of variable objects. As a response, deep learning techniques are increasingly being explored to effectively classify, analyze, and interpret these large datasets. In this paper we use two-dimensional histograms to represent Optical Gravitational Lensing Experiment (OGLE) phasefolded light curves as images. We use a Convolutional Neural Network (CNN) to classify variable objects within eight different categories (from now on labels): Classical Cepheid (CEP), RR Lyrae (RR), Long Period Variable (LPV), Miras (M), Ellipsoidal Binary (ELL), Delta Scuti (DST), Eclipsing Binary (E), and spurious class with Incorrect Periods (Rndm). We set up different training sets to train the same CNN architecture in order to characterize the impact of the training. The training sets were built from the same source of labels but different filters and balancing techniques were applied. Namely: Undersampling (U), Data Augmentation (DA), and Batch Balancing (BB). The best performance was achieved with the BB approach and a training sample size of $\sim$370000 stars. Regarding computational performance, the image representation production rate is of $\sim$76 images per core per second, and the time to predict is $\sim$ 60$\, \mu\text{s}$ per star. The accuracy of the classification improves from $\sim$ 92%, when based only on the CNN, to $\sim$ 98% when the results of the CNN are combined with the period and amplitude features in a two step approach. This methodology achieves comparable results with previous studies but with two main advantages: the identification of miscalculated periods and the improvement in computational time cost.

The vector meson-baryon interaction in a coupled channel scheme is revisited within the correlation function framework. As illustrative cases to reveal the important role played by the coupled channels, we focus on the $\phi$p and $\rho^0$p systems given their complex dynamics and the presence of quasi-bound states or resonances in the vicinity of their thresholds. We show that the $\phi$p femtoscopic data provide novel information about a $N^*$ state present in the experimental region and anticipate the relevance of a future $\rho^0$p correlation function measurement in order to pin down the $S=0, Q=+1$ vector meson-baryon interaction as well as to disclose the characterizing features of the $N^*(1700)$ state.

We have applied an operator-overloading forward-mode algorithmic differentiation tool to the Monte-Carlo particle simulation toolkit Geant4. Our differentiated version of Geant4 allows computing mean pathwise derivatives of user-defined outputs of Geant4 applications with respect to user-defined inputs. This constitutes a major step towards enabling gradient-based optimization techniques in high-energy physics, as well as other application domains of Geant4. This is a preliminary report on the technical aspects of applying operator-overloading AD to Geant4, as well as a first analysis of some results obtained by our differentiated Geant4 prototype. We plan to follow up with a more refined analysis.

Context. The vast majority of close binaries containing a compact object, including the progenitors of supernovae Ia and at least a substantial fraction of all accreting black holes in the Galaxy, form through common-envelope (CE) evolution. Despite this importance, we struggle to even understand the energy budget of CE evolution. For decades, observed long-period post-CE binaries have been interpreted as evidence of additional energies contributing during CE evolution. We have recently shown that this argument is based on simplified assumptions for all long-period post-CE binaries containing massive white dwarfs (WDs). The only remaining post-CE binary star that has been claimed to require contributions from additional energy sources to understand its formation is KOI 3278.
Aims: Here, we address in detail the potential evolutionary history of KOI 3278. In particular, we investigate whether extra energy sources, such as recombination energy, are indeed required to explain its existence.
Methods: We used the 1D stellar evolution code MESA to carry out binary evolution simulations and searched for potential formation pathways for KOI 3278 that are able to explain its observed properties.
Results: We find that KOI 3278 can be explained if the WD progenitor filled its Roche lobe during a helium shell flash. In this case, the orbital period of KOI 3278 can be reproduced if the CE binding energy is calculated taking into account gravitational energy and thermodynamic internal energy. While the CE evolution that led to the formation of KOI 3278 must have been efficient - that is, most of the available orbital energy must have been used to unbind the CE - recombination energy is not required.
Conclusions: We conclude that currently not a single observed post-CE binary requires one to assume that energy sources other than gravitational and thermodynamic energy are contributing to CE evolution. KOI 3278, however, remains an intriguing post-CE binary as, unlike its siblings, understanding its existence requires highly efficient CE ejection.

Context. Galaxy clusters are the largest gravitating structures in the universe, and their mass assembly is sensitive to the underlying cosmology. Their mass function, baryon fraction, and mass distribution have been used to infer cosmological parameters despite the presence of systematics. However, the complexity of the scaling relations among galaxy cluster properties has never been fully exploited, limiting their potential as a cosmological probe.
Aims: We propose the first machine learning (ML) method using galaxy cluster properties from hydrodynamical simulations in different cosmologies to predict cosmological parameters combining a series of canonical cluster observables, such as gas mass, gas bolometric luminosity, gas temperature, stellar mass, cluster radius, total mass, and velocity dispersion at different redshifts.
Methods: The ML model was trained on mock "measurements" of these observable quantities from Magneticum multi-cosmology simulations to derive unbiased constraints on a set of cosmological parameters. These include the mass density parameter, Ω_m, the power spectrum normalization, σ₈, the baryonic density parameter, Ω_b, and the reduced Hubble constant, h₀.
Results: We tested the ML model on catalogs of a few hundred clusters taken, in turn, from each simulation and found that the ML model can correctly predict the cosmology from where they have been picked. The cumulative accuracy depends on the cosmology, ranging from 21% to 75%. We demonstrate that this is sufficient to derive unbiased constraints on the main cosmological parameters with errors on the order of ~14% for Ω_m, ~8% for σ₈, ~6% for Ω_b, and ~3% for h₀.
Conclusions: This proof-of-concept analysis, though based on a limited variety of multi-cosmology simulations, shows that ML can efficiently map the correlations in the multidimensional space of the observed quantities to the cosmological parameter space and narrow down the probability that a given sample belongs to a given cosmological parameter combination. More large-volume, mid-resolution, multi-cosmology hydro-simulations need to be produced to expand the applicability to a wider cosmological parameter range. However, this first test is exceptionally promising, as it shows that these ML tools can be applied to cluster samples from multiwavelength observations from surveys such as Rubin/LSST, CSST, Euclid, and Roman in optical and near-infrared bands, and eROSITA in X-rays, to the constrain cosmology and effect of baryonic feedback.

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec Integral Field Unit (IFU) spectroscopy of the young Galactic supernova remnant Cassiopeia A (Cas A) to probe the physical conditions for molecular CO formation and destruction in supernova ejecta. We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the re-formation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3–5.5 μm) were obtained toward two representative knots in the NE and S fields that show very different nucleosynthesis characteristics. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V], and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudocontinuum underneath, which is attributed to the high-velocity widths of CO lines. Our results with LTE modeling of CO emission indicate a temperature of ∼1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe.

Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100−1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of 10¹⁰ M_⊙ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs (M_⋆ ≥ 10¹⁰ M_⊙), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems (M_⋆ < 10¹⁰ M_⊙) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings. Our results seem to be supported by observations, and we briefly discuss how this new suite of simulations can help our understanding of the secular evolution of main sequence disc galaxies.

We present the serendipitous discovery of (1) a large double radio relic associated with the galaxy cluster PSZ2 G277.93 + 12.34 and (2) a new odd radio circle, ORC J1027-4422, both found in the same deep MeerKAT 1.3 GHz wide-band radio continuum image. The angular separation of the two arc-shaped cluster relics is ~16 arcmin or ~2.6 Mpc for a cluster redshift of z ≈ 0.158. The thin southern relic, which shows several ridges/shocks including one possibly moving inwards, has a linear extent of ~1.64 Mpc. In contrast, the northern relic is about twice as wide, twice as bright, but only has a largest linear size of ~0.66 Mpc. Complementary SRG/eROSITA X-ray images reveal extended emission from hot intracluster gas between the two relics and around the narrow-angle tail (NAT) radio galaxy PMN J1033-4335 (z ≈ 0.153) located just east of the northern relic. The radio morphologies of the NAT galaxy and the northern relic, which are also detected with the Australian Square Kilometer Array Pathfinder (ASKAP) at 888 MHz, suggest both are moving in the same outward direction. The discovery of ORC J1027-4422 in a different part of the same MeerKAT image makes it the fourth known single ORC. It has a diameter of ~90 arcsec corresponding to 400 kpc at a tentative redshift of z ≈ 0.3 and remains undetected in X-ray emission. Supported by simulations, we discuss similarities between outward moving galaxy and cluster merger shocks as the formation mechanisms for ORCs and radio relics, respectively.

Context. The Kepler high-precision planetary sample has revealed a `radius valley' separating compact super-Earths from sub-Neptunes with lower densities. Super-Earths are generally assumed to be rocky planets that were probably born in situ, while the composition and formation of sub-Neptunes remains debated. Numerous statistical studies have explored planetary and stellar properties and their correlations to provide observational clues. However, no conclusive result on the origin of the radius valley or the composition of sub-Neptunes has been derived to date.
Aims: To provide more constraints, our aim is to investigate the distributions of the orbital spacing of sub-Neptunes and super-Earth planets in Kepler systems and compare their distributions with theoretical predictions of planet pairs of different formation pathways and compositions in synthetic planetary systems.
Methods: Based on the Kepler planetary sample, we derived the distributions of period ratios of sub-Neptune and super-Earth planet pairs. Using synthetic planetary systems generated by the Generation III Bern Model, we also obtained theoretical predictions of period ratio distributions of planet pairs of different compositions and origins.
Results: We find that Kepler sub-Neptune pairs show a significant preference to be near first-order mean motion resonances by a factor of 1.7<sub>−0.3</sub><sup>+0.3</sup>. This is smaller than the model predictions for `water-rich' pairs but larger than that of `water-poor' pairs by confidence levels of ~2σ. Kepler super-Earth pairs show no significant preference for mean motion resonances from a random distribution. The derived normalised fraction of near first-order resonances of actual Kepler super-Earth pairs is consistent with that of synthetic water-poor planet pairs but significantly (≳3σ) smaller than that of synthetic water-rich planet pairs.
Conclusions: The orbital migration has been more important for sub-Neptunes than for super-Earths, suggesting a partial ex situ formation of the former and an origin of the radius valley caused in part by distinct formation pathways. However, the model comparisons also show that sub-Neptunes in Kepler multiple systems are not likely (~2σ) to all be water-rich planets born ex situ but a mixture of the two (in situ and ex situ) pathways. Whereas, Kepler super-Earth planets are predominantly composed of water-poor planets that were born inside the ice line, likely through a series of giant impacts without large-scale migration.

Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust.
Aims: We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales.
Methods: We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits.
Results: We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk.
Conclusions: Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

We describe an attempt of string theoretic derivation of the Gibbons-Hawking entropy. Despite not admitting a de Sitter vacuum, the string theory, by the power of open-close correspondence, captures the Gibbons-Hawking entropy as the entropy of Chan-Paton species on a de Sitter-like state obtained via $D$-branes. Moreover, this derivation sheds a new light at the origin of the area-form, since the equality takes place for a critical 't Hooft coupling for which the species entropy of open strings saturates the area-law unitarity bound.

Observations and high-resolution hydrodynamical simulations indicate that massive star clusters assemble hierarchically from subclusters with a universal power-law cluster mass function. We study the consequences of such assembly for the formation of intermediate-mass black holes (IMBHs) at low metallicities (Z = 0.01 Z_⊙) with our updated N-body code BIFROST based on the hierarchical fourth-order forward integrator. BIFROST integrates few-body systems using secular and regularized techniques including post-Newtonian equations of motion up to order PN3.5 and gravitational-wave recoil kicks for BHs. Single stellar evolution is treated using the fast population synthesis code SEVN. We evolve three cluster assembly regions with N_tot = 1.70-2.35 × 10⁶ stars following a realistic IMF in ~1000 subclusters for t = 50 Myr. IMBHs with masses up to m_• ~ 2200 M_⊙ form rapidly mainly via the collapse of very massive stars (VMSs) assembled through repeated collisions of massive stars followed by growth through tidal disruption events and BH mergers. No IMBHs originate from the stars in the initially most massive clusters. We explain this by suppression of hard massive star binary formation at high velocity dispersions and the competition between core collapse and massive star lifetimes. Later the IMBHs form subsystems resulting in gravitational-wave BH-BH, IMBH-BH, and IMBH-IMBH mergers with a m_• ~ 1000 M_⊙ gravitational-wave detection being the observable prediction. Our simulations indicate that the hierarchical formation of massive star clusters in metal poor environments naturally results in formation of potential seeds for supermassive black holes.

We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\, h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping on to the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall on to its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures.

Lambda- Cold Dark Matter (LambdaCDM) has been successful at explaining the large-scale structures in the universe but faces severe issues on smaller scales when compared to observations. Introducing self-interactions between dark matter particles claims to provide a solution to the small-scale issues in the LambdaCDM simulations while being consistent with the observations at large scales. The existence of the energy region in which these self-interactions between dark matter particles come close to saturating the S-wave unitarity bound can result in the formation of dark matter bound states called darkonium. In this scenario, all the low energy scattering properties are determined by a single parameter, the inverse scattering length gamma. In this work, we set bounds on gamma by studying the impact of darkonium on the observations at direct detection experiments using data from CRESST-III and XENON1T. The exclusion limits on gamma are then subsequently converted to exclusion limits on the self-interaction cross-section and compared with the constraints from astrophysics and N-body simulations.

RR Lyrae stars are excellent tracers of stellar populations for old, metal-poor components in the Milky Way Galaxy and the Local Group. Their luminosities have a metallicity-dependence, but determining spectroscopic [Fe/H] metallicities for RR Lyrae stars, especially at distances outside the solar neighbourhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(II) and Fe(I) abundances, we verify the calibration between the [Fe/H] of RR Lyrae stars from the Calcium triplet. Our calibration is applied to all RR Lyrae stars with Gaia RVS spectra in Gaia DR3 as well as to 80 stars in the inner Galaxy from the BRAVA-RR survey. The co-added Gaia RVS RR Lyrae spectra provide RR Lyrae metallicities with an uncertainty of 0.25~dex, which is a factor of two improvement over the Gaia photometric RR Lyrae metallicities. Within our Galactic bulge RR Lyrae star sample, we find a dominant fraction with low energies without a prominent rotating component. Due to the large fraction of such stars, we interpret these stars as belonging to the $in-situ$ metal-poor Galactic bulge component, although we can not rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles.