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.

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.

This review describes the duality between color and kinematics and its applications, with the aim of gaining a deeper understanding of the perturbative structure of gauge and gravity theories. We emphasize, in particular, applications to loop-level calculations, the broad web of theories linked by the duality and the associated double-copy structure, and the issue of extending the duality and double copy beyond scattering amplitudes. The review is aimed at doctoral students and junior researchers both inside and outside the field of amplitudes and is accompanied by various exercises.

RR Lyrae stars (RRLs) are excellent tracers of stellar populations for old, metal-poor components in the the Milky Way and the Local Group. Their luminosities have a metallicity dependence, but determining spectroscopic [Fe/H] metallicities for RRLs, especially at distances outside the solar neighborhood, 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 RRLs from the calcium triplet. Our calibration is applied to all RRLs with Gaia Radial Velocity Spectrometer (RVS) spectra in Gaia DR3 and to 80 stars in the inner Galaxy from the BRAVA-RR survey. The coadded Gaia RVS RRL spectra provide RRL metallicities with an uncertainty of 0.25 dex, which is a factor of two improvement over the Gaia photometric RRL metallicities. Within our Galactic bulge RRL 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 cannot rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles.

A central requirement for the faithful implementation of large-scale lattice gauge theories (LGTs) on quantum simulators is the protection of the underlying gauge symmetry. Recent advancements in the experimental realizations of large-scale LGTs have been impressive, albeit mostly restricted to Abelian gauge groups. Guided by this requirement for gauge protection, we propose an experimentally feasible approach to implement large-scale non-Abelian <inline-formula><mml:math display="inline" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>SU</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>N</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math display="inline" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi mathvariant="normal">U</mml:mi></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>N</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> LGTs with dynamical matter in <inline-formula><mml:math display="inline" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>d</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mrow><mml:mrow><mml:mi mathvariant="normal">D</mml:mi></mml:mrow></mml:mrow></mml:math></inline-formula>, enabled by two-body spin-exchange interactions realizing local emergent gauge-symmetry stabilizer terms. We present two concrete proposals for <inline-formula><mml:math display="inline" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mn>2</mml:mn><mml:mo>+</mml:mo><mml:mn>1</mml:mn><mml:mrow><mml:mrow><mml:mi mathvariant="normal">D</mml:mi></mml:mrow></mml:mrow><mml:mspace width="0.2em"></mml:mspace><mml:mi>SU</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>2</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math display="inline" overflow="scroll" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi mathvariant="normal">U</mml:mi></mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>2</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> LGTs, including dynamical bosonic matter and induced plaquette terms, that can be readily implemented in current ultracold-molecule and next-generation ultracold-atom platforms. We provide numerical benchmarks showcasing experimentally accessible dynamics, and demonstrate the stability of the underlying non-Abelian gauge invariance. We develop a method to obtain the effective gauge-invariant model featuring the relevant magnetic plaquette and minimal gauge-matter coupling terms. Our approach paves the way towards near-term realizations of large-scale non-Abelian quantum link models in analog quantum simulators.

We present phase-corrected photometric measurements of 88 Cepheid variables in the core of the Small Magellanic Cloud (SMC), the first sample obtained with the Hubble Space Telescope's (HST) Wide Field Camera 3, in the same homogeneous photometric system as past measurements of all Cepheids on the SH0ES distance ladder. We limit the sample to the inner core and model the geometry to reduce errors in prior studies due to the nontrivial depth of this cloud. Without crowding present in ground-based studies, we obtain an unprecedentedly low dispersion of 0.102 mag for a period–luminosity (P–L) relation in the SMC, approaching the width of the Cepheid instability strip. The new geometric distance to 15 late-type detached eclipsing binaries in the SMC offers a rare opportunity to improve the foundation of the distance ladder, increasing the number of calibrating galaxies from three to four. With the SMC as the only anchor, we find H ₀ = 74.1 ± 2.1 km s^‑1 Mpc^‑1. Combining these four geometric distances with our HST photometry of SMC Cepheids, we obtain H ₀ = 73.17 ± 0.86 km s^‑1 Mpc^‑1. By including the SMC in the distance ladder, we also double the range where the metallicity ([Fe/H]) dependence of the Cepheid P–L relation can be calibrated, and we find γ = ‑0.234 ± 0.052 mag dex^‑1. Our local measurement of H ₀ based on Cepheids and Type Ia supernovae shows a 5.8σ tension with the value inferred from the cosmic microwave background assuming a Lambda cold dark matter (ΛCDM) cosmology, reinforcing the possibility of physics beyond ΛCDM.

