Source author record

Francesco A. Evangelista

Francesco A. Evangelista appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Automatic derivation of fermionic many-body theories based on general Fermi vacua

This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory (DSRG-MRPT2).

preprint2022arXiv

Leveraging small scale quantum computers with unitarily downfolded Hamiltonians

In this work, we propose a quantum unitary downfolding formalism based on the driven similarity renormalization group (QDSRG) that may be combined with quantum algorithms for both noisy and fault-tolerant hardware. The QDSRG is a classical polynomially-scaling downfolding method that avoids the evaluation of costly three- and higher-body reduced density matrices while retaining the accuracy of classical multireference many-body theories. We calibrate and test the QDSRG on several challenging chemical problems and propose a strategy for avoiding classical exponential-scaling steps in the QDSRG scheme. We report QDSRG computations of two chemical systems using the variational quantum eigensolver on IBM quantum devices: i) the dissociation curve of H$_2$ using a quintuple-$ζ$ basis and ii) the bicyclobutane isomerization reaction to $trans$-butadiene, demonstrating the reduction of problems that require several hundred qubits to a single qubit. Our work shows that the QDSRG is a viable approach to leverage near-term quantum devices for the accurate estimation of molecular properties.

preprint2020arXiv

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing a representative example on the driven similarity renormalization group (DSRG). Starting from the least approximate method - the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)] - we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-$n$ and -$n'$ ($n=1,2,3,4$) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive quadratic commutators [qDSRG(2)]. The resulting qDSRG(2)+(T) approach matches the accuracy of the "gold-standard" coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

preprint2020arXiv

Correlating AGP on a quantum computer

For variational algorithms on the near term quantum computing hardware, it is highly desirable to use very accurate ansatze with low implementation cost. Recent studies have shown that the antisymmetrized geminal power (AGP) wavefunction can be an excellent starting point for ansatze describing systems with strong pairing correlations, as those occurring in superconductors. In this work, we show how AGP can be efficiently implemented on a quantum computer with circuit depth, number of CNOTs, and number of measurements being linear in system size. Using AGP as the initial reference, we propose and implement a unitary correlator on AGP and benchmark it on the ground state of the pairing Hamiltonian. The results show highly accurate ground state energies in all correlation regimes of this model Hamiltonian.

preprint2020arXiv

Exploring Hilbert space on a budget: Novel benchmark set and performance metric for testing electronic structure methods in the regime of strong correlation

This work explores the ability of classical electronic structure methods to efficiently represent (compress) the information content of full configuration interaction (FCI) wave functions. We introduce a benchmark set of four hydrogen model systems of different dimensionality and distinctive electronic structures: a 1D chain, a 1D ring, a 2D triangular lattice, and a 3D close-packed pyramid. To assess the ability of a computational method to produce accurate and compact wave functions, we introduce the accuracy volume, a metric that measures the number of variational parameters necessary to achieve a target energy error. Using this metric and the hydrogen models, we examine the performance of three classical deterministic methods: i) selected configuration interaction (sCI) realized both via an a posteriori and variational selection of the most important determinants, ii) rank-reduced FCI, obtained by an a posteriori singular value decomposition of the FCI tensor (SVD-FCI), and iii) the matrix product state representation obtained via the density matrix renormalization group (DMRG). We find that DMRG generally gives the most efficient wave function representation for all systems, particularly in the 1D chain with a localized basis. For the 2D and 3D systems, all methods perform best with a delocalized basis, and the efficiency of sCI is closer to that of DMRG, with the former having and accuracy volume approximately twice as large in the strong correlation regime. Compared to sCI, the SVD-FCI scheme is generally found to require a slightly larger number of parameters to achieve the same energy accuracy.

preprint2019arXiv

Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory

A formal analysis is conducted on the exactness of various forms of unitary coupled cluster (UCC) theory based on particle-hole excitation and de-excitation operators. Both the conventional single exponential UCC parameterization and a disentangled (factorized) version are considered. We formulate a differential cluster analysis to determine the UCC amplitudes corresponding to a general quantum state. The exactness of conventional UCC (ability to represent any state) is explored numerically and it is formally shown to be determined by the structure of the critical points of the UCC exponential mapping. A family of disentangled UCC wave functions are shown to exactly parameterize any state, thus showing how to construct Trotter-error-free parameterizations of UCC for applications in quantum computing. From these results, we derive an exact disentangled UCC parameterization that employs an infinite sequence of particle-hole or general one- and two-body substitution operators.

preprint2019arXiv

Time Dependent Adaptive Configuration Interaction Applied to Attosecond Charge Migration

In this work, we present a time-dependent (TD) selected configuration interaction method based on our recently-introduced adaptive configuration interaction (ACI). We show that ACI, in either its ground or excited state formalisms, is capable of building a compact basis for use in real-time propagation of wave functions for computing electron dynamics. TD-ACI uses an iteratively selected basis of determinants in real-time propagation capable of capturing strong correlation effects in both ground and excited states, all with an accuracy---and associated cost---tunable by the user. We apply TD-ACI to study attosecond-scale migration of charge following ionization in small molecules. We first compute attosecond charge dynamics in a benzene model to benchmark and understand the utility of TD-ACI with respect to an exact solution. Finally, we use TD-ACI to reproduce experimentally determined ultrafast charge migration dynamics in iodoacetylene. TD-ACI is shown to be a valuable benchmark theory for electron dynamics, and it represents an important step towards accurate and affordable time-dependent multireference methods.

preprint2016arXiv

