Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Matrix”

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 1,099 records · Page 61

How many invariant polynomials are needed to decide local unitary equivalence of qubit states?

Given L-qubit states with the fixed spectra of reduced one-qubit density matrices, we find a formula for the minimal number of invariant polynomials needed for solving local unitary (LU) equivalence problem, that is, problem of deciding if two states can be connected by local unitary operations. Interestingly, this number is not the same for every collection of the spectra. Some spectra require less polynomials to solve LU equivalence problem than others. The result is obtained using geometric methods, i.e., by calculating the dimensions of reduced spaces, stemming from the symplectic reduction procedure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Communication: Reduced density matrices in molecular systems: Grand-canonical electron states

Grand-canonical like descriptions of many electron atomic and molecular open systems which are characterized by a non-integer number of electrons are presented. Their associated reduced density matrices (RDMs) are obtained by introducing the contracting mapping for this type of distributions. It is shown that there is loss of information when connecting RDMs of different order by partial contractions. The energy convexity property of these systems simplifies the description. Consequently, this formulation opens the possibility to a new look for chemical descriptors such as chemical potential and reactivity among others. Examples are presented to discuss the theoretical aspects of this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Generalization of Lieb's variational principle to Bogoliubov–Hartree–Fock theory

In its original formulation, Lieb's variational principle holds for fermion systems with purely repulsive pair interactions. As a generalization we prove for both fermion and boson systems with semi-bounded Hamiltonian that the infimum of the energy over quasifree states coincides with the infimum over pure quasifree states. In particular, the Hamiltonian is not assumed to preserve the number of particles. To shed light on the relation between our result and the usual formulation of Lieb's variational principle in terms of one-particle density matrices, we also include a characterization of pure quasifree states by means of their generalized one-particle density matrices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How accurate is the strongly orthogonal geminal theory in predicting excitation energies? Comparison of the extended random phase approximation and the linear response theory approaches

Performance of the antisymmetrized product of strongly orthogonal geminal (APSG) ansatz in describing ground states of molecules has been extensively explored in the recent years. Not much is known, however, about possibilities of obtaining excitation energies from methods that would rely on the APSG ansatz. In the paper we investigate the recently proposed extended random phase approximations, ERPA and ERPA2, that employ APSG reduced density matrices. We also propose a time-dependent linear response APSG method (TD-APSG). Its relation to the recently proposed phase including natural orbital theory is elucidated. The methods are applied to Li{sub 2}, BH, H{sub 2}O, and CH{sub 2}O molecules at equilibrium geometries and in the dissociating limits. It is shown that ERPA2 and TD-APSG perform better in describing double excitations than ERPA due to inclusion of the so-called diagonal double elements. Analysis of the potential energy curves of Li{sub 2}, BH, and H{sub 2}O reveals that ERPA2 and TD-APSG describe correctly excitation energies of dissociating molecules if orbitals involved in breaking bonds are involved. For single excitations of molecules at equilibrium geometries the accuracy of the APSG-based methods approaches that of the time-dependent Hartree-Fock method with the increase of the system size. A possibility of improving the accuracy of the TD-APSG method for single excitations by splitting the electron-electron interaction operator into the long- and short-range terms and employing density functionals to treat the latter is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

Exciton delocalization, charge transfer, and electronic coupling for singlet excitation energy transfer between stacked nucleobases in DNA: An MS-CASPT2 study

Exciton delocalization and singlet excitation energy transfer have been systematically studied for the complete set of 16 DNA nucleobase dimers in their ideal, single-strand stacked B-DNA conformation, at the MS-CASPT2 level of theory. The extent of exciton delocalization in the two lowest (π,π{sup *}) states of the dimers is determined using the symmetrized one-electron transition density matrices between the ground and excited states, and the electronic coupling is calculated using the delocalization measure and the energy splitting between the states [see F. Plasser, A. J. A. Aquino, W. L. Hase, and H. Lischka, J. Phys. Chem. A 116, 11151–11160 (2012)]. The calculated couplings lie between 0.05 eV and 0.14 eV. In the B-DNA conformation, where the interchromophoric distance is 3.38 Å, our couplings deviate significantly from those calculated with the transition charges, showing the importance of orbital overlap components for the couplings in this conformation. The calculation of the couplings is based on a two-state model for exciton delocalization. However, in three stacks with a purine in the 5{sup ′} position and a pyrimidine in the 3{sup ′} one (AT, GC, and GT), there is an energetically favored charge transfer state that mixes with the two lowest excited states. In these dimers we have applied a three-state model that considers the two locally excited diabatic states and the charge transfer state. Using the delocalization and charge transfer descriptors, we obtain all couplings between these three states. Our results are important in the context of DNA photophysics, since the calculated couplings can be used to parametrize effective Hamiltonians to model extended DNA stacks. Our calculations also suggest that the 5{sup ′}-purine-pyrimidine-3{sup ′} sequence favors the formation of charge transfer excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗

