Source author record

Y. Meurice

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

36 published item(s)

preprint2020arXiv

Quantum Joule Expansion of One-Dimensional Systems

We investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time $t = 0$. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to nonnegative temperatures in the thermodynamic limit for bosons, but is a true thermodynamic effect for fermions, which is confirmed by finite temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

preprint2019arXiv

$B$-meson semileptonic form factors on (2+1+1)-flavor HISQ ensembles

We report updates to an ongoing lattice-QCD calculation of the form factors for the semileptonic decays $B \to π\ell ν$, $B_s \to K \ell ν$, $B \to π\ell^+ \ell^-$, and $B \to K \ell^+ \ell^-$. The tree-level decays $B_{(s)} \to π(K) \ell ν$ enable precise determinations of the CKM matrix element $|V_{ub}|$, while the flavor-changing neutral-current interactions $B \to π(K) \ell^+ \ell^-$ are sensitive to contributions from new physics. This work uses MILC's (2+1+1)-flavor HISQ ensembles at approximate lattice spacings between $0.057$ and $0.15$ fm, with physical sea-quark masses on four out of the seven ensembles. The valence sector is comprised of a clover $b$ quark (in the Fermilab interpretation) and HISQ light and $s$ quarks. We present preliminary results for the form factors $f_0$, $f_+$, and $f_T$, including studies of systematic errors.

preprint2016arXiv

$B\to Kl^+l^-$ decay form factors from three-flavor lattice QCD

We compute the form factors for the $B \to Kl^+l^-$ semileptonic decay process in lattice QCD using gauge-field ensembles with 2+1 flavors of sea quark, generated by the MILC Collaboration. The ensembles span lattice spacings from 0.12 to 0.045 fm and have multiple sea-quark masses to help control the chiral extrapolation. The asqtad improved staggered action is used for the light valence and sea quarks, and the clover action with the Fermilab interpretation is used for the heavy $b$ quark. We present results for the form factors $f_+(q^2)$, $f_0(q^2)$, and $f_T(q^2)$, where $q^2$ is the momentum transfer, together with a comprehensive examination of systematic errors. Lattice QCD determines the form factors for a limited range of $q^2$, and we use the model-independent $z$ expansion to cover the whole kinematically allowed range. We present our final form-factor results as coefficients of the $z$ expansion and the correlations between them, where the errors on the coefficients include statistical and all systematic uncertainties. We use this complete description of the form factors to test QCD predictions of the form factors at high and low $q^2$. We also compare a Standard-Model calculation of the branching ratio for $B \to Kl^+l^-$ with experimental data.

preprint2016arXiv

Tensor RG calculations and quantum simulations near criticality

We discuss the reformulation of the O(2) model with a chemical potential and the Abelian Higgs model on a 1+1 dimensional space-time lattice using the Tensor Renormalization Group (TRG) method. The TRG allows exact blocking and connects smoothly the classical Lagrangian approach to the quantum Hamiltonian approach. We calculate the entanglement entropy in the superfluid phase of the O(2) model and show that it approximately obeys the logarithmic Calabrese-Cardy scaling obtained from Conformal Field Theory (CFT). We calculate the Polyakov loop in the Abelian Higgs model and discuss the possibility of a deconfinement transition at finite volume. We propose Bose-Hubbard Hamiltonians implementable on optical lattices as quantum simulators for CFT models.

preprint2015arXiv

$|V_{ub}|$ from $B\toπ\ellν$ decays and (2+1)-flavor lattice QCD

