Source author record

Hideo Matsufuru

Hideo Matsufuru 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

17works
5topics
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

17 published item(s)

preprint2022arXiv

Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

preprint2021arXiv

General purpose lattice QCD code set Bridge++ 2.0 for high performance computing

Bridge++ is a general-purpose code set for a numerical simulation of lattice QCD aiming at a readable, extensible, and portable code while keeping practically high performance. The previous version of Bridge++ is implemented in double precision with a fixed data layout. To exploit the high arithmetic capability of new processor architecture, we extend the Bridge++ code so that optimized code is available as a new branch, i.e., an alternative to the original code. This paper explains our strategy of implementation and displays application examples to the following architectures and systems: Intel AVX-512 on Xeon Phi Knights Landing, Arm A64FX-SVE on Fujitsu A64FX (Fugaku), NEC SX-Aurora TSUBASA, and GPU cluster with NVIDIA V100.

preprint2016arXiv

Gauge cooling for the singular-drift problem in the complex Langevin method --- an application to finite density QCD

We study full QCD at finite density and low temperature with light quark mass using the complex Langevin method. Since the singular drift problem turns out to be mild on a $4^3 \times 8$ lattice we use, the gauge cooling is performed only to control the unitarity norm in this exploratory study. We report on our preliminary data obtained from the complex Langevin simulation up to certain Langevin time. While the data are still noisy due to lack of statistics, the onset of the baryon number density seems to occur at larger $μ$ than half the pion mass, which is the value for the phase quenched QCD. The validity of our simulation is tested by the recently proposed criterion based on the probability distribution of the drift term.

preprint2014arXiv

Lattice simulation of SU(2) gauge theory with chirally symmetric fermions

We numerically study the SU(2) gauge theory with two dynamical flavors of the domain-wall fermions in fundamental representation. The meson spectra and the residual mass are measured on three lattice volumes and at two values of gauge coupling so as to investigate the finite volume effect. On generated configurations, eigenvalues of the overlap fermion operator are determined and compared to the random matrix theory. To quantify the effect of violation of the exact chiral symmetry, we measure the correlation between the eigenvectors of the domain-wall and the overlap operators.

preprint2013arXiv

Finite temperature study of the axial U(1) symmetry on the lattice with overlap fermion formulation

We examine the axial U(1) symmetry near and above the finite temperature phase transition in two-flavor QCD using lattice QCD simulations. Although the axial U(1) symmetry is always violated by quantization, (i.e.) the chiral anomaly, the correlation functions may manifest effective restoration of the symmetry in the high temperature phase. We explicitly study this possibility by calculating the meson correlators as well as the Dirac operator spectral density near the critical point. Our numerical simulations are performed on a $16^3\times 8$ lattice with two flavors of dynamical quarks represented by the overlap fermion formalism. Chiral symmetry and its violation due to the axial anomaly is manifestly realized with this formulation, which is a prerequisite for the study of the effective restoration of the axial U(1) symmetry. In order to avoid discontinuity in the gauge configuration space, which occurs for the exactly chiral lattice fermions, the simulation is confined in a fixed topological sector. It induces finite volume effect, which is well described by a formula based on the Fourier transform from the $θ$-vacua. We confirm this formula at finite temperature by calculating the topological susceptibility in the quenched theory. Our two flavor simulations show degeneracy of the meson correlators and a gap in the Dirac operator spectral density, which implies that the axial U(1) symmetry is effectively restored in the chirally symmetric phase.

preprint2012arXiv

Topological susceptibility and axial symmetry at finite temperature

We consider the simulation of finite temperature QCD with two flavors of dynamical overlap fermions in order to study the suppression of the axial U(1) symmetry breaking at the chiral phase transition point. As a preliminary study, pure gauge simulations are performed to investigate how fixing the topology affects physical quantities like the topological susceptibility, $χ_t$, at finite temperature, showing that it is possible to reconstruct known results from the fixed topology sector. First results on the degeneracy of meson correlators in the high temperature QGP sector are shown.

preprint2011arXiv

The Infrared behavior of SU(3) Nf=12 gauge theory -about the existence of conformal fixed point-

Incorporated with twisted boundary condition, Polyakov loop correlators can give a definition of the renormalized coupling. We employ this scheme for the step scaling method (with step size s = 2) in the search of conformal fixed point of SU(3) gauge theory with 12 massless flavors. Staggered fermion and plaquette gauge action are used in the lattice simulation with six different lattice sizes, L/a = 20, 16, 12, 10, 8 and 6. For the largest lattice size, L/a = 20, we used a large number of Graphics Processing Units (GPUs) and accumulated 3,000,000 trajectories in total. We found that the step scaling function sigma (u) is consistent with u in the low-energy region. This means the existence of conformal fixed point. Some details of our analysis and simulations will also be presented.

preprint2010arXiv

A new scheme for the running coupling constant in gauge theories using Wilson loops

We propose a new renormalization scheme of the running coupling constant in general gauge theories using the Wilson loops. The renormalized coupling constant is obtained from the Creutz ratio in lattice simulations and the corresponding perturbative coefficient at the leading order. The latter can be calculated by adopting the zeta-function resummation techniques. We perform a benchmark test of our scheme in quenched QCD with the plaquette gauge action. The running of the coupling constant is determined by applying the step-scaling procedure. Using several methods to improve the statistical accuracy, we show that the running coupling constant can be determined in a wide range of energy scales with relatively small number of gauge configurations.

