Researcher profile

Andrei Alexandru

Andrei Alexandru contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

hep-lat nucl-th hep-ph hep-th quant-ph cond-mat.dis-nn nucl-ex astro-ph.HE cond-mat.quant-gas cond-mat.str-el hep-ex Machine Learning Performance

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

30works

0followers

13topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2023arXiv

Nucleon Electric Dipole Moment from the $θ$ Term with Lattice Chiral Fermions

We calculate the nucleon electric dipole moment (EDM) from the $θ$ term with overlap fermions on three domain wall lattices with different sea pion masses at lattice spacing 0.11 fm. Due to the chiral symmetry conserved by the overlap fermions, we have well defined topological charge and chiral limit for the EDM. Thus, the chiral extrapolation can be carried out reliably at nonzero lattice spacings. We use three to four different partially quenched valence pion masses for each sea pion mass and find that the EDM dependence on the valence and sea pion masses behaves oppositely, which can be described by partially quenched chiral perturbation theory. With the help of the cluster decomposition error reduction (CDER) technique, we determine the neutron and proton EDM at the physical pion mass to be $d_{n}=-0.00148\left(14\right)\left(31\right)\barθ$ e$\cdot$fm and $d_{p}=0.0038\left(11\right)\left(8\right)\barθ$ e$\cdot$fm. This work is a clear demonstration of the advantages of using chiral fermions in the nucleon EDM calculation and paves the road to future precise studies of the strong $CP$ violation effects.

preprint2022arXiv

Anderson Metal-to-Critical Transition in QCD

A picture of thermal QCD phase change based on the analogy with metal-to-insulator transition of Anderson type was proposed in the past. In this picture, a low-$T$ thermal state is akin to a metal with deeply infrared (IR) Dirac modes abundant and extended, while a high-$T$ state is akin to an insulator with IR modes depleted and localized below a mobility edge $λ_{\text A} > 0$. Here we argue that, while $λ_{\text A}$ exists in QCD, a high-$T$ state is not an insulator in such an analogy. Rather, it is a critical state arising due to a new singular mobility edge at $λ_{\text IR}=0$. This new mobility edge appears upon the transition into the recently proposed IR phase. As a key part of such a metal-to-critical scenario, we present evidence using pure-glue QCD that deeply infrared Dirac modes in the IR phase extend to arbitrarily long distances. This is consistent with our previous suggestion that the IR phase supports scale invariance in the infrared. We discuss the role of Anderson-like aspects in this thermal regime and emphasize that the combination of gauge field topology and disorder plays a key role in shaping its IR physics. Our conclusions are conveyed by the structure of Dirac spectral non-analyticities.

preprint2022arXiv

Applications of Machine Learning to Lattice Quantum Field Theory

There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment of this approach in the future.

preprint2022arXiv

Pole position of the $a_1(1260)$ resonance in a three-body unitary framework

Masses, widths, and branching ratios of hadronic resonances are quantified by their pole positions and residues with respect to transition amplitudes on the Riemann sheets of the complex energy-plane. In this study we discuss the analytic structure in the physical energy region of three-body scattering amplitudes on such manifolds. As an application, we determine the pole position of the $a_1(1260)$ meson from the ALEPH experiment by allowing for $πρ$ coupled channels in S- and D-wave. We find it to be $\sqrt{s_0}=(1232^{+15+9}_{-0-11}-i266^{+0+15}_{-22-27})~\text{MeV}$.

preprint2021arXiv

The spectrum of qubitized QCD: glueballs in a $S(1080)$ gauge theory

Quantum simulations of QCD require digitization of the infinite-dimensional gluon field. Schemes for doing this with the minimum amount of qubits are desirable. We present a practical digitization for $SU(3)$ gauge theories via its discrete subgroup $S(1080)$. Using a modified action that allows classical simulations down to $a\approx 0.08$ fm, the low-lying glueball spectrum is computed with percent-level precision at multiple lattice spacings and shown to extrapolate to the continuum limit $SU(3)$ results. This suggests that this digitization scheme is sufficient for precision quantum simulations of QCD.

