Hana Saito

Phase structure of finite temperature QCD in the heavy quark region

We study the quark mass dependence of the finite temperature QCD phase transition in the heavy quark region using an effective potential defined through the probability distribution function of the average plaquette. Performing a simulation of SU(3) pure gauge theory, we first confirm that the distribution function has two peaks indicating that the phase transition is of first order in the heavy quark limit, while the first order transition turns into a crossover as the quark mass decreases from infinity, where the mass dependence of the distribution function is evaluated by the reweighting method combined with the hopping parameter expansion. We determine the endpoint of the first order transition region for Nf = 1, 2, 3 and 2 + 1 cases. The quark mass dependence of the latent heat is also evaluated in the first order transition region.

Thermal evolution of the Schwinger model with matrix product operators

We demonstrate the suitability of tensor network techniques for describing the thermal evolution of lattice gauge theories. As a benchmark case, we have studied the temperature dependence of the chiral condensate in the Schwinger model, using matrix product operators to approximate the thermal equilibrium states for finite system sizes with non-zero lattice spacings. We show how these techniques allow for reliable extrapolations in bond dimension, step width, system size and lattice spacing, and for a systematic estimation and control of all error sources involved in the calculation. The reached values of the lattice spacing are small enough to capture the most challenging region of high temperatures and the final results are consistent with the analytical prediction by Sachs and Wipf over a broad temperature range.

Histograms in heavy-quark QCD at finite temperature and density

We study the phase structure of lattice QCD with heavy quarks at finite temperature and density by a histogram method. We determine the location of the critical point at which the first-order deconfining transition in the heavy-quark limit turns into a crossover at intermediate quark masses through a change of the shape of the histogram under variation of coupling parameters. We estimate the effect of the complex phase factor which causes the sign problem at finite density, and show that, in heavy-quark QCD, the effect is small around the critical point. We determine the critical surface in 2+1 flavor QCD in the heavy-quark region at all values of the chemical potential mu including mu=infty.

Chiral condensate in the Schwinger model with matrix product operators

Tensor network (TN) methods, in particular the Matrix Product States (MPS) ansatz, have proven to be a useful tool in analyzing the properties of lattice gauge theories. They allow for a very good precision, much better than standard Monte Carlo (MC) techniques for the models that have been studied so far, due to the possibility of reaching much smaller lattice spacings. The real reason for the interest in the TN approach, however, is its ability, shown so far in several condensed matter models, to deal with theories which exhibit the notorious sign problem in MC simulations. This makes it prospective for dealing with the non-zero chemical potential in QCD and other lattice gauge theories, as well as with real-time simulations. In this paper, using matrix product operators, we extend our analysis of the Schwinger model at zero temperature to show the feasibility of this approach also at finite temperature. This is an important step on the way to deal with the sign problem of QCD. We analyze in detail the chiral symmetry breaking in the massless and massive cases and show that the method works very well and gives good control over a broad range of temperatures, essentially from zero to infinite temperature.

The temperature dependence of the chiral condensate in the Schwinger model with Matrix Product States

We present our recent results for the tensor network (TN) approach to lattice gauge theories. TN methods provide an efficient approximation for quantum many-body states. We employ TN for one dimensional systems, Matrix Product States, to investigate the 1-flavour Schwinger model. In this study, we compute the chiral condensate at finite temperature. From the continuum extrapolation, we obtain the chiral condensate in the high temperature region consistent with the analytical calculation by Sachs and Wipf.

Numerical study of QCD phase diagram at high temperature and density by a histogram method

We study the QCD phase structure at high temperature and density adopting a histogram method. Because the quark determinant is complex at finite density, the Monte-Carlo method cannot be applied directly. We use a reweighting method and try to solve the problems which arise in the reweighting method, i.e. the sign problem and the overlap problem. We discuss the chemical potential dependence of the probability distribution function in the heavy quark mass region and examine the applicability of the approach in the light quark region.

Finite density QCD phase transition in the heavy quark mass region

We extend our previous study of the QCD phase structure in the heavy quark region to non-zero chemical potentials. To identify the critical point where the first order deconfining transition terminates, we study an effective potential defined by the probability distribution function of the plaquette and the Polyakov loop. The reweighting technique is shown to be powerful in evaluating the effective potential in a wide range of the plaquette and Polyakov loop expectation values. We adopt the cumulant expansion to overcome the sign problem in the calculation of complex phase of the quark determinant. We find that the method provides us with an intuitive and powerful way to study the phase structure. We estimate the location of the critical point at finite chemical potential in the heavy quark region.

Towards overcoming the Monte Carlo sign problem with tensor networks

The study of lattice gauge theories with Monte Carlo simulations is hindered by the infamous sign problem that appears under certain circumstances, in particular at non-zero chemical potential. So far, there is no universal method to overcome this problem. However, recent years brought a new class of non-perturbative Hamiltonian techniques named tensor networks, where the sign problem is absent. In previous work, we have demonstrated that this approach, in particular matrix product states in 1+1 dimensions, can be used to perform precise calculations in a lattice gauge theory, the massless and massive Schwinger model. We have computed the mass spectrum of this theory, its thermal properties and real-time dynamics. In this work, we review these results and we extend our calculations to the case of two flavours and non-zero chemical potential. We are able to reliably reproduce known analytical results for this model, thus demonstrating that tensor networks can tackle the sign problem of a lattice gauge theory at finite density.

An application of the variational analysis to calculate the meson spectral functions

We present a new method to calculate meson spectral functions (SPFs) on the lattice based on a variational method. Because, on a finite volume lattice, the meson SPFs have discrete spectra only, a suitable way to extract such discrete signals is needed. Using a variational method, we can calculate several discrete quantities such as the position and the area of spectral peaks for low-lying states. Moreover data accuracy can be improved by increasing the number of basis functions. In this report, we first confirm our method in the free quark case and show that our method works well. Then, we apply the method to a quenched lattice QCD simulation and calculate the charmonium SPFs for S and P-waves at zero temperature. Our results for the ground state are well consistent with the position and the area of the lowest peaks of charmonium SPFs calculated by the conventional maximum entropy method. For first excited states, the signals may be reliablly extracted with our method because the charmonium mass converges to a value close to the experimental one when the number of basis functions is increased. We also investigate the SPFs for S-wave charmonia at below and above

The multi-flavor Schwinger model with chemical potential - Overcoming the sign problem with Matrix Product States

T_c