A projection method for statics and dynamics of lattice spin systems

A method based on Monte Carlo sampling of the probability flows projected onto the subspace of one or more slow variables is proposed for investigation of dynamic and static properties of lattice spin systems. We illustrate the method by applying it, with projection onto the order-parameter subspace, to the three-dimensional 3-state Potts model in equilibrium and to metastable decay in a three-dimensional 3-state kinetic Potts model.Comment: 4 pages, 3 figures, RevTex, final version to appear in Phys. Rev. Let

Chiral Limit of Strongly Coupled Lattice Gauge Theories

We construct a new and efficient cluster algorithm for updating strongly coupled U(N) lattice gauge theories with staggered fermions in the chiral limit. The algorithm uses the constrained monomer-dimer representation of the theory and should also be of interest to researchers working on other models with similar constraints. Using the new algorithm we address questions related to the chiral limit of strongly coupled U(N) gauge theories beyond the mean field approximation. We show that the infinite volume chiral condensate is non-zero in three and four dimensions. However, on a square lattice of size $L$ we find $\sum_x \sim L^{2-\eta}$ for large $L$ where $\eta = 0.420(3)/N + 0.078(4)/N^2$. These results differ from an earlier conclusion obtained using a different algorithm. Here we argue that the earlier calculations were misleading due to uncontrolled autocorrelation times encountered by the previous algorithm.Comment: 36 Pages, 9 figures, aps revtex forma

Algebraic Bethe ansatz approach for the one-dimensional Hubbard model

