Solvation of transition metal complexes: thermochemical approaches

Conclusion and extensions We hope that this Review has made readers more aware of solvation of inorganic complexes, and of the importance of such knowledge in understanding their chemistryperhaps particularly their reactivity. The approach just set out for inorganic complexes should be of considerable value in the field or organometallic chemistry. In particular, informed use of solvation characteristics should help in optimising conditions for organometallic reactions and in homogeneous catalysis. Unfortunately, solvation data on reactants are too sparse (the subject index ofComprehensive Organometallic Chemistry contains justthree entries under solubility!) for serious examination of reactivity trends in terms of initial state and transition state contributions to be possible in almost all areas. Moreover, there are some fundamental problems over transfer parameters. Thus, a favourite electrochemical assumption is that the ferrocene/ferrocinium redox potential is independent of solvent. Yet, the dependence of rate constants on medium for outer-sphere electron transfer in the ferrocene/ferrocinium system can only be understood⁽⁶⁶⁾ in terms of specific solvation effects which are incompatible with the parallel solvation changes of these two substrates implicit in the redox potential assumption. The solvation of organometallic species should prove a most rewarding area for continued study, but it will be some time before the overall picture becomes as clear as in the more limited area of classical transition metal complexes considered in the present Review.     

Solvation structure and transport of acidic protons in ionic liquids: a first-principles simulation study

Ab initio simulations of a single molecule of HCl in liquid dimethyl imidazolium chloride [dmim][Cl] show that the acidic proton exists as a symmetric, linear ClHCl(-) species. Details of the solvation structure around this molecule are given. The proton-transfer process was investigated by applying a force along the antisymmetric stretch coordinate until the molecule broke. Changes in the free energy and local solvation structure during this process were investigated. In the reaction mechanism identified, a free chloride approaches the proton from the side. As the original ClHCl(-) distorts and the incoming chloride forms a new bond to the proton, one of the original chlorine atoms is expelled and a new linear molecule is formed.     

On the structure of the benzene-chlorine complex: a CNDO study

An oblique geometry in which the chlorine molecule is perpendicular to the benzene ring and centered over a carbon-carbon bond has been found to have a calculated CNDO binding energy more stable than those previously obtained for other geometrical models. This geometry is equivalent to the one preferred on the basis of infrared spectral studies of matrix-isolated benzene-chlorine complexes.     

Solvation states of HCl in mixed ether:acid crystals: a computational study

Acid solvation states are investigated in the recently discovered mixed ether:acid crystalline solids. The solids are simulated using on-the-fly molecular dynamics as implemented in the density functional code QUICKSTEP employing Gaussian basis sets. The solids are shown to display a remarkably broad range of acid solvation states, depending on the ether:acid ratio, including proton sharing in the 1:1 case, proton transfer to the ether in 1:2, and perturbed molecular acid in 1:6. The observed variation of the infrared spectra with the composition is accounted for qualitatively with the help of the calculations.     

Chelate effect and thermodynamics of metal complex formation in solution: a quantum chemical study

