Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

2D-FFTLog: efficient computation of real-space covariance matrices for galaxy clustering and weak lensing

ABSTRACT Accurate covariance matrices for two-point functions are critical for inferring cosmological parameters in likelihood analyses of large-scale structure surveys. Among various approaches to obtaining the covariance, analytic computation is much faster and less noisy than estimation from data or simulations. However, the transform of covariances from Fourier space to real space involves integrals with two Bessel integrals, which are numerically slow and easily affected by numerical uncertainties. Inaccurate covariances may lead to significant errors in the inference of the cosmological parameters. In this paper, we introduce a 2D-FFTLog algorithm for efficient, accurate, and numerically stable computation of non-Gaussian real-space covariances for both 3D and projected statistics. The 2D-FFTLog algorithm is easily extended to perform real-space bin-averaging. We apply the algorithm to the covariances for galaxy clustering and weak lensing for a Dark Energy Survey Year 3-like and a Rubin Observatory’s Legacy Survey of Space and Time Year 1-like survey, and demonstrate that for both surveys, our algorithm can produce numerically stable angular bin-averaged covariances with the flat sky approximation, which are sufficiently accurate for inferring cosmological parameters. The code CosmoCov for computing the real-space covariances with or without the flat-sky approximation is released along with this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Intruder configurations in 29 Ne at the transition into the island of inversion: Detailed structure study of 28 Ne

Detailed γ-ray spectroscopy of the exotic neon isotope 28 Ne has been performed for the first time using the one-neutron removal reaction from 29 Ne on a liquid hydrogen target at 240 MeV/nucleon. Based on an analysis of parallel momentum distributions, a level scheme with spin-parity assignments has been constructed for 28 Ne and the negative-parity states are identified for the first time. The measured partial cross sections and momentum distributions reveal a significant intruder p-wave strength providing evidence of the breakdown of the N=20 and N=28 shell gaps. Only a weak, possible f-wave strength was observed to bound final states. Large-scale shell-model calculations with different effective interactions do not reproduce the large p-wave and small f-wave strength observed experimentally, indicating an ongoing challenge for a complete theoretical description of the transition into the island of inversion along the Ne isotopic chain.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Misaligned Wind‐Waves Behind Atmospheric Cold Fronts

Abstract Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes.

Oceanography↗
DOE OSTI2025

Detailed characterization of runaway electron driven whistler waves in low-density DIII-D discharges

RE-driven whistler waves during quiescent DIII-D shots have been investigated further. The waves are confirmed to be mostly perpendicularly propagating and are observed for the first time with frequencies up to 700 MHz. Phase-spectral analysis has been used to infer their toroidal mode numbers, n, which are expected to scale with the wavenumber, k, of the mode. Though we derive a theoretical scaling of k ≈ 4n, the measured mode numbers are found to exhibit a very weak dependence on k. In addition, increases in synchrotron emission have been found to consistently lag whistler wave bursts by roughly 3–5 ms, suggesting the waves are causing pitch-angle scattering, since the emitted synchrotron radiation is a strong function of the REs' perpendicular energy. The stronger the wave bursts, the greater the subsequent increase in synchrotron emission. A predator-prey model is used to describe these nonlinear wave-particle interactions, from which the wave damping rates and the loss parameter can be inferred. The damping rates are found to be of the order of (1.6 ± 0.8) × 10 4 /s, and the unitless loss parameter is found to be approximately 2, suggesting that the loss mechanism is diffusive. These observations will serve to validate models of RE-driven waves in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

AbacusHOD : a highly efficient extended multitracer HOD framework and its application to BOSS and eBOSS data

ABSTRACT We introduce the AbacusHOD model and present two applications of AbacusHOD and the AbacusSummit simulations to observations. AbacusHOD is a Halo Occupation Distribution (HOD) framework written in Python that is particle-based, multitracer, highly generalized, and highly efficient. It is designed specifically with multitracer/cosmology analyses for next-generation large-scale structure surveys in mind, and takes advantage of the volume and precision offered by the new state-of-the-art AbacusSummit cosmological simulations. The model is also highly customizable and should be broadly applicable to any upcoming surveys and a diverse range of cosmological analyses. In this paper, we demonstrate the capabilities of the AbacusHOD framework through two example applications. The first example demonstrates the high efficiency and the large HOD extension feature set through an analysis of full-shape redshift-space clustering of BOSS galaxies at intermediate to small scales ($\lt 30\, h^{-1}$ Mpc), assessing the necessity of introducing secondary galaxy biases (assembly bias). We find strong evidence for using halo environment instead of concentration to trace secondary galaxy bias, a result which also leads to a moderate reduction in the ‘lensing is low’ tension. The second example demonstrates the multitracer capabilities of the AbacusHOD package through an analysis of the extended Baryon Oscillation Spectroscopic Survey cross-correlation measurements between three different galaxy tracers: luminous red galaxies, emission-line galaxies, and quasi-stellar objects. We expect the AbacusHOD framework, in combination with the AbacusSummit simulation suite, to play an important role in a simulation-based analysis of the upcoming Dark Energy Spectroscopic Instrument data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Light scalars at the cosmological collider

