Density functional theory investigation of hydrogen bonding effects on the oxygen, nitrogen and hydrogen electric field gradient and chemical shielding tensors of anhydrous chitosan crystalline structure

A systematic computational investigation was carried out to characterize the 17O, 14N and 2H electric field gradient, EFG, as well as 17O, 15N, 13C and 1H chemical shielding tensors in the anhydrous chitosan crystalline structure. To include the hydrogen-bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a hexameric cluster. The computations were performed with the B3LYP method and 6-311++G(d,p) and 6-31++G(d,p) standard basis sets using the Gaussian 98 suite of programs. Calculated EFG and chemical shielding tensors were used to evaluate the 17O, 14N and 2H nuclear quadrupole resonance, NQR, and 17O, 15N, 13C and 1H nuclear magnetic resonance, NMR, parameters in the hexameric cluster, which are in good agreement with the available experimental data. The difference between the calculated NQR and NMR parameters of the monomer and hexamer cluster shows how much hydrogen bonding interactions affect the EFG and chemical shielding tensors of each nucleus. These results indicate that both O(3)-H(33)...O(5-3) and N-H(22)...O(6-4) hydrogen bonding have a major influence on NQR and NMR parameters. Also, the quantum chemical calculations indicate that the intra- and intermolecular hydrogen bonding interactions play an essential role in determining the relative orientation of EFG and chemical shielding principal components in the molecular frame axes.     

NQR and NMR study of hydrogen bonding interactions in anhydrous and monohydrated guanine cluster model: A computational study

In this paper extensive systematic computational study has been carried out to justify hydrogen bonding interactions and their influence on the oxygen, nitrogen and hydrogen NQR and NMR parameters of the anhydrous and monohydrated guanine crystal structures at two different levels, B3LYP and MP2, using 6-311++G** and D95** basis sets. These theoretical data have been compared with experimental NMR and NQR measurements. For further investigation, results of cluster calculation have been compared with that of a single molecule. Our theoretical NQR and NMR parameters of ¹⁷O, ¹⁵N and ²H atoms of anhydrous and monohydrated guanine exhibited extreme sensitivity to electron distribution around mentioned nuclei caused by cooperative influences of various types of hydrogen bonding interactions. Fortunately, our calculated isotropic shielding values and CS tensors for the ¹⁷O and ¹⁵N nuclei as well as obtained ¹⁴N-NQR parameters are in excellent agreement with experimental data. Therefore, we can undoubtedly conclude that for anhydrous and monohydrated guanine tetrameric clusters including intermolecular interactions, our theoretical estimates are in better agreement with observed experimental values than those in which these interactions have been ignored.     

Phonon-driven proton transfer in 3,5-pyridine dicarboxylic acid studied by 2H, 14N, and 17O nuclear quadrupole resonance

Hydrogen bonding in crystalline 3,5-pyridine dicarboxylic acid has been studied by (2)H, (14)N, and (17)O nuclear quadrupole resonance. The (2)H and (17)O data show the presence of two distinct hydrogen bonds, a "normal" O-H···O bond and a short, strong N···H···O bond, with significantly different NQR parameters. In the latter, the temperature variation of the (14)N nuclear quadrupole resonance (NQR) parameters is related to the phonon-driven proton transfer in the N···H···O hydrogen bond. The temperature dependence of the N···H and H···O distances in the N···H···O hydrogen bond is extracted from the (14)N NQR data.     

Investigation of the effects of ionic hydrogen bonds (COH–C and N–HOC) in crystalline dl-proline by ab initio and DFT calculated NQR parameters

In this paper, we have calculated the nuclear quadrupole resonance (NQR) parameters of the quadrupole nuclei involved in the hydrogen bonds (CO⁻H–C and ⁺N–H⁻OC) in the monomer and pentameric cluster of dl-proline by HF and B3LYP methods and basis sets of 6-311+G^* and 6-311++G^**. These computations are performed on the basis of X-ray diffraction structural data of dl-proline. The results indicate that the calculations including hydrogen-bonding (HB) interactions (in pentamer) are in better agreement with the experimental data than those in which these interactions are neglected (in monomer). The quantum chemical calculations show that the intermolecular hydrogen-bonding interactions play an important role in determination of the NQR parameters of ¹⁴N, ²H of group and ¹⁷O.     

A computational NQR study on the hydrogen-bonded lattice of cytosine-5-acetic acid

