SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

Surface Energy and Surface Proton Order of Ice Ih

Ice Ih is comprised of orientationally disordered water molecules giving rise to positional disorder of the hydrogen atoms in the hydrogen bonded network of the lattice. Here we arrive at a first principles determination of the surface energy of ice Ih and suggest that the surface of ice is significantly more proton ordered than the bulk. We predict that the proton order-disorder transition, which occurs in the bulk at similar to 72 K, will not occur at the surface at any temperature below surface melting. An order parameter which defines the surface energy of ice Ih surfaces is also identified.

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.

Characterisation of the Interaction of Lactate Dehydrogenase With Tween-20 Using Isothermal Titration Calorimetry, Interfacial Rheometry and Surface Tension Measurements

In this study the nature of the interaction between Tween-20 and lactate dehydrogenase (LDH) was investigated using isothermal titration calorimetry (ITC). In addition the effects of the protein and surfactant on the interfacial properties were followed with interfacial rheology and surface tension measurements in order to understand the mechanism by which the surfactant prevents protein adsorption to the air– water interface. Comparisons were made with Tween-40 and Tween-80 in order to further investigate the mechanism. ITC measurements indicated a weak, probably hydrophobic, interaction between Tween-20 and LDH. Prevention of LDH adsorption to the air–water interface by the Tween surfactants was correlated with surface energy rather than surfactant CMC. While surface pressure appears to be the main driving force for the displacement of LDH from the air–water interface by Tween-20 a solubilisation mechanism may exist for other protein molecules. More generally the results of this study highlight the value of the use of ITC and interfacial measurements in characterising the surface behaviour of mixed surfactant and protein systems.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.

Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Polymer nanotubes by wetting of ordered porous templates

We have developed a simple technique for the fabrication of polymer nanotubes with a monodisperse size distribution and uniform orientation. When either a polymer melt or solution is placed on a substrate with high surface energy, it will spread to form a thin film, known as a precursor film, similar to the behavior of low molar mass liquids. Similar wetting phenomena occur if porous templates are brought into contact with polymer solutions or melts: A thin surface film will cover the pore walls in the initial stages of wetting. This is because the cohesive driving forces for complete filling are much weaker than the adhesive forces. Wall wetting and complete filling of the pores thus take place on different time scales. The latter is prevented by thermal quenching in the case of melts or by solvent evaporation in the case of solutions, thus preserving a nanotube structure. If the template is of monodisperse size distribution, aligned or ordered, so are the nanotubes, and ordered polymer nanotube arrays can be obtained if the template is removed. Any melt-processible polymer, such as polytetrafluoroethylene (PTFE), blends, or multicomponent solutions can be formed into nanotubes with a wall thickness of a few tens of nanometers. Owing to its versatility, this approach should be a promising route toward functionalized polymer nanotubes.

Imaging domains in BaTiO3 single crystal nanostructures: comparing information from transmission electron microscopy and piezo-force microscopy

This article compares and contrasts information
obtained, using transmission electron microscopy (TEM)
and piezo-force microscopy (PFM), on domain configurations
adopted in single crystal lamellae of BaTiO3, that had
been cut directly from bulk using a focused ion beam
microscope with top and bottom surfaces parallel to
{100}pseudocubic. Both forms of imaging reveal domain
walls parallel to {110}pseudocubic, consistent with sets of 90
domains with dipoles oriented parallel to the two
\001[pseudocubic directions in the plane of the lamellae.
However, the domain width was observed to be dramatically
larger using PFM than it was using TEM. This suggests
significant differences in the surface energy densities
that drive the domain formation in the first place, that could
relate to differences in the boundary conditions in the two
modes of imaging (TEM samples are imaged under high
vacuum, whereas PFM imaging was performed in air).
Attempts were made to map local dipole orientations
directly, using a form of ‘vector’ PFM. However, information
inferred was largely inconsistent with the known
crystallography of the samples, raising concern about the
levels of care needed for accurate interpretation of PFM
images.

Model calculations for copper clusters on Au(111) surfaces

Using the semi-empirical embedded-atom method, the structure of small copper clusters on Au(111) surfaces has been investigated both by static and dynamic calculations. By varying the size of roughly circular clusters, the edge energy per atom is obtained; it agrees quite well with estimates based on experimental results. Small three-dimensional clusters tend to have the shape of a pyramid, whose sides are oriented in the directions of small surface energy. The presence of a cluster is found to distort the underlying lattice of adsorbed copper atoms. (C) 2002 Published by Elsevier Science B.V.

Compressed Metal Powders that Remain Superhydrophobic After Abrasion

Superhydrophobic “lotus effect” materials are typically not sufficiently robust for most real world applications because their small surface features are both easily damaged and vulnerable to fouling. Here, a method for preparing a new type of superhydrophobic (? > 162°) composite material by compression of superhydrophobic metal particles is reported. This material, which has no natural analogue, has low-surface-energy microstructures extending throughout its whole volume. Removing its outer layer by abrasion or cutting deep into it does not result in loss of superhydrophobicity because it merely exposes a fresh portion of the underlying superhydrophobic material. The high contact angle is therefore retained even after accidental damage, and vigorous abrasion can be used to restore hydrophobicity after fouling.

Electrochemical versus gas-phase oxidation of ru single-crystal surfaces

The electrochemical uptake of oxygen on a Ru(0001) electrode was investigated by electron diffraction, Auger spectroscopy, and cyclic voltammetry. An ordered (2 × 2)-O overlayer forms at a potential close to the hydrogen region. At +0.42 and +1.12 V vs Ag/AgCl, a (3 × 1) phase and a (1 × 1)-O phase, respectively, emerge. When the Ru electrode potential is maintained at +1.12 V for 2 min, RuO2 grows epitaxially with its (100) plane parallel to the Ru(0001) surface. In contrast to the RuO domains, the non-oxidized regions of the Ru electrode surface are flat. If, however, the electrode potential is increased to +1.98 V for 2 min, the remaining non-oxidized Ru area also becomes rough. These findings are compared with O overlayers and oxides on the Ru(0001) and Ru(101¯1) surfaces created by exposure to gaseous O under UHV conditions. On the other hand, gas-phase oxidation of the Ru(101¯0) surface leads to the formation of RuO with a (100) orientation. It is concluded that the difference in surface energy between RuO(110) and RuO(100) is quite small. RuO again grows epitaxially on Ru(0001), but with the (110) face oriented parallel to the Ru(0001) surface. The electrochemical oxidation of the Ru(0001) electrode surface proceeds via a 3-dimensional growth mechanism with a mean cluster size of 1.6 nm, whereas under UHV conditions, a 2-dimensional oxide film (1-2 nm thick) is epitaxially formed with an average domain size of 20 µm. © 2000 American Chemical Society.

Physical origin of the high reactivity of subsurface hydrogen in catalytic hydrogenation

In the catalytic hydrogenation of hydrocarbons, subsurface hydrogen is known experimentally to be much more reactive than surface hydrogen. We use density functional theory to identify low-energy pathways for the hydrogenation of methyl adsorbed on Ni(111) by surface and subsurface hydrogen. The metastability of subsurface hydrogen with respect to chemisorbed hydrogen is mainly responsible for the low activation barrier for subsurface reactions. (C) 1999 American Institute of Physics.