DFT analysis of the interaction between Hg2+ and monodentate neutral ligands using NBO, EDA, and QTAIM.

We report thermodynamic, geometric, and electronic parameters for the interaction between neutral ligands and the [Hg(H2O)]2+ dication, using the B3LYP/6-311 + G(d,p) approach. Gibbs free energies for the interaction were employed to rank the affinity order of the several neutral ligands. To identify the parameters that characterize the affinity between the two fragments, the metal-ligand interaction was analyzed according to the EDA, NBO, and QTAIM decomposition schemes. The phosphine oxide showed the highest affinity for the Hg(H2O)2+ dication, mainly due to the P=O bond polarization. Ligands containing the sulfur atom, characterized by a high covalent component for the metal-ligand interaction, are the following in the interaction order. According to the Gibbs free energy for substitution of one water molecule in the [Hg(H2O)2]2+ complex, the sequence for the affinity order is: phosphine oxide > thioketone > thioesther > lactam > amide > amine > carboxylic acid > thiophene > ketone > esther > thiol > thiocyanate > ammonia > disulfide > aldehyde > ether > haloydrin > alcohol > enol > azide. Graphical abstract Synopsis The interaction between the Hg2+ cation and monodentate ligands containing S, O, or N atoms was evaluated in terms of energetic (bond strength, electrostatic and covalent interactions, donation energy), geometric (metal-ligand distance), electronic (atomic charges, orbital overlap, orbital hybridization) and topologic parameters.

A DFT study of the interaction between [Cd(H2O)3]2+ and monodentate O-, N-, and S-donor ligands: bond interaction analysis.

A series of B3LYP/6-311+G(d,p) calculations of the affinity of monodentate ligands for [Cd(H2O)3]2+ are performed. Three types of ligands containing O (phosphine oxide, lactam, amide, carboxylic acid, ester, ketone, aldehyde, ether, halohydrin, enol, furan), N (thiocyanate, amine, ammonia, azide), and S (thioester, thioketone, thiol, thiophene, disulfide) interacting atoms are investigated. The results show that phosphine oxide has the largest affinity for the cadmium cation due to the polarization of the P=O bond. As the P atom has a large atomic radius, the O atom can polarize the electronic cloud enhancing its amount of electronic charge and favoring the interaction with Cd2+. The affinity order found is phosphine oxide > thioester > lactam > amide > carboxylic acid > ester > thioketone > ketone > thiocyanate > amine > ammonia > aldehyde > ether > thiol > thiophene > enol > halohydrin > disulfide > azide > furan ligands. These results were also corroborated by the functional M06-2X. The electronic effects (resonance and induction) of neighboring groups of the interacting atom modulate the strength of metal-ligand binding. For almost all the O-donor ligands the electrostatic component has the same magnitude as the covalent term, while for the N- and S-donor ligands the covalent term is predominant. The polarization term accounts for twice the exchange term as part of the covalent component. The dispersion term varies less than 2 kcal mol-1 for the complexes analyzed. The Pauli repulsion term is correlated with the metal ligand distance, increasing in the compounds with decreased metal-ligand bond length. The charge between the interacting atoms is also strongly correlated with both the interacting strength and the electrostatic interaction component. The natural bond orbital analysis highlights correlations of the bond order, and S and P contributions of the interacting metal-ligand orbital with the coordination strength. Graphical abstract The affinity of 20 monodentate ligands with different functional groups for the [Cd(H2O)3]2+ cation is calculated based on the interaction enthalpy and Gibbs free energy for the substitution of one water molecule from the fully hydrated cation. The affinity is correlated with geometric, electronic, and energetic parameters of the ligands and the complexes as well as with energy decomposition and natural bond order analyses results.

A density functional theory investigation of the interaction of the tetraaqua calcium cation with bidentate carbonyl ligands.

