Local and Remote Conformational Switching in 2-Fluoro-4-Hydroxy Benzoic Acid

: In this work, 2-F-4-OH benzoic acid was isolated in Ar matrices and conformational changes were induced by near-IR irradiating the sample. Upon deposition, three conformers could be observed in the matrix, denoted as A1 , A2 , and D1 , respectively. A1 and A2 are trans carboxylic acids, i.e., there is an intramolecular H bond between the H and the carbonyl O atoms in the COOH group, whereas D1 is a cis carboxylic acid with an intramolecular H bond between the F atom and the H atom in the COOH group, which otherwise has the same structure as A1 . The difference between A1 and A2 is in the orientation of the carbonyl O atom with regard to the F atom, i.e., whether they are on the opposite or on the same side of the molecule, respectively. All three conformers have their H atom in their 4-OH group, facing the opposite direction with regard to the F atom. The stretching overtones of the 4-OH and the carboxylic OH groups were selectively excited in the case of each conformer. Unlike A2 , which did not show any response to irradiation, A1 could be converted to the higher energy form D1 . The D1 conformer spontaneously converts back to A1 via tunneling; however, the conversion rate could be signiﬁcantly increased by selectively exciting the OH vibrational overtones of D1 . Quantum efﬁciencies have been determined for the ‘local’ or ‘remote’ excitations, i.e., when the carboxylic OH or the 4-OH group is excited in order to induce the rotamerization of the carboxylic OH group. Both ‘local’ and ‘remote’ conformational switching are induced by the same type of vibration, which allows for a direct comparison of how much energy is lost by energy dissipation during the two processes. The experimental ﬁndings indicate that the ‘local’ excitation is only marginally more efﬁcient than the ‘remote’ one.

Introduction
To fully understand chemical reactions occurring in the environment, an investigation of the excited states of molecules and their relaxation is of the utmost importance. A vibrationally excited molecule may emit its excess energy via intra-or intermolecular processes. One example of the latter is when an excited molecule collides with a nearby one, transferring its excess energy in forms of vibrational, rotational, or translational energy. In contrast to this, during intramolecular vibrational energy transfer or intramolecular vibrational energy relaxation (IVR), the energy flows from one excited vibrational mode to another while the energy difference dissipates into the surroundings. This also means that in certain cases the excitation of a vibrational mode on one side of the molecule may result in a change in the geometry of a functional group on another, remote side of the molecule, called remote conformational switching.
The matrix isolation (MI) technique has been successfully used to study the conformational switching of molecules isolated in inert, nearly interaction-free matrices using mostly solidified noble or other inert gases (such as Ne, Ar, Kr, Xe, N 2 , etc.) as matrix material. The intermolecular relaxation of matrix-isolated molecules is greatly limited, thereby allowing for the examination of the IVR processes. Furthermore, another primary advantage of the technique is that the resulting IR spectrum contains sharp absorption peaks, meaning that they do not overlap; thus, the different conformers can be distinguished from each other, and the induced changes in their population can be monitored. It is important to note, however, that in order to achieve this vibrationally induced conformational change, selective excitation of the vibrational modes is required, which necessitates the use of monochromatic laser irradiation sources.
Of greater interest are those cases when the excitation and the rotamerization occurs in different parts of the molecule. Compared to its short-range counterpart, far fewer matrix-isolation studies have been devoted to these long-range IVR processes. The first one was cysteine, which showed the rotamerization of the thiol (-SH) group upon exciting the νOH stretching overtone (νOH), although it occurred along with other conformational changes in the molecule, thus rendering it not a selective method [21]. Nevertheless, this work pointed to the feasibility of remote switching using molecular antennas. The next step was to find a molecular system representing the first example for selective conformational control. Accordingly, the rotamerization of the thiol group in 2-thiocytosine was found to be reversibly carried out by pumping the amino (-NH 2 ) group [26]. The first remote switching using an NH group as a molecular antenna was done in the case of 6-methoxyindole, where the methoxy group could be selectively induced to rotate around the C-O bond [27]. Another molecule that greatly exemplifies this phenomenon is kojic acid, whose hydroxymethyl (-CH 2 OH) moiety goes through rotamerization upon exciting the -OH group on the other side of the ring of the molecule [28]. The OH or NH 2 groups of 3-hydroxy-2-formyl-2H-azirine and 3-amino-2-formyl-2H-azirine, respectively, also serve as molecular antennas that facilitate conformational control over the heavy aldehyde moiety three bonds away [29]. The current 'record holder' system, investigated by matrix-isolation IR spectroscopy, is E-glutaconic acid, in which long-range IVR could act over eight covalent bonds to successfully induce the conformational change of an OH group [30].
In this work, we aimed to investigate the selective conformational switching of a matrix-isolated 2-fluoro-4-hydroxy derivative of benzoic acid (2-F-4-OH benzoic acid in short) achieved by narrowband near-IR laser irradiation. This molecule serves as a great example, in which both 'local' and 'remote' switching can be done by exciting the same type of vibration (i.e., OH overtone) selectively, and their efficiencies can be directly compared with each other. Preliminary experiments showed that upon the excitation of the 2ν(OH) stretching vibrational overtone in a hydroxy carboxylic acid, change in the geometry of the COOH group in another part of the molecule can be efficiently induced [31]. In other words, higher energy cis carboxylic acid can be generated from the more stable trans form by remote switching (using the notation of Pettersson et al. for the two types of carboxylic acids) [1]. Moreover, the cis isomer is assumed to slowly convert back to trans over time via tunneling, which can also be monitored by IR spectrometry. The experiments were supplemented by quantum-chemical computations, which made the spectral analysis easier by allowing for simple comparison between the simulated spectra and the experimental ones. It is important to emphasize that effective remote switching that acts over multiple bonds has been observed in only a handful of molecular systems. As such, in order to fully understand its mechanism, it is necessary to find more examples of this phenomenon and study its overall efficiency.

