Search NASA⌕ Search

Engineering topics

Gaztanaga, Enrique

Publications and source records attributed to Gaztanaga, Enrique.

The DESI N -body Simulation Project – II. Suppressing sample variance with fast simulations

ABSTRACT Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim 20\, h^{-3}\, \mathrm{Gpc}^{3}$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. AbacusSummit is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\, h^{-3}\, \mathrm{Gpc}^{3}$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-N-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k \lt 0.3\, h\, \mathrm{Mpc}^{-1}$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in AbacusSummit suite to extend the effective volume ∼20 times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

cosmology: theory↗

The PAU Survey: narrow-band photometric redshifts using Gaussian processes

Here, we study the performance of the hybrid template machine learning photometric redshift (photo- z ) algorithm delight , which uses Gaussian processes, on a subset of the early data release of the Physics of the Accelerating Universe Survey (PAUS). We calibrate the fluxes of the 40 PAUS narrow bands with six broad-band fluxes ( uBVriz ) in the Cosmic Evolution Survey (COSMOS) field using three different methods, including a new method that utilizes the correlation between the apparent size and overall flux of the galaxy. We use a rich set of empirically derived galaxy spectral templates as guides to train the Gaussian process, and we show that our results are competitive with other standard photometric redshift algorithms. delight achieves a photo- z 68th percentile error of σ 68 = 0.0081(1 + z ) without any quality cut for galaxies with i auto < 22.5 as compared to 0.0089(1 + z ) and 0.0202(1 + z ) for the bpz and annz 2 codes, respectively. delight is also shown to produce more accurate probability distribution functions for individual redshift estimates than bpz and annz 2. Common photo- z outliers of delight and bcnz 2 (previously applied to PAUS) are found to be primarily caused by outliers in the narrow-band fluxes, with a small number of cases potentially indicating spectroscopic redshift failures in the reference sample. In the process, we introduce performance metrics derived from the results of bcnz 2 and delight , allowing us to achieve a photo- z quality of σ 68 < 0.0035(1 + z ) at a magnitude of i auto < 22.5 while keeping 50 per cent objects of the galaxy sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

The PAU Survey: an improved photo- z sample in the COSMOS field

Here we present – and make publicly available – accurate and precise photometric redshifts in the ACS footprint from the COSMOS field for objects with i AB ≤ 23. The redshifts are computed using a combination of narrow-band photometry from PAUS, a survey with 40 narrow bands spaced at $100\,\mathring{\rm A}$ intervals covering the range from 4500 to $8500\,\mathring{\rm A}$, and 26 broad, intermediate, and narrow bands covering the UV, visible and near-infrared spectrum from the COSMOS2015 catalogue. We introduce a new method that models the spectral energy distributions as a linear combination of continuum and emission-line templates and computes its Bayes evidence, integrating over the linear combinations. The correlation between the UV luminosity and the $\mathrm{O\,{\small II}}$ line is measured using the 66 available bands with the zCOSMOS spectroscopic sample, and used as a prior which constrains the relative flux between continuum and emission-line templates. The flux ratios between the $\mathrm{O\,{\small II}}$ line and H α , H β and $\mathrm{O\,{\small III}}$ are similarly measured and used to generate the emission-line templates. Comparing to public spectroscopic surveys via the quantity Δ z ≡ (z photo – z spec )/(1 + z spec ), we find the photometric redshifts to be more precise than previous estimates, with σ 68 (Δ z ) ≈ (0.003, 0.009) for galaxies at magnitude i AB ~ 18 and i AB ~ 23, respectively, which is three times and 1.66 times tighter than COSMOS2015. Additionally, we find the redshifts to be very accurate on average, yielding a median of the Δ z distribution compatible with |median(Δ z )| ≤ 0.001 at all redshifts and magnitudes considered. Both the added PAUS data and new methodology contribute significantly to the improved results. The catalogue produced with the technique presented here is expected to provide a robust redshift calibration for current and future lensing surveys, and allows one to probe galaxy formation physics in an unexplored luminosity-redshift regime, thanks to its combination of depth, completeness, and excellent redshift precision and accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Installation of the Dark Energy Spectroscopic Instrument at the Mayall 4-meter telescope

