Extensive ab inito computations have been carried out to study the equilibrium structure, infrared spectra, and bonding characteristics of a variety of hydrated NPO2(CO3)(m)(q-) complexes by considering the solvent as a polarizable dielectric continuum as well as the corresponding anhydrate complexes in the gas phase. The computed structural parameters and vibrational results at the MP2 level in aqueous solution are in good agreement with Clark et al.'s experiments and provide realistic pictures of the neptunyl complexes in an aqueous environment. Our computed hydration energies reveal that the complex with water molecules directly bound to it yields the best results. Our analysis of the nature of the bonding of neptunyl complexes provides insight into the nature of 6d and 5f bonding in actinide complexes.