Durotun M. S. Z. Raihan P. H.

Sepuluh Nopember Institute of Technology

Bintoro A. Subagyo

Sepuluh Nopember Institute of Technology

Agus Purwanto

Sepuluh Nopember Institute of Technology

DOI: https://doi.org/10.19184/jid.v27i2.60013

ABSTRACT 

We extend the Monte Carlo PNJL framework by introducing nearest-neighbor interactions between fluctuating sub-volumes, resulting in a spatially structured lattice model. A two-dimensional implementation is developed and sampled using the Metropolis algorithm, together with several computational optimizations to improve simulation efficiency. The resulting framework reproduces the expected qualitative behavior of PNJL thermodynamics, including the deconfinement crossover in the Polyakov loop and the emergence of spatial correlations in the chiral field induced by the lattice interaction. In addition to validating the implementation, the Metropolis acceptance rate is found to exhibit a reproducible two-peak structure that brackets the crossover region, suggesting that algorithmic observables may contain information about the underlying thermodynamics. This framework provides a computational foundation for studying spatially correlated fluctuations beyond the mean-field approximation in effective QCD models.

Keywords: PNJL model, Monte Carlo simulation, field fluctuations, lattice model..