Search

We present Car-Parrinello molecular dynamics simulations of the liquid imidazolium chloride/AlCl3 by inserting one pair of [C2C1im]Cl into liquid AlCl3, forming an acidic mixture. Two different starting conditions lead to two trajectories from which we harvest structural data. For both simulations, we find large anions within the equilibrium phase: In both trajectories at longer simulation time, the anion size converges to four monomer units, i.e., to Al4Cl13-. The cluster size fluctuations indicate that Grotthus diffusion must play a role. We discuss one possible mechanism of such a reaction changing the anionic species. This process involves many steps of chlorine rattling, bond breaking, and bond forming. With the aid of the electron localization function, a probable rationale for the formation of larger anions is determined: Large anionic species are formed simply to account for the "lack of electrons" present in the acidic melt.

Theoretical bioinorganic and biomimetic chemistry involves the careful description of the electronic structure: for example, 'valence bond reading' of broken-symmetry density functional theory computations gives insight into the structure and bonding of metal-radical systems with complex electronic structures. Exploring the reactivities of such systems leads to the design of novel compounds with better reactivities. Combined quantum-mechanics/molecular-mechanics (QM/MM), where the QM part is a sophisticated ab initio method, aids in understanding nature's most complicated reaction mechanisms in atomic detail. First principles molecular dynamics simulations (Car-Parrinello simulations) open up exciting new avenues for studying transition metal centers and enable several questions to be addressed that cannot be resolved with either standard quantum chemical or traditional force-field methods.

A vibrational analysis of 2-fold hydrogen bonds between an isophthalic amide donor and different acceptors is presented. These systems can be considered as mimetics for the hydrogen-binding situation of numerous supramolecular compounds such as rotaxanes, catenanes, knotanes, and anion receptors. We calculated pronounced red-shifts up to 65 cm(-1) for the stretching modes of the acceptor carbonyl as well as for the donor NH2 groups, whereas we observe a blue shift for the NH2 bending modes and an additional weak hydrogen bond between the acceptor and the middle C-H group of the donor. The red and blue shifts observed for different modes in various complexes have been correlated with the binding energy of the complexes, independently. In comparison with comparable single hydrogen bonds, we find for the 2-fold hydrogen bonds smaller red shifts for the N-H stretch modes of the donor but larger red shifts for the C=O stretch mode of the acceptor. Furthermore, our results indicate that the pronounced blue shift of the C-H stretch mode is basically caused by the fact that the acceptor is fixed directly above this group due to the 2-fold hydrogen bond.

Protonated betaine bis(trifluoromethylsulfonyl)imide is an ionic liquid with the ability to dissolve large quantities of metal oxides. This metal-solubilizing power is selective. Soluble are oxides of the trivalent rare earths, uranium(VI) oxide, zinc(II) oxide, cadmium(II) oxide, mercury(II) oxide, nickel(II) oxide, copper(II) oxide, palladium(II) oxide, lead(II) oxide, manganese(II) oxide, and silver(I) oxide. Insoluble or very poorly soluble are iron(III), manganese(IV), and cobalt oxides, as well as aluminum oxide and silicon dioxide. The metals can be stripped from the ionic liquid by treatment of the ionic liquid with an acidic aqueous solution. After transfer of the metal ions to the aqueous phase, the ionic liquid can be recycled for reuse. Betainium bis(trifluoromethylsulfonyl)imide forms one phase with water at high temperatures, whereas phase separation occurs below 55.5 degrees C (temperature switch behavior). The mixtures of the ionic liquid with water also show a pH-dependent phase behavior: two phases occur at low pH, whereas one phase is present under neutral or alkaline conditions. The structures, the energetics, and the charge distribution of the betaine cation and the bis(trifluoromethylsulfonyl)imide anion, as well as the cation-anion pairs, were studied by density functional theory calculations.

In this paper we analyse the electronic properties of gas-phase 1-butyl-3-methylimidazolium Cl ion pairs, [C(4)C(1)im]Cl, in order to deepen our understanding of ionic liquids in general. Examination of charge densities, natural bond orbitals (NBO), and delocalised molecular orbitals computed at the B3LYP and MP2/6-31(++)G(d,p) levels have enabled us to explain a number of experimental phenomena: the relative acidity of different sites on the imidazolium ring, variations in hydrogen-bond donor and acceptor abilities, the apparent contradiction of the hydrogen-bond-donor parameters for different types of solute, the low probability of finding a Cl(-) anion at the rear of the imidazolium ring and the expansion of the imidazolium ring in the presence of a strong hydrogen-bond acceptor. The unreactive but coordinating environment and large electrochemical window have also been accounted for, as has the strong electron-donating character of the carbon atoms to the rear of the ring in associated imidazolylidenes. The electronic structure of the [C(4)C(1)im](+) cation is best described by a C(4)==C(5) double bond at the rear, and a delocalised three-centre 4 e(-) component across the front (N(1)-C(2)-N(3)) of the imidazolium ring; delocalisation between these regions is also significant. Hydrogen-bond formation is driven by Coulombic stabilisation, which compensates for an associated destabilisation of the electronic part of the system. Interactions are dominated by a large positive charge at C(2) and the build up of pi-electron density above and below the ring, particularly that associated with the double bond between C(4) and C(5). The NBO partial charges have been computed and compared with those used in a number of classical simulations.

