Showing 1–41 of 41 results for author: Xie, Z Y

The Green's function Monte Carlo combined with projected entangled pair state approach to the frustrated $J_1$-$J_2$ Heisenberg model

Authors: He-Yu Lin, Yibin Guo, Rong-Qiang He, Z. Y. Xie, Zhong-Yi Lu

Abstract: The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enha… ▽ More The tensor network algorithm, a family of prevalent numerical methods for quantum many-body problems, aptly captures the entanglement properties intrinsic to quantum systems, enabling precise representation of quantum states. However, its computational cost is notably high, particularly in calculating physical observables like correlation functions. To surmount the computational challenge and enhance efficiency, we propose integrating the Green's function Monte Carlo (GFMC) method with the projected entangled pair state (PEPS) ansatz. This approach combines the high-efficiency characteristics of Monte Carlo with the sign-free nature of tensor network states and proves effective in addressing the computational bottleneck. To showcase its prowess, we apply this hybrid approach to investigate the antiferromagnetic $J_1$-$J_2$ Heisenberg model on the square lattice, a model notorious for its sign problem in quantum Monte Carlo simulations. Our results reveal a substantial improvement in the accuracy of ground-state energy when utilizing a preliminary PEPS as the guiding wave function for GFMC. By calculating the structure factor and spin-spin correlation functions, we further characterize the phase diagram, identifying a possible columnar valence-bond state phase within the intermediate parameter range of $0.52 < J_2/J_1 < 0.58$. This comprehensive study underscores the efficacy of our combined approach, demonstrating its ability to accurately simulate frustrated quantum spin systems while ensuring computational efficiency. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 11 pages, 15 figures

Journal ref: Phys. Rev. B 109, 235133 (2024)

Susceptibility indicator for chiral topological orders emergent from correlated fermions

Authors: Rui Wang, Tao Yang, Z. Y. Xie, Baigeng Wang, X. C. Xie

Abstract: Chiral topological orders formed in correlated fermion systems have been widely explored. However, the mechanism on how they emerge from interacting fermions is still unclear. Here, we propose a susceptibility condition. Under this condition, we show that chiral topological orders can spontaneously take place in correlated fermion systems. The condition leads to a low-energy effective theory of bo… ▽ More Chiral topological orders formed in correlated fermion systems have been widely explored. However, the mechanism on how they emerge from interacting fermions is still unclear. Here, we propose a susceptibility condition. Under this condition, we show that chiral topological orders can spontaneously take place in correlated fermion systems. The condition leads to a low-energy effective theory of bosons with strong frustration, mimicking the flat band systems. The frustration then melts the long-range orders and results in topological orders with time-reversal symmetry breaking. We apply the theory to strongly-correlated semiconductors doped to the metallic phase. A novel excitonic topological order with semionic excitations and chiral excitonic edge state is revealed. We also discuss the application to frustrated magnets. The theory predicts a chiral spin liquid state, which is numerically confirmed by our tensor network calculations. These results demonstrate an unprecedented indicator for chiral topological orders, which bridges the existing gap between interacting fermions and correlated topological matter. △ Less

Submitted 23 June, 2024; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: 22 pages, 6 figures

Journal ref: Phys.Rev.B.109, L241113 (2024)

Structural, electronic, magnetic properties of Cu-doped lead-apatite Pb$_{10-x}$Cu$_x$(PO$_4$)$_6$O

Authors: Jianfeng Zhang, Huazhen Li, Ming Yan, Miao Gao, Fengjie Ma, Xun-Wang Yan, Z. Y. Xie

Abstract: The recent report of superconductivity in the Cu-doped PbPO compound stimulates the extensive researches on its physical properties. Herein, the detailed atomic and electronic structures of this compound are investigated, which are the necessary information to explain the physical properties, including possible superconductivity. By the first-principles electronic structure calculations, we find t… ▽ More The recent report of superconductivity in the Cu-doped PbPO compound stimulates the extensive researches on its physical properties. Herein, the detailed atomic and electronic structures of this compound are investigated, which are the necessary information to explain the physical properties, including possible superconductivity. By the first-principles electronic structure calculations, we find that the partial replacement of Pb at $4f$ site by Cu atom, instead of Pb at $6h$ site, plays a crucial role in dominating the electronic state at Fermi energy. The $3d$ electronic orbitals of Cu atom emerge near the Fermi energy and exhibit strong spin-polarization, resulting in the local moment around the doped Cu atom. Particularly, the ground state of Pb$_{10-x}$Cu$_x$(PO$_4$)$_6$O (x = 1) is determined to be a semiconducting phase, in good agreement with the experimental measurements. △ Less

Submitted 8 August, 2023; originally announced August 2023.

High-temperature ferromagnetism and strong $π$-conjugation feature in two-dimensional manganese tetranitride

Authors: Ming Yan, Z. Y. Xie, Miao Gao

Abstract: Two-dimensional (2D) magnetic materials have attracted tremendous research interest because of the promising application in the next-generation microelectronic devices. Here, by the first-principles calculations, we propose a two-dimensional ferromagnetic material with high Curie temperature, manganese tetranitride MnN$_4$ monolayer, which is a square-planar lattice made up of only one layer of at… ▽ More Two-dimensional (2D) magnetic materials have attracted tremendous research interest because of the promising application in the next-generation microelectronic devices. Here, by the first-principles calculations, we propose a two-dimensional ferromagnetic material with high Curie temperature, manganese tetranitride MnN$_4$ monolayer, which is a square-planar lattice made up of only one layer of atoms. The structure is demonstrated to be stable by the phonon spectra and the molecular dynamic simulations, and the stability is ascribed to the $π$-d conjugation between $π$ orbital of N=N bond and Mn $d$ orbital. More interestingly, the MnN$_4$ monolayer displays robust 2D ferromagnetism, which originates from the strong exchange couplings between Mn atoms due to the $π$-d conjugation. The high critical temperature of 247 K is determined by solving the Heisenberg model with the Monte Carlo method. △ Less

Submitted 27 November, 2022; originally announced November 2022.

Journal ref: Chin. Phys. B 32, 037104 (2023)

Long-range spin-orbital order in the spin-orbital SU(2)$\times$SU(2)$\times$U(1) model

Authors: Yang Liu, Z. Y. Xie, Hong-Gang Luo, Jize Zhao

Abstract: By using the tensor-network state algorithm, we study a spin-orbital model with SU(2)$\times$SU(2)$\times$U(1) symmetry on the triangular lattice. This model was proposed to describe some triangular $d^1$ materials and was argued to host a spin-orbital liquid ground state. In our work the trial wavefunction of its ground state is approximated by an infinite projected entangled simplex state and op… ▽ More By using the tensor-network state algorithm, we study a spin-orbital model with SU(2)$\times$SU(2)$\times$U(1) symmetry on the triangular lattice. This model was proposed to describe some triangular $d^1$ materials and was argued to host a spin-orbital liquid ground state. In our work the trial wavefunction of its ground state is approximated by an infinite projected entangled simplex state and optimized by the imaginary-time evolution. Contrary to the previous conjecture, we find that the two SU(2) symmetries are broken, resulting in a stripe spin-orbital order with the same magnitude $m=0.085(10)$. This value is about half of that in the spin-1/2 triangular Heisenberg antiferromagnet. Our result demonstrates that although the long-sought spin-orbital liquid is absent in this model the spin-orbital order is significantly reduced due to the enhanced quantum fluctuation. This suggests that high-symmetry spin-orbital models are promising in searching for exotic states of matter in condensed-matter physics. △ Less

