Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Spinon continuum in the Heisenberg quantum chain compound Sr 2 ⁢V 3 ⁢O 9

Magnetic excitations in the spin chain candidate Sr 2 ⁢V 3 ⁢O 9 have been investigated by inelastic neutron scattering on a single crystal sample. A spinon continuum with a bandwidth of ~22 meV is observed along the chain formed by alternating magnetic V 4+ and nonmagnetic V 5+ ions, which reveals the importance of the orbital degree of freedom in determining the chain axis as identified by prior electronic structure calculations. Incipient magnetic Bragg peaks due to weak ferromagnetic interchain couplings emerge when approaching the magnetic transition at T N ~ 5.3 K, while the excitations remain gapless within the instrumental resolution. Comparisons to the Bethe ansatz, density matrix renormalization group calculations, and effective field theories confirm Sr 2 ⁢V 3 ⁢O 9 as a host of weakly coupled S = 1/2 chains dominated by antiferromagnetic intrachain interactions of ~7.1(1) meV.

1-dimensional spin chains↗
DOE OSTI2024

Extending the Takagi-Taupin equations for x-ray nanobeam Bragg coherent diffraction

We present an approach for simulating x-ray nanobeam Bragg coherent diffraction patterns based on the Takagi-Taupin equations. Compared to conventional methods, the current approach can be universally applied to any weakly strained system including semi-infinite crystals that diffract dynamically. It addresses issues such as the curved wave front and redivergence of the focused incident beam. In conclusion, we show excellent agreement with experimental data for a strained La 0.7 ⁢Sr 0.3⁢ MnO 3 thin film on a SrTiO 3 substrate and a path to extracting physical information using automatic differentiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical Study of a Dual Representation of the Integer Quantum Hall Transition

Here, we study the critical properties of the noninteracting integer quantum Hall to insulator transition (IQHIT) in a “dual” composite-fermion (CF) representation. A key advantage of the CF representation over electron coordinates is that at criticality CF states are delocalized at all energies. The CF approach thus enables us to study the transition from a new vantage point. Using a lattice representation of CF mean-field theory, we compute the critical and multifractal exponents of the IQHIT. We obtain ν = 2.56 ± 0.02 and η = 0.51 ± 0.01, both of which are consistent with the predictions of the Chalker-Coddington network model formulated in the electron representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Magnetization Process of Atacamite: A Case of Weakly Coupled S = 1/2 Sawtooth Chains

We present a combined experimental and theoretical study of the mineral atacamite Cu 2 Cl(OH) 3 . Density-functional theory yields a Hamiltonian describing anisotropic sawtooth chains with weak 3D connections. Experimentally, we fully characterize the antiferromagnetically ordered state. Here, magnetic order shows a complex evolution with the magnetic field, while, starting at 31.5 T, we observe a plateaulike magnetization at about M sat /2. Based on complementary theoretical approaches, we show that the latter is unrelated to the known magnetization plateau of a sawtooth chain. Instead, we provide evidence that the magnetization process in atacamite is a field-driven canting of a 3D network of weakly coupled sawtooth chains that form giant moments.

1-dimensional spin chains↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Dual Mechanism for Transient Capacitance Anomaly in Improper Ferroelectrics

Negative capacitance (NC) effects in ferroelectrics can potentially break fundamental limits of power dissipation known as “Boltzmann tyranny.” However, the origin of transient NC of ferroelectrics, which is attributed to two different mechanisms involving free-energy landscape and nucleation, is under intense debate. Here, we report the coexistence of transient NC and an S-shaped anomaly during the switching of ferroelectric hexagonal ferrites capacitor in an RC circuit. The early-stage NC arises from the nucleation process, while the late-stage S-shaped anomaly corresponds to a nascent NC associated with the free-energy landscape. The entire waveform can be reproduced using a hybrid model that simultaneously incorporates these two mechanisms. These results highlight the multivariable free-energy landscape of hexagonal ferrites that enables an abrupt change of the internal field and demonstrate that the two mechanisms are not mutually exclusive, resolving the long-standing debate. In conclusion, the behavior of the S-shaped anomaly also provides a pathway to extract parameters of free-energy landscape and switching dynamics.

Capacitance↗
DOE OSTI2022

Transitions in the morphology and critical stresses of gliding dislocations in multiprincipal element alloys

