Structural Relevance of Intramolecular H-Bonding in Ortho-Hydroxyaryl Schiff Bases: The Case of 3-(5-bromo-2-hydroxybenzylideneamino) Phenol

A new Schiff base compound, 3-(5-bromo-2-hydroxybenzylideneamino)phenol (abbreviated as BHAP) was synthesized and characterized by 1H- and 13C- nuclear magnetic resonance and infrared spectroscopies. DFT/B3LYP/6-311++G(d,p) calculations were undertaken in order to explore the conformational space of both the E- and Z- geometrical isomers of the enol-imine and keto-amine tautomers of the compound. Optimized geometries and relative energies were obtained, and it was shown that the most stable species is the E-enol-imine form, which may exist in four low-energy intramolecularly hydrogen-bonded forms (I, II, V, and VI) that are almost isoenergetic. These conformers were concluded to exist in the gas phase equilibrium with nearly equal populations. On the other hand, the infrared spectra of the compound isolated in a cryogenic argon matrix (10 K) are compatible with the presence in the matrix of only two of these conformers (conformers II and V), while conformers I and VI convert to these ones by quantum mechanical tunneling through the barrier associated with the rotation of the OH phenolic group around the C–O bond. The matrix isolation infrared spectrum was then assigned and interpreted with help of the DFT(B3LYP)/6-311++G(d,p) calculated infrared spectra for conformers II and V. In addition, natural bond orbital (NBO) analysis was performed on the most stable conformer of the experimentally relevant isomeric form (E-enol-imino conformer V) to shed light on details of its electronic structure. This investigation stresses the fundamental structural relevance of the O–H···N intramolecular H-bond in o-hydroxyaryl Schiff base compounds.