MI-IR Experiments
A commercial 2-F-4-OH benzoic acid (≤100%) sample obtained from Sigma-Aldrich was used without further purification. The compound was inserted into a Knudsen cell with a heatable cartridge in a small quartz sample holder directly connected to the vacuum chamber. The adsorbed water and other volatiles had been removed the day before the experiments by slightly heating the sample compartment to 355 K and keeping it at that temperature for a couple of hours. A closed-cycle helium refrigeration system (Janis CCS-350R cold head cooled by a CTI Cryogenics 22 refrigerator) was used to cool down the CsI optical window to 13 K, which serves as the deposition temperature. The temperature of the cold window was measured by a silicon diode sensor connected to a LakeShore 321 digital temperature controller; the base pressure in the cooled chamber is usually found in the high 10 −9 mbar region. A sample sublimation temperature of 375 ± 2 K was used during the experiment whereupon the vapors were mixed and co-deposited in a 1:1000 ratio with Ar (Messer, 99.9999%) on the window. The orientation of the window (at a relative angle of 45 • with regard to the sample oven and the spectrometer beam) allows for a simultaneous deposition and spectral collection. For the latter, a Bruker IFS 55 FT-IR spectrometer equipped with an MCT detector cooled with liquid nitrogen was used. The transmission mid-IR spectra were taken during and after the deposition by averaging 16 and 128 scans, respectively, in the 2000-600 cm −1 region with 1 cm −1 resolution using KBr as a beam splitter and a low-pass filter with a cutoff wavenumber of 1830 cm −1 . In order to define the 2ν(OH) stretching overtones of each conformer necessary for the laser irradiation experiments, preliminary near-IR spectra were collected with the same instrument, without the cutoff filter, utilizing a W lamp and a CaF 2 beam splitter in the spectral region of the 8000-2500 cm −1 region at a resolution of 1 cm −1 .

Near-IR Laser Irradiation
The deposited matrices were in situ irradiated through the outer KBr window of the cryostat, applying a tunable narrowband laser light provided by the idler beam of an Optical Parametric Oscillator (OPO; GWU/Spectra-Physics VersaScan MB 240, fwhm ≈ 5 cm −1 ). The OPO was pumped by a pulsed Nd:YAG laser (Spectra-Physics Quanta Ray Lab 150, p ≈ 2.1-2.2 W, λ = 355 nm, f = 10 Hz, duration = 2-3 ns). The laser coming from the OPO device is perpendicular to the beam of the IR spectrometer, whereas the CsI cold window has a relative angle of 45 • to both of them, allowing for an online spectral collection during laser irradiation. The wavelength of the laser beam was selected so that the 2ν(OH) modes of the conformers of the matrix-isolated molecule were selectively excited. Accordingly, the sample was irradiated at 6952.0 cm −1 , 6994.9 cm −1 , 7077.1 cm −1 , and at 7093.2 cm −1 for approximately 60 min in each case; the photon fluxes were measured to be (3.1 ± 0.6) × 10 17 cm −2 s −1 (6952.0 cm −1 ), (3.2 ± 0.1) × 10 17 cm −2 s −1 (6994.9 cm −1 ), (3.3 ± 0.1) × 10 17 cm −2 s −1 (7077.1 cm −1 ), and (3.6 ± 0.1) × 10 17 cm −2 s −1 (7093.2 cm −1 ), respectively. Mid-IR spectra were collected online in the meantime, averaging 16 scans (during the 6952.0 cm −1 irradiation) and 12 scans (during the rest of the irradiations), respectively; 128 scans were taken before and after irradiation. After the irradiation experiments, the sample was left in the dark overnight while continuously collecting mid-IR spectra, averaging 128 scans every 15 min.

