NMR crystallography of zeolites: refinement of an NMR-solved crystal structure using ab initio calculations of 29Si chemical shift tensors

An NMR structure refinement method for the NMR crystallography of zeolites is presented and demonstrated to give an NMR-determined crystal structure for the zeolite Sigma-2 that is in very good agreement with the single-crystal X-ray diffraction structure. The Si coordinates of the zeolite framework were solved from 29Si double-quantum NMR data obtained at a low magnetic field strength (7.0 T) and the Si and O coordinates were subsequently refined using the principal components of 29Si chemical shift tensors experimentally measured at an ultrahigh-field (21.1 T) and calculated using ab initio quantum chemical methods.     

Organic peracids: a structural puzzle for 17O NMR and ab initio chemical shift calculations

We have applied (17)O NMR spectroscopy to investigate the structure of the organic peracids formed by reaction of acetic acid (AA) or lactic acid (LA) with aqueous hydrogen peroxide (HP), which are used in several "green chemistry" applications. The interpretation of the experimental spectra has been supported by ab initio calculations of the (17)O chemical shifts for several possible species, using a continuum representation of the solvent. The combined use of these tools has also allowed us to discuss the decomposition mechanism of LA/HP solutions. The calculated electric field gradients for water, HP, and CO(2) (a decomposition product of LA) correlate well with the experimental (17)O line widths.     

Combined application of 2D NMR correlation methods and ab initio chemical shift calculations to the structure determination of new heterocyclic compounds

The combined use of 2D NMR correlation methods and ab initio chemical shift calculations is efficient and, in some cases, virtually the only way to determine the structures of new organic compounds. This approach enabled us to establish the structure of the major unusual product of the three-component reaction of imidazo[1,5-a]quinoxalin-4-one, bis(2-chloroethyl)amine hydrochloride, and potassium carbonate in DMF. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2172–2179, December, 2006. Dedicated to Professor A. V. Il’yasov on the occasion of his 70th birthday.     

Correlations between 31P chemical shift anisotropy and molecular structure in polycrystalline O,O'-dialkyldithiophosphate zinc(II) and nickel(II) complexes: 31P CP/MAS NMR and ab initio quantum mechanical calculation studies

Different potassium salts and zinc(II) and nickel(II) O,O'-dialkyldithiophosphate complexes were studied by solid-state 31P CP/MAS and static NMR and ab initio quantum mechanical calculations. Spectra were obtained at different spinning frequencies, and the intensities of the spinning sidebands were used to estimate the chemical shift anisotropy parameters. Useful correlations between the shapes of the 31P chemical shift tensor and the type of ligand were found: terminal ligands have negative values of the skew kappa, while bridging and ionic ligands have positive values for this parameter. The experimental results were compared with known X-ray diffraction structures for some of these complexes as well as with ab initio quantum mechanical calculations, and a useful correlation between the delta22 component of the 31P chemical shift tensor and the S-P-S bond angle in the O,O'-dialkyldithiophoshate zinc(II) and nickel(II) complexes was found: delta22 increases more than 50 ppm with the increase of S-P-S bond angle from ca. 100 degrees to 120 degrees , while the other two principal values of the tensor, delta11 and delta33, are almost conserved. This eventually leads to the change in sign for kappa in the bridging type of ligand, which generally has a larger S-P-S bond angle than the terminally bound O,O'-dialkyldithiophosphate group forming chelating four-membered P(ss)Me heterocycles.     

Alkaline earth chloride hydrates: chlorine quadrupolar and chemical shift tensors by solid-state NMR spectroscopy and plane wave pseudopotential calculations