preprint2021arXiv

Universality of a truncated sigma-model

Bosonic quantum field theories, even when regularized using a finite lattice, possess an infinite dimensional Hilbert space and, therefore, cannot be simulated in quantum computers with a finite number of qubits. A truncation of the Hilbert space is then needed and the physical results are obtained after a double limit: one to remove the truncation and another to remove the regulator (the continuum limit). A simpler alternative is to find a model with a finite dimensional Hilbert space belonging to the same universality class as the continuum model (a "qubitization"), so only the space continuum limit is required. A qubitization of the $1+1$ dimensional asymptotically free $O(3)$ nonlinear $σ$-model based on ideas of non-commutative geometry was previously proposed arXiv:1903.06577 and, in this paper, we provide evidence that it reproduces the physics of the $σ$-model both in the infrared and the ultraviolet regimes.

preprint2020arXiv

Complex Paths Around The Sign Problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of Physics, from QCD/nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half-filling, dynamical and non-equilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called "sign-problem". We review a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it. The mathematical ideas underpinning this approach, as well as the algorithms so far developed, are described together with non-trivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

preprint2020arXiv

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors $G^c_{E,M}(Q^2)$ in the momentum transfer range $0\leq Q^2 \leq 1.4$ $\rm GeV^2$. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment $μ^c_M=-0.00127(38)_{\rm stat}(5)_{\rm sys}$, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero $G^c_{E}(Q^2)$ indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the $[c(x)-\bar{c}(x)]$ distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

preprint2020arXiv

Finite-volume energy spectrum of the $K^-K^-K^-$ system

The dynamics of multi-kaon systems are of relevance for several areas of nuclear physics. However, even the simplest systems, two and three kaons, are hard to prepare and study experimentally. Here we show how to extract this information using first-principle lattice QCD results. We (1) extend the relativistic three-body quantization condition to the strangeness sector, predicting for the first time the excited level finite-volume spectrum of three kaon systems at maximal isospin, and (2) present a first lattice QCD calculation of the excited levels of this system in a finite box. We compare our predictions with the lattice results reported here and with previous ground state calculations and find very good agreement.

preprint2020arXiv

Quantum algorithms for disordered physics

We show how a quantum computer may efficiently simulate a disordered Hamiltonian, by incorporating a pseudo-random number generator directly into the time evolution circuit. This technique is applied to quantum simulation of few-body disordered systems in the large volume limit; in particular, Anderson localization. The method requires a number of (error corrected) qubits proportional to the logarithm of the volume of the system, and each time evolution step requires a number of gates polylogarithmic in the volume. We simulate the method to observe the metal-insulator transition on a three-dimensional lattice. Additionally, we demonstrate the algorithm on a one-dimensional lattice, using physical quantum processors.

preprint2020arXiv

Roper State from Overlap Fermions

The Roper state is extracted with valence overlap fermions on a $2+1$-flavor domain-wall fermion lattice (spacing $a = 0.114$ fm and $m_π = 330$ MeV) using both the Sequential Empirical Bayes (SEB) method and the variational method. The results are consistent, provided that a large smearing-size interpolation operator is included in the variational calculation to have better overlap with the lowest radial excitation. Similar calculations carried out for an anisotropic clover lattice with similar parameters find the Roper $\approx 280$ MeV higher than that of the overlap fermion. The fact that the prediction of the Roper state by overlap fermions is consistently lower than those of clover fermions, chirally improved fermions, and twisted-mass fermions over a wide range of pion masses has been dubbed a "Roper puzzle." To understand the origin of this difference, we study the hairpin $Z$-diagram in the isovector scalar meson ($a_0$) correlator in the quenched approximation. Comparing the $a_0$ correlators for clover and overlap fermions, at a pion mass of 290 MeV, we find that the spectral weight of the ghost state with clover fermions is smaller than that of the overlap at $a = 0.12$ fm and $0.09$ fm, whereas the whole $a_0$ correlators of clover and overlap at $a = 0.06$ fm coincide within errors. This suggests that chiral symmetry is restored for clover at $a \le 0.06$ fm and that the Roper should come down at and below this $a$. We conclude that this work supports a resolution of the "Roper puzzle" due to $Z$-graph type chiral dynamics. This entails coupling to higher components in the Fock space (e.g. $Nπ$, $Nππ$ states) to induce the effective flavor-spin interaction between quarks as prescribed in the chiral quark model, resulting in the parity-reversal pattern as observed in the experimental excited states of $N, Δ$ and $Λ$.

