Revealing the hidden structure of disordered materials by parameterizing their local structural manifold.

Durable interest in developing a framework for the detailed structure of glassy materials has produced numerous structural descriptors that trade off between general applicability and interpretability. However, none approach the combination of simplicity and wide-ranging predictive power of the lattice-grain-defect framework for crystalline materials. Working from the hypothesis that the local atomic environments of a glassy material are constrained by enthalpy minimization to a low-dimensional manifold in atomic coordinate space, we develop a generalized distance function, the Gaussian Integral Inner Product (GIIP) distance, in connection with agglomerative clustering and diffusion maps, to parameterize that manifold. Applying this approach to a two-dimensional model crystal and a three-dimensional binary model metallic glass results in parameters interpretable as coordination number, composition, volumetric strain, and local symmetry. In particular, we show that a more slowly quenched glass has a higher degree of local tetrahedral symmetry at the expense of cyclic symmetry. While these descriptors require post-hoc interpretation, they minimize bias rooted in crystalline materials science and illuminate a range of structural trends that might otherwise be missed.

Water adsorption on MoS2 under realistic atmosphere conditions and impacts on tribology.

Molybdenum disulfide (MoS2) is a 2D material widely used as a dry lubricant. However, exposure to water and oxygen is known to reduce its effectiveness, and therefore an understanding of the uptake of water is important information for mitigating these effects. Here we use grand canonical Monte Carlo simulations to rigorously study water adsorption on MoS2 surfaces and edges with different concentrations of defects under realistic atmospheric conditions (i.e. various temperatures and humidity levels). We find that the amount of water adsorbed depends strongly on the number of defects. Simulations indicate that defect sites are generally saturated with water even at low ppm levels of humidity. Water binds strongly to S vacancies on interlamellar surfaces, but generally only one water molecule can fit on each of these sites. Defects on surfaces or edges of lamellae also strongly attract water molecules that then nucleate small clusters of water bonded via hydrogen bonding. We demonstrate that water preferentially binds to surface defects, but once those are saturated at a critical humidity level of about 500-1000 ppm water, water binds to edge sites where it negatively impacts the tribological performance of MoS2.

Interactions of Water with Pristine and Defective MoS2.

Molybdenum disulfide (MoS2) is a lamellar solid lubricant often used in aerospace applications because of its extremely low friction coefficient (∼0.01) in inert environments. The lubrication performance of MoS2 is significantly impaired by exposure to even small amounts of water and oxygen, and the mechanisms behind this remain poorly understood. We use density functional theory calculations to study the binding of water on MoS2 sheets with and without defects. In general, we find that pristine MoS2 is slightly hydrophilic but that defects greatly increase the binding affinity for water. Intercalated water disrupts the crystal structure of bulk MoS2 due to the limited space between lamellae (∼3.4 Å), and this leads to generally unfavorable adsorption, except in the cases where water molecules are located on the sites of sulfur vacancies. We also find that water adsorption is more favorable directly below a surface layer of MoS2 compared to in the bulk.

Formation of Coherent 1H-1T Heterostructures in Single-Layer MoS2 on Au(111).

Heterojunctions of semiconductors and metals are the fundamental building blocks of modern electronics. Coherent heterostructures between dissimilar materials can be achieved by composition, doping, or heteroepitaxy of chemically different elements. Here, we report the formation of coherent single-layer 1H-1T MoS2 heterostructures by mechanical exfoliation on Au(111), which are chemically homogeneous with matched lattices but show electronically distinct semiconducting (1H phase) and metallic (1T phase) character, with the formation of these heterojunctions attributed to a combination of lattice strain and charge transfer. The exfoliation approach employed is free of tape residues usually found in many exfoliation methods and yields single-layer MoS2 with millimeter (mm) size on the Au surface. Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), atomic force microscopy (AFM), scanning tunneling microscopy (STM), and scanning tunneling spectroscopy (STS) have collectively been employed to elucidate the structural and electronic properties of MoS2 monolayers on Au substrates. Bubbles in the MoS2 formed by the trapping of ambient adsorbates beneath the single layer during deposition, have also been observed and characterized. Our work here provides a basis to produce two-dimensional heterostructures which represent potential candidates for future electronic devices.

Evidence of Inverse Hall-Petch Behavior and Low Friction and Wear in High Entropy Alloys.

