Poly (O-Aminophenol) Produced by Plasma Polymerization Has IR Spectrum Consistent with a Mixture of Quinoid & Keto Structures

A vibrational analysis of various poly(o-aminophenol) structures has been undertaken using first principles methods. It is shown that a mixture of quinoid and keto forms of poly(o-aminophenol) gives rise to a simulated spectrum that replicates the experimental infrared spectra of plasma-produced poly(o-aminophenol) better than either the quinoid or keto poly(o-aminophenol) spectra alone. An unassigned peak in the spectrum is attributed to hydrogen bonding to the silica substrate.


Introduction
Conductive polymers, such as polyaniline (PAN) and its derivatives, are widely studied due to their promising electric and optical properties, tunable configurations, and relatively low cost. PAN is commonly utilized in energy storage devices, such as rechargeable batteries, and chemical sensor applications [1,2]. Aminophenols are derivatives of anilines that benefit from containing two reaction points, amino (-NH 2 ) and hydroxyl (-OH) groups, that can be oxidized. Thus, aminophenols yield electrochemical properties comparable to those of both anilines and phenols [3]. Specifically, poly(o-aminophenol) (PoAP) finds use in solar cell technologies [4] and as a chemical sensor for glucose [5,6] and uric acid [7]. It has recently been demonstrated by Chen et al. [8] that PoAP can be easily synthesized by exposing the ortho-aminophenol (oAP) monomer to dielectric barrier discharge (DBD) plasma. Their process produces polymerized o-aminophenol, as confirmed by Fourier Transform Infrared (FTIR) and ultraviolet-visible (UV-Vis) spectroscopy [8].
It is commonly known that the structure of PoAP can exhibit a ladder-like or openchain form [9]. Upon oxidation, the quinoid [10] and keto [11,12] forms are the predominant structures of PoAP. These non-oxidated and oxidated structures are all represented in Figure 1.
Plasma treatment is an oxidative process, and thus it is reasonable to deduce that DBD plasma-produced PoAP would not feature un-oxidized ladder-like and open-chain structures (See Figure 1), as is the conclusion of Chen et al. [8]. The researchers identified their PoAP as having a quinoid structure, based on FTIR spectra featuring C=O stretching, C-O-C stretching, C=N stretching, and C-N stretching, while lacking N-H and O-H stretching frequencies. Plasma treatment is also commonly understood to be a highly reactive process, however, due to the ionized nature of the plasma gas. In concert with aminophenol's reactive amino and hydroxyl groups, we propose that DBD plasmaproduced PoAP may exhibit both oxidized quinoid and keto structures. It is not uncommon for more than one polymeric structure to be formed due to the multiple reaction points on the monomer. The works of Al-Hossainy et al. [4] and Bicak et al. [5] both report synthesizing PoAP structures containing both quinoid and benzoid rings after oxidation via ferric chloride/hydrochloric acid and copper bromide, respectively. As previously mentioned, it is reasonable to conclude that the PoAP product synthesized by Chen et al. [8] would Plasma 2022, 5 197 feature oxidized structures, though the relative abundance of the quinoid and keto forms (See Figure 1) remains unresolved. In this paper, we report the application of theoretical modeling of the infrared (IR) spectra of PoAP to explore the relative abundance of the quinoid and keto forms.
PoAP may exhibit both oxidized quinoid and keto structures. It is not uncommon for more than one polymeric structure to be formed due to the multiple reaction points on the monomer. The works of Al-Hossainy et al. [4] and Bicak et al. [5] both report synthesizing PoAP structures containing both quinoid and benzoid rings after oxidation via ferric chloride/hydrochloric acid and copper bromide, respectively. As previously mentioned, it is reasonable to conclude that the PoAP product synthesized by Chen et al. [8] would feature oxidized structures, though the relative abundance of the quinoid and keto forms (See Figure 1) remains unresolved. In this paper, we report the application of theoretical modeling of the infrared (IR) spectra of PoAP to explore the relative abundance of the quinoid and keto forms.

