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 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗
NASA NTRS2020

Computationally-efficient Structural Models for Analysis of Woven Composites

The paper presents a novel approach to model woven composite using the computationally efficient one-dimensional models. The framework is built within the scheme of the Carrera Unified Formulation (CUF), a generalized hierarchical formulation that generates variable kinematic structural theories. Various components of the woven composite unit cell are modeled using a combination of straight and curved one-dimensional CUF models. By employing a component-wise approach, a modeling technique within CUF, the complex geometry of the woven composite components is modeled precisely. The ability of CUF models to accurately resolve stress and strain fields are exploited to capture complex deformation within a woven composite unit cell. Numerical results include analyses of a non-crimped textile composite, a curved tow under tension, and a dry woven textile unit cell.

Kaleel, Ibrahim↗
DOE OSTI2022

Investigation of methods for targeted search of dominant higher modes in subcritical systems

The reactor physics modeling of subcritical systems requires a substantially different approach from that of a critical reactor due to the presence of the so-called higher modes. Earlier investigations showed that the dominant higher modes lie in the inner part of the eigenvalue spectrum, making the conventional eigenvalue searches unfeasible and calling for a targeted search. This paper investigates the possibility of targeted and multitudinous eigenvalue calculation by testing the Krylov-Schur and the Dynamic Modes Decomposition (DMD) method for simple, analytically tractable problems. While the DMD method clearly showed its efficiency in a one-group homogeneous problem, modal analysis of a three-dimensional reflected reactor in two-groups approximation demonstrated that the increasing complexity poses challenges for both methods. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Search for efficient formulations for Hamiltonian simulation of non-Abelian lattice gauge theories

Hamiltonian formulation of lattice gauge theories (LGTs) is the most natural framework for the purpose of quantum simulation, an area of research that is growing with advances in quantum-computing algorithms and hardware. It, therefore, remains an important task to identify the most accurate, while computationally economic, Hamiltonian formulation(s) in such theories, considering the necessary truncation imposed on the Hilbert space of gauge bosons with any finite computing resources. This paper is a first step toward addressing this question in the case of non-Abelian LGTs, which further require the imposition of non-Abelian Gauss’s laws on the Hilbert space, introducing additional computational complexity. Focusing on the case of SU(2) LGT in 1+1 dimensions coupled to one flavor of fermionic matter, a number of different formulations of the original Kogut-Susskind framework are analyzed with regard to the dependence of the dimension of the physical Hilbert space on boundary conditions, system’s size, and the cutoff on the excitations of gauge bosons. The impact of such dependencies on the accuracy of the spectrum and dynamics obtained from a Hamiltonian computation is examined, and the (classical) computational-resource requirements given these considerations are studied. Besides the well-known angular-momentum formulation of the theory, the cases of purely fermionic and purely bosonic formulations (with open boundary conditions), and the Loop-String-Hadron formulation are analyzed, along with a brief discussion of a Quantum-Link-Model formulation of the same theory. Clear advantages are found in working with the Loop-String-Hadron framework which implements non-Abelian Gauss’s laws a priori using a complete set of gauge-invariant operators. Although small lattices are studied in the numerical analysis of this work, and only the simplest algorithms are considered, a range of conclusions will be applicable to larger systems and potentially to higher dimensions. Future studies will extend this investigation to the analysis of resource requirements for quantum simulations of non-Abelian LGT, with the goal of shedding light on the most efficient Hamiltonian formulation of gauge theories of relevance in nature

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effects of Aluminum and Molybdenum on the Phase Stability of Iron-Chromium Alloys: A First-Principles Study

