Mg(II) and Ca(II) Microsolvation by Ammonia: Born-Oppenheimer Molecular Dynamics Studies

We report the structural and energetic features of the Mg²⁺ and Ca²⁺ cations in ammonia microsolvation environments. Born-Oppenhemier molecular dynamics studies are carried out for [Mg(NH₃)]²⁺ and [Ca(NH₃)]²⁺ clusters with n = 2, 3, 4, 6, 8, 20, and 27 at 300 K based on hybrid density functional theory calculations. We determine binding energies per ammonia molecule and the metal cation solvation patterns as a function of the number of molecules. The general trend for Mg²⁺ is that the Mg-N distances increase as a function of n until the first solvation shell is populated by six ammonia molecules, and then the distances slightly decrease while CN = 6 does not change. For Ca²⁺, the first solvation shell at room temperature is populated by eight ammonia molecules for clusters with more than one solvation shell, leading to a different structure from that of [Ca(NH₃)₆]²⁺ hexamine. The evaporation of NH₃ molecules was found at 300 K only for Mg²⁺ clusters with n ≥ 10; this was not the case for Ca²⁺ clusters. Vibrational spectra are obtained for all of the clusters, and the evolution of the main features is discussed. EXAFS spectra are also presented for the [Mg(NH₃)₂₇(NH₃)₂₇]²⁺ and [Ca(NH₃)₂₇]²⁺ clusters, which yield valuable data to be compared with experimental data in the liquid phase, as previously done for the aqueous solvation of these dications.