We present a lattice-QCD calculation of the $B\toπ\ellν$ semileptonic form factors and a new determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad 2+1-flavor lattice configurations at four lattice spacings and light-quark masses down to 1/20 of the physical strange-quark mass. We extrapolate the lattice form factors to the continuum using staggered chiral perturbation theory in the hard-pion and SU(2) limits. We employ a model-independent $z$ parameterization to extrapolate our lattice form factors from large-recoil momentum to the full kinematic range. We introduce a new functional method to propagate information from the chiral-continuum extrapolation to the $z$ expansion. We present our results together with a complete systematic error budget, including a covariance matrix to enable the combination of our form factors with other lattice-QCD and experimental results. To obtain $|V_{ub}|$, we simultaneously fit the experimental data for the $B\toπ\ellν$ differential decay rate obtained by the BaBar and Belle collaborations together with our lattice form-factor results. We find $|V_{ub}|=(3.72\pm 0.16)\times 10^{-3}$ where the error is from the combined fit to lattice plus experiments and includes all sources of uncertainty. Our form-factor results bring the QCD error on $|V_{ub}|$ to the same level as the experimental error. We also provide results for the $B\toπ\ellν$ vector and scalar form factors obtained from the combined lattice and experiment fit, which are more precisely-determined than from our lattice-QCD calculation alone. These results can be used in other phenomenological applications and to test other approaches to QCD.

preprint2015arXiv

$B\toπ\ell\ell$ form factors for new-physics searches from lattice QCD

The rare decay $B\toπ\ell^+\ell^-$ arises from $b\to d$ flavor-changing neutral currents and could be sensitive to physics beyond the Standard Model. Here, we present the first $ab$-$initio$ QCD calculation of the $B\toπ$ tensor form factor $f_T$. Together with the vector and scalar form factors $f_+$ and $f_0$ from our companion work [J. A. Bailey $et~al.$, Phys. Rev. D 92, 014024 (2015)], these parameterize the hadronic contribution to $B\toπ$ semileptonic decays in any extension of the Standard Model. We obtain the total branching ratio ${\text{BR}}(B^+\toπ^+μ^+μ^-)=20.4(2.1)\times10^{-9}$ in the Standard Model, which is the most precise theoretical determination to date, and agrees with the recent measurement from the LHCb experiment [R. Aaij $et~al.$, JHEP 1212, 125 (2012)]. Note added: after this paper was submitted for publication, LHCb announced a new measurement of the differential decay rate for this process [T. Tekampe, talk at DPF 2015], which we now compare to the shape and normalization of the Standard-Model prediction.

preprint2015arXiv

The fine structure of the entanglement entropy in the classical XY model

We compare two calculations of the particle density in the superfluid phase of the classical XY model with a chemical potential $μ$ in 1+1 dimensions.The first relies on exact blocking formulas from the Tensor Renormalization Group (TRG) formulation of the transfer matrix. The second is a worm algorithm. We show that the particle number distributions obtained with the two methods agree well. We use the TRG method to calculate the thermal entropy and the entanglement entropy. We describe the particle density, the two entropies and the topology of the world lines as we increase $μ$ to go across the superfluid phase between the first two Mott insulating phases. For a sufficiently large temporal size, this process reveals an interesting fine structure: the average particle number and the winding number of most of the world lines in the Euclidean time direction increase by one unit at a time. At each step, the thermal entropy develops a peak and the entanglement entropy increases until we reach half-filling and then decreases in a way that approximately mirror the ascent. This suggests an approximate fermionic picture.

preprint2014arXiv

$B\toπ\ellν$ semileptonic form factors from unquenched lattice QCD and determination of $|V_{ub}|$

We compute the $B\toπ\ellν$ semileptonic form factors and update the determination of the CKM matrix element $|V_{ub}|$. We use the MILC asqtad ensembles with $N_f=2+1$ sea quarks at four different lattice spacings in the range $a \approx 0.045$~fm to $0.12$~fm. The lattice form factors are extrapolated to the continuum limit using SU(2) staggered chiral perturbation theory in the hard pion limit, followed by an extrapolation in $q^2$ to the full kinematic range using a functional $z$-parameterization. The extrapolation is combined with the experimental measurements of the partial branching fraction to extract $|V_{ub}|$. Our preliminary result is $|V_{ub}|=(3.72\pm 0.14)\times 10^{-3}$, where the error reflects both the lattice and experimental uncertainties, which are now on par with each other.

preprint2014arXiv

Progress towards quantum simulating the classical O(2) model