Submitted 13 November, 2022; originally announced November 2022.

Comments: 4 pages, 5 figures

Journal ref: Phys. Rev. B 107, L041106 (2023)

Exploring explicit coarse-grained structure in artificial neural networks

Authors: Xi-Ci Yang, Z. Y. Xie, Xiao-Tao Yang

Abstract: We propose to employ the hierarchical coarse-grained structure in the artificial neural networks explicitly to improve the interpretability without degrading performance. The idea has been applied in two situations. One is a neural network called TaylorNet, which aims to approximate the general mapping from input data to output result in terms of Taylor series directly, without resorting to any ma… ▽ More We propose to employ the hierarchical coarse-grained structure in the artificial neural networks explicitly to improve the interpretability without degrading performance. The idea has been applied in two situations. One is a neural network called TaylorNet, which aims to approximate the general mapping from input data to output result in terms of Taylor series directly, without resorting to any magic nonlinear activations. The other is a new setup for data distillation, which can perform multi-level abstraction of the input dataset and generate new data that possesses the relevant features of the original dataset and can be used as references for classification. In both cases, the coarse-grained structure plays an important role in simplifying the network and improving both the interpretability and efficiency. The validity has been demonstrated on MNIST and CIFAR-10 datasets. Further improvement and some open questions related are also discussed. △ Less

Submitted 4 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Journal ref: Chin. Phys. Lett. 40, 020501 (2023)

Efficient calculation of three-dimensional tensor networks

Authors: Li-Ping Yang, Y. F. Fu, Z. Y. Xie, T. Xiang

Abstract: We have proposed an efficient algorithm to calculate physical quantities in the translational invariant three-dimensional tensor networks, which is particularly relevant to the study of the three-dimensional classical statistical models and the (2+1)-dimensional quantum lattice models. In the context of a classical model, we determine the partition function by solving the dominant eigenvalue probl… ▽ More We have proposed an efficient algorithm to calculate physical quantities in the translational invariant three-dimensional tensor networks, which is particularly relevant to the study of the three-dimensional classical statistical models and the (2+1)-dimensional quantum lattice models. In the context of a classical model, we determine the partition function by solving the dominant eigenvalue problem of the transfer matrix, whose left and right dominant eigenvectors are represented by two projected entangled simplex states. These two projected entangled simplex states are not Hermitian conjugate to each other but are appropriately arranged so that their inner product can be computed much more efficiently than in the usual prescription. For the three-dimensional Ising model, the calculated internal energy and spontaneous magnetization agree with the published results in the literature. The possible improvement and extension to other models are also discussed. △ Less

Submitted 14 April, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 8 pages, 10 figures

Journal ref: Phys. Rev. B 107, 165127 (2023)

Exploring quantum phase transitions by the cross derivative of the ground state energy

Authors: H. Y. Wu, Yu-Chin Tzeng, Z. Y. Xie, K. Ji, J. F. Yu

Abstract: In this work, the cross derivative of the Gibbs free energy, initially proposed for phase transitions in classical spin models [Phys. Rev. B 101, 165123 (2020)], is extended for quantum systems. We take the spin-1 XXZ chain with anisotropies as an example to demonstrate its effectiveness and convenience for the Gaussian-type quantum phase transitions therein. These higher-order transitions are ver… ▽ More In this work, the cross derivative of the Gibbs free energy, initially proposed for phase transitions in classical spin models [Phys. Rev. B 101, 165123 (2020)], is extended for quantum systems. We take the spin-1 XXZ chain with anisotropies as an example to demonstrate its effectiveness and convenience for the Gaussian-type quantum phase transitions therein. These higher-order transitions are very challenging to determine by conventional methods. From the cross derivative with respect to the two anisotropic strengths, a single valley structure is observed clearly in each system size. The finite-size extrapolation of the valley depth shows a perfect logarithmic divergence, signaling the onset of a phase transition. Meanwhile, the critical point and the critical exponent for the correlation length are obtained by a power-law fitting of the valley location in each size. The results are well consistent with the best estimations in the literature. Its application to other quantum systems with continuous phase transitions is also discussed briefly. △ Less

Submitted 11 April, 2023; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 7 pages, 7 figures

Journal ref: New J. Phys. 25 (2023) 043006

Two-dimensional binary transition metal nitride $M$N$_4$ ($M$ = V, Cr, Mn, Fe, Co) with a graphene-like structure and strong magnetic properties

Authors: Shuo Zhang, Panjun Feng, Dapeng Liu, Hongfen Wu, Miao Gao, Tongshuai Xu, Xun-Wang Yan, Z. Y. Xie

Abstract: Binary transition metal nitride with a graphene-like structure and strong magnetic properties is rare. Based on the first-principles calculations, we design two kinds of $M$N$_4$ ($M$ =transition metal) monolayers, which are transition metal nitrides with a planar structure, made up of $M$N$_4$ units aligned in the rhombic and square patterns. The two structural lattices have robust stability and… ▽ More Binary transition metal nitride with a graphene-like structure and strong magnetic properties is rare. Based on the first-principles calculations, we design two kinds of $M$N$_4$ ($M$ =transition metal) monolayers, which are transition metal nitrides with a planar structure, made up of $M$N$_4$ units aligned in the rhombic and square patterns. The two structural lattices have robust stability and good compatibility with different metal atoms, and the underlying mechanism is the combination of $sp^2$ hybridization, coordinate bond, and $π$ conjugation. With the metal atom changing from V, Cr, Mn, Fe to Co, the total charge of $M$N$_4$ system increases by one electron in turn, which results in continuous adjustability of the electronic and magnetic properties. The planar ligand field is another feature of the two $M$N$_4$ lattices, which brings about the special splitting of five suborbitals of 3$d$ metal atom and gives rise to strong magnetism. Moreover, room-temperature ferromagnetism in square-CoN$_4$ monolayer with the Curie temperatures of 321 K is determined by solving the Heisenberg model combined with Monte Carlo method. △ Less

Submitted 3 October, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Journal ref: Phys. Rev. B 106, 235402 (2022)

High-temperature ferromagnetism in two-dimensional material MnSn originated from interlayer coupling

Authors: Panjun Feng, Xiaohui Zhang, Dapeng Liu, Shuo Zhang, Xun-Wang Yan, Z. Y. Xie

Abstract: MnSn monolayer synthesized recently is a novel two-dimensional ferromagnetic material with a hexagonal lattice, in which three Mn atom come together to form a trimer, making it remarkably different from other magnetic two-dimensional materials. Most impressively, there happens a sharp increase of Curie temperature from 54 K to 225 K when the number of layers increase from 1 to 3. However, no quant… ▽ More MnSn monolayer synthesized recently is a novel two-dimensional ferromagnetic material with a hexagonal lattice, in which three Mn atom come together to form a trimer, making it remarkably different from other magnetic two-dimensional materials. Most impressively, there happens a sharp increase of Curie temperature from 54 K to 225 K when the number of layers increase from 1 to 3. However, no quantitative explanation is reported in previous studies. Herein, by means of first-principle calculations method and Monte carlo method, we demonstrate that strong interlayer ferromagnetic coupling is the essential role in enhancing its critical temperature, which act as a magnetic field to stabilize the ferromagnetism in the MnSn multilayers. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Journal ref: ACS Omega 7, 43316 (2022)

Variational corner transfer matrix renormalization group method for classical statistical models

