Atomistic Scale Analysis of the Carbonization Process for C/H/O/N-Based Polymers with the ReaxFF Reactive Force Field

Porous carbon nitride fullerenes: a novel family of porous cage molecules

A novel family of cage molecules, porous carbon nitride fullerenes (PCNFs), is designed, proposed and studied theoretically. PCNFs can be considered the zero-dimensional counterparts of two-dimensional porous graphitic carbon nitrides, in analogy with icosahedral fullerenes, being the zero-dimensional counterparts of graphene. The study is focused on two representative members of the PCNF family: icosahedral C60N60 and C120N60, which are the first members of the two main sub-families of these structures. Given the advanced potential of two-dimensional graphitic carbon nitrides for several interesting applications, it is reasonable to expect that this potential extends to their zero-dimensional counterparts. The present study demonstrates the electronic, vibrational, and thermal stability of the two representative PCNFs utilizing density functional theory and molecular dynamics simulations with ReaxFF potentials. In addition, their structural, vibrational, and electronic properties are revealed.

In situ construction of multifunctional femtosecond laser-induced graphene on arbitrary substrates

The construction of laser-induced graphene (LIG) on various substrates is important for expanding new applications. However, current LIG transfer technologies are hampered by limited substrates, complicated processes, induced graphene defects, and single function. Herein, a facile laser processing method is proposed for in situ construction of multifunctional femtosecond laser-induced graphene (FsLIG) on arbitrary substrates utilizing femtosecond laser acting on polyimide tape. Unlike previous LIG transfer research, the proposed method is applicable to any substrates without introducing additional graphene defects, while also exhibiting multifunctionality. Raman spectra results confirm successful fabrication of FsLIG on various substrates involving paper, aluminum, ceramic, and silicon. Taking paper for example, FsLIG demonstrates multifunctional characteristics including high water contact angle (∼153.4°), large absorptance (∼98.8%), low sheet resistance (∼82.0 Ω sq-1), and reliable temperature sensing (∼-0.089% °C-1) properties. Our study provides a reliable pathway for fabricating multifunctional FsLIG on arbitrary substrates.

Synergistic mechanism of polymer and asphaltene in supercritical water oxidation treatment of polymer-containing oily sludge: Degradation pathway and nitrogen transfer

Asphaltenes and polymers constitute critical components governing the environmentally sound treatment of polymer-containing oily sludge, yet their synergistic mechanisms during supercritical water oxidation (SCWO) remain poorly understood. There are no SCWO studies addressing the synergistic mechanisms between multiple components. This study systematically investigates the synergistic oxidation mechanism of the asphaltene/polyacrylamide (PAM) system in supercritical water using ReaxFF molecular dynamics simulations, with a focus on product distribution, nitrogen migration pathways, degradation mechanisms, and free radical reactions. The result shows that the synergy emerges through two mechanisms: first, PAM reacts earlier than asphaltene, creating uneven free radical distribution; second, PAM fragments actively attack asphaltene structures. The presence of PAM, on the one hand, preferentially promoted the conversion of -COOH to CO2 rather than CO; on the other hand, the intermediate fragment produced by the early cleavage of PAM shortened and reconfigured the path of the ring-opening degradation of asphaltenes. Elevated temperatures reduced the temporal difference in reaction initiation between PAM and asphaltene, weakened the inhomogeneous distribution of free radicals, and strengthened the reaction competition between the two. Higher PAM/asphaltene ratios enhanced degradation efficiency through intensified fragment-driven attack mechanisms.

Enhancement of Fracture Toughness of Inner Liner Material for Type IV Hydrogen Storage Cylinders Based on Molecular Dynamics Method

To develop liner materials with improved toughness, this study combines molecular dynamics simulations and experimental testing to investigate the effect of different mass ratios (10/0, 7/3, 6/4, 4/6, 3/7, and 0/10) of high-density polyethylene (HDPE)/polyamide 6 (PA6) on their fracture toughness of the composites. The fracture toughness was quantitatively assessed using the J-integral method, while the material's behavior in terms of crack propagation during tensile deformation was examined at the molecular level. The results reveal that as the HDPE mass ratio increases, the fracture toughness of the composites also gradually improves. Furthermore, the fracture toughness of four materials (PA6, 4HDPE/6PA6, 7HDPE/3PA6, and HDPE) was tested using the essential work of the fracture method. The trend observed in the simulation results was in agreement with the experimental results, validating the reliability of the molecular dynamics simulation.

Comparison between molecular dynamics potentials for simulation of graphene-based nanomaterials for biomedical applications

OBJECTIVE

This article provides a substantial review of recent research and comparison on molecular dynamics potentials to determine which are most suitable for simulating the phenomena in graphene-based nanomaterials (GBNs).

SIGNIFICANCE

GBNs gain significant attention due to their remarkable properties and potential applications, notably in nanomedicine. However, the physical and chemical characteristics toward macromolecules that justify their nanomedical applications are not yet fully understood. The molecular interaction through molecular dynamic simulation offers the benefits for simulating inorganic molecules like GBNs, with necessary adjustments to account for physical and chemical interactions, or thermodynamic conditions.

METHOD

In this review, we explore various molecular dynamics potentials (force fields) used to simulate interactions and phenomena in graphene-based nanomaterials. Additionally, we offer a brief overview of the benefits and drawbacks of each force fields that available for analysis to assess which one is suitable to study the molecular interaction of graphene-based nanomaterials.

RESULT

We identify and compare various molecular dynamics potentials that available for analyzing GBNs, providing insights into their suitability for simulating specific phenomena in graphene-based nanomaterials. The specification of each force fields and its purpose can be used for further application of molecular dynamics simulation on GBNs.

CONCLUSION

GBNs hold significant promise for applications like nanomedicine, but their physical and chemical properties must be thoroughly studied for safe clinical use. Molecular dynamics simulations, using either reactive or non-reactive MD potentials depending on the expected chemical changes, are essential for accurately modeling these properties, requiring careful selection based on the specific application.

Noncatalytic Reduction of Nitrogen Oxide and Soot Inhibition during Cocombustion with Biotar: The Molecular Dynamics Modeling Approach

Cocombustion with biomass tar is a potential method for NO reduction during fossil fuel combustion. In this work, the molecular dynamic method based on the reactive force field was used to study the NO reduction by phenol, which is a typical tar model compound. Results indicate that phenol undergoes significant decomposition at 3000 K, resulting in the formation of small molecular fragments accompanied by the generation of large molecular, network-structured soot particles. At higher temperatures (3500 K), the morphology of the soot produced from phenol decomposition undergoes a certain degree of change. It evolves from a two-dimensional network structure to a three-dimensional particulate structure. Soot particles are significant products of the thermal decomposition process of phenol. NO acts as an oxidant during the phenol decomposition process, significantly inhibiting the formation of soot particles during this process. Phenol also promotes the reduction of NO. The corresponding NO reduction ratios of NO are 52%, 76%, and 89% for 1500, 2000, and 2500 K respectively. CO, H2O, and N2 were the three most important reaction products during phenol-NO interaction. HNO is an important intermediate during the reduction of NO by phenol. The combination of the H radical and NO is the main route for HNO. It can be seen that HNO is quickly produced at the initial stage. The calculated activated energies for NO reduction and phenol oxidation are 30.47 and 50.85 kJ/mol, respectively.