A computational study at the level of density functional theory (DFT) employing 6-311++G** standard basis set was carried out to evaluate nuclear quadrupole resonance (NQR) spectroscopy parameters in cytosine-5-acetic acid (C5AA). Since the electric field gradient (EFG) tensors are very sensitive to the electrostatic environment at the sites of quadruple nuclei, the most possible interacting molecules with the target one were considered in a five-molecule model system of C5AA using X-ray coordinates transforming. The hydrogen atoms positions were optimized and two model systems of original and H-optimized C5AA were considered in NQR calculations. The calculated EFG tensors at the sites of (17)O, (14)N, and (2)H nuclei were converted to their experimentally measurable parameters, quadrupole coupling constants and asymmetry parameters. The evaluated NQR parameters reveal that the nuclei in original and H-optimized systems contribute to different hydrogen bonding (HB) interaction. The comparison of calculated parameters between optimized isolated gas-phase and crystalline monomer also shows the relationship between the structural deformation and NQR parameters in C5AA. The basis set superposition error (BSSE) calculations yielded no significant errors for employed basis set in the evaluation of NQR parameters. All the calculations were performed by Gaussian 98 package of program.     

A study on the intermolecular hydrogen bonds of alpha-glycylglycine in its actual crystalline phase using ab initio calculated 14N and 2H nuclear quadrupole coupling constants

alpha-Glycylglycine in its actual crystalline phase is studied by ab initio calculated nuclear quadrupole coupling constants. These physical quantities are computed for 2H and 14N in the hydrogen bonds. The type of hydrogen bond is the N-H...O type. The computations are performed with the RHF and B3LYP methods and 6-31++G** and 6-311++G** basis sets using the Gaussian 98 program. Values of the calculated nuclear quadrupole coupling constants are shown in Tables 1-3. The aim of this work is the study of 2H and 14N quadrupole coupling constants which contribute in the CON2H...O=CN2H type of hydrogen bond. The computed nuclear quadrupole coupling constants of 2H nuclei meet the related experimental values. In addition, the computed chi value of 14N belonging to the -CO-14NH- group agrees well with values obtained experimentally. However, there are some discrepancies between calculated 14N chi values of the N+H3 residue and experiments. Also, the values of these physical parameters are calculated for >C2H2 of alpha-glycylglycine in its crystalline phase. Calculations for these parameters are carried out in a single molecule using X-ray diffraction coordinates, too.     

Hydrogen bond interactions in sulfamerazine: DFT study of the O-17, N-14, and H-2 electric field gradient tensors

Hydrogen bond (HB) interactions are studied in the real crystalline structure of sulfamerazine by density functional theory (DFT) calculations of the electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (single) and four-molecule (cluster) models of sulfamerazine are created by available crystal coordinates and the EFG tensors are calculated in both models to indicate the influence of HB interactions on the tensors. Directly relate to the experiments, the calculated EFG tensors are converted to the experimentally measurable nuclear quadrupole resonance (NQR) parameters, quadrupole coupling constant (qcc) and asymmetry parameter (η_Q). The evaluated NQR parameters reveal that due to contribution of the target molecule to N–HN and N–HO types of HB interactions, the EFG tensors at the sites of various nuclei are influenced from single model to the target molecule in cluster. Additionally, O2, N4, and H2 nuclei of the target molecule are significantly influenced by HB interactions, consequently, they have the major contributions to HB interactions in cluster model of sulfamerazine. The calculations are performed employing B3LYP method and 6-311++G** basis set using GAUSSIAN 98 suite of program.     

Influence of N-H...O and C-H...O hydrogen bonds on the 17O NMR tensors in crystalline uracil: computational study

We report a computational study for the 17O NMR tensors (electric field gradient and chemical shielding tensors) in crystalline uracil. We found that N-H...O and C-H...O hydrogen bonds around the uracil molecule in the crystal lattice have quite different influences on the 17O NMR tensors for the two C=O groups. The computed 17O NMR tensors on O4, which is involved in two strong N-H...O hydrogen bonds, show remarkable sensitivity toward the choice of cluster model, whereas the 17O NMR tensors on O2, which is involved in two weak C-H...O hydrogen bonds, show much smaller improvement when the cluster model includes the C-H...O hydrogen bonds. Our results demonstrate that it is important to have accurate hydrogen atom positions in the molecular models used for 17O NMR tensor calculations. In the absence of low-temperature neutron diffraction data, an effective way to generate reliable hydrogen atom positions in the molecular cluster model is to employ partial geometry optimization for hydrogen atom positions using a cluster model that includes all neighboring hydrogen-bonded molecules. Using an optimized seven-molecule model (a total of 84 atoms), we were able to reproduce the experimental 17O NMR tensors to a reasonably good degree of accuracy. However, we also found that the accuracy for the calculated 17O NMR tensors at O2 is not as good as that found for the corresponding tensors at O4. In particular, at the B3LYP/6-311++G(d,p) level of theory, the individual 17O chemical shielding tensor components differ by less than 10 and 30 ppm from the experimental values for O4 and O2, respectively. For the 17O quadrupole coupling constant, the calculated values differ by 0.30 and 0.87 MHz from the experimental values for O4 and O2, respectively.     

