Search

Gauge theories with matter fields in various representations play an important role in different branches of physics. Recently, it was proposed that several aspects of the interesting pseudogap phase of cuprate superconductors near optimal doping may be explained by an emergent $SU(2)$ gauge symmetry. Around the transition with positive hole-doping, one can construct a $(2+1)-$dimensional $SU(2)$ gauge theory coupled to four adjoint scalar fields which gives rise to a rich phase diagram with a myriad of phases having different broken symmetries. We study the phase diagram of this model on the Euclidean lattice using the Hybrid Monte Carlo algorithm. We find the existence of multiple broken phases as predicted by previous mean field studies. Depending on the quartic couplings, the $SU(2)$ gauge symmetry is broken down either to $U(1)$ or $\mathbb{Z}_2$ in the perturbative description of the model. We further study the confinement-deconfinement transition in this theory, and find that both the broken phases are deconfining in the range of volumes that we studied. However, there exists a marked difference in the behavior of the Polyakov loop between the two phases., Comment: 11 pages, 11 figures, minor changes to results and conclusions

Using tensor network methods, we simulate the real-time evolution of the lattice Thirring model quenched out of equilibrium in both the critical and massive phases, and study the appearance of dynamical quantum phase transitions, as non-analyticities in the Loschmidt rate. Whereas the presence of a dynamical quantum phase transition in the model does not correspond to quenches across the critical line of the equilibrium phase diagram at zero temperature, we identify a threshold in the energy density of the initial state, necessary for a dynamical quantum phase transition to be present. Moreover, in the case of the gapped quench Hamiltonian, we unveil a connection of this threshold to a transition between different regions in the finite temperature phase diagram.

Electronic liquid crystal (ELC) phases are spontaneous symmetry breaking states believed to arise from strong electron correlation in quantum materials such as cuprates and iron pnictides. Here, we report a direct observation of ELC phases in a Dirac nodal line (DNL) semimetal GdSbxTe2-x. Electronic nanostructures consisting of incommensurate smectic charge modulation and intense local nematic order are visualized by using spectroscopic imaging - scanning tunneling microscopy. As topological materials with symmetry protected Dirac or Weyl fermions are mostly weakly correlated, the discovery of such ELC phases are anomalous and raise questions on the origin of their emergence. Specifically, we demonstrate how chemical substitution generates these symmetry breaking phases before the system undergoes a charge density wave - orthorhombic structural transition. We further show how dopants can induce nematicity via quasiparticle scattering interference. Our results highlight the importance of impurities in realizing ELC phases and present a new material platform for exploring the interplay among quenched disorder, topology and electron correlation.

We introduce a tensor network library designed for classical and quantum physics simulations called Cytnx (pronounced as sci-tens). This library provides almost an identical interface and syntax for both C++ and Python, allowing users to effortlessly switch between two languages. Aiming at a quick learning process for new users of tensor network algorithms, the interfaces resemble the popular Python scientific libraries like NumPy, Scipy, and PyTorch. Not only multiple global Abelian symmetries can be easily defined and implemented, Cytnx also provides a new tool called Network that allows users to store large tensor networks and perform tensor network contractions in an optimal order automatically. With the integration of cuQuantum, tensor calculations can also be executed efficiently on GPUs. We present benchmark results for tensor operations on both devices, CPU and GPU. We also discuss features and higher-level interfaces to be added in the future.

We study a 3-dimensional SU(2) gauge theory with 4 Higgs fields which transform under the adjoint representation of the gauge group, that has been recently proposed by Sachdev et al. to explain the physics of cuprate superconductors near optimal doping. The symmetric confining phase of the theory corresponds to the usual Fermi-liquid phase while the broken (Higgs) phase is associated with the interesting pseudogap phase of cuprates. We employ the Hybrid Monte-Carlo algorithm to study the phase diagram of the theory. We find the existence of a variety of broken phases in qualitative accordance with earlier mean-field predictions and discuss their role in cuprates. In addition, we investigate the behavior of Polyakov loop to probe the confinement/deconfinement phase transition, and find that the Higgs phase hosts a stable deconfining phase consistent with previous studies., Comment: 7 pages, 4 figures, Proceedings of the 40th International Symposium on Lattice Field Theory (LATTICE2023), July 31st - August 4th, 2023, Fermi National Accelerator Laboratory

Topology plays a crucial role in many physical systems, leading to interesting states at the surface. The paradigmatic example is the Chern number defined in the Brillouin zone that leads to the robust gapless edge states. Here we introduce the reduced Chern number, defined in subregions of the Brillouin zone (BZ), and construct a family of Chern dartboard insulators (CDIs) with quantized reduced Chern numbers in the sub-BZ (sBZ) but with trivial bulk topology. CDIs are protected by mirror symmetries and exhibit distinct pseudospin textures, including (anti)skyrmions, inside the sBZ. These CDIs host exotic gapless edge states, such as M\"{o}bius fermions and midgap corner states, and can be realized in photonic crystals. Our work opens up new possibilities for exploring sBZ topology and nontrivial surface responses in topological systems.

