BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrei Alexandru

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

ResearcherUnclaimed source record
hep-latnucl-thhep-phhep-thquant-phcond-mat.str-elcond-mat.dis-nnnucl-exastro-ph.HEcond-mat.quant-gashep-exMachine LearningPerformance

Catalog footprint

What is connected

64works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

64 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 $Λ$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

$π$N and strangeness sigma terms at the physical point with chiral fermions

Lattice QCD calculations with chiral fermions of the $π$N sigma term $σ_{πN}$ and strangeness sigma term $σ_{sN}$ including chiral interpolation with continuum and volume corrections are provided in this work, with the excited-state contaminations subtracted properly. We calculate the scalar matrix element for the light/strange quark directly and find $σ_{πN}=45.9(7.4)(2.8)$ MeV, with the disconnected insertion part contributing 20(12)(4)\%, and $σ_{sN}=40.2(11.7)(3.5)$ MeV, which is somewhat smaller than $σ_{πN}$. The ratio of the strange/light scalar matrix elements is $y$ = 0.09(3)(1).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A study of symmetry breaking in a relativistic Bose gas using the contraction algorithm

A relativistic Bose gas at finite density suffers from a sign problem that makes direct numerical simulations not feasible. One possible solution to the sign problem is to re-express the path integral in terms of Lefschetz thimbles. Using this approach we study the relativistic Bose gas both in the symmetric phase (low-density) and the spontaneously broken phase (high-density). In the high-density phase we break explicitly the symmetry and determine the dependence of the order parameter on the breaking. We study the relative contributions of the dominant and sub-dominant thimbles in this phase. We find that the sub-dominant thimble only contributes substantially when the explicit symmetry breaking is small, a regime that is dominated by finite volume effects. In the regime relevant for the thermodynamic limit, this contribution is negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fast Estimator of jacobians in Monte Carlo Integration on Lefschetz Thimbles

A solution to the sign problem is the so-called "Lefschetz thimble approach" where the domain of integration for field variables in the path integral is deformed from the real axis to a sub-manifold in the complex space. For properly chosen sub-manifolds ("thimbles") the sign problem disappears or is drastically alleviated. The parametrization of the thimble by real coordinates require the calculation of a jacobian with a computational cost of order O(V^3), where V is proportional to the spacetime volume. In this note we propose two estimators for this jacobian with a computational cost of order O(V). We discuss analytically the regimes where we expect the estimator to work and show numerical examples in two different models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Glue Spin of the Proton

We report the progress on the lattice QCD calculation of the glue spin contribution to proton spin. This calculation is carried out with valence overlap fermion on 2+1 flavor DWF gauge configurations at two lattice spacings with the momentum of the frame in the kinematic range $0\leq p^2 \leq 2$ GeV$^2$. A mild frame dependence is observed. The matching and mixing with large-momentum effective field theory are in progress. The unrenormalized result at $p^2=4$ GeV$^2$ with $O(a^2)$ correction gives $S_G$ = 0.13(3).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Monte Carlo study of real time dynamics

Monte Carlo studies involving real time dynamics are severely restricted by the sign problem that emerges from highly oscillatory phase of the path integral. In this letter, we present a new method to compute real time quantities on the lattice using the Schwinger-Keldysh formalism via Monte Carlo simulations. The key idea is to deform the path integration domain to a complex manifold where the phase oscillations are mild and the sign problem is manageable. We use the previously introduced "contraction algorithm" to create a Markov chain on this alternative manifold. We substantiate our approach by analyzing the quantum mechanical anharmonic oscillator. Our results are in agreement with the exact ones obtained by diagonalization of the Hamiltonian. The method we introduce is generic and in principle applicable to quantum field theory albeit very slow. We discuss some possible improvements that should speed up the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rho resonance parameters from lattice QCD

We perform a high-precision calculation of the phase shifts for $π$-$π$ scattering in the I = 1, J = 1 channel in the elastic region using elongated lattices with two mass-degenerate quark favors ($N_f = 2$). We extract the $ρ$ resonance parameters using a Breit-Wigner fit at two different quark masses, corresponding to $m_π = 226$MeV and $m_π = 315$MeV, and perform an extrapolation to the physical point. The extrapolation is based on a unitarized chiral perturbation theory model that describes well the phase-shifts around the resonance for both quark masses. We find that the extrapolated value, $m_ρ = 720(1)(15)$MeV, is significantly lower that the physical rho mass and we argue that this shift could be due to the absence of the strange quark in our calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sign problem and Monte Carlo calculations beyond Lefschetz thimbles

