Comparison of ReaxFF, DFTB, and DFT for Phenolic Pyrolysis. 1. Molecular Dynamics Simulations

Efficient Parameterization of Density Functional Tight-Binding for 5f-Elements: A Th-O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th-O system, which was subsequently applied to the study of ThO2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO2. Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

Simulations on coal water slurry gasification by molecular dynamics method with ReaxFF

CONTEXT

Coal water slurry (CWS) is a new type of liquid coal product with low pollution, which is mainly used in the chemical industry to produce syngas (CO + H2). It is of great significance to study the microscopic mechanism of CWS gasification reaction for improving the efficiency of coal gasification. In this paper, the method of molecular dynamics based on reaction force fields (ReaxFF-MD) is used to study the gasification process of CWS/O2 system at different temperatures. The results show that, in the range of 1600-2400 K, the macromolecular network structure of lignite is decomposed into a large number of small molecular structures and a small number of light tar free radical fragments, and the types and quantities of reaction products increased rapidly. At 2400-4000 K, the free radical fragments of light tar are further decomposed and reacted with gasification agents, but the types and quantities of reaction products have little change. At 3600 K, a full gasification reaction occurred in the system, and the content of syngas is the highest.

METHODS

The model was established and optimized by Materials Studio (MS) software. Based on ReaxFF-MD method, Lammps software was used to simulate the gasification process of CWS/O2 system, and the reaction force field files containing C, H, O, N, and S element were used. By calculating the activation energy of gasification reaction, the rationality of the model and calculation method was illustrated. The post-processing of the results was implemented using OVITO, ChemDraw software, and self-programmed Python scripts.

Quantum chemical study on the catalytic debromination mechanism of brominated epoxy resins

The study employed Density Functional Theory (DFT) to investigate the catalytic debromination mechanism of brominated epoxy resins (BERs) by iron (Fe) and copper (Cu) catalysts. By introducing electric field (EF), intramolecular electron transfer and polarization effects on BERs debromination were explored and experimentally validated. Results indicated that the bond dissociation energy (BDE) of the C-Br bond was 312.27 kJ/mol without catalysis, while with Fe, Cu, and EF, it was 114.47 kJ/mol, 94.85 kJ/mol, and 292.59 kJ/mol, respectively, enhancing reactivity. EF parallel to the C-Br bond and oriented toward the C atom, altered electrostatic potential and dipole moment around C-Br bond, leading to 68.60% and 50.19% increment in electronic contribution difference and molecule polarity, respectively, thereby reducing the C-Br BDE. Fe and Cu facilitated electron transfers with BERs, inducing reactions between their negative electrostatic potentials and Br's positive potential, changing electron sharing, resulting in 19.87% and 12.11% increase in polarity, respectively, and further BDE reduction. Structural modifications by the EF and catalysts also intensified van der Waals forces with bromine atoms and decreased spatial hindrance, collectively making C-Br bond breakage easier. Experiments revealed the EF enhanced BERs' debromination efficiency but hindered Fe/Cu's catalysis at lower temperatures.

The evolution of atomistic structure and mechanical property of coke in the gasification process with CO2 and H2O at different temperatures: A ReaxFF molecular dynamics study

CONTEXT

An atomistic coke carbon model was constructed to simulate the structural evolution in the gasification and stretching process. The coke model was placed in a box with different CO2/H2O content to investigate the evolution of the atomistic structure of coke during the gasification. It was found that different atmospheric concentrations had different effects on the structure and reaction sites of the coke model. The CO2 molecules tended to dissolve on the surface of coke and disrupt its surface structure, while H2O molecules were more likely to enter the coke model to disrupt the internal structure. For tensile simulation, it was found that CO2 and H2O had different effects on the tensile resistance of the coke model. Controlling the composition content of the reaction gas can effectively influence the tensile strength of the coke model. By revealing the behavior of coke model at the micro scale, it provides a theoretical basis for the industrial coke application process.

METHODS

Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) was used to conduct the molecular dynamics using the reactive force field (ReaxFF). The atomistic model of coke carbon was constructed using the well-known annealing and quenching method, and its composition is determined according to the element analysis of industrial coke. The structural evolution in the gasification with CO2/H2O and the stretching process were analyzed in detail. Molecular dynamics simulations with reactive force field (ReaxFF-MD) were used to simulate the coke dissolution reaction under CO2/H2O atmosphere and the coke stretching process. The atmosphere ratio was modified to investigate the changes in coke structure under different atmosphere conditions. The Packmol software was used to place gas and coke models into the same box. During the reaction process, the Ovito software was used to perform corresponding visualization analysis on the changes in the atomic structure of coke.