We introduce a block encoding method for mapping discrete subgroups to qubits on a quantum computer. This method is applicable to general discrete groups, including crystal-like subgroups such as <inline-formula><mml:math display="inline"><mml:mrow><mml:mi mathvariant="double-struck">BI</mml:mi></mml:mrow></mml:math></inline-formula> of <inline-formula><mml:math display="inline"><mml:mi>S</mml:mi><mml:mi>U</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>2</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:mi mathvariant="double-struck">V</mml:mi></mml:math></inline-formula> of <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>. We detail the construction of primitive gates—the inversion gate, the group multiplication gate, the trace gate, and the group Fourier gate—utilizing this encoding method for <inline-formula><mml:math display="inline"><mml:mrow><mml:mi mathvariant="double-struck">BT</mml:mi></mml:mrow></mml:math></inline-formula> and for the first time <inline-formula><mml:math display="inline"><mml:mrow><mml:mi mathvariant="double-struck">BI</mml:mi></mml:mrow></mml:math></inline-formula> group. We also provide resource estimations to extract the gluon viscosity. The inversion gates for <inline-formula><mml:math display="inline"><mml:mrow><mml:mi mathvariant="double-struck">BT</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:mrow><mml:mi mathvariant="double-struck">BI</mml:mi></mml:mrow></mml:math></inline-formula> are benchmarked on the Baiwang quantum computer with estimated fidelities of <inline-formula><mml:math display="inline"><mml:msubsup><mml:mn>40</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>4</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msubsup><mml:mo>%</mml:mo></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:msubsup><mml:mn>4</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>3</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msubsup><mml:mo>%</mml:mo></mml:math></inline-formula>, respectively.

When a star is described as a spectral class G2V, we know its approximate mass, temperature, age, and size. At more than 5,700 exoplanets discovered, it is a natural developmental step to establish a classification for them, such as for example, the Harvard classification for stars. This exoplanet classification has to be easily interpreted and present the most relevant information about them and divides them into groups based on certain characteristics. We propose an exoplanet classification, which using an easily readable code, may inform you about a exoplanet's main characteristics. The suggested classification code contains four parameters by which we can quickly determine the range of temperature, mass, density and their eccentricity. The first parameter concerns the mass of an exoplanet in the form of the units of the mass of other known planets, where e.g. M represents the mass of Mercury, E that of Earth, N Neptune, or J Jupiter. The second parameter is the mean Dyson temperature of the extoplanet's orbit, for which we established four main classes: F represents the Frozen class, W the Water class, G the Gaseous class, and R the Roaster class. The third parameter is eccentricity and the fourth parameter is surface attribute which is defined as the bulk density of the exoplanet, where g represents a gaseous planet, w - water planet, t - terrestrial planet, i - iron planet and s - super dense planet. The classification code for Venus, could be EG0t (E - mass in the range of the mass of the Earth, G - Gaseous class, temperature in the range from 450 to 1000 K, 0 - circular or nearly circular orbit, t - terrestrial surface), for Earth it could be EW0t (W - Water class - a possible Habitable zone). This classification is very helpful in, for example, quickly delimiting if a planet can be found in the Habitable zone; if it is terrestrial or not.

