Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “body wave”

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 73 records · Page 4

Second Linear Response Theory and the Analytic Calculation of Excited-State Properties

We present a method based on second linear response time-dependent density func- tional theory (TDDFT) to calculate permanent and transition multipoles of excited states, which are required to compute excited-state absorption/emission spectra, multi- photon optical processes, among others. In previous work, we examined computations based on second linear response theory in which linear response TDDFT was employed twice. In contrast, the present methodology only requires information from a single linear response calculation to compute the excited-state properties. These are eval- uated analytically through various algebraic operations involving electron repulsion integrals and excitation vectors. The present derivation focuses on full many-body wave functions instead of single orbitals, as in our previous approach. Here, we test the proposed method by applying it to several diatomic and triatomic molecules. This shows that the computed excited-state dipoles are consistent with respect to reference equation-of-motion coupled cluster calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Seeking Repeating Anthropogenic Seismic Sources: Implications for Seismic Velocity Monitoring at Fault Zones

Abstract Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high‐frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring. In this work, we propose a systematic workflow to seek such powerful seismic sources in a rapid and straightforward manner. We tackle the problem from a statistical point of view, considering that persistent, powerful seismic sources yield highly coherent correlation functions (CFs) between pairs of seismic sensors. The algorithm is tested in California and Japan. Multiple sites close to fault zones show high‐frequency CFs stable for an extended period of time. These findings have great potential for monitoring fault zones, including the San Jacinto Fault and the Ridgecrest area in Southern California, Napa in Northern California, and faults in central Japan. However, extra steps, such as beamforming or polarization analysis, are required to determine the dominant seismic sources and study the source characteristics, which are crucial to interpreting the velocity monitoring results. Train tremors identified by the present approach have been successfully used for seismic velocity monitoring of the San Jacinto Fault in previous studies.

58 GEOSCIENCES↗

Binding and excitations in Si x H y molecular systems using quantum Monte Carlo

In this work, we present high-accuracy correlated calculations of small SixHy molecular systems in both the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within ≈0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between 1% and 3.5% with absolute values at most ≈0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PyQMC : An all-Python real-space quantum Monte Carlo module in PySCF

We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for a simple comparison between QMC calculations and other many-body wave function techniques, as well as access to high accuracy trial wave functions.

Wheeler, William A. (ORCID:0000000283910530)↗

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗

Many-Body Majorana Braiding without an Exponential Hilbert Space

Qubits built out of Majorana zero modes constitute the primary path toward topologically protected quantum computing. Simulating the braiding process of multiple Majorana zero modes corresponds to the quantum dynamics of a superconducting many-body system. It is crucial to study the Majorana dynamics both in the presence of all other quasiparticles and for reasonably large system sizes. We present a method to calculate arbitrary many-body wave functions as well as their expectation values, correlators, and overlaps from time evolved single-particle states of a superconductor, allowing for significantly larger system sizes. We calculate the fidelity, transition probabilities, and joint parities of Majorana pairs to track the quality of the braiding process. Here, we show how the braiding success depends on the speed of the braid. Moreover, we demonstrate the topological CNOT two-qubit gate as an example of two-qubit entanglement. Our Letter opens the path to test and analyze the many theoretical implementations of Majorana qubits. Moreover, this method can be used to study the dynamics of any noninteracting superconductor.

Majorana bound states↗

Quantum Criticality Under Imperfect Teleportation

Entanglement, measurement, and classical communication together enable teleportation of quantum states between distant parties, in principle, with perfect fidelity. To what extent do correlations and entanglement of a many-body wave function transfer under teleportation protocols? We address this question for the case of an imperfectly teleported quantum critical wave function, focusing on the ground state of a critical Ising chain. We demonstrate that imperfections, e.g., in the entangling gate adopted for a given protocol, effectively manifest as weak measurements acting on the otherwise pristinely teleported critical state. Armed with this perspective, we leverage and further develop the theory of measurement-altered quantum criticality to quantify the resilience of critical-state teleportation. We identify classes of teleportation protocols for which imperfection (i) preserves both the universal long-range entanglement and correlations of the original quantum critical state, (ii) weakly modifies these quantities away from their universal values, and (iii) obliterates long-range entanglement altogether while preserving power-law correlations, albeit with a new set of exponents. We also show that mixed states describing the average over a series of sequential imperfect teleportation events retain pristine power-law correlations due to a “built-in” decoding algorithm, though their entanglement structure measured by the negativity depends on errors similarly to individual protocol runs. These results may allow one to design teleportation protocols that optimize against errors—highlighting a potential practical application of measurement-altered criticality. Published by the American Physical Society 2024

Physics↗

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗

UMass 2-Body WEC Techno-Economic Assessment

The University of Massachusetts (UMass) is developing a 2-body wave energy converter (WEC) device that is converting mechanical power into electricity using a mechanical motion rectifier that allows the system to couple to a flywheel. UMass has completed numerical modeling, wave tank testing, and PTO sub-system testing and needed assistance in developing a techno-economic model to enable optimization of their topology, comparison to a generic heaving point absorber topology, and guide the next steps in their development efforts. The core objective was to develop a techno-economic approach and modeling tool that allows benchmarking of the two topologies across a wide range of scales to evaluate their respective competitiveness in different application spaces. This data includes the final report as well as a supporting spreadsheet containing the data produced for this report.

16 TIDAL AND WAVE POWER↗

On Magnitudes and Inferred Yields of the 39 Underground Nuclear Test Explosions at the Novaya Zemlya Test Site

