Prediction of a Novel Electromechanical Response in Polar Polymers with Rigid Backbones: Contrasting Furan-Derived Nanothreads to Poly(Vinylidene Fluoride).

Syn furan nanothreads have all oxygen atoms arranged on one side of the thread backbone; these polar threads present intriguing opportunities in electromechanical response owing to their rigid ladder-like backbone. We retrained a C/H/O reactive force field to simulate their response to external electric field for both end-anchored individual threads and bulk nanothread crystals, contrasting the results to those for poly(vinylidene fluoride) (PVDF) polymer. Whereas the field induces a length-independent torque in PVDF through backbone rotation about σ bonds, furan-derived nanothreads generate a length-dependent torque by progressively twisting their rigid backbone. This mode of response couples the rotational history of the electric field to axial tension in the anchored thread. In simulations of densely packed syn furan nanothread crystals without anchors, the crystals pole in a field (∼3 GV/m at 300 K) similar to that seen in simulations of PVDF, suggesting that crystals of polar nanothreads can be ferroelectric.

Watching (De)Intercalation of 2D Metals in Epitaxial Graphene: Insight into the Role of Defects.

Intercalation forms heterostructures, and over 25 elements and compounds are intercalated into graphene, but the mechanism for this process is not well understood. Here, the de-intercalation of 2D Ag and Ga metals sandwiched between bilayer graphene and SiC are followed using photoemission electron microscopy (PEEM) and atomistic-scale reactive molecular dynamics simulations. By PEEM, de-intercalation "windows" (or defects) are observed in both systems, but the processes follow distinctly different dynamics. Reversible de- and re-intercalation of Ag is observed through a circular defect where the intercalation velocity front is 0.5 nm s-1 ± 0.2 nm s.-1 In contrast, the de-intercalation of Ga is irreversible with faster kinetics that are influenced by the non-circular shape of the defect. Molecular dynamics simulations support these pronounced differences and complexities between the two Ag and Ga systems. In the de-intercalating Ga model, Ga atoms first pile up between graphene layers until ultimately moving to the graphene surface. The simulations, supported by density functional theory, indicate that the Ga atoms exhibit larger binding strength to graphene, which agrees with the faster and irreversible diffusion kinetics observed. Thus, both the thermophysical properties of the metal intercalant and its interaction with defective graphene play a key role in intercalation.

Advancing DBD Plasma Chemistry: Insights into Reactive Nitrogen Species such as NO2, N2O5, and N2O Optimization and Species Reactivity through Experiments and MD Simulations.

This study aims to fine-tune the plasma composition with a particular emphasis on reactive nitrogen species (RNS) including nitrogen dioxide (NO2), dinitrogen pentoxide (N2O5), and nitrous oxide (N2O), produced by a self-constructed cylindrical dielectric barrier discharge (CDBD). We demonstrated the effective manipulation of the plasma chemical profile by optimizing electrical properties, including the applied voltage and frequency, and by adjusting the nitrogen and oxygen ratios in the gas mixture. Additionally, quantification of these active species was achieved using Fourier transform infrared spectroscopy. The study further extends to exploring the aerosol polymerization of acrylamide (AM) into polyacrylamide (PAM), serving as a model reaction to evaluate the reactivity of different plasma-generated species, highlighting the significant role of NO2 in achieving high polymerization yields. Complementing our experimental data, molecular dynamics (MD) simulations, based on the ReaxFF reactive force field potential, explored the interactions between reactive oxygen species, specifically hydroxyl radicals (OH) and hydrogen peroxide (H2O2), with water molecules. Understanding these interactions, combined with the optimization of plasma chemistry, is crucial for enhancing the effectiveness of DBD plasma in environmental applications like air purification and water treatment.

Exploring Gas-Phase Chemical Reactions in NH3/B2H6 Systems for Chemical Vapor Deposition Using Reactive Molecular Dynamics.

Hexagonal boron nitride (hBN) has attracted significant attention as a two-dimensional (2D) material due to its unique structure and properties. In this paper, we investigated the gas-phase reactions between B2H6 and NH3 and the reaction pathways potentially leading to the synthesis of hBN by using ReaxFF-based reactive molecular dynamics (MD) simulations. From the hundreds of chemical reaction pathways observed in these MD simulations, we extracted a highly reduced chemical kinetic model to describe the gas-phase mixture evolution in chemical vapor deposition of hBN using B2H6 and NH3 as precursors. The intent is to integrate this chemical model into future computational fluid dynamics (CFD) simulations of actual hBN deposition testing and production reactors to provide enhanced insights for experimental synthesis processes and reactor optimization. The chemical model of this study will serve as a stepping stone for large-scale simulations.

