Monte Carlo (MC) simulations are essential computational approaches with widespread use throughout all areas of science. We present a method for accelerating lattice MC simulations using fully-connected and convolutional artificial neural networks that are trained to perform local and global moves in configuration space, respectively. Both networks take local spacetime MC configurations as input features and can, therefore, be trained using samples generated by conventional MC runs on smaller lattices before being utilized for simulations on larger systems. This new approach is benchmarked for the case of determinant quantum Monte Carlo (DQMC) studies of the twodimensional Holstein model. We find that both artificial neural networks are capable of learning an unspecified effective model that accurately reproduces the MC configuration weights of the original Hamiltonian and achieve an order of magnitude speedup over the conventional DQMC algorithm. Our approach is broadly applicable to many classical and quantum lattice MC algorithms.
We study the temperature-filling phase diagram of the single-band Holstein model in two dimensions using the self-consistent Migdal approximation, where both the electron and phonon selfenergies are treated on an equal footing. By employing an efficient numerical algorithm utilizing fast Fourier transforms to evaluate momentum and Matsubara frequency summations, we determine the charge-density-wave (CDW) and superconducting transition temperatures in the thermodynamic limit using lattice sizes that are sufficient to eliminate significant finite size effects present at lower temperatures. We obtain the temperature-filling phase diagrams for a range of coupling strengths and phonon frequencies for the model defined on a square lattice with and without next-nearest neighbor hopping. We find the appearance of a superconducting dome with a critical temperature that decreases before reaching the qmax = ( , ) CDW phase boundary. For very low phonon frequencies, we also find an incommensurate CDW phase with the ordering vector qmax ≈ ( , ) appearing between the commensurate CDW and superconducting phases. Our numerical implementation can be easily extended to treat momentum-dependent electron-phonon coupling, as well as dispersive phonon branches, and has been made available to the public. arXiv:1811.03676v1 [cond-mat.supr-con]
The electron-phonon (e-ph) interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge-density-waves (CDW), and polarons. We study the two-dimensional Holstein model for e-ph coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless e-ph coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate e-ph coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.
Determining the range of validity of Migdal's approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal's approximation holds. Our results demonstrate that questions regarding the validity of Migdal's approximation go hand in hand with questions of the validity of a linear e-ph interaction.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.