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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Extended nonergodic regime and spin subdiffusion in disordered SU(2)-symmetric Floquet systems

In this work, we explore thermalization and quantum dynamics in a one-dimensional disordered SU(2)-symmetric Floquet model, where a many-body localized phase is prohibited by the non-Abelian symmetry. Despite the absence of localization, we find an extended nonergodic regime at strong disorder where the system exhibits nonthermal behaviors. In the strong disorder regime, the level spacing statistics exhibit neither a Wigner-Dyson nor a Poisson distribution, and the spectral form factor does not show a linear-in-time growth at early times characteristic of random matrix theory. The average entanglement entropy of the Floquet eigenstates is subthermal, although violating an area-law scaling with system sizes. We further compute the expectation value of local observables and find strong deviations from the eigenstate thermalization hypothesis. The infinite-temperature spin autocorrelation function decays at long times as $t^{–β}$ with $\textit{β}$ < 0.5, indicating subdiffusive transport at strong disorders.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Femtojoule optical nonlinearity for deep learning with incoherent illumination

Optical neural networks (ONNs) are a promising computational alternative for deep learning due to their inherent massive parallelism for linear operations. However, the development of energy-efficient and highly parallel optical nonlinearities, a critical component in ONNs, remains an outstanding challenge. Here, we introduce a nonlinear optical microdevice array (NOMA) compatible with incoherent illumination by integrating the liquid crystal cell with silicon photodiodes at the single-pixel level. We fabricate NOMA with more than half a million pixels, each functioning as an optical analog of the rectified linear unit at ultralow switching energy down to 100 femtojoules per pixel. With NOMA, we demonstrate an optical multilayer neural network. Our work holds promise for large-scale and low-power deep ONNs, computer vision, and real-time optical image processing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

50 Years of quantum chromodynamics

Quantum Chromodynamics, the theory of quarks and gluons, whose interactions can be described by a local SU(3) gauge symmetry with charges called “color quantum numbers”, is reviewed; the goal of this review is to provide advanced Ph.D. students a comprehensive handbook, helpful for their research. When QCD was “discovered” 50 years ago, the idea that quarks could exist, but not be observed, left most physicists unconvinced. Then, with the discovery of charmonium in 1974 and the explanation of its excited states using the Cornell potential, consisting of the sum of a Coulomb-like attraction and a long range linear confining potential, the theory was suddenly widely accepted. This paradigm shift is now referred to as the November revolution. It had been anticipated by the observation of scaling in deep inelastic scattering, and was followed by the discovery of gluons in three-jet events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Collaborative Research: The Evolution of Magnetic Skeletons During 3D Reconnection. Final report

We investigated how the plasmoid instability affects the evolution of topological features of a magnetic field and the process of magnetic reconnection with three sub-projects. (i) We extensively studied the criterion for plasmoid-mediated current sheet disruption and the onset of fast reconnection. Scaling relations with respect to plasma parameters have been established, using both direct numerical simulations and analytic theories, for the disruption time, the current sheet width, the instability linear growth rate, and the dominant mode wavenumber. (ii) We applied the theories of plasmoid-mediated current sheet disruption to magnetohydrodynamic turbulence and predicted a steepening of the turbulence energy spectrum. This steepening has been confirmed by direct numerical simulations. (iii) We performed numerical simulations of ultraviolet burst events in the solar atmosphere with a spine-dome magnetic topology. The plasmoid instability leads to unsteady reconnection, which is consistent with solar observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Random Circuits in the Black Hole Interior

In this paper, we present a quantitative holographic relation between a microscopic measure of randomness and the geometric length of the wormhole in the black hole interior. To this end, we perturb an AdS black hole with Brownian semiclassical sources, implementing the continuous version of a random quantum circuit for the black hole. We use the random circuit to prepare ensembles of states of the black hole whose semiclassical duals contain Einstein-Rosen (ER) caterpillars: long cylindrical wormholes with large numbers of matter inhomogeneities, of linearly growing length with the circuit time. In this setup, we show semiclassically that the ensemble of ER caterpillars of average length $k\ell_Δ$ and matter correlation scale $\ell_Δ$ forms an approximate quantum state $k$-design of the black hole. At exponentially long circuit times, the ensemble of ER caterpillars becomes polynomial-copy indistinguishable from a collection of random states of the black hole. We comment on the implications of these results for holographic circuit complexity and for the holographic description of the black hole interior.

FOS: Physical sciences↗
DOE OSTI2025

Development of mesopores in superfine grain graphite neutron-irradiated at high fluence

