Investigation of rotameric conformations of substituted imidazo-[1,2-a]pyrazine: experimental and theoretical approaches

The different rotameric conformations of imidazo-[1,2-a]pyrazine have been synthesized and characterized by means of different experimental techniques, such as NMR, FTIR, and absorption spectroscopy and quantum chemical calculations. The different conformations were stabilized by hydrogen bonds, such as OH⋯N, ArH⋯N and ArH⋯ArH. The ground state optimizations and potential energy surface (PES) scanning profiles produced using density functional theory (DFT) show two stable rotameric forms for each molecule. The relative population of the conformations is affected by the strength of the hydrogen bonds. The calculated absorption spectra and isotopic shielding constants were acquired by time-dependent density functional theory (TD-DFT) and gauge invariant atomic orbitals (GIAO)-DFT, respectively. The strength of the hydrogen bonding interactions that resulted in the different conformations was studied by quantum theory of atoms in molecules (QTAIM).


Introduction
The heterocyclic compounds are essential candidates in the chemical, biological, agricultural, and veterinary elds. 1,2 Among the various heterocycles, imidazo [1,2-a]pyrazine is an aza-heterocyclic compound that has applications as an anticancer agent, 3 cardiac stimulating agent and uterine relaxant, 4 antihyperglycemic agent, 5 antihypertensive agent, 6 antiulcer agent, 7 antibronchospastic agent 8 etc., [9][10][11] The properties of compounds as an application candidate depend upon the molecular structure of the compound and can be tuned by altering the molecular modication, and orientation of the compound. 12 The spatial and chemical constraints in molecular systems give rise to different spatial arrangements or conformations. These conformations can be governed by non-covalent interactions, such as van der Waals interactions, hydrogen bonding, and steric hindrance signicantly. [12][13][14][15] These non-covalent interactions are the deciding forces in building supramolecular assemblies, and thus, this is a fascinating area in structural chemistry. These interactions connect together the synthons and building blocks of supramolecular structures. 12 Thus, the implication of such interactions in conformational analysis would result in new conformations of compounds in the solid, solution and gas phase. These interactions lock the molecular congurations through hydrogen bonding, and thus, permit the formation of different conformations. [16][17][18][19][20][21][22][23][24] Thus, the exact structure determination of a molecule is important and a big subject of research itself in science. Structural determination requires tough conformational analysis by means of experimental and theoretical calculations.
For decades, correct structural information has been very helpful in describing many chemical, physical and biological processes, such as reaction mechanisms and photophysical processes, and in designing molecular machines and understanding their functioning, etc. 25,26 The problem of exact structure determination is solved by experimental (e.g. spectroscopic techniques), single crystal X-ray diffraction (SCXRD) and theoretical (quantum chemical calculations) tools. Spectroscopic techniques, such as FTIR, NMR, absorption spectroscopy etc., provide very specic information about functional groups and their orientation, and thus, they are promising tools in solving molecular structures. Also, SCXRD provides the exact molecular conguration following good structural renement, but it requires high-quality crystals and growing good quality crystals is not always an easy task. Therefore, conformer assignment requires experimental data from different techniques, and data analysis exercises, and thus is always a challenging task through experimental techniques. Therefore, the use of computational tools is the best alternative to predict an accurate structure via different quantum chemical methods. Quantum chemical calculations can be used to predict the spectroscopic values for different possible structures. The correlation of experimental and theoretical results is also a good parameter to conrm the exact molecular structure. During a literature survey, we found some reports on the conformational examination of organic compounds by mean of experimental and theoretical techniques, but the conformation of imidazo-[1,2-a]pyrazine has not been determined. [27][28][29][30][31][32] In this work, we have synthesized a series of imidazo-[1,2-a]pyrazine derivatives that exhibit different rotameric conformations. These rotameric conformations were further puried by column chromatography. Moreover, the rotamers were characterized through NMR, UV and quantum chemical calculations. Furthermore, the experimental data have been correlated with the calculated values to gain insight into the rotameric structures.

