Hybrid Monte Carlo with Fourier Acceleration simulation package for the N=2 Principal Chiral model.
Project description
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
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.