The vibrational spectra of water cluster molecules on Pt(111) surface at 20 K

Enhanced oxygen reduction reaction on caffeine-modified platinum single-crystal electrodes

Enhancing the activity of the oxygen reduction reaction (ORR) is crucial for fuel cell development, and hydrophobic species are known to increase the ORR activity. This paper reports that caffeine enhanced the specific ORR activity of Pt(111) 11-fold compared to that without caffeine in a 0.1 M HClO4 aqueous solution. Moreover, caffeine increased the ORR activity of Pt(110) 2.5-fold; however, the activity of Pt(100) was unaffected. The infrared (IR) band of PtOH (blocking species of the ORR) decreased for all the surfaces. Caffeine was adsorbed with its molecular plane perpendicular to the Pt(111) and Pt(110) surfaces and tilted relative to the Pt(100) surface. Thus, the effects of caffeine on the ORR activity can be rationalized by a decrease in PtOH coverage and the difference in adsorption geometry of caffeine.

Electrochemical hydrogen evolution on Pt-based catalysts from a theoretical perspective

Hydrogen evolution reaction (HER) by splitting water is a key technology toward a clean energy society, where Pt-based catalysts were long known to have the highest activity under acidic electrochemical conditions but suffer from high cost and poor stability. Here, we overview the current status of Pt-catalyzed HER from a theoretical perspective, focusing on the methodology development of electrochemistry simulation, catalytic mechanism, and catalyst stability. Recent developments in theoretical methods for studying electrochemistry are introduced, elaborating on how they describe solid-liquid interface reactions under electrochemical potentials. The HER mechanism, the reaction kinetics, and the reaction sites on Pt are then summarized, which provides an atomic-level picture of Pt catalyst surface dynamics under reaction conditions. Finally, state-of-the-art experimental solutions to improve catalyst stability are also introduced, which illustrates the significance of fundamental understandings in the new catalyst design.

Ordered Water Layer on the Macroscopically Hydrophobic Fluorinated Polymer Surface and Its Ultrafast Vibrational Dynamics

Hydrophobic-like water monolayers have been predicted at the metal and some polar surfaces by theoretical simulations. However, direct experimental evidence for the presence of this water layer at surfaces, particularly at biomolecule and polymer surfaces, is yet to be validated at room temperature. Here we observe experimentally that an ordered molecular water layer is present at the hydrophobic fluorinated polymer such as polytetrafluoroethylene (PTFE) surface by using sum frequency generation vibrational spectroscopy. The macroscopic hydrophobicity of PTFE surface is actually hydrophilic at the molecular level. The macroscopically hydrophobic character of PTFE is indeed resulting from the hydrophobicity of the ordered two-dimension (2D) water layer, in which cyclic water tetramer structure is found. The water layer at humidity of ≤40% has a vibrational relaxation time of 550 ± 60 fs. The vibrational relaxation time in the frequency range of 3200-3400 cm^-1 shows remarkable difference from the interfacial water at the air/H₂O interface and the lipid/H₂O interface. No discernible frequency dependence of the vibrational relaxation time is observed, indicating the homogeneous dynamics of OH groups in the water layer. These insights into the water layer at the macroscopically hydrophobic surface may contribute to a better understanding of the hydrophobic interaction and interfacial water dynamics.

Orientational ordering in heteroepitaxial water ice on metal surfaces

Sum frequency generation spectroscopy uncovers the orientational ordering in crystalline ice films of water grown on Pt(111) and Rh(111).

Vibrational spectroscopic study of pH dependent solvation at a Ge(100)-water interface during an electrode potential triggered surface termination transition

The charge-dependent structure of interfacial water at the n-Ge(100)-aqueous perchlorate interface was studied by controlling the electrode potential. Specifically, a joint attenuated total reflection infrared spectroscopy and electrochemical experiment was used in 0.1M NaClO₄ at pH ≈ 1-10. The germanium surface transformation to an H-terminated surface followed the thermodynamic Nernstian pH dependence and was observed throughout the entire pH range. A singular value decomposition-based spectra deconvolution technique coupled to a sigmoidal transition model for the potential dependence of the main components in the spectra shows the surface transformation to be a two-stage process. The first stage was observed together with the first appearance of Ge-H stretching modes in the spectra and is attributed to the formation of a mixed surface termination. This transition was reversible. The second stage occurs at potentials ≈0.1-0.3 V negative of the first one, shows a hysteresis in potential, and is attributed to the formation of a surface with maximum Ge-H coverage. During the surface transformation, the surface becomes hydrophobic, and an effective desolvation layer, a "hydrophobic gap," developed with a thickness ≈1-3 Å. The largest thickness was observed near neutral pH. Interfacial water IR spectra show a loss of strongly hydrogen-bound water molecules compared to bulk water after the surface transformation, and the appearance of "free," non-hydrogen bound OH groups, throughout the entire pH range. Near neutral pH at negative electrode potentials, large changes at wavenumbers below 1000 cm^-1 were observed. Librational modes of water contribute to the observed changes, indicating large changes in the water structure.