N nuclear quadrupole resonance of picolinic,nicotinic, isonicotinic and dinicotinic acids

The ¹⁴N nuclear quadrupole resonance (NQR) quadrupole coupling tensors of picolinic, nicotinic, isonicotinic and dinicotinic acids have been determined. Two different ¹⁴N quadrupole coupling constants 1007 kHz and 4159 kHz have been observed for picolinic acid demonstrating the presence of both protonated and non-protonated nitrogen atoms in this system in the solid. Only one set of non-protonated ¹⁴N NQR lines has been observed in other pyridinecarboxylic acids demonstrating the absence of the protonated zwitter ion forms observed in picolinic acid. The non-protonated ¹⁴N quadrupole coupling constant is the highest for the non-protonated nitrogen in picolinic acid and decreases to 3774 kHz in nicotinic acid and 3570 kHz in isonicotinic acid. It is the lowest in dinicotinic acid where the corresponding ¹⁴N quadrupole coupling constant is 2794 kHz. The observed anomalous decrease in the ¹⁴N quadrupole coupling constant of dinicotinic acid with decreasing temperature is tentatively explained as reflecting the increase in the residence time of the N–H?O bonded proton in the potential well close to the nitrogen.     

DFT calculations of <Superscript>14</Superscript>N, <Superscript>17</Superscript>O,and <Superscript>2</Superscript>H NQR parameters: Investigation of hydrogen bond interactions in sulfapyridine crystalline structure

DFT calculations of electric field gradient (EFG) tensors at the sites of ¹⁴N, ¹⁷O, and ²H nuclei are carried out to characterize the hydrogen bond (HB) interactions in the sulfapyridine crystal structure. One-molecule (monomer) and hydrogen-bonded hexameric cluster models of sulfapyridine are constructed according to available X-ray coordinates where the proton positions are optimized. Then, EFG tensors are calculated for both monomer and target molecule in the hexameric cluster of sulfapyridine to show the effect of HB interactions on the tensors. The calculated EFG tensors are converted to the experimentally measurable nuclear quadrupole resonance (NQR) parameters: quadrupole coupling constant (C _Q ) and asymmetry parameter (η _Q ). The results reveal different contribution of various nuclei to N-H⋯N and N-H⋯O HB interactions in the cluster where the N2 and O1 have major contributions. The computations are performed with B3LYP and B3PW91 functionals DFT method and 6-311+G* and 6-311++G** standard basis sets using the Gaussian 98 package.     

Combined first-principles computational and experimental multinuclear solid-state NMR investigation of amino acids

13C, 14N, 15N, 17O, and 35Cl NMR parameters, including chemical shift tensors and quadrupolar tensors for 14N, 17O, and 35Cl, are calculated for the crystalline forms of various amino acids under periodic boundary conditions and complemented by experiment where necessary. The 13C shift tensors and 14N electric field gradient (EFG) tensors are in excellent agreement with experiment. Similarly, static 17O NMR spectra could be precisely simulated using the calculation of the full chemical shift (CS) tensors and their relative orientation with the EFG tensors. This study allows correlations to be found between hydrogen bonding in the crystal structures and the 17O NMR shielding parameters and the 35Cl quadrupolar parameters, respectively. Calculations using the two experimental structures for L-alanine have shown that, while the calculated isotropic chemical shift values of 13C and 15N are relatively insensitive to small differences in the experimental structure, the 17O shift is markedly affected.     

Unusual electron charge density in carboxylic acid. 17O quadrupole coupling in cis-cyclobutane-1,2-dicarboxylic acid