Authors: X. F. Liu, Y. F. Fu, W. Q. Yu, J. F. Yu, Z. Y. Xie

Abstract: In the context of tensor network states, we for the first time reformulate the corner transfer matrix renormalization group (CTMRG) method into a variational bilevel optimization algorithm. The solution of the optimization problem corresponds to the fixed-point environment pursued in the conventional CTMRG method, from which the partition function of a classical statistical model, represented by a… ▽ More In the context of tensor network states, we for the first time reformulate the corner transfer matrix renormalization group (CTMRG) method into a variational bilevel optimization algorithm. The solution of the optimization problem corresponds to the fixed-point environment pursued in the conventional CTMRG method, from which the partition function of a classical statistical model, represented by an infinite tensor network, can be efficiently evaluated. The validity of this variational idea is demonstrated by the high-precision calculation of the residual entropy of the dimer model, and is further verified by investigating several typical phase transitions in classical spin models, where the obtained critical points and critical exponents all agree with the best known results in literature. Its extension to three-dimensional tensor networks or quantum lattice models is straightforward, as also discussed briefly. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: 7 pages, 7 figures

Journal ref: Chin. Phys. Lett. 39, 067502 (2022)

Prediction of single-atom-thick transition metal nitride CrN$_4$ with a square-planar network and high-temperature ferromagnetism

Authors: Dapeng Liu, Panjun Feng, Shuo Zhang, Miao Gao, Fengjie Ma, Xun-Wang Yan, Z. Y. Xie

Abstract: Single-atom-thick two-dimensional materials such as graphene usually have a hexagonal lattice while the square-planar lattice is uncommon in the family of two-dimensional materials. Here, we demonstrate that single-atom-thick transition metal nitride CrN$_4$ monolayer is a stable free-standing layer with a square-planar network. The stability of square-planar geometry is ascribed to the combinat… ▽ More Single-atom-thick two-dimensional materials such as graphene usually have a hexagonal lattice while the square-planar lattice is uncommon in the family of two-dimensional materials. Here, we demonstrate that single-atom-thick transition metal nitride CrN$_4$ monolayer is a stable free-standing layer with a square-planar network. The stability of square-planar geometry is ascribed to the combination of N=N double bond, Cr-N coordination bond, and $π$-d conjugation, in which the double $π$-d conjugation is rarely reported in previous studies. This mechanism is entirely different from that of the reported two-dimensional materials, leading to lower formation energy and more robust stability compared to the synthesized g-C$_3$N$_4$ monolayer. On the other hand, CrN$_4$ layer has a ferromagnetic ground state, in which the ferromagnetic coupling between two Cr atoms is mediated by electrons of the half-filled large $π$ orbitals from $π$-d conjugation. The high-temperature ferromagnetism in CrN$_4$ monolayer is confirmed by solving the Heisenberg model with Monte Carlo method. △ Less

Submitted 22 September, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: https://link.aps.org/doi/10.1103/PhysRevB.106.125421

Journal ref: Physical Review B 106, 125421 (2022)

Achieving high-temperature ferromagnetism by means of magnetic ions dimerization

Authors: Panjun Feng, Shuo Zhang, Dapeng Liu, Miao Gao, Fengjie Ma, Xun-Wang Yan, Z. Y. Xie

Abstract: Magnetic two-dimensional materials have potential application in next-generation electronic devices and have stimulated extensive interest in condensed matter physics and material fields. However, how to realize high-temperature ferromagnetism in two-dimensional materials remains a great challenge in physics. Herein, we propose an effective approach that the dimerization of magnetic ions in two-… ▽ More Magnetic two-dimensional materials have potential application in next-generation electronic devices and have stimulated extensive interest in condensed matter physics and material fields. However, how to realize high-temperature ferromagnetism in two-dimensional materials remains a great challenge in physics. Herein, we propose an effective approach that the dimerization of magnetic ions in two-dimensional materials can enhance the exchange coupling and stabilize the ferromagnetism. Manganese carbonitride Mn$_2$N$_6$C$_6$ with a planar monolayer structure is taken as an example to clarify the method, in which two Mn atoms are gathered together to form a ferromagnetic dimer of Mn atoms and further these dimers are coupled together to form the overall ferromagnetism of the two-dimensional material. In Mn$_2$N$_6$C$_6$ monolayer, the near-room-temperature ferromagnetism with the Curie temperature of 272.3 K is determined by solving Heisenberg model using Monte Carlo simulations method. △ Less

Submitted 12 June, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: Journal of Physical Chemistry C Cite This: https://doi.org/10.1021/acs.jpcc.2c02593

Journal ref: J. Phys. Chem. C 126, 10139 (2022)

Laboratory observation of plasmoid-dominated magnetic reconnection in hybrid collisional-collisionless regime

Authors: Z. H. Zhao, H. H. An, Y. Xie, Z. Lei, W. P. Yao, W. Q. Yuan, J. Xiong, C. Wang, J. J. Ye, Z. Y. Xie, Z. H. Fang, A. L. Lei, W. B. Pei, X. T. He, W. M. Zhou, W. Wang, S. P. Zhu, B. Qiao

Abstract: Magnetic reconnection, breaking and reorganization of magnetic field topology, is a fundamental process for rapid release of magnetic energy into plasma particles that occurs pervasively throughout the universe. In most natural circumstances, the plasma properties on either side of the reconnection layer are asymmetric, in particular for the collision rates that are associated with a combination o… ▽ More Magnetic reconnection, breaking and reorganization of magnetic field topology, is a fundamental process for rapid release of magnetic energy into plasma particles that occurs pervasively throughout the universe. In most natural circumstances, the plasma properties on either side of the reconnection layer are asymmetric, in particular for the collision rates that are associated with a combination of density and temperature and critically determine the reconnection mechanism. To date, all laboratory experiments on magnetic reconnections have been limited to purely collisional or collisionless regimes. Here, we report a well-designed experimental investigation on asymmetric magnetic reconnections in a novel hybrid collisional-collisionless regime by interactions between laser-ablated Cu and CH plasmas. We show that the growth rate of the tearing instability in such a hybrid regime is still extremely large, resulting in rapid formation of multiple plasmoids, lower than that in the purely collisionless regime but much higher than the collisional case. In addition, we, for the first time, directly observe the topology evolutions of the whole process of plasmoid-dominated magnetic reconnections by using highly-resolved proton radiography. △ Less

Submitted 24 February, 2022; originally announced February 2022.

First-order transition between the plaquette valence bond solid and antiferromagnetic phases of the Shastry-Sutherland model

Authors: Ning Xi, Hongyu Chen, Z. Y. Xie, Rong Yu

Abstract: We study the ground state phase diagram of the Shastry-Sutherland model by using the variational optimization of the infinite tensor network states, and find a weakly first-order transition between the plaquette and the antiferromagnetic states. The full plaquette state strongly competes with the empty plaquette ground state, with an energy difference less than $10^{-4}J$. We show a staggered ring… ▽ More We study the ground state phase diagram of the Shastry-Sutherland model by using the variational optimization of the infinite tensor network states, and find a weakly first-order transition between the plaquette and the antiferromagnetic states. The full plaquette state strongly competes with the empty plaquette ground state, with an energy difference less than $10^{-4}J$. We show a staggered ring exchange interaction that preserves the Shastry-Sutherland lattice symmetry can stabilize the full plaquette ground state. In light of this, we propose the triple point where the full plaquette, empty plaquette, and antiferromagnetic phases meet as a deconfined quantum critical point. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: 5.5+2 pages, 4+3 figures

