Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution 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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fermionic mean-field theory as a tool for studying spin Hamiltonians

The Jordan–Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Furthermore, Jordan–Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J 1 –J 2 Heisenberg models, as well as to the pairing or reduced Bardeen–Cooper–Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Spiral Galaxies Flat Rotational Velocity Curve Explained by the Constant Group Velocity of a Nonlinear Density Wave

The rotation velocity curves of stars in galaxies, the motions of pairs of galaxies, and the behavior of galaxies in clusters and super-clusters all indicate that there is a lack of mass on different scales in the universe. In this paper, we derive the expression for rotational velocity using the nonlinear density wave theory considering only stellar components and we show that such theory can support the observed flat rotational velocity curve due to the main property of the soliton wave, which is a constant group velocity of the wave. The surface mass density (SMD) function, used to derive gravitational potential gradient and rotational velocity, is not assumed but rather derived as a solution of the nonlinear Srödinger equation, on the contrary to the widely used, in the literature, exponential disk approximation. Three parameters relevant to the curve shape are the intensities of equilibrium SMD, the amplitude of the wave, and total angular velocity or differential rotation, equivalently. Since the shape of the rotational velocity is highly sensitive to the mentioned parameters, this theory eventually provides a method for a very accurate estimation of galaxy mass and angular velocity as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Selective oxidation and nickel enrichment hinders the repassivation kinetics of multi-principal element alloy surfaces

Robust and sustained corrosion resistance in multi-principal element alloys (MPEA) requires rapid repassivation, i.e. regrowth of the passive layer once it is damaged or destroyed at the surface. In this study, we show that the repassivation of Al 0.1 CrCoFeNi MPEA in 0.6 M NaCl solution is hindered at pH of ~ 2.4 - 6.8 due to the formation of a Ni-enriched subsurface layer as a result of selective oxidation and dissolution of several principal elements, which can be fully restored at pH of ~ 14 from the oxidation of all principal elements. Specifically, surface characterization via X-ray photoelectron spectroscopy (XPS), high-angle annular dark-field (HAADF) imaging in scanning transmission electron microscopy (STEM), and atom probe tomography (APT), are coupled with density functional theory (DFT) calculations to determine surface composition, oxidation state, and electron work function to uncover the structural origin of the pH-dependent repassivation mechanisms. It was found that selective oxidation of Cr, Co, and Fe in the acidic to neutral solutions altered the surface composition to be significantly enriched in Ni as compared to the bulk. Once the original passive film is destroyed locally by either pitting or tribocorrosion, this altered surface composition exhibited a much poorer repassivation capability due to the increased electron work function and reduced surface reactivity at higher Ni concentration. Finally, these understandings could shed light on the future compositional design of non-equiatomic MPEAs towards sustained repassivation and corrosion resistance over a wide pH range.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Designing order–disorder transformation in high-entropy ferritic steels

Order–disorder transformations hold an essential place in chemically complex high-entropy ferritic steels (HEFSs) due to their critical technological application. The chemical inhomogeneity arising from mixing of multi-principal elements of varying chemistry can drive property altering changes at the atomic scale, in particular short-range order. Using density-functional theory-based linear-response theory, we predict the effect of compositional tuning on the order–disorder transformation in ferritic steels—focusing on Cr–Ni–Al–Ti–Fe HEFSs. We show that Ti content in Cr–Ni–Al–Ti–Fe solid solutions can be tuned to modify short-range order that changes the order–disorder path from BCC-B2 (Ti atomic-fraction = 0) to BCC-B2-L2 1 (Ti atomic-fraction > 0) consistent with existing experiments. Our study suggests that tuning degree of SRO through compositional variation can be used as an effective means to optimize phase selection in technologically useful alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

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

Effects of Sn addition on precipitation of Al–Zn–Mg alloy at early stage of natural aging

The effect of Sn addition on the natural aging behavior of Al–Zn–Mg alloy was investigated using atom probe tomography (APT) experiments and density functional theory (DFT) calculations. The results demonstrate that Sn retards the natural aging behavior by facilitating the preferential formation of Mg–Sn clusters, which suppress the formation of Guinier-Preston zones (I) during the earliest aging stage. Quantitative solute analysis conducted by APT has suggested increasing amount of Mg–Sn clusters in Sn-contained alloy due to their high binding energy, supported by DFT calculations on the Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Use of FLOSIC for understanding anion-solvent interactions

An Achille’s heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3−, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth’s atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi–Löwdin–Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree–Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Chemistry↗
DOE OSTI2021

Lie group analysis of the point-reactor neutron kinetics equations for various reactivity models

