Coulomb parameters and photoemission for the molecular metal TTF-TCNQ

We employ density-functional theory to calculate realistic parameters for an extended Hubbard model of the molecular metal TTF-TCNQ. Considering both intra- and intermolecular screening in the crystal, we find significant longer-range Coulomb interactions along the molecular stacks, as well as inter-stack coupling. We show that the long-range Coulomb term of the extended Hubbard model leads to a broadening of the spectral density, likely resolving the problems with the interpretation of photoemission experiments using a simple Hubbard model only.Comment: 4 pages, 2 figure

Massively parallel exact diagonalization of strongly correlated systems

The physics of strongly correlated materials poses one of the most challenging problems in condensed-matter sciences. Standard approximations applicable to wide classes of materials such as the local density approximation fail, due to the importance of the Coulomb repulsion between localized electrons. Instead, we resort to non-perturbative many-body methods. The calculations are, however, only feasible for rather small model systems. The full Hamiltonian of a real material is approximated by a model Hamiltonian comprising only the most important electronic degrees of freedom, while the effect of all other electrons is included in an average way by renormalizing the parameters. Realistic calculations of strongly correlated materials need to include sufficiently many of these electronic degrees of freedom. The new generation of massively parallel supercomputers allows for these realistic calculations. However, exploiting their computational power requires newly devised algorithms. As a solver we use the Lanczos method, which needs the full many-body state of the correlated system. It is thus limited by the available main memory. The foremost problem for a distributed-memory implementation is that the multiplication of the Hamiltonian to the many-body state leads to highly non-local memory access patterns. A simple yet important observation leads to an efficient solution: in the kinetic term of the Hamiltonian the electron-spin is conserved. Thus, writing the many-body vector as a matrix v(iup,idn), where the indices label spin-configurations, we find that the hopping term only connects vector elements that differ in one index. Hence, storing entire slices v(iup,:) on one node, the kinetic term for the spin-down electrons is local to that process. After transposing v, the same is true for the hopping of the spin-up electrons...