Antiferromagnetic Ising model in an imaginary magnetic field

We study the two-dimensional antiferromagnetic Ising model with a purely imaginary magnetic field, which can be thought of as a toy model for the usual $\theta$ physics. Our motivation is to have a benchmark calculation in a system which suffers from a strong sign problem, so that our results can be used to test Monte Carlo methods developed to tackle such problems. We analyze here this model by means of analytical techniques, computing exactly the first eight cumulants of the expansion of the effective Hamiltonian in powers of the inverse temperature, and calculating physical observables for a large number of degrees of freedom with the help of standard multi-precision algorithms. We report accurate results for the free energy density, internal energy, standard and staggered magnetization, and the position and nature of the critical line, which confirm the mean-field qualitative picture, and which should be quantitatively reliable, at least in the high-temperature regime, including the entire critical line

Massive Schwinger model at finite θ\theta

Using the approach developed in [V. Azcoiti, G. Di Carlo, A. Galante, V. Laliena, \textit{Phys. Lett.} \textbf{B563}, (2003) 117], we are able to reconstruct the behavior of the massive 1-flavor Schwinger model with a $\theta$ term and a quantized topological charge. We calculate the full dependence of the order parameter with $\theta$. Our results at $\theta = \pi$ are compatible with Coleman's conjecture on the phase diagram of this model.Comment: 8 pages, 8 figure

Frustration in Lattice Gauge Theory

We introduce a U(1) lattice gauge theory which incorporates explicit frustration in $d>2$. We show, by identifying an appropiate order parameter and through computer simulations, the existence of a frustrated region in the phase diagram of the model. We study this phase diagram and the nature of the transition lines.Comment: 3 pages, Latex, 5 figures, espcrc2.sty needed, Talk presented by E. Follana at LATTICE 9

Critical behaviour of the O(3) nonlinear sigma model with topological term at theta=pi from numerical simulations

We investigate the critical behaviour at theta=pi of the two-dimensional O(3) nonlinear sigma model with topological term on the lattice. Our method is based on numerical simulations at imaginary values of theta, and on scaling transformations that allow a controlled analytic continuation to real values of theta. Our results are compatible with a second order phase transition, with the critical exponent of the SU(2)_1 Wess-Zumino-Novikov-Witten model, for sufficiently small values of the coupling.Comment: Revised version. 24 pages, 7 figure

Non-perturbative physics in lattice gauge theories

A few decades have passed since quantum chromodynamics (QCD) was established as the theory describing strong interactions. It is broadly accepted as one of the most successful theories in modern physics, and it has been extensively tested, both from the theoretical and the experimental perspectives.At high energies, QCD is asymptotically free, which means that its fundamental constituents, quarks and gluons, interact with a strength that decreases as the energy scale reaches higher values. In this regime, it is feasible to use perturbation theory to resolve short distance interactions. On the other hand, for not-so-high energies, the strong interaction cannot be reduced to a converging series of Feynman diagrams. In fact, one of the characteristic properties of QCD is the so-called color-confinement. In this purely non-perturbative regime, there are few techniques that can analyze the theory successfully. Probably the most well-established of them is lattice QCD. Since the foundational work of Wilson in 1974, the success of the lattice approach has been growing consistently over time. Many milestones have already been reached, including precise simulations that account for the effects of virtual quark loops, the determination of the light hadron spectrum with fully controlled systematics or, more recently, the computation of the isospin splittings with great agreement with the experimental data.For the above reasons, QCD is believed to be the correct theory describing strong interactions, both for high and low energies, and lattice QCD is recognized by the community as a trustworthy ab initio approach that has an useful interaction with experiment, paraphrasing Wilson. However, there are some fundamental topics that still constitute open questions. At least two problems share this status: the behavior of matter at finite baryonic density and the studies involving topological effects in QCD. The main difficulty behind the modest progress achieved in both areas is the same: the action of the theory is complex, and there is no known reformulation that can avoid the appearance of a severe sign problem (SSP).In this context, the main part of this thesis has been devoted to study models which suffer from a SSP, such as the two-dimensional Ising model within an imaginary magnetic field or the massive 1-flavor Schwinger model with a theta term. In the first case, we study the well-known model by means of analytical techniques, exploring a region of the parameter space somewhat unattended by the literature, possibly due to the difficulty of applying either analytical or numerical techniques. Secondly, and with the aim of engaging with QCD-like systems with a theta term and to develop further the methods dealing with the SSP, we have studied the massive 1-flavor Schwinger model with a $\theta$ term, which corresponds to QED in $1+1$ dimensions, and is in fact broadly used as its toy model. Moreover, defining the topological charge on the lattice is almost trivial in this model, in contrast with any of the usual definitions of this observable in lattice QCD, which are much more involved. As a byproduct of this line of work, and driven by the necessity of optimizing further our previous algorithms, we have also analysed the 2-flavor version of the Schwinger model. In this case, we have bypassed the computation of the full fermionic determinant by following an approach based on the use of pseudofermions.Finally, beyond the study of systems afflicted by a SSP, another topic within lattice QCD has been treated during the development of this thesis: the strong running coupling alpha_S. Its dependence with the momentum transfer, which encodes the underlying interactions of quarks and gluons in the QCD framework, constitutes a very active field of research, that includes a large variety of approaches. At large momenta, where perturbative QCD can be applied, both experimental and theoretical methods try to provide the most accurate approximation. In this context, lattice-based strategies have been capable of delivering results both in the infrared region and in the high energy regime, where in fact they provide the most precise determination of the coupling constant. Our work can be framed precisely into these approaches that come from lattice QCD, and it relies upon a ghost-gluon vertex computation.<br /

D→π,lνD \rightarrow \pi, l \nu Semileptonic Decays, ∣Vcd∣|V_{cd}| and 2nd^{nd} Row Unitarity from Lattice QCD

We present a new calculation of the $D \rightarrow \pi, l \nu$ semileptonic form factor $f^{D \rightarrow \pi}_+(q^2)$ at $q^2 = 0$ based on HISQ charm and light valence quarks on MILC $N_f = 2 +1$ lattices. Using methods developed recently for HPQCD's study of $D \rightarrow K, l \nu$ decays, we find $f^{D \rightarrow \pi}_+(0) = 0.666(29)$. This signifies a better than factor of two improvement in errors for this quantity compared to previous calculations. Combining the new result with CLEO-c branching fraction data, we extract the CKM matrix element $|V_{cd}| = 0.225(6)_{exp.}(10)_{lat.}$, where the first error comes from experiment and the second from theory. With a total error of $\sim5.3$\% the accuracy of direct determination of $|V_{cd}|$ from $D$ semileptonic decays has become comparable to (and in good agreement with) that from neutrino scattering. We also check for second row unitarity using this new $|V_{cd}|$, HPQCD's earlier $|V_{cs}|$ and $|V_{cb}|$ from the Fermilab Lattice \& MILC collaborations. We find $|V_{cd}|^2 + |V_{cs}|^2 + |V_{cb}|^2 = 0.976(50)$, improving on the current PDG2010 value.Comment: 7 pages, 7 figures, and 4 table