Materials and Methods
Electronic structure calculations and geometric optimization by energy minimization were carried out on the monomer and oligomers of oAP using density functional theory (DFT) with the hybrid B3LYP [13,14] functional at the 6-31G(d) [15,16] level of theory using the General Atomic and Molecular Electronic Structure System (GAMESS) suite of codes [17]. The threshold for geometric convergence was taken to be 0.0001 Hartree/Bohr. The optimized geometries were subsequently subjected to analytic Hessian calculations. All vibrational frequencies were scaled by a factor of 0.9614, which has been recommended by Scott and Radom [18] as optimal for the DFT/B3LYP/6-31G(d) level of theory. The vibrational modes were visualized & animated using Facio graphical user interface (GUI) software [19], to identify those modes corresponding to each prominent functional group (C=O, C=N stretching, etc.). The GUI was also used to generate theoretically predicted spectra from (frequency, intensity) data obtained from the electronic structure calculations. The experimental IR spectrum of Chen et al. [8], as well as the spectra predicted with the Facio [19] software, were digitized using WebPlotDigitizer [20] to facilitate a comparison.

Procedure
First, monomer structures for quinoid and keto o-aminophenol were optimized. From these optimized monomers, initial structural models of quinoid and keto trimers, (taken here to be representative of the full polymer) were constructed and subjected to structural optimization. Fully optimized structures for quinoid and keto trimers are shown in Figure 2A,B, respectively. Subsequently, harmonic normal mode vibrational frequencies were computed. Earlier computational work on oligomers of 3-phenyl-1-ureidonitrile [21] showed C=N stretching frequencies to be well converged with respect to the degree of oligomerization in the trimer, supporting the choice of trimer to represent the full polymer. Convergence tests on the present keto and quinoid o-aminophenol systems confirmed this choice. It was found that peak intensities for the dominant peaks used in

Materials and Methods
Electronic structure calculations and geometric optimization by energy minimization were carried out on the monomer and oligomers of oAP using density functional theory (DFT) with the hybrid B3LYP [13,14] functional at the 6-31G(d) [15,16] level of theory using the General Atomic and Molecular Electronic Structure System (GAMESS) suite of codes [17]. The threshold for geometric convergence was taken to be 0.0001 Hartree/Bohr. The optimized geometries were subsequently subjected to analytic Hessian calculations. All vibrational frequencies were scaled by a factor of 0.9614, which has been recommended by Scott and Radom [18] as optimal for the DFT/B3LYP/6-31G(d) level of theory. The vibrational modes were visualized & animated using Facio graphical user interface (GUI) software [19], to identify those modes corresponding to each prominent functional group (C=O, C=N stretching, etc.). The GUI was also used to generate theoretically predicted spectra from (frequency, intensity) data obtained from the electronic structure calculations. The experimental IR spectrum of Chen et al. [8], as well as the spectra predicted with the Facio [19] software, were digitized using WebPlotDigitizer [20] to facilitate a comparison.

Procedure
First, monomer structures for quinoid and keto o-aminophenol were optimized. From these optimized monomers, initial structural models of quinoid and keto trimers, (taken here to be representative of the full polymer) were constructed and subjected to structural optimization. Fully optimized structures for quinoid and keto trimers are shown in Figure 2A,B, respectively. Subsequently, harmonic normal mode vibrational frequencies were computed. Earlier computational work on oligomers of 3-phenyl-1-ureidonitrile [21] showed C=N stretching frequencies to be well converged with respect to the degree of oligomerization in the trimer, supporting the choice of trimer to represent the full polymer. Convergence tests on the present keto and quinoid o-aminophenol systems confirmed this choice. It was found that peak intensities for the dominant peaks used in the fitting are well converged at the trimer, with a mean unsigned difference between the trimer and tetramer intensities of 12%.
Plasma 2022, 5, FOR PEER REVIEW 3 the fitting are well converged at the trimer, with a mean unsigned difference between the trimer and tetramer intensities of 12%. Once predicted vibrational frequencies were obtained for the quinoid and keto PoAP trimer structures, which are shown in Figure 3B,C, respectively, a GUI software [19] was used to animate each vibrational mode. Thus, the C=O stretching, C=N stretching, C=C stretching, C-N stretching of aromatic rings, and C-O-C stretching modes were identified by visual inspection. A simulated spectrum for the polymer, shown in Figure 3D, was generated by finding the weighted average of the predicted quinoid and keto spectra that produces intensities for the dominant peaks in the best agreement with the experimental spectrum of Chen et al. [8]. This procedure consists of two steps: digitizing and optimization.

