Electrochemical SEIRAS Analysis of Imidazole-Ring-Functionalized Self-Assembled Monolayers

An essential amino acid, histidine, has a vital role in the secondary structure and catalytic activity of proteins because of the diverse interactions its side chain imidazole (Im) ring can take part in. Among these interactions, hydrogen donating and accepting bonding are often found to operate at the charged interfaces. However, despite the great biological significance, hydrogen-bond interactions are difficult to investigate at electrochemical interfaces due to the lack of appropriate experimental methods. Here, we present a surface-enhanced infrared absorption spectroscopy (SEIRAS) and density functional theory (DFT) study addressing this issue. To probe the hydrogen-bond interactions of the Im at the electrified organic layer/water interface, we constructed Au-adsorbed self-assembled monolayers (SAMs) that are functionalized with the Im group. As the prerequisite for spectroelectrochemical investigations, we first analyzed the formation of the monolayer and the relationship between the chemical composition of SAM and its structure. Infrared absorption markers that are sensitive to hydrogen-bonding interactions were identified. We found that negative electrode polarization effectively reduced hydrogen-bonding strength at the Im ring at the organic layer–water interface. The possible mechanism governing such a decrease in hydrogen-bonding interaction strength is discussed.


Introduction
Histidine involvement in protein architecture and catalytic activity is primarily determined by the diverse range of interactions that its side chain imidazole (Im) can take part in [1]. Im acts as an aromatic system, participates in hydrogen-bond (H-bond) interactions, and coordinates metal (II) cations through the nitrogen's lone electron pair [2]. Im has a pKa value of~5.9; therefore, in physiological conditions, only one nitrogen is protonated, which results in two possible tautomeric forms. Tautomer-I (T-I) is protonated at N1 atom (N1-H, N3); Tautomer-II (T-II) is protonated at N3 (N1, N3-H). Tautomeric forms of imidazole ring are shown in Figure 1. Because of the slightly more energetically favorable form, T-I dominates neutral pH solutions at room temperature [3]. The tautomeric equilibrium can change depending on the chemical Im group environment, temperature, intermolecular interactions, and other external conditions. For example, a clear T-I to T-II shift in imidazole-containing compounds was found to be induced by the chemisorption on a metal surface [4].
Hydrogen bonding is among the most important interactions that govern protein folding, molecular conformation, and packing. Notable is the Im ring, which in neutral form can participate in both hydrogen accepting and donating interactions. However, studying these H-interactions at the electrified interface is a challenging task due to the lack of appropriate experimental techniques. Vibrational spectroscopy is highly sensitive Hydrogen bonding is among the most important interactions that govern protein folding, molecular conformation, and packing. Notable is the Im ring, which in neutral form can participate in both hydrogen accepting and donating interactions. However, studying these H-interactions at the electrified interface is a challenging task due to the lack of appropriate experimental techniques. Vibrational spectroscopy is highly sensitive to the local environment and structure of probed molecules, intermolecular forces, and electric fields. In particular, Raman spectroscopy provides detailed information on the Im structure and the interactions with divalent metal cations, both in natural biomolecules and model compounds [3][4][5][6][7][8][9][10][11][12]. Fourier-transform infrared absorption spectroscopy (FTIR) exceeds Raman by a many-orders-of-magnitude-larger process cross section, which is especially beneficial in dynamic studies and in the cases of a low number of molecules under investigation.
Surface-enhanced infrared absorption spectroscopy (SEIRAS), pioneered by Osawa [13], has become a valuable tool in the analysis of Au or Ag surface-adsorbed supramolecular systems in situ [14]. In recent years, the method became technologically simpler, which has even led to its commercialization [15,16]. Various techniques of thin Au-layer preparation were devised, such as electroless, vapor, and electrochemical deposition [17], and also hybrid ones [18], which offer increased signal-to-noise ratio, Au-film rigidness, and even enable multiple reuses of the same film. The most readily available method today remains electroless deposition since it is relatively fast and requires less investment in equipment. With the right amount of care taken due to harmful chemicals used, the method is not complicated and is even time-efficient [19].
In this study, we used an electrochemical-SEIRAS method to attain a deeper understanding of the Im ring's interactions with the local environment at the organic layerwater interface on an electrically polarized electrode. The Im-ring-functionalized self-assembled monolayers were formed on the Au surface from the N-(2-(1H-imidazol-4yl)ethyl)-6-mercaptohexanamide (IMHA) compound. The SEIRAS technique provided insights into the adsorption and reorientation dynamics of the SAM molecules with an emphasis on the changes in the functional Im ring group. Spectroscopic markers for the Hbond interactions at the Im group were identified and the analysis of their alterations at the electrified interface was provided. Spectroscopic findings were backed by quantum chemistry calculations, which are prominent when studying hydrogen-bonding interactions [20]. Our findings will shed some light on the interfacial behavior of the amino acid, histidine, which has an essential role in protein structure and activity and will be beneficial in surface chemistry applications. Surface-enhanced infrared absorption spectroscopy (SEIRAS), pioneered by Osawa [13], has become a valuable tool in the analysis of Au or Ag surface-adsorbed supramolecular systems in situ [14]. In recent years, the method became technologically simpler, which has even led to its commercialization [15,16]. Various techniques of thin Au-layer preparation were devised, such as electroless, vapor, and electrochemical deposition [17], and also hybrid ones [18], which offer increased signal-to-noise ratio, Au-film rigidness, and even enable multiple reuses of the same film. The most readily available method today remains electroless deposition since it is relatively fast and requires less investment in equipment. With the right amount of care taken due to harmful chemicals used, the method is not complicated and is even time-efficient [19].
In this study, we used an electrochemical-SEIRAS method to attain a deeper understanding of the Im ring's interactions with the local environment at the organic layer-water interface on an electrically polarized electrode. The Im-ring-functionalized self-assembled monolayers were formed on the Au surface from the N-(2-(1H-imidazol-4-yl)ethyl)-6mercaptohexanamide (IMHA) compound. The SEIRAS technique provided insights into the adsorption and reorientation dynamics of the SAM molecules with an emphasis on the changes in the functional Im ring group. Spectroscopic markers for the H-bond interactions at the Im group were identified and the analysis of their alterations at the electrified interface was provided. Spectroscopic findings were backed by quantum chemistry calculations, which are prominent when studying hydrogen-bonding interactions [20]. Our findings will shed some light on the interfacial behavior of the amino acid, histidine, which has an essential role in protein structure and activity and will be beneficial in surface chemistry applications.