This work applies Lie Group Theory to the nonlinear point-reactor neutron kinetics equations for several reactivity feedback and reactivity insertion models. The goal of this work is to connect known analytical solutions to their corresponding Lie groups as well as to obtain new analytical solutions and the circumstances from which they arise. The reactivity models we consider are the Nordheim-Fuchs Model, an arbitrary-order polynomial reactivity insertion model, and the Fuchs Ramp-Insertion Model. In all three models, we generalize the analysis when possible and obtain analytical solutions in every case. For the first two models, we succeed in finding the Lie groups corresponding to classical solutions and find additional solutions that were otherwise unknown. Finally, for the final model, we obtain a reduced-order nonlinear kinetics equation which is not analytically solvable, but by applying a second-order truncation, we find a novel closed-form solution and demonstrate its validity with a numerical example.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

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

Isotropic-to-Nematic Transition in Salt-Free Polyelectrolyte Coacervates from Coarse-Grained Simulations

Recent interest in complex coacervation between oppositely charged polyelectrolytes (PEs) has been fueled by its relevance to biology in the context of membraneless organelle formation within living cells. For PEs with limited flexibility (such as double-stranded DNA), theoretical treatments and recent experiments have reported the emergence of liquid crystalline order (LCO) within the resulting coacervate phases. In this work, we study the underlying physics of this phenomenon using coarse-grained molecular dynamics simulations of symmetric semiflexible-semiflexible and asymmetric semiflexible-flexible coacervates. By comparing coacervates with the corresponding semidilute solutions of neutral polymers, we demonstrate that the presence of Coulomb interactions in coacervates facilitates orientational ordering, in agreement with theoretical predictions. Quantitative comparisons between our simulations and theory indicate that, for asymmetric nematic coacervates, the strong orientational ordering of stiff polyanions induces a weak ordering of the flexible polycations-an effect that was not anticipated by available theoretical studies. Simulations reveal that, for nematic coacervates, the preferred orientation of the PE chains at the liquid-liquid coacervate-supernatant interface is parallel, and the alignment of semiflexible PEs is homogeneous. The results presented here provide new molecular-level insights into the intracoacervate LCO and will help motivate further experimental and theoretical activities in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembly of Metallo-Supramolecules under Kinetic or Thermodynamic Control: Characterization of Positional Isomers Using Scanning Tunneling Spectroscopy

Coordination-driven self-assembly has been extensively employed to construct a variety of discrete structures as a bottom-up strategy. However, mechanistic understanding regarding whether self-assembly is under kinetic or thermodynamic control is less explored. To date, such mechanistic investigation has been limited to distinct, assembled structures. It still remains a formidable challenge to study the kinetic and thermodynamic behavior of self-assembly systems with multiple assembled isomers due to the lack of characterization methods. Herein, we use a stepwise strategy which combined self-recognition and self-assembly processes to construct giant metallo-supramolecules with 8 positional isomers in solution. With the help of ultrahigh-vacuum, low-temperature scanning tunneling microscopy and scanning tunneling spectroscopy, we were able to unambiguously differentiate 14 isomers on the substrate which correspond to 8 isomers in solution. Through measurement of 162 structures, the experimental probability of each isomer was obtained and compared with the theoretical probability. Finally, such a comparison along with density functional theory (DFT) calculation suggested that although both kinetic and thermodynamic control existed in this self-assembly, the increased experimental probabilities of isomers compared to theoretical probabilities should be attributed to thermodynamic control.

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

Uranium-Mediated Peroxide Activation and a Precursor toward an Elusive Uranium cis -Dioxo Fleeting Intermediate

The activation of chalcogen–chalcogen bonds using organometallic uranium complexes has been well documented for S–S, Se–Se, and Te–Te bonds. In stark contrast, reports concerning the ability of a uranium complex to activate the O–O bond of an organic peroxide are exceedingly rare. Herein, we describe the peroxide O–O bond cleavage of 9,10-diphenylanthracene-9,10-endoperoxide in nonaqueous media, mediated by a uranium(III) precursor [(( Me,Ad ArO) 3 N)U III (dme)] to generate a stable uranium(V) bis-alkoxide complex, namely, [(( Me,Ad ArO) 3 N)U V (DPAP)]. This reaction proceeds via an isolable, alkoxide-bridged diuranium(IV/IV) species, implying that the oxidative addition occurs in two sequential, single-electron oxidations of the metal center, including rebound of a terminal oxygen radical. Furthermore, this uranium(V) bis-alkoxide can then be reduced with KC 8 to form a uranium(IV) complex, which upon exposure to UV light, in solution, releases 9,10-diphenylanthracene to generate a cyclic uranyl trimer through formal two-electron photooxidation. Analysis of the mechanism of this photochemical oxidation via density functional theory (DFT) calculations indicates that the formation of this uranyl trimer occurs through a fleeting uranium cis-dioxo intermediate. At room temperature, this cis-configured dioxo species rapidly isomerizes to a more stable trans configuration through the release of one of the alkoxide ligands from the complex, which then goes on to form the isolated uranyl trimer complex.

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

