Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Here we investigate twisted bilayer graphene (TBG) at filling ν = -3 in the presence of realistic heterostrain. Strain amplifies the band dispersion and drives the system beyond the strong-coupling regime of previous theoretical studies. We use DMRG to conduct an unbiased, large-scale numerical calculations that include all spin and valley degrees of freedom, up to bond dimension χ = 24576. We establish a global phase diagram that unifies a number of theoretical and experimental results. Near zero strain we find an intervalley-coherent quantized anomalous Hall (QAH-IVC) state, a competitive strong-coupling order that evaded past numerical studies. A tiny strain around 0.05% drives a transition into an incommensurate Kekulé spiral (IKS) phase, supporting the mean-field prediction in [Kwan, Phys. Rev. X 11, 041063 (2021)2160-330810.1103/PhysRevX.11.041063]. Even higher strains above 0.2% favor a flavor-symmetric metallic order, which may explain metals found at ν = -3 in many experiments.
Not provided.
Not Available
Not Available
Not Available
The overarching goal of this research program is to quantitatively predict out-of-equilibrium quantum dynamics and transport with electrons and spin degrees of freedom coupled to phonons, photons, defects and liquid environments, specifically targeting spin-selective chemistry. Spin polarized carriers, which can be generated at room temperature by circularly-polarized light, provide unique opportunities for creating new reaction pathways and increasing selectivity towards specific pathways. If realized, this can substantially increase the level of reaction control possible in the field of heterogeneous catalysis. Spin, a quantum mechanical entity without a classical counterpart, requires quantum dynamics. Capturing both coherent and incoherent dynamics and accurate correlation of spins, electrons and phonon environments are essential to quantitative description of generation, relaxation and transport, as well as chemical reaction of spin-polarized carriers in molecules and materials promising for spin-photocatalysis.
The effect of correlations between an absorber atom and perturbers in the binary-collision approximation are applied to degenerate atomic systems. A generalized absorption profile which specifies the final state of the atom after an absorption event is related to the total intensities of Rayleigh scattering and fluorescence from the atom. It is suggested that additional dynamical information to that obtainable from ordinary absorption experiments is required in order to describe redistributed atomic radiation. The scattering of monochromatic radiation by a degenerate atom is computed in a binary-collision approximation; an equation of motion is derived for the correlation function which is valid outside the quantum-regression regime. Solutions are given for the weak-field conditions in terms of generalized absorption and emission profiles that depend on the indices of the atomic multipoles.
The susceptibilities for a three-level system with arbitrary pump and signal field strengths are derived for arbitrary longitudinal and transverse relaxation rates. The results are of interest in connection with the calculation of the Raman gain in systems where resonance enhancement plays a dominant role.
The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.
Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.
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.
The region of validity of conventional rate equation formulations used to describe the excitation of flame gases for spectroscopic diagnostic purposes is examined by direct comparison with the density matrix formulation. The density matrix equations are formulated for an ensemble of atoms or molecules in the gas state. These equations are reduced to the steady-state limit for a two-level system under the effect of single-mode and multimode optical fields. The steady-state equations are shown to be valid only when the laser pulse rises slowly with respect to the collisional dephasing rate for single-mode excitation, and for certain special cases under multimode excitation.
We study the time evolution of the density matrix of a high energy quark in the presence of a dense QCD background that is modeled as a stochastic Gaussian color field. At late times, we find that only the color singlet component of the quark’s reduced density matrix survives the in-medium evolution and that the density matrix becomes asymptotically diagonal in both transverse position and momentum spaces. In addition, we observe an accelerated entropy growth due to the larger phase space being explored by the quark and that the quantum and classical quark entropies converge at late times. We further observe that the quark state loses all memory of the initial condition. Combined with the fact that the reduced density matrix satisfies Boltzmann-diffusion transport, we conclude that the quark reduced density matrix can be interpreted as a classical phase space distribution.
Simulating the nonequilibrium dynamics of thermal states is a fundamental problem across scales from high-energy to condensed-matter physics. Quantum computers may provide a way to solve this problem efficiently. Preparing a thermal state on a quantum computer is challenging, but there exist methods to circumvent this by computing a weighted sum of time-dependent matrix elements in a convenient basis. Further, while the number of basis states can be large, in this paper we show that it can be reduced by simulating only the largest density matrix elements by weight, capturing the density matrix to a specified precision. Leveraging Hamiltonian symmetries enables further reductions. This approach paves the way to more accurate thermal-state dynamics simulations on near-term quantum hardware.