We propose a simple and generic construction of the variational tensor network operators to study the quantum spin systems by the synergy of ideas from the imaginary-time evolution and variational optimization of trial wave functions. By applying these operators to simple initial states, accurate variational ground state wave functions with extremely few parameters can be obtained. Furthermore, the framework can be applied to study spontaneously symmetry breaking, symmetry protected topological, and intrinsic topologically ordered phases, and we show that symmetries of the local tensors associated with these phases can emerge directly after the optimization without any gauge fixing. This provides a universal way to identify quantum phase transitions without prior knowledge of the system., Comment: 16 pages, 12 figures

Bulk-boundary correspondence serves as an important feature of the strong topological insulators, including Chern insulators and $Z_2$ topological insulators. Under nontrivial band topology, the protected gapless edge states correspond to the Wannier obstruction or Wilson-loop winding in the bulk. Recent studies show that the bulk topological features may not imply the existence of protected gapless edge states. Here we address the opposite question: Does the existence of protected gapless edge states necessarily imply the Wannier obstruction or Wilson-loop winding? We provide an example where the protected gapless edge states arise without the aforementioned bulk topological features. This trivialized topological insulator belongs to a new class of systems with non-delta-like Wannier functions. Interestingly, the gapless edge states are not protected by the crystalline symmetry; instead the protection originates from the mirror antisymmetry, a combination of chiral and mirror symmetries. Although the protected gapless edge states cannot be captured by the bulk topological features, they can be characterized by the spectral flow in the entanglement spectrum., Comment: 7 pages, 4 figures

We numerically study the phase structure of the CP(1) model in the presence of a topological $\theta$-term, a regime afflicted by the sign problem for conventional lattice Monte Carlo simulations. Using a bond-weighted tensor renormalization group method, we compute the free energy for inverse couplings ranging from $0\leq \beta \leq 1.1$ and find a CP-violating, first-order phase transition at $\theta=\pi$. In contrast to previous findings, our numerical results provide no evidence for a critical coupling $\beta_c<1.1$ above which a second-order phase transition emerges at $\theta=\pi$ and/or the first-order transition line bifurcates at $\theta\neq\pi$. If such a critical coupling exists, as suggested by Haldane's conjecture, our study indicates that is larger than $\beta_c>1.1$., Comment: 10 pages, 9 figures, v2: updated to match journal version

We study the excitation spectrum of the spin-1 Kitaev model using the symmetric tensor network. By evaluating the virtual order parameters defined on the virtual Hilbert space in the tensor network formalism, we confirm the ground state is in a $\mathbb{Z}_2$ spin liquid phase. Using the correspondence between the transfer matrix spectrum and low-lying excitations, we find that contrary to the dispersive Majorana excitation in the spin-1/2 case, the isotropic spin-1 Kitaev model has a dispersive charge anyon excitation. Bottom of the gapped single-particle charge excitations are found at $\mathbf{K}, \mathbf{K}'=(\pm2\pi/3, \mp 2\pi/3)$, with a corresponding correlation length of $\xi \approx 6.7$ unit cells. The lower edge of the two-particle continuum, which is closely related to the dynamical structure factor measured in inelastic neutron scattering experiments, is obtained by extracting the excitations in the vacuum superselection sector in the anyon theory language, Comment: 7 pages, 8 figures

We develop two novel numerical schemes to study the conductance of the two-wire junction of inequivalent Tomonaga-Luttinger Liquids. In the first scheme we use the static current-current correlation function across the junction to extract the linear conductance through a relation that is derived via the bosonization method. In the second scheme we apply a bias and evaluate the time-dependent current across the junction to obtain the current-voltage characteristic. The conductance is then extracted from the small bias result within the linear response regime. Both schemes are based on the infinite size matrix product state to minimize the finite-size effects. Due to the lack of the translational invariance, we focus on a finite-size window containing the junction. For time-independent calculations, we use infinite boundary conditions to evaluate the correlations within the window. For time-dependent calculations, we use the window technique to evaluate the local currents within the window. The numerical results obtained by both schemes show excellent agreement with the analytical predictions., Comment: 11 pages, 10 figures

