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

Yuzhi Liu

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

ResearcherUnclaimed source record
hep-lathep-thcond-mat.stat-mechhep-phcond-mat.str-elquant-phcond-mat.quant-gascond-mat.supr-conhep-exphysics.comp-ph

Catalog footprint

What is connected

40works
10topics
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

40 published item(s)

preprint2022arXiv

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

The origin of the pseudogap behavior, found in many high-$T_c$ superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a "gap-filling" rather than "gap-closing" behavior, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solving quantum rotor model with different Monte Carlo techniques

We systematically test the performance of several Monte Carlo update schemes for the $(2+1)$d XY phase transition of quantum rotor model. By comparing the local Metropolis (LM), LM plus over-relaxation (OR), Wolff-cluster (WC), hybrid Monte Carlo (HM), hybrid Monte Carlo with Fourier acceleration (FA) scheme, it is clear that among the five different update schemes, at the quantum critical point, the WC and FA schemes acquire the smallest autocorrelation time and cost the least amount of CPU hours in achieving the same level of relative error, and FA enjoys a further advantage of easily implementable for more complicated interactions such as the long-range ones. These results bestow one with the necessary knowledge of extending the quantum rotor model, which plays the role of ferromagnetic/antiferromagnetic critical bosons or Z$_2$ topological order, to more realistic and yet challenging models such as Fermi surface Yukawa-coupled to quantum rotor models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The dynamical exponent of a quantum critical itinerant ferromagnet: a Monte Carlo study

We consider the effect of the coupling between 2D quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions.A common belief is that near a $q=0$ ferromagnetic transition, fermions induce an $Ω/q$ Landau damping of rotors (i.e., the dynamical critical exponent is $z=3$) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy $Σ\propto ω^{2/3}$. This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies.Here we show that for the XY case the behavior is different.We report the results of large scale quantum Monte Carlo simulations,which show that at small frequencies $z=2$ and $Σ\propto ω^{1/2}$. We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces $1/q$ in the Landau damping by a constant. We discuss the implication of these results to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fermion enhanced first-order phase transition and chiral Gross-Neveu tricritical point

The fluctuations of massless Dirac fermion can not only turn a first-order bosonic phase transition (in the Landau sense) to a quantum critical point, but also work reversely to enhance the first-order transition itself, depending on the implementation of finite size effects in the coupling corrections. Here, we report a case study of the latter by employing quantum Monte Carlo simulation upon a lattice model in which the bosonic part featuring the Landau-Devonshire first-order phase transition and Yukawa coupled to the Dirac fermions. We find that the parameter range for the first-order phase transition becomes larger as the Yukawa coupling increases and the microscopic mechanism of this phenomena is revealed, at a quantitative level, as the interplay between the critical fluctuations and the finite-size effects. Moreover, the scaling behavior at the separation point between the first-order and the continuous phase transitions is found to belong to the chiral tricritical Gross-Neveu universality. Our result demonstrates that the interplay of massless Dirac fermions, critical fluctuations and the finite size effects could trigger a plethora of interesting phenomena and therefore great care is called for when making generalizations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Automatic Differentiation for Second Renormalization of Tensor Networks

Tensor renormalization group (TRG) constitutes an important methodology for accurate simulations of strongly correlated lattice models. Facilitated by the automatic differentiation technique widely used in deep learning, we propose a uniform framework of differentiable TRG ($\partial$TRG) that can be applied to improve various TRG methods, in an automatic fashion. Essentially, $\partial$TRG systematically extends the concept of second renormalization [PRL 103, 160601 (2009)] where the tensor environment is computed recursively in the backward iteration, in the sense that given the forward process of TRG, $\partial$TRG automatically finds the gradient through backpropagation, with which one can deeply "train" the tensor networks. We benchmark $\partial$TRG in solving the square-lattice Ising model, and demonstrate its power by simulating one- and two-dimensional quantum systems at finite temperature. The deep optimization as well as GPU acceleration renders $\partial$TRG manybody simulations with high efficiency and accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Designer Monte Carlo Simulation for Gross-Neveu Transition

