Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density matrix methods”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Intramolecular hydrogen bonding analysis

The quasi-atomic orbital (QUAO) bonding analysis is used to study intramolecular hydrogen bonding (IMHB) in salicylic acid and an intermediate that is crucial to the synthesis of aspirin. The bonding analysis rigorously explores IMHB through directly accessing information that is intrinsic to the molecular wave function, thereby bypassing the need for intrinsically biased methods. The variables that affect the strength of IMHB are determined using kinetic bond orders, QUAO populations, and QUAO hybridizations. Important properties include both the interatomic distance between hydrogen and oxygen participating in the IMHB and the hybridization on the oxygen. The bonding analysis further shows that each intramolecular hydrogen bond is a four-electron three-center bond. Here, the bonding analysis is used to understand how aromatic reactivity changes due to the effect of functional groups on the aromatic ring.

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

First-Principles Phonon Quasiparticle Theory Applied to a Strongly Anharmonic Halide Perovskite

Understanding and predicting lattice dynamics in strongly anharmonic crystals is one of the long-standing challenges in condensed matter physics. Here we propose a first-principles method that gives accurate quasiparticle (QP) peaks of the phonon spectrum with strong anharmonic broadening. On top of the conventional first-order self-consistent phonon (SC1) dynamical matrix, the proposed method incorporates frequency renormalization effects by the bubble self-energy within the QP approximation. We apply the developed methodology to the strongly anharmonic α-CsPbBr 3 that displays phonon instability within the harmonic approximation in the whole Brillouin zone. While the SC1 theory significantly underestimates the cubic-to-tetragonal phase transition temperature (T c ) by more than 50%, we show that our approach yields T c = 404–423 K, in excellent agreement with the experimental value of 403 K. Furthermore, we also demonstrate that an accurate determination of QP peaks is paramount for quantitative prediction and elucidation of phonon linewidth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

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

Extracting off-diagonal order from diagonal basis measurements

Quantum gas microscopy has developed into a powerful tool to explore strongly correlated quantum systems. However, discerning phases with topological or off-diagonal long range order requires the ability to extract these correlations from site-resolved measurements. Here, we show that a multiscale complexity measure can pinpoint the transition to and from the bond ordered wave phase of the one-dimensional extended Hubbard model with an off-diagonal order parameter, sandwiched between diagonal charge and spin density wave phases, using only diagonal descriptors. We study the model directly in the thermodynamic limit using the recently developed variational uniform matrix product states algorithm, and draw our samples from degenerate ground states related by global spin rotations, emulating the projective measurements that are accessible in experiments. Our results will have important implications for the study of exotic phases using optical lattice experiments. Published by the American Physical Society 2024

1-dimensional systems↗
DOE OSTI2020

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

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

Fiber Sorbents – A Versatile Platform for Sorption-Based Gas Separations

Increasing demand for high-purity fine chemicals and a drive for process intensification of large-scale separations have driven significant work on the development of highly engineered porous materials with promise for sorption-based separations. While sorptive separations in porous materials offer energy-efficient alternatives to longstanding thermal-based methods, the particulate nature of many of these sorbents has sometimes limited their large-scale deployment in high-throughput applications such as gas separations, for which the necessary high feed flow rates and gas velocities accrue prohibitive operational costs. These processability limitations have been historically addressed through powder shaping methods aimed at the fabrication of structured sorbent contactors based on pellets, beads or monoliths, commonly obtained as extrudates. These structures overcome limitations such as elevated pressure drops commonly recorded across powder adsorption beds but often accrue thermal limitations arising from elevated particle density and aggregation, which ultimately cap their maximum separation performance. Furthermore, the harsh mechanical strain to which powder particles are subjected during contactor fabrication, in the form of extrusion/compression forces, can result in partial pore occlusion and framework degradation, further limiting their performance. Here, we present the development of porous fiber sorbents as an alternative sorbent contactor design capable of addressing sorbent processability limitations while enabling an array of performance-maximizing heat integration capabilities. This new sorbent form factor leverages pre-existing know-how from hollow fiber spinning to produce fiber-shaped sorbent contactors through the phase inversion of known polymers in a process known as dry-jet/wet quenching. The process of phase inversion allows microporous sorbent particles to be latched onto a macroporous polymer matrix under mild processing conditions, thus making it compatible with soft porous materials prone to amorphization under traditional pelletization conditions. Sorbent fibers can be created with different geometries through control of the spinning apparatus and process, offering the possibility to produce monolithic and hollow fibers alike, the latter of which can be integrated with thermalization fluid flows. In this Account, we summarize our progress in the field of fiber sorbents from both design and application standpoints. We further guide the reader through the evolution of this field from the early inceptive work on zeolite hollow fibers to recent developments on MOF fibers. We highlight the versatile nature of fiber sorbents, both from the composition, fabrication and structure points of view, and further demonstrate how fiber sorbents offer alternative paths in tackling new and challenging chemical separation challenges like direct air capture (DAC), with a final perspective on the future of the field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A High-Throughput Solver for Marginalized Graph Kernels on GPU