Chemical crystallography by serial femtosecond X-ray diffraction

Abstract Inorganic–organic hybrid materials represent a large share of newly reported structures, owing to their simple synthetic routes and customizable properties 1 . This proliferation has led to a characterization bottleneck: many hybrid materials are obligate microcrystals with low symmetry and severe radiation sensitivity, interfering with the standard techniques of single-crystal X-ray diffraction 2,3 and electron microdiffraction 4–11 . Here we demonstrate small-molecule serial femtosecond X-ray crystallography (smSFX) for the determination of material crystal structures from microcrystals. We subjected microcrystalline suspensions to X-ray free-electron laser radiation 12,13 and obtained thousands of randomly oriented diffraction patterns. We determined unit cells by aggregating spot-finding results into high-resolution powder diffractograms. After indexing the sparse serial patterns by a graph theory approach 14 , the resulting datasets can be solved and refined using standard tools for single-crystal diffraction data 15–17 . We describe the ab initio structure solutions of mithrene (AgSePh) 18–20 , thiorene (AgSPh) and tethrene (AgTePh), of which the latter two were previously unknown structures. In thiorene, we identify a geometric change in the silver–silver bonding network that is linked to its divergent optoelectronic properties 20 . We demonstrate that smSFX can be applied as a general technique for structure determination of beam-sensitive microcrystalline materials at near-ambient temperature and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Wetting behavior of polyelectrolyte complex coacervates on solid surfaces

The wetting behavior of complex coacervates underpins their use in many emerging applications of surface science, particularly wet adhesives and coatings. Many factors dictate if a coacervate phase will condense on a solid surface, including solution conditions, the nature of the polymer–substrate interaction, and the underlying supernatant–coacervate bulk phase behavior. In this work, we use a simple inhomogeneous mean-field theory to study the wetting behavior of complex coacervates on solid surfaces both off-coexistence (wetting transitions) and on-coexistence (contact angles). Here we focus on the effects of salt concentration, the polycation/polyanion surface affinity, and the applied electrostatic potential on the wettability. We find that the coacervate generally wets the surface via a first order wetting transition with second order transitions possible above a surface critical point. Applying an electrostatic potential to a solid surface always improves the surface wettability when the polycation/polyanion–substrate interaction is symmetric. For asymmetric surface affinity, the wettability has a nonmonotonic dependence with the applied potential. We use simple scaling and thermodynamic arguments to explain our results.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Nuclear effective field theory: Status and perspectives

The nuclear physics landscape has been redesigned as a sequence of effective field theories (EFTs) connected to the Standard Model through symmetries and lattice simulations of Quantum Chromodynamics (QCD). EFTs in this sequence are expansions around different low-energy limits of QCD, each with its own characteristics, scales, and ranges of applicability regarding energy and number of nucleons. We review each of the three main nuclear EFTs---Chiral, Pionless, Halo/Cluster---highlighting their similarities, differences, and connections. Here, in doing so, we survey the structural properties and reactions of nuclei that have been derived from the \abinitio solution of the few- and many-body problem built upon EFT input.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Solvable theory of a strange metal at the breakdown of a heavy Fermi liquid

In this work, we introduce an effective theory for quantum critical points (QCPs) in heavy-fermion systems, involving a change in carrier density without symmetry breaking. Our theory captures a strongly coupled metallic QCP, leading to robust marginal Fermi-liquid transport phenomenology, and associated linear in temperature (T) "strange metal" resistivity, all within a controlled large-N limit. In the parameter regime of strong damping of emergent bosonic excitations, the QCP also displays a near-universal "Planckian" transport lifetime τ tr ~ℏ/(k B T). This is contrasted with the conventional so-called "slave boson" theory of the Kondo breakdown, where the large-N limit describes a weak coupling fixed point and nontrivial transport behavior may only be obtained through uncontrolled 1/N corrections. We also compute the weak-field Hall coefficient within the effective model as the system is tuned across the transition. We then find that, between the two plateaus reflecting the different carrier densities in the two Fermi-liquid phases, the Hall coefficient can develop a peak in the critical crossover regime, like in recent experimental findings, in the parameter regime of weak boson damping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