We present evidence of inverse Hall-Petch behavior for a single-phase high entropy alloy (CoCrFeMnNi) in ultra-high vacuum and show that it is associated with low friction coefficients (~0.3). Grain size measurements by STEM validate a recently proposed dynamic amorphization model that accurately predicts grain size-dependent shear strength in the inverse Hall-Petch regime. Wear rates in the initially soft (coarse grained) material were shown to be remarkably low (~10-6 mm3/N-m), the lowest for any HEA tested in an inert environment where oxidation and the formation of mixed metal-oxide films is mitigated. The combined high wear resistance and low friction are linked to the formation of an ultra-nanocrystalline near-surface layer. The dynamic amorphization model was also used to predict an average high angle grain boundary energy (0.87 J/m2). This value was used to explain cavitation-induced nanoporosity found in the highly deformed surface layer, a phenomenon that has been linked to superplasticity.

Ultimate Strength of Metals.

We present a theoretical model that predicts the peak strength of polycrystalline metals based on the activation energy (or stress) required to cause deformation via amorphization. Building on extensive earlier work, this model is based purely on materials properties, requires no adjustable parameters, and is shown to accurately predict the strength of four exemplar metals (fcc, bcc, and hcp, and an alloy). This framework reveals new routes for design of more complex high-strength materials systems, such as compositionally complex alloys, multiphase systems, nonmetals, and composite structures.

New nanoscale toughening mechanisms mitigate embrittlement in binary nanocrystalline alloys.

Nanocrystalline metals offer significant improvements in structural performance over conventional alloys. However, their performance is limited by grain boundary instability and limited ductility. Solute segregation has been proposed as a stabilization mechanism, however the solute atoms can embrittle grain boundaries and further degrade the toughness. In the present study, we confirm the embrittling effect of solute segregation in Pt-Au alloys. However, more importantly, we show that inhomogeneous chemical segregation to the grain boundary can lead to a new toughening mechanism termed compositional crack arrest. Energy dissipation is facilitated by the formation of nanocrack networks formed when cracks arrested at regions of the grain boundaries that were starved in the embrittling element. This mechanism, in concert with triple junction crack arrest, provides pathways to optimize both thermal stability and energy dissipation. A combination of in situ tensile deformation experiments and molecular dynamics simulations elucidate both the embrittling and toughening processes that can occur as a function of solute content.

Achieving Ultralow Wear with Stable Nanocrystalline Metals.

Recent work suggests that thermally stable nanocrystallinity in metals is achievable in several binary alloys by modifying grain boundary energies via solute segregation. The remarkable thermal stability of these alloys has been demonstrated in recent reports, with many alloys exhibiting negligible grain growth during prolonged exposure to near-melting temperatures. Pt-Au, a proposed stable alloy consisting of two noble metals, is shown to exhibit extraordinary resistance to wear. Ultralow wear rates, less than a monolayer of material removed per sliding pass, are measured for Pt-Au thin films at a maximum Hertz contact stress of up to 1.1 GPa. This is the first instance of an all-metallic material exhibiting a specific wear rate on the order of 10-9 mm3 N-1 m-1 , comparable to diamond-like carbon (DLC) and sapphire. Remarkably, the wear rate of sapphire and silicon nitride probes used in wear experiments are either higher or comparable to that of the Pt-Au alloy, despite the substantially higher hardness of the ceramic probe materials. High-resolution microscopy shows negligible surface microstructural evolution in the wear tracks after 100k sliding passes. Mitigation of fatigue-driven delamination enables a transition to wear by atomic attrition, a regime previously limited to highly wear-resistant materials such as DLC.

Impact of Microstructure on MoS2 Oxidation and Friction.

This work demonstrates the role of microstructure in the friction and oxidation behavior of the lamellar solid lubricant molybdenum disulfide (MoS2). We report on systematic investigations of oxidation and friction for two MoS2 films with distinctively different microstructures-amorphous and planar/highly-ordered-before and after exposure to atomic oxygen (AO) and high-temperature (250 °C) molecular oxygen. A combination of experimental tribology, molecular dynamics simulations, X-ray photoelectron spectroscopy (XPS), and high-sensitivity low-energy ion scattering (HS-LEIS) was used to reveal new insights about the links between structure and properties of these widely utilized low-friction materials. Initially, ordered MoS2 films showed a surprising resistance to both atomic and molecular oxygens (even at elevated temperature), retaining characteristic low friction after exposure to extreme oxidative environments. XPS shows comparable oxidation of both coatings via AO; however, monolayer resolved compositional depth profiles from HS-LEIS reveal that the microstructure of the ordered coatings limits oxidation to the first atomic layer.

Molecular scale modeling of polymer imprint nanolithography.