We connect explicitly the classical $O(2)$ model in 1+1 dimensions, a model sharing important features with $U(1)$ lattice gauge theory, to physical models potentially implementable on optical lattices and evolving at physical time. Using the tensor renormalization group formulation, we take the time continuum limit and check that finite dimensional projections used in recent proposals for quantum simulators provide controllable approximations of the original model. We propose two-species Bose-Hubbard models corresponding to these finite dimensional projections at strong coupling and discuss their possible implementations on optical lattices using a $^{87}$Rb and $^{41}$K Bose-Bose mixture.

preprint2013arXiv

$B\toπ\ellν$ and $B\toπ\ell^+\ell^-$ semileptonic form factors from unquenched lattice QCD

We update the lattice calculation of the $B\toπ$ semileptonic form factors, which have important applications to the CKM matrix element $|V_{ub}|$ and the $B\toπ\ell^+\ell^-$ rare decay. We use MILC asqtad ensembles with $N_f=2+1$ sea quarks and over a range of lattice spacings $a \approx 0.045$--$0.12$ fm. We perform a combined chiral and continuum extrapolation of our lattice data using SU(2) staggered chiral perturbation theory in the hard pion limit. To extend the results for the form factors to the full kinematic range, we take a functional approach to parameterize the form factors using the Bourrely-Caprini-Lellouch formalism in a model-independent way. Our analysis is still blinded with an unknown off-set factor which will be disclosed when we present the final results.

preprint2013arXiv

Accurate exponents from approximate tensor renormalizations

We explain the recent numerical successes obtained by Tao Xiang's group, who developed and applied Tensor Renormalization Group methods for the Ising model on square and cubic lattices, by the fact that their new truncation method sharply singles out a surprisingly small subspace of dimension two. We show that in the two-state approximation, their transformation can be handled analytically yielding a value 0.964 for the critical exponent nu much closer to the exact value 1 than 1.338 obtained in the Migdal-Kadanoff approximation. We propose two alternative blocking procedures that preserve the isotropy and improve the accuracy to nu=0.987 and 0.993 respectively. We discuss applications to other classical lattice models, including models with fermions, and suggest that it could become a competitor for Monte Carlo methods suitable to calculate accurately critical exponents, take continuum limits and study near-conformal systems in arbitrarily large volumes.

preprint2013arXiv

Controlling sign problems in spin models using tensor renormalization

We consider the sign problem for classical spin models at complex $β=1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$β$ than the reweighting Monte Carlo method. For the Ising model with complex $β$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined precisely with the TRG method. We check the convergence of the TRG method for the O(2) model on $L\times L$ lattices when the number of states $D_s$ increases. We show that the finite size scaling of the calculated Fisher zeros agrees very well with the Kosterlitz-Thouless transition assumption and predict the locations for larger volume. The location of these zeros agree with Monte Carlo reweighting calculation for small volume. The application of the method for the O(2) model with a chemical potential is briefly discussed.

preprint2013arXiv

Exact blocking formulas for spin and gauge models

Using the example of the two-dimensional (2D) Ising model, we show that in contrast to what can be done in configuration space, the tensor renormalization group (TRG) formulation allows one to write exact, compact, and manifestly local blocking formulas and exact coarse grained expressions for the partition function. We argue that similar results should hold for most models studied by lattice gauge theorists. We provide exact blocking formulas for several 2D spin models (the O(2) and O(3) sigma models and the SU(2) principal chiral model) and for the 3D gauge theories with groups Z_2, U(1) and SU(2). We briefly discuss generalizations to other groups, higher dimensions and practical implementations.

preprint2013arXiv

Heavy-meson semileptonic decays for the Standard Model and beyond

We calculate the form factors for the semileptonic decays $B_s\to K\ellν$ and $B\to K\ell\ell$ with lattice QCD. We work at several lattice spacings and a range of light quark masses, using the MILC 2+1-flavor asqtad ensembles. We use the Fermilab method for the $b$ quark. We obtain chiral-continuum extrapolations for $E_K$ up to $\sim1.2$ GeV and then extend to the entire kinematic range with the model-independent $z$ expansion.

preprint2013arXiv

Tensor renormalization group study of classical XY model on the square lattice

