The formation energies of point defects and the interaction energies of various defect pairs in NiAl are calculated from first principles within an order N, locally self-consistent Green's-function method in conjunction with multipole electrostatic corrections to the atomic sphere approximation. The theory correctly reproduces the ground state for the off-stoichiometric NiAl alloys. The constitutional defects ͑antisite Ni atoms and Ni vacancies in Ni-rich and Al-rich NiAl, respectively͒ are shown to form ordered structures in the ground state, in which they tend to avoid each other at the shortest distance on their sublattice. The dominant thermal defects in Ni-rich and stoichiometric NiAl are calculated to be triple defects. In Al-rich alloys another type of thermal defect dominates, where two Ni vacancies are replaced by one antisite Al atom. As a result, the vacancy concentration decreases with temperature in this region. The effective defect formation enthalpies for different concentration regions of NiAl are also obtained.