preprint2019arXiv

A cross-channel study of pion scattering from lattice QCD

We use a chiral model for pion interactions, in the inverse amplitude formalism, to perform a simultaneous analysis of lattice QCD results for pion-pion scattering in all three isospin channels. The input is the finite-volume two-pion spectrum computed using lattice QCD from six ensembles on lattices elongated in one of the spatial dimensions. A two-flavor dynamical lattice QCD action is used with two quark masses corresponding to a pion mass of 315 MeV and 224 MeV. The spectrum in the elastic region is subjected to a global fit which takes into account full correlations across isospin, pion mass and decay constant. The parameters from the fit are used to perform a chiral extrapolation to the physical point. The cross-channel fit results in a more precise determination of the parameters of the model when compared with single channel fits. We obtain $m_πa_0^{I=0}=0.2132(9)$, and $m_πa_0^{I=2}=0.0433(2)$ as well as $m_σ=443(3)-i221(6)$ MeV and $m_ρ=724(4)-i67(1)$ MeV. Several aspects of scale setting and consistency with previous analyses of lattice QCD results are discussed as well.

preprint2019arXiv

Structure Factors of The Unitary Gas Under Supernova Conditions

We compute with lattice field theory the vector and axial static structure factors of the unitary gas for arbitrary temperature above the superfluid transition and for fugacities 0.1 < z < 1.0. Using the lattice formulation, we calculate beyond the validity of the virial expansion, a commonly used technique in many-body physics. We find qualitative differences in the behavior of the structure factors at high fugacity compared to the predictions of the virial expansion. Due to the large scattering length of neutrons, we expect the unitary gas structure factors to approximate the structure factors of hot neutron gases, and we therefore expect our calculations to be useful in supernova simulations, where neutron gas structure factors are needed to compute in-medium neutrino-neutron scattering rates.

preprint2013arXiv

Dynamical Local Chirality and Chiral Symmetry Breaking

We present some of the reasoning and results substantiating the notion that spontaneous chiral symmetry breaking (SChSB) in QCD is encoded in local chiral properties of Dirac eigenmodes. Such association is possible when viewing chirality as a dynamical effect, measured with respect to the benchmark of statistically independent left-right components. Following this rationale leads to describing local chiral behavior by a taylor-made correlation, namely the recently introduced correlation coefficient of polarization C_A. In this language, correlated modes (C_A>0) show dynamical preference for local chirality while anti-correlated modes (C_A<0) favor anti-chirality. Our conclusion is that SChSB in QCD can be viewed as dominance of low-energy correlation (chirality) over anti-correlation (anti-chirality) of Dirac sea. The spectral range of local chirality, chiral polarization scale Lambda_ch, is a dynamically generated scale in the theory associated with SChSB. One implication of these findings is briefly discussed.

preprint2012arXiv

Chiral Polarization Scale at Finite Temperature

We study the chiral polarization properties of low-lying Dirac eigenmodes at finite temperature using the overlap operator. Results for pure gauge theory on both sides of deconfinement phase transition are presented. We find that the polarization scale decreases as we increase the temperature, but it remains non-zero as we cross in the deconfined phase and vanishes only when $T\approx 1.4 T_c$. This is caused by the presence of near-zero modes which, we find, are chirally polarized.

preprint2012arXiv

Sea Contributions to Hadron Electric Polarizabilities through Reweighting