Electroless Au Plating
The reflecting plane of the Si crystal was polished by hand using a polishing pad (QATM, Mammelzen, Germany) and water-based diamond suspensions starting with 3 µm (5 min), then 1 µm (5 min), and 0.25 µm (20 min). Polishing was performed by shaping the "∞" symbol and turning the crystal by 90 degrees every 15 s, with thorough water rinsing in between each polishing step. Then, the mirror-polished Si crystal was ultrasonicated with acetone and water twice, for 10 min each time, and dried under a nitrogen stream. The prism sides were protected with 0.1 mm Teflon ribbon to prevent coatings from damage. Fresh Au-plating mixture was prepared in a Teflon container and stirred continuously under a nitrogen blanket. For that, equal volumes of (1) salt solution (0.15 M Na 2 SO 3 , 0.05 M Na 2 S 2 O 3 , 0.05 M NH 4 Cl), (2) NH 4 F (20 wt%), (3) HF (2 wt%), and (4) NaAuCl 4 (0.03 M) were mixed to at least 2 mL in total. The (4) solution was freshly prepared for every plating. To hydrogenate the Si surface, the polished crystal was incubated with 1 mL HF (2 wt%) for 3 min, after which it was removed and 1 mL Auplating mixture was added for another 3-4 min in room-temperature conditions. During Au plating, the crystal surface attained a bright yellow color. Reaction was stopped by transferring the crystal directly into a water bath for 15 min. SEM image in Figure 2 shows the formation of homogeneous nanostructures of approximately 50 nm in diameter.

Electroless Au Plating
The reflecting plane of the Si crystal was polished by hand using a polishing pad (QATM, Germany) and water-based diamond suspensions starting with 3 µ m (5 min), then 1 µ m (5 min), and 0.25 µ m (20 min). Polishing was performed by shaping the "∞" symbol and turning the crystal by 90 degrees every 15 s, with thorough water rinsing in between each polishing step. Then, the mirror-polished Si crystal was ultrasonicated with acetone and water twice, for 10 min each time, and dried under a nitrogen stream. The prism sides were protected with 0.1 mm Teflon ribbon to prevent coatings from damage. Fresh Au-plating mixture was prepared in a Teflon container and stirred continuously under a nitrogen blanket. For that, equal volumes of (1) salt solution (0.15 M Na2SO3, 0.05 M Na2S2O3, 0.05 M NH4Cl), (2) NH4F (20 wt%), (3) HF (2 wt%), and (4) NaAuCl4 (0.03 M) were mixed to at least 2 mL in total. The (4) solution was freshly prepared for every plating. To hydrogenate the Si surface, the polished crystal was incubated with 1 mL HF (2 wt%) for 3 min, after which it was removed and 1 mL Au-plating mixture was added for another 3-4 min in room-temperature conditions. During Au plating, the crystal surface attained a bright yellow color. Reaction was stopped by transferring the crystal directly into a water bath for 15 min. SEM image in Figure 2 shows the formation of homogeneous nanostructures of approximately 50 nm in diameter.

