Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1508.06683

An efficient variational principle for the direct optimization of excited states

Luning Zhao, Eric Neuscamman

Published Aug 26, 2015• 2 citations
arXivHTMLPDF

Abstract

We present a variational function that targets excited states directly based on their position in the energy spectrum, along with a Monte Carlo method for its evaluation and minimization whose cost scales polynomially for a wide class of approximate wave functions. Being compatible with both real and Fock space and open and periodic boundary conditions, the method has the potential to impact many areas of chemistry, physics, and materials science. Initial tests on doubly excited states show that using this method, the Hilbert space Jastrow antisymmetric geminal power ansatz can deliver order-of-magnitude improvements in accuracy relative to equation of motion coupled cluster theory, while a very modest real space multi-Slater Jastrow expansion can achieve accuracies within 0.1 eV of the best theoretical benchmarks for the carbon dimer.

DOI
10.48550/arxiv.1508.06683

Related papers

An Efficient Variational Principle for the Direct Optimization of Excited States

Luning Zhao, Eric Neuscamman

We present a variational principle that enables systematically improvable predictions for individual excited states through an efficient Monte Carlo evaluation. We demonstrate its compatibility with different ansatzes and with both real space and Fock space sampling and discuss its potential for use in the solid state. In numerical demonstrations for challenging molecular excitations, the method rivals or surpasses the accuracy of very high level methods using drastically more compact wave function approximations.

Journal of Chemical Theory and Computation•Jul 5, 2016•66 citations

Communication: Variation after response in quantum Monte Carlo

Eric Neuscamman

We present a new method for modeling electronically excited states that overcomes a key failing of linear response theory by allowing the underlying ground state ansatz to relax in the presence of an excitation. The method is variational, has a cost similar to ground state variational Monte Carlo, and admits both open and periodic boundary conditions. We present preliminary numerical results showing that, when paired with the Jastrow antisymmetric geminal power ansatz, the variation-after-response formalism delivers accuracies for valence and charge transfer single excitations on par with equation of motion coupled cluster, while surpassing coupled cluster's accuracy for excitations with significant doubly excited character.

arXiv:1605.09438•The Journal of Chemical Physics•Aug 25, 2016•19 citations

A Blocked Linear Method for Optimizing Large Parameter Sets in Variational Monte Carlo

Luning Zhao, Eric Neuscamman

We present a modification to variational Monte Carlo's linear method optimization scheme that addresses a critical memory bottleneck while maintaining compatibility with both the traditional ground state variational principle and our recently introduced variational principle for excited states. For wave function ansatzes with tens of thousands of variables, our modification reduces the required memory per parallel process from tens of gigabytes to hundreds of megabytes, making the methodology a much better fit for modern supercomputer architectures in which data communication and per-process memory consumption are primary concerns. We verify the efficacy of the new optimization scheme in small molecule tests involving both the Hilbert space Jastrow antisymmetric geminal power ansatz and real space multi-Slater Jastrow expansions. Satisfied with its performance, we have added the optimizer to the QMCPACK software package, with which we test a systematically convergent, nonperturbative approach to excitation energies on the example of a Mott-insulating hydrogen ring.

arXiv:1702.01481•Journal of Chemical Theory and Computation•May 17, 2017•26 citations

Reduced scaling Hilbert space variational Monte Carlo

Haochuan Wei, Eric Neuscamman

We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through various tensor decomposition methods, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons when using a minimal basis. Further study will be needed to assess where the crossover occurs in more compact molecular geometries and larger basis sets and to explore how in that context the crossover can be accelerated.

arXiv:1806.08778•The Journal of Chemical Physics•Nov 13, 2018•12 citations

Improved Optimization for the Cluster Jastrow Antisymmetric Geminal Power and Tests on Triple-Bond Dissociations

Eric Neuscamman

We present a novel specialization of the variational Monte Carlo linear method for the optimization of the recently introduced cluster Jastrow antisymmetric geminal power ansatz, achieving a lower-order polynomial cost scaling than would be possible with a naive application of the linear method and greatly improving optimization performance relative to that of the previously employed quasi-Newton approach. We test the methodology on highly multireference triple-bond stretches, achieving accuracy superior to those of the traditional coupled cluster theory and multireference perturbation theory in both the typical example of N2 and the transition-metal-oxide example of [ScO](+).

arXiv:1603.06605•Journal of Chemical Theory and Computation•Jun 9, 2016•34 citations

Optimization of quantum Monte Carlo wave functions by energy minimization

Julien Toulouse, Cyrus J. Umrigar