As part of our ongoing lattice study of the electric polarizabilities of hadrons using the background field approach, we use reweighting to examine the effect of the field on the sea quarks. As with other reweighting studies, the chief difficulty lies in the construction of a stochastic estimate of the ratio of the fermion determinants. In contrast to the case of reweighting in the quark mass, these estimators converge extremely slowly, and are resistant to common variance-reduction techniques such as low-mode subtraction. However, it is possible to construct an alternate estimator, taking advantage of the fact that we are interested in only perturbatively small fields; this estimator is susceptible to a variance-reduction technique based on a hopping parameter expansion.

preprint2011arXiv

$ρ$ meson decay on asymmetrical lattices

We present a lattice QCD calculation of the characteristics of the $ρ$ meson decay. The study is carried out on spatially asymmetric boxes using nHYP-smeared clover fermions in the quenched approximation. The resonance mass and coupling constant are calculate using the P-wave scattering phaseshifts, of the isospin I=1 two-pion system. We use pion masses m_π= 418 MeV and m_π=312 MeV. In both cases, the $ρ$ decay is kinematically feasible. We work on lattice sizes N_z X 24^2 X 48 with lattice spacing a=0.1 fm and N_z=24,30,34,48.

preprint2011arXiv

Absolute X-distribution and self-duality

Various models of QCD vacuum predict that it is dominated by excitations that are predominantly self-dual or anti-self-dual. In this work we look at the tendency for self-duality in the case of pure-glue SU(3) gauge theory using the overlap-based definition of the field-strength tensor. To gauge this property, we use the absolute X-distribution method which is designed to quantify the dynamical tendency for polarization for arbitrary random variables that can be decomposed in a pair of orthogonal subspaces.

preprint2011arXiv

Critical point of $N_f = 3$ QCD from lattice simulations in the canonical ensemble

A canonical ensemble algorithm is employed to study the phase diagram of $N_f = 3$ QCD using lattice simulations. We lock in the desired quark number sector using an exact Fourier transform of the fermion determinant. We scan the phase space below $T_c$ and look for an S-shape structure in the chemical potential, which signals the coexistence phase of a first order phase transition in finite volume. Applying Maxwell construction, we determine the boundaries of the coexistence phase at three temperatures and extrapolate them to locate the critical point. Using an improved gauge action and improved Wilson fermions on lattices with a spatial extent of $1.8 \fm$ and quark masses close to that of the strange, we find the critical point at $T_E = 0.925(5) T_c$ and baryon chemical potential $μ_B^E = 2.60(8) T_c$.

preprint2011arXiv

Efficient implementation of the overlap operator on multi-GPUs

Lattice QCD calculations were one of the first applications to show the potential of GPUs in the area of high performance computing. Our interest is to find ways to effectively use GPUs for lattice calculations using the overlap operator. The large memory footprint of these codes requires the use of multiple GPUs in parallel. In this paper we show the methods we used to implement this operator efficiently. We run our codes both on a GPU cluster and a CPU cluster with similar interconnects. We find that to match performance the CPU cluster requires 20-30 times more CPU cores than GPUs.

preprint2011arXiv

Electric polarizability of hadrons with overlap fermions on multi-GPUs

Electric polarizability is an important parameter for the internal structure of hadrons. Previous studies of polarizabilities have been done at relatively heavy pion masses, leaving the chiral region largely unexplored. In this report, we use overlap fermions which are known to be computationally demanding to properly capture the chiral dynamics. We present an implementation strategy to construct overlap on multi-GPUs. We find that our GPU code has an equivalent of \sim30 CPU cores to 1 GPU. We also present preliminary results for the polarizability of the neutral pion.

preprint2011arXiv

The Analysis of Space-Time Structure in QCD Vacuum II: Dynamics of Polarization and Absolute X-Distribution

