Adsorption of Water onto SrTiO3 from Periodic Møller-Plesset Second-Order Perturbation Theory.

Adsorption of water onto metal oxide surfaces is a long-standing problem motivated by relevance to many promising technological applications. In this work, we compute the adsorption energy of water on SrTiO3 using periodic Møller-Plesset second-order perturbation theory (MP2). We compare our MP2 results to density functional and hybrid density functional theory calculations with and without the widely used D3 dispersion correction. The MP2 ground-state adsorption energy of water on SrTiO3 (001) at one monolayer coverage is 0.9 eV on the TiO2 termination in the molecular configuration and 0.6 eV in the dissociative configuration, the corresponding results on the SrO termination being 0.9 eV for both modes of adsorption. These results are reproduced well by the PBE and PBE0 exchange-correlation functionals. Correcting for dispersion effects through the D3 dispersion correction leads to significantly overestimated adsorption energies for both PBE and PBE0 with respect to MP2. The D3 correction also fails to capture the difference in electron correlation between the molecular and dissociative adsorption states, similarly to the optB86b van der Waals density functional.

Diacetylene polymerization on a bulk insulator surface.

Molecular electronics has great potential to surpass known limitations in conventional silicon-based technologies. The development of molecular electronics devices requires reliable strategies for connecting functional molecules by wire-like structures. To this end, diacetylene polymerization has been discussed as a very promising approach for contacting single molecules with a conductive polymer chain. A major challenge for future device fabrication is transferring this method to bulk insulator surfaces, which are mandatory to decouple the electronic structure of the functional molecules from the support surface. Here, we provide experimental evidence for diacetylene polymerization of 3,3'-(1,3-butadiyne-1,4-diyl)bisbenzoic acid precursors on the (10.4) surface of calcite, a bulk insulator with a band gap of around 6 eV. When deposited on the surface held at room temperature, ordered islands with a (1 × 3) superstructure are observed using dynamic atomic force microscopy. A distinct change is revealed upon heating the substrate to 485 K. After heating, molecular stripes with a characteristic inner structure are formed that excellently match the expected diacetylene polymer chains in appearance and repeat distance. The corresponding density functional theory computations reveal molecular-level bonding patterns of both the (1 × 3) superstructure and the formed striped structure, confirming the assignment of on-surface diacetylene polymerization. Transferring the concept of using diacetylene polymerization for creating conductive connections to bulk insulator surfaces paves the way towards application-relevant systems for future molecular electronic devices.

The interface of SrTiO3 and H2O from density functional theory molecular dynamics.

We use dispersion-corrected density functional theory molecular dynamics simulations to predict the ionic, electronic and vibrational properties of the SrTiO3/H2O solid-liquid interface. Approximately 50% of surface oxygens on the planar SrO termination are hydroxylated at all studied levels of water coverage, the corresponding number being 15% for the planar TiO2 termination and 5% on the stepped TiO2-terminated surface. The lateral ordering of the hydration structure is largely controlled by covalent-like surface cation to H2O bonding and surface corrugation. We find a featureless electronic density of states in and around the band gap energy region at the solid-liquid interface. The vibrational spectrum indicates redshifting of the O-H stretching band due to surface-to-liquid hydrogen bonding and blueshifting due to high-frequency stretching vibrations of OH fragments within the liquid, as well as strong suppression of the OH stretching band on the stepped surface. We find highly varying rates of proton transfer above different SrTiO3 surfaces, owing to differences in hydrogen bond strength and the degree of dissociation of incident water. Trends in proton dynamics and the mode of H2O adsorption among studied surfaces can be explained by the differential ionicity of the Ti-O and Sr-O bonds in the SrTiO3 crystal.

Boundary lubrication with a liquid crystal monolayer.

We study boundary lubrication characteristics of a liquid crystal (LC) monolayer sheared between two crystalline surfaces by nonequilibrium molecular dynamics simulations, using a simplified rigid bead-necklace model of the LC molecules. We consider LC monolayers confined by surfaces with three different atomic structures, subject to different shearing velocities, thus approximating a wide variety of materials and driving conditions. The time dependence of the friction force is studied and correlated with that of the orientational order exhibited by the LC molecules, arising from the competition between the effect of the structure of the confining surfaces and that of the imposed sliding direction. We show that the observed stick-slip events for low shear rates involve order-disorder transitions, and that the LC monolayer no longer has enough time to reorder at high shear rates, resulting in a smooth sliding regime. An irregular stick-slip phase between the regular stick-slip and smooth sliding is observed for intermediate shear rates regardless of the surface structure.

Non-contact atomic force microscopy study of hydroxyl groups on the spinel MgAl2O4(100) surface.