Journal ref: Phys. Rev. B 107, L220408 (2023)

Enabling Lightweight Fine-tuning for Pre-trained Language Model Compression based on Matrix Product Operators

Authors: Peiyu Liu, Ze-Feng Gao, Wayne Xin Zhao, Z. Y. Xie, Zhong-Yi Lu, Ji-Rong Wen

Abstract: This paper presents a novel pre-trained language models (PLM) compression approach based on the matrix product operator (short as MPO) from quantum many-body physics. It can decompose an original matrix into central tensors (containing the core information) and auxiliary tensors (with only a small proportion of parameters). With the decomposed MPO structure, we propose a novel fine-tuning strategy… ▽ More This paper presents a novel pre-trained language models (PLM) compression approach based on the matrix product operator (short as MPO) from quantum many-body physics. It can decompose an original matrix into central tensors (containing the core information) and auxiliary tensors (with only a small proportion of parameters). With the decomposed MPO structure, we propose a novel fine-tuning strategy by only updating the parameters from the auxiliary tensors, and design an optimization algorithm for MPO-based approximation over stacked network architectures. Our approach can be applied to the original or the compressed PLMs in a general way, which derives a lighter network and significantly reduces the parameters to be fine-tuned. Extensive experiments have demonstrated the effectiveness of the proposed approach in model compression, especially the reduction in finetuning parameters (91% reduction on average). △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: Accepted by ACL 2021 main conference

Theory of competing Chern-Simons orders and emergent phase transitions

Authors: Rui Wang, Z. Y. Xie, Baigeng Wang, Tigran Sedrakyan

Abstract: Based on the Chern-Simons fermionization of spin-1/2 operators, we propose a systematic framework to investigate the competition between emergent phases in frustrated two-dimensional XY quantum magnets. Application of the method to the antiferromagnetic honeycomb spin-1/2 $J_1$-$J_2$ XY model reveals an unconventional phase transition between two Chern-Simons orders: the Chern-Simons superconducto… ▽ More Based on the Chern-Simons fermionization of spin-1/2 operators, we propose a systematic framework to investigate the competition between emergent phases in frustrated two-dimensional XY quantum magnets. Application of the method to the antiferromagnetic honeycomb spin-1/2 $J_1$-$J_2$ XY model reveals an unconventional phase transition between two Chern-Simons orders: the Chern-Simons superconductor and the exciton insulator of Chern-Simons fermions. We show that in the spin language, this transition translates to the transition from the planar Néel state to the non-uniform chiral spin-liquid that was proposed earlier in the literature. Namely, the Chern-Simons superconductor describes the planar Néel state, while the Chern-Simons exciton insulator corresponds to the non-uniform chiral spin-liquid. These results are further confirmed by our high-precision tensor network calculation, which provides the first numerical evidence for the transition from Néel order to a non-uniform chiral spin-liquid. We argue that the developed method can be applied to other frustrated quantum magnets of XXZ type and can detect unconventional phase transitions. △ Less

Submitted 20 January, 2021; v1 submitted 12 January, 2021; originally announced January 2021.

Comments: 15 pages, 8 figures, supplemental material included

Journal ref: Phys. Rev. B 106, L121117 (2022)

Magnetization of the spin-1/2 Heisenberg antiferromagnet on the triangular lattice

Authors: Qian Li, Hong Li, Jize Zhao, Hong-Gang Luo, Z. Y. Xie

Abstract: After decades of debate, now there is a rough consensus that at zero temperature the spin-$1/2$ Heisenberg antiferromagnet on the triangular lattice is three-sublattice $120^\circ$ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work… ▽ More After decades of debate, now there is a rough consensus that at zero temperature the spin-$1/2$ Heisenberg antiferromagnet on the triangular lattice is three-sublattice $120^\circ$ magnetically ordered, in contrast to a quantum spin liquid as originally proposed. However, there remains considerable discrepancy in the magnetization reported among various methods. To resolve this issue, in this work we revisit this model by the tensor-network state algorithm. The ground-state energy per bond $E_b$ and magnetization per spin $M_0$ in the thermodynamic limit are obtained with high precision. The former is estimated to be $E_b = -0.18334(10)$. This value agrees well with that from the series expansion. The three-sublattice magnetic order is firmly confirmed and the magnetization is determined as $M_0 = 0.161(5)$. It is about $32\%$ of its classical value and slightly below the lower bound from the series expansion. In comparison with the best estimated value by Monte Carlo and density-matrix renormalization group, our result is about $20\%$ smaller. This magnetic order is consistent with further analysis of the three-body correlation. Our work thus provides new benchmark results for this prototypical model. △ Less

Submitted 24 February, 2022; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: 7 pages, 6 figures, 2 tables, new data and references added

Journal ref: Phys. Rev. B 105, 184418 (2022)

Critical properties of the two-dimensional $q$-state clock model

Authors: Zi-Qian Li, Li-Ping Yang, Z. Y. Xie, Hong-Hao Tu, Hai-Jun Liao, T. Xiang

Abstract: We perform the state-of-the-art tensor network simulations directly in the thermodynamic limit to clarify the critical properties of the $q$-state clock model on the square lattice. We determine accurately the two phase transition temperatures through the singularity of the classical analog of the entanglement entropy, and provide extensive numerical evidences to show that both transitions are of… ▽ More We perform the state-of-the-art tensor network simulations directly in the thermodynamic limit to clarify the critical properties of the $q$-state clock model on the square lattice. We determine accurately the two phase transition temperatures through the singularity of the classical analog of the entanglement entropy, and provide extensive numerical evidences to show that both transitions are of the Berezinskii-Kosterlitz-Thouless (BKT) type for $q\ge 5$ and that the low-energy physics of this model is well described by the $\mathbb{Z}_q$-deformed sine-Gordon theory. We also determine the characteristic conformal parameters, especially the compactification radius, that govern the critical properties of the intermediate BKT phase. △ Less

Submitted 25 June, 2020; v1 submitted 24 December, 2019; originally announced December 2019.

Comments: 5 pages, 5 figures plus supplemental material; accepted for publication in Phys. Rev. E (Rapid Com.), minor typo corrections

Journal ref: Phys. Rev. E 101, 060105 (2020)

Automatic Differentiation for Second Renormalization of Tensor Networks

Authors: Bin-Bin Chen, Yuan Gao, Yi-Bin Guo, Yuzhi Liu, Hui-Hai Zhao, Hai-Jun Liao, Lei Wang, Tao Xiang, Wei Li, Z. Y. Xie

Abstract: Tensor renormalization group (TRG) constitutes an important methodology for accurate simulations of strongly correlated lattice models. Facilitated by the automatic differentiation technique widely used in deep learning, we propose a uniform framework of differentiable TRG ($\partial$TRG) that can be applied to improve various TRG methods, in an automatic fashion. Essentially, $\partial$TRG system… ▽ More Tensor renormalization group (TRG) constitutes an important methodology for accurate simulations of strongly correlated lattice models. Facilitated by the automatic differentiation technique widely used in deep learning, we propose a uniform framework of differentiable TRG ($\partial$TRG) that can be applied to improve various TRG methods, in an automatic fashion. Essentially, $\partial$TRG systematically extends the concept of second renormalization [PRL 103, 160601 (2009)] where the tensor environment is computed recursively in the backward iteration, in the sense that given the forward process of TRG, $\partial$TRG automatically finds the gradient through backpropagation, with which one can deeply "train" the tensor networks. We benchmark $\partial$TRG in solving the square-lattice Ising model, and demonstrate its power by simulating one- and two-dimensional quantum systems at finite temperature. The deep optimization as well as GPU acceleration renders $\partial$TRG manybody simulations with high efficiency and accuracy. △ Less

