Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity theory”

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 OSTI2023

Understanding the effects of fabrication process on BaZr 0.9 Y 0.1 O 3– δ grain-boundary chemistry using atom probe tomography

Proton-conducting oxides, such as BaZr 0.9 Y 0.1 O 3–δ (BZY10), are receiving significant scientific attention for application in a variety of electrochemical applications. It has been shown that the synthesis route has a significant effect on the ionic behavior of these materials, especially the grain boundaries (GB). In this study, laser-pulsed atom probe tomography (APT) is used to characterize the GB chemistry of BZY10 samples prepared by four different methods: (1) spark plasma sintering (SPS), (2) conventional sintering followed by high-temperature annealing (HT), (3) conventional sintering with NiO used as a sintering aid (SSR-Ni), and (4) solid-state reactive sintering of oxide precursor powders with NiO as a reactive sintering aid (SSRS-Ni). Oxygen depletion (corresponding to oxygen vacancy accumulation) was observed at every GB in this study, consistent with the positive space-charge GB region commonly assigned to these materials. In contrast to the consistent trends associated with oxygen depletion, cation constituents showed varying segregation/depletion behaviors. The samples that used NiO revealed Ni preferentially segregating to the GB. The most common impurities in BZY10 are Al, Fe, Mg, Si, and Sr, all of which generally accumulate at GBs. In order to quantify these chemical trends, the APT method was optimized for this material supplemented with understanding from density functional theory (DFT). Finally, this APT study reveals the complexity of BZY10 GB chemistry and shows that even adjacent GBs in the same material can sometimes show surprisingly different GB chemistry.

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

Highly transferable atomistic machine-learning potentials from curated and compact datasets across the periodic table

Machine learning atomistic potentials trained using density functional theory (DFT) datasets allow for the modeling of complex material properties with near-DFT accuracy while imposing a fraction of its computational cost. The curation of the DFT datasets can be extensive in size and time-consuming to train and refine. In this study, we focus on addressing these barriers by developing minimalistic and flexible datasets for many elements in the periodic table regardless of their mass, electronic configuration, and ground state lattice. These DFT datasets have on average, ~4000 different structures and 27 atoms per structure, which we found sufficient to maintain the predictive accuracy of DFT properties and notably with high transferability. We envision these highly curated training sets as starting points for the community to expand, modify, or use with other machine learning atomistic potential models, whatever may suit individual needs, further accelerating the utilization of machine learning as a tool for material design and discovery.

42 ENGINEERING↗
DOE OSTI2020

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

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

Downward quantum learning from element 118: Automated generation of Fermi–Löwdin orbitals for all atoms

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

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

Scaling functions for eigenstate entanglement crossovers in harmonic lattices

For quantum matter, eigenstate entanglement entropies obey an area law or log-area law at low energies and small subsystem sizes and cross over to volume laws for high energies and large subsystems. This transition is captured by crossover functions, which assume a universal scaling form in quantum critical regimes. We demonstrate this for the harmonic lattice model, which describes quantized lattice vibrations and is a regularization for free scalar field theories, modeling, e.g., spin-0 bosonic particles. In one dimension, the ground-state entanglement obeys a log-area law. For dimensions d ≥ 2, it displays area laws, even at criticality. The distribution of excited-state entanglement entropies is found to be sharply peaked around subsystem entropies of corresponding thermodynamic ensembles in accordance with the eigenstate thermalization hypothesis. Numerically, we determine crossover scaling functions for the quantum critical regime of the model and do a large-deviation analysis. We show how infrared singularities of the system can be handled and how to access the thermodynamic limit using a perturbative trick for the covariance matrix. Eigenstates for quasi-free bosonic systems are not Gaussian. Here, we resolve this problem by considering appropriate squeezed states instead. For these, entanglement entropies can be evaluated efficiently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simultaneous stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Furthermore, solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Accurate and efficient description of interacting carriers in quantum nanostructures by selected configuration interaction and perturbation theory

In this paper, we present a method to accurately and efficiently calculate many-body states of interacting carriers in quantum nanostructures based on a combination of iterative selection of configurations and perturbation theory. This method enables investigations of large excitonic complexes and multielectron systems with near full configuration interaction accuracy, even though only a small subspace of the full many-body Hilbert space is sampled, thus saving orders of magnitudes in computational resources. Important advantages of this method are that the convergence is controlled by a single parameter, the threshold, and that ground and excited states can be treated on an equal footing. On the example of InAsP nanowire quantum dots described using a million-atom tight-binding approach, we demonstrate the extreme efficiency of the method by numerical studies of large excitonic complexes filling up to the fourth electronic shell. We find that the method generally converges fast as a function of the threshold, profiting from a significant acceleration due to the perturbative corrections. The role of the choice of single-particle basis states is discussed. It is found that the algorithm converges faster in the Hartree-Fock basis only for highly charged systems, where Coulomb repulsion dominates. Finally, based on the observation that second-order perturbative energy corrections only depend on off-diagonal elements of the many-body Hamiltonian, we present a way to accurately calculate many-body states that requires only a relatively small number of Coulomb matrix elements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Orthorhombic distortion drives orbital ordering in the antiferromagnetic 3$d$ 1 Mott insulator PrTiO 3

