两腿XXZ自旋梯子模型的量子相变与量子临界现象

对于无限大尺寸两腿自旋1/2的XXZ自旋梯子模型,通过运用基于随机行走的张量网络(TN)算法数值模拟出基态波函数,首次尝试研究自旋梯子模型的约化保真度、普适序参量、纠缠熵等物理观测量,并系统研究基态保真度的三维挤点与二维分叉、约化保真度的分叉、局域序参量、普适序参量、纠缠熵和量子相变之间存在的关联关系.基于张量网络表示的算法在任意随机选择初始状态时,可以得到两腿XXZ量子自旋梯子系统简并的对称破缺基态波函数,该基态波函数是由于Z2对称破缺引起的.本文期望所提供的方法可为进一步研究凝聚态物质中热力学极限下的强关联电子量子晶格自旋梯子系统的量子相变和量子临界现象提供一种更有效的强大的工具.     

Ground-state phase diagram of the dimerizedspin-1/2 two-leg ladder

Dimerized spin-1/2 ladders exhibit a variety of phase structures, which depend on the intra-chain and inter-chain spin exchange energies as well as on the dimerization pattern of the ladder. Using the density matrix renormalization group (DMRG) algorithm, we study critical properties of the bond-alternating two-leg Heisenberg spin ladder with diagonal interaction J_×. Two types of spin systems, staggered dimerized antiferromagnetic ladder and columnar dimerized ferro-antiferromagnetic couplings ladder, are investigated. To clarify the phase transition behaviors, we simultaneously analyze the string order parameter (SOP), the twisted order parameter (TOP), as well as a measurement of the quantum information analysis. Based on measuring this different observables, we establish the phase diagram accurately and give the fitting functions of the phase boundaries. In addition, the phase transition of cross-coupled spin ladder (in the absence of intrinsic dimerization) is also discussed.     

0- and 2/3-magnetization plateaus in three-leg antiferromagnetic Heisenberg spin-1/2 ladders with leg-dimerization

Magnetic properties of three-leg antiferromagnetic Heisenberg spin-1/2 ladders with different dimerization patterns have been studied using the bond mean-field theory. Our results show that rung-columnar ladders are thermodynamically stable states for large rung-to-leg coupling ratios. Magnetization curves of leg-columnar and leg-staggered ladders always exhibit 0- and 2/3-magnetization plateaus, which do not appear in rung-columnar and rung-staggered ladders. In leg-dimerized ladders, the formation of spin dimers in the three legs results in the appearance of the 0- and 2/3-magnetization plateaus. Spin configuration in the 2/3-magnetization plateau can be understood from the mean-field bond parameters.     

一维自旋1键交替XXZ链中的量子纠缠和临界指数

利用基于张量网络表示的矩阵乘积态算法以及无限虚时间演化块抽取方法,本文研究了一维无限格点自旋1的键交替反铁磁XXZ海森伯模型中的量子相变.分别计算了系统的von Neumann熵、单位格点保真度和序参量,从而得到了系统随键交替强度的变化从拓扑有序Néel相到局域有序二聚化相的量子相变点.我们用矩阵乘积态方法拟合出了相变的中心荷c?0.5,表明此相变属于二维经典的Ising普适类.另外,通过对拓扑Néel序的数值拟合,我们得到了相变点处的特征临界指数β′=0.236和γ′=0.838.     

Bipartite entanglement in the spin-1/2 Ising-Heisenberg planar lattice constituted of identical trigonal bipyramidal plaquettes

The bipartite entanglement is rigorously examined in the spin-1/2 Ising-Heisenberg planar lattice composed of identical inter-connected bipyramidal plaquettes at zero and finite temperatures using the quantity called concurrence. It is shown that the Heisenberg spins of the same plaquette are twice stronger entangled in the two-fold degenerate quantum ground state than in the macroscopically degenerate quantum chiral one. The bipartite entanglement with chiral features completely disappears below or exactly at the critical temperature of the model, while that with no chirality may survive even above the critical temperature of the model. Non-monotonous temperature variations of the concurrence clearly evidence the activation of the entangled Heisenberg states also above classical ground state as well as their re-appearance above the critical temperature of the model.     