We propose a framework for quantitative evaluation of dynamical tendency for polarization in arbitrary random variable that can be decomposed into a pair of orthogonal subspaces. The method uses measures based on comparisons of given dynamics to its counterpart with statistically independent components. The formalism of previously considered X-distributions is used to express the aforementioned comparisons, in effect putting the former approach on solid footing. Our analysis leads to definition of a suitable correlation coefficient with clear statistical meaning. We apply the method to the dynamics induced by pure-glue lattice QCD in local left-right components of overlap Dirac eigenmodes. It is found that, in finite physical volume, there exists a non-zero physical scale in the spectrum of eigenvalues such that eigenmodes at smaller (fixed) eigenvalues exhibit convex X-distribution (positive correlation), while at larger eigenvalues the distribution is concave (negative correlation). This chiral polarization scale thus separates a regime where dynamics enhances chirality relative to statistical independence from a regime where it suppresses it, and gives an objective definition to the notion of "low" and "high" Dirac eigenmode. We propose to investigate whether the polarization scale remains non-zero in the infinite volume limit, in which case it would represent a new kind of low energy scale in QCD.

preprint2010arXiv

Finite density phase transition of QCD with $N_f=4$ and $N_f=2$ using canonical ensemble method

In a progress toward searching for the QCD critical point, we study the finite density phase transition of $N_f = 4$ and 2 lattice QCD at finite temperature with the canonical ensemble approach. We develop a winding number expansion method to accurately project out the particle number from the fermion determinant which greatly extends the applicable range of baryon number sectors to make the study feasible. Our lattice simulation was carried out with the clover fermions and improved gauge action. For a given temperature, we calculate the baryon chemical potential from the canonical approach to look for the mixed phase as a signal for the first order phase transition. In the case of $N_f=4$, we observe an "S-shape" structure in the chemical potential-density plane due to the surface tension of the mixed phase in a finite volume which is a signal for the first order phase transition. We use the Maxwell construction to determine the phase boundaries for three temperatures below $T_c$. The intersecting point of the two extrapolated boundaries turns out to be at the expected first order transition point at $T_c$ with $μ= 0$. This serves as a check for our method of identifying the critical point. We also studied the $N_f =2$ case, but do not see a signal of the mixed phase for temperature as low as 0.83 $T_c$.

preprint2010arXiv

Hadron electric polarizability -- finite volume corrections

We use the background field method to extract the polarizability for the neutral "pion". In our previous study we found that the polarizability for this system is negative which is believed to be a finite volume artifact. To address this issue, we carry out simulations for different lattice sizes and we also look at the influence of the boundary conditions on these results. We find that for pion masses lower than 400 MeV the polarizability remains negative even on larger lattices. An infinite volume extrapolation is attempted, but the results are not conclusive due mainly to a lack of an analytical form for the finite volume corrections for this system.

preprint2010arXiv

Magnetic Moments of Negative-Parity Baryons from Lattice QCD

We report preliminary results for the magnetic moments of negative-parity baryons extracted from mass shifts in the presence of static external magnetic fields. The calculations are done on $24^3\times 48$ quenched lattices using standard Wilson actions, with $β$=6.0 and pion mass down to about 520 MeV, and 1000 configurations. Reasonable signals for the negative-parity states are observed and the sign of their magnetic moments is established. The results are compared to model calculations.

preprint2010arXiv

QCD at non-zero density and canonical partition functions with Wilson fermions

We present a reduction method for Wilson Dirac fermions with non-zero chemical potential which generates a dimensionally reduced fermion matrix. The size of the reduced fermion matrix is independent of the temporal lattice extent and the dependence on the chemical potential is factored out. As a consequence the reduced matrix allows a simple evaluation of the Wilson fermion determinant for any value of the chemical potential and hence the exact projection to the canonical partition functions.

preprint2009arXiv

Neutron electric polarizability

We use the background field method to extract the "connected" piece of the neutron electric polarizability. We present results for quenched simulations using both clover and Wilson fermions and discuss our experience in extracting the mass shifts and the challenges we encountered when we lowered the quark mass. For the neutron we find that as the pion mass is lowered below $500\MeV$, the polarizability starts rising in agreement with predictions from chiral perturbation theory. For our lowest pion mass, $m_π=320\MeV$, we find that $α_n = 3.8(1.3)\times 10^{-4}\fm^3$, which is still only one third of the experimental value. We also present results for the neutral pion; we find that its polarizability turns negative for pion masses smaller than $500\MeV$ which is puzzling.