We propose a unified scheme to identify phase transitions out of the $\mathbb{Z}_2$ Abelian topological order, including the transition to a non-Abelian chiral spin liquid. Using loop gas and and string gas states [H.-Y. Lee, R. Kaneko, T. Okubo, N. Kawashima, Phys. Rev. Lett. 123, 087203 (2019)] on the star lattice Kitaev model as an example, we compute the overlap of minimally entangled states through transfer matrices. We demonstrate that, similar to the anyon condensation, continuous deformation of a $\mathbb{Z}_2$-injective projected entangled-pair state (PEPS) also allows us to study the transition between Abelian and non-Abelian topological orders. We show that the charge and flux anyons defined in the Abelian phase transmute into the $\sigma$ anyon in the non-Abelian topological order. Furthermore, we show that contrary to the claim in [Phys. Rev. B 101, 035140 (2020)], both the LG and SG states have infinite correlation length in the non-Abelian regime, consistent with the no-go theorem that a chiral PEPS has a gapless parent Hamiltonian., Comment: 9 pages, 9 figures, published version

We present concluding results from our study for zero-temperature phase structure of the massive Thirring model in 1+1 dimensions with staggered regularisation. Employing the method of matrix product states, several quantities, including two types of correlators, are investigated, leading to numerical evidence of a Berezinskii-Kosterlitz-Thouless phase transition. Exploratory results for real-time dynamics pertaining to this transition, obtained using the approaches of variational uniform matrix product state and time-dependent variational principle, are also discussed., Comment: 7 pages, 4 figures, contribution to the 37th International Symposium on Lattice Field Theory (LATTICE2019), 16-22 June 2019, Wuhan, China

Quantum entanglement is fragile to thermal fluctuations, which raises the question whether finite temperature phase transitions support long-range entanglement similar to their zero temperature counterparts. Here we use quantum Monte Carlo simulations to study the third Renyi negativity, a generalization of entanglement negativity, as a proxy of mixed-state entanglement in the 2D transverse field Ising model across its finite temperature phase transition. We find that the area-law coefficient of the Renyi negativity is singular across the transition, while its subleading constant is zero within the statistical error. This indicates that the entanglement is short-ranged at the critical point despite a divergent correlation length. Renyi negativity in several exactly solvable models also shows qualitative similarities to that in the 2D transverse field Ising model., Comment: 5 pages, 3 figures, 4 page appendix

We study the two-dimensional kagome-ice model derived from a pyrochlore lattice with second- and third-neighbor interactions. The canted moments align along the local $\langle 111 \rangle$ axes of the pyrochlore and respond to both in-plane and out-of-plane external fields. We find that the combination of further-neighbor interactions together with the external fields introduces a rich phase diagram with different spin textures. Close to the phase boundaries, metastable $\textit{"snake"}$ domains emerge with extremely long relaxation time. Our kinetic Monte Carlo analysis of the magnetic-field quench process from saturated state shows unusually slow dynamics. Despite that the interior spins are almost frozen in snake domains, the spins on the edge are free to fluctuate locally, leading to frequent creation and annihilation of monopole-anti-monopole bound states. Once the domains are formed, these excitations are localized and can hardly propagate due to the energy barrier of snakes. The emergence of such snake domains may shed light on the experimental observation of dipolar spin ice under tilted fields, and provide a new strategy to manipulate both spin and charge textures in artificial spin ice., Comment: 8 pages, 6 figures

Employing matrix product states as an ansatz, we study the non-thermal phase structure of the (1+1)-dimensional massive Thirring model in the sector of vanishing total fermion number with staggered regularization. In this paper, details of the implementation for this project are described. To depict the phase diagram of the model, we examine the entanglement entropy, the fermion bilinear condensate and two types of correlation functions. Our investigation shows the existence of two phases, with one of them being critical and the other gapped. An interesting feature of the phase structure is that the theory with non-zero fermion mass can be conformal. We also find clear numerical evidence that these phases are separated by a transition of the Berezinskii-Kosterlitz-Thouless type. Results presented in this paper establish the possibility of using the matrix product states for probing this type of phase transition in quantum field theories. They can provide information for further exploration of scaling behaviour, and serve as an important ingredient for controlling the continuum extrapolation of the model., Comment: 31 pages, 18 figures; minor changes to the text, typos corrected, references added; version published in Physical Review D

We present preliminary results of a study on the non-thermal phase structure of the (1+1) dimensional massive Thirring model, employing the method of matrix product states. Through investigating the entanglement entropy, the fermion correlators and the chiral condensate, it is found that this approach enables us to observe numerical evidence of a Kosterlitz-Thouless phase transition in the model., Comment: 7 pages, 4 figures; contribution to the proceedings of Lattice 2018 conference