Optimal control of fast and high-fidelity quantum state transfer in spin-1/2 chains

Spin chains are promising candidates for quantum communication and computation. Using quantum optimal control (OC) theory based on the Krotov method, we present a protocol to perform quantum state transfer with fast and high fidelity by only manipulating the boundary spins in a quantum spin-1/2 chain. The achieved speed is about one order of magnitude faster than that is possible in the Lyapunov control case for comparable fidelities. Additionally, it has a fundamental limit for OC beyond which optimization is not possible. The controls are exerted only on the couplings between the boundary spins and their neighbors, so that the scheme has good scalability. We also demonstrate that the resulting OC scheme is robust against disorder in the chain.     

Genuine tripartite entanglement and quantum phase transition

A new simplified formula is presented to characterize genuine tripartite entanglement of （2 2 n）-dimensional quantum pure states. The formula turns out equivalent to that given in （Quant. Inf. Comp. 7（7） 584 （2007））, hence it also shows that the genuine tripartite entanglement can be described only on the basis of the local （2 2）-dimensional reduced density matrix. In particular, the two exactly solvable models of spin system studied by Yang （Phys. Rev. A 71 030302（R） （2005）） are reconsidered by employing the formula. The results show that a discontinuity in the first derivative of the formula or in the formula itself of the ground state just corresponds to the existence of quantum phase transition, which is obviously different from the concurrence.     

Identification of intrinsic gapped behavior in spin-1/2 ladder with staggered dimerization

The zero- and low-temperature behaviors of spin-1/2 two-leg ladder with staggered dimerization are investigated by the Green’s function theory. At zero temperature, the ground state phase diagram is explored, wherein the leg-dimer and rung-singlet phases are revealed, which reflect two different intrinsic gapped behaviors. The former is attributed to the bond alternation along the legs, while the latter is due to the strong rung coupling. It is found that the quantum phase transition from one to another is of the first order, which can be clearly signaled by the rung entanglement entropy. At finite temperatures, the temperature dependence of thermodynamic quantities such as the magnetic susceptibility, specific heat, thermal Drude weight and rung entanglement entropy are calculated to characterize the corresponding quantum phases. It is shown that the magnetic behaviors clearly manifest a typical antiferromagnetism at low temperature, which is in accordance with the experimental results. It is also found that the intrinsic gapped low-lying excitations are responsible for the observed thermodynamic behaviors.     

Bipartite entanglement in spin-1/2 Heisenberg model

The bipartite entanglement of the two- and three-spin Heisenberg model was investigated by using the concept of negativity. It is found that for the ground-state entanglement of the two-spin model, the negativity always decreases as B increases if Δ<γ－1, and it may keep a steady value of 0.5 in the region of B2－γ²]^1/2 if Δ>γ－1, while for that of the three-spin model, the negativity exhibits square wave structures if γ=0 or Δ=0. For thermal states, there are two areas showing entanglement, namely, the main region and the sub-region. The main region exists only when Δ>Δ_c (Δ_c=γ－1 and (γ²－1)/2 for the 2- and 3-spin model respectively) and extends in terms of B and T as Δ increases, while the sub-region survives only when γ≠0 and shrinks in terms of B and T as Δ increases.     

Ground-state entanglement in a three-spin transverse Ising model with energy current

The ground-state entanglement associated with a three-spin transverse Ising model is studied. By introducing an energy current into the system, a quantum phase transition to energy-current phase may be presented with the variation of external magnetic field; and the ground-state entanglement varies suddenly at the critical point of quantum phase transition. In our model, the introduction of energy current makes the entanglement between any two qubits become maximally robust.     

Thermodynamic properties of the spin-1/2 antiferromagnetic ladder under magnetic field