Electrode potential dependent desolvation and resolvation of germanium(100) in contact with aqueous perchlorate electrolytes

The electrode potential dependence of the hydration layer on an n-Ge(100) surface was studied by a combination of in situ and operando electrochemical attenuated total reflection infrared (ATR-IR) spectroscopy and real space density functional theory (DFT) calculations. Constant-potential DFT calculations were coupled to a modified generalised Poisson-Boltzmann ion distribution model and applied within an ab initio molecular dynamics (AIMD) scheme. As a result, potential-dependent vibrational spectra of surface species and surface water were obtained, both experimentally and by simulations. The experimental spectra show increasing absorbance from the Ge-H stretching modes at negative potentials, which is associated with an increased negative difference absorbance of water-related OH modes. When the termination transition of germanium from OH to H termination occurs, the surface switches from hydrophilic to hydrophobic. This transition is fully reversible. During the switching, the interface water molecules are displaced from the surface forming a "hydrophobic gap". The gap thickness was experimentally estimated by a continuum electrodynamic model to be ≈2 Å. The calculations showed a shift in the centre of mass of the interface water by ≈0.9 Å due to the surface transformation. The resulting IR spectra of the interfacial water in contact with the hydrophobic Ge-H show an increased absorbance of free OH groups, and a decreased absorbance of strongly hydrogen bound water. Consequently, the surface transformation to a Ge-H terminated surface leads to a surface which is weakening the H-bond network of the interfacial water in contact.

Long-range influence of steps on water adsorption on clean and D-covered Pt surfaces

Water wets the D-covered Pt(111) surface (right), while it clusters at steps of D-covered Pt(533), (755), and (977) (left).

Theoretical vibrational Raman and surface-enhanced Raman scattering spectra of water interacting with silver clusters

In this study, we analyzed the Raman spectrum of a water molecule adsorbed on a cluster of 20 silver atoms, and the plasmonic electromagnetic effect of the silver surface was also considered to give a theoretical prediction of the surface-enhanced Raman scattering spectrum. The calculations were performed at the density functional theory (DFT) level by using both frozen and unfrozen silver clusters. Two different models were used to consider the plasmonic enhancement; one of them was a modified classical (dipole) model and the other was the coupled perturbed Hartree-Fock method with excitation frequencies obtained from time-dependent DFT calculations and with proper detuning of these frequencies. The importance of small geometrical distortions of the silver surface in the orientation of the adsorbed water was shown. Moreover, it was shown how the symmetry of the transition dipole moment and the symmetry of the vibrational modes influence the Raman intensities of the SERS spectrum.

Adsorption of water dimer on platinum(111): identification of the -OH···Pt hydrogen bond

The fundamental structure of an isolated water dimer on Pt(111) was determined by means of a spectroscopic method using scanning tunneling microscopy (STM) and density functional theory (DFT) calculations. Two water molecules on adjacent atop sites form a dimer through a hydrogen bond, and they rotate even at a substrate temperature of 5 K. Action spectroscopy using STM (STM-AS) for water dimer hopping allows us to obtain the vibrational spectrum of a single water dimer on Pt(111). Comparisons between the experiments and theory show that one of the OH groups of the acceptor water molecule points toward the surface to form an -OH···Pt hydrogen bond.

Material functions of liquid n-hexadecane under steady shear via nonequilibrium molecular dynamics simulations: temperature, pressure, and density effects

Computer experiments of rheology regarding the effects of temperature (T), pressure (P), and density (rho) on steady shear flow material functions, which include viscosity (eta) and first and second normal stress coefficients (psi(1) and psi(2)) depending on shear rate (gamma), have been conducted via nonequilibrium molecular dynamics simulations for liquid n-hexadecane. Straightforwardly, using both characteristic values of a zero-shear-rate viscosity and critical shear rate, eta-gamma flow curves are well normalized to achieve the temperature-, pressure-, and density-invariant master curves, which can be formulary described by the Carreau-Yasuda rheological constitutive equation. Variations in the rate of shear thinning, obviously exhibiting in eta-gamma, psi(1)-gamma, and -psi(2)-gamma relationships, under different T, P, and rho values, are concretely revealed through the power-law model's exponent. More importantly, at low shear rates, the fluid explicitly possesses Newtonian fluidic characteristics according to both manifestations; first and second normal stress differences decay to near zero, while nonequilibrium states are close to equilibrium ones. Significantly, the tendency to vary of the degree of shear thinning in rheology is qualitatively contrary to that of shear dilatancy in thermodynamics. In addition, a convergent transition point is evidently observed in the -psi(2)/psi(1)-gamma curves undergoing dramatic variations, which should be associated with shear dilatancy, as addressed analytically.

