Published papers

We use Monte Carlo simulations and tensor network methods to study a classical monomer-dimer model on the square lattice with a hole (monomer) fugacity z, an aligning dimer-dimer interaction u that favors columnar order, and an attractive dimer-dimer interaction v between two adjacent dimers that lie on the same principal axis of the lattice. The Monte Carlo simulations of finite size systems rely on our grand-canonical generalization of the dimer worm algorithm, while the tensor network computations are based on a uniform matrix product ansatz for the eigenvector of the row-to-row transfer matrix that work directly in the thermodynamic limit. The phase diagram has nematic, columnar order and fluid phases, and a nonzero temperature multicritical point at which all three meet. For any fixed v/u<∞, we argue that this multicritical point continues to be located at a nonzero hole fugacity zmc(v/u)>0; our numerical results confirm this theoretical expectation but find that zmc(v/u)→0 very rapidly as v/u→∞. Our numerical results also confirm the theoretical expectation that the corresponding multicritical behavior is in the universality class of the four-state Potts multicritical point on critical line of the two-dimensional Ashkin-Teller model.

We investigate a quantum phase transition between a Néel phase and a valence bond solid (VBS) phase, each of which breaks a different Z2 symmetry, in a spin-1/2 two-leg XXZ ladder with a four-spin interaction. The model can be viewed as a one-dimensional variant of the celebrated J-Q model on a square lattice. By means of variational uniform matrix product state calculations and an effective field theory, we determine the phase diagram of the model and present evidence that the Néel-VBS transition is continuous and belongs to the Gaussian universality class with the central charge c=1. In particular, the critical exponents β,η, and, ν are found to satisfy the constraints expected for a Gaussian transition within numerical accuracy. These exponents do not detectably change along the phase boundary while they are in general allowed to do so for the Gaussian class.

A tensor network renormalization algorithm with global optimization based on the corner transfer matrix is proposed. Since the environment is updated by the corner transfer matrix renormalization group method, the forward-backward iteration is unnecessary, which is a time-consuming part of other methods with global optimization. In addition, a further approximation reducing the order of the computational cost of contraction for the calculation of the coarse-grained tensor is proposed. The computational time of our algorithm in two dimensions scales as the sixth power of the bond dimension, while the higher-order tensor renormalization group and the higher-order second renormalization group methods have the seventh power. We perform benchmark calculations in the Ising model on the square lattice and show that the time-to-solution of the proposed algorithm is faster than that of other methods.

We study the ground-state phase diagram of a spin-12 frustrated XXZ ladder, in which two antiferromagnetic chains are coupled by competing rung and diagonal interactions, J⊥ and J×. Previous studies on the isotropic model have revealed that a fluctuation-induced effective dimer attraction between the legs stabilizes the columnar dimer (CD) phase in the highly frustrated regime J⊥≈2J×, especially for ferromagnetic J⊥,×<0. By means of effective field theory and numerical analyses, we extend this analysis to the XXZ model, and we obtain a rich phase diagram. The diagram includes four gapped featureless phases with no symmetry breaking: the rung singlet (RS) and Haldane phases as well as their twisted variants, the RS∗ and Haldane∗ phases, which are all distinct in the presence of certain symmetries. Significantly, the Haldane-CD transition point in the isotropic model turns out to be a crossing point of two transition lines in the XXZ model, and the stripe Néel and RS∗ phases appear between these lines. This indicates a nontrivial interplay between the effective dimer attraction and the exchange anisotropy. In the easy-plane regime, the four featureless phases and two critical phases are found to compete in a complex manner depending on the signs of J⊥,×.

We propose a class of tight-binding systems of interacting bosons with a flat band, which are exactly solvable in the sense that one can explicitly write down the unique ground state. The ground state is expressed in terms of local creation operators and apparently resembles that of a Mott insulator. Based on an exact representation in terms of a classical loop-gas model, we conjecture that the ground state may exhibit quasi Bose-Einstein condensation (BEC) and genuine BEC in dimension 2 and in dimension 3 or higher, respectively, still keeping Mott-insulator-like character. Our Monte Carlo simulation of the loop-gas model strongly supports this conjecture, i.e., the ground state undergoes a Kosterlitz-Thouless transition and exhibits quasi-BEC in two dimensions.

