Effects of surface chemistry on the mechanochemical decomposition of tricresyl phosphate

The growth of protective tribofilms from lubricant antiwear additives on rubbing surfaces is initiated by mechanochemically promoted dissociation reactions. These processes are not well understood at the molecular scale for many important additives, such as tricresyl phosphate (TCP). One aspect that needs further clarification is the extent to which the surface properties affect the mechanochemical decomposition. Here, we use nonequilibrium molecular dynamics (NEMD) simulations with a reactive force field (ReaxFF) to study the decomposition of TCP molecules confined and pressurised between sliding ferrous surfaces at a range of temperatures. We compare the decomposition of TCP on native iron, iron carbide, and iron oxide surfaces. We show that the decomposition rate of TCP molecules on all the surfaces increases exponentially with temperature and shear stress, implying that this is a stress-augmented thermally activated (SATA) process. The presence of base oil molecules in the NEMD simulations decreases the shear stress, which in turn reduces the rate constant for TCP decomposition. The decomposition is much faster on iron surfaces than iron carbide, and particularly iron oxide. The activation energy, activation volume, and pre-exponential factor from the Bell model are similar on iron and iron carbide surfaces, but significantly differ for iron oxide surfaces. These findings provide new insights into the mechanochemical decomposition of TCP and have important implications for the design of novel lubricant additives for use in high-temperature and high-pressure environments.

Effect of Fe-O ReaxFF on Liquid Iron Oxide Properties Derived from Reactive Molecular Dynamics

As iron powder nowadays attracts research attention as a carbon-free, circular energy carrier, molecular dynamics (MD) simulations can be used to better understand the mechanisms of liquid iron oxidation at elevated temperatures. However, prudence must be practiced in the selection of a reactive force field. This work investigates the influence of currently available reactive force fields (ReaxFFs) on a number of properties of the liquid iron-oxygen (Fe-O) system derived (or resulting) from MD simulations. Liquid Fe-O systems are considered over a range of oxidation degrees ZO, which represents the molar ratio of O/(O + Fe), with 0 < ZO < 0.6 and at a constant temperature of 2000 K, which is representative of the combustion temperature of micrometric iron particles burning in air. The investigated properties include the minimum energy path, system structure, (im)miscibility, transport properties, and the mass and thermal accommodation coefficients. The properties are compared to experimental values and thermodynamic calculation results if available. Results show that there are significant differences in the properties obtained with MD using the various ReaxFF parameter sets. Based on the available experimental data and equilibrium calculation results, an improved ReaxFF is required to better capture the properties of a liquid Fe-O system.

Fe3 C/nanocarbon-Enabled Lithium Dendrite Mitigation in Lithium-Sulfur batteries

Lithium dendrite-induced short circuits and material loss are two major obstacles to the commercialization of lithium-sulfur (Li-S) batteries. Here, a nanocarbon composite consisting of cotton-derived Fe3 C-encapsulated multiwalled carbon nanotubes (Fe3 C-MWCNTs) and graphene effectively traps polysulfides to suppress lithium dendrite growth is reported. Machine learning combined with molecular dynamics (MD) simulations unveils a new polysulfide-induced lithium dendrite formation mechanism: the migration of polysulfides away from the anode drags out lithium protrusions through localized lattice distortion of the lithium anode and traps lithium ions in the surrounding electrolyte, leading to lithium dendrite formation. The Li-S battery, constructed using the composite of cotton-derived Fe3 C-MWCNTs and graphene that serves as both the sulfur host and the anode interlayer, exhibits exceptional cycling stability, impressive capacity retention, and effective mitigation of lithium dendrite formation. The findings offer valuable strategies to prevent lithium dendrite formation and enhance understanding of lithium dendrite growth in Li-S batteries.

Evaluation of autophagy inhibition to combat cancer: (vanadium complex)-protein interactions, parameterization, and validation of a new force field