We present the results of large-scale molecular dynamics simulations of two different nanolithographic processes, step-flash imprint lithography (SFIL), and hot embossing. We insert rigid stamps into an entangled bead-spring polymer melt above the glass transition temperature. After equilibration, the polymer is then hardened in one of two ways, depending on the specific process to be modeled. For SFIL, we cross-link the polymer chains by introducing bonds between neighboring beads. To model hot embossing, we instead cool the melt to below the glass transition temperature. We then study the ability of these methods to retain features by removing the stamps, both with a zero-stress removal process in which stamp atoms are instantaneously deleted from the system as well as a more physical process in which the stamp is pulled from the hardened polymer at fixed velocity. We find that it is necessary to coat the stamp with an antifriction coating to achieve clean removal of the stamp. We further find that a high density of cross-links is necessary for good feature retention in the SFIL process. The hot embossing process results in good feature retention at all length scales studied as long as coated, low surface energy stamps are used.

Forces between functionalized silica nanoparticles in solution.

To prevent the flocculation and phase separation of nanoparticles in solution, nanoparticles are often functionalized with short chain surfactants. Here we present fully atomistic molecular dynamics simulations which characterize how these functional coatings affect the interactions between nanoparticles and with the surrounding solvent. For 5-nm-diameter silica nanoparticles coated with poly(ethylene oxide) (PEO) oligomers in water, we determined the hydrodynamic drag on two approaching nanoparticles moving through solvent and on a single nanoparticle as it approaches a planar surface. In most circumstances, macroscale fluid theory accurately predicts the drag on these nanoscale particles. Good agreement is seen with Brenner's analytical solutions for wall separations larger than the soft nanoparticle radius. For two approaching coated nanoparticles, the solvent-mediated (velocity independent) and lubrication (velocity-dependent) forces are purely repulsive and do not exhibit force oscillations that are typical of uncoated rigid spheres.

Molecular dynamics simulations of water confined between matched pairs of hydrophobic and hydrophilic self-assembled monolayers.

We have conducted a molecular dynamics (MD) simulation study of water confined between methyl-terminated and carboxyl-terminated alkylsilane self-assembled monolayers (SAMs) on amorphous silica substrates. In doing so, we have investigated the dynamic and structural behavior of the water molecules when compressed to loads ranging from 20 to 950 MPa for two different amounts of water (27 and 58 water molecules/nm2). Within the studied range of loads, we observe that no water molecules penetrate the hydrophobic region of the carboxyl-terminated SAMs. However, we observe that at loads larger than 150 MPa water molecules penetrate the methyl-terminated SAMs and form hydrogen-bonded chains that connect to the bulk water. The diffusion coefficient of the water molecules decreases as the water film becomes thinner and pressure increases. When compared to bulk diffusion coefficients of water molecules at the various loads, we found that the diffusion coefficients for the systems with 27 water molecules/nm2 are reduced by a factor of 20 at low loads and by a factor of 40 at high loads, while the diffusion coefficients for the systems with 58 water molecules/nm2 are reduced by a factor of 25 at all loads.

Water in nanoconfinement between hydrophilic self-assembled monolayers.

Molecular dynamics (MD) simulations of water confined to subnanometer thicknesses between carboxyl-terminated alkanethiol self-assembled monolayers (SAMs) on gold were performed to address conflicts in the literature on the structure and response of water in confinement. The amount of water was varied to yield submonolayer to bilayer structures. The orientation of the water is affected by the confinement, especially in the submonolayer case. We find that the diffusion coefficient decreases as the film becomes thinner and at higher pressures. However, in all cases studied, liquid diffusion is always found. At maximal suppression, the diffusion constant is 2 orders of magnitude smaller than the bulk value.

Water penetration of damaged self-assembled monolayers.

The interaction of water with self-assembled monolayers (SAMs) on amorphous silica is investigated using classical molecular dynamics simulation. Damage is induced through shear simulations with model atomic force microscopy (AFM) tips and separately with controlled extraction. We find that SAM coatings that have been slightly damaged (by normal loads close to 10 nN from a 10-nm-diameter AFM tip) are susceptible to water penetration and migration to the underlying hydrophilic substrate. The controlled damage studies indicate that the presence of water tends to heal damage below a threshold radius and exploits and magnifies damage above this threshold. For the systems studied here, Si(OH)3(CH2)10CH3 alkylsilane chains on amorphous silica, this threshold radius is between 0.5 and 1.0 nm.

Simulations of nanotribology with realistic probe tip models.