Theoretical Computations
The quantum chemical geometry optimizations, as well as the harmonic and anharmonic frequency computations, were carried out in Gaussian 09 (Rev. D01) with Becke's three-parameter hybrid functional B3LYP, using the non-local and local exchange correlation functionals as described by Lee-Yang-Parr and Vosko-Wilk-Nusair III, respectively [32][33][34][35]. Dunning's correlation-consistent cc-pVTZ basis set was used [36]. The anharmonic frequencies were used to simulate the vibrational spectra by convoluting them with Lorentzian functions with an FWHM of 0.5 cm −1 in the Synspec software [37]. The isomerization barriers were estimated from the optimized transition states (TSs) using the Berny algorithm, as implemented in Gaussian [38]. This was followed by the computation of intrinsic reaction paths (IRPs) in the Cartesian coordinates [39]. The tunneling rates were estimated using the Wentzel-Kramers-Brillouin (WKB) model described in detail elsewhere [40].

Structure and Energy of the Conformers
The conformational structure of the molecule is primarily defined by its two functional groups, the 4-hydroxyl (-OH) and the carboxylic group (-COOH). The H atom on both moieties may have two orientations (the respective torsional angles are denoted by ϕ 1 and ϕ 3 ), whereas -COOH has two forms based on the alignment of its carbonyl (C=O) group (ϕ 2 ). This means 2 3 = 8 different conformers, which are depicted in Scheme 1, in which the torsional angles most important to the molecular structure are also marked. The conformers with ϕ 3 close to 180 • are named A (ϕ 1 ≈ 180 • ) or B (ϕ 1 ≈ 0 • ), respectively. The same applies to structures C and D, which can have a ϕ 3 value of 0 • with ϕ 1 ≈ 0 • or 180 • , respectively, whereas the labels 1 or 2 after them differentiate the ones with ϕ 2 ≈ 180 • and 0 • , respectively. This means that the following notations for the eight conformers are used throughout the manuscript: A1, A2, B1, B2, C1, C2, D1 and D2; their optimized geometries are listed Tables S1-S20 in the Supplementary Material.
According to the orientation of the groups discussed above, the conformers have different stabilities. The most stable conformer is A1, which is closely followed by B1; the only difference between them is the alignment of the 4-OH group. Apart from this, they both have a trans-COOH carboxylic group (i.e., there is a weak intramolecular H bond between the H atom of the OH group with the C=O oxygen atom), whereas the O atom in the C=O group is on the opposite side with regard to the F atom. It should be noted that the energy difference is predicted to be 0.8 kJ mol −1 , which is within the accuracy of the computational level. The next pair of conformers have slightly higher relative energies (by ca. 3 kJ mol −1 ) denoted with A2 and B2. They, similarly to A1 and B1, only differ from each other in the orientation of the 4-OH hydrogen atom as well, but, unlike A1 and B1, their C=O oxygen atoms are on the same side as the F atom. The next four conformers have cis-COOH groups, meaning inherently higher relative energies due to breakage of the abovementioned intramolecular H bond, which is counterbalanced to some extent in the case of the C1 and D1 by the formation of a less energetic H bond with the F atom. Their relative energy (with regard to the A1 form) is approximately 8 kJ mol −1 ; this value should more or less reflect on the difference of the H atom bond strength between the cases when the H atom creates the bond with the C=O oxygen and when it bonds with the F atom. Nevertheless, C1 and D1 only differ in the orientation of the 4-OH hydrogen atom, whereas their C=O oxygen atoms (similarly to A1 and B1) are on the opposite side of the molecule with regard to the F atom, thus allowing for H bonding between the carboxylic H atom and the F atom. Lastly, C2 and D2 are unique in the sense that, due to steric effects, their -COOH group is not in the same plane as the aromatic ring, which means that the ϕ 2 torsional angle is closer to 30 • instead of 0 • and can be both positive and negative (i.e., the C=O oxygen atom and the carboxylic OH group are slightly above or below the plane of the aromatic ring). This results in the fact that both C2 and D2 are actually chiral and have two enantiomeric forms differentiated by an asterisk in Table 1 (i.e., C2, C2*, D2, D2*), and the only difference between C2 and D2 is the orientation of the 4-OH group. Apart from this, their C=O oxygen atoms, similarly to A2 and B2, is on the same side of the molecule with regard to the F atom and have cis-COOH groups (similarly to C1 and D1). As a matter of fact, this causes distortion of the molecule as the COOH hydrogen atom becomes relatively close to the H atom on the C atom in position six of the ring. Due to the unfavorable effects listed above, there is no intramolecular H bond in them, and therefore, they are much less stable than the other conformers; their relative energy lies more than 30 kJ mol −1 above the most stable A1 form. This also implies that their Boltzmann population at the sample inlet temperature (375 ± 2 K) is expected to be significantly less than 1%, which prohibits their IR detection in the deposited matrix.