Material and methods
All the chemicals used for synthesis were purchased from Sigma-Aldrich Chemical Ltd., Loba Chemie etc. depending upon their availability. All the solvents used were of spectroscopic grade and purchased from Spectrochem and Rankem Ltd. All chemicals and solvents were used without further purication. The progress of the chemical reaction was monitored by thin-layer chromatography (TLC). Melting points were recorded using the open capillary tube method and were uncorrected. 1 H NMR and 13 C NMR spectra were recorded on a JEOL ECS-400 MHz spectrometer at ambient temperature in CDCl 3 with TMS as an internal reference. All chemical shis were reported in ppm relative to the reference. Mass spectra of the synthesised compounds were recorded using a Waters Micromass Q-Tof micro mass spectrometer. The absorption spectra were recorded on a SHIMADZU-2600 spectrophotometer using quartz cuvettes of 1 cm path length.

Computational details
All theoretical calculations were performed using the GAUSSIAN-03W program. 35 The geometry optimizations for all molecules were carried out at the DFT/B3LYP level of theory using the 6-31++G(d,p) basis set. [36][37][38] Frequency calculations were performed on the obtained optimized structures at the same level of theory. The potential energy surfaces for the dihedral angles were acquired using DFT at the B3LYP/6-31+g* level of theory to identify the local minima. The theoretical absorption spectra were calculated from the vertical excitation energies using the time-dependent density functional theory (TD-DFT) method. The NMR shielding tensors were calculated using the GAIO method. 39 The conformational optimization was performed in different solvents (CH 3 CN and CHCl 3 ) for the optimized structures using a self-consistent reaction eld (SCRF) approach coupled with the integral equation formation of polarizable continuum model (IEFPCM) solvation method. Topological analysis of the electron distribution was conducted on selected conformations according to Bader's "Quantum Theory of Atoms in Molecule (QTAIM)". 40 QTAIM calculations were carried out using the Multiwfn program 41 to evaluate the nature and strength of the hydrogen bonds present in the different conformations.

Results and discussion
3.1. Synthesis of 2-amino-3,5-dibromopyrazine (4) 33 N-Bromosuccinamide (14.95 g, 83.99 mmol) was added over 50 min to a mixture of 2-aminopyrazine (3.80 g, 40 mmol) in 80 ml of DMSO and 2 ml of H 2 O below 15 C. The mixture was then stirred for 6 h at room temperature. Aer completion of the reaction, the mixture was extracted with water and ethyl acetate. The ethyl acetate layer was dried over Na 2 SO 4 and concentrated under vacuum. The crude product was puried by column chromatography using hexane : ethyl acetate (9 : 1) as an eluent. White solid; yield: 90%; mp 115-116 C; 1  To 2-amino-3,5-dibromopyrazine (5.0 g, 19.8 mmol) in 100 ml of isopropyl alcohol (IPA), a 50% aqueous solution of chloroacetaldehyde (99 mmol) was added dropwise. The reaction mixture was reuxed at 110 C for 24 h. Aer the completion of the reaction, the reaction mixture was cooled to room temperature and then extracted with water and chloroform. The chloroform layer was dried over sodium sulphate and concentrated under vacuum to get the crude product. The product was puried by column chromatography using hexane : ethyl acetate (6 : 4) as an eluent to obtain a pure compound as a white

Synthesis of compound 1
A vial equipped with a stirring bar was charged with 6,8dibromo-imidazo[1,2-a]pyrazine (5) (0.5 g, 1.8 mmol), Cs 2 CO 3 (0.6 g, 1.8 mmol) and boronic acid (1.8 mmol), dissolved in CH 3 CN : H 2 O (9 : 1) at 100 C under an inert atmosphere. Then, 5 mol% of Pd(PPh 3 ) 4 was added and the vial was capped. The reaction mixture was reuxed for 7-12 h. Aer the completion of the reaction (monitored by TLC), the reaction mixture was cooled, and then extracted with water and chloroform. The organic layer was dried over Na 2 SO 4 , ltered and concentrated under vacuum to obtain the crude product. The residue was puried by silica gel (60-120 mesh) column chromatography using hexane : ethyl acetate as an eluent. Two different coloured solid products, 1A and 1B, were obtained. 1A