Submitted 4 July, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: 13 pages, 9 figures

Journal ref: Phys. Rev. B 101, 220409 (2020)

Cross derivative: a universal and efficient method for phase transitions in classical spin models

Authors: Y. Chen, K. Ji, Z. Y. Xie, J. F. Yu

Abstract: With an auxiliary weak external magnetic field, we reexamine the fundamental thermodynamic function, Gibbs free energy F(T, h), to study the phase transitions in the classical spin lattice models. A cross derivative, i.e. the second-order partial derivative of F(T, h) with respect to both temperature and field, is calculated to precisely locate the critical temperature, which also reveals the natu… ▽ More With an auxiliary weak external magnetic field, we reexamine the fundamental thermodynamic function, Gibbs free energy F(T, h), to study the phase transitions in the classical spin lattice models. A cross derivative, i.e. the second-order partial derivative of F(T, h) with respect to both temperature and field, is calculated to precisely locate the critical temperature, which also reveals the nature of a transition. The strategy is efficient and universal, as exemplified by the 5-state clock model, 2-dimensional (2D) and 3D Ising models, and the XY model, no matter a transition is trivial or exotic with complex excitations. More importantly, other conjugate pairs could also be integrated into a similar cross derivative if necessary, which would greatly enrich our vision and means to investigate phase transitions both theoretically and experimentally. △ Less

Submitted 14 October, 2019; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. B 101, 165123 (2020)

Compressing deep neural networks by matrix product operators

Authors: Ze-Feng Gao, Song Cheng, Rong-Qiang He, Z. Y. Xie, Hui-Hai Zhao, Zhong-Yi Lu, Tao Xiang

Abstract: A deep neural network is a parametrization of a multilayer mapping of signals in terms of many alternatively arranged linear and nonlinear transformations. The linear transformations, which are generally used in the fully connected as well as convolutional layers, contain most of the variational parameters that are trained and stored. Compressing a deep neural network to reduce its number of varia… ▽ More A deep neural network is a parametrization of a multilayer mapping of signals in terms of many alternatively arranged linear and nonlinear transformations. The linear transformations, which are generally used in the fully connected as well as convolutional layers, contain most of the variational parameters that are trained and stored. Compressing a deep neural network to reduce its number of variational parameters but not its prediction power is an important but challenging problem toward the establishment of an optimized scheme in training efficiently these parameters and in lowering the risk of overfitting. Here we show that this problem can be effectively solved by representing linear transformations with matrix product operators (MPOs), which is a tensor network originally proposed in physics to characterize the short-range entanglement in one-dimensional quantum states. We have tested this approach in five typical neural networks, including FC2, LeNet-5, VGG, ResNet, and DenseNet on two widely used data sets, namely, MNIST and CIFAR-10, and found that this MPO representation indeed sets up a faithful and efficient mapping between input and output signals, which can keep or even improve the prediction accuracy with a dramatically reduced number of parameters. Our method greatly simplifies the representations in deep learning, and opens a possible route toward establishing a framework of modern neural networks which might be simpler and cheaper, but more efficient. △ Less

Submitted 9 June, 2020; v1 submitted 11 April, 2019; originally announced April 2019.

Comments: 8+9 pages, 3+7 figures, 2+11 tables

Journal ref: Phys. Rev. Research 2, 023300 (2020)

Phase transitions of the five-state clock model on the square lattice

Authors: Y. Chen, Z. Y. Xie, J. F. Yu

Abstract: Using the tensor renormalization group method based on the higher-order singular value decomposition,we have studied the phase transitions of the five-state clock model on the square lattice.The temperature dependence of the specific heat indicates the system has two phase transitions, as verified clearly in the correlation function. By investigating the magnetic susceptibility, we can only obtain… ▽ More Using the tensor renormalization group method based on the higher-order singular value decomposition,we have studied the phase transitions of the five-state clock model on the square lattice.The temperature dependence of the specific heat indicates the system has two phase transitions, as verified clearly in the correlation function. By investigating the magnetic susceptibility, we can only obtain the upper critical temperature as Tc2 = 0.9565(7). From the fixed-point tensor, we locate the transition temperatures at Tc1 = 0.9029(1) and Tc2 = 0.9520(1), consistent with the MC and the DMRG results. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 5 pages, 7 figures

Journal ref: Chin. Phys. B 27, 080503 (2018)

Reorthonormalization of Chebyshev matrix product states for dynamical correlation functions

Authors: H. D. Xie, R. Z. Huang, X. J. Han, X. Yan, H. H. Zhao, Z. Y. Xie, H. J. Liao, T. Xiang

Abstract: The Chebyshev expansion offers a numerically efficient and easy-implement algorithm for evaluating dynamic correlation functions using matrix product states (MPS). In this approach, each recursively generated Chebyshev vector is approximately represented by an MPS. However, the recurrence relations of Chebyshev polynomials are broken by the approximation, leading to an error which is accumulated w… ▽ More The Chebyshev expansion offers a numerically efficient and easy-implement algorithm for evaluating dynamic correlation functions using matrix product states (MPS). In this approach, each recursively generated Chebyshev vector is approximately represented by an MPS. However, the recurrence relations of Chebyshev polynomials are broken by the approximation, leading to an error which is accumulated with the increase of the order of expansion. Here we propose a reorthonormalization approach to remove this error introduced in the loss of orthogonality of the Chebyshev polynomials. Our approach, as illustrated by comparison with the exact results for the one-dimensional XY and Heisenberg models, improves significantly the accuracy in the calculation of dynamical correlation functions. △ Less

Submitted 25 November, 2017; originally announced November 2017.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. B 97, 075111 (2018)

Optimized contraction scheme for tensor-network states

Authors: Z. Y. Xie, H. J. Liao, R. Z. Huang, H. D. Xie, J. Chen, Z. Y. Liu, T. Xiang

Abstract: In the tensor-network framework, the expectation values of two-dimensional quantum states are evaluated by contracting a double-layer tensor network constructed from initial and final tensor-network states. The computational cost of carrying out this contraction is generally very high, which limits the largest bond dimension of tensor-network states that can be accurately studied to a relatively s… ▽ More In the tensor-network framework, the expectation values of two-dimensional quantum states are evaluated by contracting a double-layer tensor network constructed from initial and final tensor-network states. The computational cost of carrying out this contraction is generally very high, which limits the largest bond dimension of tensor-network states that can be accurately studied to a relatively small value. We propose an optimized contraction scheme to solve this problem by mapping the double-layer tensor network onto an intersected single-layer tensor network. This reduces greatly the bond dimensions of local tensors to be contracted and improves dramatically the efficiency and accuracy of the evaluation of expectation values of tensor-network states. It almost doubles the largest bond dimension of tensor-network states whose physical properties can be efficiently and reliably calculated, and it extends significantly the application scope of tensor-network methods. △ Less

Submitted 21 July, 2017; v1 submitted 23 May, 2017; originally announced May 2017.

Comments: 7 pages, 6 figures, accepted version for Phys. Rev. B. The title is changed for easier understanding

Journal ref: Phys. Rev. B 96, 045128 (2017)