Structure and Energy of the Conformers
The conformational structure of the molecule is primarily defined by its two functional groups, the 4-hydroxyl (-OH) and the carboxylic group (-COOH). The H atom on both moieties may have two orientations (the respective torsional angles are denoted by φ1 and φ3), whereas -COOH has two forms based on the alignment of its carbonyl (C=O) group (φ2). This means 2 3 = 8 different conformers, which are depicted in Scheme 1, in which the torsional angles most important to the molecular structure are also marked. The conformers with φ3 close to 180° are named A (φ1 ≈ 180°) or B (φ1 ≈ 0°), respectively. The same applies to structures C and D, which can have a φ3 value of 0° with φ1 ≈ 0° or 180°, respectively, whereas the labels 1 or 2 after them differentiate the ones with φ2 ≈ 180° and 0°, respectively. This means that the following notations for the eight conformers are used throughout the manuscript: A1, A2, B1, B2, C1, C2, D1 and D2; their optimized geometries are listed Tables S1-S20 in the Supplementary Material. According to the orientation of the groups discussed above, the conformers have different stabilities. The most stable conformer is A1, which is closely followed by B1; the only difference between them is the alignment of the 4-OH group. Apart from this, they both have a trans-COOH carboxylic group (i.e., there is a weak intramolecular H bond between the H atom of the OH group with the C=O oxygen atom), whereas the O atom in the C=O group is on the opposite side with regard to the F atom. It should be noted that Scheme 1. Structures of the eight conformers of 2-F-4-OH benzoic acid, namely A1, A2, B1, B2, C1, C2, D1, and D2 as denoted by the labels in bold and in parentheses below the structures. The three torsional angles ϕ 1 , ϕ 2 and ϕ 3 are defined in structure A1. In case of C2 and D2, only one of the two enantiomeric forms are depicted. Table 1. Torsional angles (in degrees) and anharmonic zero-point energy-corrected relative energy values (kJ mol -1 ) of the 2-F-4-OH-benzoic acid conformers. The relative energies of the TSs, along with their optimized geometries, are summarized in Tables S21-S49 in the Supplementary Material. Based on the values listed there, the rotation of the 4-OH group along the ϕ 1 torsional angle has a barrier of approximately 20 kJ mol −1 (A1↔B1, A2↔B2, C1↔D1, C2↔D2), and the same holds true for the rotation along ϕ 2 in the case of A1↔A2 and B1↔B2, but not for C1↔C2 and D1↔D2 where the barriers of the 'rightward' direction (C1→C2 and D1→D2, ≈30 kJ mol −1 ) are much higher than those in the reverse one (C1←C2, D1←D2, <10 kJ mol −1 ) due to the energy difference of the minima A1, B1 and A2, B2, respectively. Unsurprisingly, the barrier of the rotation around ϕ 3 is the highest one with values lying between 40-50 kJ mol −1 (A1↔D1, B1↔C1), except if one of the conformers is much more stable than the other one. In these cases (A2↔D2, B2↔C2), the barrier heights are about 50 and 20 kJ mol −1 , respectively, depending on the direction of the transition. It should be noted that all the rotations around the three torsional angles may go either clockwise or counterclockwise, but the TSs produced by the two different rotational directions are mirror images of each other and, therefore, have identical energy and vibrational spectra. It is also worth mentioning that the TSs between the enantiomers C2↔C2* and D2↔D2* have a planar structure, and their energy is only slightly above the minima (≈5 kJ mol −1 , approximately 400 in cm −1 ), which allows for their rapid interconversion upon exposure to the spectrometer beam source.