Here, we present the design and optimization of a solver for efficient and high-throughput computation of the marginalized graph kernel on General Purpose GPUs. The graph kernel is computed using the conjugate gradient method to solve a generalized Laplacian of the tensor product between a pair of graphs. To cope with the large gap between the instruction throughput and the memory bandwidth of the GPUs, our solver forms the graph tensor product on-the-fly without storing it in memory. This is achieved by using threads in a warp cooperatively to stream the adjacency and edge label matrices of individual graphs by small square matrix blocks called tiles, which are then staged in registers and the shared memory for later reuse. Warps across a thread block can further share tiles via the shared memory to increase data reuse. We exploit the sparsity of the graphs hierarchically by storing only non-empty tiles using a coordinate format and nonzero elements within each tile using bitmaps. We propose a new partition-based reordering algorithm for aggregating nonzero elements of the graphs into fewer but denser tiles to further exploit sparsity. We carry out extensive theoretical analyses on the graph tensor product primitives for tiles of various density and evaluate their performance on synthetic and real-world datasets. Our solver delivers three to four orders of magnitude speedup over existing CPU-based solvers such as GraKeL and GraphKernels. The capability of the solver enables kernel-based learning tasks at unprecedented scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

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

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

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

First-principles ionized-impurity scattering and charge transport in doped materials

Scattering of carriers with ionized impurities governs charge transport in doped semiconductors. However, electron interactions with ionized impurities cannot be fully described with quantitative first-principles calculations, so their understanding relies primarily on simplified models. Here we show an ab initio approach to compute the interactions between electrons and ionized impurities or other charged defects. It includes the short- and long-range electron-defect (e-d) interactions on equal footing and allows for efficient interpolation of the e-d matrix elements. Here we combine the e-d and electron-phonon interactions in the Boltzmann transport equation to compute the carrier mobilities in doped silicon over a wide range of temperature and doping concentrations, seamlessly spanning the defect- and phonon-limited transport regimes. The individual contributions of the defect- and phonon-scattering mechanisms to the carrier relaxation times and mean-free paths are analyzed. Our method provides a powerful tool to study electronic interactions in doped materials. It broadens the scope of first-principles transport calculations, enabling studies of a wide range of doped semiconductors and oxides with application to electronics, energy, and quantum technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the speciation of actinide salts in the presence of contaminants related to molten salt reactors

Molten salt reactors (MSR) are up-and-coming Generation IV nuclear reactors with either coolant and/or radioactive fuel in molten salt form. Due to their improved safety, efficiency, affordability and convenient waste processing system, these reactors are considered a superior alternative to conventional ones. However, prior to commercial implementation, it is necessary to have an extensive understanding of the reactions occurring in the MSR. The safety parameters of the MSR are determined by the rheological properties such as density and viscosity of the molten salt and these properties are governed by the local structure of the molten salts. Due to the high operational temperatures of MSR's, corrosion plays an important role. Because of their influence, the container corrosion products become part of the fuel salt matrix. In this work, we are investigating how these corrosion products, e.g., Mo, Ni, Cr, would impact the speciation of the actinide metal centers. For the initial studies, lanthanides are used as a surrogate for actinides. Lanthanide chloride salts (LnCl3) are mixed with transition metals and excess of alkali or alkali earth metal salts or salt mixtures (LiCl, NaCl, KCl). These are added to an alumina crucible, placed in a quartz tube, and sealed under a vacuum. This is heated up to around 1000 °C and slowly cooled to a temperature above the melting point of the salt and at that temperature, the reaction is taken out of the furnace. The excess salt is removed, and the resulting product is analyzed using X-ray diffraction techniques, scanning electron microscopy, and spectroscopic methods such as Raman, UV-Visible, and FTIR spectroscopy.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Extracting Higher Central Charge from a Single Wave Function

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_ is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_ to gapping out the edge. Here, in this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν = 1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1) 2 and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2021

Anisotropic quantum Hall states in the presence of interactions with fourfold rotational symmetry

We study the effects of anisotropic interactions in the quantum Hall effect in the presence of a fourfold discrete rotational (C 4 ) symmetry. Employing the density matrix renormalization group technique on an infinite cylinder geometry (iDMRG), we calculate the anisotropy response of the Laughlin state at ν = 1/3 and the composite-Fermi liquid (CFL) state at ν = 1/2. We find that the anisotropy transferred from the interaction potential to the ν = 1/3 state is stronger when compared with the complementary case of an anisotropic band. Further, the strength of anisotropy reduces as the interaction is made shorter ranged. Quite surprisingly, at ν = 1/2, the deformation in the CF Fermi surface changes sign as the interaction range is reduced. Our results imply that the short-distance and long-distance parts of the interaction potential have opposite effects on the CFL state in the presence of C 4 -symmetric anisotropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Interpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices

Due to the presence of non-linear equations, iterative Green’s function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge resonance and covalent characters of these solutions. When applied within the unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, and biquadratic spin–spin interactions. These Hamiltonian parameterizations are compared to previous wave-function calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