Autophagy has drawn attention from the scientific community, mainly because of its significant advantages over chemotherapeutic processes. One of these advantages is its direct action on cancer cells, avoiding possible side effects, unlike chemotherapy, which reaches tumor cells and affects healthy cells in the body, leading to a great loss in the quality of life of patients. In this way, it is known that vanadium complex (VC) [VO(oda)(phen)] has proven inhibition effect on autophagy process in pancreatic cancer cells. Keeping that in mind, molecular dynamics (MD) simulations can be considered excellent strategies to investigate the interaction of metal complexes and their biological targets. However, simulations of this type are strongly dependent on the appropriate choice of force field (FF). Therefore, this work proposes the development of AMBER FF parameters for VC, having a minimum energy structure as a starting point, obtained through DFT calculations with B3LYP/def2-TZVP level of theory plus ECP for the vanadium atom. An MD simulation in vacuum was performed to validate the developed FF. From the structural analyses, satisfying values of VC bond lengths and angles were obtained, where a good agreement with the experimental data and the quantum reference was found. The RMSD analysis showed an average of only 0.3%. Finally, we performed docking and MD (120 ns) simulations with explicit solvent between VC and PI3K. Overall, our findings encourage new parameterizations of metal complexes with significant biological applications, as well as allow to contribute to the elucidation of the complex process of autophagy.

Hydrophobic Nanoconfinement Enhances CO2 Conversion to H2CO3

Understanding the formation of H₂CO₃ in water from CO₂ is important in environmental and industrial processes. Although numerous investigations have studied this reaction, the conversion of CO₂ to H₂CO₃ in nanopores, and how it differs from that in bulk water, has not been understood. We use ReaxFF metadynamics molecular simulations to demonstrate striking differences in the free energy of CO₂ conversion to H₂CO₃ in bulk and nanoconfined aqueous environments. We find that nanoconfinement not only reduces the energy barrier but also reverses the reaction from endothermic in bulk water to exothermic in nanoconfined water. Also, charged intermediates are observed more often under nanoconfinement than in bulk water. Stronger solvation and more favorable proton transfer with increasing nanoconfinement enhance the thermodynamics and kinetics of the reaction. Our results provide a detailed mechanistic understanding of an important step in the carbonation process, which depends intricately on confinement, surface chemistry, and CO₂ concentration.

Effects of Iron Addition on the Collision of Polycyclic Aromatic Hydrocarbon Clusters: A Molecular Dynamics Study

In this work, soot particle size distributions in iron-doped premixed ethylene flames are examined using scanning mobility particle sizer measurements. It is found that iron addition promotes the growth in soot particle size, and the enhanced particle coagulation is inferred to be an important reason. To support that, the influence of iron addition on the coagulation of polycyclic aromatic hydrocarbon (PAH) clusters, the analogue of incipient soot particles, is further investigated using molecular dynamics simulations. Based on the results of hundreds of binary head-on collision simulations, the collision between two coronene-Fe-coronene dimers is found to have a significantly higher coagulation efficiency than that between two coronene dimers. However, this enhancement effect weakens as the size of the PAH monomer increases. Although the coagulation efficiency can be increased by iron addition, the collision frequency is almost unaffected, as revealed from the binary off-central collision simulations. Moreover, the simulation results of coronene cluster growth via coagulation show that iron addition promotes coronene cluster growth, leading to larger cluster size, which may explain the larger soot particle size observed in iron-doped flames.

Initial stage of carbonization of iron during hydrocarbons dissociation: a molecular dynamics study

The carbonization of iron is a very important early phenomenon in the field of heterogeneous catalysis and the petrochemical industry, but the mechanism is still controversial. In this work, the carbonization mechanism and carbonization structure of iron nanoparticles by different carbon sources (CH₄, C₂H₆, C₂H₄, C₂H₂) were systematically investigated using the reactive molecular dynamics method. The results show that saturated alkanes are dehydrogenated while adsorbed, but unsaturated olefins and alkynes undergo bond-breaking while adsorbed. The C-H bond is more likely to break than the C-C bond. Hydrocarbons with high carbon content have a strong ability to carbonize Fe nanoparticles under the same conditions. For C₂H₄ and C₂H₂, the C atoms generated from dissociation form a large number of long carbon chains intertwined with branched chains and multiple carbon rings. The C2 species formed by C₂H₂ after complete dehydrogenation diffuse rapidly to the interior of the nanoparticles, releasing the surface active sites and accelerating the carbonization process. Carbon-rich iron carbides (FeC_x) with different Fe/C ratios were obtained by carbonization with different carbon sources. In addition, the Fe(110) surface exhibits the strongest carburizing ability. These findings provide systematic insights into the initial stages of metal Fe carburization.

