Search NASA⌕ Search

SEARCH · Search NASA

Results for “small bodies”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Concentration-Driven Slowdown of Chain-Exchange Kinetics in Block Copolymer Micelle Solutions

Here, we investigate the dependence of chain-exchange kinetics on polymer concentration for spherical micelles prepared from a polystyrene-b-poly(ethylene-alt-propylene) (SEP) diblock copolymer ( M n = 75 kDa, Ð = 1.05, and f PS = 0.18) in the PEP-selective solvent squalane. Solutions with concentrations of 15 vol% or more adopt a body-centered-cubic (BCC) packing of micelles, or a dense liquid-like packing (LLP), as confirmed by small-angle X-ray scattering. Concurrently, the mean aggregation number of the micelles grows steadily from about 60 in dilute solution to nearly 200 by 50 vol%. Time-resolved small-angle neutron scattering (TR-SANS) analysis reveals a monotonic and substantial (i.e., more than four orders of magnitude) decrease in the rate of chain exchange as the polymer concentration increases from 1 to 50 vol%. The TR-SANS relaxation functions are analyzed, taking into account both changes in contrast (due to chain exchange) and evolution in the structure factor (due to refinement of micelle packing). The activation energy for chain pullout is extracted using an established model and can be expressed as a linear function of polymer concentration, indicating that, in addition to the enthalpic barrier for extracting the core block, the evolution of micelle characteristics with increasing concentration imposes further constraints on chain exchange. Consistent with a theory of Halperin, these results suggest that the dominant factor reflects steric penalties during core block pullout arising from a reduced interfacial area per chain and associated increased corona crowding, due to the increase in aggregation number. The rate of exchange is independent of whether the micelles adopt a BCC or LLP arrangement, confirming that the rate-limiting step is extraction of the core block.

chemical structure↗
DOE OSTI2024

Topography Controls Variability in Circumpolar Permafrost Thaw Pond Expansion

Abstract One of the most conspicuous signals of climate change in high‐latitude tundra is the expansion of ice wedge thermokarst pools. These small but abundant water features form rapidly in depressions caused by the melting of ice wedges (i.e., meter‐scale bodies of ice embedded within the top of the permafrost). Pool expansion impacts subsequent thaw rates through a series of complex positive and negative feedbacks which play out over timescales of decades and may accelerate carbon release from the underlying sediments. Although many local observations of ice wedge thermokarst pool expansion have been documented, analyses at continental to pan‐Arctic scales have been rare, hindering efforts to project how strongly this process may impact the global carbon cycle. Here we present one of the most geographically extensive and temporally dense records yet compiled of recent pool expansion, in which changes to pool area from 2008 to 2020 were quantified through satellite‐image analysis at 27 survey areas (measuring 10–35 km 2 each, or 400 km 2 in total) dispersed throughout the circumpolar tundra. The results revealed instances of rapid expansion at 44% (15%) of survey areas. Considered alone, the extent of departures from historical mean air temperatures did not account for between site variation in rates of change to pool area. Pool growth was most clearly associated with upland (i.e., hilly) terrain and elevated silt content at soil depths greater than one meter. These findings suggest that, at short time scales, pedologic and geomorphologic conditions may exert greater control on pool dynamics in the warming Arctic than spatial variability in the rate of air temperature increases.

Abolt, C. J.↗
DOE OSTI2022

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗
DOE OSTI2020

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

$S^5$: Tidal Disruption in Crater 2 and Formation of Diffuse Dwarf Galaxies in the Local Group