The properties of isolated AlCl3 clusters and the bulk system are investigated by means of static and dynamic electronic structure methods. We find important structural motifs with the edge connectivity dominant in a dimer and the corner connectivity dominant in a trimer. Furthermore, the trimer cluster exhibits an interesting ring structure with large cooperative effects relative to the dimer. Comparing the found structural motifs in isolated molecule calculations with the structure of the liquid allows us to determine the dominance of edge connectivity in the liquid. The size of the clusters present in the liquid indicates indeed that the dimer is the most abundant species, but there are also trimers, tetramers, and pentamers present. From the local dipole analysis both for the isolated clusters as well as for the liquid, further proof for the edge connectivity is given. However, all results point to the fact that there is also some small percentage of corner connectivity present that might be attributed to the most stable corner-connected cluster, namely the trimer. Importantly, we find that energetic considerations of isolated (static) clusters only do not represent the findings in liquid phase. Instead, a quantum cluster equilibrium approach or simulations are needed.

Cooperativity in ionic liquids is investigated by means of static quantum chemical calculations. Larger clusters of the dimethylimidazolium cation paired with a chloride anion are calculated within density functional theory combined with gradient corrected functionals. Tests of the monomer unit show that density functional theory performs reasonably well. Linear chain and ring aggregates have been considered and geometries are found to be comparable with liquid phase structures. Cooperative effects occur when the total energy of the oligomer differs from a simple sum of monomer energies. Cooperative effects have been found in the structural motifs examined. A systematic study of linear chains of increasing length (up to nine monomer units) has shown that cooperativity plays a more important role than expected and is stronger than in water. The Cl...H distance of the chloride to the most acidic proton increases with an increasing number of monomer units. The average bond distance approaches 218.9 pm asymptotically. The dipole moment grows almost linearly and the dipole moment per monomer unit reaches the asymptotic value of 16.3 D. The charge on the chloride atoms decreases with an increasing chain length. In order to detect local hydrogen bonding in the clusters a new parametrization of the shared-electron number method is introduced. We find decreasing hydrogen bond energies with an increasing cluster size for both the first hydrogen bond to the most acidic proton and the average hydrogen bond.

In this Article we extend the idea of detecting a hydrogen bond solely on one single quantum chemically determined descriptor. We present an improvement of the method introduced by Reiher et al. (Theor. Chim. Acta 2001, 106, 379), who mapped the strength of the hydrogen bond onto an easily accessible quantity, namely, the two-center shared-electron number sigma(HA). First, we show that the linear dependence between the interaction energy from the supermolecular approach and sigma(HA) is valid for a test set of about 120 hydrogen-bonded complexes. Furthermore, we demonstrate that a classification according to acceptor atoms of the hydrogen-bonded complexes can give more accurate results. We thus recommend to detect hydrogen bonds with a specific acceptor atom according to our subset linear regression analysis. Case studies on alcohols and isolated base pairs and trimers from RNA and DNA show the utility of the detection criterion. The shared-electron number method yields that the strength of the N1...N3 hydrogen bond is in the range of 30 kJ/mol. Furthermore the A-U pair is indeed stronger bound than the A-T complex if environmental effects are incorporated in the calculations.

Density functional theory results on key steps of Schrock's catalytic cycle are presented. The quantum chemical modeling of the dinitrogen-reducing reaction steps is conducted without simplifying the bulky HIPT [HIPT = 3,5-(2,4,6-iPr(3)C(6)H(2))(2)C(6)H(3)] substituents at the triamidoamine ligand.

Chiral, enantiopure metallo-supramolecular rhombs self-assemble in solution through coordination of bis-pyridyl-substituted ligands with (en)M(NO3)2 (en = ethylenediamine, M = Pd(II), Pt(II)). Characterization by NMR and CD spectroscopy in solution and by ESI-FT-ICR mass spectrometry in the gas phase suggests that an equilibrium exists in water/methanol of a major 2:2 complex and a minor 3:3 complex of ligands and metal corners. In the gas phase, doubly charged 2:2 complexes fragment into two identical singly charged halves followed by metal-mediated C-H and C-C bond activation reactions within the ethylenediamine ligands. Electrochemical scanning tunneling microscopy (EC-STM) provides in situ imaging of the complexes even with submolecular resolution. Flat-lying rhombs are deposited under potential control from an aqueous electrolyte on a Cu(100) electrode surface precovered by a tetragonal pattern of chloride anions from the supporting electrolyte. Chirality induces the formation of only one domain orientation. Density functional calculations help to interpret the STM images.