Attenuated Total Reflection (ATR)-SEIRAS Setup and Au-Film Activation
Spectral measurements were carried out by using Vertex 80v spectrometer (Bruker, Ettlingen, Germany) equipped with a liquid-nitrogen-cooled narrow-band MCT (HgCdTe) detector working at a 40 kHz scanner speed. The resolution was set to 4 cm −1 and aperture to 2 mm; sample and background scans were collected from 50 and 100 iterations, accordingly. Freshly prepared nitrogen-blow-dried Au/Si crystal was assembled into a VeeMax III variable angle accessory with Jackfish cell J1F (Pike Technologies, Fitchburg, WI, USA). The angle of the ATR unit was set to 63 degrees. After assembly, the spectrometer was purged with dry air overnight to eliminate residual water vapor from the spectrometer chambers.
Gold layer of the Si crystal was electropolished and additionally activated in pH 5. by PGSTAT101 potentiostat (Methrom, Riverview, FL, USA). A three-electrode system comprised the gold layer as working electrodes, reference Ag/AgCl, and counter Pt electrodes. A cell filled with sodium acetate solution (0.1 M) was purged with nitrogen for 30 min and CV was performed, starting in the range of ±200 mV from open circuit potential (OCP) at 20 mV/s speed. After 3 full cycles, potential range was increased by 100 mV towards the anodic direction. The procedure continued with 3 time repeats for each 100 mV increment until gold oxidation became detectable. SEIRAS measurements at designated CV steps were completed to confirm surface cleaning. Reference (sample) spectra were collected at -100 mV (600 mV) to identify the marker band of the gold surface-adsorbed acetate. The υν(COO -) mode near 1400 cm −1 gained intensity with increasing surface cleanness at later steps of the CV treatment. The couple of milli-absorption intensities of the ν(COO -) mode is sought after since longer CV treatment increases the risk of gold-layer delamination.

Self-Assembled Monolayer (SAM) Formation
Gold layer was washed several times with H 2 O before collecting the reference SEIRAS spectrum for later use. Then, the cell was thoroughly rinsed with ethanol without letting the gold surface to dry. Approximately 0.2 mL of IMHA, Frag, or IMHA/Frag ethanol mixture ( Figure 1) was injected into a cell to a total 0.5 mM concentration and incubated for 90 min. Monolayers were formed from 100, 50, 20, and 0 mol% IMHA solutions and denoted as IMHA, IMHA/Frag (1/1), IMHA/Frag (1/5), and Frag. Each experiment was repeated in the same order by using the deuterated solvents (D 2 O and ethanol-d 6 ).

Attenuated Total Reflection Spectroscopy (ATR-FTIR) Measurements
ATR-FTIR spectra of IMHA and its fragment were measured using Alpha spectrometer (Bruker, Ettlingen, Germany) equipped with a diamond ATR accessory and DTGS detector working at a 7.5 kHz scan speed. The spectral resolution was set to 4 cm −1 and 200 scans were co-added for sample and background channels.

Scanning Electron Microscopy (SEM)
The SEM images were acquired by JSM-IT200 InTouchScope TM (Jeol, Kyoto, Japan) microscope working at 20 kV accelerating voltage and a secondary electron detector.

Density Functional Theory (DFT) Modeling
Theoretical modeling was performed using Gaussian 09W [22]. Geometry optimization and frequency calculations were completed by using the B3LYP-functional and the 6-311++G(2d,p)-basis sets without anharmonic correction. Calculated frequencies and intensities were scaled according to the procedure described elsewhere [23]. In short, wavenumber axis was scaled according to where ν and ν' are original and corrected wavenumbers, respectively. Parameters were α F = 0.97, ν F = 4000, and ν 0 = 600. No imaginary wavenumbers were obtained in the calculated spectra.