The applications of this type of Schiff base are mostly determined by specific interactions involving their -CH=N− and phenolic o-hydroxyl substituent moieties, in particular the characteristic O-H···N intramolecular H-bond interaction between these two groups The intramolecular H-bonding in o-hydroxyaryl Schiff bases has been shown very much sensitive to the characteristics of the substituents present in the phenolic azomethine groups, with bulky substituents forcing these groups to lost co-planarity leading to a remarkable lengthening and weakening of the H-bond [19][20][21]. The intr lecular proton transfer in this type of compounds depends also on solvent polarity temperature [22,23], while optical excitation of the intramolecularly H-bonded Schiff can lead to excited state intramolecular proton transfer (ESIPT) [24][25][26][27], which may be accompanied by isomerization around the azomethine bond [28][29][30].
One of the most interesting properties that is in general exhibited by o-hydrox Schiff base derivatives is solvatochromism, which ultimately is a manifestation of th ative populations of the enol-imine and keto-amine tautomers present in the differen dia [14][15][16][17]. Most of times, solvatochromism is investigated by room temperature tronic absorption and fluorescence emission spectroscopies, and the simultaneous ence of several isomers and conformers of each tautomer in the samples has been co ered a factor leading to an increased complexity in the analysis of the experimental [15][16][17]. Hence, creating conditions where the number of conformers might be red appears of interest to undertake a detailed structural characterization of the tautom forms. On the other hand, the understanding of the main characteristics of the confo tional space of this type of compounds is instrumental for the understanding of their p erties. For this purpose, minimization of the intermolecular interactions is of fundam importance.
Matrix isolation is an experimental sampling technique that suits all the requirem stated above. In one side, being a low temperature technique (the work temperatu The keto-amine tautomer exists in the Z and E isomers as shown in the left and right bottom images, respectively. However, the geometric isomers of this tautomer do not have a direct correspondence to the geometric isomers of the enol-imine tautomer. To avoid confusion, in this article the keto-amine conformers will be grouped in two sets, one related to the E-enol-imine tautomer and other to the Z-enol-imine tautomer, which will be designated as "keto-amine forms of the E-enol-imine tautomer" and "ketoamine forms of the Z-enol-imine tautomer", respectively. In both sets, there are Z-keto-amine and E-keto-amine conformers. The intramolecular H-bonding in o-hydroxyaryl Schiff bases has been shown to be very much sensitive to the characteristics of the substituents present in the phenolic and azomethine groups, with bulky substituents forcing these groups to lost co-planarity and leading to a remarkable lengthening and weakening of the H-bond [19][20][21]. The intramolecular proton transfer in this type of compounds depends also on solvent polarity and temperature [22,23], while optical excitation of the intramolecularly H-bonded Schiff base can lead to excited state intramolecular proton transfer (ESIPT) [24][25][26][27], which may also be accompanied by isomerization around the azomethine bond [28][29][30].
One of the most interesting properties that is in general exhibited by o-hydroxyaryl Schiff base derivatives is solvatochromism, which ultimately is a manifestation of the relative populations of the enol-imine and keto-amine tautomers present in the different media [14][15][16][17]. Most of times, solvatochromism is investigated by room temperature electronic absorption and fluorescence emission spectroscopies, and the simultaneous presence of several isomers and conformers of each tautomer in the samples has been considered a factor leading to an increased complexity in the analysis of the experimental data [15][16][17]. Hence, creating conditions where the number of conformers might be reduced appears of interest to undertake a detailed structural characterization of the tautomeric forms. On the other hand, the understanding of the main characteristics of the conformational space of this type of compounds is instrumental for the understanding of their properties. For this purpose, minimization of the intermolecular interactions is of fundamental importance.
Matrix isolation is an experimental sampling technique that suits all the requirements stated above. In one side, being a low temperature technique (the work temperature is usually of a few degrees Kelvin), conditions are created that allow reducing the number of conformers present in the samples, and, on the other side, intermolecular interactions between the matrix-isolated molecules and the inert matrix medium (in general a frozen noble gas, e.g., argon) are minimal [31][32][33][34]. Matrix isolation has also many practical advantages compared with other possible approaches (e.g., molecular beams), in particular when the probing method is a spectroscopic technique [31,32]. Moreover, the concerted use of matrix isolation infrared spectroscopy and quantum chemical electronic structure calculations has been proven to be a very powerful tool for structural studies on conformationally flexible and multi-isomeric complex systems [31][32][33][34][35][36][37].
The study reported herein focused on 3-(5-bromo-2-hydroxybenzylideneamino)phenol (abbreviated as BHAP, Figure 1), which is a conformationally flexible o-hydroxyaryl Schiff base that exhibits tautomerism and geometric isomerism (E/Z). The conformational spaces of the two main tautomers of the molecule were investigated using electronic structure quantum chemical calculations undertaken at the DFT/B3LYP/6-311++G(d,p) level of theory, and the compound was experimentally studied under matrix isolation conditions by infrared spectroscopy. As it will be described in details in the next sections, it was concluded that the most stable species of BHAP is its E-enol-imine form, which has 4 lowenergy intramolecularly hydrogen-bonded forms (I, II, V and VI) which are of experimental relevance. These four conformers are concluded to exist in the gas phase with nearly equal populations, while only two of them (II and V) could be efficiently trapped in a cryogenic (10 K) argon matrix. Conformers I and VI convert to conformer II and V, respectively, by quantum mechanical tunneling through the barrier associated with the rotation of the OH phenolic group around the C-O bond. Natural bond orbital (NBO) analysis was performed on the most stable conformer of the experimentally relevant isomeric form (E-enol-imino conformer V) and used to shed light on details of its electronic structure. As a whole, this investigation is a comprehensive structural study of BHAP, and stresses the relevance of the O-H···N intramolecular interaction in determining the structure and properties of o-hydroxyaryl Schiff base compounds. of conformers present in the samples, and, on the other side, intermolecular interactions between the matrix-isolated molecules and the inert matrix medium (in general a frozen noble gas, e.g., argon) are minimal [31,32,33,34]. Matrix isolation has also many practical advantages compared with other possible approaches (e.g., molecular beams), in particular when the probing method is a spectroscopic technique [31,32]. Moreover, the concerted use of matrix isolation infrared spectroscopy and quantum chemical electronic structure calculations has been proven to be a very powerful tool for structural studies on conformationally flexible and multi-isomeric complex systems [31,32,33,34,35,36,37].
The study reported herein focused on 3-(5-bromo-2hydroxybenzylideneamino)phenol (abbreviated as BHAP, Figure 1), which is a conformationally flexible o-hydroxyaryl Schiff base that exhibits tautomerism and geometric isomerism (E/Z). The conformational spaces of the two main tautomers of the molecule were investigated using electronic structure quantum chemical calculations undertaken at the DFT/B3LYP/6-311++G(d,p) level of theory, and the compound was experimentally studied under matrix isolation conditions by infrared spectroscopy. As it will be described in details in the next sections, it was concluded that the most stable species of BHAP is its E-enol-imine form, which has 4 low-energy intramolecularly hydrogen-bonded forms (I, II, V and VI) which are of experimental relevance. These four conformers are concluded to exist in the gas phase with nearly equal populations, while only two of them (II and V) could be efficiently trapped in a cryogenic (10 K) argon matrix. Conformers I and VI convert to conformer II and V, respectively, by quantum mechanical tunneling through the barrier associated with the rotation of the OH phenolic group around the C-O bond. Natural bond orbital (NBO) analysis was performed on the most stable conformer of the experimentally relevant isomeric form (E-enol-imino conformer V) and used to shed light on details of its electronic structure. As a whole, this investigation is a comprehensive structural study of BHAP, and stresses the relevance of the O-H … N intramolecular interaction in determining the structure and properties of ohydroxyaryl Schiff base compounds.

Synthesis
3-(5-Bromo-2-hydroxybenzylideneamino)phenol was synthesized by mixing the equivalent amounts of the substituted amine (3-aminophenol) and substituted aldehyde (5-bromo-2-hydroxybenzaldehyde) precursors in methanol, at ca. 45 °C (see Figure 1) [15,16,17]. 3-Aminophenol (1.091 g; 0.01 mol) and 5-bromo-2-hydroxybenzaldehyde (2.000 g; 0.01 mol) were first dissolved in 25 mL of methanol by heating, and the solutions were then added. After 1 h at ca. 45 °C, under slow stirring, the reaction completed and the precipitated product was filtered, purified by recrystallization from methanol, and dried in a vacuum desiccator at room temperature. The IR spectrum of the purified crystalline material (in a KBr pellet at room temperature), and the 1 H-and 13

Synthesis
3-(5-Bromo-2-hydroxybenzylideneamino)phenol was synthesized by mixing the equivalent amounts of the substituted amine (3-aminophenol) and substituted aldehyde (5-bromo-2-hydroxybenzaldehyde) precursors in methanol, at ca. 45 • C (see Figure 1) [15][16][17]. 3-Aminophenol (1.091 g; 0.01 mol) and 5-bromo-2-hydroxybenzaldehyde (2.000 g; 0.01 mol) were first dissolved in 25 mL of methanol by heating, and the solutions were then added. After 1 h at ca. 45 • C, under slow stirring, the reaction completed and the precipitated product was filtered, purified by recrystallization from methanol, and dried in a vacuum desiccator at room temperature. The IR spectrum of the purified crystalline material (in a KBr pellet at room temperature), and the 1 H-and 13 C-NMR spectra in dimethylsulfoxide (DMSO-d 6 ) solution were found to be compatible with the desired product, and are provided in Supplementary Figures S1-S3. The IR spectrum was obtained in a Perkin Elmer FT-IR 100 spectrometer, and the 1 H-and 13 C-NMR spectra were obtained using a Bruker Biospin Ultrashield TM 300 MHz NMR spectrometer (Bruker, Karlsruhe, Germany), at room temperature. The melting point (m. p.) of the compound was measured using a Gallenkamp Sanyo Heater.

Matrix-Isolation Experiments
To prepare the cryogenic matrices, BHAP was placed in a homemade Knudsen cell assembled inside the vacuum chamber of the cryostat (APD Cryogenics closed-cycle helium refrigeration system with a DE202 expander), and was sublimated through thermoelectrical heating with help of a DC Power Supply VITEECOM (model 75-HY5003). The vapor of the compound was then co-deposited with argon (Air Liquide, N60), in a~1:1000 solute:host molar ratio, onto the optical CsI substrate of the cryostat, kept at 10.0 ± 0.1 K.
The IR spectra of the matrices were obtained using a Nicolet 6700 FTIR spectrometer, equipped with a mercury cadmium telluride (MCT) detector and a KBr beam splitter. The IR spectra were recorded in the 4000-500 cm −1 range, with 0.5 cm −1 resolution. The instrument was purged by a stream of dry/CO 2 -filtered air in order to avoid interference from atmospheric H 2 O and CO 2 vapors.

Theoretical Calculations
The geometries and IR spectra (within the harmonic approximation) of the conformers of E-enol-imine and Z-enol-imine isomers of BHAP, and the corresponding keto-amine forms (Figures 2-4) were fully optimized at the DFT level of theory, with the B3LYP functional and the 6-311++G(d,p) basis set [38][39][40][41][42][43]. The harmonic vibrational frequencies were scaled by the factors 0.994 and 0.945, below and above 1800 cm −1 , respectively, to approximately correct them for the neglected anharmonicity and method/basis set limitations. These scaling factors were used throughout the study, being applied to the calculated spectra of all structures considered. Relaxed potential energy profiles for the relevant internal rotations in the BHAP molecule were calculated at the same theory level. All calculations were performed with GAUSSIAN 09 (revision A.02) [44]. Assignment of the vibrational modes was carried out with the help of the animation module of GaussView (version 5.0) program [45]. In the simulated spectra presented in the figures, the IR bands were broadened by Lorentzian profiles (fwhm = 10 cm −1 ) centered at the calculated (scaled) frequencies, using the program ChemCraft [46].    Table  S1. The conformationally relevant dihedral angles α, β, ϒ, δ are indicated over the structure of conformer I. Note that all conformers have a symmetry related structure.  Table  S1. The conformationally relevant dihedral angles α, β, ϒ, δ are indicated over the structure of conformer I. Note that all conformers have a symmetry related structure.   Table  S1. The conformationally relevant dihedral angles α, β, ϒ, δ are indicated over the structure of con former I. Note that all conformers have a symmetry related structure.  Table  S1. The conformationally relevant dihedral angles α, β, ϒ, δ are indicated over the structure of con former I. Note that all conformers have a symmetry related structure.  Figure 3. Conformers of the BHAP Z-enol-imine form. The B3LYP/6-311++G(d,p) calculated optimized Cartesian coordinates of each conformer are provided in the Supporting Information Table S1. The conformationally relevant dihedral angles α, β, γ, δ are indicated over the structure of conformer I. Note that all conformers have a symmetry related structure.   Table S1. The conformationally relevant dihedral angles α, β, ϒ are indicated over the structure of conformer I (A). Note that all conformers have a symmetry related structure.

Structures and Conformational Landscapes of BHAP Isomeric Forms
The two geometrical enol-imine isomers of the BHAP differ in the placement of the substituents around the azomethine C=N bond, which leads to isomers E and Z. These isomers have four conformationally relevant internal axes of rotation: two of these (α and δ in Figures 2 and 3) correspond to the rotation around the two phenolic C-O bonds and define the orientation of the phenolic OH hydrogen atoms, while the remaining two axes of internal rotation (β and ϒ in Figures 2 and 3) are associated with the C-N and C-C exocyclic bonds of the bridge that connects the two phenolic rings, thus determining the relative orientation of the two aromatic rings. The performed B3LYP/6-311++G(d,p) calculations predicted the existence of 16 distinct conformers for each enol-imine geometrical isomer, whose optimized structures are depicted in Figures 2 and 3. The conformationally   I  II  III  IV   V  VI  VII  VIII   I  II  III  IV   V  VI  VII Table S1. The conformationally relevant dihedral angles α, β, γ are indicated over the structure of conformer I (A). Note that all conformers have a symmetry related structure.

Structures and Conformational Landscapes of BHAP Isomeric Forms
The two geometrical enol-imine isomers of the BHAP differ in the placement of the substituents around the azomethine C=N bond, which leads to isomers E and Z. These isomers have four conformationally relevant internal axes of rotation: two of these (α and δ in Figures 2 and 3) correspond to the rotation around the two phenolic C-O bonds and define the orientation of the phenolic OH hydrogen atoms, while the remaining two axes of internal rotation (β and γ in Figures 2 and 3) are associated with the C-N and C-C exocyclic bonds of the bridge that connects the two phenolic rings, thus determining the relative orientation of the two aromatic rings. The performed B3LYP/6-311++G(d,p) calculations predicted the existence of 16 distinct conformers for each enol-imine geometrical isomer, whose optimized structures are depicted in Figures 2 and 3. The conformationally relevant dihedral angles for these conformers are given in Table 1, while the corresponding electronic energies (with and without zero-point correction), Gibbs standard energies (at 298.15 K), and dipole moments are given in Table 2.
As shown in Table 2, the B3LYP/6-311++G(d,p) calculations predicted E-enol-imine conformer V as the lowest energy species among all possible structures. All the conformers of the Z-enol-imine isomer have much higher energies than the E-enol-amine V form, the lowest energy Z-enol-imine conformer being conformer Z-enol-imine V, which stays 62 kJ mol −1 above E-enol-imine V. Notably, conformer Z-enol-imine V is in fact higher in energy than all E-enol-imine forms. The main reason for the high energy of the Z-enolimine conformers is the repulsion between the two phenolic rings, which in the Z-isomer stay in the same side of the azomethine C=N bond (see Figure 3). Table 1. B3LYP/6-311++G(d,p) optimized conformationally-relevant dihedral angles of the E-enol-imine and Z-enol-imine conformers of BHAP a . In the case of the conformers of Z-enol-imine, the first set of energy values are relative to the most stable form, E-enol-imine V, and the second set are relative to the most stable conformer of the Z-enol-imine isomer.

Dihedral Angles
The four most stable conformers of the E-enol-imine isomer (conformers I, II, V and VI) are more than 38 kJ mol −1 lower in energy than the next conformer in the increasing order of energy (conformer XV). These conformers exhibit an intramolecular hydrogen-bond interaction of the O-H···N type (see Figure 2), which leads to the formation of a nearly planar pseudoaromatic six-membered ring where the dihedral angles γ and δ are both~0 • , and that is responsible for their stabilization [47]. Within the two pairs of conformers (I, II) and (V, VI), the main structural difference is the orientation of the non-hydrogen bonded phenolic OH group (the α dihedral angle being near 0 or 180 • ), while the orientation of the two aromatic rings is the same (i.e., the dihedrals β and γ are similar in each pair of conformers). The two pairs differ essentially in the value of the dihedral angle β, which is ca. 37 • in I and II and about −144 • in V and VI. Note that in all the low-energy intramolecularly H-bonded E-enol-imine conformers the bromo-substituted phenolic ring is roughly planar with the azomethine moiety (γ is~0 • ) due to the involvement of these two molecular fragments in the hydrogen bond, while the non-substituted phenol ring deviates considerably from the plane formed by the bromo-substituted phenolic ring and azomethine group, mostly because of the repulsive interaction between the azomethine hydrogen atom and the closest located hydrogen atom of the unsubstituted phenolic moiety (see Figure 2).
The four hydrogen bonded E-enol-imine conformers have predicted energies within less than 2 kJ mol −1 , which after zero-point and thermal corrections lead to standard Gibbs energies (at 298.15 K) that are within less than 1.5 kJ mol −1 , i.e., in practical terms the calculations predict the four conformers as being nearly isoenergetical. The electronic energy of conformer V was predicted to be the lowest (followed by those of conformers II, VI and I, which are higher than that of conformer V by 0.10, 0.84, and 1.69 kJ mol −1 , respectively), and inclusion of the zero-point correction does not change the order of stability of the conformers nor their relative energies significantly (see Table 2). However, conformers I and II have slightly larger entropy than forms V and VI, so that their relative standard Gibbs energies at 298.15 K comparatively reduce. In particular, this makes conformer II to become the lowest Gibbs energy form at the considered temperature. Nevertheless, all differences in the energies and entropy of the four intramolecularly Hbonded, low-energy conformers of the E-enol-imine isomer of the studied Schiff base are very small and, as mentioned above, in practical terms the four conformers shall be considered nearly isoenergetical and expected to be approximately equally populated in the gas phase equilibrium (i.e., each conformer shall have a population of approximately 25%, since all the other conformers have much higher energy and can be expected to have negligible populations).
The intramolecular H-bond in the four experimentally relevant E-enol-imine conformers (I, II, V, VI) is predicted by the calculations to be of nearly the same strength, taking into account the associated geometric parameters: in all forms, the O-H bond length is 0.993 Å (considerably longer than that of the free phenolic OH group that is 0.962 Å long), and the H···N and O···N distances are within the narrow ranges of 1.742-1.746 and 2.632-2.635 Å, respectively. The O-H···N angle is also very much similar in all the four conformers, being within the 146.9-147.1 • range. In each conformer, the two C-O bonds are predicted to be considerably different from each other, with that associated with the H-bonded phenolic group being considerably shorter (1.339 Å in all the conformers) than that associated with the free OH group (1.368 Å in I and II and 1.367 Å in V and VI). The shortening of both C-O and O-H bonds in the H-bonded phenolic moiety in comparison to the second phenol group amounts to~0.03 Å, respectively.
Conformers III, IV, VII and VIII of the E-enol-imine tautomer are higher in energy than the most stable form by over 55 kJ mol −1 , due to the breaking of the intramolecular O-H···N bond and its replacement by the repulsive interaction between the nitrogen and oxygen lone electron pairs. Among the high-energy conformers, forms XI, XII, XV and XVI are the less energetic (relative energies between 38 and 40 kJ mol −1 ; see Table 2) because they bear a C-H···O(H) attractive interaction between the CH group of the azomethine bridge and the phenolic oxygen atom (of the bromo-substituted phenolic ring), while conformers IX, X, XIII and XIV possess a C-H···H-O repulsive interaction and are intermediate in energy, with relative energies between 44 and 47 kJ mol −1 . The C-C-O-H (δ) and C-C-C=N (γ) dihedral angles in forms IX, X, XIII and XIV clearly reveal the relevance of the C-H···H-O repulsive interaction: in these conformers, (i) the OH group is tilted from the plane of the aromatic ring (δ is in the range 6-9 • ), while it is deviated by less than 1 • from the ring plane in all the other E-enol-imine conformers, and (ii) the dihedral angle γ, which is a measure of the coplanarity between the azomethine group and the bromo-substituted phenol ring, have values corresponding to geometries where the two groups deviate from by 7-9 • from coplanarity, while for all the other E-enol-imine conformers the deviation is in the range of 0-4 • (see Table 1).
In the case of the Z-enol-imine isomer, none of the sixteen distinct conformers have any strong intramolecular H-bond, this being one of the main reasons for their much higher relative energies (63-75 kJ mol −1 ) in comparison with the four lowest energy conformers of the E-enol-imine isomer (together with the already mentioned repulsions between the two phenolic rings). The most stable Z-enol-imine conformers (all with energies~63 kJ mol −1 relatively to E-enol-imine V) are forms I, II, V and VI, and possess a weak O-H···π intramolecular interaction (see Figure 3 and Table 2), which justifies their lower energy compared with the remaining Z-enol-imine conformers. Conformers XI, XI, XV, and XVI are the next in order of growing energy (~67 kJ mol −1 ). Like the related conformers of E-enol-imine, these conformers bear a weak C-H···O(H) attractive interaction between the CH group of the azomethine bridge and the phenolic oxygen atom of the bromo-substituted phenol ring. Also as their analogues in E-enol-imine, conformers Z-enol-imine IX, X, XIII and XIV (~70 kJ mol −1 ) possess a C-H···H-O repulsive interaction, which likewise reflects in the values of the C-C-O-H (δ) and C-C-C=N (γ) dihedral angles that evidentiate the tilt of the OH group out of the plane of the aromatic ring (δ is in the range 7-10 • ) and the strong deviation from coplanarity (by 25-31 • ) of the azomethine and bromo-substituted phenol moieties (see Table 1). It is interesting to note that in the case of the Z-enol-imine conformers, the most out of the plane tilted OH group, as measured by the dihedral angle δ, are conformers I, II, V, and VI, due to the fact that in these conformers the OH group has to adapt its geometry to better fulfil the O-H···π stabilizing intramolecular interaction (see Figure 3). This fact also explains the strong non-coplanarity of the azomethine and bromo-substituted phenol moieties in these conformers, as measured by the dihedral angle γ (see Table 1). The highest energy Z-enol-imine conformers (73-75 kJ mol −1 ) are forms III, IV, VII, and VIII, mostly due to effect of the repulsive interaction between the lone electron pairs of the azomethine nitrogen and bromo-substituted phenolic oxygen atoms (see Figure 3). It is in these conformers that the dihedral angle γ assumes the values (~47 • ) that correlate with a maximum deviation from coplanarity of the azomethine and bromo-substituted phenolic fragments.
Considering only the Z-enol-imine conformers, the energies are within 12 kJ mol −1 , i.e., they span a much narrower range of values compared to what happens for E-enol-imine. The main reason for this is that in Z-enol-imine no conformers exist that are stabilized by a strong intramolecular H-bond (if we do not count with those, the energies of the E-enol-imine conformers differ only by 19 kJ mol −1 , nevertheless still a somewhat larger value compared to what is observed for Z-enol-imine).
A final note shall be done regarding the geometries around the N=C double bond. For both Eand Z-enol-imine isomers, all conformers show the azomethine bridge slightly non-planar. The C-C=N-C dihedral stays in the range~5-8 • in the case of the Z-enol-imine conformers, where the repulsions between the two phenol moieties dominate, and betweeñ 2-4 • in the case of the E-enol-imine conformers. This structural feature has already been pointed out before for other aryl Schiff bases [26,[28][29][30]47,48].
The keto-amine tautomer of BHAP can also exist in two geometrical isomeric species (E and Z) differing in the placement of the substituents around the C=C bond of the bridge. Since the geometric isomers of this tautomer do not have a direct structural correspondence to the geometric isomers of the enol-imine tautomer, in this article the keto-amine conformers of BHAP have been grouped in two sets, one structurally related to the E-enol-imine tautomer and the other to the Z-enol-imine tautomer, which will be designated as "keto-amine forms structurally related to the E-enol-imine tautomer" and "keto-amine forms structurally related to the Z-enol-imine tautomer", respectively. In both sets, there are Z-keto-amine and E-keto-amine conformers. In the keto-amine forms, there are three conformationally relevant internal axes of rotation: one of these (α in Figure 4) corresponds to the rotation around the phenolic C-O bond and define the orientation of the phenolic hydrogen atom, while the others (β and γ in Figure 4) are associated with the C-N bonds of the bridge (the first connecting the N atom to the phenol ring, and the second being the central bond of the bridge), and determine the relative orientation of the two aromatic rings. The performed B3LYP/6-311++G(d,p) calculations predicted the existence of 8 distinct conformers for each keto-amine geometrical isomer, which have been grouped in Figure 4 according to their structural relation with the two enol-imine isomers, as stated above. The conformationally relevant dihedral angles of the keto-amine conformers are provided in Table 3, and the corresponding electronic energies (with and without zero-point correction), Gibbs standard energies (at 298.15 K), and dipole moments are given in Table 4. Table 3. B3LYP/6-311++G(d,p) optimized conformationally-relevant dihedral angles of the ketoamine forms of the Eand Z-enol-imine isomers of BHAP a .  Figure 4. Note that all conformers have a symmetry related structure.  Figure 4 for the structures of the conformers. The first set of energy values is relative to the most stable form, E-enol-imine V, and the second set is relative to the most stable conformer of each keto-amine isomer. Values in italic style are for E-keto-amine forms and values in round style for E-keto-amine forms.

Dihedral Angles
The DFT performed calculations indicate that all keto-amine conformers structurally related to the E-enol-imino isomer have lower energies than those related with the Z-enolimino isomer (see Table 4). Conformers I, II, V, and VI of the first group are considerably lower in energy (by over~40 kJ mol −1 , compared to the remaining E-enol-imino related keto-amine conformers, and by over 80 kJ mol −1 , compared to all the Z-enol-imine related keto-amine conformers). This could be anticipated, since these conformers have a strongly stabilizing N-H···O=C intramolecular H-bond interaction.
The N-H···O=C H-bond interaction is indeed strong, as reflected in the H-bond structural parameters. The N-H bond length in all the keto-amine conformers having this interaction is 1.042 Å, and the H···O distance stays in the range 1.682-1.700 Å, which is considerably shorter than the sum of H+O van der Walls radii (1.20 + 1.52 = 2.72 Å). The O···N distance is in the 2.579-2.588 Å range, being shorter (as the H···O distance) in conformer V, and longer in conformer II, while the N-H···O angles are~140 • in all the four conformers.
As seen in Table 3, conformers I, II, V and VI are almost planar, whereas conformers III, IV, VII and VIII are considerably non-planar, exhibiting a dihedral angle β equal to 166.5, 169.9, −8.7, and −9.0 • , respectively, because of the repulsive interaction between the hydrogen linked to the carbon atom of the bridge and the closely located hydrogen of the phenol group (see Figure 4). All the keto-amino conformers related to the Z-enol-imine isomer are strongly non-planar due to the close proximity of the two phenolic moieties.
Looking now to the energy differences within the conformers of the E-keto-amine and Z-keto-amine forms (these values are given in the last columns of Table 4 and distinguished using italic type and round type respectively), one can see that in each case the four conformers related to the E-enol-imine form are considerably more stable than those related to the Z-enol-imine isomer. In the case of the E-keto-amine isomer, the energies of the conformers span through a range of~37 kJ mol −1 , the four most stable conformers being within a~3 kJ mol −1 energy range. For the Z-keto-amino isomer, the energies of the conformers span for a range of~66 kJ mol −1 , the four most stable conformers being again within a~3 kJ mol −1 energy range.
A general conclusion that can be extracted from the theoretical calculations is that, when the Eand Zisomeric forms of the two tautomers (enol-imine and keto-amine) are considered as independent species, the gas phase equilibrium populations of the corresponding most stable four conformers shall be roughly equal to~25% each. As described below, this is an important piece of information for the analysis of the obtained matrix isolation experimental data.

Natural Bond Orbital Analysis
In order to look into more detail to the electronic structure characteristic of BHAP, natural bond orbital (NBO) analysis [49] was performed using the DFT/B3LYP/6-311++G(d,p) data. The lowest energy conformer E-enol-imine V was used as target for these calculations.
The calculated NBO atomic charges are listed in Table 5. As expected, the NBO atomic charges show that the hydrogen-bonded phenolic OH group (O15-H26, according to the adopted atom numbering scheme, shown in the insert of Table 5) is considerably more polarized than the free phenolic OH group (O17-H27), due to the involvement of the first in the hydrogen bond: the positive charge of the H-bonded H27 atom is considerably larger than that of H26 (+0.507 vs. +0.468 e), but the negative charge of O15 is also slightly more negative than that of O17 (-0.677 vs. -0.670 e).
The charges of the aromatic hydrogen atoms vary in the range +0.207 to +0.224 e, the more positive ones belonging to the hydrogen atoms of the bromo-substituted phenolic ring, while that of the hydrogen atom of the azomethinic bridging group is +0.170 e. The aromatic carbon atoms have calculated negative charges within the range −0.171 to −0.278 e, except that connected to the Br atom, which has a less negative charge (−0.138 e), and those connected to the oxygen and nitrogen atoms, which are positive. As expected, the carbon atom connected to the more negative O15 atom has a larger positive charge (+0.376 e) than that connected to O17 (+0.325 e), while the carbon atoms linked to the nitrogen atom (whose charge is −0.522 e, i.e., less negative than those of the oxygen atoms) have positive charges, +0.154 e (C5) and +0.162 e (C8). The Br atom is predicted to have a slightly positive charge, pointing to a more pronounced π-donation to the ring, compared to the σ-withdrawn. to the adopted atom numbering scheme, shown in the insert of Table 5) is considerably more polarized than the free phenolic OH group (O17-H27), due to the involvement of the first in the hydrogen bond: the positive charge of the H-bonded H27 atom is considerably larger than that of H26 (+0.507 vs. +0.468 e), but the negative charge of O15 is also slightly more negative than that of O17 (-0.677 vs. -0.670 e). The charges of the aromatic hydrogen atoms vary in the range +0.207 to +0.224 e, the more positive ones belonging to the hydrogen atoms of the bromo-substituted phenolic ring, while that of the hydrogen atom of the azomethinic bridging group is +0.170 e. The aromatic carbon atoms have calculated negative charges within the range −0.171 to −0.278 e, except that connected to the Br atom, which has a less negative charge (−0.138 e), and those connected to the oxygen and nitrogen atoms, which are positive. As expected, the carbon atom connected to the more negative O15 atom has a larger positive charge (+0.376 e) than that connected to O17 (+0.325 e), while the carbon atoms linked to the nitrogen atom (whose charge is −0.522 e, i.e., less negative than those of the oxygen atoms) have positive charges, +0.154 e (C5) and +0.162 e (C8). The Br atom is predicted to have a slightly positive charge, pointing to a more pronounced π-donation to the ring, compared to the σ-withdrawn.
The nature of the NBOs (composition in terms of atomic orbitals) provides an elegant way to define the hybridization adopted by the different atoms in the molecule. Table 6 shows the description of the occupied bonding NBOs and those associated with lone electron pairs in terms of the atomic contributions and the corresponding occupancies. The descriptions allow to extract information about the hybridization assumed by the different atoms in the molecule.
All carbon atoms except C1 exhibit the expected hybridization pattern, using a p orbital to form the π component of the double bonds, and typical sp 2  The nature of the NBOs (composition in terms of atomic orbitals) provides an elegant way to define the hybridization adopted by the different atoms in the molecule. Table 6 shows the description of the occupied bonding NBOs and those associated with lone electron pairs in terms of the atomic contributions and the corresponding occupancies. The descriptions allow to extract information about the hybridization assumed by the different atoms in the molecule.
All carbon atoms except C1 exhibit the expected hybridization pattern, using a p orbital to form the π component of the double bonds, and typical sp 2 hybrid orbitals to form the σ bonds (including the σ component of the double bonds). In all cases, the sp 2 hybrids used to establish a bond with a hydrogen atom have slightly larger p character, while those used to form the σ bonds with other carbon atoms have a larger s contribution. The p contribution to the hybrid orbitals are large when the carbon is connected to the more electronegative Br, N, and O atoms, the extreme case corresponding to the bond between C1 and O17, where the carbon atom uses an essentially pure p orbital. The orbital mixing calculated for C1 is in fact unusual, since it uses a p orbital to establish the linkage with O17, two hybrids orbitals with substantial s character to form the σ bonds with the vicinal ring carbon atoms (C2 and C6) and an orbital with large p character, but still significant s contribution, to form the π bond in which it is involved. The pattern of hybridization adopted by the two oxygen atoms is also interesting, though both can be considered to be formally hybridized sp 2 . O15 has its lone pairs occupying a p and a sp 2 hybrid orbital (the latter with an increased s character), the two additional hybrid orbitals being used to form the bonds with C14 and H26, the last one having an increased p character (approaching the typical p composition of a sp 3 hybrid orbital). In turn, O17 has its lone pairs occupying sp 2 hybrid orbitals, one of them with a large s character, uses the third hybrid orbital, which has a large p character, to establish the bond with the hydrogen atom, and binds to the ring carbon atom (C1) using a p orbital. The difference revealed by the NBO analysis in the electronic structures of the two oxygen atoms is most probably a result of the involvement of O15 in the strong intramolecular O-H···N hydrogen bond, which also has implications on the electronic structure around the N atom. Indeed, the N atom has a hybridization rather close to sp, with two essentially pure p orbitals being used to accommodate the lone pair and form the π bond with C8, and the two hybrid orbitals to form the σ bonds with C5 and C8. a Occupancy is given with an exaggerated accuracy, as in the Gaussian output file. b The A and B values correspond to the contributions of the atomic orbitals of the two atoms forming a bond, by the order indicated in the first column. Numbering of atoms as in the in-set figure shown in Table 5.
Finally, the bromine atom is not hybridized, using one of its p orbitals to form the bond with C11, while the three lone pairs occupy the s and the remaining two p valence orbitals.
The bond orders (BO) were evaluated using the natural atomic orbitals Wiberg bond index [50]. The results are given in Table 7. The bond orders for the carbon-carbon and C-H bonds of the aromatic rings are in the range of 1.27-1.47 and equal to 0.92, respectively. In the case of the carbon-carbon bonds, the bonds with the smallest bond orders in the bromosubstituted phenol ring are the C9-C14, C9-C10, and C13-C14 bonds, while those with the largest bond orders are the C12-C13 and C10-C11 bonds. The first involve the carbon atoms bearing the OH and azomethine substituents, and participate in the 6-membered pseudo-ring containing the intramolecular hydrogen-bond; as expected these bonds are the longest bonds in the bromo-substituted phenol ring (1.419, 1407, and 1.400 Å, respectively). In turn, the largest bond order bonds are vicinal to the previous ones and the shortest bonds in the ring (1.386 and 1.381 Å, respectively). a Numbering of atoms as in the in-set figure shown in Table 5.
A similar trend is observed for the second phenol ring, with the bonds involving C5 and C1, which bear the substituents, having the smallest bond orders and being the longest bonds. The bond order of the exocyclic C8-C9 bond is somewhat larger than the unity (1.14), expressing the π delocalization between the bromo-substituted phenol ring and the azomethine bridge, which are co-planar. This is also reflected in the calculated bond order for the central N7=C8 bond (1.70), which is smaller than the formal value for a non-delocalized double bond. In turn, the bond order of the C5-N7 bond is close to 1, since the phenol ring and the azomethine group are considerably tilted in relation to each other (see Table 1 and discussion on Section 3.1) and delocalization is then precluded. The two C-O bonds have bond orders close to the unity, but that of the C1-O17 bond is somewhat higher than the one of the C14-O15 bond, in consonance with their relative lengths (1.367 vs. 1.340 Å). On the other hand, the bond orders and corresponding bond lengths of the two O-H bonds are very much different, in result of the participation of the O15-H26 group in the intramolecular hydrogen bond. While for this group the bond order is only 0.64 and the bond length is 1.993 Å, for the O17-H27 group the bond order is 0.77 and the bond length is 1.963 Å. The calculated bond order for the N···H hydrogen bond is 0.09, which is a value similar to those found for other molecules bearing a strong intramolecular hydrogen bond [51].

