DFT calculation