Microstructural changes induced by neutron irradiation of superfine grain graphite G347A (Tokai Carbon, Japan) were examined by nitrogen adsorption at 77K and by three microscopy techniques (SEM, TEM and FIB-SEM tomography). The specimens were irradiated at doses of up to 30 dpa, covering stages before and after the turnaround fluence at three temperatures (300, 450, 750 °C) of their irradiation envelope. The initial graphite densification at low fluences did not produce any detectable effect in the pore size range (<350 nm) measured by gas adsorption. However, graphite irradiated at high fluences, after turnaround, showed severe structural changes. At all three temperatures and high irradiation fluences, gas adsorption revealed significant increase of the volume of narrow mesopores (<5–20nm) and up to five times increase of BET surface area, both in linear relationship with the relative volume expansion. Analysis of microscopy images showed multiplication of fine macropores (>50 nm) at high irradiation fluences and more structural changes on multiple scales, from nanometers to microns. This work demonstrates the unique ability of gas adsorption techniques to analyze open pores with sizes between sub-nanometer and sub-micron in bulk nuclear graphite, with supporting microscopy results.

I. Contescu, Cristian↗
DOE OSTI2026

Scalable Quantum Monte Carlo Method for Polariton Chemistry via Mixed Block Sparsity and Tensor Hypercontraction Method

We present a reduced-scaling auxiliary-field quantum Monte Carlo (AFQMC) framework designed for large molecular systems and ensembles, with or without coupling to optical cavities. Our approach leverages the natural block sparsity of the Cholesky decomposition (CD) of electron repulsion integrals in molecular ensembles and employs tensor hypercontraction (THC) to efficiently compress low-rank Cholesky blocks. By representing the Cholesky vectors in a mixed format, keeping high-rank blocks in block-sparse form and compressing low-rank blocks with THC, we reduce the scaling of exchange-energy evaluation from quartic to robust cubic in the number of molecular orbitals N, while lowering memory from cubic toward quadratic. Benchmark analyses on one-, two-, and three-dimensional molecular ensembles (up to ∼1,200 orbitals) show that (a) the number of nonzeros in Cholesky tensors grows linearly with system size across dimensions; (b) the average numerical rank increases sublinearly and does not saturate at these sizes; and (c) rank heterogeneity─some blocks nearly full rank and many low rank, naturally motivates the proposed mixed block sparsity and THC scheme for efficient calculation of exchange energy. In conclusion, we demonstrate that the mixed scheme yields cubic wall-time scaling with favorable prefactors and preserves AFQMC accuracy.

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

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗
DOE OSTI2025

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constraining the baryonic feedback with cosmic shear using the DES Year-3 small-scale measurements

ABSTRACT We use the small scales of the Dark Energy Survey (DES) Year-3 cosmic shear measurements, which are excluded from the DES Year-3 cosmological analysis, to constrain the baryonic feedback. To model the baryonic feedback, we adopt a baryonic correction model and use the numerical package baccoemu to accelerate the evaluation of the baryonic non-linear matter power spectrum. We design our analysis pipeline to focus on the constraints of the baryonic suppression effects, utilizing the implication given by a principal component analysis on the Fisher forecasts. Our constraint on the baryonic effects can then be used to better model and ameliorate the effects of baryons in producing cosmological constraints from the next-generation large-scale structure surveys. We detect the baryonic suppression on the cosmic shear measurements with a ∼2σ significance. The characteristic halo mass for which half of the gas is ejected by baryonic feedback is constrained to be $M_c \gt 10^{13.2} \, h^{-1} \, \mathrm{M}_{\odot }$ (95 per cent C.L.). The best-fitting baryonic suppression is $\sim 5{{\ \rm per\ cent}}$ at $k=1.0 \, {\rm Mpc}\ h^{-1}$ and $\sim 15{{\ \rm per\ cent}}$ at $k=5.0 \, {\rm Mpc} \ h^{-1}$. Our findings are robust with respect to the assumptions about the cosmological parameters, specifics of the baryonic model, and intrinsic alignments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Superdiffusive transport of energy in one-dimensional metals

Metals in one spatial dimension are described at the lowest energy scales by the Luttinger liquid theory. It is well understood that this free theory, and even interacting integrable models, can support ballistic transport of conserved quantities including energy. In contrast, realistic one-dimensional metals, even without disorder, contain integrability-breaking interactions that are expected to lead to thermalization and conventional diffusive linear response. We argue that the expansion of energy when such a nonintegrable Luttinger liquid is locally heated above its ground state shows superdiffusive behavior (i.e., spreading of energy that is intermediate between diffusion and ballistic propagation), by combining an analytical anomalous diffusion model with numerical matrix-product–state calculations on a specific perturbed spinless fermion chain. Different metals will have different scaling exponents and shapes in their energy spreading, but the superdiffusive behavior is stable and should be visible in time-resolved experiments.