Changes upon Near-IR Irradiation, Vibrational Analysis
Some selected regions of the near-and mid-IR spectra of the sample molecule deposited in an Ar matrix are visualized in Figures 1a, 2a, 3a and 4a. Furthermore, by having a look at the mid-IR spectra in Figures 1b, 2c, 3c and 4c, one can see that irradiating some 2ν(OH) stretching overtone bands results in a selective conformational switch. The spectral features can be divided into three different groups based on their general behavior upon irradiation: (a) the ones that increase, (b) the ones that decrease, and (c) the ones that do not change, which should belong to three different conformers or groups of conformers. One should also bear in mind that the enantiomers C2, C2* and D2, D2* cannot be differentiated by their IR spectrum. Furthermore, based on theoretical results, the conformers that differ only in the orientation of their 4-OH groups have almost identical vibrational frequencies; thus, they cannot be distinguished either (Tables S11-S20). This leads us to the conclusion that only four band groups should be detected in the matrix belonging to the conformer pairs (1) A1/B1, (2) A2/B2, (3) C1/D1, and (4) C2/D2 (or C2*/D2*). However, as mentioned above, only three groups of bands could be distinguished based on their response to laser irradiation. As such, the conformer pairs with the highest energy (i.e., group (4)) could be excluded based on the following considerations: First, they are 30 kJ mol −1 above the most stable members of group (1), meaning that their Boltzmann population at the deposition temperature is expected to be 0.003%. For comparison, the values for the other groups are 72% for group (1), 23% for group (2), and 5% for group (3). It is important to note that the abundance of conformer groups in the freshly deposited matrix is very similar to these theoretically obtained ones and can be determined to be 61% for group (1), 33% for group (2), and 6% for group (3). Consequently, group (4) is not expected in the freshly deposited matrix. Secondly, no new bands arise during any of the near-IR excitation (only the relative intensities of the bands that are already present change), suggesting that members of this group are not generated upon irradiation. Furthermore, it can be assumed that only the most stable members of each of the three groups (i.e., A1, A2, and D1) are expected to be present in the matrix, since the higher energy members can be converted to them via tunneling. This hypothesis will be justified in Section 3.5 where an estimate on the tunneling rates are given. The bands belonging to A1 all decrease when its 2ν(OH) stretching overtone bands, found at 6952.0 cm −1 and 7077.1 cm −1 , are excited. This occurs with the simultaneous increase in the bands belonging to D1, i.e., the A1→D1 conversion is induced. Similarly, when the matrix is irradiated at 6994.9 cm −1 and at 7093.2 cm −1 , the opposite can be observed, i.e., the D1→A1 conversion; therefore these bands must belong to the 2ν(OH) stretching overtones of the former conformer. It is important to note the findings described above show that no matter which OH group is actually excited (i.e., the 4-OH or the carboxylic one), all near-IR irradiation results in the same rotamerization process, which is the cis-trans rotamerization of the carboxylic moiety. As a consequence, the irradiation of the carboxylic OH group (6952.0 cm −1 for A1 and 6994.9 cm −1 for D1) will be called 'local' excitation, whereas that of the 4-OH group (7077.1 cm −1 for A1 and 7093.2 cm −1 for D1) will be called 'remote' excitation hereafter.  In contrast to the results discussed above, the bands of A2 are not affected by any of the irradiations, even if an excitation wavelength of 6978.8 cm −1 is used, which is the frequency of its local carboxylic 2ν(OH) stretching overtone. Accordingly, Table 2 sums up the general response of these three conformers to the excitation of various O-H stretching overtone bands. The different response to near-IR laser irradiation, accompanied by a comparison of the experimental and theoretical spectra, allows for assignment of the vibrational bands; Tables 3-5 lists their spectral assignment for part of the mid-IR spectral range. overtone bands. The different response to near-IR laser irradia comparison of the experimental and theoretical spectra, allows f brational bands; Tables 3-5 lists their spectral assignment for pa range.       a Normalized experimental band areas obtained by multiplying each integrated band area by the sum of the theoretical infrared intensities divided by the sum of the experimental band areas; b ν: stretching, β: in-plane bending, γ: out-of-plane bending, δ: in-plane deformation, τ: out-of-plane deformation, ?: unassigned peak.  a Normalized experimental band areas obtained by multiplying each integrated band area by the sum of the theoretical infrared intensities divided by the sum of the experimental band areas; b ν: stretching, β: in-plane bending, γ: out-of-plane bending, δ: in-plane deformation, τ: out-of-plane deformation, ?: unassigned peak. a Normalized experimental band areas obtained by multiplying each integrated band area by the sum of the theoretical infrared intensities divided by the sum of the experimental band areas; b ν: stretching, β: in-plane bending, γ: out-of-plane bending, δ: in-plane deformation, τ: out-of-plane deformation, ?: unassigned peak.

