Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Sequential optical response suppression for chemical mixture characterization

The characterization of mixtures of non-interacting, spectroscopically similar quantum components has important applications in chemistry, biology, and materials science. We introduce an approach based on quantum tracking control that allows for determining the relative concentrations of constituents in a quantum mixture, using a single pulse which enhances the distinguishability of components of the mixture and has a length that scales linearly with the number of mixture constituents. To illustrate the method, we consider two very distinct model systems: mixtures of diatomic molecules in the gas phase, as well as solid-state materials composed of a mixture of components. A set of numerical analyses are presented, showing strong performance in both settings.

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

An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations

We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the Illustris-TNG suite. We create an HOD mock with Illustris-like correlations using our method, constraining the associated IA model parameters, with the $\mathcal{X}$$^{2}_{dof}$ between our model’s correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation–position and orientation–orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to 0.1 $\mathcal{h}$ –1 Mpc. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Coupled cluster theory on modern heterogeneous supercomputers

This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory—a linear-scaling, massively parallel framework—as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.

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

CLUE: A Fast Parallel Clustering Algorithm for High Granularity Calorimeters in High-Energy Physics

One of the challenges of high granularity calorimeters, such as that to be built to cover the endcap region in the CMS Phase-2 Upgrade for HL-LHC, is that the large number of channels causes a surge in the computing load when clustering numerous digitized energy deposits (hits) in the reconstruction stage. In this article, we propose a fast and fully parallelizable density-based clustering algorithm, optimized for high-occupancy scenarios, where the number of clusters is much larger than the average number of hits in a cluster. The algorithm uses a grid spatial index for fast querying of neighbors and its timing scales linearly with the number of hits within the range considered. We also show a comparison of the performance on CPU and GPU implementations, demonstrating the power of algorithmic parallelization in the coming era of heterogeneous computing in high-energy physics.

Rovere, Marco↗
DOE OSTI2022

Aqueous titanium redox flow batteries—State-of-the-art and future potential

Market-driven deployment of inexpensive (but intermittent) renewable energy sources, such as wind and solar, in the electric power grid necessitates grid-stabilization through energy storage systems Redox flow batteries (RFBs), with their rated power and energy decoupled (resulting in a sub-linear scaling of cost), are an inexpensive solution for the efficient electrochemical storage of large amounts of electrical energy. Titanium-based RFBs, first developed by NASA in the 1970s, are an interesting albeit less examined chemistry and are the focus of the present review. Ti, constituting 0.6% of the Earth’s crust and an ingredient in inexpensive white paints, is amongst the few elements (V and Mn being some others) which exhibit multiple soluble oxidation states in aqueous electrolytes. Further, the very high (approaching 10 M) solubility of Ti in low pH solutions suggests the possibility of developing exceptionally high energy density aqueous Redox Flow Batteries systems. With these advantages in mind, we present the state-of the-art in Ti-RFBs with a focus on Ti/Mn, Ti/Fe and Ti/Ce couples and systems that use Ti as an additive (such as Ti/V/Mn). The inherent advantages of inexpensive Ti actives and relatively high energy density is contrasted with potential side-reactions resulting in reduced energy efficiency. Technological pathways are presented with a view to overcoming critical bottlenecks and a vision is presented for the future development of Ti-RFBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pulse Length Monitor for Breakdown Diagnostics in THz and Mm-Wave Accelerators

The development of novel high-gradient accelerating structures operating at THz frequencies is critical for future free-electron lasers and TeV scale linear colliders. To reach high energies with reasonable length requires high accelerating gradients of ~100 MV/m. The main limitation to reaching these high-energy gradients is the vacuum RF breakdown phenomenon, which disrupts normal accelerator operations. For stable operations and to understand the breakdown microscopic dynamics, a new device capable of detecting the breakdown occurrences is required. In this paper, we provide the design of a pulse length monitor based on an analog to digital converter for fast signal digitization without the need to use high-speed digitizers to be used in a commercial mm-wave heterodyne spectrometer.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

A Fast, Two-dimensional Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization

Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of data points has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions where, typically, the second dimension is small. This method scales linearly with the total number of data points when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Radius-to-frequency Mapping and FRB Frequency Drifts

We build a model of radius-to-frequency mapping in magnetospheres of neutron stars and apply it to frequency drifts observed in fast radio bursts (FRBs). We assume that an emission patch propagates along the dipolar magnetic field lines, producing coherent emission with frequency, direction, and polarization defined by the local magnetic field. The observed temporal evolution of the frequency depends on the relativistic effects of time contraction and the curvature of the magnetic field lines. The model generically produces linear scaling of the drift rate, $\dot{\omega}$ ∝-${\omega}$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dot{\omega}$ is also possible. Fast rotating magnetospheres produce higher drifts rates for similar viewing parameters than the slowly rotating ones. In the case of repeaters, the same source may show variable drift patterns depending on the observing phase. We expect rotational of polarization position angle through a burst, though by smaller amount than in radio pulsars. All of these findings compare favorably with properties of FBRs, strengthening their possible loci in the magnetospheres of neutron stars.