A preponderance of evidence suggests that the ground state of the nearest-neighbor $S = 1/2$ antiferromagnetic Heisenberg model on the kagome lattice is a gapless spin liquid. Many candidate materials for the realization of this model possess in addition a Dzyaloshinskii-Moriya (DM) interaction. We study this system by tensor-network methods and deduce that a weak but finite DM interaction is required to destabilize the gapless spin-liquid state. The critical magnitude, $D_c/J \simeq 0.012(2)$, lies well below the DM strength proposed in the kagome material herbertsmithite, indicating a need to reassess the apparent spin-liquid behavior reported in this system., Comment: 9 pages, 9 figures, to appear in Phys. Rev. B

Quantum effect is expected to dictate the behaviour of physical systems at low temperature. For quantum magnets with geometrical frustration, quantum fluctuation usually lifts the macroscopic classical degeneracy, and exotic quantum states emerge. However, how different types of quantum processes entangle wave functions in a constrained Hilbert space is not well understood. Here, we study the topological entanglement entropy (TEE) and the thermal entropy of a quantum ice model on a geometrically frustrated kagome lattice. We find that the system does not show a $Z_2$ topological order down to extremely low temperature, yet continues to behave like a classical kagome ice with finite residual entropy. Our theoretical analysis indicates an intricate competition of off-diagonal and diagonal quantum processes leading to the quasi-degeneracy of states and effectively, the classical degeneracy is restored., Comment: 14 pages, 16 figures, added ED results for the effective model

An impurity in a Tomonaga-Luttinger liquid leads to a crossover between short- and long-distance regime which describes many physical phenomena. However, calculation of the entire crossover of correlation functions over different length scales has been difficult. We develop a powerful numerical method based on infinite DMRG utilizing a finite system with infinite boundary conditions, which can be applied to correlation functions near an impurity. For the $S=1/2$ chain, we demonstrate that the full crossover can be precisely obtained, and that their limiting behaviors show a good agreement with field-theory predictions.

Tensor Network methods have been established as a powerful technique for simulating low dimensional strongly-correlated systems for over two decades. Employing the formalism of Matrix Product States, we investigate the phase diagram of the massive Thirring model. We also show the possibility of studying soliton dynamics and topological phase transition via the Thirring model., Comment: 8 pages, 3 figures; contribution to the proceedings of the 35th International Symposium on Lattice Field Theory, 18 - 24 June 2017, Granada, Spain

The spin-1 Haldane chain is an example of the symmetry-protected-topological (SPT) phase in one dimension. Experimental realization of the spin chain materials usually involves both the uniaxial-type, $D(S^z)^2$, and the rhombic-type, $E[(S^x)^2-(S^y)^2]$, single-ion anisotropies. Here, we provide a precise ground-state phase diagram for spin-1 Haldane chain with these single-ion anisotropies. Using quantum numbers, we find that the $\mathbb{Z}_2$ symmetry breaking phase can be characterized by double degeneracy in the entanglement spectrum. Topological quantum phase transitions take place on particular paths in the phase diagram, from the Haldane phase to the Large-$E_x$, Large-$E_y$, or Large-$D$ phases. The topological critical points are determined by the level spectroscopy method with a newly developed parity technique in the density matrix renormalization group [Phys. Rev. B 86, 024403 (2012)], and the Haldane-Large-$D$ critical point is obtained with an unprecedentedly precision, $(D/J)_c=0.9684713(1)$. Close to this critical point, a small rhombic single-ion anisotropy $|E|/J\ll1$ can destroy the Haldane phase and bring the system into a $y$-N\'eel phase. We propose that the compound [Ni(HF$_2$)(3-Clpy)$_4$]BF$_4$ is a candidate system to search for the $y$-N\'eel phase., Comment: Supplemental Material is included

Directly in the thermodynamic limit, we show how to combine imaginary and real time evolution of tensor networks to efficiently and accurately find the nonequilibrium steady states (NESS) of one-dimensional dissipative quantum lattices governed by the Lindblad master equation. The imaginary time evolution first bypasses any highly correlated portions of the real-time evolution trajectory by directly converging to the weakly correlated subspace of the NESS, after which real time evolution completes the convergence to the NESS with high accuracy. We demonstrate the power of the method with the dissipative transverse field quantum Ising chain. We show that a crossover of an order parameter shown to be smooth in previous finite-size studies remains smooth in the thermodynamic limit., Comment: 5+3 pages, 5 figures, 2 tables

