Hydrogen-Bond and Solvent Dynamics in Transition Metal Complexes: A Combined Simulation and NMR-Investigation

Spectroscopy and dynamics of double proton transfer in formic acid dimer

We present the isolated gas phase infrared spectra of formic acid dimer, (HCOOH)2, and its deuterated counterpart formic-d acid, (DCOOH)2, at room temperature. The formic acid dimer spectrum was obtained by spectral subtraction of a spectrum of formic acid vapor recorded at low pressure from that recorded at a higher pressure. The spectra of formic acid vapor contain features from both formic acid monomer and formic acid dimer, but at low and high pressures of formic acid, the equilibrium is pushed towards the monomer and dimer, respectively. A similar approach was used for the formic-d acid dimer. Building on the previous development of the Molecular Mechanics with Proton Transfer (MMPT) force field for simulating proton transfer reactions, molecular dynamics (MD) simulations were carried out to interpret the experimental spectra in the OH-stretching region. Within the framework of MMPT, a combination of symmetric single and double minimum potential energy surfaces (PESs) provides a good description of the double proton transfer PES. In a next step, potential morphing together with electronic structure calculations at the B3LYP and MP2 level of theory was used to align the computed and experimentally observed spectral features in the OH-stretching region. From this analysis, a barrier for double proton transfer between 5 and 7 kcal mol(-1) was derived, which compares with a CCSD(T)/aug-cc-pVTZ calculated barrier of 7.9 kcal mol(-1). Such a combination of experimental and computational techniques for estimating barriers for proton transfer in gas phase systems is generic and holds promise for further improved PESs and energetics of these important systems. Additional MD simulations at the semi-empirical DFTB level of theory agree quite well for the center band position but underestimate the width of the OH-stretching band.

Reproducing kernel potential energy surfaces in biomolecular simulations: Nitric oxide binding to myoglobin

Multidimensional potential energy surfaces based on reproducing kernel-interpolation are employed to explore the energetics and dynamics of free and bound nitric oxide in myoglobin (Mb). Combining a force field description for the majority of degrees of freedom and the higher-accuracy representation for the NO ligand and the Fe out-of-plane motion allows for a simulation approach akin to a mixed quantum mechanics/molecular mechanics treatment. However, the kernel-representation can be evaluated at conventional force-field speed. With the explicit inclusion of the Fe-out-of-plane (Fe-oop) coordinate, the dynamics and structural equilibrium after photodissociation of the ligand are correctly described compared to experiment. Experimentally, the Fe-oop coordinate plays an important role for the ligand dynamics. This is also found here where the isomerization dynamics between the Fe-ON and Fe-NO state is significantly affected whether or not this co-ordinate is explicitly included. Although the Fe-ON conformation is metastable when considering only the bound (2)A state, it may disappear once the (4)A state is included. This explains the absence of the Fe-ON state in previous experimental investigations of MbNO.

The performance of methallyl nickel complexes and boron adducts in the catalytic activation of ethylene: a conceptual DFT perspective

In this work, global and local descriptors of chemical reactivity and selectivity are used to explain the differences in reactivities toward ethylene of methallyl nickel complexes and their B(C6F5)3 and BF3 adducts. DFT calculations were used to explain why nickel complexes alone are inactive in ethylene polymerization while their boron adducts can activate it. It is shown that chemical potential, hardness, electrophilicity and molecular electrostatic potential surfaces describe fairly well the reactivity and selectivity of these organometallic systems toward ethylene. Experimental data indicates that addition of a borane molecule to nickel complexes changes dramatically their reactivity-behavior that is confirmed computationally. Our results show that bare complexes are unable to activate ethylene-a Lewis base-because they also behave as Lewis bases. The addition of the co-catalyst-a Lewis acid-turns the adducts into Lewis acids, making them active towards ethylene.

Full kinetic analysis of a rhodium-catalyzed hydroformylation: beyond the rate-limiting step picture

A complete dynamic kinetic analysis beyond the steady state approximation of the rhodium-catalyzed hydroformylation with the 6-DPPon ligand is presented. The results show that not one single step but several transition states and intermediates control the selectivity and activity of the catalysis.

Molecular insights into the damping mechanism of poly(vinyl acetate)/hindered phenol hybrids by a combination of experiment and molecular dynamics simulation

By combining experiment and MD simulation, the relationship between hydrogen bond evolution and damping property variation of PVAc was revealed.

Toward force fields for atomistic simulations of iridium-containing complexes

The structural and energetic characterization of metal complexes is important in catalysis and photochemical applications. Unraveling their modes-of-action can be greatly assisted by computation, which typically is restricted to computationally demanding methods including electronic structure calculations with density functional theory. Here, we present an empirical force field based on valence bond theory applicable to a range of octahedral Ir(III) complexes with different coordinating ligands, including iridium complexes with a chiral P,N ligand. Using an approach applicable to metal-containing complexes in general, it is shown that with one common parametrization 85% of the 116 diastereomers--all within 21 kcal/mol of the lowest energy conformation of each series--can be correctly ranked. For neutral complexes, all diastereomers are ranked correctly. This helps to identify the most relevant diastereomers which, if necessary, can be further investigated by more demanding computational methods. Furthermore, if one specific complex is considered, the root mean square deviation between reference data from electronic structure calculations and the force field is ≈1 kcal/mol. This, together with the possibility to carry out explicit simulations in solution paves the way for an atomistic understanding of iridium-containing complexes in catalysis.