Infrared Spectrum of BHAP in an Ar Matrix
As described in the Section 2.2 vapors of BHAP were co-deposited with a large excess of argon onto a cold (10 K) CsI substrate mounted at the cold tip of the cryostat. The purposes of this experiment were: (i) to verify the nature of the species trapped in the matrix; (ii) to perform the vibrational study of these species. A few considerations shall be made a priori. Firstly, the four isomeric species discussed in the previous sections (E-and Z-enol-imine and Eand Z-keto-amine isomers) cannot interconvert in the gas phase nor in the matrices, since the energy barriers associated with these processes are very high (over a few hundreds of kJ mol −1 [52], because they imply breaking of a double bond, N=C or C=C). Secondly, the IR data obtained for the crystalline material in a KBr (room temperature; Supplementary Figure  S1) indicate that in the crystal an enol-imino form is present, since no band ascribable to the stretching vibration of the carbonyl group of the keto-amine forms is observed. Because the E-enol-imine isomer is considerably more stable that the Z-enol-imine one, the first form could be expected to be the one resulting from the synthesis. Nevertheless, this hypothesis should be verified. Finally, as already mentioned, the conformational population in the gas phase equilibrium, after sublimation of the compound, shall comprehend conformers I, II, V, and VI in nearly equal amounts. This applies for the E-enol-imine isomer, but also for the Z-enol-imine isomer in the case this higher energy species reveals to be experimentally relevant.
It is also worth mentioning that, in a matrix isolation experiment, higher-energy conformers may convert into lower energy forms during deposition (conformational cooling effect), if the energy barriers separating the conformers are small enough (a few kJ mol −1 ) [33,34,[53][54][55][56]. In addition, moderate size barriers (up to ca. 20-25 kJ mol −1 ) associated with transformations involving only movement of a hydrogen atom may also allow conformational decay by quantum mechanical tunneling through the barrier, this phenomenon being common for molecules bearing the hydroxyl moiety, in particular for phenols and naphtols and for carboxylic acids [57][58][59][60][61][62]. In order to estimate the energy barriers separating the low-energy conformers of BHAP, relaxed potential energy scans were performed (at the B3LYP/6-311++G(d,p) level) along the relevant coordinates. In the studied molecule, the four low-energy E-enol-imine conformers (I, II, V and VI) can be interconverted by two ways: (i) by rotation of the phenolic OH around the C-O bond; and (ii) by rotation of the whole phenol group around the C-N bond connecting to phenol and azomethine fragments. The corresponding potential energy scan profiles are shown in Figure 5. still large enough to preclude the conformational cooling to take place during deposition (at least in a significant way) [53]. So, taking into account only this conformational path, in principle, all 4 low-energy conformers should be trapped in the cryogenic matrix. On the other hand, the internal rotation of the OH group implies to overcome a barrier of ~15 kJ mol -1 ( Figure 5A). This barrier is enough high to preclude the OH group to rotate at the cryogenic temperature of the matrix in an over-the-barrier process. However, as mentioned above, hydrogen tunneling through a barrier of such height is highly probable, so that conformers I and VI can be expected to quickly convert into conformers II and V, respectively. This has been observed to happen in phenols and other similar compounds [59,60]. As a whole, the calculated data on the energy barriers then support the presence of only conformers II and V of the E-enol-imine isomer in the matrix in nearly equal amounts.  The conformational isomerization energy barriers associated with the internal rotation around the C5-N bond are~7-9 kJ mol −1 ( Figure 5B). These barriers are not large, but still large enough to preclude the conformational cooling to take place during deposition (at least in a significant way) [53]. So, taking into account only this conformational path, in principle, all 4 low-energy conformers should be trapped in the cryogenic matrix. On the other hand, the internal rotation of the OH group implies to overcome a barrier of 15 kJ mol −1 ( Figure 5A). This barrier is enough high to preclude the OH group to rotate at the cryogenic temperature of the matrix in an over-the-barrier process. However, as mentioned above, hydrogen tunneling through a barrier of such height is highly probable, so that conformers I and VI can be expected to quickly convert into conformers II and V, respectively. This has been observed to happen in phenols and other similar compounds [59,60]. As a whole, the calculated data on the energy barriers then support the presence of only conformers II and V of the E-enol-imine isomer in the matrix in nearly equal amounts. Figure 6 shows the infrared spectrum of BHAP isolated in an argon matrix, at 10 K. In the figure, the sum spectrum of the B3LYP/6-311++G(d,p) calculated spectra for E-enolimine conformers II and V is also presented. Figure S4, in the Supporting Information, shows also the sum spectra of the B3LYP/6-311++G(d,p) calculated infrared spectra for Z-enol-imine conformers II and V, and also those for Eand Z-keto-amine isomers built based on the calculated spectra of their hypothetically experimentally relevant conformers, for comparison.  Figure 6 shows the infrared spectrum of BHAP isolated in an argon matrix, at 10 K. In the figure, the sum spectrum of the B3LYP/6-311++G(d,p) calculated spectra for E-enolimine conformers II and V is also presented. Figure S4, in the Supporting Information, shows also the sum spectra of the B3LYP/6-311++G(d,p) calculated infrared spectra for Zenol-imine conformers II and V, and also those for E-and Z-keto-amine isomers built based on the calculated spectra of their hypothetically experimentally relevant conformers, for comparison.
As it can be seen in Figure 6 and Supplementary Figure S4, the data clearly shows that in the matrix only the E-enol-imine isomer is present, the spectra of the remaining isomeric species strongly differing from the experimental spectrum. The good reproduction of the experimental spectrum by the simulated one based on conformers II and V of the E-enol-imino isomer indicated also that the considerations above were most probably correct. However, because the infrared spectra of the four low-energy conformers of the E-enol-imine isomer do not differ significantly, it is not possible, based only on the comparison of the experimental and simulated spectra, to conclude with all certainty that conformers I and VI are not also contributing to the observed spectrum. Nevertheless, while the frequencies are very similar all among the four conformers, intensities are somehow different in several spectral regions. Thus, annealing of the matrix should, in principle, allow for isomerization of the putative conformers I and VI into conformers V and II, respectively, through internal rotation around the C5-N bond (see Figure 5B). Nonetheless, annealing of the matrix up to 25 K did not lead to any spectral changes, strongly supporting the sole presence of conformers II and V in the matrix.
. Figure 6. Infrared spectrum of BHAP in an Ar matrix (10 K) in the 1800-500 cm −1 range (bottom), and sum spectrum of the B3LYP/6-311++G(d,p) calculated spectra for E-enol-imine conformers II and V. As it can be seen in Figure 6 and Supplementary Figure S4, the data clearly shows that in the matrix only the E-enol-imine isomer is present, the spectra of the remaining isomeric species strongly differing from the experimental spectrum. The good reproduction of the experimental spectrum by the simulated one based on conformers II and V of the E-enolimino isomer indicated also that the considerations above were most probably correct. However, because the infrared spectra of the four low-energy conformers of the E-enolimine isomer do not differ significantly, it is not possible, based only on the comparison of the experimental and simulated spectra, to conclude with all certainty that conformers I and VI are not also contributing to the observed spectrum. Nevertheless, while the frequencies are very similar all among the four conformers, intensities are somehow different in several spectral regions. Thus, annealing of the matrix should, in principle, allow for isomerization of the putative conformers I and VI into conformers V and II, respectively, through internal rotation around the C5-N bond (see Figure 5B). Nonetheless, annealing of the matrix up to 25 K did not lead to any spectral changes, strongly supporting the sole presence of conformers II and V in the matrix.
The molecule of BHAP has 27 atoms, i.e., 75 fundamental vibrations, all active in IR, with 47 modes predicted by the calculations to appear in the 1800-500 cm −1 region. The proposed assignments are presented in Table 8. Below, some of the assignments are discussed briefly.  For full calculated spectra of the conformers see Table S2, in the Supporting Information. b Values in italic correspond to absorptions tentatively assigned to conformer II; values in round style are assigned to both conformers or to conformer V only (when a corresponding band ascribed to conformer II is indicated). c Approximate descriptions were obtained by visual inspection of the animated vibrations using GaussView 5.0. ν, stretching; δ, bending; γ, rocking/torsion; n.o., not observed. Subscript "p" refers to the phenol ring, "bp" to the bromo-substituted phenol ring, "am" to the azomethine bridging group.
Above 1800 cm −1 , only the OH and CH stretching vibrations are expected to absorb. The aromatic CH stretchings were found to span the~3100-2900 cm −1 range, while the azomethine CH stretching mode is observed as a pair of bands at 2864 and 2846 cm −1 , which are tentatively assigned to conformers V (calc. 2878 cm −1 ) and II (calc. 2877 cm −1 ), respectively. The stretching mode of the free OH group is assigned to the pair of bands at 3642 (II; calc. 2625 cm −1 ) and 3640 cm −1 (V; calc. 3623 cm −1 ), while the intramolecularly hydrogen bonded OH group gives rise to a broad and structured feature spreading over the very wide 3573-3035 cm −1 region, as it is normally observed for hydrogen bonded hydroxyl groups [63].
In the 1800-500 cm −1 spectral region, the experimental spectrum shows intense bands in the 1650-1590 and 1200-1150 cm −1 ranges, and around 1500, 1320, 960, 840, 760 and 680 cm −1 . This is in full agreement with the results of the calculations, as can be seen in Figure 6 and also in Table 8.
The 1650-1590 cm −1 region comprehends the bands due to the νN=C stretching mode (exp.: 1630 and 1626 cm −1 ; calc.: II, 1658 cm −1 and II, 1656 cm −1 ) and the highest frequency νCC stretching modes of the aromatic rings, some of them significantly mixed with the bending vibration (δOH) of the corresponding phenolic OH group. Indeed, according to the calculations, the intense band observed at 1593 cm −1 corresponds to the mode with a larger contribution of the bending mode of the H-bonded phenolic OH group. As expected, this mode appears considerably shifted to higher frequencies when compared to the equivalent vibration of the free OH group, which is assigned to the bands at 1198 and 1181 cm −1 , in conformer II (calc: 1195 cm −;1 ) and V (calc: 1187 cm −1 ), respectively.
The group of observed bands around 1500 cm −1 correspond essentially to mixed ring νCC stretching and δCH bending modes, in some cases also mixed with the δOH bending modes, which contribute to the intensification of the bands. The δCH bending vibration of the azomethine group is assigned to the experimental bands at 1368 and 1358 cm −1 , tentatively ascribed to conformer II and V, respectively (calc: 1374 and 1372 cm −1 ).
Among the vibrations contributing to the bands observed in the experimental spectrum at around 1320 cm −1 , the νCO stretching modes of the phenolic groups shall be highlighted. As expected, the higher polarized (but with a higher bond order) C-O bond of the H-bonded phenolic group gives rise to a considerably more intense band (observed at 1314 cm −1 ) than the C-O bond of the free phenolic group (observed at 1308 cm −1 for V and 1285 cm −1 for II).
In the 1200-1150 cm −1 range, besides the bands assigned to the δOH bending mode of the free phenol OH group in the two conformers, another intense band is observed at 1158 cm −1 , which, according to the calculations, is ascribable to a mixed δCH/νCO/νC-N mode (Table 8).
Finally, the intense experimental bands at ca. 960, 840, 760, and 680 cm −1 are due to a ring δCC/δCH deformation mode with a relevant contribution from the νCO streching of the free phenolic moiety, the γOH torsion of the hydrogen bonded OH group, the ring allin-phase out-of-the-plane γCH rocking modes, and ring γCC torsional modes, respectively. As expected [64], the γOH torsion of the hydrogen bonded OH group appears at a much higher frequency than the equivalent mode in the free OH group, which is predicted by the calculations at~320 cm −1 .

Conclusions
In this study, a detailed structural investigation of a new o-hydroxyaryl Schiff base was undertaken. Quantum chemical electronic structure calculations, performed at the DFT(B3LYP)/6-311++G(d,p) level, were used to study the conformational space of the two geometrical isomers (E and Z) of the enol-imine and keto-amine tautomers of the compound. These studies allowed to conclude that the most stable isomeric form is the E-enol-imino form, and that this form has four low-energy, intramolecularly H-bonded conformers which are nearly degenerate.
The infrared spectrum of the crystalline material (in a KBr pellet at room temperature) showed that the synthesized solid material was an enol-imino form of the compound, a result that was also in consonance with the 1 H-and 13 C-NMR spectroscopy data obtained in DMSO-d 6 solution. In turn, the low temperature matrix isolation infrared spectra demonstrated that the matrix isolated molecules, trapped from the room temperature Schiff base vapor, correspond to the E-enol-imine isomer in two specific conformers, II and V, which exist in the matrix in nearly equal amounts. The conjugation of these experimental results with simple chemical considerations on the structural relations between the enolimine and keto-amine isomers permitted to conclude that the synthesized crystalline material corresponds to the E-enol-imine tautomer of the studied compound, which upon sublimation gives rise to a vapor where the four low-energy conformers (I, II, V and VI) bearing a strong O-H···N intramolecular hydrogen bond coexist with identical populations in the equilibrium. Deposition of the vapor of the compound with a large excess of argon onto a substrate kept at 10 K gives rise to a solid argon matrix of the compound, where only conformers II and V subsist, conformers I and VI spontaneously converting into forms II and V, respectively, by the quantum mechanical tunneling mechanism. These conformational decays, involve the conformational change of the free OH phenolic group and takes place by hydrogen atom tunneling through a barrier for internal rotation about the C-O bond that is of ca. 15 kJ mol −1 .
Additional details of the electronic structure of the experimentally relevant isomeric species were obtained using the natural bond orbital (NBO) theory. This analysis targeted the most stable conformer of the E-enol-imino isomer (conformer V), and allowed to get insights on the hybridization states of the different atoms in the molecule, and their charges, as well as on the polarization and bond orders of the different bonds in the molecule, including the intramolecular hydrogen-bond. From all the different analyses performed, the structural relevance of the O-H···N intramolecular H-bond in the studied o-hydroxyaryl Schiff base was highlighted. Additionally, the matrix isolation infrared spectrum of the studied compound was fully assigned and interpreted with help of the DFT(B3LYP)6-311++G(d,p) calculated infrared spectra for the E-enol-imine conformers II and V.
Supplementary Materials: The following are available online. Figures S1-S3, with the IR spectrum of the compound (in KBr pellet; room temperature) and 1H and 13C NMR spectra (in DMSO; room temperature); Figure S4, with calculated IR spectra of the isomeric forms of the enol-imine and keto-amine tautomers of BHAP vs. the experimental matrix isolation IR spectrum of the compound; Table S1, with the Cartesian coordinates of the optimized calculated conformers of the Eand Zisomers of the enol-imine and keto-amine tautomers of BHAP; Table S2, with the calculated IR spectra of the 4 most stable conformers of the experimentally relevant E-enol-imine isomer of BHAP.