Magnetars↗
DOE OSTI2024

The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters

We present cosmological constraints from a gravitational lensing mass map covering 9400 deg2 reconstructed from measurements of the cosmic microwave background (CMB) made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with measurements of baryon acoustic oscillations and big bang nucleosynthesis, we obtain the clustering amplitude σ8 = 0.819 ± 0.015 at 1.8% precision, S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.840 ± 0.028, and the Hubble constant H 0 = (68.3 ± 1.1) km s -1 Mpc -1 at 1.6% precision. A joint constraint with Planck CMB lensing yields σ 8 = 0.812 ± 0.013, S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.831 ± 0.023 , and H0 = (68.1 ± 1.0) km s -1 Mpc -1 . These measurements agree with ΛCDM extrapolations from the CMB anisotropies measured by Planck. We revisit constraints from the KiDS, DES, and HSC galaxy surveys with a uniform set of assumptions and find that S8 from all three are lower than that from ACT+Planck lensing by levels ranging from 1.7σ to 2.1σ. This motivates further measurements and comparison, not just between the CMB anisotropies and galaxy lensing but also between CMB lensing probing z ~ 0.5–5 on mostly linear scales and galaxy lensing at z 0.5 on smaller scales. We combine with CMB anisotropies to constrain extensions of ΛCDM, limiting neutrino masses to Σm ν < 0.13 eV (95% c.l.), for example. We describe the mass map and related data products that will enable a wide array of cross-correlation science. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the ΛCDM model, while paving a promising path for neutrino physics with lensing from upcoming ground-based CMB surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The Atacama Cosmology Telescope: Cosmology from Cross-correlations of unWISE Galaxies and ACT DR6 CMB Lensing

We present tomographic measurements of structure growth using cross-correlations of Atacama Cosmology Telescope (ACT) DR6 and Planck cosmic microwave background (CMB) lensing maps with the unWISE Blue and Green galaxy samples, which span the redshift ranges 0.2 ≲ z ≲ 1.1 and 0.3 ≲ z ≲ 1.8, respectively. We improve on prior unWISE cross-correlations not just by making use of the new, high-precision ACT DR6 lensing maps, but also by including additional spectroscopic data for redshift calibration and by analyzing our measurements with a more flexible theoretical model. We determine the amplitude of matter fluctuations at low redshifts (z ≃ 0.2–1.6), finding S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.813 ± 0.021 using the ACT cross-correlation alone and S 8 = 0.810 ± 0.015 with a combination of Planck and ACT cross-correlations; these measurements are fully consistent with the predictions from primary CMB measurements assuming standard structure growth. The addition of baryon acoustic oscillation data breaks the degeneracy between σ 8 and Ω m , allowing us to measure σ 8 = 0.813 ± 0.020 from the cross-correlation of unWISE with ACT and σ 8 = 0.813 ± 0.015 from the combination of cross-correlations with ACT and Planck. These results also agree with the expectations from primary CMB extrapolations in ΛCDM cosmology; the consistency of σ 8 derived from our two redshift samples at z ∼ 0.6 and 1.1 provides a further check of our cosmological model. Our results suggest that structure formation on linear scales is well described by ΛCDM even down to low redshifts z ≲ 1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Rapid-cycling ReBCO dipole magnet concept for muon acceleration

We present a concept of the rapid-cycling 1.7 T ReBCO magnet for 1 to 10 kT/s ramping range as required for muon acceleration in a future muon collider presently under study. This approach is based on the 6 kA CORC-like cable constructed with 12 superconducting ReBCO tapes of 2 mm width. Based on theoretical prediction there is a linear scaling of the ReBCO cable hysteresis loss with the crossing magnetic field possibly generating power loss independent of the magnetic field ramping speed. This feature makes this ReBCO cable suitable for the construction of rapid-cycling magnets. In this work we outline the design of a dipole magnet for 1.7 T gap magnetic field and discuss the helium coolant parameters for the minimal electric power.