Bulchandani, Vir B.↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2021

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. Additionally, we compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Resolving the Inner Parsec of the Blazar J1924–2914 with the Event Horizon Telescope

The blazar J1924–2914 is a primary Event Horizon Telescope (EHT) calibrator for the Galactic center’s black hole Sagittarius A*. Here we present the first total and linearly polarized intensity images of this source obtained with the unprecedented 20 μas resolution of the EHT. J1924–2914 is a very compact flat-spectrum radio source with strong optical variability and polarization. In April 2017 the source was observed quasi-simultaneously with the EHT (April 5–11), the Global Millimeter VLBI Array (April 3), and the Very Long Baseline Array (April 28), giving a novel view of the source at four observing frequencies, 230, 86, 8.7, and 2.3 GHz. These observations probe jet properties from the subparsec to 100 pc scales. We combine the multifrequency images of J1924–2914 to study the source morphology. We find that the jet exhibits a characteristic bending, with a gradual clockwise rotation of the jet projected position angle of about 90° between 2.3 and 230 GHz. Linearly polarized intensity images of J1924–2914 with the extremely fine resolution of the EHT provide evidence for ordered toroidal magnetic fields in the blazar compact core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Binary operations on neuromorphic hardware with application to linear algebraic operations and stochastic equations

Abstract Non-von Neumann computational hardware, based on neuron-inspired, non-linear elements connected via linear, weighted synapses—so-called neuromorphic systems—is a viable computational substrate. Since neuromorphic systems have been shown to use less power than CPUs for many applications, they are of potential use in autonomous systems such as robots, drones, and satellites, for which power resources are at a premium. The power used by neuromorphic systems is approximately proportional to the number of spiking events produced by neurons on-chip. However, typical information encoding on these chips is in the form of firing rates that unarily encode information. That is, the number of spikes generated by a neuron is meant to be proportional to an encoded value used in a computation or algorithm. Unary encoding is less efficient (produces more spikes) than binary encoding. For this reason, here we present neuromorphic computational mechanisms for implementing binary two’s complement operations. We use the mechanisms to construct a neuromorphic, binary matrix multiplication algorithm that may be used as a primitive for linear differential equation integration, deep networks, and other standard calculations. We also construct a random walk circuit and apply it in Brownian motion simulations. We study how both algorithms scale in circuit size and iteration time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantum time dynamics employing the Yang-Baxter equation for circuit compression

Quantum time dynamics (QTD) is considered a promising problem for quantum supremacy on near-term quantum computers. However, QTD quantum circuits grow with increasing time simulations. This study focuses on simulating the time dynamics of one-dimensional (1D) integrable spin chains with nearest-neighbor interactions. We have proved the existence of a reflection symmetry in the quantum circuit employed for simulating the time evolution of certain classes of 1D Heisenberg model Hamiltonians by virtue of the quantum Yang-Baxter equation, and how this symmetry can be exploited to compress and produce a shallow quantum circuit. With this compression scheme, the depth of the quantum circuit becomes independent of step size and only depends on the number of spins. We show that the depth of the compressed circuit is rigorously a linear function of the system size for the studied Heisenberg model Hamiltonians in the present work. As a consequence, the number of CNOT gates in the compressed circuit only scales quadratically with the system size, which allows for the simulations of time dynamics of very large 1D spin chains. We derive the compressed circuit representations for different special cases of the Heisenberg Hamiltonian. We compare and demonstrate the effectiveness of this approach by performing simulations on quantum computers.

1-dimensional spin chains↗
DOE OSTI2022

Observation of nonlinear planar Hall effect in magnetic-insulator–topological-insulator heterostructures

Interfacing topological insulators (TIs) with magnetic insulators (MIs) have been widely used to study the interaction between topological surface states and magnetism. Previous transport studies typically interpret the suppression of weak antilocalization or appearance of the anomalous Hall effect as signatures of the magnetic proximity effect (MPE) imposed to TIs. Here, in this work, we report the observation of the nonlinear planar Hall effect (NPHE) in Bi 2 ⁢Se 3 films grown on MI thulium and yttrium-iron-garnet (TmIG and YIG) substrates, which is an order of magnitude larger than that in Bi 2 ⁢Se 3 grown on nonmagnetic gadolinium-gallium-garnet (GGG) substrate. The nonlinear Hall resistance in TmIG/Bi 2 ⁢Se 3 depends linearly on the external magnetic field, while that in YIG/Bi 2 ⁢Se 3 exhibits an extra hysteresis loop around zero field. The magnitude of the NPHE is found to scale inversely with carrier density. We speculate that the observed NPHE is related to the MPE-induced exchange gap opening and out-of-plane spin textures in the TI surface states, which may be used as an alternative transport signature of the MPE in MI/TI heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