Ying-Jer Kao
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
Research Interests
Ying-Jer Kao has a total of 205 co-authors. Here are the top co-authors by number of shared publications:
- Pochung Chen
National Tsing Hua University
15 shared publications K. Levin
James Franck Institute and Department of Physics, University of Chicago, Chicago, Illinois 60637, USA
8 shared publications- Yu-Ping Lin
University of California
7 shared publications - C. -J. David Lin
National yang Ming Chiao Tung University
6 shared publications - Qijin Chen
University of Science and Technology of China
6 shared publications Yun-Da Hsieh
Department of Physics, National Taiwan University, Taipei, 10607, Taiwan
6 shared publications
Latest Publications
- Open Access
Improved Ising Meson Spectroscopy Simulation on a Noisy Digital Quantum Device
The transverse-field Ising model serves as a paradigm for studying confinement and excitation spectra, particularly the emergence of $E_8$ symmetry near criticality. However, experimentally resolving the Ising meson spectroscopy required to verify these symmetries is challenging on near-term quantum hardware due to the depth of circuits required for real-time evolution. Here, we demonstrate improved spectroscopy of confined excitations using two distinct error-resilient circuit construction techniques on the IBM Torino device: first-order Trotter decomposition utilizing native fractional gates, and a tensor-network-based circuit compression via Riemannian optimization. By analyzing the Fourier spectrum of error-mitigated time-series data, we successfully identify key signatures of $E_8$ symmetry despite hardware noise. These results validate the viability of both circuit compression and hardware-efficient compilation for probing complex topological phenomena on NISQ devices.Open AccessImproved Ising Meson Spectroscopy Simulation on a Noisy Digital Quantum Device
The transverse-field Ising model serves as a paradigm for studying confinement and excitation spectra, particularly the emergence of $E_8$ symmetry near criticality. However, experimentally resolving the Ising meson spectroscopy required to verify these symmetries is challenging on near-term quantum hardware due to the depth of circuits required for real-time evolution. Here, we demonstrate improved spectroscopy of confined excitations using two distinct error-resilient circuit construction techniques on the IBM Torino device: first-order Trotter decomposition utilizing native fractional gates, and a tensor-network-based circuit compression via Riemannian optimization. By analyzing the Fourier spectrum of error-mitigated time-series data, we successfully identify key signatures of $E_8$ symmetry despite hardware noise. These results validate the viability of both circuit compression and hardware-efficient compilation for probing complex topological phenomena on NISQ devices. - Open Access
- 27 October 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (10) , 1-36
Lattice investigation of custodial two-Higgs-doublet model at weak quartic couplings
A bstract: The SU(2)–gauged custodial two-Higgs-doublet model, which shares the same global-symmetry properties with the standard model, is studied non-perturbatively on the lattice. The additional Higgs doublet enlarges the scalar spectrum and opens the possibility for spontaneous breaking of the global symmetry. In this work we start by showing the occurrence of spontaneous breaking of the custodial symmetry in a region of the parameter space of the model. Following this, both the spectrum and the running of the gauge coupling of are examined at weak quartic couplings in the presence of the custodial symmetry. The calculations are performed with energy cutoffs ranging from 300 to 600 GeV on a line of constant standard model physics, obtained by tuning bare couplings to fix the ratio between the masses of the Higgs and the W bosons, as well as the value of the renormalized gauge coupling at the scale of the W boson mass. The realizable masses for the additional scalar states are explored. For the choice of bare quartic couplings in this work, the estimated lower bound of these masses is found to be well below the W boson mass, and independent of the cutoff. We also study the finite temperature electroweak transition along this line of constant standard model physics, revealing properties of a smooth crossover behavior.
- 27 October 2025
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2025 (10) , 1-36
Open AccessLattice investigation of custodial two-Higgs-doublet model at weak quartic couplings
A bstract: The SU(2)–gauged custodial two-Higgs-doublet model, which shares the same global-symmetry properties with the standard model, is studied non-perturbatively on the lattice. The additional Higgs doublet enlarges the scalar spectrum and opens the possibility for spontaneous breaking of the global symmetry. In this work we start by showing the occurrence of spontaneous breaking of the custodial symmetry in a region of the parameter space of the model. Following this, both the spectrum and the running of the gauge coupling of are examined at weak quartic couplings in the presence of the custodial symmetry. The calculations are performed with energy cutoffs ranging from 300 to 600 GeV on a line of constant standard model physics, obtained by tuning bare couplings to fix the ratio between the masses of the Higgs and the W bosons, as well as the value of the renormalized gauge coupling at the scale of the W boson mass. The realizable masses for the additional scalar states are explored. For the choice of bare quartic couplings in this work, the estimated lower bound of these masses is found to be well below the W boson mass, and independent of the cutoff. We also study the finite temperature electroweak transition along this line of constant standard model physics, revealing properties of a smooth crossover behavior.
- Open Access
- 27 October 2025
- Published by IOP Publishing in Machine Learning: Science and Technology
- Vol. 6 (4) , 045017
Learning phases with quantum Monte Carlo simulation cell
We propose the use of the ‘spin-opstring’, derived from Stochastic Series Expansion quantum Monte Carlo (QMC) simulations as machine learning (ML) input data. It offers a compact, memory-efficient representation of QMC simulation cells, combining the initial state with an operator string that encodes the state’s evolution through imaginary time. Using supervised ML, we demonstrate the input’s effectiveness in capturing both conventional and topological phase transitions, and in a regression task to predict non-local observables. We also demonstrate the capability of spin-opstring data in transfer learning by training models on one quantum system and successfully predicting on another, as well as showing that models trained on smaller system sizes generalize well to larger ones. Importantly, we illustrate a clear advantage of spin-opstring over conventional spin configurations in the accurate prediction of a quantum phase transition. Finally, we show how the inherent structure of spin-opstring provides an elegant framework for the interpretability of ML predictions. Using two state-of-the-art interpretability techniques, Layer-wise Relevance Propagation and SHapley Additive exPlanations, we show that the ML models learn and rely on physically meaningful features from the input data. Together, these findings establish the spin-opstring as a broadly-applicable and interpretable input format for ML in quantum many-body physics.
- 27 October 2025
- Published by IOP Publishing in Machine Learning: Science and Technology
- Vol. 6 (4) , 045017
Open AccessLearning phases with quantum Monte Carlo simulation cell
We propose the use of the ‘spin-opstring’, derived from Stochastic Series Expansion quantum Monte Carlo (QMC) simulations as machine learning (ML) input data. It offers a compact, memory-efficient representation of QMC simulation cells, combining the initial state with an operator string that encodes the state’s evolution through imaginary time. Using supervised ML, we demonstrate the input’s effectiveness in capturing both conventional and topological phase transitions, and in a regression task to predict non-local observables. We also demonstrate the capability of spin-opstring data in transfer learning by training models on one quantum system and successfully predicting on another, as well as showing that models trained on smaller system sizes generalize well to larger ones. Importantly, we illustrate a clear advantage of spin-opstring over conventional spin configurations in the accurate prediction of a quantum phase transition. Finally, we show how the inherent structure of spin-opstring provides an elegant framework for the interpretability of ML predictions. Using two state-of-the-art interpretability techniques, Layer-wise Relevance Propagation and SHapley Additive exPlanations, we show that the ML models learn and rely on physically meaningful features from the input data. Together, these findings establish the spin-opstring as a broadly-applicable and interpretable input format for ML in quantum many-body physics.
- Open Access
- 7 June 2025
- Published by Springer Nature in npj Quantum Materials
- Vol. 10 (1) , 1-8
Direct visualization of a disorder driven electronic smectic phase in nonsymmorphic square-net semimetal GdSbTe
- 7 June 2025
- Published by Springer Nature in npj Quantum Materials
- Vol. 10 (1) , 1-8
Open AccessDirect visualization of a disorder driven electronic smectic phase in nonsymmorphic square-net semimetal GdSbTe
- Open Access
- 28 May 2025
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 7 (2)
Dynamical quantum phase transition and thermal equilibrium in the lattice Thirring model
Using tensor network methods, we simulate the real-time evolution of the lattice Thirring model quenched out of equilibrium in both the critical and massive phases and study the appearance of dynamical quantum phase transitions, as nonanalyticities in the Loschmidt rate. Although the presence of a dynamical quantum phase transition in the model does not correspond to quenches across the critical line of the equilibrium phase diagram at zero temperature, we identify a threshold in the energy density of the initial state, necessary for a dynamical quantum phase transition to be present. Moreover, in the case of the gapped quench Hamiltonian, we unveil a connection of this threshold to a transition between different regions in the finite-temperature phase diagram. Published by the American Physical Society 2025
- 28 May 2025
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 7 (2)
Open AccessDynamical quantum phase transition and thermal equilibrium in the lattice Thirring model
Using tensor network methods, we simulate the real-time evolution of the lattice Thirring model quenched out of equilibrium in both the critical and massive phases and study the appearance of dynamical quantum phase transitions, as nonanalyticities in the Loschmidt rate. Although the presence of a dynamical quantum phase transition in the model does not correspond to quenches across the critical line of the equilibrium phase diagram at zero temperature, we identify a threshold in the energy density of the initial state, necessary for a dynamical quantum phase transition to be present. Moreover, in the case of the gapped quench Hamiltonian, we unveil a connection of this threshold to a transition between different regions in the finite-temperature phase diagram. Published by the American Physical Society 2025
- Open Access
- 13 March 2025
- Published by Stichting SciPost in SciPost Physics Codebases
The Cytnx library for tensor networks
We introduce a tensor network library designed for classical and quantum physics simulations called Cytnx (pronounced as sci-tens). This library provides an almost identical interface and syntax for both C++ and Python, allowing users to effortlessly switch between the two languages. Aiming at a quick learning process for new users of tensor network algorithms, the interfaces resemble the popular Python scientific libraries like NumPy, Scipy, and PyTorch. Not only multiple global Abelian symmetries can be easily defined and implemented, Cytnx also provides a new tool called Network that allows users to store large tensor networks and perform tensor network contractions in an optimal order automatically. With the integration of cuQuantum, tensor calculations can also be executed efficiently on GPUs. We present benchmark results for tensor operations on both devices, CPU and GPU. We also discuss features and higher-level interfaces to be added in the future.
- 13 March 2025
- Published by Stichting SciPost in SciPost Physics Codebases
Open AccessThe Cytnx library for tensor networks
We introduce a tensor network library designed for classical and quantum physics simulations called Cytnx (pronounced as sci-tens). This library provides an almost identical interface and syntax for both C++ and Python, allowing users to effortlessly switch between the two languages. Aiming at a quick learning process for new users of tensor network algorithms, the interfaces resemble the popular Python scientific libraries like NumPy, Scipy, and PyTorch. Not only multiple global Abelian symmetries can be easily defined and implemented, Cytnx also provides a new tool called Network that allows users to store large tensor networks and perform tensor network contractions in an optimal order automatically. With the integration of cuQuantum, tensor calculations can also be executed efficiently on GPUs. We present benchmark results for tensor operations on both devices, CPU and GPU. We also discuss features and higher-level interfaces to be added in the future.
- Open Access
- 19 November 2024
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 6 (4)
Lattice study of SU(2) gauge theory coupled to four adjoint Higgs fields
Gauge theories with matter fields in various representations play an important role in different branches of physics. Recently, it was proposed that several aspects of the pseudogap phase of cuprate superconductors near optimal doping may be explained by an emergent SU ( 2 ) gauge symmetry. Around the transition with positive hole doping, one can construct a ( 2 + 1 ) − dimensional SU(2) gauge theory coupled to four adjoint scalar fields which gives rise to a rich phase diagram with a myriad of phases having different broken symmetries. We study the phase diagram of this model on the Euclidean lattice using the hybrid Monte Carlo algorithm. We find the existence of multiple broken phases as predicted by previous mean-field studies. Depending on the quartic couplings, the SU(2) gauge symmetry is broken down either to U(1) or Z 2 in the perturbative description of the model. We further study the confinement-deconfinement transition in this theory, and find that both the broken phases are deconfining in the range of volumes that we studied. However, there exists a marked difference in the behavior of the Polyakov loop between the two phases. Published by the American Physical Society 2024
- 19 November 2024
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 6 (4)
Open AccessLattice study of SU(2) gauge theory coupled to four adjoint Higgs fields
Gauge theories with matter fields in various representations play an important role in different branches of physics. Recently, it was proposed that several aspects of the pseudogap phase of cuprate superconductors near optimal doping may be explained by an emergent SU ( 2 ) gauge symmetry. Around the transition with positive hole doping, one can construct a ( 2 + 1 ) − dimensional SU(2) gauge theory coupled to four adjoint scalar fields which gives rise to a rich phase diagram with a myriad of phases having different broken symmetries. We study the phase diagram of this model on the Euclidean lattice using the hybrid Monte Carlo algorithm. We find the existence of multiple broken phases as predicted by previous mean-field studies. Depending on the quartic couplings, the SU(2) gauge symmetry is broken down either to U(1) or Z 2 in the perturbative description of the model. We further study the confinement-deconfinement transition in this theory, and find that both the broken phases are deconfining in the range of volumes that we studied. However, there exists a marked difference in the behavior of the Polyakov loop between the two phases. Published by the American Physical Society 2024
- Open Access
- 15 January 2024
- Published by Springer Nature in Communications Physics
- Vol. 7 (1) , 1-8
Chern dartboard insulator: sub-Brillouin zone topology and skyrmion multipoles
Topology plays a crucial role in many physical systems, leading to interesting states at the surface. A paradigmatic example is the Chern number defined in the Brillouin zone that leads to robust gapless edge states. Here we introduce the reduced Chern number, defined in the subregions of Brillouin zone, and construct a family of Chern dartboard insulators with quantized reduced Chern numbers but with trivial bulk topology. Chern dartboard insulators are protected by the mirror symmetries and exhibit distinct pseudospin textures, including (anti)skyrmions, inside the sub-Brillouin zone. These Chern dartboard insulators host exotic gapless edge states, such as Möbius fermions and midgap corner states, and can be realized in the photonic crystals. Our work opens up new possibilities for exploring sub-Brillouin zone topology and nontrivial surface responses in topological systems.
- 15 January 2024
- Published by Springer Nature in Communications Physics
- Vol. 7 (1) , 1-8
Open AccessChern dartboard insulator: sub-Brillouin zone topology and skyrmion multipoles
Topology plays a crucial role in many physical systems, leading to interesting states at the surface. A paradigmatic example is the Chern number defined in the Brillouin zone that leads to robust gapless edge states. Here we introduce the reduced Chern number, defined in the subregions of Brillouin zone, and construct a family of Chern dartboard insulators with quantized reduced Chern numbers but with trivial bulk topology. Chern dartboard insulators are protected by the mirror symmetries and exhibit distinct pseudospin textures, including (anti)skyrmions, inside the sub-Brillouin zone. These Chern dartboard insulators host exotic gapless edge states, such as Möbius fermions and midgap corner states, and can be realized in the photonic crystals. Our work opens up new possibilities for exploring sub-Brillouin zone topology and nontrivial surface responses in topological systems.
- 11 May 2023
- Published by American Physical Society (APS) in Physical Review B
- Vol. 107 (20) , 205123
Tensor network based finite-size scaling for two-dimensional Ising model
We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density, which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension , the absolute errors of the critical temperature and the critical exponent , are at the order of , , respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cutoff in bond dimension and elucidate its physical meaning in this context.- 11 May 2023
- Published by American Physical Society (APS) in Physical Review B
- Vol. 107 (20) , 205123
Tensor network based finite-size scaling for two-dimensional Ising model
We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density, which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension , the absolute errors of the critical temperature and the critical exponent , are at the order of , , respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cutoff in bond dimension and elucidate its physical meaning in this context.- Open Access
Tensor Network Based Finite-Size Scaling for Two-Dimensional Classical Models
We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension $D=70$, the absolute errors of the critical temperature $T_c$ and the critical exponent $ν$, $β$ are at the order of $10^{-7}, 10^{-5}$, $10^{-4}$ respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cut-off in bond dimension and elucidate its physical meaning in this context.Open AccessTensor Network Based Finite-Size Scaling for Two-Dimensional Classical Models
We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension $D=70$, the absolute errors of the critical temperature $T_c$ and the critical exponent $ν$, $β$ are at the order of $10^{-7}, 10^{-5}$, $10^{-4}$ respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cut-off in bond dimension and elucidate its physical meaning in this context.