The interaction between solute atoms is critical to the thermodynamic behavior of Fe-Cr alloys, but the effects of non-dilute Al and Mo on the Fe-Cr phase stability and vacancy formation energy are not clearly understood. In this study, density functional theory, cluster expansion, and Monte Carlo simulation are used to predict the effects of ternary solute elements on the thermodynamic properties in multicomponent Fe-Cr alloys. The machine learning regression approach is applied to train and construct energy models that accurately describe the complex chemical interactions. The computational outcomes include the prediction of the partial ternary phase diagram, mixing enthalpy, and vacancy formation energy for different compositions. The phase boundary calculation predicts a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from a’ precipitates. The mixing enthalpy calculation shows that Mo may also reduce the Cr solubility in bcc Fe. Additionally, the simulations show that non-dilute Cr decreases the vacancy formation energy in bcc Fe, while adding Al results in a less significant effect. Finally, the results demonstrate important applications of using machine learning energy models to study model or commercial alloys with multicomponent solute species and point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗
DOE OSTI2022

Reaction processes at step edges on S-decorated Cu(111) and Ag(111) surfaces: MD analysis utilizing machine learning derived potentials

We report a variety of complexation, reconstruction, and sulfide formation processes can occur at step edges on the {111} surfaces of coinage metals (M) in the presence of adsorbed S under ultra-high vacuum conditions. Given the cooperative many-atom nature of these reaction processes, Molecular Dynamics (MD) simulation of the associated dynamics is instructive. However, only quite restricted Density Functional Theory (DFT)-level ab initio MD is viable. Thus, for M = Ag and Cu, we instead utilize the DeePMD framework to develop machine-learning derived potentials, retaining near-DFT accuracy for the M–S systems, which should have broad applicability. These potentials are validated by comparison with DFT predictions for various key quantities related to the energetics of S on M(111) surfaces. The potentials are then utilized to perform extensive MD simulations elucidating the above diverse restructuring and reaction processes at step edges. Key observations from MD simulations include the formation of small metal–sulfur complexes, especially MS 2 ; development of a local reconstruction at A-steps featuring an S-decorated {100} motif; and 3D sulfide formation. Additional analysis yields further information on the kinetics for metal–sulfur complex formation, where these complexes can strongly enhance surface mass transport, and on the propensity for sulfide formation.

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

Evidence of antiferromagnetism in ultrathin metallic (111)-oriented LaNiO 3 films

Antiferromagnets with exotic spin textures promise low-power spintronic devices with extremely high operating frequencies and resistance to external perturbations. In particular, the combination of highly tunable correlated electron physics, as in complex oxides, with metallicity and antiferromagnetism is desirable but exceedingly rare. LaNiO 3 , the lone example of a perovskite nickelate which is metallic across all temperatures, has long been a promising candidate, but the antiferromagnetic metallic state has remained elusive. Here, we demonstrate the emergence of this state in ultrathin films of (111)-oriented LaNiO 3 using a combination of polarized neutron reflectometry, low-energy muon spectroscopy and anomalous Hall effect measurements, and density functional theory calculations with a Hubbard U term. Further, we find a highly strained symmetry-breaking interfacial region which may support canting of the AFM moments leading to a vanishingly small net magnetization at the film-substrate interface, providing a convenient route toward control of the Néel order. Evidence of antiferromagnetic metallic behavior in (111)LNO films highlights the role of crystal symmetry in tuning the novel quantum states in complex oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comments on the quantum field theory of the Coulomb gas formalism

The holomorphic Coulomb gas formalism, as developed by Feigin-Fuchs, Dotsenko-Fateev and Felder, is a set of rules for computing minimal model observables using free field techniques. We attempt to derive and clarify these rules using standard techniques of quantum field theory. We begin with a careful examination of the timelike linear dilaton. Although the background charge of the model breaks the scalar field’s continuous shift symmetry, the exponential of the action remains invariant under a discrete shift because the background charge is imaginary. Gauging this symmetry makes the dilaton compact and introduces winding modes into the spectrum. One of these winding operators corresponds to the anti-holomorphic completion of the BRST current first introduced by Felder, and the full left/right cohomology of this BRST charge isolates the irreducible representations of the Virasoro algebra within the degenerate Fock space of the linear dilaton. The “supertrace” in the BRST complex reproduces the minimal model partition function and exhibits delicate cancellations between states with both momentum and winding. The model at the radius $R=\sqrt{pp^{\prime }}$ has two marginal operators corresponding to the Dotsenko-Fateev “screening charges”. Deforming by them, we obtain a model that might be called a “BRST quotiented compact timelike Liouville theory”. The Hamiltonian of the zero-mode quantum mechanics of this model is not Hermitian, but it is PT -symmetric and exactly solvable. Its eigenfunctions have support on an infinite number of plane waves, suggesting an infinite reduction in the number of independent states in the full quantum field theory. Applying conformal perturbation theory to the exponential interactions reproduces the Coulomb gas calculations of minimal model correlation functions. In contrast to spacelike Liouville, these “resonance correlators” are finite because the zero mode is compact. We comment on subtleties regarding the reflection operator identification, as well as naive violations of truncation in correlators with multiple reflection operators inserted. This work is part of an attempt to understand the relationship between the JT model of two dimen- sional gravity and the worldsheet description of the (2 , p ) minimal string as suggested by Seiberg and Stanford.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An Adsorptive Membrane Platform for Precision Ion Separation: Membrane Design and First‐Principles Studies