We present UV–optical–near-infrared observations and modeling of supernova (SN) 2024ggi, a type II supernova (SN II) located in NGC 3621 at 7.2 Mpc. Early-time ("flash") spectroscopy of SN 2024ggi within +0.8 days of discovery shows emission lines of H I, He I, C III, and N III with a narrow core and broad, symmetric wings (i.e., "IIn-like") arising from the photoionized, optically thick, unshocked circumstellar material (CSM) that surrounded the progenitor star at shock breakout (SBO). By the next spectral epoch at +1.5 days, SN 2024ggi showed a rise in ionization as emission lines of He II, C IV, N IV/V, and O V became visible. This phenomenon is temporally consistent with a blueward shift in the UV–optical colors, both likely the result of SBO in an extended, dense CSM. The IIn-like features in SN 2024ggi persist on a timescale of t _IIn = 3.8 ± 1.6 days, at which time a reduction in CSM density allows the detection of Doppler-broadened features from the fastest SN material. SN 2024ggi has peak UV–optical absolute magnitudes of M _w2 = ‑18.7 mag and M _g = ‑18.1 mag, respectively, that are consistent with the known population of CSM-interacting SNe II. Comparison of SN 2024ggi with a grid of radiation hydrodynamics and non–local thermodynamic equilibrium radiative-transfer simulations suggests a progenitor mass-loss rate of <inline-formula> <mml:math overflow="scroll"><mml:mover accent="true"><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>̇</mml:mo></mml:mrow></mml:mover><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mspace width="0.25em"></mml:mspace><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>⊙</mml:mo></mml:mrow></mml:msub></mml:math> </inline-formula> yr^‑1 (v _w = 50 km s^‑1), confined to a distance of r < 5 × 10¹⁴ cm. Assuming a wind velocity of v _w = 50 km s^‑1, the progenitor star underwent an enhanced mass-loss episode in the last ∼3 yr before explosion.

We study the contribution of large scalar perturbations sourced by a sharp feature during cosmic inflation to the stochastic gravitational wave background (SGWB), extending our previous work to include the SGWB sourced during the inflationary era. We focus in particular on three-field inflation, since the third dynamical field is the first not privileged by the perturbations' equations of motion and allows a more direct generalization to $N$-field inflation. For the first time, we study the three-field isocurvature perturbations sourced during the feature and include the effects of isocurvature masses. In addition to a two-field limit, we find that the third field's dynamics during the feature can source large isocurvature transients which then later decay, leaving an inflationary-era-sourced SGWB as their only observable signature. We find that the inflationary-era signal shape near the peak is largely independent of the number of dynamical fields and has a greatly enhanced amplitude sourced by the large isocurvature transient, suppressing the radiation-era contribution and opening a new window of detectable parameter space with small adiabatic enhancement. The largest enhancements we study could easily violate backreaction constraints, but much of parameter space remains under perturbative control. These SGWBs could be visible in LISA and other gravitational wave experiments, leaving an almost universal signature of sharp features during multi-field inflation, even when the sourcing isocurvature decays to unobservability shortly afterwards.

