SEARCH · Engineering Papers

Results for “complex systems 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Chain trajectories, domain shapes, and terminal boundaries in block copolymers

The packing geometry of macromolecules in complex mesophases is of key importance to self-organization in synthetic and biological soft materials. While approximate or heuristic models rely on often-untested assumptions about how flexible molecules “fit in” to distinct locations of complex assemblies, physical assemblies derive from ensembles of fluctuating conformations, obscuring the connection between mesophase geometry and the underlying arrangements. Here, we present an approach to extract and analyze features of molecular packing in diblock block copolymer (BCP) melts, a prototypical soft matter system, based on the statistical description of chain conformations in self-consistent field (SCF) theory. We show how average BCP chain trajectories in ordered morphologies can be analyzed from the SCF-derived orientational order parameter of chain segments. We use these extracted trajectories to analyze the features of local packing geometry, including chain bending and tilt, as well as the terminal boundaries that delineate distinct domains in ordered BCP morphologies. In conclusion, we illustrate this analysis by focusing on measurable features of packing frustration in 2D (columnar) and 3D (spherical and bicontinuous) morphologies, notably establishing an explicit link between chain conformations in complex morphologies and their medial geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

X-ray induced Coulomb explosion imaging of transient excited-state structural rearrangements in CS2

Abstract Structural imaging of transient excited-state species is a key goal of molecular physics, promising to unveil rich information about the dynamics underpinning photochemical transformations. However, separating the electronic and nuclear contributions to the spectroscopic observables is challenging, and typically requires the application of high-level theory. Here, we employ site-selective ionisation via ultrashort soft X-ray pulses and time-resolved Coulomb explosion imaging to interrogate structural dynamics of the ultraviolet photochemistry of carbon disulfide. This prototypical system exhibits the complex motifs of polyatomic photochemistry, including strong non-adiabatic couplings, vibrational mode couplings, and intersystem crossing. Immediately following photoexcitation, we observe Coulomb explosion signatures of highly bent and stretched excited-state geometries involved in the photodissociation. Aided by a model to interpret such changes, we build a comprehensive picture of the photoinduced nuclear dynamics that follows initial bending and stretching motions, as the reaction proceeds towards photodissociation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Principles of two-dimensional terahertz spectroscopy of collective excitations: The case of Josephson plasmons in layered superconductors

Two-dimensional terahertz spectroscopy (2DTS), a terahertz analog of nuclear magnetic resonance, is a new technique poised to address many open questions in complex condensed matter systems. The conventional theoretical framework used ubiquitously for interpreting multidimensional spectra of discrete quantum level systems is, however, insufficient for the continua of collective excitations in strongly correlated materials. Here, in this work, we develop a theory for 2DTS of a model collective excitation, the Josephson plasma resonance in layered superconductors. Starting from a mean-field approach at temperatures well below the superconducting phase transition, we obtain expressions for the multidimensional nonlinear responses that are amenable to intuition derived from the conventional single-mode scenario. We then consider temperatures near the superconducting critical temperature T c , where dynamics beyond mean-field become important and conventional intuition fails. As fluctuations proliferate near T c , the dominant contribution to nonlinear response comes from an optical parametric drive of counterpropagating Josephson plasmons, which gives rise to 2D spectra that are qualitatively different from the mean-field predictions. As such, and in contrast to one-dimensional spectroscopy techniques, such as third harmonic generation, 2DTS can be used to directly probe thermally excited finite-momentum plasmons and their interactions. Our theory can readily be tested in cuprates, and we discuss implications beyond the present context of Josephson plasmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Nonequilibrium thermodynamics of erasure with superconducting flux logic

We implement a thermal-fluctuation-driven logical bit reset on a superconducting flux logic cell. We show that the logical state of the system can be continuously monitored with only a small perturbation to the thermally activated dynamics at 500 mK. We use the trajectory information to derive a single-shot estimate of the work performed on the system per logical cycle. We acquire a sample of 10 5 erasure trajectories per protocol and show that the work histograms agree with both microscopic theory and global fluctuation theorems. The results demonstrate how to design and diagnose complex, high-speed, and thermodynamically efficient computing using superconducting technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Exploring the nexus of many-body theories through neural network techniques: the tangent model

Abstract In this paper, we present a physically informed neural network (NN) representation of the effective interactions associated with coupled-cluster downfolding models to describe chemical systems and processes. The NN representation not only allows us to evaluate the effective interactions efficiently for various geometrical configurations of chemical systems corresponding to various levels of complexity of the underlying wave functions, but also reveals that the bare and effective interactions are related by a tangent function of some latent variables. We refer to this characterization of the effective interaction as a tangent model. We discuss the connection between this tangent model for the effective interaction with the previously developed theoretical analysis that examines the difference between the bare and effective Hamiltonians in the corresponding active spaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Understanding the Applicability of Gas/Particle Partitioning Theory to SOA using Novel Chemically-Speciated Partitioning Measurements (Final Report)

The main objective of this research grant project was to investigate the physical-chemical processes and properties of particles and partitioning gases that influence gas/particle (G/P) partitioning and timescales for equilibration, to directly determine the conditions under which equilibrium G/P partitioning theory can be applied to secondary organic aerosol (SOA). SOA constitutes a major component of submicron particles in the atmosphere, which have important implications for air quality, human and ecosystem health, radiation balance, and climate; however SOA formation, evolution and losses are not well understood. Use of novel on-line measurements of chemically-speciated partitioning measurements of both gas-phase and particle-phase components was a central component to advancing these scientific goals. The key approach was to conduct environmental chamber studies with systems of increasing complexity to determine G/P partitioning rates and equilibrium states for chemically-speciated SOA components. The project builds on our previous work on gas-wall partitioning that must be understood in order to study G/P partitioning in chambers, and to help extrapolate the results to be atmospherically relevant. The goal was to employ measurements and modeling to systematically determine when G/P partitioning theory accurately describes the distribution of oxidized organic compounds in the gas and aerosol phases, and when and why it does not (e.g., aerosol phase-state kinetic limitations, non-ideal mixing, phase separation, oligomer formation,), for chemical and physical conditions typically found throughout the continental boundary layer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Nickel(II) Metal Complexes as Optically Addressable Qubit Candidates

