Showing 1–50 of 56 results for author: Appleby, S

Probing cosmology via the clustering of critical points

Authors: Junsup Shim, Christophe Pichon, Dmitri Pogosyan, Stephen Appleby, Corentin Cadiou, Juhan Kim, Katarina Kraljic, Changbom Park

Abstract: Exclusion zones in the cross-correlations between critical points (peak-void, peak-wall, filament-wall, filament-void) of the density field define quasi-standard rulers that can be used to constrain dark matter and dark energy cosmological parameters. The average size of the exclusion zone is found to scale linearly with the typical distance between extrema. The latter changes as a function of the… ▽ More Exclusion zones in the cross-correlations between critical points (peak-void, peak-wall, filament-wall, filament-void) of the density field define quasi-standard rulers that can be used to constrain dark matter and dark energy cosmological parameters. The average size of the exclusion zone is found to scale linearly with the typical distance between extrema. The latter changes as a function of the matter content of the universe in a predictable manner, but its comoving size remains essentially constant in the linear regime of structure growth on large scales, unless the incorrect cosmology is assumed in the redshift-distance relation. This can be used to constrain the dark energy parameters when considering a survey that scans a range of redshifts. The precision of the parameter estimation is assessed using a set of cosmological simulations, and is found to be a 4$σ$ detection of a change in matter content of 5%, or about 3.8$σ$ detection of 50% shift in the dark energy parameter using a full sky survey up to redshift 0.5. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 12 pages (including a 3-page appendix, 9 figures, 2 tables. Submitted to MNRAS

The effects of non-linearity on the growth rate constraint from velocity correlation functions

Authors: Motonari Tonegawa, Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: The two-point statistics of the cosmic velocity field, measured from galaxy peculiar velocity (PV) surveys, can be used as a dynamical probe to constrain the growth rate of large-scale structures in the universe. Most works use the statistics on scales down to a few tens of Megaparsecs, while using a theoretical template based on the linear theory. In addition, while the cosmic velocity is volume-… ▽ More The two-point statistics of the cosmic velocity field, measured from galaxy peculiar velocity (PV) surveys, can be used as a dynamical probe to constrain the growth rate of large-scale structures in the universe. Most works use the statistics on scales down to a few tens of Megaparsecs, while using a theoretical template based on the linear theory. In addition, while the cosmic velocity is volume-weighted, the observable line-of-sight velocity two-point correlation is density-weighted, as sampled by galaxies, and therefore the density-velocity correlation term also contributes, which has often been neglected. These effects are fourth order in powers of the linear density fluctuation $δ_{\rm L}^4$, compared to $δ_{\rm L}^2$ of the linear velocity correlation function, and have the opposite sign. We present these terms up to $δ_{\rm L}^4$ in real space based on the standard perturbation theory, and investigate the effect of non-linearity and the density-velocity contribution on the inferred growth rate $fσ_8$, using $N$-body simulations. We find that for a next-generation PV survey of volume $\sim {\cal O}(500 \, h^{-1} \, {\rm Mpc})^3$, these effects amount to a shift of $fσ_8$ by $\sim 10$ per cent and is comparable to the forecasted statistical error when the minimum scale used for parameter estimation is $r_{\rm min} = 20 \, h^{-1} \, {\rm Mpc}$. △ Less

Submitted 6 March, 2024; v1 submitted 25 September, 2023; originally announced September 2023.

Comments: 17 pages, 6 figures. Minor corrections to non-linear calculation in section 2.2, conclusions unchanged. Accepted for publication in MNRAS

The One-hundred-deg^2 DECam Imaging in Narrowbands (ODIN): Survey Design and Science Goals

Authors: Kyoung-Soo Lee, Eric Gawiser, Changbom Park, Yujin Yang, Francisco Valdes, Dustin Lang, Vandana Ramakrishnan, Byeongha Moon, Nicole Firestone, Stephen Appleby, Maria Celeste Artale, Moira Andrews, Franz E. Bauer, Barbara Benda, Adam Broussard, Yi-Kuan Chiang, Robin Ciardullo, Arjun Dey, Rameen Farooq, Caryl Gronwall, Lucia Guaita, Yun Huang, Ho Seong Hwang, Sanghyeok Im, Woong-Seob Jeong , et al. (17 additional authors not shown)

Abstract: We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.… ▽ More We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.2, 7.4, and 9.8 nm, corresponding to Lya at z=2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broad-band data from Hyper Suprime-Cam, DECam, and in the future, LSST, the ODIN narrow-band images will enable the selection of over 100,000 Lya-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrow-band images enable detection of highly extended Lya blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The two narrow-band filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z=0.34, Lya and He II 1640 at z=3.1, and Lyman continuum plus Lya at z=4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby star-forming galaxies resolved into numerous [O III] and [S II] emitting regions. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 21 pages, 8 figures, submitted to Astrophysical Journal

Tadpole Cosmology: Milne Solution as a Cosmological Constant Hideout

Authors: Stephen Appleby, Reginald Christian Bernardo

Abstract: Dynamical cancellation frameworks present a potential means of mitigating the effect of a large vacuum energy, that would otherwise ruin the late-time, low energy dynamics of the Universe. Certain models in the literature, such as the Fab Four and Well Tempering, realize this idea by introducing some degeneracy in the dynamical equations. In this paper, we introduce a third potential route to self… ▽ More Dynamical cancellation frameworks present a potential means of mitigating the effect of a large vacuum energy, that would otherwise ruin the late-time, low energy dynamics of the Universe. Certain models in the literature, such as the Fab Four and Well Tempering, realize this idea by introducing some degeneracy in the dynamical equations. In this paper, we introduce a third potential route to self-tuning, and infer the existence of a new, exact Milne solution in the simplest tadpole plus cubic-Galileon scalar-tensor theory. We study the dynamics of the scalar field and metric in the vicinity of the Milne coordinate singularity, and find that the vacuum solution belongs to a more general family of Milne-like metrics. By numerically evolving the field equations for a range of initial conditions, we show that the Milne solution is not an attractor, and varying the initial scalar field data can lead to completely different asymptotic states; exponential growth of the scale factor, a static non-spatially flat metric or a severe finite-time instability in the scalar field and metric. We generalise the Milne solution to a class of FLRW spacetimes, finding that the tadpole-cubic Galileon model admits perfect-fluid-like solutions in the presence of matter. Finally, we present a second Horndeski model which also admits an exact Milne solution, hinting at the existence of a larger undiscovered model space containing vacuum-energy-screened solutions. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 20 pages + refs, 4 figures, comments welcome

Journal ref: JCAP 12 (2023) 003

Cosmological Parameter Constraints from the SDSS Density and Momentum Power Spectra

Authors: Stephen Appleby, Motonari Tonegawa, Changbom Park, Sungwook E. Hong, Juhan Kim, Yongmin Yoon

Abstract: We extract the galaxy density and momentum power spectra from a subset of early-type galaxies in the SDSS DR7 main galaxy catalog. Using galaxy distance information inferred from the improved fundamental plane described in \citet{Yoon_2020}, we reconstruct the peculiar velocities of the galaxies and generate number density and density-weighted velocity fields, from which we extract the galaxy dens… ▽ More We extract the galaxy density and momentum power spectra from a subset of early-type galaxies in the SDSS DR7 main galaxy catalog. Using galaxy distance information inferred from the improved fundamental plane described in \citet{Yoon_2020}, we reconstruct the peculiar velocities of the galaxies and generate number density and density-weighted velocity fields, from which we extract the galaxy density and momentum power spectra. We compare the measured values to the theoretical expectation of the same statistics, assuming an input $Λ$CDM model and using a third-order perturbative expansion. After validating our analysis pipeline with a series of mock data sets, we apply our methodology to the SDSS data and arrive at constraints $fσ_{8} = 0.471_{-0.080}^{+0.077}$ and $b_{1}σ_{8} = 0.920_{-0.070}^{+0.070}$ at a mean redshift $\bar{z} = 0.04$. Our result is consistent with the Planck cosmological best fit parameters for the $Λ$CDM model. The momentum power spectrum is found to be strongly contaminated by small scale velocity dispersion, which suppresses power by $\sim {\cal O}(30\%)$ on intermediate scales $k \sim 0.05 \, h \, {\rm Mpc}^{-1}$. △ Less

Submitted 2 October, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 15 figures, 4 tables, accepted for ApJ

Tomographic Alcock-Paczynski Test with Redshift-Space Correlation Function: Evidence for the Dark Energy Equation of State Parameter w>-1

Authors: Fuyu Dong, Changbom Park, Sungwook E. Hong, Juhan Kim, Ho Seong Hwang, Hyunbae Park, Stephen Appleby

Abstract: The apparent shape of galaxy clustering depends on the adopted cosmology used to convert observed redshift to comoving distance, the $r(z)$ relation, as it changes the line elements along and across the line of sight differently. The Alcock-Paczyński (AP) test exploits this property to constrain the expansion history of the universe. We present an extensive review of past studies on the AP test. W… ▽ More The apparent shape of galaxy clustering depends on the adopted cosmology used to convert observed redshift to comoving distance, the $r(z)$ relation, as it changes the line elements along and across the line of sight differently. The Alcock-Paczyński (AP) test exploits this property to constrain the expansion history of the universe. We present an extensive review of past studies on the AP test. We adopt an extended AP test method introduced by Park et al. (2019), which uses the full shape of redshift-space two-point correlation function (CF) as the standard shape, and apply it to the SDSS DR7, BOSS, and eBOSS LRG samples covering the redshift range up to $z=0.8$.We calibrate the test against the nonlinear cosmology-dependent systematic evolution of the CF shape using the Multiverse simulations. We focus on examining whether or not the flat $Λ$CDM `concordance' model is consistent with observation. We constrain the flat $w$CDM model to have $w=-0.892_{-0.050}^{+0.045}$ and $Ω_m=0.282_{-0.023}^{+0.024}$ from our AP test alone, which is significantly tighter than the constraints from the BAO or SNe I$a$ methods by a factor of 3 - 6. When the AP test result is combined with the recent BAO and SNe I$a$ results, we obtain $w=-0.903_{-0.023}^{+0.023}$ and $Ω_m=0.285_{-0.009}^{+0.014}$. This puts a strong tension with the flat $Λ$CDM model with $w=-1$ at $4.2σ$ level. Consistency with $w=-1$ is obtained only when the Planck CMB observation is combined. It remains to see if this tension between observations of galaxy distribution at low redshifts and CMB anisotropy at the decoupling epoch becomes greater in the future studies and leads us to a new paradigm of cosmology. △ Less

Submitted 29 April, 2023; originally announced May 2023.

Comments: 21 pages, 11 figures, accepted by ApJ

Mapping Circumgalactic Medium Observations to Theory Using Machine Learning

Authors: Sarah Appleby, Romeel Davé, Daniele Sorini, Christopher Lovell, Kevin Lo

Abstract: We present a random forest framework for predicting circumgalactic medium (CGM) physical conditions from quasar absorption line observables, trained on a sample of Voigt profile-fit synthetic absorbers from the Simba cosmological simulation. Traditionally, extracting physical conditions from CGM absorber observations involves simplifying assumptions such as uniform single-phase clouds, but by usin… ▽ More We present a random forest framework for predicting circumgalactic medium (CGM) physical conditions from quasar absorption line observables, trained on a sample of Voigt profile-fit synthetic absorbers from the Simba cosmological simulation. Traditionally, extracting physical conditions from CGM absorber observations involves simplifying assumptions such as uniform single-phase clouds, but by using a cosmological simulation we bypass such assumptions to better capture the complex relationship between CGM observables and underlying gas conditions. We train random forest models on synthetic spectra for HI and selected metal lines around galaxies across a range of star formation rates, stellar masses, and impact parameters, to predict absorber overdensities, temperatures, and metallicities. The models reproduce the true values from Simba well, with normalised transverse standard deviations of $0.50-0.54$ dex in overdensity, $0.32-0.54$ dex in temperature, and $0.49-0.53$ dex in metallicity predicted from metal lines (not HI), across all ions. Examining the feature importance, the random forest indicates that the overdensity is most informed by the absorber column density, the temperature is driven by the line width, and the metallicity is most sensitive to the specific star formation rate. Alternatively examining feature importance by removing one observable at a time, the overdensity and metallicity appear to be more driven by the impact parameter. We introduce a normalising flow approach in order to ensure the scatter in the true physical conditions is accurately spanned by the network. The trained models are available online. △ Less

Submitted 23 July, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 16 pages, 14 figures. Accepted for publication in MNRAS

Minkowski Tensors in Redshift Space -- Beyond the Plane Parallel Approximation

Authors: Stephen Appleby, Joby P. Kochappan, Pravabati Chingangbam, Changbom Park

Abstract: The Minkowski tensors (MTs) can be used to probe anisotropic signals in a field, and are well suited for measuring the redshift space distortion (RSD) signal in large scale structure catalogs. We consider how the linear RSD signal can be extracted from a field without resorting to the plane parallel approximation. A spherically redshift space distorted field is both anisotropic and inhomogeneous.… ▽ More The Minkowski tensors (MTs) can be used to probe anisotropic signals in a field, and are well suited for measuring the redshift space distortion (RSD) signal in large scale structure catalogs. We consider how the linear RSD signal can be extracted from a field without resorting to the plane parallel approximation. A spherically redshift space distorted field is both anisotropic and inhomogeneous. We derive expressions for the two point correlation functions that elucidate the inhomogeneity, and then explain how the breakdown of homogeneity impacts the volume and ensemble averages of the tensor Minkowski functionals. We construct the ensemble average of these quantities in curvilinear coordinates and show that the ensemble and volume averages can be approximately equated, but this depends on our choice of definition of the volume average of a tensor and the radial distance between the observer and field. We then extract the tensor Minkowski functionals from spherically redshift space distorted, Gaussian random fields and gravitationally evolved dark matter density fields at $z=0$ to test if we can successfully measure the Kaiser RSD signal. For the dark matter field we find a significant, $\sim 10\%$ anomalous signal in the MT component parallel to the line of sight that is present even on large scales $R_{\rm G} \gtrsim 15 \, {\rm Mpc}$, in addition to the Kaiser effect. This is due to the line of sight component of the MT being significantly contaminated by the Finger of God effect, which can be approximately modelled by an additional damping term in the cumulants. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 35 pages, 12 figures

The Physical Nature of Circumgalactic Medium Absorbers in Simba

Authors: Sarah Appleby, Romeel Davé, Daniele Sorini, Weiguang Cui, Jacob Christiansen

Abstract: We study the nature of the low-redshift CGM in the Simba cosmological simulations as traced by ultraviolet absorption lines around galaxies in bins of stellar mass ($M_\star>10^{10}M_\odot$) for star-forming, green valley and quenched galaxies at impact parameters $r_\perp\leq 1.25r_{200}$. We generate synthetic spectra for HI, MgII, CII, SiIII, CIV, and OVI, fit Voigt profiles to obtain line prop… ▽ More We study the nature of the low-redshift CGM in the Simba cosmological simulations as traced by ultraviolet absorption lines around galaxies in bins of stellar mass ($M_\star>10^{10}M_\odot$) for star-forming, green valley and quenched galaxies at impact parameters $r_\perp\leq 1.25r_{200}$. We generate synthetic spectra for HI, MgII, CII, SiIII, CIV, and OVI, fit Voigt profiles to obtain line properties, and estimate the density, temperature, and metallicity of the absorbing gas. We find that CGM absorbers are most abundant around star forming galaxies with $M_\star < 10^{11}M_\odot$, while the abundance of green valley galaxies show similar behaviour to those of quenched galaxies, suggesting that the CGM "quenches" before star formation ceases. HI absorbing gas exists across a broad range of cosmic phases (condensed gas, diffuse gas, hot halo gas and Warm-Hot Intergalactic Medium), while essentially all low-ionisation metal absorption arises from condensed gas. OVI absorbers are split between hot halo gas and the WHIM. The fraction of collisionally ionised CGM absorbers is $\sim 25-55\%$ for CIV and $\sim 80-95\%$ for OVI, depending on stellar mass and impact parameter. In general, the highest column density absorption features for each ion arise from dense gas. Satellite gas, defined as that within $10r_{1/2,\star},$ contributes $\sim 3\%$ of overall HI absorption but $\sim 30\%$ of MgII absorption, with the fraction from satellites decreasing with increasing ion excitation energy. △ Less

Submitted 5 January, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 24 pages, 15 figures; accepted for publication in MNRAS

Tadpole Cosmology: Self Tuning Without Degeneracy

Authors: Stephen Appleby, Reginald Christian Bernardo

Abstract: Degeneracy is a method to accommodate exact, low energy vacuum states in scalar-tensor gravitational models despite the presence of an arbitrarily large vacuum energy. However, this approach requires very particular combinations of scalar field and metric couplings in the Lagrangian. In this work we study departures from the restrictive degeneracy condition -- starting from a fiducial model contai… ▽ More Degeneracy is a method to accommodate exact, low energy vacuum states in scalar-tensor gravitational models despite the presence of an arbitrarily large vacuum energy. However, this approach requires very particular combinations of scalar field and metric couplings in the Lagrangian. In this work we study departures from the restrictive degeneracy condition -- starting from a fiducial model containing an exact Minkowski space solution, we break the degeneracy condition in numerous simple ways to test if the resulting models maintain certain key features -- specifically the dynamical cancellation of a large vacuum energy by the scalar field and the existence of a low energy vacuum state. We highlight the role the tadpole plays in eliminating the fixed points of the dynamical system, generically rendering both the scalar field and metric time dependent. Our results indicate that when violating the degeneracy condition but preserving shift symmetry, the metric maintains an asymptotic Minkowski state, irrespective of the presence of the cosmological constant. In contrast, when shift symmetry is also broken the asymptotic behaviour can radically alter. Regardless, the non-degenerate models in this work share an attractive quality; harboring low energy, late-time asymptotic states that are independent of the vacuum energy. The tadpole allows for a broader class of non-degenerate, self-tuning models than was previously realized. △ Less

Submitted 14 June, 2022; v1 submitted 17 February, 2022; originally announced February 2022.

Comments: 19 pages, 4 figures, 1 table, v2: to appear in JCAP

Journal ref: JCAP 07 (2022) 035

How baryons affect halos and large-scale structure: a unified picture from the Simba simulation

Authors: Daniele Sorini, Romeel Dave, Weiguang Cui, Sarah Appleby

Abstract: Using the state-of-the-art suite of hydrodynamic simulations Simba, as well as its dark-matter-only counterpart, we study the impact of the presence of baryons and of different stellar/AGN feedback mechanisms on large-scale structure, halo density profiles, and on the abundance of different baryonic phases within halos and in the intergalactic medium (IGM). The unified picture that emerges from ou… ▽ More Using the state-of-the-art suite of hydrodynamic simulations Simba, as well as its dark-matter-only counterpart, we study the impact of the presence of baryons and of different stellar/AGN feedback mechanisms on large-scale structure, halo density profiles, and on the abundance of different baryonic phases within halos and in the intergalactic medium (IGM). The unified picture that emerges from our analysis is that the main physical drivers shaping the distribution of matter at all scales are star formation-driven galactic outflows at $z>2$ for lower mass halos and AGN jets at $z<2$ in higher mass halos. Feedback suppresses the baryon mass function with time relative to the halo mass function, and it even impacts the halo mass function itself at the ~20% level, particularly evacuating the centres and enhancing dark matter just outside halos. At early epochs baryons pile up in the centres of halos, but by late epochs and particularly in massive systems gas has mostly been evacuated from within the inner halo. AGN jets are so efficient at such evacuation that at low redshifts the baryon fraction within $\sim 10^{12}-10^{13} \, \rm M_{\odot}$ halos is only 25% of the cosmic baryon fraction, mostly in stars. The baryon fraction enclosed in a sphere around such halos approaches the cosmic value $Ω_{\rm b}/Ω_{\rm m}$ only at 10-20 virial radii. As a result, 87% of the baryonic mass in the Universe lies in the IGM at $z=0$, with 67% being in the form of warm-hot IGM ($T>10^5 \, \rm K$). △ Less

Submitted 16 March, 2023; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: Published in the Monthly Notices of the Royal Astronomical Society (MNRAS). This is a pre-copyedited, author-produced version of the article accepted for publication in MNRAS following peer review. The version of record (MNRAS, Vol. 516, Issue 1, pp. 883-906) is available online at: DOI 10.1093/mnras/stac2214 . Minor updates to text and some figures wrt previous version. Conclusions unchanged

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 1, pp.883-906, October 2022

Minkowski Functionals of SDSS-III BOSS : Hints of Possible Anisotropy in the Density Field?

Authors: Stephen Appleby, Changbom Park, Pratyush Pranav, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim, Thomas Buchert

Abstract: We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data we generate a s… ▽ More We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data we generate a set of four density fields in three dimensions corresponding to the northern and southern skies of LOWZ and CMASS catalogs, smoothing over large scales such that the field is perturbatively non-Gaussian. We extract the Minkowski functionals from each data set separately, and measure their shapes and amplitudes by fitting a Hermite polynomial expansion. For the shape parameter of the Minkowski functional curves $a_0$, that is related to the bispectrum of the field, we find that the LOWZ-South data presents a systematically lower value of $a_0 = -0.080 \pm 0.040$ than its northern sky counterpart $a_0 = 0.032 \pm 0.024$. Although the significance of this discrepancy is low, it potentially indicates some systematics in the data or that the matter density field exhibits anisotropy at low redshift. By assuming a standard isotropic flat $Λ$CDM cosmology, the amplitudes of Minkowski functionals from the combination of northern and southern sky data give the constraints $Ω_{\rm c} h^2 n_{\rm s} = 0.110 \pm 0.006$ and $0.111 \pm 0.008$ for CMASS and LOWZ, respectively, which is in agreement with the Planck $Λ$CDM best-fit $Ω_{\rm c}h^{2} n_{\rm s} = 0.116 \pm 0.001$. △ Less

Submitted 12 October, 2021; originally announced October 2021.

The geometrical meaning of statistical isotropy of smooth random fields in two dimensions

Authors: Pravabati Chingangbam, Priya Goyal, K. P. Yogendran, Stephen Appleby

Abstract: We revisit the geometrical meaning of statistical isotropy that is manifest in excursion sets of smooth random fields in two dimensions. Using the contour Minkowski tensor, $\W_1$, as our basic tool we first examine geometrical properties of single structures. For simple closed curves in two dimensions we show that $\W_1$ is proportional to the identity matrix if the curve has $m$-fold symmetry, w… ▽ More We revisit the geometrical meaning of statistical isotropy that is manifest in excursion sets of smooth random fields in two dimensions. Using the contour Minkowski tensor, $\W_1$, as our basic tool we first examine geometrical properties of single structures. For simple closed curves in two dimensions we show that $\W_1$ is proportional to the identity matrix if the curve has $m$-fold symmetry, with $m\ge 3$. Then we elaborate on how $\W_1$ maps any arbitrary shaped simple closed curve to an ellipse that is unique up to translations of its centroid. We also carry out a comparison of the shape parameters, $α$ and $β$, defined using $\W_1$, with the filamentarity parameter defined using two scalar Minkowski functionals - area and contour length. We show that they contain complementary shape information, with $\W_1$ containing additional information of orientation of structures. Next, we apply our method to boundaries of excursion sets of random fields and examine what statistical isotropy means for the geometry of the excursion sets. Focusing on Gaussian isotropic fields, and using a semi-numerical approach we quantify the effect of finite sampling of the field on the geometry of the excursion sets. In doing so we obtain an analytic expression for $α$ which takes into account the effect of finite sampling. Finally we derive an analytic expression for the ensemble expectation of $\W_1$ for Gaussian anisotropic random fields. Our results provide insights that are useful for designing tests of statistical isotropy using cosmological data. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 14 pages

Well-Tempered Minkowski Solutions in Teleparallel Horndeski Theory

Authors: Reginald Christian Bernardo, Jackson Levi Said, Maria Caruana, Stephen Appleby

Abstract: Well-tempering stands among the few classical methods of screening vacuum energy to deliver a late-time, low energy vacuum state. We build on the class of Horndeski models that admit a Minkowski vacuum state despite the presence of an arbitrarily large vacuum energy to obtain a much larger family of models in teleparallel Horndeski theory. We set up the routine for obtaining these models and prese… ▽ More Well-tempering stands among the few classical methods of screening vacuum energy to deliver a late-time, low energy vacuum state. We build on the class of Horndeski models that admit a Minkowski vacuum state despite the presence of an arbitrarily large vacuum energy to obtain a much larger family of models in teleparallel Horndeski theory. We set up the routine for obtaining these models and present a variety of cases, all of which are able to screen a natural particle physics scale vacuum energy using degeneracy in the field equations. We establish that well-tempering is the unique method of utilizing degeneracy in Horndeski scalar-tensor gravity -- and its teleparallel generalisation -- that can accommodate self-tuned flat Minkowski solutions, when the explicit scalar field dependence in the action is minimal (a tadpole and a conformal coupling to the Ricci scalar). Finally, we study the dynamics of the well-tempered teleparallel Galileon. We generate its phase portraits and assess the attractor nature of the Minkowski vacuum under linear perturbations and through a phase transition of vacuum energy. △ Less

Submitted 4 November, 2021; v1 submitted 5 August, 2021; originally announced August 2021.

Comments: 34 pages, 5 figures, v2: minor revisions, template changed, to appear in CQG

Journal ref: 2022 Class. Quantum Grav. 39 015013

Well-tempered teleparallel Horndeski cosmology: a teleparallel variation to the cosmological constant problem

Authors: Reginald Christian Bernardo, Jackson Levi Said, Maria Caruana, Stephen Appleby

Abstract: Well-tempering is a promising classical method of dynamically screening an arbitrarily large vacuum energy and generating a late-time, low energy de Sitter vacuum state. In this paper, we study for the first time self-tuning in teleparallel gravity and obtain well-tempered cosmological models in the teleparallel gravity analogue of Horndeski theory. This broadens the scope of well-tempered cosmolo… ▽ More Well-tempering is a promising classical method of dynamically screening an arbitrarily large vacuum energy and generating a late-time, low energy de Sitter vacuum state. In this paper, we study for the first time self-tuning in teleparallel gravity and obtain well-tempered cosmological models in the teleparallel gravity analogue of Horndeski theory. This broadens the scope of well-tempered cosmology and teases the potentially far richer cosmological dynamics that could be anchored on teleperallel gravity. We expand the well-tempered recipe to its most general form so far and use it to search for the first well-tempered cosmologies in teleparallel gravity. We also study the cosmological dynamics in a well-tempered model and demonstrate the dynamical stability of the vacuum, the compatibility with a matter era, and the stability of the vacuum through a phase transition. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 35 pages, 5 figures, 1 table, v3: some typos fixed

Journal ref: JCAP 10 (2021) 078

The Low Redshift Circumgalactic Medium in Simba

Authors: Sarah Appleby, Romeel Davé, Daniele Sorini, Kate Storey-Fisher, Britton Smith

Abstract: We examine the properties of the low-redshift circumgalactic medium (CGM) around star-forming and quenched galaxies in the Simba cosmological hydrodynamic simulations, focusing on comparing HI and metal line absorption to observations from the COS-Halos and COS-Dwarfs surveys. Halo baryon fractions are generally $\lesssim 50\%$ of the cosmic fraction due to stellar feedback at low masses, and jet-… ▽ More We examine the properties of the low-redshift circumgalactic medium (CGM) around star-forming and quenched galaxies in the Simba cosmological hydrodynamic simulations, focusing on comparing HI and metal line absorption to observations from the COS-Halos and COS-Dwarfs surveys. Halo baryon fractions are generally $\lesssim 50\%$ of the cosmic fraction due to stellar feedback at low masses, and jet-mode AGN feedback at high masses. Baryons and metals in the CGM of quenched galaxies are $\gtrsim 90\%$ hot gas, while the CGM of star-forming galaxies is more multi-phase. Hot CGM gas has low metallicity, while warm and cool CGM gas have metallicity close to that of galactic gas. Equivalent widths, covering fractions and total path absorption of HI and selected metal lines (MgII, SiIII, CIV and OVI) around a matched sample of Simba star-forming galaxies are mostly consistent with COS-Halos and COS-Dwarfs observations to $\lesssim 0.4$~dex, depending on ion and assumed ionising background. Around matched quenched galaxies, absorption in all ions is lower, with HI absorption significantly under-predicted. Metal-line absorption is sensitive to choice of photo-ionising background; assuming recent backgrounds, Simba matches OVI but under-predicts low ions, while an older background matches low ions but under-predicts OVI. Simba reproduces the observed dichotomy of OVI absorption around star forming and quenched galaxies. CGM metals primarily come from stellar feedback, while jet-mode AGN feedback reduces absorption particularly for lower ions. △ Less

Submitted 6 September, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Comments: 24 pages, 12 figures, accepted for publication in MNRAS

Cosmological Parameter Estimation from the Two-Dimensional Genus Topology -- Measuring the Expansion History using the Genus Amplitude as a Standard Ruler

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim, Motonari Tonegawa

Abstract: We measure the genus of the galaxy distribution in two-dimensional slices of the SDSS-III BOSS catalog to constrain the cosmological parameters governing the expansion history of the Universe. The BOSS catalogs are divided into twelve concentric shells over the redshift range $0.25 < z < 0.6$ and we repeatedly measure the genus from the two-dimensional galaxy density fields, each time varying the… ▽ More We measure the genus of the galaxy distribution in two-dimensional slices of the SDSS-III BOSS catalog to constrain the cosmological parameters governing the expansion history of the Universe. The BOSS catalogs are divided into twelve concentric shells over the redshift range $0.25 < z < 0.6$ and we repeatedly measure the genus from the two-dimensional galaxy density fields, each time varying the cosmological parameters used to infer the distance-redshift relation to the shells. We also indirectly reconstruct the two-dimensional genus amplitude using the three-dimensional genus measured from SDSS Main Galaxy Sample with galaxies at low redshift $z < 0.12$. We combine the low- and high-redshift measurements, finding the cosmological model which minimizes the redshift evolution of the genus amplitude, using the fact that this quantity should be conserved. Being a distance measure, the test is sensitive to the matter density parameter ($Ω_{\rm m}$) and equation of state of dark energy ($w_{\rm de}$). We find a constraint of $w_{\rm de} = -1.05^{+0.13}_{-0.12}$, $Ω_{\rm m} = 0.303 \pm 0.036$ after combining the high- and low-redshift measurements and combining with Planck CMB data. Higher redshift data and combining data sets at low redshift will allow for stronger constraints. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: Accepted for publication in the Astrophysical Journal. 19 pages, 13 figures

Prospects of constraining reionization model parameters using Minkowski tensors and Betti numbers

Authors: Akanksha Kapahtia, Pravabati Chingangbam, Raghunath Ghara, Stephen Appleby, Tirthankar Roy Choudhury

Abstract: We explore the possibility of constraining model parameters of the Epoch of Reionization (EoR) from 21cm brightness temperature maps, using a combination of morphological descriptors constructed from the eigenvalues of the Contour Minkowski Tensor (CMT), Betti numbers (count of connected regions $n_{con}$ and holes $n_{hole}$) and the area of structures in the excursion set of the field. We use a… ▽ More We explore the possibility of constraining model parameters of the Epoch of Reionization (EoR) from 21cm brightness temperature maps, using a combination of morphological descriptors constructed from the eigenvalues of the Contour Minkowski Tensor (CMT), Betti numbers (count of connected regions $n_{con}$ and holes $n_{hole}$) and the area of structures in the excursion set of the field. We use a three parameter model of EoR simulated using 21\textrm{cmFAST}, namely the ionizing efficiency of sources $ζ$, the minimum virial temperature $T_{vir}$ required for collapse into a halo and the maximum radius for ionizing radiation described by $R_{mfp}$. We performed a Bayesian analysis to recover model parameters for a mock 21cm image from SKA phase I at a redshift of $z=7.4$ corresponding to a mean neutral hydrogen fraction of $\mathrm{\bar x}_{HI} \simeq 0.5$. We find that in the absence of noise the average size of structures in the field with $x_{HI} \lesssim 0.5$ is smaller than regions with $x_{HI} \gtrsim 0.5$ and the structures are equally isotropic when $\mathrm{\bar x}_{HI}=0.5$ . We also find that in order to recover the input model to within $1-σ$ accuracy for a mock noisy image at a single frequency channel of $1~\mathrm{MHz}$, for an observation time $t_{obs}<2000~\mathrm{hrs}$, the noisy $δT_b$ map needs to be smoothed at a scale $R_s>9.5~\mathrm{Mpc}$. Finally we show that the systematic behaviour of the statistic as ionization progresses, enables us to obtain stringent constraints on $\mathrm{\bar x}_{HI}$ (with a coefficient of variation $\sim 0.05$ as compared to $\sim 0.1-0.2$ for model parameter constraints), thereby making these descriptors a promising statistic for constraining EoR model parameters and the ionization history of the universe. △ Less

Submitted 22 April, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in JCAP

An Expansion of Well Tempered Gravity

Authors: Eric V. Linder, Stephen Appleby

Abstract: When faced with two nigh intractable problems in cosmology -- how to remove the original cosmological constant problem and how to parametrize modified gravity to explain current cosmic acceleration -- we can make progress by counterposing them. The well tempered solution to the cosmological constant through degenerate scalar field dynamics also relates disparate Horndeski gravity terms, making the… ▽ More When faced with two nigh intractable problems in cosmology -- how to remove the original cosmological constant problem and how to parametrize modified gravity to explain current cosmic acceleration -- we can make progress by counterposing them. The well tempered solution to the cosmological constant through degenerate scalar field dynamics also relates disparate Horndeski gravity terms, making them contrapuntal. We derive the connection between the kinetic term $K$ and braiding term $G_3$ for shift symmetric theories (including the running Planck mass $G_4$), extending previous work on monomial or binomial dependence to polynomials of arbitrary finite degree. We also exhibit an example for an infinite series expansion. This contrapuntal condition greatly reduces the number of parameters needed to test modified gravity against cosmological observations, for these "golden" theories of gravity. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: 7 pages

Journal ref: JCAP 2103, 074 (2021)

The Well-Tempered Cosmological Constant: Fugue in B$^\flat$

Authors: Stephen Appleby, Eric V. Linder

Abstract: Zero point fluctuations of quantum fields should generate a large cosmological constant energy density in any spacetime. How then can we have anything other than de Sitter space without fine tuning? Well tempering -- dynamical cancellation of the cosmological constant using degeneracy within the field equations -- can replace a large cosmological constant with a much lower energy state. Here we gi… ▽ More Zero point fluctuations of quantum fields should generate a large cosmological constant energy density in any spacetime. How then can we have anything other than de Sitter space without fine tuning? Well tempering -- dynamical cancellation of the cosmological constant using degeneracy within the field equations -- can replace a large cosmological constant with a much lower energy state. Here we give an explicit mechanism to obtain a Minkowski solution, replacing the cosmological constant with zero, and testing its attractor nature and persistence through a vacuum phase transition. We derive the general conditions that Horndeski scalar-tensor gravity must possess, and evolve in a fugue of functions, to deliver nothing and make the universe be flat. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 15 pages, 3 figures

Report number: APCTP Pre2020-023

Journal ref: JCAP 2012, 037 (2020)

The Well-Tempered Cosmological Constant: The Horndeski Variations

Authors: Stephen Appleby, Eric V. Linder

Abstract: Well tempering is one of the few classical field theory methods for solving the original cosmological constant problem, dynamically canceling a large (possibly Planck scale) vacuum energy and leaving the matter component intact, while providing a viable cosmology with late time cosmic acceleration and an end de Sitter state. We present the general constraints that variations of Horndeski gravity m… ▽ More Well tempering is one of the few classical field theory methods for solving the original cosmological constant problem, dynamically canceling a large (possibly Planck scale) vacuum energy and leaving the matter component intact, while providing a viable cosmology with late time cosmic acceleration and an end de Sitter state. We present the general constraints that variations of Horndeski gravity models with different combinations of terms must satisfy to admit an exact de Sitter spacetime that does not respond to an arbitrarily large cosmological constant. We explicitly derive several specific scalar-tensor models that well temper and can deliver a standard cosmic history including current cosmic acceleration. Stability criteria, attractor behavior of the de Sitter state, and the response of the models to pressureless matter are considered. The well tempered conditions can be used to focus on particular models of modified gravity that have special interest -- not only removing the original cosmological constant problem but providing relations between the free Horndeski functions and reducing them to a couple of parameters, suitable for testing gravity and cosmological data analysis. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 25 pages, 3 figures

Report number: APCTP Pre2020-022

Journal ref: JCAP 2012, 036 (2020)

Identification of MHD compressible modes in interstellar plasma with synchrotron emission polarization

Authors: Alexey Chepurnov, Reinaldo Santos de Lima, Sarah Appleby

Abstract: We provide a procedure for identification of dominating compressible and Alfvenic MHD modes or isotropic turbulence in synchrotron emission polarization maps of Galactic objects. The results for the region of North Galactic Pole, Orion molecular cloud complex and the star-forming complex Cygnus X are presented. We provide a procedure for identification of dominating compressible and Alfvenic MHD modes or isotropic turbulence in synchrotron emission polarization maps of Galactic objects. The results for the region of North Galactic Pole, Orion molecular cloud complex and the star-forming complex Cygnus X are presented. △ Less

Submitted 17 June, 2022; v1 submitted 7 April, 2020; originally announced April 2020.

Cosmological Parameter Estimation from the Two-Dimensional Genus Topology -- Measuring the Shape of the Matter Power Spectrum

Authors: Stephen A. Appleby, Changbom Park, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim

Abstract: We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian… ▽ More We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian expectation values, exploiting the fact that this quantity is relatively insensitive to non-linear gravitational collapse. The genus amplitude provides a measure of the shape of the linear matter power spectrum, and is principally sensitive to $Ω_{\rm c}h^{2}$ and scalar spectral index $n_{\rm s}$. A strong negative degeneracy between $Ω_{\rm c}h^{2}$ and $n_{\rm s}$ is observed, as both can increase small scale power by shifting the peak and tilting the power spectrum respectively. We place a constraint on the particular combination $n_{\rm s}^{3/2} Ω_{\rm c}h^{2}$ -- we find $n_{\rm s}^{3/2} Ω_{\rm c}h^{2} = 0.1121 \pm 0.0043$ after combining the LOWZ and CMASS data sets, assuming a flat $Λ$CDM cosmology. This result is practically insensitive to reasonable variations of the power spectrum amplitude and linear galaxy bias. Our results are consistent with the Planck best fit $n_{\rm s}^{3/2}Ω_{\rm c}h^{2} = 0.1139 \pm 0.0009$. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: ApJ submitted, comments welcome

Morphology of CMB fields -- effect of weak gravitational lensing

Authors: Priya Goyal, Pravabati Chingangbam, Stephen Appleby

Abstract: We study the morphology of the cosmic microwave background temperature and polarization fields using the shape and alignment parameters, $β$ and $α$, that are constructed from the contour Minkowski tensor. The primary goal of our paper is to understand the effect of weak gravitational lensing on the morphology of the CMB fields. In order to isolate different physical effects that can be potentiall… ▽ More We study the morphology of the cosmic microwave background temperature and polarization fields using the shape and alignment parameters, $β$ and $α$, that are constructed from the contour Minkowski tensor. The primary goal of our paper is to understand the effect of weak gravitational lensing on the morphology of the CMB fields. In order to isolate different physical effects that can be potentially confused with the effect of lensing, we first study the effect of varying the cosmology on $α$ and $β$, and show that they are relatively insensitive to variation of cosmological parameters. Next we analyze the signatures of hemispherical anisotropy, and show that information of such anisotropy in $α$ gets washed out at small angular scales and become pronounced only at large angular scales. For $β$ we find characteristic distortions which vary with the field threshold. We then study the effect of weak gravitational lensing using simulations of lensed temperature and $E$ and $B$ modes. We quantify the distortion induced in the fields across different angular scales. We find that lensing makes structures of all fields increasingly more anisotropic as we probe down to smaller scales. We find distinct behaviour of morphological distortions as a function of threshold for the different fields. The effect is small for temperature and $E$ mode, while it is significantly large for $B$ mode. Further, we find that lensing does not induce statistical anisotropy, as expected from the isotropic distribution of large scale structure of matter. We expect that the results obtained in this work will provide insights on the reconstruction of the lensing potential. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: 26+1 pages

The impact of quenching on galaxy profiles in the Simba simulation

Authors: Sarah Appleby, Romeel Davé, Katarina Kraljic, Daniel Anglés-Alcázar, Desika Narayanan

Abstract: We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to num… ▽ More We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to numerical heating. We compare Simba $z=0$ sSFR radial profiles for star forming and green valley galaxies to observations. Simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in green valley galaxies and massive star-forming systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, Simba yields green valley galaxies with strongly dropping sSFR profiles beyond $\gtrsim R_{\rm half}$, regardless of AGN feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. Simba's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At $z=2$ Simba does not show central depressions in massive star-forming galaxies, suggesting Simba's X-ray feedback should be more active at high-$z$. High resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in $z=0$ green valley galaxies represents a unique success of Simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms. △ Less

Submitted 20 May, 2020; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: 21 pages, 11 figures, accepted by MNRAS

Ensemble Average of Three-Dimensional Minkowski Tensors of a Gaussian Random Field in Redshift Space

Authors: Stephen Appleby, Joby P. K., Pravabati Chingangbam, Changbom Park

Abstract: We present the ensemble expectation values for the translation invariant, rank-2 Minkowski tensors in three-dimensions, for a linearly redshift space distorted Gaussian random field. The Minkowski tensors $W^{0,2}_{1}$, $W^{0,2}_{2}$ are sensitive to global anisotropic signals present within a field, and by extracting these statistics from the low redshift matter density one can place constraints… ▽ More We present the ensemble expectation values for the translation invariant, rank-2 Minkowski tensors in three-dimensions, for a linearly redshift space distorted Gaussian random field. The Minkowski tensors $W^{0,2}_{1}$, $W^{0,2}_{2}$ are sensitive to global anisotropic signals present within a field, and by extracting these statistics from the low redshift matter density one can place constraints on the redshift space distortion parameter $β= f/b$. We begin by reviewing the calculation of the ensemble expectation values $\langle W^{0,2}_{1} \rangle$, $\langle W^{0,2}_{2} \rangle $ for isotropic, Gaussian random fields, then consider how these results are modified by the presence of a linearly anisotropic signal. Under the assumption that all fields remain Gaussian, we calculate the anisotropic correction due to redshift space distortion in a coordinate system aligned with the line of sight, finding inequality between the diagonal elements of $\langle W^{0,2}_{1} \rangle $, $\langle W^{0,2}_{2} \rangle $. The ratio of diagonal elements of these matrices provides a set of statistics that are sensitive only to the redshift space distortion parameter $β$. We estimate the Fisher information that can be extracted from the Minkowski tensors, and find $W^{0,2}_{1}$ is more sensitive to $β$ than $W^{0,2}_{2}$, and a measurement of $W^{0,2}_{1}$ accurate to $\sim 1\%$ can yield a $\sim 4\%$ constraint on $β$. Finally, we discuss the difference between using the matrix elements of the Minkowski tensors directly against measuring the eigenvalues. For the purposes of cosmological parameter estimation we advocate the use of the matrix elements, to avoid spurious anisotropic signals that can be generated by the eigenvalue decomposition. △ Less

Submitted 7 August, 2019; originally announced August 2019.

Morphology of 21cm brightness temperature during the Epoch of Reioinization using Contour Minkowski Tensor

Authors: Akanksha Kapahtia, Pravabati Chingangbam, Stephen Appleby

Abstract: We use morphological descriptors, Betti numbers and Contour Minkowski Tensor (CMT) on 21cm brightness temperature excursion sets, to study the ionization and heating history of the intergalactic medium (IGM) during and before the Epoch of Reionization (EoR). The ratio of eigenvalues of the CMT denoted by $β$, gives shape information while it's trace gives the contour length of holes and connected… ▽ More We use morphological descriptors, Betti numbers and Contour Minkowski Tensor (CMT) on 21cm brightness temperature excursion sets, to study the ionization and heating history of the intergalactic medium (IGM) during and before the Epoch of Reionization (EoR). The ratio of eigenvalues of the CMT denoted by $β$, gives shape information while it's trace gives the contour length of holes and connected regions. We simulate the matter density, neutral hydrogen fraction, spin temperature and brightness temperature field using the publicly available code 21cmFAST in a redshift range of $z=20.22$ to $z=6$. We study the redshift evolution of three quantities - the Betti number counts $N_{con,hole}$, the characteristic size $r^{ch}_{con,hole}$ and shape anisotropy parameter $β^{ch}_{con,hole}$ of connected regions and holes for these fields and investigate the different physical origins of their evolution. We make a qualitative comparison of different models of heating and ionization during the EoR. We obtain different regimes of morphological evolution of brightness temperature, depending upon how the shapes and sizes of connected regions and holes change with redshift for different astrophysical settings affecting the ionization and heating history of the IGM during and before the EoR. We find that the morphology of the brightness temperature field traces the morphology of ionized regions below a certain redshift value depending upon the model, where $Δr^{ch}_{hole}<10 \%$ and $Δβ^{ch}_{hole}<1 \%$ relative to the $x_{HI}$ field. This difference decreases with redshift. Therefore, the ionization history of the IGM can be reconstructed using the morphological description of $δT_b$ in real space. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 31 Pages, 15 figures

The VANDELS survey: the star-formation histories of massive quiescent galaxies at 1.0 < z < 1.3

Authors: A. C. Carnall, R. J. McLure, J. S. Dunlop, F. Cullen, D. J. McLeod, V. Wild, B. D. Johnson, S. Appleby, R. Davé, R. Amorin, M. Bolzonella, M. Castellano, A. Cimatti, O. Cucciati, A. Gargiulo, B. Garilli, F. Marchi, L. Pentericci, L. Pozzetti, C. Schreiber, M. Talia, G. Zamorani

Abstract: We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints… ▽ More We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity. △ Less

Submitted 24 October, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: 25 pages, 13 figures, MNRAS accepted

Simba: Cosmological Simulations with Black Hole Growth and Feedback

Authors: Romeel Davé, Daniel Anglés-Alcázar, Desika Narayanan, Qi Li, Mika H. Rafieferantsoa, Sarah Appleby

Abstract: We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Anglés-Alcázar et al. (2017) from cold gas and Bondi accretion from… ▽ More We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Anglés-Alcázar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies. △ Less

Submitted 2 April, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: 25 pages, MNRAS accepted, some new figures

Identification of plasma modes in Galactic turbulence with synchrotron polarization

Authors: Heshou Zhang, Alexey Chepurnov, Huirong Yan, Kirit Makwana, Reinaldo Santos-Lima, Sarah Appleby

Abstract: Magneto-hydrodynamic (MHD) turbulence is ubiquitous and a fundamental ingredient underlying many astrophysical phenomena. The multiphase nature of interstellar medium and diversity of driving mechanisms give rise to spatial variation of turbulence properties, particularly their plasma properties. There has been no observational diagnosis of the plasma modes so far beyond the solar system. Here we… ▽ More Magneto-hydrodynamic (MHD) turbulence is ubiquitous and a fundamental ingredient underlying many astrophysical phenomena. The multiphase nature of interstellar medium and diversity of driving mechanisms give rise to spatial variation of turbulence properties, particularly their plasma properties. There has been no observational diagnosis of the plasma modes so far beyond the solar system. Here we report the identification of different plasma modes in various Galactic environments, including active star forming zones and supernova remnants, based on our synchrotron polarization analysis. The observed high degree of consistency between the $γ-$ray excess in Cygnus cocoon and the location of magnetosonic modes provides strong observational evidence for the long-advocated theory that magnetosonic modes dominate the cosmic ray (CR) scattering and acceleration. Our results open up a new avenue for the study of interstellar turbulence and demonstrate the indispensability to account for their plasma property in all the relevant processes including CR transport and star formation. △ Less

Submitted 21 February, 2021; v1 submitted 3 August, 2018; originally announced August 2018.

Comments: 17 pages (including Methods), 5 figures, Nature Astronomy published

Journal ref: Nature Astronomy (2020), Volume 4, p. 1001-1008

Minkowski Tensors in Three Dimensions - Probing the Anisotropy Generated by Redshift Space Distortion

Authors: Stephen Appleby, Pravabati Chingangbam, Changbom Park, K. P. Yogendran, P. K. Joby

Abstract: We apply the Minkowski tensor statistics to three dimensional Gaussian random fields. Minkowski tensors contain information regarding the orientation and shape of excursion sets, that is not present in the scalar Minkowski functionals. They can be used to quantify globally preferred directions, and additionally provide information on the mean shape of subsets of a field. This makes them ideal stat… ▽ More We apply the Minkowski tensor statistics to three dimensional Gaussian random fields. Minkowski tensors contain information regarding the orientation and shape of excursion sets, that is not present in the scalar Minkowski functionals. They can be used to quantify globally preferred directions, and additionally provide information on the mean shape of subsets of a field. This makes them ideal statistics to measure the anisotropic signal generated by redshift space distortion in the low redshift matter density field. We review the definition of the Minkowski tensor statistics in three dimensions, focusing on two coordinate invariant quantities $W^{0,2}_{1}$ and $W^{0,2}_{2}$. We calculate the ensemble average of these $3 \times 3$ matrices for an isotropic Gaussian random field, finding that they are proportional to products of the identity matrix and a corresponding scalar Minkowski functional. We show how to numerically reconstruct $W^{0,2}_{1}$ and $W^{0,2}_{2}$ from discretely sampled fields and apply our algorithm to isotropic Gaussian fields generated from a linear $Λ$CDM matter power spectrum. We then introduce anisotropy by applying a linear redshift space distortion operator to the matter density field, and find that both $W^{0,2}_{1}$ and $W^{0,2}_{2}$ exhibit a distinct signal characterised by inequality between their diagonal components. We discuss the physical origin of this signal and how it can be used to constrain the redshift space distortion parameter $Υ\equiv f/b$. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 21 pages

The Well-Tempered Cosmological Constant

Authors: Stephen Appleby, Eric V. Linder

Abstract: Self tuning is one of the few methods for dynamically cancelling a large cosmological constant and yet giving an accelerating universe. Its drawback is that it tends to screen all sources of energy density, including matter. We develop a model that tempers the self tuning so the dynamical scalar field still cancels an arbitrary cosmological constant, including the vacuum energy through any high en… ▽ More Self tuning is one of the few methods for dynamically cancelling a large cosmological constant and yet giving an accelerating universe. Its drawback is that it tends to screen all sources of energy density, including matter. We develop a model that tempers the self tuning so the dynamical scalar field still cancels an arbitrary cosmological constant, including the vacuum energy through any high energy phase transitions, without affecting the matter fields. The scalar-tensor gravitational action is simple, related to cubic Horndeski gravity, with a nonlinear derivative interaction plus a tadpole term. Applying shift symmetry and using the property of degeneracy of the field equations we find families of functions that admit de Sitter solutions with expansion rates that are independent of the magnitude of the cosmological constant and preserve radiation and matter dominated phases. That is, the method can deliver a standard cosmic history including current acceleration, despite the presence of a Planck scale cosmological constant. △ Less

Submitted 30 April, 2018; originally announced May 2018.

Comments: 19 pages, 4 figures

Journal ref: JCAP 1807, 034 (2018)

Cosmological Parameter Estimation Using the Genus Amplitude - Application to Mock Galaxy Catalogs

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering… ▽ More We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering both its redshift evolution and magnitude. The constancy of the genus amplitude with redshift can be used as a standard population, from which we derive constraints on the equation of state of dark energy $w_{\rm de}$ - by measuring $A$ at $z \sim 0.1$ and $z \sim 1$, we can place an order $Δw_{\rm de} \sim {\cal O}(15\%)$ constraint on $w_{\rm de}$. By comparing $A$ to its Gaussian expectation value we can potentially derive an additional stringent constraint on the matter density $ΔΩ_{\rm mat} \sim 0.01$. We discuss the primary sources of contamination associated with the two measurements - redshift space distortion and shot noise. With accurate knowledge of galaxy bias, we can successfully remove the effect of redshift space distortion, and the combined effect of shot noise and non-linear gravitational evolution is suppressed by smoothing over suitably large scales $R_{\rm G} \ge 15 {\rm Mpc}/h$. Without knowledge of the bias, we discuss how joint measurements of the two and three dimensional genus can be used to constrain the growth factor $β= f/b$. The method can be applied optimally to redshift slices of a galaxy distribution generated using the drop-off technique. △ Less

Submitted 19 January, 2018; originally announced January 2018.

Comments: Published 01/2018. 15 pages, 10 figures

Journal ref: ApJ 853 17 (2018)

A novel probe of bubble size statistics and time scales of the epoch of reionization using the contour Minkowski Tensor

Authors: Akanksha Kapahtia, Pravabati Chingangbam, Stephen Appleby, Changbom Park

Abstract: We employ the rank-2 {\em contour} Minkowski Tensor in two dimensions to probe length and time scales of ionized bubbles during the epoch of reionization. We demonstrate that the eigenvalues of this tensor provide excellent probes of the distribution of the sizes of ionized bubbles, and from it the characteristic bubble sizes, at different redshifts. We show that ionized bubbles are not circular,… ▽ More We employ the rank-2 {\em contour} Minkowski Tensor in two dimensions to probe length and time scales of ionized bubbles during the epoch of reionization. We demonstrate that the eigenvalues of this tensor provide excellent probes of the distribution of the sizes of ionized bubbles, and from it the characteristic bubble sizes, at different redshifts. We show that ionized bubbles are not circular, and hence not spherical in three dimensions, as is often assumed for simplified analytic arguments. We quantify their shape anisotropy by using the ratio of the two eigenvalues. The shape parameter provides the characteristic time epochs when bubble mergers begin and end. Our method will be very useful to reconstruct the reionization history using data of the brightness temperature field. △ Less

Submitted 19 November, 2020; v1 submitted 26 December, 2017; originally announced December 2017.

Comments: 22 pages, matches published version, factor of 2 typo in definition of CMT corrected

Journal ref: JCAP Issue 10, 011 (2018)

Minkowski Tensors in Two Dimensions - Probing the Morphology and Isotropy of the Matter and Galaxy Density Fields

Authors: Stephen Appleby, Pravabati Chingangbam, Changbom Park, Sungwook E. Hong, Juhan Kim, Vidhya Ganesan

Abstract: We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimension… ▽ More We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($Λ_{1},Λ_{2}$) and for individual connected regions and holes within the set ($λ_{1},λ_{2}$). The ratio of eigenvalues $Λ_{2}/Λ_{1}$ informs us of the presence of global anisotropies in the data, and $\langle λ_{2}/λ_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $Λ_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $Λ_{1,2}$, suggesting that they can be used to probe the large-scale velocity field. △ Less

Submitted 22 May, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: 17 pages, accepted for publication in APJ

Tensor Minkowski Functionals for random fields on the sphere

Authors: Pravabati Chingangbam, K P Yogendran, Joby P. K., Vidhya Ganesan, Stephen Appleby, Changbom Park

Abstract: We generalize the translation invariant tensor-valued Minkowski Functionals which are defined on two-dimensional flat space to the unit sphere. We apply them to level sets of random fields. The contours enclosing boundaries of level sets of random fields give a spatial distribution of random smooth closed curves. We obtain analytic expressions for the ensemble expectation values for the matrix ele… ▽ More We generalize the translation invariant tensor-valued Minkowski Functionals which are defined on two-dimensional flat space to the unit sphere. We apply them to level sets of random fields. The contours enclosing boundaries of level sets of random fields give a spatial distribution of random smooth closed curves. We obtain analytic expressions for the ensemble expectation values for the matrix elements of the tensor-valued Minkowski Functionals for isotropic Gaussian and Rayleigh fields. We elucidate the way in which the elements of the tensor Minkowski Functionals encode information about the nature and statistical isotropy (or departure from isotropy) of the field. We then implement our method to compute the tensor-valued Minkowski Functionals numerically and demonstrate how they encode statistical anisotropy and departure from Gaussianity by applying the method to maps of the Galactic foreground emissions from the PLANCK data. △ Less

Submitted 29 November, 2017; v1 submitted 14 July, 2017; originally announced July 2017.

Comments: 1+23 pages, 5 figures, Significantly expanded from version 1. To appear in JCAP

Topology of Large-Scale Structures of Galaxies in Two Dimensions - Systematic Effects

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic a… ▽ More We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic and numerical artifacts such as shot noise, galaxy bias, and finite pixel effects. We model the systematics using a Hermite polynomial expansion and perform a comprehensive analysis of known effects on the two-dimensional genus, with a view toward using the statistic for cosmological parameter estimation. We find that the finite pixel effect is dominated by an amplitude drop and can be made less than $1\%$ by adopting pixels smaller than $1/3$ of the angular smoothing length. Nonlinear gravitational evolution introduces time-dependent coefficients of the zeroth, first, and second Hermite polynomials, but the genus amplitude changes by less than $1\%$ between $z=1$ and $z=0$ for smoothing scales $R_{\rm G} > 9 {\rm Mpc/h}$. Non-zero terms are measured up to third order in the Hermite polynomial expansion when studying RSD. Differences in shapes of the genus curves in real and redshift space are small when we adopt thick redshift shells, but the amplitude change remains a significant $\sim {\cal O}(10\%)$ effect. The combined effects of galaxy biasing and shot noise produce systematic effects up to the second Hermite polynomial. It is shown that, when sampling, the use of galaxy mass cuts significantly reduces the effect of shot noise relative to random sampling. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 836 45 (2017)

Cosmology and Fundamental Physics with the Euclid Satellite

Authors: Luca Amendola, Stephen Appleby, Anastasios Avgoustidis, David Bacon, Tessa Baker, Marco Baldi, Nicola Bartolo, Alain Blanchard, Camille Bonvin, Stefano Borgani, Enzo Branchini, Clare Burrage, Stefano Camera, Carmelita Carbone, Luciano Casarini, Mark Cropper, Claudia de Rham, Joerg P. Dietrich, Cinzia Di Porto, Ruth Durrer, Anne Ealet, Pedro G. Ferreira, Fabio Finelli, Juan Garcia-Bellido, Tommaso Giannantonio , et al. (47 additional authors not shown)

Abstract: Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of c… ▽ More Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: This article provides an update of arXiv:1206.1225, with different authors. Forecasts are not updated in this version

Journal ref: Living Reviews in Relativity 21 (2018)

Direct search for features in the primordial bispectrum

Authors: Stephen Appleby, Jinn-Ouk Gong, Dhiraj Kumar Hazra, Arman Shafieloo, Spyros Sypsas

Abstract: We study features in the bispectrum of the primordial curvature perturbation correlated with the reconstructed primordial power spectrum from the observed cosmic microwave background temperature data. We first show how the bispectrum can be completely specified in terms of the power spectrum and its first two derivatives, valid for any configuration of interest. Then using a model-independent reco… ▽ More We study features in the bispectrum of the primordial curvature perturbation correlated with the reconstructed primordial power spectrum from the observed cosmic microwave background temperature data. We first show how the bispectrum can be completely specified in terms of the power spectrum and its first two derivatives, valid for any configuration of interest. Then using a model-independent reconstruction of the primordial power spectrum from the Planck angular power spectrum of temperature anisotropies, we compute the bispectrum in different triangular configurations. We find that in the squeezed limit at k ~ 0.06/Mpc and k ~ 0.014/Mpc there are marginal 2sigma deviations from the standard featureless bispectrum, which meanwhile is consistent with the reconstructed bispectrum in the equilateral configuration. △ Less

Submitted 15 July, 2016; v1 submitted 30 December, 2015; originally announced December 2015.

Comments: (v1) 6 pages, 4 figures; (v2) improved results but conclusions not changed, to appear in Physics Letters B

Report number: APCTP-Pre2015-027

Journal ref: Phys. Lett. B 760 (2016) 297-301

Self Tuning Scalar Fields in Spherically Symmetric Spacetimes

Authors: Stephen Appleby

Abstract: We search for self tuning solutions to the Einstein-scalar field equations for the simplest class of `Fab-Four' models with constant potentials. We first review the conditions under which self tuning occurs in a cosmological spacetime, and by introducing a small modification to the original theory - introducing the second and third Galileon terms - show how one can obtain de Sitter states where th… ▽ More We search for self tuning solutions to the Einstein-scalar field equations for the simplest class of `Fab-Four' models with constant potentials. We first review the conditions under which self tuning occurs in a cosmological spacetime, and by introducing a small modification to the original theory - introducing the second and third Galileon terms - show how one can obtain de Sitter states where the expansion rate is independent of the vacuum energy. We then consider whether the same self tuning mechanism can persist in a spherically symmetric inhomogeneous spacetime. We show that there are no asymptotically flat solutions to the field equations in which the vacuum energy is screened, other than the trivial one (Minkowski space). We then consider the possibility of constructing Schwarzschild de Sitter spacetimes for the modified Fab Four plus Galileon theory. We argue that the only model that can successfully screen the vacuum energy in both an FLRW and Schwarzschild de Sitter spacetime is one containing `John' $\sim G^μ{}_ν \partial_μφ\partial^νφ$ and a canonical kinetic term $\sim \partial_αφ\partial^αφ$. This behaviour was first observed in (Babichev&Charmousis,2013). The screening mechanism, which requires redundancy of the scalar field equation in the `vacuum', fails for the `Paul' term in an inhomogeneous spacetime. △ Less

Submitted 18 March, 2015; originally announced March 2015.

Probing bulk flow with nearby SNe Ia data

Authors: Stephen Appleby, Arman Shafieloo, Andrew Johnson

Abstract: We test the isotropy of the local Universe using low redshift supernova data from various catalogs and the non-parametric method of smoothed residuals. Using a recently developed catalog which combines supernova data from various surveys, we show that the isotropic hypothesis of a Universe with zero velocity perturbation can be rejected with moderate significance, with $p$-value $\sim 0.07$ out to… ▽ More We test the isotropy of the local Universe using low redshift supernova data from various catalogs and the non-parametric method of smoothed residuals. Using a recently developed catalog which combines supernova data from various surveys, we show that the isotropic hypothesis of a Universe with zero velocity perturbation can be rejected with moderate significance, with $p$-value $\sim 0.07$ out to redshift $z < 0.045$. We estimate the direction of maximal anisotropy on the sky for various pre-existing catalogs and show that it remains relatively unaffected by the light curve fitting procedure. However the recovered direction is biased by the underlying distribution of data points on the sky. We estimate both the uncertainty and bias in the direction by creating mock data containing a randomly oriented bulk flow and using our method to reconstruct its direction. We conclude that the inhomogeneous nature of the data introduces a directional bias in galactic latitude of approximately $|Δb_{\rm max}| \sim 18^{\circ}$ for the supernova catalog considered in this work, and after correcting for this effect we infer the direction of maximum anisotropy as $(b,\ell) = (20^{\circ}, 276^{\circ}) \pm (12^{\circ},29^{\circ})$ in galactic coordinates. Finally we compare the anisotropic signal in the data to mock realisations in which large scale velocity perturbations are consistently accounted for at the level of linear perturbation theory. We show that including the effect of the velocity perturbation in our mock catalogs degrades the significance of the anisotropy considerably, with $p$-value increasing to $\sim 0.29$. One can conclude from our analysis that there is a moderate deviation from isotropy in the supernova data, but the signal is consistent with a large scale bulk velocity expected within $Λ$CDM. △ Less

Submitted 11 June, 2015; v1 submitted 21 October, 2014; originally announced October 2014.

Comments: 13 pages, 3 figures, 4 tables, discussions extended, main results unchanged, matches published version in ApJ

Journal ref: Astrophys. J. 801 (2015) 2, 76

Testing Isotropy in the Local Universe

Authors: Stephen Appleby, Arman Shafieloo

Abstract: We test the isotropy of the local distribution of galaxies using the 2MASS extended source catalogue. By decomposing the full sky survey into distinct patches and using a combination of photometric and spectroscopic redshift data, we use both parametric and non-parametric methods to obtain the shape of the luminosity function in each patch. We use the shape of the luminosity function to test the s… ▽ More We test the isotropy of the local distribution of galaxies using the 2MASS extended source catalogue. By decomposing the full sky survey into distinct patches and using a combination of photometric and spectroscopic redshift data, we use both parametric and non-parametric methods to obtain the shape of the luminosity function in each patch. We use the shape of the luminosity function to test the statistical isotropy of the underlying galaxy distribution. The parametric estimator shows some evidence of a hemispherical asymmetry in the north/south Galactic plane. However the non-parametric estimator exhibits no significant anisotropy, with the galaxy distribution being consistent with the assumption of isotropy in all regions considered. The parametric asymmetry is attributed to the relatively poor fit of the functional form to the underlying data. When using the non-parametric estimator, we do find a dipole in the shape of the luminosity function, with maximal deviation from isotropy at galactic coordinate $(b,l)=(30^{\circ},315^{\circ})$. However we can ascribe no strong statistical significance to this observation. △ Less

Submitted 28 October, 2014; v1 submitted 19 May, 2014; originally announced May 2014.

Comments: 12 pages, 8 figures, 5 tables. Paper significantly modified, multiple estimators of the LF now used and a new method of selecting maximally anisotropic direction presented and implemented. Matches final version published in JCAP

Journal ref: JCAP 10 (2014) 070

Testing Local Anisotropy Using the Method of Smoothed Residuals I - Methodology

Authors: Stephen Appleby, Arman Shafieloo

Abstract: We discuss some details regarding the method of smoothed residuals, which has recently been used to search for anisotropic signals in low-redshift distance measurements (Supernovae). In this short note we focus on some details regarding the implementation of the method, particularly the issue of effectively detecting signals in data that are inhomogeneously distributed on the sky. Using simulated… ▽ More We discuss some details regarding the method of smoothed residuals, which has recently been used to search for anisotropic signals in low-redshift distance measurements (Supernovae). In this short note we focus on some details regarding the implementation of the method, particularly the issue of effectively detecting signals in data that are inhomogeneously distributed on the sky. Using simulated data, we argue that the original method proposed in Colin et al. [1] will not detect spurious signals due to incomplete sky coverage, and that introducing additional Gaussian weighting to the statistic as in [2] can hinder its ability to detect a signal. Issues related to the width of the Gaussian smoothing are also discussed. △ Less

Submitted 7 May, 2014; v1 submitted 12 December, 2013; originally announced December 2013.

Comments: 8 pages, 2 figures, 1 table, references added, matches final version published in JCAP

Journal ref: JCAP 03 (2014) 007

Spherical collapse and halo mass function in f(R) theories

Authors: Michael Kopp, Stephen A. Appleby, Ixandra Achitouv, Jochen Weller

Abstract: We compute the critical density of collapse for spherically symmetric overdensities in a class of f(R) modified gravity models. For the first time we evolve the Einstein, scalar field and non-linear fluid equations, making the minimal simplifying assumptions that the metric potentials and scalar field remain quasi-static throughout the collapse. Initially evolving a top hat profile, we find that t… ▽ More We compute the critical density of collapse for spherically symmetric overdensities in a class of f(R) modified gravity models. For the first time we evolve the Einstein, scalar field and non-linear fluid equations, making the minimal simplifying assumptions that the metric potentials and scalar field remain quasi-static throughout the collapse. Initially evolving a top hat profile, we find that the density threshold for collapse depends significantly on the initial conditions imposed, specifically the choice of size and shape. By imposing `natural' initial conditions, we obtain a fitting function for the spherical collapse delta_c as a function of collapse redshift, mass of the overdensity and f_{R0}, the background scalar field value at z=0. By extending delta_c into drifting and diffusing barrier within the context of excursion set theory, we obtain a realistic mass function that might be used to confront this class of scalar-tensor models with observations of dark matter halos. The proposed analytic formula for the halo mass function was tested against Monte Carlo random walks for a wide class of moving barriers and can therefore be applied to other modified gravity theories. △ Less

Submitted 23 November, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: 24 pages, 13 figures, corresponds to published version

Journal ref: Phys. Rev. D 88, 084015 (2013)

Cluster Probes of Dark Energy Clustering

Authors: Stephen A. Appleby, Eric V. Linder, Jochen Weller

Abstract: Cluster abundances are oddly insensitive to canonical early dark energy. Early dark energy with sound speed equal to the speed of light cannot be distinguished from a quintessence model with the equivalent expansion history for $z<2$ but negligible early dark energy density, despite the different early growth rate. However, cold early dark energy, with a sound speed much smaller than the speed of… ▽ More Cluster abundances are oddly insensitive to canonical early dark energy. Early dark energy with sound speed equal to the speed of light cannot be distinguished from a quintessence model with the equivalent expansion history for $z<2$ but negligible early dark energy density, despite the different early growth rate. However, cold early dark energy, with a sound speed much smaller than the speed of light, can give a detectable signature. Combining cluster abundances with cosmic microwave background power spectra can determine the early dark energy fraction to 0.3 % and distinguish a true sound speed of 0.1 from 1 at 99 % confidence. We project constraints on early dark energy from the Euclid cluster survey, as well as the Dark Energy Survey, using both current and projected Planck CMB data, and assess the impact of cluster mass systematics. We also quantify the importance of dark energy perturbations, and the role of sound speed during a crossing of $w=-1$. △ Less

Submitted 29 May, 2013; originally announced May 2013.

Probing Dark Energy Anisotropy

Authors: Stephen A. Appleby, Eric V. Linder

Abstract: Wide area cosmological surveys enable investigation of whether dark energy properties are the same in different directions on the sky. Cosmic microwave background observations strongly restrict any dynamical effects from anisotropy, in an integrated sense. For more local constraints we compute limits from simulated distance measurements for various distributions of survey fields in a Bianchi I ani… ▽ More Wide area cosmological surveys enable investigation of whether dark energy properties are the same in different directions on the sky. Cosmic microwave background observations strongly restrict any dynamical effects from anisotropy, in an integrated sense. For more local constraints we compute limits from simulated distance measurements for various distributions of survey fields in a Bianchi I anisotropic universe. We then consider the effects of fitting for line of sight properties where isotropic dynamics is assumed (testing the accuracy through simulations) and compare sensitivities of observational probes for anisotropies, from astrophysical systematics as well as dark energy. We also point out some interesting features of anisotropic expansion in Bianchi I cosmology. △ Less

Submitted 4 February, 2013; v1 submitted 30 October, 2012; originally announced October 2012.

Comments: 9 pages, 5 figures; v2 corrected equations, conclusions unchanged, matches PRD version

Journal ref: Phys. Rev. D 87, 023532 (2013)

Fab 5: Noncanonical Kinetic Gravity, Self Tuning, and Cosmic Acceleration

Authors: Stephen A. Appleby, Antonio De Felice, Eric V. Linder

Abstract: We investigate circumstances under which one can generalize Horndeski's most general scalar-tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic gravity terms can give rise to an accelerating universe without the addition of extra propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological constan… ▽ More We investigate circumstances under which one can generalize Horndeski's most general scalar-tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic gravity terms can give rise to an accelerating universe without the addition of extra propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological constant under control. This nonlinear approach leads to new properties that may be instructive for exploring the behaviors of gravity. △ Less

Submitted 20 August, 2012; originally announced August 2012.

Comments: 16 pages, 7 figures

Journal ref: JCAP 1210 (2012) 060

Cosmology and fundamental physics with the Euclid satellite

Authors: Luca Amendola, Stephen Appleby, David Bacon, Tessa Baker, Marco Baldi, Nicola Bartolo, Alain Blanchard, Camille Bonvin, Stefano Borgani, Enzo Branchini, Clare Burrage, Stefano Camera, Carmelita Carbone, Luciano Casarini, Mark Cropper, Claudia deRham, Cinzia di Porto, Anne Ealet, Pedro G. Ferreira, Fabio Finelli, Juan Garcia-Bellido, Tommaso Giannantonio, Luigi Guzzo, Alan Heavens, Lavinia Heisenberg , et al. (38 additional authors not shown)

Abstract: Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of… ▽ More Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission. △ Less

Submitted 22 October, 2015; v1 submitted 6 June, 2012; originally announced June 2012.

Comments: 236 pages, minor edits to match the journal version 2013

Journal ref: Living Rev. Relativity 16, (2013), 6

Galileons on Trial

Authors: Stephen A. Appleby, Eric V. Linder

Abstract: Galileon gravity is a robust theoretical alternative to general relativity with a cosmological constant for explaining cosmic acceleration, with interesting properties such as having second order field equations and a shift symmetry. While either its predictions for the cosmic expansion or growth histories can approach standard ΛCDM, we demonstrate the incompatibility of both doing so simultaneous… ▽ More Galileon gravity is a robust theoretical alternative to general relativity with a cosmological constant for explaining cosmic acceleration, with interesting properties such as having second order field equations and a shift symmetry. While either its predictions for the cosmic expansion or growth histories can approach standard ΛCDM, we demonstrate the incompatibility of both doing so simultaneously. Already current observational constraints can severely disfavor an entire class of Galileon gravity models that do not couple directly to matter, ruling them out as an alternative to ΛCDM. △ Less

Submitted 28 August, 2012; v1 submitted 19 April, 2012; originally announced April 2012.

Comments: v2 matches JCAP version

Journal ref: JCAP08(2012)026

The Paths of Gravity in Galileon Cosmology

Authors: Stephen A. Appleby, Eric V. Linder

Abstract: Galileon gravity offers a robust gravitational theory for explaining cosmic acceleration, having a rich phenomenology of testable behaviors. We explore three classes of Galileon models -- standard uncoupled, and linearly or derivatively coupled to matter -- investigating the expansion history with particular attention to early time and late time attractors, as well as the linear perturbations. Fro… ▽ More Galileon gravity offers a robust gravitational theory for explaining cosmic acceleration, having a rich phenomenology of testable behaviors. We explore three classes of Galileon models -- standard uncoupled, and linearly or derivatively coupled to matter -- investigating the expansion history with particular attention to early time and late time attractors, as well as the linear perturbations. From the relativistic and nonrelativistic Poisson equations we calculate the generalizations of the gravitational strength (Newton's constant), deriving its early and late time behavior. By scanning through the parameters we derive distributions of the gravitational strength at various epochs and trace the paths of gravity in its evolution. Using ghost-free and stability criteria we restrict the allowed parameter space, finding in particular that the linear and derivative coupled models are severely constrained by classical instabilities in the early universe. △ Less

Submitted 19 April, 2012; v1 submitted 8 December, 2011; originally announced December 2011.

Comments: References added, minor modifications to the text (results unchanged). Accepted for publication in JCAP

Journal ref: JCAP 1203:043,2012