Synthesis of compound 2
A vial equipped with a stirring bar was charged with 6-bromo-8substituted-imidazo[1,2-a]pyrazine (1A or 1B) (0.1 g, 0.323 mmol), Cs 2 CO 3 (0.105 g, 0.323 mmol) and 2-hydroxyphenylboronic acid (0.323 mmol), dissolved in CH 3 CN : H 2 O (9 : 1) at 100 C under an inert atmosphere. Then, 5 mol% of Pd(PPh 3 ) 4 was added, the vial was capped and the mixture was reuxed for 6-12 h. Aer the completion of the reaction (monitored by TLC), the reaction mixture was cooled, and then extracted with water and chloroform. The organic layer was dried over Na 2 SO 4 , ltered and concentrated under vacuum to obtain the crude product. The residue was puried by silica gel (100-200 mesh) column chromatography using hexane : ethyl acetate as an eluent. Interestingly, the molecular mass from ESI-MS was found to be the same for compounds 1A and 1B (Fig. S3, and S7 †), and likewise for 2A and 2B (Fig. S14, and S19 †), but with different NMR signals (Fig. S1, S8, S5, S6, S12, S13, S17 and S18 †). Furthermore, all the synthesized compounds have different absorption maxima ( Fig. S10 and S21 †), and characteristic FTIR signals (Fig. S4, S8, S15 and S20 †), which clearly indicates that the obtained compounds 1A and 1B are isomeric forms of each other, as are 2A and 2B. The obtained compounds have a rigid architecture, except for rotation around the single bond connecting the phenyl ring to the imidazo-[1,2-a]pyrazine moiety. Therefore, the only possibility is for the isomeric forms to be different rotameric conformations. Thus, in order to determine the accurate structures of the possible conformers of compounds 1 and 2, theoretical calculations were performed.

