Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates

Numerical simulations in cosmology require trade-offs between volume, resolution and run-time that limit the volume of the Universe that can be simulated, leading to sample variance in predictions of ensemble-average quantities such as the power spectrum or correlation function(s). Sample variance is particularly acute at large scales, which is also where analytic techniques can be highly reliable. This provides an opportunity to combine analytic and numerical techniques in a principled way to improve the dynamic range and reliability of predictions for clustering statistics. In this paper we extend the technique of Zel'dovich control variates, previously demonstrated for 2-point functions in real space, to reduce the sample variance in measurements of 2-point statistics of biased tracers in redshift space. We demonstrate that with this technique, we can reduce the sample variance of these statistics down to their shot-noise limit out to $k \sim 0.2\, h\rm Mpc^{-1}$. This allows a better matching with perturbative models and improved predictions for the clustering of e.g.~quasars, galaxies and neutral Hydrogen measured in spectroscopic redshift surveys at very modest computational expense. We discuss the implementation of ZCV, give some examples and provide forecasts for the efficacy of the method under various conditions.Comment: 17 pages main text, 5 pages of appendices, 9 figure

Cosmological Analysis of Three-Dimensional BOSS Galaxy Clustering and Planck CMB Lensing Cross Correlations via Lagrangian Perturbation Theory

We present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain $\Omega_\mathrm{m} = 0.303\pm 0.008$, $H_0 = 69.21\pm 0.78$ and $\sigma_8 = 0.743\pm 0.043$. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to $\Omega_m = 0.300 \pm 0.008$, $H_0 = 69.21 \pm 0.77$ and $\sigma_8 = 0.707 \pm 0.035$, in tension with CMB and cosmic shear constraints. The combination of $\Omega_m$ and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $\sigma_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($z\simeq 0.38$) and at large angular scales, which show a $20\%$ deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.Comment: 46 pages, 15 figures, updated to match version accepted by JCA

The cosmology dependence of galaxy clustering and lensing from a hybrid NN-body-perturbation theory model

We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from $N$-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of $N$-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these nonlinear spectra from redshifts $z = 0$ to $z=2$. We quantify the accuracy of our emulation procedure, which is sub-per cent at $k=1\, h {\rm Mpc}^{-1}$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift $z\simeq 0.4$ to scales of $k_{\rm max} \approx 0.6\, h\mathrm{Mpc}^{-1}$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy--galaxy lensing analyses with data drawn from an independent $N$-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming datasets from galaxy surveys.Comment: 19 pages, 17 figures. Updated to reflect the journal version. Code available at https://github.com/kokron/anz

Priors on red galaxy stochasticity from hybrid effective field theory

We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use Halo Occupation Distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with datasets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of hybrid EFT methods to characterize models of the galaxy--halo connection at the field level, revealing new connections between once-disparate modelling frameworks.Comment: 16 pages, 10 figures. Revised version accepted to MNRAS. Revisions include a new appendix on the covariance of the field-level bias estimato

Hepatic Homeostasis of Metal Ions Following Acute Repeated Stress Exposure in Rats.

Essential metals such as copper, iron, and zinc are cofactors in various biological processes including oxygen utilisation, cell growth, and biomolecular synthesis. The homeostasis of these essential metals is carefully controlled through a system of protein transporters involved in the uptake, storage, and secretion. Some metal ions can be transformed by processes including reduction/oxidation (redox) reactions, and correspondingly, the breakdown of metal ion homeostasis can lead to formation of reactive oxygen and nitrogen species. We have previously demonstrated rapid biochemical responses to stress involving alterations in the redox state to generate free radicals and the resultant oxidative stress. However, the effects of stress on redox-active metals including iron and copper and redox-inert zinc have not been well characterised. Therefore, this study aims to examine the changes in these essential metals following exposure to short-term repeated stress, and to further elucidate the alterations in metal homeostasis through expression analysis of different metal transporters. Outbred male Wistar rats were exposed to unrestrained (control), 1 day, or 3 days of 6 h restraint stress (n = 8 per group). After the respective stress treatment, blood and liver samples were collected for the analysis of biometal concentrations and relative gene expression of metal transporter and binding proteins. Exposure to repeated restraint stress was highly effective in causing hepatic redox imbalance. Stress was also shown to induce hepatic metal redistribution, while modulating the mRNA levels of key metal transporters. Overall, this study is the first to characterise the gene expression profile of metal homeostasis following stress and provide insight into the changes occurring prior to the onset of chronic stress conditions