The accuracy of quantum chemical predictions of structures and thermodynamic data for metal complexes depends both on the quantum chemical methods and the chemical models used. A thermodynamic analogue of the Eigen-Wilkins mechanism for ligand substitution reactions (Model A) turns out to be sufficiently simple to catch the essential chemistry of complex formation reactions and allows quantum chemical calculations at the ab initio level of thermodynamic quantities both in gas phase and solution; the latter by using the conductor-like polarizable continuum (CPCM) model. Model A describes the complex formation as a two-step reaction: 1. [M(H2O)x](aq) + L(aq) <==>[M(H2O)x], L(aq); 2. [M(H2O)x], L(aq) <==>[M(H2O)(x-1)L],(H2O)(aq). The first step, the formation of an outer-sphere complex is described using the Fuoss equation and the second, the intramolecular exchange between an entering ligand from the second and water in the first coordination shell, using quantum chemical methods. The thermodynamic quantities for this model were compared to those for the reaction: [M(H2O)x](aq) + L(aq) <==>[M(H2O)(x-1)L](aq) + (H2O)(aq) (Model B), as calculated for each reactant and product separately. The models were tested using complex formation between Zn(2+) and ammonia, methylamine, and ethylenediamine, and complex formation and chelate ring closure reactions in binary and ternary UO(2)(2+)-oxalate systems. The results show that the Gibbs energy of reaction for Model A are not strongly dependent on the number of water ligands and the structure of the second coordination sphere; it provides a much more precise estimate of the thermodynamics of complex formation reactions in solution than that obtained from Model B. The agreement between the experimental and calculated data for the formation of Zn(NH(3))(2+)(aq) and Zn(NH(3))(2)(2+)(aq) is better than 8 kJ/mol for the former, as compared to 30 kJ/mol or larger, for the latter. The Gibbs energy of reaction obtained for the UO(2)(2+) oxalate systems using model B differs between 80 and 130 kJ/mol from the experimental results, whereas the agreement with Model A is better. The errors in the quantum chemical estimates of the entropy and enthalpy of reaction are somewhat larger than those for the Gibbs energy, but still in fair agreement with experiments; adding water molecules in the second coordination sphere improves the agreement significantly. Reasons for the different performance of the two models are discussed. The quantum chemical data were used to discuss the microscopic basis of experimental enthalpy and entropy data, to determine the enthalpy and entropy contributions in chelate ring closure reactions and to discuss the origin of the so-called "chelate effect". Contrary to many earlier suggestions, this is not even in the gas phase, a result of changes in translation entropy contributions. There is no simple explanation of the high stability of chelate complexes; it is a result of both enthalpy and entropy contributions that vary from one system to the other.     

Solvation of transmembrane proteins by isotropic membrane mimetics: a molecular dynamics study

Mixtures of organic solvents are often used as membrane mimetics in structure determination of transmembrane proteins by solution NMR; however, the mechanism through which these isotropic solvents mimic the anisotropic environment of cell membranes is not known. Here, we use molecular dynamics simulations to study the solvation thermodynamics of the c-subunit of Escherichia coli F1F0 ATP synthase in membrane mimetic mixtures of methanol, chloroform, and water with varying fractions of components as well as in lipid bilayers. We show that the protein induces a local phase separation of the solvent components into hydrophobic and hydrophilic layers, which provides the anisotropic solvation environment to stabilize the amphiphilic peptide. The extent of this effect varies with solvent composition and is most pronounced in the ternary methanol-chloroform-water mixtures. Analysis of the solvent structure, including the local mole fraction, density profiles, and pair distribution functions, reveals considerable variation among solvent mixtures in the solvation environment surrounding the hydrophobic transmembrane region of the protein. Hydrogen bond analysis indicates that this is primarily driven by the hydrogen-bonding propensity of the essential Asp(61) residue. The impact of the latter on the conformational stability of the solvated protein is discussed. Comparison with the simulations in explicit all-atom models of lipid bilayer indicates a higher flexibility and reduced structural integrity of the membrane mimetic solvated c-subunit. This was particularly true for the deprotonated form of the protein and found to be linked to solvent stabilization of the charged Asp(61).     

Solvation of the Whelk-O1 chiral stationary phase: a molecular dynamics study

Density functional theory calculations and molecular dynamics simulations are employed to explore the solvation of the Whelk-O1 chiral stationary phase. First, a semi-flexible representation of the Whelk-O1 selective molecule is extracted from an extensive series of B3LYP/6-311+ G(2d,p) calculations. The resulting model is used to build a chiral surface, including end-caps, for molecular dynamics study of the interface between solvent and Whelk-O1. Three solvent environments in common use for Whelk-O1 HPLC have been examined: a normal-phase solvent of n-hexane/2-propanol; a reversed-phase solvent of water/methanol; and a supercritical solvent of CO(2) and methanol. In each case, we analyze the interface with an emphasis on solvent composition and solvent hydrogen bonding to the Whelk-O1 selector.     

Solvation free energies and hydration structure of N-methyl-p-nitroaniline