Geometry optimization and reactivity of compound 1
Compound 1 was optimized, and was further subjected to Fukui indices analysis, which is widely used to predict the reactivity indices in order to observe the reactive sites (electrophilic and nucleophilic sites) present on a molecule. 42 Fukui functions are the partial derivatives of electron density with respect to the number of electrons present. where r(r) denotes the electron density, and N is the number of electrons present in the molecular system. Generally, the indices are calculated for the HOMO and LUMO, and represent the nucleophilicity and electrophilicity index, respectively. Reactive sites in the molecular system have larger Fukui index values ( Fig. 1). This journal is © The Royal Society of Chemistry 2018 3.5.1. Conformational analysis of compound 1. The aim of nding the most stable conformation has been fullled using a simple general strategy.
1. Compound 1 has rotational exibility with respect to the single bond joining the hydroxy phenyl ring and the rigid imidazo-[1,2-a]pyrazine ring. Therefore, a conformational search was carried out using the rotational potential energy surface (PES) produced by DFT at the B3LYP/6-31+G* level. All possible conformations of compound 1 at local minima are shown in Scheme 2.
2. The obtained structures at local minima were analyzed to determine their relative Boltzmann contributions. From the Boltzmann distribution principle, any conformers with a free energy in the region of 5 kcal mol À1 can make some contribution to the Boltzmann population at room temperature (RT). Thus, the contributing conformers have been further optimized using DFT at the B3LYP/6-31++G(d,p) level and conrmed as having no imaginary frequencies.
3. The absorption spectra and NMR chemical shis were computed for each of the optimized conformers using the TD-DFT and GIAO methods.
4. The experimental and calculated data for the absorption spectra, FTIR signals, and NMR shielding tensors were compared and evaluated to determine the structure of the compounds.
3.5.2. Structural analysis of compound 1. At rst, compound 1 was optimized in the form represented by 1A. Then, the potential energy surface (PES) of the dihedral angle for rotation of the single bond connecting the phenyl ring and imidazo-[1,2-a]pyrazine moiety was acquired for 1A. The PES resulted in a new local minimum, with the conformation represented by 1B. Both conformers were found to be stabilized by hydrogen bonding, primarily examined by looking at the bond distance of OH and NH, and it was found that the conformer at the local minimum 1B has an energy close to that of 1A. On the other hand, compound 1 was also optimized as conformer 1C to determine the impact of hydrogen bonding in the stabilization of conformers 1A and 1B. Also, the dihedral angle PES was analyzed to check for other local minima. The new local minimum observed was assigned as conformer 1D. The relative energy diagram for all conformations is shown in Fig. 2. It is found that both 1C and 1D have high energy relative to 1A and 1B. The relative free energies of the optimized conformations, along with their Boltzmann distributions, are summarized in Table 1. From the Boltzmann distribution principle, conformers with a free energy in the region of 5 kcal mol À1 can make some contribution to the Boltzmann population at room temperature (RT). Furthermore, conformers 1A and 1B have been optimized at the B3LYP/6-31++G(d,p) level. Both conformers, 1A and 1B, were found to exhibit a planar geometry and be stabilized by six-membered and seven-membered hydrogen bonding, respectively. Conformer 1A is the major product with an abundance of about $92.5%, while 1B, which is close in energy to conformation 1A, makes a relatively low ($7.5%) contribution. No other conformations make a signicant contribution. The optimized structures of the conformations are shown in Fig. 3, along with important structural parameters, such as bond lengths and bond angles. The stabilization of the conformations is discussed below: 3.5.3. Conformer 1A. According to the Boltzmann distribution and optimized geometry, conformer 1A is the key conformation in the gas phase. The stability of this conformation can be explained by the relative charges and hydrogen bond distances.
The calculated non-bonded distance between N7 and H24 bonded to O23 [d (N7/H24) ¼ 1.678 A] is less than the sum of the van der Waals radii of hydrogen and nitrogen. The smaller Scheme 2 Possible conformations for compound 1.   ))] atoms establish an intramolecular dipolar interaction between a phenolic hydrogen (O23-H24) and the pyrrole nitrogen (N4).
The non-bonded distance between the phenolic hydrogen and pyrrole nitrogen [d (H24/N4) ¼ 1.592 A], which is less than the sum of the van der Waals radii of hydrogen and nitrogen, supports hydrogen bonding. The smaller distance and columbic interaction with opposite charges on H24 and N4 conrmed the seven-membered intramolecular hydrogen bonding.
In addition to the N/H-O hydrogen bonding interactions in conformers 1A and 1B, there are ArH/N interactions, where the aromatic hydrogen of the phenolic moiety interacts with the pyrrole nitrogen for conformer 1A and the imidazole nitrogen for conformer 1B. These interactions are weak in nature, as conrmed by the interaction distances and topological parameters. Thus, conformer 1B exhibits seven-membered intramolecular hydrogen bonding, which is less stable than the sixmembered hydrogen bonding exhibited by conformer 1A, and this could be the reason for the low Boltzmann contribution for conformer 1B.
Additionally, both conformers 1A and 1B were free from interelectronic repulsion, which is present in conformers 1C and 1D. Furthermore, the hydrogen bonding interaction was also estimated by QTAIM calculations using the Multiwfn program (Table 2).
3.5.5. Topology parameters. In order to determine the nature of the intramolecular hydrogen bonds (HBs) existing within the key conformers, they were studied by means of Bader's quantum theory of atoms in molecules (QTAIM) using  Table 2). The positive value of the Laplacian of the electron density [V 2 r(r c )] for both conformers 1A and 1B, indicated electrostatic closed shell interactions and the value was found to be within Popelier's hydrogen bonding criteria. The Laplacian of the electron density for N/HO hydrogen bonding in conformer 1A is higher than that in 1B (Table 2). On the other hand, ArH/N has more Laplacian electron density in conformer 1B. Furthermore, the degree of covalency and strength of the interactions were characterized according to Rozas' rules: (i) V 2 r(r c ) < 0 and H(r) < 0 for strong H-bonds of covalent nature, (ii) V 2 r(r c ) > 0 and H(r) < 0 for medium H-bonds of partially covalent nature, and (iii) V 2 r(r c ) > 0 and H(r) > 0 for weak H-bonds of electrostatic character. 30 All of the bond critical points in Table 2 (V 2 r(r c ) > 0 and H(r) < 0) are designated as medium range interactions with partial covalent character for both conformers. Moreover, the hydrogen bond energy is determined using the Espinosa equation, 44 which states that the interaction energy of a H/X contact is dened as E HX ¼ V(r)/2 at the BCP, where V(r) is the potential electron density at the bond critical point. Conformer 1A exhibited two interactions, N/HO and ArH/N, which contributed À14.18 and À4.11 kcal mol À1 of stabilization. On the other hand, conformer 1B also exhibited N/HO and ArH/N interactions, which contributed À18. 22 and À4.39 kcal mol À1 of stabilization, respectively. The ArH/N interactions make nearly the same contribution for both conformers 1A and 1B, while the N/HO hydrogen bonds were found to be weaker in conformer 1A than 1B. This is due to the presence of an electron withdrawing bromine atom, which attracts the electron density, and therefore weakens the hydrogen bonding. On the other hand, the lone pair on the nitrogen atom could easily interact with a polarized hydrogen atom in conformer 1B, forming strong hydrogen bonding. The bond critical points and bond paths of conformers 1A and 1B are shown in Fig. 4. 3.5.6. MBO and LBO analysis for compound 1. In addition, the Laplacian bond order (LBO) and Mayer bond order (MBO) were analyzed for stable conformers 1A and 1B. LBO has a direct correlation with the bond polarity and vibrational frequency. Bonds with a bond order value smaller than 1.0 are polarized in nature, and the lower the value, the greater the polarization. The Mayer bond order and Laplacian bond order are summarized in Table 3. The MBO for C15-O23 in conformer 1A is slightly lower than that for conformer 1B. On the other hand, a smaller value of MBO and LBO for O23-H24 reects more bond polarity for conformer 1B than conformer 1A. The polarization of the OH bond causes a difference in bond strength, which results in a shi in the vibrational frequency of the two different conformations and validates the intramolecular hydrogen bonding.
3.5.7. Frequency analysis. FTIR is a good tool for functional group analysis. The FTIR signals of functional groups including OH, C]N, C]C etc. were observed for both conformers 1A and 1B. All computed band positions are in good agreement with the experimental values for both conformers. However, the OH stretching vibrations, which are usually environmentally dependent, suggestively distinguish both conformations. The The computed frequency results showed the O-H stretching band at 3230.25 cm À1 in FTIR for conformer 1A. Experimentally a broad O-H stretching band is observed at 3158 cm À1 for conformer 1A (Fig. S4 †). The computed and experimental FTIR values for conformer 1B were at 2999.47 cm À1 and 3069 cm À1 , respectively (Fig. S8 †). The broadness in the band also conrmed the presence of intramolecular hydrogen bonding. The difference in the FTIR values of the O-H bands for 1A and 1B revealed that the OH bond was slightly weak in conformer 1A relative to 1B. As discussed earlier, the more polarized nature of the O-H bond as indicated by LBO bond order analysis and the strong hydrogen bonding indicated by QTAIM analysis for conformer 1B further weaken the O-H bond strength for 1B. Therefore, the FTIR values were in good agreement with the observations from topology analysis, hydrogen bond strength, and LBOs.
3.5.8. NMR analysis. The 1 H and 13 C NMR shielding tensors were calculated using the GIAO method for the optimized structures, and the effect of the solvent was included using the IEFPCM model provided in Gaussian 03 W. The experimental chemical shis of 1 H NMR and 13 C NMR were located between 7.07 and 13.12 ppm, and 114.65 and 160.84 ppm, respectively. The comparison between the computed and experimental 1 H NMR and 13 C NMR results gave correlation values of R 2 ¼ 0.9888 and R 2 ¼ 0.9331, respectively, for conformer 1A (Fig. 5). The deviation in the experimental chemical shi of C11 from the computed one is only due to the attachment of the hydroxyl group to the phenyl ring. In a similar way, the experimental NMR of conformer 1B has been compared with the calculated data. Although the correlation was not good, the combined information from FTIR and NMR could distinguish the two conformers. The correlation of the experimental and calculated NMR spectra for conformer 1B is shown in Fig. S9. † Also, it has been observed that the hydrogen atoms involved in hydrogen bonding resonate at different strengths for both 1A and 1B. The aromatic proton at C6 0 Table 2 Topology parameters, including electron density (r), Laplacian electron density (V 2 r), potential energy density [V(r)], total energy density [H(r)], and hydrogen bonding energy [E HB (kcal mol À1 )] at the bond critical point of non-covalent interactions (D/HA) for conformers 1A and 1B at the B3LYP/6-31++G** level  involved in hydrogen bonding resonates at 9.54 ppm for compound 1A, while it resonates at 7.99 in the case of compound 1B (Table S1 †). 3.5.9. Absorption analysis. The vertical excitations were calculated for the optimized structures of conformers 1A and 1B in the gaseous and acetonitrile solution phases using TD-DFT. The obtained absorption value for conformation 1A at 380, 361, and 318 nm has a contribution from the rst three Franck-Condon excited states [(S 0 -S 1 ), (S 0 -S 2 ), and (S 0 -S 3 )], respectively. The rst transition results primarily from the HOMO to the LUMO, and has the nature of charge transfer from the HOMO to LUMO (from the phenolic moiety to the imidazole moiety), while the second and third transitions are kind of p-p* in nature and have major contributions from (HOMO-1) to LUMO and (HOMO-2) to LUMO transitions, respectively. On the other hand, the experimental absorption maximum value for conformer 1B was observed to be 290 nm, found due to the third Franck-Condon transition. The transition occurred from (HOMO-2) to LUMO due to charge transition in the opposite direction. The calculated and experimental absorption values for conformers 1A and 1B are summarized in Table S2 † and the contributing Frontier molecular orbitals are shown in Fig. 6 and S11, † respectively.
Thus, the synthesized conformers have been optimized in the form of 1A and 1B. Both conformers 1A and 1B have a planar geometry and were found to be stabilised by six-and sevenmembered hydrogen bonding, respectively. The stability and strength of the hydrogen bonds were established by QTAIM. The observed and calculated spectral data for NMR and FTIR were found to be in good correlation for each conformer.

