Chaotic thermalization in Yang-Mills-Higgs theory on a spacial lattice

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

REDUCTION OF TODA LATTICE HIERARCHY TO GENERALIZED KDV HIERARCHIES AND THE 2-MATRIX MODEL

Toda lattice hierarchy and the associated matrix formulation of the 2M-boson KP hierarchies provide a framework for the Drinfeld-Sokolov reduction scheme realized through Hamiltonian action within the second KP Poisson bracket. By working with free currents, which Abelianize the second KP Hamiltonian structure, we are able to obtain a unified formalism for the reduced SL(M + 1, M - k) KdV hierarchies interpolating between the ordinary KP and KdV hierarchies. The corresponding Lax operators are given as superdeterminants of graded SL(M + 1, M - k) matrices in the diagonal gauge and we describe their bracket structure and field content. In particular, we provide explicit free field representations of the associated W(M, M - k) Poisson bracket algebras generalising the familiar nonlinear W-M+1 algebra. Discrete Backlund transformations for SL(M + 1, M - k) KdV are generated naturally from lattice translations in the underlying Toda-like hierarchy. As an application we demonstrate the equivalence of the two-matrix string model to the SL(M + 1, 1) KdV hierarchy.

Spinodal decomposition in pure-gauge QCD

Spinodal decomposition in a model of pure-gauge SU(2) theory that incorporates a deconfinement phase transition is investigated by means of real-time lattice simulations of the fully nonlinear Ginzburg-Landau equation. Results are compared with a Glauber dynamical evolution using Monte Carlo simulations of pure-gauge lattice QCD. © 2005 American Institute of Physics.

Charmed-meson scattering on nucleons in a QCD Coulomb gauge quark model

The scattering of charmed mesons on nucleons is investigated within a chiral quark model inspired on the QCD Hamiltonian in Coulomb gauge. The microscopic model incorporates a longitudinal Coulomb confining interaction derived from a self-consistent quasi-particle approximation to the QCD vacuum, and a traverse hyperfine interaction motivated from lattice simulations of QCD in Coulomb gauge. From the microscopic interactions at the quark level, effective meson-baryon interactions are derived using a mapping formalism that leads to quark-Born diagrams. As an application, the total cross-section of heavy-light D-mesons scattering on nucleons is estimated.

Fractional fermion charges induced by axial-vector and vector gauge potentials and parity anomaly in planar graphenelike structures

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

No-pole condition in Landau gauge: Properties of the Gribov ghost form factor and a constraint on the 2d gluon propagator

We study general properties of the Landau-gauge Gribov ghost form factor sigma(p(2)) for SU(N-c) Yang-Mills theories in the d-dimensional case. We find a qualitatively different behavior for d = 3, 4 with respect to the d = 2 case. In particular, considering any (sufficiently regular) gluon propagator D(p(2)) and the one-loop-corrected ghost propagator, we prove in the 2d case that the function sigma(p(2)) blows up in the infrared limit p -> 0 as -D(0) ln(p(2)). Thus, for d = 2, the no-pole condition sigma(p(2)) < 1 (for p(2) > 0) can be satisfied only if the gluon propagator vanishes at zero momentum, that is, D(0) = 0. On the contrary, in d = 3 and 4, sigma(p(2)) is finite also if D(0) > 0. The same results are obtained by evaluating the ghost propagator G(p(2)) explicitly at one loop, using fitting forms for D(p(2)) that describe well the numerical data of the gluon propagator in two, three and four space-time dimensions in the SU(2) case. These evaluations also show that, if one considers the coupling constant g(2) as a free parameter, the ghost propagator admits a one-parameter family of behaviors (labeled by g(2)), in agreement with previous works by Boucaud et al. In this case the condition sigma(0) <= 1 implies g(2) <= g(c)(2), where g(c)(2) is a "critical" value. Moreover, a freelike ghost propagator in the infrared limit is obtained for any value of g(2) smaller than g(c)(2), while for g(2) = g(c)(2) one finds an infrared-enhanced ghost propagator. Finally, we analyze the Dyson-Schwinger equation for sigma(p(2)) and show that, for infrared-finite ghost-gluon vertices, one can bound the ghost form factor sigma(p(2)). Using these bounds we find again that only in the d = 2 case does one need to impose D(0) = 0 in order to satisfy the no-pole condition. The d = 2 result is also supported by an analysis of the Dyson-Schwinger equation using a spectral representation for the ghost propagator. Thus, if the no-pole condition is imposed, solving the d = 2 Dyson-Schwinger equations cannot lead to a massive behavior for the gluon propagator. These results apply to any Gribov copy inside the so-called first Gribov horizon; i.e., the 2d result D(0) = 0 is not affected by Gribov noise. These findings are also in agreement with lattice data.

