Hydroxyl radical (HO·) plays an important role in the initial pyrolysis of hydroxyl-containing polymers, such as phenolic resin (PR). In this study, the reaction mechanism between HO· and bisphenol-F (BPF) or tetra-methyl substituted BPFs, which were taken as the model molecules of PR, was studied with the density functional theory approach. The results based on the Fukui function and reduced density gradient function showed that, both the hydroxyls and the carbon atoms in the phenolic groups are the reactive sites for HO· attack. The hydroxyls are most likely to be attacked by HO· owing to the strong electrostatic potential around the hydroxyls and the low reaction barriers, especially for cis-o-o' type BPF. The strong pAp (CHAp and OHAp) interaction between the phenolic rings in BPF leads to decreased conjugative effect of the phenolic rings, which further lead to decreased addition barriers and reaction rate constant. K E Y W O R D S bisphenol-F, conformation, density functional theory, hydroxyl radical, reaction mechanism 1 | I N T R O D U C T I O N Bis(hydroxyphenyl)methanes, also named as bisphenol-Fs or BPFs, are a group of chemical compounds with two hydroxyphenyl functionalities bridged by one methylene group ( Figure 1). Owing to the structurally similarity to bisphenol-A, BPFs are desirable materials with many attractive applications paralleling bisphenol-A in use, such as the primary raw materials of epoxy, polycarbonate, or phenolic resins, but considerably lower toxicity than bisphenol-A. [1] However, the widespread application of BPFs also leads to considerable environmental problems. [2] A detailed investigation on the thermal decomposition mechanism of BPFs will be conducive to guiding the application and degradation of BPFs in reasonable ways. In particular, BPFs can also be regarded as the basic structural units or model compounds for phenolic resins (PRs), which are generally made from phenol and formaldehyde. [3] Hence, the elucidation of the thermal decomposition mechanism of BPF will be helpful for understanding the relationship of structural-thermal stability of PRs.The constitutional and conformational isomerism of BPFs lead to multiple structural characteristics of BPFs. There are three constitutional BPF isomers, 2,2 0 -bis(hydroxyphenyl)methane, 2,4 0 -bis(hydroxyphenyl)methane, and 4,4 0 -bis(hydroxyphenyl)methane (Figure 1), and this is also the primary reason for the rather complex crosslinked structure of PR. Wang et al.'s computational results [4] showed that there are more than 24 stable conformers of o-o' type BPF. Pillsbury et al. [5] stated that owing to the hydrogen bonds and OHAp interaction formation, the conformers with hydroxyls paralleled or perpendicular to each other are the dominant structures. Gruber et al. [6] proved that a p-p 0 type BPF conformer with the C 2 symmetry has the lowest energy, leading to a strong crystallinity of p-p 0 type BPF. In the highly crosslinked network of cured PR, most of the BPF-like structures are four methylene substituted BPFs, this fu...