Investigating the structure-property correlations of pyrolyzed phenolic resin as a function of degree of carbonization

Carbon-carbon (C/C) composites are attractive materials for high-speed flights and terrestrial atmospheric reentry applications due to their insulating thermal properties, thermal resistance, and high strength-to-weight ratio. It is important to understand the evolving structure-property correlations in these materials during pyrolysis, but the extreme laboratory conditions required to produce C/C composites make it difficult to quantify the properties in situ. This work presents an atomistic modeling methodology to pyrolyze a crosslinked phenolic resin network and track the evolving thermomechanical properties of the skeletal matrix during simulated pyrolysis. First, the crosslinked resin is pyrolyzed and the resulting char yield and mass density are verified to match experimental values, establishing the model's powerful predictive capabilities. Young's modulus, yield stress, Poisson's ratio, and thermal conductivity are calculated for the polymerized structure, intermediate pyrolyzed structures, and fully pyrolyzed structure to reveal structure-property correlations, and the evolution of properties are linked to observed structural features. It is determined that reduction in fractional free volume and densification of the resin during pyrolysis contribute significantly to the increase in thermomechanical properties of the skeletal phenolic matrix. A complex interplay of the formation of six-membered carbon rings at the expense of five and seven-membered carbon rings is revealed to affect thermal conductivity. Increased anisotropy was observed in the latter stages of pyrolysis due to the development of aligned aromatic structures. Experimentally validated predictive atomistic models are a key first step to multiscale process modeling of C/C composites to optimize next-generation materials.

Computational Insights into Tunable Reversible Network Materials: Accelerated ReaxFF Kinetics of Furan-Maleimide Diels-Alder Reactions for Self-Healing and Recyclability

In this study, ReaxFF molecular dynamics simulations were benchmarked and used to study the relative kinetics of the retro Diels-Alder reaction between furan and N-methylmaleimide. This reaction is very important for the creation of polymer networks with self-healing and recyclable properties, since they can be used as reversible linkers in the network. So far, the reversible Diels-Alder reaction has not yet been studied by using reactive molecular dynamics simulations. This work is, thus, the first step in simulating a covalent adaptable network (CAN) using Diels-Alder reactions as reversible linkers. For both endo and exo, the bond breaking in 40 product molecules was simulated using the bond boost method and the endo/exo ratio was evaluated. This ratio was benchmarked against density functional theory (DFT) and experimental results for a changing set of bond boost parameters. Given their importance to understand how the CAN performs, the effect of the addition of a polymer backbone and the effect of temperature were successfully simulated using our newly parametrized reactive force field.

A Study of the Methane Oxidation Mechanism and Reaction Pathways Using Reactive Molecular Simulation and Nonlinear Manifold Learning

Methane, as the primary component of natural gas, is a vital energy resource extensively utilized through oxidation reactions. These reactions yield diverse radicals and molecules via varying intermediate reaction routes, contingent upon the oxidation conditions. In this study, we employ reactive molecular dynamics simulations to investigate the early-stage mechanism of methane oxidation across different temperatures and methane/oxygen conditions. Our analysis reveals distinct variations in species count, initial reaction times, and the spectrum of the main reactions/molecules under diverse conditions. Notably, both full oxidation of methane (FOM) and partial oxidation of methane (POM) are observed in all simulations, with FOM favored under high-temperature and fuel-lean conditions, while POM prevails in low-temperature and fuel-rich environments. Furthermore, we utilize nonlinear manifold learning techniques to extract a 2D manifold from the reaction state space, identifying two collective variables governing the reaction pathways. This research provides a systematic understanding of the initial stage mechanisms of methane oxidation under varying conditions, offering useful insights into chemical science and fuel engineering.

Degradation pathways and product formation mechanisms of asphaltene in supercritical water

Asphaltene is the compound with the most complex structure and the most difficult degradation in oily sludge, which is the key to limit the efficiency of supercritical water oxidation treatment of oily sludge. In this paper, the supercritical water oxidation process of asphaltene was investigated in terms of free radical reaction, degradation pathway, and product generation mechanism using ReaxFF molecular dynamics simulation method. The results showed that increasing temperature, increasing O2, and increasing H2O have different effects on HO2·generation. Benzene rings undergo fusion and condensation through hydrogenation abstraction and oxygen addition reactions, subsequently breaking down into long-chain alkanes. Increasing O2 can effectively promote the ring-opening of nitrogen-containing heterocycles. -COOH is the most important intermediate fragment for CO and CO2 generation, and there is a reaction competition with -CHO3 and -CO3. When the number of oxygen molecules increases from 300 to 700, the reaction frequency of -CHO3 and -CO3 to generate CO and CO2 increases by 17.14 % and 12.77 %·H2O determines the production of H2 by controlling the number of H·radicals present. As the amount of H2O increases from 500 to 1500, the product ratio of H2 increases from 12.73 % to 21.31 %. ENVIRONMENTAL IMPLICATION: Asphaltene is the most structurally complex organic matter in oily sludge, and its presence makes it difficult for oily sludge to be completely degraded by conventional treatment methods such as pyrolysis and incineration. Polycyclic aromatic hydrocarbons (PAHs) represented by asphaltene increase the carcinogenicity and mutagenicity of oily sludge, and even irreversibly pollute soil and groundwater. Supercritical water oxidation, as an efficient organic waste treatment technology, can realize harmlessness in a green and efficient way. So the study on the mechanism of supercritical water oxidation of asphaltene is of great significance for environmental protection.

Mechanistic Insights into the Evolution of Graphitic Carbon from Sulfur Containing Polar Aromatic Feedstock

In the realm of carbon fiber research, a variety of structural configurations is noted, comprising crystalline, noncrystalline, and semicrystalline forms. Recent investigations into this domain have revealed an array of intriguing phases of carbon, among which amorphous graphite is the most notable for its unique mechanical, thermal, and electrical properties that arise from its inherent topological disorders. In this study, we utilized the ReaxFF molecular dynamics (MD) simulations to investigate the carbonization and graphitization processes involved in the production of amorphous graphite from benzothiophene, a sulfur-containing polar aromatic precursor. We developed C/H/S ReaxFF force field parameters to describe the high-temperature chemistry of benzothiophene. Our investigation reveals the reaction mechanisms, providing critical insights into the underlying chemical processes toward the formation of amorphous graphite and the structural characteristics of the end products. The formation of volatile gaseous molecules and their continuous elimination led to the development of noncontinuous layered graphite structures analogous to amorphous graphite consisting of pentagons, hexagons, and heptagons. These findings offer unprecedented insights into the carbonization and graphitization processes of sulfur-containing heavy-end aromatic feedstock. This knowledge lays the groundwork for advancing synthesis methods and developing amorphous graphite materials with specific properties.