We present results of a spectroscopic campaign around the diffuse dwarf galaxy Crater 2 (Cra2) and its tidal tails as part of the Southern Stellar Stream Spectroscopic Survey ($S^5$). Cra2 is a Milky Way dwarf spheroidal satellite with extremely cold kinematics, but a huge size similar to the Small Magellanic Cloud, which may be difficult to explain within collisionless cold dark matter. We identify 143 Cra2 members, of which 114 belong to the galaxy's main body and 29 are deemed part of its stellar stream. We confirm that Cra2 is dynamically cold (central velocity dispersion $2.51^{+0.33}_{-0.30}\,{\rm km\,s^{-1}}$) and also discover a $\approx$7$σ$ velocity gradient consistent with its tidal debris track. We separately estimate the stream velocity dispersion to be $5.74^{+0.98}_{-0.83}\,{\rm km\,s^{-1}}$. We develop a suite of $N$-body simulations with both cuspy and cored density profiles on a realistic Cra2 orbit to compare with $S^5$ observations. We find that the velocity dispersion ratio between Cra2 stream and galaxy ($2.30^{+0.41}_{-0.35}$) is difficult to reconcile with a cuspy halo with fiducial concentration and an initial mass predicted by standard stellar mass$-$halo mass relationships. Instead, either a cored halo with relatively small core radius or a low-concentration cuspy model can reproduce this ratio. Despite tidal mass loss, Cra2 is metal-poor ($\langle \rm[Fe/H]\rangle=-2.16\pm0.04$) compared to the stellar mass$-$metallicity relation for its luminosity. Other diffuse dwarf galaxies similar to Cra2 in the Local Group (Antlia 2 and Andromeda 19) also challenge galaxy formation models. Finally, we discuss possible formation scenarios for Cra2, including ram-pressure stripping of a gas-rich progenitor combined with tides.

Limberg, Guilherme [Chicago U., KICP; Chicago U.] ↗
DOE OSTI2023

Constraints on sub-terrestrial free-floating planets from Subaru microlensing observations

ABSTRACT The abundance of protoplanetary bodies ejected from their parent star system is presently poorly constrained. With only two existing optical observations of interstellar objects in the 108–1010 kg mass range and a small number of robust microlensing observations of free-floating planets (FFPs) in the 1024–1025 kg mass range, there is a large range of masses for which there are no existing measurements of the unbound population. The three primary microlensing surveys currently searching for FFPs operate at a cadence greater than 15 min, which limits their ability to observe events associated with bodies with a mass much below an Earth mass. We demonstrate that existing high-cadence observations of M31 with the Subaru Hyper Suprime-Cam place constraints on the abundance of unbound objects at sub-terrestrial masses, with peak sensitivity at 10−4 M⊕ for Milky Way lenses and 10−1 M⊕ for lenses in M31. For a fiducial $\frac{dn}{dM}\propto M^{-2}$ mass distribution, we find that the abundance of unbound objects is constrained to $n_\text{unbound} \lt 1.4 \times 10^{7} ~\rm {pc}^{-3}$ for masses within 1 dex of 10−4 M⊕. Additionally, we compute limits on an artificial ‘monochromatic’ distribution of unbound objects and compare to existing literature, demonstrating that the assumed spatial distribution of lenses has very significant consequences for the sensitivity of microlensing surveys. While the observations ultimately do not probe abundances suggested by current models of planetary formation, our limits place direct observational constraints on the unbound population in the sub-terrestrial mass range and motivate new observational strategies for microlensing surveys.

Astronomy & Astrophysics↗
DOE OSTI2021

The Dutch–Belgian Registry of Stereotactic Body Radiation Therapy for Liver Metastases: Clinical Outcomes of 515 Patients and 668 Metastases

Although various studies have reported that stereotactic body radiation therapy (SBRT) for liver metastases has high local control rates and relatively low toxicity, most series included a small number of patients. We aimed to validate these outcomes in a large multi-institution patient cohort treated in accordance with a common protocol.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

On the Detection of Habitable Trojan Planets in the Kepler Circumbinary Systems