preprint2010arXiv

Search for the IR fixed point in the Twisted Polyakov Loop scheme (II)

We measure the renormalized coupling in the Twisted Polyakov loop scheme for SU(3) gauge theory coupled with $N_f=12$ fundamental fermions. To find the infrared fixed point of this theory, we focus on the step scaling function for the renormalized coupling. We take the continuum limit using the linear function of $(a/L)^2$ and a constant fit function. We find that there is a sizeable systematic error due to the choice of the continuum extrapolation function in the low energy region. We will give some directions to reduce the systematic errors.

preprint2010arXiv

Study of the scaling properties in SU(2) gauge theory with eight flavors

We present our preliminary study of the SU(2) gauge theory with 8 flavors of fermions in fundamental representation. This theory could be a candidate of the gauge theory with conformal fixed point. By using Wilson/Polyakov loop in a finite volume with twisted boundary conditions, we study the renormalization group flow of the gauge coupling constant. Our calculation gives consistent result with the perturbative prediction of the running coupling in the weak coupling region. We investigate a possible signal for conformal behavior in the strong coupling region.

preprint2009arXiv

Search for the IR fixed point in the twisted Polyakov loop scheme

We present a non-perturbative study of the running coupling constant in the Twisted Polyakov Loop (TPL) scheme. We investigate how the systematic and statistical errors can be controlled {\it via} a feasibility study in SU(3) pure Yang-Mills theory. We show that our method reproduces the perturbative determination of the running coupling in the UV. In addition, our numerical result agrees with the theoretical prediction of this coupling constant in the IR. We also present our preliminary results for $N_f=12$ QCD, where an IR fixed point may be present.

preprint2008arXiv

A New Method of Calculating the Running Coupling Constant

We propose a new method to compute the running coupling constant of gauge theories on the lattice. We first give the definition of the running coupling in the new scheme using the Wilson loops in a finite volume, and explain how the running of the coupling constant is extracted from the measurement of the volume dependence. The perturbative calculation of the renormalization constant to define the scheme is also given at the leading order. As a benchmark test of the new scheme we apply the method the case of the quenched QCD. We show the preliminary result from our numerical simulations which are carried out with plaquette gauge action for various lattice sizes and bare lattice couplings. With techniques to improve the statistical accuracy, we show that we can determine the non-perturbative running of the coupling constant in a wide range of the energy scale with relatively small number of gauge configurations in our scheme. We compare our lattice data of the running coupling constant with perturbative renormalization group evolution at one- and two-loop order, and confirm the consistency between them at high energy.

preprint2004arXiv

Contribution of disconnected diagrams to the hyperfine splitting of charmonium

We discuss the contribution of OZI-suppressed diagrams to the hyperfine splitting of charmonium in lattice QCD. We study valence quark mass regions from strange to charm quark masses. No contribution of the disconnected diagram is seen in the vector meson channel. In the pseudo-scalar channel and for valence quark masses around the strange quark, the disconnected contribution induces a considerable increase of the meson mass. This contribution quickly decreases as the quark mass increases. For charmonium the effect is very small although a decrease of the pseudoscalar mass induced by the disconnected contribution cannot be ruled out.

preprint2003arXiv

Charmonium at finite temperature in quenched lattice QCD

We study charmonium correlators in pseudoscalar and vector channels at finite temperature using lattice QCD simulation in the quenched approximation. Anisotropic lattices are used in order to have sufficient numbers of degrees of freedom in the Euclidean temporal direction. We focus on the low energy structure of the spectral function, corresponding to the ground state in the hadron phase, by applying the smearing technique to enhance the contribution to the correlator from this region. We employ two analysis procedures: the maximum entropy method (MEM) for the extraction of the spectral function without assuming a specific form, to estimate the shape of the spectral function, and the standard $χ^2$ fit analysis using typical forms in accordance with the result of MEM, for a more quantitative evaluation. To verify the applicability of the procedures, we first analyze the smeared correlators as well as the point correlators at zero temperature. We find that by shortening the $t$-interval used for the analysis (a situation inevitable at $T>0$) the reliability of MEM for point correlators is lost, while it subsists for smeared correlators. Then the smeared correlators at $T\simeq 0.9 T_c$ and $1.1 T_c$ are analyzed. At $T\simeq 0.9 T_c$, the spectral function exhibits a strong peak, well approximated by a delta function corresponding to the ground state with almost the same mass as at T=0. At $T\simeq 1.1 T_c$, we find that the strong peak structure still persists at almost the same place as below $T_c$, but with a finite width of a few hundred MeV. This result indicates that the correlators possess a nontrivial structure even in the deconfined phase.

preprint2001arXiv

O(a)-improved quark action on anisotropic lattices and perturbative renormalization of heavy-light currents

We investigate the Symanzik improvement of the Wilson quark action on anisotropic lattices. Taking first a general action with nearest-neighbor and clover interactions, we study the mass dependence of the ratio of the hopping parameters, the clover coefficients, and an improvement coefficient for heavy-light vector and axial vector currents. We show how tree-level improvement can be achieved. For a particular choice of the spatial Wilson coupling, the results simplify, and $O(m_0a_τ)$ improvement is possible. (Here $m_0$ is the bare quark mass and $a_τ$ the temporal lattice spacing.) With this choice we calculate the renormalization factors of heavy-light bilinear operators at one-loop order of perturbation theory employing the standard plaquette gauge action.