Research on the pyrolysis characteristics and mechanisms of waste printed circuit boards at fast and slow heating rates

The pyrolysis treatment of waste printed circuit boards (WPCBs) shows great potential for sustainable treatment and hazard reduction. In this work, based on thermogravimetry (TG), pyrolysis-gas chromatography/mass spectrometry (Py-GC/MS), and density functional theory (DFT), the thermal weight loss, product distribution, and kinetics of WPCBs pyrolysis were studied by single-step and multi-step pyrolysis at fast (600 °C/min) and slow (10 °C/min) heating rates. The heating rates of TG and Py-GC/MS were the same for each group of experiments. In addition, the bond dissociation energy (BDE) of WPCBs polymer monomers was calculated by DFT method. Compared with slow pyrolysis, the final weight loss of fast pyrolysis is reduced by 0.76 wt%. The kinetic analysis indicates that the activation energies of main pyrolysis stages range from 98.29 kJ/mol to 177.59 kJ/mol. The volatile products of fast pyrolysis are mainly phenols and aromatics. With the increase of multi-step pyrolysis temperature, the order of the escaping volatiles is phenols, hydrocarbyl phenols, aromatics, and benzene (or diphenyl phenol). The pyrolysis residue of WPCBs may contains phenolics and polymers. Based on the free radical reactions, the mechanism and reaction pathways of WPCBs pyrolysis were deduced by the DFT. Moreover, a large amount of benzene is produced by pyrolysis, and its formation mechanism was elaborated.

Structural stability and initial decomposition mechanisms of BTF crystal induced by vacancy defects: a computational study

Density functional tight binding (DFTB) and DFTB-based molecular dynamics (DFTB-MD) were used to study the effects of vacancy defects on the structure, stability, and initial decomposition mechanisms of condensed phase...

Structure, intermolecular interactions, and dynamic properties of NTO crystals with impurity defects: a computational study

Frontier orbitals distribute in the position of impurity molecules, whose adjacent NTO molecules begin to decompose first.

Multiscale Simulation on Product Distribution from Pyrolysis of Styrene-Butadiene Rubber

Pyrolysis of styrene-butadiene rubber receives renewed attention due to its application in tackling the waste tire disposal problem while allowing energy recovery. The density functional theory calculation (DFT) and ReaxFF molecular dynamics simulation (MD) are adopted to study the pyrolysis process with the variation of temperature and pressure. The bond dissociation energies of intramonomer and intermonomer bonds in trimers with different linking methods are calculated by DFT, where the bond with low energy tends to break during the pyrolysis process. The following MD simulation shows the pyrolysis product distribution of chain segments in styrene-butadiene rubber, where bond breaking positions in MD agree well with corresponding results in DFT and experiment. The next nearest neighbor bonds (single bonds) connected with double bond or benzene usually have lower dissociation energies than other single bonds and prone to break during the pyrolysis process. And thus, the intermonomer bonds tend to break at relatively low temperatures (around 650 K in experiment) prior to intramonomer bonds, which result in the emergence of monomers. With the temperature increase, intramonomer bonds are broken and thus large fragments are further pyrolyzed into small ones (e.g., C₂ and C). Besides, the pressure strongly influences the product distribution, where high pressures promote the occurrence of secondary reactions.

Effect of Ambient Chemistry on Friction at the Basal Plane of Graphite

Graphite is widely used as a solid lubricant due to its layered structure, which enables ultralow friction. However, the lubricity of graphite is affected by ambient conditions and previous studies have shown a sharp contrast between frictional behavior in vacuum or dry environments compared to humid air. Here, we studied the effect of organic gaseous species in the environment, specifically comparing the adsorption of phenol and pentanol vapor. Atomic force microscopy experiments and reactive molecular dynamics simulations showed that friction was larger with phenol than with pentanol. The simulation results were analyzed to test multiple hypotheses to explain the friction difference, and it was found that mechanically driven chemical bonding between the tip and phenol molecules plays a critical role. Bonding increases the number of phenol molecules in the contact, which increases the adhesion as well as the number of atoms in registry with the topmost graphene layer acting as a pinning site to resist sliding. The findings of this research provide insight into how the chemistry of the operating environment can affect the frictional behavior of graphite and layered materials more generally.

Generating Converged Accurate Free Energy Surfaces for Chemical Reactions with a Force-Matched Semiempirical Model