We formulate in terms of the quantum inverse scattering method the algebraic Bethe ansatz solution of the one-dimensional Hubbard model. The method developed is based on a new set of commutation relations which encodes a hidden symmetry of 6-vertex type.Comment: appendix additioned with Boltzmann weigths and R-matrix. Version to be published in J.Phys.A:math.Gen. (1997

SO(4) Symmetry of the Transfer Matrix for the One-Dimensional Hubbard Model

The SO(4) invariance of the transfer matrix for the one-dimensional Hubbard model is clarified from the QISM (quantum inverse scattering method) point of view. We demonstrate the SO(4) symmetry by means of the fermionic R-matrix, which satisfy the graded Yang-Baxter relation. The transformation law of the fermionic L-operator under the SO(4) rotation is identified with a kind of gauge transformation, which determines the corresponding transformation of the fermionic creation and annihilation operators under the SO(4) rotation. The transfer matrix is confirmed to be invariant under the SO(4) rotation, which ensures the SO(4) invariance of the conserved currents including the Hamiltonian. Furthermore, we show that the representation of the higher conserved currents in terms of the Clifford algebra gives manifestly SO(4) invariant forms.Comment: 20 pages, LaTeX file using citesort.st

Direct Observation of Field-Induced Incommensurate Fluctuations in a One-Dimensional S=1/2 Antiferromagnet

Neutron scattering from copper benzoate, Cu(C6D5COO)2 3D2O, provides the first direct experimental evidence for field-dependent incommensurate low energy modes in a one-dimensional spin S = 1/2 antiferromagnet. Soft modes occur for wavevectors q=\pi +- dq(H) where dq(H) ~ 2 \pi M(H)/g\mu_B as predicted by Bethe ansatz and spinon descriptions of the S = 1/2 chain. Unexpected was a field-induced energy gap $\Delta(H) \propto H^\alpha$, where $\alpha = 0.65(3)$ as determined from specific heat measurements. At H = 7 T (g\mu_B H/J = 0.52), the magnitude of the gap varies from 0.06 - 0.3 J depending on the orientation of the applied field.Comment: 11 pages, 5 postscript figures, LaTeX, Submitted to PRL 3/31/97, e-mail comments to [email protected]

Dynamical Structure Factors of the S=1/2 Bond-Alternating Spin Chain with a Next-Nearest-Neighbor Interaction in Magnetic Fields

The dynamical structure factor of the S=1/2 bond-alternating spin chain with a next-nearest-neighbor interaction in magnetic field is investigated using the continued fraction method based on the Lanczos algorithm. When the plateau exists on the magnetization curve, the longitudinal dynamical structure factor shows a large intensity with a periodic dispersion relation, while the transverse one shows a large intensity with an almost dispersionless mode. The periodicity and the amplitude of the dispersion relation in the longitudinal dynamical structure factor are sensitive to the coupling constants. The dynamical structure factor of the S=1/2 two-leg ladder in magnetic field is also calculated in the strong interchain-coupling regime. The dynamical structure factor shows gapless or gapful behavior depending on the wave vector along the rung.Comment: 8 pages, 4 figures, to appear in Journal of the Physical Society of Japan, vol. 69, no. 10, (2000

Energy level statistics of the two-dimensional Hubbard model at low filling

The energy level statistics of the Hubbard model for $L \times L$ square lattices (L=3,4,5,6) at low filling (four electrons) is studied numerically for a wide range of the coupling strength. All known symmetries of the model (space, spin and pseudospin symmetry) have been taken into account explicitly from the beginning of the calculation by projecting into symmetry invariant subspaces. The details of this group theoretical treatment are presented with special attention to the nongeneric case of L=4, where a particular complicated space group appears. For all the lattices studied, a significant amount of levels within each symmetry invariant subspaces remains degenerated, but except for L=4 the ground state is nondegenerate. We explain the remaining degeneracies, which occur only for very specific interaction independent states, and we disregard these states in the statistical spectral analysis. The intricate structure of the Hubbard spectra necessitates a careful unfolding procedure, which is thoroughly discussed. Finally, we present our results for the level spacing distribution, the number variance $\Sigma^2$, and the spectral rigidity $\Delta_3$, which essentially all are close to the corresponding statistics for random matrices of the Gaussian ensemble independent of the lattice size and the coupling strength. Even very small coupling strengths approaching the integrable zero coupling limit lead to the Gaussian ensemble statistics stressing the nonperturbative nature of the Hubbard model.Comment: 31 pages (1 Revtex file and 10 postscript figures

Antiferromagnetism in the Exact Ground State of the Half Filled Hubbard Model on the Complete-Bipartite Graph

As a prototype model of antiferromagnetism, we propose a repulsive Hubbard Hamiltonian defined on a graph \L={\cal A}\cup{\cal B} with ${\cal A}\cap {\cal B}=\emptyset$ and bonds connecting any element of ${\cal A}$ with all the elements of ${\cal B}$. Since all the hopping matrix elements associated with each bond are equal, the model is invariant under an arbitrary permutation of the ${\cal A}$-sites and/or of the ${\cal B}$-sites. This is the Hubbard model defined on the so called $(N_{A},N_{B})$-complete-bipartite graph, $N_{A}$ ($N_{B}$) being the number of elements in ${\cal A}$ (${\cal B}$). In this paper we analytically find the {\it exact} ground state for $N_{A}=N_{B}=N$ at half filling for any $N$; the repulsion has a maximum at a critical $N$-dependent value of the on-site Hubbard $U$. The wave function and the energy of the unique, singlet ground state assume a particularly elegant form for N \ra \inf. We also calculate the spin-spin correlation function and show that the ground state exhibits an antiferromagnetic order for any non-zero $U$ even in the thermodynamic limit. We are aware of no previous explicit analytic example of an antiferromagnetic ground state in a Hubbard-like model of itinerant electrons. The kinetic term induces non-trivial correlations among the particles and an antiparallel spin configuration in the two sublattices comes to be energetically favoured at zero Temperature. On the other hand, if the thermodynamic limit is taken and then zero Temperature is approached, a paramagnetic behavior results. The thermodynamic limit does not commute with the zero-Temperature limit, and this fact can be made explicit by the analytic solutions.Comment: 19 pages, 5 figures .ep

Complete Pseudohole and Heavy-Pseudoparticle Operator Representation for the Hubbard Chain

We introduce the pseudohole and heavy-pseudoparticle operator algebra that generates all Hubbard-chain eigenstates from a single reference vacuum. In addition to the pseudoholes already introduced for the description of the low-energy physics, this involves the heavy pseudoparticles associated with Hamiltonian eigenstates whose energy spectrum has a gap relatively to the many-electron ground state. We introduce a generalized pseudoparticle perturbation theory which describes the relevant finite-energy ground state transitions. In the present basis these excitations refer to a small density of excited pseudoparticles. Our operator basis goes beyond the Bethe-ansatz solution and it is the suitable and correct starting point for the study of the finite-frequency properties, which are of great relevance for the understanding of the unusual spectral properties detected in low-dimensional novel materials.Comment: LaTeX, 32 pages, no Figures. To be published in Phys. Rev. B (15th of August 1997