Development of a ReaxFF Reactive Force Field for Pt/Cl Systems with Application to Platinum Metal Etching with Chlorine and Hydrogen Chloride Gases.

In this study, we present the development of a ReaxFF Pt/Cl/H reactive force field designed to elucidate the etching process by Cl for Pt surfaces. The ReaxFF force field parameters were optimized based on a quantum mechanical training set, which included adsorption energies of Cl and dissociation of HCl on Pt(100) and Pt(111) surfaces, energy/volume relations of PtCl2 crystals, and Cl diffusion on Pt(100) and Pt(111) surfaces. The predictive capability of the force field was further established through molecular dynamics simulations, which investigated the interactions of Cl2 and HCl molecules with the (100) and (111) surfaces of c-Pt crystalline solid slabs. A comparative analysis revealed that the Pt (100) surface exhibited higher susceptibility to chlorination and etching, leading to a more dominant removal of surface Pt atoms, whereas the Pt (111) surface showed greater resistance to these processes. This resistance impeded the access of Cl atoms to the Pt surface, resulting in a slower formation of PtxCly molecules. The etching ratios between HCl and Cl2 were compared with experimental results, yielding satisfactory agreement. This indicates that the developed ReaxFF protocol serves as a valuable tool for studying atomistic-scale details of the etching process in platinum metal systems.

Development and validation of a general-purpose ReaxFF reactive force field for earth material modeling.

ReaxFF reactive force field bridges the gap between nonreactive molecular simulations and quantum mechanical calculations and has been widely applied during the past two decades. However, its application to earth materials, especially those under high T-P conditions relevant to Earth's interior, is still limited due to the lack of available parameters. Here, we present the development and validation of a ReaxFF force field containing several of the most common elements in Earth's crust, i.e., Si/Al/O/H/Na/K. The force field was trained against a large data set obtained from density functional theory (DFT) calculations, including charges, bond/angle distortion curves, equation of states, ion migration energy profiles, and condensation reaction energies. Different coordination environments were considered in the training set. The fitting results showed that the current force field can well reproduce the DFT data (the Pearson correlation coefficient, Rp, is 0.95). We validated the force field on mineral-water interfaces, hydrous melts/supercritical geofluids, and bulk crystals. It was found that the current force field performed excellently in predicting the structural, thermodynamic, and transport properties of various systems (Rp = 0.95). Moreover, possible applications and future development have been discussed. The results obtained in this study suggest that the current force field holds good promise to model a wide range of processes and thus open opportunities to advance the application of ReaxFF in earth material modeling.

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.

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics.

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

Stable metal anodes enabled by a labile organic molecule bonded to a reduced graphene oxide aerogel.

Metallic anodes (lithium, sodium, and zinc) are attractive for rechargeable battery technologies but are plagued by an unfavorable metal-electrolyte interface that leads to nonuniform metal deposition and an unstable solid-electrolyte interphase (SEI). Here we report the use of electrochemically labile molecules to regulate the electrochemical interface and guide even lithium deposition and a stable SEI. The molecule, benzenesulfonyl fluoride, was bonded to the surface of a reduced graphene oxide aerogel. During metal deposition, this labile molecule not only generates a metal-coordinating benzenesulfonate anion that guides homogeneous metal deposition but also contributes lithium fluoride to the SEI to improve Li surface passivation. Consequently, high-efficiency lithium deposition with a low nucleation overpotential was achieved at a high current density of 6.0 mA cm-2 A Li|LiCoO2 cell had a capacity retention of 85.3% after 400 cycles, and the cell also tolerated low-temperature (-10 °C) operation without additional capacity fading. This strategy was applied to sodium and zinc anodes as well.

Molecular Alignment of a Meta-Aramid on Carbon Nanotubes by In Situ Interfacial Polymerization.

Molecularly organized nanocomposites of polymers and carbon nanotubes (CNTs) have great promise as high-performance materials; in particular, conformal deposition of polymers can control interfacial properties for mechanical load transfer, electrical or thermal transport, or electro/chemical transduction. However, controllability of polymer-CNT interaction remains a challenge with common processing methods that combine CNTs and polymers in melt or in solution, often leading to nonuniform polymer distribution and CNT aggregation. Here, we demonstrate CNTs within net-shape sheets can be controllably coated with a conformal coating of meta-aramid by simultaneous capillary infiltration and interfacial polymerization. We determine that π-interaction between the polymer and CNTs results in chain alignment parallel to the CNT outer wall. Subsequent nucleation and growth of the precipitated aramid forms a smooth continuous layered sheath around the CNTs. These findings motivate future investigation of mechanical properties of the resulting composites, and adaptation of the in situ polymerization method to other substrates.

Effects of interlayer confinement and hydration on capacitive charge storage in birnessite.