We present numerical studies of dipolar spin ice in the presence of a magnetic field slightly tilted away from the [111] axis. We find a first-order transition from a kagome ice to a $\mathbf{q}=$X state when the external field is tilted toward the $[11\bar{2}]$ direction. This is consistent with the anomalous critical scattering previously observed in the neutron scattering experiment on the spin ice material $\textrm{Ho}_{2}\textrm{Ti}_{2}\textrm{O}_{7}$ in a tilted field [Nat. Phys. 3, 566 (2007)]. We show that this ordering originates from the antiferromagnetic alignment of spin chains on the kagome planes. The residual entropy of the kagome ice is fully recovered. Our result captures the features observed in the experiments and points to the importance of the dipolar interaction in determining ordered states in the spin ice materials. We place our results in the context of recent susceptibility measurements on $\textrm{Dy}_{2}\textrm{Ti}_{2}\textrm{O}_{7}$, showing two features for a [111] field., Comment: 5+1 pages; 5+2 figures

Inverse melting refers to the rare thermodynamic phenomenon in which a solid melts into a liquid upon cooling, a transition that can occur only when the ordered (solid) phase has more entropy than the disordered (liquid) phase, and that has so far only been observed in a handful of systems. Here we report the first experimental observation for the inverse melting of an electronic liquid crystalline order in strontium-doped lanthanum nickelate, a compound isostructural with the superconducting cuprates, with a hole doping concentration of 1/3. Using x-ray scattering, we demonstrate that the isotropic charge modulation is driven to nematic order by fluctuating spins and shows an inverse melting transition. Using a phenomenological Landau theory, we show that this inverse melting transition is due to the interlayer coupling between the charge and spin orders. This discovery points to the importance of the interlayer correlations in the system, and provides a new perspective to study the intricate nature of the electronic liquid crystal phases in strongly correlated electronic systems, including possibly the Cu- and Fe-based high-Tc superconductors., Comment: 7 pages, 7 figures

We propose an efficient numerical method, which combines the advantages of recently developed tensor-network based methods and standard trial wave functions, to study the ground state properties of quantum many-body systems. In this approach, we apply a projector in the form of a tensor-product operator to an input wave function, such as a Jastrow-type or Hartree-Fock wave function, and optimize the tensor elements via variational Monte Carlo. The entanglement already contained in the input wave function can considerably reduce the bond dimensions compared to the regular tensor-product state representation. In particular, this allows us to also represent states that do not obey the area law of entanglement entropy. In addition, for fermionic systems, the fermion sign structure can be encoded in the input wave function. We show that the optimized states provide good approximations of the ground-state energy and correlation functions in the cases of two-dimensional bosonic and fermonic systems., Comment: 7 pages, 5 figures, published version

We study numerically the universal conductance of Luttinger liquids wire with a single impurity via the Muti-scale Entanglement Renormalization Ansatz (MERA). The scale invariant MERA provides an efficient way to extract scaling operators and scaling dimensions for both the bulk and the boundary conformal field theories. By utilizing the key relationship between the conductance tensor and ground-state correlation function, the universal conductance can be evaluated within the framework of the boundary MERA. We construct the boundary MERA to compute the correlation functions and scaling dimensions for the Kane-Fisher fixed points by modeling the single impurity as a junction (weak link) of two interacting wires. We show that the universal behavior of the junction can be easily identified within the MERA and argue that the boundary MERA framework has tremendous potential to classify the fixed points in general multi-wire junctions., Comment: 14 pages, 18 figures

By using a different quantum-to-classical mapping from the Trotter-Suzuki decomposition, we identify the entanglement structure of the maximal eigenvectors for the associated quantum transfer matrix. This observation provides a deeper insight into the problem of linear growth of the entanglement entropy in time evolution using conventional methods. Based on this observation, we propose a general method for arbitrary temperatures using the biorthonormal transfer-matrix renormalization group. Our method exhibits a competitive accuracy with a much cheaper computational cost in comparison with two recent proposed methods for long-time dynamics based on a folding algorithm [Phys. Rev. Lett. 102, 240603 (2009)] and a modified time-dependent density-matrix renormalization group [Phys. Rev. Lett. 108, 227206 (2012)].

We report here numerical results of the low-temperature behavior of a dipolar spin ice in a magnetic field along the [100] direction. Tuning the magnetic field, the system exhibit a half-magnetization plateau at low temperature. This half-polarized phase should correspond to a quantum solid phase in an effective 2D quantum boson model, and the transition from the Coulomb phase with a power-law correlation to this state can be regarded as a superfluid to a quantum solid transition. We discuss possible experimental signatures of this half-polarized state., Comment: 7 pages, 9 figures; Published version

Symmetry-protected topological (SPT) phases are short-range entangled quantum phases with symmetry, which have gapped excitations in the bulk and gapless modes at the edge. In this paper, we study the SPT phases in the spin-1 Heisenberg chain with a single-ion anisotropy D, which has a quantum phase transition between a Haldane phase and a large-D phase. Using symmetric multiscale entanglement renormalization ansatz (MERA) tensor networks, we study the nonlocal order parameters for time-reversal and inversion symmetry. For the inversion symmetric MERA, we propose a brick-and-rope representation that gives a geometrical interpretation of inversion symmetric tensors. Finally, we study the symmetric renormalization group (RG) flow of the inversion symmetric string-order parameter, and show that entanglement renormalization with symmetric tensors produces proper behavior of the RG fixed-points., Comment: 11 pages, 18 figures

