SU2xSU2 2

Hybrid Monte Carlo with Fourier Acceleration simulation package for the N=2 Principal Chiral model.

This python package offers efficient simulation and data analysis routines for the $SU(2) \times SU(2)$ Principal Chiral model. The key feature offered is the integration of Fourier Acceleration into the Hybrid Monte Carlo algorithm which leads to a significant reduction in the degree of critical slowing down.

The simulation is supported for cubic lattices with even side length in arbitary dimensions.

Installation

To install SU2xSU2 using pip run:

pip install SU2xSU2

Its is recommended to work in a virtual environment. The package comes with a custom style sheet which is used by default.

Documentation

Read the docs here.

Example

A basic example showing how to set up a simulation using Fourier accelerated HMC to measure the wall-to-wall correlation function. Further examples can be found here.

from SU2xSU2.SU2xSU2 import SU2xSU2

# define model and lattice parameters 
model_paras = {'D':2, L':40, 'a':1, 'ell':5, 'eps':1/5, 'beta':0.6}
model = SU2xSU2(**model_paras)
# define simulation parameters and measurements
sim_paras = {'M':500, 'burnin_frac':0.5, 'accel':True, 'measurements':[model.ww_correlation_func], 'chain_paths':['corfunc_chain.npy']}
model.run_HMC(**sim_paras) 

Attribution

Please cite the following paper if you found this code useful in your research:

@article{WackFAHMC2024,
    author = {Roger Horsley and Brian Pendleton and Julian Wack},
    title = {Hybrid Monte Carlo simulation with Fourier acceleration of the N = 2 principal chiral model in two dimensions},
    journal = {Physics Letters B},
    volume = {849},
    pages = {138429},
    year = {2024},
    issn = {0370-2693},
    doi = {https://doi.org/10.1016/j.physletb.2023.138429},
    eprint={2308.14628v2},
    archivePrefix={arXiv},
    primaryClass={hep-lat}
}

Licence

SU2xSU2 is free software made available under the MIT License. For details see the LICENSE file.

To Do

• Runtime warning in correlations l.64
• implement weak coupling expansion for all D