Gapless spin-liquid ground state in the $S = 1/2$ kagome antiferromagnet

Authors: H. J. Liao, Z. Y. Xie, J. Chen, Z. Y. Liu, H. D. Xie, R. Z. Huang, B. Normand, T. Xiang

Abstract: Frustrated quantum magnetism has moved to the forefront of physics research, posing fundamental questions concerning quantum disordered states, entanglement, topology and the nature of the quantum wavefunction. The defining problem in the field is one of the simplest, the ground state of the nearest-neighbour $S = 1/2$ antiferromagnetic Heisenberg model on the kagome lattice, but has defied all th… ▽ More Frustrated quantum magnetism has moved to the forefront of physics research, posing fundamental questions concerning quantum disordered states, entanglement, topology and the nature of the quantum wavefunction. The defining problem in the field is one of the simplest, the ground state of the nearest-neighbour $S = 1/2$ antiferromagnetic Heisenberg model on the kagome lattice, but has defied all theoretical and numerical methods employed to date. We apply the formalism of tensor-network states (TNS), specifically the method of projected entangled simplex states (PESS), whose combination of a correct accounting for multipartite entanglement and infinite system size provides qualitatively new insight. By studying the ground-state energy, the staggered magnetization we find at all finite tensor bond dimensions and the effects of a second-neighbour coupling, we demonstrate that the ground state is a gapless spin liquid. We discuss the comparison with other numerical studies and the physical interpretation of the gapless ground state. △ Less

Submitted 28 February, 2017; v1 submitted 15 October, 2016; originally announced October 2016.

Comments: 8 pages, 4 figures, revised version accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 118, 137202 (2017)

Charge dynamics of the antiferromagnetically ordered Mott insulator

Authors: X. J. Han, Y. Liu, Z. Y. Liu, X. Li, J. Chen, Z. Y. Xie, B. Normand, T. Xiang

Abstract: We study the charge dynamics of the half-filled Hubbard model on the square lattice at zero temperature. We employ a slave-fermion formulation in which the charge degrees of freedom are represented by fermionic holons and doublons and the bosonic spin degrees of freedom are assumed to be described by the antiferromagnetic Heisenberg model. The ground state has long-ranged magnetic (Néel) order and… ▽ More We study the charge dynamics of the half-filled Hubbard model on the square lattice at zero temperature. We employ a slave-fermion formulation in which the charge degrees of freedom are represented by fermionic holons and doublons and the bosonic spin degrees of freedom are assumed to be described by the antiferromagnetic Heisenberg model. The ground state has long-ranged magnetic (Néel) order and its Mott-insulating characteristics are the consequence of holon-doublon bound-state formation. Within this framework, we calculate the average double occupancy, the electron density of states, and the spectral function in the self-consistent Born approximation. The lower and upper Hubbard bands are reproduced, there is spectral-weight transfer into a coherent quasiparticle band at their lower edges, and the Mott gap, associated with holon-doublon binding, is renormalized due to the interactions of the holons and doublons with the magnons. We compute the zeros of the Green function at the chemical potential to discuss the Luttinger volume and show that the poles of the self-energy reflect the underlying quasiparticle dispersion with a spin-renormalized hopping parameter. Optical conductivity results relate the optical gap to the Mott gap of the insulating system. We conclude that a self-consistent treatment of the spin fluctuation effects on the charge degrees of freedom captures all the essential physics of the antiferromagnetic Mott-Hubbard insulator. △ Less

Submitted 14 December, 2016; v1 submitted 3 February, 2016; originally announced February 2016.

Journal ref: New J. Phys. 18, 103004 (2016)

Heisenberg antiferromagnet on the Husimi lattice

Authors: H. J. Liao, Z. Y. Xie, J. Chen, X. J. Han, H. D. Xie, B. Normand, T. Xiang

Abstract: We perform a systematic study of the antiferromagnetic Heisenberg model on the Husimi lattice using numerical tensor-network methods based on Projected Entangled Simplex States (PESS). The nature of the ground state varies strongly with the spin quantum number, $S$. For $S = 1/2$, it is an algebraic (gapless) quantum spin liquid. For $S = 1$, it is a gapped, non-magnetic state with spontaneous bre… ▽ More We perform a systematic study of the antiferromagnetic Heisenberg model on the Husimi lattice using numerical tensor-network methods based on Projected Entangled Simplex States (PESS). The nature of the ground state varies strongly with the spin quantum number, $S$. For $S = 1/2$, it is an algebraic (gapless) quantum spin liquid. For $S = 1$, it is a gapped, non-magnetic state with spontaneous breaking of triangle symmetry (a trimerized simplex-solid state). For $S = 2$, it is a simplex-solid state with a spin gap and no symmetry-breaking; both integer-spin simplex-solid states are characterized by specific degeneracies in the entanglement spectrum. For $S = 3/2$, and indeed for all spin values $S \ge 5/2$, the ground states have $120$-degree antiferromagnetic order. In a finite magnetic field, we find that, irrespective of the value of $S$, there is always a plateau in the magnetization at $m = 1/3$. △ Less

Submitted 7 March, 2016; v1 submitted 29 October, 2015; originally announced October 2015.

Comments: 18 pages, 25 figures;minor changes, refs updated

Journal ref: Phys. Rev. B 93, 075154 (2016)

The fine structure of the entanglement entropy in the classical XY model

Authors: Li-Ping Yang, Yuzhi Liu, Haiyuan Zou, Z. Y. Xie, Y. Meurice

Abstract: We compare two calculations of the particle density in the superfluid phase of the classical XY model with a chemical potential $μ$ in 1+1 dimensions.The first relies on exact blocking formulas from the Tensor Renormalization Group (TRG) formulation of the transfer matrix. The second is a worm algorithm. We show that the particle number distributions obtained with the two methods agree well. We us… ▽ More We compare two calculations of the particle density in the superfluid phase of the classical XY model with a chemical potential $μ$ in 1+1 dimensions.The first relies on exact blocking formulas from the Tensor Renormalization Group (TRG) formulation of the transfer matrix. The second is a worm algorithm. We show that the particle number distributions obtained with the two methods agree well. We use the TRG method to calculate the thermal entropy and the entanglement entropy. We describe the particle density, the two entropies and the topology of the world lines as we increase $μ$ to go across the superfluid phase between the first two Mott insulating phases. For a sufficiently large temporal size, this process reveals an interesting fine structure: the average particle number and the winding number of most of the world lines in the Euclidean time direction increase by one unit at a time. At each step, the thermal entropy develops a peak and the entanglement entropy increases until we reach half-filling and then decreases in a way that approximately mirror the ascent. This suggests an approximate fermionic picture. △ Less

Submitted 12 July, 2015; v1 submitted 6 July, 2015; originally announced July 2015.

Comments: 12 pages, 12 figures

Journal ref: Phys. Rev. E 93, 012138 (2016)

Partial long-range order in antiferromagnetic Potts models

Authors: M. P. Qin, Q. N. Chen, Z. Y. Xie, J. Chen, J. F. Yu, H. H. Zhao, B. Normand, T. Xiang

