Single hole dynamics in the Kondo Necklace and Bilayer Heisenberg models on a square lattice

We study single hole dynamics in the bilayer Heisenberg and Kondo Necklace models. Those models exhibit a magnetic order-disorder quantum phase transition as a function of the interlayer coupling J_perp. At strong coupling in the disordered phase, both models have a single-hole dispersion relation with band maximum at p = (\pi,\pi) and an effective mass at this p-point which scales as the hopping matrix element t. In the Kondo Necklace model, we show that the effective mass at p = (\pi,\pi) remains finite for all considered values of J_perp such that the strong coupling features of the dispersion relation are apparent down to weak coupling. In contrast, in the bilayer Heisenberg model, the effective mass diverges at a finite value of J_perp. This divergence of the effective mass is unrelated to the magnetic quantum phase transition and at weak coupling the dispersion relation maps onto that of a single hole doped in a planar antiferromagnet with band maximum at p = (\pi/2,\pi/2). We equally study the behavior of the quasiparticle residue in the vicinity of the magnetic quantum phase transition both for a mobile and static hole. In contrast to analytical approaches, our numerical results do not unambiguously support the fact that the quasiparticle residue of the static hole vanishes in the vicinity of the critical point. The above results are obtained with a generalized version of the loop algorithm to include single hole dynamics on lattice sizes up to 20 X 20.Comment: 12 pages, 13 Fig

Comment on "Quantum Monte Carlo Evidence for Superconductivity in the Three-Band Hubbard Model in Two Dimensions"

In a recent Letter, Kuroki and Aoki [Phys. Rev. Lett. 76, 440 (1996)] presented quantum Monte-Carlo (QMC) results for pairing correlations in the three-band Hubbard model, which describes the Cu-d_{x^2-y^2} and O-p_{x,y} orbitals present in the CuO_2 planes of high-T_c materials. In this comment we argue that (i) the used parameter set is not appropriate for the description of high-T_c materials since it does not satisfy the minimal requirement of a charge-transfer gap at half-filling, and (ii) the observed increase in the d_{x^2-y^2} channel is dominantly produced by the pair-field correlations without the vertex part. Hence, the claim of evidence of ODLRO is not justified.Comment: 1 page latex and 2 eps-figures, uses epsfig, submitted to PR

Dynamic response of trapped ultracold bosons on optical lattices

We study the dynamic response of ultracold bosons trapped in one-dimensional optical lattices using Quantum Monte Carlo simulations of the boson Hubbard model with a confining potential. The dynamic structure factor reveals the inhomogeneous nature of the low temperature state, which contains coexisting Mott insulator and superfluid regions. We present new evidence for local quantum criticality and shed new light on the experimental excitation spectrum of 87Rb atoms confined in one dimension.Comment: 4 pages, 5 figure

Impact of Climate Change on the Production of Wheat and Rice in India

India is a world agricultural power, the share of agriculture in GDP is 17.8%, agricultural employment represents 54.87% of total employment and the country is considered the second largest producer wheat and rice. But this advantage is facing several obstacles because India is likely to be severely affected by climate change. India is one of the country's most vulnerable to disasters around the world and many of its 1.2 billion people live in areas vulnerable to hazards such as floods, cyclones and droughts. In this work, we propose a model that highlights the impacts of climate changes (changes in temperatures and precipitations) on the production of wheat and rice. The results confirm the hypothesis according which the impact of climate changes in India are Important

Dynamic Exponent of t-J and t-J-W Model

