Quantum Chemical Topology Study of the Water-Platinum(II) Interaction

Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions

We present ab initio simulations based on subsystem DFT of group 10 aqua ions accurately compared against experimental data on hydration structure. Our simulations provide insights into the molecular structures and dynamics of hydration shells, offering recalibrated interpretations of experimental results. We observe a soft, but distinct second hydration shell in Palladium (Pd) due to a balance between thermal fluctuations, metal-water interactions, and hydrogen bonding. Nickel (Ni) and platinum (Pt) exhibit more rigid hydration shells. Notably, our simulations align with experimental findings for Pd, showing axial hydration marked by a broad peak at about 3 Å in the Pd-O radial distribution function, revising the previously sharp "mesoshell" prediction. We introduce the "hydrogen bond dome" concept to describe a resilient network of hydrogen-bonded water molecules around the metal, which plays a critical role in the axial hydration dynamics.

Square-Planar Pt(II) and Ir(I) Complexes as the Lewis Bases: Donor-Acceptor Adducts with Group 13 Trihalides and Trihydrides

The stability and properties of donor-acceptor adducts of square-planar Pt(II) and Ir(I) complexes (designated as PtX, IrX, or generally MX complexes) with trihydrides and trihalides of group 13 elements of general formula YZ₃ (Y = B, Al, Ga; Z = H, F, Cl, Br) were studied theoretically using DFT methodology in the gas phase. MX complexes were represented by wide range of the ligand environment which included model complexes [Ir(NH₃)₃X]⁰ and cis-[Pt(NH₃)₂X₂]⁰ (X = H, CH₃, F, Cl, Br) and isoelectronic complexes [Ir(NNN)(CH₃)]⁰ and [Pt(NCN)(CH₃)]⁰ with tridentate NNN and NCN pincer ligands. MX complexes acted as the Lewis bases donating electron density from the doubly occupied 5d _z² atomic orbital of the metal M atom to the empty valence p _z orbital of Y whose evidence was clearly provided by the natural atomic orbital (NAO) analysis. This charge transfer led to the formation of pentacoordinated square pyramidal MX·(YZ₃) adducts with M·Y dative bond. Binding energies were -44.7 and -75.2 kcal/mol for interaction of GaF₃ as the strongest acid with PtNCN and IrNNN pincer ligands complexes. Only M·B bonds had covalent character although MX·BZ₃ adducts were the least stable due to large values of Pauli repulsion and deformation energies. The highest degree of covalent character was found for adducts of BH₃ in all series of structures studied. Al and Ga adducts showed remarkably similar behavior with respect to geometry and binding energies.

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.

On understanding the chemical origin of band gaps

Conceptual DFT and quantum chemical topology provide two different approaches based on the electron density to grasp chemical concepts. We present a model merging both approaches, in order to obtain physical properties from chemically meaningful fragments (bonds, lone pairs) in the solid. One way to do so is to use an energetic model that includes chemical quantities explicitly, so that the properties provided by conceptual DFT are directly related to the inherent organization of electrons within the regions provided by topological analysis. An example of such energy model is the bond charge model (BCM) by Parr and collaborators. Bonds within an ELF-BCM coupled approach present very stable chemical features, with a bond length of ca. 1 Å and 2[Formula: see text]. Whereas the 2[Formula: see text] corroborate classical views of chemical bonding, the fact that bonds always expand along 1 Å introduces the concept of geometrical transferability and enables estimating crystalline cell parameters. Moreover, combining these results with conceptual DFT enables deriving a model for the band gap where the chemical hardness of a solid is given by the bond properties, charge, length, and a Madelung factor, where the latter plays the major role. In short, the fundamental gap of zinc-blende solids can be understood as given by a 2[Formula: see text] bond particle asymmetrically located on a 1 Å length box and electrostatically interacting with other bonds and with a core matrix. This description is able to provide semi-quantitative insight into the band gap of zinc-blende semiconductors and insulators on equal footing, as well as a relationship between band gap and compressibility. In other words, merging these different approaches to bonding enables to connect measurable macroscopic behavior with microscopic electronic structure properties and to obtain microscopic insight into the chemical origin of band gaps, whose prediction is still nowadays a difficult task. Graphical Abstract Conceptual DFT couples to quatum chemcial topology to explain the band gap of zinc-blende solids.

Prediction of strong O-H/M hydrogen bonding between water and square-planar Ir and Rh complexes

Intermolecular O-H/M interactions, between a water molecule and square-planar acac complexes ([M(acac)L₂]), with different types of L ligands (en, H₂O, CO, CN^-, and OH^-) and different types of metal atoms (Ir(i), Rh(i), Pt(ii), and Pd(ii)) were studied by high level ab initio calculations. Among the studied neutral complexes, the [Pd(acac)(CN)(CO)] complex forms the weakest interaction, -0.62 kcal mol^-1, while the [Ir(acac)(en)] complex forms the strongest interaction, -9.83 kcal mol^-1, which is remarkably stronger than the conventional hydrogen bond between two water molecules (-4.84 kcal mol^-1).

Interactions of the "piano-stool" [ruthenium(II)(η(6) -arene)(quinolone)Cl](+) complexes with water; DFT computational study

Full optimizations of stationary points along the reaction coordinate for the hydration of several quinolone Ru(II) half-sandwich complexes were performed in water environment using the B3PW91/6-31+G(d)/PCM/UAKS method. The role of diffuse functions (especially on oxygen) was found crucial for correct geometries along the reaction coordinate. Single-point (SP) calculations were performed at the B3LYP/6-311++G(2df,2pd)/DPCM/saled-UAKS level. In the first part, two possible reaction mechanisms-associative and dissociative were compared. It was found that the dissociative mechanism of the hydration process is kinetically slightly preferred. Another important conclusion concerns the reaction channels. It was found that substitution of chloride ligand (abbreviated in the text as dechlorination reaction) represents energetically and kinetically the most feasible pathway. In the second part the same hydration reaction was explored for reactivity comparison of the Ru(II)-complexes with several derivatives of nalidixic acid: cinoxacin, ofloxacin, and (thio)nalidixic acid. The hydration process is about four orders of magnitude faster in a basic solution compared to neutral/acidic environment with cinoxacin and nalidixic acid as the most reactive complexes in the former and latter environments, respectively. The explored hydration reaction is in all cases endergonic; nevertheless the endergonicity is substantially lower (by ∼6 kcal/mol) in basic environment. © 2016 Wiley Periodicals, Inc.

3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One

Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.

On the nature of hydrogen bonds to platinum(II)--which interaction can predict their strength?

The interaction between hydrogen bond donors and platinum has been analysed. Our results point to an interaction that can be entirely predicted from the dz(2) orbital energy of the platinum centre indicating strong charge transfer, with significant dispersion contribution to the bonding, very different from classical hydrogen bonds.