The orbital, which represents the shape of the electron cloud, very often strongly influences the manifestation of various exotic phenomena, e.g., magnetism, metal-insulator transition, colossal magnetoresistance, unconventional superconductivity, etc. in solid-state systems. The observation of the antiferromagnetism in $RE$TiO 3 ($RE$ = rare-earth) series has been puzzling since the celebrated Kugel-Khomskii model of spin-orbital superexchange predicts ferromagnetism in an orbitally degenerate $d$ 1 system. Further, the existence of the orbitally ordered vs. orbital liquid phase in both antiferromagnetic and paramagnetic phase have been unsettled issues thus far. To address these longstanding questions, we investigate single crystalline film of PrTiO 3 . Our synchrotron x-ray diffraction measurements confirm the retention of bulklike orthorhombic ($D_{2h}$) symmetry in the thin film geometry. We observe similar x-ray linear dichroism signal in both paramagnetic and antiferromagnetic phase, which can be accounted by ferro-orbital ordering (FOO). While the presence of $D_{2h}$ crystal field does not guarantee lifting of orbital degeneracy always, we find it to be strong enough in these rare-earth titanates, leading to the FOO state. Thus, our work demonstrates the orthorhombic distortion is the driving force for the orbital ordering of antiferromagnetic $RE$TiO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structural and electronic properties of the first iridium containing mixed B-site spinel oxide: Cu [ Ir 1.5 Cu 0.5 ] O 4

Geometrically frustrated systems populated with large spin-orbit coupled ions are an ideal setting for the exploration of novel exotic states of matter. In this study we present an example of iridium on a mixed B-site spinel oxide structure: Cu [ Ir 1.498 ( 2 ) Cu 0.502 ( 2 ) ] O 4 . Synchrotron XRD refinements reveal a face-centered-cubic structure with space group F d 3 ¯ m and mixed Cu-Ir site disorder within the B 2 O 4 rocksalt substructure. Electrical properties reveal a metallic state within the 50–600-K range with a Kondo effect at T < 50 K . X-ray absorption spectroscopy (XAS) measurements show a mixed Cu 1 + / 2 + and Ir 3 + / 4 + charge partitioned picture, which suggests a metallic/band description with reduced on-site Coulomb interactions. Spin-glass-like freezing is seen at T g = 49 K , and the hysteresis behavior for T > T g resembles that of a strongly frustrated magnet. DFT calculations show sizable hybridization between the Cu 3 d and Ir 5 d states with an effective mixed Ir 3 + / 4 + charge partitioned picture, supporting the electronic and XAS results.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying uncertainty in machine learning for nuclear binding energy

Techniques from artificial intelligence and machine learning are increasingly employed in nuclear theory; however, the uncertainties that arise from the complex parameter manifold encoded by the neural networks are often overlooked. Epistemic uncertainties arising from training the same network multiple times for an ensemble of initial weight sets offer a first insight into the confidence of machine learning predictions, but they often come with a high computational cost. Instead, we apply a single-model uncertainty quantification method called Δ-UQ that gives epistemic uncertainties with one-time training. Here, we demonstrate our approach on a two-feature model of nuclear binding energies per nucleon with proton and neutron number pairs as inputs. We show that Δ-UQ can produce reliable and self-consistent epistemic uncertainty estimates and can be used to assess the degree of confidence in predictions made with deep neural networks.