Kinetics of the Near-IR Induced Rotamerization, Quantum Efficiencies
The temporal evolution of the spectrum can be monitored throughout the near-IR irradiation, which provides us with a quantitative means of obtaining the conversion rate. The kinetic curves are plotted in Figure 5 for the 1742.5 and 1628.2 cm −1 bands for conformers A1 and D1, respectively. A single exponential function was used to fit the decay profiles: The following equation could be used to fit the growth profiles: Moreover, the quantum efficiencies can be estimated using the form [3,6,41] The following equation could be used to fit the growth profiles: Moreover, the quantum efficiencies can be estimated using the formula given below [3,6,41] φ(i) = k(ṽ)/σ i (ṽ)I(ṽ) (3) where φ(i) is the quantum yield of the process when exciting the ith vibrational mode at the overtone wavenumberṽ; k(ṽ) is the rotamerization rate in s −1 , whereas σ i (ṽ) is the absorption cross-section of that particular vibrational mode in cm 2 ; and I(ṽ) is the average photon intensity of the laser beam at that wavelength in cm 2 s −1 . We can obtain σ i (ṽ) in cm 2 when the absorbance of the excited overtone band (A, dimensionless value) is divided by the calculated column density of the decaying conformer at the beginning of irradiation (N, in cm −2 ). Furthermore, the latter can be deduced from the integrated peak area of one absorption band of the conformer (A int , in cm) and by knowing the absorption coefficient (α, in cm) and the area of irradiation (S = 1 cm 2 ), respectively: Even though α is not known for our system, a rough estimate can be given based on the computed anharmonic IR intensities (Tables 3 and 4). In order to determine N, the strong and well-resolved ν(C=O) stretching vibrational modes were used. Moreover, I(ṽ) can be obtained if the output power of the produced near-IR laser light (P, in W) is divided by the photon energy (E photon , in J) times the surface area S. It should be noted that P was measured without the KBr window in the beam path; thus, P and I(ṽ) both represent an upper estimate, which also means that the quantum efficiencies deduced here are lower estimates. Table 6 lists all the derived values. Table 6. Rotamerization rates (k X , in sec −1 , where X = A1 or D1, respectively) of the conformers upon near-IR irradiation. Positive values show growth, whereas negative ones indicate decay. Computed quantum yields for the rotamerization of the carboxylic OH group when exciting the OH stretching overtones.
As an alternative approach, φ(i) can also be estimated as follows: φ(i) = N iso /N abs (5) where N iso denotes the number of molecules converted, and N abs is equal to the number of photons absorbed per time unit (both in s −1 ).
The averaged value φ(i) average varies between 1.4 × 10 −1 and 7 × 10 −2 (Table 6), which agree well with, for instance, those of the formic acid (1.7 × 10 −1 , 7 × 10 −2 ) [3], acetic acid (2.2 × 10 −2 ) [6], and propionic acid (1.4 × 10 −2 ) [41], respectively, but they are significantly higher than those of amino acids, such as glycine (8 × 10 −4 ) [18] and alanine (5 × 10 −4 and 1 × 10 −3 ), respectively [19]. It is also worth noting that the relative uncertainty of φ(i) average is found to be around 30-40%, which is also in agreement with previous findings [6,41]. By having a look at the φ(i) average values belonging to the different excitation frequencies, it can be deduced that those of the 'local' processes do not differ significantly from each other. Furthermore, the φ(i) average values of the 'remote' processes are comparable to the 'local' ones, which is somewhat unexpected, since it would be straightforward to think that energy dissipation is not negligible when exciting a distant functional group that results in lower efficiencies for the 'remote' excitations. Nevertheless, this does not necessarily hold true because, for instance, thioacetamide, which also shows a remote molecular switching property upon near-IR irradiation, has φ(i) average values of 3.7 × 10 −2 to 7.2 × 10 −2 [42]; these are not significantly lower than the quantum efficiencies listed above. It has to be kept in mind, however, that the IVR process in thioacetamide acts through only four bonds instead of the eight in 2-F-4-OH benzoic acid, thus the current findings might still be surprising. A possible explanation for the minor difference between the quantum efficiencies of the 'local' and 'remote' excitations might be the rapid redistribution of the vibrational energy within the molecule, which has a much higher rate than that of the energy dissipation into the surrounding matrix. In this case, the latter would be the one that determines the overall efficiency, i.e., it would matter less which vibrational mode of the molecule is actually excited.

