DFT and AIMD insights into heterogeneous dissociation of 2-chlorothiophenol on CuO(111) surface: Impact of H2O and OH.

Copper oxides are vital catalysts in facilitating the formation of polychlorinated thianthrenes/dibenzothiophenes (PCTA/DTs) through heterogeneous reactions in high-temperature industrial processes. Chlorothiophenols (CTPs) are the most crucial precursors for PCTA/DT formation. The initial step in this process is the metal-catalyzed production of chlorothiophenoxy radicals (CTPRs) from CTPs via dissociation reactions. This work combines density functional theory (DFT) calculations with ab initio molecular dynamics (AIMD) simulations to explore the formation mechanism of the adsorbed 2-CTPR from 2-CTP, with the assistance of CuO(111). Our study demonstrates that flat adsorption configurations of 2-CTP on the CuO(111) surface are more stable than vertical configurations. The CuO(111) surface acts as a strong catalyst, facilitating the dissociation of 2-CTP into the adsorbed 2-CTPR. Surface oxygen vacancies enhance the adsorption of 2-CTP on the CuO(111) surface, while moderately suppressing the dissociation of 2-CTP. More importantly, water molecules and surface hydroxyl groups actively promote the dissociation of 2-CTP. Specifically, water directly participates in the reaction through "water bridge", enabling a barrier-free process. This research provides molecular-level insights into the heterogeneous generation of dioxins with the catalysis of metal oxides in fly ash from static and dynamic aspects, providing novel approaches for reducing dioxin emissions and establishing dioxin control strategies.

Formation of Pre-PCTA/DT Intermediates from 2-Chlorothiophenol on Silica Clusters: A Quantum Mechanical Study.

Silica (SiO2), accounting for the main component of fly ash, plays a vital role in the heterogeneous formation of polychlorinated thianthrenes/dibenzothiophenes (PCTA/DTs) in high-temperature industrial processes. Silica clusters, as the basic units of silica, provide reasonable models to understand the general trends of complex surface reactions. Chlorothiophenols (CTPs) are the most crucial precursors for PCTA/DT formation. By employing density functional theory, this study examined the formation of 2-chlorothiophenolate from 2-CTP adsorbed on the dehydrated silica cluster ((SiO2)3) and the hydroxylated silica cluster ((SiO2)3O2H4). Additionally, this study investigated the formation of pre-PCTA/DTs, the crucial intermediates involved in PCTA/DT formation, from the coupling of two adsorbed 2-chlorothiophenolates via the Langmuir-Hinshelwood (L-H) mechanism and the coupling of adsorbed 2-chlorothiophenolate with gas-phase 2-CTP via the Eley-Rideal (E-R) mechanism on silica clusters. Moreover, the rate constants for the main elementary steps were calculated over the temperature range of 600-1200 K. Our study demonstrates that the 2-CTP is more likely to adsorb on the termination of the dehydrated silica cluster, which exhibits more effective catalysis in the formation of 2-chlorothiophenolate compared with the hydroxylated silica cluster. Moreover, the E-R mechanism mainly contributes to the formation of pre-PCTAs, whereas the L-H mechanism is prone to the formation of pre-PCDTs on dehydrated and hydroxylated silica clusters. Silica can act as a relatively mild catalyst in facilitating the heterogeneous formation of pre-PCTA/DTs from 2-CTP. This research provides new insights into the surface-mediated generation of PCTA/DTs, further providing theoretical foundations to reduce dioxin emission and establish dioxin control strategies.

Periodic DFT calculations for the heterogeneous formation of 2-chlorothiophenoxy radical from 2-chlorothiophenol on Cu(111) surface in fly ash.