Piekarz, Henryk [Fermilab] (ORCID:0000000194225662↗
DOE OSTI2025

A new technique to retrieve aerosol vertical profiles using micropulse lidar and ground-based aerosol measurements

Accurately characterizing the vertical distribution of aerosols and their cloud-forming properties is crucial for understanding aerosol-cloud interactions and their impact on climate. This study presents a novel technique for retrieving vertical profiles of aerosols, cloud condensation nuclei (CCN), and ice nucleating particles (INP) by combining micropulse lidar, radiosonde, and ground-based aerosol measurements. Herein, the technique is applied to data collected by our team at Texas A&M University during the Tracking Aerosol Convection Interactions ExpeRiment (TRACER) campaign. Ground-based aerosol size distribution and CCN counter data are used to estimate the value of the aerosol hygroscopicity parameter, κ. The derived κ, together with Mie scattering theory and the relative humidity profile from the radiosonde, is used to estimate aerosol size growth and the associated increase in backscatter at each altitude. We then correct the lidar backscatter to dry conditions to produce the dry aerosol backscatter coefficient profile. The dry aerosol backscatter coefficient profile is linearly scaled to collocated surface measurements of aerosols, CCN, and INP to produce corresponding vertical profiles. Combining lidar backscatter profiles with aerosol and cloud nucleation measurements leads to a more realistic representation of vertical distributions of aerosol properties. The method could be readily applied to lidar measurements in future field campaigns.

Chen, Bo [Texas A & M Univ., College Station, TX (↗
DOE OSTI2023

Methane, ethane, and propane production in Greenland ice core samples and a first isotopic characterization of excess methane

Air trapped in polar ice provides unique records of the past atmospheric composition ranging from key greenhouse gases such as methane (CH 4 ) to short-lived trace gases like ethane (C 2 H 6 ) and propane (C 3 H 8 ). Recently, the comparison of CH 4 records obtained using different extraction methods revealed disagreements in the CH 4 concentration for the last glacial in Greenland ice. Elevated methane levels were detected in dust-rich ice core sections measured discretely, pointing to a process sensitive to the melt extraction technique. To shed light on the underlying mechanism, we performed targeted experiments and analyzed samples for methane and the short-chain alkanes ethane and propane covering the time interval from 12 to 42 kyr. Here, we report our findings of these elevated alkane concentrations, which scale linearly with the amount of mineral dust within the ice samples. The alkane production happens during the melt extraction step of the classic wet-extraction technique and reaches 14 to 91 ppb of CH 4 excess in dusty ice samples. We document for the first time a co-production of excess methane, ethane, and propane, with the observed concentrations for ethane and propane exceeding their past atmospheric background at least by a factor of 10. Independent of the produced amounts, excess alkanes were produced in a fixed molar ratio of approximately 14 : 2 : 1, indicating a shared origin. The measured carbon isotopic signature of excess methane is (-47.0 ± 2.9) ‰ and its deuterium isotopic signature is (-326 ± 57) ‰. With the co-production ratios of excess alkanes and the isotopic composition of excess methane we established a fingerprint that allows us to constrain potential formation processes. This fingerprint is not in line with a microbial origin. Moreover, an adsorption–desorption process of thermogenic gas on dust particles transported to Greenland does not appear very likely. Instead, the alkane pattern appears to be indicative of abiotic decomposition of organic matter as found in soils and plant leaves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DMTN-221: Periodicity Analysis in Alert Production

The baselined timeseries features to be computed in Alert Production include a Lomb-Scargle periodogram ran on two classes of variable systems: RR Lyrae and Eclipsing Binaries. Based on a simulated LSST-like cadence light curves taken from the Extended LSST Astronomical Time-series Classification Challenge (ELAsTiCC) we perform an end-to-end test to characterize the periodicity recovery on the Alert Production multi-band light curves. In both variable classes, we found that a single-band Lomb-Scargle implementation yields to a low fraction of recovered periods, with a significant preference on the simple periodic phenomena such as RR Lyrae. We also investigated the results from a multi-band Lomb-Scargle and found an increased fraction of recovered periodicities above 15% for the eclipsing binaries, and over 80$\%$ for the RR Lyrae stars. Our findings suggest that a multi-band Lomb-Scargle should be implemented for searching periodic phenomena through AP. We also asses the computational and scientific performance of several configurations on simulated alert data and find that our current configuration scales linearly with the number of detections while assuming an heuristic frequency grid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2021

Real time evolution for ultracompact Hamiltonian eigenstates on quantum hardware

In this work we present a detailed analysis of variational quantum phase estimation (VQPE), a method based on real-time evolution for ground and excited state estimation on near-term hardware. We derive the theoretical ground on which the approach stands, and demonstrate that it provides one of the most compact variational expansions to date for solving strongly correlated Hamiltonians. At the center of VQPE lies a set of equations, with a simple geometrical interpretation, which provides conditions for the time evolution grid in order to decouple eigenstates out of the set of time evolved expansion states, and connects the method to the classical filter diagonalization algorithm. Further, we introduce what we call the unitary formulation of VQPE, in which the number of matrix elements that need to be measured scales linearly with the number of expansion states, and we provide an analysis of the effects of noise which substantially improves previous considerations. The unitary formulation allows for a direct comparison to iterative phase estimation. Our results mark VQPE as both a natural and highly efficient quantum algorithm for ground and excited state calculations of general many-body systems. We demonstrate a hardware implementation of VQPE for the transverse field Ising model. Further, we illustrate its power on a paradigmatic example of strong correlation (Cr2 in the SVP basis set), and show that it is possible to reach chemical accuracy with as few as ~50 timesteps.

Klymko, Katherine↗