Random lattice models

This thesis deals with three different physical models, where each model involves a random component which is linked to a cubic lattice. First, a model is studied, which is used in numerical calculations of Quantum Chromodynamics.In these calculations random gauge-fields are distributed on the bonds of the lattice. The formulation of the model is fitted into the mathematical framework of ergodic operator families. We prove, that for small coupling constants, the ergodicity of the underlying probability measure is indeed ensured and that the integrated density of states of the Wilson-Dirac operator exists. The physical situations treated in the next two chapters are more similar to one another. In both cases the principle idea is to study a fermion system in a cubic crystal with impurities, that are modeled by a random potential located at the lattice sites. In the second model we apply the Hartree-Fock approximation to such a system. For the case of reduced Hartree-Fock theory at positive temperatures and a fixed chemical potential we consider the limit of an infinite system. In that case we show the existence and uniqueness of minimizers of the Hartree-Fock functional. In the third model we formulate the fermion system algebraically via C*-algebras. The question imposed here is to calculate the heat production of the system under the influence of an outer electromagnetic field. We show that the heat production corresponds exactly to what is empirically predicted by Joule's law in the regime of linear response.

Hadronic correlation functions with quark-disconnected contributions in lattice QCD

One of the fundamental interactions in the Standard Model of particle physicsrnis the strong force, which can be formulated as a non-abelian gauge theoryrncalled Quantum Chromodynamics (QCD). rnIn the low-energy regime, where the QCD coupling becomes strong and quarksrnand gluons are confined to hadrons, a perturbativernexpansion in the coupling constant is not possible.rnHowever, the introduction of a four-dimensional Euclidean space-timernlattice allows for an textit{ab initio} treatment of QCD and provides arnpowerful tool to study the low-energy dynamics of hadrons.rnSome hadronic matrix elements of interest receive contributionsrnfrom diagrams including quark-disconnected loops, i.e. disconnected quarkrnlines from one lattice point back to the same point. The calculation of suchrnquark loops is computationally very demanding, because it requires knowledge ofrnthe all-to-all propagator. In this thesis we use stochastic sources and arnhopping parameter expansion to estimate such propagators.rnWe apply this technique to study two problems which relay crucially on therncalculation of quark-disconnected diagrams, namely the scalar form factor ofrnthe pion and the hadronic vacuum polarization contribution to the anomalousrnmagnet moment of the muon.rnThe scalar form factor of the pion describes the coupling of a charged pion torna scalar particle. We calculate the connected and the disconnected contributionrnto the scalar form factor for three different momentum transfers. The scalarrnradius of the pion is extracted from the momentum dependence of the form factor.rnThe use ofrnseveral different pion masses and lattice spacings allows for an extrapolationrnto the physical point. The chiral extrapolation is done using chiralrnperturbation theory ($chi$PT). We find that our pion mass dependence of thernscalar radius is consistent with $chi$PT at next-to-leading order.rnAdditionally, we are able to extract the low energy constant $ell_4$ from thernextrapolation, and ourrnresult is in agreement with results from other lattice determinations.rnFurthermore, our result for the scalar pion radius at the physical point isrnconsistent with a value that was extracted from $pipi$-scattering data. rnThe hadronic vacuum polarization (HVP) is the leading-order hadronicrncontribution to the anomalous magnetic moment $a_mu$ of the muon. The HVP canrnbe estimated from the correlation of two vector currents in the time-momentumrnrepresentation. We explicitly calculate the corresponding disconnectedrncontribution to the vector correlator. We find that the disconnectedrncontribution is consistent with zero within its statistical errors. This resultrncan be converted into an upper limit for the maximum contribution of therndisconnected diagram to $a_mu$ by using the expected time-dependence of therncorrelator and comparing it to the corresponding connected contribution. Wernfind the disconnected contribution to be smaller than $approx5%$ of thernconnected one. This value can be used as an estimate for a systematic errorrnthat arises from neglecting the disconnected contribution.rn

Up, down, strange and charm quark masses with Nf = 2+1+1 twisted mass lattice QCD

