Cuprate parameters from numerical Wannier functions

We have developed a constrained-occupation local-density-functional technique to calculate individual model Hamiltonian parameters for states of arbitrary localization. The method generates numerical Wannier functions for the relevant states, and by decoupling these functions, permits direct calculation of the individual parameters. We have tested this method for the cuprate La₂{CuO}₄$, and obtained values for the Cu(3${d}ₓ²$-${y}²$) and O(2${p}_{{{σ}}}$) parameters in general agreement with previous empirical and theoretical estimates. These calculations provide even-handed treatment of both Cu(3d) and the more extended O(2p) states, and yield the important difference ε_p(σ)-εd(x²-y²) consistent with expectations based on the superexchange frequency. We also report parameter values for the Cu(3d3z²-r²) states, and for compatible O(2p) states, which may be of use in theoretical attempts to confirm the symmetry of the doped-in carriers.