How Polytetrafluoroethylene Lubricates Iron: An Atomistic View by Reactive Molecular Dynamics

The tribochemistry and transfer film formation at the metal/polymer interface plays an essential role in surface protection, wear reduction, and lubrication. Although the topic has been studied for decades, challenges persist in clarifying the nanoscale mechanism and dynamic evolution of tribochemical reactions. To investigate the tribochemistry between iron and polytetrafluoroethylene (PTFE) in ambient and cryogenic environments, we have trained and expanded a ReaxFF reactive force field to describe iron-oxygen-water-PTFE systems (C/H/O/F/Fe). Using ReaxFF molecular dynamics simulations, we find that mechanical shearing of single asperity induced the degradation of PTFE molecules and radicals, showing subsequent oxidation and hydroxylation reactions of the radicals initiated by C-C bond cleavage, in agreement with previous experimental observations. Furthermore, we studied mechanisms of interfacial tribochemical reactions and formation of transfer films. We found that tribochemical wear and Fe-C and Fe-F bonding networks are important mechanisms for anchoring molecular chains to form a transfer film on the iron countersurface. Hydroxyl groups can dehydrogenate to form short and strong chelation bonds with the Fe₂O₃ countersurface. A friction-induced oriented molecular layer plays a key role in reducing friction, which is responsible for the excellent lubrication property. By varying temperatures in the range of 10-300 K, we found a nonmonotonic change in friction with a maxima at 100 K. At cryogenic temperatures, the molecular mobility was obviously suppressed, while the chain rigidity was enhanced, resulting in the less oriented interface and brittle-like shear interface, which is responsible for nonmonotonic friction. This work elucidates mechanisms of tribochemical reactions and transfer film formation between iron and PTFE at the atomistic level, facilitating design and development of self-lubricating materials, especially under harsh conditions.

A Reactive Molecular Dynamics Study of Bi-modal Particle Size Distribution in Binder-Jetting Additive Manufacturing using Stainless-Steel Powders

In order to improve a final part’s density and achieve desired mechanical properties, binder-jetting additive manufacturing usually requires lengthy post-processing steps such as curing, sintering, and infiltration. The role of...

Modelling of high-temperature order-disorder phase transitions of non-stoichiometric Mo2C and Ti2C from first principles

High-temperature order-disorder phase transitions play an important role in determining the structure and physical and chemical properties of non-stoichiometric transition metal carbides. Due to the large number of possible carbon vacancy arrangements, it is difficult to study these systems with first-principles calculations. Here, we construct a simple atomistic potential capable of accurately reproducing the energetics of the carbon vacancy arrangements in cubic Mo₂C and Ti₂C obtained from density functional theory calculations. We show that this potential can be applied to correctly predict the transition temperatures between the ordered and disordered states in Monte Carlo simulations on large supercells and reveal the extend of local order in the disordered phases of Mo₂C and Ti₂C that show interesting physical and chemical properties. We find that even the high-temperature disordered phase exhibit a relatively high degree of local order as indicated by the relatively small change in the root mean square number of C atom neighbours of Mo/Ti compared to the ordered phase (from 3.0 to 3.1-3.2). This atomistic potential enables the study of how the structure of these carbides can be tuned through the synthesis temperature to control the properties of carbide materials that are related to the degree of disorder in the system such as catalytic activity and electrical conductivity and play an important role in applications of these carbides. Fundamentally, the successful modelling of these carbides suggests that despite the presence of metallic, covalent and ionic interactions, bonding in carbides can be modelled by simple and physically intuitive interatomic potentials.