The Hubble parameter, H ₀, is not an univocally defined quantity: It relates redshifts to distances in the near Universe, but it is also a key parameter of the ΛCDM standard cosmological model. As such, H ₀ affects several physical processes at different cosmic epochs and multiple observables. We have counted more than a dozen H ₀s that are expected to agree if (a) there are no significant systematics in the data and their interpretation and (b) the adopted cosmological model is correct. With few exceptions (proverbially confirming the rule), these determinations do not agree at high statistical significance; their values cluster around two camps: the low (68 km s¹ Mpc¹) and high (73 km s¹ Mpc¹) camps. It appears to be a matter of anchors. The shape of the Universe expansion history agrees with the model; it is the normalizations that disagree. Beyond systematics in the data/analysis, if the model is incorrect, there are only two viable ways to "fix" it: by changing the early time (z ≳ 1,100) physics and, thus, the early time normalization or by a global modification, possibly touching the model's fundamental assumptions (e.g., homogeneity, isotropy, gravity). None of these three options has the consensus of the community. The research community has been actively looking for deviations from ΛCDM for two decades; the one we might have found makes us wish we could put the genie back in the bottle.

Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be further improved with time-stepping schemes that use input from cosmological perturbation theory. Specifically, we introduce a class of time-stepping schemes derived by matching the particle trajectories in a single leapfrog/Verlet drift-kick-drift step to those predicted by Lagrangian perturbation theory (LPT). As a corollary, these schemes exactly yield the analytic Zel'dovich solution in 1D in the pre-shell-crossing regime (i.e. before particle trajectories cross). One representative of this class is the popular 'FastPM' scheme by Feng et al. 2016[1], which we take as our baseline. We then construct more powerful LPT-inspired integrators and show that they outperform FastPM and standard integrators in fast simulations in two and three dimensions with <mml:math altimg="si1.svg"><mml:mi mathvariant="script">O</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo linebreak="badbreak" linebreakstyle="after">‑</mml:mo><mml:mn>100</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:math> timesteps, requiring fewer steps to accurately reproduce the power spectrum and bispectrum of the density field. Furthermore, we demonstrate analytically and numerically that, for any integrator, convergence is limited in the post-shell-crossing regime (to order <mml:math altimg="si2.svg"><mml:mfrac bevelled="true"><mml:mrow><mml:mn>3</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mfrac></mml:math> for planar-wave collapse), owing to the lacking regularity of the acceleration field, which makes the use of high-order integrators in this regime futile. Also, we study the impact of the timestep spacing and of a decaying mode present in the initial conditions. Importantly, we find that symplecticity of the integrator plays a minor role for fast approximate simulations with a small number of timesteps.

We combine stellar orbits with the abundances of the heavy, r-process element europium and the light, $\alpha$-element, silicon to separate in situ and accreted populations in the Milky Way (MW) across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in situ, the levels of [Eu/Si] are lower. These systematically different levels of [Eu/Si] in the MW and the accreted halo imply that the scatter in [Eu/$\alpha$] within a single galaxy is smaller than previously thought. At the lowest metallicities, we find that both accreted and in situ populations trend down in [Eu/Si], consistent with enrichment via neutron star mergers. Through compiling a large data set of abundances for 54 globular clusters (GCs), we show that differences in [Eu/Si] extend to populations of in situ/accreted GCs. We interpret this consistency as evidence that in r-process elements GCs trace the star formation history of their hosts, motivating their use as sub-Gyr timers of galactic evolution. Furthermore, fitting the trends in [Eu/Si] using a simple galactic chemical evolution model, we find that differences in [Eu/Si] between accreted and in situ MW field stars cannot be explained through star formation efficiency alone. Finally, we show that the use of [Eu/Si] as a chemical tag between GCs and their host galaxies extends beyond the Local Group, to the halo of M31 - potentially offering the opportunity to do Galactic Archaeology in an external galaxy.

While LCDM cosmology is the most successful cosmological model at our disposal today, being able to explain most of the observed phenomena, it has been challenged by more and more tensions. One of the greatest, both in terms of numerical tension and of the importance of the parameter measured, is the infamous Hubble tension. This refers to the disagreement between measurements of the Hubble constant, which describes the rate of expansion of the Universe, a cornerstone of our cosmological understanding. In recent years the methods for measuring H0 have grown in number and sophistication, and yet, as the uncertainties of the measurements have decreased, the tension has not been solved; in fact, it has increased.