In this manuscript, we study quantum criticality of Dirac fermions via large-scale numerical simulations, focusing on the Gross-Neveu-Yukawa(GNY) chiral-Ising quantum critical point with critical bosonic modes coupled with Dirac fermions. We show that finite-size effects at this quantum critical point can be efficiently minimized via model design, which maximizes the ultraviolet cutoff and at the same time places the bare control parameters closer to the nontrivial fixed point to better expose the critical region. Combined with the efficient self-learning quantum Monte Carlo algorithm, which enables non-local update of the bosonic field, we find that moderately-large system size (up to $16\times 16$) is already sufficient to produce robust scaling behavior and critical exponents.The conductance of the Dirac fermions is also calculated and its frequency dependence is found to be consistent with the scaling behavior predicted by the conformal field theory. The methods and model-design principles developed for this study can be generalized to other fermionic QCPs, and thus provide a promising direction for controlled studies of strongly-correlated itinerant systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment from four-flavor lattice QCD

We calculate the contribution to the muon anomalous magnetic moment hadronic vacuum polarization from {the} connected diagrams of up and down quarks, omitting electromagnetism. We employ QCD gauge-field configurations with dynamical $u$, $d$, $s$, and $c$ quarks and the physical pion mass, and analyze five ensembles with lattice spacings ranging from $a \approx 0.06$ to~0.15~fm. The up- and down-quark masses in our simulations have equal masses $m_l$. We obtain, in this world where all pions have the mass of the $π^0$, $10^{10} a_μ^{ll}({\rm conn.}) = 637.8\,(8.8)$, in agreement with independent lattice-QCD calculations. We then combine this value with published lattice-QCD results for the connected contributions from strange, charm, and bottom quarks, and an estimate of the uncertainty due to the fact that our calculation does not include strong-isospin breaking, electromagnetism, or contributions from quark-disconnected diagrams. Our final result for the total $\mathcal{O}(α^2)$ hadronic vacuum polarization to the muon's anomalous magnetic moment is~$10^{10}a_μ^{\rm HVP,LO} = 699(15)_{u,d}(1)_{s,c,b}$, where the errors are from the light-quark and heavy-quark contributions, respectively. Our result agrees with both {\it ab-initio} lattice-QCD calculations and phenomenological determinations from experimental $e^+e^-$-scattering data. It is $1.3σ$ below the "no new physics" value of the hadronic-vacuum-polarization contribution inferred from combining the BNL E821 measurement of $a_μ$ with theoretical calculations of the other contributions.

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

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Lattice study of large $N_c$ QCD

We present a lattice simulation study of large $N_c$ regularities of meson and baryon spectroscopy in $SU(N_c)$ gauge theory with two flavors of dynamical fundamental representation fermions. Systems investigated include $N_c=2$, 3, 4, and 5, over a range of fermion masses parametrized by a squared pseudoscalar to vector meson mass ratio between about 0.2 to 0.7. Good agreement with large $N_c$ scaling is observed in the static potential, in meson masses and decay constants, and in baryon spectroscopy. This is an update of the published version of the paper (Phys. Rev. D94 (2016) 034506).

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Approaching conformality with the Tensor Renormalization Group method

We discuss the Tensor Renormalization Group (TRG) method for the O(2) model with a chemical potential in 1+1 dimensions with emphasis on near gapless/conformal situations. We emphasize the role played by the late Leo Kadanoff in the development of this theoretical framework. We describe the entanglement entropy in the superfluid phase (see arXiv:1507.01471 for details). We present recent progress on optimized truncation methods.

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The renormalization group step scaling function of the 2-flavor SU(3) sextet model