Conditions for uniqueness of product representations for separable quantum channels and separable quantum states

We give a sufficient condition that an operator sum representation of a separable quantum channel in terms of product operators is the unique product representation for that channel, and then provide examples of such channels for any number of parties. This result has implications for efforts to determine whether or not a given separable channel can be exactly implemented by local operations and classical communication. By the Choi-Jamiolkowski isomorphism, it also translates to a condition for the uniqueness of product state ensembles representing a given quantum state. These ideas follow from considerations concerning whether or not a subspace spanned by a given set of product operators contains at least one additional product operator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Resource cost results for one-way entanglement distillation and state merging of compound and arbitrarily varying quantum sources

We consider one-way quantum state merging and entanglement distillation under compound and arbitrarily varying source models. Regarding quantum compound sources, where the source is memoryless, but the source state an unknown member of a certain set of density matrices, we continue investigations begun in the work of Bjelaković et al. [“Universal quantum state merging,” J. Math. Phys. 54, 032204 (2013)] and determine the classical as well as entanglement cost of state merging. We further investigate quantum state merging and entanglement distillation protocols for arbitrarily varying quantum sources (AVQS). In the AVQS model, the source state is assumed to vary in an arbitrary manner for each source output due to environmental fluctuations or adversarial manipulation. We determine the one-way entanglement distillation capacity for AVQS, where we invoke the famous robustification and elimination techniques introduced by Ahlswede. Regarding quantum state merging for AVQS we show by example that the robustification and elimination based approach generally leads to suboptimal entanglement as well as classical communication rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unbiased reduced density matrices and electronic properties from full configuration interaction quantum Monte Carlo

Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the full configuration interaction quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself) and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Semi-stochastic full configuration interaction quantum Monte Carlo: Developments and application

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The multi-configurational time-dependent Hartree approach revisited

The multi-configurational time-dependent Hartree (MCTDH) approach facilitates accurate high-dimensional quantum dynamics simulations. In the approach, the wavefunction is expanded in a direct product of self-adapting time-dependent single-particle functions (SPFs). The equations of motion for the expansion coefficients and the SPFs are obtained via the Dirac-Frenkel variational principle. While this derivation yields well-defined differential equations for the motion of occupied SPFs, singularities in the working equations resulting from unoccupied SPFs have to be removed by a regularization procedure. Here, an alternative derivation of the MCTDH equations of motion is presented. It employs an analysis of the time-dependence of the single-particle density matrices up to second order. While the analysis of the first order terms yields the known equations of motion for the occupied SPFs, the analysis of the second order terms provides new equations which allow one to identify optimal choices for the unoccupied SPFs. The effect of the optimal choice of the unoccupied SPFs on the structure of the MCTDH equations of motion and their regularization is discussed. Generalized equations applicable in the multi-layer MCTDH framework are presented. Finally, the effects resulting from the initial choice of the unoccupied SPFs are illustrated by a simple numerical example.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent second-order Green’s function perturbation theory for periodic systems

Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green’s function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green’s function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice — a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction

Over the past few years, pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled-cluster based basis sets for valence correlation calculations

Novel basis sets are generated that target the description of valence correlation in atoms H through Ar. The new contraction coefficients are obtained according to the Atomic Natural Orbital (ANO) procedure from CCSD(T) (coupled-cluster singles and doubles with perturbative triples correction) density matrices starting from the primitive functions of Dunning et al. [J. Chem. Phys. 90, 1007 (1989); ibid. 98, 1358 (1993); ibid. 100, 2975 (1993)] (correlation consistent polarized valence X-tuple zeta, cc-pVXZ). The exponents of the primitive Gaussian functions are subject to uniform scaling in order to ensure satisfaction of the virial theorem for the corresponding atoms. These new sets, named ANO-VT-XZ (Atomic Natural Orbital Virial Theorem X-tuple Zeta), have the same number of contracted functions as their cc-pVXZ counterparts in each subshell. The performance of these basis sets is assessed by the evaluation of the contraction errors in four distinct computations: correlation energies in atoms, probing the density in different regions of space via 〈r{sup n}〉 (−3 ≤ n ≤ 3) in atoms, correlation energies in diatomic molecules, and the quality of fitting potential energy curves as measured by spectroscopic constants. All energy calculations with ANO-VT-QZ have contraction errors within “chemical accuracy” of 1 kcal/mol, which is not true for cc-pVQZ, suggesting some improvement compared to the correlation consistent series of Dunning and co-workers.

