General variational many-body theory with complete self-consistency for trapped bosonic systems

In this work we develop a complete variational many-body theory for a system of $N$ trapped bosons interacting via a general two-body potential. In this theory both the many-body basis functions {\em and} the respective expansion coefficients are treated as variational parameters. The optimal variational parameters are obtained {\em self-consistently} by solving a coupled system of non-eigenvalue -- generally integro-differential -- equations to get the one-particle functions and by diagonalizing the secular matrix problem to find the expansion coefficients. We call this theory multi-configurational Hartree for bosons or MCHB(M), where M specifies explicitly the number of one-particle functions used to construct the configurations. General rules for evaluating the matrix elements of one- and two-particle operators are derived and applied to construct the secular Hamiltonian matrix. We discuss properties of the derived equations. It is demonstrated that for any practical computation where the configurational space is restricted, the description of trapped bosonic systems strongly depends on the choice of the many-body basis set used, i.e., self-consistency is of great relevance. As illustrative examples we consider bosonic systems trapped in one- and two-dimensional symmetric and asymmetric double-well potentials. We demonstrate that self-consistency has great impact on the predicted physical properties of the ground and excited states and show that the lack of self-consistency may lead to physically wrong predictions. The convergence of the general MCHB(M) scheme with a growing number M is validated in a specific case of two bosons trapped in a symmetric double-well.Comment: 53 pages, 8 figure

Understanding the Heavy Fermion Phenomenology from Microscopic Model

We solve the 3D periodic Anderson model via two impurity DMFT. We obtain the temperature v.s. hybridization phase diagram. In approaching the quantum critical point (QCP) both the Neel and lattice Kondo temperatures decrease and they do not cross at the lowest temperature we reached. While strong ferromagnetic spin fluctuation on the Kondo side is observed, our result indicates the critical static spin susceptibility is local in space at the QCP. We observe in the crossover region logarithmic temperature dependence in the specific heat coefficient and spin susceptibility

Crossed spin-1/2 Heisenberg chains as a quantum impurity problem

Using equivalencies between different models we reduce the model of two spin-1/2 Heisenberg chains crossed at one point to the model of free fermions. The spin-spin correlation function is calculated by summing the perturbation series in the interchain coupling. The result reveals a power law decay with a nonuniversal exponent.Comment: 3 pages, the background information is adde

Lattice susceptibility for 2D Hubbard Model within dual fermion method

In this paper, we present details of the dual fermion (DF) method to study the non-local correction to single site DMFT. The DMFT two-particle Green's function is calculated using continuous time quantum monte carlo (CT-QMC) method. The momentum dependence of the vertex function is analyzed and its renormalization based on the Bethe-Salpeter equation is performed in particle-hole channel. We found a magnetic instability in both the dual and the lattice fermions. The lattice fermion susceptibility is calculated at finite temperature in this method and also in another recently proposed method, namely dynamical vertex approximation (D$\Gamma$A). The comparison between these two methods are presented in both weak and strong coupling region. Compared to the susceptibility from quantum monte carlo (QMC) simulation, both of them gave satisfied results.Comment: 10 pages, 11 figure

Neutrino Scattering in Heterogeneous Supernova Plasmas

Neutrinos in core collapse supernovae are likely trapped by neutrino-nucleus elastic scattering. Using molecular dynamics simulations, we calculate neutrino mean free paths and ion-ion correlation functions for heterogeneous plasmas. Mean free paths are systematically shorter in plasmas containing a mixture of ions compared to a plasma composed of a single ion species. This is because neutrinos can scatter from concentration fluctuations. The dynamical response function of a heterogeneous plasma is found to have an extra peak at low energies describing the diffusion of concentration fluctuations. Our exact molecular dynamics results for the static structure factor reduce to the Debye Huckel approximation, but only in the limit of very low momentum transfers.Comment: 11 pages, 13 figure

London's limit for the lattice superconductor

A stability problem for the current state of the strong coupling superconductor has been considered within the lattice Ginzburg-Landau model. The critical current problem for a thin superconductor film is solved within the London limit taking into account the crystal lattice symmetry. The current dependence on the order parameter modulus is computed for the superconductor film for various coupling parameter magnitudes. The field penetration problem is shown to be described in this case by the one-dimensional sine-Gordon equation. The field distribution around the vortex is described at the same time by the two-dimensional elliptic sine-Gordon equation.Comment: 7 pages, 3 figures, Revtex4, mostly technical correction; extended abstrac

Fragmented and Single Condensate Ground States of Spin-1 Bose Gas

We show that the ground state of a spin-1 Bose gas with an antiferro- magnetic interaction is a fragmented condensate in uniform magnetic fields. The number fluctuations in each spin component change rapidly from being enormous (order $N$) to exceedingly small (order 1) as the magnetization of the system increases. A fragmented condensate can be turned into a single condensate state by magnetic field gradients. The conditions for existence and the method of detecting fragmented states are presented.Comment: 4 pages, no figure

Magnetically Stabilized Nematic Order I: Three-Dimensional Bipartite Optical Lattices

We study magnetically stabilized nematic order for spin-one bosons in optical lattices. We show that the Zeeman field-driven quantum transitions between non-nematic Mott states and quantum spin nematic states in the weak hopping limit are in the universality class of the ferromagnetic XXZ (S=1/2) spin model. We further discuss these transitions as condensation of interacting magnons. The development of O(2) nematic order when external fields are applied corresponds to condensation of magnons, which breaks a U(1) symmetry. Microscopically, this results from a coherent superposition of two non-nematic states at each individual site. Nematic order and spin wave excitations around critical points are studied and critical behaviors are obtained in a dilute gas approximation. We also find that spin singlet states are unstable with respect to quadratic Zeeman effects and Ising nematic order appears in the presence of any finite quadratic Zeeman coupling. All discussions are carried out for states in three dimensional bipartite lattices.Comment: 16 pages, 3 figure

Quantum phase transition in an atomic Bose gas with a Feshbach resonance

We show that in an atomic Bose gas near a Feshbach resonance a quantum phase transition occurs between a phase with only a molecular Bose-Einstein condensate and a phase with both an atomic and a molecular Bose-Einstein condensate. We show that the transition is characterized by an Ising order parameter. We also determine the phase diagram of the gas as a function of magnetic field and temperature: the quantum critical point extends into a line of finite temperature Ising transitions.Comment: 4 pages, 2 figure