One of the key challenges in separation science is the lack of precise ion separation methods and mechanistic understanding crucial for efficiently recovering critical materials from complex aqueous matrices. Herein, first-principles electronic structure calculations and in situ Raman spectroscopy are studied to elucidate the factors governing ion discrimination in an adsorptive membrane specifically designed for transition metal ion separation. Density functional theory calculations and in situ Raman data jointly reveal the thermodynamically favorable binding preferences and detailed adsorption mechanisms for competing ions. How membrane binding preferences correlate with the electronic properties of ligands is explored, such as orbital hybridization and electron localization. The findings underscore the importance of the phenolate group in oxime ligands for achieving high selectivity among competing transition metal ions. In-depth understanding on which specific atomistic site within the microenvironment of metal-ligand binding pockets governs the ion discrimination behaviors of the host will build a solid foundation to guide the rational design of next-generation materials for precision separation essential for energy technologies and environment remediation. In tandem, synthetic controllability is demonstrated to transform 3D micrometer-scale crystals to a 2D crystalline selective layer in membranes, paving the way for more precise and sustainable advances in separation science.

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

One-loop central-emission vertex for two gluons in $\mathcal{N}$ = 4 super Yang-Mills theory

A necessary ingredient for extending the BFKL equation to next-to-next-to-leading logarithmic (NNLL) accuracy is the one-loop central emission vertex (CEV) for two gluons which are not strongly ordered in rapidity. Here we consider the one-loop six-gluon amplitude in $\mathcal{N}$ = 4 super Yang-Mills (SYM) theory in a central next-to-multi-Regge kinematic (NMRK) limit, we show that its dispersive part factorises in terms of the two-gluon CEV, and we use it to extract the one-loop two-gluon CEV for any helicity configuration within this theory. This is a component of the two-gluon CEV in QCD. Although computed in the NMRK limit, both the colour structure and the kinematic dependence of the two-gluon CEV capture much of the complexity of the six-gluon amplitudes in general kinematics. In fact, the transcendental functions of the latter can be conveniently written in terms of impact factors, trajectories, single-emission CEVs and a remainder, which is a function of the conformally invariant cross ratios which characterise the six-gluon amplitudes in planar $\mathcal{N}$ = 4 SYM. Finally, as expected, in the MRK limit the two-gluon CEV neatly factorises in terms of two single-emission CEVs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computational Design of Metal-Exchanged MFI Zeolites for Coupling CO 2 –Ethylene to Form Acrylic Acid

The production of acrylic acid by coupling CO 2 and ethylene is highly desirable, as it uses a C 1 molecule to introduce the carboxylic group into a large commodity chemical, forming an unsaturated carboxylic acid. Molecular complexes have been used as catalysts in the literature for this reaction; however, these homogeneous studies were all limited by low activity and stability, and efficient heterogeneous catalysts remain to be explored. Here, we report density functional theory calculations, through which we investigated metal-exchanged MFI zeolites with metals from 3d, 4d, and 5d series for this reaction. Screening of reaction thermodynamics has shown that Sc, Y, Cr, Mo, W, Mn, Tc, and Re dispersed into MFI zeolites are potential catalysts. We further calculated the activation barriers for the β–H transfer, a key elementary step that normally determines the overall rates, and found the lowest barrier over Sc- and Y-MFI; such a feasible kinetics is attributed to the low electronegativity of Sc and Y that strongly affects the nature of the transition states and the charge redistribution during the β–H transfer. In conclusion, this work thus provides guidelines for designing a zeolite-based catalyst for this valuable reaction and identifies promising metal centers to be used.

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