Atom-resolved non-contact atomic force microscopy (NC-AFM) studies of the magnesium aluminate (MgAl(2)O(4)) surface have revealed that, contrary to expectations, the (100) surface is terminated by an aluminum and oxygen layer. Theoretical studies have suggested that hydrogen plays a strong role in stabilizing this surface through the formation of surface hydroxyl groups, but the previous studies did not discuss in depth the possible H configurations, the diffusion behaviour of hydrogen atoms and how the signature of adsorbed H is reflected in atom-resolved NC-AFM images. In this work, we combine first principles calculations with simulated and experimental NC-AFM images to investigate the role of hydrogen on the MgAl(2)O(4)(100) surface. By means of surface energy calculations based on density functional theory, we show that the presence of hydrogen adsorbed on the surface as hydroxyl groups is strongly predicted by surface stability considerations at all relevant partial pressures of H(2) and O(2). We then address the question of how such adsorbed hydrogen atoms are reflected in simulated NC-AFM images for the most stable surface hydroxyl groups, and compare with experimental atom-resolved NC-AFM data. In the appendices we provide details of the methods used to simulate NC-AFM using first principles methods and a virtual AFM.

Topography and work function measurements of thin MgO(001) films on Ag(001) by nc-AFM and KPFM.

The surface topography and local surface work function of ultrathin MgO(001) films on Ag(001) have been studied by noncontact atomic force microscopy (nc-AFM) and Kelvin probe force microscopy (KPFM). First principles calculations have been used to explain the contrast formation of nc-AFM images. In agreement with literature, thin MgO films grow in islands with a quasi rectangular shape. Contrary to alkali halide films supported on metal surfaces, where the island heights can be correctly measured, small MgO islands are either imaged as depressions or elevations depending on the electrostatic potential of the tip apex. Correct island heights therefore cannot be given without knowing the precise contrast formation discussed in this paper. KPFM shows a silver work function which is reduced by the MgO islands. The values for the work function differences for one and two layer thin films are -1.1 and -1.4 eV, respectively, in good agreement with recent calculations and experiments.

Embedding transition-metal atoms in graphene: structure, bonding, and magnetism.

We present a density-functional-theory study of transition-metal atoms (Sc-Zn, Pt, and Au) embedded in single and double vacancies (SV and DV) in a graphene sheet. We show that for most metals, the bonding is strong and the metal-vacancy complexes exhibit interesting magnetic behavior. In particular, an Fe atom on a SV is not magnetic, while the Fe@DV complex has a high magnetic moment. Surprisingly, Au and Cu atoms at SV are magnetic. Both bond strengths and magnetic moments can be understood within a simple local-orbital picture, involving carbon sp(2) hybrids and the metal spd orbitals. We further calculate the barriers for impurity-atom migration, and they agree well with available experimental data. We discuss the experimental realization of such systems in the context of spintronics and nanocatalysis.

Finite-element implementation for electron transport in nanostructures.

We have modeled transport properties of nanostructures using Green's-function method within the framework of the density-functional theory. The scheme is computationally demanding, so numerical methods have to be chosen carefully. A typical solution to the numerical burden is to use a special basis-function set, which is tailored to the problem in question, for example, the atomic-orbital basis. In this paper we present our solution to the problem. We have used the finite-element method with a hierarchical high-order polynomial basis, the so-called p elements. This method allows the discretation error to be controlled in a systematic way. The p elements work so efficiently that they can be used to solve interesting nanosystems described by nonlocal pseudopotentials. We demonstrate the potential of the implementation with two different systems. As a test system a simple Na-atom chain between two leads is modeled and the results are compared with several previous calculations. Secondly, we consider a thin hafnium dioxide (HfO2) layer on a silicon surface as a model for a gate structure of the next generation of microelectronics.

Reactions and clustering of water with silica surface.

The interaction between silica surface and water is an important topic in geophysics and materials science, yet little is known about the reaction process. In this study we use first-principles molecular dynamics to simulate the hydrolysis process of silica surface using large cluster models. We find that a single water molecule is stable near the surface but can easily dissociate at three-coordinated silicon atom defect sites in the presence of other water molecules. These extra molecules provide a mechanism for hydrogen transfer from the original water molecule, hence catalyzing the reaction. The two-coordinated silicon atom is inert to the water molecule, and water clusters up to pentamer could be stably adsorbed at this site at room temperature.

Probing organic layers on the TiO2(110) surface.