Using the tensor renormalization group method based on the higher-order singular value decom- position, we have studied the thermodynamic properties of the continuous XY model on the square lattice. The temperature dependence of the free energy, the internal energy and the specific heat agree with the Monte Carlo calculations. From the field dependence of the magnetic susceptibility, we find the Kosterlitz-Thouless transition temperature to be 0.8921 \pm 0.0019, consistent with the Monte Carlo as well as the high temperature series expansion results. At the transition temperature, the critical exponent δis estimated as 14.5, close to the analytic value by Kosterlitz.

preprint2012arXiv

B_s->D_s/B->D Semileptonic Form-Factor Ratios and Their Application to BR(B^0_s->μ^+μ^-)

We calculate form-factor ratios between the semileptonic decays \bar{B}->D^+\ell^-\barν and \bar{B}_s->D_s^+\ell^-\barν with lattice QCD. These ratios are a key theoretical input in a new strategy to determine the fragmentation fractions of the neutral B decays, which are needed for measurements of BR(B^0_s-> μ^+μ^-). We use the MILC ensembles of gauge configurations with 2+1 flavors of sea quarks at two lattice spacings of approximately 0.12 fm and 0.09 fm. We use the model-independent z-parametrization to extrapolate our simulation results at small recoil toward maximum recoil. Our results for the form-factor ratios are $f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_K) =1.046(44)_{stat.}(15)_{syst.}$ and $f_0^{(s)}(M^2_π)/f_0^{(d)}(M^2_π)=1.054(47)_{stat.}(17)_{syst.}$. In contrast to a QCD sum-rule calculation, no significant departure from U-spin (d<->s) symmetry is observed.

preprint2012arXiv

Comparison of Improved Perturbative Methods

In many cases of interest, the perturbative series based on conventional Feynman diagrams have a zero radius of convergence. Series with a finite radius of convergence can be obtained by either introducing a large field cutoff or by replacing the exponential of the perturbation by a sequence of approximants as recently proposed by Pollet, Prokof'ev, and Svistunov. We compare these two methods for integrals and quantum mechanical problems. The two methods perform well in complementary regime (strong coupling for the large field cutoff and intermediate coupling for the other method). We briefly discuss potential applications for lattice gauge theory with compact groups (which have a build-in large field cutoff).

preprint2012arXiv

Density of states and Fisher's zeros in compact U(1) pure gauge theory