A series of alkaline earth chloride hydrates has been studied by solid-state (35/37)Cl NMR spectroscopy in order to characterize the chlorine electric field gradient (EFG) and chemical shift (CS) tensors and to relate these observables to the structure around the chloride ions. Chlorine-35/37 NMR spectra of solid powdered samples of pseudopolymorphs (hydrates) of magnesium chloride (MgCl(2).6H(2)O), calcium chloride (CaCl(2).2H(2)O), strontium chloride (SrCl(2), SrCl(2).2H(2)O, and SrCl(2).6H(2)O), and barium chloride (BaCl(2).2H(2)O) have been acquired under stationary and magic-angle spinning conditions in magnetic fields of 11.75 and 21.1 T. Powder X-ray diffraction was used as an additional tool to confirm the purity and identity of the samples. Chlorine-35 quadrupolar coupling constants (C(Q)) range from essentially zero in cubic anhydrous SrCl(2) to 4.26+/-0.03 MHz in calcium chloride dihydrate. CS tensor spans, Omega, are between 40 and 72 ppm, for example, Omega= 45+/-20 ppm for SrCl(2).6H(2)O. Plane wave-pseudopotential density functional theory, as implemented in the CASTEP program, was employed to model the extended solid lattices of these materials for the calculation of their chlorine EFG and nuclear magnetic shielding tensors, and allowed for the assignment of the two-site chlorine NMR spectra of barium chloride dihydrate. This work builds upon our current understanding of the relationship between chlorine NMR interaction tensors and the local molecular and electronic structure, and highlights the particular sensitivity of quadrupolar nucleus solid-state NMR spectroscopy to the differences between various pseudopolymorphic structures in the case of strontium chloride.     

31P chemical shift tensors for canonical and non-canonical conformations of nucleic acids: a DFT study and NMR implications

31P chemical shift anisotropy (CSA) tensors have been calculated for a set of selected DNA and RNA backbone conformations using density functional theory. The set includes canonical A-RNA, A-DNA, BI-DNA, BII-DNA, ZI-DNA, and ZII-DNA as well as four A-RNA-type, seven non-A-RNA-type, and three non-canonical DNA conformations. Hexahydrated dimethyl phosphate has been employed as a model. The 31P chemical shift tensors obtained are discussed in terms of similarities in the behavior observed for gauche-gauche (gg) and gauche-trans (gt) conformations around the P-O bonds. We show that torsion angles alpha and zeta are major determinants of the isotropic chemical shift deltaiso and of the deltaCSA11 component of the traceless chemical shift tensor, which is revealed in separate ranges of both deltaiso and deltaCSA11 for gg- and gt-conformers, respectively. A clear distinction between the two conformation types has not been found for the deltaCSA22 and deltaCSA33 components, which is attributed to their different directional properties. The 31P CSA tensors exhibit considerable variations resulting in large spans of approximately 16 ppm for deltaCSA11 and approximately 22 ppm for deltaCSA22 and deltaCSA33. We examine the consequences of the CSA variations for predicting the chemical shift changes upon partial alignment deltacsa and for the values of CSA order parameters extracted from the analysis of 31P NMR relaxation data. The theoretical 31P CSA tensors as well as the experimental 31P CSA tensor of barium diethyl phosphate (BDEP) are used to calculate deltacsa for two eclipsed orientations of the CSA and molecular alignment tensors. Percentage differences between the CSA order parameters obtained using the theoretical 31P CSA tensors and the experimental 31P CSA tensor of BDEP, respectively, are also determined.     

Cluster models and ab initio calculations of (19)F NMR isotropic chemical shifts for inorganic fluorides

(19)F NMR isotropic chemical shift (delta(iso)) calculations are performed in crystallized compounds using the GIAO method with the B3LYP hybrid functional at DFT level. Clusters centered on the studied fluorine atoms mimic the crystalline structures. The 6-311+G(d) basis set is chosen for the central fluorine atom, and the LanL2DZ basis set for the others. The metal atoms are described by the 3-21G(2d) basis set or, when not available, by the CRENBL basis set with the corresponding ECP, and augmented with 2d polarization functions when existing. First, for high-symmetry systems (MF, MF(2), and MF(3) compounds), a systematization of the cluster building up from coordination spheres is proposed, generalized to fluoroperovskites and fluoroaluminates KAlF(4) and RbAlF(4). When applied to rather low symmetry systems such as barium fluorometalates BaMgF(4), BaZnF(4), and Ba(2)ZnF(6), the definition of the coordination spheres is far from easy. Then, for structures built up from a MF(6) octahedron network, we may define different "starting clusters": [FM(2)F(8)] for the shared fluorine atoms, [FMF(4)] for the unshared ones, and [FBa(4)](7+) for the "free" ones. Analogous "starting clusters" are then tested on compounds from the NaF-AlF(3), BaF(2)-AlF(3), and CaF(2)-AlF(3) binary systems and for alpha-BaCaAlF(7) that are also built up from a MF(6) octahedron network. For each of these corresponding fluorine sites, delta(iso) values are calculated with the "starting clusters" and several larger clusters and compared to the experimental delta(iso) values. For the barium-containing clusters, the RMS deviation is equal to 51 ppm. It is suggested that this result may be related to the poor quality of the barium basis sets for which no polarization functions are available for the moment. In total, chemical shifts were calculated for 122 fluorine sites, in a various range of compounds. For the clusters without barium, the ab initio method leads to a RMS equal to 22 ppm, which is a quite nice result keeping in mind that the (19)F chemical shift range is larger than 200 ppm.     

