Consistent and simultaneous modelling of galaxy clustering and galaxy-galaxy lensing with Subhalo Abundance Matching by Sergio Contreras et al. on Tuesday 22 November

The spatial distribution of galaxies and their gravitational lensing signal

offer complementary tests of galaxy formation physics and cosmology. However,

their synergy can only be fully exploited if both probes are modelled

accurately and consistently. In this paper, we demonstrate that this can be

achieved using an extension of Subhalo Abundance Matching, dubbed SHAMe.

Specifically, we use mock catalogues built from the TNG300 hydrodynamical

simulation to show that SHAMe can simultaneously model the multipoles of the

redshift-space galaxy correlation function and galaxy-galaxy lensing, without

noticeable bias within the statistical sampling uncertainties of a SDSS volume

and on scales r = [0.6-30] Mpc/h. Modelling the baryonic processes in

galaxy-galaxy lensing with a baryonification scheme allows SHAMe's range of

validity to be extended to r = [0.1-30] Mpc/h. Remarkably, our model achieves

this level of precision with just five free parameters beyond those describing

the baryonification model. At fixed cosmology, we find that galaxy-galaxy

lensing provides a general consistency test but little additional information

on galaxy modelling parameters beyond that encoded in the redshift-space

multipoles. It does, however, improve constraints if only the projected

correlation function is available, as in surveys with only photometric

redshifts. We expect SHAMe to have a higher fidelity across a wider range of

scales than more traditional methods such as Halo Occupation Distribution

modelling. Thus it should provide a significantly more powerful and more robust

tool for analysing next-generation large-scale surveys.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.11745v1

Low-Energy Electron-Track Imaging for a Liquid Argon Time-Projection-Chamber Telescope Concept using Probabilistic Deep Learning by M. Buuck et al. on Monday 21 November

The GammaTPC is an MeV-scale single-phase liquid argon

time-projection-chamber gamma-ray telescope concept with a novel dual-scale

pixel-based charge-readout system. It promises to enable a significant

improvement in sensitivity to MeV-scale gamma-rays over previous telescopes.

The novel pixel-based charge readout allows for imaging of the tracks of

electrons scattered by Compton interactions of incident gamma-rays. The two

primary contributors to the accuracy of a Compton telescope in reconstructing

an incident gamma-ray's original direction are its energy and position

resolution. In this work, we focus on using deep learning to optimize the

reconstruction of the initial position and direction of electrons scattered in

Compton interactions, including using probabilistic models to estimate

predictive uncertainty. We show that the deep learning models are able to

predict locations of Compton scatters of MeV-scale gamma-rays from simulated

pixel-based data to better than 0.6 mm RMS error, and are sensitive to the

initial direction of the scattered electron. We compare and contrast different

deep learning uncertainty estimation algorithms for reconstruction

applications. Additionally, we show that event-by-event estimates of the

uncertainty of the locations of the Compton scatters can be used to select

those events that were reconstructed most accurately, leading to improvement in

locating the origin of gamma-ray sources on the sky.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.07805v2

Large Field Polynomial Inflation: Parameter Space, Predictions and Double Eternal Nature by Manuel Drees et al. on Monday 21 November

Simple monomial inflationary scenarios have been ruled out by recent

observations. In this work we revisit the next simplest scenario, a

single--field model where the scalar potential is a polynomial of degree four

which features a concave ``almost'' saddle point. We focus on trans--Planckian

field values. We reparametrize the potential, which greatly simplifies the

procedure for finding acceptbale model parameters. This allows for the first

comprehensive scan of parameter space consistent with recent Planck and

BICEP/Keck 2018 measurements. Even for trans--Planckian field values the

tensor--to--scalar ratio $r$ can be as small as $\mathcal{O}(10^{-8})$, but the

model can also saturate the current upper bound. In contrast to the

small--field version of this model, radiative stability does not lead to strong

constraints on the parameters of the inflaton potential. For very large field

values the potential can be approximated by the quartic term; as well known,

