Exploring structure and dynamics of solvated Ca(II) in liquid ammonia: A quantum mechanical charge field (QMCF) molecular dynamics simulation

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₃)_n]²⁺ and [Ca(NH₃)_n]²⁺ 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.

Cu2+ in liquid ammonia-The impact of solvent flexibility and electron correlation in ab initio quantum mechanical charge field molecular dynamics

The impact of solvent flexibility and electron correlation on the simulation results of Cu²⁺ in liquid ammonia has been investigated via an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach. To achieve this, three different simulation systems were considered in this study, namely Cu²⁺ in rigid and flexible ammonia at Hartree-Fock (HF) level of theory, as well as resolution of identity second order Møller-Plesset (MP2) perturbation theory in the rigid body case. In all cases, a stable octahedral [Cu(NH₃ )₆ ]²⁺ complex subject to dynamic Jahn-Teller distortions without the occurrence of ligand exchange was observed. The Cu²⁺  - NH₃ distance in the first shell agrees well with the experimental and other theoretical data. In all three cases, the structural data shows that the rigid-body ammonia model in conjunction with the HF level of theory provides accurate data for the first solvation shell, while at the same time, the computational demand and thus the achievable simulation time are much more beneficial. The vibrational analysis of the Cu²⁺  - NH₃ interaction yields similar force constants in the three investigated systems indicating that there is no distinct difference on the dynamical properties of the first solvation shell. In addition to the QMCF MD simulations, a number of natural bond orbital (NBO) analyses were carried out, confirming the strong electrostatic character of the Cu²⁺  - NH₃ interaction.

The Jahn-Teller effect in mixed aqueous solution: the solvation of Cu2+ in 18.6% aqueous ammonia obtained from ab initio quantum mechanical charge field molecular dynamics

The solvation structure and dynamics of Cu2+ in 18.6 % aqueous ammonia have been investigated using an ab initio quantum mechanical charge field molecular dynamics (QMCF MD) simulation approach at the Hartree–Fock (HF) level of theory applying the LANL2DZ ECP and Dunning DZP basis sets for Cu2+, ammonia and water, respectively. During a simulation time of 20 ps, only NH3 molecules are observed within the first solvation shell of Cu2+, resulting in the formation of an octahedral [Cu(NH3)6]2+ complex. While no exchange of these ligands with the second solvation shell are observed along the simulation, the monitoring of the associated N-Ntrans distances highlight the dynamics of the associated Jahn-Teller distortions, showing on average 2 elongated axial (2.19 Å) and 4 equatorial Cu–N bonds (2.39 Å). The observed structural properties are found in excellent agreement with experimental studies. In addition, an NBO analysis was carried out, confirming the strong electrostatic character of the Cu2+–NH3 interaction.

Square planar or octahedral after all? The indistinct solvation of platinum(ii)

The solvation structures of Pd(ii) and Pt(ii) are typically reduced to the well-known square-planar structural motif, although it has been shown, in both experimental and theoretical investigations, that these solutes demonstrate the affinity to bind ligand molecules at elongated distance in axial coordination sites.