Hideo Matsufuru
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Hideo Matsufuru has a total of 111 co-authors. Here are the top co-authors by number of shared publications:
- Tetsuya Onogi
Yukawa Institute for Theoretical Physics Kyoto University 606-8502 Kyoto Japan
10 shared publications - Asato Tsuchiya
Graduate School of Science and Technology Shizuoka University 836 Ohya Suruga-ku 422-8529 Shizuoka Japan
9 shared publications - Jun Nishimura
High Energy Accelerator Research Organization
9 shared publications - Yuta Ito
National Institute of Technology Tokuyama College Gakuendai 745-8585 Shunan Yamaguchi Japan
9 shared publications - Issaku Kanamori
RIKEN Center for Computational Science
7 shared publications - Shoji Hashimoto
KEK Theory Center, High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801, Japan
6 shared publications
Latest Publications
- Top Cited Papers
- 11 May 2026
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 136 (19) , 191002
Bifurcated Impact of Neutrino Fast Flavor Conversion on Core-Collapse Supernovae Informed by Multiangle Neutrino Radiation Hydrodynamics
In this Letter, we present a compelling and robust argument for the roles of neutrino fast flavor conversion (FFC) in the explosion mechanism of core-collapse supernova (CCSN), combining the multiangle FFC subgrid model rooted in quantum kinetic theory with the multidimensional four-species Boltzmann neutrino radiation hydrodynamics. Employing various progenitor masses and the nuclear equations of states, we find that the effect of FFC on CCSN explosion is bifurcated depending on the progenitors. For the lowest-mass progenitor, FFC facilitates the shock revival and enhances the explosion energy, whereas for higher-mass progenitors its impact is inhibitory. We identify the mass accretion rate as the key determinant governing this bifurcation. When the mass accretion rate is low (high), the contribution of FFC to neutrino heating becomes positive (negative) because the heating efficiency enhancement via FFC-driven spectral hardening of electron-type neutrinos dominates over (is outweighed by) the concurrent reduction in neutrino luminosity. Our results further highlight the limitations of approximate neutrino transport and demonstrate that a multiangle treatment is essential for accurately capturing FFC effects; otherwise, FFCs are missed and even generated spuriously.
- 11 May 2026
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 136 (19) , 191002
Top Cited PapersBifurcated Impact of Neutrino Fast Flavor Conversion on Core-Collapse Supernovae Informed by Multiangle Neutrino Radiation Hydrodynamics
In this Letter, we present a compelling and robust argument for the roles of neutrino fast flavor conversion (FFC) in the explosion mechanism of core-collapse supernova (CCSN), combining the multiangle FFC subgrid model rooted in quantum kinetic theory with the multidimensional four-species Boltzmann neutrino radiation hydrodynamics. Employing various progenitor masses and the nuclear equations of states, we find that the effect of FFC on CCSN explosion is bifurcated depending on the progenitors. For the lowest-mass progenitor, FFC facilitates the shock revival and enhances the explosion energy, whereas for higher-mass progenitors its impact is inhibitory. We identify the mass accretion rate as the key determinant governing this bifurcation. When the mass accretion rate is low (high), the contribution of FFC to neutrino heating becomes positive (negative) because the heating efficiency enhancement via FFC-driven spectral hardening of electron-type neutrinos dominates over (is outweighed by) the concurrent reduction in neutrino luminosity. Our results further highlight the limitations of approximate neutrino transport and demonstrate that a multiangle treatment is essential for accurately capturing FFC effects; otherwise, FFCs are missed and even generated spuriously.
- Open Access
Accelerating iterative linear equation solver using modified domain-wall fermion matrix in lattice QCD simulations
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations of lattice QCD is typically solving a linear equation for the quark matrix. In particular, a discretized quark formulation called the domain-wall fermion operator requires a high numerical cost, while retaining the lattice version of the chiral symmetry to good precision. The domain-wall operator is defined on a five-dimensional (5D) space extending the four-dimensional (4D) spacetime with an extra fifth coordinate. After solving the linear equation in 5D space, the result vector is projected onto the original 4D space. There is a variant of the domain-wall operator that improves the convergence of the 5D linear equation while unchanging the 4D solution vector. In this paper, we examine how this variant of the domain-wall operator accelerates the iterative linear equation solver in practical setups. We also measure the eigenvalues of the operator and compare the condition number with the convergence of the solver. We use a generic lattice QCD code set Bridge++ that is planned to be released including the improved form of the domain-wall operator examined in this work with code for the GPU.
Open AccessAccelerating iterative linear equation solver using modified domain-wall fermion matrix in lattice QCD simulations
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations of lattice QCD is typically solving a linear equation for the quark matrix. In particular, a discretized quark formulation called the domain-wall fermion operator requires a high numerical cost, while retaining the lattice version of the chiral symmetry to good precision. The domain-wall operator is defined on a five-dimensional (5D) space extending the four-dimensional (4D) spacetime with an extra fifth coordinate. After solving the linear equation in 5D space, the result vector is projected onto the original 4D space. There is a variant of the domain-wall operator that improves the convergence of the 5D linear equation while unchanging the 4D solution vector. In this paper, we examine how this variant of the domain-wall operator accelerates the iterative linear equation solver in practical setups. We also measure the eigenvalues of the operator and compare the condition number with the convergence of the solver. We use a generic lattice QCD code set Bridge++ that is planned to be released including the improved form of the domain-wall operator examined in this work with code for the GPU.
- Open Access
- 25 January 2026
- Published by Association for Computing Machinery (ACM)
- p. 133-144
Solving large-scale eigen problem in quantum few-body system on massive parallel computer
- 25 January 2026
- Published by Association for Computing Machinery (ACM)
- p. 133-144
Open AccessSolving large-scale eigen problem in quantum few-body system on massive parallel computer
- Open Access
- 25 January 2026
- Published by Association for Computing Machinery (ACM)
- p. 145-155
Mixed precision solvers with half-precision floating point numbers for Lattice QCD on A64FX processor
- 25 January 2026
- Published by Association for Computing Machinery (ACM)
- p. 145-155
Open AccessMixed precision solvers with half-precision floating point numbers for Lattice QCD on A64FX processor
- Open Access
- 14 October 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (10) , 1-24
On the validity of the complex Langevin method near the deconfining phase transition in QCD at finite density
A bstract: In our previous paper [1], we found that the complex Langevin (CL) method works for QCD at finite density on the 16 3 × 32 lattice in the low-temperature high-density regime within the range μ / T = 1.6 – 9.6 with μ and T being the quark chemical potential and the temperature, which enabled us to see a clear trend towards the formation of the Fermi sphere. Here we investigate the validity of the CL method on the 24 3 × 12 lattice in the deconfined phase near the deconfinement phase transition. As before, we use four-flavor staggered fermions and judge the validity using the criterion based on the probability distribution of the drift term. The spatial extent is L = (1.3 − 2.7 fm) > $$ {\Lambda}_{\textrm{LQCD}}^{-1} $$ ∼ 1 fm, in contrast to our previous study with L < $$ {\Lambda}_{\textrm{LQCD}}^{-1} $$ . We find that the CL method works in a broad region up to μ / T = 4.8, while it starts to fail as we approach the phase boundary due to the singular drift problem, which can be understood qualitatively by extending the Banks-Casher relation to the case at finite density.- 14 October 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (10) , 1-24
Open AccessOn the validity of the complex Langevin method near the deconfining phase transition in QCD at finite density
A bstract: In our previous paper [1], we found that the complex Langevin (CL) method works for QCD at finite density on the 16 3 × 32 lattice in the low-temperature high-density regime within the range μ / T = 1.6 – 9.6 with μ and T being the quark chemical potential and the temperature, which enabled us to see a clear trend towards the formation of the Fermi sphere. Here we investigate the validity of the CL method on the 24 3 × 12 lattice in the deconfined phase near the deconfinement phase transition. As before, we use four-flavor staggered fermions and judge the validity using the criterion based on the probability distribution of the drift term. The spatial extent is L = (1.3 − 2.7 fm) > $$ {\Lambda}_{\textrm{LQCD}}^{-1} $$ ∼ 1 fm, in contrast to our previous study with L < $$ {\Lambda}_{\textrm{LQCD}}^{-1} $$ . We find that the CL method works in a broad region up to μ / T = 4.8, while it starts to fail as we approach the phase boundary due to the singular drift problem, which can be understood qualitatively by extending the Banks-Casher relation to the case at finite density. - 28 June 2025
- Published by Springer Nature
- p. 54-69
Accelerating Iterative Linear Equation Solver Using Modified Domain-Wall Fermion Matrix in Lattice QCD Simulations
- 28 June 2025
- Published by Springer Nature
- p. 54-69
Accelerating Iterative Linear Equation Solver Using Modified Domain-Wall Fermion Matrix in Lattice QCD Simulations
- Open Access
- 24 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 414
Spectrum of Preconditioned Moebius Domain-wall Operators
The convergence property of iterative solvers strongly depends on the spectrum of the Dirac operator. For most standard algorithms to work, the real part of the spectrum should be positive. The domain-wall operator does not satisfy this condition, which is one of the reasons for the difficulty in applying the multi-grid algorithms. In this presentation, we examine several preconditioning operators for the Moebius domain-wall operator and investigate their spectra.
- 24 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 414
Open AccessSpectrum of Preconditioned Moebius Domain-wall Operators
The convergence property of iterative solvers strongly depends on the spectrum of the Dirac operator. For most standard algorithms to work, the real part of the spectrum should be positive. The domain-wall operator does not satisfy this condition, which is one of the reasons for the difficulty in applying the multi-grid algorithms. In this presentation, we examine several preconditioning operators for the Moebius domain-wall operator and investigate their spectra.
- Open Access
- 19 February 2025
- Published by Association for Computing Machinery (ACM)
- p. 80-89
Lattice QCD code on GPUs: Implementation and performancecomparison with OpenACC and CUDA
- 19 February 2025
- Published by Association for Computing Machinery (ACM)
- p. 80-89
Open AccessLattice QCD code on GPUs: Implementation and performancecomparison with OpenACC and CUDA
- Open Access
- 13 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 411
International Lattice Data Grid 2.0: Status and Progress
In this proceeding contribution we discuss the status and progress towards a modernised and extended International Lattice Data Grid (ILDG), which has seen major developments, updates, and upgrades over the last year. In particular, metadata and file schemata have been extended. Moreover, the registration and authentication services have been modernised, and the file and metadata catalogues re-implemented.
- 13 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 411
Open AccessInternational Lattice Data Grid 2.0: Status and Progress
In this proceeding contribution we discuss the status and progress towards a modernised and extended International Lattice Data Grid (ILDG), which has seen major developments, updates, and upgrades over the last year. In particular, metadata and file schemata have been extended. Moreover, the registration and authentication services have been modernised, and the file and metadata catalogues re-implemented.