Drude weight of optical conductivity is calculated at zero temperature by exact diagonalization for the two-dimensional t-J model with the two-particle term, $W$. For the ordinary t-J model with $W$=0, the scaling of the Drude weight $D \propto \delta^2$ for small doping concentration $\delta$ is obtained, which indicates anomalous dynamic exponent $z$=4 of the Mott transition. When $W$ is switched on, the dynamic exponent recovers its conventional value $z$=2. This corresponds to an incoherent-to-coherent transition associated with the switching of the two-particle transfer.Comment: LaTeX, JPSJ-style, 4 pages, 5 eps files, to appear in J. Phys. Soc. Jpn. vol.67, No.6 (1998

Doping induced metal-insulator transition in two-dimensional Hubbard, t−Ut-U, and extended Hubbard, t−U−Wt-U-W, models

We show numerically that the nature of the doping induced metal-insulator transition in the two-dimensional Hubbard model is radically altered by the inclusion of a term, $W$, which depends upon a square of a single-particle nearest-neighbor hopping. This result is reached by computing the localization length, $\xi_l$, in the insulating state. At finite values of $W$ we find results consistent with $\xi_l \sim | \mu - \mu_c|^{- 1/2}$ where $\mu_c$ is the critical chemical potential. In contrast, $\xi_l \sim | \mu - \mu_c|^{-1/4}$ for the Hubbard model. At finite values of $W$, the presented numerical results imply that doping the antiferromagnetic Mott insulator leads to a $d_{x^2 - y ^2}$ superconductor.Comment: 19 pages (latex) including 7 figures in encapsulated postscript format. Submitted for publication in Phys. Rev.

Insulator-Metal Transition in the One and Two-Dimensional Hubbard Models

We use Quantum Monte Carlo methods to determine $T=0$ Green functions, $G(\vec{r}, \omega)$, on lattices up to $16 \times 16$ for the 2D Hubbard model at $U/t =4$. For chemical potentials, $\mu$, within the Hubbard gap, $|\mu | < \mu_c$, and at {\it long} distances, $\vec{r}$, $G(\vec{r}, \omega = \mu) \sim e^{ -|\vec{r}|/\xi_l}$ with critical behavior: $\xi_l \sim | \mu - \mu_c |^{-\nu}$, $\nu = 0.26 \pm 0.05$. This result stands in agreement with the assumption of hyperscaling with correlation exponent $\nu = 1/4$ and dynamical exponent $z = 4$. In contrast, the generic band insulator as well as the metal-insulator transition in the 1D Hubbard model are characterized by $\nu = 1/2$ and $z = 2$.Comment: 9 pages (latex) and 5 postscript figures. Submitted for publication in Phys. Rev. Let

Quantum Monte Carlo study of confined fermions in one-dimensional optical lattices

Using quantum Monte Carlo (QMC) simulations we study the ground-state properties of the one-dimensional fermionic Hubbard model in traps with an underlying lattice. Since due to the confining potential the density is space dependent, Mott-insulating domains always coexist with metallic regions, such that global quantities are not appropriate to describe the system. We define a local compressibility that characterizes the Mott-insulating regions and analyze other local quantities. It is shown that the momentum distribution function, a quantity that is commonly considered in experiments, fails in giving a clear signal of the Mott-insulator transition. Furthermore, we analyze a mean-field approach to these systems and compare it with the numerically exact QMC results. Finally, we determine a generic form for the phase diagram that allows us to predict the phases to be observed in the experiments.Comment: RevTex file, 13 pages, 19 figures, published versio

Enhancement of Pairing Correlation and Spin Gap through Suppression of Single-Particle Dispersion in One-Dimensional Models

We investigate the effects of suppression of single-particle dispersion near the Fermi level on the spin gap and the singlet-pairing correlation by using the exact diagonalization method for finite-size systems. We consider strongly correlated one-dimensional models, which have flat band dispersions near wave number k=\pi/2, if the interactions are switched off. Our results for strongly correlated models show that the spin gap region expands as the single-particle dispersion becomes flatter. The region where the singlet-pairing correlation is the most dominant also expands in models with flatter band dispersions. Based on our numerical results, we propose a pairing mechanism induced by the flat-band dispersion.Comment: 5 pages, including 5 eps figures, to appear in J.Phys.Soc.Jpn Vol.69 No.