We propose a multi-impurity method for the bond-weighted tensor renormalization group (BWTRG) to compute the higher-order moment of physical quantities in a two-dimensional system. The replacement of the bond weight with an impurity matrix in a bond-weighted triad tensor network represents a physical quantity such as the magnetization and the energy. We demonstrate that the accuracy of the proposed method is much higher than the conventional tensor renormalization group for the Ising model and the five-state Potts model. Furthermore, we perform the finite-size scaling analysis and observe that the dimensionless quantity characterizing the structure of the fixed point tensor satisfies the same scaling relation in the critical region as the Binder parameter. The estimated critical temperature dependence on the bond dimension indicates that the exponent relating the correlation length to the bond dimension varies continuously with respect to the BWTRG hyperparameter. We find that BWTRG with the optimal hyperparameter is more efficient in terms of computational time than alternative approaches based on the matrix product state in estimating the critical temperature.

We study the magnetization process of the S=1 Heisenberg model on a two-leg ladder with farther neighbor spin-exchange interaction. We consider the interaction that couples up to the next-nearest neighbor rungs and find an exactly solvable regime where the ground states become product states. The next-nearest neighbor interaction tends to stabilize magnetization plateaus at multiples of 1/6. In most of the exactly solvable regime, a single magnetization curve shows two series of plateaus with different periodicities.

Quantum mechanical control of the position of a particle by using a trapping potential well is an important problem for the manipulation of a quantum particle. We study the probability of successful conveyance of the particle trapped in a potential well for a given length within a given fixed time, i.e., survival probability after the motion. For the actual motion of conveyance, we need to accelerate the particle to move and then decelerate it to stop at the destination. Acceleration and deceleration cause a dropoff of a particle from the trapping potential well. The relaxation of the survival probability in a process with a constant acceleration rate is studied in detail. First, the process is studied as the relaxation of the survival probability of the trapped particle by direct numerical calculations. The survival probability obeys an exponential decay in a long time, which is analyzed from a viewpoint of eigenvalue problem. The value of survival probability is also estimated by the Wentzel-Kramers-Brillouin method with connection formulas using the Airy function and the Weber function. The value is further estimated by a method of the resonance states. We emphasize the fact that an important source of dropoff comes from a nonanalytic change of velocity at the starting point. When the rested particle begins to move, the ground state of the rest frame is redistributed to eigenstates of the moving frame, and then each eigenstate of the moving frame evolves in time. The dephasing of wave functions of the distributed populations reduces the probability of successful conveyance. In general, a smooth start gives a small initial disturbance but it requires a large acceleration during the process to reach the destination in the fixed time which causes a larger dropoff in the process. Considering these conflicting facts, we study the survival probability in concrete conveyance schemes, i.e., protocols with (1) a sudden change of velocity to a constant velocity, (2) a sudden change of acceleration rate to a constant acceleration, and (3) a smooth change of acceleration by studying the real-time change of populations of the adiabatic (instantaneous) eigenstates. We observe the time evolution of the trapped probability and the population distribution during the conveyance process. In cases that the potential well has several bound states, we propose a method to select the particle trapped at the ground state by making use of the difference of survival probabilities of bound states.

We study a spin-12 XXZ model with a four-spin interaction on a two-leg ladder. By means of effective field theory and matrix product state calculations, we obtain rich ground-state phase diagrams that consist of eight distinct gapped phases. Four of them exhibit spontaneous symmetry breaking with either a magnetic or valence-bond-solid (VBS) long-range order. The other four are featureless, i.e., the bulk ground state is unique and does not break any symmetry. The featureless phases include the rung singlet (RS) and Haldane phases as well as their variants, the RS∗ and Haldane∗ phases, in which twisted singlet pairs (|↑↓)+|↓↑))/2 are formed between the two legs. We argue and demonstrate that Gaussian transitions with the central charge c=1 occur between the featureless phases and between the ordered phases while Ising transitions with c=1/2 occur between the featureless and ordered phases. The two types of transition lines cross at the SU(2)-symmetric point, where the criticality is described by the SU(2)2 Wess-Zumino-Witten theory with c=3/2. The RS-Haldane∗ and RS∗-Haldane transitions give examples of topological phase transitions. Interestingly, the RS∗ and Haldane∗ phases, which have highly anisotropic nature, appear even in the vicinity of the isotropic case. We demonstrate that all the four featureless phases are distinguished by topological indices in the presence of certain symmetries.