Such methods can be roughly divided between "early" and "late" probes of H0, approximately referring to the time of origin of the phenomenon observed. While "early" probes, based for example on the cosmic microwave background, are strongly dependent on the assumed cosmology, "late" probes are generally model-independent but are more susceptible to systematic errors in the measurements. In this context, the time delay cosmographic method is a "late" time probe which can measure H0 directly, without requiring any calibration. This analysis is based on the well-tested general relativity phenomenon of strong gravitational lensing. Given a background variable source and a foreground strong gravitational lens, the time delay between the multiple lensed images can be measured by monitoring and analysing their luminosity over time. A separate modelling analysis of the system can then constrain the mass profile of the lens. The two combined information can then be used to constrain the Hubble constant. In this work, I implemented this analysis based on Hubble Space Telescope archival data and a dedicated observational campaign from the 2.1-meter telescope at Wendelstein. I employed the space-based data by taking advantage of the multiple filters available and their higher resolution to model the lens mass, obtaining a result with 3% precision on the Fermat potential.

I instead used the data from the Wendelstein observational campaign to produce the lightcurves of the image and analyse them in order to constrain the time delay, which was obtained with a precision ranging from 8% to 15% depending on the image pair.

I then combined the results following a Bayesian approach, reaching a constraint on H0 of 71.3+5.0 -4.5 km/(s*Mpc) with a precision ~6.7% considering random uncertainty.

Notably, this work has been mostly independent of major collaborations, such as TDCOSMO, thus providing an unbiased validation of the methodology. Furthermore, the result is proof of the capabilities of the Wendelstein observatory, which should be considered a reliable asset for time delay cosmography or similar projects that require high-sampling, high-quality data.

This thesis reports a determination of the branching fraction (B) and CP-violating charge asymmetry (ACP) of the three-body decay B 0 -> K + π − π 0 at the Belle II experiment. In addition to the inclusive B and ACP, i.e. for B0 -> K+ π− π0 decays, we measure B and ACP exclusively for individual two-body resonances appearing in the K+ π− π0 system. To this end, we employ a model dependent Dalitz plot analysis, including the seven dominant intermediate resonances and a non-resonant contribution. The analyzed data were recorded between 2019 and 2022 and correspond to an integrated luminosity of 362 fb^−1 produced in e+ e− collisions at the Y(4S) resonance by the SuperKEKB collider containing 387 × 106 pairs of bottom-antibottom mesons. The analysis is developed on simulated data and control mode data solely. The branching fractions and CP-asymmetries are extracted in a four-dimensional extended maximum likelihood fit. As the analysis is still under Belle II internal review, we blind central values and state uncertainties only. We measure the branching fraction and CP -violating charge asymmetry inclusively as well as exclusively for the channels B0 -> K∗(892)+ π− , B0 -> K∗(892)0 π0 , B0 -> ρ(770)− K+ , B0 -> (Kπ)∗+ π-, B0 -> (Kπ)0 π0, B0 -> ρ(1450)- K+, B0 -> ρ(1700)- K+ und B0 -> K+ π- π0 non-resonant.

This thesis presents the first model dependent Dalitz plot analysis at Belle II. We achieve uncertainties on par with known determinations. The B 0 -> K ∗ (892)π modes will serve as inputs for a sum rule based on isospin to probe the Standard Model.

Alongside the search of a potential origin of life on Earth, it is clear, that all processes leading towards life's first molecules had to be compatible with the geo-physical circumstances provided by Earth's prebiotic environment. This implies that life not only had to emerge from simple chemical building blocks but also from the highly diluted solutions posed by the prebiotic oceans. As all living systems require information storing molecules, nucleic acids are assumed to be the starting point for molecular evolution. However, the mechanisms of how molecules could accumulate in early geological settings are poorly understood and especially the formation of nucleic acids primary building blocks poses one of the biggest challenges in the research field.[...]

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.

While it is possible to estimate the dark matter density at the Sun distance from the galactic center, this does not give information on actual dark matter density in the Solar system. There can be considerable local enhancement of dark matter density in the vicinity of gravitating centers, including the Sun, the Earth, as well as other planets in the solar system. Generic mechanisms for the formation of such halos were recently elucidated. In this work, we studies the possible halo dark matter overdensities and corresponding dark matter masses allowed for various objects in the solar system. We explore spacecraft missions to detect such halos with instruments such as quantum clocks, atomic and molecular spectrometers designed to search for fast (tens of hertz to gigahertz) oscillations of fundamental constants, highly sensitive comagnetometers, and other quantum sensors and sensor networks.