Huang, Mengyao [Lawrence Livermore National Labora↗
DOE OSTI2022

Leaf water potential measurements using the pressure chamber: Synthetic testing of assumptions towards best practices for precision and accuracy

Leaf water potential (ψ leaf ), typically measured using the pressure chamber, is the most important metric of plant water status, providing high theoretical value and information content for multiple applications in quantifying critical physiological processes including drought responses. Pressure chamber measurements of ψ leaf (ψ leafPC ) are most typical, yet, the practical complexity of the technique and of the underlying theory has led to ambiguous understanding of the conditions to optimize measurements. Consequently, specific techniques and precautions diversified across the global research community, raising questions of reliability and repeatability. Here, we surveyed specific methods of ψ leafPC from multiple laboratories, and synthesized experiments testing common assumptions and practices in ψ leafPC for diverse species: (i) the need for equilibration of previously transpiring leaves; (ii) leaf storage before measurement; (iii) the equilibration of ψ leaf for leaves on bagged branches of a range of dehydration; (iv) the equilibration of ψ leaf across the lamina for bagged leaves, and the accuracy of measuring leaves with artificially ‘elongated petioles’; (v) the need in ψ leaf measurements for bagging leaves and high humidity within the chamber; (vi) the need to avoid liquid water on leaf surfaces; (vii) the use of ‘pulse’ pressurization versus gradual pressurization; and (viii) variation among experimenters in ψ leafPC determination. Based on our findings we provide a best practice protocol to maximise accuracy, and provide recommendations for ongoing species-specific tests of important assumptions in future studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spectral Equivalence of Low-Order Discretizations for High-Order H(curl) and H(div) Spaces

In this paper, we present spectral equivalence results for high-order tensor product edge- and face-based finite elements for the H(curl) and H(div) function spaces. Specifically, we show for certain choices of shape functions that the mass and stiffness matrices of the high-order elements are spectrally equivalent to those for an assembly of low-order elements on the associated Gauss--Lobatto--Legendre mesh. Based on this equivalence, efficient preconditioners can be designed with favorable computational complexity. Numerical results are presented which confirm the theory and demonstrate the benefits of the equivalence results for overlapping Schwarz preconditioners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spectral Equivalence Properties of Higher-Order Tensor Product Finite Elements and Applications to Preconditioning

In this study, we present spectral equivalence results for higher-order tensor product edge-, face- and interior-based finite elements. Specifically, we show for certain choices of shape functions that the mass and stiffness matrices of the higher-order elements are spectrally equivalent to those for an assembly of lowest-order elements on the associated Gauss-Lobatto-Legendre mesh. Based on this equivalence, efficient preconditioners can be designed with favorable computational complexity. Numerical results are presented which confirm the theory and demonstrate the benefits of the equivalence results for overlapping Schwarz preconditioners.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Approach for Inferring Full-Scope Human Reliability Data Based on Simplified Simulator Data

This paper proposes a method for inferring full-scope human reliability data based on the Simplified Human Error Experimental Program (SHEEP) data. It mainly focuses on the human errors observed when using simulators with different complexity levels. In the proposed method, the manner in which human error probabilities (HEPs) change as a result of increasing simulator complexity and how simulator complexity levels are quantified represent key information for inferring full-scope data. In the present study, SHEEP error data pertaining to actual professional operators using Rancor Microworld (Rancor) (i.e., a more simplified simulator) and Compact Nuclear Simulator (CNS) (i.e., a less simplified simulator) were compared with the HuREX error data. An approach to quantifying simulator complexity levels was then proposed based on information theory and acquired eye-tracker data.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2020

Application of Quantum Computing to Biochemical Systems: A Look to the Future

Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.

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

A comparison of Gaussian processes and neural networks for computer model emulation and calibration

The Department of Energy relies on complex physics simulations for prediction in domains like cosmology, nuclear theory, and materials science. These simulations are often extremely computationally intensive, with some requiring days or weeks for a single simulation. In order to assure their accuracy, these models are calibrated against observational data in order to estimate inputs and systematic biases. Because of their great computational complexity, this process typically requires the construction of an emulator, a fast approximation to the simulation. In this paper, two emulator approaches are compared: Gaussian process regression and neural networks. Their emulation accuracy and calibration performance on three real problems of Department of Energy interest is considered. On these problems, the Gaussian process emulator tends to be more accurate with narrower, but still well-calibrated uncertainty estimates. The neural network emulator is accurate, but tends to have large uncertainty on its predictions. Finally, as a result, calibration with the Gaussian process emulator produces more constrained posteriors that still perform well in prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Influence of Ga doping on magnetic properties, magnetocaloric effect, and electronic structure of pseudobinary GdZn 1 - x Ga x ( x = 0 - 0.1 )

Here, we explore the impact of introducing IIIA-group element Ga in place of IIB-group element Zn in binary intermetallic GdZn on its magnetic and magnetocaloric properties, as well as explicate the modified electronic band structure of the compound. The magnetic transition temperature of the compound decreases with the increase of Ga concentration in GdZn 1-x ⁢Ga x ⁢(x = 0–0.1) while the crystal structure (CsCl prototype) and lattice parameters remain unchanged. Our detailed analysis of magnetization and magnetocaloric data conclusively proves that long-ranged magnetic ordering exists in the sample, despite the magnetic interaction considerably weakening with the increase of Ga. The experimental data are rationalized using both a theoretical machine-learning model and first-principle density functional theory. The electronic band structure of GdZn is manifested with some unusual complex features which gradually diminish with Ga doping and the conventional sinusoidal feature of Ruderman-Kittel-Kasuya-Yosida-type interactions also disappears. A mean-field theory model is developed and can successfully describe the overall magnetocaloric behavior of the GdZn 1-x⁢ Ga x series of samples.

36 MATERIALS SCIENCE↗