Carbon-based Flame Retardants for Polymers: A Bottom-up Review

This state-of-the-art review is geared toward elucidating the molecular understanding of the carbon-based flame-retardant mechanisms for polymers via holistic characterization combining detailed analytical assessments and computational material science. The use of carbon-based flame retardants, which include graphite, graphene, carbon nanotubes (CNTs), carbon dots (CDs), and fullerenes, in their pure and functionalized forms are initially reviewed to evaluate their flame retardancy performance and to determine their elevation of the flammability resistance on various types of polymers. The early transition metal carbides such as MXenes, regarded as next-generation carbon-based flame retardants, are discussed with respect to their superior flame retardancy and multifunctional applications. At the core of this review is the utilization of cutting-edge molecular dynamics (MD) simulations which sets a precedence of an alternative bottom-up approach to fill the knowledge gap through insights into the thermal resisting process of the carbon-based flame retardants, such as the formation of carbonaceous char and intermediate chemical reactions offered by the unique carbon bonding arrangements and microscopic in-situ architectures. Combining MD simulations with detailed experimental assessments and characterization, a more targeted development as well as a systematic material synthesis framework can be realized for the future development of advanced flame-retardant polymers.

Controllable preparation of graphene glass fiber fabric towards mass production and its application in self-adaptive thermal management

Direct synthesis of graphene on nonmetallic substrates via chemical vapor deposition (CVD) has become a frontier research realm targeting transfer-free applications of CVD graphene. However, the stable mass production of graphene with a favorable growth rate and quality remains a grand challenge. Herein, graphene glass fiber fabric (GGFF) was successfully developed through the controllable growth of graphene on non-catalytic glass fiber fabric, employing a synergistic binary-precursor CVD strategy to alleviate the dilemma between growth rate and quality. The binary precursors consisted of acetylene and acetone, where acetylene with high decomposition efficiency fed rapid graphene growth while oxygen-containing acetone was adopted for improving the layer uniformity and quality. Notably, the bifurcating introducing-confluent premixing (BI-CP) system was self-built for the controllable introduction of gas and liquid precursors, enabling the stable production of GGFF. GGFF features solar absorption and infrared emission properties, based on which the self-adaptive dual-mode thermal management film was developed. This film can automatically switch between heating and cooling modes by spontaneously perceiving the temperature, achieving excellent thermal management performances with heating and cooling power of ∼501.2 and ∼108.6 W m-2, respectively. These findings unlock a new strategy for the large-scale batch production of graphene materials and inspire advanced possibilities for further applications.

Size-dependent shock response mechanisms in nanogranular RDX: a reactive molecular dynamics study

Understanding the shock initiation mechanisms of explosives is pivotal for advancing physicochemical theories and enhancing experimental methodologies. This study delves into the size-dependent shock responses of nanogranular hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) through nonequilibrium reactive molecular dynamics simulations. Utilizing the ReaxFF-lg force field, we examine the influence of the particle size on the decomposition dynamics of RDX under varying shock velocities. Our findings reveal that larger particles promote more significant RDX decomposition at lower velocities due to fluid jet formation and gas compression during void collapse. Conversely, smaller particles exhibit a higher average temperature and a faster decomposition rate under high-velocity shocks, attributed to their increased specific surface area. Detailed chemical reaction pathways are analyzed to elucidate the growth and initiation of reactions during shock waves. The results contribute to resolving the discrepancies observed in experimental studies of shocked granular explosives and provide a deeper understanding of the underlying mechanisms governing their behavior. This research offers valuable insights into the design and control of nano- and submicron-sized explosives with tailored sensitivity to external stimuli.

The impact of molecular configuration on the bond breaking rates of hydrocarbons: a computational study

The dissociation of hydrocarbon bonds plays a pivotal role in their utilization, whether through fuel combustion or the thermo-cracking of large hydrocarbons in petroleum refinement. Previous studies have primarily focused on the effects of temperature, pressure, and chemical environment on hydrocarbon reactions. However, the influence of molecular configuration on bond breaking rates has not been thoroughly explored. In this study, we propose an approach to compute bond dissociation rates, and apply it to the reactive molecular dynamics simulation (ReaxFF) trajectories of three molecules: n-tridecane, n-pentane, and 1,3-propanediol. Our results reveal that the bond dissociation rate depends not only on the bond position in the chain, but also on the molecular configuration. Stretched configurations exhibit higher dissociation rates, particularly favoring the breaking of central bonds. Conversely, when the molecule is coiled, resulting in a reduced size, terminal bonds exhibit higher dissociation rates. This research contributes to a deeper understanding of molecular dissociation properties in the oxidation of hydrocarbons, and provides a way to explore the bond breaking properties of other molecules.

Insights into hydro thermal gasification process of microplastic polyethylene via reactive molecular dynamics simulations

Microplastics have become a pressing environmental issue due to their widespread presence in our ecosystems. Among various plastic components, polyethylene (PE) is a prevalent and persistent contaminant. Hydrothermal gasification (HTG), a promising technology for converting PE into syngas, holds great promise for mitigating the microplastic problem. In this study, we employ ReaxFF molecular dynamics simulations to investigate the HTG process of PE, shedding light on the intricate relationships between temperature, water content, carbon conversion efficiency, and product distributions. The results reveal that hydrothermal gasification of PE is a complex process involving multiple reaction pathways. Consistently with experimental findings, the calculations indicate that the gas phase exhibits a substantial hydrogen fraction, reaching up to 70%. Interestingly, our simulations reveal a dual role of water content in the HTG process. On one hand, water enhances hydrogen production by promoting the gas formation. On the other hand, it elevates the activation energy required for PE decomposition. Depending on the water content, the calculated activation energies range from 176 to 268 kJ/mol, which are significantly lower than those reported for thermal gasification (TG). This suggests that HTG may be a more efficient route for PE conversion. Furthermore, this study highlights the importance of optimizing both temperature and water content in HTG systems to achieve high yields of hydrogen-rich syngas. The results obtained from our ReaxFF MD simulations demonstrate the robustness of this computational methodology in elucidating complex chemical reactions under extreme conditions. Our findings offer critical insights into the design of advanced waste management strategies for microplastics and contribute to the development of sustainable practices for resource recovery. This work underscores the potential of HTG as a key technology for addressing the global challenge of plastic pollution.

Biomimetic Multifunctional Graphene-Based Coating for Thermal Management, Solar De-Icing, and Fire Safety: Inspired from the Antireflection Nanostructure of Compound Eyes

