Development of a novel ReaxFF reactive potential for organochloride molecules.

This article presents a new reactive potential in the ReaxFF formalism. It aims to include the chlorine element and opens up the fields of use of ReaxFF to the whole class of organochloride compounds including conjugated or aromatic groups. Numerous compounds in this family raise global awareness due to their environmental impact, and such a reactive potential will help investigate their degradation pathways. The new force field, named CHONCl-2022_weak, belongs to the aqueous branch. The force field parameters were fitted against high-level quantum chemistry calculations, including complete active space self-consistent field/NEVPT2 calculations and density functional theory calculations, and their accuracy was evaluated using a validation set. The root means square deviation against quantum mechanics energies is 0.38 eV (8.91 kcal mol-1). From a structural point of view, the root means square deviation is about 0.06 Å for the bond lengths, 11.86° for the angles, and 4.12° for the dihedral angles. With CHONCl-2022_weak new force field, we successfully investigated the regioselectivity for nucleophilic or electrophilic attacks on polychlorinated biphenyls, which are toxic and permanent pollutants. The rotation barriers along the bond linking the two benzene rings, which is crucial in the toxicity of these compounds, are well reproduced by CHONCl-2022_weak. Then, our new reactive potential is used to investigate the chlorobenzene reactivity in the presence of hydroxyl radicals in atmospheric condition or in aqueous solution. The reaction pathways computed with ReaxFF agree with the quantum mechanics results. We showed that, in the presence of dioxygen molecules, in atmospheric condition, the oxidation of chlorobenzene likely leads to the formation of highly oxygenated compounds after the abstraction of hydrogen radicals. In water, the addition of a hydroxyl radical leads to the formation of chlorophenol or phenol molecules, as already predicted from plasma-induced degradation experiments.

Gas-to-Particle Partitioning of Cyclohexene- and α-Pinene-Derived Highly Oxygenated Dimers Evaluated Using COSMOtherm.

Oxidized organic compounds are expected to contribute to secondary organic aerosol (SOA) if they have sufficiently low volatilities. We estimated saturation vapor pressures and activity coefficients (at infinite dilution in water and a model water-insoluble organic phase) of cyclohexene- and α-pinene-derived accretion products, "dimers", using the COSMOtherm19 program. We found that these two property estimates correlate with the number of hydrogen bond-donating functional groups and oxygen atoms in the compound. In contrast, when the number of H-bond donors is fixed, no clear differences are seen either between functional group types (e.g., OH or OOH as H-bond donors) or the formation mechanisms (e.g., gas-phase radical recombination vs liquid-phase closed-shell esterification). For the cyclohexene-derived dimers studied here, COSMOtherm19 predicts lower vapor pressures than the SIMPOL.1 group-contribution method in contrast to previous COSMOtherm estimates using older parameterizations and nonsystematic conformer sampling. The studied dimers can be classified as low, extremely low, or ultra-low-volatility organic compounds based on their estimated saturation mass concentrations. In the presence of aqueous and organic aerosol particles, all of the studied dimers are likely to partition into the particle phase and thereby contribute to SOA formation.