skip to main content


Title: Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

 
more » « less
Award ID(s):
1931258
NSF-PAR ID:
10440032
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Institute of Physics
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
157
Issue:
19
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The azimuthal ($$\Delta \varphi $$Δφ) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p–Pb collisions at$$\sqrt{s_{\mathrm{{NN}}}} = 5.02$$sNN=5.02TeV. Results are reported for electrons with transverse momentum$$44<pT<16$$\textrm{GeV}/c$$GeV/c and pseudorapidity$$|\eta |<0.6$$|η|<0.6. The associated charged particles are selected with transverse momentum$$11<pT<7$$\textrm{GeV}/c$$GeV/c, and relative pseudorapidity separation with the leading electron$$|\Delta \eta | < 1$$|Δη|<1. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p–Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The$$\Delta \varphi $$Δφdistribution and the peak observables in pp and p–Pb collisions are compared with calculations from various Monte Carlo event generators.

     
    more » « less
  2. In a recent paper (JCTC, 16, 6098 (2020)), we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces, called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excitonic character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons (PAH) with respect to the extrapolated result (i.e. near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 biexcitonic states. 
    more » « less
  3. Abstract Purpose

    To develop an independent log file–based intensity‐modulated radiation therapy (IMRT) quality assurance (QA) tool for the 0.35 T magnetic resonance‐linac (MR‐linac) and investigate the ability of various IMRT plan complexity metrics to predict the QA results. Complexity metrics related to tissue heterogeneity were also introduced.

    Methods

    The tool for particle simulation (TOPAS) Monte Carlo code was utilized with a previously validated linac head model. A cohort of 29 treatment plans was selected for IMRT QA using the developed QA tool and the vendor‐supplied adaptive QA (AQA) tool. For 27 independent patient cases, various IMRT plan complexity metrics were calculated to assess the deliverability of these plans. A correlation between the gamma pass rates (GPRs) from the AQA results and calculated IMRT complexity metrics was determined using the Pearson correlation coefficients. Tissue heterogeneity complexity metrics were calculated based on the gradient of the Hounsfield units.

    Results

    The median and interquartile range for the TOPAS GPRs (3%/3 mm criteria) were 97.24% and 3.75%, respectively, and were 99.54% and 0.36% for the AQA tool, respectively. The computational time for TOPAS ranged from 4 to 8 h to achieve a statistical uncertainty of <1.5%, whereas the AQA tool had an average calculation time of a few minutes. Of the 23 calculated IMRT plan complexity metrics, the AQA GPRs had correlations with 7 out of 23 of the calculated metrics. Strong correlations (|r| > 0.7) were found between the GPRs and the heterogeneity complexity metrics introduced in this work.

    Conclusions

    An independent MC and log file–based IMRT QA tool was successfully developed and can be clinically deployed for offline QA. The complexity metrics will supplement QA reports and provide information regarding plan complexity.

     
    more » « less
  4. This work supports the existence of extended nonergodic states in the intermediate region between the chaotic (thermal) and the many‐body localized phases. These states are identified through an extensive analysis of static and dynamical properties of a finite one‐dimensional system with onsite random disorder. The long‐time dynamics is particularly sensitive to changes in the spectrum and in the structures of the eigenstates. The study of the evolution of the survival probability, Shannon information entropy, and von Neumann entanglement entropy enables the distinction between the chaotic and the intermediate region.image

     
    more » « less
  5. Trace inequalities are general techniques with many applications in quantum information theory, often replacing the classical functional calculus in noncommutative settings. The physics of quantum field theory and holography, however, motivates entropy inequalities in type III von Neumann algebras that lack a semifinite trace. The Haagerup and Kosaki Lp spaces enable re-expressing trace inequalities in non-tracial von Neumann algebras. In particular, we show this for the generalized Araki–Lieb–Thirring and Golden–Thompson inequalities from the work of Sutter et al. [Commun. Math. Phys. 352(1), 37 (2017)]. Then, using the Haagerup approximation method, we prove a general von Neumann algebra version of universal recovery map corrections to the data processing inequality for relative entropy. We also show subharmonicity of a logarithmic p-fidelity of recovery. Furthermore, we prove that the non-decrease of relative entropy is equivalent to the existence of an L1-isometry implementing the channel on both input states. 
    more » « less