In this work, we use first principles simulations to provide features of the dynamic scanning force microscopy imaging of adsorbed organic layers on insulating surfaces. We consider monolayers of formic (HCOOH) and acetic (CH(3)COOH) acid and a mixed layer of acetic and trifluoroacetic acids (CF(3)COOH) on the TiO(2)(110) surface and study their interaction with a silicon dangling bond tip. The results demonstrate that the silicon tip interacts more strongly with the substrate and the COO(-) group than the adsorbed acid headgroups, and, therefore, molecules would appear dark in images. The pattern of contrast and apparent height of molecules is determined by the repulsion between the tip and the molecular headgroups and by significant deformation of the monolayer and individual molecules. The height of the molecule on the surface and the size of the headgroup play a large role in determining access of the tip to the substrate and, hence, the contrast in images. Direct imaging of the molecules themselves could be obtained by providing a functionalized tip with attraction to the molecular headgroups, for example, a positive potential tip.

Irradiation-induced magnetism in graphite: a density functional study.

Recent experiments indicate that proton irradiation triggers ferromagnetism in originally nonmagnetic graphite samples while He ion bombardment has a much smaller effect. To understand the origin of irradiation-induced magnetism, we have performed spin-polarized density functional theory calculations of the magnetic properties of the defects which are most likely to appear under irradiation vacancies and vacancy-hydrogen complexes. Both defects are magnetic, but as for the latter we find that H adsorption on one of the vacancy dangling bonds gives rise to a magnetic moment double that of the naked vacancy. We show that for small irradiation doses vacancy-hydrogen complexes result in a macroscopic magnetic signal which agrees well with the experimental values.

Adsorption of acetic and trifluoroacetic acid on the TiO2(110) surface.

We use the first-principles static and dynamic simulations to study the adsorption of acetic (CH(3)COOH) and trifluoroacetic (CF(3)COOH) acid on the TiO(2)(110) surface. The most favorable adsorption for both molecules is a dissociative process, which results in the two oxygens of the carboxylate ion bonding to in-plane titanium atoms in the surface. The remaining proton then bonds to a bridging oxygen site, forming a hydroxyl group. We further show that, by comparing the calculated dipoles of the molecules on the surface, it is possible to understand the difference in contrast over the acetate and trifluoroacetate molecules in the atomically resolved noncontact atomic force microscopy images.

Interaction of silicon dangling bonds with insulating surfaces.

We use first principles density functional theory calculations to study the interaction of a model dangling bond silicon tip with the surfaces of CaF2, Al2O3, TiO2, and MgO. In each case the strongest interaction is with the highest anions in the surface. We show that this is due to the onset of chemical bonding with the surface anions, which can be controlled by an electric field across the system. Combining our results and previous studies on semiconductor surfaces suggests that using dangling bond Si tips can provide immediate identification of surface species in atomically resolved noncontact atomic force microscopy and facilitate selective measurements of short-range interactions with surface sites.

Magnetic properties and diffusion of adatoms on a graphene sheet.

We use ab initio methods to calculate the properties of adatom defects on a graphite surface. By applying a full spin-polarized description to the system we demonstrate that these defects have a magnetic moment of about 0.5micro(B) and also calculate its role in diffusion over the surface. The magnetic nature of these intrinsic carbon defects suggests that it is important to understand their role in the recently observed magnetism in pure carbon systems.

Mechanism of interstitial oxygen diffusion in hafnia.

We have performed density functional calculations of oxygen incorporation and diffusion in monoclinic hafnia (HfO2) for a range of oxygen charge states. The calculations demonstrate that oxygen favors atomic incorporation and that O2- is the most stable species. We find that oxygen interstitials diffuse via exchange with lattice oxygen sites in hafnia, and that O- species have the smallest diffusion barrier.

Unambiguous interpretation of atomically resolved force microscopy images of an insulator.

The (111) surface of CaF2 was imaged with dynamic mode scanning force microscopy and modeled using atomistic simulation. Both experiment and theory showed a clear triangular contrast pattern in images, and theory demonstrated that the contrast pattern is due to the interaction of a positive electrostatic potential tip with fluorine ions in the two topmost surface layers. We find a good agreement of position and relative height of scan line features between theory and experiment and thus establish for the first time an unambiguous identification of sublattices of an insulator imaged by force microscopy.

Reading achievement by learning disabled students in resource and regular classes.

K-TEA Comprehensive Reading scores of 34 elementary boys in either resource rooms or regular settings were compared. The boys were identified as learning disabled in reading. They were pretested at the beginning of the school year and posttested at the end. Treatment was one year of daily instruction in reading provided by six teachers in resource setting and six teachers in regular settings. K-TEA Reading Decoding and Reading Comprehension scores, separately compared in 2 x 2 repeated-measures analysis of variance, were not significantly different.