this allows eternal inflation even for field energy well below the reduced

Planck mass $M_{\rm Pl}$, with Hubble parameter $H \sim 10^{-2} M_{\rm Pl}$.

More interestingly, we find a region of parameter space that even supports {\em

two phases of eternal inflation}. The second epoch only occurs if the slope at

the would--be saddle point is very small, and has $H \sim 10^{-5} M_{\rm Pl}$;

it can only be realized if $r \sim 10^{-2}$, within the sensitivity range of

next--generation CMB observations.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2209.07545v2

Three-dimensional numerical simulations of structured GRB jets by Gerardo Urrutia et al. on Monday 21 November

After the detection of GRB 170817A, the first unambiguous off-axis gamma-ray

burst (GRB), several studies tried to understand the structure of GRB jets. The

initial jet structure (directly produced by the central engine) can be

partially preserved, or can be completely modified by the interaction with the

environment. In this study, we perform three-dimensional, special relativistic

hydrodynamics simulations of long GRB jets evolving through a massive

progenitor star. Different jet scenarios were considered: Top-hat, Gaussian

jets dominated by pressure or by kinetic energy, as well as a model of a

supernova (SN) plus a jet both propagating through the progenitor. We found

that, while propagating inside the progenitor star, jets with different initial

structures are nearly indistinguishable. Kinetic dominated jets are faster and

more collimated than pressure dominated jets. The dynamics of jets inside the

progenitor star strongly depends on the presence of an associated SN, which can

substantially decelerate the jet propagation. We show that the initial

structure of GRB jets is preserved, or not, mainly depending on the jet

collimation. The initial structure is preserved in uncollimated jets, i.e. jets

which move through low density environments. Meanwhile, jets which move through

dense environments are shaped by the interaction with the medium and remain

collimated.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.07925v2

Impact of Rastall gravity on mass, radius and sound speed of the pulsar PSR J0740+6620 by Waleed El Hanafy. on Monday 21 November

Millisecond pulsars are perfect laboratories to test possible matter-geometry

coupling and its physical implications in light of recent Neutron Star Interior

Composition Explorer (NICER) observations. We apply Rastall field equations of

gravity, where matter and geometry are nonminimally coupled, to Krori-Barua

interior spacetime whereas the matter source is assumed to be anisotropic

fluid. We show that all physical quantities inside the star can be expressed in

terms of Rastall, $\epsilon$, and compactness, $C=2GM/Rc^2$, parameters. Using

NICER and X-ray Multi-Mirror X-ray observational constraints on the mass and

radius of the pulsar PSR J0740+6620 we determine Rastall parameter to be at

most $\epsilon=0.041$ in the positive range. The obtained solution provides a

stable compact object; in addition the squared sound speed does not violate the

conjectured sound speed $c_s^2\leq c^2/3$ unlike the general relativistic

treatment. We note that no equations of state are assumed; the model however

fits well with linear patterns with bag constants. In general, for

$\epsilon>0$, the theory predicts a slightly larger size star in comparison to

general relativity for the same mass. This has been explained as an additional

force, due to matter-geometry coupling, in the hydrodynamic equilibrium

equation, which contributes to partially diminish the gravitational force

effect. Consequently, we calculate the maximal compactness as allowed by the

strong energy condition to be $C = 0.735$ which is $\sim 2\%$ higher than

general relativity prediction. Moreover, for the surface density at saturation

nuclear density $\rho_{\text{nuc}} = 2.7\times 10^{14}$ g/cm$^3$ we estimate

the maximum mass $M=4 M_\odot$ at radius $R=16$ km.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2209.10287v2

Accuracy and precision of triaxial orbit models II: Viewing angles, shape and orbital structure by Stefano de Nicola et al. on Monday 21 November

We explore the potential of our novel triaxial modeling machinery in

recovering the viewing angles, the shape and the orbit distribution of galaxies

by using a high-resolution $N$-body merger simulation. Our modelling technique

