We present release 1.3 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an international open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. Development is centered on common XML and binary data formats, on libraries to simplify and speed up code development, and on full-featured simulation programs. The programs enable non-experts to start carrying out numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), as well as the density matrix renormalization group (DMRG). Changes in the new release include a DMRG program for interacting models, support for translation symmetries in the diagonalization programs, the ability to define custom measurement operators, and support for inhomogeneous systems, such as lattice models with traps. The software is available from our web server at http://alps.comp-phys.org/.

We report high-resolution high-energy photoemission spectra together with parameter-free LDA + DMFT (local density approximation + dynamical mean-field theory) results for Sr1-xCaxVO3, a prototype 3d(1) system. In contrast to earlier investigations the bulk spectra are found to be insensitive to x. The good agreement between experiment and theory confirms the bulk sensitivity of the high-energy photoemission spectra.

Ab initio analyses of A(2)IrO(4) (A=Sr,Ba) are presented. Effective Hubbard-type models for Ir 5d t(2g) manifolds downfolded from the global band structure are solved based on the dynamical mean-field theory. The results for A=Sr and Ba correctly reproduce paramagnetic metals undergoing continuous transitions to insulators below the Néel temperature T(N). These compounds are classified not into Mott insulators but into Slater insulators. However, the insulating gap opens by a synergy of the Néel order and significant band renormalization, which is also manifested by a 2D bad metallic behavior in the paramagnetic phase near the quantum criticality.

We propose a computational scheme for the ab initio calculation of Wannier functions (WFs) for correlated electronic materials. The full-orbital HamiltonianĤ is projected into the WF subspace defined by the physically most relevant partially filled bands. The HamiltonianĤ W F obtained in this way, with interaction parameters calculated by constrained LDA for the Wannier orbitals, is used as an ab initio setup of the correlation problem, which can then be solved by many-body techniques, e.g., dynamical mean-field theory (DMFT). In such calculations the self-energy operator Σ(ε) is defined in WF basis which then can be converted back into the full-orbital Hilbert space to compute the full-orbital interacting Green function G(r, r ′ , ε). Using G(r, r ′ , ε) one can evaluate the charge density, modified by correlations, together with a new set of WFs, thus defining a fully selfconsistent scheme. The Green function can also be used for the calculation of spectral, magnetic and electronic properties of the system. Here we report the results obtained with this method for SrVO 3 and V 2 O 3 . Comparisons are made with previous results obtained by the LDA+DMFT approach where the LDA DOS was used as input, and with new bulk-sensitive experimental spectra.

We present parameter-free LDA+DMFT (local density approximation + dynamical mean field theory) results for the many-body spectra of cubic SrVO 3 and orthorhombic CaVO 3 . Both systems are found to be strongly correlated metals, but not on the verge of a metal-insulator transition.In spite of the considerably smaller V-O-V bond angle in CaVO 3 the LDA+DMFT spectra of the two systems for energies E < E F are very similar, their quasiparticle parts being almost identical.The calculated spectrum for E > E F shows more pronounced, albeit still small, differences. This is in contrast to earlier theoretical and experimental conclusions, but in good agreement with recent bulk-sensitive photoemission and x-ray absorption experiments.

LDA+DMFT (Local Density Approximation combined with Dynamical Mean-Field Theory) computation scheme has been used to calculate spectral properties of LaFeAsO -the parent compound for new high-T c iron oxypnictides. Coulomb repulsion U and Hund's exchange J parameters for iron 3d electrons were calculated using first principles constrained density functional theory scheme in Wannier functions formalism. Resulting values strongly depend on the number of states taken into account in calculations: when full set of O-2p, As-4p, and Fe-3d orbitals with corresponding bands are included, computation results in U =3÷4 eV and J=0.8 eV. In contrast to that when the basis set is restricted to Fe-3d orbitals and bands only, computation gives much smaller parameter values F 0 =0.8 eV, J=0.5 eV. However, DMFT calculations with both parameter sets and corresponding to them choice of basis functions result in weakly correlated electronic structure that is in agreement with experimental X-ray and photoemission spectra.

Electronic correlations in the Fe-pnictide BaFe2As2 are explored within LDA+DMFT, the combination of density functional theory with dynamical mean-field theory. While the correlated band structure is substantially renormalized there is only little transfer of spectral weight. The computed k-integrated and k-resolved spectral functions are in good agreement with photoemission spectroscopy (PES) and angular resolved PES experiments. Making use of a general classification scheme for the strength of electronic correlations we conclude that BaFe2As2 is a moderately correlated system.

Starting from a single band Hubbard model in the Wannier function basis, we revisit the problem of the ligand contribution to exchange and derive explicit formulae for the exchange integrals in metal oxide compounds in terms of atomic parameters that can be calculated with constrained LDA and LDA+U. The analysis is applied to the investigation of the isotropic exchange interactions of LiCu2O2, a compound where the Cu-O-Cu angle of the dominant exchange path is close to 90 • . Our results show that the magnetic moments are localized in Wannier orbitals which have strong contribution from oxygen atomic orbitals, leading to exchange integrals that considerably differ from the estimates based on kinetic exchange only. Using LSDA+U approach, we also perform a direct ab-initio determination of the exchange integrals LiCu2O2. The results agree well with those obtained from the Wannier function approach, a clear indication that this modelization captures the essential physics of exchange. A comparison with experimental results is also included, with the conclusion that a very precise determination of the Wannier function is crucial to reach quantitative estimates.