Min-map bias Monte Carlo for chain molecules: biased Monte Carlo sampling based on bijective minimum-to-minimum mapping

A novel Monte Carlo (MC) simulation scheme based on Theodorou's bijective mapping strategy [D. N. Theodorou, J. Chem. Phys. 124, 034109 (2006)] is introduced. This min-map bias Monte Carlo acts in combination with any other proper, bare MC. It carries over the bare MC move from the original configuration space Omega(0), where trial move acceptance may be low, to a different configuration space, Omega(1), where acceptance is higher. The bare MC move is then performed in Omega(1) and the resulting configuration is finally mapped back to Omega(0). Mappings between Omega(0) and Omega(1) entail weighted selection of trial configurations, the bias of which is subsequently removed in the overall acceptance criterion. The new method is applied, in conjunction with continuum configurational bias as bare MC scheme, to the simulation of explicit hydrogen linear alkanes in the canonical ensemble. Min-map bias MC is found to alleviate the pervasive problem of very low acceptance rates encountered when using an explicit molecular description.

Structure and internal dynamics of a side chain liquid crystalline polymer in various phases by molecular dynamics simulations: a step towards coarse graining

Side chain liquid crystalline polymer with relatively long spacer was modeled on a semiatomistic level and studied in different liquid crystalline phases with the aid of molecular dynamics simulations. Well equilibrated isotropic, polydomain smectic and monodomain smectic phases were studied for their structural and dynamic properties. Particular emphasis was given to the analysis on a coarse-grained level, where backbones, side chains, and mesogens were considered in terms of their equivalent ellipsoids. The authors found that the liquid crystalline phase had a minor influence on the metrics of these objects but affected essentially their translational and orientational order. In the monodomain smectic phase, mesogens, backbones, and side chains are confined spatially. Their diffusion and shape dynamics are frozen along the mesogen director (the one-dimensional solidification) and the reorientation times increase by one to one-and-half orders of magnitude. In this phase, besides obvious orientational order of mesogens and side chains, a stable detectable order of the backbones was also observed. The backbone director is confined in the plane perpendicular to the mesogen director and constantly changes its orientation within this plane. The backbone diffusion in these planes is of the same range as in the polydomain smectic phase at the same temperature. A detailed analysis of the process of field-induced growth of the smectic phase was performed. The study revealed properties of liquid crystalline polymers that may enable their future fully coarse-grained modeling.

Dynamics in atomistic simulations of phospholipid membranes: Nuclear magnetic resonance relaxation rates and lateral diffusion

It is shown that a long, near microsecond, atomistic simulation can shed some light upon the dynamical processes occurring in a lipid bilayer. The analysis focuses on reorientational dynamics of the chains and lateral diffusion of lipids. It is shown that the reorientational correlation functions exhibits an algebraic decay (rather than exponential) for several orders of magnitude in time. The calculated nuclear magnetic resonance relaxation rates agree with experiments for carbons at the C7 position while there are some differences for C3. Lateral diffusion can be divided into two stages. In a first stage occurring at short times, t<5 ns, the center of mass of the lipid moves due to conformational changes of the chains while the headgroup position remains relatively fixed. In this stage, the center of mass can move up to approximately 0.8 nm. The fitted short-time diffusion coefficient is D(1)=13 x 10(-7) cm(2) s(-1) On a longer time scale, the diffusion coefficient becomes D(2)=0.79 x 10(-7) cm(2) s(-1).

Constrained dynamics and extraction of normal modes fromab initiomolecular dynamics: Application to ammonia

We present a methodology for extracting phonon data from ab initio Born-Oppenheimer molecular dynamics calculations of molecular crystals. Conventional ab initio phonon methods based on perturbations are difficult to apply to lattice modes because the perturbation energy is dominated by intramolecular modes. We use constrained molecular dynamics to eliminate the effect of bond bends and stretches and then show how trajectories can be used to isolate and define in particular, the eigenvalues and eigenvectors of modes irrespective of their symmetry or wave vector. This is done by k-point and frequency filtering and projection onto plane wave states. The method is applied to crystalline ammonia: the constrained molecular dynamics allows a significant speed-up without affecting structural or vibrational modes. All Gamma point lattice modes are isolated: the frequencies are in agreement with previous studies; however, the mode assignments are different.