Abstract: The Potts model plays an essential role in classical statistical mechanics, illustrating many fundamental phenomena. One example is the existence of partially long-range-ordered states, in which some degrees of freedom remain disordered. This situation may arise from frustration of the interactions, but also from an irregular but unfrustrated lattice structure. We study partial long-range order in… ▽ More The Potts model plays an essential role in classical statistical mechanics, illustrating many fundamental phenomena. One example is the existence of partially long-range-ordered states, in which some degrees of freedom remain disordered. This situation may arise from frustration of the interactions, but also from an irregular but unfrustrated lattice structure. We study partial long-range order in a range of antiferromagnetic $q$-state Potts models on different two-dimensional lattices and for all relevant values of $q$. We exploit the power of tensor-based numerical methods to evaluate the partition function of these models and hence to extract the key thermodynamic properties -- entropy, specific heat, magnetization, and susceptibility -- giving deep insight into the phase transitions and ordered states of each system. Our calculations reveal a range of phenomena related to partial ordering, including different types of entropy-driven phase transition, the role of lattice irregularity, very large values of the critical $q_c$, and double phase transitions. △ Less

Submitted 23 October, 2014; originally announced October 2014.

Comments: 23 pages, 27 figures

Journal ref: Phys. Rev. B 90, 144424(2014)

Phase Transitions of Ferromagnetic Potts Models on the Simple Cubic Lattice

Authors: S. Wang, Z. Y. Xie, J. Chen, B. Normand, T. Xiang

Abstract: We investigate the 2- and 3-state ferromagnetic Potts models on the simple cubic lattice using the tensor renormalization group method with higher-order singular value decomposition (HOTRG). HOTRG works in the thermodynamic limit, where we use the $Z_q$ symmetry of the model, combined with a new measure for detecting the transition, to improve the accuracy of the critical point for the 2-state mod… ▽ More We investigate the 2- and 3-state ferromagnetic Potts models on the simple cubic lattice using the tensor renormalization group method with higher-order singular value decomposition (HOTRG). HOTRG works in the thermodynamic limit, where we use the $Z_q$ symmetry of the model, combined with a new measure for detecting the transition, to improve the accuracy of the critical point for the 2-state model by two orders of magnitude, obtaining $T_c = 4.51152469(1)$. The 3-state model is far more complex, and we improve the overall understanding of this case by calculating its thermodynamic quantities with high accuracy. Our results verify the first-order nature of the phase transition and the HOTRG transition temperature benchmarks the most recent Monte Carlo result. △ Less

Submitted 6 May, 2014; originally announced May 2014.

Comments: 6 pages, 5 figures, 2 tables

Journal ref: Chin. Phys. Lett. 31, 070503 (2014)

Progress towards quantum simulating the classical O(2) model

Authors: Haiyuan Zou, Yuzhi Liu, Chen-Yen Lai, J. Unmuth-Yockey, Li-Ping Yang, A. Bazavov, Z. Y. Xie, T. Xiang, S. Chandrasekharan, S. -W. Tsai, Y. Meurice

Abstract: We connect explicitly the classical $O(2)$ model in 1+1 dimensions, a model sharing important features with $U(1)$ lattice gauge theory, to physical models potentially implementable on optical lattices and evolving at physical time. Using the tensor renormalization group formulation, we take the time continuum limit and check that finite dimensional projections used in recent proposals for quantum… ▽ More We connect explicitly the classical $O(2)$ model in 1+1 dimensions, a model sharing important features with $U(1)$ lattice gauge theory, to physical models potentially implementable on optical lattices and evolving at physical time. Using the tensor renormalization group formulation, we take the time continuum limit and check that finite dimensional projections used in recent proposals for quantum simulators provide controllable approximations of the original model. We propose two-species Bose-Hubbard models corresponding to these finite dimensional projections at strong coupling and discuss their possible implementations on optical lattices using a $^{87}$Rb and $^{41}$K Bose-Bose mixture. △ Less

Submitted 25 November, 2014; v1 submitted 20 March, 2014; originally announced March 2014.

Comments: 7 pages, 6 figures, uses revtex, new material and one author added, as to appear in Phys. Rev. A

Journal ref: Phys. Rev. A 90, 063603 (2014)

Controlling sign problems in spin models using tensor renormalization

Authors: Alan Denbleyker, Yuzhi Liu, Y. Meurice, M. P. Qin, T. Xiang, Z. Y. Xie, J. F. Yu, Haiyuan Zou

Abstract: We consider the sign problem for classical spin models at complex $β=1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$β$ than the reweighting Monte Carlo method. For the Ising model with complex $β$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined prec… ▽ More We consider the sign problem for classical spin models at complex $β=1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$β$ than the reweighting Monte Carlo method. For the Ising model with complex $β$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined precisely with the TRG method. We check the convergence of the TRG method for the O(2) model on $L\times L$ lattices when the number of states $D_s$ increases. We show that the finite size scaling of the calculated Fisher zeros agrees very well with the Kosterlitz-Thouless transition assumption and predict the locations for larger volume. The location of these zeros agree with Monte Carlo reweighting calculation for small volume. The application of the method for the O(2) model with a chemical potential is briefly discussed. △ Less

Submitted 25 September, 2013; originally announced September 2013.

Comments: 8 pages, 9 figures

Report number: FERMILAB-PUB-13-394-T

Journal ref: Phys. Rev. D 89, 016008 (2014)

Tensor renormalization group study of classical XY model on the square lattice

Authors: J. F. Yu, Z. Y. Xie, Y. Meurice, Yuzhi Liu, A. Denbleyker, Haiyuan Zou, M. P. Qin, J. Chen, T. Xiang

Abstract: Using the tensor renormalization group method based on the higher-order singular value decom- position, we have studied the thermodynamic properties of the continuous XY model on the square lattice. The temperature dependence of the free energy, the internal energy and the specific heat agree with the Monte Carlo calculations. From the field dependence of the magnetic susceptibility, we find the K… ▽ More Using the tensor renormalization group method based on the higher-order singular value decom- position, we have studied the thermodynamic properties of the continuous XY model on the square lattice. The temperature dependence of the free energy, the internal energy and the specific heat agree with the Monte Carlo calculations. From the field dependence of the magnetic susceptibility, we find the Kosterlitz-Thouless transition temperature to be 0.8921 \pm 0.0019, consistent with the Monte Carlo as well as the high temperature series expansion results. At the transition temperature, the critical exponent δis estimated as 14.5, close to the analytic value by Kosterlitz. △ Less

Submitted 19 September, 2013; originally announced September 2013.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. E 89, 013308 (2014)

Exact blocking formulas for spin and gauge models

Authors: Yuzhi Liu, Y. Meurice, M. P. Qin, J. Unmuth-Yockey, T. Xiang, Z. Y. Xie, J. F. Yu, Haiyuan Zou

Abstract: Using the example of the two-dimensional (2D) Ising model, we show that in contrast to what can be done in configuration space, the tensor renormalization group (TRG) formulation allows one to write exact, compact, and manifestly local blocking formulas and exact coarse grained expressions for the partition function. We argue that similar results should hold for most models studied by lattice gaug… ▽ More Using the example of the two-dimensional (2D) Ising model, we show that in contrast to what can be done in configuration space, the tensor renormalization group (TRG) formulation allows one to write exact, compact, and manifestly local blocking formulas and exact coarse grained expressions for the partition function. We argue that similar results should hold for most models studied by lattice gauge theorists. We provide exact blocking formulas for several 2D spin models (the O(2) and O(3) sigma models and the SU(2) principal chiral model) and for the 3D gauge theories with groups Z_2, U(1) and SU(2). We briefly discuss generalizations to other groups, higher dimensions and practical implementations. △ Less

Submitted 24 July, 2013; originally announced July 2013.

Comments: 12 pages, 10 figures