We demonstrate the capability of creating robust density functional tight binding (DFTB) models for chemical reactivity in prebiotic mixtures through force matching to short time scale quantum free energy estimates. Molecular dynamics using density functional theory (DFT) is a highly accurate approach to generate free energy surfaces for chemical reactions, but the extreme computational cost often limits the time scales and range of thermodynamic states that can feasibly be studied. In contrast, DFTB is a semiempirical quantum method that affords up to a thousandfold reduction in cost and can recover DFT-level accuracy. Here, we show that a force-matched DFTB model for aqueous glycine condensation reactions yields free energy surfaces that are consistent with experimental observations of reaction energetics. Convergence analysis reveals that multiple nanoseconds of combined trajectory are needed to reach a steady-fluctuating free energy estimate for glycine condensation. Predictive accuracy of force-matched DFTB is demonstrated by direct comparison to DFT, with the two approaches yielding surfaces with large regions that differ by only a few kcal mol^-1.

Direct observation of realistic-temperature fuel combustion mechanisms in atomistic simulations

Advanced accelerated molecular dynamics simulations provide a detailed atomic-level picture of combustion at realistic temperatures and pressures.

Atomistic simulations can in principle provide an unbiased description of all mechanisms, intermediates, and products of complex chemical processes. However, due to the severe time scale limitation of conventional simulation techniques, unrealistically high simulation temperatures are usually applied, which are a poor approximation of most practically relevant low-temperature applications. In this work, we demonstrate the direct observation at the atomic scale of the pyrolysis and oxidation of n-dodecane at temperatures as low as 700 K through the use of a novel simulation technique, collective variable-driven hyperdynamics (CVHD). A simulated timescale of up to 39 seconds is reached. Product compositions and dominant mechanisms are found to be strongly temperature-dependent, and are consistent with experiments and kinetic models. These simulations provide a first atomic-level look at the full dynamics of the complicated fuel combustion process at industrially relevant temperatures and time scales, unattainable by conventional molecular dynamics simulations.

Thermal stability of phenolic resin: new insights based on bond dissociation energy and reactivity of functional groups

Based on bisphenol-F-like model molecules, the bond dissociation energies and Fukui function were calculated to interpret the relationship between the atomistic structure and thermal properties of the phenolic resin.

Quantum Chemical Molecular Dynamics Simulations of 1,3-Dichloropropene Combustion

Oxidative decomposition of 1,3-dichloropropene was investigated using quantum chemical molecular dynamics (QM/MD) at 1500 and 3000 K. Thermal oxidation of 1,3-dichloropropene was initiated by (1) abstraction of allylic H/Cl by O2 and (2) intra-annular C-Cl bond scission and elimination of allylic Cl. A kinetic analysis shows that (2) is the more dominant initiation pathway, in agreement with QM/MD results. These QM/MD simulations reveal new routes to the formation of major products (H2O, CO, HCl, CO2), which are propagated primarily by the chloroperoxy (ClO2), OH, and 1,3-dichloropropene derived radicals. In particular, intra-annular C-C/C-H bond dissociation reactions of intermediate aldehydes/ketones are shown to play a dominant role in the formation of CO and CO2. Our simulations demonstrate that both combustion temperature and radical concentration can influence the product yield, however not the combustion mechanism.

Investigation on the Gas-phase decomposition of trichlorfon by GC-MS and theoretical calculation

The gas phase pyrolysis of trichlorfon was investigated by the on-line gas chromatography – mass spectrometry (GC-MS) pyrolysis and theoretical calculations. Two reaction channels were proposed in the pyrolytic reaction, by analyzing the detected pyrolytic products in the total ion chromatography, including 2,2,2-trichloroacetaldehyde, dimethyl phosphite, and dichlorvos. Theoretical calculations showed that there is an intramolecular hydrogen bond between the hydroxyl group and the phosphate O atom in trichlorfon, through which the hydroxyl H atom can be easily transferred to phosphate O atom to trigger two pyrolytic channels. In path-a, migration of H atom results in direct decomposition of trichlorfon to give 2,2,2-trichloroacetaldehyde and dimethyl phosphite in one step. In path-b, migration of H atom in trichlorfon is combined with formation of the O-P bond to give an intermediate, followed by HCl elimination to afford dichlorvos. Path-a is kinetically more favorable than path-b, which is consistent with the GC-MS results.

Dipropargyl ether bisphenol A based boron-containing polymer: synthesis, characterization and molecular dynamics simulations of the resulting pyrolysis and carbonization

The time evolution of major pyrolysis products including small-molecule species of a dipropargyl ether bisphenol A based novel boron-containing polymer was examined via ReaxFF-MD simulation (Color code: C, grey; O, red; H, white; B, yellow).

Theoretical investigation of atomic oxygen erosion mechanisms of 1,3-didecyl cyclopentane, 1,3-dioctyldodecyl cyclopentane and alkylated cyclopentane

Atomic oxygen erosion mechanisms of 1,3-didecyl cyclopentane, 1,3-dioctyldodecyl cyclopentane and alkylated cyclopentane.