Aemulus ν\nu: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos

We present the Aemulus $\nu$ simulations: a suite of 150 $(1.05 h^{-1}\rm Gpc)^3$ $N$-body simulations with a mass resolution of $3.51\times 10^{10} \frac{\Omega_{cb}}{0.3} ~ h^{-1} M_{\odot}$ in a $w\nu$CDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in $\sigma_8$ to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to $0.5\, \rm eV$, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of $10-10^5$ depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve $\le 1\%$ accuracy for $k\le 1\, h\,\rm Mpc^{-1}$ and $0\le z \le 3$, and $\le 2\%$ accuracy for $k\le 4\, h\,\rm Mpc^{-1}$ for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come.Comment: 37 pages, 15 figures, matching version accepted by JCA

Restraint Stress Alters Expression of Glucocorticoid Bioavailability Mediators, Suppresses Nrf2, and Promotes Oxidative Stress in Liver Tissue.

Hepatic glutathione synthesis and antioxidant protection are critically important for efficient detoxification processes in response to metabolic challenges. However, this biosynthetic pathway, regulated by nuclear factor (erythroid-derived 2)-like 2 (Nrf2), previously demonstrated paradoxical repression following exposure to glucocorticoid stress hormones in cultured hepatic cells. Therefore, the present study used an in vivo model of sub-acute psychological stress to investigate the relationship between hepatic corticosteroid regulation and antioxidant systems. Male Wistar rats were kept under control conditions or subjected to six hours of restraint stress applied for 1 or 3 days (n = 8 per group) after which the liver was isolated for assays of oxidative/nitrosative status and expression of corticosteroid regulatory and Nrf2-antioxidant response element pathway members. A single stress exposure produced a significant increase in the expression of corticosterone reactivator, 11-beta-hydroxysteroid dehydrogenase 1 (11β-Hsd1), while the 11β-Hsd2 isozyme and corticosteroid-binding globulin were down-regulated following stress, indicative of an elevated availability of active corticosterone. Exposure to restraint significantly decreased hepatic concentrations of total cysteine thiols and the antioxidant reduced glutathione on Day 1 and increased 3-nitrotyrosinated and carbonylated proteins on Day 3, suggestive of oxidative/nitrosative stress in the liver following stress exposure. Conversely, there was a sustained down-regulation of Nrf2 mRNA and protein in addition to significant reductions in downstream glutamate-cysteine ligase catalytic subunit (Gclc), the rate-limiting enzyme in glutathione synthesis, on Day 1 and 3 of stress treatment. Interestingly, other antioxidant genes including superoxide dismutase 1 and 2, and glutathione peroxidase 4 were significantly up-regulated following an episode of restraint stress. In conclusion, the results of the present study indicate that increased expression of 11β-Hsd1, indicative of elevated tissue glucocorticoid concentrations, may impair the Nrf2-dependent antioxidant response

Impacts of the Tropical Pacific/Indian Oceans on the Seasonal Cycle of the West African Monsoon

The current consensus is that drought has developed in the Sahel during the second half of the twentieth century as a result of remote effects of oceanic anomalies amplified by local land–atmosphere interactions. This paper focuses on the impacts of oceanic anomalies upon West African climate and specifically aims to identify those from SST anomalies in the Pacific/Indian Oceans during spring and summer seasons, when they were significant. Idealized sensitivity experiments are performed with four atmospheric general circulation models (AGCMs). The prescribed SST patterns used in the AGCMs are based on the leading mode of covariability between SST anomalies over the Pacific/Indian Oceans and summer rainfall over West Africa. The results show that such oceanic anomalies in the Pacific/Indian Ocean lead to a northward shift of an anomalous dry belt from the Gulf of Guinea to the Sahel as the season advances. In the Sahel, the magnitude of rainfall anomalies is comparable to that obtained by other authors using SST anomalies confined to the proximity of the Atlantic Ocean. The mechanism connecting the Pacific/Indian SST anomalies with West African rainfall has a strong seasonal cycle. In spring (May and June), anomalous subsidence develops over both the Maritime Continent and the equatorial Atlantic in response to the enhanced equatorial heating. Precipitation increases over continental West Africa in association with stronger zonal convergence of moisture. In addition, precipitation decreases over the Gulf of Guinea. During the monsoon peak (July and August), the SST anomalies move westward over the equatorial Pacific and the two regions where subsidence occurred earlier in the seasons merge over West Africa. The monsoon weakens and rainfall decreases over the Sahel, especially in August.Peer reviewe