On the analysis of conformational dynamics in polymers with several rotational isomers

The ability of different correlation functions to shed some light onto the conformational dynamics of an amorphous polymer has been analyzed. The study has been performed on a polyethylene model polymer, which has been simulated at decreasing temperatures towards its glass transition, via the molecular dynamics technique. Three rotational isomers are allowed by the considered torsional potential. The correlation times associated with the evaluated transition rates have shown to be Arrhenius in nature, with activation energies resulting basically from internal rotation barriers. Overall torsional autocorrelation functions have been calculated. We have observed that they are dominated by slow events. Alternatively, a set of torsional autocorrelation functions associated with every isomeric state has been evaluated. Stretched exponential fits lead to correlation times that display Vogel-Fulcher temperature dependence.

Free energy of a trans-membrane pore calculated from atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations of a lipid bilayer were performed to calculate the free energy of a trans-membrane pore as a function of its radius. The free energy was calculated as a function of a reaction coordinate using a potential of mean constraint force. The pore radius was then calculated from the reaction coordinate using Monte Carlo particle insertions. The main characteristics of the free energy that comes out of the simulations are a quadratic shape for a radius less than about 0.3 nm, a linear shape for larger radii than this, and a rather abrupt change without local minima or maxima between the two regions. In the outer region, a line tension can be calculated, which is consistent with the experimentally measured values. Further, this line tension can be rationalized and understood in terms of the energetic cost for deforming a part of the lipid bilayer into a hydrophilic pore. The region with small radii can be described and understood in terms of statistical mechanics of density fluctuations. In the region of crossover between a quadratic and linear free energy there was some hysteresis associated with filling and evacuation of the pore with water. The metastable prepore state hypothesized to interpret the experiments was not observed in this region.

Free-energy analysis of solubilization in micelle

A statistical-mechanical treatment of the solubilization in micelle is presented in combination with molecular simulation. The micellar solution is viewed as an inhomogeneous and partially finite, mixed solvent system, and the method of energy representation is employed to evaluate the free-energy change for insertion of a solute into the micelle inside with a realistic set of potential functions. Methane, benzene, and ethylbenzene are adopted as model hydrophobic solutes to analyze the solubilization in sodium dodecyl sulfate micelle. It is shown that these solutes are more favorably located within the micelle than in bulk water and that the affinity to the micelle inside is stronger for benzene and ethylbenzene than for methane. The micellar system is then divided into the hydrophobic core, the head-group region in contact with water, and the aqueous region outside the micelle to assess the relative importance of each region in the solubilization. In support of the pseudophase model, the aqueous region is found to be unimportant to determine the extent of solubilization. The contribution from the hydrophobic-core region is shown to be dominant for benzene and ethylbenzene, while an appreciable contribution from the head-group region is observed for methane. The methodology presented is not restricted to the binding of a molecule to micelle, and will be useful in treating the binding to such nanoscale structures as protein and membrane.

Molecular dynamics study of carbon nanotube oscillators revisited

We performed molecular dynamics simulation of double walled carbon nanotube (DWCNT) oscillators under constant energy and constant temperatures with various commensurations and nanotube lengths. We clarify and resolve questions and differences raised by previous simulation results of similar systems. At constant energy, sustained oscillation is available for a wide range of initial temperatures. But low initial temperature is advantageous for DWCNTs to sustain oscillation under constant energy. We observed sustained oscillation at constant energy for both commensurate and incommensurate DWCNTs. On the other hand, under constant temperatures, both high and low temperatures are disadvantageous to sustain DWCNT oscillations. At constant low temperature, neither commensurate nor incommensurate DWCNTs can maintain oscillation. At appropriate constant temperatures, the oscillatory behavior of incommensurate nanotubes is much more sustained than that of commensurate tubes. The oscillatory frequency of DWCNTs depends significantly on the length of tubes. The initial oscillatory frequency is inversely proportional to the DWCNT lengths. The oscillation frequency of DWCNTs is insensitive to the initial temperatures at constant energy, but slightly dependent on the temperature at constant temperatures.

The structure and crystallization of thin water films on Pt(111)