We investigate the discrete $β$ function of the 2-flavor SU(3) sextet model using the finite volume gradient flow scheme. Our results, using clover improved nHYP smeared Wilson fermions, follow the (non-universal) 4-loop $\overline{\textrm{MS}}$ perturbative predictions closely up to $g^2 \approx 5.5$, the strongest coupling reached in our simulation. At strong couplings the results are in tension with a recently published work using the same gradient flow renormalization scheme with staggered fermions. Since these calculations define the discrete $β$ function in the same continuum renormalization scheme, they should lead to the same continuum predictions, irrespective of the lattice fermion action. In order to test systematic effects in our computation we compare two different lattice operators, three different flow definitions, and two volume extrapolations. We find agreement among these different approaches in the continuum limit when the gradient flow parameter $c\gtrsim0.35$. Considering the potential phenomenological impact of this model, it is important to understand the origin of the disagreement between our work and the staggered fermion results.

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exploring the phase structure of 12-flavor $SU(3)$

We are studying the $SU(3)$ gauge theory with 12 staggered fermions, searching for the endpoint of the line of first-order phase transitions in the mass--beta plane. This endpoint plays an important role in our understanding of the phase diagram of this model. Having found this endpoint with high statistics on a small lattice using unimproved staggered fermions, we are working to find it on larger lattices and with improved actions. For an action improved with nHYP-smeared staggered fermions, we discuss the effect of slowly turning off the improvement on the broken shift symmetry phase.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Improving the continuum limit of gradient flow step scaling

We introduce a non-perturbative improvement for the renormalization group step scaling function based on the gradient flow running coupling, which may be applied to any lattice gauge theory of interest. Considering first SU(3) gauge theory with $N_f = 4$ massless staggered fermions, we demonstrate that this improvement can remove $O(a^2)$ lattice artifacts, and thereby increases our control over the continuum extrapolation. Turning to the 12-flavor system, we observe an infrared fixed point in the infinite-volume continuum limit. Applying our proposed improvement reinforces this conclusion by removing all observable $O(a^2)$ effects. For the finite-volume gradient flow renormalization scheme defined by $c = \sqrt{8t} / L = 0.2$, we find the continuum conformal fixed point to be located at $g_\star^2 = 6.2(2)$

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sampling versus Blocking

The idea of blocking in configuration space has played an important role in the development of the RG ideas. However, despite being half a century old and having had a huge intellectual impact, generic numerical methods to perform blocking for lattice models have progressed more slowly than sampling methods. Blocking may be essential to deal with near conformal situations. Typically, blocking methods have smaller statistical errors but larger systematic errors than sampling methods. This situation is evolving with recent developments based on the Tensor RG (TRG) method. We report recent results for spin and gauge lattice models obtained with this new method regarding searches for fixed points, calculations of critical exponents and resolutions of sign problems. An interesting model for comparison is the 2-dimensional O(2) model with a chemical potential which has a sign problem with conventional Monte Carlo but allows sampling with the worm algorithm and blocking with various TRG formulations. We compare the efficiency and accuracy of these two methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectroscopy of SU(4) lattice gauge theory with fermions in the two index anti-symmetric representation

We present a study of spectroscopy of SU(4) lattice gauge theory coupled to two flavors of Dirac fermions in the anti-symmetric two index representation. The fermion representation is real, and the pattern of chiral symmetry breaking is SU(2Nf) -> SO(2Nf) with Nf flavors of Dirac fermions. It is an interesting generalization of QCD, for several reasons: it allows direct exploration of an alternate large Nc expansion, it can be simulated at non-zero chemical potential with no sign problem, and several UV completions of composite Higgs systems are built on it. We present preliminary results on the baryon and meson spectra of the theory and compare them with SU(3) results and with expectations for large Nc scaling.

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparing Tensor Renormalization Group and Monte Carlo calculations for spin and gauge models