Calcium complexes with bidentate carbonyl ligands are important in biological systems, medicine and industry, where the concentration of Ca2+ is controlled using chelating ligands. The exchange of two water molecules of [Ca(H2O)6]2+ for one bidentate monosubstituted and homo disubstituted dicarbonyl ligand was investigated using the B3LYP/6-311++G(d,p) method. The ligand substituents NH2, OCH3, OH, CH3, H, F, Cl, CN and NO2 are functional groups with distinct electron-donating and -withdrawing effects that bond directly to the sp2 C atom of the carbonyl group. The geometry, charge and energy characteristics of the complexes were analyzed to help understand the effects of substituents, spacer length and chelation. Coordination strength was quantified in terms of the enthalpy and free energy of the exchange reaction. The most negative enthalpies were calculated for the coordination of bidentate ligands containing three to five methylene group spacers between carbonyls. The chelate effect contribution was analyzed based on the thermochemistry. The electronic character of the substituent modulates the strength of binding to the metal cation, as ligands containing electron-donor substituents coordinate stronger than those with electron-acceptor substituents. This is reflected in the geometric (bond length and chelating angle), electronic (atomic charges) and energetic (components of the total interacting energy) characteristics of the complexes. Energy decomposition analysis (EDA)-an approach for partitioning of the energy into its chemical origins-shows that the electrostatic component of the coordination is predominant, and yields relevant contribution of the covalent term, especially for the electron-withdrawing substituted ligands. The chelate effect of the bidentate ligands was noticeable when compared with substitution by two monodentate ligands. Graphical abstract The affinity of 18 bidentate carbonyl ligands toward the [Ca(H2O)4]2+ cation is evaluated in terms of energetic, geometric and electronic parameters of the isolated ligands and the substituted aqua complexes. The electronic effects-inductive and mesomeric-intrinsic to the molecular structure of each ligand are found to modulate the strength of the metal-ligand interaction. The effects of polysubstitution, chelation and the length of the alkyl spacers between the anchor points of the ligand are also analyzed.

DFT studies of imino and thiocarbonyl ligands with the pentaaqua Mg²âº cation: affinity and associated parameters.

The affinity of the pentaaqua Mg(2+) cation for a set of para-substituted imino [HN = CHC6H4(R)] and thiocarbonyl [S = CHC6H4(R)] ligands (R = H, F, Cl, Br, OH, OCH3, CH3, CN, NH2 and NO2) was analyzed with DFT (B3LYP/6-31+G(d)) and semi-empirical (PM6-DH2) methods. The interaction enthalpy was calculated to quantify the affinity of the Mg(2+) cation for the ligands. Additionally, geometric and electronic parameters were correlated with the intensity of the metal-ligand interaction. The imino ligands have stronger interaction with the pentaaqua Mg(2+) cation than the thiocarbonyl derivatives. The electronic nature of the substituent is the main parameter that determines the interaction enthalpy. Ligands with electron donor substituents have more exothermic interaction enthalpies than those with electron withdrawing groups. The HSBA analysis showed that the interaction between the Mg(2+) cation and hard bases (imino ligands) is stronger than with soft bases (thiocarbonyl derivatives). The EDA analysis showed that the electrostatic, covalent and repulsion components of the interaction are the most affected by the substituent, whereas the dispersion and exchange components are almost constant.

Interaction between alkaline earth cations and oxo-ligands. DFT study of the affinity of the Ca2+ cation for carbonyl ligands.

The affinity of the Ca(2+) ion for a set of substituted carbonyl ligands was analyzed with both the DFT (B3LYP/6-31+G(d)) and semi-empirical (PM6) methods. Two types of ligands were studied: a set of monosubstituted [O=CH(R)] and a set of disubstituted ligands [O=C(R)(2)] (R=H, F, Cl, Br, OH, OCH(3), CH(3), CN, NH(2) and NO(2)), with R either directly bound to the carbonyl carbon atom or to the para position of a phenyl ring. The interaction energy was calculated to quantify the affinity of the Ca(2+) cation for the ligands. Geometric and electronic parameters were correlated with the intensity of the metal-ligand interaction. The electronic nature of the substituent is the main parameter that determines the interaction energy. Donor groups make the interaction energy more negative (stabilizing the complex formed), while acceptor groups make the interaction energy less negative (destabilizing the complex formed).

Interactions between alkaline earth cations and oxo ligands. DFT study of the affinity of the Mg²+ cation for phosphoryl ligands.

DFT (B3LYP/6-31+G(d)) calculations of Mg(2+) affinities for a set of phosphoryl ligands were performed. Two types of ligands were studied: a set of trivalent [O = P(R)] and a set of pentavalent phosphoryl ligands [O = P(R)(3)] (R = H, F, Cl, Br, OH, OCH(3), CH(3), CN, NH(2) and NO(2)), with R either bound directly to the phosphorus atom or to the para position of a phenyl ring. The affinity of the Mg(2+) cation for the ligands was quantified by means of the enthalpy for the substitution of one water molecule in the [Mg(H(2)O)(6)](2+) complex for a ligand. The enthalpy of substitution was correlated with electronic and geometric parameters. Electron-donor groups increase the interaction between the cation and the ligand, while electron-acceptor groups decrease the interaction enthalpy.