Ab initio study of (13)C(alpha) chemical shift anisotropy tensors in peptides

This study reports magnitudes and the orientation of the (13)C(alpha) chemical shift anisotropy (CSA) tensors of peptides obtained using quantum chemical calculations. The dependency of the CSA tensor parameters on the energy optimization of hydrogen atom positions and hydrogen bonding effects and the use of zwitterionic peptides in the calculations are examined. Our results indicate that the energy optimization of the hydrogen atom positions in crystal structures is necessary to obtain accurate CSA tensors. The inclusion of intermolecular effects such as hydrogen bonding in the calculations provided better agreement between the calculated and experimental values; however, the use of zwitterionic peptides in calculations, with or without the inclusion of hydrogen bonding, did not improve the results. In addition, our calculated values are in good agreement with tensor values obtained from solid-state NMR experiments on glycine-containing tripeptides. In the case of peptides containing an aromatic residue, calculations on an isolated peptide yielded more accurate isotropic shift values than the calculations on extended structures of the peptide. The calculations also suggested that the presence of an aromatic ring in the extended crystal peptide structure influences the magnitude of the delta(22) which the present level of ab initio calculations are unable to reproduce.     

Solid-state (13)C NMR chemical shift anisotropy tensors of polypeptides.

Carbon-13 chemical shift anisotropy (CSA) tensors for various carbon sites of polypeptides, and for carbon sites in alpha-helical and beta-sheet conformations of poly-L-alanine, and polyglycine, are presented. The carbonyl (13)C CSA tensors were determined from one-dimensional CPMAS spectra obtained at a slow spinning speed, whereas the CSA tensors of C(alpha) and other carbons in side chains of peptides were determined using 2D PASS experiments on powder samples. The results suggest that the spans of (13)Carbonyl CSA tensors of alanine and glycine residues in various peptides are similar, even though the magnitude of individual components of the CSA tensor and the isotropic chemical shift are different. In addition, the delta(22) element is the only component of the (13)Carbonyl CSA tensor that significantly depends on the CO.HN hydrogen-bond length. Solid-state NMR experimental results also suggest that (13)Carbonyl and (13)C(alpha) CSA tensors are similar for alpha-helical and beta-sheet conformations of poly-L-alanine, which is in agreement with the reported quantum chemical calculation studies and previous solid-state NMR experimental studies on other systems. On the other hand, the (13)C(alpha) CSA tensor of the first alanine residue is entirely different from that of the second or later alanine residues of the peptide. While no clear trends in terms of the span and the anisotropic parameter were predicted for (13)C(beta) CSA tensors of alanine, they mainly depend on the conformation and dynamics of the side chain as well as on the packing interactions in the solid state of peptides.     

Determination of chemical shift anisotropy tensors of carbonyl nuclei in proteins through cross-correlated relaxation in NMR.

The principal components and orientations of the chemical shift anisotropy (CSA) tensors of nearly all 13C carbonyl nuclei in a small protein have been determined in isotropic solution by a combination of three complementary cross-correlation measurements.     

DFT calculations of 31P spin-spin coupling constants and chemical shift in dioxaphosphorinanes

One-bond heteronuclear spin-spin coupling constants (1)J(PX) (X=H, O, S, Se, C and N) between the phosphorus atom and axial and equatorial substituents in dioxaphosphorinanes are computed using density functional theory (DFT). The experimental values of these coupling constants for a variety of substituents can be applied to identify different diastereoisomers. The DFT calculations confirm the systematic trend observed in experiment, and indicate that the computed (1)J(PX) coupling constants are related to the length of the axial and equatorial bonds. A similar relation between the phosphorus chemical shift and the R(PX) bond length appears to be valid, with the exception of selenium substituents.     