Astronomical and cosmological observations indicate that dark matter should interact very weakly with the electromagnetic radiation. Nevertheless, the existence of such interactions is not precluded by observations nor by theoretical considerations. A promising approach to probe the dark matter electromagnetic properties is through the search of photon-mediated dark matter-nucleus interactions in direct detection experiments. In this paper we present a simple methodology to calculate the scattering rate in a direct detection experiment for given values of the dark matter electric charge, charge radius, electric- and magnetic- dipole moments and anapole moment. In our work we include contributions to the scattering from nuclear recoils and from the Migdal effect. We finally apply this formalism to determine exclusion limits on the five electromagnetic interactions using data from XENON1T, LZ, PICO-60 and DS50 experiments, and we discuss the implications for a simplified dark matter model with t-channel mediators.

The photon parton distribution function (PDF) of the proton is crucial for precise comparisons of LHC cross sections with theoretical predictions. However, it was previously affected by very large uncertainties of around ${\cal O}(100\%)$ or dependent upon phenomenologically inspired models. In the paper~\cite{Manohar:2016nzj}, we demonstrated how the photon PDF could be determined using the proton structure functions $F_2$ and $F_L$ measured in electron--proton scattering experiments. We provided an explicit formula for the PDF, which can be systematically improved order by order in perturbation theory. We obtained a photon PDF with errors $\lesssim 2$\% for $10^{-4} < x < 0.1$. Here, we recall the underlying idea and method used to obtain this result, as well as the progress made since then.

The Gaia mission has provided highly accurate observations that have significantly reduced the scatter in the colour-magnitude diagrams of open clusters. As a result of the improved isochrone sequence of the open cluster M67, we have created new stellar models that avoid commonly used simplifications in 1D stellar modelling, such as mass-independent core overshooting and a constant mixing length parameter. This has enabled us to deliver a precise isochrone specifically designed for M67, available for download. We follow a commonly used qualitative approach to adjust the input physics to match the well-defined colour-magnitude sequence, and we test the model-predicted masses against a known eclipsing binary system at the main sequence turnoff of the cluster. Despite using improvements in photometry and stellar physics we cannot match the masses of both binary components with the same theoretical isochrone. A $\chi ^{2}$-based isochrone fitting approach using our preferred input physics results in a cluster age of $3.95^{+ 0.16}_{- 0.15}$ Gyr.

If the gradient of a probability distribution on a landscape of vacua aligns with the variation of some fundamental parameter, the parameter may be likely to take some non-generic value. Such non-generic values can be associated to critical boundaries, where qualitative changes of the landscape properties happen, or an anthropic bound is located. Assuming the standard volume-weighted and the local probability measures, we discuss ordered landscapes which can produce several types of the aligned probability gradients. The resulting values of the gradients are defined by the "closeness" of a given vacuum to the highest- or the lowest-energy vacuum. Using these ingredients we construct a landscape scanning independently the Higgs mass and the cosmological constant (CC). The probability gradient pushes the Higgs mass to its observed value, where a structural change of the landscape takes place, while the CC is chosen anthropically.

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.

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.

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.

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 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.

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.

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.

Context. Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipation, which affects only very small scales.
Aims: We want to test and quantify the reconnection hypothesis in decaying hydromagnetic turbulence.
Methods: We performed high-resolution numerical simulations with zero net magnetic helicity using the PENCIL CODE with up to 2048³ mesh points and relate the decay time to the Alfvén time for different resistivities and viscosities.
Results: The decay time is found to be longer than the Alfvén time by a factor that increases with increasing Lundquist number to the 1/4 power. The decay exponent is as expected from the conservation of the Hosking integral, but a timescale dependence on resistivity is unusual for developed turbulence and not found for hydrodynamic turbulence. In two dimensions, the Lundquist number dependence is shown to be leveling off above values of ≈25 000, independently of the value of the viscosity.
Conclusions: Our numerical results suggest that resistivity effects have been overestimated in earlier work. Instead of reconnection, it may be the magnetic helicity density in smaller patches that is responsible for the resistively slow decay. The leveling off at large Lundquist number cannot currently be confirmed in three dimensions.

