Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Reactions of One-Electron-Oxidized Methionine with Oxygen: An ab initio study

A one-electron oxidation of a methionine residue is thought to be a key step in the neurotoxicity of the beta amyloid peptide of Alzheimer's disease. The chemistry of the radical cation of N-formylmethioninamide (11+) and two model systems, dimethyl sulfide (1+) and ethyl methyl sulfide (6+), in the presence of oxygen have been studied by B3LYP/6-31G(d) and CBS-RAD calculations. The stable form of 11+ has a three-electron bond between the sulfur radical cation and the carbonyl oxygen atom of the i - 1 residue. The radical cation may lose a proton from the methyl or methylene groups flanking the oxidized sulfur. Both 11+ and the resultant C-centered radicals may add oxygen to form peroxy radicals. The calculations indicate that unlike C-centered radicals the sulfur radical cation does not form a covalent bond to oxygen but rather forms a loose ion-induced dipole complex with an S-O separation of about 2.7 Å, and is bound by about 13 kJ mol-1 (on the basis of 1+ + O2). Direct intramolecular abstraction of an H atom from the C site is unlikely. It is endothermic by more than 20 kJ mol-1 and involves a high barrier (G = 79 kJ mol-1). The -to-S C-centered radicals will add oxygen to form peroxy radicals. The OH BDEs of the parent hydroperoxides are in the range of 352-355 kJ mol-1, similar to SH BDEs (360 kJ mol-1) and C-H BDEs (345-350 kJ mol-1). Thus, the peroxy radicals are oxidizing species comparable in strength to thiyl radicals and peptide backbone C-centered radicals. Each peroxy radical can abstract a hydrogen atom from the backbone C site of the Met residue to yield the corresponding C-centered radical/hydroperoxide in a weakly exothermic process with modest barriers in the range of 64-92 kJ mol-1.

Raman and ab initio studies of simple and binary 1-Alkyl-3-methylimidazolium ionic liquids

Raman spectra of the ionic liquids, 1-butyl-3-methylimidazolium hexafluorophosphate ([C(4)mim][PF(6)]), 1-hexyl-3-methylimidazolium chloride ([C(6)mim]Cl), and 1-hexyl-3-methylimidazolium hexafluorophosphate ([C(6)mim][PF(6)]), and binary mixtures thereof, have been assigned using ab initio MP2 calculations. The previously reported anti and gauche forms of the [C(4)mim](+) cation have been observed, and this study reveals this to be a general feature of the long-chain I-alkyl derivatives. Analysis of mixtures Of [C(6)mim]Cl and [C(6)mim][PF(6)] has provided information on the nature of the hydrogen bonding between the imidazolium headgroup and the anions, and the invariance of the essentially 50:50 mixture of the predominant conformers informs on the nature of glass formation in these systems.

Clusters, liquids and crystals of dialkylimidazolium salts. A combined perspective from ab initio computer simulations

We summarize results obtained by a combination of ab initio and classical computer simulations of dialkylimidazolium ionic liquids in different states of aggregation, from crystals to liquids and clusters. Unusual features arising from the competition between electrostatic, dispersion, and hydrogen-bonding interactions are identified at the origin of observed structural patterns. We also discuss the way Brønsted acids interact with ionic liquids leading to the formation of hydrogen-bonded anions.

Softened C-H modes of adsorbed methyl and their implications for dehydrogenation: An ab initio study

To investigate the softening of CH vibrational frequencies and their implications for dehydrogenation of adsorbed hydrocarbons, an issue of scientific and technological importance, density functional theory calculations have been performed on the chemisorption and dehydrogenation of CH3 on Cu(111) and Pt(111) surfaces. By comparing these results with those of Ni(111) we find that the CH bonds of the adsorbate, when close enough, interact with metal atoms of the surface. It is this interaction and its associated lengthening and weakening of CH bonds that is the physical origin of mode softening. We rule out the possibility of a relationship between the mere presence of mode softening and dehydrogenation. We do show, however, that there is a clear relationship between the extent to which a surface can induce mode softening and the activation energy to dehydrogenation. In addition, periodic trends concerning the extent of mode softening are reproduced. (C) 2001 American Institute of Physics.

Ultrafast Laser-Driven Excitation Dynamics in Ne: An Ab Initio Time-Dependent R-Matrix Approach

An attosecond pump-probe scheme that combines the use of a free-electron laser pulse with an ultrashort pulse is applied in order to explore the ultrafast excitation dynamics in Ne. We describe the multielectron dynamics using a new nonperturbative time-dependent R-matrix theory. This theory enables the interaction of ultrashort light fields with multielectron atoms and atomic ions to be determined from first principles. By probing the emission of an inner 2s electron from Ne we are also able to study the bound state population dynamics during the free-electron laser pulse.

Crystal structure and ab initio calculations of a cyano-carbamimidic acid ethyl ester

The structure of one tautomer (amine form) of cyano-carbamimidic acid ethyl ester or (amino-ethoxy-methylidene)aminoformonitrile (CAS: 13947-84-7) was determined by single crystal X-ray diffraction. Ab initio quantum chemical calculations at the B3LYP, MP2 and G3 levels were performed to investigate the stability and the formation of the different tautomers and conformers. The calculations indicate that the amine form is the more stable tautomer, showing a high degree of election conjugation. The most stable amine conformer located by the calculations corresponds to the crystallized structure. On the contrary, in the less stable imine form, the conjugation is separated by a N2-C2 single bond. (C) 2007 Elsevier B.V. All rights reserved.

Prediction of the formation and stabilities of energetic salts and ionic liquids based on ab initio electronic structure calculations

A computational approach to predict the thermodynamics for forming a variety of imidazolium-based salts and ionic liquids from typical starting materials is described. The gas-phase proton and methyl cation acidities of several protonating and methylating agents, as well as the proton and methyl cation affinities of many important methyl-, nitro-, and cyano- substituted imidazoles, have been calculated reliably by using the computationally feasible DFT (B3LYP) and MP2 (extrapolated to the complete basis set limit) methods. These accurately calculated proton and methyl cation affinities of neutrals and anions are used in conjunction with an empirical approach based on molecular volumes to estimate the lattice enthalpies and entropies of ionic liquids, organic solids, and organic liquids. These quantities were used to construct a thermodynamic cycle for salt formation to reliably predict the ability to synthesize a variety of salts including ones with potentially high energetic densities. An adjustment of the gas phase thermodynamic cycle to account for solid- and liquid-phase chemistries provides the best overall assessment of salt formation and stability. This has been applied to imidazoles (the cation to be formed) with alkyl, nitro, and cyano substituents. The proton and methyl cation donors studied were as follows: HCl, HBr, HI, (HO)(2)SO2, HSO3CF3 (TfOH), and HSO3(C6H4)CH3 (TsOH); CH3Cl, CH3Br, CH3I, (CH3O)(2)SO2, CH3SO3CF3 (TfOCH3) and CH3SO3(C6H4)CH3 (TsOCH3). As substitution of the cation with electron-withdrawing groups increases, the triflate reagents appear to be the best overall choice as protonating and methylating agents. Even stronger alkylating agents should be considered to enhance the chances of synthetic success. When using the enthalpies of reaction for the gas-phase reactants (eq 6) to form a salt, a cutoff value of - 13 kcal mol(-1) or lower (more negative) should be used as the minimum value for predicting whether a salt can be synthesized.