The Dark Energy Spectroscopic Instrument (DESI) is a Stage IV ground-based dark energy experiment that will measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 square degrees will be measured during the life of the experiment. We describe the installation of the major elements of the instrument at the Mayall 4m telescope, completed in late 2019. The previous prime focus corrector, spider vanes, and upper rings were removed from the Mayall’s Serrurier truss and replaced with the newlyconstructed DESI ring, vanes, cage, hexapod, and optical corrector. The new corrector was optically aligned with the primary mirror using a laser tracker system. The DESI focal plane system was integrated to the corrector, with each of its ten 500-fiber-positioner petal segments installed using custom installation hardware and the laser tracker. Ten DESI spectrographs with 30 cryostats were installed in a newly assembled clean room in the Large Coude Room. The ten cables carrying 5000 optical fibers from the positioners in the focal plane were routed down the telescope through cable wraps at the declination and hour angle axes, and their integral slitheads were integrated with the ten spectrographs. The fiber view camera assembly was installed to the Mayall’s primary mirror cell. Servers for the instrument control system replaced existing computer equipment. The fully integrated instrument has been commissioned and is ready to start its operations phase.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Performance of Kitt Peak’s Mayall 4-meter telescope during DESI commissioning

In preparation for the Dark Energy Spectroscopic Instrument (DESI), a new top end was installed on the Mayall 4-meter telescope at Kitt Peak National Observatory. The refurbished telescope and the DESI instrument were successfully commissioned on sky between 2019 October and 2020 March. Here we describe the pointing, tracking and imaging performance of the Mayall telescope equipped with its new DESI prime focus corrector, as measured by six guider cameras sampling the outer edge of DESI's focal plane. Analyzing ~500,000 guider images acquired during commissioning, we find a median delivered image FWHM of 1.1 arcseconds (in the r-band at 650 nm), with the distribution extending to a best-case value of ~0.6 arcseconds. The point spread function is well characterized by a Moffat profile with a power-law index of β ~ 3.5 and little dependence of β on FWHM. The shape and size of the PSF delivered by the new corrector at a field angle of 1.57 degrees are very similar to those measured with the old Mayall corrector on axis. Here, we also find that the Mayall achieves excellent pointing accuracy (several arcseconds RMS) and minimal open-loop tracking drift (< 1 milliarcsecond per second), improvements on the telecope's pre-DESI performance. In the future, employing DESI's active focus adjustment capabilities will likely further improve the Mayall/DESI delivered image quality.

47 OTHER INSTRUMENTATION↗

Discovery of a Candidate Binary Supermassive Black Hole in a Periodic Quasar from Circumbinary Accretion Variability

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ~10 –3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67–002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302–102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.

79 ASTRONOMY AND ASTROPHYSICS↗

Candidate Periodically Variable Quasars from the Dark Energy Survey and the Sloan Digital Sky Survey

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg 2 overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-year long multi-color ( griz ) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to r ~23.5 mag) powered by less massive black holes (with masses ≳ 10 8.5 M ⊙ ) at high redshift for the first time. We find five candidates with significant (at >99.74% single-frequency significance in at least two bands with a global p-value of ~7 × 10 -4 –3× 10 -3 accounting for the look-elsewhere effect) periodicity with observed periods of ~3–5 years (i.e., 1–2 years in rest frame) having ~4–6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of |${\sim }0.8^{+0.5}_{-0.3}$| % or |${\sim }1.1^{+0.7}_{-0.5}$| quasar per deg 2 . Our detection rate is 4–80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. Furthermore, we discuss implications on the future direct detection of low-frequency gravitational waves. Additionally, continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Energy Survey Identification of A Low-Mass Active Galactic Nucleus at Redshift 0.823 from Optical Variability

We report the identification of a low-mass active galactic nucleus (AGN), DES J0218−0430, in a redshift z = 0.823 galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218−0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 yr. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg ii and broad H β lines, confirming its nature as a broad-line AGN. Archival XMM–Newton X-ray observations suggest an intrinsic hard X-ray luminosity of |$L_{{\rm 2-12\, keV}}\approx 7.6\pm 0.4\times 10^{43}$| erg s^−1, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H β from SDSS spectrum, we estimate a virial black hole (BH) mass of M_• ≈ 10^6.43–10^6.72 M_⊙ (with the error denoting the systematic uncertainty from different calibrations), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift > 0.4 detected in optical. We estimate the host galaxy stellar mass to be M_* ≈ 10^10.5 ± 0.3 M_⊙ based on modelling the multiwavelength spectral energy distribution. DES J0218−0430 extends the M_•–M_* relation observed in luminous AGNs at z ∼ 1 to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗

