The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein-Boltzmann solvers CLASS and CAMB in the context of Euclid. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of other forecasting methods and tools. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collaboration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein-Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov Chains with MontePython while using the exact same mock likelihoods. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models.
Using the relativistic Hartree-Bogoliubov framework with separable pairing force coupled with the latest covariant density functionals, i.e., PC-L3R, PC-X, DD-PCX, and DD-MEX, we systematically explore the ground-state properties of all isotopes of Z=8-110. These properties consist of the binding energies, one- and two-neutron separation energies (Sn and S2n), root-mean-square radius of matter, of neutron, of proton, and of charge distributions, Fermi surfaces, ground-state spins and parities. We then predict the edges of nuclear landscape and bound nuclei for the isotopic chains from oxygen (Z=8) to darmstadtium (Z=110) based on these latest covariant density functionals. The number of bound nuclei predicted by PC-L3R, PC-X, DD-PCX, and DD-MEX, are 9004, 9162, 6799, and 7112, respectively. The root-mean-square deviations of Sn (S2n) yielded from PC-L3R, PCX, DD-PCX, and DD-MEX are 0.962 (1.300) MeV, 0.920 (1.483) MeV, 0.993 (1.753) MeV, and 1.010 (1.544) MeV, respectively. The root-mean-square deviations of charge radius distributions of comparing the available experimental values with the theoretical counterparts resulted from PC-L3R, PC-X, DD-PCX, and DD-MEX are 0.035 fm, 0.037 fm, 0.035 fm, and 0.034 fm, respectively. We notice pronounced differences between the empirical and theoretical root-mean-square radii of neutron at nuclei near the neutron drip line of the Mg, Ca, and Kr isotopic chains, suggesting the possible existence of the halo or giant halo phenomena.
A relativistic microscopic optical model potential for nucleon-nucleus scattering is developed based on the \emph{ab initio} relativistic Brueckner-Hartree-Fock (RBHF) theory with the improved local density approximation, which is abbreviated as the RBOM potential. Both real and imaginary parts of the single-particle potentials in symmetric and asymmetric nuclear matter at various densities are determined uniquely in the full Dirac space. The density distributions of the target nuclei are calculated by the covariant energy density functional theory with the density functional PC-PK1. The central and spin-orbit terms of the optical potentials are quantitatively consistent with the relativistic phenomenological optical potentials. The performance of the RBOM potential is evaluated by considering proton scattering with incident energy E≤200 MeV on five target nuclei, $\prescript{208}{}{\text{Pb}}$, $\prescript{120}{}{\text{Sn}}$, $\prescript{90}{}{\text{Zr}}$, $\prescript{48}{}{\text{Ca}}$, and $\prescript{40}{}{\text{Ca}}$. Scattering observables including the elastic scattering angular distributions, analyzing powers, spin rotation functions, and reaction cross sections are analyzed. Theoretical predictions show good agreements with the experimental data and the results derived from phenomenological optical potentials. We anticipate that the RBOM potential can provide reference for other phenomenological and microscopic optical model potentials, as well as reliable descriptions for nucleon scattering on exotic nuclei in the era of rare-isotope beams.
Self-assembly is a fundamental concept in biology and of significant interest to nanotechnology. Significant progress has been made in characterizing and controlling the properties of the resulting structures, both experimentally and theoretically. However, much less is known about kinetic constraints and determinants of dynamical properties like time efficiency, although these constraints can become severe limiting factors of self-assembly processes. Here we investigate how the time efficiency and other dynamical properties of reversible self-assembly depend on the morphology (shape) of the building blocks for systems in which the binding energy between the constituents is large. As paradigmatic examples, we stochastically simulate the self-assembly of constituents with triangular, square, and hexagonal morphology into two-dimensional structures of a specified size. We find that the constituents' morphology critically determines the assembly time and how it scales with the size of the target structure. Our analysis reveals three key structural parameters defined by the morphology: The nucleation size and attachment order, which describe the effective order of the chemical reactions by which clusters nucleate and grow, respectively, and the growth exponent, which determines how the growth rate of an emerging structure scales with its size. Using this characterization, we formulate an effective theory of the self-assembly kinetics, which we show exhibits an inherent scale invariance. This allows us to identify general scaling laws that describe the minimal assembly time as a function of the size of the target structure. We show how these insights on the kinetics of self-assembly processes can be used to design assembly schemes that could significantly increase the time efficiency and robustness of artificial self-assembly processes.
We present and validate the catalog of Lyman-α forest fluctuations for 3D analyses using the Early Data Release (EDR) from the Dark Energy Spectroscopic Instrument (DESI) survey. We used 88,511 quasars collected from DESI Survey Validation (SV) data and the first two months of the main survey (M2). We present several improvements to the method used to extract the Lyman-α absorption fluctuations performed in previous analyses from the Sloan Digital Sky Survey (SDSS). In particular, we modify the weighting scheme and show that it can improve the precision of the correlation function measurement by more than 20%. This catalog can be downloaded from this https URL and it will be used in the near future for the first DESI measurements of the 3D correlations in the Lyman-α forest.
At z≲1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T∼104 K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with T>105K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, T=T0(ρ/ρ¯)γ−1. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T0,γ], based on the b-NHI distribution of the Lyman-α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T0, γ] can still be reliably measured from the b-NHI distribution at z=0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z=0.1. The results suggest that the underlying [T0,γ] of both simulations can be recovered with biases as low as |Δlog(T0/K)|≲0.05 dex, |Δγ|≲0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the b-NHI distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.
We present the angular diameter distance measurement obtained with the Baryonic Acoustic Oscillation feature from galaxy clustering in the completed Dark Energy Survey, consisting of six years (Y6) of observations. We use the Y6 BAO galaxy sample, optimized for BAO science in the redshift range 0.6<z<1.2, with an effective redshift at zeff=0.85 and split into six tomographic bins. The sample has nearly 16 million galaxies over 4,273 square degrees. Our consensus measurement constrains the ratio of the angular distance to sound horizon scale to DM(zeff)/rd = 19.51±0.41 (at 68.3% confidence interval), resulting from comparing the BAO position in our data to that predicted by Planck ΛCDM via the BAO shift parameter α=(DM/rd)/(DM/rd)Planck. To achieve this, the BAO shift is measured with three different methods, Angular Correlation Function (ACF), Angular Power Spectrum (APS), and Projected Correlation Function (PCF) obtaining α= 0.952±0.023, 0.962±0.022, and 0.955±0.020, respectively, which we combine to α= 0.957±0.020, including systematic errors. When compared with the ΛCDM model that best fits Planck data, this measurement is found to be 4.3% and 2.1σ below the angular BAO scale predicted. To date, it represents the most precise angular BAO measurement at z>0.75 from any survey and the most precise measurement at any redshift from photometric surveys. The analysis was performed blinded to the BAO position and it is shown to be robust against analysis choices, data removal, redshift calibrations and observational systematics.
In this paper we present and validate the galaxy sample used for the analysis of the baryon acoustic oscillation (BAO) signal in the Dark Energy Survey (DES) Y6 data. The definition is based on a color and redshift-dependent magnitude cut optimized to select galaxies at redshifts higher than 0.6, while ensuring a high-quality photo-z determination. The optimization is performed using a Fisher forecast algorithm, finding the optimal i-magnitude cut to be given by i<19.64+2.894zph. For the optimal sample, we forecast an increase in precision in the BAO measurement of ∼25% with respect to the Y3 analysis. Our BAO sample has a total of 15,937,556 galaxies in the redshift range 0.6<zph<1.2, and its angular mask covers 4,273.42 deg2 to a depth of i=22.5. We validate its redshift distributions with three different methods: directional neighborhood fitting algorithm (DNF), which is our primary photo-z estimation; direct calibration with spectroscopic redshifts from VIPERS; and clustering redshift using SDSS galaxies. The fiducial redshift distribution is a combination of these three techniques performed by modifying the mean and width of the DNF distributions to match those of VIPERS and clustering redshift. In this paper we also describe the methodology used to mitigate the effect of observational systematics, which is analogous to the one used in the Y3 analysis. This paper is one of the two dedicated to the analysis of the BAO signal in DES Y6. In its companion paper, we present the angular diameter distance constraints obtained through the fitting to the BAO scale.
We describe a family of twisted partition functions for the relativistic spinning particle models. For suitable choices of fugacities this computes a refined Euler characteristics that counts the dimension of the physical states for arbitrary picture and, furthermore, encodes the complete BV-spectrum of the effective space-time gauge theory originating from this model upon second quantization. The relation between twisted world-line partition functions and the spectrum of the space-time theory is most easily seen on-shell but we will give an off-shell description as well. Finally we discuss the construction of a space-time action in terms of the world-line fields in analogy to string field theory.
The orbital distribution of the S-star cluster surrounding the supermassive black hole in the center of the Milky Way is analyzed. A tight dependence of the pericenter distance r p on orbital eccentricity e ⋆ is found, $\mathrm{log}({r}_{{\rm{p}}})\sim (1-{e}_{\star })$ , which cannot be explained simply by a random distribution of semimajor axis and eccentricities. No stars are found in the region with high e ⋆ and large $\mathrm{log}({r}_{{\rm{p}}})$ or in the region with low e ⋆ and small $\mathrm{log}({r}_{{\rm{p}}})$ . Although the sample is still small, the G-clouds show a very similar distribution. The likelihood $P(\mathrm{log}({r}_{{\rm{p}}}),(1-{e}_{\star }))$ to determine the orbital parameters of S-stars is determined. P is very small for stars with large e ⋆ and large $\mathrm{log}({r}_{{\rm{p}}})$ . S-stars might exist in this region. To determine their orbital parameters, one however needs observations over a longer time period. On the other hand, if stars would exist in the region of low $\mathrm{log}({r}_{{\rm{p}}})$ and small e ⋆, their orbital parameters should by now have been determined. That this region is unpopulated therefore indicates that no S-stars exist with these orbital characteristics, providing constraints for their formation. We call this region, defined by $\mathrm{log}({r}_{{\rm{p}}}/\mathrm{AU})\lt 1.57+2.6(1-{e}_{\star })$ , the zone of avoidance. Finally, it is shown that the observed frequency of eccentricities and pericenter distances is consistent with a random sampling of $\mathrm{log}({r}_{{\rm{p}}})$ and e ⋆ if one takes into account the fact that no stars exist in the zone of avoidance and that orbital parameters cannot yet be determined for stars with large r p and large e ⋆.
Nuclear double-beta decays are a unique probe to search for new physics beyond the standard model. Hypothesized particles, non-standard interactions, or the violation of fundamental symmetries would affect the decay kinematics, creating detectable and characteristic experimental signatures. In particular, the energy distribution of the electrons emitted in the decay gives an insight into the decay mechanism and has been studied in several isotopes and experiments. No deviations from the prediction of the standard model have been reported yet. However, several new experiments are underway or in preparation and will soon increase the sensitivity of these beyond-the-standard-model physics searches, exploring uncharted parts of the parameter space. This review brings together phenomenological and experimental aspects related to new-physics searches in double-beta decay experiments, focusing on the testable models, the most-sensitive detection techniques, and the discovery opportunities of this field.
The holographic principle suggests that regions of space contain fewer physical degrees of freedom than would be implied by conventional quantum field theory. Meanwhile, in Hilbert spaces of large dimension $2^n$, it is possible to define $N \gg n$ Pauli algebras that are nearly anti-commuting (but not quite) and which can be thought of as "overlapping degrees of freedom". We propose to model the phenomenology of holographic theories by allowing field-theory modes to be overlapping, and derive potential observational consequences. In particular, we build a Fermionic quantum field whose effective degrees of freedom approximately obey area scaling and satisfy a cosmic Bekenstein bound, and compare predictions of that model to cosmic neutrino observations. Our implementation of holography implies a finite lifetime of plane waves, which depends on the overall UV cutoff of the theory. To allow for neutrino flux from blazar TXS 0506+056 to be observable, our model needs to have a cutoff $k_{\mathrm{UV}} \lesssim 500\, k_{\mathrm{LHC}}\,$. This is broadly consistent with current bounds on the energy spectrum of cosmic neutrinos from IceCube, but high energy neutrinos are a potential challenge for our model of holography. We motivate our construction via quantum mereology, \ie using the idea that EFT degrees of freedom should emerge from an abstract theory of quantum gravity by finding quasi-classical Hilbert space decompositions. We also discuss how to extend the framework to Bosons. Finally, using results from random matrix theory we derive an analytical understanding of the energy spectrum of our theory. The numerical tools used in this work are publicly available within the GPUniverse package, github.com/OliverFHD/GPUniverse .
The origin of obscuration in active galactic nuclei (AGNs) is still an open debate. In particular, it is unclear what drives the relative contributions to the line-of-sight column densities from galaxy-scale and torus-linked obscuration. The latter source is expected to play a significant role in Unification Models, while the former is thought to be relevant in both Unification and Evolutionary models. In this work, we make use of a combination of cosmological semi-analytic models and semi-empirical prescriptions for the properties of galaxies and AGN, to study AGN obscuration. We consider a detailed object-by-object modelling of AGN evolution, including different AGN light curves (LCs), gas density profiles, and also AGN feedback-induced gas cavities. Irrespective of our assumptions on specific AGN LC or galaxy gas fractions, we find that, on the strict assumption of an exponential profile for the gas component, galaxy-scale obscuration alone can hardly reproduce the fraction of log (NH/cm-2) ≥ 24 sources at least at z ≲ 3. This requires an additional torus component with a thickness that decreases with luminosity to match the data. The torus should be present in all evolutionary stages of a visible AGN to be effective, although galaxy-scale gas obscuration may be sufficient to reproduce the obscured fraction with 22 < log (NH/cm-2) < 24 (Compton-thin, CTN) if we assume extremely compact gas disc components. The claimed drop of CTN fractions with increasing luminosity does not appear to be a consequence of AGN feedback, but rather of gas reservoirs becoming more compact with decreasing stellar mass.
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 study the 10 Myr evolution of parsec-scale stellar discs with initial masses of Mdisc = 1.0-$7.5 \times 10^4\, \mathrm{M}_\odot$ and eccentricities einit = 0.1-0.9 around supermassive black holes (SMBHs). Our disc models are embedded in a spherical background potential and have top-heavy single and binary star initial mass functions (IMF) with slopes of 0.25-1.7. The systems are evolved with the N-body code BIFROST, including post-Newtonian (PN) equations of motion and simplified stellar evolution. All discs are unstable and evolve on Myr time-scales towards similar eccentricity distributions peaking at e⋆ ~ 0.3-0.4. Models with high einit also develop a very eccentric (e⋆ ≳ 0.9) stellar population. For higher disc masses Mdisc ≳ 3 × 104 M⊙, the disc disruption dynamics is more complex than the standard secular eccentric disc instability with opposite precession directions at different disc radii - a precession direction instability. We present an analytical model describing this behaviour. A milliparsec population of N ~ 10-100 stars forms around the SMBH in all models. For low einit, stars migrate inward while for einit ≳ 0.6 stars are captured by the Hills mechanism. Without PN, after 6 Myr, the captured stars have a sub-thermal eccentricity distribution. We show that including PN effects prevents this thermalization by suppressing resonant relaxation effects and cannot be ignored. The number of tidally disrupted stars is similar or larger than the number of milliparsec stars. None of the simulated models can simultaneously reproduce the kinematic and stellar population properties of the Milky Way centre clockwise disc and the S-cluster.
Disc winds and planet-disc interactions are two crucial mechanisms that define the structure, evolution and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. We study the interplay between an X-ray photoevaporative disc wind and the substructures generated due to planet-disc interactions to determine how their mutual interactions affect the disc's and the planet's evolution. We perform three-dimensional hydrodynamic simulations of viscous ($\alpha = 6.9\cdot10^{-4}$) discs that host a Jupiter-like planet and undergo X-ray photoevaporation. We trace the gas flows within the disc and wind and measure the accretion rate onto the planet, as well as the gravitational torque that is acting on it. Our results show that the planetary gap takes away the wind's pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower, and the flow of mass across the gap in both directions is significantly increased, as well as the planet's mass-accretion rate (by factors $\approx 5$ and $\approx 2$, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes of the inner disc, contrary to the expectation in a planet-induced photoevaporation (PIPE) scenario that has been proposed in the past.
High-velocity stellar collisions driven by a supermassive black hole (BH) or BH-driven disruptive collisions in dense, nuclear clusters can rival the energetics of supergiant star explosions following the gravitational collapse of their iron core. Starting from a sample of red-giant star collisions simulated with the hydrodynamics code AREPO, we generated photometric and spectroscopic observables using the nonlocal thermodynamic equilibrium time-dependent radiative transfer code CMFGEN. Collisions from more extended giants or more violent collisions (with higher velocities or smaller impact parameters) yield bolometric luminosities on the order of 1043 erg s−1 at 1 d, evolving on a timescale of a week to a bright plateau at ∼1041 erg s−1 before plunging precipitously after 20-40 d at the end of the optically thick phase. This luminosity falls primarily in the UV in the first few days, thus when it is at its maximum, and shifts to the optical thereafter. Collisions at lower velocities or from less extended stars produce ejecta that are fainter but can remain optically thick for up to 40 d if they have a low expansion rate. This collision debris shows a similar spectral evolution as that observed or modeled for Type II supernovae from blue-supergiant star explosions, differing only in the more rapid transition to the nebular phase. Such BH-driven disruptive collisions should be detectable by high-cadence surveys in the UV such as ULTRASAT.
Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2". Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher-resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 days and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ~2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (> 1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually nonexistent outside of interacting binaries.
Dark photons, aside from constituting non-relativistic dark matter, can also be generated relativistically through the decay or annihilation of other dark matter candidates, contributing to a galactic dark photon background. The production of dark photons tends to favor specific polarization modes, determined by the microscopic coupling between dark matter and dark photons. We leverage data obtained from previous searches for dark photon dark matter using a superconducting radio-frequency cavity to explore galactic dark photon fluxes. The interplay of anisotropic directions and Earth's rotation introduces a diurnal modulation of signals within the cavities, manifesting distinct variation patterns for longitudinal and transverse modes. Our findings highlight the efficacy of superconducting radio-frequency cavities, characterized by significantly high-quality factors, as powerful telescopes for detecting galactic dark photons, unveiling a novel avenue in the indirect search for dark matter through multi-messenger astronomy.
The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4) line emission at the spatial resolution of 50 au. We detect extended SO emission, peaking towards the A and B protostars. Emission blue-shifted down to 6.6 km s$^{-1}$ reveals for the first time a long ($\sim$ 2000 au) accelerating streamer plausibly feeding the VLA1623 B protostar. Using SO, we derive for the first time an estimate of the excitation temperature of an accreting streamer: 33$\pm$9 K. The SO column density is $\sim$ 10$^{14}$ cm$^{-2}$, and the SO/H$_2$ abundance ratio is $\sim$ 10$^{-8}$. The total mass of the streamer is 3 $\times$ 10$^{-3}$ $Msun$, while its accretion rate is 3--5 $\times$ 10$^{-7}$ Msun yr$^{-1}$. This is close to the mass accretion rate of VLA1623 B, in the 0.6--3 $\times$ 10$^{-7}$ Msun yr$^{-1}$ range, showing the importance of the streamer in contributing to the mass of protostellar disks. The highest blue- and red-shifted SO velocities behave as the SiO(5--4) emission, the latter species detected for the first time in VLA1623-2417: the emission is compact (100-200 au), and associated only with the B protostar. The SO excitation temperature is $\sim$ 100 K, supporting the occurrence of shocks associated with the jet, traced by SiO.
The general epidemic process (GEP), also known as susceptible-infected-recovered model (SIR), describes how an epidemic spreads within a population of susceptible individuals who acquire permanent immunization upon recovery. This model exhibits a second-order absorbing state phase transition, commonly studied assuming immobile healthy individuals. We investigate the impact of mobility on disease spreading near the extinction threshold by introducing two generalizations of GEP, where the mobility of susceptible and recovered individuals is examined independently. In both cases, including mobility violates GEP's rapidity reversal symmetry and alters the number of absorbing states. The critical dynamics of the models are analyzed through a perturbative renormalization group approach and large-scale stochastic simulations using a Gillespie algorithm. The renormalization group analysis predicts both models to belong to the same novel universality class describing the critical dynamics of epidemic spreading when the infected individuals interact with a diffusive species and gain immunization upon recovery. At the associated renormalization group fixed point, the immobile species decouples from the dynamics of the infected species, dominated by the coupling with the diffusive species. Numerical simulations in two dimensions affirm our renormalization group results by identifying the same set of critical exponents for both models. Violation of the rapidity reversal symmetry is confirmed by breaking the associated hyperscaling relation. Our study underscores the significance of mobility in shaping population spreading dynamics near the extinction threshold.
Aims: The BL Lac 1ES 2344+514 is known for temporary extreme properties characterised by a shift of the synchrotron spectral energy distribution (SED) peak energy νsynch, p above 1 keV. While those extreme states have only been observed during high flux levels thus far, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to very high energy (VHE) performed so far, focussing on a systematic characterisation of the intermittent extreme states.
Methods: We organised a monitoring campaign covering a 3-year period from 2019 to 2021. More than ten instruments participated in the observations in order to cover the emission from radio to VHE. In particular, sensitive X-ray measurements by XMM-Newton, NuSTAR, and AstroSat took place simultaneously with multi-hour MAGIC observations, providing an unprecedented constraint of the two SED components for this blazar.
Results: While our results confirm that 1ES 2344+514 typically exhibits νsynch, p > 1 keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase in the electron acceleration efficiency without a change in the electron injection luminosity. On the other hand, we also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of νsynch, p. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. By combining Swift-XRT and Swift-UVOT measurements during a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions significantly contribute to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3−2 keV band. Using time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model.
Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star.
Aims: We investigated here how strong MB torques need to be to reproduce the observationally inferred relative numbers of white dwarf plus M dwarf post-common-envelope binaries under the assumption of magnetic saturation.
Methods: We carried out binary population simulations with the BSE code adopting empirically derived inter-correlated main-sequence binary distributions as initial binary populations and compared the simulation outcomes with observations.
Results: We found that the dearth of extreme mass ratio binaries in the inter-correlated initial distributions is key to reproduce the large fraction of post-common-envelope binaries hosting low-mass M dwarfs (∼0.1 − 0.2 M⊙). In addition, orbital angular momentum loss rates due to MB should be high for M dwarfs with radiative cores and orders of magnitude smaller for fully convective stars to explain the observed dramatic change of the fraction of short-period binaries at the fully convective boundary.
Conclusions: We conclude that saturated but disrupted, that is, dropping drastically at the fully convective boundary, MB can explain the observations of both close main-sequence binaries containing M and K dwarfs and post-common-envelope binaries. Whether a similar prescription can explain the spin down rates of single stars and of binaries containing more massive stars needs to be tested.
Massive ($\gtrsim9$ M$_\odot$) and very massive ($\gtrsim100$ M$_\odot$) stars are expected to release large amounts of energy and metal-enriched material throughout their relatively short lives. Their stellar winds may play an important role in the metal-enrichment during the formation of star clusters. With novel high-resolution hydrodynamical GRIFFIN-project simulations, we investigate the rapid recycling of stellar wind-material during the formation of massive star clusters up to $M_\mathrm{cluster}\sim2\times10^5$ M$_\odot$ in a low-metallicity dwarf galaxy starburst. The simulation realises new stars from a stellar initial mass function (IMF) between $0.08$ M$_\odot$ and $\sim400$ M$_\odot$and follows stellar winds, radiation and supernova-feedback of single massive stars with evolution tracks. Star clusters form on timescales less than $\sim5$ Myr, and their supernova-material is very inefficiently recycled. Stellar wind-material, however, is trapped in massive clusters resulting in the formation of stars self-enriched in Na, Al, and N within only a few Myr. Wind-enriched (second population) stars can be centrally concentrated in the most massive clusters ($\gtrsim10^4$ M$_\odot$) and the locked wind-material increases approximately as $M_\mathrm{cluster}^{2}$. These trends resemble the characteristics of observed second population stars in globular clusters. We fit scaling relations to the log-normal distributed wind-mass fractions and extrapolate to possible globular cluster progenitors of $M_\mathrm{cluster}=10^7$ M$_\odot$ to investigate whether a dominant second population could form. This can only happen if the IMF is well sampled, single massive stars produce at least a factor of a few more enriched winds e.g. through a top-heavy IMF, and a significant fraction of the first population (unenriched) stars is lost during cluster evolution.
Extragalactic and galactic cosmic rays scatter with the cosmic neutrino background during propagation to Earth, yielding a flux of relic neutrinos boosted to larger energies. If an overdensity of relic neutrinos is present in galaxies, and neutrinos are massive enough, this flux might be detectable by high-energy neutrino experiments. For a lightest neutrino of mass $m_{\nu} \sim 0.1$ eV, we find an upper limit on the local relic neutrino overdensity of $\sim 10^{13}$ and an upper limit on the relic neutrino overdensity at TXS 0506+056 of $\sim 10^{10}$. Future experiments like GRAND or IceCube-Gen2 could improve these bounds by orders of magnitude.
Our Sun lies within 300 pc of the 2.7-kpc-long sinusoidal chain of dense gas clouds known as the Radcliffe Wave. The structure's wave-like shape was discovered using 3D dust mapping, but initial kinematic searches for oscillatory motion were inconclusive. Here we present evidence that the Radcliffe Wave is oscillating through the Galactic plane while also drifting radially away from the Galactic Center. We use measurements of line-of-sight velocity for 12CO and 3D velocities of young stellar clusters to show that the most massive star-forming regions spatially associated with the Radcliffe Wave (including Orion, Cepheus, North America, and Cygnus X) move as if they are part of an oscillating wave driven by the gravitational acceleration of the Galactic potential. By treating the Radcliffe Wave as a coherently oscillating structure, we can derive its motion independently of the local Galactic mass distribution, and directly measure local properties of the Galactic potential as well as the Sun's vertical oscillation period. In addition, the measured drift of the Radcliffe Wave radially outward from the Galactic Center suggests that the cluster whose supernovae ultimately created today's expanding Local Bubble may have been born in the Radcliffe Wave.
The eROSITA telescope array aboard the Spektrum Roentgen Gamma (SRG) satellite began surveying the sky in December 2019, with the aim of producing all-sky X-ray source lists and sky maps of an unprecedented depth. Here we present catalogues of both point-like and extended sources using the data acquired in the first six months of survey operations (eRASS1; completed June 2020) over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process, the data analysis pipelines, and the characteristics of the X-ray sources. With nearly 930 000 entries detected in the most sensitive 0.2-2.3 keV energy range, the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%, and provides a comprehensive inventory of all classes of X-ray celestial objects, covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.3-5 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We present methods to identify and flag potential spurious sources in the catalogues, which we applied for this work, and we tested and validated the astrometric accuracy via cross-comparison with other X-ray and multi-wavelength catalogues. We show that the number counts of X-ray sources in eRASSl are consistent with those derived over narrower fields by past X-ray surveys of a similar depth, and we explore the number counts variation as a function of the location in the sky. Adopting a uniform all-sky flux limit (at 50% completeness) of F05-2 keV > 5 × 10−14 erg s−1 cm−2, we estimate that the eROSITA all-sky survey resolves into individual sources about 20% of the cosmic X-ray background in the 1-2 keV range. The catalogues presented here form part of the first data release (DR1) of the SRG/eROSITA all-sky survey. Beyond the X-ray catalogues, DR1 contains all detected and calibrated event files, source products (light curves and spectra), and all-sky maps. Illustrative examples of these are provided.
The catalogue is 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/682/A34
For decades, the boundary of cosmic filaments have 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 z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience 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.
Recent observations of galaxy mergers inside galaxy cluster environments report high star formation rates in the ejected tidal tails, which point towards currently developing tidal dwarf galaxies. We test whether these dwarf objects could get stripped from the galaxy potential by the galaxy cluster and thus populate it with dwarf galaxies. To this end, we perform three high-resolution hydrodynamical simulations of mergers between spiral galaxies in a cluster environment, varying the initial orbit of the infalling galaxies with respect to the cluster center. We demonstrate that cluster environments are indeed capable of stripping tidal dwarf galaxies from the host potential in all tested setups. In the three orbit scenarios, we find 3, 7, and 8 tidal dwarf galaxies per merger, respectively, which survive longer than 1 Gyr after the merger event. Exposed to ram pressure, these gas dominated dwarf galaxies exhibit high star formation rates while also losing gas to the environment. Experiencing a strong headwind due to their motion through the intracluster medium, they quickly lose momentum and start spiraling towards the cluster center, reaching distances on the order of ~Mpc from their progenitor. About 4 Gyr after the merger event, we still find several intact dwarf galaxies, demonstrating that such objects can prevail for a significant fraction of the Hubble time. Comparing their contribution to the observed galaxy mass function in clusters, our results indicate that ~30% of dwarf galaxies in clusters could have been formed by stripping from galaxy mergers.
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 half-decay time. 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 BBN and CMB 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^9$g which enter the memory burden stage before BBN and are still present today as DM.
The phenomenon of multiple stellar populations is exacerbated in massive globular clusters, with the fraction of first-population (1P) stars a decreasing function of the cluster present-day mass. We decipher this relation in far greater detail than has been done so far. We assume (i) a fixed stellar mass threshold for the formation of second-population (2P) stars, (ii) a power-law scaling $F_{1P} \propto m_{ecl}^{-1}$ between the mass $m_{ecl}$ of newly-formed clusters and their 1P-star fraction $F_{1P}$, and (iii) a constant $F_{1P}$ over time. The $F_{1P}(m_{ecl})$ relation is then evolved up to an age of 12Gyr for tidal field strengths representative of the entire Galactic halo. The 12Gyr-old model tracks cover extremely well the present-day distribution of Galactic globular clusters in (mass,$F_{1P}$) space. The distribution is curtailed on its top-right side by the scarcity of clusters at large Galactocentric distances, and on its bottom-left side by the initial scarcity of very high-mass clusters, and dynamical friction. Given their distinct dissolution rates, "inner" and "outer" model clusters are offset from each other, as observed. The locus of Magellanic Clouds clusters in (mass,$F_{1P}$) space is as expected for intermediate-age clusters evolving in a gentle tidal field. Given the assumed constancy of $F_{1P}$, we conclude that 2P-stars do not necessarily form centrally-concentrated. We infer a minimum mass of $4 \cdot 10^5~M_{\odot}$ for multiple-populations clusters at secular evolution onset. This high-mass threshold severely limits the amount of 2P-stars lost from evolving clusters, thereby fitting the low 2P-star fraction of the Galactic halo field.
This is the second paper in a series presenting the results from a 500 $h^{-1}$Mpc large constrained hydro-dynamical simulation of the local Universe (SLOW). The initial conditions are based on peculiar velocities derived from the CosmicFlows-2 catalogue. The inclusion of galaxy formation treatment, allows to directly predict observable properties of the Intra-Cluster Medium (ICM) within galaxy clusters. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts, enables us to assess the effectiveness of the initial condition constraints in accurately replicating the non-linear properties of the largest, collapsed objects within the simulation. Based on the combination of several, publicly available surveys, we identified 45 local Universe galaxy clusters in SLOW, including the 13 most massive from the Planck SZ catalog and 70% of those with $M_{500} > 2\times 10^{14}$ M$_{\odot}$. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature and Compton-y by comparing it to a random selection. In relation to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross identified sample allows to independently deduce a hydrostatic mass bias of $(1-b)\approx0.87$. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of the galaxy clusters within our local neighbourhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM.
We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combines time-delay cosmography in lens galaxy clusters with pure samples of `cosmic chronometers' by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits the most massive and passive cluster member galaxies to measure the differential time evolution of the Universe. We applied two different statistical techniques, adopting realistic errors on the measured quantities, to assess the accuracy and the gain in precision on the values of the cosmological parameters. We demonstrate that the proposed combined method allows for a robust and accurate measurement of the value of the Hubble constant. In addition, this provides valuable information on the other cosmological parameters thanks to the complementarity between the two different probes in breaking parameter degeneracies. Finally, we showcased the immediate observational feasibility of the proposed joint method by taking advantage of the existing high-quality spectro-photometric data for several lens galaxy clusters.
We present high-definition observations with the James Webb Space Telescope (JWST) of >1000 Cepheids in a geometric anchor of the distance ladder, NGC 4258, and in five hosts of eight Type Ia supernovae, a far greater sample than previous studies with JWST. These galaxies individually contain the largest samples of Cepheids, an average of >150 each, producing the strongest statistical comparison to those previously measured with the Hubble Space Telescope (HST) in the near-infrared (NIR). They also span the distance range of those used to determine the Hubble constant with HST, allowing us to search for a distance-dependent bias in HST measurements. The superior resolution of JWST negates crowding noise, the largest source of variance in the NIR Cepheid period–luminosity relations (Leavitt laws) measured with HST. Together with the use of two epochs to constrain Cepheid phases and three filters to remove reddening, we reduce the dispersion in the Cepheid P–L relations by a factor of 2.5. We find no significant difference in the mean distance measurements determined from HST and JWST, with a formal difference of ‑0.01 ± 0.03 mag. This result is independent of zero-points and analysis variants including metallicity dependence, local crowding, choice of filters, and slope of the relations. We can reject the hypothesis of unrecognized crowding of Cepheid photometry from HST that grows with distance as the cause of the "Hubble tension" at 8.2σ, i.e., greater confidence than that of the Hubble tension itself. We conclude that errors in photometric measurements of Cepheids across the distance ladder do not significantly contribute to the tension.
We investigate the main tensions within the current standard model of cosmology from the perspective of the main statistics of cosmic voids, using the final BOSS DR12 data set. For this purpose, we present the first estimate of the S8 ≡ σ8 Ωm/0.3 and H0 parameters obtained from void number counts and shape distortions. To analyze void counts we relied on an extension of the popular volume-conserving model for the void size function, tailored to the application on data, including geometric and dynamic distortions. We calibrated the two nuisance parameters of this model with the official BOSS Collaboration mock catalogs and propagated their uncertainty through the statistical analysis of the BOSS void number counts. The constraints from void shapes come from the study of the geometric distortions of the stacked void-galaxy cross-correlation function. In this work we focus our analysis on the Ωm − σ8 and Ωm − H0 parameter planes and derive the marginalized constraints S8 = 0.813−0.068+0.093 and H0 = 67.3−9.1+10.0 km s−1 Mpc−1, which are fully compatible with constraints from the literature. These results are expected to notably improve in precision when analyzed jointly with independent probes and will open a new viewing angle on the rising cosmological tensions in the near future.
The deposition of energy and momentum by supernova explosions has been subject to numerous studies in the past few decades. However, while there has been some work that focused on the transition from the adiabatic to the radiative stage of a supernova remnant (SNR), the late radiative stage and merging with the interstellar medium (ISM) have received little attention. Here, we use three-dimensional, hydrodynamic simulations, focusing on the evolution of SNRs during the radiative phase, considering a wide range of physical explosion parameters ($n_{\text{H, ISM}} \in \left[0.1, 100\right] \text{cm}^{-3}$ and $E_{\text{SN}} \in \left[1, 14\right]\times 10^{51} \text{erg}$). We find that the radiative phase can be subdivided in four stages: A pressure driven snowplow phase during which the hot overpressurized bubble gas is evacuated and pushed into the cold shell, a momentum conserving snowplow phase which is accompanied by a broadening of the shell, an implosion phase where cold material from the back of the shell is flooding the central vacuum and a final cloud phase, during which the imploding gas is settling as a central, compact overdensity. The launching timescale for the implosion ranges from a few 100 kyr to a few Myr, while the cloud formation timescale ranges from a few to about 10 Myr. The highly chemically enriched clouds can become massive ($M_{\text{cl}} \sim 10^3 - 10^4 \, \text{M}_{\odot}$) and self-gravitating within a few Myr after their formation, providing an attractive, novel pathway for supernova induced star and planet formation in the ISM.
We present a comprehensive catalog of 2824 RR Lyrae stars (RRLs) residing in 115 Galactic globular clusters (GCs). Our catalog includes 1594 fundamental-mode (RRab), 824 first-overtone (RRc), and 28 double-mode (RRd) RRLs, as well as 378 RRLs of an unknown pulsation mode. We cross-matched 481349 RRLs reported in the third data release (DR3) of the ESA mission Gaia and the literature to 170 known GCs. Membership probabilities were computed as the products of a position and shape-dependent prior and a likelihood was computed using parallaxes, proper motions, and, where available, radial velocities from Gaia. Membership likelihoods of RRLs were computed by comparing cluster average parameters based on known member stars and the cross-matched RRLs. We determined empirical RRL instability strip (IS) boundaries based on our catalog and detected three new cluster RRLs inside this region via their excess Gaia G-band photometric uncertainties. We find that 77% of RRLs in GCs are included in the Gaia DR3 Specific Object Study, and 82% were classified as RRLs by the Gaia DR3 classifier, with the majority of the missing sources being located at the crowded GC centers. Surprisingly, we find that 25% of cluster member stars located within the empirical IS are not RRLs and appear to be non-variable. Additionally, we find that 80% of RRab, 84% of RRc, and 100% of the RRd stars are located within theoretical IS boundaries predicted using MESA models with Z = 0.0003, M = 0.7 (M_\odot), and Y = 0.290. Unexpectedly, a higher Y = 0.357 is required to fully match the location of RRc stars, and lower Y = 0.220 is needed to match the location of RRab stars. Lastly, our catalog does not exhibit an Oosterhoff dichotomy, with at least 22 GCs located inside the Oosterhoff "gap," which is close to the mode of the distribution of mean RRL periods in GCs.
To shed light on the origin of magnetic symbiotic stars, we investigated the system FN Sgr in detail. We searched for a reasonable formation pathway to explain its stellar and binary parameters including the magnetic field of the accreting white dwarf. We used the MESA code to carry out pre-CE and post-CE binary evolution and determined the outcome of CE evolution assuming the energy formalism. For the origin and evolution of the white dwarf magnetic field, we adopted the crystallization scenario. We found that FN Sgr can be explained as follows. First, a non-magnetic white dwarf is formed through CE evolution. Later, during post-CE evolution, the white dwarf starts to crystallize and a weak magnetic field is generated. After a few hundred Myr, the magnetic field penetrates the white dwarf surface and becomes detectable. Meanwhile, its companion evolves and becomes an evolved red giant. Subsequently, the white dwarf accretes part of the angular momentum from the red giant stellar winds. As a result, the white dwarf spin period decreases and its magnetic field reaches super-equipartition, getting amplified due to a rotation- and crystallization-driven dynamo. The binary then evolves into a symbiotic star, with a magnetic white dwarf accreting from an evolved red giant through atmospheric Roche-lobe overflow. We conclude that the rotation- and crystallization-driven dynamo scenario, or any age-dependent scenario, can explain the origin of magnetic symbiotic stars reasonably well. This adds another piece to the pile of evidence supporting this scenario. If our formation channel is correct, our findings suggest that white dwarfs in most symbiotic stars formed through CE evolution might be magnetic, provided that the red giant has spent >3 Gyr as a main-sequence star.
The paradigm-changing possibility of collective neutrino-antineutrino oscillations was recently advanced in analogy to collective flavor oscillations. However, the amplitude for the backward scattering process $\nu_{\mathbf{p}_1}\overline\nu_{\mathbf{p}_2}\to\nu_{\mathbf{p}_2}\overline\nu_{\mathbf{p}_1}$ is helicity-suppressed and vanishes for massless neutrinos, implying that there is no off-diagonal refractive index between $\nu$ and $\overline\nu$ of a single flavor of massless neutrinos. For a nonvanishing mass, collective helicity oscillations are possible, representing de-facto $\nu$--$\overline\nu$ oscillations in the Majorana case. However, such phenomena are suppressed by the smallness of neutrino masses as discussed in the previous literature.
We develop a neural network based pipeline to estimate masses of galaxy clusters with a known redshift directly from photon information in X-rays. Our neural networks are trained using supervised learning on simulations of eROSITA observations, focusing in this paper on the Final Equatorial Depth Survey (eFEDS). We use convolutional neural networks which are modified to include additional information of the cluster, in particular its redshift. In contrast to existing work, we utilize simulations including background and point sources to develop a tool which is usable directly on observational eROSITA data for an extended mass range from group size halos to massive clusters with masses in between 1013M⊙<M<1015M⊙. Using this method, we are able to provide for the first time neural network mass estimation for the observed eFEDS cluster sample from Spectrum-Roentgen-Gamma/eROSITA observations and we find consistent performance with weak lensing calibrated masses. In this measurement, we do not use weak lensing information and we only use previous cluster mass information which was used to calibrate the cluster properties in the simulations. When compared to simulated data, we observe a reduced scatter with respect to luminosity and count-rate based scaling relations.
We comment on the application for other upcoming eROSITA All-Sky Survey observations.
The last few years have seen the development of a promising theoretical framework for statistics of the cosmic large-scale structure -- the theory of large deviations (LDT) for modelling weak-lensing one-point statistics in the mildly non-linear regime. The goal of this series of papers is to make the leap and lay out the steps to perform an actual data analysis with this theoretical tool. Building upon the LDT framework, in this work (Paper I) we demonstrate how to accurately model the Probability Distribution Function (PDF) of a reconstructed Kaiser-Squires convergence field under a realistic mask, that of the third data release of the Dark Energy Survey (DES). We also present how weak lensing systematics and higher-order lensing corrections due to intrinsic alignments, shear biases, photo-z errors and baryonic feedback can be incorporated in the modelling of the reconstructed convergence PDF. In an upcoming work (Paper II) we will then demonstrate the robustness of our modelling through simulated likelihood analyses, the final step required before applying our method to actual data.
In this series of papers, we present an emulator-based halo model for the non-linear clustering of galaxies in modified gravity cosmologies. In the first paper, we present emulators for the following halo properties: the halo mass function, concentration-mass relation and halo-matter cross-correlation function. The emulators are trained on data extracted from the FORGE and BRIDGE suites of N-body simulations, respectively, for two modified gravity (MG) theories: f(R) gravity, and the DGP model, varying three standard cosmological parameters Ωm0, H0, σ8, and one MG parameter, either $\bar{f}_{R0}$ or rc. Our halo property emulators achieve an accuracy of ${\lesssim}1\ \hbox{per cent}$ on independent test data sets. We demonstrate that the emulators can be combined with a galaxy-halo connection prescription to accurately predict the galaxy-galaxy and galaxy-matter correlation functions using the halo model framework.
Accurate redshift calibration is required to obtain unbiased cosmological information from large-scale galaxy surveys. In a forward modelling approach, the redshift distribution n(z) of a galaxy sample is measured using a parametric galaxy population model constrained by observations. We use a model that captures the redshift evolution of the galaxy luminosity functions, colours, and morphology, for red and blue samples. We constrain this model via simulation-based inference, using factorized Approximate Bayesian Computation (ABC) at the image level. We apply this framework to HSC deep field images, complemented with photometric redshifts from COSMOS2020. The simulated telescope images include realistic observational and instrumental effects. By applying the same processing and selection to real data and simulations, we obtain a sample of n(z) distributions from the ABC posterior. The photometric properties of the simulated galaxies are in good agreement with those from the real data, including magnitude, colour and redshift joint distributions. We compare the posterior n(z) from our simulations to the COSMOS2020 redshift distributions obtained via template fitting photometric data spanning the wavelength range from UV to IR. We mitigate sample variance in COSMOS by applying a reweighting technique. We thus obtain a good agreement between the simulated and observed redshift distributions, with a difference in the mean at the 1σ level up to a magnitude of 24 in the i band. We discuss how our forward model can be applied to current and future surveys and be further extended. The ABC posterior and further material will be made publicly available at this https URL.
We present a new high-precision strong lensing model of PLCK G287.0+32.9, a massive lens galaxy cluster at z=0.383, with the aim to get an accurate estimation of its effective Einstein radius and total mass distribution. We also present a spectroscopic catalog containing accurate redshift measurements for 490 objects, including multiply-lensed sources and cluster member galaxies. We exploit high-quality spectroscopic data from three pointings of the VLT Multi Unit Spectroscopic Explorer, covering a central 3 arcmin2 region of the cluster. We complete the spectroscopic catalog by including redshift measurements from VLT-VIMOS and KECK-DEIMOS. We identify 129 spectroscopic cluster member galaxies, with redshift values 0.360≤z≤0.405 and mF160W≤21, and 24 photometric ones identified with a Convolutional Neural Network from ancillary HST imaging. We also identify 114 multiple images from 28 background sources, of which 84 images from 16 sources are new and the remaining ones were identified in previous work. The best-fitting lens model shows a root mean square separation value between the predicted and observed positions of the multiple images of 0.75″, corresponding to an improvement in reconstructing the observed positions of the multiple images of a factor of 2.5 with respect to previous models. Using the predictive power of our new lens model we find 3 new multiple images and we confirm the configuration of three systems of multiple images that were not used for the optimization of the model. The derived total mass distribution confirms this cluster to be a very prominent gravitational lens with an effective Einstein θE=43.4″±0.1″, that is in agreement with previous estimates and corresponds to a total mass enclosed in the critical curve of ME=3.33+0.02−0.07×1014M⊙.
High-resolution 3D maps of interstellar dust are critical for probing the underlying physics shaping the structure of the interstellar medium, and for foreground correction of astrophysical observations affected by dust. We aim to construct a new 3D map of the spatial distribution of interstellar dust extinction out to a distance of 1.25 kpc from the Sun. We leveraged distance and extinction estimates to 54 million nearby stars derived from the Gaia BP/RP spectra. Using the stellar distance and extinction information, we inferred the spatial distribution of dust extinction. We modeled the logarithmic dust extinction with a Gaussian process in a spherical coordinate system via iterative charted refinement and a correlation kernel inferred in previous work. In total, our posterior has over 661 million degrees of freedom. We probed the posterior distribution using the variational inference method MGVI. Our 3D dust map has an angular resolution of up to 14' (Nside = 256), and we achieve parsec-scale distance resolution, sampling the dust in 516 logarithmically spaced distance bins spanning 69 pc to 1250 pc. We generated 12 samples from the variational posterior of the 3D dust distribution and release the samples alongside the mean 3D dust map and its corresponding uncertainty. Our map resolves the internal structure of hundreds of molecular clouds in the solar neighborhood and will be broadly useful for studies of star formation, Galactic structure, and young stellar populations. It is available for download in a variety of coordinate systems online and can also be queried via the publicly available dustmaps Python package.
The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.
The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg2 over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg2 using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg2 program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval z<1.1, 0.39% over the redshift interval 1.1<z<1.9, and 0.46% over the redshift interval 1.9<z<3.5.
We investigate the potential of low-redshift Lyman alpha (Lyα) forest for constraining active galactic nuclei (AGN) feedback models by analyzing the Illustris and IllustrisTNG simulation at z=0.1. These simulations are ideal for studying the impact of AGN feedback on the intergalactic medium (IGM) as they share initial conditions with significant differences in the feedback prescriptions. Both simulations reveal that the IGM is significantly impacted by AGN feedback. Specifically, feedback is stronger in Illustris and results in reducing cool baryon fraction to 23% relative to 39% in IllustrisTNG. However, when comparing various statistics of Lyα forest such as 2D and marginalized distributions of Doppler widths and H I column density, line density, and flux power spectrum with real data, we find that most of these statistics are largely insensitive to the differences in feedback models. This lack of sensitivity arises because of the fundamental degeneracy between the fraction of cool baryons and the H I photoionization rate (ΓHI) as their product determines the optical depth of the Lyα forest. Since the ΓHI cannot be precisely predicted from first principles, it needs to be treated as a nuisance parameter adjusted to match the observed Lyα line density. After adjusting ΓHI, the distinctions in the considered statistics essentially fade away. Only the Lyα flux power spectrum at small spatial scales exhibits potentially observable differences, although this may be specific to the relatively extreme feedback model employed in Illustris. Without independent constraints on either ΓHI or cool baryon fraction, constraining AGN feedback with low-redshift Lyα forest will be very challenging.
Neutrinos in dense environments like core-collapse supernovae (CCSNe) and neutron star mergers (NSMs) can undergo fast flavor conversions (FFCs) once the angular distribution of neutrino lepton number crosses zero along a certain direction. Recent advancements have demonstrated the effectiveness of machine learning (ML) in detecting these crossings. In this study, we enhance prior research in two significant ways. First, we utilize realistic data from CCSN simulations, where neutrino transport is solved using the full Boltzmann equation. We evaluate the ML methods' adaptability in a real-world context, enhancing their robustness. In particular, we demonstrate that when working with artificial data, simpler models outperform their more complex counterparts, a noteworthy illustration of the bias-variance tradeoff in the context of ML. We also explore methods to improve artificial datasets for ML training. In addition, we extend our ML techniques to detect the crossings in the heavy-leptonic channels, accommodating scenarios where νx and ν¯x may differ. Our research highlights the extensive versatility and effectiveness of ML techniques, presenting an unparalleled opportunity to evaluate the occurrence of FFCs in CCSN and NSM simulations.
HD 54879 is the most recently discovered magnetic O-type star. Previous studies ruled out a rotation period shorter than 7 years, implying that HD 54879 is the second most slowly-rotating known magnetic O-type star. We report new high-resolution spectropolarimetric measurements of HD 54879, which confirm that a full stellar rotation cycle has been observed, and derive a rotation period of $P = 2562^{+63}_{-58}$ d (about 7.02 yr). The radial velocity of HD 54879 has been stable over the last decade of observations. We explore equivalent widths and longitudinal magnetic fields calculated from lines of different elements, and conclude the atmosphere of HD 54879 is likely chemically homogeneous, with no strong evidence for chemical stratification or lateral abundance nonuniformities. We present the first detailed magnetic map of the star, with an average surface magnetic field strength of 2954 G, and a strength for the dipole component of 3939 G. There is a significant amount of magnetic energy in the quadrupole components of the field (23 percent). Thus, we find HD 54879 has a strong magnetic field with a significantly complex topology.
The core-collapse of massive stars and the merger of neutron star binaries are among the most promising candidate sites for the production of high-energy cosmic neutrinos. We demonstrate that the high-energy neutrinos produced in such extreme environments can experience efficient flavor conversions on scales much shorter than those expected in vacuum, due to their coherent forward scatterings with the bath of decohered low-energy neutrinos emitted from the central engine. These low-energy neutrinos, which exist as mass eigenstates, provide a very special and peculiar dominant background for the propagation of the high-energy ones. We point out that the high-energy neutrino flavor ratio is modified to a value independent of neutrinos energies, which is distinct from the conventional prediction with the matter effect. We also suggest that the signals can be used as a novel probe of new neutrino interactions beyond the Standard Model. This is yet another context where neutrino-neutrino interactions can play a crucial role in their flavor evolution.
When QCD is described by a nonrelativistic effective field theory, operators consisting of gluonic correlators of two chromoelectric or -magnetic fields will often appear in descriptions of quarkonium physics. At zero T, these correlators give the masses of gluelumps and the moments of these correlators can be used to understand the inclusive P-wave decay of quarkonium. At finite T these correlators define the diffusion of the heavy quarkonium. However, these correlators come with a divergent term in lattice spacing which needs to be taken care of. We inspect these correlators in pure gauge theory with gradient flow smearing, which should allow us to reduce and remove the divergence more carefully. In these proceedings, we focus on the effect of gradient flow to these correlators and the reduction of this divergence.
We introduce our novel Bayesian parton density determination code, PartonDensity.jl. The motivation for this new code, the framework and its validation are described. As we show, PartonDensity.jl provides both a flexible environment for the determination of parton densities and a wealth of information concerning the knowledge update provided by the analyzed data set.
Recent observations of protostellar cores suggest that most of the material in the protostellar phase is accreted along streamers. Streamers in this context are defined as velocity coherent funnels of denser material potentially connecting the large scale environment to the small scales of the forming accretion disk. Using simulations which simultaneously resolve the driving of turbulence on the filament scale as well as the collapse of the core down to protostellar disk scales, we aim to understand the effect of the turbulent velocity field on the formation of overdensities in the accretion flow. We perform a three-dimensional numerical study on a core collapse within a turbulent filament using the RAMSES code and analyse the properties of overdensities in the accretion flow. We find that overdensities are formed naturally by the initial turbulent velocity field inherited from the filament and subsequent gravitational collimation. This leads to streams which are not really filamentary but show a sheet-like morphology. Moreover, they have the same radial infall velocities as the low density material. As a main consequence of the turbulent initial condition, the mass accretion onto the disk does not follow the predictions for solid body rotation. Instead, most of the mass is funneled by the overdensities to intermediate disk radii.
The streaming instability, a promising mechanism to drive planetesimal formation in dusty protoplanetary discs, relies on aerodynamic drag naturally induced by the background radial pressure gradient. This gradient should vary in disks, but its effect on the streaming instability has not been sufficiently explored. For this purpose, we use numerical simulations of an unstratified disc to study the non-linear saturation of the streaming instability with mono-disperse dust particles and survey a wide range of gradients for two distinct combinations of the particle stopping time and the dust-to-gas mass ratio. As the gradient increases, we find most kinematic and morphological properties increase but not always in linear proportion. The density distributions of tightly-coupled particles are insensitive to the gradient whereas marginally-coupled particles tend to concentrate by more than an order of magnitude as the gradient decreases. Moreover, dust-gas vortices for tightly-coupled particles shrink as the gradient decreases, and we note higher resolutions are required to trigger the instability in this case. In addition, we find various properties at saturation that depend on the gradient may be observable and may help reconstruct models of observed discs dominated by streaming turbulence. In general, increased dust diffusion from stronger gradients can lower the concentration of dust filaments and can explain the higher solid abundances needed to trigger strong particle clumping and the reduced planetesimal formation efficiency previously found in vertically-stratified simulations.
We discuss F-theory backgrounds associated to flat torus bundles over Ricci-flat manifolds. In this setting the F-theory background can be understood as a IIB orientifold with a large radius limit described by a supersymmetric compactification of IIB supergravity on a smooth, Ricci flat, but in general non-spin geometry. When compactified on an additional circle these backgrounds are T-dual to IIA compactifications on smooth non-orientable manifolds with a Pin− structure.
Prevailing $N$-body planet formation models typically start with lunar-mass embryos and show a general trend of rapid migration of massive planetary cores to the inner Solar System in the absence of a migration trap. This setup cannot capture the evolution from a planetesimal to embryo, which is crucial to the final architecture of the system. We aim to model planet formation with planet migration starting with planetesimals of $\sim10^{-6}$ -- $10^{-4}M_\oplus$ and reproduce the giant planets of the Solar System. We simulated a population of 1,000 -- 5,000 planetesimals in a smooth protoplanetary disc, which was evolved under the effects of their mutual gravity, pebble accretion, gas accretion, and planet migration, employing the parallelized $N$-body code SyMBAp. We find that the dynamical interactions among growing planetesimals are vigorous and can halt pebble accretion for excited bodies. While a set of results without planet migration produces one to two gas giants and one to two ice giants beyond 6 au, massive planetary cores readily move to the inner Solar System once planet migration is in effect. Dynamical heating is important in a planetesimal disc and the reduced pebble encounter time should be considered in similar models. Planet migration remains a challenge to form cold giant planets in a smooth protoplanetary disc, which suggests an alternative mechanism is required to stop them at wide orbits.
We analyze TYPHOON long-slit-absorption line spectra of the starburst barred spiral galaxy NGC 1365 obtained with the Progressive Integral Step Method covering an area of 15 kpc2. Applying a population synthesis technique, we determine the spatial distribution of ages and metallicities of the young and old stellar populations together with star formation rates, reddening, extinction, and the ratio R V of extinction to reddening. We detect a clear indication of inside-out growth of the stellar disk beyond 3 kpc characterized by an outward increasing luminosity fraction of the young stellar population, a decreasing average age, and a history of mass growth, which was finished 2 Gyr later in the outermost disk. The metallicity of the young stellar population is clearly super solar but decreases toward larger galactocentric radii with a gradient of -0.02 dex kpc-1. On the other hand, the metal content of the old population does not show a gradient and stays constant at a level roughly 0.4 dex lower than that of the young population. In the center of NGC 1365, we find a confined region where the metallicity of the young population drops dramatically and becomes lower than that of the old population. We attribute this to the infall of metal-poor gas, and additionally, to interrupted chemical evolution where star formation is stopped by active galactic nuclei and supernova feedback and then after several gigayears resumes with gas ejected by stellar winds from earlier generations of stars. We provide a simple model calculation as support for the latter.
The existence of a nucleon-ϕ (N-ϕ) bound state has been subject of theoretical and experimental investigations for decades. In this letter, indication of a p-ϕ bound state is found, using for the first time two-particle correlation functions as alternative to invariant mass spectra. Newly available lattice calculations for the spin 3/2 N-ϕ interaction by the HAL QCD collaboration are used to constrain the spin 1/2 counterpart from the fit of the experimental p-ϕ correlation function measured by ALICE. The corresponding scattering length and effective range are f0(1/2) = (-1.54-0.53+0.53 (stat .)-0.09+0.16 (syst .) + i ⋅0.00-0.00+0.35 (stat .)-0.00+0.16 (syst .)) fm and d0(1/2) = (0.39-0.09+0.09 (stat .)-0.03+0.02 (syst .) + i ⋅ 0.00-0.04+0.00 (stat .)-0.02+0.00 (syst .)) fm, respectively. The results imply the appearance of a p-ϕ bound state with an estimated binding energy in the range of 12.8 - 56.1 MeV.
The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.
Flavor-dependent neutrino transport is described by a well-known kinetic equation for occupation-number matrices in flavor space. However, as an overlooked theoretical problem, we show that in the inhomogeneous case, neutrino-neutrino refractive energy is not conserved. We derive the missing gradient terms in the fast flavor limit (vanishing neutrino masses), and prove that the missing refractive energy is traded with the huge reservoir of neutrino kinetic energy through gradients of the weak interaction potential. Even small changes of the kinetic energy accommodate the refractive energy gained or lost. Flavor evolution alone is negligibly affected by the new terms.
There is currently no evidence for a baryon asymmetry in our Universe. Instead, cosmological observations have only demonstrated the existence of a quark-antiquark asymmetry, which does not necessarily imply a baryon asymmetric Universe, since the baryon number of the dark sector particles is unknown. In this paper we discuss a framework where the total baryon number of the Universe is equal to zero, and where the observed quark-antiquark asymmetry arises from neutron portal interactions with a dark sector fermion N that carries baryon number. In order to render a baryon symmetric universe throughout the whole cosmological history, we introduce a complex scalar χ, with opposite baryon number and with the same initial abundance as N. Notably, due to the baryon number conservation, χ is absolutely stable and could have an abundance today equal to the observed dark matter abundance. Therefore, in this simple framework, the existence of a quark-antiquark asymmetry is intimately related to the existence (and the stability) of dark matter.
We present initial results from a JWST survey of the youngest Galactic core-collapse supernova remnant Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI/MRS IFU spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include: 1) a web-like network of unshocked ejecta filaments resolved to 0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, neutrino and radioactively heated high-entropy matter, 2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (approximately one arcsecond) round holes formed by knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks, 3) dozens of light echoes with angular sizes between 0.1 arcsecond to 1 arcminute reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.
We present forecasts for constraints on the Hu \& Sawicki $f(R)$ modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the $f(R)$ gravity parameter $\log_{10}|f_{R0}|$. In this analysis we vary the cosmological parameters $[\Omega_{\rm m}, \Omega_\nu h^2, h^2, A_s, n_s, \log_{10}|f_{R0}|]$, which allows us to account for possible degeneracies between cosmological parameters and $f(R)$ modified gravity. The analysis accounts for $f(R)$ gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT-3G$\times$ngWL is expected to obtain an upper limit of $\log_{10}|f_{R0}| < -5.95$ at $95\,\%$ credibility, which significantly improves upon the current best bounds. The CMB-S4$\times$ngWL dataset is expected to improve this even further to $\log_{10}|f_{R0}| < -6.23$. Furthermore, $f(R)$ gravity models with $\log_{10}|f_{R0}| \geq -6$, which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between $\log_{10}|f_{R0}|$ and $\Omega_{\mathrm{m}}$ as well as $\sigma_8$ for $\log_{10}|f_{R0}| > -6$ and $\log_{10}|f_{R0}| > -5$ respectively. Our forecasts indicate that future cluster abundance studies of $f(R)$ gravity will enable substantially improved constraints that are competitive with other cosmological probes.
We review the current status of astrophysical bounds on QCD axions, primarily based on the observational effects of nonstandard energy losses on stars, including black-hole superradiance. Over the past few years, many of the traditional arguments have been reexamined both theoretically and using modern data and new ideas have been put forth. This compact review updates similar Lecture Notes written by one of us in 2006 [Lect. Notes Phys. 741 (2008) 51-71].
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), which determines the kicks of NSs born from stars near the low-mass end of SN progenitors. In more massive progenitors aspherical accretion downflows around the PNS can also lead to anisotropic neutrino emission (absorption) with a neutrino-induced NS kick roughly opposite to (aligned with) the kick by asymmetric mass ejection. We estimate upper bounds for the final neutrino kicks of 150-260 km/s, whereas the hydrodynamic kicks can reach up to more than 1300 km/s. Therefore the hydrodynamic kicks dominate for NSs from explosions of higher-mass progenitors, whereas the neutrino kicks dominate in the case of NSs from the lowest-mass progenitors. Our models suggest that the Crab pulsar as a representative of the latter category could have received its velocity of about 160 km/s by a neutrino kick due to the LESA asymmetry. Such neutrino kicks of 100-200 km/s define a nearly ubiquitous floor value, which may shed new light on the origin of pulsars in globular clusters. 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.
Biomolecular condensates help organize the cell cytoplasm and nucleoplasm into spatial compartments with different chemical compositions. A key feature of such compositional patterning is the local enrichment of enzymatically active biomolecules which, after transient binding via molecular interactions, catalyze reactions among their substrates. Thereby, biomolecular condensates provide a spatial template for non-uniform concentration profiles of substrates. In turn, the concentration profiles of substrates, and their molecular interactions with enzymes, drive enzyme fluxes which can enable novel non-equilibrium dynamics. To analyze this generic class of systems, with a current focus on self-propelled droplet motion, we here develop a self-consistent sharp interface theory. In our theory, we diverge from the usual bottom-up approach, which involves calculating the dynamics of concentration profiles based on a given chemical potential gradient. Instead, reminiscent of control theory, we take the reverse approach by deriving the chemical potential profile and enzyme fluxes required to maintain a desired condensate form and dynamics. The chemical potential profile and currents of enzymes come with a corresponding power dissipation rate, which allows us to derive a thermodynamic consistency criterion for the passive part of the system (here, reciprocal enzyme-enzyme interactions). As a first use case of our theory, we study the role of reciprocal interactions, where the transport of substrates due to reactions and diffusion is, in part, compensated by redistribution due to molecular interactions. More generally, our theory applies to mass-conserved active matter systems with moving phase boundaries.
The direct detection of core-collapse supernova (SN) progenitor stars is a powerful way of probing the last stages of stellar evolution. However, detections in archival Hubble Space Telescope images are limited to about one detection per year. Here, we explore whether we can increase the detection rate by using data from ground-based wide-field surveys. Due to crowding and atmospheric blurring, progenitor stars can typically not be identified in preexplosion images alone. Instead, we combine many pre-SN and late-time images to search for the disappearance of the progenitor star. As a proof of concept, we implement our search of ZTF data. For a few hundred images, we achieve limiting magnitudes of ~23 mag in the g and r bands. However, no progenitor stars or long-lived outbursts are detected for 29 SNe within z ≤ 0.01, and the ZTF limits are typically several magnitudes less constraining than detected progenitors in the literature. Next, we estimate progenitor detection rates for the Legacy Survey of Space and Time (LSST) with the Vera C. Rubin telescope by simulating a population of nearby SNe. The background from bright host galaxies reduces the nominal LSST sensitivity by, on average, 0.4 mag. Over the 10 yr survey, we expect the detection of ~50 red supergiant progenitors and several yellow and blue supergiants. The progenitors of Type Ib and Ic SNe will be detectable if they are brighter than -4.7 or -4.0 mag in the LSST i band, respectively. In addition, we expect the detection of hundreds of pre-SN outbursts depending on their brightness and duration.
We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including the Dark Energy Survey (DES), DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of z = 2.35 to a depth of m r = 27. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group-galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to further improve the color assignment model and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the DES data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys.
The J-region Asymptotic Giant Branch (JAGB) method is a standard candle that leverages the constant luminosities of color-selected, carbon-rich AGB stars, measured in the near-infrared at 1.2 μm. The Chicago-Carnegie Hubble Program has obtained JWST imaging of the SN Ia host galaxies NGC 7250, NGC 4536, and NGC 3972. With these observations, the JAGB method can be studied for the first time using JWST. Lee et al. demonstrated the JAGB magnitude is optimally measured in the outer disks of galaxies, because in the inner regions the JAGB magnitude can vary significantly due to a confluence of reddening, blending, and crowding effects. However, determining where the "outer disk" lies can be subjective. Therefore, we introduce a novel method for systematically selecting the outer disk. In a given galaxy, the JAGB magnitude is first separately measured in concentric regions, and the "outer disk" is then defined as the first radial bin where the JAGB magnitude stabilizes to a few hundredths of a magnitude. After successfully employing this method in our JWST galaxy sample, we find the JAGB stars are well segregated from other stellar populations in color-magnitude space, and have observed dispersions about their individual F115W modes of σ N7250 = 0.32 mag, σ N4536 = 0.34 mag, and σ N3972 = 0.35 mag. These measured dispersions are similar to the scatter measured for the JAGB stars in the LMC using 2MASS data (σ = 0.33 mag). In conclusion, the JAGB stars as observed with JWST clearly demonstrate their considerable power both as high-precision extragalactic distance indicators and as SN Ia supernova calibrators.
In dense neutrino environments like core-collapse supernovae (CCSNe) and neutron star mergers (NSMs), neutrinos can undergo fast flavor conversions (FFC) when their angular distribution of neutrino electron lepton number ($\nu$ELN) crosses zero along some directions. While previous studies have demonstrated the detection of axisymmetric $\nu$ELN crossings in these extreme environments, non-axisymmetric crossings have remained elusive, mostly due to the absence of models for their angular distributions. In this study, we present a pioneering analysis of the detection of non-axisymmetric $\nu$ELN crossings using machine learning (ML) techniques. Our ML models are trained on data from two CCSN simulations, one with rotation and one without, where non-axisymmetric features in neutrino angular distributions play a crucial role. We demonstrate that our ML models achieve detection accuracies exceeding 90\%. This is an important improvement, especially considering that a significant portion of $\nu$ELN crossings in these models eluded detection by earlier methods.
Dark matter particles captured in neutron stars deposit their energy as heat. This DM heating effect can be observed only if it dominates over other internal heating effects in neutron stars. In this work, as an example of such an internal heating source, we consider the frictional heating caused by the creep motion of neutron superfluid vortex lines in the neutron star crust. The luminosity of this heating effect is controlled by the strength of the interaction between the vortex lines and nuclei in the crust, which can be estimated from the many-body calculation of a high-density nuclear system as well as through the temperature observation of old neutron stars. We show that both the temperature observation and theoretical calculation suggest that the vortex creep heating dominates over the DM heating. The vortex-nuclei interaction must be smaller than the estimated values by several orders of magnitude to overturn this.
When hypothetical neutrino secret interactions (ν SI ) are large, they form a fluid in a supernova (SN) core, flow out with sonic speed, and stream away as a fireball. For the first time, we tackle the complete dynamical problem and solve all steps, systematically using relativistic hydrodynamics. The impact on SN physics and the neutrino signal is remarkably small. For complete thermalization within the fireball, the observable spectrum changes in a way that is independent of the coupling strength. One potentially large effect beyond our study is quick deleptonization if ν SI violate lepton number. By present evidence, however, SN physics leaves open a large region in parameter space, where laboratory searches and future high-energy neutrino telescopes will probe ν SI .
Neutrino-neutrino scattering could have a large secret component that would turn neutrinos within a supernova (SN) core into a self-coupled fluid. Neutrino transport within the SN core, emission from its surface, expansion into space, and the flux spectrum and time structure at Earth might all be affected. We examine these questions from first principles. First, diffusive transport differs only by a modified spectral average of the interaction rate. We next study the fluid energy transfer between a hot and a cold blackbody surface in plane-parallel and spherical geometry. The key element is the decoupling process within the radiating bodies, which themselves are taken to be isothermal. For a zero-temperature cold plate, mimicking radiation into free space by the hot plate, the energy flux is 3%-4% smaller than the usual Stefan-Boltzmann law. The fluid energy density just outside the hot plate is numerically 0.70 of the standard case, the outflow velocity is the speed of sound vs=c /√{3 } , conspiring to a nearly unchanged energy flux. Our results provide the crucial boundary condition for the expansion of the self-interacting fluid into space, assuming an isothermal neutrino sphere. We also derive a dynamical solution, assuming the emission suddenly begins at some instant. A neutrino front expands in space with luminal speed, whereas the outflow velocity at the radiating surface asymptotically approaches vs from above. Asymptotically, one thus recovers the steady-state emission found in the two-plate model. A sudden end to neutrino emission leads to a fireball with constant thickness equal to the duration of neutrino emission.
Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave (GW) sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the DRAGON-II simulation database, a suite of 19 direct N-body models representing dense star clusters with up to 106 stars and $<33~{{\%}}$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30~{{\%}}$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around 20 % of DRAGON-II mergers are eccentric in the LISA band and 5 % in the LIGO band. We infer a mean cosmic merger rate of $\mathcal {R}\sim 30(4.4)(1.2)$ yr-1 Gpc-3 for BBHs, NS-BH, and WD-BH binary mergers, respectively, and discuss the prospects for multimessenger detection of WD-BH binaries with LISA. We model the rate of pair-instability supernovae (PISNe) in star clusters and find that surveys with a limiting magnitude mbol = 25 can detect ~1 - 15 yr-1 PISNe. Comparing these estimates with future observations could help to pin down the impact of massive star evolution on the mass spectrum of compact stellar objects in star clusters.
We present the data release of the Uchuu-SDSS galaxies: a set of 32 high-fidelity galaxy lightcones constructed from the large Uchuu 2.1 trillion particles N-body simulation using Planck cosmology. We adopt subhalo abundance matching to populate the Uchuu-box halo catalogues with SDSS galaxy luminosities. These box catalogues generated at several redshifts are combined to create a set of lightcones with redshift-evolving galaxy properties. The Uchuu-SDSS galaxy lightcones are built to reproduce the footprint and statistical properties of the SDSS main galaxy survey, along with stellar masses and star formation rates. This facilitates a direct comparison of the observed SDSS and simulated Uchuu-SDSS data. Our lightcones reproduce a large number of observational results, such as the distribution of galaxy properties, galaxy clustering, stellar mass functions, and halo occupation distributions. Using simulated and real data we select samples of bright red galaxies at zeff = 0.15 to explore Redshift Space Distortions and Baryon Acoustic Oscillations (BAO) by fitting the full two-point correlation function and the BAO peak. We create a set of 5100 galaxy lightcones using GLAM N-body simulations to compute covariance errors. We report a $\sim 30~{{\%}}$ precision increase on fσ8 and the pre-reconstruction BAO scale, due to our better estimate of the covariance matrix. From our BAO-inferred α∥ and α⊥ parameters, we obtain the first SDSS measurements of the Hubble and angular diameter distances $D_\mathrm{H}(z=0.15) / r_d = 27.9^{+3.1}_{-2.7}$, $D_\mathrm{M}(z=0.15) / r_d = 5.1^{+0.4}_{-0.4}$. Overall, we conclude that the Planck Λ CDM cosmology nicely explains the observed large-scale structure statistics of SDSS. All data sets are made publicly available.
Long-lived radioactive by-products of nucleosynthesis provide an opportunity to trace the flow of ejecta away from its sources for times beyond where ejecta can be seen otherwise. Gamma rays from such radioactive decay in interstellar space can be measured with space-borne telescopes. A prominent useful example is 26Al with a radioactive decay time of one My. Such observations have revealed that typical surroundings of massive stars are composed of large cavities, extending to kpc sizes. Implications are that material recycling into new stars is twofold: rather direct as parental clouds are hosts to new star formation triggered by feedback, and more indirect as these large cavities merge with ambient interstellar gas after some delay. Kinematic measurements of hot interstellar gas carrying such ejecta promises important measurements complementing stellar and dense gas kinematics.
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.
The application of Graph Neural Networks (GNN) in track reconstruction is a promising approach to cope with the challenges arising at the High-Luminosity upgrade of the Large Hadron Collider (HL-LHC). GNNs show good track-finding performance in high-multiplicity scenarios and are naturally parallelizable on heterogeneous compute architectures. Typical high-energy-physics detectors have high resolution in the innermost layers to support vertex reconstruction but lower resolution in the outer parts. GNNs mainly rely on 3D space-point information, which can cause reduced track-finding performance in the outer regions. In this contribution, we present a novel combination of GNN-based track finding with the classical Combinatorial Kalman Filter (CKF) algorithm to circumvent this issue: The GNN resolves the track candidates in the inner pixel region, where 3D space points can represent measurements very well. These candidates are then picked up by the CKF in the outer regions, where the CKF performs well even for 1D measurements. Using the ACTS infrastructure, we present a proof of concept based on truth tracking in the pixels as well as a dedicated GNN pipeline trained on $t\bar{t}$ events with pile-up 200 in the OpenDataDetector.
We compute the static force on the lattice in the quenched case directly through generalized Wilson loops. We modify the Wilson loop by inserting an $E$-field component on one of the temporal Wilson lines. However, chromo-field components prevent us from performing the continuum limit properly, hence, we use gradient flow to renormalize the field insertion. As a result, we obtain continuum results and compare them to perturbative expression to extract $\Lambda_0$, and we predict the value $\sqrt{8t_0} \Lambda_{\overline{\textrm{MS}}}^{n_f=0} =0.629^{+22}_{-26}$. This work serves as preparation for similar operators with field insertions required in nonrelativistic effective field theories.
Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms.
Aims: We aim to characterise the young binary system HD 34700 A, which shows a wealth of structures.
Methods: Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and of ALMA observations, we analyse this system at multiple wavelengths. We study the morphology of the rings and spiral arms and the scattering properties of the dust. We discuss the possible causes of all the observed features.
Results: We detect for the first time, in the Hα band, a ring extending from ~65 au to ~120 au, inside the ring which is already known from recent studies. These two have different physical and geometrical properties. Based on the scattering properties, the outer ring may consist of grains with a typical size of aout ≥ 4 µm, while the inner ring has a smaller typical size of ain ≤ 0.4 µm. Two extended logarithmic spiral arms stem from opposite sides of the disk. The outer ring appears as a spiral arm itself, with a variable radial distance from the centre and extended substructures. ALMA data confirm the presence of a millimetric dust substructure centred just outside the outer ring, and detect misaligned gas rotation patterns for HD 34700 A and B.
Conclusions: The complexity of HD 34700 A, revealed by the variety of observed features, suggests the existence of one or more disk-shaping physical mechanisms. Our findings are compatible with the presence inside the disk of an as of yet undetected planet of several Jupiter masses and the system interaction with the surroundings, by means of gas cloudlet capture or flybys. Further observations with JWST/MIRI or ALMA (gas kinematics) could shed more light on them.
When a galaxy falls into a cluster, its outermost parts are the most affected by the environment. In this paper, we are interested in studying the influence of a dense environment on different galaxy's components to better understand how this affects the evolution of galaxies. We use, as laboratory for this study, the Hydra cluster which is close to virialization; yet it still shows evidence of substructures. We present a multiwavelength bulge-disc decomposition performed simultaneously in 12 bands from S-PLUS (Southern Photometric Local Universe Survey) data for 52 galaxies brighter than mr = 16. We model the galaxies with a Sérsic profile for the bulge and an exponential profile for the disc. We find that the smaller, more compact, and bulge-dominated galaxies tend to exhibit a redder colour at a fixed stellar mass. This suggests that the same mechanisms (ram-pressure and tidal stripping) that are causing the compaction in these galaxies are also causing them to stop forming stars. The bulge size is unrelated to the galaxy's stellar mass, while the disc size increases with greater stellar mass, indicating the dominant role of the disc in the overall galaxy mass-size relation found. Furthermore, our analysis of the environment unveils that quenched galaxies are prevalent in regions likely associated with substructures. However, these areas also harbour a minority of star-forming galaxies, primarily resulting from galaxy interactions. Lastly, we find that ~37 per cent of the galaxies exhibit bulges that are bluer than their discs, indicative of an outside-in quenching process in this type of dense environments.
To date, the hot Jupiter WASP-12 b has been the only planet with confirmed orbital decay. The late F-type host star has been hypothesized to be surrounded by a large structure of circumstellar material evaporated from the planet. We obtained two high-resolution spectral transit time series with CARMENES and extensively searched for absorption signals by the atomic species Na, H, Ca, and He using transmission spectroscopy, thereby covering the He I triplet with high resolution for the first time. We apply SYSREM for atomic line transmission spectroscopy, introduce the technique of signal protection to improve the results for individual absorption lines, and compare the outcomes to those of established methods. No transmission signals were detected and the most stringent upper limits as of yet were derived for the individual indicators. Nonetheless, we found variation in the stellar Halpha and He I lines, the origin of which remains uncertain but is unlikely to be activity. To constrain the enigmatic activity state of WASP-12, we analyzed XMM-Newton X-ray data and found the star to be moderately active at most. We deduced an upper limit for the X-ray luminosity and the irradiating X-ray and extreme ultraviolet (XUV) flux of WASP-12 b. Based on the XUV flux upper limit and the lack of the He I signal, our hydrodynamic models slightly favor a moderately irradiated planet with a thermospheric temperature of <= 12,000 K, and a conservative upper limit of <= 4e12 g/s on the mass-loss rate. Our study does not provide evidence for an extended planetary atmosphere or absorption by circumstellar material close to the planetary orbit.
Life continuously transduces energy to perform critical functions using energy stored in reactive molecules like ATP or NADH. ATP dynamically phosphorylates active sites on proteins and thereby regulates their function. Inspired by such machinery, regulating supramolecular functions using energy stored in reactive molecules has gained traction. Enzyme-free, synthetic systems that use dynamic phosphorylation to regulate supramolecular processes do not exist.
We present an enzyme-free reaction cycle that consumes phosphorylating agents by transiently phosphorylating amino acids. The phosphorylated amino acids are labile and deactivate through hydrolysis. The cycle exhibits versatility and tunability, allowing for the dynamic phosphorylation of multiple precursors with a tunable half-life. Notably, we show the resulting phosphorylated products can regulate the peptide’s phase separation, leading to active droplets that require the continuous conversion of fuel to sustain. Our new reaction cycle will be valuable as a model for biological phosphorylation but can also offer insights into protocell formation.
Unobscured quasars (QSOs) are predicted to be the final stage in the evolutionary sequence from gas-rich mergers to gas-depleted, quenched galaxies. Studies of this population, however, find a high incidence of far-infrared-luminous sources -suggesting significant dust-obscured star formation-but direct observations of the cold molecular gas fuelling this star formation are still necessary. We present a NOEMA study of CO(2-1) emission, tracing the cold molecular gas, in ten lensed z=1-1.5 unobscured QSOs. We detected CO(2-1) in seven of our targets, four of which also show continuum emission (\lambda_rest = 1.3mm). After subtracting the foreground galaxy contribution to the photometry, spectral energy distribution fitting yielded stellar masses of 10^9-11 M_\odot, with star formation rates of 25-160 M_\odot yr^-1 for the host galaxies. These QSOs have lower $L'_\mathrm{CO}$ than star-forming galaxies with the same L_IR, and show depletion times spanning a large range (50-900 Myr), but with a median of just 90 Myr. We find molecular gas masses in the range 2-40 x 10^9(alpha_CO/4) M_\odot, which suggest gas fractions above ~50% for most of the targets. Despite the presence of an unobscured QSO, the host galaxies are able to retain significant amounts of cold gas. However, with a median depletion time of ~90 Myr, the intense burst of star formation taking place in these targets will quickly deplete their molecular gas reservoirs in the absence of gas replenishment, resulting in a quiescent host galaxy. The non-detected QSOs are three of the four radio-loud QSOs in the sample, and their properties indicate that they are likely already transitioning into quiescence. Recent cosmological simulations tend to overestimate the depletion times expected for these z~1 QSO-host galaxies, which is likely linked to their difficulty producing starbursts across the general high-redshift galaxy population.
We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree.
In this first of a series of papers related to long-period post-common-envelope (CE) binaries, we investigated whether extra energy is required or not to explain the currently known post-CE binaries with sufficiently long orbital periods consisting of oxygen-neon white dwarfs with AFGK-type main-sequence star companions. We carried out binary population simulations with the BSE code and searched for their formation pathways. Unlike what has been claimed for a long time, we show that all such post-CE binaries can be explained by assuming inefficient CE evolution, which is consistent with results achieved for the remaining post-CE binaries. There is therefore no need for an extra energy source. We also found that for CE efficiency close to 100%, post-CE binaries hosting oxygen-neon white dwarfs with orbital periods as long as a thousand days can be explained. For all known systems we found formation pathways consisting of CE evolution triggered when a highly evolved (i.e. the envelope mass being comparable to the core mass) thermally-pulsing asymptotic giant branch star fills its Roche lobe at an orbital period of several thousand days. Due to the sufficiently low envelope mass and sufficiently long orbital period, the resulting post-CE orbital period can easily be several tens of days. We conclude that the known post-CE binaries with oxygen-neon white dwarfs and AFGK-type main-sequence stars can be explained without invoking any energy source other than orbital and thermal energy. Our results strengthen the idea that the most common formation pathway of the overall population of post-CE binaries hosting white dwarfs is through inefficient CE evolution.
We present a Lagrangian approach to counting degrees of freedom in first-order field theories. The emphasis is on the systematic attainment of a complete set of constraints. In particular, we provide the first comprehensive procedure to ensure the functional independence of all constraints and discuss in detail the possible closures of the constraint algorithm. We argue degrees of freedom can but need not correspond to physical modes. The appendix comprises fully worked out, physically relevant examples of varying complexity.
J191213.72 - 441045.1 is a binary system composed of a white dwarf and an M-dwarf in a 4.03-h orbit. It shows emission in radio, optical, and X-ray, all modulated at the white dwarf spin period of 5.3 min, as well as various orbital sideband frequencies. Like in the prototype of the class of radio-pulsing white dwarfs, AR Scorpii, the observed pulsed emission seems to be driven by the binary interaction. In this work, we present an analysis of far-ultraviolet spectra obtained with the Cosmic Origins Spectrograph at the Hubble Space Telescope, in which we directly detect the white dwarf in J191213.72 - 441045.1. We find that the white dwarf has a temperature of Teff = 11485 ± 90 K and mass of 0.59 ± 0.05 M⊙. We place a tentative upper limit on the magnetic field of ≈50 MG. If the white dwarf is in thermal equilibrium, its physical parameters would imply that crystallization has not started in the core of the white dwarf. Alternatively, the effective temperature could have been affected by compressional heating, indicating a past phase of accretion. The relatively low upper limit to the magnetic field and potential lack of crystallization that could generate a strong field pose challenges to pulsar-like models for the system and give preference to propeller models with a low magnetic field. We also develop a geometric model of the binary interaction which explains many salient features of the system.
Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs. Aims. In order to ascertain the infalling nature of observed streamer-like structures and then systematically characterize their dynamics, we developed the code TIPSY (Trajectory of Infalling Particles in Streamers around Young stars). Methods. Using TIPSY, the streamer molecular line emission is first isolated from the disk emission. Then the streamer emission, which is effectively a point cloud in three-dimensional (3D) position-position-velocity space, is simplified to a curve-like representation. The observed streamer curve is then compared to the theoretical trajectories of infalling material. The best-fit trajectories are used to constrain streamer features, such as the specific energy, the specific angular momenta, the infall timescale, and the 3D morphology. Results. We used TIPSY to fit molecular-line ALMA observations of streamers around a Class II binary system, S CrA, and a Class I/II protostar, HL Tau. Our results indicate that both of the streamers are consistent with infalling motion. TIPSY results and mass estimates suggest that S CrA and HL Tau are accreting material at a rate of $\gtrsim27$ M$_{jupiter}$ Myr$^{-1}$ and $\gtrsim5$ M$_{jupiter}$ Myr$^{-1}$, respectively, which can significantly increase the mass budget available to form planets. Conclusions. TIPSY can be used to assess whether the morphology and kinematics of observed streamers are consistent with infalling motion and to characterize their dynamics, which is crucial for quantifying their impact on the protostellar systems.
We compute the planar three-loop Quantum Chromodynamics (QCD) corrections to the helicity amplitudes involving a vector boson V = Z ,W± ,γ⁎, two quarks and a gluon. These amplitudes are relevant to vector-boson-plus-jet production at hadron colliders and other precision QCD observables. The planar corrections encompass the leading colour factors N3, N2Nf, NNf2 and Nf3 . We provide the finite remainders of the independent helicity amplitudes in terms of multiple polylogarithms, continued to all kinematic regions and in a form which is compact and lends itself to efficient numerical evaluation. The presented amplitude respects the conjectured symbol-adjacency constraints for amplitudes with three massless and one massive leg.
The vast majority of close binaries containing a compact object form through common-envelope (CE) evolution. Despite this importance, we struggle to even understand the energy budget of CE evolution. For decades, observed long-period post-CE binaries have been interpreted as evidence for additional energies to contribute during CE evolution. We have recently shown that this argument is based on simplified assumptions for all long-period post-CE binaries containing massive white dwarfs. The only remaining post-CE binary star that has been claimed to require contributions from additional energy sources to understand its formation is KOI 3278. Here we address in detail the potential evolutionary history of KOI 3278. In particular, we investigated whether extra energy sources, such as recombination energy, are indeed required to explain its existence. We used the 1D stellar evolution code MESA to carry out binary evolution simulations and searched for potential formation pathways for KOI 3278 that are able to explain its observed properties. We found that KOI 3278 can be explained if the white dwarf progenitor filled its Roche lobe during a helium shell flash. In this case, the orbital period of KOI 3278 can be reproduced if the CE binding energy is calculated taking into account gravitational energy and thermodynamic internal energy. While the CE evolution that led to the formation of KOI 3278 must have been efficient, that is, most of the available orbital energy must have been used to unbind the CE, recombination energy is not required. We conclude that currently not a single observed post-CE binary requires to assume energy sources other than gravitational and thermodynamic energy to contribute to CE evolution. KOI 3278, however, remains an intriguing post-CE binary as, unlike its siblings, understanding its existence requires highly efficient CE ejection.
The eROSITA instrument aboard the Spectrum Roentgen Gamma (SRG) satellite has performed its first all-sky survey between December 2019 and June 2020. This paper presents the resulting hard X-ray (2.3-5 keV) sample, the first created from an all-sky imaging survey in the 2-8 keV band, for sources within western galactic sky. The 5466 hard X-ray selected sources detected with eROSITA are presented and discussed. The Bayesian statistics-based code NWAY is used to identify the counterparts for the X-ray sources. These sources are classified based on their multiwavelength properties, and the literature is searched to identify spectroscopic redshifts, which further inform the source classification. A total of 2547 sources are found to have good-quality counterparts, and 111 of these are detected only in the hard band. Comparing with other hard X-ray selected surveys, the eROSITA hard sample covers a larger redshift range and probes dimmer sources, providing a complementary and expanded sample as compared to Swift-BAT. Examining the column density distribution of missed and detected eROSITA sources present in the follow-up catalog of Swift BAT 70 month sources, it is demonstrated that eROSITA can detect obscured sources with column densities $>10^{24}$ cm$^{-2}$, but that the completeness drops rapidly after $10^{23}$ cm$^{-2}$. A sample of hard-only sources, many of which are likely to be heavily obscured AGN, is also presented and discussed. X-ray spectral fitting reveals that these sources have extremely faint soft X-ray emission and their optical images suggest that they are found in more edge-on galaxies with lower b/a. The resulting X-ray catalog is demonstrated to be a powerful tool for understanding AGN, in particular heavily obscured AGN found in the hard-only sample.
Cold dark matter axions produced in the post-inflationary Peccei-Quinn symmetry breaking scenario serve as clear targets for their experimental detection, since it is in principle possible to give a sharp prediction for their mass once we understand precisely how they are produced from the decay of global cosmic strings in the early Universe. In this paper, we perform a dedicated analysis of the spectrum of axions radiated from strings based on large scale numerical simulations of the cosmological evolution of the Peccei-Quinn field on a static lattice. Making full use of the massively parallel code and computing resources, we executed the simulations with up to $11264^3$ lattice sites, which allows us to improve our understanding of the dependence on the parameter controlling the string tension and thus give a more accurate extrapolation of the numerical results. We found that there are several systematic effects that have been overlooked in previous works, such as the dependence on the initial conditions, contaminations due to oscillations in the spectrum, and discretisation effects, some of which could explain the discrepancy in the literature. We confirmed the trend that the spectral index of the axion emission spectrum increases with the string tension, but did not find a clear evidence of whether it continues to increase or saturates to a constant at larger values of the string tension due to the severe discretisation effects. Taking this uncertainty into account and performing the extrapolation with a simple power law assumption on the spectrum, we find that the dark matter mass is predicted in the range of $m_a \approx 95$-$450\,\mu\mathrm{eV}$.
Context. About 30% - 40% of the baryons in the local Universe remain unobserved. Many of these "missing" baryons are expected to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments that connect clusters of galaxies. SRG/eROSITA performance verification (PV) observations covered about 15 square degrees of the A3391/95 system and have revealed a ~15 Mpc continuous soft emission connecting several galaxy groups and clusters.
Aims: We aim to characterize the gas properties in the cluster outskirts (R500 < r < R200) and in the detected inter-cluster filaments (> R200) and to compare them to predictions.
Methods: We performed X-ray image and spectral analyses using the eROSITA PV data in order to assess the gas morphology and properties in the outskirts and the filaments in the directions of the previously detected Northern and Southern Filament of the A3391/95 system. We constructed surface brightness profiles using particle-induced background-subtracted, exposure- and Galactic absorption-corrected eROSITA products in the soft band (0.3-2.0 keV). We constrained the temperatures, metallicities, and electron densities through X-ray spectral fitting and compared them with the expected properties of the WHIM. We took particular care of the foreground.
Results: In the filament-facing outskirts of A3391 and the Northern Clump, we find higher temperatures than typical cluster outskirts profiles, with a significance of between 1.6 and 2.8σ, suggesting heating due to their connections with the filaments. We confirm surface brightness excess in the profiles of the Northern, Eastern, and Southern Filaments. From spectral analysis, we detect hot gas of 0.96−0.14+0.17 keV and 1.09−0.06+0.09 for the Northern and Southern Filament, respectively, which are close to the upper WHIM temperature limit. The filament metallicities are below 10% solar metallicity and the electron densities are ranging between 2.6 and 6.3 × 10−5 cm−3. The characteristic properties of the Little Southern Clump (LSC), which is located at a distance of ~1.5R200 from A3395S in the Southern Filament, suggest that it is a small galaxy group. Excluding the LSC from the analysis of the Southern Filament does not significantly change the temperature or metallicity of the gas, but it decreases the gas density by 30%. This shows the importance of taking into account any clumps in order to avoid overestimation of the gas measurement in the outskirts and filament regions.
Conclusions: We present measurements of morphology, temperature, metallicity, and density of individual warm-hot filaments. The electron densities of the filaments are consistent with the WHIM properties as predicted by cosmological simulations, but the temperatures are higher. As both filaments are short (1.8 and 2.7 Mpc) and located in a denser environment, stronger gravitational heating may be responsible for this temperature enhancement. The metallicities are low, but still within the expected range from the simulations.
Image that is displayed in Fig. 1 is 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/681/A108
We present homogeneous multiband (grizJHKs) time-series observations of 78 Cepheids including 49 fundamental mode variables and 29 first-overtone mode variables. These observations were collected simultaneously using the ROS2 and REMIR instruments at the Rapid Eye Mount telescope. The Cepheid sample covers a large range of distances (0.5 - 19.7 kpc) with varying precision of parallaxes, and thus astrometry-based luminosity fits were used to derive PL and PW relations in optical Sloan (griz) and near-infrared (JHKs) filters. These empirically calibrated relations exhibit large scatter primarily due to larger uncertainties in parallaxes of distant Cepheids, but their slopes agree well with those previously determined in the literature. Using homogeneous high-resolution spectroscopic metallicities of 61 Cepheids covering -1.1 < [Fe/H] < 0.6 dex, we quantified the metallicity dependence of PL and PW relations which varies between $-0.30\pm0.11$ (in Ks) and $-0.55\pm0.12$ (in z) mag/dex in grizJHKs bands. However, the metallicity dependence in the residuals of the PL and PW relations is predominantly seen for metal-poor stars ([Fe/H] < -0.3 dex), which also have larger parallax uncertainties. The modest sample size precludes us from separating the contribution to the residuals due to parallax uncertainties, metallicity effects, and reddening errors. While this Cepheid sample is not optimal for calibrating the Leavitt law, upcoming photometric and spectroscopic datasets of the C-MetaLL survey will allow the accurate derivation of PL and PW relations in the Sloan and near-infrared bandpasses, which will be useful for the distance measurements in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time and upcoming extremely large telescopes.
Neutrinos can undergo fast flavor conversions (FFCs) within extremely dense astrophysical environments such as core-collapse supernovae (CCSNe) and neutron star mergers (NSMs). In this study, we explore FFCs in a \emph{multi-energy} neutrino gas, revealing that when the FFC growth rate significantly exceeds that of the vacuum Hamiltonian, all neutrinos (regardless of energy) share a common survival probability dictated by the energy-integrated neutrino spectrum. We then employ physics-informed neural networks (PINNs) to predict the asymptotic outcomes of FFCs within such a multi-energy neutrino gas. These predictions are based on the first two moments of neutrino angular distributions for each energy bin, typically available in state-of-the-art CCSN and NSM simulations. Our PINNs achieve errors as low as $\lesssim6\%$ and $\lesssim 18\%$ for predicting the number of neutrinos in the electron channel and the relative absolute error in the neutrino moments, respectively.
Multiparton interactions are a fascinating phenomenon that occur in almost every high-energy hadron-hadron collision yet are remarkably difficult to study quantitatively. In this Letter, we present a strategy to optimally disentangle multiparton interactions from the primary scattering in a collision. That strategy enables probes of multiparton interactions that are significantly beyond the state of the art, including their characteristic momentum scale, the interconnection between primary and secondary scatters, and the pattern of three and potentially even more simultaneous hard scatterings. This opens a path to powerful new constraints on multiparton interactions for LHC phenomenology and to the investigation of their rich field-theoretical structure.
In this work we demonstrate that significant gains in performance and data efficiency can be achieved in High Energy Physics (HEP) by moving beyond the standard paradigm of sequential optimization or reconstruction and analysis components. We conceptually connect HEP reconstruction and analysis to modern machine learning workflows such as pretraining, finetuning, domain adaptation and high-dimensional embedding spaces and quantify the gains in the example usecase of searches of heavy resonances decaying via an intermediate di-Higgs system to four $b$-jets.
We investigate the impact of AGN feedback, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms. We analyze deeper eROSITA observations of 1178 galaxy groups detected in eRASS1. We divide the sample into 271 subsamples and extract average thermodynamic properties, including electron density, temperature, and entropy at three characteristic radii along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. We present the tightest constraints on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium temperatures, while for the low mass groups, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. The comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL simulations) after the applications of the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Our measurements will pave the way for more realistic AGN feedback implementations in simulations.