Nanostructured birnessite exhibits high specific capacitance and nearly ideal capacitive behaviour in aqueous electrolytes, rendering it an important electrode material for low-cost, high-power energy storage devices. The mechanism of electrochemical capacitance in birnessite has been described as both Faradaic (involving redox) and non-Faradaic (involving only electrostatic interactions). To clarify the capacitive mechanism, we characterized birnessite's response to applied potential using ex situ X-ray diffraction, electrochemical quartz crystal microbalance, in situ Raman spectroscopy and operando atomic force microscope dilatometry to provide a holistic understanding of its structural, gravimetric and mechanical responses. These observations are supported by atomic-scale simulations using density functional theory for the cation-intercalated structure of birnessite, ReaxFF reactive force field-based molecular dynamics and ReaxFF-based grand canonical Monte Carlo simulations on the dynamics at the birnessite-water-electrolyte interface. We show that capacitive charge storage in birnessite is governed by interlayer cation intercalation. We conclude that the intercalation appears capacitive due to the presence of nanoconfined interlayer structural water, which mediates the interaction between the intercalated cation and the birnessite host and leads to minimal structural changes.

Cellulose Nanocrystals: Tensile Strength and Failure Mechanisms Revealed Using Reactive Molecular Dynamics.

Cellulose nanocrystals (CNCs) offer excellent mechanical properties. However, measuring the strength by performing reliable experiments at the nanoscale is challenging. In this paper, we model Iß crystalline cellulose using reactive molecular dynamics simulations. Taking the fibril twist into account, structural changes and hydrogen-bonding characteristics of CNCs during the tensile test are inspected and the failure mechanism of CNCs is analyzed down to the scale of individual bonds. The C4-O4 glycosidic bond is found to be responsible for the failure of CNCs. Finally, the effect of strain rate on ultimate properties is analyzed and a nonlinear model is used to predict the ultimate strength of 9.2 GPa and ultimate strain of 8.5% at a 1 s-1 strain rate. This study sheds light on the applications of cellulose in nanocomposites and further modeling of cellulose nanofibres.

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions.

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na+, Ca2+, and Mg2+ cations and the CO32- anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me-O-C (Me = metal) three-body valence angle parameters and Me-C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na+ predominantly forms carbonates due to weaker angular strain, while Ca2+ and Mg2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. The nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na+ and Mg2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

A reactive molecular dynamics study of bi-modal particle size distribution in binder-jetting additive manufacturing using stainless-steel powders.

Binder-jetting is a powder-bed-based additive manufacturing (AM) process that is uniquely different from other powder-bed "fusion" metal AM technologies because it is a binder-based consolidation process similar to powder metallurgy "green" part and offers a larger selection of materials and part design capabilities. In order to improve a final part's density and achieve desired mechanical properties, binder-jetting usually requires lengthy post-processing steps such as curing, sintering, and infiltration. The role of particle size distribution in this process has been demonstrated to have a major impact. When comparing different distributions such as mono- and bi-modal sizes, the latter, consisting of a mix between coarser and finer particles, has shown to increase packing density and decrease porosity for a printed part. In this present work, we employ ReaxFF reactive force-field-based molecular dynamics (MD) simulations to study the atomistic level mechanism of binder-jetting using a bi-modal austenitic stainless-steel powder mixture. In addition, we compare the fracture process of the bi-modal powder mixture system with that of a system with mono-modal particle size, aiming to understand how the finer particles in the bi-modal powder mixture contribute to raising rupture strength. The MD simulation results show that the energy barriers after curing and sintering in the bi-modal particle system increase by 42.9% and 40.9%, respectively than in the mono-modal particle system. Moreover, the analysis of chemical composition and microstructure shows that iron is dominantly oxidized by oxygen atoms rather than hydroxyl radicals. Besides, the finer particle is subject to internal oxidation during sintering because its iron core melts. In contrast, the iron core of the coarser particle remains crystalline. Additionally, the statistical analysis of bonding oxygen atoms for each reference iron atom indicates that both particles have a small ratio of iron oxidized to Fe(II) but only slowly oxidizes to Fe(III) in the binder-jetting process. The coarser particle has a lot of non-oxidized iron atoms, while the majority of iron atoms in the finer particle bond with one oxygen atom during the time scale of our MD simulations. Furthermore, de-hydroxylation and oxygen inward diffusion lead to the reduction of chromium cations throughout sintering. The original findings of this study provide a nanoscale explanation for the mechanical property improvement using a bi-modal powder mixture. Moreover, the study of chemical composition and microstructure also contributes to improving the chemical properties of binder-jetting products.

Role of tilt grain boundaries on the structural integrity of WSe2 monolayers.

Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe2) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe2, of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe2 monolayer with tilted GBs. WSe2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures.