We study the self-energies of weakly interacting scalar fields in de Sitter space with one field much lighter than the Hubble scale. We argue that self-energies drastically simplify in this light limit. We illustrate this in theories with two scalar fields, one heavy and one light, interacting with one another through either cubic or quartic interactions. To regulate infrared divergences, we compute these self-energies in Euclidean de Sitter space and then carefully analytically continue to Lorentzian signature. In particular, we do this for the most general renormalizable theory of two scalar fields with even interactions to leading order in the coupling and the mass of the light field. These self-energies are determined by de Sitter sunset diagrams, whose analytic structure and UV divergences we derive. Even at very weak couplings, the light field can substantially change how the heavy field propagates over long distances. The light field’s existence may then be inferred from how it modifies the heavy field’s oscillatory contribution to the primordial bispectrum in the squeezed limit, i.e. its cosmological collider signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Future changes in extreme precipitation over the San Francisco Bay Area: Dependence on atmospheric river and extratropical cyclone events

Extreme precipitation poses a major challenge for local governments, including the City and County of San Francisco, California, as flooding can damage and destroy infrastructure and property. As the climate continues to warm, reliable future precipitation projections are needed to provide the best possible information to decision makers. However, future changes in the magnitude of extreme precipitation are uncertain, as current state-of-the-art global climate models are typically run at relatively coarse horizontal resolutions that require the use of convective parameterization and have difficulty simulating observed extreme rainfall rates. Here, we performed ensembles of convection-permitting regional climate model simulations to investigate how five historically impactful extreme precipitation events over the San Francisco Bay Area could change if similar events occurred in future climates. We found that changes in storm-total precipitation depend strongly on storm type. Precipitation associated with an atmospheric river (AR) accompanied by an extratropical cyclone (ETC) is projected to increase at a rate exceeding (by up to 1.5 times) the theoretical Clausius Clapeyron scaling of 6–7% per °C warming. On the other hand, future precipitation changes are weak or negative for events characterized by an AR only, despite increases in precipitable water and integrated vapor transport that are similar to those of the co-occurring AR and ETC events. The differences in the sign of future precipitation change between AR-only events and co-occurring AR and ETC events is instead linked with changes in mid-tropospheric vertical velocity. Given that the majority of observed ARs are associated with an ETC, this research has important implications for future precipitation impacts over the Bay Area, as it indicates that storm-total precipitation associated with the most common type of storm event may increase by up to 26–37% in 2100 relative to historical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Broadband unidirectional visible imaging using wafer-scale nano-fabrication of multi-layer diffractive optical processors

We present a broadband and polarization-insensitive unidirectional imager that operates at the visible part of the spectrum, where image formation occurs in one direction, while in the opposite direction, it is blocked. This approach is enabled by deep learning-driven diffractive optical design with wafer-scale nano-fabrication using high-purity fused silica to ensure optical transparency and thermal stability. Our design achieves unidirectional imaging across three visible wavelengths (covering red, green, and blue parts of the spectrum), and we experimentally validated this broadband unidirectional imager by creating high-fidelity images in the forward direction and generating weak, distorted output patterns in the backward direction, in alignment with our numerical simulations. This work demonstrates wafer-scale production of diffractive optical processors, featuring 16 levels of nanoscale phase features distributed across two axially aligned diffractive layers for visible unidirectional imaging. This approach facilitates mass-scale production of ~0.5 billion nanoscale phase features per wafer, supporting high-throughput manufacturing of hundreds to thousands of multi-layer diffractive processors suitable for large apertures and parallel processing of multiple tasks. Beyond broadband unidirectional imaging in the visible spectrum, this study establishes a pathway for artificial-intelligence-enabled diffractive optics with versatile applications, signaling a new era in optical device functionality with industrial-level, massively scalable fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Compound Drought and Heatwaves at a Global Scale: The Role of Natural Climate Variability-Associated Synoptic Patterns and Land-Surface Energy Budget Anomalies