The nature of dark matter (DM) and its interaction with the Standard Model (SM) is one of the biggest open questions in physics nowadays. The vast majority of theoretically motivated ultralight-DM (ULDM) models predict that ULDM couples dominantly to the SM strong/nuclear sector. This coupling leads to oscillations of nuclear parameters that are detectable by comparing clocks with different sensitivities to these nature's constants. Vibrational transitions of molecular clocks are more sensitive to a change in the nuclear parameters than the electronic transitions of atomic clocks. Here, we propose the iodine molecular ion, I₂⁺, as a sensitive detector for such a class of ULDM models. The iodine's dense spectrum allows us to match its transition frequency to that of an optical atomic clock (Ca⁺) and perform correlation spectroscopy between the two clock species. With this technique, we project a few-orders-of-magnitude improvement over the most sensitive clock comparisons performed to date. We also briefly consider the robustness of the corresponding "Earth-bound" under modifications of the Z_N -QCD axion model.

For non-relativistic thermal dark matter, close-to-threshold effects largely dominate the evolution of the number density for most of the times after thermal freeze-out, and hence affect the cosmological relic density. A precise evaluation of the relevant interaction rates in a thermal medium representing the early universe includes accounting for the relative motion of the dark matter particles and the thermal medium. We consider a model of dark fermions interacting with a plasma of dark gauge bosons, which is equivalent to thermal QED. The temperature is taken to be smaller than the dark fermion mass and the inverse of the typical size of the dark fermion-antifermion bound states, which allows for the use of non-relativistic effective field theories. For the annihilation cross section, bound-state formation cross section, bound-state dissociation width and bound-state transition width of dark matter fermion-antifermion pairs, we compute the leading recoil effects in the reference frame of both the plasma and the center-of-mass of the fermion-antifermion pair. We explicitly verify the Lorentz transformations among these quantities. We evaluate the impact of the recoil corrections on the dark matter energy density. Our results can be directly applied to account for the relative motion of quarkonia in the quark-gluon plasma formed in heavy-ion collisions. They may be also used to precisely assess thermal effects in atomic clocks based on atomic transitions; the present work provides a first field theory derivation of time dilation for these processes in vacuum and in a medium.

We use functional methods to match the two-Higgs-doublet model with heavy scalars in the nondecoupling regime to the appropriate nonlinear effective field theory, which takes the form of an electroweak chiral Lagrangian (HEFT). The effective Lagrangian is derived to leading order in the chiral counting. This includes the loop induced h →γ γ and h →Z γ local terms, which enter at the same chiral order as their counterparts in the Standard Model. An algorithm is presented that allows us to compute the coefficient functions to all orders in h . Some of the all-orders results are given in closed form. The parameter regimes for decoupling, nondecoupling, and alignment scenarios in the effective field theory context and some phenomenological implications are briefly discussed.

High-frequency gravitational waves (HFGWs) are predicted in various exotic scenarios involving both cosmological and astrophysical sources. These elusive signals have recently sparked the interest of a diverse community of researchers, due to the possibility of HFGW detection in the laboratory through graviton-photon conversion in strong magnetic fields. Notable examples include the redesign of the resonant cavities currently under development to detect the cosmic axion. In this work, we derive the sensitivities of some existing and planned resonant cavities to detect a HFGW background. As a concrete scenario, we consider the collective signals that originate from the merging of compact objects, such as two primordial black holes (PBHs) in the asteroid mass window. Our findings improve over existing work by explicitly discussing and quantifying the loss in the experimental reach due to the actual coherence of the source. We elucidate on the approach we adopt in relation with recent literature on the topic. Most notably, we give a recipe for the estimate of the stochastic background that focuses on the presence of the signal in the cavity at all times and showing that, in the relevant PBH mass region, the signal is dominated by coherent binary mergers.