A deterministic projector configuration interaction approach for the ground state of quantum many-body systems

In this work we propose a novel approach to solve the Schrödinger equation which combines projection onto the ground state with a path-filtering truncation scheme. The resulting projector configuration interaction (PCI) approach realizes a deterministic version of the full configuration interaction quantum Monte Carlo (FCIQMC) method [Booth, G. H.; Thom, A. J. W.; Alavi, A. J. Chem. Phys. 2009, 131, 054106]. To improve upon the linearized imaginary-time propagator, we develop an optimal projector scheme based on an exponential Chebyshev expansion in the limit of an infinite imaginary time step. After writing the exact projector as a path integral in determinant space, we introduce a path filtering procedure that truncates the size of the determinantal basis and approximates the Hamiltonian. The path filtering procedure is controlled by one real threshold that determines the accuracy of the PCI energy and is not biased towards any determinant. Therefore, the PCI approach can equally well describe static and dynamic electron correlation. This point is illustrated in benchmark computation on N$_2$ at both equilibrium and stretched geometries. In both cases, the PCI achieves chemical accuracy with wave functions that contain less than 0.5% of the full CI space. We also report computations on the ground state of C$_2$ with up to quaduple-$ζ$ basis sets and wave functions as large as 200 million determinants, which allow a direct comparison of the PCI, FCIQMC, and density matrix renormalization group (DMRG) methods. The size of the PCI wave function grows modestly with the number of unoccupied orbitals and its accuracy may be tuned to match that of FCIQMC and DMRG.

preprint2016arXiv

An adaptive configuration interaction approach for strongly correlated electrons with tunable accuracy

We introduce a new procedure for iterative selection of determinant spaces capable of describing highly correlated systems. This adaptive configuration interaction (ACI) determines an optimal basis by an iterative procedure in which the determinant space is expanded and coarse grained until self consistency. Two importance criteria control the selection process and tune the ACI to a user-defined level of accuracy. The ACI is shown to yield potential energy curves of N$_2$ with nearly constant errors, and it predicts singlet-triplet splittings of acenes up to decacene that are in good agreement with the density matrix renormalization group.

preprint2016arXiv

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 ($Δ_\text{ST}$) of the naphthyne isomers strongly depend on the equilibrium structures. For a consistent set of geometries, the $Δ_\text{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.

preprint2016arXiv

Towards numerically robust multireference theories: The driven similarity renormalization group truncated to one- and two-body operators

The first nonperturbative version of the multireference driven similarity renormalization group (MR-DSRG) theory [C. Li and F. A. Evangelista, J. Chem. Theory Comput. $\mathbf{11}$, 2097 (2015)] is introduced. The renormalization group structure of the MR-DSRG equations ensures numerical robustness and avoidance of the intruder state problem, while the connected nature of the amplitude and energy equations guarantees size consistency and extensivity. We approximate the MR-DSRG equations by keeping only one- and two-body operators and using a linearized recursive commutator approximation of the Baker--Campbell--Hausdorff expansion [T. Yanai and G. K.-L. Chan, J. Chem. Phys. $\mathbf{124}$, 194106 (2006)]. The resulting MR-LDSRG(2) equations contain only 39 terms and scales as ${\cal O}(N^2 N_{\rm P}^2 N_{\rm H}^2)$ where $N_{\rm H}$, $N_{\rm P}$, and $N$ correspond to the number of hole, particle, and total orbitals, respectively. Benchmark MR-LDSRG(2) computations on the hydrogen fluoride and molecular nitrogen binding curves and the singlet-triplet splitting of $p$-benzyne yield results comparable in accuracy to those from multireference configuration interaction, Mukherjee multireference coupled cluster theory, and internally-contracted multireference coupled cluster theory.

preprint2014arXiv

A driven similarity renormalization group approach to quantum many-body problems

Applications of the similarity renormalization group (SRG) approach [F. Wegner, Ann. Phys. 506, 77 (1994), S. D. Głazek and K. G. Wilson, Phys. Rev. D 49, 4214 (1994)] to the formulation of useful many-body theories of electron correlation are considered. In addition to presenting a production-level implementation of the SRG based on a single-reference formalism, a novel integral version of the SRG is reported, in which the flow of the Hamiltonian is driven by a source operator. It is shown that this driven SRG (DSRG) produces a Hamiltonian flow that is analogous to that of the SRG. Compared to the SRG, which requires propagating a set of ordinary differential equations, the DSRG is computationally advantageous since it consists of a set of polynomial equations. The equilibrium distances, harmonic vibrational frequencies, and vibrational anharmonicities of a series of diatomic molecules computed with the SRG and DSRG approximated with one- and two-body normal ordered operators are in good agreement with benchmark values from coupled cluster with singles, doubles, and perturbative triples [CCSD(T)]. Particularly surprising results are found when the SRG and DSRG methods are applied to C$_2$ and F$_2$. In the former case both methods fail to converge, while in the latter case an unbound potential energy curve is obtained. A modified commutator approximation is shown to correct these problems in the case of the DSRG method.

preprint2014arXiv

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff $Λ$. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than $Λ$. The resulting $Λ$-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization ($Λ$+SD-CI), which is based on a small $Λ$-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build $Λ$-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second consists of a breadth-first search coupled with pruning. The $Λ$-CI and $Λ$+SD-CI approaches are used to compute the dissociation curve of N$_2$ and the potential energy curves for the first three singlet states of C$_2$. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the $Λ$-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu$_2$O$_2^{2+}$ core, illustrates an alternative use of the $Λ$-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.