QCDF90: Lattice QCD with Fortran 90

We have used Fortran 90 to implement lattice QCD. We have designed a set of machine independent modules that define fields (gauge, fermions, scalars, etc...) and overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write high-level efficient programs for QCD simulations. To increase performances our modules also implements assignments that do not require temporaries, and a machine independent precision definition. We have also created a useful compression procedure for storing the lattice configurations, and a parallel implementation of the random generators. We have widely tested our program and modules on several parallel and single processor supercomputers obtaining excellent performances.Comment: LaTeX file, 8 pages, no figures. More information available at: http://hep.bu.edu/~leviar/qcdf90.htm

Are People Ashamed of Paying with Food Stamps?

As is amply documented, there appears to be a large discrepancy between the marginal propensity to purchase food out of cash income and that out of food stamps. In this paper, we have examined both the formal and empirical bases of the claim that marginal welfare stigma explains this puzzling empirical regularity.SOCIAL WELFARE

Electronic structure and exchange interactions of the ladder vanadates CaV2O5 and MgV2O5

We have performed ab-initio calculations of the electronic structure and exchange couplings in the layered vanadates CaV2O5 and MgV2O5. Based on our results we provide a possible explanation of the unusual magnetic properties of these materials, in particular the large difference in the spin gap between CaV2O5 and MgV2O5

The Taming of QCD by Fortran 90

We implement lattice QCD using the Fortran 90 language. We have designed machine independent modules that define fields (gauge, fermions, scalars, etc...) and have defined overloaded operators for all possible operations between fields, matrices and numbers. With these modules it is very simple to write QCD programs. We have also created a useful compression standard for storing the lattice configurations, a parallel implementation of the random generators, an assignment that does not require temporaries, and a machine independent precision definition. We have tested our program on parallel and single processor supercomputers obtaining excellent performances.Comment: Talk presented at LATTICE96 (algorithms) 3 pages, no figures, LATEX file with ESPCRC2 style. More information available at: http://hep.bu.edu/~leviar/qcdf90.htm

Third-Generation TB-LMTO

We describe the screened Korringa-Kohn-Rostoker (KKR) method and the third-generation linear muffin-tin orbital (LMTO) method for solving the single-particle Schroedinger equation for a MT potential. The simple and popular formalism which previously resulted from the atomic-spheres approximation (ASA) now holds in general, that is, it includes downfolding and the combined correction. Downfolding to few-orbital, possibly short-ranged, low-energy, and possibly orthonormal Hamiltonians now works exceedingly well, as is demonstrated for a high-temperature superconductor. First-principles sp3 and sp3d5 TB Hamiltonians for the valence and lowest conduction bands of silicon are derived. Finally, we prove that the new method treats overlap of the potential wells correctly to leading order and we demonstrate how this can be exploited to get rid of the empty spheres in the diamond structure.Comment: latex2e, 32 printed pages, Postscript figs, to be published in: Tight-Binding Approach to Computational Materials Science, MRS Symposia Proceedings No. 491 (MRS, Pittsburgh, 1998

Effect of spin orbit coupling and Hubbard UU on the electronic structure of IrO2_2

We have studied in detail the electronic structure of IrO$_2$ including spin-orbit coupling (SOC) and electron-electron interaction, both within the GGA+U and GGA+DMFT approximations. Our calculations reveal that the Ir t$_{2g}$ states at the Fermi level largely retain the J$_{\rm eff}$ = $\frac{1}{2}$ character, suggesting that this complex spin-orbit entangled state may be robust even in metallic IrO$_2$. We have calculated the phase diagram for the ground state of IrO$_2$ as a function of $U$ and find a metal insulator transition that coincides with a magnetic phase change, where the effect of SOC is only to reduce the critical values of $U$ necessary for the transition. We also find that dynamic correlations, as given by the GGA+DMFT calculations, tend to suppress the spin-splitting, yielding a Pauli paramagnetic metal for moderate values of the Hubbard $U$. Our calculated optical spectra and photoemission spectra including SOC are in good agreement with experiment demonstrating the importance of SOC in IrO$_2$

NiS - An unusual self-doped, nearly compensated antiferromagnetic metal

NiS, exhibiting a text-book example of a first-order transition with many unusual properties at low temperatures, has been variously described in terms of conflicting descriptions of its ground state during the past several decades. We calculate these physical properties within first-principle approaches based on the density functional theory and conclusively establish that all experimental data can be understood in terms of a rather unusual ground state of NiS that is best described as a self-doped, nearly compensated, antiferromagnetic metal, resolving the age-old controversy. We trace the origin of this novel ground state to the specific details of the crystal structure, band dispersions and a sizable Coulomb interaction strength that is still sub-critical to drive the system in to an insulating state. We also show how the specific antiferromagnetic structure is a consequence of the less-discussed 90 degree and less than 90 degree superexchange interactions built in to such crystal structures