Compound drought and heatwave (CDHW) events have garnered much attention in recent studies. However, thus far, the identification of such events is oversimplified, and their association with natural climate variability is not fully explored. Here, we derive anomalies in the weekly self-calibrated Palmer Drought Severity Index (sc_PDSI) and daily maximum temperatures to identify CDHW events from 1982 to 2016 over 26 climate regions across the globe. Using a Poisson Generalized Linear Model (GLM), we analyze yearly occurrences of seasonal CDHW events and their association with the warm and cold phases of El Nino Southern Oscillation (ENSO), Pacific Decadal Oscillation (PDO), and North Atlantic Oscillation (NAO). ENSO exhibits robust association with CDHW events over the Southern Hemisphere during the austral summer and fall, while PDO influences their occurrences over the Western North America in the Northern Hemisphere during the boreal summer, which is supported by the composites of anomalies in the atmospheric circulations and surface energy budget. However, NAO association with CDHW events is relatively weak. The CDHW occurrence over other regions is driven by a combination of these large-scale natural forcing. Our analyses also highlight that the cooccurrence of weekly to submonthly scale anomalies in the observed temperature and precipitation may not be always aligned between the observations and the reanalysis. Therefore, caution must be exercised while explaining such observed anomalies on the basis of reanalysis-based circulations and surface energy budget. Lastly, our analyses provide a new insight towards concurrent extremes and should help foster research efforts in this area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Cluster cosmology with anisotropic boosts: validation of a novel forward modelling analysis and application on SDSS redMaPPer clusters

We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large-scale structure, a consequence of the limited discrimination between line-of-sight interlopers and true cluster members offered by photometric surveys. We validate our analysis via a blind cosmology challenge on mocks, and find that we can obtain tight and unbiased cosmological constraints without informative priors or external calibrations on any of our model parameters. We then apply our analysis on the SDSS redMaPPer clusters, and find results favoring low Ωm and high σ8, combining to yield the lensing strength constraint $S_8 = 0.715_{-0.021}^{+0.024}$. We investigate potential drivers behind these results through a series of post-unblinding tests, noting that our results are consistent with existing cluster cosmology constraints but clearly inconsistent with other CMB/LSS based cosmology results. From these tests, we find hints that a suppression in the cluster lensing signal may be driving our results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius, R B = 10 3 R g . The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outside R B and suppress BH accretion to 1.5% of the Bondi rate, M ̇ B . Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of M ̇ B reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% M ̇ B . Because the BAD and RAD states tangle up the jets and destroy them well inside R B , they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

Astronomy & Astrophysics↗
DOE OSTI2022

Signatures of Light Massive Relics on non-linear structure formation

ABSTRACT Cosmologies with Light Massive Relics (LiMRs) as a subdominant component of the dark sector are well-motivated from a particle physics perspective, and can also have implications for the σ8 tension between early and late time probes of clustering. The effects of LiMRs on the cosmic microwave background (CMB) and structure formation on large (linear) scales have been investigated extensively. In this paper, we initiate a systematic study of the effects of LiMRs on smaller, non-linear scales using cosmological N-body simulations; focusing on quantities relevant for photometric galaxy surveys. For most of our study, we use a particular model of non-thermal LiMRs but the methods developed generalizing to a large class of LiMR models – we explicitly demonstrate this by considering the Dodelson–Widrow velocity distribution. We find that, in general, the effects of LiMR on small scales are distinct from those of a ΛCDM universe, even when the value of σ8 is matched between the models. We show that weak lensing measurements around massive clusters, between ∼0.1 h−1Mpc and ∼10 h−1Mpc, should have sufficient signal-to-noise in future surveys to distinguish between ΛCDM and LiMR models that are tuned to fit both CMB data and linear scale clustering data at late times. Furthermore, we find that different LiMR cosmologies indistinguishable by conventional linear probes can be distinguished by non-linear probes if their velocity distributions are sufficiently different. LiMR models can, therefore, be best tested by jointly analyzing the CMB and late-time structure formation on both large and small scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

High-Performance Near-Infrared Quantum Emission from Color Centers in hBN

