Contributor info: Dr Steven Johnston

James Neuhaus, Nathan S. Nichols, Debshikha Banerjee, Benjamin Cohen-Stead, Thomas A. Maier, Adrian Del Maestro, Steven Johnston

SciPost Phys. Codebases 39 (2024) · published 12 November 2024 | Toggle abstract · pdf

We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.

James Neuhaus, Nathan S. Nichols, Debshikha Banerjee, Benjamin Cohen-Stead, Thomas A. Maier, Adrian Del Maestro, Steven Johnston

SciPost Phys. Codebases 39-r1.1 (2024) · published 12 November 2024 | Toggle abstract · src

We introduce the SmoQyDEAC.jl package, a Julia implementation of the Differential Evolution Analytic Continuation (DEAC) algorithm [N. S. Nichols et al., Phys. Rev. E 106, 025312 (2022)] for analytically continuing noisy imaginary time correlation functions to the real frequency axis. Our implementation supports fermionic and bosonic correlation functions on either the imaginary time or Matsubara frequency axes, and treatment of the covariance error in the input data. This paper presents an overview of the DEAC algorithm and the features implemented in the SmoQyDEAC.jl package. It also provides detailed benchmarks of the package's output against the popular maximum entropy and stochastic analytic continuation methods. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDEAC.jl or installed using the Julia package manager. The online documentation, including examples, can be accessed at https://smoqysuite.github.io/SmoQyDEAC.jl/stable/.

Benjamin Cohen-Stead, Sohan Malkaruge Costa, James Neuhaus, Andy Tanjaroon Ly, Yutan Zhang, Richard Scalettar, Kipton Barros, Steven Johnston

SciPost Phys. Codebases 29 (2024) · published 28 May 2024 | Toggle abstract · pdf

We introduce the SmoQyDQMC.jl package, a Julia implementation of the determinant quantum Monte Carlo algorithm. SmoQyDQMC.jl supports generalized tight-binding Hamiltonians with on-site Hubbard and generalized electron-phonon ($e$-ph) interactions, including non-linear $e$-ph coupling and anharmonic lattice potentials. Our implementation uses hybrid Monte Carlo methods with exact forces for sampling the phonon fields, enabling efficient simulation of low-energy phonon branches, including acoustic phonons. The SmoQyDQMC.jl package also uses a flexible scripting interface, allowing users to adapt it to different workflows and interface with other software packages in the Julia ecosystem. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDQMC.jl or installed using the Julia package manager. The online documentation, including examples, can be obtained from our document page at https://smoqysuite.github.io/SmoQyDQMC.jl/stable/.

Benjamin Cohen-Stead, Sohan Malkaruge Costa, James Neuhaus, Andy Tanjaroon Ly, Yutan Zhang, Richard Scalettar, Kipton Barros, Steven Johnston

SciPost Phys. Codebases 29-r0.3 (2024) · published 28 May 2024 | Toggle abstract · src

We introduce the SmoQyDQMC.jl package, a Julia implementation of the determinant quantum Monte Carlo algorithm. SmoQyDQMC.jl supports generalized tight-binding Hamiltonians with on-site Hubbard and generalized electron-phonon ($e$-ph) interactions, including non-linear $e$-ph coupling and anharmonic lattice potentials. Our implementation uses hybrid Monte Carlo methods with exact forces for sampling the phonon fields, enabling efficient simulation of low-energy phonon branches, including acoustic phonons. The SmoQyDQMC.jl package also uses a flexible scripting interface, allowing users to adapt it to different workflows and interface with other software packages in the Julia ecosystem. The code for this package can be downloaded from our GitHub repository at https://github.com/SmoQySuite/SmoQyDQMC.jl or installed using the Julia package manager. The online documentation, including examples, can be obtained from our document page at https://smoqysuite.github.io/SmoQyDQMC.jl/stable/.

Krzysztof Bieniasz, Steven Johnston, Mona Berciu

SciPost Phys. 11, 062 (2021) · published 20 September 2021 | Toggle abstract · pdf

Resonant inelastic X-ray scattering (RIXS) is used increasingly for characterizing low-energy collective excitations in materials. RIXS is a powerful probe, which often requires sophisticated theoretical descriptions to interpret the data. In particular, the need for accurate theories describing the influence of electron-phonon ($e$-p) coupling on RIXS spectra is becoming timely, as instrument resolution improves and this energy regime is rapidly becoming accessible. To date, only rather exploratory theoretical work has been carried out for such problems. We begin to bridge this gap by proposing a versatile variational approximation for calculating RIXS spectra in weakly doped materials, for a variety of models with diverse $e$-p couplings. Here, we illustrate some of its potential by studying the role of electron mobility, which is completely neglected in the widely used local approximation based on Lang-Firsov theory. Assuming that the electron-phonon coupling is of the simplest, Holstein type, we discuss the regimes where the local approximation fails, and demonstrate that its improper use may grossly \textit{underestimate} the $e$-p coupling strength.