Tunneling Decay Kinetics
After the irradiation experiments, the matrix was irradiated once more with the 6952.0 cm −1 laser light to bleach conformer A1 and generate as much D1 as possible. Then, the sample was left in the dark overnight while continuously collecting the mid-IR spectra. Furthermore, it is important to recall that a low-pass filter with a cutoff wavelength of 1830 cm −1 was installed between the spectrometer and the sample during the experiment. This was done in order to prevent the processes induced by the broadband IR light originating from the spectrometer beam source, which is known to facilitate unwanted rotamerization processes [43][44][45]. It should be noted that the photon energy at the filter cutoff wavelength is 21.9 kJ mol −1 , which is significantly below the A1←D1 rotamerization barrier (36.7 kJ mol −1 , Table S49). However, this cutoff energy is comparable to, or even higher than, the barrier of processes involving the change of other torsional angles. For instance, the 4-OH group almost rotates freely around, and the rotation of the carboxylic group should also be made possible upon exposure to the IR beam. It is also important to note that the 36.7 kJ mol −1 barrier (which equals 3068 cm −1 ) should be easily overcome without the presence of the filter by exciting the vibrational modes above this threshold, such as the O-H stretching vibrations; this finding justifies the use of the filter. Figure 6 shows the kinetic growth/decay of conformers A1 and D1, respectively, when the already irradiated matrix was left in the dark overnight. To fit the experimental data, the same single exponential functions defined by equations E1 and E2 were used as in Section 3.3. The k values, as expected, are significantly (by roughly 2 orders of magnitude) lower than those obtained for the vibrationally induced processes and have a fair agreement with each other. The half-life of D1 was found to be around 115,000 ± 4000 s, which is some 32 h. It is interesting to note, however, that the second-order exponentials could also be fitted on these data providing a somewhat improved fit. The same holds true for the irradiation-induced processes, in most cases. This finding implies that there may be matrix sites with different vibrationally induced rotamerization and tunneling rates; thus, they can be classified as 'slow' and 'fast', which is a phenomenon that has been previously described [10,23,46].

Estimating the Tunneling Half-Life
The tunneling rates of the various feasible pathways can be estimated based on computational considerations, which can be used to interpret the experimental findings. The intrinsic reaction paths can be obtained after running IRC computations in Gaussian, which allows for estimation of the width and the height of the barriers that fundamentally determine the tunneling half-lives of the processes. Figures S1-S4 show the computed IRC profiles between the isomers, which also contain the barrier heights and widths used in the calculations. The WKB method mentioned in Section 2.3 was used to predict the tunneling rates using the equation below: where P(E) is the probability of the tunneling process, m stands for the particle mass (1.68 × 10 -27 kg for H), V0 is the barrier height, E is the particle energy (both in J), h is Planck's constant, and w is the barrier width (in m). The tunneling rate k (in s -1 ) can be obtained by multiplying P(E) by the frequency of attempts (n, in s -1 ), which in the latter can be estimated using the anharmonic frequency of the vibrational mode that takes place in the process (ν, in m -1 ) by means of the formula ν = νc, where c is the speed of light (299,792,458 m s -1 ). The half-life t1/2 (in s -1 ) can be derived from k using the equation t1/2 = ln2/k. It is important to note that, due to the inherent uncertainty of the computational results, it is not extraordinary to have a difference of orders of magnitude between the computational and experimental values. What really matters is how the calculated values compare to

Estimating the Tunneling Half-Life
The tunneling rates of the various feasible pathways can be estimated based on computational considerations, which can be used to interpret the experimental findings. The intrinsic reaction paths can be obtained after running IRC computations in Gaussian, which allows for estimation of the width and the height of the barriers that fundamentally determine the tunneling half-lives of the processes. Figures S1-S4 show the computed IRC profiles between the isomers, which also contain the barrier heights and widths used in the calculations. The WKB method mentioned in Section 2.3 was used to predict the tunneling rates using the equation below: where P(E) is the probability of the tunneling process, m stands for the particle mass (1.68 × 10 -27 kg for H), V 0 is the barrier height, E is the particle energy (both in J), h is Planck's constant, and w is the barrier width (in m). The tunneling rate k (in s -1 ) can be obtained by multiplying P(E) by the frequency of attempts (n, in s -1 ), which in the latter can be estimated using the anharmonic frequency of the vibrational mode that takes place in the process (ṽ, in m -1 ) by means of the formula ν =ṽc, where c is the speed of light (299,792,458 m s -1 ). The half-life t 1/2 (in s -1 ) can be derived from k using the equation t 1/2 = ln2/k. It is important to note that, due to the inherent uncertainty of the computational results, it is not extraordinary to have a difference of orders of magnitude between the computational and experimental values. What really matters is how the calculated values compare to each other, and this is why the theoretically obtained values can be scaled up with a uniform factor, so that they will be equal to their experimental counterparts [47]. Here, only the D1→A1 process that can be observed, so the theoretical value was multiplied by a scaling factor of 10,500, and the same was done for all of the other processes as well. Scheme 2 visualizes the tunneling rates of all possible processes, including the rotation of the 4-OH group, the OH group in the COOH moiety, and that of the whole COOH group. each other, and this is why the theoretically obtained values can be scaled up with a uniform factor, so that they will be equal to their experimental counterparts [47]. Here, only the D1→A1 process that can be observed, so the theoretical value was multiplied by a scaling factor of 10,500, and the same was done for all of the other processes as well. Scheme 2 visualizes the tunneling rates of all possible processes, including the rotation of the 4-OH group, the OH group in the COOH moiety, and that of the whole COOH group.

