A new method for Monte Carlo simulation of theories with Grassmann variables

TL;DR

This paper introduces a novel Monte Carlo simulation algorithm for theories with fermions, utilizing a new configuration update method based on solving a specific matrix equation to improve simulation efficiency.

Contribution

The paper presents a new algorithm for simulating theories with fermions, enhancing current methods by solving a matrix equation for configuration updates.

Findings

Algorithm accelerates fermionic theory simulations

Effective in generating new configurations close to previous ones

Potential to improve computational efficiency in lattice QCD

Abstract

A new algorithm for simulation of theories with dynamical fermions is presented. The algorithm is based on obtaining the new configuration U' from the old one U by solving the equation M(U')\eta= \omega M(U)\eta, where M is fermionic operator, \eta is random Gaussian vector, and \omega is random real number close to unity. This algorithm can be used for acceleration of current simulations in theories with fermions.