Report number: FERMILAB-PUB-13-280-T}

Journal ref: Phys. Rev. D 88, 056005 (2013)

Tensor renormalization of quantum many-body systems using projected entangled simplex states

Authors: Z. Y. Xie, J. Chen, J. F. Yu, X. Kong, B. Normand, T. Xiang

Abstract: We propose a new class of tensor-network states, which we name projected entangled simplex states (PESS), for studying the ground-state properties of quantum lattice models. These states extend the pair-correlation basis of projected entangled pair states (PEPS) to a simplex. PESS are an exact representation of the simplex solid states and provide an efficient trial wave function that satisfies th… ▽ More We propose a new class of tensor-network states, which we name projected entangled simplex states (PESS), for studying the ground-state properties of quantum lattice models. These states extend the pair-correlation basis of projected entangled pair states (PEPS) to a simplex. PESS are an exact representation of the simplex solid states and provide an efficient trial wave function that satisfies the area law of entanglement entropy. We introduce a simple update method for evaluating the PESS wave function based on imaginary-time evolution and the higher-order singular-value decomposition of tensors. By applying this method to the spin-1/2 antiferromagnetic Heisenberg model on the kagome lattice, we obtain accurate and systematic results for the ground-state energy, which approach the lowest upper bounds yet estimated for this quantity. △ Less

Submitted 1 January, 2014; v1 submitted 22 July, 2013; originally announced July 2013.

Comments: 12 pages, 8 figures

Journal ref: Phys. Rev. X 4, 011025 (2014)

Partial order in Potts models on the generalized decorated square lattice

Authors: M. P. Qin, J. Chen, Q. N. Chen, Z. Y. Xie, X. Kong, H. H. Zhao, B. Normand, T. Xiang

Abstract: We explore the Potts model on the generalized decorated square lattice, with both nearest (J1) and next-neighbor (J2) interactions. Using the tensor renormalization-group method augmented by higher-order singular value decompositions, we calculate the spontaneous magnetization of the Potts model with q = 2, 3, and 4. The results for q = 2 allow us to benchmark our numerics using the exact solution… ▽ More We explore the Potts model on the generalized decorated square lattice, with both nearest (J1) and next-neighbor (J2) interactions. Using the tensor renormalization-group method augmented by higher-order singular value decompositions, we calculate the spontaneous magnetization of the Potts model with q = 2, 3, and 4. The results for q = 2 allow us to benchmark our numerics using the exact solution. For q = 3, we find a highly degenerate ground state with partial order on a single sublattice, but with vanishing entropy per site, and we obtain the phase diagram as a function of the ratio J2=J1. There is no finite-temperature transition for the q = 4 case when J1 = J2, but the magnetic susceptibility diverges as the temperature goes to zero, showing that the model is critical at T = 0. △ Less

Submitted 26 April, 2013; originally announced April 2013.

Comments: 5 pages, 6 figures

Journal ref: Chin. Phys. Lett. 30, 076402 (2013)

Coarse-graining renormalization by higher-order singular value decomposition

Authors: Z. Y. Xie, J. Chen, M. P. Qin, J. W. Zhu, L. P. Yang, T. Xiang

Abstract: We propose a novel coarse graining tensor renormalization group method based on the higher-order singular value decomposition. This method provides an accurate but low computational cost technique for studying both classical and quantum lattice models in two- or three-dimensions. We have demonstrated this method using the Ising model on the square and cubic lattices. By keeping up to 16 bond basis… ▽ More We propose a novel coarse graining tensor renormalization group method based on the higher-order singular value decomposition. This method provides an accurate but low computational cost technique for studying both classical and quantum lattice models in two- or three-dimensions. We have demonstrated this method using the Ising model on the square and cubic lattices. By keeping up to 16 bond basis states, we obtain by far the most accurate numerical renormalization group results for the 3D Ising model. We have also applied the method to study the ground state as well as finite temperature properties for the two-dimensional quantum transverse Ising model and obtain the results which are consistent with published data. △ Less

Submitted 7 August, 2012; v1 submitted 5 January, 2012; originally announced January 2012.

Comments: 11 pages, 16 figures

Journal ref: Phys. Rev. B 86, 045139 (2012)

Translation invariant tensor product states in a finite lattice system

Authors: J. W. Cai, Q. N. Chen, H. H. Zhao, Z. Y. Xie, M. P. Qin, Z. C. Wei, T. Xiang

Abstract: We show that the matrix (or more generally tensor) product states in a finite translation invariant system can be accurately constructed from the same set of local matrices (or tensors) that are determined from an infinite lattice system in one or higher dimensions. This provides an efficient approach for studying translation invariant tensor product states in finite lattice systems. Two methods a… ▽ More We show that the matrix (or more generally tensor) product states in a finite translation invariant system can be accurately constructed from the same set of local matrices (or tensors) that are determined from an infinite lattice system in one or higher dimensions. This provides an efficient approach for studying translation invariant tensor product states in finite lattice systems. Two methods are introduced to determine these size-independent local tensors. △ Less

Submitted 9 November, 2010; originally announced November 2010.

Journal ref: Chin. Phys. Lett. 27, 076402 (2010)

Renormalization of tensor-network states

Authors: H. H. Zhao, Z. Y. Xie, Q. N. Chen, Z. C. Wei, J. W. Cai, T. Xiang

Abstract: We have discussed the tensor-network representation of classical statistical or interacting quantum lattice models, and given a comprehensive introduction to the numerical methods we recently proposed for studying the tensor-network states/models in two dimensions. A second renormalization scheme is introduced to take into account the environment contribution in the calculation of the partition… ▽ More We have discussed the tensor-network representation of classical statistical or interacting quantum lattice models, and given a comprehensive introduction to the numerical methods we recently proposed for studying the tensor-network states/models in two dimensions. A second renormalization scheme is introduced to take into account the environment contribution in the calculation of the partition function of classical tensor network models or the expectation values of quantum tensor network states. It improves significantly the accuracy of the coarse grained tensor renormalization group method. In the study of the quantum tensor-network states, we point out that the renormalization effect of the environment can be efficiently and accurately described by the bond vector. This, combined with the imaginary time evolution of the wavefunction, provides an accurate projection method to determine the tensor-network wavfunction. It reduces significantly the truncation error and enable a tensor-network state with a large bond dimension, which is difficult to be accessed by other methods, to be accurately determined. △ Less

Submitted 26 April, 2010; v1 submitted 6 February, 2010; originally announced February 2010.

Comments: 18 pages 23 figures, minor changes, references added

Journal ref: Physical Review B 81, 174411 (2010)

Second Renormalization of Tensor-Network States

Authors: Z. Y. Xie, H. C. Jiang, Q. N. Chen, Z. Y. Weng, T. Xiang

Abstract: We propose a second renormalization group method to handle the tensor-network states or models. This method reduces dramatically the truncation error of the tensor renormalization group. It allows physical quantities of classical tensor-network models or tensor-network ground states of quantum systems to be accurately and efficiently determined. We propose a second renormalization group method to handle the tensor-network states or models. This method reduces dramatically the truncation error of the tensor renormalization group. It allows physical quantities of classical tensor-network models or tensor-network ground states of quantum systems to be accurately and efficiently determined. △ Less

Submitted 14 October, 2009; v1 submitted 1 September, 2008; originally announced September 2008.

Comments: 5 figures, 4 pages

Journal ref: Phys.Rev.Lett.103, 160601 (2009)