We point out that Monte Carlo simulations of theories with severe sign problems can be profitably performed over manifolds in complex space different from the one with fixed imaginary part of the action. We describe a family of such manifolds that interpolate between the tangent space at one critical point, where the sign problem is milder compared to the real plane but in some cases still severe, and the union of relevant thimbles, where the sign problem is mild but a multimodal distribution function complicates the Monte Carlo sampling. We exemplify this approach using a simple 0 + 1 dimensional fermion model previously used on sign problem studies and show that it can solve the model for some parameter values where a solution using Lefshetz thimbles was elusive.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic method with low mode substitution for nucleon isovector matrix elements

We introduce a stochastic sandwich method with low-mode substitution to evaluate the connected three-point functions. The isovector matrix elements of the nucleon for the axial-vector coupling $g_A^3$, scalar couplings $g_S^3$ and the quark momentum fraction $\langle x\rangle_{u -d}$ are calculated with overlap fermion on 2+1 flavor domain-wall configurations on a $24^3 \times 64$ lattice at $m_π = 330$ MeV with lattice spacing $a = 0.114$ fm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Monte Carlo algorithm for simulating fermions on Lefschetz thimbles

A possible solution of the notorious sign problem preventing direct Monte Carlo calculations for systems with non-zero chemical potential is to deform the integration region in the complex plane to a Lefschetz thimble. We investigate this approach for a simple fermionic model. We introduce an easy to implement Monte Carlo algorithm to sample the dominant thimble. Our algorithm relies only on the integration of the gradient flow in the numerically stable direction, which gives it a distinct advantage over the other proposed algorithms. We demonstrate the stability and efficiency of the algorithm by applying it to an exactly solvable fermionic model and compare our results with the analytical ones. We report a very good agreement for a certain region in the parameter space where the dominant contribution comes from a single thimble, including a region where standard methods suffer from a severe sign problem. However, we find that there are also regions in the parameter space where the contribution from multiple thimbles is important, even in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Charm and strange quark masses and $f_{D_s}$ from overlap fermions

We use overlap fermions as valence quarks to calculate meson masses in a wide quark mass range on the $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaborations. The well-defined quark masses in the overlap fermion formalism and the clear valence quark mass dependence of meson masses observed from the calculation facilitate a direct derivation of physical current quark masses through a global fit to the lattice data, which incorporates $O(a^2)$ and $O(m_c^4a^4)$ corrections, chiral extrapolation, and quark mass interpolation. Using the physical masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs, Sommer's scale parameter $r_0$ and the masses of charm quark and strange quark in the $\overline{\rm MS}$ scheme are determined to be $r_0=0.465(4)(9)$ fm, $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.118(6)(24)$ GeV (or $m_c^{\overline{\rm MS}}(m_c)=1.304(5)(20)$ GeV), and $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.101(3)(6)\,{\rm GeV}$, respectively. Furthermore, we observe that the mass difference of the vector meson and the pseudoscalar meson with the same valence quark content is proportional to the reciprocal of the square root of the valence quark masses. The hyperfine splitting of charmonium, $M_{J/ψ}-M_{η_c}$, is determined to be 119(2)(7) MeV, which is in good agreement with the experimental value. We also predict the decay constant of $D_s$ to be $f_{D_s}=254(2)(4)$ MeV. The masses of charmonium $P$-wave states $χ_{c0}, χ_{c1}$ and $h_c$ are also in good agreement with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chiral Symmetry Breaking and Chiral Polarization: Tests for Finite Temperature and Many Flavors

It was recently conjectured that, in SU(3) gauge theories with fundamental quarks, valence spontaneous chiral symmetry breaking is equivalent to condensation of local dynamical chirality and appearance of chiral polarization scale $Λ_{ch}$. Here we consider more general association involving the low-energy layer of chirally polarized modes which, in addition to its width ($Λ_{ch}$), is also characterized by volume density of participating modes ($Ω$) and the volume density of total chirality ($Ω_{ch}$). Few possible forms of the correspondence are discussed, paying particular attention to singular cases where $Ω$ emerges as the most versatile characteristic. The notion of finite-volume "order parameter", capturing the nature of these connections, is proposed. We study the effects of temperature (in N$_f$=0 QCD) and light quarks (in N$_f$=12), both in the regime of possible symmetry restoration, and find agreement with these ideas. In N$_f$=0 QCD, results from several volumes indicate that, at the lattice cutoff studied, the deconfinement temperature $T_c$ is strictly smaller than the overlap-valence chiral transition temperature $T_{ch}$ in real Polyakov line vacuum. Somewhat similar intermediate phase (in quark mass) is also seen in N$_f$=12. It is suggested that deconfinement in N$_f$=0 is related to indefinite convexity of absolute X-distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Classifying the Phases of Gauge Theories by Spectral Density of Probing Chiral Quarks

