Abstract We extend a semiclassical numerical method, bosonic auxiliary-field Monte Carlo, to quantum spin systems. This method breaks the lattice into clusters, solves each cluster precisely and couples them with classical auxiliary fields through classical Monte Carlo simulation. We test the method with antiferromagnetic spin models in one-dimensional chains, square lattices and triangular lattices, and obtain reasonable results at finite temperatures. This algorithm builds a bridge between classical Monte Carlo method and quantum methods. The algorithm can be improved with either progress in classical Monte Carlo sampling or the development of quantum solvers, and can also be further applied to systems with different lattices or interactions.Abstract We extend a semiclassical numerical method, bosonic auxiliary-field Monte Carlo, to quantum spin systems. This method breaks the lattice into clusters, solves each cluster precisely and couples them with classical auxiliary fields through classical Monte Carlo simulation. We test the method with antiferromagnetic spin models in one-dimensional chains, square lattices and triangular lattices, and obtain reasonable results at finite temperatures. This algorithm builds a bridge between classical Monte Carlo method and quantum methods. The algorithm can be improved with either progress in classical Monte Carlo sampling or the development of quantum solvers, and can also be further applied to systems with different lattices or interactions.

Chinese Physics BJan 1, 20250 citations

Accelerating and Converging Stochastic Quantum Chemistry

Verena A. Neufeld

In the last decade, stochastic versions of quantum chemistry methods such as coupled cluster Monte Carlo (CCMC) or full configuration interaction quantum Monte Carlo (FCIQMC) have made highly accurate energy calculations possible that are not accessible to the corresponding deterministic methods (full configuration interaction and coupled cluster) at the same accuracy. CCMC and FCIQMC parallelize well and exploit the sparsity in the wavefunction which decreases memory costs and makes calculations in larger systems tractable. With CCMC it is straightforward to set up high order coupled cluster calculations, such as CCSDTQ5, which includes quintuple excitations explicitly. In this thesis, the convergence of the energy accuracy with the coupled cluster levels up to CCSDTQ5 was tested on the uniform electron gas, a model solid system, for various degrees of electron correlation. This gave information on what coupled cluster level is needed to reach sufficient accuracy when modelling a solid system. Before large solid systems can be modelled, the CCMC and FCIQMC algorithms need to be optimised. The efficiency in one of the crucial steps in these algorithms, the $spawn$ step, was improved, keeping computational and memory costs as low as possible. Furthermore, the convergence of CCMC and FCIQMC was accelerated by employing a quasi-Newton propagation. Using the model system information of what coupled cluster level is needed and having made great progress towards accelerating these methods, the computation of highly accurate energies in solid or large molecular systems should be more feasible in the future.In the last decade, stochastic versions of quantum chemistry methods such as coupled cluster Monte Carlo (CCMC) or full configuration interaction quantum Monte Carlo (FCIQMC) have made highly accurate energy calculations possible that are not accessible to the corresponding deterministic methods (full configuration interaction and coupled cluster) at the same accuracy. CCMC and FCIQMC parallelize well and exploit the sparsity in the wavefunction which decreases memory costs and makes calculations in larger systems tractable. With CCMC it is straightforward to set up high order coupled cluster calculations, such as CCSDTQ5, which includes quintuple excitations explicitly. In this thesis, the convergence of the energy accuracy with the coupled cluster levels up to CCSDTQ5 was tested on the uniform electron gas, a model solid system, for various degrees of electron correlation. This gave information on what coupled cluster level is needed to reach sufficient accuracy when modelling a solid system. Before large solid systems can be modelled, the CCMC and FCIQMC algorithms need to be optimised. The efficiency in one of the crucial steps in these algorithms, the $spawn$ step, was improved, keeping computational and memory costs as low as possible. Furthermore, the convergence of CCMC and FCIQMC was accelerated by employing a quasi-Newton propagation. Using the model system information of what coupled cluster level is needed and having made great progress towards accelerating these methods, the computation of highly accurate energies in solid or large molecular systems should be more feasible in the future.

Dec 20, 20190 citations

METHODS IN QUANTUM MONTE CARLO

Malvin H. Kalos, G.A. Mikhailov

De Gruyter eBooksDec 31, 19870 citations

The Continuous Time Quantum Monte Carlo method as a cluster solver in the Dynamical Cluster Approximation

Karlis Mikelsons, Alexandru Macridin, Mark Jarrell +4 more

Bulletin of the American Physical SocietyMar 13, 20080 citations

spawn

step, was improved, keeping computational and memory costs as low as possible. Furthermore, the convergence of CCMC and FCIQMC was accelerated by employing a quasi-Newton propagation. Using the model system information of what coupled cluster level is needed and having made great progress towards accelerating these methods, the computation of highly accurate energies in solid or large molecular systems should be more feasible in the future.