We present the results of a study of the prospect of detecting habitable Trojan planets in the Kepler Habitable Zone circumbinary planetary systems (Kepler-16, -47, -453, -1647, and -1661). We integrated the orbits of 10,000 separate N-body systems (N = 4, 6), each with a one Earth-mass body in a randomly selected orbit near the L {sub 4} and L {sub 5} Lagrangian points of the host HZ circumbinary planet. We find that stable Trojan planets are restricted to a narrow range of semimajor axes in all five systems and limited to small eccentricities in Kepler-16, -47, and -1661. To assess the prospect of the detection of these habitable Trojan planets, we calculated the amplitudes of the variations they cause in the transit timing of their host bodies. Results show that the mean amplitudes of the transit timing variations (TTVs) correlate with the mass of the transiting planet and range from 70 minutes for Kepler-16b to 390 minutes for Kepler-47c. Our analysis indicates that the TTVs of the circumbinary planets caused by these Trojan bodies fall within the detectable range of timing precision obtained from the Kepler telescope’s long-cadence data. The latter points to Kepler data as a viable source to search for habitable Trojan planets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Additionally, Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bayesian constraints on the origin and geology of exoplanetary material using a population of externally polluted white dwarfs

ABSTRACT White dwarfs that have accreted planetary bodies are a powerful probe of the bulk composition of exoplanetary material. In this paper, we present a Bayesian model to explain the abundances observed in the atmospheres of 202 DZ white dwarfs by considering the heating, geochemical differentiation, and collisional processes experienced by the planetary bodies accreted, as well as gravitational sinking. The majority (>60 per cent) of systems are consistent with the accretion of primitive material. We attribute the small spread in refractory abundances observed to a similar spread in the initial planet-forming material, as seen in the compositions of nearby stars. A range in Na abundances in the pollutant material is attributed to a range in formation temperatures from below 1000 K to higher than 1400 K, suggesting that pollutant material arrives in white dwarf atmospheres from a variety of radial locations. We also find that Solar System-like differentiation is common place in exoplanetary systems. Extreme siderophile (Fe, Ni, or Cr) abundances in eight systems require the accretion of a core-rich fragment of a larger differentiated body to at least a 3σ significance, whilst one system shows evidence that it accreted a crust-rich fragment. In systems where the abundances suggest that accretion has finished (13/202), the total mass accreted can be calculated. The 13 systems are estimated to have accreted masses ranging from the mass of the Moon to half that of Vesta. Our analysis suggests that accretion continues for 11 Myrs on average.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Infinitely rugged intra-cage potential energy landscape in metallic glasses caused by many-body interaction

The absence of translational symmetry in glassy materials poses a significant challenge in establishing effective structure-property relationships in real space. Consequently, the potential energy landscape (PEL) in phase space is widely utilized to comprehend the complex phenomena in glasses. The classical PEL features a two-scale profile comprising mega-basins and sub-basins, corresponding to α-relaxations (e.g. glass transition) and β-relaxations (e.g. local cage-breaking atomic rearrangements), respectively. Recent studies, however, reveal that sub-basins are not smooth and contain finer structures, the origins of which remain elusive. Here we probe the smoothness of sub-basin bottoms in glasses' PEL by introducing small intra-cage cyclic loading and then measuring the net changes in atomic-level stresses. Compared to glasses with pair interaction, glasses with many-body interaction exhibit orders-of-magnitude larger and loading-dependent stress changes even before the first cage-breaking event takes place, which reflect much more feature-rich sub-basins. We further demonstrate this stark contrast stems from the spatial distribution of individual atom's constraining force field. Specifically, at vanishing perturbations, many-body interactions disrupt the positive-definite synchrony in energy variations of the perturbed atom and the whole system, causing inherently less confined atomic responses and infinitely rugged sub-basins. The implications of these findings for the selective addition or removal of fine structures in the PEL and the subsequent tuning of glassy materials' responses to external stimuli are also explored.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Universal scaling laws and density slopes for dark matter haloes