When water is adsorbed on Pt(111) above 135 K several different ice structures crystallize, depending on the thickness of the ice layer. At low coverage water forms extended islands of ice with a (square root(37) x square root(37))R25(o) unit cell, which compresses as the monolayer saturates to form a (square root(39) x square root(39))R16(o) structure. The square root(39) low-energy electron diffraction (LEED) pattern becomes more intense as the second layer grows, remaining bright for films up of 10-15 layers and then fading and disappearing for films more than ca. 40 layers thick. The ice multilayer consists of an ordered square root(39) wetting layer, on which ice grows as a crystalline film which progressively loses its registry to the wetting layer. Ice films more than ca. 50 layers thick develop a hexagonal LEED pattern, the entire film and wetting layer reorienting to form an incommensurate bulk ice. These changes are reflected in the vibrational spectra which show changes in line shape and intensity associated with the different ice structures. Thin amorphous solid water films crystallize to form the same phases observed during growth, implying that these structures are thermodynamically stable and not kinetic phases formed during growth. The change from a square root(39) registry to incommensurate bulk ice at ca. 50 layers is associated with a change in crystallization kinetics from nucleation at the Pt(111) interface in thin films to nucleation of incommensurate bulk ice in amorphous solid water films more than 50 layers thick.

Water Adsorption on Rh(111) at 20 K: From Monomer to Bulk Amorphous Ice

The adsorption of water (D(2)O) molecules on Rh(111) at 20 K was investigated using infrared reflection absorption spectroscopy (IRAS). At the initial stage of adsorption, water molecules exist as monomers on Rh(111). With increasing water coverage, monomers aggregate into dimers, larger clusters (n = 3-6), and two-dimensional (2D) islands. Further exposure of water molecules leads to the formation of three-dimensional (3D) water islands and finally to a bulk amorphous ice layer. Upon heating, the monomer and dimer species thermally migrate on the surface and aggregate to form larger clusters and 2D islands. Based on the temperature dependence of OD stretching peaks, we succeeded in distinguishing water molecules inside 2D islands from those at the edge of 2D islands. From the comparison with the previous vibrational spectra of water clusters on other metal surfaces, we conclude that the number of water molecules at the edge of 2D islands is comparable with that of water molecules inside 2D islands on the Rh(111) surface at 20 K. This indicates that the surface migration of water molecules on Rh(111) is hindered as compared with the cases on Pt(111) and Ni(111) and thus the size of 2D islands on Rh(111) is relatively small.

Different surface chemistries of water on Ru[0001]: from monomer adsorption to partially dissociated bilayers

Density functional theory has been used to perform a comparative theoretical study of the adsorption and dissociation of H(2)O monomers and icelike bilayers on Ru[0001]. H(2)O monomers bind preferentially at atop sites with an adsorption energy of approximately 0.4 eV/H(2)O. The main bonding interaction is through the H(2)O 1b(1) molecular orbital which mixes with Ru d(z)2 states. The lower-lying set of H(2)O molecules in an intact H(2)O bilayer bond in a similar fashion; the high-lying H(2)O molecules, however, do not bond directly with the surface, rather they are held in place through H bonding. The H(2)O adsorption energy in intact bilayers is approximately 0.6 eV/H(2)O and we estimate that H bonding accounts for approximately 70% of this. In agreement with Feibelman (Science 2002, 295, 99) we find that a partially dissociated OH + H(2)O overlayer is energetically favored over pure intact H(2)O bilayers on the surface. The barrier for the dissociation of a chemisorbed H(2)O monomer is 0.8 eV, whereas the barrier to dissociate a H(2)O incorporated in a bilayer is just 0.5 eV.

Structure and bonding of water on Pt(111)

We address the adsorption of water on Pt(111) using x-ray absorption, x-ray emission, and x-ray photoelectron spectroscopy along with calculations in the framework of density functional theory. Using the direct relationship between the electronic structure and adsorbate geometry, we show that in the first layer all the molecules bind directly to the surface and to each other through the in-layer H bonds without dissociation, creating a nearly flat overlayer. The water molecules are adsorbed through alternating metal-oxygen (M-O) and metal-hydrogen (M-HO) bonds.

Water diffusion and clustering on Pd(111)

The adsorption, diffusion, and clustering of water molecules on a Pd(111) surface were studied by scanning tunneling microscopy. At 40 kelvin, low-coverage water adsorbs in the form of isolated molecules, which diffuse by hopping to nearest neighbor sites. Upon collision, they form first dimers, then trimers, tetramers, and so on. The mobility of these species increased by several orders of magnitude when dimers, trimers, and tetramers formed, and decreased again when the cluster contained five or more molecules. Cyclic hexamers were found to be particularly stable. They grow with further exposure to form a commensurate hexagonal honeycomb structure relative to the Pd(111) substrate. These observations illustrate the change in relative strength between intermolecular hydrogen bonds and molecule-substrate bonds as a function of water cluster size, the key property that determines the wetting properties of materials.