Solvation Gibbs energies of N-methyl-p-nitroaniline (MNA) in water and 1-octanol are calculated using the expanded ensemble molecular dynamics method with a force field taken from the literature. The accuracy of the free energy calculations is verified with the experimental Gibbs free energy data and found to reproduce the experimental 1-octanol∕water partition coefficient to within ±0.1 in log unit. To investigate the hydration structure around N-methyl-p-nitroaniline, an independent NVT molecular dynamics simulation was performed at ambient conditions. The local organization of water molecules around the solute MNA molecule was investigated using the radial distribution function (RDF), the coordination number, and the extent of hydrogen bonding. The spatial distribution functions (SDFs) show that the water molecules are distributed above and below the nitrogen atoms parallel to the plane of aromatic ring for both the methylamino and nitro functional groups. It is found that these groups have a significant effect on the hydration of MNA with water molecules forming two weak hydrogen bonds with both the methylamino and nitro groups. The hydration structures around the functional groups in MNA in water are different from those that have been found for methylamine, nitrobenzene, and benzene in aqueous solutions, and these differences together with weak hydrogen bonds explain the lower solubility of MNA in water. The RDFs together with SDFs provide a tool for the understanding the hydration of MNA (and other molecules) and therefore their solubility.     

石杉碱甲-AChE复合物中石杉碱甲的结构特征: 量子化学研究

在运用量子化学从头计算方法(HF/4-31G)结合点电荷模型方法对AChE-HupA复合物活性位点的410个原子和1929个点电荷进行理论计算的基础上, 比较了石杉碱甲分子在形成复合物前后的结构变化特征。发现复合物中石杉碱甲分子构象并非能量最低构象, 它的能量比HF/4-31G全优化得到的构象的能量高91.8kj/mol。和单分子状态相比, 形成复合物后季铵基和内酰胺基的N-H, C=O键的键长变长、键强减弱, 其总原子净电荷也发生了明显的变化。且这些基团的空间取向都有不同程度的改变, C(8)-N(21)键的旋转达20ⅲ。这些信息将有益于设计新的AChE抑制剂。     

Polymeric structure of a coproporphyrin I ruthenium(II) complex: a powder diffraction study

Porphyrin complexes of ruthenium are widely used as models for the heme protein system, for modelling naturally occurring iron–porphyrin systems and as catalysts in epoxidation reactions. The structural diversity of ruthenium complexes offers an opportunity to use them in the design of multifunctional supramolecular assemblies. Coproporphyrins and metallocoproporphyrins are used as sensors in bioassay and the potential use of derivatives as multiparametric sensors for oxygen and H⁺ is one of the main factors driving a growing interest in the synthesis of new porphyrin derivatives. In the coproporphyrin I Ru^II complex catena‐poly[[carbonylruthenium(II)]‐μ‐2,7,12,17‐tetrakis[2‐(ethoxycarbonyl)ethyl]‐3,8,13,18‐tetramethylporphyrinato‐κ⁵N ,N ′,N ′′,N ′′′:O ], [Ru(C₄₄H₅₂N₄O₈)(CO)]_n , the Ru^II centre is coordinated by four N atoms in the basal plane, and by axial C (carbonyl ligand) and O (ethoxycarbonylethyl arm from a neighbouring complex) atoms. The complex adopts a distorted octahedral geometry. Self‐assembly of the molecules during crystallization from a methylene chloride–ethanol (1:10 v /v ) solution at room temperature gives one‐dimensional polymeric chains.     

Synthesis, structure, and reactivity of a dinuclear metal complex with linear M-H-M bonding

Reaction of two equiv of 1-adamantylzinc bromide with (dippm)NiBr2 (dippm = bis(di-isopropylphosphino)methane) led to a dinuclear metal complex containing a unique linear bridging hydride ligand. The hydride was characterized by neutron diffraction methods, which confirmed a linear bonding mode. Preliminary reactivity studies of this unusual dimer are reported.     

Crystal structure evolution of complex metal aluminum hydrides upon hydrogen release

Complex aluminum hydrides have been widely studied as potential hydrogen storage materials but also,for some time now, for electrochemical applications. This review summarizes the crystal structures of alkali and alkaline earth aluminum hydrides and correlates structure properties with physical and chemical properties of the hydride compounds. The crystal structures of the alkali metal aluminum hydrides change significantly during the stepwise dehydrogenation. The general pathway follows a transformation of structures built of isolated [AlH4]~- tetrahedra to structures built of isolated [Al H6]~(3-) octahedra.The crystal structure relations in the group of alkaline earth metal aluminum hydrides are much more complicated than those of the alkali metal aluminum hydrides. The structures of the alkaline earth metal aluminum hydrides consist of isolated tetrahedra but the intermediate structures exhibit chains of cornershared octahedra. The coordination numbers within the alkali metal group increase with cation sizes which goes along with an increase of the decomposition temperatures of the primary hydrides. Alkaline earth metal hydrides have higher coordination numbers but decompose at slightly lower temperatures than their alkali metal counterparts. The decomposition pathways of alkaline metal aluminum hydrides have not been studied in all cases and require future research.     