Atomic-scale investigation of the effect of surface carbon coatings on the oxidation and mechanical properties of iron nanowires

The understanding of the complex atomistic-scale mechanisms of the oxidation process of carbon (C) coated iron nanowires (Fe NW) and also the resulting modulation of mechanical properties is a highly challenging task.

Interaction and thermal stability of carboxymethyl cellulose on α-Fe2O3(001) surface: ReaxFF molecular dynamics simulations study

Reactive (ReaxFF) molecular dynamic simulations were performed to elucidate the nature of interaction between hematite and carboxymethyl cellulose (CMC) considering effect of moisture and temperature. Simulations showed that the presence of moisture prohibited CMC to interact directly with hematite surface. However, the moisture helps to disperse CMC along the hematite surface thus maximizing the interaction. In dry condition, it was found that the negatively charged oxygen present in the functional groups and in CMC backbone (-CH2OCH2-) take part to form chemical bond with the positively charged surface Fe. The bonding interaction with polymer back bone is a new finding from this simulation. The binding energy estimated as -56.2 kcal/mol clearly indicates chemisorption with bond length ∼2.00 Å. Upon heating up the complex in presence of atmospheric oxygen it was observed that CMC decomposed at high temperature rather than desorbed from the surface. This study clearly shows that CMC is a suitable binder for iron ore pelletization and will pave the way for more rationale design of organic binders.

Reactive Molecular Dynamics Simulations of Thermal Film Growth from Di- tert-butyl Disulfide on an Fe(100) surface

Iron sulfide films are present in many applications, including lubricated interfaces where protective films are formed through the reactions of lubricant additive molecules with steel surfaces during operation. Such films are critical to the efficiency and useful lifetime of moving components. However, the mechanisms by which films form are still poorly understood because the reactions occur between two surfaces and so cannot be directly probed experimentally. To address this, we explore the thermal contribution to film formation of di- tert-butyl disulfide-an important extreme pressure additive-on an Fe(100) surface using reactive molecular dynamics simulations, where the reactive potential parameters are validated by comparison to ab initio calculations. The reaction pathway leading to the formation of iron sulfide surfaces is characterized using the reactive simulations. Then, the film formation process is mimicked by simulations where di- tert-butyl disulfide molecules are cyclically added to the surface and subjected to temperatures comparable to those expected due to frictional heating. The use of a reactive empirical potential is a novel approach to modeling the iterative nature of thermal film growth with realistic lubricant additive molecules.

Simulating the effects of carbon nanotube continuity and interfacial bonding on composite strength and stiffness

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.

Development of a reactive force field for the Fe-C interaction to investigate the carburization of iron

The approach of molecular dynamics with Reactive Force Field (ReaxFF) is a promising way to investigate the carburization of iron which is pivotal in the preparation of desired iron-based materials and catalysts. However, it is a challenge to develop a reliable ReaxFF to describe the Fe-C interaction, especially when it involves bond rearrangement. In this work, we develop an exclusive set of Reactive Force Field (ReaxFF) parameters, denoted RPOIC-2017, to describe the diffusion behavior of carbon atoms in the α-Fe system. It inherited some partial parameters in 2012 (ReaxFF-2012) which are suitable for hydrogen adsorption and dissociation. This set of parameters is trained against data from first-principles calculations, including the equations of state of α-Fe, the crystal constant of Fe₃C and Fe₄C, a variety of periodic surface structures with varying carbon coverages, as well as the barriers of carbon diffusion in the α-Fe bulk and on diverse surfaces. The success in predicting the carbon diffusion coefficient and the diffusion barrier using the developed RPOIC-2017 potential demonstrates that the performance is superior to that of the traditional MEAM potential. The new ReaxFF for the Fe-C interaction developed in this work is not only essential for the design of novel iron based materials, but could also help understand atomic arrangements and the interfacial structure of iron carbides.

ReaxFF based molecular dynamics simulations of ignition front propagation in hydrocarbon/oxygen mixtures under high temperature and pressure conditions

This work investigates the ignition front speed of hydrocarbon fuels at atomistic level for the first time using the ReaxFF reactive force field method.