The Curious Case of PHL 293B: A Long-lived Transient in a Metal-poor Blue Compact Dwarf Galaxy

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy’s unusual spectral features (P Cygni-like profiles with ~800 km s -1 blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of ~0.1 mag in the g band offset by 100 ± 21 pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or nonterminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Furthermore, our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at the Vera C. Rubin Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral variability of a sample of extreme variability quasars and implications for the MgII broad-line region

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed bymore than 1.5mag in the g band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey epochs (separated by a few years in the quasar rest frame). These EVQs are selected from quasars in the SDSS Stripe 82 region, covering a redshift range of 0.5 < z < 2.1. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the g band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on time-scales of only months. To increase sample statistics, we use a supplemental sample of 33 EVQs with multi-epoch spectra from SDSS that cover the broad Mg II lambda 2798 line. Leveraging on the large dynamic range in continuum variability between the multi-epoch spectra, we explore the associated variations in the broad MgII line, whose variability properties have not been well studied before. The broad MgII flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the full width at half-maximum (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad MgII to estimate the black hole mass with single epoch spectra therefore introduces a luminosity-dependent bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Bias and high-order galaxy correlation functions in the APM galaxy survey

On large scales, the higher order moments of the mass distribution, S(sub J) = bar-zeta(sub J)/bar-zeta(sup J-1)(sub 2), e.g., the skewness S(sub 3) and kurtosis S(sub 4), can be predicted using nonlinear perturbation theory. Comparison of these predictions with moments of the observed galaxy distribution probes the bias between galaxies and mass. Applying this method to models with initially Gaussian fluctuations and power spectra P(k) similar to that of galaxies in the Automatic Plate Measuring (APM) survey, we find that the predicted higher order moments S(sub J)(R) are in good agreement with those directly inferred from the APM survey in the absence of bias. We use this result to place limits on the linear and nonlinear bias parameters. Models in which the extra power observed on large scales (with respect to the standard cold dark matter (CDM) model) is produced by scale-dependent bias match the APM higher order amplitudes only if nonlinear bias (rather than nonlinear gravity) generates the observed higher order moments. When normalized to Cosmic Background Explorer Differential Microwave Radiometer (COBE DMR), these models are siginificantly ruled out by the S(sub 3) observations. The cold plus hot dark matter model normalized to COBE can reproduce the APM higher order correlations if one introduces nonlinear bias terms, while the low-density CDM model with a cosmological constant does not require any bias to fit the large-scale amplitudes.

Gaztanaga, Enrique↗

The three-point function as a probe of models for large-scale structure

We analyze the consequences of models of structure formation for higher order (n-point) galaxy correlation functions in the mildly nonlinear regime. Several variations of the standard Omega = 1 cold dark matter model with scale-invariant primordial perturbations have recently been introduced to obtain more power on large scales, R(sub p) is approximately 20/h Mpc, e.g., low matter-density (nonzero cosmological constant) models, 'tilted' primordial spectra, and scenarios with a mixture of cold and hot dark matter. They also include models with an effective scale-dependent bias, such as the cooperative galaxy formation scenario of Bower et al. We show that higher-order (n-point) galaxy correlation functions can provide a useful test of such models and can discriminate between models with true large-scale power in the density field and those where the galaxy power arises from scale-dependent bias: a bias with rapid scale dependence leads to a dramatic decrease of the the hierarchical amplitudes Q(sub J) at large scales, r is greater than or approximately R(sub p). Current observational constraints on the three-point amplitudes Q(sub 3) and S(sub 3) can place limits on the bias parameter(s) and appear to disfavor, but not yet rule out, the hypothesis that scale-dependent bias is responsible for the extra power observed on large scales.

Frieman, Joshua A.↗

Redshift distortions of galaxy correlation functions

To examine how peculiar velocities can affect the two-, three-, and four-point redshift correlation functions, we evaluate volume-average correlations for configurations that emphasize and minimize redshift distortions for four different volume-limited samples from each of the CfA, SSRS, and IRAS redshift catalogs. We present the results as the correlation length r(sub 0) and power index gamma of the two-point correlations, bar-xi(sub 0) = (r(sub 0)/r)(exp gamma), and as the hierarchical amplitudes of the three- and four-point functions, S(sub 3) = bar-xi(sub 3)/bar-xi(exp 2)(sub 2) and S(sub 4) = bar-xi(sub 4)/bar-xi(exp 3)(sub 2). We find a characteristic distortion for bar-xi(sub 2), the slope gamma is flatter and the correlation length is larger in redshift space than in real space; that is, redshift distortions 'move' correlations from small to large scales. At the largest scales (up to 12 Mpc), the extra power in the redshift distribution is compatible with Omega(exp 4/7)/b approximately equal to 1. We estimate Omega(exp 4/7)/b to be 0.53 +/- 0.15, 1.10 +/- 0.16, and 0.84 +/- 0.45 for the CfA, SSRS, and IRAS catalogs. Higher order correlations bar-xi(sub 3) and bar-xi(sub 4) suffer similar redshift distortions but in such a way that, within the accuracy of our ananlysis, the normalized amplitudes S(sub 3) and S(sub 4) are insensitive to this effect. The hierarchical amplitudes S(sub 3) and S(sub 4) are constant as a function of scale between 1 and 12 Mpc and have similar values in all samples and catalogs, S(sub 3) approximately equal to 2 and S(sub 4) approximately equal to 6, despite the fact that bar-xi(sub 2), bar-xi(sub 3), and bar-xi(sub 4) differ from one sample to another by large factors (up to a factor of 4 in bar-xi(sub 2), 8 for bar-xi(sub 3), and 12 for bar-xi(sub 4)). The agreement between the independent estimations of S(sub 3) and S(sub 4) is remarkable given the different criteria in the selection of galaxies and also the difference in the resulting range of densities, luminosities, and locations between samples.

Fry, J. N.↗

Biasing and hierarchical statistics in large-scale structure

In this paper we consider the consequences for galaxy formation of an arbitrary, effectively local biasing transformation of a hierarchical underlying matter distribution. We show that a general form of such a transformation preserves the hierarchical properties and the shape of the dispersion in the limit of small fluctuations, i.e., on large scales, although the values of the hierarchical amplitudes may change arbitrarily. We present expressions for the induced hierarchical amplitudes S(g,j) of the galaxy distribution in terms of the matter amplitudes S(j) and biasing parameters for j = 3-7. For higher order correlations, j greater than 2, restricting to a linear bias is not a consistent approximation even at very large scales. To draw any conclusions from the galaxy distribution about matter correlations of order j, properties of biasing must be specified completely to order j - 1.

Fry, J. N.↗

The three-point function as a probe of models for large-scale structure

The consequences of models of structure formation for higher-order (n-point) galaxy correlation functions in the mildly non-linear regime are analyzed. Several variations of the standard Omega = 1 cold dark matter model with scale-invariant primordial perturbations were recently introduced to obtain more power on large scales, R(sub p) is approximately 20 h(sup -1) Mpc, e.g., low-matter-density (non-zero cosmological constant) models, 'tilted' primordial spectra, and scenarios with a mixture of cold and hot dark matter. They also include models with an effective scale-dependent bias, such as the cooperative galaxy formation scenario of Bower, etal. It is shown that higher-order (n-point) galaxy correlation functions can provide a useful test of such models and can discriminate between models with true large-scale power in the density field and those where the galaxy power arises from scale-dependent bias: a bias with rapid scale-dependence leads to a dramatic decrease of the hierarchical amplitudes Q(sub J) at large scales, r is approximately greater than R(sub p). Current observational constraints on the three-point amplitudes Q(sub 3) and S(sub 3) can place limits on the bias parameter(s) and appear to disfavor, but not yet rule out, the hypothesis that scale-dependent bias is responsible for the extra power observed on large scales.

Frieman, Joshua A.↗

Probing the statistics of primordial fluctuations and their evolution

The statistical distribution of fluctuations on various scales is analyzed in terms of the counts in cells of smoothed density fields, using volume-limited samples of galaxy redshift catalogs. It is shown that the distribution on large scales, with volume average of the two-point correlation function of the smoothed field less than about 0.05, is consistent with Gaussian. Statistics are shown to agree remarkably well with the negative binomial distribution, which has hierarchial correlations and a Gaussian behavior at large scales. If these observed properties correspond to the matter distribution, they suggest that our universe started with Gaussian fluctuations and evolved keeping hierarchial form.

Gaztanaga, Enrique↗

N-point correlation functions in the CfA and SSRS redshift distribution of galaxies

Using counts in cells, we estimate the volume-average N-point galaxy correlation functions for N = 2, 3, and 4, in redshift samples of the CfA and SSRS catalogs. Volume-limited samples of different sizes are used to study the uncertainties at different scales, the shot noise, and the problem with the boundaries. The hierarchical constants S3 and S4 agree well in all samples in CfA and SSRS, with average S3 = 194 +/- 0.07 and S4 = 4.56 +/- 0.53. We compare these results with estimates obtained from angular catalogs and recent analysis over IRAS samples. The amplitudes SJ seem larger in real space than in redshift space, although the values from the angular analysis correspond to smaller scales, where we might expect larger nonperturbative effects. It is also found that S3 and S4 are smaller for IRAS than for optical galaxies. This, together with the fact that IRAS galaxies have smaller amplitude for the above correlation functions, indicates that the density fluctuations of IRAS galaxies cannot be simply proportional to the density fluctuations of optical galaxies, i.e., biasing has to be nonlinear between them.

Gaztanaga, Enrique↗