Copper plays a crucial role in the heterogenous dissociation of chlorothiophenols (CTPs) to form chlorothiophenoxy radicals (CTPRs), which is the initial and critical step in the formation of polychlorinated thianthrenes/dibenzothiophenes (PCTA/DTs). Here, first-principles calculations were performed to investigate the activity of Cu(111) surface towards the formation of adsorbed 2-CTPR from 2-CTP. The interaction between 2-CTP and Cu(111) surface was explored to find stable adsorption configurations. Besides, the decomposition routes of 2-CTP on the Cu(111) surface were further explored. Moreover, the effects of water on the formation of absorbed 2-CTPR on the Cu(111) surface were examined. Our results demonstrate that the flat adsorption of 2-CTP on the surface with adsorption energy in the range of -33.21 kcal/mol to -28.37 kcal/mol is more stable than the vertical adsorption with adsorption energy ranging from -23.53 kcal/mol to -13.38 kcal/mol. The Cu(111) surface catalyzes the conversion of 2-CTP into the adsorbed 2-CTPR with a modest energy barrier of 9.46 kcal/mol. Furthermore, water molecules exhibit stronger catalytic activity in this process with a decreased energy barrier of 5.87 kcal/mol through "water bridge" and hydrogen bonding. Specifically, the water accepts the hydrogen atom from 2-CTP and donates another hydrogen to the surface via "water bridge". This research provides a molecular-level understanding of the heterogeneous formation of PCTA/DTs by fly ash, suggesting novel approaches for control strategy and legislation of dioxin analogues.

The effects of the gas-liquid interface and gas phase on Cl/ClO radical interaction with water molecules.

In the marine boundary layer (MBL), chlorine (Cl) and chlorine monoxide (ClO) are powerful oxidants with high concentrations. The gas-liquid interface is also ubiquitous in the MBL as a favorable site for atmospheric reactions. Understanding the role of water in Cl/ClO radical chemistry is essential for predicting their behavior in the atmosphere and developing effective strategies for mitigating their harmful effects. However, the research studies on the system of Cl/ClO radicals on the surface of water droplets are still insufficient. In previous studies, we have found unique results related to the hydroxyl radical at the interface using ab initio molecular dynamics (AIMD). In this work, we have used AIMD to investigate interactions between Cl/ClO radicals and water molecules at the gas-liquid interface. Radical mobility, radial distribution functions, coordination, and population analyses were conducted to investigate the surface preference, bonding pattern, and track Cl/ClO radicals in the water droplets. In addition, density functional theory (DFT) analysis was conducted to compare the results at the gas-liquid interface with those in the gas phase. We found that Cl/ClO radicals tend to remain near the gas-liquid interface in water droplet systems and outside of water clusters in gas phase systems. The ClO radical can form O*-H and Cl-O bonds with water molecules; however, neither the O*-O hemibond nor the Cl-H bond was detected in all systems. Different dominant structures were obtained for ClO in the interface and gas phase. The ClO radical can be bonded to one water molecule from its oxygen side, (H2O)0-Cl-O*-(H2O)1 at the interface, or to two water molecules from the chlorine and oxygen sides, (H2O)1-Cl-O*-(H2O)1 in the gas phase. Meanwhile, the Cl radical can only form a dominant structure like Cl*-(H2O)1 at the gas-liquid interface by making a Cl*-O hemibond. Providing a thorough explanation of the Cl/ClO radical behavior at the gas-liquid interface, this study will improve our understanding of the MBL's oxidizing capacity and pollution causes.

Theoretical insights into the degradation mechanisms, kinetics and eco-toxicity of oxcarbazepine initiated by OH radicals in aqueous environments.