Specific heat (C_V) measurements in the spin-1/2 Cu₂(C₂H₁₂N₂)₂Cl₄ system under a magnetic field up to H =8.25 T are reported and compared to the results of numerical calculations based on the 2-leg antiferromagnetic Heisenberg ladder. While the temperature dependences of both the susceptibility and the low-field specific heat are accurately reproduced by this model, deviations are observed above the critical field H_C1 at which the spin gap closes. In this Quantum High Field phase, the contribution of the low-energy quantum fluctuations are stronger than in the Heisenberg ladder model. We argue that this enhancement can be attributed to dynamical lattice fluctuations. Finally, we show that such a Heisenberg ladder, for H > H C1, is unstable, when coupled to the 3D lattice, against a lattice distortion. These results provide an alternative explanation for the observed low temperature ( K-0.8 K) phase (previously interpreted as a 3D magnetic ordering) as a new type of incommensurate gapped state. Received: 23 July 1998 / Accepted: 24 August 1998     

交错Dzyaloshinskii-Moriya相互作用对反铁磁Heisenberg链纠缠的影响

主要研究了具有交错Dzyaloshinskii-Moriya(DM)相互作用的反铁磁Heisenberg链的纠缠.基于 density-matrix renormalization group(DMRG)的数值计算表明,交错DM相互作用消除了系统在外磁场H=2处的二级量子相变,从而量子纠缠反常行为也随之消失;同时纠缠范围的发散也被消除,意味着该模型因子化点的消失.交错DM相互作用导致系统在任意强场下也不会达到铁磁饱和状态,从而保持着自旋纠缠.交错DM相互作用有利于通过外场调控纠缠程度和纠缠范 关键词： Dzyaloshinskii-Moriya相互作用 量子纠缠 量子相变 纠缠范围     

Bipartite entanglement in spin-1/2Heisenberg model

The bipartite entanglement of the two-and three-spin Heisenberg model was investigated by using the concept of negativity.It is found that for the ground-state entanglement of the two-spin model,the negativity always decreases as B increases if A Δ＜y-1,and it may keep a steady value of 0.5in the region of B＜J[(Δ+1)2-y2]1/2if Δ＞y-1,while for that of the three-spin model,the negativity exhibits square wave structures if y=0 or Δ=0.For thermal states,there are two areas showing entanglement,namely,the main region and the sub-region.The main region exists only when Δ＞Δc(Δc1=and(y2-1)/2for the 2-and 3-spin model respectively)and extends in terms of B and T as Δ increases,while the sub-region survives only when y≠0 and shrinks in terms of B and T as Δ increases.     

二维3态Potts量子系统的保真度与序参量

二维无限正方格子上的量子3态Potts模型是发生一级相变还是二级相变?通过运用无限纠缠投影对态(i PEPS)算法,在进行数值模拟时任意选取初态,能得到二维无限正方格子上的3态Potts模型的三个不同的简并基态波函数,这些简并的情况是由自发对称性破缺引起的.首先,揭示了在二维系统中自发对称性破缺引起的相变可以运用单点基态保真度的分叉来研究,也反映了在二维系统中约化保真度同样有一个分叉行为;再者,还提出了二维系统的普适序参量以及多分量的复数局域序参量的行为来尝试研究二维3态Potts模型,共同确定系统发生的量子相变的临界点及其类型.即基于i PEPS算法,从单点基态保真度、约化保真度、普适序参量以及局域序参量的角度,来研究3态Potts模型的量子相变,其为一级相变.     

二维３态Potts量子系统的保真度与序参量

二维无限正方格子上的量子3态Potts模型是发生一级相变还是二级相变?通过运用无限纠缠投影对态 (iPEPS) 算法,在进行数值模拟时任意选取初态,能得到二维无限正方格子上的3态Potts模型的三个不同的简并基态波函数,这些简并的情况是由自发对称性破缺引起的.首先,揭示了在二维系统中自发对称性破缺引起的相变可以运用单点基态保真度的分叉来研究,也反映了在二维系统中约化保真度同样有一个分叉行为;再者,还提出了二维系统的普适序参量以及多分量的复数局域序参量的行为来尝试研究二维3态Potts模型,共同确定系统发生的量子相变的临界点及其类型.即基于iPEPS算法,从单点基态保真度、约化保真度、普适序参量以及局域序参量的角度,来研究3态Potts模型的量子相变,其为一级相变.     

Concurrence and entanglement entropy in a dimerised spin-1/2 two-leg ladder

