Source author record

Piotr Korcyl

Piotr Korcyl appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2026arXiv

Sampling two-dimensional spin systems with transformers

Autoregressive Neural Networks based on dense or convolutional layers have recently been shown to be a viable strategy for generating classical spin systems. Unlike these methods, sampling with transformers is commonly considered to be computationally inefficient. In this work, we propose a novel approach to transformer-based neural samplers in which we generate not a single spin per step but groups of spins. As an additional improvement, we construct a model of approximated probabilities, further improving the efficiency of the algorithm. Despite our approach being computationally heavier than dense networks or CNN-based approaches, we were able to sample larger systems of up to $180 \times 180$ spins in case of the Ising model. The Effective Sample Size of our sampler is $\sim 20$ times larger than that of the previous state-of-the-art neural sampler when trained for the $128 \times 128$ Ising model at critical temperature. Finally, we also test our algorithm on the 2D Edwards-Anderson model, where we train $64\times 64$ spin systems.

preprint2026arXiv

Variational Autoregressive Networks with probability priors

Monte Carlo methods are essential across diverse scientific fields, yet their efficiency is frequently hampered by critical slowing down-a sharp increase in autocorrelation times near phase transitions. Although deep learning approaches, such as neural-network-based samplers, have been proposed to alleviate this issue, they face another serious problem: the difficulty of training the models. This difficulty partially stems from the overly general nature of original machine-learning architectures, which often ignore underlying physical symmetries and force networks to relearn them from scratch. In this paper, we demonstrate that incorporating physical priors into the model significantly enhances performance. Building upon existing strategies that integrate spin-spin interactions, we propose a framework that utilizes a prior probability distribution as a starting point for training. Our results for the Ising model, as well as for the Edwards-Anderson spin glass model, suggest that moving away from `blank slate' models in favor of physics-informed priors reduces the training burden and facilitates the simulation of larger system sizes in discrete spin models.

preprint2023arXiv

Analysis of autocorrelation times in Neural Markov Chain Monte Carlo simulations

We provide a deepened study of autocorrelations in Neural Markov Chain Monte Carlo (NMCMC) simulations, a version of the traditional Metropolis algorithm which employs neural networks to provide independent proposals. We illustrate our ideas using the two-dimensional Ising model. We discuss several estimates of autocorrelation times in the context of NMCMC, some inspired by analytical results derived for the Metropolized Independent Sampler (MIS). We check their reliability by estimating them on a small system where analytical results can also be obtained. Based on the analytical results for MIS we propose a new loss function and study its impact on the autocorelation times. Although, this function's performance is a bit inferior to the traditional Kullback-Leibler divergence, it offers two training algorithms which in some situations may be beneficial. By studying a small, $4 \times 4$, system we gain access to the dynamics of the training process which we visualize using several observables. Furthermore, we quantitatively investigate the impact of imposing global discrete symmetries of the system in the neural network training process on the autocorrelation times. Eventually, we propose a scheme which incorporates partial heat-bath updates which considerably improves the quality of the training. The impact of the above enhancements is discussed for a $16 \times 16$ spin system. The summary of our findings may serve as a guidance to the implementation of Neural Markov Chain Monte Carlo simulations for more complicated models.

preprint2022arXiv

Gradient estimators for normalising flows

Recently a machine learning approach to Monte-Carlo simulations called Neural Markov Chain Monte-Carlo (NMCMC) is gaining traction. In its most popular form it uses neural networks to construct normalizing flows which are then trained to approximate the desired target distribution. In this contribution we present new gradient estimator for Stochastic Gradient Descent algorithm (and the corresponding \texttt{PyTorch} implementation) and show that it leads to better training results for $ϕ^4$ model. For this model our estimator achieves the same precision in approximately half of the time needed in standard approach and ultimately provides better estimates of the free energy. We attribute this effect to the lower variance of the new estimator. In contrary to the standard learning algorithm our approach does not require estimation of the action gradient with respect to the fields, thus has potential of further speeding up the training for models with more complicated actions.

preprint2021arXiv

Finite volume effects in the McLerran-Venugopalan initial condition for the JIMWLK equation

We revisit the numerical construction of the initial condition for the dipole amplitude from the McLerran-Venugopalan model in the context of the JIMWLK evolution equation. We observe large finite volume effects induced by the Poisson equation formulated on a torus. We show that the situation can be partially cured by introducing an infrared regularization. We propose a procedure that has negligible finite volume corrections. The control of the finite volume and finite lattice spacings effects is crucial when considering the numerical solutions of the JIMWLK evolution equation with the collinear improvement.

preprint2020arXiv

Running coupling constant from position-space current-current correlation functions in three-flavor lattice QCD

In this Letter, we provide a determination of the coupling constant in three-flavor quantum chromodynamics (QCD), $α^{\overline{\mathrm{MS}}}_s(μ)$, for $\overline{\mathrm{MS}}$ renormalization scales $μ\in (1,\,2)$ GeV. The computation uses gauge field configuration ensembles with $\mathcal{O}(a)$-improved Wilson-clover fermions generated by the Coordinated Lattice Simulations (CLS) consortium. Our approach is based on current-current correlation functions and has never been applied before in this context. We convert the results perturbatively to the QCD $Λ$-parameter and obtain $Λ_{\overline{\mathrm{MS}}}^{N_f=3} = 342 \pm 17$ MeV, which agrees with the world average published by the Particle Data Group and has competing precision. The latter was made possible by a unique combination of state-of-the-art CLS ensembles with very fine lattice spacings, further reduction of discretization effects from a dedicated numerical stochastic perturbation theory simulation, combining data from vector and axial-vector channels and matching to high-order perturbation theory.

preprint2016arXiv

Non-perturbative running of renormalization constants from correlators in coordinate space using step scaling

Working in a quenched setup with Wilson twisted mass valence fermions, we explore the possibility to compute non-perturbatively the step scaling function using the coordinate (X-space) renormalization scheme. This scheme has the advantage of being on-shell and gauge invariant. The step scaling method allows us to calculate the running of the renormalization constants of quark bilinear operators. We describe here the details of this calculation. The aim of this exploratory study is to identify the feasibility of the X-space scheme when used in small volume simulations required by the step scaling technique. Eventually, we translate our final results to the continuum MSbar scheme and compare against four-loop analytic formulae finding satisfactory agreement.

preprint2016arXiv

Step scaling in coordinate space: running of the quark mass

We perform a benchmark study of the step scaling procedure for the ratios of renormalization constants extracted from position space correlation functions. We work in the quenched approximation and consider the pseudoscalar, scalar, vector and axial vector bilinears. The pseudoscalar/scalar cases allow us to obtain the non-perturbative running of the quark mass over a wide range of energy scales - from around 17 GeV to below 1.5 GeV - which agrees well with the 4-loop prediction of continuum perturbation theory. We find that step scaling is feasible in X-space and we discuss its advantages and potential problems.

preprint2015arXiv

Simulation of QCD with N_f=2+1 flavors of non-perturbatively improved Wilson fermions

We describe a new set of gauge configurations generated within the CLS effort. These ensembles have N_f=2+1 flavors of non-perturbatively improved Wilson fermions in the sea with the Luescher-Weisz action used for the gluons. Open boundary conditions in time are used to address the problem of topological freezing at small lattice spacings and twisted-mass reweighting for improved stability of the simulations. We give the bare parameters at which the ensembles have been generated and how these parameters have been chosen. Details of the algorithmic setup and its performance are presented as well as measurements of the pion and kaon masses alongside the scale parameter t_0.

preprint2014arXiv

Numerical investigations of Supersymmetric Yang-Mills Quantum Mechanics with 4 supercharges

We report on our non-perturbative investigations of supersymmetric Yang-Mills quantum mechanics with 4 supercharges. We employ two independent numerical methods. First of them is the cut Fock space method whose numerical implementation was recently generalized to include the SU(N) gauge group. It allowed us to calculate for the first time the spectrum of the model with SU(3) symmetry in all fermionic sectors. Independently, we implemented the Rational Hybrid Monte Carlo algorithm and reproduced the accessible part of the low-energy spectrum of the model with SU(2) gauge symmetry. We argue that by simulating at imaginary chemical potential one can get access to observables defined in sectors of Hilbert space with a given quark number.

preprint2014arXiv

On the extraction of spectral quantities with open boundary conditions

We discuss methods to extract decay constants, meson masses and gluonic observables in the presence of open boundary conditions. The ensembles have been generated by the CLS effort and have 2+1 flavors of O(a)-improved Wilson fermions with a small twisted-mass term as proposed by Lüscher and Palombi. We analyse the effect of the associated reweighting factors on the computation of different observables.

preprint2014arXiv

The CLS 2+1 flavor simulations

We report on the status of large volume simulations with 2+1 dynamical fermions which are being performed by the CLS initiative. The algorithmic details include: open boundary conditions, twisted mass reweighting and RHMC, whereas the main feature of the simulation strategy is the approach to the physical point along a trajectory of constant trace of the mass matrix. We comment on the practical side of the above issues using as examples some of the newly generated ensembles, which presently cover lattice spacings between 0.05 fm and 0.11 fm and pion masses between 150 MeV and 415 MeV.

preprint2013arXiv

Fixing the parameters of Lattice HQET including 1/m_B terms

The study of the CKM matrix elements with increasing precision requires a reliable evaluation of hadronic matrix elements of axial and vector currents which can be done with Lattice QCD. The tiniest entry, |V_{ub}|, can be estimated independently from B to tau nu and B to pi l nu decays. The ALPHA collaboration has undertaken the effort to evaluate non-perturbatively the decay constant f_B and the f^+(q) form factor for q^2 close to q_{max}^2 entering these determinations. Since for the b quark m_b >> a^{-1} for available lattice sizes, an effective description of the b quark is necessary. HQET provides an example of the latter. As any effective theory, HQET is predictive only when a set of parameters have been determined through a process called matching. The non-perturbative matching procedure applied by the ALPHA collaboration consists of 19 matching conditions needed to fix all the relevant parameters at order 1/m_B of the HQET action and the axial and vector currents. We present a study of one-loop corrections to two representative matching conditions. Our results enable us to quantify the quality of the observables used in the matching procedure.

preprint2013arXiv

On non-trivial spectra of trivial gauge theories

In this Letter we point out that the analytic solution of the two dimensional U(1) gauge theory, on a finite lattice, reveals in the continuum limit the renowned Manton's spectrum of topological electric fluxes together with their effective hamiltonian and wave functions. We extend this result for the system with strings and external charges providing also a novel interpretation of the theta parameter. Some further generalizations are also outlined.

preprint2012arXiv

Non-perturbative renormalization in coordinate space for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action

We present results of a lattice QCD application of a coordinate space renormalization scheme for the extraction of renormalization constants for flavour non-singlet bilinear quark operators. The method consists in the analysis of the small-distance behaviour of correlation functions in Euclidean space and has several theoretical and practical advantages, in particular: it is gauge invariant, easy to implement and has relatively low computational cost. The values of renormalization constants in the X-space scheme can be converted to the MSbar scheme via 4-loop continuum perturbative formulae. Our results for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action are compared to the ones from the RI-MOM scheme and show full agreement with this method.

preprint2012arXiv

Screening in two-dimensional gauge theories

We analyze the problem of screening in 1+1 dimensional gauge theories. Using QED2 as a warm-up for the non-abelian models we show the mechanism of the string breaking, in particular the vanishing overlap of the Wilson loops to the broken-string ground state that has been conjectured in higher-dimensional analyses. We attempt to extend our analysis to non-integer charges in the quenched and unquenched cases, in pursuit of the numerical check of a renowned result for the string tension between arbitrarily-charged fermions in the massive Schwinger model.

preprint2011arXiv

Preliminary study of two-dimensional SU(N) Yang-Mills theory with adjoint matter by Hybrid Monte Carlo approach

Two-dimensional non-abelian quantum field models provide a useful laboratory for analytic and numerical investigations of quantum theories with gauge symmetry. They can exhibit various features, such as charge confinement, which are known from D=4 theories like QCD. Several analytic predictions concerning the spectra of two-dimensional systems with adjoint matter were postulated and numerical results were obtained using Discrete Light Cone Quantization techniques, however none of them has been checked via Monte Carlo simulations. In this Letter we present two such models which are particularly interesting from the physical point of view and discuss first numerical results.

preprint2009arXiv

Exact solutions to D=2 Supersymmetric Yang-Mills Quantum Mechanics with SU(3) gauge group

In this article we present the cut Fock space approach to the D=d+1=2, Supersymmetric Yang-Mills Quantum Mechanics (SYMQM). We start by briefly introducing the main features of the framework. We concentrate on those properties of the method which make it a convenient set up not only for numerical calculations but also for analytic computations. In the main part of the article a sample of results are discussed, namely, analytic and numerical analysis of the D=2, SYMQM systems with SU(2) and SU(3) gauge symmetry.

preprint2009arXiv

Recursive approach to supersymmetric quantum mechanics for arbitrary fermion occupation number

We present in details a numerical approach for solving supersymmetric quantum mechanical systems with a gauge symmetry valid in all fermionic sectors. The method uses a recursive algorithm to calculate matrix elements of any gauge invariant operator in the Fock basis, in particular of the Hamiltonian operator, and can be used for any gauge group. We describe its application to a supersymmetric anharmonic oscillator model with discrete spectrum.