Due to the ubiquitous and inexhaustible solar source, photothermal materials have gained considerable attention for their potential in heating and de-icing. Nevertheless, traditional photothermal materials, exemplified by graphene, frequently encounter challenges emanating from their elevated reflectance. Inspired by ocular structures, this study uses the Fresnel equation to enhance the photo-thermal conversion efficiency of graphene by introducing a polydimethylsiloxane (PDMS)/silicon dioxide (SiO2) coating, which reduces the light reflectance (≈20%) through destructive interference. The designed coating achieves an equilibrium temperature of ≈77 °C at one sun and a quick de-icing in ≈65 s, all with a thickness of 5 µm. Simulations demonstrate that applying this coating to high-rise buildings results in energy savings of ≈31% in winter heating. Furthermore, the combination of PDMS/SiO2 and graphene confers a notable enhancement in thermal stability through a synergistic flame-retardant mechanism, effectively safeguarding polyurethane against high temperatures and conflagrations, leading to marked reduction of 58% and 28% in heat release rate and total heat release. This innovative design enhances the photo-thermal conversion, de-icing function, and flame retardancy of graphene, thereby advancing its applications in outdoor equipment, high-rise buildings, and aerospace vessels.

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole-dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. This study enables the ability to predict the structure-property relations for PAN and analogues with sustainable ab initio accuracy across scales.

LUNAR: Automated Input Generation and Analysis for Reactive LAMMPS Simulations

Generating simulation-ready molecular models for the LAMMPS molecular dynamics (MD) simulation software package is a difficult task and impedes the more widespread and efficient use of MD in materials design and development. Fixed-bond force fields generally require manual assignment of atom types, bonded interactions, charges, and simulation domain sizes. A new LAMMPS pre- and postprocessing toolkit (LUNAR) is presented that efficiently builds molecular systems for LAMMPS. LUNAR automatically assigns atom types, generates bonded interactions, assigns charges, and provides initial configuration methods to generate large molecular systems. LUNAR can also incorporate chemical reactivity into simulations by facilitating the use of the REACTER protocol. Additionally, LUNAR provides postprocessing for free volume calculations, cure characterization calculations, and property predictions from LAMMPS thermodynamic outputs. LUNAR has been validated via building and simulation of pure epoxy and cyanate ester polymer systems with a comparison of the corresponding predicted structures and properties to benchmark values, including experimental results from the literature. LUNAR provides the tools for the computationally driven development of next-generation composite materials in the Integrated Computational Materials Engineering (ICME) and Materials Genome Initiative (MGI) frameworks. LUNAR is written in Python with the usage of NumPy and can be used via a graphical user interface, a command line interface, or an integrated design environment. LUNAR is freely available via GitHub.

Fabricating Ultrastrong Carbon Nanotube Fibers via a Microwave Welding Interface

Liquid crystal wet-spun carbon nanotube fibers (CNTFs) offer notable advantages, such as precise alignment and scalability. However, these CNTFs usually suffer premature failure through intertube slippage due to the weak interfacial interactions between adjacent shells and bundles. Herein, we present a microwave (MW) welding strategy to enhance intertube interactions by partially carbonizing interstitial heterocyclic aramid polymers. The resulting CNTFs exhibit ultrahigh static tensile strength (6.74 ± 0.34 GPa) and dynamic tensile strength (9.52 ± 1.31 GPa), outperforming other traditional high-performance fibers. This work provides a straightforward yet effective approach to strengthening CNTFs for advanced engineering applications.

Dynamic Insights into the Growth Mechanisms of 2D Covalent Organic Frameworks on Graphene Surfaces

Producing high-quality two-dimensional (2D) covalent organic frameworks (COFs) is crucial for industrial applications. However, this remains significantly challenging with current synthetic techniques. A deep understanding of the intermolecular interactions, reaction temperature, and oligomers is essential to facilitate the growth of highly crystalline COF films. Herein, molecular dynamics simulations were employed to explore the growth of 2D COFs from monomer assemblies on graphene. Our results showed that chain growth reactions dominated the COF surface growth and that van der Waals (vdW) interactions were important in enhancing the crystallinity through monomer preorganization. Moreover, appropriately tuning the reaction temperature improved the COF crystallinity and minimized the effects of amorphous oligomers. Additionally, the strength of the interface between the COF and the graphene substrate indicated that the adhesion force was proportional to the crystallinity of the COF. This work reveals the mechanisms for nucleation and growth of COFs on surfaces and provides theoretical guidance for fabricating high-quality 2D polymer-based crystalline nanomaterials.

Effect of Hyperthermal Hybrid Gas Composition on the Interfacial Oxidation and Nitridation Mechanisms of Graphene Sheet

Graphene is one of the most promising thermal protection materials for high-speed aircraft due to its lightweight and excellent thermophysical properties. At high Mach numbers, the extremely high postshock temperature would dissociate the surrounding air into a mixture of atomic and molecular components in a highly thermochemical nonequilibrium state, which greatly affects the subsequent thermal chemical reactions of the graphene interface. Through establishing a reactive gas-solid interface model, the reactive molecular dynamics method is employed in this study to reveal the influences of the thermochemical nonequilibrium gas mixture on the thermal oxidation and nitridation mechanisms of graphene sheet. The results show that three distinctive stages can be recognized during bombardment of various nonequilibrium gas components toward the graphene sheet: (i) collision and surface adsorption stage, (ii) gas-solid heterogeneous reaction stage, and (iii) gas phase homogeneous reaction stage. The surface catalysis effect is found to be dominant during the first two stages, which can influence the following ablation behavior of graphene significantly at high-temperature conditions. Moreover, surface catalysis, oxidation, nitridation, and ablation mechanisms between nonequilibrium gas and graphene interface are revealed, which is of high relevance for future interfacial design and application of graphene as a thermal protection material.

Simulation and experimental evaluation of laser-induced graphene on the cellulose and lignin substrates

In this article, the formation of laser-induced graphene on the two natural polymers, cellulose, and lignin, as precursors was investigated with molecular dynamics simulations and some experiments. These eco-friendly polymers provide significant industrial advantages due to their low cost, biodegradability, and recyclable aspects. It was discovered during the simulation that LIG has numerous defects and a porous structure. Carbon monoxide, H2, and water vapor are gases released by cellulose and lignin substrates. H2O and CO are released when the polymer transforms into an amorphous structure. Later on, as the amorphous structure changes into an ordered graphitic structure, H2 is released continuously. Since cellulose monomer has a higher mass proportion of oxygen (49%) than lignin monomer (29%), it emits more CO. The LIG structure contains many 5- and 7-carbon rings, which cause the structure to have bends and undulations that go out of the plane. In addition, to verify the molecular dynamics simulation results with experimental tests, we used a carbon dioxide laser to transform filter paper, as a cellulose material, and coconut shell, as a lignin material, into graphene. Surprisingly, empirical experiments confirmed the simulation results.