Smalls scale challenges suggest some missing pieces in our current understanding of dark matter. In this work, a cascade theory for dark matter is proposed to provide extra insights, similar to the cascade phenomenon in hydrodynamic turbulence. The kinetic energy is cascaded in dark matter from small to large scales involves a constant rate $ε_u$ (≈–4.6×10 –7 m 2 / s 3 ). Confirmed by N-body simulations, the energy cascade leads to a two-thirds law for kinetic energy $v^2_r$ on scale $\textit{r}$ such that $v^2_r∝(ε_ur)^{2/3}$. Equivalently, a four-thirds law can be established for mean halo density $ρ_s$ enclosed in the scale radius $r_s$ such that $ρ_s ∝ε^{2/3}_uG^{–1}r^{–4/3}_s$, which was confirmed by galaxy rotation curves. Critical properties of dark matter might be obtained by identifying key constants on relevant scales. First, the largest halo scale $r_l$ can be determined by $–u^3_0/ε_u$, where u0 is the velocity dispersion. Second, the smallest scale rη is dependent on the nature of dark matter. For collisionless dark matter, $r_η∝(–Gℏ/ε_u)^{1/3}$ ≈ 10 –13 m , where ℏ is the Planck constant. An uncertainty principle for momentum and acceleration fluctuations is also postulated. For self-interacting dark matter, $r_η∝ε^2_uG^{–3}(σ/m)^3$, where σ/m is the cross-section of interaction. On halo scale, the energy cascade leads to an asymptotic density slope γ = –4/3 for fully virialized haloes with a vanishing radial flow, which might explain the nearly universal halo density. Based on the continuity equation, halo density is analytically shown to be closely dependent on the radial flow and mass accretion, such that simulated haloes can have different limiting slopes. A modified Einasto density profile is proposed accordingly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The cross-section measurement for the 3 H($e, e'K^+)nnΛ$ reaction

The small binding energy of the hypertriton leads to predictions of the non-existence of bound hypernuclei for isotriplet three-body systems such as nnΛ. However, invariant mass spectroscopy at GSI has reported events that may be interpreted as the bound nnΛ state. The nnΛ state was sought by missing-mass spectroscopy via the (e, e'K + ) reaction at Jefferson Lab’s experimental Hall A. The present experiment has higher sensitivity to the nnΛ-state investigation in terms of better precision by a factor of about three. The analysis shown in this article focuses on the derivation of the reaction cross-section for the 3 H(γ*, K + )X reaction. Events that were detected in an acceptance, where a Monte Carlo simulation could reproduce the data well (|δp/p| < 4%), were analyzed to minimize the systematic uncertainty. No significant structures were observed with the acceptance cuts, and the upper limits of the production cross-section of the nnΛ state were obtained to be 21 and $31 nbsr -1 at the 90% confidence level when theoretical predictions of (-B Λ , Γ) = (0.25, 0.8) MeV and (0.55, 4.7) MeV, respectively, were assumed. The cross-section result provides valuable information for examining the existence of nnΛ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Feasibility of an Accelerometer-Based Structural Health Monitoring System for the LANL Blast Tube

A modeling- and simulation-based study was conducted on the feasibility of implementing an accelerometer-based SHM system on the Los Alamos National Laboratory blast tube. A blast tube experiment was modeled using the Abaqus explicit finite element solver. A custom user subroutine was written to apply test-like pressure loading to the inside surface of the blast tube. The subroutine applies analytically defined pressure loads derived from tracer output taken from a Compressible Flow Computational Fluid Dynamics Solver model of the blast tube. Five unique versions of the model were created: an undamaged reference model at 65°F was used as the baseline and compared to equivalent models at 10°F and 100°F. These three models were compared to models with small damage at the reference temperature. The two types of damage considered were a radial (circumferential) crack in the main tube body and a longitudinal crack in the supports. Acceleration outputs were extracted from accelerometer bodies included in the model and were post processed using a variety of standard SHM techniques. Different potential features signaling failure were extracted and compared using statistical methods in the time and frequency domains. A method was identified that clearly shows that differences in structural response resulting from the modeled damage can be differentiated from the structural response resulting from changing environmental conditions. However, the amount of damage applied to create observable differences in the accelerometer data was so large that simpler methods of damage detection would be more cost effective in locating damage.

42 ENGINEERING↗
DOE OSTI2024

Dynamical instability in multi-orbiter systems with gas friction

ABSTRACT Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

Li, Jiaru (ORCID:0000000155507421)↗