ABSTRACT We review the previously compiled seismic magnitudes for underground nuclear explosions at the Union of the Soviet Socialist Republics Novaya Zemlya Test Site, and analyze these data in the context of published details about individual tests and official Russian reports of annual cumulative explosive energy release at the test site. We derive relationships between teleseismic body wave magnitude mb and explosive yield Y by predicting the yield of each test using mb and minimizing the difference between predicted and reported annual cumulative yield. A linear relationship between mb and logY with a slope of 0.82 provides a good fit to the annual total yield, with a root mean square (rms) misfit of 0.095 logY units. A spectral magnitude relationship based on the explosion source model of Mueller and Murphy (1971) with a single-free-parameter provides a slightly better fit. For the 10 tests for which individual yields can be directly deduced from the official reports, the rms misfit of reported yields is less than 0.1 logY units or about 25%. More complex relationships that differentiate between the northern and southern test areas, or take account of variations in the scaled depth of burial, do not provide substantially better fits to the reported yields.

Geochemistry & Geophysics↗

Nucleon momentum distributions for local chiral interactions [Slides]

We use quantum Monte Carlo methods to calculate single- and two-nucleon momentum distributions in 4 He, 12 C, and 16 O. We use correlated many-body wave functions optimized for local chiral interactions up to next-to-next-to-leading order (N 2 L0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗

Neural-network quantum states for the nuclear many-body problem

A long-standing goal of nuclear theory is to explain how the structure and dynamics of atomic nuclei and neutron-star matter emerge from the underlying interactions among protons and neutrons. Achieving this goal requires solving the nuclear quantum many-body problem with high accuracy across a wide range of length scales and density regimes. In this review, we discuss how artificial neural network representations of the nuclear many-body wave function have significantly extended the capabilities of continuum quantum Monte Carlo methods. In particular, neural network quantum states enable calculations of larger systems than were previously accessible and provide a flexible framework for capturing phenomena that challenge conventional approaches, including the emergence of nuclear clusters and superfluid phases in dense matter. We highlight recent applications to finite nuclei, infinite nuclear and neutron matter, and dynamical processes relevant to lepton-nucleus and nucleus-nucleus scattering. We also discuss conceptual and methodological connections with condensed matter physics, emphasizing developments in neural network quantum states that bridge strongly correlated systems across disciplines. Together, these developments demonstrate how neural-network methods open new avenues toward unified and accurate descriptions of nuclear structure, matter, and reactions.

Lovato, Alessandro [Argonne; TIFPA-INFN, Trento; V↗

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗

RidgeAlloy: A high pressure die casting alloy that captures the approaching wave of automotive body sheet scrap

Oak Ridge National Laboratory (ORNL) has recently launched an effort to economically expand the North American aluminum supply chain by designing and developing a new family of Al-Mg-Si-Fe-Mn structural die cast aluminum alloys that can be made from up to 100% mixed 5xxx and 6xxx series automotive post-consumer sheet scrap and do not require heat treatment. The goal is to develop an alloy capable of capturing the approaching scrap wave of aluminum sheet intensive vehicles (in the 2030’s) into integrated structural high pressure die castings, rather than down-cycling these high-grade sheet alloys into non-structural castings. ORNL discovered a significant gap in the commercial thermodynamic databases that were unable to predict the trends of a key primary intermetallic phase as a function of composition. Experiments were performed to correct this gap and to create a unique, corrected quaternary database for Al-Mg-Si-Fe. High throughput computational methods, including solidification simulations, were used to identify a composition volume capable of avoiding this embrittling primary intermetallic phase, even at higher Fe + Si contents. Laboratory scale castings validated the thermodynamic and solidification predictions and resulted in lab-cast alloys with little primary intermetallic, and which met or exceeded structural casting properties requirements for yield strength and ductility. An accelerated demonstration of an HPDC automotive part was conducted in collaboration with two U.S. small businesses. Ingots made from 100% recycled 5xxx and 6xxx series scrap (plus Fe) were used to die cast one variant of this new family of Al-Mg-Si-Fe-Mn alloys (with high Si + Fe content) into a mid-sized HPDC structural-type automotive component. The HPDC alloy demonstration part showed good properties (especially ductility) and good castability for a complex component.

Plotkowski, Alex [ORNL] (ORCID:0000000154718681)↗

Charge ordering in Ir dimers in the ground state of Ba 5 AlIr 2 O 11

It has been well established experimentally that the interplay of electronic correlations and spin-orbit interactions in Ir 4+ and Ir 5+ oxides results in insulating J eff = 1/2 and J eff = 0 ground states, respectively. However, in compounds where the structural dimerization of iridium ions is favorable, the direct Ir d-d hybridization can be significant and takes a key role. Here, in this study, we investigate the effects of direct Ir d-d hybridization in comparison with electronic correlations and spin-orbit coupling in Ba 5 AlIr 2 O 11 , a compound with Ir dimers. Using a combination of ab initio many-body wave-function quantum chemistry calculations and resonant inelastic x-ray scattering experiments, we elucidate the electronic structure of Ba 5 AlIr 2 O 11 . We find excellent agreement between the calculated and the measured spin-orbit excitations. Contrary to expectations, the analysis of the many-body wave function shows that the two Ir (Ir 4+ and Ir 5+ ) ions in the Ir 2 O 9 dimer unit in this compound preserve their local J eff character close to 1/2 and 0, respectively. The local point group symmetry at each of the Ir ions plays an important role, significantly limiting the direct d-d hybridization. Our results emphasize that minute details in the local crystal field environment can lead to dramatic differences in the electronic states in iridates and 5d oxides in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