The (17)O NQR frequencies have been measured in cis-cyclobutane-1,2-dicarboxylic acid and the quadrupole coupling tensors have been determined at various temperatures. Two O···H oxygen positions and two O-H oxygen positions are observed, showing the presence of two different types of O-H···O hydrogen bonds in the unit cell. The quadrupole coupling constants at the O-H oxygen positions are approximately 30% lower than the lowest quadrupole coupling constants experimentally observed at the C-O-H positions in other carboxylic acids with either ordered or disordered hydrogen bonds. The O-H distances have been calculated from the (17)O-(1)H dipole-dipole interaction at the O-H oxygen positions. The obtained values are longer than the O-H distances usually found in O-H···O hydrogen bonds with comparable O···O distance, in agreement with the proposed proton exchange O-H···O ? O···H-O, which partially averages the dipole-dipole interaction. The energy difference of the two proton configurations, O-H···O and O···H-O, is calculated from the O-H distances determined by NQR. The temperature dependence of the (17)O quadrupole coupling tensors at the (17)O···H-O oxygen positions is analyzed in the model of proton exchange and the energy differences of the two proton configurations obtained by this analysis agree with the values obtained from the O-H distances. The quadrupole coupling tensors are analyzed in a model based on the Townes and Dailey model. The model shows that the population of an oxygen lone pair orbital is at this oxygen position reduced from 2 to approximately 1.3. The electron electric charge is most probably transferred to the oxygen σ and π electron orbitals. This may be associated with the structure of the cyclobutane ring, where the X-ray data show the presence of two unusually short C-C bonds.     

Study of hydrogen bonds in N-methylacetamide by DFT calculations of oxygen, nitrogen, and hydrogen solid-state NMR parameters

Solid-state nuclear magnetic resonance (NMR) parameters of ¹⁷O, ¹⁴N/¹⁵N, and ²H/¹H nuclei were evaluated in two available neutron crystalline structures of N-methylacetamide (NMA) at 250 and 276 K, NMA-I and NMA-II, respectively. Density functional theory calculations were performed by B3LYP method and 6-311++G** and IGLO-II type basis sets to calculate the electric field gradient (EFG) and chemical shielding (CS) tensors at the sites of mentioned nuclei. In order to investigate hydrogen bonds (HBs) effects on NMR tensors, calculations were performed on four-model systems of NMA: an optimized isolated gas-phase, crystalline monomers, crystalline dimers, and crystalline trimers. Comparing the calculated results reveal the influence of N–H···O=C and C–H···O=C HB types on the NMR tensors which are observable by the evaluated parameters including quadrupole coupling constant, C _Q, and isotropic CS, σ _iso. Furthermore, the results demonstrate more influence of HB on the EFG and CS tensors of NMA at 276 K rather than that of 250 K.     

B24N24 nanocages: a GIAO density functional theory study of 14N and 11B nuclear magnetic shielding and electric field gradient tensors

Abstract  

Density functional theory (DFT) calculations were performed to determine boron-11 and nitrogen-14 nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) spectroscopy parameters in the three most stable B₂₄N₂₄ fullerenes for the first time. The considered samples were first allowed to relax entirely, and then the NMR and NQR calculations were performed on the geometrically optimized models. The calculations of the ¹¹B and ¹⁴N nuclear magnetic shielding tensors and electric field gradient tensors employed the Gaussian 98 software implementation of the gauge-including atomic orbital (GIAO) method using the Becke3, Lee-Yang-Parr (B3LYP) DFT level and 6-311G** and 6-311++G** standard basis sets in each of the three optimized forms, and converted the results to experimentally measurable NMR parameters.The calculated NMR chemical shieldings of the three cages show significant differences, providing a way to identify these clusters. The evaluated NQR parameters of the ¹¹B and ¹⁴N nuclei in the clusters are also reported and discussed.     

Calculation of the 14N and 2H quadrupole coupling tensor in trimeric imidazole

The ¹⁴N and ²H quadrupole coupling tensors at both nitrogen and at all hydrogen sites, have been calculated by ab initio SCF MO methods, in a trimeric imidazole molecule having the geometry found in the solid state. The ¹⁴N changes produced at N(1)H by hydrogen bonding are correctly reproduced to a noteworthy degree of precision, as is the ²H quadrupole interaction at the same site. The changes in the values of the three principal components of the tensor, which have been inferred from NQR measurements when imidazole passes from the gas to the solid state, are supported by the calculations, both for N(1)H and N(3).     

Theoretical 14N and 17O nuclear quadrupole resonance parameters for tirapazamine and related metabolites