Atomistically informed hierarchical modeling for revisiting the constituent structures from heredity and nano-micro mechanics of sheath-core carbon fiber

To better understand the heterogeneous anisotropic nanocomposite features and provide reliable underlying constitutive parameters of carbon fiber for continuum-level simulations, hierarchical modeling approaches combining quantum chemistry, molecular dynamics, numerical and analytical micromechanics are employed for studying the structure-performance relationships of the precursor-inherited sheath-core carbon fiber layers. A robust debonding force field is derived from energy matching protocols, including bond dissociation enthalpy calculations and rigid-constraint potential energy surface scan. Logistic long range bond stretching curves with exponential parameters and shifted force vdW curves are designed to diminish energy perturbations. The pseudo-crystalline microstructure is proposed and validated using virtual wide angle X-ray diffraction patterns and bond-orientational order parameters. The distribution or alignment features of the nanocomposite microstructures are collected from quantum chemical topology analysis and normal vector extractions. Non-equilibrium tensile loading simulation predicts the decomposed strain energy contributions, principal-axis modulus, strength limit, localized stress, and fracture morphologies of the model. Finally, an atomistically-informed stiffness prediction model combining numerical homogenization and analytical self-consistent Eshelby-Mori-Tanaka-type effective mean field micromechanics theory is proposed, giving a successful estimation of the overall stiffness matrix of the sheath-core carbon fiber system. The hierarchical models in combination with the carbonization reaction template will help in providing efficient and feasible schemes for the synergistic process-performance control of distinct types of carbon fiber.

Reactive Molecular Dynamics Simulations of Polystyrene Pyrolysis

Polymers' controlled pyrolysis is an economical and environmentally friendly solution to prepare activated carbon. However, due to the experimental difficulty in measuring the dependence between microstructure and pyrolysis parameters at high temperatures, the unknown pyrolysis mechanism hinders access to the target products with desirable morphologies and performances. In this study, we investigate the pyrolysis process of polystyrene (PS) under different heating rates and temperatures employing reactive molecular dynamics (ReaxFF-MD) simulations. A clear profile of the generation of pyrolysis products determined by the temperature and heating rate is constructed. It is found that the heating rate affects the type and amount of pyrolysis intermediates and their timing, and that low-rate heating helps yield more diverse pyrolysis intermediates. While the temperature affects the pyrolytic structure of the final equilibrium products, either too low or too high a target temperature is detrimental to generating large areas of the graphitized structure. The reduced time plots (RTPs) with simulation results predict a PS pyrolytic activation energy of 159.74 kJ/mol. The established theoretical evolution process matches experiments well, thus, contributing to preparing target activated carbons by referring to the regulatory mechanism of pyrolytic microstructure.

Artificial Neural Networks for Predicting Mechanical Properties of Crystalline Polyamide12 via Molecular Dynamics Simulations

Predicting material properties of 3D printed polymer products is a challenge in additive manufacturing due to the highly localized and complex manufacturing process. The microstructure of such products is fundamentally different from the ones obtained by using conventional manufacturing methods, which makes the task even more difficult. As the first step of a systematic multiscale approach, in this work, we have developed an artificial neural network (ANN) to predict the mechanical properties of the crystalline form of Polyamide12 (PA12) based on data collected from molecular dynamics (MD) simulations. Using the machine learning approach, we are able to predict the stress-strain relations of PA12 once the macroscale deformation gradient is provided as an input to the ANN. We have shown that this is an efficient and accurate approach, which can provide a three-dimensional molecular-level anisotropic stress-strain relation of PA12 for any macroscale mechanics model, such as finite element modeling at arbitrary quadrature points. This work lays the foundation for a multiscale finite element method for simulating semicrystalline polymers, which will be published as a separate study.

Evolution of Glassy Carbon Derived from Pyrolysis of Furan Resin

Glassy carbon (GC) material derived from pyrolyzed furan resin was modeled by using reactive molecular dynamics (MD) simulations. The MD polymerization simulation protocols to cure the furan resin precursor material are validated via comparison of the predicted density and Young's modulus with experimental values. The MD pyrolysis simulations protocols to pyrolyze the furan resin precursor is validated by comparison of calculated density, Young's modulus, carbon content, sp2 carbon content, the in-plane crystallite size, out-of-plane crystallite stacking height, and interplanar crystallite spacing with experimental results from the literature for furan resin derived GC. The modeling methodology established in this work can provide a powerful tool for the modeling-driven design of next-generation carbon-carbon composite precursor chemistries for thermal protection systems and other high-temperature applications.

Effect of Methyl Groups on Formation of Ordered or Layered Graphitic Materials from Aromatic Molecules

Developing functionally complex carbon materials from small aromatic molecules requires an understanding of how the chemistry and structure of its constituent molecules evolve and crosslink, to achieve a tailorable set of functional properties. Here, molecular dynamics (MD) simulations are used to isolate the effect of methyl groups on condensation reactions during the oxidative process and evaluate the impact on elastic modulus by considering three monodisperse pyrene-based systems with increasing methyl group fraction. A parameter to quantify the reaction progression is designed by computing the number of new covalent bonds formed. Utilizing the previously developed MD framework, it is found that increasing methylation leads to an almost doubling of bond formation, a larger fraction of the new bonds oriented in the direction of tensile stress, and a higher basal plane alignment of the precursor molecules along the direction of tensile stress, resulting in enhanced tensile modulus. Additionally, via experiments, it is demonstrated that precursors with a higher fraction of methyl groups result in a higher alignment of molecules. Moreover, increased methylation results in the lower spread of single molecule alignment which may lead to smaller variations in tensile modulus and more consistent properties in carbon materials derived from methyl-rich precursors.

Joint Theoretical and Experimental Study of Stress Graphitization in Aligned Carbon Nanotube/Carbon Matrix Composites

Stress graphitization is a unique phenomenon at the carbon nanotube (CNT)-matrix interfaces in CNT/carbon matrix (CNT/C) composites. A lack of fundamental atomistic understanding of its evolution mechanisms and a gap between the theoretical and experimental research have hindered the pursuit of utilizing this phenomenon for producing ultrahigh-performance CNT/C composites. Here, we performed reactive molecular dynamics simulations along with an experimental study to explore stress graphitization mechanisms of a CNT/polyacrylonitrile (PAN)-based carbon matrix composite. Different CNT contents in the composite were considered, while the nanotube alignment was controlled in one direction in the simulations. We observe that the system with a higher CNT content exhibits higher localized stress concentration in the periphery of CNTs, causing alignment of the nitrile groups in the PAN matrix along the CNTs, which subsequently results in preferential dehydrogenation and clustering of carbon rings and eventually graphitization of the PAN matrix when carbonized at 1500 K. These simulation results have been validated by experimentally produced CNT/PAN-based carbon matrix composite films, with transmission electron microscopy images showing the formation of additional graphitic layers converted by the PAN matrix around CNTs, where 82 and 144% improvements of the tensile strength and Young's modulus are achieved, respectively. The presented atomistic details of stress graphitization can provide guidance for further optimizing CNT-matrix interfaces in a more predictive and controllable way for the development of novel CNT/C composites with high performance.

