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

Synergy Effects in Heavy Metal Ion Chelation with Aryl- and Aroyl-Substituted Thiourea Derivatives

Detection and removal of metal ion contaminants have attracted great interest due to the health risks that they represent for humans and wildlife. Among the proposed compounds developed for these purposes, thiourea derivatives have been shown as quite efficient chelating agents of metal cations and have been proposed for heavy metal ion removal and for components of high-selectivity sensors. Understanding the nature of metal-ionophore activity for these compounds is thus of high relevance. We present a theoretical study on the interaction between substituted thioureas and metal cations, namely, Cd²⁺, Hg²⁺, and Pb²⁺. Two substituent groups have been chosen: 2-furoyl and m-trifluoromethylphenyl. Combining density functional theory simulations with wave function analysis techniques, we study the nature of the metal-thiourea interaction and characterize the bonding properties. Here, it is shown how the N,N'-disubstituted derivative has a strong affinity for Hg²⁺, through cation-hydrogen interactions, due to its greater oxidizing capacity.

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(H₂O)]²⁺ 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(H₂O)²⁺ 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(H₂O)₂]²⁺ 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 Hg²⁺ 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(H₂O)₃]²⁺ 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 Cd²⁺. 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(H₂O)₃]²⁺ 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.

An efficient ab initio DFT and PCM assessment of the potentiometric selectivity of a salophen type Schiff base

As a neutral carrier component for the preparation of a potentiometric membrane sensor, the affinity and selectivity of the salophen type Schiff base ligand obtained by 1:2 condensation of 2.3-diaminopyridine with salicylaldehyde toward a series of common cations has been fully examined by DFT/B3LYP and integral equation formalism polarizable continum model (IEF-PCM or only given with PCM as default input in the computations) in combination with the experimental data. Both the potentiometric measurements and DFT calculations have exhibited that the ionophore shows appreciable selectivity for Cu(2+) ion over other cations. Four different approaches where the last three are the modified version of each other have been evaluated and compared with potentiometric data. Based upon the results of comparison among the approaches suggested to verify the selective behavior of ionophore toward Cu(2+), PCM implemented approach having a whole computational groundwork has given well-matched results with the observed data and with the method augmented with experimental hydration energies. The foremost interferences were detected by determining potentiometric selectivity coefficients for each metal ion relative to Cu(2+) and compared to the results obtained by the DFT calculations.

Theoretical strategy to build structural models of microhydrated inorganic systems for the knowledge of their vibrational properties: the case of the hydrated nitrate aerosols

This study provides theoretical anharmonic calculations for microhydrated NaNO3-labeled (NaNO3, nH2O)x with a water-to-solute ratio (n) ranging from 1 to 3. A representative geometrical model of these forms was first investigated by simulating the molecular clusters as (NaNO3,1H2O)x with x = 1 to 4. The comparison between the calculated time independent anharmonic frequencies using the B3LYP-D/6-311+G(d,p) method and their experimental counterparts led to the choice of a supercluster model. The most probable structures of (NaNO3,nH2O)3 molecular system were investigated by using our global search algorithm we developed recently (GSAM code) both at the B3LYP/6-311+G(d,p) and the B3LYP-D/6-311+G(d,p) levels of theory. The quality of the structural model is illustrated by comparing the B3LYP/6-311+G(d,p) and B3LYP-D/6-311+G(d,p) anharmonic vibrational signatures with those obtained from IR experiments. While an average deviation of 16 cm(-1) is observed in the case of the B3LYP computations, the deviation is reduced to 7 cm(-1) for the B3LYP-D computations.

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.