High-Energy Hybridized States Enable Long-Lived Hot Electrons in Cobaloxime-Silicon Nanocrystal System

Strong electronic coupling is achieved between the molecular catalyst cobaloxime ([Co]) and silicon nanocrystals (Si NCs) bridged by an ethylenepyridine group derived from vinylpyridine (vpy) covalently bound to the Si NC surface (Si-vpy-[Co]). The ethylenepyridine tether in Si-vpy-[Co] is key to dramatic changes to the system’s physical properties which are not observed in the corresponding formylpyridine (fpy) system (Si-fpy-[Co]) consistent with strong electronic coupling previously observed only in dark electrochemical systems. UV−vis absorption spectroscopy reveals new [Co]-centered electronic states in Si-vpy-[Co], and transient absorption spectroscopy finds a strong absorption feature appearing within 250 fs and persisting for at least 5 ns. Astoundingly, spectroelectrochemical measurements reveal that this absorption feature is consistent with both the singly reduced [Co] − and doubly reduced [Co] 2− complexes, leading to the conclusion that these long-lived charges are derived from high-energy “hot” electrons residing in [Co]-centered states. Detailed analysis using cyclic voltammetry, spectroelectrochemistry, electron paramagnetic resonance spectroscopy, and density functional theory (DFT) calculations provides insight into the unique electronic structure created in Si-vpy-[Co]. DFT reveals that the new electronic states arise from hybridization between deep Si NC band states and high-energy molecular orbitals of the ethylenepyridine tether and the [Co] catalyst and are facilitated by σ-bonding character at the ethylenepyridine linkage. This study demonstrates that strong electronic coupling achieved through precise molecular chemistry can change the paradigm of otherwise fixed energy levels in hybrid photoelectrochemical systems for artificial photosynthesis and related applications.

14 SOLAR ENERGY↗
DOE OSTI2022

Resolving atomistic structure and oxygen evolution activity in nickel antimonates

The oxygen evolution reaction (OER) requires electrodes that are not only catalytically active, but also stable under harsh electrochemical environments to enable efficient, durable technologies. Our recent report of a stable amorphous Ni 0.5 Sb 0.5 O z OER photoanode established Ni–Sb–O as an important system for computational understanding of both the structural and catalytic behavior of these complex oxides. In the present work we show that Ni x Sb 1-x O z with x > 0.33 crystallizes into a previously unknown phase. Guided by experimental X-ray diffraction, here, we use density functional theory calculations to perform a prototype phase search to identify a broad family of stable and metastable mixed rutile and hexagonal-like phases for x = 0.33, 0.50, and 0.66 compositions. For the identified phases, we predict favorable oxygen vacancy formation energies for Ni-rich compositions under the reducing synthesis conditions which match measured Ni K-edge X-ray absorption spectra. The calculated overpotential for the most active site decreases with increasing Ni content, from 0.91 V (x = 0.33) to 0.49 V (x = 0.66), which captures the experimentally observed trend. We find the active site changes from the Ni–O–Sb bridge to a Ni–O–Ni bridge at increasing Ni concentrations, rather than the commonly studied singly under-coordinated sites. Finally, detailed Pourbaix analysis of the identified phases show excellent electrochemical stability, consistent with experimentally measured low metal ion concentrations in the electrolyte of photoelectrochemical cells. Collectively, our consideration of an ensemble of structures enables identification of the most catalytically prolific structural motifs, aiding the understanding of crystalline and amorphous catalysts and elucidating the co-optimization of activity and durability in nickel antimonates.

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

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