Formation of nanoribbons by carbon atoms confined in a single-walled carbon nanotube-A molecular dynamics study

Carbon nanotubes can serve as one-dimensional nanoreactors for the in-tube synthesis of various nanostructures. Experimental observations have shown that chains, inner tubes, or nanoribbons can grow by the thermal decomposition of organic/organometallic molecules encapsulated in carbon nanotubes. The result of the process depends on the temperature, the diameter of the nanotube, and the type and amount of material introduced inside the tube. Nanoribbons are particularly promising materials for nanoelectronics. Motivated by recent experimental results observing the formation of carbon nanoribbons inside carbon nanotubes, molecular dynamics calculations were performed with the open source LAMMPS code to investigate the reactions between carbon atoms confined within a single-walled carbon nanotube. Our results show that the interatomic potentials behave differently in quasi-one-dimensional simulations of nanotube-confined space than in three-dimensional simulations. In particular, the Tersoff potential performs better than the widely used Reactive Force Field potential in describing the formation of carbon nanoribbons inside nanotubes. We also found a temperature window where the nanoribbons were formed with the fewest defects, i.e., with the largest flatness and the most hexagons, which is in agreement with the experimental temperature range.

Insights into the carbonization mechanism of PAN-derived carbon precursor fibers and establishment of a kinetics-driven accelerated reaction template for atomistic simulation

To better understand the chemistry behind the carbonization process of the polyacrylonitrile (PAN)-based precursor fibers and provide a more authentic virtual counterpart of the process-inherited model for process optimization and rational performance design, we develop arrow-pushing reaction routes for primary exhaust gas product (H₂O/H₂/HCN/N₂/tar vapor) formation and a pragmatic kinetics-driven accelerated reaction template for atomistic simulation of the carbonization process overcoming traditional challenges in time scale discrepancy of the reaction-diffusion system. The results of enthalpy barriers from hybrid first principles calculations validate the rationality and sequence of conjectured reactions during the two-stage carbonization process. Conversion rates of the rate-determining steps under 300 s carbonization are also estimated based on Eyring's transition state theory realizing kinetics equivalency of the reaction extent. Process-control measurements are further demonstrated corresponding to the proposed mechanism. The iterative densified crosslinking scheme specially designed for the surface layer is implanted into the topological reaction molecular dynamics template and a series of highly devisable structural models during the whole evolutionary process from the pre-oxidized fiber to the pristine carbon fiber surface are successfully predicted. The ultimate structure of the model presents excellent similarity in carbon yield and elemental composition with the type II high strength carbon fiber surface.

Depolymerization of plastics by means of electrified spatiotemporal heating

Depolymerization is a promising strategy for recycling waste plastic into constituent monomers for subsequent repolymerization1. However, many commodity plastics cannot be selectively depolymerized using conventional thermochemical approaches, as it is difficult to control the reaction progress and pathway. Although catalysts can improve the selectivity, they are susceptible to performance degradation2. Here we present a catalyst-free, far-from-equilibrium thermochemical depolymerization method that can generate monomers from commodity plastics (polypropylene (PP) and poly(ethylene terephthalate) (PET)) by means of pyrolysis. This selective depolymerization process is realized by two features: (1) a spatial temperature gradient and (2) a temporal heating profile. The spatial temperature gradient is achieved using a bilayer structure of porous carbon felt, in which the top electrically heated layer generates and conducts heat down to the underlying reactor layer and plastic. The resulting temperature gradient promotes continuous melting, wicking, vaporization and reaction of the plastic as it encounters the increasing temperature traversing the bilayer, enabling a high degree of depolymerization. Meanwhile, pulsing the electrical current through the top heater layer generates a temporal heating profile that features periodic high peak temperatures (for example, about 600 °C) to enable depolymerization, yet the transient heating duration (for example, 0.11 s) can suppress unwanted side reactions. Using this approach, we depolymerized PP and PET to their monomers with yields of about 36% and about 43%, respectively. Overall, this electrified spatiotemporal heating (STH) approach potentially offers a solution to the global plastic waste problem.

ReaxFF molecular dynamics simulation of the thermal decomposition reaction of bio-based polyester materials

Bio-based polyester elastomers have been widely studied by researchers in recent years because of their comprehensive sources of monomers and environmentally friendly characteristics. However, compared with traditional petroleum-based elastomers, the thermal decomposition temperature of bio-based polyester elastomers is generally low, limiting the application of bio-based elastomers. An effective strategy to increase the intrinsic thermal decomposition temperature (T_d) of bio-based elastomers is to increase the length of the monomer carbon chain in the bio-based elastomers. In this work, the content of dodecanedioic acid (DDA) in a bio-based polyester elastomer composed of butanediol (BDO) and succinic acid (SUA) was increased to improve the T_d of the bio-based polyester elastomer through the reaction force-field molecular dynamics (ReaxFF-MD) simulations. And the thermal decomposition mechanism of the bio-based polyester was analyzed in detail. By calculating the change rate of the molecular chain mean square displacement (MSD), it was determined that when the content of DDA was 50%, the T_d of the bio-based elastomer was up to 718 K. By calculating the activation energy of thermal decomposition and further analyzing the thermal decomposition process, it is found that the thermal decomposition of the bio-based polyester elastomer is mainly through breaking the C-O bond in the backbone. This work is expected to provide theoretical guidance for designing and fabricating highly heat-resistant bio-based elastomers by systematically exploring the thermal decomposition mechanism of bio-based polyester elastomers.

Sustainable valorization of asphaltenes via flash joule heating

The refining process of petroleum crude oil generates asphaltenes, which poses complicated problems during the production of cleaner fuels. Following refining, asphaltenes are typically combusted for reuse as fuel or discarded into tailing ponds and landfills, leading to economic and environmental disruption. Here, we show that low-value asphaltenes can be converted into a high-value carbon allotrope, asphaltene-derived flash graphene (AFG), via the flash joule heating (FJH) process. After successful conversion, we develop nanocomposites by dispersing AFG into a polymer effectively, which have superior mechanical, thermal, and corrosion-resistant properties compared to the bare polymer. In addition, the life cycle and technoeconomic analysis show that the FJH process leads to reduced environmental impact compared to the traditional processing of asphaltene and lower production cost compared to other FJH precursors. Thus, our work suggests an alternative pathway to the existing asphaltene processing that directs toward a higher value stream while sequestering downstream emissions from the processing.

Probing Friction Properties of Hydrogen-Free DLC Films in a Nitrogen Environment Based on ReaxFF Molecular Dynamics