Severe early onset preeclampsia: short and long term clinical, psychosocial and biochemical aspects

Preeclampsia is a pregnancy specific disorder commonly defined as de novo hypertension and proteinuria after 20 weeks gestational age. It occurs in approximately 3-5% of pregnancies and it is still a major cause of both foetal and maternal morbidity and mortality worldwide1. As extensive research has not yet elucidated the aetiology of preeclampsia, there are no rational preventive or therapeutic interventions available. The only rational treatment is delivery, which benefits the mother but is not in the interest of the foetus, if remote from term. Early onset preeclampsia (<32 weeks’ gestational age) occurs in less than 1% of pregnancies. It is, however often associated with maternal morbidity as the risk of progression to severe maternal disease is inversely related with gestational age at onset2. Resulting prematurity is therefore the main cause of neonatal mortality and morbidity in patients with severe preeclampsia3. Although the discussion is ongoing, perinatal survival is suggested to be increased in patients with preterm preeclampsia by expectant, non-interventional management. This temporising treatment option to lengthen pregnancy includes the use of antihypertensive medication to control hypertension, magnesium sulphate to prevent eclampsia and corticosteroids to enhance foetal lung maturity4. With optimal maternal haemodynamic status and reassuring foetal condition this results on average in an extension of 2 weeks. Prolongation of these pregnancies is a great challenge for clinicians to balance between potential maternal risks on one the eve hand and possible foetal benefits on the other. Clinical controversies regarding prolongation of preterm preeclamptic pregnancies still exist – also taking into account that preeclampsia is the leading cause of maternal mortality in the Netherlands5 - a debate which is even more pronounced in very preterm pregnancies with questionable foetal viability6-9. Do maternal risks of prolongation of these very early pregnancies outweigh the chances of neonatal survival? Counselling of women with very early onset preeclampsia not only comprises of knowledge of the outcome of those particular pregnancies, but also knowledge of outcomes of future pregnancies of these women is of major clinical importance. This thesis opens with a review of the literature on identifiable risk factors of preeclampsia

Measurements of the pp → ZZ production cross section and the Z → 4ℓ branching fraction, and constraints on anomalous triple gauge couplings at √s = 13 TeV

Four-lepton production in proton-proton collisions, pp -> (Z/gamma*)(Z/gamma*) -> 4l, where l = e or mu, is studied at a center-of-mass energy of 13 TeV with the CMS detector at the LHC. The data sample corresponds to an integrated luminosity of 35.9 fb(-1). The ZZ production cross section, sigma(pp -> ZZ) = 17.2 +/- 0.5 (stat) +/- 0.7 (syst) +/- 0.4 (theo) +/- 0.4 (lumi) pb, measured using events with two opposite-sign, same-flavor lepton pairs produced in the mass region 60 4l) = 4.83(-0.22)(+0.23) (stat)(-0.29)(+0.32) (syst) +/- 0.08 (theo) +/- 0.12(lumi) x 10(-6) for events with a four-lepton invariant mass in the range 80 4GeV for all opposite-sign, same-flavor lepton pairs. The results agree with standard model predictions. The invariant mass distribution of the four-lepton system is used to set limits on anomalous ZZZ and ZZ. couplings at 95% confidence level: -0.0012 < f(4)(Z) < 0.0010, -0.0010 < f(5)(Z) < 0.0013, -0.0012 < f(4)(gamma) < 0.0013, -0.0012 < f(5)(gamma) < 0.0013