Solvation effects on equilibria: Triazoles and N-methyl piperidinol

We have performed calculations of the solvation effects on a number of equilibrium constants in water using a recently proposed hybrid quantum classical scheme in which the liquid environment is modelled using classical solvent molecules and the solute electronic structure is computed using modern quantum chemical methods. The liquid phase space is sampled from a fully classical simulation. We find that solvation effects on both triazole tautomeric equilibrium constants and piperidinol conformational equilibrium constants can be interpreted in terms of subtle differences in the local environment which can be seen in probability densities and radial distribution functions. Lower level calculations were performed for comparison and we conclude that the solvation thermodynamics can be predicted from a good classical model of solvent and solute molecules, but the implicit models that we tried are less successful.

Relevance of heterometallic binding energy for metal underpotential deposition

We present first-principles calculations for a number of metals adsorbed on several different metallic substrates. Some of these systems are very relevant in electrochemistry, especially in the field of underpotential deposition phenomena. The present studies reveal the existence of a relationship between the excess binding energy and the surface energy difference between substrate and adsorbate. Comparisons with experimental underpotential shifts show that excess binding energies are systematically underestimated. By analyzing experimental information on different systems, we conclude that this discrepancy between our vacuum calculations and experiments carried out in an electrolytic solution is likely to be due to anion adsorption and/or solvent effects.

Quantum electronic transport in a configuration interaction basis

An overview of a many-body approach to calculation of electronic transport in molecular systems is given. The physics required to describe electronic transport through a molecule at the many-body level, without relying on commonly made assumptions such as the Landauer formalism or linear response theory, is discussed. Physically, our method relies on the incorporation of scattering boundary conditions into a many-body wavefunction and application of the maximum entropy principle to the transport region. Mathematically, this simple physical model translates into a constrained nonlinear optimization problem. A strategy for solving the constrained optimization problem is given. (C) 2004 Wiley Periodicals, Inc.

Monte Carlo simulation of cluster growth in surface defects induced by the tip of a scanning tunnelling microscope

First steps are taken to model the electrochemical deposition of metals in nanometer-sized cavities. In the present work, the electrochemical deposition of Cu atoms in nanometer-sized holes dug on Au(111) is investigated through Monte Carlo simulations using the embedded atom method to represent particle interactions. By sweeping the chemical potential of Cu, a cluster is allowed to grow within the hole rising four atomic layers above the surface. Its lateral extension remains confined to the area defined by the borders of the original defect. (C) 2004 Elsevier B.V. All rights reserved.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

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.

Monte Carlo simulation of properties of monolayers and metal islands adsorbed on metallic (111) surfaces

To obtain the surface stress changes due to the adsorption of metal monolayers onto metallic surfaces, a new model derived from thermodynamic considerations is presented. Such a model is based on continuum Monte Carlo simulations with embedded atom method potentials in the canonical ensemble, and it is extended to consider the behavior on different islands adsorbed onto (111) substrate surfaces. Homoepitaxial and heteroepitaxial systems are studied. Pseudomorphic growth is not observed for small metal islands with considerable positive misfit with the substrate. Instead, the islands become compressed upon increase of their size. A simple model is proposed to interpolate between the misfits of atoms in small islands and the pseudomorphic behavior of the monolayer.

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.

A combined experimental and theoretical study of the generation of palladium clusters on Au(111) with a scanning tunnelling microscope

Palladium clusters have been deposited on the surface of a Au(111) electrode with the tip of a scanning tunnelling microscope. The distance over which the tip was moved towards the surface has a decisive influence on the properties of the clusters: the larger this distance, the larger the generated clusters, and the more stable they are. These findings are supported by computer simulations, which further suggest that the larger clusters contain a sizable amount of gold, which enhances their stability. Dissolution of the clusters occurs from the edges rather than layer by layer.

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

Quantitative analysis of MRP-8 in gingival crevicular fluid in periodontal health and disease using microbore HPLC.

BACKGROUND:

The protein components of GCF can be separated by reverse-phase microbore HPLC on a C18 column with detection on the basis of 214 nm absorbance. A single major symmetrical protein peak eluting with a retention time of 26 min (50% acetonitrile) was evident in gingival crevicular fluid (GCF) from periodontitis patients but not in healthy GCF. This protein was identified as human MRP-8 by N-terminal amino acid sequencing and liquid chromatography quadropole mass spectrometry.

AIMS:

To quantify the amount of MRP-8 detectable in GCF from individual healthy, gingivitis and periodontitis affected sites and to study the relationship, if any, between the levels of this responsive protein and periodontal health and disease.

METHODS:

GCF was sampled (30 s) from healthy, gingivitis, and periodontitis sites in peridontitis subjects (n=15) and from controls (n=5) with clinically healthy gingiva and no periodontitis. Purified MRP-8 was sequenced by Edmann degradation and the phenylthiohydantoin (PTH) amino acid yield determined (by comparison of peak area with external PTH amino acid standards). This value was subsequently used to calculate the relative amount of protein in the peak eluting with a retention time of 26.0 min (MRP-8) in individual GCF chromatograms.

RESULTS:

Higher levels of MRP-8 were detected in inflammatory sites: periodontitis 457.0 (281.0) ng; gingivitis 413.5 (394.5) ng compared with periodontally healthy sites in diseased subjects 14.6 (14.3) ng and in controls 18.6 (18.5) ng, p=0.003. There was at least 20-fold more MRP-8 in the inflammatory compared with the healthy sites studied.

CONCLUSIONS:

The preliminary data indicate that MRP-8 is present in GCF, with significantly greater amounts present at diseased than healthy sites. A systematic study of the relationship of this protein to periodontal disease could prove useful in further clarifying whether MRP-8 could be a reliable GCF biomarker of gingivitis and periodontitis.