We show that the Tensor Renormalization Group (TRG) method can be applied to O(N) spin models, principal chiral models and pure gauge theories (Z2, U(1) and SU(2)) on (hyper) cubic lattices. We explain that contrarily to some common belief, it is very difficult to write compact formulas expressing the blockspinning of lattice models. We show that in contrast to other approaches, the TRG formulation allows us to write exact blocking formulas with numerically controllable truncations. The basic reason is that the TRG blocking separates neatly the degrees of freedom inside the block and which are integrated over, from those kept to communicate with the neighboring blocks. We argue that the TRG is a method that can handle large volumes, which is crucial to approach quasi-conformal systems. The method can also get rid of some sign problems. We discuss recent results regarding the critical properties of the 2D O(2) nonlinear sigma model with complex beta and chemical potential. As some of these results appeared in a recently published paper (PRD 88, 056005) and two recent preprints (arXiv:1309.4963 and arXiv:1309.6623), these proceedings rather emphasize the conceptual aspects of our ongoing effort.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite size scaling of conformal theories in the presence of a near-marginal operator

The slowly evolving gauge coupling of gauge-fermion systems near the conformal window makes numerical investigations of these models challenging. We consider finite size scaling and show that this often used technique leads to inconsistent results if the leading order scaling corrections are neglected. When the corrections are included the results become consistent not only between different operators but even when data obtained at different gauge couplings or with different lattice actions are combined. Our results indicate that the SU(3) 12-fermion system is conformal with mass anomalous dimension $γ_m=0.235(15)$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fisher zeros and RG flows for $SU(3)$ with $N_f$ flavors

We calculate the Fisher zeros for $SU(3)$ gauge theory with different $N_f$ flavors of staggered fermions for various values of the fermion mass. We discuss the finite-size scaling near the end point of the line of discontinuity of $\barψ ψ$ in the beta-mass plane and in the larger beta-lower mass region. We discuss possible interpretations of these results in terms of Wilsonian RG flows and their possible relevance to construct composite Higgs models.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fisher zeros and conformality in lattice models

Fisher zeros are the zeros of the partition function in the complex beta=2N_c/g^2 plane. When they pinch the real axis, finite size scaling allows one to distinguish between first and second order transition and to estimate exponents. On the other hand, a gap signals confinement and the method can be used to explore the boundary of the conformal window. We present recent numerical results for 2D O(N) sigma models, 4D U(1) and SU(2) pure gauge and SU(3) gauge theory with N_f=4 and 12 flavors. We discuss attempts to understand some of these results using analytical methods. We discuss the 2-lattice matching and qualitative aspects of the renormalization group (RG) flows in the Migdal-Kadanoff approximation, in particular how RG flows starting at large beta seem to move around regions where bulk transitions occur. We consider the effects of the boundary conditions on the nonperturbative part of the average energy and on the Fisher zeros for the 1D O(2) model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local gauge symmetry on optical lattices?

The versatile technology of cold atoms confined in optical lattices allows the creation of a vast number of lattice geometries and interactions, providing a promising platform for emulating various lattice models. This opens the possibility of letting nature take care of sign problems and real time evolution in carefully prepared situations. Up to now, experimentalists have succeeded to implement several types of Hubbard models considered by condensed matter theorists. In this proceeding, we discuss the possibility of extending this effort to lattice gauge theory. We report recent efforts to establish the strong coupling equivalence between the Fermi Hubbard model and SU(2) pure gauge theory in 2+1 dimensions by standard determinantal methods developed by Robert Sugar and collaborators. We discuss the possibility of using dipolar molecules and external fields to build models where the equivalence holds beyond the leading order in the strong coupling expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fisher's zeros, complex RG flows and confinement in LGT models

The zeros of the partition function in the complex beta plane (Fisher's zeros) play an important role in our understanding of phase transitions and RG flows. Recently, we argued that they act as gates or separatrices for complex RG flows. Using histogram reweighting to construct the density of states, we calculate the Fisher's zeros for pure gauge SU(2) and U(1) on L^4 lattices. For SU(2), these zeros appear to move almost horizontally when the volume increases. They stay away from the real axis which indicates a confining theory at zero temperature. We discuss the effect of an adjoint term on these results. In contrast, using recent multicanonical simulations for the U(1) model for L up to 8 we find that the zeros pinch the real axis near beta =1.0113. Preliminary results concerning U(1) at larger volumes, SU(3) with 3 light flavors and plans to delimit the boundary of the conformal window are briefly discussed.

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access