Flavor-dependent neutrino transport is described by a well-known kinetic equation for occupation-number matrices in flavor space. However, in the context of fast flavor conversion, we identify an unforeseen predicament: the pivotal self-induced exponential growth of small inhomogeneities strongly violates conservation of neutrino-neutrino refractive energy. We prove that it is traded with the huge reservoir of neutrino kinetic energy through gradients of neutrino flavor coherence (the off-diagonal piece of the flavor density matrix) and derive the missing gradient terms. The usual equations remain sufficient to describe flavor evolution, at the cost of renouncing energy conservation, which cannot play any role in explaining the numerically observed final state.

We propose a simple model that can alleviate the H₀ tension while remaining consistent with big bang nucleosynthesis (BBN). It is based on a dark sector described by a standard Lagrangian featuring a S U (N ) gauge symmetry with N ≥3 and a massive scalar field with a quartic coupling. The scalar acts as a dark Higgs leading to spontaneous symmetry breaking S U (N )→S U (N -1 ) via a first-order phase transition à la Coleman-Weinberg. This setup naturally realizes previously proposed scenarios featuring strongly interacting dark radiation (SIDR) with a mass threshold within hot new early dark energy. For a wide range of reasonable model parameters, the phase transition occurs between the BBN and recombination epochs and releases a sufficient amount of latent heat such that the model easily respects bounds on extra radiation during BBN while featuring a sufficient SIDR density around recombination for increasing the value of H₀ inferred from the cosmic microwave background. Our model can be summarized as a natural mechanism providing two successive increases in the effective number of relativistic degrees of freedom after BBN but before recombination Δ N_BBN→Δ N_NEDE→Δ N_IR alleviating the Hubble tension. The first step is related to the phase transition, and the second is related to the dark Higgs becoming nonrelativistic. This setup predicts further signatures, including a stochastic gravitational wave background and features in the matter power spectrum that can be searched for with future pulsar timing and Lyman-α forest measurements.

We examine the physical significance of torsion co-cycles in the cohomology of a projective Calabi-Yau three-fold for the (2,2) superconformal field theory (SCFT) associated to the non-linear sigma model with such a manifold as a target space. There are two independent torsion subgroups in the cohomology. While one is associated to an orbifold construction of the SCFT, the other encodes the possibility of turning on a topologically non-trivial flat gerbe for the NS-NS B-field. Inclusion of these data enriches mirror symmetry by providing a refinement of the familiar structures and points to a generalization of the duality symmetry, where the topology of the flat gerbe enters on the same footing as the topology of the underlying manifold.

In this paper we present a novel method to extract information on hadron-hadron interactions using for the first time femtoscopic data to constrain the low-energy constants of a QCD effective Lagrangian. This method offers a new way to investigate the nonperturbative regime of QCD in sectors where scattering experiments are not feasible, such as the multistrange and charm ones. As an example of its application, we use the very precise K^-Λ correlation function data, recently measured in p p collisions at LHC, to constrain the strangeness S =-2 meson-baryon interaction. The model obtained delivers new insights on the molecular nature of the Ξ (1620 ) and Ξ (1690 ) states.

We present a novel method for identifying transients suitable for both strong signal-dominated and background-dominated objects. By employing the unsupervised machine learning algorithm known as expectation maximization, we achieve computing time reductions of over 10⁴ on a single CPU compared to conventional brute-force methods. Furthermore, this approach can be readily extended to analyze multiple flares. We illustrate the algorithm's application by fitting the IceCube neutrino flare of TXS 0506+056.