In several recent studies Schrock and collaborators demonstrated for the first time how molecular dinitrogen can be catalytically transformed under mild and ambient conditions to ammonia by a molybdenum triamidoamine complex. In this work, we investigate the geometrical and electronic structures involved in this process of dinitrogen activation with quantum chemical methods. Density functional theory (DFT) has been employed to calculate the coordination energies of ammonia and dinitrogen relevant for the dissociation/association step in which ammonia is substituted by dinitrogen. In the DFT calculations the triamidoamine chelate ligand has been modeled by a systematic hierarchy of increasingly complex substituents at the amide nitrogen atoms. The most complex ligand considered is an experimentally known ligand with an HMT = 3,5-(2,4,6-Me3C6H2)2C6H3 substituent. Several assumptions by Schrock and collaborators on key reaction steps are confirmed by our calculations. Additional information is provided on many species not yet observed experimentally. Particular attention is paid to the role of the charge of the complexes. The investigation demonstrates that dinitrogen coordination is enhanced for the negatively charged metal fragment, that is, coordination is more favorable for the anionic metal fragment than for the neutral species. Coordination of N2 is least favorable for the cationic metal fragment. Furthermore, ammonia abstraction from the cationic complex is energetically unfavorable, while NH3 abstraction is less difficult from the neutral and easily feasible from the anionic low-spin complex.

We implemented the quantum cluster equilibrium theory in our postprocessing program PEACEMAKER. This program may be run in conjunction with the very efficient vibrational frequency analysis code SNF and can therefore provide access to all electronic structure programs combined with this program. We applied the quantum cluster equilibrium theory in order to investigate the influence of a wide range of electronic structure models on the description of the liquid state. This investigation revealed much about the relevance of approximations in modern simulations of associated liquids such as water. While it is often claimed that the use of density-functional theory in condensed matter is leading to gravely erroneous results, we found that, contrary to these assertions, the exact exchange functional B3LYP and the gradient-corrected functional BP perform very well in combination with sizable basis sets as compared to second-order Moller-Plesset perturbation theory employing the same basis set. The use of density-functional theory with smaller basis sets does, in fact, lead to better results in the liquid state than the use of second-order Moller-Plesset perturbation theory in combination with these small basis sets. Most importantly, the neglect of cooperative effects disturbs a good description much more evenly if we apply second-order Moller-Plesset perturbation theory in combination with large basis sets than density-functional theory including cooperativity with smaller basis sets or Hartree-Fock using a very small basis set.

A biomimetic complex which mimics the arginine-phosphonate diester interaction of the arginine fork is investigated with respect to structure and energetics of stable configurations. Within this work, we provide knowledge on local minima of the isolated system obtained from first-principles calculations. Non-negligible solvation effects are studied in a microsolvation approach. The interactions which govern the structural patterns of molecular recognition in this tweezer-guest complex can be significantly modulated by the action of hydrogen bond accepting and donating solvent molecules, such as dimethyl sulfoxide or water, which were present in experimental investigations on this system. Different tweezer-guest structures are evaluated with respect to their temperature-dependent thermodynamical properties as products of the first association reaction step of the bisphosphonate tweezer template and the guanidinium moiety.

We have studied reduction reactions for nitrogen fixation at Sellmann-type model complexes with Car-Parrinello simulation techniques. These dinuclear complexes are especially designed to emulate the so-called open-side FeMoco model. The main result of this work shows that in order to obtain the reduced species several side reactions have to be suppressed. These involve partial dissociation of the chelate ligands and hydrogen atom transfer to the metal center. Working at low temperature turns out to be one necessary pre-requisite in carrying out successful events. The successful events cannot be described by simple reaction coordinates. Complicated processes are involved during the initiation of the reaction. Our theoretical study emphasizes two experimental strategies which are likely to inhibit the side reactions. Clamping of the two metal fragments by a chelating phosphane ligand should prevent dissociation of the complex. Furthermore, introduction of tert-butyl substituents could improve the solubility and should thus allow usage of a wider range of (mild) acids, reductants, and reaction conditions.