We apply the biorthonormal transfer-matrix renormalization group (BTMRG) [Phys. Rev. E 83, 036702 (2011)] to study low-temperature properties of quantum spin chains. Simulation on isotropic Heisenberg spin-1/2 chain demonstrates that the BTMRG outperforms the conventional transfer-matrix renormalization group (TMRG) by successfully accessing far lower temperature unreachable by conventional TMRG, while retaining the same level of accuracy. The power of the method is further illustrated by the calculation of the low-temperature specific heat for a frustrated spin chain., Comment: 5 pages, 4 figures

In a class of frustrated magnets known as spin ice, magnetic monopoles emerge as classical defects and interact via the magnetic Coulomb law. With quantum-mechanical interactions, these magnetic charges are carried by fractionalised bosonic quasi-particles, spinons, which can undergo Bose-Einstein condensation through a first-order transition via the Higgs mechanism. Here, we report evidence of a Higgs transition from a magnetic Coulomb liquid to a ferromagnet in single-crystal Yb2Ti2O7. Polarised neutron-scattering experiments show that the diffuse [111]-rod scattering and pinch-point features which develop on cooling are suddenly suppressed below T_C ~ 0.21 K, where magnetic Bragg peaks and a full depolarisation of the neutron spins are observed with thermal hysteresis, indicating a first-order ferromagnetic transition. Our results are explained on the basis of a quantum spin-ice model, whose high-temperature phase is effectively described as a magnetic Coulomb liquid, while the ground state shows a nearly collinear ferromagnetism with gapped spin excitations., Comment: 36 pages, including 3 figures and supplement with 4 more figures

We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E, 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg J1-J2 model on a L*L square lattice with L = 8, 16 and 32. By treating tensor elements as variational parameters, we obtain the ground states for different J2/J1 values, and investigate the staggered magnetic order parameters, the nearest-neighbor spin-spin correlations and the plaquette order parameters. In addition to the well-known N\'eel-order and collinear-order at low and high J2/J1, we observe a plaquette order at intermediate J2/J1. A continuous transition between the N\'eel order and the plaquette order near J2 = 0.40J1 is observed. The collinear order emerges at J2 = 0.62J1 through a first-order phase transition., Comment: 6 pages, 5 figures total; added one figure; Revised argument in section 3, results unchanged; added note and reference wherein

We study the electronic transport in a molecular junction in which each site is coupled to a local phonon bath using the non-equilibrium Green's function method. We observe the length period of the oscillatory conductance in odd-numbered chains depends strongly on the applied bias, and the oscillatory behavior is smeared out for the bias voltage near the phonon energy. In addition, a crossover from tunneling to thermally activated hopping transport as the length of the molecule increases is found for the phonon-free case. In the presence of electron-phonon interaction, hopping transport is dominant and a transition from the thermally suppressed to assisted conduction is observed., Comment: 9 pages, 10 figures

We discuss variationally optimized matrix-product states for the transverse-field Ising chain, using D*D matrices with small D=2-10. For finite system size N there are energy minimums for symmetric as well as symmetry-broken states, which cross each other at a field value hc(N,D); thus the transition is first-order. A continuous transition develops as N->infinity. The asymptotic critical behavior is then always of mean-field type (the magnetization exponent beta=1/2), but a window of field strengths where true Ising scaling holds (beta=1/8) emerges with increasing D. We also demonstrate asymptotic mean-field behavior for infinite-size two-dimensional tensor-product (iPEPS) states with small tensors., Comment: 4 pages, 5 figures

We study the fermionic spectral function in a holographic superconductor model. At zero temperature, the black hole has zero horizon and hence the entropy of the system is zero after the back reaction of the condensate is taken into account. We find the system exhibits the famous peak-dip-hump lineshape with a sharp low-energy peak followed by a dip then a hump at higher energies. This feature is widely observed in the spectrum of several high-T_c superconductors. We also find a linear relation between the gap in the fermionic spectrum and the condensate, indicating the condensate is formed by fermion pairing., Comment: 4 pages, revtex4

We present a method for contracting a square-lattice tensor network in two dimensions, based on auxiliary tensors accomplishing successive truncations (renormalization) of 8-index tensors for 2 by 2 plaquettes into 4-index tensors. The scheme is variational, and thus the tensors can be optimized by minimizing the energy. Test results for the quantum phase transition of the transverse-field Ising model confirm that even the smallest possible tensors (two values for each tensor index at each renormalization level) produce much better results than the simple product (mean-field) state., Comment: 4 pages, 7 figures