We study three wave function optimization methods based on energy minimization in a variational Monte Carlo framework: the Newton, linear, and perturbative methods. In the Newton method, the parameter variations are calculated from the energy gradient and Hessian, using a reduced variance statistical estimator for the latter. In the linear method, the parameter variations are found by diagonalizing a nonsymmetric estimator of the Hamiltonian matrix in the space spanned by the wave function and its derivatives with respect to the parameters, making use of a strong zero-variance principle. In the less computationally expensive perturbative method, the parameter variations are calculated by approximately solving the generalized eigenvalue equation of the linear method by a nonorthogonal perturbation theory. These general methods are illustrated here by the optimization of wave functions consisting of a Jastrow factor multiplied by an expansion in configuration state functions (CSFs) for the C2 molecule, including both valence and core electrons in the calculation. The Newton and linear methods are very efficient for the optimization of the Jastrow, CSF, and orbital parameters. The perturbative method is a good alternative for the optimization of just the CSF and orbital parameters. Although the optimization is performed at the variational Monte Carlo level, we observe for the C2 molecule studied here, and for other systems we have studied, that as more parameters in the trial wave functions are optimized, the diffusion Monte Carlo total energy improves monotonically, implying that the nodal hypersurface also improves monotonically.

arXiv:physics/0701039•The Journal of Chemical Physics•Feb 22, 2007•290 citations
We present a variational function that targets excited states directly based on their position in the energy spectrum, along with a Monte Carlo method for its evaluation and minimization whose cost scales polynomially for a wide class of approximate wave functions. Being compatible with both real and Fock space and open and periodic boundary conditions, the method has the potential to impact many areas of chemistry, physics, and materials science. Initial tests on doubly excited states show that using this method, the Hilbert space Jastrow antisymmetric geminal power ansatz can deliver order-of-magnitude improvements in accuracy relative to equation of motion coupled cluster theory, while a very modest real space multi-Slater Jastrow expansion can achieve accuracies within 0.1 eV of the best theoretical benchmarks for the carbon dimer.
We present a variational principle that enables systematically improvable predictions for individual excited states through an efficient Monte Carlo evaluation. We demonstrate its compatibility with different ansatzes and with both real space and Fock space sampling and discuss its potential for use in the solid state. In numerical demonstrations for challenging molecular excitations, the method rivals or surpasses the accuracy of very high level methods using drastically more compact wave function approximations.
We present a new method for modeling electronically excited states that overcomes a key failing of linear response theory by allowing the underlying ground state ansatz to relax in the presence of an excitation. The method is variational, has a cost similar to ground state variational Monte Carlo, and admits both open and periodic boundary conditions. We present preliminary numerical results showing that, when paired with the Jastrow antisymmetric geminal power ansatz, the variation-after-response formalism delivers accuracies for valence and charge transfer single excitations on par with equation of motion coupled cluster, while surpassing coupled cluster's accuracy for excitations with significant doubly excited character.
We present a modification to variational Monte Carlo's linear method optimization scheme that addresses a critical memory bottleneck while maintaining compatibility with both the traditional ground state variational principle and our recently introduced variational principle for excited states. For wave function ansatzes with tens of thousands of variables, our modification reduces the required memory per parallel process from tens of gigabytes to hundreds of megabytes, making the methodology a much better fit for modern supercomputer architectures in which data communication and per-process memory consumption are primary concerns. We verify the efficacy of the new optimization scheme in small molecule tests involving both the Hilbert space Jastrow antisymmetric geminal power ansatz and real space multi-Slater Jastrow expansions. Satisfied with its performance, we have added the optimizer to the QMCPACK software package, with which we test a systematically convergent, nonperturbative approach to excitation energies on the example of a Mott-insulating hydrogen ring.
We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through various tensor decomposition methods, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons when using a minimal basis. Further study will be needed to assess where the crossover occurs in more compact molecular geometries and larger basis sets and to explore how in that context the crossover can be accelerated.
We present a novel specialization of the variational Monte Carlo linear method for the optimization of the recently introduced cluster Jastrow antisymmetric geminal power ansatz, achieving a lower-order polynomial cost scaling than would be possible with a naive application of the linear method and greatly improving optimization performance relative to that of the previously employed quasi-Newton approach. We test the methodology on highly multireference triple-bond stretches, achieving accuracy superior to those of the traditional coupled cluster theory and multireference perturbation theory in both the typical example of N2 and the transition-metal-oxide example of [ScO](+).
We study three wave function optimization methods based on energy minimization in a variational Monte Carlo framework: the Newton, linear, and perturbative methods. In the Newton method, the parameter variations are calculated from the energy gradient and Hessian, using a reduced variance statistical estimator for the latter. In the linear method, the parameter variations are found by diagonalizing a nonsymmetric estimator of the Hamiltonian matrix in the space spanned by the wave function and its derivatives with respect to the parameters, making use of a strong zero-variance principle. In the less computationally expensive perturbative method, the parameter variations are calculated by approximately solving the generalized eigenvalue equation of the linear method by a nonorthogonal perturbation theory. These general methods are illustrated here by the optimization of wave functions consisting of a Jastrow factor multiplied by an expansion in configuration state functions (CSFs) for the C2 molecule, including both valence and core electrons in the calculation. The Newton and linear methods are very efficient for the optimization of the Jastrow, CSF, and orbital parameters. The perturbative method is a good alternative for the optimization of just the CSF and orbital parameters. Although the optimization is performed at the variational Monte Carlo level, we observe for the C2 molecule studied here, and for other systems we have studied, that as more parameters in the trial wave functions are optimized, the diffusion Monte Carlo total energy improves monotonically, implying that the nodal hypersurface also improves monotonically.