Assignment of Spectral Bands of IMHA and Molecular Fragment
The IMHA comprises four structure elements: (1) the thiol group (-SH) that interacts with Au surface, (2) the alkane chain ((-CH 2 -) 5 ) spacer, (3) the amide group that increases SAM stability [24], and (4) the functional Im ring. In order to have control over the imidazole ring concentration at the interface, a surface backfiller (molecular fragment, Frag) was introduced into the monolayers. Molecular fragment has a similar chemical structure to IMHA, except for the Im ring, which is substituted with the methyl group ( Figure 1). Figure 3 shows ATR-FTIR spectra of IMHA, Frag, and an IMHA molecule recrystallized from ethanol-d 6 (IMHA-d). Assignments of the Im ring and amide spectral modes, based on the literature and our density functional theory (DFT) calculations, are given in Table 1 [4,[6][7][8]11,[25][26][27][28][29].
The two most pronounced bands in the ATR-FTIR spectrum of solid IMHA are assigned to the Amide-I (Am-I) vibration at 1635 cm −1 , which is mainly C=O coupled with the out-of-phase C-N stretching and C-C-N bending motion, and Amide-II (Am-II) at 1573 cm −1 , which is associated with the NH bending coupled with the C-N stretching [30,31]. Due to recrystallization in ethanol-d 6 , the mobile N-hydrogens exchange with D from the solvent (H/D exchange), which decreases the vibrational frequencies of the Am-I and Am-II modes to 1622 and 1467 cm −1 , respectively. Deuterium substitution also breaks the vibrational coupling between ND bending and C-N stretching and effectively changes the direction of the transition dipole moment (TDM) vector of the Am-II vibration [32]. The aliphatic chains of IMHA and Frag give rise to modes near 2852 and 2926 cm −1 ascribed to symmetric and asymmetric methylene-stretching motion and also scissoring deformation of CH 2 at 1451-1463 cm −1 . We identify symmetric deformation of the terminal methyl group, δ s (CH 3 ), of the molecular fragment for which the frequency at 1407 cm −1 is affected by the induction effect from the neighboring nitrogen of the amide group [33]. The complex N-H stretching pattern above 3000 cm −1 upon H/D exchange shifts to 2240-2350 cm −1 region, where three well-defined bands emerge. Medium intensity features at 2352 and 2288 cm −1 are linked with Amide-A' and Amide-B' [34,35]. The 2242 cm −1 mode in the IMHA spectrum is assigned to the vibration of the imidazole ring. The H/D exchange allowed for the identification of weak imidazole ring features otherwise obscured by strong infrared absorption bands. Modes at 3138 and 3117 cm −1 in the IMHA-d spectrum are associated with ν(C5-H) and ν(C2-H) vibrations, respectively-2242 cm −1 with N1-D stretching and 1560 cm −1 with ν(C4=C5). In the IMHA spectrum, the mode at 1494 cm −1 is easily recognized as C2-N3 stretch coupled with in-plane C2H bending of the Im ring, which in the IMHA-d spectrum appears as a shoulder at 1480 cm −1 (1484 cm −1 in DFT for IMHA) [27,29]. Figure 4A shows the time-dependent adsorption SEIRAS spectra of the IMHA monolayer formation in the ethanol incubation solution. The planar amide bond, -CO-NH-, gives rise to Am-I and Am-II vibrations, for which the transition dipole moments (TDM) are orientated perpendicularly to each other ( Figure S1) [25]. Metal surface selection rule states that the intensities of spectral modes scale with respect to the TDM vector projection to the surface normal [36]; therefore, any changes in the amide plane orientation will be reflected in the altered Am-I and Am-II intensity patterns, which, conveniently, could be expressed as the intensity ratio Am-I/Am-II. After 15 s of incubation, the SEIRAS spectrum already exhibited clear features of Am-I (1626 cm −1 and shoulder at 1649 cm −1 ) and C-N stretching of the imidazole ring (1487 cm −1 ). At increasingly longer times, the Am-II mode at 1558 cm −1 emerged and subsequently became the predominant band. In previous studies, the Am-II was found to be the most intense infrared absorption mode for neat SAMs on atomically smooth Au surfaces formed from alkanethiol molecules with the intrachain amide group [37][38][39][40]. Strong Am-II is related to the amide C-N bond oriented perpendicular to the surface.  The wavenumbers of the amide spectral modes are telltale of the H-bond strength. The Am-I is predominantly sensitive to the H-bonding state of the amide's carbonyl group, whereas the Am-II mode is mostly sensitive to the hydrogen-bonding interactions at the N-H site [41]. The weakening of the H-bonding interactions at C=O···H results in upshifted Am-I frequencies, while the weakening at the N-H···O results in downshifted Am-II frequencies [25,40,42]. Based on the changes in the wavenumbers of amide modes, whereas the Am-II mode is mostly sensitive to the hydrogen-bonding interactions at the N-H site [41]. The weakening of the H-bonding interactions at C=O···H results in upshifted Am-I frequencies, while the weakening at the N-H···O results in downshifted Am-II frequencies [25,40,42]. Based on the changes in the wavenumbers of amide modes, we deduce that the formation of SAM adversely affects the intrachain packing of IMHA molecules at the amide group. The bulk IMHA molecules in the ATR-FTIR spectrum exhibited an Am-II absorption maximum at 1573 cm −1 , which, upon adsorption on the Au surface, shifted to 1552 cm −1 (at 15 s). In the next five minutes, frequency upshifted to 1559 cm −1 but then slowly downshifted by several cm −1 in a time course of 90 min ( Figure S2). Such a difference between the ATR-FTIR and SEIRAS modes (δ = 14-21 cm −1 ) shows a drastic reduction in H-bonding strength. Am-I, on the other hand, showed only a small-scale steady frequency decrease in the first five minutes, after which the changes remained minimal. Figure 4B clearly shows that, within the first 10 min of adsorption, the Am-I/Am-II ratio sharply decreased from 5.2 to 0.8 for the IMHA SAM, meaning that reorientation into more vertically aligned amide planes for molecules occurs concurrently with the adsorption on the Au surface. Such an alignment is necessary for the rigid lateral H-bond network formation of molecules in the monolayer.

SEIRAS Analysis of the Formation of Self-Assembled Monolayer
To allow control over the IMHA molecule concentration on the surface and to enhance the motional flexibility of the Im group, we have synthesized and introduced molecular fragments into the monolayer composition. The increasingly higher molecular-fragment concentration clearly affected the amide group reorientation dynamics, such that it became increasingly faster and the final Am-I/Am-II ratios for mature SAMs were lower. For example, the ratio asymptotically approached progressively lower values, specifically 0.4, 0.2, 0.2, and 0.1 for the SAMs formed from 100, 50, 20, and 0 mol% of IMHA. The reduction of the surface IMHA concentration accelerated the narrowing of the Am-II band, expressed as the full width at half maximum (FWHM) as well ( Figure 4C), which suggests that the increased surface Frag concentration and incubation time synergistically acted on the formation of more homogeneous H-bond strength distribution in the SAM. Overall, the presented data suggest that the relatively bulky imidazole ring group of the IMHA molecule affects the packing ability of the molecules at the amide group and the formation of an extended H-bond network.
In our recent study based on the temperature-Raman and theoretical modeling, the Im ring ν(C2-N3) + β(C2H) mode was found to depend on the H-bond strength at the N1/N3 atoms [21]. DFT predicted that interaction of the Im ring with explicit water molecules resulted in a vibrational frequency shift by 2-6 cm −1 to higher frequencies, while an addition of two water molecules shifted the vibrational mode by 9 cm −1 . These data linked higher ν(C2-N3) vibrational frequency with the stronger H-bonding interactions of the Im ring. Here, we observe a wavenumber decrease from 1494 cm −1 for bulk IMHA (ATR-FTIR) to 1487 cm −1 for IMHA in monolayer (SEIRAS) at initial adsorption times. Such a decrease might be related to reduced H-bonding strength in ethanol-diluted IMHA molecules (0.5 mM) in the incubation solution. As the molecules chaotically adsorbed on a surface, they retained the melted-state character, which, upon reorientation into more concerted and solid-like structures with a gradual SAM development, was lost and resulted in upshifted wavenumbers ( Figure 5). The increasing wavenumbers suggest the incrementally stronger H-bond interactions at the Im ring occur at longer incubation times. Notably, the dynamic curve in Figure 5 reaches a plateau in approximately 60 min, which is in contrast to ca. 30 min found for the development of Am-I/Am-II intensities of the IMHA SAM. Thus, the Im ring properties at the interface directly depend on the structure and packing of the interchain amides.
in upshifted wavenumbers ( Figure 5). The increasing wavenumbers suggest the incrementally stronger H-bond interactions at the Im ring occur at longer incubation times. Notably, the dynamic curve in Figure 5 reaches a plateau in approximately 60 min, which is in contrast to ca. 30 min found for the development of Am-I/Am-II intensities of the IMHA SAM. Thus, the Im ring properties at the interface directly depend on the structure and packing of the interchain amides.

SEIRAS Analysis of the Structure of IMHA SAMs
To study in detail the way the self-assembled monolayer composition affects its structure, we compared the midrange spectra of various compositions of mature SAMs in ethanol and also the high-frequency spectra of the same SAMs in D2O ( Figure 6). We performed an ethanol/D2O exchange to extract spectral information from the aliphatic and aromatic C-H stretching region, otherwise obscured by strong ethanol-induced aberrations and complex N-H-and O-H-stretching patterns. Analysis of the Am-I/Am-II ratio shows the amide group tendency to become increasingly more perpendicularly oriented with the surface for SAMs, which are composed of higher molecular fragment concentrations (lower concentrations of IMHA), presumably because of reduced sterical interference from the rather bulky neighboring Im ring. Such an assumption is confirmed by an 8.5 cm −1 frequency upshift of Am-II mode, which is linked with stronger H-bonding interaction through the N-H group. The decreasing bandwidths from 44 to 41 cm −1 for Am-II

SEIRAS Analysis of the Structure of IMHA SAMs
To study in detail the way the self-assembled monolayer composition affects its structure, we compared the midrange spectra of various compositions of mature SAMs in ethanol and also the high-frequency spectra of the same SAMs in D 2 O ( Figure 6). We performed an ethanol/D 2 O exchange to extract spectral information from the aliphatic and aromatic C-H stretching region, otherwise obscured by strong ethanol-induced aberrations and complex N-H-and O-H-stretching patterns. Analysis of the Am-I/Am-II ratio shows the amide group tendency to become increasingly more perpendicularly oriented with the surface for SAMs, which are composed of higher molecular fragment concentrations (lower concentrations of IMHA), presumably because of reduced sterical interference from the rather bulky neighboring Im ring. Such an assumption is confirmed by an 8.5 cm −1 frequency upshift of Am-II mode, which is linked with stronger H-bonding interaction through the N-H group. The decreasing bandwidths from 44 to 41 cm −1 for Am-II and from 48 to 35 cm −1 for Am-I show that the molecular fragment favors the more uniformly arranged H-bonding network.
The weak mode at 1493 cm −1 ascribed to ν(C2-N3) + β(C2H) motion is detected only for the monolayers of the highest IMHA concentrations. Two other Im-related features at 3115 and 3152 cm −1 are assigned to ν(C2-H) and ν(C5-H) in Figure 6B. These bands are absent from the molecular-fragment spectrum, confirming the assignment to the =C-H stretching motion of the Im ring. Unexpectedly, the relative intensity of this doublet in the D 2 O environment remains largely unaffected by the IMHA's concentration in the SAM. This could only be possible if the Im ring group adopts an increasingly more perpendicular surface orientation as its concentration diminishes. Notably, a more perpendicular alignment is observed in D 2 O rather than ethanol-d 6 solutions ( Figure S3), as evidenced by the intensity of the high-frequency C-H stretching doublet.
In the high-frequency region, 2943-2947 cm −1 mode is assigned to the ν s (CH 3 ) motion of the molecular fragment. The relatively high mode's position is related to the induction effect produced by the neighboring nitrogen atom, which effectively increases the stretching frequency of the methyl group [33]. The vibrational modes at 2856 and 2927 cm −1 are assigned to ν s (CH 2 ) and ν as (CH 2 ), respectively. Their intensities depend on the carbon chain alignment with respect to the surface. There is an acute decrease in the intensity of the 2856 cm −1 mode when bicomponent 1/1 SAMs is compared to 1/5 ( Figure 6B). This could be explained by alkane chains adopting a vertical surface orientation so that the TDMs of their methylene groups are parallel with the surface and no longer SEIRAS-active. We also observe an increased intensity of the molecular-fragment-related symmetric CH 3 stretching mode. In fact, the absence of its asymmetric counterpart, ν as (CH 3 ), suggests that methyl groups are being oriented orthogonally with the surface. It is known that the ν s (CH 2 ) and ν as (CH 2 ) are highly sensitive to the ordering of the alkane chains so that densely packed and well-ordered SAM absorbs infrared at 2850 and 2917 cm −1 , whereas those of less-ordered arrangements have absorptions at higher frequencies [43][44][45][46]. Methylene bands are found to be upshifted by ca. 10 cm −1 in SEIRAS spectra, suggesting a melted structure of the alkane chains in the SAM, despite the fact that, from the intensity pattern, we found an almost vertical chain arrangement. Obviously, the alkane chain packing is strongly affected by the relatively bulky amide and imidazole ring groups. For comparison, molecules in the bulk phase exhibited wavenumbers nearly the same as the SAM, e.g., 2852 and 2926-2928 cm −1 .  The weak mode at 1493 cm −1 ascribed to ν(C2-N3) + β(C2H) motion is detected only for the monolayers of the highest IMHA concentrations. Two other Im-related features at 3115 and 3152 cm −1 are assigned to ν(C2-H) and ν(C5-H) in Figure 6B. These bands are absent from the molecular-fragment spectrum, confirming the assignment to the =C-H stretching motion of the Im ring. Unexpectedly, the relative intensity of this doublet in the D2O environment remains largely unaffected by the IMHA's concentration in the SAM. This could only be possible if the Im ring group adopts an increasingly more perpendicular surface orientation as its concentration diminishes. Notably, a more perpendicular alignment is observed in D2O rather than ethanol-d6 solutions ( Figure S3), as evidenced by the intensity of the high-frequency C-H stretching doublet.
In the high-frequency region, 2943-2947 cm −1 mode is assigned to the νs(CH3) motion of the molecular fragment. The relatively high mode's position is related to the induction effect produced by the neighboring nitrogen atom, which effectively increases the stretching frequency of the methyl group [33]. The vibrational modes at 2856 and 2927 cm −1 are assigned to νs(CH2) and νas(CH2), respectively. Their intensities depend on the carbon chain alignment with respect to the surface. There is an acute decrease in the intensity of the 2856 cm −1 mode when bicomponent 1/1 SAMs is compared to 1/5 ( Figure 6B). This could be explained by alkane chains adopting a vertical surface orientation so that the TDMs of their methylene groups are parallel with the surface and no longer SEIRAS-active. We also observe an increased intensity of the molecular-fragment-related symmetric CH3 stretching mode. In fact, the absence of its asymmetric counterpart, νas(CH3), suggests that methyl groups are being oriented orthogonally with the surface. It is known that the νs(CH2) and νas(CH2) are highly sensitive to the ordering of the alkane chains so that densely packed and well-ordered SAM absorbs infrared at 2850 and 2917 cm −1 , whereas We deduce from the SEIRAS spectra of varied-composition SAMs that the decrease in the relative IMHA concentration (i) affected the imidazole ring, amide group, and the alkane chain orientation, prompting them to adopt a more perpendicular orientation to the surface and (ii) strengthened the H-bonding interaction at the N-H group of the amide moiety. These observations imply stiffer head and tail group segments of the mixed SAMs, including an increased flexibility and accessibility to the interactions with the functional imidazole ring group.

Hydrogen-Bonding Interactions Revealed by Potential-Dependent SEIRAS
In order to understand the influence that electric potential has on the structure of mature SAMs and to investigate the H-bonding interactions at Im, series of potential-dependent SEIRAS measurements were performed for various composition SAMs. Figure 7A,B show the spectra of IMHA SAM at electrode polarizations with 0.3, −0.1, and −0.5 V and at 0.3 V after the excursion to negative potentials. Interestingly, both Am-II and Am-I shift to lower wavenumbers by 5-6 cm −1 as potential is stepped from 0.3 to −0.5 V, which is related to a decreasing H-bonding strength at N-H site and an increase at the C=O group. We can only speculate possible reasons for such a wavenumber shift discrepancy. One explanation could be that the decreased H-bonding interaction strength between the interchain N-H···O=C groups is related to an increased interaction between the carbonyl and the water molecules (C=O···HOH) at a negatively electrified interface. The Am-II wavenumber dependency on potential is similar for all tested SAM compositions and is related to the gradually weaker H-bonding strength towards the negatively polarized electrode ( Figure 7C). The weakening is related to reorientation of the amide planes into a less perpendicular arrangement with the surface (Figure 7D), supported by a sharp Am-I/Am-II ratio increase at -0.5 V for all but the IMHA monolayer. For the IMHA monolayer ( Figure 7D red line), strong intermolecular interactions at the Im ring groups bring rigidness to the molecular structure that suppresses amide group flexibility. Other SAM compositions showed a much stronger potential dependency of Am-I/Am-II intensity ratio and an altered Am-I band shape ( Figure S4). An unexpected increase in the low-frequency component at 1628-1630 cm −1 of the Am-I band was found for each SAM-containing molecular fragment at the −0.5 V bias. Such low wavenumbers in proteins were linked with a rigid and highly ordered β-sheet secondary structure element [30]; however, here we ascribe changes with the reorientation of C=O dipoles more closely along the external field lines (thus the increased Am-I intensity) and the establishment of H-bonding interactions with H 2 O molecules in a liquid solution. The SEIRAS information on the Im ring could be accessed from the 1492 cm −1 mode of the monolayer submerged in H2O and from a duplet at 3117 and 3151 cm −1 in D2O. While the former mode is poorly resolved, the duplet intensity suffices to provide a deeper understanding of the biased electric-potential-induced changes in the imidazole ring group. In the case of the T-I form, our DFT calculations and the literature ascribe the lower frequency component at 3117 cm −1 to =C2-H stretching and the 3151 cm −1 mode to the =C5-H [6,47]. For the molecule in the T-II state, these vibrations become coupled and engage in out-of-phase and in-phase vibrational motions, respectively. Surface adsorption induces the transition from the T-I to the T-II form; therefore, we consider the T-II (N1, N3-D) as the main imidazole tautomeric form in our system [4,21]. The frequency tuning rate of −10.4 and −6.3 cm −1 V −1 for the low and high-frequency modes was identified when electrode polarization was changed from 0.3 to −0.5 V ( Figure 8B  The SEIRAS information on the Im ring could be accessed from the 1492 cm −1 mode of the monolayer submerged in H 2 O and from a duplet at 3117 and 3151 cm −1 in D 2 O. While the former mode is poorly resolved, the duplet intensity suffices to provide a deeper understanding of the biased electric-potential-induced changes in the imidazole ring group. In the case of the T-I form, our DFT calculations and the literature ascribe the lower frequency component at 3117 cm −1 to =C2-H stretching and the 3151 cm −1 mode to the =C5-H [6,47]. For the molecule in the T-II state, these vibrations become coupled and engage in outof-phase and in-phase vibrational motions, respectively. Surface adsorption induces the transition from the T-I to the T-II form; therefore, we consider the T-II (N1, N3-D) as the main imidazole tautomeric form in our system [4,21]. The frequency tuning rate of −10.4 and −6.3 cm −1 V −1 for the low and high-frequency modes was identified when electrode polarization was changed from 0.3 to −0.5 V ( Figure 8B). Bicomponent monolayers exhibited similar tuning rates of −8.  We conducted a DFT modeling of 4-ethyl-1-imidazole in Tautomer-II form (N1, N3-D) with zero, one, and two explicit D2O molecules in order to explain the effects that Hbonding interaction has on the vibrational frequency of ν(=C-H) modes ( Figure 9). We found that the N1 atom in T-II is more sensitive to hydrogen bonding than the N3-D deuteron. For example, the interaction at N3-D group (N3-D···D2O) has little effect on the vibrational frequency of the modes (δ = −3 cm −1 ), whereas the interaction at the N1 atom increases the frequencies by 8 cm −1 . An almost exact tendency was found for a compound in T-I form. A considerably shorter H-bond at N1 compared to N3-D by 8.5 pm shows the greater stabilization over the ImN1···D2O complex. Indeed, a supersonic jet FTIR spectroscopy study suggested that Im prefers accepting the H-bond at N···H−O rather than donating it at N−H···O [48]. However, even stronger Im stabilization was attained with an addition of a second D2O molecule to the Im-D2O complex. The coordination with two D2O molecules resulted in the shortening of both H-bonds by 1.89 and 2.69 pm for the N1 and N3-D sites. The infrared absorption spectroscopic measurements of the model compound imidazole-4-methanol (Im-CH2-OH) in solvents of different polarity ( Figure S5) clearly shows decreased frequencies of the ν(=C-H) modes in methanol-d4, where a less-strong hydrogen-bonding interaction is expected compared to that in D2O. We conducted a DFT modeling of 4-ethyl-1-imidazole in Tautomer-II form (N1, N3-D) with zero, one, and two explicit D 2 O molecules in order to explain the effects that Hbonding interaction has on the vibrational frequency of ν(=C-H) modes ( Figure 9). We found that the N1 atom in T-II is more sensitive to hydrogen bonding than the N3-D deuteron. For example, the interaction at N3-D group (N3-D···D 2 O) has little effect on the vibrational frequency of the modes (δ = −3 cm −1 ), whereas the interaction at the N1 atom increases the frequencies by 8 cm −1 . An almost exact tendency was found for a compound in T-I form. A considerably shorter H-bond at N1 compared to N3-D by 8.5 pm shows the greater stabilization over the ImN1···D 2 O complex. Indeed, a supersonic jet FTIR spectroscopy study suggested that Im prefers accepting the H-bond at N···H−O rather than donating it at N−H···O [48]. However, even stronger Im stabilization was attained with an addition of a second D 2 O molecule to the Im-D 2 O complex. The coordination with two D 2 O molecules resulted in the shortening of both H-bonds by 1.89 and 2.69 pm for the N1 and N3-D sites. The infrared absorption spectroscopic measurements of the model compound imidazole-4-methanol (Im-CH 2 -OH) in solvents of different polarity ( Figure S5) clearly shows decreased frequencies of the ν(=C-H) modes in methanol-d 4 , where a less-strong hydrogen-bonding interaction is expected compared to that in D 2 O.
Based on our theoretical and experimental analysis, we linked the =C-H frequency shift towards lower wavenumbers with decreased imidazole H-bond strength. The potential dependency in Figure 8 clearly shows decreasing wavenumbers with negative electrode polarization indicating the reversible reduction of imidazole's H-bond strength with its immediate surroundings. One can associate such a wavenumber shift with the vibrational Stark effect; however, no significant difference in the frequency tuning rates was found for Im-SAM in the environments of various ionic strengths, as shown in Figure S6. The weakening of the H-bond may be explained by the structural changes of interfacial water. The dielectric constant of water is reduced near the surfaces or under confinement, largely because of the reorientation of water molecule dipoles [49]. The bias electric potential polarizes medium by acting on molecular dipoles and competes against the isotropic H-bond network. Substantially negative polarization orients O-D dipoles to point away from the bulk to the surface [50]. Overall, the arrangement of water molecules at the electrified interface affects the hydrogen interaction strength with exposed Im groups. Based on our theoretical and experimental analysis, we linked the =C-H frequency shift towards lower wavenumbers with decreased imidazole H-bond strength. The potential dependency in Figure 8 clearly shows decreasing wavenumbers with negative electrode polarization indicating the reversible reduction of imidazole's H-bond strength with its immediate surroundings. One can associate such a wavenumber shift with the vibrational Stark effect; however, no significant difference in the frequency tuning rates was found for Im-SAM in the environments of various ionic strengths, as shown in Figure S6. The weakening of the H-bond may be explained by the structural changes of interfacial water. The dielectric constant of water is reduced near the surfaces or under confinement, largely because of the reorientation of water molecule dipoles [49]. The bias electric potential polarizes medium by acting on molecular dipoles and competes against the isotropic H-bond network. Substantially negative polarization orients O-D dipoles to point away from the bulk to the surface [50]. Overall, the arrangement of water molecules at the electrified interface affects the hydrogen interaction strength with exposed Im groups.
The relative intensities of the 3150 and 3115 cm −1 modes require some explanation as well. The calculations predict that both displacement vectors of the T-II form of the Im ring lay within the Im plane. The displacement vector of the 3115 cm −1 vibrational mode is parallel and points along the C2-H bond direction, whereas the displacement vector of the 3150 cm −1 vibrational mode is almost parallel to the N3-D bond direction ( Figure S8). Thus, changes in the Im ring surface orientation are reflected in the relative intensities of the 3150 and 3115 cm −1 modes. The transition to the negative potentials decreases the I3150/I3115 intensity ratio, showing that the imidazole tilts towards the position where C2-H is slightly more exposed to the solution. In contrast to the two component SAMs, we found that the potential-induced intensity alterations were slightly dampened for the monolayer composed of 100 mol% IMHA, most likely due to the stronger intermolecular interactions between the Im ring groups. The relative intensities of the 3150 and 3115 cm −1 modes require some explanation as well. The calculations predict that both displacement vectors of the T-II form of the Im ring lay within the Im plane. The displacement vector of the 3115 cm −1 vibrational mode is parallel and points along the C2-H bond direction, whereas the displacement vector of the 3150 cm −1 vibrational mode is almost parallel to the N3-D bond direction ( Figure S8). Thus, changes in the Im ring surface orientation are reflected in the relative intensities of the 3150 and 3115 cm −1 modes. The transition to the negative potentials decreases the I 3150 /I 3115 intensity ratio, showing that the imidazole tilts towards the position where C2-H is slightly more exposed to the solution. In contrast to the two component SAMs, we found that the potential-induced intensity alterations were slightly dampened for the monolayer composed of 100 mol% IMHA, most likely due to the stronger intermolecular interactions between the Im ring groups.

Conclusions
In this work, we studied the imidazole (Im) ring, terminated at alkanethiol molecule self-assembled monolayers (SAMs), formed from IMHA and the surface-dilutor molecular fragment. SAM formation dynamics and interfacial Im ring structure were probed by H/D substitution SEIRAS. We found that despite the SAM composition, the alkane chain-nested amide groups adopt optimal surface orientation within the first 30 min of incubation. Based on the frequency shifting of the Im structure marker band near 1490 cm −1 , changes in the SAM terminal functional group did not cease until the 60 min of incubation. These changes were linked with the gradual increase of the hydrogen-bonding strength at the Im group. The IMHA SAM formed a structurally rigid organic layer with somewhat disturbed hydrogen-bonding network at amide groups. The introduction of a surfacedilutor molecular fragment affected the Im orientation in such a way that it became more vertically oriented, and the amide groups formed a stronger H-bonding network. Negative electric potentials reversibly decreased the vibrational frequency of Im's =C-H stretching modes at 3115 and 3150 cm −1 (frequency tuning rates were −10.4 and −6.3 cm −1 V −1 ). We assign a frequency decrease to the weakening of the H-bond interaction of the Im ring with the immediate water at the negatively electrified interface.