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 415 records · Page 23

Benchmark calculations on the lowest-energy singlet, triplet, and quintet states of the four-electron harmonium atom

For a wide range of confinement strengths ω, explicitly-correlated calculations afford approximate energies E(ω) of the ground and low-lying excited states of the four-electron harmonium atom that are within few μhartree of the exact values, the errors in the respective energy components being only slightly higher. This level of accuracy constitutes an improvement of several orders of magnitude over the previously published data, establishing a set of benchmarks for stringent calibration and testing of approximate electronic structure methods. Its usefulness is further enhanced by the construction of differentiable approximants that allow for accurate computation of E(ω) and its components for arbitrary values of ω. The diversity of the electronic states in question, which involve both single- and multideterminantal first-order wavefunctions, and the availability of the relevant natural spinorbitals and their occupation numbers make the present results particularly useful in research on approximate density-matrix functionals. The four-electron harmonium atom is found to possess the {sup 3}P{sub +} triplet ground state at strong confinements and the {sup 5}S{sub −} quintet ground state at the weak ones, the energy crossing occurring at ω ≈ 0.0240919.

74 ATOMIC AND MOLECULAR PHYSICS↗

Quantum tunneling splittings from path-integral molecular dynamics

We illustrate how path-integral molecular dynamics can be used to calculate ground-state tunnelling splittings in molecules or clusters. The method obtains the splittings from ratios of density matrix elements between the degenerate wells connected by the tunnelling. We propose a simple thermodynamic integration scheme for evaluating these elements. Numerical tests on fully dimensional malonaldehyde yield tunnelling splittings in good overall agreement with the results of diffusion Monte Carlo calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of the finite linewidth of the laser radiation spectrum on the shape of the coherent population trapping resonance line in an optically dense medium with a buffer gas

The theory of coherent population trapping resonance is developed for the finite linewidth of the laser radiation spectrum in an optically dense medium of Λ atoms in a cell with a buffer gas. Equations are derived for the atomic density matrix and laser emission spectrum transfer in a cell with working and buffer gases at a finite temperature. The dependence of the quality factor of coherent population trapping resonance on the linewidth of the laser radiation spectrum is studied by measuring transmitted radiation and fluorescence signals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum and classical correlations in electron-nuclear spin echo

The quantum properties of dynamic correlations in a system of an electron spin surrounded by nuclear spins under the conditions of free induction decay and electron spin echo have been studied. Analytical results for the time evolution of mutual information, classical part of correlations, and quantum part characterized by quantum discord have been obtained within the central-spin model in the high-temperature approximation. The same formulas describe discord in both free induction decay and spin echo although the time and magnetic field dependences are different because of difference in the parameters entering into the formulas. Changes in discord in the presence of the nuclear polarization β{sub I} in addition to the electron polarization β{sub S} have been calculated. It has been shown that the method of reduction of the density matrix to a two-spin electron-nuclear system provides a qualitatively correct description of pair correlations playing the main role at β{sub S} ≈ β{sub I} and small times. At large times, such correlations decay and multispin correlations ensuring nonzero mutual information and zero quantum discord become dominant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

To the Theory of Inhomogeneous Electron Gas

Nonrigorous character of the density-functional theory for inhomogeneous electron gas based on the hypothesis assuming the existence of a universal density functional is demonstrated. A single-particle density matrix must be determined to calculate the ground-state energy of a finite system with a finite number of electrons. A single-particle Green function can be used to unambiguously determine the ground-state energy of an inhomogeneous electron system that satisfies the thermodynamic limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Charge order and superconductivity in a two-band model for infinite-layer nickelates

The recent discovery of superconductivity in infinite-layer nickelates has drawn considerable attention; however, a consensus on the fundamental building blocks and common ingredients necessary to understand and describe their ground states and emergent properties is lacking. A series of experimental and theoretical studies have suggested that an effective two-band Hubbard model with Ni 3d x 2 -y 2 and rare-earth (R) 5d character may describe the low-energy physics. In this report we study the ground state properties of this two-band model on four-leg cylinders using the density-matrix renormalization group (DMRG) technique to better grasp whether such a simple model can embody the essential physics. A key difference compared to single-band Hubbard materials is that the system is self-doped: even at overall half-filling, the R-band acts as an electron reservoir, hole-doping the Ni-layer and fundamentally altering the physics expected from an undoped antiferromagnet. On the four-leg cylinder, the ground state is consistent with a Luttinger liquid, with anti-phase modulations of the charge density in the Ni- and R-layers having corresponding wavevectors that lock together. Light hole doping away from 1/2 filling releases the locking between the Ni and the R charge modulations, as the electron density in the R-band decreases and eventually becomes exhausted at a hole doping concentration that depends sensitively on the effective splitting between the Ni and the R orbitals. The ground state of the doped system is consistent with a Luther-Emery liquid, possessing quasi-long-range superconducting correlations in the Ni layer, similar to the single-band Hubbard model. Our results are consistent with experimental observations and may help to reveal the microscopic mechanism for pairing and other emergent properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Moduli Space of Wigner Quasiprobability Distributions for N-Dimensional Quantum Systems