includes several recent advancements. (i) Our new triaxial deprojection

algorithm SHAPE3D is able to significantly shrink the range of possible

orientations of a triaxial galaxy and therefore to constrain its shape relying

only on photometric information. It also allows to probe degeneracies, i.e. to

recover different deprojections at the same assumed orientation. With this

method we can constrain the intrinsic shape of the $N$-body simulation, i.e.

the axis ratios $p=b/a$ and $q=c/a$, with $\Delta p$ and $\Delta q$ $\lesssim$

0.1 using only photometric information. The typical accuracy of the viewing

angles reconstruction is 15-20$^\circ$. (ii) Our new triaxial Schwarzschild

code SMART exploits the full kinematic information contained in the entire

non-parametric line-of-sight velocity distributions (LOSVDs) along with a 5D

orbital sampling in phase space. (iii) We use a new generalised information

criterion AIC$_p$ to optimise the smoothing and to select the best-fit model,

avoiding potential biases in purely $\chi^2$-based approaches. With our

deprojected densities, we recover the correct orbital structure and anisotropy

parameter $\beta$ with $\Delta \beta$ $\lesssim$ 0.1. These results are valid

regardless of the tested orientation of the simulation and suggest that even

despite the known intrinsic photometric and kinematic degeneracies the above

described advanced methods make it possible to recover the shape and the

orbital structure of triaxial bodies with unprecedented accuracy.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2210.01893v2

A measure of cosmological distance using the civ Baldwin effect in quasars by L. Huang et al. on Monday 21 November

We use the anticorrelation between the equivalent width (EW) of the

C\,\textsc{iv} 1549 {\AA} emission line and the continuum luminosity in the

quasars rest frame (Baldwin effect) to measure their luminosity distance as

well as estimate cosmological parameters. We obtain a sample of 523 Type I

quasars with the UV/optical spectra and EW (C\,\textsc{iv}) measurements in the

redshift range of $2.3< z< 7.1$ including 38 objects at $5 < z < 7.1$, which

can be used to investigate the C\,\textsc{iv} Baldwin effect and determine

cosmological luminosity distance. The relation $EW(C\,\textsc{iv}) \propto

{(\lambda {L_\lambda })^\gamma }$ can be applied to check the inverse

correlation between the C\,\textsc{iv} EW and ${L_\lambda }$ of quasars and

give their distance. The data suggest that the EW of C\,\textsc{iv} is

inversely correlated with continuum monochromatic luminosities. On the other

hand, we also consider dividing the Type I quasars sample into various redshift

bins, which can be used to check if the C\,\textsc{iv} EW-luminosity relation

depends on the redshift. Finally, we apply a combination of Type I quasars and

SNIa Pantheon to test the property of dark energy concerning whether or not its

density deviates from the constant, and give the statistical results.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.09456v3

DEMNUni: comparing nonlinear power spectra prescriptions in the presence of massive neutrinos and dynamical dark energy by G. Parimbelli et al. on Monday 21 November

We provide an accurate comparison, against large cosmological $N$-body

simulations, of different prescriptions for modelling nonlinear matter power

spectra in the presence of massive neutrinos and dynamical dark energy. We test

the current most widely used approaches: fitting functions (HALOFIT and

HMcode), the halo-model reaction (ReACT) and emulators (baccoemu and

EuclidEmulator2). Focussing on redshifts $z\leq2$ and scales $k\lesssim 1 \

h/$Mpc (where the simulation mass resolution provides $\sim 1\%$ accuracy), we

find that HMcode and ReACT considerably improve over the HALOFIT prescriptions

of Smith and Takahashi (both combined with the Bird correction), with an

overall agreement of 2\% for all the cosmological scenarios considered.

Concerning emulators, we find that, especially at low redshifts,

EuclidEmulator2 remarkably agrees with the simulated spectra at $\lesssim 1\%$

level in scenarios with dynamical dark energy and massless neutrinos, reaching