Interaction of Cd and Zn with biologically important ligands characterized using solid-state NMR and ab initio calculations

For the first time, coordination geometry and structure of metal binding sites in biologically relevant systems are studied using chemical shift parameters obtained from solid-state NMR experiments and quantum chemical calculations. It is also the first extensive report looking at metal-imidazole interaction in the solid state. The principal values of the (113)Cd chemical shift anisotropy (CSA) tensor in crystalline cadmium histidinate and two different cadmium formates (hydrate and anhydrate) were experimentally measured to understand the effect of coordination number and geometry on (113)Cd CSA. Further, (13)C and (15)N chemical shifts have also been experimentally determined to examine the influence of cadmium on the chemical shifts of (15)N and (13)C nuclei present near the metal site in the cadmium-histidine complex. These values were then compared with the chemical shift values obtained from the isostructural bis(histidinato)zinc(II) complex as well as from the unbound histidine. The results show that the isotropic chemical shift values of the carboxyl carbons shift downfield and those of amino and imidazolic nitrogens shift upfield in the metal (Zn,Cd)-histidine complexes relative to the values of the unbound histidine sample. These shifts are in correspondence with the anticipated values based on the crystal structure. Ab initio calculations on the cadmium histidinate molecule show good agreement with the (113)Cd CSA tensors determined from solid-state NMR experiments on powder samples. (15)N chemical shifts for other model complexes, namely, zinc glycinate and zinc hexaimidazole chloride, are also considered to comprehend the effect of zinc binding on (15)N chemical shifts.     

Understanding solvation in hydrofluoroalkanes: ab initio calculations and chemical force microscopy

Understanding solvation in hydrofluoroalkane (HFA) propellants is of great importance for the development of novel pressurized metered-dose inhaler (pMDI) formulations. HFA-based pMDIs are not only the most widely used inhalation therapy devices for treating lung diseases, but they also hold promise as vehicles for the systemic delivery of biomolecules to and through the lungs. In this work we propose a combined microscopic experimental and computational approach to quantitatively relate the chemistry of moieties to their HFA-philicity. Binding energy calculations are used to determine the degree of interaction between a propellant HFA and candidate fragments. We define a new quantity, the enhancement factor E, which also takes into account fragment-fragment interactions. This quantity is expected to correlate well with the solubility and the ability of the moieties of interest to impart stability to colloidal dispersions in HFAs. We use a methyl-based (CH) segment and its fluorinated analog (CF) to test our approach. CH is an important baseline case since it represents the tails of surfactants in FDA-approved pMDIs. CF was chosen due to the improved solubility and ability of this chemistry to stabilize aqueous dispersions in HFAs. Adhesion force from Chemical Force Microscopy (CFM) is used as an experimental analog to the binding energy calculations. The force of interaction between a chemically modified AFM tip and substrate is measured in a model HFA, which is a liquid at ambient conditions. Silanes with the same chemistry as the fragments used in the ab initio calculations allow for direct comparison between the two techniques. The CFM results provide an absolute scale for HFA-philicity. Single molecule (pair) forces calculated from the CFM experiments are shown to be in very good agreement to the E determined from the ab initio calculations. The ab initio calculations and CFM are corroborated by previous experimental studies where propellants HFAs are seen to better solvate the CF functionality.     

Relationships between 31P chemical shift tensors and conformation of nucleic acid backbone: a DFT study

Density functional theory (DFT) has been applied to study the conformational dependence of 31P chemical shift tensors in B-DNA. The gg and gt conformations of backbone phosphate groups representing BI- and BII-DNA have been examined. Calculations have been carried out on static models of dimethyl phosphate (dmp) and dinucleoside-3',5'-monophosphate with bases replaced by hydrogen atoms in vacuo as well as in an explicit solvent. Trends in 31P chemical shift anisotropy (CSA) tensors with respect to the backbone torsion angles alpha, zeta, beta, and epsilon are presented. Although these trends do not change qualitatively upon solvation, quantitative changes result in the reduction of the chemical shift anisotropy. For alpha and zeta in the range from 270 degrees to 330 degrees and from 240 degrees to 300 degrees , respectively, the delta22 and delta33 principal components vary within as much as 30 ppm, showing a marked dependence on backbone conformation. The calculated 31P chemical shift tensor principal axes deviate from the axes of O-P-O bond angles by at most 5 degrees . For solvent models, our results are in a good agreement with experimental estimates of relative gg and gt isotropic chemical shifts. Solvation also brings the theoretical deltaiso of the gg conformation closer to the experimental gg data of barium diethyl phosphate.     