We present the results of massively parallel molecular dynamics simulations aimed at understanding the nanotribological properties of alkylsilane self-assembled monolayers (SAMs) on amorphous silica. In contrast to studies with opposing flat plates, as found in the bulk of the simulation literature, we use a model system with a realistic AFM tip (radius of curvature ranging from 3 to 30 nm) in contact with a SAM-coated silica substrate. We compare the differences in response between systems in which chains are fully physisorbed, fully chemisorbed, and systems with a mixture of the two. Our results demonstrate that the ubiquitous JKR and DMT models do not accurately describe the contact mechanics of these systems. In shear simulations, we find that the chain length has minimal effects on both the friction force and coefficient. The tip radius affects the friction force only (i.e., the coefficient is unchanged) by a constant shift in magnitude due to the increase in pull-off force with increasing radius. We also find that at extremely low loads, on the order of 10 nN, shearing from the tip causes damage to the physisorbed monolayers by removal of molecules.

Density functional theory and DFT+U study of transition metal porphines adsorbed on Au(111) surfaces and effects of applied electric fields.

We apply density functional theory (DFT) and the DFT+U technique to study the adsorption of transition metal porphine molecules on atomistically flat Au(111) surfaces. DFT calculations using the Perdew-Burke-Ernzerhof exchange correlation functional correctly predict the palladium porphine (PdP) low-spin ground state. PdP is found to adsorb preferentially on gold in a flat geometry, not in an edgewise geometry, in qualitative agreement with experiments on substituted porphyrins. It exhibits no covalent bonding to Au(111), and the binding energy is a small fraction of an electronvolt. The DFT+U technique, parametrized to B3LYP-predicted spin state ordering of the Mn d-electrons, is found to be crucial for reproducing the correct magnetic moment and geometry of the isolated manganese porphine (MnP) molecule. Adsorption of Mn(II)P on Au(111) substantially alters the Mn ion spin state. Its interaction with the gold substrate is stronger and more site-specific than that of PdP. The binding can be partially reversed by applying an electric potential, which leads to significant changes in the electronic and magnetic properties of adsorbed MnP and approximately 0.1 A changes in the Mn-nitrogen distances within the porphine macrocycle. We conjecture that this DFT+U approach may be a useful general method for modeling first-row transition metal ion complexes in a condensed-matter setting.

Tribological properties of alkylsilane self-assembled monolayers.

In this study, we perform molecular dynamics simulations of adhesive contact and friction between alkylsilane Si(OH)(3)(CX(2))(10)CX(3) and alkoxylsilane Si(OH)(2)(CX(2))(10)CX(3) (where X = H or F) self-assembled monolayers (SAMs) on an amorphous silica substrate. The alkylsilane SAMs are primarily hydrogen-bonded or physisorbed to the surface. The alkoxylsilane SAMs are covalently bonded or chemisorbed to the surface. Previously, we studied the chemisorbed systems. In this work, we study the physisorbed systems and compare the tribological properties with the chemisorbed systems. Furthermore, we examine how water at the interface of the SAMs and substrate affects the tribological properties of the physisorbed systems. When less than a third of a monolayer is present, very little difference in the microscopic friction coefficient mu or shear stresses is observed. For increasing amounts of water, the values of mu and the shear stresses decrease; this effect is somewhat more pronounced for fluorocarbon alkylsilane SAMs than for the hydrocarbon SAMs. The observed decrease in friction is a consequence of a slip plane that occurs in the water as the amount of water is increased. We studied the frictional behavior using relative shear velocities ranging from v = 2 cm/s to 2 m/s. Similar to previously reported results for alkoxylsilane SAMs, the values of the measured stress and mu for the alkylsilane SAM systems decrease monotonically with v.

Frictional dynamics of fluorine-terminated alkanethiol self-assembled monolayers.

The frictional dynamics of fluorine-terminated alkanethiol (S(CH2)8CF3) self-assembled monolayers (SAMs) on gold are studied using molecular dynamics simulations. The simulations treat the interactions between two SAMs on flat surfaces. The structure and frictional behavior are investigated as a function of applied pressure (200 MPa to 1 GPa) for a shear velocity of 2 m/s and compared to methyl-terminated alkanethiol SAMs. The maximum adhesive pressure between the SAMs is 220 MPa for both end groups. In agreement with experiments on the molecular scale, the shear stress and the coefficient of friction for CF3-terminated alkanethiols are larger than for CH3-terminated alkanethiols. The main source for the difference is primarily the tighter packing of the fluorinated terminal group resulting in a higher degree of order. The molecular scale coefficient of friction is correlated with the degree of order among all the systems.