Using the example of Zn-doped La2CuO4, we demonstrate that a spinless impurity doped into a non-frustrated antiferromagnet can induce substantial frustrating interactions among the spins surrounding it. This counterintuitive result is the key to resolving discrepancies between experimental data and earlier theories. Analytic and quantum Monte Carlo studies of the impurity-induced frustration are in a close accord with each other and experiments. The mechanism proposed here should be common to other correlated oxides as well., Comment: 4 pages, updated figures, accepted version

We study the ground state phase diagram of the hard-core extended boson Hubbard model on the square lattice with both nearest- (nn) and next-nearest-neighbor (nnn) hopping and repulsion, using Gutzwiller mean field theory and quantum Monte Carlo simulations. We observe the formation of supersolid states with checkerboard, striped, and quarter-filled crystal structures, when the system is doped away from commensurate fillings. In the striped supersolid phase, a strong anisotropy in the superfluid density is obtained from the simulations; however, the transverse component remains finite, indicating a true two-dimensional superflow. We find that upon doping, the striped supersolid transitions directly into the supersolid with quarter-filled crystal structure, via a first-order stripe melting transition., Comment: Revtex 4, 6 pages, 9 figures

We study effects of Fermi surface fluctuations on a single-particle life time near the diagonal electronic nematic phase on a two-dimensional square lattice. It has been shown that there exists a quantum critical point (QCP) between the diagonal nematic and isotropic phases. We study the longitudinal fluctuations of the order parameter near the critical point, where the singular forward scattering leads to a non-Fermi liquid behavior over the whole Fermi surface except along the k_x- and k_y-directions. We will also discuss the temperature and chemical potential dependence of the single-particle decay rate., Comment: 4 pages, 3 figures, revtex 4

We study spin and charge susceptibilities in the d-wave superconducting state whose underlying electronic dispersion is anisotropic due to the formation of the electronic nematic order. We show that the amplitude of the incommensurate peaks in the spin susceptibility near (\pi,\pi) reveals a pronounced anisotropy in the momentum space. The relevance of our findings to the magnetic scattering pattern observed in a recent neutron scattering measurement on untwinned YBa_2Cu_3O_{6+x} is discussed. In the charge channel, we identify a well-defined collective mode at small momentum transfer with strong anisotropic amplitude depending on the direction of momentum transfer, which is associated with the broken symmetry due to nematic ordering., Comment: Revtex 4, 7 pages, 4 figures, references added

We study possible quantum ground states of the Sp(N) generalized Heisenberg model on a cubic lattice with nearest-neighbor and next-nearest-neighbor exchange interactions. The phase diagram is obtained in the large-N limit and fluctuation effects are considered via appropriate gauge theories. In particular, we find three U(1) spin liquid phases with different short-range magnetic correlations. These phases are characterized by deconfined gapped spinons, gapped monopoles, and gapless ``photons''. As N becomes smaller, a confinement transition from these phases to valence bond solids (VBS) may occur. This transition is studied by using duality and analyzing the resulting theory of monopoles coupled to a non-compact dual gauge field; the condensation of the monopoles leads to VBS phases. We determine the resulting VBS phases emerging from two of the three spin liquid states. On the other hand, the spin liquid state near J_1 \approx J_2 appears to be more stable against monopole condensation and could be a promising candidate for a spin liquid state in real systems., Comment: revtex file 12 pages, 17 figures

We study the sensitivity of $T_c$ and the pseudogap onset temperature, $T^*$, to low fields, $H$, for cuprate superconductors, using a BCS-based approach extended to arbitrary coupling. We find that $T^*$ and $T_c$, which are of the same superconducting origin, have very different $H$ dependences. The small coherence length makes $T^*$ rather insensitive to the field. However, the presence of the pseudogap at $T_c$ makes $T_c$ more sensitive to $H$. Our results for the coherence length $\xi$ fit well with existing experiments. We predict that very near the insulator $\xi$ will rapidly increase., Comment: 4 pages, 1 figure, contribution to the PPHMF-IV conference, Oct. 2001

We demonstrate that the observed dependences of $T_c$ and $T^*$ on small magnetic fields can be readily understood in a precursor superconductivity approach to the pseudogap phase. In this approach, the presence of a pseudogap at $T_c$ (but not at $T^*$) and the associated suppression of the density of states lead to very different sensitivities to pair-breaking perturbations for the two temperatures. Our semi-quantitative results address the puzzling experimental observation that the coherence length $\xi$ is weakly dependent on hole concentration $x$ throughout most of the phase diagram. We present our results in a form which can be compared with the recent experiments of Shibauchi et al, and argue that orbital effects contribute in an important way to the $H$ dependence of $T^*$., Comment: 6 pages, 1 figure, elsart.cls included. Submitted to the proceeding of SNS 2001, Chicago