On phosphorus chemical shift anisotropy in metal (IV) phosphates and phosphonates: A 31P NMR experimental study

The phosphorus chemical shift anisotropies, ³¹PΔcs, and asymmetry parameters η were measured by the ³¹P{¹H} NMR experiments in static and low-frequency spinning samples of the zirconium phosphates and phosphonates and also in the mixed Zr (IV)/Sn (IV) phosphate/phosphonate material. The data obtained have shown a 111 connectivity in the HPO₄ and PO₃ groups, which does not change at modification and intercalation of the materials. The ³¹PΔcs values of the phosphonate groups (43–49 ppm) significantly surpass the values characterizing the HPO₄ groups (23–37 ppm). The ³¹P Δcs values obtained for the metal (IV) phosphates were discussed in terms of P-O distances. The ³¹P chemical shift anisotropy parameters can help at elucidation of local structures in phosphate and phosphonate materials.     

Determination of 31P chemical shielding tensors from NMR studies on single crystals of chlorotris(triphenylphosphine) rhodium(I)

The ³¹P chemical shielding tensors, and their direction cosines, of the three chemically different phosphorus nuclei in chlorotris(triphenylphosphine) rhodium(I) were determined in single-crystal studies using high-resolution cross-polarization NMR with proton decoupling. The tensors are not axially symmetric and the directions of the most shielded principal values for each of the three phosphorus nuclei are almost parallel to the corresponding PRh bonds.     

Hydrogen bonding between phenols and fatty acid esters: (1)H NMR study and ab initio calculations

[structure: see text] (1)H NMR measurements and ab initio calculations were used to study the interactions between hindered/nonhindered phenols and carboxylic acid esters. The dihedral angle (phi) between the OH group and a plane of the aromatic ring is close to 0 degrees in the hydrogen-bonded nonhindered phenols, whereas for 2,6-di-tert-butyl-4-methylphenol the OH group is completely twisted out of the aromatic plane (phi approximately 90 degrees ).     

Photoionization of C(4) molecular beam: ab initio calculations

Large computations are performed on the C(4) (+) cation in order to characterize its stable isomers and its lowest electronic excited states using configuration interaction methods and large basis sets. Several stable isomers are found including a linear C(4) (+)(l-C(4) (+)), a rhombic C(4) (+)(r-C(4) (+)) (or cyclic), and a branched (d-C(4) (+)) structure. Our calculations show a high density of electronic states for all of these isomers favoring their interactions. By combining the present ab initio data and those on neutral C(4), the l-C(4)(X)+hnu-->l-C(4) (+)(X(+))+e(-), d-C(4)(X)+hnu-->d-C(4) (+)(X(+))+e(-), and r-C(4)(X)+hnu-->r-C(4) (+)(X(+))+e(-) vertical photoionization transition energies are computed at 10.87, 10.92, and 10.77 eV, respectively. Photoionizing a C(4) molecular beam results on an onset at 10.4-10.5 eV and then to a linear increase of the signal due to the opening of several ionization channels involving most of the C(4) and C(4) (+) isomers and electronic states.     

Germanium(II)-mediated reductive cross-aldol reaction of bromoaldehydes with aldehydes: NMR studies and ab initio calculations

A highly practical reductive cross-aldol reaction of alpha-bromoaldehydes with various aldehydes has been developed using Ge(II)Cl 2 to produce aldehyde germanium(IV) aldolates, which were directly transformed to various multifunctionalized compounds. A remarkable change in stereoselectivity depended on the alpha-bromoaldehydes employed; secondary alpha-bromoaldehydes gave syn selectivities, while tertiary alpha-bromoaldehydes accomplished the synthesis of anti-selective aldol products with a quaternary carbon center. NMR studies and X-ray analysis strongly suggested the formation of germanium enolate in the reaction of alpha-bromoaldehyde 2h with GeCl 2-dioxane. Detailed mechanistic studies, including NMR analysis and ab initio calculations, revealed the generation of stable germanium aldolates, which was due to the remarkably low Lewis acidity of the germanium(IV).