Surface plasmon polariton propagation across a gentle silver step

Surface plasmon polaritons (SPPs) are excited with light of wavelength lambda (1) = 632.8 nm on or near a gentle Ag/Ag step structure using focused beam, prism coupling and detected using a bare, sharpened fibre tip. The tip-sample separation is controlled by means of an evanescent optical field at wavelength lambda (2) = 543.5 nm in a photon scanning tunnelling microscope (PSTM). The SPP propagation properties are first characterised on both the thin and thick sections of the Ag film structure either side of the step, both macroscopically, using attenuated total reflection, and microscopically from the PSTM images; the two techniques yield very good agreement. It is found that the SPP propagation length is similar to 10-11 mum across the step in each direction (thick to thin and vice versa) as observed in the PSTM images. Thus, with reference to the propagation lengths of 14.2 and 11.7 mum for the thick and thin planar parts of the Ag film respectively, it is concluded that the SPPs negotiate the step reasonably successfully. Importantly, also, it is shown that images may be produced, displaying SPPs with either an artificially enhanced (similar to 15-20 mum) or truncated (5-8 mum) propagation length across the step. Consideration of such images leads us to suggest the possibility that the photon tunnelling occurs in a local water environment. (C) 2001 Elsevier Science B.V. All rights reserved.

Measurement and interpretation of the mid-infrared properties of single crystal and polycrystalline gold

The contribution of electron-phonon scattering and grain boundary scattering to the mid-IR (lambda = 3.392 mum) properties of An has been assessed by examining both bulk, single crystal samples-Au(1 1 1) and Au(1 1 0)-and thin film, polycrystalline An samples at 300 K and 100 K by means of surface plasmon polariton excitation. The investigation constitutes a stringent test for the in-vacuo Otto-configuration prism coupler used to perform the measurements, illustrating its strengths and limitations. Analysis of the optical response is guided by a physically based interpretation of the Drude model. Relative to the reference case of single crystal Au at 100 K (epsilon = - 568 + i17.5), raising the temperature to 300 K causes increased electron-phonon scattering that accounts for a reduction of similar to40 nm in the electron mean free path. Comparison of a polycrystalline sample to the reference case determines a mean free path due to grain boundary scattering of similar to 17 nm, corresponding to about half the mean grain size as determined from atomic force microscopy and indicating a high reflectance coefficient for the An grain boundaries. An analysis combining consideration of grain boundary scattering and the inclusion of a small percentage of voids in the polycrystalline film by means of an effective medium model indicates a value for the grain boundary reflection coefficient in the range 0.55-0.71. (C) 2005 Elsevier B.V. All rights reserved.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

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.