A parsec-scale dusty torus is thought to be the cause of active galactic nuclei (AGN) dichotomy in the 1/2 types, narrow/broad emission lines. In a previous work, on the basis of parsec-scale resolution infrared/optical dust maps, it was found that dust filaments, few parsecs wide and several hundred parsecs long, were ubiquitous features crossing the centre of type 2 AGN, their optical thickness being sufficient to fully obscure the optical nucleus. This work presents the complementary view for type 1 and intermediate-type AGN. The same type of narrow, collimated, dust filaments are equally found at the centre of these AGN. The difference now resides in their location with respect to the nucleus, next to it but not crossing it, as it is the case in type 2, and their reduced optical thickness towards the centre, $A_V \lesssim 1.5\, \rm {mag}$, insufficient to obscure at ultraviolet nucleus wavelengths. It is concluded that large-scale, hundred parsecs to kiloparsecs long, dust filaments and lanes, reminiscent of those seen in the Milky Way, are a common ingredient to the central parsec of galaxies. Their optical thickness changes along their structure in type 2 reaching optical depths high enough to obscure the nucleus in full. Their location with respect to the nucleus and increasing gradient in optical depth towards the centre could naturally lead to the canonical type 1/2 AGN classification, making these filaments to play the role of the torus. Dust filaments and lanes show equivalent morphologies in molecular gas. Gas kinematic in the filaments indicates mass inflows at rates ${\lt}1 \, \mathrm{M}_{\odot }~ \mathrm{yr}^{-1}$.
The Sunyaev-Zel'dolvich (SZ) effect is expected to be instrumental in measuring velocities of distant clusters in near future telescope surveys. We simplify the calculation of peculiar velocities of galaxy clusters using deep learning frameworks trained on numerical simulations to avoid the independent estimation of the optical depth. Images of distorted photon backgrounds are generated for idealized observations using one of the largest cosmological hydrodynamical simulations, the Magneticum simulations. The model is tested to determine its ability of estimating peculiar velocities from future kinetic SZ observations under different noise conditions. The deep learning algorithm displays robustness in estimating peculiar velocities from kinetic SZ effect by an improvement in accuracy of about 17 per cent compared to the analytical approach.
In this paper we present the complete expressions of the lepton and neutron electric dipole moments (EDMs) in the Standard Model Effective Field Theory (SMEFT), up to 1-loop and dimension-6 level and including both RG running contributions and finite corrections. The latter play a fundamental role in the cases of operators that do not renormalize the dipoles, but there are also classes of operators for which they provide an important fraction, $10-20\%$, of the total 1-loop contribution, if the new physics scale is around $\Lambda=5$ TeV. We present the full set of bounds on each individual Wilson coefficient contributing to the EDMs using both the current experimental constraints, as well as those from future experiments, which are expected to improve by at least an order of magnitude.
Precision studies at electron-positron colliders with centre-of-mass energies in the charm-tau region and below have strongly contributed to our understanding of light-meson interactions at low energies. We focus on the processes involving two or three light mesons with invariant masses below nucleon-antinucleon threshold. A prominent role is given to the interactions of the nine lightest pseudoscalar mesons (pions, kaons, η, and η′) and the two narrow neutral isoscalar vector mesons ω and ϕ. Experimental methods used to produce the mesons are reviewed as well as theory tools to extract properties of the meson-meson interactions. Examples of recent results from the DA ΦNE, BEPCII, and VEPP-2000 colliders are presented. In the outlook we briefly discuss prospects for further studies at future super-charm-tau factories.
We present the integrated three-point shear correlation function iζ± - a higher order statistic of the cosmic shear field - which can be directly estimated in wide-area weak lensing surveys without measuring the full three-point shear correlation function, making this a practical and complementary tool to two-point statistics for weak lensing cosmology. We define it as the one-point aperture mass statistic Map measured at different locations on the shear field correlated with the corresponding local two-point shear correlation function ξ±. Building upon existing work on the integrated bispectrum of the weak lensing convergence field, we present a theoretical framework for computing the integrated three-point function in real space for any projected field within the flat-sky approximation and apply it to cosmic shear. Using analytical formulae for the non-linear matter power spectrum and bispectrum, we model iζ± and validate it on N-body simulations within the uncertainties expected from the sixth year cosmic shear data of the Dark Energy Survey. We also explore the Fisher information content of iζ± and perform a joint analysis with ξ± for two tomographic source redshift bins with realistic shape noise to analyse its power in constraining cosmological parameters. We find that the joint analysis of ξ± and iζ± has the potential to considerably improve parameter constraints from ξ± alone, and can be particularly useful in improving the figure of merit of the dynamical dark energy equation of state parameters from cosmic shear data.
We present observations of a giant Lyα blob (LAB) in the SSA22 protocluster at z = 3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array. Dust continuum, along with [C II] 158 μm and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z = 3.0987-3.1016 in the surroundings are identified in [C II] and dust continuum emission, with two of them potential companions or tidal structures associated with the most massive galaxies. Spatially resolved [C II] and infrared luminosity ratios for the widely distributed media (L[Cɪɪ]/LIR ≍ 10-2-10-3) suggest that the observed extended interstellar media are likely to have originated from star formation activity and the contribution from shocked gas is probably not dominant. LAB1 is found to harbor a total molecular gas mass Mmol = (8.7 ± 2.0) × 1010 M⊙, concentrated in the core region of the Lyα-emitting area. While (primarily obscured) star formation activity in the LAB1 core is one of the most plausible power sources for the Lyα emission, multiple major mergers found in the core may also play a role in making LAB1 exceptionally bright and extended in Lyα as a result of cooling radiation induced by gravitational interactions.
We test the adequacy of ultraviolet (UV) spectra for characterizing the outer structure of Type Ia supernova (SN) ejecta. For this purpose, we perform spectroscopic analysis for ASASSN-14lp, a normal SN Ia showing low continuum in the mid-UV regime. To explain the strong UV suppression, two possible origins have been investigated by mapping the chemical profiles over a significant part of their ejecta. We fit the spectral time series with mid-UV coverage obtained before and around maximum light by HST, supplemented with ground-based optical observations for the earliest epochs. The synthetic spectra are calculated with the one-dimensional MC radiative transfer code TARDIS from self-consistent ejecta models. Among several physical parameters, we constrain the abundance profiles of nine chemical elements. We find that a distribution of 56Ni (and other iron-group elements) that extends towards the highest velocities reproduces the observed UV flux well. The presence of radioactive material in the outer layers of the ejecta, if confirmed, implies strong constraints on the possible explosion scenarios. We investigate the impact of the inferred 56Ni distribution on the early light curves with the radiative transfer code TURTLS, and confront the results with the observed light curves of ASASSN-14lp. The inferred abundances are not in conflict with the observed photometry. We also test whether the UV suppression can be reproduced if the radiation at the photosphere is significantly lower in the UV regime than the pure Planck function. In this case, solar metallicity might be sufficient enough at the highest velocities to reproduce the UV suppression.
We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$ which, at 4.8 per cent precision, provides one of the tightest constraints on S8 from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four w cold dark matter parameters (Ωm, σ8, h, w0). We also determine the joint covariance matrix with over 1000 realizations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalizing over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3 per cent precision on S8 for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.
We investigate strongly gravitationally lensed type II supernovae (LSNe II) for time-delay cosmography, incorporating microlensing effects; this expands on previous microlensing studies of type Ia supernovae (SNe Ia). We use the radiative-transfer code TARDIS to recreate five spectra of the prototypical SN 1999em at different times within the plateau phase of the light curve. The microlensing-induced deformations of the spectra and light curves are calculated by placing the SN into magnification maps generated with the code GERLUMPH. We study the impact of microlensing on the color curves and find that there is no strong influence on them during the investigated time interval of the plateau phase. The color curves are only weakly affected by microlensing due to the almost achromatic behavior of the intensity profiles. However, the lack of nonlinear structure in the color curves during the plateau phase of type II-plateau supernovae makes time-delay measurements more challenging compared to SN Ia color curves, given the possible presence of differential dust extinction. Therefore, we further investigate SN phase inference through spectral absorption lines under the influence of microlensing and Gaussian noise. As the spectral features shift to longer wavelengths with progressing time after explosion, the measured wavelength of a specific absorption line provides information on the epoch of the SN. The comparison between retrieved epochs of two observed lensing images then gives the time delay of the images. We find that the phase retrieval method that uses spectral features yields accurate delays with uncertainties of ≲2 days, making it a promising approach.
Time-delay strong lensing (TDSL) is a powerful probe of the current expansion rate of the Universe. However, in light of the discrepancies between early and late-time cosmological studies, efforts revolve around the characterisation of systematic uncertainties in the methods. Here, we focus on the mass-sheet degeneracy (MSD), which is considered a significant source of systematics in TDSL, and aim to assess the constraining power provided by IFU stellar kinematics. We approximate the MSD with a cored, two-parameter extension to the lensing mass profiles (with core radius $r_{\rm c}$ and mass-sheet parameter $\lambda_{\rm int}$). In addition, we utilise mock IFU stellar kinematics of time-delay strong lenses, given the prospects of obtaining such data with JWST. We construct joint strong lensing and stellar dynamical models, where the time delays, mock imaging and IFU observations are used to constrain the mass profile of lens galaxies, and yield joint constraints on the time-delay distance ($D_{\Delta t}$) and angular diameter distance ($D_{\rm d}$) to the lens. We find that mock JWST-like stellar kinematics constrain the internal mass sheet and limit its contribution to the uncertainties of $D_{\Delta t}$ and $D_{\rm d}$, each at the < 4% level, without assumptions on the background cosmological model. These distance constraints would translate to a < 4% precision measurement on $H_{\rm 0}$ in flat $\Lambda CDM$ for a single lens. Our study shows that IFU stellar kinematics of time-delay strong lenses will be key in lifting the MSD on a per lens basis, assuming reasonable core sizes. However, even in the limit of infinite $r_{\rm c}$, where $D_{\Delta t}$ is degenerate with $\lambda_{\rm int}$, stellar kinematics of the deflector, time delays and imaging data will provide powerful constraints on $D_{\rm d}$, which becomes the dominant source of information in the cosmological inference.
The impact of new and highly precise neutron β decay data is reviewed. We focus on recent results from neutron lifetime, β asymmetry, and electron-neutrino correlation experiments. From these results, weak interaction parameters are extracted with unprecedented precision, which is possible also because of progress in effective field theory and lattice QCD. Limits on New Physics beyond the Standard Model derived from neutron decay data are sharper than those derived from high-energy experiments, except for processes involving right-handed neutrinos.
The strong interaction among hadrons has been measured in the past by scattering experiments. Although this technique has been extremely successful in providing information about the nucleon-nucleon and pion-nucleon interactions, when unstable hadrons are considered the experiments become more challenging. In the last few years, the analysis of correlations in the momentum space for pairs of stable and unstable hadrons measured in pp and p+Pb collisions by the ALICE Collaboration at the LHC has provided a new method to investigate the strong interaction among hadrons. In this article, we review the numerous results recently achieved for hyperon-nucleon, hyperon-hyperon, and kaon-nucleon pairs, which show that this new method opens the possibility of measuring the residual strong interaction of any hadron pair.
SU(2) gauge fields coupled to an axion field can acquire an isotropic background solution during inflation. We study homogeneous but anisotropic inflationary solutions in the presence of such (massless) gauge fields. A gauge field in the cosmological background may pose a threat to spatial isotropy. We show, however, that such models generally isotropize in Bianchi type-I geometry, and the isotropic solution is the attractor. Restricting the setup by adding an axial symmetry, we revisited the numerical analysis presented in [1]. We find that the reported numerical breakdown in the previous analysis is an artifact of parametrization singularity. We use a new parametrization that is well-defined all over the phase space. We show that the system respects the cosmic no-hair conjecture and the anisotropies always dilute away within a few e-folds.
Ghost-free bimetric theory describes two nonlinearly interacting spin-2 fields, one massive and one massless, thus extending general relativity. We confront bimetric theory with observations of Supernovae type 1a, Baryon Acoustic Oscillations and the Cosmic Microwave Background in a statistical analysis, utilising the recently proposed physical parametrisation. This directly constrains the physical parameters of the theory, such as the mass of the spin-2 field and its coupling to matter. We find that all models under consideration are in agreement with the data. Next, we compare these results to bounds from local tests of gravity. Our analysis reveals that all two- and three parameter models are observationally consistent with both cosmological and local tests of gravity. The minimal bimetric model (only β1) is ruled out by our combined analysis.
The statistical models used to derive the results of experimental analyses are of incredible scientific value and are essential information for analysis preservation and reuse. In this paper, we make the scientific case for systematically publishing the full statistical models and discuss the technical developments that make this practical. By means of a variety of physics cases -- including parton distribution functions, Higgs boson measurements, effective field theory interpretations, direct searches for new physics, heavy flavor physics, direct dark matter detection, world averages, and beyond the Standard Model global fits -- we illustrate how detailed information on the statistical modelling can enhance the short- and long-term impact of experimental results.
We provide a simple computation in order to estimate the probability of a given hierarchy between two scales. In particular, we work in a model provided with a gauge symmetry, with two scalar doublets. We start from a scale-invariant classical Lagrangian, but by taking into account the Coleman-Weinberg mechanism, we obtain masses for the gauge bosons and the scalars. This approach typically provides a light (L) and a heavy (H) sector related to the two different vacuum expectation values of the two scalars. We compute the size of the hyper-volume of the parameter space of the model associated with an interval of mass ratios between these two sectors. We define the probability as proportional to this size and conclude that probabilities of very large hierarchies are not negligible in the type of models studied in this letter.
A fraction of the dark matter in the solar neighborhood might be composed of non-galactic particles with speeds larger than the escape velocity of the Milky Way. The non-galactic dark matter flux would enhance the sensitivity of direct detection experiments, due to the larger momentum transfer to the target. In this note, we calculate the impact of the dark matter flux from the Local Group and the Virgo Supercluster diffuse components in nuclear and electron recoil experiments. The enhancement in the signal rate can be very significant, especially for experiments searching for dark matter induced electron recoils.
We consider gauged linear sigma models with gauge group U(1) that exhibit a geometric as well as a Landau Ginzburg phase. We construct defects that implement the transport of D-branes from the Landau-Ginzburg phase to the geometric phase. Through their fusion with boundary conditions these defects in particular provide functors between the respective D-brane categories. The latter map (equivariant) matrix factorizations to coherent sheaves and can be formulated explicitly in terms of complexes of matrix factorizations.
We find the complete set of conditions satisfied by the forward 2 →2 scattering amplitude in unitary and causal theories. These are based on an infinite set of energy dependent quantities (the arcs) which are dispersively expressed as moments of a positive measure defined at (arbitrarily) higher energies. We identify optimal finite subsets of constraints, suitable to bound effective field theories (EFTs), at any finite order in the energy expansion. At tree level arcs are in a one to one correspondence with Wilson coefficients. We establish under which conditions this approximation applies, identifying seemingly viable EFTs where it never does. In all cases, we discuss the range of validity in both energy and couplings, where the latter has to satisfy two-sided bounds. We also extend our results to the case of small but finite t . A consequence of our study is that EFTs in which the scattering amplitude in some regime grows in energy faster than E6 cannot be UV completed.
The differential cross section for the quasi-free photoproduction reaction $\gamma n\rightarrow K^0\Sigma^0$ was measured at BGOOD at ELSA from threshold to a center-of-mass energy of 2400 MeV. An increase in the cross section is observed at forward angles above 2000 MeV. The available statistics prevent an accurate description of this behavior, however it appears consistent with models describing a resonance of dynamically generated vector meson-baryon states, where an equivalent model predicted the $P_C$ states observed at LHCb. If proven correct, this could indicate parallels between the strange and charmed quark sectors.
We set out to determine stellar labels from low-resolution survey spectra of hot, OBA stars with effective temperature (Teff) higher than 7500K. This fills a gap in the scientific analysis of large spectroscopic stellar surveys such as LAMOST, which offers spectra for millions of stars at R=1800. We first explore the theoretical information content of such spectra for determining stellar labels, via the Cramér-Rao bound. We show that in the limit of perfect model spectra and observed spectra with S/N of 100, precise estimates are possible for a wide range of stellar labels: not only the effective temperature Teff, surface gravity logg, and projected rotation velocity vsini, but also the micro-turbulence velocity, Helium abundance and the elemental abundances [C/H], [N/H], [O/H], [Si/H], [S/H], and [Fe/H]. Our analysis illustrates that the temperature regime of around 9500K is challenging, as the dominant Balmer and Paschen line strength vary little with Teff. We implement the simultaneous fitting of these 11 stellar labels to LAMOST hot-star spectra using the Payne approach, drawing on Kurucz's ATLAS12/SYNTHE LTE spectra as the underlying models. We then obtain stellar parameter estimates for a sample of about 330,000 hot stars with LAMOST spectra, an increase by about two orders of magnitude in sample size. Among them, about 260,000 have good Gaia parallaxes (S/N>5), and more than 95 percent of them are luminous stars, mostly on the main sequence; the rest reflects lower luminosity evolved stars, such as hot subdwarfs and white dwarfs. We show that the fidelity of the abundance estimates is limited by the systematics of the underlying models, as they do not account for NLTE effects. Finally, we show the detailed distribution of vsini of stars with 8000-15,000K, illustrating that it extends to a sharp cut-off at the critical rotation velocity, across a wide range of temperatures.
Using 324 numerically modelled galaxy clusters as provided by THE THREE HUNDRED project, we study the evolution of the kinematic properties of the stellar component of haloes on first infall. We selected objects with Mstar > 5 × 1010 h−1 M⊙ within 3R200 of the main cluster halo at z = 0 and followed their progenitors. We find that although haloes are stripped of their dark matter and gas after entering the main cluster halo, there is practically no change in their stellar kinematics. For the vast majority of our `galaxies' - defined as the central stellar component found within the haloes that form our sample - their kinematic properties, as described by the fraction of ordered rotation, and their position in the specific stellar angular momentum−stellar mass plane jstar − Mstar are mostly unchanged by the influence of the central host cluster. However, for a small number of infalling galaxies, stellar mergers and encounters with remnant stellar cores close to the centre of the main cluster, particularly during pericentre passage, are able to spin up their stellar component by z = 0.
Catalytic nucleic acids, such as ribozymes, are central to a variety of origin-of-life scenarios. Typically, they require elevated magnesium concentrations for folding and activity, but their function can be inhibited by high concentrations of monovalent salts. Here we show that geologically plausible high-sodium, low-magnesium solutions derived from leaching basalt (rock and remelted glass) inhibit ribozyme catalysis, but that this activity can be rescued by selective magnesium up-concentration by heat flow across rock fissures. In contrast to up-concentration by dehydration or freezing, this system is so far from equilibrium that it can actively alter the Mg:Na salt ratio to an extent that enables key ribozyme activities, such as self-replication and RNA extension, in otherwise challenging solution conditions. The principle demonstrated here is applicable to a broad range of salt concentrations and compositions, and, as such, highly relevant to various origin-of-life scenarios.
We investigate the impact of different assumptions in the modeling of one-loop galaxy bias on the recovery of cosmological parameters, as a follow-up of the analysis done in the first paper of the series at fixed cosmology. To carry out these tests we focus on the real-space galaxy-power spectrum from a set of three different synthetic galaxy samples whose clustering properties are meant to match the ones of the CMASS and LOWZ catalogs of BOSS and the SDSS Main Galaxy Sample. We investigate the relevance of allowing for either short range nonlocality or scale-dependent stochasticity by fitting the real-space galaxy autopower spectrum or the combination of galaxy-galaxy and galaxy-matter power spectrum. From a comparison among the goodness of fit (χ2 ), unbiasedness of cosmological parameters (FoB), and figure of merit (FoM) of the model, we find that a simple four-parameter model (linear, quadratic, cubic nonlocal bias, and constant shot noise) with fixed quadratic tidal bias provides a robust modeling choice for the autopower spectrum of the three galaxy samples, up to kmax=0.3 h Mpc-1 and for an effective volume of 6 h-3 Gpc3. Instead, a joint analysis of the two observables fails at larger scales, and a model extension with either higher derivatives or scale-dependent shot noise is necessary to reach a similar kmax, with the latter providing the most accurate and stable results. Throughout the majority of the paper, we fix the description of the nonlinear matter evolution using a hybrid perturbative-N-body approach, RESPRESSO, that was found in the first paper to be the closest performing to the measured matter spectrum. We also test the impact of different modeling assumptions based on perturbative approaches, such as galilean-invariant Renormalised Perturbation Theory (gRPT) and effective field theory (EFT). In all cases, we find the inclusion of scale-dependent shot noise to increase the range of validity of the model in terms of FoB and χ2. Interestingly, these model extensions with additional free parameters do not necessarily lead to an increase in the maximally achievable FoM for the cosmological parameters (h ,Ωch2,As), which are generally consistent with those of the simpler model at smaller kmax.
We study the transition widths of ϒ (10753 ) and ϒ (11020 ) into standard bottomonium under the hypothesis that they correspond to the two lowest laying 1-- hybrid bottomonium states. We employ weakly coupled potential NRQCD an effective field theory incorporating the heavy-quark and multipole expansions. We consider the transitions generated by the leading order and next-to-leading order singlet-octet operators. In the multipole expansion the heavy-quark matrix elements factorize from the production of light-quark mesons by gluonic operators. For the leading order operator we compute the widths with a single π0, η or η' in the final state and for the next-to-leading operator for π+π- or K+K-. The hadronization of the gluonic operators is obtained, in the first case, from the axial anomaly and a standard π0-η -η' mixing scheme and, in the second case, we employ a coupled-channel dispersive representation matched to chiral perturbation theory for both the S - and D -wave pieces of the gluonic operator. We compare with experimental values and semi-inclusive widths. Our results strongly suggest that ϒ (11020 ) is indeed a hybrid bottomonium state.
The BV-Laplacian $\Delta$ in quantum field theory is singular, by construction, but can be regularized by deforming the classical BV-action. Taking inspiration from string theory we describe a non-local deformation of the latter by adding stubs to the interaction vertices while keeping classical BV-invariance manifest. This is achieved using a version of homotopy transfer resulting in a non-polynomial action for which the quantum master equation is now well defined and will be satisfied by adding additional vertices at loop level. The latter can be defined with the help of standard regularization schemes and is independent of the definition of $\Delta$. In particular, the determination of anomalies reduces to the standard text-book calculation. Finally, we describe how the deformed (quantum) action can be obtained as a canonical transformation. As an example, we illustrate this procedure for quantum electrodynamics.
Manual fits to spectral times series of Type Ia supernovae have provided a method of reconstructing the explosion from a parametric model but due to lack of information about model uncertainties or parameter degeneracies direct comparison between theory and observation is difficult. In order to mitigate this important problem we present a new way to probabilistically reconstruct the outer ejecta of the normal Type Ia supernova SN 2002bo. A single epoch spectrum, taken 10 days before maximum light, is fit by a 13-parameter model describing the elemental composition of the ejecta and the explosion physics (density, temperature, velocity, and explosion epoch). Model evaluation is performed through the application of a novel rapid spectral synthesis technique in which the radiative transfer code, TARDIS, is accelerated by a machine-learning framework. Analysis of the posterior distribution reveals a complex and degenerate parameter space and allows direct comparison to various hydrodynamic models. Our analysis favors detonation over deflagration scenarios and we find that our technique offers a novel way to compare simulation to observation.
The SuperKEKB accelerator in Tsukuba, Japan is providing e$^+$e$^-$ beams for the Belle II experiment since March 2019. To deal with the aimed peak luminosity being forty times higher than the one recorded at Belle, a pixel detector based on DEPFET technology has been installed. It features a long integration time of 20 $\mu$s resulting in an expected data rate of 20 GByte/s (160 GBit/s) at a maximum occupancy of 3 %. To deal with this high amount of data, the data handling hub (DHH) has been developed. It contains all necessary functionality for the control and readout of the detector. In this paper we describe the architecture and features of the DHH system. Further we will show the key performance characteristics after one year of operation.
A free-floating planet is a planetary-mass object that orbits around a non-stellar massive object (e.g. a brown dwarf) or around the Galactic Center. The presence of exomoons orbiting free-floating planets has been theoretically predicted by several models. Under specific conditions, these moons are able to retain an atmosphere capable of ensuring the long-term thermal stability of liquid water on their surface. We model this environment with a one-dimensional radiative-convective code coupled to a gas-phase chemical network including cosmic rays and ion-neutral reactions. We find that, under specific conditions and assuming stable orbital parameters over time, liquid water can be formed on the surface of the exomoon. The final amount of water for an Earth-mass exomonoon is smaller than the amount of water in Earth oceans, but enough to host the potential development of primordial life. The chemical equilibrium time-scale is controlled by cosmic rays, the main ionization driver in our model of the exomoon atmosphere.
We present an empirical model for the number of globular clusters (GCs) in galaxies based on recent data showing a tight relationship between dark matter halo virial masses and GC numbers. While a simple base model forming GCs in low-mass haloes reproduces this relation, we show that a second formation pathway for GCs is needed to account for observed younger GC populations. We confirm previous works that reported the observed linear correlation as being a consequence of hierarchical merging and its insensitivity to the exact GC formation processes at higher virial masses, even for a dual formation scenario. We find that the scatter of the linear relation is strongly correlated with the relative amount of smooth accretion: the more dark matter is smoothly accreted, the fewer GCs a halo has compared to other haloes of the same mass. This scatter is smaller than that introduced by halo mass measurements, indicating that the number of GCs in a galaxy is a good tracer for its dark matter mass. Smooth accretion is also the reason for a lower average dark matter mass per GC in low-mass haloes. Finally, we successfully reproduce the observed general trend of GCs being old and the tendency of more massive haloes hosting older GC systems. Including the second GC formation mechanism through gas-rich mergers leads to a more realistic variety of GC age distributions and also introduces an age inversion in the halo virial mass range log Mvir/M⊙ = 11-13.
We present a set of nonlocal thermodynamic equilibrium steady-state calculations of radiative transfer for one-year-old Type II supernovae (SNe) starting from state-of-the-art explosion models computed with detailed nucleosynthesis. This grid covers single-star progenitors with initial masses between 9 and 29 M⊙, all evolved with the code KEPLER at solar metallicity and ignoring rotation. The [O I] λλ 6300, 6364 line flux generally grows with progenitor mass, and Hα exhibits an equally strong and opposite trend. The [Ca II] λλ 7291, 7323 strength increases at low 56Ni mass, at low explosion energy, or with clumping. This Ca II doublet, which forms primarily in the explosively produced Si/S zones, depends little on the progenitor mass but may strengthen if Ca+ dominates in the H-rich emitting zones or if Ca is abundant in the O-rich zones. Indeed, Si-O shell merging prior to core collapse may boost the Ca II doublet at the expense of the O I doublet, and may thus mimic the metal line strengths of a lower-mass progenitor. We find that the 56Ni bubble effect has a weak impact, probably because it is too weak to induce much of an ionization shift in the various emitting zones. Our simulations compare favorably to observed SNe II, including SN 2008bk (e.g., the 9 M⊙ model), SN 2012aw (12 M⊙ model), SN 1987A (15 M⊙ model), or SN 2015bs (25 M⊙ model with no Si-O shell merging). SNe II with narrow lines and a low 56Ni mass are well matched by the weak explosion of 9-11 M⊙ progenitors. The nebular-phase spectra of standard SNe II can be explained with progenitors in the mass range 12-15 M⊙, with one notable exception for SN 2015bs. In the intermediate mass range, these mass estimates may increase by a few M⊙, with allowance for clumping of the O-rich material or CO molecular cooling.
The temperatures of red supergiants (RSGs) are expected to depend on metallicity (Z) in such a way that lower Z RSGs are warmer. In this work, we investigate the Z-dependence of the Hayashi limit by analysing RSGs in the low-Z galaxy Wolf-Lundmark-Mellote, and compare with the RSGs in the higher Z environments of the Small Magellanic Cloud and Large Magellanic Cloud. We determine the effective temperature (Teff) of each star by fitting their spectral energy distributions, as observed by VLT + SHOOTER, with MARCS model atmospheres. We find average temperatures of $T_{\textrm {eff}_{\textrm {WLM}}}=4400\pm 202$ K, $T_{\textrm {eff}_{\textrm {SMC}}}=4130\pm 103$ K, and $T_{\textrm {eff}_{\textrm {LMC}}}=4140\pm 148$ K. From population synthesis analysis, we find that although the Geneva evolutionary models reproduce this trend qualitatively, the RSGs in these models are systematically too cool. We speculate that our results can be explained by the inapplicability of the standard solar mixing length to RSGs.
We present non-radiative, cosmological zoom-simulations of galaxy cluster formation with magnetic fields and (anisotropic) thermal conduction of one very massive galaxy cluster with a mass at redshift zero that corresponds to $M_\mathrm{vir} \sim 2 \times 10^{15} M_{\odot}$. We run the cluster on three resolution levels (1X, 10X, 25X), starting with an effective mass resolution of $2 \times 10^8M_{\odot}$, subsequently increasing the particle number to reach $4 \times 10^6M_{\odot}$. The maximum spatial resolution obtained in the simulations is limited by the gravitational softening reaching $\epsilon=1.0$ kpc at the highest resolution level, allowing to resolve the hierarchical assembly of the structures in very fine detail. All simulations presented, have been carried out with the SPMHD-code Gadget-3 with a heavily updated SPMHD prescription. The primary focus is to investigate magnetic field amplification in the Intracluster Medium (ICM). We show that the main amplification mechanism is the small scale-turbulent-dynamo in the limit of reconnection diffusion. In our two highest resolution models we start to resolve the magnetic field amplification driven by this process and we explicitly quantify this with the magnetic power-spectra and the magnetic tension that limits the bending of the magnetic field lines consistent with dynamo theory. Furthermore, we investigate the $\nabla \cdot \mathbf{B}=0$ constraint within our simulations and show that we achieve comparable results to state-of-the-art AMR or moving-mesh techniques, used in codes such as Enzo and Arepo. Our results show for the first time in a fully cosmological simulation of a galaxy cluster that dynamo action can be resolved in the framework of a modern Lagrangian magnetohydrodynamic (MHD) method, a study that is currently missing in the literature.
Within the transport model evaluation project (TMEP) of simulations for heavy-ion collisions, the mean-field response is examined here. Specifically, zero-sound propagation is considered for neutron-proton symmetric matter enclosed in a periodic box, at zero temperature and around normal density. The results of several transport codes belonging to two families (BUU-like and QMD-like) are compared among each other and to exact calculations. For BUU-like codes, employing the test particle method, the results depend on the combination of the number of test particles and the spread of the profile functions that weight integration over space. These parameters can be properly adapted to give a good reproduction of the analytical zero-sound features. QMD-like codes, using molecular dynamics methods, are characterized by large damping effects, attributable to the fluctuations inherent in their phase-space representation. Moreover, for a given nuclear effective interaction, they generally lead to slower density oscillations, as compared to BUU-like codes. The latter problem is mitigated in the more recent lattice formulation of some of the QMD codes. The significance of these results for the description of real heavy-ion collisions is discussed.
The baryon acoustic oscillation feature can be used as a standard cosmological ruler. In practice, for subpercent level accuracy on the distance scale, it must be standardized. The physical reason why is understood, so we use this to develop an algorithm which improves the estimated scale. The algorithm exploits the fact that, over the range of scales where the initial correlation function is well fit by a polynomial, the leading order effects which distort the length of the ruler can be accounted for analytically. Tests of the method in numerical simulations show that it provides simple and fast reconstruction of the full shape of the BAO feature, as well as subpercent determination of the linear point in the correlation function of biased tracers with minimal assumptions about the underlying cosmological model or the nature of the observed tracers. Our results also suggest that, for least squares estimators of the correlation function, half-integer generalized Laguerre functions are a particularly useful choice.
We present Large Millimeter Telescope (LMT)/AzTEC 1.1 mm observations of ~100 luminous high-redshift dusty star-forming galaxy candidates from the $\sim 600\,$ sq.deg Herschel-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500\, \mu \rm m}=35-80$ mJy. With an effective $\theta _{\rm FWHM}\approx 9.5\,$arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with signal-to-noise ratio ≥4 members. The fraction of multiple systems increases to ~23 per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended Herschel photometry, we derive photometric redshifts, infrared luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar (zmed ≍ 3.6), the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources, ~85 per cent lie at zphot > 3 while ~33 per cent are at zphot > 4. This corresponds to a lower limit on the space density of ultrared sources at 4 < z < 6 of $\sim 3\times 10^{-7}\, \textrm {Mpc}^{-3}$ with a contribution to the obscured star formation of $\gtrsim 8\times 10^{-4}\, \textrm {M}_\odot \, \textrm {yr}^{-1} \, \textrm {Mpc}^{-3}$. Some of the multiple systems have members with photometric redshifts consistent among them suggesting possible physical associations. Given their angular separations, these systems are most likely galaxy over-densities and/or early-stage pre-coalescence mergers. Finally, we present 3 mm LMT/RSR spectroscopic redshifts of six red-Herschel galaxies at zspec = 3.85-6.03, two of them (at z ~ 4.7) representing new redshift confirmations. Here, we release the AzTEC and deblended Herschel photometry as well as catalogues of the most promising interacting systems and z > 4 galaxies.
A search for charginos and neutralinos at the Large Hadron Collider using fully hadronic final states and missing transverse momentum is reported. Pair-produced charginos or neutralinos are explored, each decaying into a high-<math display="inline"><msub><mi>p</mi><mi mathvariant="normal">T</mi></msub></math> Standard Model weak boson. Fully hadronic final states are studied to exploit the advantage of the large branching ratio, and the efficient rejection of backgrounds by identifying the high-<math display="inline"><msub><mi>p</mi><mi mathvariant="normal">T</mi></msub></math> bosons using large-radius jets and jet substructure information. An integrated luminosity of <math display="inline"><mrow><mn>139</mn><mtext> </mtext><mtext> </mtext><msup><mrow><mi>fb</mi></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow></math> of proton-proton collision data collected by the ATLAS detector at a center-of-mass energy of 13 TeV is used. No significant excess is found beyond the Standard Model expectation. Exclusion limits at the 95% confidence level are set on wino or higgsino production with various assumptions about the decay branching ratios and the type of lightest supersymmetric particle. A wino (higgsino) mass up to 1060 (900) GeV is excluded when the lightest supersymmetry particle mass is below 400 (240) GeV and the mass splitting is larger than 400 (450) GeV. The sensitivity to high-mass winos and higgsinos is significantly extended relative to previous LHC searches using other final states.
In this paper we investigate the potential of current and upcoming cosmological surveys to constrain the mass and abundance of ultra-light axion (ULA) cosmologies with galaxy cluster number counts. ULAs, sometimes also referred to as Fuzzy Dark Matter, are well-motivated in many theories beyond the Standard Model and could potentially solve the ΛCDM small-scale crisis. Galaxy cluster counts provide a robust probe of the formation of structures in the Universe. Their distribution in mass and redshift is strongly sensitive to the underlying linear matter perturbations. In this forecast paper we explore two scenarios, firstly an exclusion limit on axion mass given a no-axion model and secondly constraints on an axion model. With this we obtain lower limits on the ULA mass on the order of ma ≳ 10-24 eV. However, this result depends heavily on the mass of the smallest reliably observable clusters for a given survey. Cluster counts, like many other cosmological probes, display an approximate degeneracy in the ULA mass vs. abundance parameter space, which is dependent on the characteristics of the probe. These degeneracies are different for other cosmological probes. Hence galaxy cluster counts might provide a complementary window on the properties of ultra-light axions.
We present non-radiative, cosmological zoom-in simulations of galaxy-cluster formation with magnetic fields and (anisotropic) thermal conduction of one massive galaxy cluster with M
$_{vir}$ ∼ 2 × 10$^{15}$
M
$_{⊙}$ at z ∼ 0. We run the cluster on three resolution levels (1×, 10×, 25×), starting with an effective mass resolution of 2 × 10$^{8}$
M
$_{⊙}$, subsequently increasing the particle number to reach 4 × 10$^{6}$
M
$_{⊙}$. The maximum spatial resolution obtained in the simulations is limited by the gravitational softening reaching ϵ = 1.0 kpc at the highest resolution level, allowing one to resolve the hierarchical assembly of the structures in fine detail. All simulations presented are carried out with the SPMHD code gadget3 with an updated SPMHD prescription. The primary focus of this paper is to investigate magnetic field amplification in the intracluster medium. We show that the main amplification mechanism is the small-scale turbulent dynamo in the limit of reconnection diffusion. In our two highest resolution models we start to resolve the magnetic field amplification driven by the dynamo and we explicitly quantify this with the magnetic power spectra and the curvature of the magnetic field lines, consistent with dynamo theory. Furthermore, we investigate the ∇ ·
B
= 0 constraint within our simulations and show that we achieve comparable results to state-of-the-art AMR or moving-mesh techniques, used in codes such as enzo and arepo. Our results show for the first time in a cosmological simulation of a galaxy cluster that dynamo action can be resolved with modern numerical Lagrangian magnetohydrodynamic methods, a study that is currently missing in the literature.
Time-delay cosmography with gravitationally lensed quasars plays an important role in anchoring the absolute distance scale and hence measuring the Hubble constant, H0, independent of traditional distance ladder methodology. A current potential limitation of time-delay distance measurements is the mass-sheet transformation (MST), which leaves the lensed imaging unchanged but changes the distance measurements and the derived value of H0. In this work we show that the standard method of addressing the MST in time-delay cosmography, through a combination of high-resolution imaging and the measurement of the stellar velocity dispersion of the lensing galaxy, depends on the assumption that the ratio, Ds/Dds, of angular diameter distances to the background quasar and between the lensing galaxy and the quasar can be constrained. This is typically achieved through the assumption of a particular cosmological model. Previous work (TDCOSMO IV) addressed the mass-sheet degeneracy and derived H0 under the assumption of the ΛCDM model. In this paper we show that the mass-sheet degeneracy can be broken without relying on a specific cosmological model by combining lensing with relative distance indicators such as supernovae Type Ia and baryon acoustic oscillations, which constrain the shape of the expansion history and hence Ds/Dds. With this approach, we demonstrate that the mass-sheet degeneracy can be constrained in a cosmological model-independent way. Hence model-independent distance measurements in time-delay cosmography under MSTs can be obtained.
We use hydrodynamical separate universe simulations with the IllustrisTNG model to predict the local primordial non-Gaussianity (PNG) bias parameters b
$_{ϕ}$ and b
$_{ϕδ}$, which enter at leading order in the galaxy power spectrum and bispectrum. This is the first time that b
$_{ϕδ}$ is measured from either gravity-only or galaxy formation simulations. For dark matter halos, the popular assumption of universality overpredicts the b
$_{ϕδ}$(b
$_{1}$) relation in the range 1 ≲ b
$_{1}$ ≲ 3 by up to Δ b
$_{ϕδ}$ ∼ 3 (b
$_{1}$ is the linear density bias). The adequacy of the universality relation is worse for the simulated galaxies, with the relations b
$_{ϕ}$(b
$_{1}$) and b
$_{ϕδ}$(b
$_{1}$) being generically redshift-dependent and very sensitive to how galaxies are selected (we test total, stellar and black hole mass, black hole mass accretion rate and color). The uncertainties on b
$_{ϕ}$ and b
$_{ϕδ}$ have a direct, often overlooked impact on the constraints of the local PNG parameter f
$_{NL}$, which we study and discuss. For a survey with V = 100 Gpc$^{3}$/h$^{3}$ at z=1, uncertainties Δ b
$_{ϕ}$ ≲ 1 and Δ b
$_{ϕδ}$ ≲ 5 around values close to the fiducial can yield relatively unbiased constraints on f
$_{NL}$ using power spectrum and bispectrum data. We also show why priors on galaxy bias are useful even in analyses that fit for products f
$_{NL}$
b
$_{ϕ}$ and f
$_{NL}$
b
$_{ϕδ}$. The strategies we discuss to deal with galaxy bias uncertainties can be straightforwardly implemented in existing f
$_{NL}$ constraint analyses (we provide fits for some of the bias relations). Our results motivate more works with galaxy formation simulations to refine our understanding of b
$_{ϕ}$ and b
$_{ϕδ}$ towards improved constraints on f
$_{NL}$.
The properties of quasar-host galaxies might be determined by the growth and feedback of their supermassive black holes (SMBHs, 10^8−10 M_⊙). We investigate such connection with a suite of cosmological simulations of massive (halo mass ≈10^12 M_⊙) galaxies at z ≃ 6 that include a detailed subgrid multiphase gas and accretion model. BH seeds of initial mass 10^5 M_⊙ grow mostly by gas accretion, and become SMBH by z = 6 setting on the observed M_BH−M_⋆ relation without the need for a boost factor. Although quasar feedback crucially controls the SMBH growth, its impact on the properties of the host galaxy at z = 6 is negligible. In our model, quasar activity can both quench (via gas heating) or enhance (by interstellar medium overpressurization) star formation. However, we find that the star formation history is insensitive to such modulation as it is largely dominated, at least at z > 6, by cold gas accretion from the environment that cannot be hindered by the quasar energy deposition. Although quasar-driven outflows can achieve velocities |$\gt 1000~\rm km~s^{-1}$|, only ≈4 per cent of the outflowing gas mass can actually escape from the host galaxy. These findings are only loosely constrained by available data, but can guide observational campaigns searching for signatures of quasar feedback in early galaxies.
The emergence of evermore complex entities from prebiotic building blocks is a key aspect of origins of life research. The RNA-world hypothesis posits that RNA oligomers known as ribozymes acted as the first self-replicating entities. However, the mechanisms governing the self-assembly of complex informational polymers from the shortest prebiotic building blocks were unclear. One open issue concerns the relation between concentration and oligonucleotide length, usually assumed to be exponentially decreasing. Here, we show that a competition of timescales in the self-assembly of informational polymers by templated ligation generically leads to nonmonotonic strand-length distributions with two distinct length scales. The first length scale characterizes the onset of a strongly nonequilibrium regime and is visible as a local minimum. Dynamically, this regime is governed by extension cascades, where the elongation of a "primer" with a short building block is more likely than its dehybridization. The second length scale appears as a local concentration maximum and reflects a balance between degradation and dehybridization of completely hybridized double strands in a heterocatalytic extension-reassembly process. Analytical arguments and extensive numerical simulations within a sequence-independent model allowed us to predict and control these emergent length scales. Nonmonotonic strand-length distributions confirming our theory were obtained in thermocycler experiments using random DNA sequences from a binary alphabet. Our work emphasizes the role of structure-forming processes already for the earliest stages of prebiotic evolution. The accumulation of strands with a typical length reveals a possible starting point for higher-order self-organization events that ultimately lead to a self-replicating, evolving system.
Simulating a survey of fluxes and redshifts (distances) from an astrophysical population is a routine task. \texttt{popsynth} provides a generic, object-oriented framework to produce synthetic surveys from various distributions and luminosity functions, apply selection functions to the observed variables and store them in a portable (HDF5) format. Population synthesis routines can be constructed either using classes or from a serializable YAML format allowing flexibility and portability. Users can not only sample the luminosity and distance of the populations, but they can create auxiliary distributions for parameters which can have arbitrarily complex dependencies on one another. Thus, users can simulate complex astrophysical populations which can be used to calibrate analysis frameworks or quickly test ideas.
Molecular outflows contributing to the matter cycle of star-forming galaxies are now observed in small and large systems at low and high redshift. Their physical origin is still unclear. In most theoretical studies, only warm ionized/neutral and hot gas outflowing from the interstellar medium is generated by star formation. We investigate an in situ H2 formation scenario in the outflow using high-resolution simulations, including non-equilibrium chemistry and self-gravity, of turbulent, warm, and atomic clouds with densities 0.1, 0.5, and $1\, \mathrm{cm}^{-3}$ exposed to a magnetized hot wind. For cloud densities $\gtrsim 0.5\, \mathrm{cm}^{-3}$, a magnetized wind triggers H2 formation before cloud dispersal. Up to 3 per cent of the initial cloud mass can become molecular on $\sim \! 10\, \mathrm{Myr}$ time-scales. The effect is stronger for winds with perpendicular B-fields and intermediate density clouds ($n_\mathrm{c}\sim 0.5\, \mathrm{cm}^{-3}$). Here, H2 formation can be boosted by up to one order of magnitude compared to isolated cooling clouds independent of self-gravity. Self-gravity preserves the densest clouds well past their $\sim \! 15\, \mathrm{Myr}$ cloud crushing time-scales. This model could provide a plausible in situ origin for the observed molecular gas. All simulations form warm ionized gas, which represents an important observable phase. The amount of warm ionized gas is almost independent of the cloud density but solely depends on the magnetic field configuration in the wind. For low-density clouds ($0.1\, \mathrm{cm}^{-3}$), up to 60 per cent of the initially atomic cloud mass can become warm and ionized.
Hydroxyl ($\rm OH$) is known to form efficiently in cold gas (T ~ 100 K) along with the molecule $\rm H_2$ and can be used as an efficient tracer of the diffuse molecular gas in the interstellar medium (ISM). Using a simple formalism describing the $\rm H\, I/H_2$ transition and a reduced network of major chemical reactions, we present a semi-analytical prescription to estimate the abundances of O-bearing molecules in the diffuse ISM. We show that predictions based on our prescription are in good agreement with the estimates obtained using the MEUDON PDR code which utilizes the full reaction network. We investigate the dependence of the relative abundances of $\rm OH/H\, I$ and $\rm OH/H_2$ on the variations of physical conditions i.e. the metallicity, number density (n), cosmic ray ionization rate (ζ), and strength of UV field (χ) in the medium. We find that the $\rm OH/H\, I$ abundances observed in the Galactic ISM can be reproduced by models with n ~ 50 cm-3, χ ~ 1 (Mathis field), and ζ ~ 3 × 10-17 s-1, with a variation of about 1 dex allowed around these values. Using the constrained $\rm H_2$ column density distribution function at z ~ 3, we estimate the $\rm OH$ column density distribution function and discuss future prospects with the upcoming large radio absorption line surveys.
To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 1051 erg (=1 bethe=1 B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 M⊙ progenitors. We present the first such model, considering a nonrotating, solar-metallicity 18.88 M⊙ progenitor, whose final 7 minutes of convective oxygen-shell burning were simulated in 3D and showed a violent oxygen-neon shell merger prior to collapse. A large set of 3D SN models was computed with the PROMETHEUS-VERTEX code, whose improved convergence of the two-moment equations with Boltzmann closure allows now to fully exploit the implicit neutrino-transport treatment. Nuclear burning is treated with a 23-species network. We vary the angular grid resolution and consider different nuclear equations of state and muon formation in the proto-neutron star (PNS), which requires six-species transport with coupling of all neutrino flavors across all energy-momentum groups. Elaborate neutrino transport was applied until ~2 s after bounce. In one case, the simulation was continued to >7 s with an approximate treatment of neutrino effects that allows for seamless continuation without transients. A spherically symmetric neutrino-driven wind does not develop. Instead, accretion downflows to the PNS and outflows of neutrino-heated matter establish a monotonic rise of the explosion energy until ~7 s post-bounce, when the outgoing shock reaches ~50,000 km and enters the He layer. The converged value of the explosion energy at infinity (with overburden subtracted) is ~1 B and the ejected 56Ni mass ≲0.087 M⊙, both within a few 10% of the SN 1987A values. The final NS mass and kick are ~1.65 M⊙ and >450 km s-1, respectively.
Blazars are among the most studied sources in high-energy astrophysics as they form the largest fraction of extragalactic gamma-ray sources and are considered prime candidates for being the counterparts of high-energy astrophysical neutrinos. Their reliable identification amid the many faint radio sources is a crucial step for multimessenger counterpart associations. As the astronomical community prepares for the coming of a number of new facilities able to survey the non-thermal sky at unprecedented depths, from radio to gamma-rays, machine-learning techniques for fast and reliable source identification are ever more relevant. The purpose of this work was to develop a deep learning architecture to identify Blazar within a population of active galactic nucleus (AGN) based solely on non-contemporaneous spectral energy distribution information, collected from publicly available multifrequency catalogues. This study uses an unprecedented amount of data, with spectral energy distributions (SEDs) for ≍14 000 sources collected with the Open Universe VOU-Blazars tool. It uses a convolutional long short-term memory neural network purposefully built for the problem of SED classification, which we describe in detail and validate. The network was able to distinguish Blazars from other types of active galactic nuclei (AGNs) to a satisfying degree (achieving a receiver operating characteristic area under curve of 0.98), even when trained on a reduced subset of the whole sample. This initial study does not attempt to classify Blazars among their different sub-classes, or quantify the likelihood of any multifrequency or multimessenger association, but is presented as a step towards these more practically oriented applications.
The structure of protostellar cores can often be approximated by isothermal Bonnor-Ebert spheres (BES), which are stabilized by an external pressure. For the typical pressure of 104kB K cm-3 to 105kB K cm-3 found in molecular clouds, cores with masses below 1.5 M⊙ are stable against gravitational collapse. In this paper, we analyze the efficiency of triggering gravitational collapse with a nearby stellar wind, which represents an interesting scenario for triggering low-mass star formation. We analytically derive a new stability criterion for a BES compressed by a stellar wind, which depends on its initial nondimensional radius ${\xi }_{\max }$ . If the stability limit is violated the wind triggers a core collapse. Otherwise, the core is destroyed by the wind. We estimate its validity range to $2.5\lt {\xi }_{\max }\lt 4.2$ and confirm this in simulations with the SPH-Code GADGET-3. The efficiency of triggering a gravitational collapse strongly decreases for ${\xi }_{\max }\lt 2.5$ since in this case destruction and acceleration of the whole sphere begin to dominate. We were unable to trigger a collapse for ${\xi }_{\max }\lt 2$ , which leads to the conclusion that a stellar wind can move the smallest unstable stellar mass to 0.5 M⊙ and that destabilizing even smaller cores would require external pressure larger than 105kB K cm-3. For ${\xi }_{\max }\gt 4.2$ the expected wind strength according to our criterion is small enough that the compression is slower than the sound speed of the BES and sound waves can be triggered. In this case our criterion somewhat underestimates the onset of collapse and detailed numerical analyses are required.
Orbit superposition models are a non-parametric dynamical modelling technique to determine the mass of a galaxy's central supermassive black hole (SMBH), its stars, or its dark matter halo. One of the main problems is how to decide which model out of a large pool of trial models based on different assumed mass distributions represents the true structure of an observed galaxy best. We show that the traditional approach to judge models solely by their goodness-of-fit can lead to substantial biases in estimated galaxy properties caused by varying model flexibilities. We demonstrate how the flexibility of the models can be estimated using bootstrap iterations and present a model selection framework that removes these biases by taking the variable flexibility into account in the model evaluation. We extend the model selection approach to optimize the degree of regularization directly from the data. Altogether, this leads to a significant improvement of the constraining power of the modelling technique. We show with simulations that one can reconstruct the mass, anisotropy, and viewing angle of an axisymmetric galaxy with a few per cent accuracy from realistic observational data with fully resolved line-of-sight velocity distributions (LOSVDs). In a first application, we reproduce a photometric estimate of the inclination of the disc galaxy NGC 3368 to within 5° accuracy from kinematic data that cover only a few sphere-of-influence radii around the galaxy's SMBH. This demonstrates the constraining power that can be achieved with orbit models based on fully resolved LOSVDs and a model selection framework.
We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on the European Southern Observatory Very Large Telescope. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided Multi-Unit Spectroscopic Explorer (MUSE) observations of fields selected from the Galaxy and Mass Assembly (GAMA) survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive ( ${>}7 × 10^{10} {M}_\odot$ ) central galaxies at $0.25 < z <0.35$ in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics ( $0.6-0.8$ arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, for example, SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at $0.25 < z <0.35$ , as well as hundreds of emission-line sources at $z < 6$ . This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at $z≈ 0.3$ . MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory.
We introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.
Our work presents an independent calibration of the J-region Asymptotic Giant Branch (JAGB) method using Infrared Survey Facility photometric data and a custom luminosity function profile to determine JAGB mean magnitudes for nine galaxies. We determine a mean absolute magnitude of carbon stars of MLMC = -6.212 ± 0.010 (stat.) ±0.030 (syst.) mag. We then use near-infrared photometry of a number of nearby galaxies, originally obtained by our group to determine their distances from Cepheids using the Leavitt law, in order to independently determine their distances with the JAGB method. We compare the JAGB distances obtained in this work with the Cepheid distances resulting from the same photometry and find very good agreement between the results from the two methods. The mean difference is 0.01 mag with an rms scatter of 0.06 mag after taking into account seven out of the eight analyzed galaxies that had their distances determined using Cepheids. The very accurate distance to the Small Magellanic Cloud based on detached eclipsing binaries is also in very good agreement with the distance obtained from carbon stars.
Self-interacting dark matter (SIDM) models have the potential to solve the small-scale problems that arise in the cold dark matter paradigm. Simulations are a powerful tool for studying SIDM in the context of astrophysics, but it is numerically challenging to study differential cross-sections that favour small-angle scattering (as in light-mediator models). Here, we present a novel approach to model frequent scattering based on an effective drag force, which we have implemented into the N-body code GADGET-3. In a range of test problems, we demonstrate that our implementation accurately models frequent scattering. Our implementation can be used to study differences between SIDM models that predict rare and frequent scattering. We simulate core formation in isolated dark matter haloes, as well as major mergers of galaxy clusters and find that SIDM models with rare and frequent interactions make different predictions. In particular, frequent interactions are able to produce larger offsets between the distribution of galaxies and dark matter in equal-mass mergers.
We propose the construction of LEGEND-1000, the ton-scale Large Enriched Germanium Experiment for Neutrinoless $\beta \beta$ Decay. This international experiment is designed to answer one of the highest priority questions in fundamental physics. It consists of 1000 kg of Ge detectors enriched to more than 90% in the $^{76}$Ge isotope operated in a liquid argon active shield at a deep underground laboratory. By combining the lowest background levels with the best energy resolution in the field, LEGEND-1000 will perform a quasi-background-free search and can make an unambiguous discovery of neutrinoless double-beta decay with just a handful of counts at the decay $Q$ value. The experiment is designed to probe this decay with a 99.7%-CL discovery sensitivity in the $^{76}$Ge half-life of $1.3\times10^{28}$ years, corresponding to an effective Majorana mass upper limit in the range of 9-21 meV, to cover the inverted-ordering neutrino mass scale with 10 yr of live time.
We present a method to derive conservative upper limits on the coupling constants of the effective theory of dark matter-nucleon interactions, taking into account the interference among operators. The method can be applied in any basis, and can be easily particularized to any UV complete model. To illustrate our method, we use the IceCube constraints on an exotic neutrino flux from dark matter annihilations in the Sun to derive conservative upper limits on the dark matter-nucleon coupling constants of the effective theory, as well as to derive conservative upper limits on the dark matter-proton and dark matter-neutron scattering cross-sections.
In the last decade, quantum simulators, and in particular cold atoms in optical lattices, have emerged as a valuable tool to study strongly correlated quantum matter. These experiments are now reaching regimes that are numerically difficult or impossible to access. In particular they have started to fulfill a promise which has contributed significantly to defining and shaping the field of cold atom quantum simulations, namely the exploration of doped and frustrated quantum magnets and the search for the origins of high-temperature superconductivity in the fermionic Hubbard model. Despite many future challenges lying ahead, such as the need to further lower the experimentally accessible temperatures, remarkable studies have already been conducted. Among them, spin-charge separation in one-dimensional systems has been demonstrated, extended-range antiferromagnetism in two-dimensional systems has been observed, connections to modern day large-scale numerical simulations were made, and unprecedented comparisons with microscopic trial wavefunctions have been carried out at finite doping. In many regards, the field has acquired new realms, putting old ideas to a new test and producing new insights and inspiration for the next generation of physicists. In the first part of this paper, we review the results achieved in cold atom realizations of the Fermi–Hubbard model in recent years. We put special emphasis on the new probes available in quantum gas microscopes, such as higher-order correlation functions, full counting statistics, the ability to study far-from-equilibrium dynamics, machine learning and pattern recognition of instantaneous snapshots of the many-body wavefunction, and access to non-local correlators. Our review is written from a theoretical perspective, but aims to provide basic understanding of the experimental procedures. We cover one-dimensional systems, where the phenomenon of spin-charge separation is ubiquitous, and two-dimensional systems where we distinguish between situations with and without doping. Throughout, we focus on the strong coupling regime where the Hubbard interactions <math display="inline" id="d1e3228" altimg="si1.svg"><mi>U</mi></math> dominate and connections to <math display="inline" id="d1e3233" altimg="si421.svg"><mrow><mi>t</mi><mo linebreak="goodbreak" linebreakstyle="after">−</mo><mi>J</mi></mrow></math> models can be justified. In the second part of this paper, with the stage set and the current state of the field in mind, we propose a new direction for cold atoms to explore: namely mixed-dimensional bilayer systems, where the charge motion is restricted to individual layers which remain coupled through spin-exchange. These systems can be directly realized experimentally and we argue that they have a rich phase diagram, potentially including a strongly correlated BEC-to-BCS cross-over and regimes with different superconducting order parameters, as well as complex parton phases and possibly even analogs of tetraquark states. In particular, we propose a novel, strong pairing mechanism in these systems, which puts the formation of hole pairs at experimentally accessible, elevated temperatures within reach. Ultimately we propose to explore how the physics of the mixed-dimensional bilayer system can be connected to the rich phenomenology of the single-layer Hubbard model. •Comprehensive review of cold atom experiments on the Fermi–Hubbard model.•Focus on physics highlights, from a theoretical perspective.•New results for hole-pairing in bilayer and ladder systems.•Including a discussion of contemporary analysis tools: from machine-learning to ARPES.
We present a model for the halo-mass correlation function that explicitly incorporates halo exclusion and allows for a redefinition of the halo boundary in a flexible way. We assume that haloes trace mass in a way that can be described using a single scale-independent bias parameter. However, our model exhibits scale-dependent biasing due to the impact of halo-exclusion, the use of a 'soft' (i.e. not infinitely sharp) halo boundary, and differences in the one halo term contributions to ξhm and ξmm. These features naturally lead us to a redefinition of the halo boundary that lies at the 'by eye' transition radius from the one-halo to the two-halo term in the halo-mass correlation function. When adopting our proposed definition, our model succeeds in describing the halo-mass correlation function with $\approx 2{{\ \rm per\ cent}}$ residuals over the radial range 0.1 h-1 Mpc < r < 80 h-1 Mpc, and for halo masses in the range 1013 h-1 M⊙ < M < 1015 h-1 M⊙. Our proposed halo boundary is related to the splashback radius by a roughly constant multiplicative factor. Taking the 87 percentile as reference we find rt/Rsp ≍ 1.3. Surprisingly, our proposed definition results in halo abundances that are well described by the Press-Schechter mass function with δsc = 1.449 ± 0.004. The clustering bias parameter is offset from the standard background-split prediction by $\approx 10{{\ \rm per\ cent}}\!-\!15{{\ \rm per\ cent}}$. This level of agreement is comparable to that achieved with more standard halo definitions.
We study the gravitational radiation emitted during the scattering of two spinless bodies in the post-Minkowskian effective field theory approach. We derive the conserved stress-energy tensor linearly coupled to gravity and the classical probability amplitude of graviton emission at leading and next-to-leading order in the Newton's constant G . The amplitude can be expressed in compact form as one-dimensional integrals over a Feynman parameter involving Bessel functions. We use it to recover the leading-order radiated angular momentum expression. Upon expanding it in the relative velocity between the two bodies v , we compute the total four-momentum radiated into gravitational waves at leading-order in G and up to an order v, 8 finding agreement with what was recently computed using scattering amplitude methods. Our results also allow us to investigate the zero frequency limit of the emitted energy spectrum.
The evolution of young stars and disks is driven by the interplay of several processes, notably the accretion and ejection of material. These processes, critical to correctly describe the conditions of planet formation, are best probed spectroscopically. Between 2020 and 2022, about 500orbits of the Hubble Space Telescope (HST) are being devoted in to the ULLYSES public survey of about 70 low-mass (M⋆ ≤ 2 M⊙) young (age < 10 Myr) stars at UV wavelengths. Here, we present the PENELLOPE Large Program carried out with the ESO Very Large Telescope (VLT) with the aim of acquiring, contemporaneously to the HST, optical ESPRESSO/UVES high-resolution spectra for the purpose of investigating the kinematics of the emitting gas, along with UV-to-NIR X-shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time and the fully reduced spectra are being made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of 13 targets located in the Orion OB1 association and in the σ-Orionis cluster, observed in November-December 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days (≲3). The comparison of the fits to the continuum excess emission obtained with a slab model on the X-shooter spectra and the HST/STIS spectra shows a shortcoming in the X-shooter estimates of ≲10%, which is well within the assumed uncertainty. Its origin can be either due to an erroneous UV extinction curve or to the simplicity of the modeling and, thus, this question will form the basis of the investigation undertaken over the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key in attaining a better understanding of the accretion and ejection mechanisms in young stars.
Based on observations collected at the European Southern Observatory under ESO programme 106.20Z8.
We use the forward modeling approach to galaxy clustering combined with the likelihood from the effective-field theory of large-scale structure to measure assembly bias, i.e. the dependence of halo bias on properties beyond the total mass, in the linear (b1) and second order bias parameters (b2 and bK2) of dark matter halos in N-body simulations. This is the first time that assembly bias in the tidal bias parameter bK2 is measured. We focus on three standard halo properties: the concentration c, spin λ, and sphericity s, for which we find an assembly bias signal in bK2 that is opposite to that in b1. Specifically, at fixed mass, halos that get more (less) positively biased in b1, get less (more) negatively biased in bK2. We also investigate the impact of assembly bias on the b2(b1) and bK2(b1) relations, and find that while the b2(b1) relation stays roughly unchanged, assembly bias strongly impacts the bK2(b1) relation. This impact likely extends also to the corresponding relation for galaxies, which motivates future studies to design better priors on bK2(b1) for use in cosmological constraints from galaxy clustering data.
We study and model the properties of galaxy clusters in the normal-branch Dvali–Gabadadze–Porrati (nDGP) model of gravity, which is representative of a wide class of theories that exhibit the Vainshtein screening mechanism. Using the first cosmological simulations that incorporate both full baryonic physics and nDGP, we find that, despite being efficiently screened within clusters, the fifth force can raise the temperature of the intracluster gas, affecting the scaling relations between the cluster mass and three observable mass proxies: the gas temperature, the Compton Y-parameter of the Sunyaev–Zel’dovich effect, and the X-ray analogue of the Y-parameter. Therefore, unless properly accounted for, this could lead to biased measurements of the cluster mass in tests that use cluster observations, such as cluster number counts, to probe gravity. Using a suite of dark-matter-only simulations, which span a wide range of box sizes and resolutions, and which feature very different strengths of the fifth force, we also calibrate general fitting formulae that can reproduce the nDGP halo concentration at percent accuracy for 0 ≤ z ≤ 1, and halo mass function with |${\lesssim}3{{\ \rm per\ cent}}$| accuracy at 0 ≤ z ≤ 1 (increasing to |${\lesssim}5{{\ \rm per\ cent}}$| for 1 ≤ z ≤ 2), over a halo mass range spanning four orders of magnitude. Our model for the concentration can be used for converting between halo mass overdensities and predicting statistics such as the non-linear matter power spectrum. The results of this work will form part of a framework for unbiased constraints of gravity using the data from ongoing and upcoming cluster surveys.
Six binary-merger progenitors of supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak -69°202 are exploded by neutrino heating and evolved until long after shock breakout in 3D and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize 56Ni masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km s-1 compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive 56Ni mixing depends mainly on high growth factors of Rayleigh-Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak -69°202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first ~7 days. Models that mix about 0.5 M⊙ of hydrogen into the He-shell and exhibit strong outward mixing of 56Ni with maximum velocities exceeding the 3000 km s-1 observed for the bulk of ejected 56Ni have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one.
We present theoretical predictions for the free-free emission at centimeter wavelengths obtained from photoevaporation and magnetohydrodynamic (MHD) wind disk models adjusted to the case of the TW Hydrae young stellar object. For this system, disk photoevaporation with heating due to the high-energy photons from the star has been proposed as a possible mechanism to open the gap observed in the dust emission with the Atacama Large Millimeter/submillimeter Array. We show that the photoevaporation disk model predicts a radial profile for the free-free emission that is made of two main spatial components, one originated from the bound disk atmosphere at 0.5-1 au from the star, and another more extended component from the photoevaporative wind at larger disk radii. We also show that the stellar X-ray luminosity has a significant impact on both these components. The predicted radio emission from the MHD wind model has a smoother radial distribution which extends to closer distances to the star than the photoevaporation case. We also show that a future radio telescope such as the Next Generation Very Large Array would have enough sensitivity and angular resolution to spatially resolve the main structures predicted by these models.
We calculate models of stellar evolution for very massive stars and include the effects of modified gravity to investigate the influence on the physical properties of blue supergiant stars and their use as extragalactic distance indicators. With shielding and fifth force parameters in a similar range as those in previous studies of Cepheid and tip of the red giant branch (TRGB) stars, we find clear effects on stellar luminosity and flux-weighted gravity. The relationship between flux-weighted gravity, gF ≡ g/ ${T}_{\mathrm{eff}}^{4}$ , and bolometric magnitude Mbol, which has been used successfully for accurate distance determinations, is systematically affected. While the stellar evolution of flux-weighted gravity-luminosity relationships (FGLRs) show a systematic offset from the observed relation, we can use the differential shifts between models with Newtonian and modified gravity to estimate the influence on FGLR distance determinations. Modified gravity leads to an increase in distance of 0.05-0.15 magnitudes in distance modulus. These changes are comparable to the ones found for Cepheid stars. We compare observed FGLR and TRGB distances of nine galaxies to constrain the free parameters of modified gravity. Not accounting for systematic differences between TRGB and FGLR distances shielding parameters of 5 × 10-7 and 10-6 and fifth force parameters of 1/3 and 1 can be ruled out with about 90% confidence. Allowing for potential systematic offsets between TRGB and FGLR distances no determination is possible for a shielding parameter of 10-6. For 5 × 10-7 a fifth force parameter of 1 can be ruled out to 92% but 1/3 is unlikely only to 60%.
The impact of new and highly precise neutron \b{eta} decay data is reviewed. We focus on recent results from neutron lifetime, \b{eta} asymmetry, and electron-neutrino correlation experiments. From these results, weak interaction parameters are extracted with unprecedented precision, possible also due to progress in effective field theory and lattice QCD. Limits on New Physics beyond the Standard Model derived from neutron decay data are sharper than those derived from high-energy experiments, except for processes involving right-handed neutrinos.
Recent millimeter and infrared observations have shown that gap- and ring-like structures are common in both dust thermal emission and scattered light of protoplanetary disks. We investigate the impact of the so-called thermal wave instability (TWI) on the millimeter and infrared scattered light images of disks. We perform 11D simulations of the TWI and confirm that the TWI operates when the disk is optically thick enough for stellar light, i.e., small-grain-to-gas mass ratio of 0.0001. The midplane temperature varies as the waves propagate, and hence gap and ring structures can be seen in both millimeter and infrared emission. The millimeter substructures can be observed even if the disk is fully optically thick since it is induced by the temperature variation, while density-induced substructures would disappear in the optically thick regime. The fractional separation between TWI-induced ring and gap is r/r 0.20.4 at 1050 au, which is comparable to those found by the Atacama Large Millimeter/submillimeter Array. Due to the temperature variation, snow lines of volatile species move radially and multiple snow lines are observed even for a single species. The wave propagation velocity is as fast as 0.6 au yr1, which can be potentially detected with a multiepoch observation with a time separation of a few years.
Context. In 2019, the eROSITA telescope on board the Russian-German satellite Spectrum-Roentgen-Gamma (SRG) began to perform a deep all-sky X-ray survey with the aim of identifying ~100 000 clusters and groups over the course of four years. As part of its performance verification phase, a ~140 deg2 survey, called eROSITA Final Equatorial-Depth Survey (eFEDS), was performed. With a depth typical of the all-sky survey after four years, it allows tests of tools and methods as well as improved predictions for the all-sky survey.Aims. As part of this effort, a catalog of 542 X-ray selected galaxy group and cluster candidates was compiled. In this paper we present the optical follow-up, with the aim of providing redshifts and cluster confirmation for the full sample. Furthermore, we aim to provide additional information on the dynamical state, richness, and optical center of the clusters. Finally, we aim to evaluate the impact of optical cluster confirmation on the purity and completeness of the X-ray selected sample.Methods. We used optical imaging data from the Hyper Suprime-Cam Subaru Strategic Program and from the Legacy Survey to identify optical counterparts to the X-ray detected cluster candidates. We make use of the multi-component matched filter cluster confirmation tool (MCMF), as well as of the optical cluster finder CAMIRA to derive cluster redshifts and richnesses. MCMF provided the probabilities with which an optical structure would be a chance superposition with the X-ray candidate. These probabilities were used to identify the best optical counterpart as well as to confirm an X-ray candidate as a cluster. The impact of this confirmation process on catalog purity and completeness was estimated using optical to X-ray scaling relations as well as simulations. The resulting catalog was furthermore matched with public group and cluster catalogs. Optical estimators of the cluster dynamical state were constructed based on density maps of the red-sequence galaxies at the cluster redshift.Results. By providing redshift estimates for all 542 candidates, we construct an optically confirmed sample of 477 clusters and groups with a residual contamination of 6%. Of these, 470 (98.5%) are confirmed using MCMF, and 7 systems are added through cross-matching with spectroscopic group catalogs. Using observable-to-observable scaling and the applied confirmation threshold, we predict that 8 ± 2 real systems have been excluded with the MCMF cut required to build this low-contamination sample. This number agrees well with the 7 systems found through cross-matching that were not confirmed with MCMF. The predicted redshift and mass distribution of this catalog agree well with simulations. Thus, we expect that these 477 systems include >99% of all true clusters in the candidate list. Using an MCMF-independent method, we confirm that the catalog contamination of the confirmed subsample is 6 ± 3%. Application of the same method to the full candidate list yields 17 ± 3%, consistent with estimates coming from the fraction of confirmed systems of ~17% and with expectations from simulations of ~20%. We also present a sample of merging cluster candidates based on the derived estimators of the cluster dynamical state.Key words: catalogs / galaxies: clusters: general / galaxies: distances and redshifts / galaxies: clusters: intracluster medium / X-rays: galaxies: clusters★ The catalog is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/661/A4
We provide the largest and most homogeneous sample of α-element (Mg, Ca, Ti) and iron abundances for field RR Lyrae (RRLs; 162 variables) by using high-resolution spectra. The current measurements were complemented with similar abundances available in the literature for 46 field RRLs brought to our metallicity scale. We ended up with a sample of old (t ≥ 10 Gyr), low-mass stellar tracers (208 RRLs: 169 fundamental, 38 first overtone, and 1 mixed mode) covering 3 dex in iron abundance (-3.00 ≤ [Fe/H] ≤ 0.24). We found that field RRLs are ~0.3 dex more α poor than typical halo tracers in the metal-rich regime ([Fe/H] ≥ -1.2), while in the metal-poor regime ([Fe/H] ≤ -2.2) they seem to be on average ~0.1 dex more α enhanced. This is the first time that the depletion in α elements for solar iron abundances is detected on the basis of a large, homogeneous, and coeval sample of old stellar tracers. Interestingly, we also detected a close similarity in the [α/Fe] trend between α-poor, metal-rich RRLs and red giants (RGs) in the Sagittarius dwarf galaxy as well as between α-enhanced, metal-poor RRLs and RGs in ultrafaint dwarf galaxies. These results are supported by similar elemental abundances for 46 field horizontal branch stars. These stars share with RRLs the same evolutionary phase and the same progenitors. This evidence further supports the key role that old stellar tracers play in constraining the early chemical enrichment of the halo and, in particular, in investigating the impact that dwarf galaxies have had in the mass assembly of the Galaxy. * Based on observations obtained with the du Pont telescope at Las Campanas Observatory, operated by Carnegie Institution for Science. Based in part on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. Based partly on data obtained with the STELLA robotic telescopes in Tenerife, an AIP facility jointly operated by AIP and IAC. Some of the observations reported in this paper were obtained with the Southern African Large Telescope (SALT). Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere.
Context. In November 2019, eROSITA on board of the Spektrum-Roentgen-Gamma (SRG) observatory started to map the entire sky in X-rays. After the four-year survey program, it will reach a flux limit that is about 25 times deeper than ROSAT. During the SRG performance verification phase, eROSITA observed a contiguous 140 deg2 area of the sky down to the final depth of the eROSITA all-sky survey (eROSITA Final Equatorial-Depth Survey; eFEDS), with the goal of obtaining a census of the X-ray emitting populations (stars, compact objects, galaxies, clusters of galaxies, and active galactic nuclei) that will be discovered over the entire sky.Aims. This paper presents the identification of the counterparts to the point sources detected in eFEDS in the main and hard samples and their multi-wavelength properties, including redshift.Methods. To identifyy the counterparts, we combined the results from two independent methods (NWAY and ASTROMATCH), trained on the multi-wavelength properties of a sample of 23k XMM-Newton sources detected in the DESI Legacy Imaging Survey DR8. Then spectroscopic redshifts and photometry from ancillary surveys were collated to compute photometric redshifts.Results. Of the eFEDS sources, 24 774 of 27 369 have reliable counterparts (90.5%) in the main sample and 231 of 246 sourcess (93.9%) have counterparts in the hard sample, including 2514 (3) sources for which a second counterpart is equally likely. By means of reliable spectra, Gaia parallaxes, and/or multi-wavelength properties, we have classified the reliable counterparts in both samples into Galactic (2695) and extragalactic sources (22 079). For about 340 of the extragalactic sources, we cannot rule out the possibility that they are unresolved clusters or belong to clusters. Inspection of the distributions of the X-ray sources in various optical/IR colour-magnitude spaces reveal a rich variety of diverse classes of objects. The photometric redshifts are most reliable within the KiDS/VIKING area, where deep near-infrared data are also available.Conclusions. This paper accompanies the eROSITA early data release of all the observations performed during the performance and verification phase. Together with the catalogues of primary and secondary counterparts to the main and hard samples of the eFEDS survey, this paper releases their multi-wavelength properties and redshifts.Key words: methods: data analysis / X-rays: general / catalogs / surveys / galaxies: active / galaxies: distances and redshifts★ The data are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/661/A3
New neutral heavy gauge bosons (Z') are predicted within many extensions of the Standard Model. While in case they couple to quarks the LHC bounds are very stringent, leptophilic Z' bosons (even with sizable couplings) can be much lighter and therefore lead to interesting quantum effects in precision observables (like (g − 2)μ) and generate flavour violating decays of charged leptons. In particular, ℓ →ℓ 'v v ¯ decays, anomalous magnetic moments of charged leptons, ℓ → ℓ'γ and ℓ → 3ℓ' decays place stringent limits on leptophilic Z' bosons. Furthermore, in case of mixing Z' with the SM Z, Z pole observables are affected. In light of these many observables we perform a global fit to leptophilic Z' models with the main goal of finding the bounds for the Z' couplings to leptons. To this end we consider a number of scenarios for these couplings. While in generic scenarios correlations are weak, this changes once additional constraints on the couplings are imposed. In particular, if one considers an Lμ− Lτ symmetry broken only by left-handed rotations, or considers the case of τ − μ couplings only. In the latter setup, on can explain the (g − 2)μ anomaly and the hint for lepton flavour universality violation in τ →μv v ¯/τ →ev v ¯ without violating bounds from electroweak precision observables.
We study gauge theories of background fields associated to BRST quantized spinning particle models and identify background-independent algebraic structures which allow to systematically reduce the spectrum of fields and subject some of them to dynamical equations of motion. More specifically, we construct a manifestly background-independent extension of the model based on $N=2$ spinning particle. The extended system describes the on-shell spin-1 field coupled to off-shell background fields including metric and dilaton. Tensoring with a given Lie algebra results in a non-abelian extension of the model.
We present a joint likelihood analysis of the real-space power spectrum and bispectrum measured from a variety of halo and galaxy mock catalogs. A novel aspect of this work is the inclusion of nonlinear triangle configurations for the bispectrum, made possible by a complete next-to-leading order ("one-loop") description of galaxy bias, as is already common practice for the power spectrum. Based on the goodness of fit and the unbiasedness of the parameter posteriors, we accomplish a stringent validation of this model compared to the leading order ("tree-level") bispectrum. Using measurement uncertainties that correspond to an effective survey volume of 6 (Gpc /h )3 , we determine that the one-loop corrections roughly double the applicable range of scales, from ∼0.17 h /Mpc (tree level) to ∼0.3 h /Mpc . This converts into a 1.5-2x improvement on constraints of the linear bias parameter at fixed cosmology, and a 1.5-2.4x shrinkage of uncertainties on the amplitude of fluctuations As, which clearly demonstrates the benefit of extracting information from nonlinear scales despite having to marginalize over a larger number of bias parameters. Besides, our precise measurements of galaxy bias parameters up to fourth order allow for thorough comparisons to coevolution relations, showing excellent agreement for all contributions generated by the nonlocal action of gravity. Using these relations in the likelihood analysis does not compromise the model validity and is crucial for obtaining the quoted improvements on As. We also analyzed the impact of higher-derivative and scale-dependent stochastic terms, finding that for a subset of our tracers the former can boost the performance of the tree-level model with constraints on As that are only slightly degraded compared to the one-loop model.
The Method of Moments is a powerful framework to disentangle the relative contributions of amplitudes of a specific process to its various phase space regions. We apply this method to carry out a fully differential analysis of the Higgs decay channel h → 4ℓ and constrain gauge-Higgs coupling modifications parametrised by dimension-six effective operators. We find that this analysis approach provides very good constraints and minimises degeneracies in the parameter space of the effective theory. By combining the decay h → 4ℓ with Higgs-associated production processes, Wh and Zh, we obtain the strongest reported bounds on anomalous gauge-Higgs couplings.
Finite density effects can destabilize the metastable vacua in relaxion models. Focusing on stars as nucleation seeds, we derive the conditions that lead to the formation and runaway of a relaxion bubble of a lower energy minimum than in vacuum. The resulting late-time phase transition in the universe allows us to set new constraints on the parameter space of relaxion models. We also find that similar instabilities can be triggered by the large electromagnetic fields around rotating neutron stars.
We present simulations of the multiphase interstellar medium (ISM) at solar neighbourhood conditions including thermal and non-thermal ISM processes, star cluster formation, and feedback from massive stars: stellar winds, hydrogen ionizing radiation computed with the novel TREERAY radiative transfer method, supernovae (SN), and the injection of cosmic rays (CR). N-body dynamics is computed with a 4th-order Hermite integrator. We systematically investigate the impact of stellar feedback on the self-gravitating ISM with magnetic fields, CR advection and diffusion, and non-equilibrium chemical evolution. SN-only feedback results in strongly clustered star formation with very high star cluster masses, a bi-modal distribution of the ambient SN densities, and low volume-filling factors (VFF) of warm gas, typically inconsistent with local conditions. Early radiative feedback prevents an initial starburst, reduces star cluster masses and outflow rates. Furthermore, star formation rate surface densities of $\Sigma _{\dot{M}_\star } = 1.4-5.9 \times 10^{-3}$$\mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$, VFFwarm = 60-80 per cent as well as thermal, kinetic, magnetic, and cosmic ray energy densities of the model including all feedback mechanisms agree well with observational constraints. On the short, 100 Myr, time-scales investigated here, CRs only have a moderate impact on star formation and the multiphase gas structure and result in cooler outflows, if present. Our models indicate that at low gas surface densities SN-only feedback only captures some characteristics of the star-forming ISM and outflows/inflows relevant for regulating star formation. Instead, star formation is regulated on star cluster scales by radiation and winds from massive stars in clusters, whose peak masses agree with solar neighbourhood estimates.
We develop a formalism for computing inclusive production cross sections of heavy quarkonia based on the nonrelativistic QCD and the potential nonrelativistic QCD effective field theories. Our formalism applies to strongly coupled quarkonia, which include excited charmonium and bottomonium states. Analogously to heavy quarkonium decay processes, we express nonrelativistic QCD long-distance matrix elements in terms of quarkonium wavefunctions at the origin and universal gluonic correlators. Our expressions for the long-distance matrix elements are valid up to corrections of order $ 1/{N}_c^2 $. These expressions enhance the predictive power of the nonrelativistic effective field theory approach to inclusive production processes by reducing the number of nonperturbative unknowns, and make possible first-principle determinations of long-distance matrix elements once the gluonic correlators are known. Based on this formalism, we compute the production cross sections of P-wave charmonia and bottomonia at the LHC, and find good agreement with measurements.
Cosmogenic radio-nuclei are an important source of background for low-energy neutrino experiments. In Borexino, cosmogenic $^{11}$C decays outnumber solar pep and CNO neutrino events by about ten to one. In order to extract the flux of these two neutrino species, a highly efficient identification of this background is mandatory. We present here the details of the most consolidated strategy, used throughout Borexino solar neutrino measurements. It hinges upon finding the space-time correlations between $^{11}$C decays, the preceding parent muons and the accompanying neutrons. This article describes the working principles and evaluates the performance of this Three-Fold Coincidence (TFC) technique in its two current implementations: a hard-cut and a likelihood-based approach. Both show stable performances throughout Borexino Phases II (2012–2016) and III (2016–2020) data sets, with a $^{11}$C tagging efficiency of $\sim 90$ % and $\sim $ 63–66 % of the exposure surviving the tagging. We present also a novel technique that targets specifically $^{11}$C produced in high-multiplicity during major spallation events. Such $^{11}$C appear as a burst of events, whose space-time correlation can be exploited. Burst identification can be combined with the TFC to obtain about the same tagging efficiency of $\sim 90\%$ but with a higher fraction of the exposure surviving, in the range of $\sim $ 66–68 %.
On the way of a microscopic derivation of covariant density functionals, the first complete solution of the relativistic Brueckner-Hartree-Fock (RBHF) equations is presented for symmetric nuclear matter. In most of the earlier investigations, the G-matrix is calculated only in the space of positive energy solutions. On the other side, for the solution of the relativistic Hartree-Fock (RHF) equations, also the elements of this matrix connecting positive and negative energy solutions are required. So far, in the literature, these matrix elements are derived in various approximations. We discuss solutions of the Thompson equation for the full Dirac space and compare the resulting equation of state with those of earlier attempts in this direction.
We compare gravitational-wave (GW) signals from eight 3D simulations of core-collapse supernovae, using two different progenitors with zero-age main-sequence masses of 9 and 20 solar masses (M⊙). The collapse of each progenitor was simulated four times, at two different grid resolutions and with two different neutrino transport methods, using the AENUS-ALCAR code. The main goal of this study is to assess the validity of recent concerns that the so-called 'Ray-by-Ray+' (RbR+) approximation is problematic in core-collapse simulations and can adversely affect theoretical GW predictions. Therefore, signals from simulations using RbR+ are compared to signals from corresponding simulations using a fully multidimensional (FMD) transport scheme. The 9 M⊙ progenitor successfully explodes, whereas the 20 M⊙ model does not. Both the standing accretion shock instability and hot-bubble convection develop in the post-shock layer of the non-exploding models. In the exploding models, neutrino-driven convection in the post-shock flow is established around 100 ms after core bounce and lasts until the onset of shock revival. We can, therefore, judge the impact of the numerical resolution and neutrino transport under all conditions typically seen in non-rotating core-collapse simulations. We find excellent qualitative agreement in all GW features. We find minor quantitative differences between simulations, but find no systematic differences between simulations using different transport schemes. Resolution-dependent differences in the hydrodynamic behaviour of low-resolution and high-resolution models have a greater impact on the GW signals than consequences of the different transport methods. Furthermore, increasing the resolution decreases the discrepancies between models with different neutrino transport.
We consider matter density effects in theories with a false ground state. Large and dense systems, such as stars, can destabilize a metastable minimum and allow for the formation of bubbles of the true minimum. We derive the conditions under which these bubbles form, as well as the conditions under which they either remain confined to the dense region or escape to infinity. The latter case leads to a phase transition in the universe at star formation. We explore the phenomenological consequences of such seeded phase transitions.
Our work presents an independent calibration of the J-region Asymptotic Giant Branch (JAGB) method using Infrared Survey Facility (IRSF) photometric data and a custom luminosity function profile to determine JAGB mean magnitudes for nine galaxies. We determine a mean absolute magnitude of carbon stars of $M_{LMC}=-6.212 \pm 0.010$ (stat.) $\pm 0.030$ (syst.) mag. We then use near-infrared photometry of a number of nearby galaxies, originally obtained by our group to determine their distances from Cepheids using the Leavitt law, in order to independently determine their distances with the JAGB method. We compare the JAGB distances obtained in this work with the Cepheid distances resulting from the same photometry and find very good agreement between the results from the two methods. The mean difference is 0.01 mag with an rms scatter of 0.06 mag after taking into account seven out of the eight analyzed galaxies that had their distances determined using Cepheids. The very accurate distance to the Small Magellanic Cloud (SMC) based on detached eclipsing binaries (Graczyk et al. 2020) is also in very good agreement with the distance obtained from carbon stars.
Upcoming surveys will use a variety of galaxy selections to map the large-scale structure of the Universe. It is important to make accurate predictions for the properties and clustering of such galaxies, including the errors on these statistics. Here, we describe a novel technique which uses the semi-analytical model of galaxy formation GALFORM, embedded in the high-resolution N-body Planck-Millennium simulation, to populate a thousand halo catalogues generated using the Parallel-PM N-body GLAM code. Our hybrid scheme allows us to make clustering predictions on scales that cannot be modelled in the original N-body simulation. We focus on luminous red galaxies (LRGs) selected in the redshift range z = 0.6 - 1 from the GALFORM output using similar colour-magnitude cuts in the r, z, and W1 bands to those that will be applied in the Dark Energy Spectroscopic Instrument (DESI) survey, and call this illustrative sample 'DESI-like' LRGs. We find that the LRG-halo connection is non-trivial, leading to the prediction of a non-standard halo occupation distribution; in particular, the occupation of central galaxies does not reach unity for the most massive haloes, and drops with increasing mass. The GLAM catalogues reproduce the abundance and clustering of the LRGs predicted by GALFORM. We use the GLAM mocks to compute the covariance matrices for the two-point correlation function and power spectrum of the LRGs and their background dark matter density field, revealing important differences. We also make predictions for the linear-growth rate and the baryon acoustic oscillations distances at z = 0.6, 0.74, and 0.93. All 'DESI-like' LRG catalogues are made publicly available.
We investigate the algebra of vector fields on the sphere. First, we find that linear deformations of this algebra are obstructed under reasonable conditions. In particular, we show that $hs[\lambda]$, the one-parameter deformation of the algebra of area-preserving vector fields, does not extend to the entire algebra. Next, we study some non-central extensions through the embedding of $\mathfrak{v}\mathfrak{e}\mathfrak{c}\mathfrak{t}(S^2)$ into $\mathfrak{v}\mathfrak{e}\mathfrak{c}\mathfrak{t}(\mathbb{C}^*)$. For the latter, we discuss a three parameter family of non-central extensions which contains the symmetry algebra of asymptotically flat and asymptotically Friedmann spacetimes at future null infinity, admitting a simple free field realization.
Stars in the mass range from 8 M⊙ to 10 M⊙ are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes and the observables of these SNe are, therefore, expected to be similar. In this study, we explore the differences in these two types of SNe. Specifically, we investigate three different progenitor models: a solar-metallicity ECSN progenitor with an initial mass of 8.8 M⊙, a zero-metallicity progenitor with 9.6 M⊙, and a solar-metallicity progenitor with 9 M⊙, carrying out radiative transfer simulations for these progenitors. We present the resulting light curves for these models. The models exhibit very low photospheric velocity variations of about 2000 km s-1; therefore, this may serve as a convenient indicator of low-mass SNe. The ECSN has very unique light curves in broad-bands, especially the U band, and does not resemble any currently observed SN. This ECSN progenitor being part of a binary will lose its envelope for which reason the light curve becomes short and undetectable. The SN from the 9.6 M⊙ progenitor exhibits also quite an unusual light curve, explained by the absence of metals in the initial composition. The artificially iron-polluted 9.6 M⊙ model demonstrates light curves closer to normal SNe IIP. The SN from the 9 M⊙ progenitor remains the best candidate for so-called low-luminosity SNe IIP like SN 1999br and SN 2005cs.
We present the first application of the extended Fast Action Minimization method (eFAM) to a real data set, the SDSS-DR12 Combined Sample, to reconstruct galaxies orbits back-in-time, their two-point correlation function (2PCF) in real-space, and enhance the baryon acoustic oscillation (BAO) peak. For this purpose, we introduce a new implementation of eFAM that accounts for selection effects, survey footprint, and galaxy bias. We use the reconstructed BAO peak to measure the angular diameter distance, $D_\mathrm{A}(z)r^\mathrm{fid}_\mathrm{s}/r_\mathrm{s}$ , and the Hubble parameter, $H(z)r_\mathrm{s}/r^\mathrm{fid}_\mathrm{s}$ , normalized to the sound horizon scale for a fiducial cosmology $r^\mathrm{fid}_\mathrm{s}$ , at the mean redshift of the sample z = 0.38, obtaining $D_\mathrm{A}(z=0.38)r^\mathrm{fid}_\mathrm{s}/r_\mathrm{s}=1090\pm 29$ (Mpc)-1, and $H(z=0.38)r_\mathrm{s}/r^\mathrm{fid}_\mathrm{s}=83\pm 3$ (km s-1 Mpc-1), in agreement with previous measurements on the same data set. The validation tests, performed using 400 publicly available SDSS-DR12 mock catalogues, reveal that eFAM performs well in reconstructing the 2PCF down to separations of ∼25h-1Mpc, i.e. well into the non-linear regime. Besides, eFAM successfully removes the anisotropies due to redshift-space distortion (RSD) at all redshifts including that of the survey, allowing us to decrease the number of free parameters in the model and fit the full-shape of the back-in-time reconstructed 2PCF well beyond the BAO peak. Recovering the real-space 2PCF, eFAM improves the precision on the estimates of the fitting parameters. When compared with the no-reconstruction case, eFAM reduces the uncertainty of the Alcock-Paczynski distortion parameters α⊥ and α∥ of about 40 per cent and that on the non-linear damping scale Σ∥ of about 70 per cent. These results show that eFAM can be successfully applied to existing redshift galaxy catalogues and should be considered as a reconstruction tool for next-generation surveys alternative to popular methods based on the Zel'dovich approximation.
We present cosmological zoom-in hydrodynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain-size distribution, with the MUPPI star formation and feedback subresolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas, and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio, and small-to-large grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose.
We showed how to use trained neural networks to perform Bayesian reasoning in order to solve tasks outside their initial scope. Deep generative models provide prior knowledge, and classification/regression networks impose constraints. The tasks at hand were formulated as Bayesian inference problems, which we approximately solved through variational or sampling techniques. The approach built on top of already trained networks, and the addressable questions grew super-exponentially with the number of available networks. In its simplest form, the approach yielded conditional generative models. However, multiple simultaneous constraints constitute elaborate questions. We compared the approach to specifically trained generators, showed how to solve riddles, and demonstrated its compatibility with state-of-the-art architectures.
We carry out 'full-physics' hydrodynamical simulations of galaxy formation in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld model using a new modified version of the AREPO code and the IllustrisTNG galaxy formation model. We simulate two nDGP models (N5 and N1) that represent, respectively, weak and moderate departures from general relativity (GR), in boxes of sizes $62$ and $25\, h^{-1}{\rm Mpc}$ using 2 × 5123 dark matter particles and initial gas cells. This allows us to explore, for the first time, the impact of baryonic physics on galactic scales in braneworld models of modified gravity and to make predictions on the stellar content of dark matter haloes and galaxy evolution through cosmic time in these models. We find significant differences between the GR and nDGP models in the power spectra and correlation functions of gas, stars and dark matter of up to ~25 per cent on large scales. Similar to their impact in the standard cosmological model (Λ cold dark matter), baryonic effects can have a significant influence over the clustering of the overall matter distribution, with a sign that depends on scale. Studying the degeneracy between modified gravity and galactic feedback in these models, we find that these two physical effects on matter clustering can be cleanly disentangled, allowing for a method to accurately predict the matter power spectrum with baryonic effects included, without having to run hydrodynamical simulations. Depending on the braneworld model, we find differences compared with GR of up to ~15 per cent in galaxy properties such as the stellar-to-halo-mass ratio, galaxy stellar mass function, gas fraction, and star formation rate density. The amplitude of the fifth force is reduced by the presence of baryons in the very inner part of haloes, but this reduction quickly becomes negligible above ~0.1 times the halo radius.
We present here a self-consistent cosmological zoom-in simulation of a triple supermassive black hole (SMBH) system forming in a complex multiple galaxy merger. The simulation is run with an updated version of our code KETJU, which is able to follow the motion of SMBHs down to separations of tens of Schwarzschild radii while simultaneously modeling the large-scale astrophysical processes in the surrounding galaxies, such as gas cooling, star formation, and stellar and AGN feedback. Our simulation produces initially an SMBH binary system for which the hardening process is interrupted by the late arrival of a third SMBH. The KETJU code is able to accurately model the complex behavior occurring in such a triple SMBH system, including the ejection of one SMBH to a kiloparsec-scale orbit in the galaxy due to strong three-body interactions as well as Lidov-Kozai oscillations suppressed by relativistic precession when the SMBHs are in a hierarchical configuration. One pair of SMBHs merges ∼3 Gyr after the initial galaxy merger, while the remaining binary is at a parsec-scale separation when the simulation ends at redshift z = 0. We also show that KETJU can capture the effects of the SMBH binaries and triplets on the surrounding stellar population, which can affect the binary merger timescales as the stellar density in the system evolves. Our results demonstrate the importance of dynamically resolving the complex behavior of multiple SMBHs in galactic mergers, as such systems cannot be readily modeled using simple orbit-averaged semianalytic models.
Parametric and non-parametric classifiers often have to deal with real-world data, where corruptions like noise, occlusions, and blur are unavoidable - posing significant challenges. We present a probabilistic approach to classify strongly corrupted data and quantify uncertainty, despite the model only having been trained with uncorrupted data. A semi-supervised autoencoder trained on uncorrupted data is the underlying architecture. We use the decoding part as a generative model for realistic data and extend it by convolutions, masking, and additive Gaussian noise to describe imperfections. This constitutes a statistical inference task in terms of the optimal latent space activations of the underlying uncorrupted datum. We solve this problem approximately with Metric Gaussian Variational Inference (MGVI). The supervision of the autoencoder's latent space allows us to classify corrupted data directly under uncertainty with the statistically inferred latent space activations. Furthermore, we demonstrate that the model uncertainty strongly depends on whether the classification is correct or wrong, setting a basis for a statistical "lie detector" of the classification. Independent of that, we show that the generative model can optimally restore the uncorrupted datum by decoding the inferred latent space activations.
The frequency of Earth-sized planets in habitable zones appears to be higher around M-dwarfs, making these systems exciting laboratories to investigate planet formation. Observations of protoplanetary disks around very low-mass stars and brown dwarfs remain challenging and little is known about their properties. The disk around CIDA 1 (~0.1-0.2 M⊙) is one of the very few known disks that host a large cavity (20 au radius in size) around a very low-mass star. We present new ALMA observations at Band 7 (0.9 mm) and Band 4 (2.1 mm) of CIDA 1 with a resolution of ~0.05″ × 0.034″. These new ALMA observations reveal a very bright and unresolved inner disk, a shallow spectral index of the dust emission (~2), and a complex morphology of a ring located at 20 au. We also present X-shooter (VLT) observations that confirm the high accretion rate of CIDA 1 of Ṁacc = 1.4 × 10−8 M⊙ yr−1. This high value of Ṁacc, the observed inner disk, and the large cavity of 20 au exclude models of photo-evaporation to explain the observed cavity. When comparing these observations with models that combine planet-disk interaction, dust evolution, and radiative transfer, we exclude planets more massive than 0.5 MJup as the potential origin of the large cavity because with these it is difficult to maintain a long-lived and bright inner disk. Even in this planet mass regime, an additional physical process may be needed to stop the particles from migrating inwards and to maintain a bright inner disk on timescales of millions of years. Such mechanisms include a trap formed by a very close-in extra planet or the inner edge of a dead zone. The low spectral index of the disk around CIDA 1 is difficult to explain and challenges our current dust evolution models, in particular processes like fragmentation, growth, and diffusion of particles inside pressure bumps.
We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.
We report the serendipitous detection of an H2-bearing damped Lyα absorber at z = 0.576 in the spectrum of the QSO J0111-0316 in the Cosmic Ultraviolet Baryon Survey. Spectroscopic observations from Hubble Space Telescope-COS in the far-ultraviolet reveal a damped absorber with log[N(H I)/cm-2] = 20.1 ± 0.2 and log[N(H2)/cm-2] $={18.97}_{-0.06}^{+0.05}$ . The diffuse molecular gas is found in two velocity components separated by Δ ν ≍ 60 km s-1, with >99.9% of the total H2 column density concentrated in one component. At a metallicity of ≍50% of solar, there is evidence for Fe enhancement and dust depletion, with a dust-to-gas ratio κO ≍ 0.4. A galaxy redshift survey conducted with IMACS and LDSS-3C on Magellan reveals an overdensity of nine galaxies at projected distance d ≤ 600 proper kpc (pkpc) and line-of-sight velocity offset Δ νg ≤ 300 km s-1 from the absorber. The closest is a massive, early-type galaxy at d = 41 pkpc that contains ≍70% of the total stellar mass identified at d ≤ 310 pkpc of the H2 absorber. The close proximity of the H2-bearing gas to the quiescent galaxy and the Fe-enhanced chemical abundance pattern of the absorber suggest a physical connection, in contrast to a picture in which DLAs are primarily associated with gas-rich dwarfs. This case study illustrates that deep galaxy redshift surveys are needed to gain insight into the diverse environments that host dense and potentially star-forming gas. * Based on data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory and the NASA/ESA Hubble Space Telescope operated by the Space Telescope Science Institute and the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.
We use field-level forward models of galaxy clustering and the EFT likelihood formalism to study, for the first time for self-consistently simulated galaxies, the relations between the linear b_1 and second-order bias parameters b 2 and b K 2 . The forward models utilize all of the information available in the galaxy distribution up to a given order in perturbation theory, which allows us to infer these bias parameters with high signal-to-noise, even from relatively small volumes (L box = 205 Mpc/h). We consider galaxies from the simulations, and our main result is that the b 2(b 1) and b K 2 (b 1) relations obtained from gravity-only simulations for total mass selected objects are broadly preserved for simulated galaxies selected by stellar mass, star formation rate, color and black hole accretion rate. We also find good agreement between the bias relations of the simulated galaxies and a number of recent estimates for observed galaxy samples. The consistency under different galaxy selection criteria suggests that theoretical priors on these bias relations may be used to improve cosmological constraints based on observed galaxy samples. We do identify some small differences between the bias relations in the hydrodynamical and gravity-only simulations, which we show can be linked to the environmental dependence of the relation between galaxy properties and mass. We also show that the EFT likelihood recovers the value of σ8 to percent-level from various galaxy samples (including splits by color and star formation rate) and after marginalizing over 8 bias parameters. This demonstration using simulated galaxies adds to previous works based on halos as tracers, and strengthens further the potential of forward models to infer cosmology from galaxy data.
In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). We examine magneto-hydrodynamic galaxy simulations with cosmic rays (CRs), accounting for their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those inject turbulence and drive multiple small-scale dynamos, which exponentially amplify magnetic fields. They grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc, bulge, and bubble-shaped magnetically-loaded outflows of our simulated Milky Way-mass galaxy. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. (abridged)
We study the approach to scaling in axion string networks in the radiation era, through measuring the root-mean-square velocity v as well as the scaled mean string separation x . We find good evidence for a fixed point in the phase-space analysis in the variables (x ,v ), providing a strong indication that standard scaling is taking place. We show that the approach to scaling can be well described by a two parameter velocity-one-scale (VOS) model, and show that the values of the parameters are insensitive to the initial state of the network. The string length has also been commonly expressed in terms of a dimensionless string length density ζ , proportional to the number of Hubble lengths of string per Hubble volume. In simulations with initial conditions far from the fixed point ζ is still evolving after half a light-crossing time, which has been interpreted in the literature as a long-term logarithmic growth. We show that all our simulations, even those starting far from the fixed point, are accounted for by a VOS model with an asymptote of ζ*=1.20 ±0.09 (calculated from the string length in the cosmic rest frame) and v*=0.609 ±0.014 .
High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to 108-9 M⊙ dubbed 'giant clumps'. However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below 107-8 M⊙ that are themselves located inside giant complexes with masses up to 108 and sometimes 109 M⊙. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.
We use a modified version of the peak patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N -body simulations at significantly lower computational cost. We employ a spherical collapse model, and provide fitting functions for the modified barrier in the radiation era. The halo mass function at z =629 has a power-law distribution M-0.6 for masses within the range 10-15≲M ≲10-10 M⊙ , with all masses scaling as (ma/50 μ eV )-0.5 . We construct merger trees to estimate the collapse redshift and concentration mass relation, C (M ), which is well described using analytical results from the initial power spectrum and linear growth. Using the calibrated analytic results to extrapolate to z =0 , our method predicts a mean concentration C ∼O (few )×104. The low computational cost of our method makes future investigation of the statistics of rare, dense miniclusters easy to achieve.
We present a novel hierarchical formulation of the fourth-order forward symplectic integrator and its numerical implementation in the GPU-accelerated direct-summation N-body code frost. The new integrator is especially suitable for simulations with a large dynamical range due to its hierarchical nature. The strictly positive integrator sub-steps in a fourth-order symplectic integrator are made possible by computing an additional gradient term in addition to the Newtonian accelerations. All force calculations and kick operations are synchronous so the integration algorithm is manifestly momentum-conserving. We also employ a time-step symmetrization procedure to approximately restore the time-reversibility with adaptive individual time-steps. We demonstrate in a series of binary, few-body and million-body simulations that frost conserves energy to a level of |ΔE/E| ∼ 10-10 while errors in linear and angular momentum are practically negligible. For typical star cluster simulations, we find that frost scales well up to $N_\mathrm{GPU}^\mathrm{max}\sim 4\times N/10^5$ GPUs, making direct-summation N-body simulations beyond N = 106 particles possible on systems with several hundred and more GPUs. Due to the nature of hierarchical integration, the inclusion of a Kepler solver or a regularized integrator with post-Newtonian corrections for close encounters and binaries in the code is straightforward.