Digitizing
First, the experimental PoAP spectra of Chen et al. [8] were reconstructed based on a digitization of the published spectrum, as shown in Figure 3A. This was done by uploading an image of the published spectrum to WebPlotDigitizer [20] and collecting an intensity data point at each wavenumber from 500-3500 cm −1 . The major peaks were identified as those ca. = 1639, 1601, 1461, 1281, 1230, 1110, and 1050 cm −1 in the experimental spectrum. All spectral intensities were converted from percent transmittance to absorbance to facilitate intensity comparisons. The theoretical keto and quinoid trimer spectra were also digitized to facilitate a direct comparison between the theoretical intensities and the experimental PoAP spectral intensities.

Optimization
Next, a linear combination of theoretical intensities was generated to find the best match to the experimental spectrum. A simulated spectrum was produced by taking a weighted average of the theoretical intensity ( ) and ( ) values. By conservation of mass, the weights on the quinoid and keto forms must sum to 1, and thus, the weighted average of the theoretical intensities can be rewritten as: in which is a normalization factor applied to all theoretical frequencies to facilitate comparison with the experimental intensities. The relative abundances of the quinoid and keto forms are represented by and (1 − ), respectively. The optimal weighting factor was determined by minimizing the sum of squares of differences ( ) between the ( ) and experimental ( ) intensities: Once predicted vibrational frequencies were obtained for the quinoid and keto PoAP trimer structures, which are shown in Figure 3B,C, respectively, a GUI software [19] was used to animate each vibrational mode. Thus, the C=O stretching, C=N stretching, C=C stretching, C-N stretching of aromatic rings, and C-O-C stretching modes were identified by visual inspection. A simulated spectrum for the polymer, shown in Figure 3D, was generated by finding the weighted average of the predicted quinoid and keto spectra that produces intensities for the dominant peaks in the best agreement with the experimental spectrum of Chen et al. [8]. This procedure consists of two steps: digitizing and optimization.

Digitizing
First, the experimental PoAP spectra of Chen et al. [8] were reconstructed based on a digitization of the published spectrum, as shown in Figure 3A. This was done by uploading an image of the published spectrum to WebPlotDigitizer [20] and collecting an intensity data point at each wavenumber from 500-3500 cm −1 . The major peaks were identified as those ca. ν = 1639, 1601, 1461, 1281, 1230, 1110, and 1050 cm −1 in the experimental spectrum. All spectral intensities were converted from percent transmittance to absorbance to facilitate intensity comparisons. The theoretical keto and quinoid trimer spectra were also digitized to facilitate a direct comparison between the theoretical intensities and the experimental PoAP spectral intensities.

Optimization
Next, a linear combination of theoretical intensities was generated to find the best match to the experimental spectrum. A simulated spectrum was produced by taking a weighted average of the theoretical intensity I quinoid theo(ν) and I keto theo(ν) values. By conservation of mass, the weights on the quinoid and keto forms must sum to 1, and thus, the weighted average of the theoretical intensities can be rewritten as: (1) in which A is a normalization factor applied to all theoretical frequencies to facilitate comparison with the experimental intensities. The relative abundances of the quinoid and keto forms are represented by c and (1 − c), respectively. The optimal weighting factor c was determined by minimizing the sum of squares of differences (ε) between the I sim(ν) and experimental I exp (ν) intensities: by varying factors and . Here, the sum runs over the three most prominent composite peaks ( = 1605, 1406, and 1179 cm −1 ) and their experimental counterparts ( = 1601, 1461, and 1281 cm −1 ), respectively.

