Finding the key to accurate excited states in quantum Monte Carlo

Andrew Williamson, Randolph Q. Hood, Richard J. Needs +1 more

The band structure of silicon is calculated at the $\ensuremath{\Gamma}$, $X$, and $L$ wave vectors using diffusion quantum Monte Carlo (DMC) methods. Excited states are formed by promoting an electron from the valence band into the conduction band. We obtain good agreement with experiment for states around the gap region, and demonstrate that the method works equally well for direct and indirect excitations, and that one can calculate many excited states at each wave vector. This work establishes the fixed-node DMC approach as an accurate method for calculating the energies of low-lying excitations in solids.

Monika Dash, Claudia Filippi

Kristoffer Simula, Ilja Makkonen

Certain point defects in solids can efficiently be used as qubits for applications in quantum technology. They have spin states that are initializable, readable, robust, and can be manipulated optically. New theoretical methods are needed to find the best host materials and defect configurations. Most methods proposed so far rely either on cluster models or restrict the many-body treatment of the defects to a subspace of single-particle orbitals. We explore best practices and theory for the use of quantum Monte Carlo to predict the excitation spectra for spin defects, by using the negatively charged nitrogen vacancy (${\mathrm{NV}}^{\ensuremath{-}}$) center in diamond as a test system. Quantum Monte Carlo can be used to explicitly simulate electronic correlations with larger systems and sets of orbitals than previous methods due to favorable scaling with respect to system size and computing power. We consider different trial wave functions for variational and diffusion Monte Carlo methods, explore the nodal surface errors of the ground and excited state wave functions, and study whether the variational principle holds for the excited states. We compute the vertical excitation energies in different simulation cell sizes and extrapolate to infinite system size, and include backflow corrections to the extrapolated energies. The final results for vertical excitation energies are found to overestimate the experimental estimates, but the triplet-to-triplet and singlet-to-singlet transitions are accurate against experiment. Finally, we list further developments for QMC needed to address the problem of accurately predicting structural and spin properties of the solid-state defects.

D. Blume, Marius Lewerenz, K. Birgitta Whaley

Stuart Shepard, Anthony Scemama, Saverio Moroni +1 more

Quantum Monte Carlo (QMC) methods have proven to be highly accurate for computing excited states, but the choice of optimization strategies for multiple states remains an active topic of investigation. In this work, we revisit the calculation of double excitation energies in nitroxyl, glyoxal, tetrazine, and cyclopentadienone, exploring different objective functionals and their impact on the accuracy and robustness of QMC. A previous study for these systems employed a penalty functional to enforce orthogonality among the states, but the chosen prefactors did not strictly ensure convergence to the target states. Here, we confirm the reliability of previous results by comparing excitation energies obtained with different functionals and analyzing their consistency. In addition, we investigate the performance of different functionals when starting from a pre-collapsed excited state, providing insight into their ability to recover the target wave functions.

Richard J. Needs, Andrew R. Porter, Mike D. Towler