The inherent atomic level structural control of synthetic chemistry enables the creation of qubits, the base units of a quantum information science system, designed for a target application. For quantum sensing applications, enabling optical read-out of spin in tunable molecular systems, akin to defect-based systems, would be transformative. This approach would bring together molecular tunability with optical read-out technology. In theory, nickel ions in octahedral symmetry meet all the criteria for optical readout of spin. Yet, to the best of our knowledge, there are no pulse EPR studies on Ni 2+ molecules. We identified two compounds featuring highly symmetric Ni 2+ centers, thereby engendering weak zero-field splitting to enable EPR addressability:[Ni(phen) 3 ](BF 4 ) 2 (1) and [Ni(pyr 3 ) 2 ](BF 4 ) 2 (2) (phen = 1,10-phenanthroline; pyr 3 = tris-2-pyridyl-methane). Crucially, these complexes feature the requisite strong field ligands to enable emission for optical addressability. We extracted axial zero-field splitting parameters of D = +0.9 cm –1 and +2.7 cm –1 for 1 and 2, respectively, enabling pulse EPR measurements. Both compounds produce emission at λ max = 938–944 nm. Furthermore, the aggregate of these results expands the catalogue of qubit materials to Ni 2+ -based compounds and offers a future pathway for optical readout of these molecules.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Path integrals, complex probabilities and the discrete Weyl representation

Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.

Physics↗

DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Formulation of the density eigenvalue problem in neutron transport for relevant engineering applications

A new formulation of the density eigenvalue problem for the neutron transport equation is presented. This new eigenvalue, named ζ eigenvalue can be introduced freely in the transport model, acting on a selected portion of the phase space. Despite its broader applications and its connection with the nature of the transport operator, the ζ eigenvalue has been presented here mainly as a design-oriented technique for the efficient evaluation of the critical concentration for a specific nuclide (or mixture of nuclides). This new eigenvalue is particularly adequate to study the definition of the material composition in the criticality design process of a multiplying system. The method is then applied for the study of classical problems such as the critical moderation ratio and the poison concentration to control the reactor. Some results are presented in one dimensional configuration using the multi- group spherical harmonics approach. This eigenvalue formulation proves to be a convenient and useful way to attain criticality, also for complex, realistic systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

[(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2– : Foundations toward the Development of Trifluoromethylations at Unsupported Nickel

Nickel anions [(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2– were prepared by the formal addition of 3 and 4 equiv, respectively, of AgCF 3 to [(dme)NiBr 2 ] in the presence of the [PPh 4 ] + counterion. Detailed insights into the electronic properties of these new compounds were obtained through the use of density functional theory (DFT) calculations, spectroscopy-oriented configuration interaction (SORCI) calculations, X-ray absorption spectroscopy, and cyclic voltammetry. The data collectively show that trifluoromethyl complexes of nickel, even in the most common oxidation state of nickel(II), are highly covalent systems whereby a hole is distributed on the trifluoromethyl ligands, surprisingly rendering the metal to a physically more reduced state. In the cases of [(MeCN)Ni(CF 3 ) 3 ] - and [Ni(CF 3 ) 4 ] 2- , these complexes are better physically described as d 9 metal complexes. [(MeCN)Ni(CF 3 ) 3 ] - is electrophilic and reacts with other nucleophiles such as phenoxide to yield the unsupported [(PhO)Ni(CF 3 ) 3 ] 2– salt, revealing the broader potential of [(MeCN)Ni(CF 3 ) 3 ] - in the development of “ligandless” trifluoromethylations at nickel. Proof-in-principle experiments show that the reaction of [(MeCN)Ni(CF 3 ) 3 ] - with an aryl iodonium salt yields trifluoromethylated arene, presumably via a high-valent, unsupported, and formal organonickel(IV) intermediate. Evidence of the feasibility of such intermediates is provided with the structurally characterized [PPh 4 ] 2 [Ni(CF 3 ) 4 (SO 4 )], which was derived through the two-electron oxidation of [Ni(CF 3 ) 4 ] 2– .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Complex Electrical Conductivity of a Single‐Fractured Rock: Fracture‐and‐Matrix Coupling Mechanism and Aperture Size Predictions

Fractured rocks play a crucial role in myriad natural and engineered systems, including Earth's critical zone, oil/gas/geothermal reservoirs, and geological CO 2 /H 2 /waste storage systems. While complex electrical conductivity is extensively used to estimate the pore and grain sizes of conventional porous rocks and soils, it is rarely used to predict the aperture size of fractured rocks and this remains poorly understood. Here, integrating theory, simulations, and experiments, we show that under external fields, fractured rocks follow the fracture-and-matrix coupling to make the bulk complex conductivity non-linear with respect to water conductivity. We find that the relaxation time and quadrature conductivity for porous media do not apply to fractured rocks, but, instead, reasonably accurate predictions of aperture size can be made based on the true formation factor. This study unravels the fundamental mechanism governing conduction and polarization of fractured rocks and paves the way for the non-invasive investigation of global fractured rocks.

58 GEOSCIENCES↗

DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