Density functional theory (DFT) calculations were performed to realize the effects of the N–O group on the reactivity and electronic properties of 3-amino-1,2,4-benzotriazines. The electric field gradient, EFG, tensors of ¹⁴N and ¹⁷O nuclei and natural bond orbital (NBO) analysis in the tirapazamine (TPZ) and its four derivatives were calculated at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) method in the gas phase. The NBO analysis reveal that the bond strength, proton affinity and position of N–O group in the heterocyclic ring have major influence on the reactivity of considered molecules. Accordingly, we suggest that the TPZ and 4-oxide (d) structures due to having a weaker N–O bond with larger negative charge on the oxygen atom at the 4-position are more active than the other ones. Calculated ¹⁴N and ¹⁷O EFG tensors were used to evaluate nuclear quadrupole coupling tensors, χ, and asymmetry parameters, η _Q . Results showed that oxidation of a nitrogen atom at any position have significant influence on its ¹⁴N nuclear quadrupole resonance (NQR) parameters. Also, the occupancy of nitrogen lone pair plays an important role in determination of the q _zz and χ values of mentioned nuclei. It is found that the η _Q and χ are appropriate parameters to study the contribution of lone pair electrons of nitrogen atom in the formation of chemical bond or conjugation with the aromatic system. Finally, a linear correlation is observed between the χ(¹⁴N) and χ(¹⁷O) values in the N–O bond which may be associated with the reactivity of these compounds.     

A theoretical investigation of hydrogen bonding effects on oxygen and hydrogen chemical shielding tensors of aspirin

A computational investigation was carried out to characterize the ¹⁷O and ¹H chemical shielding (CS) tensors in crystalline aspirin. It was found that O–H⋯O and C–H⋯O hydrogen bonds around the aspirin molecule in the crystal lattice have a different influence on the calculated ¹⁷O and ¹H CS eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the BLYP, B3LYP, and M06 functionals employing 6-311++G(d,p) standard basis set. Calculated CS tensors were used to evaluate the ¹⁷O and ¹H chemical shift isotropy (δ_iso) and anisotropy (Δσ) in crystalline aspirin, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the CS tensors of each nucleus.     

Characterization of intermolecular interactions in crystalline aspirin: A computational NQR study

A computational study at the level of density functional theory was carried out to characterize the ¹⁷O and ²H nuclear quadrupole resonance (NQR) spectroscopy parameters in crystalline aspirin. To include O? H···O and C? H···O hydrogen bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a pentamer cluster. The NQR calculations were performed with BLYP, B3LYP, and M06 functionals employing 6‐311++G** and Jensen's polarization‐consistent pcJ‐1 basis sets. Linear correlations are observed between the calculated ¹⁷O and ²H NQR parameters and the hydrogen bond strengths, suggesting the possibility of estimating hydrogen bonding information from calculated NQR data. Different contributions of various nuclei to hydrogen bonding interactions and observed trends of calculated NQR parameters are well justified by atoms in molecules analyses at the BCPs of these interactions. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

14N NQR study of nicotinamide and related compounds

14N nuclear quadrupole resonance (NQR) frequencies have been measured in picolinamide, nicotinamide, isonicotinamide, 2,6-pyridine dicarboxamide, and acetamide by double resonance. The 14N NQR spectra in picolinamide, nicotinamide, isonicotinamide, and 2,6-pyridine dicarboxamide show the presence of two distinct nitrogen positions: the ring position with the quadrupole coupling constant about 4,5 MHz and the amide position with the quadrupole coupling constant about 2.6 MHz. The NQR data are related to the structure of the investigated compounds and to the N--H...O hydrogen bonds.     

Microsolvation of formamide: a rotational study

Microsolvated formamide clusters have been generated in a supersonic jet expansion and characterized using Fourier transform microwave spectroscopy. Three conformers of the monohydrated cluster and one of the dihydrated complex have been observed. Seven monosubstituted isotopic species have been measured for the most stable conformer of formamide...H(2)O, which adopts a closed planar ring structure stabilized by two intermolecular hydrogen bonds (N-H...O(H)-H...O=C). The two higher energy forms of formamide...H(2)O have been observed for the first time. The second most stable conformer is stabilized by a O-H...O=C and a weak C-H...O hydrogen bond, while, in the less stable form, water accepts a hydrogen bond from the anti hydrogen of the amino group. For formamide...(H(2)O)(2), the parent and nine monosubstituted isotopic species have been observed. In this cluster the two water molecules close a cycle with the amide group through three intermolecular hydrogen bonds (N-H...O(H)-H...O(H)-H...O=C), the nonbonded hydrogen atoms of water adopting an up-down configuration. Substitution (r(s)) and effective (r(0)) structures have been determined for formamide, the most stable form of formamide...H(2)O and formamide...(H(2)O)(2). The results on monohydrated formamide clusters can help to explain the observed preferences of bound water in proteins. Clear evidence of sigma-bond cooperativity effects emerges when comparing the structures of the mono- and dihydrated formamide clusters. No detectable structural changes due to pi-bond cooperativity are observed on formamide upon hydration.