In this paper, ReaxFF molecular dynamics simulations were used to look at how load and the number of nitrogen molecules affect how friction behavior in hydrogen-free diamond-like carbon (DLC). The presence of nitrogen molecules will inhibit the formation of C-C covalent bonds between the contact surfaces of the upper and lower DLC, thereby effectively suppressing the increase in friction during the initial friction phase. After the initial friction stage, the mechanical mixing of the contact surfaces brought on by the diffusion of nitrogen molecules results in considerable shear stress, which has significant impacts on the friction force. In addition, due to the existence of nitrogen molecules, the effect of graphitization of hydrogen-free DLC on friction is almost negligible.

Oxidation and hydrogenation of monolayer MoS2 with compositing agent under environmental exposure: The ReaxFF Mo/Ti/Au/O/S/H force field development and applications

An atomistic modeling tool is essential to an in-depth understanding upon surface reactions of transition metal dichalcogenides (TMDs), such as molybdenum disulfide (MoS2), with the presence of compositing agents, including Ti and Au, under different environmental exposures. We report a new ReaxFF reactive force field parameter set for Mo, Ti, Au, O, S, and H interactions. We apply the force field in a series of molecular dynamics (MD) simulations to unravel the impact of the Ti dopant on the oxidation/hydrogenation behaviors of MoS2 surface. The simulation results reveal that, in the absence of Ti clusters, the MoS2 surface is ruptured and oxidized at elevated temperatures through a process of adsorption followed by dissociation of the O2 molecules on the MoS2 surface during the temperature ramp. When the MoS2 surface is exposed to H2O molecules, surface hydrogenation is most favored, followed by oxidation, then hydroxylation. The introduction of Ti clusters to the systems mitigates the oxidation/hydrogenation of MoS2 at a low or intermediate temperature by capturing the O2/H2O molecules and locking the O/H-related radicals inside the clusters. However, OH− and H3O+ are emitted from the Ti clusters in the H2O environment as temperature rises, and the accelerating hydrogenation of MoS2 is consequently observed at an ultra-high temperature. These findings indicate an important but complex role of Ti dopants in mitigating the oxidation and hydrogenation of MoS2 under different environmental exposures. The possible mechanisms of oxidation and hydrogenation revealed by MD simulations can give an insight to the design of oxidation resistant TMDs and can be useful to the optical, electronic, magnetic, catalytic, and energy harvesting industries.

Computational study of effect of radiation induced crosslinking on the properties of flattened carbon nanotubes

Flattened carbon nanotubes (flCNTs) are a primary component of many carbon nanotube (CNT) yarn and sheet materials, which are promising reinforcements for the next generation of ultra-strong composites for aerospace applications. Significant improvements in the performance of CNT materials can be realized with improvements in the load transfer between flCNTs, which are generally oriented at different angles with respect to each other. An intriguing approach to improving the load transfer is via irradiation-induced chemical crosslinking between adjacent flCNTs. The objective of this research is to use molecular dynamics (MD) simulations to predict the behavior of flCNT junctions with 0- and 90-degree orientations and varying levels of crosslinking. The results indicate that crosslinking improves the flCNT interfacial load transfer for both orientations, but degrades the flCNT tensile response. The primary toughening mechanism at the flCNT/flCNT interface is the formation of carbon chains that provide load transfer up to the point of total rupture. Based on these results, it is clear that irradiation-induced crosslinking is beneficial in CNT-based composite systems in which interfacial load transfer between flCNTs is of primary importance, even though individual flCNTs may lose some mechanical integrity with crosslinking.

Evaluating the performance of ReaxFF potentials for sp2 carbon systems (graphene, carbon nanotubes, fullerenes) and a new ReaxFF potential

We study the performance of eleven reactive force fields (ReaxFF), which can be used to study sp2 carbon systems. Among them a new hybrid ReaxFF is proposed combining two others and introducing two different types of C atoms. The advantages of that potential are discussed. We analyze the behavior of ReaxFFs with respect to 1) the structural and mechanical properties of graphene, its response to strain and phonon dispersion relation; 2) the energetics of (n, 0) and (n, n) carbon nanotubes (CNTs), their mechanical properties and response to strain up to fracture; 3) the energetics of the icosahedral C60 fullerene and the 40 C40 fullerene isomers. Seven of them provide not very realistic predictions for graphene, which made us focusing on the remaining, which provide reasonable results for 1) the structure, energy and phonon band structure of graphene, 2) the energetics of CNTs versus their diameter and 3) the energy of C60 and the trend of the energy of the C40 fullerene isomers versus their pentagon adjacencies, in accordance with density functional theory (DFT) calculations and/or experimental data. Moreover, the predicted fracture strain, ultimate tensile strength and strain values of CNTs are inside the range of experimental values, although overestimated with respect to DFT. However, they underestimate the Young's modulus, overestimate the Poisson's ratio of both graphene and CNTs and they display anomalous behavior of the stress - strain and Poisson's ratio - strain curves, whose origin needs further investigation.

Shape Entropy of a Reconfigurable Ising Surface

Disclinations in a 2D sheet create regions of Gaussian curvature whose inversion produces a reconfigurable surface with many distinct metastable shapes, as shown by molecular dynamics of a disclinated graphene monolayer. This material has a near-Gaussian "density of shapes" and an effectively antiferromagnetic interaction between adjacent cones. A∼10  nm patch has hundreds of distinct metastable shapes with tunable stability and topography on the size scale of biomolecules. As every conical disclination provides an Ising-like degree of freedom, we call this technique "Isigami."

2D-Topology-Seeded Graphitization for Highly Thermally Conductive Carbon Fibers

Highly thermally conductive carbon fibers (CFs) have become an important material to meet the increasing demand for efficient heat dissipation. To date, high thermal conductivity has been only achieved in specific pitch-based CFs with high crystallinity. However, obtaining high graphitic crystallinity and high thermal conductivity beyond pitch-CFs remains a grand challenge. Here, a 2D-topology-seeded graphitization method is presented to mediate the topological incompatibility in graphitization by seeding 2D graphene oxide (GO) sheets into the polyacrylonitrile (PAN) precursor. Strong mechanical strength and high thermal conductivity up to 850 W m^{- 1} K^-1 are simultaneously realized, which are one order of magnitude higher in conductivity than commercial PAN-based CFs. The self-oxidation and seeded graphitization effect generate large crystallite size and high orientation to far exceed those of conventional CFs. Topologically seeded graphitization, verified in experiments and simulations, allows conversion of the non-graphitizable into graphitizable materials by incorporating 2D seeds. This method extends the preparation of highly thermally conductive CFs, which has great potential for lightweight thermal-management materials.

Microstructure Evolution and Its Correlation with Performance in Nitrogen-Containing Porous Carbon Prepared by Polypyrrole Carbonization: Insights from Hybrid Calculations