74 ATOMIC AND MOLECULAR PHYSICS↗

An integral-factorized implementation of the driven similarity renormalization group second-order multireference perturbation theory

We report an efficient implementation of a second-order multireference perturbation theory based on the driven similarity renormalization group (DSRG-MRPT2) [C. Li and F. A. Evangelista, J. Chem. Theory Comput. 11, 2097 (2015)]. Our implementation employs factorized two-electron integrals to avoid storage of large four-index intermediates. It also exploits the block structure of the reference density matrices to reduce the computational cost to that of second-order Møller–Plesset perturbation theory. Our new DSRG-MRPT2 implementation is benchmarked on ten naphthyne isomers using basis sets up to quintuple-ζ quality. We find that the singlet-triplet splittings (Δ{sub ST}) of the naphthyne isomers strongly depend on the equilibrium structures. For a consistent set of geometries, the Δ{sub ST} values predicted by the DSRG-MRPT2 are in good agreements with those computed by the reduced multireference coupled cluster theory with singles, doubles, and perturbative triples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Local unitary equivalence of quantum states and simultaneous orthogonal equivalence

The correspondence between local unitary equivalence of bipartite quantum states and simultaneous orthogonal equivalence is thoroughly investigated and strengthened. It is proved that local unitary equivalence can be studied through simultaneous similarity under projective orthogonal transformations, and four parametrization independent algorithms are proposed to judge when two density matrices on ℂ{sup d{sub 1}} ⊗ ℂ{sup d{sub 2}} are locally unitary equivalent in connection with trace identities, Kronecker pencils, Albert determinants and Smith normal forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Distillation of secret-key from a class of compound memoryless quantum sources

We consider secret-key distillation from tripartite compound classical-quantum-quantum (cqq) sources with free forward public communication under strong security criterion. We design protocols which are universally reliable and secure in this scenario. These are shown to achieve asymptotically optimal rates as long as a certain regularity condition is fulfilled by the set of its generating density matrices. We derive a multi-letter formula which describes the optimal forward secret-key capacity for all compound cqq sources being regular in this sense. We also determine the forward secret-key distillation capacity for situations where the legitimate sending party has perfect knowledge of his/her marginal state deriving from the source statistics. In this case regularity conditions can be dropped. Our results show that the capacities with and without the mentioned kind of state knowledge are equal as long as the source is generated by a regular set of density matrices. We demonstrate that regularity of cqq sources is not only a technical but also an operational issue. For this reason, we give an example of a source which has zero secret-key distillation capacity without sender knowledge, while achieving positive rates is possible if sender marginal knowledge is provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamically weighted multireference perturbation theory: Combining the advantages of multi-state and state-averaged methods

In this paper, we introduce two new approaches to compute near-degenerate electronic states based on the driven similarity renormalization group (DSRG) framework. The first approach is a unitary multi-state formalism based on the DSRG (MS-DSRG), whereby an effective Hamiltonian is built from a set of state-specific solutions. The second approach employs a dynamic weighting parameter to smoothly interpolate between the multi-state and the state-averaged DSRG schemes. The resulting dynamically weighted DSRG (DW-DSRG) theory incorporates the most desirable features of both multi-state approaches (ability to accurately treat many states) and state-averaged methods (correct description of avoided crossings and conical intersections). We formulate second-order perturbation theories (PT2) based on the MS- and DW-DSRG and study the potential energy curves of LiF, the conical intersection of the two lowest singlet states of NH 3 , and several low-lying excited states of benzene, naphthalene, and anthracene. The DW-DSRG-PT2 predicts the correct avoided crossing of LiF and avoids artifacts produced by the corresponding state-specific and multi-state theories. Excitation energies of the acenes computed with the DW-DSRG-PT2 are found to be more accurate than the corresponding state-averaged values, showing a small dependence on the number of states computed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fourth order expressions for the electronic absorption lineshape of molecular excitons

The line shape of electronic absorption spectroscopy reflects the information on quantum dynamical processes accompanying the electronic excitation, and its accurate description is an important component for validating theoretical models and assumptions. The present work provides detailed expressions for the absorption line shape of molecular excitons that are valid up to the fourth order of exciton-bath interactions within the quantum master equation approach. These expressions can serve as the basis for developing general and systematic methods to model the line shape for a broad class of molecular exciton systems and environments. For the bath model of linearly coupled harmonic oscillators, more detailed expressions employing the spectral densities of the bath are presented. These expressions are then tested for a linear aggregate of identical chromophores each coupled to the harmonic oscillator bath. Finally, calculation results for a super-Ohmic spectral density with exponential cutoff demonstrate the feasibility of calculations and also offer insights into the utility and difficulty of going beyond the second order approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