Condensation and growth of amorphous aluminosilicate nanoparticles via an aggregation process.

The precipitation of zeolite nanoparticles involves the initial formation of metastable precursors, such as amorphous entities, that crystallize through non-classical pathways. Here, using reactive force field-based simulations, we reveal how aluminosilicate oligomers grow concomitantly to the decondensation of silicate entities during the initial step of the reaction. Aluminate clusters first form in the solution, thus violating the Loewenstein rule in the first instant of the reaction, which is then followed by their connection with silicate oligomers at the terminal silanol groups before reorganization to finally diffuse within the silicate oligomers to form stable amorphous aluminosilicate nanoparticles that do obey the Loewenstein rule. Our results clearly indicate that aluminate does not serve as the nucleation center for the growth of aluminosilicates in a nucleation-like process but rather proceeds via an aggregation process. The coexistence of aluminosilicate oligomers and small silicate entities induces a phase separation that promotes the precipitation of zeolites with aging.

Atomistic level aqueous dissolution dynamics of NASICON-Type Li1+xAlxTi2-x(PO4)3 (LATP).

Advancing the atomistic level understanding of aqueous dissolution of multicomponent materials is essential. We combined ReaxFF and experiments to investigate the dissolution at the Li1+xAlxTi2-x(PO4)3-water interface. We demonstrate that surface dissolution is a sequentially dynamic process. The phosphate dissolution destabilizes the NASICON structure, which triggers a titanium-rich secondary phase formation.

ReaxFF Force Field Development for Gas-Phase hBN Nanostructure Synthesis.

Two-dimensional (2D) hexagonal boron nitride materials are isomorphs of carbon nanomaterials and hold promise for electronics applications owing to their unique properties. Despite the recent advances in synthesis, the current production capacity for boron nitride (BN) nanostructures is far behind that for carbon-based nanostructures. Understanding the growth mechanism of BN nanostructures through modeling and experiments is key to improving this situation. In the current work, we present the development of a ReaxFF-based force field capable of modeling the gas-phase chemistry important for the chemical vapor deposition (CVD) synthesis process. This force field is parameterized to model the boron nitride nanostructure (BNNS) formation in the gas phase using BN and HBNH as precursors. Our ReaxFF simulations show that BN is the best of these two precursors in terms of quality and the size of BNNSs. The BN precursors lead to the formation of closed BNNSs. However, BNNSs are replaced with complex polymeric structures at temperatures of 2500 K and higher due to entropic effects. Compared to the BN precursors, the HBNH precursors form relatively small, flat, and low-quality BNNSs, but this structure is less affected by temperature. Additives like H2 significantly affect the BNNS formation by preventing closed BNNS formation. Our results show the ReaxFF capability in predicting the BN gas-phase chemistry and BNNS formation, thus providing key insights for experimental synthesis.

Dynamics of the Chemically Driven Densification of Barium Titanate Using Molten Hydroxides.

Molten hydroxides, often used for crystal growth and nanoparticle synthesis, have recently been applied for the single step densification of several inorganic materials under moderate uniaxial pressures and 1000 °C below their usual sintering temperatures. The latter approach, termed cold sintering process (CSP), is a mechanochemically driven process that enables the densification of inorganic materials through a dissolution-precipitation creep mechanism. In this study, we report the main densification mechanisms of BaTiO3 in a NaOH-KOH eutectic mixture. A chemical insight at the atomistic level, investigated by ReaxFF molecular dynamics simulations, offers plausible ionic complex formation scenarios and reactions at the BaTiO3/molten hydroxide interface, enabling the dissolution-precipitation reactions and the subsequent cold sintering of BaTiO3.

Illuminating Invisible Grain Boundaries in Coalesced Single-Orientation WS2 Monolayer Films.

Engineering atomic-scale defects is crucial for realizing wafer-scale, single-crystalline transition metal dichalcogenide monolayers for electronic devices. However, connecting atomic-scale defects to larger morphologies poses a significant challenge. Using electron microscopy and ReaxFF reactive force field-based molecular dynamics simulations, we provide insights into WS2 crystal growth mechanisms, providing a direct link between synthetic conditions and microstructure. Dark-field TEM imaging of coalesced monolayer WS2 films illuminates defect arrays that atomic-resolution STEM imaging identifies as translational grain boundaries. Electron diffraction and high-resolution imaging reveal that the films have nearly a single orientation with imperfectly stitched domains that tilt out-of-plane when released from the substrate. Imaging and ReaxFF simulations uncover two types of translational mismatch, and we discuss their origin related to relatively fast growth rates. Statistical analysis of >1300 facets demonstrates that microstructural features are constructed from nanometer-scale building blocks, describing the system across sub-Ångstrom to multimicrometer length scales.