As an anticonvulsant, oxcarbazepine (OXC) has attracted considerable attention for its potential threat to aquatic organisms. Density functional theory has been used to study the mechanisms and kinetics of OXC degradation initiated by OH radicals in aqueous environment. A total of fourteen OH-addition pathways were investigated, and the addition to the C8 position of the right benzene ring was the most vulnerable pathway, resulting in the intermediate IM8. The H-abstraction reactions initiated by OH radicals were also explored, where the extraction site of the methylene group (C14) on the seven-member carbon heterocyclic ring was found to be the optimal path. The calculations show that the total rate constant of OXC with OH radicals is 9.47 × 109 (mol/L)-1sec-1, and the half-life time is 7.32 s at 298 K with the [·OH] of 10-11 mol/L. Moreover, the branch ratio values revealed that OH-addition (89.58%) shows more advantageous than H-abstraction (10.42%). To further understand the potential eco-toxicity of OXC and its transformation products to aquatic organisms, acute toxicity and chronic toxicity were evaluated using ECOSAR software. The toxicity assessment revealed that most degradation products such as OXC-2OH, OXC-4OH, OXC-1O-1OOH, and OXC-1OH' are innoxious to fish and daphnia. Conversely, green algae are more sensitive to these compounds. This study can provide an extensive investigation into the degradation of OXC by OH radicals and enrich the understanding of the aquatic oxidation processes of pharmaceuticals and personal care products (PPCPs).

The Homogeneous Gas-Phase Formation Mechanism of PCNs from Cross-Condensation of Phenoxy Radical with 2-CPR and 3-CPR: A Theoretical Mechanistic and Kinetic Study.

Chlorophenols (CPs) and phenol are abundant in thermal and combustion procedures, such as stack gas production, industrial incinerators, metal reclamation, etc., which are key precursors for the formation of polychlorinated naphthalenes (PCNs). CPs and phenol can react with H or OH radicals to form chlorophenoxy radicals (CPRs) and phenoxy radical (PhR). The self-condensation of CPRs or cross-condensation of PhR with CPRs is the initial and most important step for PCN formation. In this work, detailed thermodynamic and kinetic calculations were carried out to investigate the PCN formation mechanisms from PhR with 2-CPR/3-CPR. Several energetically advantageous formation pathways were obtained. The rate constants of key elementary steps were calculated over 600~1200 K using the canonical variational transition-state theory (CVT) with the small curvature tunneling (SCT) contribution method. The mechanisms were compared with the experimental observations and our previous works on the PCN formation from the self-condensation of 2-CPRs/3-CPRs. This study shows that naphthalene and 1-monochlorinated naphthalene (1-MCN) are the main PCN products from the cross-condensation of PhR with 2-CPR, and naphthalene and 2-monochlorinated naphthalene (2-MCN) are the main PCN products from the cross-condensation of PhR with 3-CPR. Pathways terminated with Cl elimination are preferred over those terminated with H elimination. PCN formation from the cross-condensation of PhR with 3-CPR can occur much easier than that from the cross-condensation of PhR with 2-CPR. This study, along with the study of PCN formation from the self-condensation 2-CPRs/3-CPRs, can provide reasonable explanations for the experimental observations that the formation potential of naphthalene is larger than that of 1-MCN using 2-CP as a precursor, and an almost equal yield of 1-MCN and 2-MCN can be produced with 3-CP as a precursor.

The mobility and solvation structure of a hydroxyl radical in a water nanodroplet: a Born-Oppenheimer molecular dynamics study.

Hydroxyl radicals (OH*) play a crucial role in atmospheric chemistry and biological processes. In this study, Born-Oppenheimer molecular dynamics simulations are performed under ambient conditions for a hydroxyl radical in a water nanodroplet containing 191 water molecules. Density functional theory calculations are performed at the BLYP-D3 level with some test calculations at the B3LYP-D3 level. In two 150 ps trajectories, either with OH* initially located in the interior region or at the surface of the water nanodroplet, the OH* radical ends up in the subsurface layer of the nanodroplet, which is different from the "surface preference" predicted from previous empirical force field simulations. The solvation structure of OH* contains fluctuating hydrogen bonds, as well as a two-center three-electron hemibond in some cases. The mobility of OH* is enhanced by hydrogen transfer, which has a free energy barrier of ∼4.6 kcal mol-1. The results presented in this study deepen our understanding of the structure and dynamics of OH* in aqueous solutions, especially around the air-water interface.