The preparation of nitrogen-containing porous carbon (NCPC) materials by controlled carbonization is an exciting topic due to their high surface area and good conductivity for use in the fields of electrochemical energy storage and conversion. However, the poor controllability of amorphous porous carbon prepared by carbonization has always been a tough problem due to the unclear carbonation mechanism, which thus makes it hard to reveal the microstructure-performance relationship. To address this, here, we comprehensively employed reactive molecular dynamics (ReaxFF-MD) simulations and first-principles calculations, together with machine learning technologies, to clarify the carbonation process of polypyrrole, including the deprotonation and formation of pore structures with temperature, as well as the relationship between microstructure, conductance, and pore size. This work constructed ring expressions for PPy thermal conversion at the atomic level. It revealed the structural factors that determine the conductivity and pore size of carbonized products. More significantly, physically interpretable machine learning models were determined to quantitatively express structure factors and performance structure-activity relationships. Our study also confirmed that deprotonation preferentially occurred by desorbing the dihydrogen atom on nitrogen atoms during the carbonization of PPy. This theoretical work clearly reproduces the microstructure evolution of polypyrrole on an atomic scale that is hard to do via experimentation, thus paving a new way to the design and development of nitrogen-containing porous carbon materials with controllable morphology and performance.

Cost-effective carbon fiber precursor selections of polyacrylonitrile-derived blend polymers: carbonization chemistry and structural characterizations

Blending polyacrylonitrile (PAN) with plastic wastes and bio-based polymers provides a convenient and inexpensive method to realize cost-effective carbon fiber (CF) precursors. In this work, PAN-based blend precursors are investigated using ReaxFF reactive molecular dynamics simulations with respect to the formation of all-carbon rings, the evolutions of oxygen-containing and nitrogen-containing species, and the migration of carbon atoms to form turbostratic graphene. From these simulations, we identify that PAN/cellulose (CL) blend manifests the highest carbon yield and the most substantial all-carbon ring formation. This ReaxFF-based finding is confirmed by Raman and TEM experiments indicating high crystallinity for PAN/CL-derived blend CFs. We trace the pathway of gasification and carbonization of PAN/CL to elaborate the mechanism of the formation of all-carbon ring networks. We discover that the acetals of CL can catalyze the cyclization of the blend precursor, allowing for the search for CL derivatives or the other kinds of bio-based polymers with similar functionalities as alternative blends. In addition, we examine the structural characteristics using the carbon-carbon (C-C) radial distribution functions, C-C bond length distributions, and sp2 C atom ratios for the four representative precursors, i.e., PAN, oxidized PAN, PAN/nylon 6,6, and PAN/CL. Our simulation results show the most extensive all-carbon ring cluster and graphitic structure growths for PAN/CL. Therefore, we propose PAN/CL as a cost-effective alternative CF precursor, since (a) CL is naturally abundant and eco-friendly for production, (b) the blend precursor PAN/CL does not require oxidation treatment, (c) PAN/CL has a high carbon yield with substantial all-carbon ring formation, and (d) PAN/CL based CFs potentially provide a mechanical property enhancement.

Molecular Dynamics Modeling of Interfacial Interactions between Flattened Carbon Nanotubes and Amorphous Carbon: Implications for Ultra-Lightweight Composites

Flattened carbon nanotubes (flCNTs) naturally form in many carbon nanotube-based materials and can exhibit mechanical properties similar to round carbon nanotubes but with tighter packing and alignment. To facilitate the design, fabrication, and testing of flCNT-based composites for aerospace structures, computational modeling can be used to efficiently and accurately predict their performance as a function of processing parameters, such as reinforcement/matrix cross-linking. In this study, molecular dynamics modeling is used to predict the load transfer characteristics of the interface region between the flat region of flCNTs (i.e., bi-layer graphene) and amorphous carbon (AC) with various levels and locations of covalent bond cross-linking and AC mass density. The results of this study show that increasing the mass density of AC at the interface improves the load transfer capability of the interface. However, a much larger improvement is observed when cross-linking is added both to the flCNT-AC interface and between the flCNT sheets. With both types of cross-linking, substantial improvements in interfacial shear strength, transverse tension strength, and transverse tension toughness are predicted. The results of this study are important for optimizing the processing of flCNT/AC composites for demanding engineering applications.

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically heterogeneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp³ fraction, achieved via generating high-density CFs. In addition, the high sp³ fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

Machine Learning-Assisted Hybrid ReaxFF Simulations

We have developed a machine learning (ML)-assisted Hybrid ReaxFF simulation method ("Hybrid/Reax"), which alternates reactive and non-reactive molecular dynamics simulations with the assistance of ML models to simulate phenomena that require longer time scales and/or larger systems than are typically accessible to ReaxFF. Hybrid/Reax uses a specialized tracking tool during the reactive simulations to further accelerate chemical reactions. Non-reactive simulations are used to equilibrate the system after the reactive simulation stage. ML models are used between reactive and non-reactive stages to predict non-reactive force field parameters of the system based on the updated bond topology. Hybrid/Reax simulation cycles can be continued until the desired chemical reactions are observed. As a case study, this method was used to study the cross-linking of a polyethylene (PE) matrix analogue (decane) with the cross-linking agent dicumyl peroxide (DCP). We were able to run relatively long simulations [>20 million molecular dynamics (MD) steps] on a small test system (4660 atoms) to simulate cross-linking reactions of PE in the presence of DCP. Starting with 80 PE molecules, more than half of them cross-linked by the end of the Hybrid/Reax cycles on a single Xeon processor in under 48 h. This simulation would take approximately 1 month if run with pure ReaxFF MD on the same machine.

A pre-oxidation strategy to improve architecture stability and electrochemical performance of Na2MnPO4F particles-embedded carbon nanofibers

The rational design of an excellent architecture for active materials combined with carbon matrix is of particular importanceto obtain flexible electrode material with high electrochemical properties. Well-designed nanofibers possess unique 3D network structure, which can significantly improve the electron/ion transportation and supplies sufficient active sites for Li⁺/Na⁺ insertion. Electrospinning-carbonization technology is a popular strategy to prepare nanofibers with active material embedded in carbon. It is found that the architecture of nanofibers tended to be wrecked and destroyed during the carbonization process without pre-oxidation treatment. In this study, we prepared Na₂MnPO₄F particles embedded in carbon nanofibers (Na₂MnPO₄F/C) using PVP as carbon source and investigated the strengthen mechanism of pre-oxidation on their architecture. The experiment and simulation results demonstrate that, without pre-oxidation, the main chain of PVP is severely ruptured during the carbonization procedure, consequently leads to fractured architecture of Na₂MnPO₄F/C nanofibers. In contrast, with pre-oxidation treatment, a long-chain and heat-resistance structured carbon matrix formed, and Na₂MnPO₄F/C nanofibers with stable architecture and improved electrochemical performance can be achieved. This study demonstrates a promising guide to construct carbon based nanofiber electrodes with stable architecture and high electrochemical performance.