a maximum difference of $\sim 2\%$ at $z=2$. baccoemu has a similar behaviour

as EuclidEmulator2, except for a couple of dark energy models. In cosmologies

with massive neutrinos, at $z=0$ all the nonlinear prescriptions improve their

agreement with respect to the massless neutrino case, except for the Bird and

TakaBird models which, however, are not tailored to $w_0$--$w_a$ models. At

$z>0$ we do not find a similar improvement when including massive neutrinos,

probably due to the lower impact of neutrino free-streaming at higher

redshifts; rather at $z=2$ EuclidEmulator2 exceeds $2\%$ agreement for some

dark energy equation of state. When considering ratios between the matter power

spectrum computed in a given cosmological model and its $\Lambda$CDM

counterpart, all the tested prescriptions agree with simulated data, at

sub-percent or percent level, depending on $z$. [ABRIDGED]

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2207.13677v2

Temporal evolution and rigidity dependence of the solar modulation lag of Galactic cosmic rays by Nicola Tomassetti et al. on Monday 21 November

When traveling in the heliosphere, Galactic cosmic rays (GCRs) are subjected

to the solar modulation effect, a quasiperiodical change of their intensity

caused by the 11-year cycle of solar activity. Here we investigate the

association of solar activity and cosmic radiation over five solar cycles, from

1965 to 2020, using a collection of multichannel data from neutron monitors,

space missions, and solar observatories. In particular, we focus on the time

lag between the monthly sunspot number and the GCR flux variations. We show

that the modulation lag is subjected to a 22-year periodical variation, ranging

from about 2 to 14 months and following the polarity cycle of the Sun's

magnetic field. We also show that the lag is remarkably decreasing with

increasing energy of the GCR particles. These results reflect the interplay of

basic physics phenomena that cause the GCR modulation effect: the drift motion

of charged particles in the interplanetary magnetic field, the latitudinal

dependence of the solar wind, the energy dependence of their residence time in

the heliosphere. Based on this interpretation, we end up with a global

effective formula for the modulation lag and testable predictions for the flux

evolution of cosmic particles and antiparticles over the solar cycle.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2210.05693v2

Decomposition of galactic X-ray emission with Phox: Contributions from hot gas and XRBs by Stephan Vladutescu-Zopp et al. on Monday 21 November

We provide a numerical framework with which spatially and spectrally accurate

representations of X-ray binary populations can be studied from hydrodynamical

cosmological simulations. We construct average spectra accounting for a hot gas

component and verify the emergence of observed scaling relations between galaxy

wide X-ray luminosity ($L_{X}$) and stellar mass ($M_{\star}$) as well as

star-formation rate (SFR). Using simulated galaxy halos extracted from the

$(48\,h^{-1} \mathrm{cMpc})^3$ volume of the Magneticum Pathfinder cosmological

simulations at $z = 0.07$ we generate mock spectra with the X-ray

photon-simulator Phox. We extend the Phox code to account for the stellar

component in the simulation and study the resulting contribution in composite

galactic spectra. Average X-ray luminosity functions are perfectly reproduced

up to the one-photon luminosity limit. Comparing our resulting

$L_{X}-\mathrm{SFR}-M_{\star}$ relation for X-ray binaries with recent

observations of field galaxies in the Virgo galaxy cluster we find significant

overlap. Invoking a metallicity dependent model for high-mass X-ray binaries

yields an anti-correlation between mass-weighted stellar metallicity and SFR

normalized luminosity. The spatial distribution of high-mass X-ray binaries

coincides with star-formation regions of simulated galaxies while low-mass

X-ray binaries follow the stellar mass surface density. X-ray binary emission

is the dominant contribution in the 2-10 keV band in the absence of an actively

accreting central super-massive black hole with 50% contribution in the 0.5-2

keV band rivaling the hot gas component. Our modelling remains consistent with

observations despite uncertainties connected to our approach. The predictive

power and easily extendable framework hold great value for future

investigations of galactic X-ray spectra.

arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2208.04975v2