A mapping between operators on the Hilbert space of an N-dimensional quantum system and Wigner quasiprobability distributions defined on the symplectic flag manifold is discussed. The Wigner quasiprobability distribution is constructed as a dual pairing between the density matrix and the Stratonovich–Weyl kernel. It is shown that the moduli space of Stratonovich–Weyl kernels is given by the intersection of the coadjoint orbit space of the group SU(N) and a unit (N − 2)-dimensional sphere. The general considerations are exemplified by a detailed description of the moduli space of 2, 3, and 4-dimensional systems.

97 MATHEMATICS AND COMPUTING↗

Donor impurity states and related terahertz range nonlinear optical response in GaN cylindrical quantum wires: Effects of external electric and magnetic fields

We report a study on the optical absorption coefficient associated to hydrogenic impurity interstate transitions in zinc-blende GaN quantum wires of cylindrical shape taking into account the effects of externally applied static electric and magnetic fields. The electron states emerge within the effective mass approximation, via the exact diagonalization of the donor-impurity Hamiltonian with parabolic confinement and external field effects. The nonlinear optical absorption is calculated using a recently derived expression for the dielectric susceptibility, obtained via a nonperturbative solution of the density-matrix Bloch equation. Our results show that this treatment eliminates not only the intensity-dependent bleaching effect but also the change in sign of the nonlinear contribution due to the combined effect of asymmetric impurity location and the applied electric field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Estimating the entropy and quantifying the impurity of a swarm of surface-hopping trajectories: A new perspective on decoherence

In this article, we consider the intrinsic entropy of Tully's fewest switches surface hopping (FSSH) algorithm (as estimated by the impurity of the density matrix) [J. Chem. Phys. 93, 1061 (1990)]. We show that, even for a closed system, the total impurity of a FSSH calculation increases in time (rather than stays constant). This apparent failure of the FSSH algorithm can be traced back to an incorrect, approximate treatment of the electronic coherence between wavepackets moving along different potential energy surfaces. This incorrect treatment of electronic coherence also prevents the FSSH algorithm from correctly describing wavepacket recoherences (which is a well established limitation of the FSSH method). Nevertheless, despite these limitations, the FSSH algorithm often predicts accurate observables because the electronic coherence density is modulated by a phase factor which varies rapidly in phase space and which often integrates to almost zero. Adding “decoherence” events on top of a FSSH calculation completely destroys the incorrect FSSH electronic coherence and effectively sets the Poincaré recurrence time for wavepacket recoherence to infinity; this modification usually increases FSSH accuracy (assuming there are no recoherences) while also offering long-time stability for trajectories. In practice, we show that introducing “decoherence” events does not change the total FSSH impurity significantly, but does lead to more accurate evaluations of the impurity of the electronic subsystem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-consistent determination of the many-body state of ultracold bosonic atoms in a one-dimensional harmonic trap

We study zero-temperature quantum fluctuations in harmonically trapped one-dimensional interacting Bose gases,using the self-consistent multiconfigurational time-dependent Hartree method. We define phase fluctuations from the full single-particle density matrix by the spatial decay exponent of off-diagonal long-range order. In a regime of mesoscopic particle numbers and moderate contact couplings, we derive the spatial dependence of the amplitude of phase fluctuations, determined from the self-consistently derived shape of the field operator orbitals and Fock space orbital occupation amplitudes. It is shown that the phase fluctuations display a peak, which in turn corresponds to a dip of the first-order correlations in position space, akin to what has previously been obtained in the Tonks–Girardeau limit of very large interactions and low densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Derivation of general analytic gradient expressions for density-fitted post-Hartree-Fock methods: An efficient implementation for the density-fitted second-order Møller–Plesset perturbation theory

General analytic gradient expressions (with the frozen-core approximation) are presented for density-fitted post-HF methods. An efficient implementation of frozen-core analytic gradients for the second-order Møller–Plesset perturbation theory (MP2) with the density-fitting (DF) approximation (applying to both reference and correlation energies), which is denoted as DF-MP2, is reported. The DF-MP2 method is applied to a set of alkanes, conjugated dienes, and noncovalent interaction complexes to compare the computational cost of single point analytic gradients with MP2 with the resolution of the identity approach (RI-MP2) [F. Weigend and M. Häser, Theor. Chem. Acc. 97, 331 (1997); R. A. Distasio, R. P. Steele, Y. M. Rhee, Y. Shao, and M. Head-Gordon, J. Comput. Chem. 28, 839 (2007)]. In the RI-MP2 method, the DF approach is used only for the correlation energy. Our results demonstrate that the DF-MP2 method substantially accelerate the RI-MP2 method for analytic gradient computations due to the reduced input/output (I/O) time. Because in the DF-MP2 method the DF approach is used for both reference and correlation energies, the storage of 4-index electron repulsion integrals (ERIs) are avoided, 3-index ERI tensors are employed instead. Further, as in case of integrals, our gradient equation is completely avoid construction or storage of the 4-index two-particle density matrix (TPDM), instead we use 2- and 3-index TPDMs. Hence, the I/O bottleneck of a gradient computation is significantly overcome. Therefore, the cost of the generalized-Fock matrix (GFM), TPDM, solution of Z-vector equations, the back transformation of TPDM, and integral derivatives are substantially reduced when the DF approach is used for the entire energy expression. Further application results show that the DF approach introduce negligible errors for closed-shell reaction energies and equilibrium bond lengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SS{sub p}G: A strongly orthogonal geminal method with relaxed strong orthogonality

