|
Figueroa, D. G., Florio, A., Torrenti, F., & Valkenburg, W. (2021). The art of simulating the early universe. Part I. Integration techniques and canonical cases. J. Cosmol. Astropart. Phys., 04(4), 035–108pp.
Abstract: We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce the basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of O(delta t(2)) accuracy like staggered leapfrog and Verlet integration, to Runge-Kutta methods up to O(delta t(4)) accuracy, and the Yoshida and Gauss-Legendre higher-order integrators, accurate up to O(delta t(10)) We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in 3+1 dimensions, including the case of 'self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: i) Interacting scalar fields, ii) Abelian U(1) gauge theories, and iii) Non-Abelian SU(2) gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from O(delta t(2)) up to O(delta t(10)) For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. We note that all our algorithms for gauge theories always respect the Gauss constraint to machine precision, including when 'self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under SU(2) x U(1). We note that the present manuscript is meant to be part of the theoretical basis for the code CosmoLattice, a multi-purpose MPI-based package for simulating the non-linear evolution of field theories in an expanding universe, publicly available at http://www.cosrnolattice.net.
|
|
|
Bennett, J. J., Buldgen, G., de Salas, P. F., Drewes, M., Gariazzo, S., Pastor, S., et al. (2021). Towards a precision calculation of the effective number of neutrinos N-eff in the Standard Model. Part II. Neutrino decoupling in the presence of flavour oscillations and finite-temperature QED. J. Cosmol. Astropart. Phys., 04(4), 073–33pp.
Abstract: We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, N-eff(SM), that quantifies the cosmological neutrinoto-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to O(e(3)), where e is the elementary electric charge, and a full evaluation of the neutrino-neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark N(eff)(SM )value, through testing the value's dependence on (i) optional approximate modelling of the weak collision integrals, (ii) measurement errors in the physical parameters of the weak sector, and (iii) numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is N-eff(SM) 3.0440 +/- 0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (vertical bar delta N-eff vertical bar similar to 10(-4)), augmented by measurement errors in the solar mixing angle sin(2) theta(12) (vertical bar delta N-eff vertical bar similar to 10(-4)).
|
|
|
Valdes-Cortez, C., Mansour, I., Rivard, M. J., Ballester, F., Mainegra-Hing, E., Thomson, R. M., et al. (2021). A study of Type B uncertainties associated with the photoelectric effect in low-energy Monte Carlo simulations. Phys. Med. Biol., 66(10), 105014–14pp.
Abstract: Purpose. To estimate Type B uncertainties in absorbed-dose calculations arising from the different implementations in current state-of-the-art Monte Carlo (MC) codes of low-energy photon cross-sections (<200 keV). Methods. MC simulations are carried out using three codes widely used in the low-energy domain: PENELOPE-2018, EGSnrc, and MCNP. Three dosimetry-relevant quantities are considered: mass energy-absorption coefficients for water, air, graphite, and their respective ratios; absorbed dose; and photon-fluence spectra. The absorbed dose and the photon-fluence spectra are scored in a spherical water phantom of 15 cm radius. Benchmark simulations using similar cross-sections have been performed. The differences observed between these quantities when different cross-sections are considered are taken to be a good estimator for the corresponding Type B uncertainties. Results. A conservative Type B uncertainty for the absorbed dose (k = 2) of 1.2%-1.7% (<50 keV), 0.6%-1.2% (50-100 keV), and 0.3% (100-200 keV) is estimated. The photon-fluence spectrum does not present clinically relevant differences that merit considering additional Type B uncertainties except for energies below 25 keV, where a Type B uncertainty of 0.5% is obtained. Below 30 keV, mass energy-absorption coefficients show Type B uncertainties (k = 2) of about 1.5% (water and air), and 2% (graphite), diminishing in all materials for larger energies and reaching values about 1% (40-50 keV) and 0.5% (50-75 keV). With respect to their ratios, the only significant Type B uncertainties are observed in the case of the water-to-graphite ratio for energies below 30 keV, being about 0.7% (k = 2). Conclusions. In contrast with the intermediate (about 500 keV) or high (about 1 MeV) energy domains, Type B uncertainties due to the different cross-sections implementation cannot be considered subdominant with respect to Type A uncertainties or even to other sources of Type B uncertainties (tally volume averaging, manufacturing tolerances, etc). Therefore, the values reported here should be accommodated within the uncertainty budget in low-energy photon dosimetry studies.
|
|
|
LHCb Collaboration(Aaij, R. et al), Henry, L., Jashal, B. K., Martinez-Vidal, F., Oyanguren, A., Remon Alepuz, C., et al. (2021). Search for the doubly heavy baryons Omega(0)(bc) and Xi(0)(bc) decaying to Lambda(+)(c)pi(-) and Xi(+)(c)pi-. Chin. Phys. C, 45(9), 093002–12pp.
Abstract: The first search for the doubly heavy Omega(0)(bc) baryon and a search for the Xi(0)(bc) baryon are performed using collision data collected via the experiment from 2016 to 2018 at a centre-of-mass energy of, corresponding to an integrated luminosity of 5.2 fb(-1). The baryons are reconstructed via their decays to Lambda(+)(-)(c)(pi) and Xi(+)(c)pi(-). No significant excess is found for invariant masses between 6700 and 7300 MeV/c(2), in a rapidity range from 2.0 to 4.5 and a transverse momentum range from 2 to 20 MeV/c. Upper limits are set on the ratio of the Omega(0)(bc) and Xi(0)(bc) production cross-section times the branching fraction to Lambda(+)(c)pi(-)(Xi(+)(c)pi(-)) relative to that of the Lambda(0)(b)(Xi(0)(b)) baryon, for different lifetime hypotheses, at 95% confidence level. The upper limits range from 0.5x10(-4) to 2.5x10(-4) for the Omega(0)(bc) -> Lambda(+)(c)pi(-) (Xi(0)(bc) -> Lambda(+)(c)pi(-)) decay, pending on the considered mass and lifetime of the Omega(0)(bc) (Xi(0)(bc)) baryon.
|
|
|
Barenboim, G., Blinov, N., & Stebbins, A. (2021). Smallest remnants of early matter domination. J. Cosmol. Astropart. Phys., 12(12), 026–50pp.
Abstract: The evolution of the universe prior to Big Bang Nucleosynthesis could have gone through a phase of early matter domination which enhanced the growth of small-scale dark matter structure. If this period was long enough, self-gravitating objects formed prior to reheating. We study the evolution of these dense early halos through reheating. At the end of early matter domination, the early halos undergo rapid expansion and eventually eject their matter. We find that this process washes out structure on scales much larger than naively expected from the size of the original halos. We compute the density profiles of the early halo remnants and use them to construct late-time power spectra that include these non-linear effects. We evolve the resulting power spectrum to estimate the properties of microhalos that would form after matter-radiation equality. Surprisingly, cosmologies with a short period of early matter domination lead to an earlier onset of microhalo formation compared to those with a long period. In either case, dark matter structure formation begins much earlier than in the standard cosmology, with most dark matter bound in microhalos in the late universe.
|
|