Scheme 2.
Calculated tunneling half-lives for all conformers. The values are scaled with a uniform factor of 10,500, so that the half-life of the D1→A1 process is identical to that of the experimentally obtained one. The fastest processes for each conformer are highlighted in bold, whereas the labels of the conformers that can be experimentally observed are outlined.
What can be deduced from Scheme 2 is that all of the higher energy conformers (C1, C2, D1, D2), even if they are in the sample in detectable amounts, quickly convert to lower energy ones, such as A2, B2, or D1. For instance, if there is any C1 in the matrix at the beginning, most of it should be converted to D1 by the end of the deposition, and the other processes are even faster than that one. B2 converts to A2 in less than a minute, whereas D1 is shown to slowly transform into A1. It is important to note that neither A2 nor B2 may convert to their counterparts A1 or B1 via tunneling over a reasonable timescale and that the conversion rate of D2→A2 is some seven orders of magnitude faster than that of D1→A1, which may explain why the former process cannot be observed experimentally. Most importantly, the scheme also nicely shows why only three forms may be expected to be present in the matrix after deposition, namely A1, A2, and D1, which is in accordance with the experimental results. What can be deduced from Scheme 2 is that all of the higher energy conformers (C1, C2, D1, D2), even if they are in the sample in detectable amounts, quickly convert to lower energy ones, such as A2, B2, or D1. For instance, if there is any C1 in the matrix at the beginning, most of it should be converted to D1 by the end of the deposition, and the other processes are even faster than that one. B2 converts to A2 in less than a minute, whereas D1 is shown to slowly transform into A1. It is important to note that neither A2 nor B2 may convert to their counterparts A1 or B1 via tunneling over a reasonable timescale and that the conversion rate of D2→A2 is some seven orders of magnitude faster than that of D1→A1, which may explain why the former process cannot be observed experimentally. Most importantly, the scheme also nicely shows why only three forms may be expected to be present in the matrix after deposition, namely A1, A2, and D1, which is in accordance with the experimental results.

Conclusions
In this study, 2-F-4-OH benzoic acid was isolated in Ar matrices at low temperatures while collecting its mid-IR spectra. The molecule has eight different conformers based on the position of three fundamental torsional angles, which are the orientation of the H atom in the 4-OH group and that of the O and H atoms in the -COOH group at the other end of the aromatic ring. The spectral assignment can be made by comparing the experimental spectrum with the theoretically predicted spectra of the conformers. Based on the vibrational analysis, only the three most stable conformers (A1, A2, and D1) are expected to be present in the matrix after deposition, which are supported by their calculated Boltzmann populations at the deposition temperature. When exciting the 2ν(OH) stretching overtones of the conformers, it can be deduced that A1 converts to D1 and vice versa, meaning that rotation of the H atom on the COOH group is induced. In contrast to this, no similar effect could be observed for conformer A2. The response of A1 or D1 depends only marginally on whether the irradiation occurred 'locally' (i.e., by the excitation of the OH of the COOH group) or 'remotely' (by the excitation of the 4-OH moiety). The behavior of the bands upon excitation further confirms the vibrational assignment. The kinetic rates of the rotamerization, as well as the quantum efficiencies of the 'local' and 'remote' excitations, were also determined. According to this, the efficiency is similar in all cases independently on the excitation wavelength. The spontaneous conversion of the higher energy conformer D1 to the more stable A1 via tunneling was also examined, and, as such, its rate was determined and found to be roughly two orders of magnitude slower than that of the vibrationally induced one.
This molecule represents a great example, where both the local and remote switching could be studied simultaneously, thus allowing for their direct comparison. Therefore, the molecular system presented here further extends our understanding of intramolecular vibrational energy transfer processes.