Results and Discussion
Vibrational assignments for the major peaks in the experimental spectrum of DBDplasma-produced PoAP are collected in Table 1 together with the corresponding theoretically computed frequencies for the quinoid and keto trimers and analogous experimental values from the literature. The theoretical frequencies are in reasonable quantitative agreement with those reported by Chen et al. [8] and with the literature values. Though the quinoid spectra do reproduce some of the major peaks seen in the experimental spectrum, the intensities of the peaks are not consistent with the experiment. The DBD plasmaproduced PoAP spectrum exhibits the most intense peak at 1281 cm −1 , corresponding to C-N stretching of aromatic rings [8]. The theoretical quinoid spectrum exhibits an absorbance at 1289 cm −1 , but its intensity (relative to the other peaks in the theoretical spectrum) is much lower than that seen in the experimental spectrum. The theoretical keto structure, however, does exhibit an intense peak at 1298 cm −1 . Because the keto structure exhibits a C-N aromatic ring stretching peak that is more consistent with the experimental data, it is reasonable to postulate that a linear combination of quinoid and keto frequencies could produce a simulated spectrum in better agreement with the experimental data than the quinoid spectra alone.

Simulated Spectrum
The predicted spectra for quinoid and keto trimers are shown in Figure 3B,C, respectively. It was found that c = 0.27 generated the lowest sum-of-squares of differences expressed in Equation (3), suggesting that the experimental spectrum represented a mixture of 27% quinoid and 73% keto PoAP. The simulated spectrum is shown in Figure 3D. The intensities of the peaks in the simulated spectrum arising from a mixture of forms are qualitatively a better match to the experimental spectrum than either the quinoid or keto spectra alone.

Packing Efficiency
One factor that introduces some uncertainty into the estimated composition is the packing efficiency. The theoretical intensities may be assumed to be proportional to the absorbance of a single molecule, but because the monomers differ in size, one can reasonably expect a higher number density of the smaller species in a physical sample. Two approaches to estimating this factor produce nearly identical results. The first approach is to compare molecular volumes [24]. The keto monomer has an estimated volume of V k = 270 Å 3  is consistent with the observation that many organic polymers have very similar material densities, typically in the range of [0.9-1.0] g/cc.) Given that density is mass per unit volume, equating densities written in terms of molecular weight (M) and molar quantity (n) gives, n k M k V = n q M q V , which may be rearranged into, n q n k = M k M q = 317 236 = 1.34. To investigate the possible effect of packing efficiency on the peak fitting analysis, a composite spectrum was generated from a linear combination of keto and quinoid intensities, in which the quinoid intensities were multiplied by the factor 1.34 before optimization of the fitting parameters. Equation (1), when modified to include this concentration factor, may be written as: Optimizing the c and A factors in Equation (3) to minimize the sum-of-squares of difference as expressed in Equation (3) produced c = 0.21. Thus, after applying the packing efficiency factor to the theoretical quinoid intensities, the experimental spectrum can be represented by a composite spectrum of 21% quinoid and 79% keto structures. This result is 6% different than the composite spectrum generated without the packing efficiency factor, which is one possible estimate of the uncertainty in the peak fitting analysis.

