RESUMEN
Boron doping of diamond-like carbon coatings has multiple effects on their tribological properties. While boron typically reduces wear in cutting applications, some B-doped coatings show poor tribological performance compared with undoped films. This is the case of the tribological tests presented in this work in which an alumina ball is placed in frictional contact with different undoped and B-doped amorphous carbon coatings in humid air. With B-doped coatings, a higher friction coefficient at a steady state with respect to their undoped counterparts was observed. Estimates of the average contact shear stress based on experimental friction coefficients, surface topographies, and Persson's contact theory suggest that the increased friction is compatible with the formation of a sparse network of interfacial ether bonds leading to a mild cold-welding friction regime, as documented in the literature. Tight binding and density functional theory simulations were performed to investigate the chemical effect of B-doping on the interfacial properties of the carbon coatings. The results reveal that OH groups that normally passivate carbon surfaces in humid environments can be activated by boron and form B-O dative bonds across the tribological interfaces, leading to a mild cold-welding friction regime. Simulations performed on different tribological pairs suggest that this mechanism could be valid for B-doped carbon surfaces in contact with a variety of materials. In general, this study highlights the impact that subtle modifications in surface and interface chemistry caused by the presence of impurities can have on macroscopic properties, such as friction and wear.
RESUMEN
Tribological interfaces between silicon-based ceramics, such as Si3N4 or SiC, are characterized by high friction and wear in unlubricated conditions. A solution to this problem is to use them in combination with a hydrogenated amorphous carbon (a-C:H) countersurface from which a passivating carbon film is transferred onto the ceramic surface. However, the mechanisms underlying a stable film transfer process and the conditions that favor it remain elusive. Here, we present friction experiments in ultrahigh vacuum in which friction coefficients lower than 0.01 are achieved by sliding Si3N4 against a-C:H with 36 at. % hydrogen but not against a-C:H with 20 at. % hydrogen. Chemical surface analyses confirm that the superlubric interface forms via the transfer of a hydrocarbon nanofilm onto the Si3N4 surface. Quantum-mechanical simulations reveal that a stable passivating a-C:H film can only be transferred if, after initial cold welding of the tribological interface, the plastic shear deformation is localized within the a-C:H coating. This occurs if the yield shear stress for plastic flow of a-C:H is lower than that of the ceramic and of the shear strength of the a-C:H-ceramic interface, i.e., if the a-C:H hydrogen content ranges between â¼30 and â¼50 at. %. While the importance of a relatively high hydrogen content to achieve an efficient passivation of a-C:H surfaces in a vacuum is well-documented, this work reveals how the hydrogen content is also crucial for obtaining a stable a-C:H transfer film. These results can be extended to glass, SiC, and steel, supporting the generality of the proposed mechanism.
RESUMEN
The Reynolds lubrication equation (RLE) is widely used to design sliding contacts in mechanical machinery. While providing an excellent description of hydrodynamic lubrication, friction in boundary lubrication regions is usually considered by empirical laws, because continuum theories are expected to fail for lubricant film heights h0 ⪠10 nm, especially in highly loaded tribosystems with normal pressures pn â« 0.1 GPa. Here, the performance of RLEs is validated by molecular dynamics simulations of pressurized (with pn = 0.2 to 1 GPa) hexadecane in a gold converging-diverging channel with minimum gap heights h0 = 1.4 to 9.7 nm. For pn ≤ 0.4 GPa and h0 ≥ 5 nm, agreement with the RLE requires accurate constitutive laws for pressure-dependent density and viscosity. An additional nonlinear wall slip law relating wall slip velocities to local shear stresses extends the RLE's validity to even the most severe loading condition pn = 1 GPa and h0 = 1.4 nm. Our results demonstrate an innovative route for continuum modeling of highly loaded tribological contacts under boundary lubrication.
RESUMEN
Friction in boundary lubrication is strongly influenced by the atomic structure of the sliding surfaces. In this work, friction between dry amorphous carbon (a-C) surfaces with chemisorbed fragments of lubricant molecules is investigated employing molecular dynamic simulations. The influence of length, grafting density and polarity of the fragments on the shear stress is studied for linear alkanes and alcohols. We find that the shear stress of chain-passivated a-C surfaces is independent of the a-C density. Among all considered chain-passivated systems, those with a high density of chains of equal length exhibit the lowest shear stress. However, shear stress in chain-passivated a-C is consistently higher than in a-C surfaces with atomic passivation. Finally, surface passivation species with OH head groups generally lead to higher friction than their non-polar analogs. Beyond these qualitative trends, the shear stress behavior for all atomic- and chain-passivated, non-polar systems can be explained semi-quantitatively by steric interactions between the two surfaces that cause resistance to the sliding motion. For polar passivation species electrostatic interactions play an additional role. A corresponding descriptor that properly captures the interlocking of the two surfaces along the sliding direction is developed based on the maximum overlap between atoms of the two contacting surfaces.
RESUMEN
We study mechanically induced phase transitions at tribological interfaces between silicon crystals using reactive molecular dynamics. The simulations reveal that the interplay between shear-driven amorphization and recrystallization results in an amorphous shear interface with constant thickness. Different shear elastic responses of the two anisotropic crystals can lead to the migration of the amorphous interface normal to the sliding plane, causing the crystal with lowest elastic energy density to grow at the expense of the other one. This triboepitaxial growth can be achieved by crystal misorientation or exploiting elastic finite-size effects, enabling the direct deposition of homoepitaxial silicon nanofilms by a crystalline tip rubbing against a substrate.
RESUMEN
Friction and wear reduction by diamond-like carbon (DLC) in automotive applications can be affected by zinc-dialkyldithiophosphate (ZDDP), which is widely used in engine oils. Our experiments show that DLC's tribological behaviour in ZDDP-additivated oils can be optimised by tailoring its stiffness, surface nano-topography and hydrogen content. An optimal combination of ultralow friction and negligible wear is achieved using hydrogen-free tetrahedral amorphous carbon (ta-C) with moderate hardness. Softer coatings exhibit similarly low wear and thin ZDDP-derived patchy tribofilms but higher friction. Conversely, harder ta-Cs undergo severe wear and sub-surface sulphur contamination. Contact-mechanics and quantum-chemical simulations reveal that shear combined with the high local contact pressure caused by the contact stiffness and average surface slope of hard ta-Cs favour ZDDP fragmentation and sulphur release. In absence of hydrogen, this is followed by local surface cold welding and sub-surface mechanical mixing of sulphur resulting in a decrease of yield stress and wear.
RESUMEN
A stable passivation of surface dangling bonds underlies the outstanding friction properties of diamond and diamond-like carbon (DLC) coatings in boundary lubrication. While hydrogen is the simplest termination of a carbon dangling bond, fluorine can also be used as a monoatomic termination, providing an even higher chemical stability. However, whether and under which conditions a substitution of hydrogen with fluorine can be beneficial to friction is still an open question. Moreover, which of the chemical differences between C-H and C-F bonds are responsible for the change in friction has not been unequivocally understood yet. In order to shed light on this problem, we develop a density functional theory-based, nonreactive force field that describes the relevant properties of hydrogen- and fluorine-terminated diamond and DLC tribological interfaces. Molecular dynamics and nudged elastic band simulations reveal that the frictional stress at such interfaces correlates with the corrugation of the contact potential energy, thus ruling out a significant role of the mass of the terminating species on friction. Furthermore, the corrugation of the contact potential energy is almost exclusively determined by steric factors, while electrostatic interactions only play a minor role. In particular, friction between atomically flat diamond surfaces is controlled by the density of terminations, by the C-H and C-F bond lengths, and by the H and F atomic radii. For sliding DLC/DLC interfaces, the intrinsic atomic-scale surface roughness plays an additional role. While surface fluorination decreases the friction of incommensurate diamond contacts, it can negatively affect the friction performance of carbon surfaces that are disordered and not atomically flat. This work provides a general framework to understand the impact of chemical structure of surfaces on friction and to generate design rules for optimally terminated low-friction systems.
RESUMEN
Superlubricity of tetrahedral amorphous carbon (ta-C) coatings under boundary lubrication with organic friction modifiers is important for industrial applications, but the underlying mechanisms remain elusive. Here, combined experiments and simulations unveil a universal tribochemical mechanism leading to superlubricity of ta-C/ta-C tribopairs. Pin-on-disc sliding experiments show that ultra- and superlow friction with negligible wear can be achieved by lubrication with unsaturated fatty acids or glycerol, but not with saturated fatty acids and hydrocarbons. Atomistic simulations reveal that, due to the simultaneous presence of two reactive centers (carboxylic group and C=C double bond), unsaturated fatty acids can concurrently chemisorb on both ta-C surfaces and bridge the tribogap. Sliding-induced mechanical strain triggers a cascade of molecular fragmentation reactions releasing passivating hydroxyl, keto, epoxy, hydrogen and olefinic groups. Similarly, glycerol's three hydroxyl groups react simultaneously with both ta-C surfaces, causing the molecule's complete mechano-chemical fragmentation and formation of aromatic passivation layers with superlow friction.
RESUMEN
Large-scale quantum molecular dynamics of water-lubricated diamond (111) surfaces in sliding contact reveals multiple friction regimes. While water starvation causes amorphization of the tribological interface, small H_{2}O traces are sufficient to preserve crystallinity. This can result in high friction due to cold welding via ether groups or in ultralow friction due to aromatic surface passivation triggered by tribo-induced Pandey reconstruction. At higher water coverage, Grotthuss-type diffusion and H_{2}O dissociation yield dense H/OH surface passivation leading to another ultralow friction regime.
RESUMEN
Despite the pronounced polarity of C-F bonds, many fluorinated carbon compounds are hydrophobic: a controversial phenomenon known as "polar hydrophobicity". Here, its underlying microscopic mechanisms are explored by ab initio calculations of fluorinated and hydrogenated diamond (111) surfaces interacting with single water molecules. Gradient- and van der Waals-corrected density functional theory simulations reveal that "polar hydrophobicity" of the fully fluorinated surfaces is caused by a negligible surface/water electrostatic interaction. The densely packed C-F surface dipoles generate a short-range electric field that decays within the core repulsion zone of the surface and hence vanishes in regions accessible by adsorbates. As a result, water physisorption on fully F-terminated surfaces is weak (adsorption energies Ead < 0.1 eV) and dominated by van der Waals interactions. Conversely, the near-surface electric field generated by loosely packed dipoles on mixed F/H-terminated surfaces has a considerably longer range, resulting in a stronger water physisorption (Ead > 0.2 eV) that is dominated by electrostatic interactions. The suppression of electrostatic interactions also holds for perfluorinated molecular carbon compounds, thus explaining the prevalent hydrophobicity of fluorocarbons. In general, densely packed polar terminations do not always lead to short-range electric fields. For example, surfaces with substantial electron density spill-out give rise to electric fields with a much slower decay. However, electronic spill-out is limited in F/H-terminated carbon materials. Therefore, our ab initio results can be reproduced and rationalized by a simple classical point-charge model. Consequently, classical force fields can be used to study the wetting of F/H-terminated diamond, revealing a pronounced correlation between adsorption energies of single H2O molecules and water contact angles.
RESUMEN
We report comparisons between energy-based quantum mechanics/molecular mechanics (QM/MM) and buffered force-based QM/MM simulations in silica. Local quantities-such as density of states, charges, forces, and geometries-calculated with both QM/MM approaches are compared to the results of full QM simulations. We find the length scale over which forces computed using a finite QM region converge to reference values obtained in full quantum-mechanical calculations is â¼10 Å rather than the â¼5 Å previously reported for covalent materials such as silicon. Electrostatic embedding of the QM region in the surrounding classical point charges gives only a minor contribution to the force convergence. While the energy-based approach provides accurate results in geometry optimizations of point defects, we find that the removal of large force errors at the QM/MM boundary provided by the buffered force-based scheme is necessary for accurate constrained geometry optimizations where Si-O bonds are elongated and for finite-temperature molecular dynamics simulations of crack propagation. Moreover, the buffered approach allows for more flexibility, since special-purpose QM/MM coupling terms that link QM and MM atoms are not required and the region that is treated at the QM level can be adaptively redefined during the course of a dynamical simulation.
RESUMEN
We present a quantum-accurate multiscale study of how hydrogen-filled discoidal "platelet" defects grow inside a silicon crystal. Dynamical simulations of a 10-nm-diameter platelet reveal that H2 molecules form at its internal surfaces, diffuse, and dissociate at its perimeter, where they both induce and stabilize the breaking up of highly stressed silicon bonds. A buildup of H2 internal pressure is neither needed for nor allowed by this stress-corrosion growth mechanism, at odds with previous models. Slow platelet growth up to micrometric sizes is predicted as a consequence, making atomically smooth crystal cleavage possible in implantation experiments.