We consider isotropic spin-1/2 two-leg ladders with dominant spatially-modulated rung exchanges. We study the effect of a uniform magnetic field on the ground state phase diagram of the model using perturbation theory and the numerical Lanczos method. The ground state phase diagram consists of two gapless Luttinger liquid (LL) and three gapped phases. Numerically, we calculate the concurrence between two spins and the entanglement entropy between legs. Numerical experiment shows that the gapless LL phases are fundamentally different. In the first LL phase, only spins on rungs are entangled, but in the second LL phase the spins on legs are long-distance entangled. Therefore, the concurrence between spins on legs can be considered as a function to distinguish the LL phases.     

Quantum phase transitions in matrix product states of one-dimensional spin-1/2 chains

For the matrix product system of a one-dimensional spin-1/2 chain, we present a new model of quantum2 phase transitions and find that in the thermodynamic limit, both sides of the critical point are respectively described by phases ｜Ψa 〉=｜1··· 1 representing all particles spin up and ｜Ψb 〉=｜0··· 0 representing all particles spin down, while the phase transition point is an isolated intermediate-coupling point where√ the two phases coexist equally, which is2 described by the so-called N-qubit maximally entangled GHZ state ｜Ψpt =√2/2（｜1··· 1 ＋｜0··· 0）. At the critical point,2the physical quantities including the entanglement are not discontinuous and the matrix product system has longrange correlation and N-qubit maximal entanglement. We believe that our work is helpful for having a comprehensive understanding of quantum phase transitions in matrix product states of one-dimensional spin chains and of potential directive significance to the preparation and control of one-dimensional spin lattice models with stable coherence and N-qubit maximal entanglement.     

Diagonal entropy in many-body systems: Volume effect and quantum phase transitions

We investigate the diagonal entropy(DE) of the ground state for quantum many-body systems, including the XY model and the Ising model with next nearest neighbor interactions. We focus on the DE of a subsystem of L continuous spins. We show that the DE in many-body systems, regardless of integrability, can be represented as a volume term plus a logarithmic correction and a constant offset. Quantum phase transition points can be explicitly identified by the three coefficients thereof. Besides, by combining entanglement entropy and the relative entropy of quantum coherence, as two celebrated representatives of quantumness, we simply obtain the DE, which naturally has the potential to reveal the information of quantumness. More importantly, the DE is concerning only the diagonal form of the ground state reduced density matrix, making it feasible to measure in real experiments, and therefore it has immediate applications in demonstrating quantum supremacy on state-of-the-art quantum simulators.     

自旋-1/2量子罗盘链的量子相与相变

基于矩阵乘积态表述的无限时间演化块算法,研究了具有x,y,z三个自旋方向的轨道自由度和轨道序竞争的量子罗盘自旋链模型.为了刻画该模型的量子相和相变,计算了基态能量、局域序参量、弦关联序参量、临界指数、冯诺依曼熵、有限纠缠标度和中心荷.结果表明:该量子基态相图由条纹反铁磁相、反铁磁相、单调奇数Haldane相和振荡奇数Haldane相构成.从条纹反铁磁相到反铁磁相,以及从单调奇数Haldane相到振荡奇数Haldane相发生了非连续相变;从振荡奇数Haldane相到条纹反铁磁相,以及从反铁磁相到单调奇数Haldane相发生了连续相变;连续相变线和非连续相变线的交点是多临界点.此外,连续相变点处的临界指数β=1/8和中心荷c=1/2表明连续相变的普适类属于Ising类.由此揭示了该模型量子基态相图的本性,对今后研究更高自旋以及更为复杂轨道序竞争的量子罗盘链模型的量子相与相变具有一定借鉴与参考意义.     

两粒子自旋系统量子纠缠的时间演化

研究了两个具有海森伯耦合的自旋为1/2的粒子在随时间变化的磁场中的运动情况.系统的哈密顿量具有SU（2）代数结构,运用代数动力学方法对此系统进行求解,得到了时间演化算子的严格解.基于严格解,求得两粒子体系随时间变化的波函数,从而计算得到两粒子体系的纠缠.对不同初始波函数,研究了系统纠缠随时间的变化情况.讨论了外场影响纠缠的条件. 关键词： 二粒子系统纠缠 代数动力学解法