Structural analysis of compound 2
As discussed above, the isomeric forms of compound 1 were found to have 1A and 1B geometry. Subsequently, the separated isomers were reacted with another boronic acid and this resulted in two new compounds. Therefore, the conformers that resulted from the second substitution could have the possible structures shown in Scheme 3. Conformer 1A, stabilized by six-membered hydrogen bonding, could further result in either form 2A or 2B. Therefore, compound 2 was rst optimized based on conformer 2A using the same protocol as for compound 1. The optimized geometry was conrmed by it having no imaginary frequencies. Furthermore, the optimized structure of conformer 2A was analysed through potential energy surface analysis for a dual dihedral scan of two phenolic moieties with the imidazo-[1,2-a]pyrazine moiety with variation from 0 to 180 at a step size of 20 degrees. The Boltzmann contributions of signicant structures with a relative energy of 5 kcal mol À1 on the potential energy surface were considered for further analysis. Four conformers represented by 2A, 2B, 2C and 2D were found to be within the Boltzmann energy limit for a signicant contribution (Fig. 7). The relative energy and Boltzmann distributions for all conformers at the minima have been calculated in Table 4. All the conformers were found to be stabilized by hydrogen bonding, established by the interaction distance primarily. Conformer 2B was found to be the most stable with the highest contribution of 91.9%, followed by conformer 2A with an $6.7% contribution. On the other hand, 2C and 2D have negligible Boltzmann distributions. The molecular structures with adequate Boltzmann distributions have been optimized and analyzed by frequency calculations. The optimized structures for the conformers are shown in Fig. 8 along with their characteristic structural parameters, such as bond lengths and bond angles. The stabilization of the conformers has been discussed on the basis of hydrogen bonding.
3.6.1. Conformer 2A. Conformer 2A has an $7% Boltzmann contribution, stabilized by two six-membered hydrogen bonds, which was further quantied by the relative charges, and hydrogen bonding distances.
The The calculated non-bonded distances between the phenolic hydrogen (H12, H13) and pyrazine/pyrrole nitrogen (N3/N2) [d (N2/H12) ¼ 1.663 A, and d (N3/H13) ¼ 1.764 A] were found to be less than the sum of the van der Waals radii of hydrogen and nitrogen. The smaller distance and columbic interaction with opposite charges on H12, H13 and N3 support six-membered intramolecular hydrogen bonding.
Additionally, both conformers 2A and 2B were free from interelectronic repulsion. Furthermore, the nature of the hydrogen bonding interaction was estimated by QTAIM calculations using the Multiwfn program.      (Fig. S15 †). On the other hand, conformer 2B has O-H bonds of almost the same strength and high stabilization; thus, the calculated FTIR signals were found at lower wavenumbers with a small difference of $70 cm À1 . Likewise, the observed FTIR signals for conformer 2B were at 3000, and 2930 cm À1 (Fig. S20 †). The observed and calculated FTIR signals for the O-H bonds in both conformers were in correlation with their strength and the polarization of the bonds.
3.6.6. NMR. The 1 H and 13 C NMR shielding tensors were calculated for conformers 2A and 2B using the GIAO method on the optimized structures, and the effect of the solvent is included using the IEFPCM model provided in Gaussian 03 W. The assessment of the computed and experimental 13 C NMR results revealed good correlation for both conformers 2A and 2B. The correlation value for 13 C NMR was found to be R 2 ¼ 0.97 and 0.99 for conformer 2A and 2B, respectively. Furthermore, the correlation value of the computed and experimental 1 H NMR results for conformer 2B was R 2 ¼ 0.97. However, the   Fig. S16. † 3.6.7. Absorption properties. The vertical Franck-Condon excitations were calculated for the optimized structures of conformers 2A and 2B, using TD-DFT, and further compared with the experimental absorption spectra. The obtained absorption values for conformer 2A at 380, 360 and 265 nm have a major contribution from three Franck-Condon excited states [(S 0 -S 1 ), (S 0 -S 5 ) and (S 0 -S 6 )], respectively. On the other hand, the experimental absorption peaks for conformer 2B were observed at 375, 315, and 260 nm. The experimentally observed transitions were found to be due to [(S 0 -S 1 ), (S 0 -S 2 ), and (S 0 -S 6 )] Franck-Condon transitions, respectively. The calculated and experimental absorption spectral data are tabulated in Table  S4. † The different contributing molecular orbitals are shown in Fig. S22. † Thus, the differences in the signals for the NMR, FTIR, and absorption peaks for the two compounds and their correlation with theoretical calculations clearly predict the specic orientation of the rotamers.

Conclusions
In summary, we have successfully synthesized a series of imidazo-[1,2-a]pyrazines, possessing different rotameric conformations, which were further puried by column chromatography. Theoretical investigation of the accurate structure of the different rotameric forms was done through PES and hydrogen bond strength analysis. The obtained rotameric forms were found to be stabilized by intramolecular hydrogen bond locking, which was further established by QTAIM analysis. The Laplacian and Mayer bond order (LBO/ MBO) revealed that differences in bond strength for the different conformations, and thus distinctions in the vibrational frequency and NMR signals, assisted in the recognition of the accurate structure for the different conformations.

Conflicts of interest
There are no conicts of interest to declare.