Dynamics of the Hubbard model: a general approach by time dependent variational principle

We describe the quantum dynamics of the Hubbard model at semi-classical level, by implementing the Time-Dependent Variational Principle (TDVP) procedure on appropriate macroscopic wavefunctions constructed in terms of su(2)-coherent states. Within the TDVP procedure, such states turn out to include a time-dependent quantum phase, part of which can be recognized as Berry's phase. We derive two new semi-classical model Hamiltonians for describing the dynamics in the paramagnetic, superconducting, antiferromagnetic and charge density wave phases and solve the corresponding canonical equations of motion in various cases. Noticeably, a vortex-like ground state phase dynamics is found to take place for U>0 away from half filling. Moreover, it appears that an oscillatory-like ground state dynamics survives at the Fermi surface at half-filling for any U. The low-energy dynamics is also exactly solved by separating fast and slow variables. The role of the time-dependent phase is shown to be particularly interesting in the ordered phases.Comment: ReVTeX file, 38 pages, to appear on Phys. Rev.

Neel Order and Electron Spectral Functions in the Two-Dimensional Hubbard Model: a Spin-Charge Rotating Frame Approach

Using recently developed quantum SU(2)xU(1) rotor approach, that provides a self-consistent treatment of the antiferromagnetic state we have performed electronic spectral function calculations for the Hubbard model on the square lattice. The collective variables for charge and spin are isolated in the form of the space-time fluctuating U(1) phase field and rotating spin quantization axis governed by the SU(2) symmetry, respectively. As a result interacting electrons appear as composite objects consisting of bare fermions with attached U(1) and SU(2) gauge fields. This allows us to write the fermion Green's function in the space-time domain as the product CP^1 propagator resulting from the SU(2) gauge fields, U(1) phase propagator and the pseudo-fermion correlation function. As a result the problem of calculating the spectral line shapes now becomes one of performing the convolution of spin, charge and pseudo-fermion Green's functions. The collective spin and charge fluctuations are governed by the effective actions that are derived from the Hubbard model for any value of the Coulomb interaction. The emergence of a sharp peak in the electron spectral function in the antiferromagnetic state indicates the decay of the electron into separate spin and charge carrying particle excitations.Comment: 16 pages, 5 figures, submitted to Phys. Rev.

The breakdown of the Nagaoka phase in the 2D t-J model

In the limit of weak exchange, J, at low hole concentration, the ground state of the 2D t-J model is believed to be ferromagnetic. We study the leading instability of this Nagaoka state, which emerges with increasing J. Both exact diagonalization of small clusters, and a semiclassical analytical calculation of larger systems show that above a certain critical value of the exchange, Nagaoka's state is unstable to phase separation. In a finite-size system a bubble of antiferromagnetic Mott insulator appears in the ground state above this threshold. The size of this bubble depends on the hole concentration and scales as a power of the system size, N

Antiferromagnetism and single-particle properties in the two-dimensional half-filled Hubbard model: a non-linear sigma model approach

We describe a low-temperature approach to the two-dimensional half-filled Hubbard model which allows us to study both antiferromagnetism and single-particle properties. This approach ignores amplitude fluctuations of the antiferromagnetic (AF) order parameter and is valid below a crossover temperature $T_X$ which marks the onset of AF short-range order. Directional fluctuations (spin waves) are described by a non-linear sigma model (NL$\sigma$M) that we derive from the Hubbard model. At zero temperature and weak coupling, our results are typical of a Slater antiferromagnet. The AF gap is exponentially small; there are well-defined Bogoliubov quasi-particles (QP's) (carrying most of the spectral weight) coexisting with a high-energy incoherent excitation background. As $U$ increases, the Slater antiferromagnet progressively becomes a Mott-Heisenberg antiferromagnet. The Bogoliubov bands evolve into Mott-Hubbard bands separated by a large AF gap. A significant fraction of spectral weight is transferred from the Bogoliubov QP's to incoherent excitations. At finite temperature, there is a metal-insulator transition between a pseudogap phase at weak coupling and a Mott-Hubbard insulator at strong coupling. Finally, we point out that our results straightforwardly translate to the half-filled attractive Hubbard model, where the ${\bf q}=(\pi,\pi)$ charge and ${\bf q}=0$ pairing fluctuations combine to form an order parameter with SO(3) symmetry.Comment: Revtex4, 19 pages, 14 figures; (v2) final version as publishe

A quantum Monte Carlo study of the one-dimensional ionic Hubbard model

Quantum Monte Carlo methods are used to study a quantum phase transition in a 1D Hubbard model with a staggered ionic potential (D). Using recently formulated methods, the electronic polarization and localization are determined directly from the correlated ground state wavefunction and compared to results of previous work using exact diagonalization and Hartree-Fock. We find that the model undergoes a thermodynamic transition from a band insulator (BI) to a broken-symmetry bond ordered (BO) phase as the ratio of U/D is increased. Since it is known that at D = 0 the usual Hubbard model is a Mott insulator (MI) with no long-range order, we have searched for a second transition to this state by (i) increasing U at fixed ionic potential (D) and (ii) decreasing D at fixed U. We find no transition from the BO to MI state, and we propose that the MI state in 1D is unstable to bond ordering under the addition of any finite ionic potential. In real 1D systems the symmetric MI phase is never stable and the transition is from a symmetric BI phase to a dimerized BO phase, with a metallic point at the transition

Surface heat flux feedback controlled impurity seeding experiments with Alcator C-Mod’s high-Z vertical target plate divertor: performance, limitations and implications for fusion power reactors

The Alcator C-Mod team has recently developed a feedback system to measure and control surface heat flux in real-time. The system uses real-time measurements of surface heat flux from surface thermocouples and a pulse-width modulated piezo valve to inject low-Z impurities (typically N2) into the private flux region. It has been used in C-Mod to mitigate peak surface heat fluxes >40 MW m−2 down to 1. While the system works quite well under relatively steady conditions, use of it during transients has revealed important limitations on feedback control of impurity seeding in conventional vertical target plate divertors. In some cases, the system is unable to avoid plasma reattachment to the divertor plate or the formation of a confinement damaging x-point MARFE. This is due to the small operational window for mitigated heat flux in the parameters of incident plasma heat flux, plasma density, and impurity density as well as the relatively slow response of the impurity gas injection system compared to plasma transients. Given the severe consequences for failure of such a system to operate reliably in a reactor, there is substantial risk that the conventional vertical target plate divertor will not provide an adequately controllable system in reactor-class devices. These considerations motivate the need to develop passively stable, highly compliant divertor configurations and experimental facilities that can test such possible solutions

Charge Asymmetries for D, D_s and Lambda_c Production in Sigma- - Nucleus Interactions at 340 GeV/c

CERN experiment WA89 has studied charmed particles produced by a Sigma^- beam at 340 GeV/c on nuclear targets. Production of particles which have light quarks in common with the beam is compared to production of those which do not. Considerable production asymmetries between D^- and D^p, D_s^ and D_s^+ and Lambda_c and Antilambda_c are observed. The results are compared with pion beam data and with theoretical calculations.Comment: LaTeX ; 16 pages including 4 ps figure