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.

Simulation of the surface structure of butylmethylimidazolium ionic liquids

Molecular dynamics simulations of the liquid/vacuum surfaces of the room temperature ionic liquids [bmim][PF6], [bmim][BF4] and [bmim][Cl] have been carried out at various temperatures. The surfaces are structured with a top monolayer containing oriented cations and anions. The butyl side chains tend to face the vacuum and the methyl side chains the liquid. However, as the butyl chains are not densely packed, both anions and rings are visible from the vacuum phase. The effects of temperature and the anion on the degree of cation orientation is small, but the potential drop from the vacuum to the interior of the liquid is greater for liquids with smaller anions. We compare the simulation results with a range of experimental observations and suggest that neutron reflection from samples with protiated butyl groups would be a sensitive probe of the structure.

On the generation of metal clusters with the electrochemical scanning tunneling microscope

In the present work we consider two aspects of the deposition of metal clusters on an electrode surface. The formation of such clusters with the tip of a scanning tunneling microscope is simulated by atom dynamics. Subsequently the stability of these clusters is investigated by Monte Carlo simulations in a grand-canonical ensemble. In particular, the following systems were considered explicitly: Pd clusters on Au(111), Cu on Au(111), Ag on Au(111), Pb on Au(111) and Cu on Ag(111). The analysis of the results obtained for the different systems leads to the conclusion that optimal systems for nanostructuring are those where the metals participating have similar cohesive energies and negative heats of alloy formation. In this respect, the system Cu-Pd(111) is predicted as a good candidate for the formation of stable clusters. (c) 2005 Elsevier B.V. All rights reserved.

Kinetic Monte Carlo study of electrochemical growth in a heteroepitaxial system

Structural and kinetic aspects of 2-D irreversible metal deposition under potentiostatic conditions are analyzed by means of dynamic Monte Carlo simulations employing embedded atom potentials for a model system. Three limiting models, all considering adatom diffusion, were employed to describe adatom deposition. The first model (A) considers adatom deposition on any free substrate site on the surface at the same rate. The second model (B) considers adatom deposition only on substrate sites which exhibit no neighboring sites occupied by adatoms. The third model (C) allows deposition at higher rates on sites presenting neighboring sites occupied by adatoms. Under the proper conditions, the coverage (theta) versus time (t) relationship for the three cases can be heuristically fitted to the functional form theta = 1 - exp(-betat(alpha)), where alpha and beta are parameters. We suggest that the value of the parameter alpha can be employed to distinguish experimentally between the three cases. While model A trivially delivers a = 1, models B and C are characterized by alpha 1, respectively.

Glory and thresholds effects in H+D-2 reactive angular scattering

We analyse H + D-2 reactive angular scattering using the S- matrix elements obtained by Aoiz et al. and Althorpe et al. Enhancement of small angle scattering in the v' = 3

Novel Porphyrin-Incorporated Hydrogels for Photoactive Intraocular Lens Biomaterials

Novel surface-modified hydrogel materials have been prepared by binding charged porphyrins TMPyP (tetrakis-(4-N-methylpyridyl)porphyrin) and TPPS (tetrakis(4-sulfonatophenyl)porphyrin) to copolymers of HEMA (2-hydroxyethyl methacrylate) with either MAA (methacrylic acid) or DEAEMA (2-(diethylamino)ethylmethacrylate). The charged hydrogels display strong electrostatic interactions with the appropriate cationic or anionic porphyrins to give materials which are intended to be used to generate cytotoxic singlet oxygen (1O2) on photoexcitation and can therefore be used to reduce postoperative infection of the intraocular hydrogel-based replacement lenses that are used in cataract surgery. The UV/vis spectra of TMPyP in MAA:HEMA copolymers showed a small shift in the Soret band and a change from single exponential (161 ���­s) triplet decay lifetime in solution to a decay that could be fitted to a biexponential fit with two approximately equal components with ���´ ) 350 and 1300 ���­s. O2 bubbling reduced the decay to a dominant (90%) component with a much reduced lifetime of 3 ���­s and a minor, longer lived (20 ���­s) component. With D2O solvent the 1O2 lifetime was measured by 1270 nm fluorescence as 35 ���­s in MAA:HEMA, compared to 67 ���­s in solution, although absorbance-matched samples showed similar yield of 1O2 in the polymers and in aqueous solution. In contrast to the minor perturbation in photophysical properties caused by binding TMPyP to MAA:HEMA, TPPS binding to DEAEMA:HEMA copolymers profoundly changed the 1O2 generating ability of the TPPS. In N2-bubbled samples, the polymer-bound TPPS behaved in a similar manner to TMPyP in its copolymer host; however, O2 bubbling had only a very small effect on the triplet lifetime and no 1O2 generation could be detected. The difference in behavior may be linked to differences in binding in the two systems. With TMPyP in MAA:HEMA, confocal fluorescence microscopy showed significant penetration of the porphyrin into the core of the polymer film samples (>150 ���­m). However, for TPPS in DEAEMA:HEMA copolymers, although the porphyrin bound much more readily to the polymer, it remained localized in the first 20 ���­m, even in heavily loaded samples. It is possible that the resulting high concentration of TPPS may have cross-linked the hydrogels to such an extent that it significantly reduced the solubility and/or diffusion rate of oxygen into the doped polymers. This effect is significant since it demonstrates that even simple electrostatic binding of charged porphyrins to hydrogels can have an unexpectedly large effect on the properties of the system as a whole. In this case it makes the apparently promising TPPS/DEAEMA:HEMA system a poor candidate for clinical application as a postoperative antibacterial treatment for intraocular lenses while the apparently equivalent cationic system TMPyP/MAA:HEMA displays all the required properties.

