skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms

Authors:
; ; ; ; ;
Publication Date:
Research Org.:
SLAC National Accelerator Lab., Menlo Park, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1328994
Report Number(s):
SLAC-PUB-16839
Journal ID: ISSN 2470-0010; arXiv:1603.04405
DOE Contract Number:
AC02-76SF00515
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physical Review D; Journal Volume: 93; Journal Issue: 10
Country of Publication:
United States
Language:
English
Subject:
Astrophysics,ASTRO, GRQC

Citation Formats

Schmittfull, Marcel, /BCCP, Berkeley /UC, Berkeley /LBNL, Berkeley, Vlah, Zvonimir, /Stanford U., ITP /KIPAC, Menlo Park /SLAC, McDonald, Patrick, and /LBNL, Berkeley. Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms. United States: N. p., 2017. Web. doi:10.1103/PhysRevD.93.103528.
Schmittfull, Marcel, /BCCP, Berkeley /UC, Berkeley /LBNL, Berkeley, Vlah, Zvonimir, /Stanford U., ITP /KIPAC, Menlo Park /SLAC, McDonald, Patrick, & /LBNL, Berkeley. Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms. United States. doi:10.1103/PhysRevD.93.103528.
Schmittfull, Marcel, /BCCP, Berkeley /UC, Berkeley /LBNL, Berkeley, Vlah, Zvonimir, /Stanford U., ITP /KIPAC, Menlo Park /SLAC, McDonald, Patrick, and /LBNL, Berkeley. Mon . "Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms". United States. doi:10.1103/PhysRevD.93.103528. https://www.osti.gov/servlets/purl/1328994.
@article{osti_1328994,
title = {Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms},
author = {Schmittfull, Marcel and /BCCP, Berkeley /UC, Berkeley /LBNL, Berkeley and Vlah, Zvonimir and /Stanford U., ITP /KIPAC, Menlo Park /SLAC and McDonald, Patrick and /LBNL, Berkeley},
abstractNote = {},
doi = {10.1103/PhysRevD.93.103528},
journal = {Physical Review D},
number = 10,
volume = 93,
place = {United States},
year = {Mon Mar 06 00:00:00 EST 2017},
month = {Mon Mar 06 00:00:00 EST 2017}
}
  • Cited by 13
  • The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional fast Fourier transforms. This is exact and a fewmore » orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. In conclusion, the method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.« less
  • Disk interleaving (disk striping) distributes fragments of a data block across a group of disks. For applications with regular I/O reference patterns and requiring large block transfers, interleaving can speed data transfers and reduce I/O times. Computing fast Fourier transforms is one such application; the algorithm assumes that the entire array to be transformed should fit in the main memory. Occasionally arrays exceed the capacity of the main memory and reside in secondary storage. The authors have used synchronous and asynchronous disk striping to compute very large three-dimensional FFTs. They present a two-pass algorithm for computing data stored on interleavedmore » disks and analyze the I/O times. A large FFT example provides a quantitative view of the advantages of disk interleaving for this application.« less
  • A change of variables is proposed with the following property: given a {ital preselected} {ital density}, when a {ital uniform} grid in the new coordinates is mapped into Cartesian coordinates, the resulting grid points are distributed according to the preselected density. The preselected density is chosen so that the grid points are concentrated in the atomic-core regions. The grids generated in this way are useful for fully numerical or basis-set-free electronic-structure calculations, in particular those using plane-wave expansions: they give a good treatment of the nuclear cusps (because of the high density of grid points in the core regions) {italmore » and} the fast Fourier transform can be used for numerical interpolation (because the grids are uniform in the new coordinates). When compared with uniform or flat grids in Cartesian coordinates, the present scheme is clearly superior: it gives numerical integration errors that are typically two orders of magnitude smaller, and it very accurately interpolates atomic and molecular densities, even for atoms like uranium, which have a very high core density.« less
  • A viscoplastic approach using the Fast Fourier Transform (FFT) method for obtaining local mechanical response is utilized to study microstructure-property relationships in composite materials. Specifically, three-dimensional, two-phase digital materials containing isotropically coarsened particles surrounded by a matrix phase, generated through a Kinetic Monte Carlo Potts model for Ostwald ripening, are used as instantiations in order to calculate the stress and strain rate fields under uniaxial tension. The effects of the morphology of the matrix phase, the volume fraction and the contiguity of particles, and the polycrystallinity of matrix phase, on the stress and strain rate fields under uniaxial tension aremore » examined. It is found that the first moments of the stress and strain rate fields have a different dependence on the particle volume fraction and the particle contiguity from their second moments. The average stresses and average strain rates of both phases and of the overall composite have rather simple relationships with the particle volume fraction whereas their standard deviations vary strongly, especially when the particle volume fraction is high, and the contiguity of particles has a noticeable effect on the mechanical response. It is also found that the shape of stress distribution in the BCC hard particle phase evolves as the volume fraction of particles in the composite varies, such that it agrees with the stress field in the BCC polycrystal as the volume of particles approaches unity. Finally, it is observed that the stress and strain rate fields in the microstructures with a polycrystalline matrix are less sensitive to changes in volume fraction and contiguity of particles.« less