Color centers hosted in hexagonal boron nitride (hBN) have emerged as a highly promising platform for single-photon emission and spin-photon technologies relevant to quantum communication and quantum networking. As a wide bandgap van der Waals material, hBN can host optically active quantum defects across a broad spectral range. Here, we demonstrate a simple and scalable oxygen-plasma process that reproducibly creates single quantum emitters in hBN with blinking-free zero-phonon lines (ZPLs) spanning near-infrared (NIR) from 700 up to 971 nm. These emitters combine MHz-level brightness, single-photon purity up to 99.9%, and ultranarrow cryogenic line widths down to 2.7 GHz under quasi-resonant excitation, placing them in a particularly attractive regime for quantum photonics. Photostability measurements further reveal resistance to photobleaching, subnanometer spectral stability over long time scales, and near-shot-noise-limited intensity fluctuations. Analysis of the phonon sidebands shows weak vibronic coupling and ZPL-dominated emission, with Debye–Waller factors approaching 50%. Control experiments together with elemental mapping support oxygen incorporation as a necessary ingredient in activating the NIR emitter population, while first-principles calculations identify O N V N and O N V N H as the leading defect candidates. These results establish a high-performance NIR quantum-emitter platform in hBN for free-space quantum networking and future integrated quantum-photonic architectures.

2D materials↗
DOE OSTI2025

Unsupervised Clustering of Microseismic Events and Focal Mechanism Analysis at the CO 2 Injection Site in Decatur, Illinois

Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Willis, Rachel Marie [Sandia National Laboratories↗
DOE OSTI2020

Planck 2018 results. VI. Cosmological parameters

In this paper, we present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5 σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ω c h 2 = 0.120 ± 0.001, baryon density Ω b h 2 = 0.0224 ± 0.0001, scalar spectral index n s = 0.965 ± 0.004, and optical depth τ = 0.054 ± 0.007 (in this abstract we quote 68% confidence regions on measured parameters and 95% on upper limits). The angular acoustic scale is measured to 0.03% precision, with 100 θ * = 1.0411 ± 0.0003. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H 0 = (67.4 ± 0.5) km s –1 Mpc –1 ; matter density parameter Ω m = 0.315 ± 0.007; and matter fluctuation amplitude σ 8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be N eff = 2.99 ± 0.17, in agreement with the Standard Model prediction N eff = 3.046, and find that the neutrino mass is tightly constrained to Σ m ν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2 σ , which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, Ω K = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w 0 = –1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r 0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6 σ , tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

High mass and halo resolution from fast low resolution simulations

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. Here we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition resolution and develop a halo finder "relaxed-FoF", where we allow different linking lengths for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also calibrate small-scale velocities of small halos to improve the power spectrum in redshift space. For the matter statistics, we incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scales. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

One-dimensional neutron diffraction from layered graphite: Reciprocal space structure and grating behavior

In this work we report observations of one-dimensional neutron diffraction from highly oriented pyrolytic graphite (HOPG), where the scattering angle varies continuously with incident angle following classical grating-like behavior. The 2D polycrystalline structure of HOPG—with highly aligned layers along the 𝑐 axis but random in-plane rotations—creates planes of scattering intensity in reciprocal space at 𝑄 𝑐 =𝑛⁢(2⁢𝜋/𝑑) where 𝑑=3.35Å is the interlayer spacing. As the Ewald sphere sweeps through reciprocal space during sample rotation, it continuously intersects these planes, producing the observed angular dispersion. We observe both first-order (𝑛=1) and second-order (𝑛=2) diffraction at conventional scattering angles (25°–70°), with peak positions that remain temperature-independent between 10 K and 294 K and follow quantitative agreement with momentum conservation 𝑄 𝑐 =𝑘⁢[sin⁡𝜓−sin⁡𝜓 𝑓 ]=𝑛⁢(2⁢𝜋/𝑑). X-ray diffraction under similar conditions shows no comparable behavior, confirming that sharp nuclear-vacuum contrast is essential. While diffraction intensities are weak (∼10 −6 of Bragg peaks), the observations demonstrate how the interplay of atomic-scale periodicity, nuclear contrast, and structural disorder enables observation of continuous diffraction curves at thermal neutron wavelengths, illustrating how HOPG's unique microstructure determines its scattering properties beyond conventional Bragg diffraction.

2-dimensional systems↗