Solvation structure of hydroxyl radical by Car-Parrinello molecular dynamics

Car-Parrinello molecular dynamics simulations of a hydroxyl radical in liquid water have been performed. Structural and dynamical properties of the solvated structure have been studied in details. The partial atom-atom radial distribution functions for the hydrated hydroxyl do not show drastic differences with the radial distribution functions for liquid water. The OH is found to be a more active hydrogen bond donor and acceptor than the water molecule, but the accepted hydrogen bonds are much weaker than for the hydroxide OH- ion. The first solvation shell of the OH is less structured than the water's one and contains a considerable fraction of water molecules that are not hydrogen bonded to the hydroxyl. Part of them are found to come closer to the solvated radical than the hydrogen bonded molecules do. The lifetime of the hydrogen bonds accepted by the hydroxyl is found to be shorter than the hydrogen bond lifetime in water. A hydrogen transfer between a water molecule and the OH radical has been observed, though it is a much rarer event than a proton transfer between water and an OH- ion. The velocity autocorrelation power spectrum of the hydroxyl hydrogen shows the properties both of the OH radical in clusters and of the OH- ion in liquid.     

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.     

Intramolecular charge-transfer state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution: RISM-SCF study

Intramolecular charge-transfer (ICT) state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution is studied by the reference interaction site model self-consistent field (RISM-SCF) method. Geometry optimizations are performed for each electronic state in solution with the complete-active-space SCF wave functions. Dynamic electron correlation effects are taken into account by using the multiconfigurational quasidegenerate perturbation theory. Two-dimensional free energy surfaces are constructed as the function of the twisting and wagging angles of the dimethylamino group for the ground and locally excited (LE) states. The calculated absorption and fluorescence energies are in good agreement with experiments. The validity of the twisted ICT (TICT) model is confirmed in explaining the dual fluorescence, and the possibility of the planar ICT model is ruled out. To examine the mechanism of the TICT state formation, a "crossing" seam between the LE and charge-transfer (CT) state surfaces is determined. The inversion of two electronic states occurs at a relatively small twisting angle. The effect of solvent reorganization is also examined. It is concluded that the intramolecular twisting coordinate is more important than the solvent fluctuation for the TICT state formation, because the energy difference between the two states is minimally dependent on the solvent configuration.     

Solvation of small hydrophobic molecules in formamide: A calorimetric study

The solubility and enthalpy of solution of benzene, cyclohexane, hexane, and heptane in formanide have been determined from titration microcalorimetric experiments at 25°C. The solution enthalpies are significantly more endothermic than in water but still the solubility is much higher. The entropy changes in formamide are small and positive and do not vary significantly with size. The enthalpies of solution of some 1-alkanols, 1-chloro- and 1,5-dichloropentane and pentane-1,5-diol were measured as functions of concentration. The solution enthalpies for 1-alkanols from methanol to decanol increase linearly with chain length. The enthalpic interaction coefficients h_xx are small and negative in formamide while they are large and positive in water. The partial molar heat capacities C _p,2 ^o for 1-propanol, 1-pentanol, benzene and cyclohexane in formamide were determined at 25°C from drop heat capacity measurements. Values of C _p,2 ^o are only slightly larger than the molar heat capacities of the liquid solutes.     

Solvation of ethylmaltol and of its iron(III) complex

Summary Solubilities of tris(ethylmaltolato)iron(III) (ethylmaltol = 3-hydroxy-2-ethyl-4-pyrone) were measured in MeOH-H₂O, t-BuOH-H₂O and diol-H₂O mixtures, and in several primary alcohols. Solvation of the ethylmaltol ligand and of two 4-pyridinone analogues has been investigated through solubility measurements in MeOH- H₂O and in t-BuOH-H₂O mixtures, and in a series of primary alcohols. The solvation characteristics of these compounds are compared with those of the parent maltol, its iron(III) complex and a number of other nonelectrolytes.