Strong orthogonality is an important constraint placed on geminal wavefunctions in order to make variational minimization tractable. However, strong orthogonality prevents certain, possibly important, excited configurations from contributing to the ground state description of chemical systems. The presented method lifts strong orthogonality constraint from geminal wavefunction by computing a perturbative-like correction to each geminal independently from the corrections to all other geminals. The method is applied to the Singlet-type Strongly orthogonal Geminals variant of the geminal wavefunction. Comparisons of this new SS{sub p}G method are made to the non-orthogonal AP1roG and the unconstrained Geminal Mean-Field Configuration Interaction method using small atomic and molecular systems. The correction is also compared to Density Matrix Renormalization Group calculations performed on long polyene chains in order to assess its scalability and applicability to large strongly correlated systems. The results of these comparisons demonstrate that although the perturbative correction is small, it may be a necessary first step in the systematic improvement of any strongly orthogonal geminal method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Numerically Exact Approach to Quantum Impurity Problems in Realistic Lattice Geometries

We review a controlled numerical approach to quantum impurity problems in realistic geometries, consisting of exactly mapping the complete lattice Hamiltonian onto an equivalent one dimensional system through a unitary transformation. The resulting dimensional and entanglement reduction allows one to study the quantum many-body problem on arbitrary d-dimensional lattices using the density matrix renormalization group (DMRG) method. The real-space resolution allows one to position the impurities at the boundary or bulk of the sample and to study screening effects due to edge or surface modes. We describe how to generalize this approach to multi-impurity problems, discuss applications and possible extensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Intertwined charge, spin, and pairing orders in doped iron ladders

Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe2S3. The model is constructed by using ab initio hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the (π,0) magnetic ordering at half filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard U coupling, other magnetic patterns, such as (0,π), are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intraorbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge density wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such a scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Magnetoelectric effects in the spin-1/2 XXZ model with Dzyaloshinskii-Moriya interaction

We study the 1D spin-1/2 XXZ chain in the presence of the Dzyaloshinskii-Moriya (D-M) interaction and with longitudinal and transverse magnetic fields. We assume the spin-current mechanism of Katsura-Nagaosa-Balatsky at play and interpret the D-M interaction as a coupling between the local electric polarization and an external electric field. We study the interplay of electric and magnetic order in the ground state using the numerical density matrix renormalization group(DMRG) method. Specifically, we investigate the dependences of the magnetization and electric polarization on the external electric and magnetic fields. We find that for transverse magnetic fields, there are two different regimes of polarization while for longitudinal magnetic fields, there are three different regimes of polarization. The different regimes can be tuned by the external magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-Hermitian wave packet approximation for coupled two-level systems in weak and intense fields

We introduce a non-Hermitian Schrödinger-type approximation of optical Bloch equations for two-level systems. This approximation provides a complete and accurate description of the coherence and decoherence dynamics in both weak and strong laser fields at the cost of losing accuracy in the description of populations. In this approach, it is sufficient to propagate the wave function of the quantum system instead of the density matrix, providing that relaxation and dephasing are taken into account via automatically adjusted time-dependent gain and decay rates. The developed formalism is applied to the problem of scattering and absorption of electromagnetic radiation by a thin layer comprised of interacting two-level emitters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Angular correlations in the inelastic-scattering process {sup 24}Mg(d, d{gamma}){sup 24}Mg at E{sub d} = 15.3 MeV

The double-differential cross sections for the reaction {sup 24}Mg(d, d{gamma}{sub 1.369}){sup 24}Mg at the projectiledeuteron energy of E{sub d} = 15.3 MeV weremeasured for deuteron emission angles in the forward hemisphere. All even spin-tensor components of the density matrix for the 2{sup +} state of the {sup 24}Mg nucleus at 1.369 MeV were reconstructed, and its orientation properties were determined. These experimental results were compared with the results of calculations based on various versions of the coupled-channel method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Communication: CDFT-CI couplings can be unreliable when there is fractional charge transfer

Constrained density functional theory with configuration interaction (CDFT-CI) is a useful, low-cost tool for the computational prediction of electronic couplings between pseudo-diabatic constrained electronic states. Such couplings are of paramount importance in electron transfer theory and transition state theory, among other areas of chemistry. Unfortunately, CDFT-CI occasionally fails significantly, predicting a coupling that does not decay exponentially with distance and/or overestimating the expected coupling by an order of magnitude or more. In this communication, we show that the eigenvalues of the difference density matrix between the two constrained states can be used as an a priori metric to determine when CDFT-CI are likely to be reliable: when the eigenvalues are near 0 or ±1, transfer of a whole electron is occurring, and CDFT-CI can be trusted. We demonstrate the utility of this metric with several illustrative examples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