Modulation of surface charge, particle size and morphological properties of chitosan-TPP nanoparticles intended for gene delivery

This work investigates the polyanion initiated gelation process in fabricating chitosan-TPP (tripolyphosphate) nanoparticles in the size range of 100-250 nm intended to be used as carriers for the delivery of gene or protein macromolecules. It demonstrates that ionic gelation of cationic chitosan molecules offers a flexible and easily controllable process for systematically and predictably manipulating particle size and surface charge which are important properties in determining gene transfection efficacy if the nanoparticles are used as non-viral vectors for gene delivery, or as delivery carriers for protein molecules. Variations in chitosan molecular weight, chitosan concentration, chitosan to TPP weight ratio and solution pH value were examined systematically for their effects on nanoparticle size, intensity of surface charge, and tendency of particle aggregation so as to enable speedy fabrication of chitosan nanoparticles with predetermined properties. The chitosan-TPP nanoparticles exhibited a high positive surface charge across a wide pH range, and the isoelectric point (IEP) of the nanoparticles was found to be at pH 9.0. Detailed imaging analysis of the particle morphology revealed that the nanoparticles possess typical shapes of polyhedrons (e.g., pentagon and hexagon), indicating a similar crystallisation mechanism during the particle formation and growth process. This study demonstrates that systematic design and modulation of the surface charge and particle size of chitosan-TPP nanoparticles can be readily achieved with the right control of critical processing parameters, especially the chitosan to TPP weight ratio. (c) 2005 Elsevier B.V. All rights reserved.

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.

Polarization Relaxation in an Ionic Liquid Confined between Electrified Walls

The response of a room temperature molten salt to an external electric field when it is confined to a nanoslit is studied by molecular dynamics simulations. The fluid is confined between two parallel and oppositely charged walls, emulating two electrified solid-liquid interfaces. Attention is focused on structural, electrostatic, and dynamical properties, which are compared with those of the nonpolarized fluid. It is found that the relaxation of the electrostatic potential, after switching the electric field off, occurs in two stages. A first, subpicosecond process accounts for 80% of the decay and is followed by a second subdiffusive process with a time constant of 8 ps. Diffusion is not involved in the relaxation, which is mostly driven by small anion translations. The relaxation of the polarization in the confined system is discussed in terms of the spectrum of charge density fluctuations in the bulk.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Simulation of interfaces between room temperature ionic liquids and other liquids

The structure and properties of the interfaces between the room temperature ionic liquid dimethylimidazolium chloride ([dmim]Cl) and different Lennard-Jones fluids and between ionic liquid and water have been studied by molecular dynamics simulations, and compared to the ionic liquid-vapour interface. Two contrasting types of interface were investigated, thermodynamically stable interfaces between ionic liquid and vapour and between ionic liquid and Lennard-Jones fluids, and diffusing interfaces between miscible phases of different compositions involving water. The density profiles of different species through the interface are presented. The cations and water molecules near the former type of interface are aligned relative to the surface, but no orientational preference was found near or in the broad diffusing interface. The ionic liquid has a negative electrostatic potential relative to vapour or Lennard-Jones fluid, but is more positive than pure water. This contrast is explained in terms of the relative importance of orientation and concentration differences in the two types of interface.