preprint2008arXiv

Dominance of Sign Geometry and the Homogeneity of the Fundamental Topological Structure

We propose and support the possibility that the shape of topological density 2-point function in pure-glue QCD is crucially, and possibly entirely, determined by the space-time folding (geometry) of the double-sheet sign-coherent structure of Ref.[1], while the distribution of topological density within individual sheets only determines the overall magnitude of the correlator at finite physical distances. A specific manifestation of this, discussed here, is that the shape of the correlation function (encoding e.g. the masses of pseudoscalar glueballs) is reproduced upon the replacement q(x) -> sgn(q(x)), i.e. by considering the double sheet of the same space-time geometry but with constant magnitude of topological density. Combined with previous results on the fundamental topological structure, this suggests that a collective degree of freedom describing topological fluctuations of QCD vacuum can be viewed as a global space-filling homogeneous double membrane. Selected possibilities for practical uses of this are discussed.

preprint2008arXiv

The background field method on the lattice

The background field method has been used successfully to determine hadron electromagnetic polarizabilities. Recently questions have been raised regarding the proper way to deal with the electric field on the lattice. In this paper, we show how the presence of a background electric field affects the quark hopping matrix. We use this formulation to carry out simulations on quenched configurations, and we present the results for neutron electric polarizability for pion masses as low as $500\MeV$. We find that the polarizability is roughly constant for the quark masses considered, $α\sim 1.5 \times 10^{-4} \fm^3$. While the polarizability is positive, it is significantly smaller than the experimental value due to the fact that the quark masses used are too large.

preprint2008arXiv

Winding number expansion for the canonical approach to finite density simulations

The canonical partition function approach was designed to avoid the overlap problem that affects the lattice simulations of nuclear matter at high density. The method employs the projections of the quark determinant on a fix quark number sector. When the quark number is large, the evaluation of the projected determinant becomes numerically unstable. In this paper a different evaluation method based on expanding the determinant in terms of loops winding around the lattice is studied. We show that this method is stable and significantly faster than our original algorithm. This greatly expands the range of quark numbers that we can simulate effectively.

Andrei Alexandru

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

30 published item(s)

Nucleon Electric Dipole Moment from the $θ$ Term with Lattice Chiral Fermions

Anderson Metal-to-Critical Transition in QCD

Applications of Machine Learning to Lattice Quantum Field Theory

Pole position of the $a_1(1260)$ resonance in a three-body unitary framework

The spectrum of qubitized QCD: glueballs in a $S(1080)$ gauge theory

Universality of a truncated sigma-model

Complex Paths Around The Sign Problem

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

Finite-volume energy spectrum of the $K^-K^-K^-$ system

Quantum algorithms for disordered physics

Roper State from Overlap Fermions

A cross-channel study of pion scattering from lattice QCD

Structure Factors of The Unitary Gas Under Supernova Conditions

Dynamical Local Chirality and Chiral Symmetry Breaking

Chiral Polarization Scale at Finite Temperature

Sea Contributions to Hadron Electric Polarizabilities through Reweighting

$ρ$ meson decay on asymmetrical lattices

Absolute X-distribution and self-duality

Critical point of $N_f = 3$ QCD from lattice simulations in the canonical ensemble

Efficient implementation of the overlap operator on multi-GPUs

Electric polarizability of hadrons with overlap fermions on multi-GPUs

The Analysis of Space-Time Structure in QCD Vacuum II: Dynamics of Polarization and Absolute X-Distribution

Finite density phase transition of QCD with $N_f=4$ and $N_f=2$ using canonical ensemble method

Hadron electric polarizability -- finite volume corrections

Magnetic Moments of Negative-Parity Baryons from Lattice QCD

QCD at non-zero density and canonical partition functions with Wilson fermions

Neutron electric polarizability

Dominance of Sign Geometry and the Homogeneity of the Fundamental Topological Structure

The background field method on the lattice

Winding number expansion for the canonical approach to finite density simulations