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.

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

Bradykinins and their precursor Cdnas from the skin of the Fire-Bellied Toad (Bombina orientalis)

Bradykinin and (Thr6)-bradykinin have been identified in the defensive skin secretion of the fire-bellied toad, Bombina orientalis. The homologous cDNAs for both peptides were cloned from a skin library using a 3'- and 5'-RACE strategy. Kininogen-1 (BOK-1) contained an open-reading frame of 167 amino acid residues encoding four repeats of bradykinin, and kininogen-2 (BOK-2) contained an open-reading frame of 161 amino acid residues encoding two repeats of (Thr6)-bradykinin. Alignment of both precursor nucleotide and amino acid sequences revealed a high degree of structural similarity. These amphibian skin kininogens/preprobradykinins are not biologically analogous to mammalian kininogens.

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.