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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data

In this work, we present a machine learning approach to calculating electronic specific heat capacities for a variety of benchmark molecular systems. Our models are based on data from density matrix quantum Monte Carlo, which is a stochastic method that can calculate the electronic energy at finite temperature. As these energies typically have noise, numerical derivatives of the energy can be challenging to find reliably. In order to circumvent this problem, we use Gaussian process regression to model the energy and use analytical derivatives to produce the specific heat capacity. From there, we also calculate the entropy by numerical integration. We compare our results to cubic splines and finite differences in a variety of molecules in which Hamiltonians can be diagonalized exactly with full configuration interaction. We finally apply this method to look at larger molecules where exact diagonalization is not possible and make comparisons with more approximate ways to calculate the specific heat capacity and entropy.

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

Efficient simulation of open quantum systems coupled to a reservoir through multiple channels

It is challenging to simulate open quantum systems that are connected to a reservoir through multiple channels. For example, vibrations may induce fluctuations in both energy gaps and electronic couplings, which represent two independent channels of system–bath couplings. Systems of this kind are ubiquitous in the processes of excited state radiationless decay. Combined with density matrix renormalization group (DMRG) and matrix product states (MPS) methods, we develop an interaction-picture chain mapping strategy for vibrational reservoirs to simulate the dynamics of these open systems, resulting in time-dependent spatially local system–bath couplings in the chain-mapped Hamiltonian. This transformation causes the entanglement generated by the system–bath interactions to be restricted within a narrow frequency window of vibrational modes, enabling efficient DMRG/MPS dynamical simulations. We demonstrate the utility of this approach by simulating singlet fission dynamics using a generalized spin-boson Hamiltonian with both diagonal and off-diagonal system–bath couplings. This approach generalizes an earlier interaction-picture chain mapping scheme, allowing for efficient and exact simulation of systems with multi-channel system–bath couplings using matrix product states, which may further our understanding of nonlocal exciton–phonon couplings in exciton transport and the non-Condon effect in energy and electron transfer.

Chemistry↗
DOE OSTI2021

Witnessing entanglement in quantum magnets using neutron scattering

We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF 3 . We apply three entanglement witnesses—one tangle, two tangle, and quantum Fisher information—to its inelastic neutron spectrum and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find that each witness provides direct access to entanglement. Of these, quantum Fisher information is the most robust experimentally and indicates the presence of at least bipartite entanglement up to at least 50 K, corresponding to around 10% of the spinon zone-boundary energy. We apply quantum Fisher information to higher spin-S Heisenberg chains and show theoretically that the witnessable entanglement gets suppressed to lower temperatures as the quantum number increases. Finally, we outline how these results can be applied to higher dimensional quantum materials to witness and quantify entanglement.

1-dimensional systems↗
DOE OSTI2021

Abelian SU(N) 1 chiral spin liquids on the square lattice

In the physics of the Fractional Quantum Hall (FQH) effect, a zoo of Abelian topological phases can be obtained by varying the magnetic field. Aiming to reach the same phenomenology in spin-like systems, in this work we propose a family of SU(N)-symmetric models in the fundamental representation, on the square lattice with short-range interactions restricted to triangular units, a natural generalization for arbitrary N of an SU(3) model studied previously where time-reversal symmetry is broken explicitly. Guided by the recent discovery of SU(2) 1 and SU(3) 1 chiral spin liquids (CSL) on similar models we search for topological SU(N) 1 CSL in some range of the Hamiltonian parameters via a combination of complementary numerical methods such as exact diagonaliza- tions (ED), infinite density matrix renormalization group (iDMRG) and infinite Projected Entangled Pair State (iPEPS). Extensive ED on small (periodic and open) clusters up to N = 10 and an innovative SU(N)-symmetric version of iDMRG to compute entanglement spectra on (infinitely-long) cylinders in all topological sectors pro- vide unambiguous signatures of the SU(N) 1 character of the chiral liquids. An SU(4)-symmetric chiral PEPS, constructed in a manner similar to its SU(2) and SU(3) analogs, is shown to give a good variational ansatz of the N = 4 ground state, with chiral edge modes originating from the PEPS holographic bulk-edge correspondence. Finally, we discuss the possible observation of such Abelian CSL in ultracold atom setups where the possibility of varying N provides a tuning parameter similar to the magnetic field in the physics of the FQH effect.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1980

Collisional redistribution of radiation. I - The density matrix

As a first step in deriving accurate, calculable expressions for collisional redistribution of radiation, valid outside the impact regime a method for calculating an accurate binary-collision density operator for an atom in the presence of a driving field and perturbers is presented. Projection-operator techniques to establish first the precise validity of the Markoff approximation for radiative relaxation are used. Using the same techniques, it is shown how the collisional-relaxation problem may be analyzed in a practical manner outside the Markoff approximation. The effect of correlations between radiative and collisional events is included in a consistent way and the physical implications of these correlation effects are demonstrated for a simple example, the two-level atom with nondegenerate levels.

Burnett, K.↗
DOE OSTI2024

Density-matrix mean-field theory

Mean-field theories have proven to be efficient tools for exploring diverse phases of matter, complementing alternative methods that are more precise but also more computationally demanding. Conventional mean-field theories often fall short in capturing quantum fluctuations, which restricts their applicability to systems with significant quantum effects. In this article, we propose an improved mean-field theory, density-matrix mean-field theory (DMMFT). DMMFT constructs effective Hamiltonians, incorporating quantum environments shaped by entanglements, quantified by the reduced density matrices. Therefore, it offers a systematic and unbiased approach to account for the effects of fluctuations and entanglements in quantum ordered phases. As demonstrative examples, we show that DMMFT can not only quantitatively evaluate the renormalization of order parameters induced by quantum fluctuations, but can also detect the topological quantum phases. Additionally, we discuss the extensions of DMMFT for systems at finite temperatures and those with disorders. Our work provides an efficient approach to explore phases exhibiting unconventional quantum orders, which can be particularly beneficial for investigating frustrated spin systems in high spatial dimensions.

