A theoretical study to understand molecular structures and coordination modes of neptunium complexes with three different substituent hydroxamates viz., acetohydroxamate (AHA), benzhydroxamate (BHA) and salicylhydroxamate (SHA) was carried out using density functional theory (DFT). The geometries for different complexes reported in literature M x L y H z (x:y:z) viz 1:1:1, 1:2:2 and 1:2:0 were optimized and interaction energies (DE) for complexation process were calculated to obtain the most probable structures of complexes. The DE value revealed, the complex stability order as NpO 2 + -SHA > NpO 2 -AHA > NpO 2 -BHA which is in agreement with the reported experimental data. Higher negative charge on AHA compared to BHA leads to its higher interaction with NpO 2 which in turn reflects in shorter Np-O eq (ligand) bond distances. In 1:1:1 complex of NpO 2 -SHA, the binding mode involving ortho hydroxyl group and -NHO À oxygen atom (seven membered ring) shows most negative interaction energy.However, for binding of second ligand in 1:2:2 NpO 2 -SHA complex, steric effects dominate and the mode involving five membered ring through oxygen atoms of -CONHO À is most stable.