Refractory multiprincipal element alloys (MPEAs) are promising material candidates for high-temperature, high-strength applications. However, their deformation mechanisms, particularly at the dislocation scale, are unusual compared to those of conventional alloys, making it challenging to understand the origin of their high strength. Here, using an atomistically informed phase-field dislocation dynamics model, we study transitions in the morphology and critical stresses of long gliding screw dislocations over extended distances in a refractory MPEA. The model MPEA crystal accounts for the atomic-scale fluctuations in chemical composition across the glide planes via spatially correlated lattice energies and for the differences in glide resistances between screw and edge dislocations of unit length. Here, we show that the dislocation moves in a stop-start motion, alternating between wavy morphology in free flight and nearly recovered straight screw orientation in full arrest. The periods of wavy glide are due to variable kink-pair formation and migration rates along the length, where portions with higher rates glide more quickly. The critical stress to initiate motion corresponds to the stress required to form and migrate a kink pair at the weakest region along the length of the dislocation. Heterogeneity in lattice energy leads to variability in the local stress-strain response and to a strain hardening-like response, in which the critical stress to reactivate glide increases with glide distance. Statistical assessment of hundreds of realizations of dislocations indicates that the amount of hardening directly scales with the dispersion in underlying lattice energy.

3-dimensional systems↗
DOE OSTI2025

Quantifying the impact of precision errors on quantum approximate optimization algorithms

The quantum approximate optimization algorithm (QAOA) is a hybrid quantum-classical algorithm that seeks to achieve approximate solutions to optimization problems by iteratively alternating between intervals of controlled quantum evolution. Here, we examine the effect of analog precision errors on QAOA performance from the perspective of both algorithmic training and performance guarantees. Leveraging cumulant expansions, we recast the faulty QAOA as a control problem in which precision errors are expressed as multiplicative control noise and derive bounds on the performance of QAOA. We show using both analytical techniques and numerical simulations that fixed precision implementations of QAOA circuits are subject to an exponential degradation in performance dependent upon the number of optimal QAOA layers and magnitude of the precision error. Despite this significant reduction, we show that it is possible to mitigate precision errors in QAOA via digitization of the variational parameters at the cost of increasing circuit depth.

quantum algorithms↗
DOE OSTI2021

Many-Body Physics in the NISQ Era: Quantum Programming a Discrete Time Crystal

Recent progress in the realm of noisy intermediate-scale quantum (NISQ) devices represents an exciting opportunity for many-body physics by introducing new laboratory platforms with unprecedented control and measurement capabilities. We explore the implications of NISQ platforms for many-body physics in a practical sense: we ask which physical phenomena, in the domain of quantum statistical mechanics, they may realize more readily than traditional experimental platforms. While a universal quantum computer can simulate any system, the eponymous noise inherent to NISQ devices practically favors certain simulation tasks over others in the near term. As a particularly well-suited target, we identify discrete time crystals (DTCs), novel nonequilibrium states of matter that break time translation symmetry. These can only be realized in the intrinsically out-of-equilibrium setting of periodically driven quantum systems stabilized by disorder-induced many-body localization. While promising precursors of the DTC have been observed across a variety of experimental platforms—ranging from trapped ions to nitrogen-vacancy centers to NMR crystals—none have all the necessary ingredients for realizing a fully fledged incarnation of this phase, and for detecting its signature long-range spatiotemporal order. We show that a new generation of quantum simulators can be programmed to realize the DTC phase and to experimentally detect its dynamical properties, a task requiring extensive capabilities for programmability, initialization, and readout. Specifically, the architecture of Google’s Sycamore processor is a remarkably close match for the task at hand. We also discuss the effects of environmental decoherence, and how they can be distinguished from ‘internal’ decoherence coming from closed-system thermalization dynamics. Already with existing technology and noise levels, we find that DTC spatiotemporal order would be observable over hundreds of periods, with parametric improvements to come as the hardware advances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗
DOE OSTI2023

Factorization Machine Learning for Disaggregation of Transmission Load Profiles with High Penetration of Behind-the-Meter Solar

The ever-growing high penetration of ubiquitously distributed energy resources, especially behind-the-meter solar (BTM) generations, has significant impacts on nodal load (i.e., net injection) profiles and consequently caused imperative operational challenges to system operators such as regional transmission organizations (RTOs). Illustrated by real-world nodal data and examples at PJM Interconnection, this paper first discusses the application and necessity of effectively extracting daily nodal load profiles in a non-intrusive manner. More importantly, a novel bi-level architecture, including Factorization Machines (FM) learning procedure has been proposed to effectively disaggregate not only one node but every node in an RTO service territory. Specifically, FM leaning is adopted to capture the interconnections between related features to better utilize the correlation between buses in the same region and between a single bus and the zonal load. The proposed bi-level technique is numerically validated using real-world, minute-level, normalized, and anonymized nodal data at PJM service territory.

behind the meter solar, load disaggregation, load ↗
DOE OSTI2020

Computationally efficient Bayesian quantum state tomography

We describe a method for Bayesian quantum state estimation combining efficient parameterization, a pseudo-likelihood, and advanced numerical sampling techniques. Examples reveal significant computational speedup, indicating the approach’s promise in practical quantum state tomography.

Lukens, Joseph↗