We present high-accuracy calculations of the density of states using multicanonical methods for lattice gauge theory with a compact gauge group U(1) on 4^4, 6^4 and 8^4 lattices. We show that the results are consistent with weak and strong coupling expansions. We present methods based on Chebyshev interpolations and Cauchy theorem to find the (Fisher's) zeros of the partition function in the complex beta=1/g^2 plane. The results are consistent with reweighting methods whenever the latter are accurate. We discuss the volume dependence of the imaginary part of the Fisher's zeros, the width and depth of the plaquette distribution at the value of beta where the two peaks have equal height. We discuss strategies to discriminate between first and second order transitions and explore them with data at larger volume but lower statistics. Higher statistics and even larger lattices are necessary to draw strong conclusions regarding the order of the transition.

preprint2012arXiv

QCD calculations with optical lattices?

By trapping cold polarizable atoms in periodic potentials created by crossed laser beams, it is now possible to experimentally create "clean" lattice systems. Experimentalists have successfully engineered local and nearest-neighbor interactions that approximately recreate Hubbard-like models on table tops. I discuss the possibility of using this new technology in the context of lattice gauge theory and in particular, relativistic dispersion relations, flavor symmetry, functional derivatives and emerging local gauge symmetry.

preprint2012arXiv

Refining new-physics searches in B -> D tau nu decay with lattice QCD

The semileptonic decay channel B -> D tau nu is sensitive to the presence of a scalar current, such as that mediated by a charged-Higgs boson. Recently the BaBar experiment reported the first observation of the exclusive semileptonic decay B -> D tau nu, finding an approximately 2-sigma disagreement with the Standard-Model prediction for the ratio R(D)=BR(B->D tau nu)/BR(B->D l nu), where l=e,mu. We compute this ratio of branching fractions using hadronic form factors computed in unquenched lattice QCD and obtain R(D) = 0.316(12)(7), where the errors are statistical and total systematic, respectively. This result is the first Standard-Model calculation of R(D) from ab initio full QCD. Its error is smaller than that of previous estimates, primarily due to the reduced uncertainty in the scalar form factor f_0(q^2). Our determination of R(D) is approximately 1-sigma higher than previous estimates and, thus, reduces the tension with experiment. We also compute R(D) in models with electrically charged scalar exchange, such as the type II two-Higgs doublet model. Once again, our result is consistent with, but approximately 1-sigma higher than, previous estimates for phenomenologically relevant values of the scalar coupling in the type II model. As a byproduct of our calculation, we also present the Standard-Model prediction for the longitudinal polarization ratio P_L (D)= 0.325(4)(3).

preprint2012arXiv

Remarks about Dyson's instability in the large-N limit

There are known examples of perturbative expansions in the 't Hooft coupling lt with a finite radius of convergence. This seems to contradict Dyson's argument suggesting that the instability at negative coupling implies a zero radius of convergence. Using the example of the linear sigma model in three dimensions, we discuss to which extent the two points of view are compatible. We show that a saddle point persists for negative values of lt until a critical value -|lt_c| is reached. A numerical study of the perturbative series for the renormalized mass confirms an expected singularity of the form (lt +|lt_c|)^1/2. However, for -|lt_c|< lt <0, the effective potential does not exist if phi^2 >phi^2_{max}(lt) and not at all if lt<-|lt_c|. We show that phi^2_{max}(lt) propto 1/|lt | for small negative lt. The finite radius of convergence can be justified if the effective theory is defined with a large field cutoff phi^2_{max}(lt) which provides a quantitative measure of the departure from the original model considered.

preprint2011arXiv

Complex RG flows for 2D nonlinear O(N) sigma models

Motivated by recent attempts to find nontrivial infrared fixed points in 4-dimensional lattice gauge theories, we discuss the extension of the renormalization group (RG) transformations to complex coupling spaces for O(N) models on LxL lattices, in the large-N limit. We explain the Riemann sheet structure and singular points of the finite L mappings between the mass gap and the 't Hooft coupling. We argue that the Fisher's zeros appear on "strings" ending approximately near these singular points. We show that for the spherical model at finite N and L, the density of states is stripwise polynomial in the complex energy plane. We compare finite volume complex flows obtained from the rescaling of the ultraviolet cutoff in the gap equation and from the two-lattice matching. In both cases, the flows are channelled through the singular points and end at the strong coupling fixed points, however strong scheme dependence appear when the Compton wavelength of the mass gap is larger than the linear size of the system. We argue that the Fisher's zeros control the global properties of the complex flows. We briefly discuss the implications for perturbation theory, proofs of confinement and searches for nontrivial infrared fixed points in models beyond the standard model.

preprint2011arXiv

Lines of Fisher's zeros as separatrices for complex renormalization group flows

We extend the renormalization group transformation based on the two-lattice matching to the complex inverse temperature plane for Dyson's hierarchical Ising model. We consider values of the dimensional parameter above, below and exactly at the critical value where the ordered low temperature phase becomes impossible for a real positive temperature. We show numerically that, as the volume increases, the Fisher's zeros appear to accumulate along lines that separate the flows ending on different fixed points. We justify these findings in terms of finite size scaling. We argue that the location of the Fisher's zeros at large volume determine the phase diagram in the complex plane. We discuss the implications for nontrivial infrared fixed points in lattice gauge theory.

preprint2011arXiv

New applications of the renormalization group method in physics -- a brief introduction

The renormalization group method developed by Ken Wilson more than four decades ago has revolutionized the way we think about problems involving a broad range of energy scales such as phase transitions, turbulence, continuum limits and bifurcations in dynamical systems. The theme issue provides articles reviewing recent progress made using the renormalization group method in atomic, condensed matter, nuclear and particle physics. In the following we introduce these articles in a way that emphasizes common themes and the universal aspects of the method.

preprint2011arXiv

Volume Effects in Discrete beta functions

We calculate discrete beta functions corresponding to the two-lattice matching for the 2D O(N) models and Dyson's hierarchical model. We describe and explain finite-size effects such as the appearance of a nontrivial infrared fixed point that goes to infinity at infinite volume or the merging of an infrared and an ultraviolet fixed point. We present extensions of the RG flows to the complex coupling plane. We discuss the possibility of constructing a continuous beta function from the discrete one by using functional conjugation methods. We briefly discuss the relevance of these findings for the search of nontrivial fixed points in multiflavor lattice gauge theory models.

preprint2010arXiv

Fisher's zeros as boundary of renormalization group flows in complex coupling spaces

We propose new methods to extend the renormalization group transformation to complex coupling spaces. We argue that the Fisher's zeros are located at the boundary of the complex basin of attraction of infra-red fixed points. We support this picture with numerical calculations at finite volume for two-dimensional O(N) models in the large-N limit and the hierarchical Ising model. We present numerical evidence that, as the volume increases, the Fisher's zeros of 4-dimensional pure gauge SU(2) lattice gauge theory with a Wilson action, stabilize at a distance larger than 0.15 from the real axis in the complex beta=4/g^2 plane. We discuss the implications for proofs of confinement and searches for nontrivial infra-red fixed points in models beyond the standard model.

preprint2010arXiv

Fisher's zeros as boundary of RG flows in complex coupling space

We discuss the possibility of extending the RG flows to complex coupling spaces. We argue that the Fisher's zeros are located at the boundary of the complex basin of attraction of IR fixed points. We support this picture with numerical calculations at finite volume for2D O(N) models in the large-N limit and the hierarchical Ising model using the two-lattice matching method. We present numerical evidence supporting the idea that, as the volume increases, the Fisher's zeros of 4-dimensional pure gauge SU(2) lattice gauge theory with a Wilson action, stabilize at a distance larger than 0.1 from the real axis in the complex beta=4/g^2 plane. We show that when a positive adjoint term is added, the zeros get closer to the real axis. We compare the situation with the U(1) case. We discuss the implications of this new framework for proofs of confinement and searches for nontrivial IR fixed points in models beyond the standard model.

preprint2009arXiv

Dyson instability for 2D nonlinear O(N) sigma models

For lattice models with compact field integration (nonlinear sigma models over compact manifolds and gauge theories with compact groups) and satisfying some discrete symmetry, the change of sign of the bare coupling g_0^2 at zero results in a mere discontinuity in the average energy rather than the catastrophic instability occurring in theories with integration over arbitrarily large fields. This indicates that the large order of perturbative series and the non-perturbative contributions should have unexpected features. Using the large-N limit of 2-dimensional nonlinear O(N) sigma model, we discuss the complex singularities of the average energy for complex 't Hooft coupling lambda= g_0^2N. A striking difference with the usual situation is the absence of cut along the negative real axis. We show that the zeros of the partition function can only be inside a clover shape region of the complex lambda plane. We calculate the density of states and use the result to verify numerically the statement about the zeros. We propose dispersive representations of the derivatives of the average energy for an approximate expression of the discontinuity. The discontinuity is purely non-perturbative and contributions at small negative coupling in one dispersive representation are essential to guarantee that the derivatives become exponentially small when lambda -> 0^+ We discuss the implications for gauge theories.

preprint2009arXiv

Dyson's Instability in Lattice Gauge Theory

We discuss Dyson's argument that the vacuum is unstable under a change g^2 -> - g^2, in the context of lattice gauge theory. For compact gauge groups, the partition function is well defined at negative g^2, but the average plaquette P has a discontinuity when g^2 changes sign. This reflects a change of vacuum rather than a loss of vacuum. In addition, P has poles in the complex g^2 plane, located at the complex zeros of the partition function (Fisher's zeros). We discuss the relevance of these singularities for lattice perturbation theory. We present new methods to locate Fisher's zeros using numerical values for the density of state in SU(2) and U(1) pure gauge theory. We briefly discuss similar issues for O(N) nonlinear sigma models where the local integrals are also over compact spaces.

