State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin–Orbit Coupling: Benchmark and Perspectives
Not Available
SEARCH · Engineering Papers
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.
Not Available
Risk-informed decision-making requires a probabilistic assessment of the likelihood of success of control action, given the system status. This paper presents a systematic state transition modeling approach integrating dynamic probabilistic risk assessment with a decision-making process using a dynamic Bayesian network (DBN) coupled with functional modeling. A functional model designed with multilevel flow modeling (MFM) technique was used to build a system state structure inferred by energy, mass, and information flow so that one can verify the developed model with respect to system functionality. The MFM model represents the causal relationship among the nodes, which captures the structure of process parameters and control units. Each node may have multiple possible states, and the DBN structured by the MFM model represents the time-domain transitions among the defined states. Furthermore, the MFM-DBN integrated state transition modeling is a white-box approach that allows one to draw the system's risk profile by updating the system states and supports the decisions probabilistically with physical inference. An example of a simple heating system has been used to illustrate this process, including decision-making support based on quantitative risk profile. For demonstrating its applicability to a complex system operational decision making, a case study of station blackout accident scenario leading to the seal loss of coolant accident in a nuclear power plant is presented. The proposed approach effectively provided the risk profile along time for each option so that the operators can make the best decision, which minimizes the plant risk.
While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of major facilitator superfamily transporters. With a conformation-selective nanobody, we determined a low-sugar affinity inward-facing Na + -bound cryoEM structure. The available outward-facing sugar-bound structures showed that the N- and C-terminal residues of the inner barrier contribute to the sugar selectivity. The inward-open conformation shows that the sugar selectivity pocket is also broken when the inner barrier is broken. Isothermal titration calorimetry measurements revealed that this inward-facing conformation trapped by this nanobody exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is supported by molecular dynamics simulations. Furthermore, the hydron/deuterium exchange mass spectrometry analyses allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.
Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.
An isolation coupler for coupling a functional element to a support structure includes a first bracket. The first bracket includes a number of first-bracket sides. The number of first-bracket sides forms a closed polygonal shape, in plan view. The isolation coupler further includes a number of isolators coupled to each one of the first-bracket sides. The isolation coupler also includes a second bracket. The second bracket includes a number of second-bracket sides. The second bracket sides are coupled to the isolators. The number of second-bracket sides is equal to the number of first-bracket sides and forms the closed polygonal shape, in plan view. The isolators separate each one of the first-bracket sides from a corresponding one of the second-bracket sides to attenuate a load transferred from the first bracket to the second bracket.
An isolation coupler for coupling a functional element to a support structure includes a first bracket. The first bracket includes a number of first-bracket sides. The number of first-bracket sides forms a closed polygonal shape, in plan view. The isolation coupler further includes a number of isolators coupled to each one of the first-bracket sides. The isolation coupler also includes a second bracket. The second bracket includes a number of second-bracket sides. The second bracket sides are coupled to the isolators. The number of second-bracket sides is equal to the number of first-bracket sides and forms the closed polygonal shape, in plan view. The isolators separate each one of the first-bracket sides from a corresponding one of the second-bracket sides to attenuate a load transferred from the first bracket to the second bracket.
Not Available
Concentrating Solar Power (CSP) is a promising solar technology for electricity generation with thermal energy storage and with the additional benefit of industrial heat production. Wind loading on CSP collector structures, such as parabolic troughs or heliostats, is one of the primary drivers of their structural design costs. In particular, dynamic wind loading is a major source of uncertainty in the collector design process, which heavily relies on wind tunnel testing. In the field, the turbulent nature of the incoming wind creates fluctuating loads (support structure loads and resulting mirror deflections) on the collectors, with impacts on fatigue lifetime and optical performance. As is well known, wind tunnel tests cannot entirely reproduce the complex turbulent wind conditions typically observed at full-scale plants. To shed light on this topic, NREL initiated a field campaign at the operational Nevada Solar One (NSO) powerplant that uses parabolic troughs as solar collectors. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs. We use the published 2-year dataset of high-resolution combined wind and structural loads measurements [1] to characterize the dynamic structural wind response. For quantifying dynamic wind loading, we apply the concept of admittance functions, which are spectral transfer functions that couple the turbulent wind to resulting structural loads (aerodynamic admittance), and to the structural response (mechanical admittance). In practice, aerodynamic admittance describes which turbulent eddy sizes are effective in creating structural loads. The mechanical admittance describes in which frequency ranges these loads are reinforced or dampened by the structure. While these functions are an established concept in civil engineering, their recent application to a single full-scale heliostat [2] proved their broader applicability to CSP collectors. Here, we present a characterization of admittance functions for full-scale parabolic trough collectors and show how wind characteristics (mean wind speed and direction, turbulent kinetic energy, turbulent length scales), the sun-tracking trough angle, and row position alter the admittance functions. Further, we study to which extent the admittance functions are universal for a specific trough geometry and how our findings compare to reported heliostat results. References [1] https://data.openei.org/submissions/5938. [2] Blume, K., Roger, M., and Pitz-Paal, R. 2023b. "Simplified analytical model to describe wind loads and wind-induced tracking deviations of heliostats." Solar Energy, 256, 96-109. https://doi.org/10.1016/j.solener.2023.03.055.
Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.
The coupling of early-time high-altitude electromagnetic pulse signals into a power generation facility was investigated. This study encompasses both electromagnetic simulation and on-site experimental efforts that are relevant to shielding effectiveness estimation with a passive technique. As an initial step in understanding the attenuation models involved, simulations were carried out as a function of frequency, building properties, and signal incidence angle. Testing was then performed at the US Department of Energy’s Oak Ridge National Laboratory campus to identify the equipment needed for signal collection and to fine-tune the experimental procedure. Finally, as an application of the models and methods presented, measurements were conducted at a power generation plant to derive the field transfer functions and coupled high-altitude electromagnetic pulse waveforms for various locations within the facility.
We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.
The recent pioneering campaigns conducted by the Lead Radius Experiment (PREX) and the Calcium Radius Experiment (CREX) Collaborations have uncovered major deficiencies in the theoretical description of some fundamental properties of atomic nuclei. Following a recent refinement to the isovector sector of covariant energy density functionals we present here additional improvements to the functional by including both tensor couplings and an isoscalar-isovector mixing term in the scalar sector. Motivated by the distinct surface properties of calcium and lead, we expect that the tensor terms that generate derivative couplings will help break the linear correlation between the neutron skin thickness of these two nuclei. Moreover, the addition of these new terms mitigates most of the problems identified by Reed et al. in describing the properties of both finite nuclei and neutron stars. While significant progress has been made in reconciling the PREX-CREX results without compromising other observables, the final resolution awaits the completion of a proper calibration for this new class of functionals. As a result, we expect that powerful reduced basis methods used recently to create efficient emulators will be essential to accomplish this task.
Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.
Although perovskite oxyhydrides (POHs) have shown promise as a support for selective hydrogenation, it is unknown if they themselves can catalyze alkyne semi-hydrogenation. Here, we use first-principles density functional theory, coupled with microkinetic modeling, to investigate acetylene semi-hydrogenation on a prototypical POH, cubic BaTiO 2.5 H 0.5 (BTOH). Two different mechanisms are examined on a representative surface of BTOH under the reaction conditions: both are based on the Horiuti–Polanyi mechanism, but the way of H 2 dissociation is different. In mechanism 1, a lattice hydride H atom and then a surface adsorbed H atom sequentially hydrogenate the adsorbed acetylene. In mechanism 2, two lattice hydride H atoms from the BTOH sequentially hydrogenate the adsorbed acetylene. In both mechanisms, the H atoms are replenished from gas phase H 2 dissociation. Using density functional theory (DFT), we have calculated pathways for these two mechanisms in hydrogenation of acetylene to ethylene, as well as further to ethane. Microkinetic modeling based on the DFT energetics indicates that at 523 K, hydrogenation to ethylene and ethane via mechanism 1 occurs at a rate up to two orders of magnitude faster than via mechanism 2. A selectivity analysis for the temperature range of 373–673 K shows that the product observed is essentially only ethylene within the more active mechanism 1. This is because at the steady state, there is a significant amount of surface anion vacancies that facilitate heterolytic H 2 dissociation and help stabilize a vinyl intermediate. In conclusion, the present findings suggest that POHs are capable of selective hydrogenation, thanks to their rich hydrogen surface-redox chemistry.
Experiments show that Au25(SR)$\overline{_{18}}$ is a potential electrocatalyst for the electrochemical reduction of CO 2 , but elucidation of the reaction mechanisms remain challenging experimentally. We evaluate key steps in the active site formation, CO 2 reduction, and H 2 evolution on the nanocluster employing density functional theory coupled with an explicit solvent model of the electrochemical interface to calculate activation barriers. The predicted preferred pathway for activating the nanocluster involves dethiolation, resulting in an exposed Au site that is both active and selective for CO 2 reduction. Finally, the alternative S site is found to be kinetically prohibitive and does not facilitate CO 2 reduction.
Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.
The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.
Seven acentric sulfides Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) were grown by a high-temperature salt flux method. The crystal structures of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) compounds were determined by single-crystal X-ray diffraction with the aid of solid-state NMR spectroscopy. The Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) compounds are isostructural and crystallize in the Ba 6 Ag 4 Sn 4 S 16 structure type. The Sn-containing compound exhibits high structural similarity to Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) with the presence of an interstitial atomic position partially occupied by Sn atoms. The chemical bonding characteristics of Ba 6 (Cu 2.9 Sn 0.4 )Sn 4 S 16 were understood with electron localization function calculations coupled with crystal orbital Hamilton population calculations. The Ba–S and Cu–S interactions are dominantly ionic, but the Sn–S interactions consist of strong covalent bonding characteristics in Ba 6 (Cu 2.9 Sn 0.4 )Sn 4 S 16 . The monovalent Cu atoms, mixed with certain metals with various oxidation states, significantly shift the optical properties of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) compounds. This results in a good balance between the second-harmonic-generation (SHG) response and laser damage threshold (LDT). Ba 6 (Cu 1.9 Zn 1.1 )Sn 4 S 16 possesses a high SHG response and a high LDT of 2.8 × AGS and 3 × AGS, respectively. Here, a density functional theory calculation revealed that CuS 4 and SnS 4 tetrahedra significantly contribute to the SHG response in Ba 6 (Cu 2 Mg)Sn 4 S 16 , which also confirmed that CuS 4 tetrahedra are crucial for the stability and optical properties of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) compounds revealed by electronic structure analysis.