Quantum many-body systems are challenging targets for computational physics due to their large number of degrees of freedom. The tensor networks, particularly Tensor Product States (TPS) and Projected Entangled Pair States (PEPS), effectively represent these systems on two-dimensional lattices. However, the technical complexity of TPS/PEPS-based coding can be challenging for many researchers to manage effectively. To reduce this problem, we developed TeNeS (Tensor Network Solver). This paper introduces TeNeS-v2, which extends TeNeS with real-time and finite temperature simulations, providing deeper insights into quantum many-body systems. We detail the new algorithms, input/output design, and application examples, demonstrating TeNeS-v2's applicability to various quantum spin and Bose models on two-dimensional lattices. New version program summary: Program Title: TeNeS CPC Library link to program files: https://doi.org/10.17632/psm26xxbvd.2 Developer's repository link: https://github.com/issp-center-dev/TeNeS Licensing provisions: GPLv3 Programming language: C++11 Journal reference of previous version: Comput. Phys. Commun. (2022) 279, 108437, doi:10.1016/j.cpc.2022.108437 Reasons for the new version: To extend the capabilities of TeNeS to include real-time and finite temperature simulations, enabling deeper insights into quantum many-body systems beyond ground states. Summary of revisions: TeNeS-v2 introduces real-time evolution and finite temperature simulation capabilities, expanding the range of quantum many-body phenomena that can be studied. Nature of problem: Quantum many-body systems are extremely difficult to simulate because of the huge dimension of the Hilbert space. Conventional methods have difficulty in accurately representing these systems, especially for properties beyond small lattices and ground states. Advanced computational techniques are required to efficiently capture the dynamics and thermal properties of these systems. Solution method: TeNeS employs tensor networks, specifically Tensor Product States (TPS) and Projected Entangled Pair States (PEPS), to efficiently represent quantum many-body states on two-dimensional lattices. The real-time evolution is handled through a straightforward extension of imaginary time evolution methods, allowing the study of dynamical properties. Finite temperature simulations are conducted using an imaginary time evolution approach starting from the infinite-temperature mixed state. These methods enable the accurate calculation of various physical properties in different quantum states. Additional comments including restrictions and unusual features: TeNeS-v2 maintains ease of use by providing flexible input file configurations and supporting a variety of two-dimensional lattice models. Increasing the bond dimension to improve accuracy requires more computational resources. Finite temperature calculations may exhibit unphysical behavior due to limitations of the tensor network approximation method. Users should be aware of these potential problems and interpret the results following the suggestions made in the text.

TeNeS (Tensor Network Solver) [1,2] is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and correlation functions. For certain predefined models and lattices, there is a tool that makes it easy to generate input files. TeNeS uses an OpenMP/MPI hybrid parallelized tensor operation library and thus can perform large-scale calculations using massively parallel machines. Program summary: Program Title: TeNeS CPC Library link to program files: https://doi.org/10.17632/psm26xxbvd.1 Code Ocean capsule: https://codeocean.com/capsule/1510058 Licensing provisions: GNU General Public License version 3 Programming language: C++ and python3 External routines/libraries: mptensor Nature of problem: TeNeS calculates the approximate ground states and their properties of user-defined two-dimensional quantum lattice models using user-friendly input files. Numerically exact solutions of such tasks generally require an exponentially diverging computational time, whereas the error in the output of TeNeS is well controlled and can be reduced with a polynomial cost. Solution method: TeNeS implements the tensor networks method based on a tensor-product-state (TPS) wavefunction and the corner transfer matrix renormalization group method. TeNeS also supports massively parallel computing using the library mptensor, which implements OpenMP/MPI hybrid parallelized tensor operations.

Using the nuclear norm regularization techniques on tensor network renormalization algorithm, we study the phase diagram, the critical behavior, and the duality property of the antiferromagnetic 6-state clock model on the Union Jack lattice. We find that this model undergoes multiple phase transitions; there is the Berezinskii-Kosterlitz-Thouless, Z6 symmetry breaking, and chiral transition with decreasing temperature. Furthermore, we provide convincing numerical evidence that its quasi-long-range order is well explained by the compactified boson conformal field theory (CFT) and the chiral transition is in perfect agreement with the Ising CFT, including central charge, scaling dimension spectrum, and operator product expansion coefficients.

The behavior of dimensionless quantities defined as ratios of partition functions is analyzed to investigate phase transitions and critical phenomena. At criticality, the universal values of these ratios can be predicted from conformal field theory (CFT) through the modular-invariant partition functions on a torus. We perform numerical calculations using the bond-weighted tensor renormalization group for three two-dimensional models belonging to different universality classes: the Ising model, the three-state Potts model, and the four-state Potts model. The partition-function ratios obey the same finite-size scaling form as the Binder parameter, and their critical values agree well with the universal values predicted by CFT. In the four-state Potts model, we observe logarithmic corrections in the system-size dependence of these ratios.