I work on non-perturbative analyses of quantum many-body systems, with a current emphasis on quantum computing as a platform to explore nonequilibrium dynamics. My research combines ideas from quantum-circuit dynamics, digital quantum simulation, and strongly correlated condensed matter, aiming to identify robust dynamical phenomena and connect them to microscopic models and experimentally relevant probes.
Figure caption: a, Two-dimensional lattice geometry of IBM’s 133-qubit device (ibm_torino), which serves as the platform for the two-dimensional quantum many-body system. b, Quantum circuit implementing the Floquet unitary designed to realize a two-dimensional discrete time crystal. c, Period-doubling oscillations of the magnetization observed on ibm_torino. The results are compared with two-dimensional tensor-network simulations (2dTNS) [from K. Shinjo, K. Seki, T. Shirakawa, R.-Y. Sun, and S. Yunoki, npj Quantim Inf. 12, 41 (2026)].
I study nonequilibrium phases and dynamical regimes emerging in quantum unitary circuits and on digital quantum computers. Key interests include:
Exotic dynamical phases and prethermal behavior in driven or Floquet-like settings
Characterizing many-body dynamics using circuit-based perspectives
Connecting circuit phenomena to condensed-matter notions (transport, relaxation, emergent structures)
K. Shinjo, K. Seki, and S. Yunoki, arXiv:2603.11910.
K. Shinjo, K. Seki, and S. Yunoki, arXiv:2510.13577.
K. Shinjo, K. Seki, T. Shirakawa, R.-Y. Sun, and S. Yunoki, npj Quantim Inf. 12, 41 (2026).
Nonequilibrium quantum many-body systems have become a central theme in modern condensed-matter physics, as they can host forms of order and universality that are fundamentally distinct from equilibrium, especially under pulse excitation or periodic driving. At the same time, real-time dynamics typically generate rapid entanglement growth, which restricts the practical reach of classical simulations based on tensor-network methods. This challenge provides a natural motivation for using quantum computers as a new computational resource. However, because today’s devices are noisy intermediate-scale quantum (NISQ) processors, increasing attention has turned to quantum-circuit dynamics themselves—including discrete-time unitary evolution and measurements—as a realistic platform for nonequilibrium physics, rather than waiting for fully error-corrected quantum algorithms.
In this context, discrete time crystals—defined by the spontaneous breaking of discrete time-translation symmetry in periodically driven (Floquet) systems—represent a prototypical example of nonequilibrium order. Their hallmark is a long-lived subharmonic response relative to the drive period, whose stability relies on mechanisms that suppress heating and stabilize the dynamics. Quantum circuits natively implement discrete-time evolution and are directly realizable on quantum hardware, providing a natural setting for connecting time-crystal physics with the NISQ era.
In our recent studies, we investigated the conditions for realizing and stabilizing discrete time crystals in quantum-circuit many-body systems, focusing on circuit design, error mitigation, and careful comparisons with classical simulations. These efforts deepen our understanding of nonequilibrium many-body dynamics while articulating a concrete strategy: quantum devices can be most impactful precisely in regimes where entanglement growth makes classical computation difficult. More broadly, our work highlights the potential of quantum circuits not only as computational tools but also as experimental platforms for synthesizing and probing novel nonequilibrium quantum phases.
Figure caption: Optical conductivity Reσ(ω) of the Hubbard model on a 6×6 square lattice, obtained by the tDMRG method. The sharp peak at the absorption edge can be interpreted as follows. When light creates a doublon–holon pair, charge motion in two dimensions inevitably disturbs the antiferromagnetic spin correlations, and this spin background provides an effective binding for the pair, producing the peak [from K. Shinjo, Y. Tamaki, S. Sota, T. Tohyama, Phys. Rev. B 104, 205123 (2021)].
Figure caption: Current-induced activation of Higgs modes in the attractive Hubbard model on a 6 × 6 square lattice, obtained using tDMRG. (a) The BCS-type Higgs mode for U=−4, activated along both the parallel and perpendicular directions relative to the pump pulse. (b) The BEC-type Higgs mode for U=−8, activated only along the direction parallel to the pump pulse [from K. Shinjo, S. Sota, S. Yunoki, and T. Tohyama, Phys. Rev. B 111, 235125 (2025)]. Here, BCS and BEC denote Bardeen–Cooper–Schrieffer superconductivity and Bose–Einstein condensation, respectively.
In parallel, I investigate nonequilibrium phenomena in strongly correlated models, especially from a pump–probe–inspired viewpoint. Topics include:
Driven dynamics and relaxation of intertwined charge, spin, and orbital degrees of freedom in Hubbard-type models
Optical/transport responses (e.g., conductivity) under nonequilibrium conditions
Building bridges between microscopic dynamics and measurable observables
K. Shinjo, S. Sota, S. Yunoki, and T. Tohyama, Phys. Rev. B 111, 235125 (2025).
K. Shinjo, S. Sota, and T. Tohyama, Phys. Rev. Res. 4, L032019 (2022).
K. Shinjo, Y. Tamaki, S. Sota, and T. Tohyama, Phys. Rev. B 104, 205123 (2021).
In strongly correlated electron systems, electrons interact with each other so strongly that they no longer behave as simple, freely moving particles. As a result, the two fundamental degrees of freedom—charge and spin—become tightly coupled, giving rise to optical responses that differ markedly from those of ordinary metals or semiconductors. In particular, in Mott insulators, which do not conduct electricity under normal conditions, light irradiation can generate doublon–holon pairs. Their motion disturbs the surrounding spin order, producing dynamics characteristic of strong electron correlations.
In this process, antiferromagnetic spin correlations provide an effective attractive interaction between the photo-generated electron and hole. This leads to a characteristic peak near the absorption edge of the optical spectrum, known as an exciton. Such a feature cannot be explained by simple interband transitions; rather, it is a direct manifestation of the cooperative motion of spin and charge. By examining the optical conductivity in detail, one can probe how electrons interact and move collectively, both in energy and in time.
In recent years, non-perturbative numerical methods such as the time-dependent density matrix renormalization group (tDMRG) have made it possible to quantitatively compute the optical response of two-dimensional Hubbard models. The resulting spectra agree well with experimental measurements in cuprate Mott insulators and related materials. This agreement provides concrete insight into how the coupling between spin and charge appears in the optical response of strongly correlated systems.
Thus, the study of optical responses in strongly correlated electron systems offers a direct way to observe and understand new quantum phenomena that emerge from collective electron behavior. It not only deepens our understanding of fundamental electronic states in solids but also opens pathways toward discovering novel optoelectronic materials and quantum functionalities.
A spin loop-current texture in the Hubbard model on a square lattice with 8×6 cluster. A next-nearest neighbor hopping t'/t=-0.266 and on-site Coulomb interaction U/t=4 (in unit of a nearest neighbor hopping t) are used. Their normalized amplitudes are shown by arrows with heatmap at each bond for hole carrier density δ=0.125 [from K. Shinjo, S. Sota, S. Yunoki, and T. Tohyama, Phys. Rev. B 108, 195118 (2023)].
Figure caption: The entanglement spectrum for the Kitaev-Heisenberg model (with a parameter α). Blue crosses represent entanglement levels and gray lines connect the spectrum belonging to the same entanglement levels. [𝑛] denotes 𝑛-fold degeneracy of the lowest entanglement level in each phase [from K. Shinjo, S. Sota and T. Tohyama, Phys. Rev. B 91, 054401 (2015)].
I am also interested in quantum spin liquids and quantum liquid-crystal phases in two-dimensional systems, where strong correlations and frustration can stabilize exotic quantum states and intertwined orders. I study how emergent excitations and collective behavior can be diagnosed using experimentally motivated probes.
K. Shinjo, S. Sota, S. Yunoki, and T. Tohyama, Phys. Rev. B 108, 195118 (2023).
K. Shinjo, S. Sota, and T. Tohyama, Phys. Rev. B 91, 054401 (2015).
In strongly correlated electron systems, strong electron–electron repulsion invalidates the single-particle picture in which electrons can be treated as independent quasiparticles. Instead, the system exhibits properties of an emergent collective quantum state. A hallmark of such systems is the appearance of novel quantum phases that are not easily characterized by conventional order parameters. Typical examples include quantum spin liquids, which lack local magnetic order yet exhibit long-range entanglement, and quantum liquid-crystal phases, which break only part of the spatial symmetry and display intermediate forms of order, such as stripe-like charge patterns or loop-current structures.
A paradigmatic example is the Kitaev model, where strong bond-directional exchange interactions generate frustration and enhance quantum fluctuations, stabilizing a spin-liquid phase. However, when additional spin interactions present in real materials are introduced as perturbations, a variety of competing magnetic orders emerge nearby in the phase diagram. Two-dimensional DMRG studies of the extended Kitaev–Heisenberg model have shown that while the spin-liquid phase exhibits characteristic pair degeneracy in the entanglement spectrum, magnetically ordered phases tend to display non-degenerate structures. Moreover, the Schmidt gap changes across phase boundaries adjacent to the spin-liquid phase, demonstrating that the structure of quantum entanglement provides a useful diagnostic for identifying Kitaev spin liquids.
When charge degrees of freedom are incorporated, the landscape becomes even richer, allowing for stripe-type charge order and superconducting states. In studies focusing on spin loop currents (sLCs), loop currents can be reinterpreted as manifestations of quantum liquid-crystal behavior. Using DMRG combined with pinning fields to visualize spatial patterns of spin currents, we investigate the conditions under which sLC correlations become long ranged. We find that such correlations grow in regimes associated with excitonic condensation or block-type antiferromagnetic order. Furthermore, in the square-lattice 𝑡-𝑡′-𝑈 model (the Hubbard model with next-nearest-neighbor hopping), at parameters relevant to high-Tc superconductors (𝑡′ / 𝑡 = − 0.25), axial sLC order emerges at hole doping 𝛿 = 0.125, coexisting with charge stripes.
Taken together, the study of novel quantum phases in strongly correlated systems aims to uncover collective quantum states that arise from the competition between strong interactions and quantum fluctuations—states that go beyond conventional “simple” orders and require genuinely quantum-mechanical descriptions.
Across these themes, I place non-perturbative numerical approaches at the core of my research strategy, as they enable the discovery and characterization of emergent phenomena beyond perturbation theory in strongly correlated quantum systems. In regimes where analytic expansions fail or become uncontrolled, such methods provide essential quantitative insight. At the same time, they clarify the parameter regimes in which quantum devices can offer genuine computational advantages. In particular, I develop and apply:
Lanczos-based exact diagonalization to resolve low-energy properties and dynamical correlation functions.
Density-matrix renormalization group (DMRG) with a focus on higher-dimensional extensions and real-time nonequilibrium dynamics
Circuit-based modeling and analysis of quantum dynamics (unitary circuits, digital quantum simulation)
Computation of response functions and nonequilibrium observables relevant to pump–probe scenarios (e.g., time-dependent optical conductivity)
Bridging unitary-circuit perspectives with microscopic correlated-electron models to assess robustness, scalability, and experimental relevance
Beyond numerical analysis, I also pursue the construction of phenomenological and/or effective models to uncover the theoretical structures underlying the observed numerical results. This parallel effort aims to extract essential mechanisms and universal features, elevating numerical findings into conceptually transparent theoretical frameworks.
Left: AI-generated concept image