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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
DOE OSTI2021

Efficient Iridium Catalysts for Formic Acid Dehydrogenation: Investigating the Electronic Effect on the Elementary β-Hydride Elimination and Hydrogen Formation Steps

We report herein a series of Cp*Ir complexes containing a rigid 8-aminoquinolinesulfonamide moiety as highly efficient catalysts for the dehydrogenation of formic acid (FA). The complex [Cp*Ir(L)Cl] (HL = N-(quinolin-8-yl)benzenesulfonamide) displayed a high turnover frequency (TOF) of 2.97 × 10 4 h –1 and a good stability (>100 h) at 60 °C. Comparative studies of [Cp*Ir(L)Cl] with the rigid ligand and [Cp*Ir(L')Cl] (HL' = N-propylpypridine-2-sulfonamide) without the rigid aminoquinoline moiety demonstrated that the 8-aminoquinoline moiety could dramatically enhance the stability of the catalyst. The electron-donating ability of the N,N'-chelating ligand was tuned by functionalizing the phenyl group of the L ligand with OMe, Cl, and CF 3 to have a systematical perturbation of the electronic structure of [Cp*Ir(L)Cl]. Experimental kinetic studies and density functional theory (DFT) calculations on this series of Cp*Ir complexes revealed that (i) the electron-donating groups enhance the hydrogen formation step while slowing down the β-hydride elimination and (ii) the electron-withdrawing groups display the opposite effect on these reaction steps, which in turn leads to lower optimum pH for catalytic activity compared to the electron-donating groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

3d-3d correspondence for mapping tori

One of the main challenges in 3d-3d correspondence is that no existent approach offers a complete description of 3d N = 2 SCFT T [ M 3 ] — or, rather, a “collection of SCFTs” as we refer to it in the paper — for all types of 3-manifolds that include, for example, a 3-torus, Brieskorn spheres, and hyperbolic surgeries on knots. The goal of this paper is to overcome this challenge by a more systematic study of 3d-3d correspondence that, first of all, does not rely heavily on any geometric structure on M 3 and, secondly, is not limited to a particular supersymmetric partition function of T [ M 3 ]. In particular, we propose to describe such “collection of SCFTs” in terms of 3d N = 2 gauge theories with “non-linear matter” fields valued in complex group manifolds. As a result, we are able to recover familiar 3-manifold invariants, such as Turaev torsion and WRT invariants, from twisted indices and half-indices of T [ M 3 ], and propose new tools to compute more recent q -series invariants Z * ( M 3 ) in the case of manifolds with b 1 > 0. Although we use genus-1 mapping tori as our “case study,” many results and techniques readily apply to more general 3-manifolds, as we illustrate throughout the paper.

Topological Field Theories↗
DOE OSTI2022

Microscopic Theory of Nuclear Fission

Nuclear fission represents the ultimate test for microscopic theories of nuclear structure and reactions. Fission is a large-amplitude, time-dependent phenomenon taking place in a self-bound, strongly interacting many-body system. It should, at least in principle, emerge from the complex interactions of nucleons within the nucleus. The goal of microscopic theories is to build a consistent and predictive theory of nuclear fission by using as only ingredients protons and neutrons, nuclear forces, and quantum many-body methods. Thanks to a constant increase in computing power, such a goal has never seemed more within reach. This chapter gives an overview both of the set of techniques used in microscopic theories to describe the fission process and of some recent successes achieved by this class of methods.

Nuclear physics and radiation physics↗
DOE OSTI2025

Can Dithiolate Ligands Report Electronic Communication in Transuranium Complexes?

The actinide elements, and in particular the transuranic elements, are some of the least studied elements on the periodic table. The degree of covalent bonding between the actinide elements and ligands in is not well understood. Dithiolene/dithiolate ligands provided a revolution in understanding the bonding in transition metals through the geometric differences the ligands exhibit when interacting with electron deficient or abundant metal centers. To test if dithiolene/dithiolate ligands can be utilized to provide an analogous understanding of transuranic element bonding, density functional theory calculations have been performed on U, Np, and Pu complexes of the formula Cp 2 AnS 2 C 2 H 2 . These calculations show that the dithiolate maintains its non-innocent redox active nature when bonding with actinides. In stark contrast with the transition metals, the direction of electron donation is reversed, in the actinide series the direction of electron flow in high fold cases is from the metal center to the ligand. In conclusion, this suggests the communication is between the π* system of the ligand interacting with the f orbitals of the metal center.

Prange, Micah P. [Pacific Northwest National Labor↗
DOE OSTI2021

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

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

Symmetry reduction of tensor networks in many-body theory: I. Automated symbolic evaluation of SU(2) algebra

Abstract The ongoing progress in (nuclear) many-body theory is accompanied by an ever-rising increase in complexity of the underlying formalisms used to solve the stationary Schrödinger equation. The associated working equations at play in state-of-the-art ab initio nuclear many-body methods can be analytically reduced with respect to angular-momentum, i.e. SU (2), quantum numbers whenever they are effectively employed in a symmetry-restricted context. The corresponding procedure constitutes a tedious and error-prone but yet an integral part of the implementation of those many-body frameworks. Indeed, this symmetry reduction is a key step to advance modern simulations to higher accuracy since the use of symmetry-adapted tensors can decrease the computational complexity by orders of magnitude. While attempts have been made in the past to automate the (anti-) commutation rules linked to Fermionic and Bosonic algebras at play in the derivation of the working equations, there is no systematic account to achieve the same goal for their symmetry reduction. In this work, the first version of an automated tool performing graph-theory-based angular-momentum reduction is presented. Taking the symmetry-unrestricted expressions of a generic tensor network as an input, the code provides their angular-momentum-reduced form in an error-safe way in a matter of seconds. Several state-of-the-art many-body methods serve as examples to demonstrate the generality of the approach and to highlight the potential impact on the many-body community.

Tichai, A.↗
DOE OSTI2022

Interlayer spacing in pillared and grafted MCM-22 type silicas: density functional theory analysis versus experiment

Pillaring of synthetic layered crystalline silicates and aluminosilicates provides a strategy to enhance their adsorption and separation performance, and can facilitate the understanding of such behavior in more complex natural clays. In this study, we perform the first-principles density functional theory calculations for the pillaring of the pure silica polymorph of an MCM-22 type molecular sieve. Starting with a precursor material MCM-22P with fully hydroxylated layers, a pillaring agent, (EtO) 3 SiR, can react with hydroxyl groups (–OH) on adjacent internal surfaces, 2(–OH) + (EtO) 3 SiR + H 2 O → (–O) 2 SiOHR + 3EtOH, to form a pillar bridging these surfaces, or with a single hydroxyl, –OH + (EtO) 3 SiR + 2H 2 O → (–O)Si(OH) 2 R + 3EtOH, grafting to one surface. For computational efficiency, we replace the experimental organic ligand, R, by a methyl group. We find that the interlayer spacing in MCM-22 is reduced by 2.66 Å relative to weakly bound layers in the precursor MCM-22P. Including (–O) 2 SiR bridges for 50% (100%) of the hydroxyl sites in MCM-22P increases the interlayer spacing relative to MCM-22 by 2.52 Å (2.46 Å). For comparison, we also analyze the system where all –OH groups in MCM-22P are replaced by non-bridging grafted (–O)Si(OH) 2 R which results in a smaller interlayer spacing expansion of 2.17 Å relative to MCM-22. Our results for the interlayer spacing in the pillared materials are compatible with experimental observations for a similar MCM-22 type material with low Al content (Si : Al = 51 : 1) of an expansion relative to MCM-22 of roughly 2.8 Å and 2.5 Å from our x-ray diffraction and scanning transmission electron microscopy analyses, respectively. The latter analysis reveals significant variation in individual layer spacings.

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

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Modeling dark- and light-induced crystal structures and single-crystal optical absorption spectra of ruthenium-based complexes that undergo SO 2 -linkage photoisomerization

A family of coordination complexes of the type [Ru(SO 2 )(NH 3 ) 4 X] m+ Y n – (m, n = 1 or 2) exhibit optical switching capabilities in their single-crystal states. This striking effect is caused by the light-induced formation of SO 2 -linkage photoisomers, which are metastable if kept at suitably cool temperatures. We modeled the dark- and light-induced states of these large crystalline complexes via plane-wave (PW)- and molecular-orbital (MO)-based density functional theory (DFT) and time-dependent DFT in order to calculate their structural and optical properties; the calculated results are compared with experimental data. We show that the PW-DFT-based periodic models replicate the structural properties of these complexes more effectively than the MO-DFT-based molecular-fragment models, observing only small deviations in key bond lengths relative to the experimentally derived crystal structures. The periodic models were also found to more effectively simulate trends seen in experimental optical absorption spectra, with optical absorbance and coverage of the visible region increasing with the formation of the photoinduced geometries. The contribution of the metastable photoisomeric species more heavily focuses on the lower-energy end of the spectra. Spectra generated from the molecular-fragment models are limited by the geometry of the fragment used and the number of excited-state roots considered in those calculations. In general, periodic models outperform the molecular-fragment models owing to their ability to better appreciate the periodic phenomena that are present in these crystalline materials as opposed to MO approaches, which are finite methods. We thus demonstrate that PW-DFT-based periodic models should be considered as a more than viable method for simulating the optical and electronic properties of these single-crystal optical switches.

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

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

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

Toward simulating quantum field theories with controlled phonon-ion dynamics: A hybrid analog-digital approach

Quantum field theories are the cornerstones of modern physics, providing relativistic and quantum mechanical descriptions of physical systems at the most fundamental level. Simulating real-time dynamics within these theories remains elusive in classical computing. This provides a unique opportunity for quantum simulators, which hold the promise of revolutionizing our simulation capabilities. Trapped-ion systems are successful quantum-simulator platforms for quantum many-body physics and can operate in digital, or gate-based, and analog modes. Inspired by the progress in proposing and realizing quantum simulations of a number of relativistic quantum field theories using trapped-ion systems, and by the hybrid analog-digital proposals for simulating interacting boson-fermion models, we propose hybrid analog-digital quantum simulations of selected quantum field theories, taking recent developments to the next level. On one hand, the semi-digital nature of this proposal offers more flexibility in engineering generic model interactions compared with a fully-analog approach. On the other hand, encoding the bosonic fields onto the phonon degrees of freedom of the trapped-ion system allows a more efficient usage of simulator resources, and a more natural implementation of intrinsic quantum operations in such platforms. This opens up new ways for simulating complex dynamics of e.g., Abelian and non-Abelian gauge theories, by combining the benefits of digital and analog schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coordination-Induced N–H Bond Weakening in a Molybdenum Pyrrolidine Complex: Isotopic Labeling Provides Insight into the Pathway for H 2 Evolution

The synthesis and characterization of a cationic molybdenum pyrrolidine complex is described that exhibits significant coordination-induced N–H bond weakening. Here, the N–H bond dissociation free energy (BDFE) of the coordinated pyrrolidine in [( Ph Tpy)(PPh 2 Me) 2 Mo(NH(pyrr))][BArF 24 ] ([1-NH(pyrr)] + ; PhTpy = 4'-Ph-2,2',6',2''-terpyridine, NH(pyrr) = pyrrolidine, ArF 24 = [C 6 H 3 -3,5-(CF 3 ) 2 ] 4 ) was determined to be between 41–51 kcal mol -1 by thermochemical analysis and supported by a density functional theory (DFT) computed value of 48 kcal mol -1 . The complex [1-NH(pyrr)] + underwent proton-coupled electron transfer (PCET) to 2,4,6-tri-tert-butylphenoxyl radical, as well as spontaneous H 2 evolution upon gentle heating to furnish the corresponding molybdenum pyrrolidide complex, [( Ph Tpy)(PPh 2 Me) 2 Mo(N(pyrr))][BArF 24 ] ([1-N(pyrr)] + ). Thermolysis of the deuterated isotopolog, [1-ND(pyrr)] + still produced H 2 with concomitant incorporation of the isotopic label into the pyrrolidide ligand in the product [(1-N(pyrr-d n )] + (n = 0–2), consistent with an H 2 evolution pathway involving intramolecular H–H bond formation followed by an intermolecular product-forming PCET step. These observations provide context for understanding H 2 evolution in the nonclassical ammine complex [( Ph Tpy)(PPh 2 Me) 2 Mo(NH 3 )][BArF 24 ] ([1-NH 3 ] + ) and are supported by DFT-computed reaction thermochemistry. Overall, these studies offer rare insight into the H 2 formation pathway in nonclassical amine complexes with N–H BDFEs below the thermodynamic threshold for H 2 evolution and inform the development of well-defined, thermody-namically potent PCET reagents.

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

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2024

Can we mitigate tropical cyclone formation using aerosols? A review of cyclogenesis and aerosol effects as a theoretical basis

Tropical cyclogenesis is a complex phenomenon, involving processes across multiple scales, with various theories proposed to explain its mechanisms. This paper provides an up-to-date summary of the established physical mechanisms and controlling factors of cyclogenesis. Understanding the nature of tropical cyclone formation is necessary to improve tropical cyclone forecasts and inform tropical cyclone mitigation strategies targeting early-stage intervention. Despite the disparities among theories, the critical role of persistent, intensified convection in tropical cyclone initiation is universally acknowledged. Therefore, our hypothesis is that interrupting convective development using aerosol injection is a potential strategy for mitigating tropical cyclone formation. Targetted injection of coarse and fine aerosols at the formation stage are two proposed potential approaches for mitigating tropical cyclones. These methods have completely different microphysical effects but each may ultimately interfere with the formation of the nascent cyclone. In conclusion, these potential interventions open new areas of further research, which are important and necessary to establish a solid scientific foundation for cyclogenesis mitigation.

Tropical cyclone formation↗