We present a lattice QCD calculation of the up, down, strange and charm quark masses performed using the gauge configurations produced by the European Twisted Mass Collaboration with Nf=2+1+1 dynamical quarks, which include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their physical values. The simulations are based on a unitary setup for the two light quarks and on a mixed action approach for the strange and charm quarks. The analysis uses data at three values of the lattice spacing and pion masses in the range 210–450 MeV, allowing for accurate continuum limit and controlled chiral extrapolation. The quark mass renormalization is carried out non-perturbatively using the RI′-MOM method. The results for the quark masses converted to the scheme are: mud(2 GeV)=3.70(17) MeV, ms(2 GeV)=99.6(4.3) MeV and mc(mc)=1.348(46) GeV. We obtain also the quark mass ratios ms/mud=26.66(32) and mc/ms=11.62(16). By studying the mass splitting between the neutral and charged kaons and using available lattice results for the electromagnetic contributions, we evaluate mu/md=0.470(56), leading to mu=2.36(24) MeV and md=5.03(26) MeV.

Interfaces, strings, and a soft mode in the square lattice quantum dimer model

The quantum dimer model on the square lattice is a U(1) gauge theory that addresses aspects of the physics of high-Tc superconductors. Using a quantum Monte Carlo method, we show that the theory exists in a confining columnar valence bond solid phase. The interfaces separating distinct columnar phases display plaquette order, which, however, is not realized as a bulk phase. Static “electric” charges are confined by flux tubes that consist of multiple strands, each carrying a fractionalized flux ¼. A soft pseudo-Goldstone mode (which becomes exactly massless at the Rokhsar-Kivelson point) extends deep into the columnar phase, with potential implications for high-Tc physics.

Standard model cross-over on the lattice

With the physical Higgs mass the standard model symmetry restoration phase transition is a smooth cross-over. We study the thermodynamics of the cross-over using numerical lattice Monte Carlo simulations of an effective SU(2)×U(1) gauge+Higgs theory, significantly improving on previously published results. We measure the Higgs field expectation value, thermodynamic quantities like pressure, energy density, speed of sound and heat capacity, and screening masses associated with the Higgs and Z fields. While the cross-over is smooth, it is very well defined with a width of only ∼5  GeV. We measure the cross-over temperature from the maximum of the susceptibility of the Higgs condensate, with the result Tc=159.5±1.5  GeV. Outside of the narrow cross-over region the perturbative results agree well with nonperturbative ones.

Nucleon and pion structure with lattice QCD simulations at physical value of the pion mass

We present results on the nucleon scalar, axial, and tensor charges as well as on the momentum fraction, and the helicity and transversity moments. The pion momentum fraction is also presented. The computation of these key observables is carried out using lattice QCD simulations at a physical value of the pion mass. The evaluation is based on gauge configurations generated with two degenerate sea quarks of twisted mass fermions with a clover term. We investigate excited states contributions with the nucleon quantum numbers by analyzing three sink-source time separations. We find that, for the scalar charge, excited states contribute significantly and to a less degree to the nucleon momentum fraction and helicity moment. Our result for the nucleon axial charge agrees with the experimental value. Furthermore, we predict a value of 1.027(62) in the MS¯¯¯¯¯ scheme at 2 GeV for the isovector nucleon tensor charge directly at the physical point. The pion momentum fraction is found to be ⟨x⟩π±u−d=0.214(15)(+12−9) in the MS¯¯¯¯¯ at 2 GeV.

An improved single-plaquette gauge action

We describe and test a nonperturbatively improved single-plaquette lattice action for 4-d SU(2) and SU(3) pure gauge theory, which suppresses large fluctuations of the plaquette, without requiring the naive continuum limit for smooth fields. We tune the action parameters based on torelon masses in moderate cubic physical volumes, and investigate the size of cut-off effects in other physical quantities, including torelon masses in asymmetric spatial volumes, the static quark potential, and gradient flow observables. In 2-d O(N) models similarly constructed nearest-neighbor actions have led to a drastic reduction of cut-off effects, down to the permille level, in a wide variety of physical quantities. In the gauge theories, we find significant reduction of lattice artifacts, and for some observables, the coarsest lattice result is very close to the continuum value. We estimate an improvement factor of 40 compared to using the Wilson gauge action to achieve the same statistical accuracy and suppression of cut-off effects.

Comparison of different lattice definitions of the topological charge

We present a comparison of different definitions of the topological charge on the lattice, using a small-volume ensemble with 2 flavours of dynamical twisted mass fermions. The investigated definitions are: index of the overlap Dirac operator, spectral projectors, spectral flow of the HermitianWilson- Dirac operator and field theoretic with different kinds of smoothing of gauge fields (HYP and APE smearings, gradient flow, cooling). We also show some results on the topological susceptibility.

Aspects of Chiral Symmetry Breaking in Lattice QCD

Thesis (Ph.D.)--University of Washington, 2016-08