Physics↗
NASA NTRS1978

Pulse-echo probe of rock permeability near oil wells

Processing method involves sequential insonifications of borehole wall at number of different frequencies. Return signals are normalized in amplitude, and root-mean-square (rms) value of each signal is determined. Values can be processed to yield information on size and number density of microfractures at various depths in rock matrix by using averaging methods developed for pulse-echo technique.

Narasimhan, K. Y.↗
DOE OSTI2023

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

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

Holographic Quantum Simulation of Strongly Correlated Electron Systems

The project aimed to demonstrate a new holographic quantum simulation approach and co‐ designed quantum hardware to tackle three specific problems that fall within the broad umbrella of unraveling the physics of strongly correlated electron systems (SCES). These tasks were: (1) holographic preparation of ground‐ and thermal‐ states of correlated magnetic and electronic systems including quasi‐2d frustrated‐spin, Fermi‐Hubbard, and fractional quantum Hall (FQH) systems, (2) holographic‐simulation of long‐time out‐of‐equilibrium dynamics and (3) holographic analogs of embedding methods such as dynamical mean‐ field theory (DMFT) and density‐matrix embedding theory (DMET) to solve systems with complex structure or long‐range interactions. These tasks are prototypes for the kinds of material simulation problems of interest to BES, such as the simulation of multiferroic materials, perovskite photovoltaics and high‐temperature superconductors, that tax the capabilities of the most powerful classical supercomputers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Quantum theory of an optical maser. VI - Transient behavior.

The transient behavior of a laser is discussed using the quantum theory as did Scully and Lamb. The formal solution of the density-matrix equation is expressed in terms of exponentially decaying eigenmodes. Some of the lower decay constants are obtained numerically. The equations for the moments of the density matrix are then derived and solved by a truncation method. The equations of motion are integrated numerically for the case where the average number of photons in a laser cavity has the realistically large value 1.3 x 100,000. An alternative Fokker-Planck-equation approach is discussed.

Wang, Y. K.↗
DOE OSTI2023

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗
NASA NTRS1991

An efficient formulation and implementation of the analytic energy gradient method to the single and double excitation coupled-cluster wave function - Application to Cl2O2

The analytic energy gradient for the single and double excitation coupled-cluster (CCSD) wave function has been reformulated and implemented in a new set of programs. The reformulated set of gradient equations have a smaller computational cost than any previously published. The iterative solution of the linear equations and the construction of the effective density matrices are fully vectorized, being based on matrix multiplications. The new method has been used to investigate the Cl2O2 molecule, which has recently been postulated as an important intermediate in the destruction of ozone in the stratosphere. In addition to reporting computational timings, the CCSD equilibrium geometries, harmonic vibrational frequencies, infrared intensities, and relative energetics of three isomers of Cl2O2 are presented.

Rendell, Alistair P.↗
DOE OSTI2021

Analytic techniques for solving the transport equations in electroweak baryogenesis

We develop an efficient method for solving transport equations, particularly in the context of electroweak baryogenesis. It provides fully-analytical results under mild approximations and can also test semi-analytical results, which are applicable in more general cases. Key elements of our method include the reduction of the second-order differential equations to first order, representing the set of coupled equations as a block matrix of the particle densities and their derivatives, identification of zero modes, and block decomposition of the matrix. We apply our method to calculate the baryon asymmetry of the Universe (BAU) in a Standard Model effective field theory framework of complex Yukawa couplings to determine the sensitivity of the resulting BAU to modifications of various model parameters and rates, and to estimate the effect of the commonly-used thin-wall approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Hybrid gausslet/Gaussian basis sets

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain.

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

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Tunable Magnetic and Optical Anisotropy in ZrO 2 ‐Co Vertically Aligned Nanocomposites

Abstract Metamaterials have gained great research interest in recent years owing to their potential for property tunability, multifunctionality, and property coupling. As a new group of self‐assembled thin films, vertically aligned nanocomposite (VAN)‐based hybrid metamaterials have been demonstrated with significant anisotropic physical properties and a broad range of property tailorability, such as optical anisotropy, magnetic anisotropy, hyperbolic dispersion, and enhanced second harmonic generation properties. Herein, self‐assembled ZrO 2 ‐Co nanocomposite films, with high epitaxial quality and ultra‐fine vertically aligned Co nanopillars (with an average diameter of ≈2 nm) embedded in a ZrO 2 matrix, are fabricated using a pulsed laser deposition (PLD) method. The Co pillar density can be effectively tuned by varying the Co concentration in the target, which results in tunable optical properties and magnetic properties. Specifically, a high saturation magnetization of 100 emu cm −3 , strong out‐of‐plane magnetic anisotropy and tailorable magnetization properties are achieved via tuning the Co nanopillar density. Coupled with hyperbolic dispersion of dielectric constant from 950 to 1500 nm in wavelength, plasmonic Co metal nanopillars, and the unique dielectric ZrO 2 matrix, this new nanoscale hybrid metamaterial shows great potential for future integrated optical and magnetic device designs.

36 MATERIALS SCIENCE↗