Unidentified Peak
One feature in the experimental spectrum of PoAP that is not captured by the simulated spectrum is the broad peak in the range of 2600-3300 cm −1 . The work of Chen et al. [8] did not report an assignment for this peak, as it was not featured in their polymerization identification process. It should also be noted that this broad peak is seen in their FTIR spectra for the monomer oAP as well as PoAP, indicating that it does not arise from a structure created by polymerization [8].
To create the most accurate linear combination of quinoid and keto structures, several attempts were made to identify the origin of this broad peak and remove it (and other associated contaminant peaks) from the experimental spectra. While C-H stretching is known to produce absorbance in the~3000 cm −1 region, the broad peak was not reproduced by the vibrational frequency calculations for either the quinoid or keto structure, nor does it appear in the computed spectra of the monomer (not shown here). The work of Bicak et al. also described a broad peak at~3200 cm −1 in their FTIR spectroscopic studies of PoAP, reportedly due to the hydroxyl stretching of free -OH groups, the presence of which they confirmed by proton NMR [5]. Thus, the first possible assignment considered for the broad peak was hydroxyl stretching from the ethanol solvent used in the experiment of Chen et al. [8]. The removal of the hypothesized ethanol absorbance from the experimental spectrum to create a new simulated spectrum consisted of three steps: digitization, subtraction, and optimization.
First, an FTIR spectrum of ethanol was obtained from the National Institute of Standards and Technology (NIST) [25] and digitized (Figure 4), as described in the Procedure section. For consistency, an ethanol spectrum having the same resolution (4 cm −1 ) as the experimental spectra of Chen et al. [8] was selected. Since the intensities of the ethanol peaks were not on the same scale as those in the experimental PoAP spectrum, the ethanol spectrum had to be normalized to facilitate an intensity comparison. This was done by identifying the most intense ethanol peak (ν =2974 cm −1 , A= 1.92) and calculating a scaling factor (S) (See Equation (4)) such that this peak would be of equal intensity to the (assumed corresponding) peak at 2974 cm −1 in Chen et al.'s spectrum of PoAP (S = 0.01859).
factor ( ) (See Equation (4)) such that this peak would be of equal intensity to the (as-  . Experimental spectra of ethanol obtained from NIST database [25]. The ethanol absorbances were all multiplied by this scaling factor , generating a normalized ethanol spectrum. The normalized ethanol intensities at each wavenumber from 500-3500 cm −1 were then subtracted from those of the experimental PoAP spectrum, generating a "decontaminated" experimental PoAP spectrum with the (assumed) contribution from ethanol removed.
After the decontaminated experimental PoAP spectrum with intensities ( ) was generated, the peak fitting analysis was completed as described in the Procedure section. Equation (2) was utilized to generate a weighted sum of the quinoid and keto spectra that would best reproduce the intensities at 1601, 1461, and 1281 cm −1 in the decontaminated experimental spectrum. The substitution of Equation (1) into Equation (2) and optimization of the and parameters in the resulting expression found that = 0.30 minimized the sum-of-squares of difference between ( ) and ( ) , suggesting that the scaled experimental spectrum with ethanol removed can be represented by a composite spectrum featuring 30% quinoid and 70% keto forms. This change in of 3% from the original fitting is smaller in magnitude than the change introduced by estimating the effect of the packing efficiency. Therefore, it is reasonable to conclude that removal of ethanol contributions does not change the predicted percentages of quinoid and keto structures at a significant level.
Further analysis suggests that there are several additional reasons to conclude that ethanol is unlikely to be responsible for the unidentified peak. First, ethanol exhibits a strong, sharp absorbance of ~1050 cm −1 , which does not appear in the spectrum of the polymer. Consequently, after subtraction of the (postulated) contribution from ethanol, the experimental spectrum exhibits a strong negative absorbance, which is unphysical. Second, the OH stretching band of ethanol is not well matched in shape or position to the broad unidentified peak in the experimental spectrum of the polymer. Third, ethanol would most likely be vaporized upon plasma treatment and thus would not be visible in the experimental spectra.
Another possible explanation for the broad unidentified peak is based on the observation that the shape of the broad peak is reminiscent of those resulting from the hydroxyl stretching of carboxylic acid dimers, which produces strong and broad peaks in the 2300-3300 cm −1 region [26]. Carboxylic acids form dimers through hydrogen bonding, which effectively red-shifts the O-H stretching adsorption compared to that of the monomer The ethanol absorbances were all multiplied by this scaling factor S, generating a normalized ethanol spectrum. The normalized ethanol intensities at each wavenumber from 500-3500 cm −1 were then subtracted from those of the experimental PoAP spectrum, generating a "decontaminated" experimental PoAP spectrum with the (assumed) contribution from ethanol removed.
After the decontaminated experimental PoAP spectrum with intensities I scaled exp (ν) was generated, the peak fitting analysis was completed as described in the Procedure section. Equation (2) was utilized to generate a weighted sum of the quinoid and keto spectra that would best reproduce the intensities at 1601, 1461, and 1281 cm −1 in the decontaminated experimental spectrum. The substitution of Equation (1) into Equation (2) and optimization of the A and c parameters in the resulting expression found that c = 0.30 minimized the sum-of-squares of difference between I sim(ν) and I scaled exp (ν) , suggesting that the scaled experimental spectrum with ethanol removed can be represented by a composite spectrum featuring 30% quinoid and 70% keto forms. This change in c of 3% from the original fitting is smaller in magnitude than the change introduced by estimating the effect of the packing efficiency. Therefore, it is reasonable to conclude that removal of ethanol contributions does not change the predicted percentages of quinoid and keto structures at a significant level.
Further analysis suggests that there are several additional reasons to conclude that ethanol is unlikely to be responsible for the unidentified peak. First, ethanol exhibits a strong, sharp absorbance of~1050 cm −1 , which does not appear in the spectrum of the polymer. Consequently, after subtraction of the (postulated) contribution from ethanol, the experimental spectrum exhibits a strong negative absorbance, which is unphysical. Second, the OH stretching band of ethanol is not well matched in shape or position to the broad unidentified peak in the experimental spectrum of the polymer. Third, ethanol would most likely be vaporized upon plasma treatment and thus would not be visible in the experimental spectra.
Another possible explanation for the broad unidentified peak is based on the observation that the shape of the broad peak is reminiscent of those resulting from the hydroxyl stretching of carboxylic acid dimers, which produces strong and broad peaks in the 2300-3300 cm −1 region [26]. Carboxylic acids form dimers through hydrogen bonding, which effectively red-shifts the O-H stretching adsorption compared to that of the monomer [27,28]. Although there is no obvious source of carboxylic acids in the polymer sample, red-shifted O-H stretching vibrations could arise from H-bonding interactions between the polymer (or monomer) and the SiO 2 surface, which was used as a support in the experimental synthesis of plasma-produced PoAP [8].
To explore this possibility, calculations were carried out to investigate the behavior of the O-H stretching frequency when the -OH moiety was engaged in hydrogen bonding to the silicon oxide surface. A cluster model of aminophenol hydrogen-bonded to a silicacontaining species is shown schematically in Figure 5.
red-shifted O-H stretching vibrations could arise from H-bonding interactions between the polymer (or monomer) and the SiO2 surface, which was used as a support in the experimental synthesis of plasma-produced PoAP [8].
To explore this possibility, calculations were carried out to investigate the behavior of the O-H stretching frequency when the -OH moiety was engaged in hydrogen bonding to the silicon oxide surface. A cluster model of aminophenol hydrogen-bonded to a silicacontaining species is shown schematically in Figure 5. First, a model of an o-aminophenol monomer was fully optimized, and the energy second derivatives were calculated to find its vibrational frequencies. A GUI software [19] was used to animate each vibrational mode, and the scaled [18] DFT/B3LYP/6-31G(d) frequency corresponding to O-H stretching was found to be 3862 cm −1 with an intensity of 3.6 (arb. units). The o-aminophenol monomer was then re-optimized with a disiloxane molecule in close proximity. The disiloxane molecule was oriented such that its oxygen atom was adjacent to the hydrogen on the oAP's OH functional group and the amino group was out of plane from the disiloxane's oxygen to attenuate potential hydrogen bonding between the oAP amino group and the disiloxane oxygen. The vibrational frequencies of the oAP and disiloxane molecule were calculated. An empirical dispersion correction [29] was specified to model hydrogen-bonding. The scaled O-H stretching frequency of AP-OH-silanol was found to be 3413 cm −1 with an intensity of 16.
The change in with the addition of disiloxane was −449 cm −1 , paired with a fourfold increase in intensity. Clearly, H-bonding to the heavy silicon atom effectively decreases/red-shifts the vibrational frequency of O-H stretching. Additionally, the increase in intensity with proximity to Si atoms would explain the relative strength of the broad unidentified peak in the experimental spectrum of Chen et al. [8]. These results can be extrapolated to suggest that the plasma-produced PoAP is exhibiting H-bonding with the SiO2 surface, effectively causing its O-H stretching modes to absorb strongly in the 2600-3300 cm −1 range. It seems most reasonable to conclude that the broad unidentified peak arises from H-bonding interactions with the SiO2 substrate. Under that assumption, no contaminant subtraction is necessary. Regardless, the predicted quinoid:keto ratio is ca. 27:73. This influence of the interaction between aminophenol and silica was simulated by calculating the vibrational frequencies for an aminophenol molecule and determining its O-H stretching frequency, then adding a disiloxane molecule near the aminophenol's hydroxyl group and re-calculating the aminophenol's O-H stretching absorption. First, a model of an o-aminophenol monomer was fully optimized, and the energy second derivatives were calculated to find its vibrational frequencies. A GUI software [19] was used to animate each vibrational mode, and the scaled [18] DFT/B3LYP/6-31G(d) frequency corresponding to O-H stretching was found to be 3862 cm −1 with an intensity of 3.6 (arb. units). The o-aminophenol monomer was then re-optimized with a disiloxane molecule in close proximity. The disiloxane molecule was oriented such that its oxygen atom was adjacent to the hydrogen on the oAP's OH functional group and the amino group was out of plane from the disiloxane's oxygen to attenuate potential hydrogen bonding between the oAP amino group and the disiloxane oxygen. The vibrational frequencies of the oAP and disiloxane molecule were calculated. An empirical dispersion correction [29] was specified to model hydrogen-bonding. The scaled O-H stretching frequency of AP-OH-silanol was found to be 3413 cm −1 with an intensity of 16.
The change in ν OH with the addition of disiloxane was −449 cm −1 , paired with a fourfold increase in intensity. Clearly, H-bonding to the heavy silicon atom effectively decreases/red-shifts the vibrational frequency of O-H stretching. Additionally, the increase in intensity with proximity to Si atoms would explain the relative strength of the broad unidentified peak in the experimental spectrum of Chen et al. [8]. These results can be extrapolated to suggest that the plasma-produced PoAP is exhibiting H-bonding with the SiO 2 surface, effectively causing its O-H stretching modes to absorb strongly in the 2600-3300 cm −1 range. It seems most reasonable to conclude that the broad unidentified peak arises from H-bonding interactions with the SiO 2 substrate. Under that assumption, no contaminant subtraction is necessary. Regardless, the predicted quinoid:keto ratio is ca. 27:73.

Conclusions
Optimized geometries of quinoid and keto PoAP trimers (taken to be representative of the polymer) were obtained using first principles methods. Harmonic vibrational frequencies were then calculated using analytic second derivatives. A least-squares fitting was used to find the weighting factors for the two forms that produced a simulated spectrum in best agreement with the experiment. It corresponds to 27% quinoid and 73% keto PoAP and is qualitatively more similar to the experimental spectrum of DBD plasma-produced PoAP than the spectra of either quinoid or keto structures alone. A second simulated spectrum was created after subtracting a (supposed) contribution of ethanol from the experimental spectra, but the experimental spectrum after subtraction exhibited a strong unphysical negative absorbance. This observation, coupled with the likely loss of ethanol during plasma treatment, leads us to dismiss ethanol as an appreciable contributor to the experimental spectrum. The vibrational analysis of the shift in OH stretching upon hydrogen bonding and literature reports of the shape of IR peaks corresponding to hydrogen-bond-shifted OH vibrations suggest that the most probable origin of a broad unidentified peak in the experimental spectrum of both the polymer and its unpolymerized monomer is the red-shifting of OH vibrations due to H-bonding interactions with the SiO 2 substrate.
Overall, the present vibrational frequency analysis reveals that the IR spectrum of poly(o-aminophenol) produced by DBD plasma treatment can be more accurately described by the presence of quinoid and keto structures than by the quinoid structure alone.