We describe our recent proposal that distinct phases of gauge theories with fundamental quarks translate into specific types of low-energy behavior in Dirac spectral density. The resulting scenario is built around new evidence substantiating the existence of a phase characterized by bimodal (anomalous) density, and corresponding to deconfined dynamics with broken valence chiral symmetry. We argue that such anomalous phase occurs quite generically in these theories, including in "real world" QCD above the crossover temperature, and in zero-temperature systems with many light flavors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phases of SU(3) Gauge Theories with Fundamental Quarks via Dirac Spectral Density

We propose that, in SU(3) gauge theories with fundamental quarks, confinement can be inferred from spectral density of the Dirac operator. This stems from the proposition that its possible behaviors are exhausted by three distinct types (Fig.1). The monotonic cases are standard and entail confinement with valence chiral symmetry breaking (A) or the lack of both (C,C'). The bimodal (anomalous) option (B) was frequently regarded as an artifact (lattice or other) in previous studies, but we show for the first time that it persists in the continuum limit, and conclude that it informs of a non-confining phase with broken valence chiral symmetry. This generalization rests on the following. $(α)$ We show that bimodality in $N_f$=0 theory past deconfinement temperature $T_c$ is stable with respect to removal of both infrared and ultraviolet cutoffs, indicating that anomalous phase is not an artifact. $(β)$ We demonstrate that transition to bimodality in $N_f$=0 is simultaneous with the loss of confinement: anomalous phase occurs for $T_c < T < T_{ch}$, where $T_{ch}$ is the valence chiral restoration temperature. $(γ)$ Evidence is presented for thermal anomalous phase in $N_f$=2+1 QCD at physical quark masses, whose onset too coincides with the conventional "crossover $T_c$''. We conclude that the anomalous regime $T_c < T < T_{ch}$ is very likely a feature of nature's strong interactions. $(δ)$ Our past studies of zero-temperature $N_f$=12 theories revealed that bimodality also arises via purely light-quark effects. As a result, we expect to encounter anomalous phase on generic paths to valence chiral restoration. We predict its existence also for $N_f$ massless flavors ($T=0$) in the range $N_f^c < N_f < N_f^{ch}$, where $N_f^c$ could be quite low. Conventional arguments would associate $N_f^{ch}$ with the onset of conformal window.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pion electric polarizability from lattice QCD

Electromagnetic polarizabilities are important parameters for understanding the interaction between photons and hadrons. For pions these quantities are poorly constrained experimentally since they can only be measured indirectly. New experiments at CERN and Jefferson Lab are planned that will measure the polarizabilities more precisely. Lattice QCD can be used to compute these quantities directly in terms of quark and gluons degrees of freedom, using the background field method. We present results for the electric polarizability for two different quark masses, light enough to connect to chiral perturbation theory. These are currently the lightest quark masses used in polarizability studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resonance Parameters for the rho-meson from Lattice QCD

We present a calculation of the phase-shift for $π$-$π$ scattering in isospin-1, spin-1 channel in the elastic region. The mass and width of the rho resonance is extracted by fitting these phaseshifts. To vary the scattering momentum we employ asymmetric boxes. We use $N_f = 2$ nHYP-smeared clover fermions and generate two sets of ensembles with pion masses about 315MeV and 227MeV. To determine the phase shifts we compute the energy spectrum both for states at rest and boosted. We employ a variational analysis with interpolating fields including several $q$-$\bar{q}$ and $π$-$π$ interpolating fields with different scattering momenta.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Broken Valence Chiral Symmetry and Chiral Polarization of Dirac Spectrum in N$_f$=12 QCD at Small Quark Mass

The validity of recently proposed equivalence between valence spontaneous chiral symmetry breaking (vSChSB) and chiral polarization of low energy Dirac spectrum (ChP) in SU(3) gauge theory, is examined for the case of twelve mass-degenerate fundamental quark flavors. We find that the vSChSB-ChP correspondence holds for regularized systems studied. Moreover, our results suggest that vSChSB occurs in two qualitatively different circumstances: there is a quark mass $m_c$ such that for $m > m_c$ the mode condensing Dirac spectrum exhibits standard monotonically increasing density, while for $m_{ch} < m < m_c$ the peak around zero separates from the bulk of the spectrum, with density showing a pronounced depletion at intermediate scales. Valence chiral symmetry restoration may occur at yet smaller masses $m < m_{ch}$, but this has not yet been seen by overlap valence probe, leaving the $m_{ch}=0$ possibility open. The latter option could place massless N$_f$=12 theory outside of conformal window. Anomalous behavior of overlap Dirac spectrum for $m_{ch} < m < m_c$ is qualitatively similar to one observed previously in zero and few-flavor theories as an effect of thermal agitation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deconfinement, Chiral Symmetry Breaking and Chiral Polarization

We examine the feasibility of the proposition that there is a temperature range T$_c$ < T < T$_{ch}$ in N$_f$=0 QCD, where real Polyakov line (deconfined) vacuum exhibits valence spontaneous chiral symetry breaking and dynamical chiral polarization of Dirac eigenmodes. Detailed finite-volume analysis convincingly demonstrates the existence of such phase at fixed cutoff (a=0.085 fm). Moreover, it is found that this behavior also takes place closer to the continuum limit (a=0.060 fm) without qualitative change in its properties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distribution of Canonical Determinants in QCD

The distribution of canonical determinants in QCD is determined by means of chiral perturbation theory. For a non-zero quark charge the canonical determinants take complex values. In the dilute pion gas approximation, we compute all moments of the magnitude of the canonical determinants, as well as the first nonvanishing moments of the real and imaginary parts. The non-trivial cancellation between the real and the imaginary parts of the canonical determinants is derived and the signal to noise ratio is discussed. The analytical distributions are compared to lattice data. The average density of the magnitude of the canonical determinants is determined as well and is shown to be given by a variant of the log-normal distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electric polarizability of neutral hadrons from dynamical lattice QCD ensembles

We present a valence calculation of the electric polarizability of the neutron, neutral pion, and neutral kaon on two dynamically generated nHYP-clover ensembles. The pion masses for these ensembles are 227(2) MeV and 306(1) MeV, which are the lowest ones used in polarizability studies. This is part of a program geared towards determining these parameters at the physical point. We carry out a high statistics calculation that allows us to: (1) perform an extrapolation of the kaon polarizability to the physical point; we find $α_K =0.269(43)\times10^{-4}$fm$^{3}$, (2) quantitatively compare our neutron polarizability results with predictions from $χ$PT, and (3) analyze the dependence on both the valence and sea quark masses. The kaon polarizability varies slowly with the light quark mass and the extrapolation can be done with high confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sea contributions to the electric polarizability of the hadrons

We present a lattice QCD calculation of the polarizability of the neutron and other neutral hadrons that includes the effects of the background field on the sea quarks. This is done by perturbatively reweighting the charges of the sea quarks to couple them to the background field. The main challenge in such a calculation is stochastic estimation of the weight factors, and we discuss the difficulties in this estimation. Here we use an extremely aggressive dilution scheme to reduce the stochastic noise to a manageable level. The pion mass in our calculation is 306 MeV and the lattice size is 3 fm. For neutron, we find that $α_E = 2.70(55) * 10^{-4}fm^3$, which is the most precise lattice QCD determination of the polarizability to date that includes sea effects.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Update on the sea contributions to hadron polarizabilities via reweighting

We have made significant progress on extending lattice QCD calculation of the polarizability of the neutron and other hadrons to include the effects of charged dynamical quarks. This is done by perturbatively reweighting the charges of the sea quarks to couple them to the background field. The dominant challenge in such a calculation is stochastic estimation of the weight factors, and we discuss the difficulties in this estimation. Here we use an extremely aggressive dilution scheme with N = 124,416 sources per configuration to reduce the stochastic noise to a manageable level. We find that α_E = 2.70(55) * 10^-4 fm^3 for the neutron on one ensemble. We show that low-mode substitution can be used in tandem with dilution to construct an even better estimator, and introduce the offdiagonal matrix element mapping technique for predicting estimator quality.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

QCD at imaginary chemical potential with Wilson fermions

We investigate the phase diagram in the temperature, imaginary chemical potential plane for QCD with three degenerate quark flavors using Wilson type fermions. While more expensive than the staggered fermions used in past studies in this area, Wilson fermions can be used safely to simulate systems with three quark flavors. In this talk, we focus on the (pseudo)critical line that extends from $μ=0$ in the imaginary chemical potential plane, trace it to the Roberge-Weiss line, and determine its location relative to the Roberge-Weiss transition point. In order to smoothly follow the (pseudo)critical line in this plane we perform a multi-histogram reweighting in both temperature and chemical potential. To perform reweighting in the chemical potential we use the compression formula to compute the determinants exactly. Our results are compatible with the standard scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spontaneous Chiral Symmetry Breaking as Condensation of Dynamical Chirality

The occurrence of spontaneous chiral symmetry breaking (SChSB) is equivalent to sufficient abundance of Dirac near-zeromodes. However, dynamical mechanism leading to breakdown of chiral symmetry should be naturally reflected in chiral properties of the modes. Here we offer such connection, presenting evidence that SChSB in QCD proceeds via the appearance of modes exhibiting dynamical tendency for local chiral polarization. These modes form a band of finite width Lambda_ch (chiral polarization scale) around the surface of otherwise anti--polarized Dirac sea, and condense. Lambda_ch characterizes the dynamics of the breaking phenomenon and can be converted to a quark mass scale, thus offering conceptual means to determine which quarks of nature are governed by broken chiral dynamics. It is proposed that, within the context of SU(3) gauge theories with fundamental Dirac quarks, mode condensation is equivalent to chiral polarization. This makes Lambda_ch an "order parameter" of SChSB, albeit without local dynamical field representation away from chiral limit. Several uses of these features, both at zero and finite temperature, are discussed. Our initial estimates are Lambda_ch~150 MeV (N_f=0), Lambda_ch~80 MeV (N_f=2+1, physical point), and that the strange quark is too heavy to be crucially influenced by broken chiral symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Update on the Sea Contributions to Hadron Electric Polarizabilities through Reweighting

We present the results of a reweighting calculation to compute the contribution of the charged quark sea to the neutron electric polarizability. The chief difficulty is the stochastic estimation of weight factors, and we present a hopping parameter expansion-based technique for reducing the stochastic noise, along with a discussion of why this particular reweighting is so difficult. We used this technique to estimate weight factors for 300 configurations of nHYP-clover fermions and compute the neutron polarizability, but the reweighting greatly inflates the overall statistical error, driven by the stochastic noise in the weight factors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Valence calculation of the electric polarizability on nHYP-Clover ensembles

We present preliminary calculations for the electric polarizability of the neutral pion and neutron on three dynamically generated nHYP-Clover ensembles. We use two different pion masses ($m_π \simeq 300$ and 220 MeV) to gauge the chiral behavior. The effects of partial quenching are analyzed by computing a string of partial quenched valence masses for each ensemble. We also analyzed the volume dependence using elongated lattices, where the elongation is in the direction of the electric field.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chiral polarization scale of QCD vacuum and spontaneous chiral symmetry breaking

It has recently been found that dynamics of pure glue QCD supports the low energy band of Dirac modes with local chiral properties qualitatively different from that of a bulk: while bulk modes suppress chirality relative to statistical independence between left and right, the band modes enhance it. The width of such chirally polarized zone - chiral polarization scale Lambda_ch - has been shown to be finite in the continuum limit at fixed physical volume. Here we present evidence that Lambda_ch remains non-zero also in the infinite volume, and is therefore a dynamical scale in the theory. Our experiments in N_f=2+1 QCD support the proposition that the same holds in the massless limit, connecting Lambda_ch to spontaneous chiral symmetry breaking. In addition, our results suggest that thermal agitation in quenched QCD destroys both chiral polarization and condensation of Dirac modes at the same temperature T_ch > T_c.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Resonance parameters of the rho-meson from asymmetrical lattices

We present a lattice QCD calculation of the parameters of the $ρ$ meson decay. The study is carried out on spatially asymmetric boxes using nHYP-smeared clover fermions with two mass-degenerate quark flavors. Our calculations are carried out at a pion mass $m_π=304(2)$ MeV on the set of lattices $V=24^2\times η24\times48$ with $η=1.0,1.25$, and 2.0 with lattice spacing $a=0.1255(7)$ fm. The resonance mass $m_ρ=827(3)(5)$ MeV and coupling constant $g_{ρππ}=6.67(42)$ are calculate using the P-wave scattering phase shifts. We construct a 2$\times$2 correlation matrix to extract the energy of the scattering states and compute the phase shifts using the finite volume formula. By varying the degree of asymmetry, we are able to compute a set of phase shifts that are evenly distributed throughout the spectral region where the $ρ$ decays.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How Self-Dual is QCD?

Vacuum characteristics quantifying dynamical tendency toward self-duality in gauge theories could be used to judge the relevance of classical solutions or the viability of classically motivated vacuum models. Here we decompose the field strength of equilibrium gauge configurations into self-dual and anti-self-dual parts, and apply absolute X-distribution method to the resulting polarization dynamics in order to construct such characteristics. Using lattice regularization and focusing on pure-glue SU(3) gauge theory at zero temperature, we find evidence for positive but very small dynamical tendency for self-duality of vacuum in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spin Polarizabilities on the Lattice

Spin polarizabilities provide information on the internal structure of hadrons in the presence of weak external electromagnetic fields, and are actively studied by Compton scattering experiments. They provide finer detail than the regular polarizabilities since they require space and time-varying fields. Using an effective action in the weak field limit, we have identified methods to isolate each of the physical quantities ($μ, α, β, γ_{E1}, γ_{M1}, γ_{E2}, γ_{M2}$) for spin-1/2 hadrons, both neutral and charged. We also perform a lattice QCD simulation to investigate the feasibility of the effective action approach.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Absolute Measure of Local Chirality and the Chiral Polarization Scale of the QCD Vacuum

The use of the absolute measure of local chirality is championed since it has a uniform distribution for randomly reshuffled chiral components so that any deviations from uniformity in the associated "X-distribution" are directly attributable to QCD-induced dynamics. We observe a transition in the qualitative behavior of this absolute X-distribution of low-lying eigenmodes which, we propose, defines a chiral polarization scale of the QCD vacuum.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Lattice QCD at finite density via a new canonical approach

We carry out a finite density calculation based on a canonical approach which is designed to address the overlap problem. Two degenerate flavor simulations are performed using Wilson gauge action and Wilson fermions on $4^4$ lattices, at temperatures close to the critical temperature $T_c\approx 170\MeV$ and large densities (5 to 20 times nuclear matter density). In this region, we find that the algorithm works well. We compare our results with those from other approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Locality and Scaling of Quenched Overlap Fermions

The overlap fermion offers the tremendous advantage of exact chiral symmetry on the lattice, but is numerically intensive. This can be made affordable while still providing large lattice volumes, by using coarse lattice spacing, given that good scaling and localization properties are established. Here, using overlap fermions on quenched Iwasaki gauge configurations, we demonstrate directly that the overlap Dirac operator's range is comfortably small in lattice units for each of the lattice spacings 0.20 fm, 0.17 fm, and 0.13 fm (and scales to zero in physical units in the continuum limit). In particular, our direct results contradict recent speculation that an inverse lattice spacing of $1 {\rm GeV}$ is too low to have satisfactory localization. Furthermore, hadronic masses (available on the two coarser lattices) scale very well.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Improved Measure of Local Chirality

It is popular to probe the structure of the QCD vacuum indirectly by studying individual fermion eigenmodes, because this provides a natural way to filter out UV fluctuations. The double-peaking in the distribution of the local chiral orientation parameter (X) has been offered as evidence, by some, in support of a particular model of the vacuum. Here we caution that the X-distribution peaking varies significantly with various versions of the definition of X. Furthermore, each distribution varies little from that resulting from a random reshuffling of the left-handed (and independently the right-handed) fields, which destroys any QCD-induced left-right correlation; that is, the double-peaking is mostly a phase-space effect. We propose a new universal definition of the X parameter whose distribution is uniform for randomly reshuffled fields. Any deviations from uniformity for actual data can then be directly attributable to QCD-induced dynamics. We find that the familiar double peak disappears.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Progress on a canonical finite density algorithm

We test the finite density algorithm in the canonical ensemble which combines the HMC update with the accept/reject step according to the ratio of the fermion number projected determinant to the unprojected one as a way of avoiding the determinant fluctuation problem. We report our preliminary results on the Polyakov loop in different baryon number sectors which exhibit deconfinement transitions on small lattices. The largest density we obtain around $T_c$ is an order of magnitude larger than that of nuclear matter. From the conserved vector current, we calculate the quark number and verify that the mixing of different baryon sectors is small.

0 citations0 reviewsNo code linkedOpen access