A general photochemical activation process of inert dinitrogen coordinated to two metal centers is presented on the basis of high-level DFT and ab initio calculations. The central feature of this activation process is the occupation of an antibonding pi* orbital upon electronic excitation from the singlet ground state S0 to the first excited singlet state S1. Populating the antibonding LUMO weakens the triple bond of dinitrogen. After a vertical excitation, the excited complex may structurally relax in the S1 state and approaches its minimum structure in the S1 state. This excited-state minimum structure features the dinitrogen bound in a diazenoid form, which exhibits a double bond and two lone pairs localized at the two nitrogen atoms, ready to be protonated. Reduction and de-excitation then yield the corresponding diazene complex; its generation represents the essential step in a nitrogen fixation and reduction protocol. The consecutive process of excitation, protonation, and reduction may be rearranged in any experimentally appropriate order. The protons needed for the reaction from dinitrogen to diazene can be provided by the ligand sphere of the complexes, which contains sulfur atoms acting as proton acceptors. These protonated thiolate functionalities bring protons close to the dinitrogen moiety. Because protonation does not change the pi*-antibonding character of the LUMO, the universal and well-directed character of the photochemical activation process makes it possible to protonate the dinitrogen complex before it is irradiated. The pi*-antibonding LUMO plays the central role in the activation process, since the diazenoid structure was obtained by excitation from various occupied orbitals as well as by a direct two-electron reduction (without photochemical activation) of the complex; that is, the important bending of N2 towards a diazenoid conformation can be achieved by populating the pi*-antibonding LUMO.

We present an efficient implementation for the calculation of maximally localized Wannier functions (MLWFs) during parallel Car-Parrinello molecular dynamics simulations. The implementation is based on a block Jacobi method. The calculation of MLWFs results in only a moderate (10%-20%) increase in computer time. Consequently it is possible to calculate MLWFs routinely during Car-Parrinello simulations. The Wannier functions are then applied to derive molecular dipole moments of dimethyl sulfoxide (DMSO) in gas phase and aqueous solution. We observe a large increase of the local dipole moment from 3.97 to 7.39 D. This large solvent effect is caused by strong hydrogen bonding at the DMSO oxygen atom and methyl groups. Decomposing the dipole moment into local contributions from the S-O bond and the methyl groups is used to understand the electrostatic response of DMSO in aqueous solution. A scheme is given to derive charges on individual atoms from the MLWFs using the D-RESP methodology. The charges also display large solvent effects and give insight into the transferability of recent force field models for DMSO.

An approach is presented that allows NMR relaxation rates to be determined for a complex mixture, and it is applied to a dimethyl sulfoxide/water solution. This approach is novel for such systems, having only been used for simple systems such as atomic liquids or atomic ions in liquids until now. It involves use of a predetermined, quantum mechanical, multidimensional property surface in a simulation. The results are used in conjunction with the simulated rotational correlation time to calculate the deuteron quadrupole coupling constant (DQCC), in an analogous approach to the one used by experimentalists, and to examine the surprising experimental findings for the composition dependence of the DQCC in the dimethyl sulfoxide/water mixture. Experiments have suggested that the DQCC for a mixture of 5% dimethyl sulfoxide in water is close to the DQCC of ice, whereas its value increases to a value close to the gas value with further dilution.(1) The results are further critically analyzed using combinations of different experimental and theoretical results from the literature.

Using Car-Parrinello molecular dynamics a structural diffusion mechanism for the simplest hydrophobic species in water, an H atom, is proposed. The hydrophobic solvation cavity is a highly dynamical aggregate that actually drives, by its own hydrogen-bond fluctuations, the diffusion of the enclosed solute. This makes possible an anomalously fast diffusion that falls only short of that of "Grotthuss structural diffusion" of H+ in water. Here, the picture of a static, i.e., "iceberglike," clathrate cage is a misleading concept. The uncovered scenario is similar to the "dynamical hole mechanism" found in a very different context, that is, large molecules moving in hot polymeric melts.

DMSO-water mixtures exhibit a marked freezing point depression, reaching close to 60 K at n(DMSO) = 0.33. The phase diagram indicates that stable DMSO-water clusters may be responsible for this phenomenon. Using time-independent quantum chemical methods, we investigate possible candidates for stable supermolecules at mole fractions n(DMSO) = 0.25 and 0.33. The model clusters are built by adding various numbers of water molecules to a single DMSO molecule. Structures and interaction energetics are discussed in the light of experimental and theoretical results from the literature. A comparison with results from molecular dynamics simulations is of particular interest. Our optimized structures are spatially very different from those previously identified through MD simulations. To identify the structural patterns characterizing the clusters, we classify them on the basis of hydrogen-acceptor interactions. These are well separated on an interaction energy scale. For the hydrophobic interactions of the methyl groups with water, attractive interactions of up to 8 kJ/mol are found. In forming clusters corresponding to a range of different mole fractions, up to four water molecules are added to each DMSO molecule. This corresponds to a rough local model of solvation. Examination of the trends in the interactions indicates that the methyl-water interaction becomes more important upon solvation. Finally, we investigate how the clusters interact and attempt to explain which role is played by the various structures and their intercluster interaction modes in the freezing behavior of DMSO-water.