In the last few years evidence has been accumulating that there are a multiplicity of energy scales which characterize superconductivity in the underdoped cuprates. In contrast to the situation in BCS superconductors, the phase coherence temperature Tc is different from the energy gap onset temperature T*. In addition, thermodynamic and tunneling spectroscopies have led to the inference that the order parameter $\Delta_{sc}$ is to be distinguished from the excitation gap $\Delta$; in this way, pseudogap effects persist below Tc. It has been argued by many in the community that the presence of these distinct energy scales demonstrates that the pseudogap is unrelated to superconductivity. In this paper we show that this inference is incorrect. We demonstrate that the difference between the order parameter and excitation gap and the contrasting dependences of T* and Tc on hole concentration $x$ and magnetic field $H$ follow from a natural generalization of BCS theory. This simple generalized form is based on a BCS-like ground state, but with self consistently determined chemical potential in the presence of arbitrary attractive coupling $g$. We have applied this mean field theory with some success to tunneling, transport, thermodynamics and magnetic field effects. We contrast the present approach with the phase fluctuation scenario and discuss key features which might distinguish our precursor superconductivity picture from that involving a competing order parameter., Comment: 4 pages, 2 EPS figures, use LaTeX package espcrc2.sty from Elsevier, submitted to SNS'01 conference proceedings

We study the sensitivity of T_c and T^* to low fields, H, within the pseudogap state using a BCS-based approach extended to arbitrary coupling. We find that T^* and T_c, which are of the same superconducting origin, have very different H dependences. This is due to the pseudogap, \Delta_{pg}, which is present at the latter, but not former temperature. Our results for the coherence length \xi fit well with existing experiments.We predict that very near the insulator \xi will rapidly increase., Comment: 4 pages, 4 figures, RevTex

We study the evolution of the d-wave neutron cross-section with variable frequency \omega and fixed T (below and above Tc) in two different cuprate families. The evolution from incommensurate to commensurate to incommensurate peaks is rather generic within an RPA-like scheme. This behavior seems to be in reasonable accord with experiments, and may help distinguish between this and the "stripe" scenario., Comment: 2 pages; submitted to Proceedings of M2S-HTSC-VI

We study the evolution of the neutron cross-section with variable frequency $\omega$ and fixed $T$ below $T_c$ in two different cuprate families. Our calculations, which predominantly probe the role of d-wave pairing, lead to generic features, independent of Fermi surface shapes. Among our findings, reasonably consistent with experiment, are (i) for $\omega $ near the gap energy $\Delta$, both optimal {LaSrCuO} and slightly underdoped YBCO exhibit (comparably) incommensurate peaks (ii) peak sharpening below $T_c$ is seen in {LaSrCuO}, (iii) quite generically, a frequency evolution from incommensurate to commensurate and then back to incommensurate structure is found with increasing $\omega$. Due to their narrow $\omega$ regime of stability, commensurate peaks in {LaSrCuO} should be extremely difficult to observe., Comment: RevTex 5pages, 4figures; Manuscript rewritten, figures revised, and direct comparisons with experiments added

In this paper we study the BCS Bose Einstein condensation (BEC) crossover scenario within the superconducting state, using a T-matrix approach which yields the ground state proposed by Leggett. Here we extend this ground state analysis to finite temperatures T and interpret the resulting physics. We find two types of bosonic-like excitations of the system: long lived, incoherent pair excitations and collective modes of the superconducting order parameter, which have different dynamics. Using a gauge invariant formalism, this paper addresses their contrasting behavior as a function of T and superconducting coupling constant g. At a more physical level, our paper emphasizes how, at finite T, BCS-BEC approaches introduce an important parameter \Delta^2_{pg} = \Delta^2 - \Delta_{sc}^2 into the description of superconductivity. This parameter is governed by the pair excitations and is associated with particle-hole asymmetry effects which are important for sufficiently large g. In the fermionic regime, \Delta_{pg}^2 represents the difference between the square of the excitation gap \Delta^2 and that of the superconducting order parameter \Delta_{sc}^2. The parameter \Delta_{pg}^2, which is necessarily zero in the BCS (mean field) limit increases monotonically with the strength of the attractive interaction g. It follows that there is a significant physical distinction between this BCS-BEC crossover approach (in which g is the essential variable which determines \Delta_{pg}) and the widely discussed (Coulomb-modulated) phase fluctuation scenario in which the plasma frequency is the tuning parameter., Comment: RevTeX3.x, 22 pages, 6 EPS figures (included)