preprint2009arXiv

Finite Size Scaling and Universality in SU(2) at Finite Temperature

We study the 4-th Binder cumulant on $4\times{N_σ}^3$ lattices for a pure SU(2) gauge theory. We use 20 data points for a sequence of ${N_σ}$ in $β$ intervals shrinking when ${N_σ}$ increases, in order to reduce the nonlinear effects. Using a log-log fit of the slope versus ${N_σ}$, we obtain the preliminary result $ν=0.637(11)$ in reasonably good agreement with the value for the 3D Ising model universality class. The corrections due to irrelevant directions appear to be dominated by a term proportional to ${N_σ}^{-2.03(4)}$ which seems compatible with the breaking of rotational symmetry.

preprint2008arXiv

Approximate forms of the density of states

We compare MC calculations of the density of states in SU(2) pure gauge theory with the weak and strong coupling expansions. Surprisingly, the range of validity of the two approximations overlap significantly, however the large order behavior of both expansions appear to be similar to the corresponding expansions of the plaquette. We discuss the implications for the calculation of the Fisher's zeros of the partition function.

preprint2008arXiv

Volume dependence of Fisher's zeros

We study the location of the partition function zeros in the complex beta plane (Fisher's Zeros) for SU(2) lattice gauge theory on L^4 lattices. We discuss recent attempts to locate complex zeros for L=4 and 6. We compare results obtained using various polynomial approximations of the logarithm of the density of states and a straightforward MC reweighting. We conclude that the method based on a combination of discrete Chebyshev orthogonality and patching plaquette distributions at different beta provides the more reliable estimates.

preprint2006arXiv

Definition and parametrization of non-perturbative effects in quenched QCD

The notion of a non-perturbative effect is ambiguous if it requires the subtraction of a perturbative part defined by a diverging series. A common procedure consists in dropping the order of minimal contribution and the higher orders. This allows us to isolate very accurately the one-instanton effect for the double-well potential. For the one plaquette gauge theory, an exact analytical expression can be written for the non-perturbative part. We report recent attempts to extend this approach to the average plaquette of quenched QCD. Our goal is to express the non-perturbative effects in terms of expressions of the form beta^B exp(-A beta) calculable semi-classically. The situation is complicated by zeroes of the partition function in the complex $β$ plane (presumably near 5.75 pm i 0.2). We discuss two methods to describe the intermediate and large order behavior of the perturbative series. One is inspired by mean field theory (logarithmic specific heat) and reproduces accurately the known perturbative series with only two free parameters. A diagrammatic interpretation of this fact is still lacking. The other is based on infra-red renormalons with a factorial growth showing up at order larger than 20 and a possible effective theory interpretation. These extrapolations are compatible with the non-perturbative part of the plaquette being proportional to a^4. We propose an exponential parametrization to the corrections to the universal part of the beta function and find results compatible with the suggestion of a^2 corrections made by C. Allton.

preprint1998arXiv

Accurate Checks of Universality for Dyson's Hierarchical Model

Using recently developed methods, we perform high-accuracy calculations of the susceptibility near beta_c for the D=3 version of Dyson's hierarchical model. Using linear fits, we estimate the leading gamma and subleading Delta exponents. Independent estimates are obtained by calculating the first two eigenvalues of the linearized renormalization group transformation. We found gamma = 1.29914073 (with an estimated error of 10^{-8}) and, Delta=0.4259469 (with an estimated error of 10^{-7}) independently of the choice of local integration measure (Ising or Landau-Ginzburg). After a suitable rescaling, the approximate fixed points for a large class of local measure coincide accurately with a fixed point constructed by Koch and Wittwer.

preprint1994arXiv

The High-Temperature Expansion of the Hierarchical Ising Model: From Poincaré Symmetry to an Algebraic Algorithm

We show that the hierarchical model at finite volume has a symmetry group which can be decomposed into rotations and translations as the familiar Poincaré groups. Using these symmetries, we show that the intricate sums appearing in the calculation of the high-temperature expansion of the magnetic susceptibility can be performed, at least up to the fourth order, using elementary algebraic manipulations which can be implemented with a computer. These symmetries appear more clearly if we use the 2-adic fractions to label the sites. We then apply the new algebraic methods to the calculation of quantities having a random walk interpretation. In particular, we show that the probability of returning at the starting point after $m$ steps has poles at $D=-2,-4,....-2m$ , where $D$ is a free parameter playing a role similar to the dimensionality in nearest neighbor models.