Calculations of chemical shieldings: Theory and applications

Confined Lewis Pairs: Investigation of the X- →Si20 Interaction in Halogen-Encapsulating Silafulleranes

A joint theoretical and experimental study on 32 endohedral silafullerane derivatives [X@Si20 Y20 ]- (X=F-I; Y=F-I, H, Me, Et) andT h ${T_h }$ -[Cl@Si20 H12 Y8 ]- (Y=F-I) is presented. First, we evaluated the structure-determining template effect of Cl- in a systematic series of concave silapolyquinane model systems. Second, we investigated the X- →Si20 interaction energy (E int ${E_{{\rm{int}}} }$ ) as a function of X- and Y and found the largestE int ${E_{{\rm{int}}} }$ values for electron-withdrawing exohedral substituents Y. Given that X- ions can be considered as Lewis bases and empty Si20 Y20 clusters as Lewis acids, we classify our inseparable host-guest complexes [X@Si20 Y20 ]- as "confined Lewis pairs". Third, 35 Cl NMR spectroscopy proved to be highly diagnostic for an experimental assessment of the Cl- →Si20 interaction as the paramagnetic shielding and, in turn, δ ${\delta }$ (35 Cl) of the endohedral Cl- ion correlate inversely withE int ${E_{{\rm{int}}} }$ . Finally, we disclose the synthesis of [PPN][Cl@Si20 Y20 ] (Y=Me, Et, Br) and provide a thorough characterization of these new silafulleranes.

Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems

We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFμ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fμ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.

DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking

Density functional theory (DFT) benchmark studies of ¹H and ¹³C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured ¹H and ¹³C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for ¹H and ¹³C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for ¹H and 0.5 to 2.9 for ¹³C. Maximum deviations were less than 0.5 and 6.5 ppm for ¹H and ¹³C, respectively.

Relativistic Density Functional NMR Tensors Analyzed with Spin-free Localized Molecular Orbitals

The implementation of fast relativistic methods based on density functional theory, in conjunction with localized molecular orbital (LMO) based analysis, allows straightforward interpretations of NMR parameters in terms of contributions from core shells, lone pairs, and bonds, for compounds containing elements from across the periodic table. We present a conceptual review of a frequently used LMO analysis of NMR parameters calculated in the presence of spin-orbit interactions and other relativistic effects. An accompanying example focuses on the ¹⁵ N shielding in a heavy metal complex.

Spatial Contributions to 1H NMR Chemical Shifts of Free-Base Porphyrinoids

A recently developed methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities is used for studying spatial contributions including ring-current contributions to 1H nuclear magnetic resonance (NMR) chemical shifts of aromatic and anti-aromatic free-base porphyrinoids. Our approach allows a visual inspection of the spatial origin of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constants. Diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions implying that not merely the tropicity of the current density determines whether the contribution has a shielding or deshielding character. Instead the shielding or deshielding contribution is determined by the direction of the current-density flux with respect to the studied nucleus.

[Cl@Si20H20]-: Parent Siladodecahedrane with Endohedral Chloride Ion

Fullerenes and diamondoids are at the core of nanoscience. Comparable monodisperse silicon analogues are scarce. Herein, we report the synthesis of the parent siladodecahedrane, which represents the largest Platonic solid. It shares its pattern of pentagonal faces with the smallest fullerene, C₂₀, and its saturated, H-terminated skeleton with diamondoids. Similar to endofullerenes, the silicon cage encapsulates a chloride ion ([Cl@Si₂₀H₂₀]^-); similar to diamondoids, its Si-H termini offer a wealth of opportunities for further functionalization. Mere treatment with chloromethanes leads to the perchlorinated cluster [Cl@Si₂₀Cl₂₀]^-. Both compounds were characterized by mass spectrometry, X-ray crystallography, NMR spectroscopy, and quantum-chemical calculations. The experimentally determined ³⁵Cl resonances of the endohedral chloride ions are particularly diagnostic to probe the Cl^- → Si₂₀ interaction strength as a function of the different surface substituents, as we have proven by high-level computational analyses.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Spatial Contributions to Nuclear Magnetic Shieldings

We develop a methodology for calculating, analyzing, and visualizing nuclear magnetic shielding densities which are calculated from the current density via the Biot-Savart relation. Atomic contributions to nuclear magnetic shielding constants can be estimated within our framework with a Becke partitioning scheme. The new features have been implemented in the GIMIC program and are applied in this work to the study of the 1H and 13C nuclear magnetic shieldings in benzene (C6H6) and cyclobutadiene (C4H4). The new methodology allows a visual inspection of the spatial origins of the positive (shielding) and negative (deshielding) contributions to the nuclear magnetic shielding constant of a single nucleus, something which has not been hitherto easily accomplished. Analysis of the shielding densities shows that diatropic and paratropic current-density fluxes yield both shielding and deshielding contributions, as the shielding or deshielding is determined by the direction of the current-density flux with respect to the studied nucleus instead of the tropicity. Becke partitioning of the magnetic shieldings shows that the magnetic shielding contributions mainly originate from the studied atom and its nearest neighbors, confirming the localized character of nuclear magnetic shieldings.

Accurate and cost-effective NMR chemical shift predictions for proteins using a molecules-in-molecules fragmentation-based method

We have developed an efficient protocol using our two-layer Molecules-in-Molecules (MIM2) fragmentation-based quantum chemical method for the prediction of NMR chemical shifts of large biomolecules. To investigate the performance of our fragmentation approach and demonstrate its applicability, MIM-NMR calculations are first calibrated on a test set of six proteins. The MIM2-NMR method yields a mean absolute deviation (MAD) from unfragmented full molecule calculations of 0.01 ppm for 1H and 0.06 ppm for 13C chemical shifts. Thus, the errors from fragmentation are only about 3% of our target accuracy of ∼0.3 ppm for 1H and 2-3 ppm for 13C chemical shifts. To compare with experimental chemical shifts, a standard protocol is first derived using two smaller proteins 2LHY (176 atoms) and 2LI1 (146 atoms) for obtaining an appropriate protein structure for NMR chemical shift calculations. The effect of the solvent environment on the calculated NMR chemical shifts is incorporated through implicit, explicit, or explicit-implicit solvation models. The expensive first solvation shell calculations are replaced by a micro-solvation model in which only the immediate interaction between the protein and the explicit solvation environment is considered. A single explicit water molecule for each amine and amide proton is found to be sufficient to yield accurate results for 1H chemical shifts. The 1H and 13C NMR chemical shifts calculated using our protocol give excellent agreement with experiments for two larger proteins, 2MC5 (the helical part with 265 atoms) and 3UMK (33 residue slice with 547 atoms). Overall, our target accuracy of ∼0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for the larger proteins. The proposed MIM-NMR method is accurate and computationally cost-effective and should be applicable to study a wide range of large proteins.

Calculating nuclear magnetic resonance chemical shifts in solvated systems

The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.

Influence of Solvent Representation on Nuclear Shielding Calculations of Protonation States of Small Biological Molecules

In this study, we assess the influence of solvation on the accuracy and reliability of isotropic nuclear magnetic shielding calculations for amino acids in comparison to experimental data. We focus particularly on the performance of solvation methods for different protonation states, as biological molecules occur almost exclusively in aqueous solution and are subject to protonation with pH. We identify significant shortcomings of current implicit solvent models and present a hybrid solvation approach that improves agreement with experimental data by taking into account the presence of direct interactions between amino acid protonation state and water molecules.

Crystal packing and crystallization tendency from the melt of 2-((2-ethylphenyl)amino)nicotinic acid

2-((2-ethylphenyl)amino)nicotinic acid (2EPNA) was synthesized and its crystal structure was determined. It was observed that alkylation of the phenyl ring with ethyl group disrupts the planar conformation of the molecule by steric repulsion, resulting in formation of an acid-pyridine heterosynthon (instead of acid-acid homosynthon) in the crystal. Crystallization tendency from the melt state of the polymorph was studied by differential scanning calorimetry (DSC). It was revealed that this compound could form a very stable amorphous phase on melt quenching and not crystallize even on re-heating. The formation of acid-pyridine hydrogen bonding in the amorphous state is believed to be responsible for its good glass forming ability.

Effect of positional isomerism and vanadium substitution on 51V magic angle spinning NMR Spectra Of Wells-Dawson polyoxotungstates

We examined the positional isomerism and vanadium substitution on the 51V magic angle spinning NMR spectra of potassium salts of vanadium-substituted polyoxotungstates of the Wells-Dawson series. NMR parameters of this class of catalytically active polyoxotungstates effect of are reported. Multiple species, indicative of differences in the local environment at the substitution sites, are observed in solid-state NMR spectra of the di- and tri- substituted complexes in contrast to solution NMR spectra, where single average chemical shift was observed. The quadrupolar and chemical shift anisotropy parameters depend strongly on the position and the degree of the vanadium substitution into the oxoanion core establishing 51V SATRAS NMR spectroscopy as a sensitive probe of the local electronic environment in these catalytically active solids.

State of the art and perspectives in the application of quantum chemical prediction of 1 H and 13 C chemical shifts and scalar couplings for structural elucidation of organic compounds

Quantum mechanical prediction of chemical shifts and scalar couplings involving ¹ H and ¹³ C nuclei is now a popular tool in structural elucidation of organic compounds. Here, we summarize the current state of the art and present to the reader present limitations and problems, mostly related to treatment of conformational flexibility, as well as future perspectives and potential applications in the field.

(1)H NMR spectra of alcohols in hydrogen bonding solvents: DFT/GIAO calculations of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in hydrogen bonding solvents have been computed on the basis of density functional theory by applying the gauge-including atomic orbital method to geometry-optimized alcohol/solvent complexes. The OH proton shifts and hydrogen bond distances for methanol or ethanol complexed with pyridine depend very much on the functional employed and very little on the basis set, provided it is sufficiently large to give the correct quasi-linear hydrogen bond geometry. The CH proton shifts are insensitive to both the functional and the basis set. NMR shifts for all protons in several alcohol/pyridine complexes are calculated at the Perdew, Burke and Ernzerhof PBE0/cc-pVTZ//PBE0/6-311 + G(d,p) level in the gas phase. The results correlate with the shifts for the pyridine-complexed alcohols, determined by analysing data from the NMR titration of alcohols against pyridine. More pragmatically, computed shifts for a wider range of alcohols correlate with experimental shifts in neat pyridine. Shifts for alcohols in dimethylsulfoxide, based on the corresponding complexes in the gas phase, correlate well with the experimental values, but the overall root mean square difference is high (0.23 ppm), shifts for the OH, CHOH and other CH protons being systematically overestimated, by averages of 0.42, 0.21 and 0.06 ppm, respectively. If the computed shifts are corrected accordingly, a very good correlation is obtained with a gradient of 1.00 ± 0.01, an intercept of 0.00 ± 0.02 ppm and a root mean square difference of 0.09 ppm. This is a modest improvement on the result of applying the CHARGE programme to a slightly different set of alcohols. Some alcohol complexes with acetone and acetonitrile were investigated both in the gas phase and in a continuum of the relevant solvent.

Solvent Dependence of (14)N Nuclear Magnetic Resonance Chemical Shielding Constants as a Test of the Accuracy of the Computed Polarization of Solute Electron Densities by the Solvent

Although continuum solvation models have now been shown to provide good quantitative accuracy for calculating free energies of solvation, questions remain about the accuracy of the perturbed solute electron densities and properties computed from them. Here we examine those questions by applying the SM8, SM8AD, SMD, and IEF-PCM continuum solvation models in combination with the M06-L density functional to compute the (14)N magnetic resonance nuclear shieldings of CH3CN, CH3NO2, CH3NCS, and CH3ONO2 in multiple solvents, and we analyze the dependence of the chemical shifts on solvent dielectric constant. We examine the dependence of the computed chemical shifts on the definition of the molecular cavity (both united-atom models and models based on superposed individual atomic spheres) and three kinds of treatments of the electrostatics, namely the generalized Born approximation with the Coulomb field approximation, the generalized Born model with asymmetric descreening, and models based on approximate numerical solution schemes for the nonhomogeneous Poisson equation. Our most systematic analyses are based on the computation of relative (14)N chemical shifts in a series of solvents, and we compare calculated shielding constants relative to those in CCl4 for various solvation models and density functionals. While differences in the overall results are found to be reasonably small for different solvation models and functionals, the SMx models SM8, and SM8AD, using the same cavity definitions (which for these models means the same atomic radii) as those employed for the calculation of free energies of solvation, exhibit the best agreement with experiment for every functional tested. This suggests that in addition to predicting accurate free energies of solvation, the SM8 and SM8AD generalized Born models also describe the solute polarization in a manner reasonably consistent with experimental (14)N nuclear magnetic resonance spectroscopy. Models based on the nonhomogeneous Poisson equation show slightly reduced accuracy. Scaling the intrinsic Coulomb radii to larger values (as has sometimes been suggested in the past) does not uniformly improve the results for any kind of solvent model; furthermore it uniformly degrades the results for generalized Born models. Use of a basis set that increases the outlying charge diminishes the accuracy of continuum models that solve the nonhomogeneous Poisson equation, which we ascribe to the inability of the numerical schemes for approximately solving the nonhomogeneous Poisson equation to fully account for the effects of electronic charge outside the solute cavity.

1H and 13C NMR scaling factors for the calculation of chemical shifts in commonly used solvents using density functional theory

Calculation of NMR chemical shifts and coupling constants using quantum mechanical calculations [density functional theory (DFT)], has become a very popular tool for the determination of conformation and the assignment of stereochemistry within a molecule. We present the scaling factors (linear regression parameters) from 10 DFT methods for 10 commonly used NMR solvents using the same set of reference compounds. The results were compared with the corresponding gas-phase calculations to assess the inclusion of the polarizable continuum model for solvent effects.

Investigation of solute-solvent interactions in phenol compounds: accurate ab initio calculations of solvent effects on 1H NMR chemical shifts

Accurate (1)H chemical shifts of the -OH groups of polyphenol compounds can be calculated, compared to experimental values, using a combination of DFT, polarizable continuum model (PCM) and discrete solute-solvent hydrogen bond interactions. The study focuses on three molecular solutes: phenol, 4-methylcatechol and the natural product genkwanin in DMSO, acetone, acetonitrile, and chloroform. Excellent linear correlation between experimental and computed chemical shifts (with the GIAO method at the DFT/B3LYP/6-311++G(2d,p) level) was obtained with minimization of the solvation complexes at the DFT/B3LYP/6-31+G(d) and DFT/B3LYP/6-311++G(d,p) level of theory with a correlation coefficient of 0.991. The use of the DFT/B3LYP/6-31+G(d) level of theory for minimization could provide an excellent means for the accurate prediction of the experimental OH chemical shift range of over 8 ppm due to: (i) strong intramolecular and solute-solvent intermolecular hydrogen bonds, (ii) flip-flop intramolecular hydrogen bonds, and (iii) conformational effects of substituents of genkwanin. The combined use of ab initio calculations and experimental (1)H chemical shifts of phenol -OH groups provides a method of primary interest in order to obtain detailed structural, conformation and electronic description of solute-solvent interactions at a molecular level.

Accurate prediction of 195Pt NMR chemical shifts for a series of Pt(ii) and Pt(iv) antitumor agents by a non-relativistic DFT computational protocol

Exhaustive benchmark DFT calculations reveal that the non-relativistic GIAO-PBE0/SARC-ZORA(Pt)∪6-31+G(d)(E) computational protocol predicts accurate ¹⁹⁵Pt NMR chemical shifts for a wide range of square planar Pt(ii) and octahedral Pt(iv) anticancer agents.

Quantum chemical calculations of amide-15N chemical shift anisotropy tensors for a membrane-bound cytochrome-b5

There is considerable interest in determining amide-(15)N chemical shift anisotropy (CSA) tensors from biomolecules and understanding their variation for structural and dynamics studies using solution and solid-state NMR spectroscopy and also by quantum chemical calculations. Due to the difficulties associated with the measurement of CSA tensors from membrane proteins, NMR-based structural studies heavily relied on the CSA tensors determined from model systems, typically single crystals of model peptides. In the present study, the principal components of backbone amide-(15)N CSA tensors have been determined using density functional theory for a 16.7 kDa membrane-bound paramagnetic heme containing protein, cytochrome-b(5) (cytb(5)). All the calculations were performed by taking residues within 5 Å distance from the backbone amide-(15)N nucleus of interest. The calculated amide-(15)N CSA spans agree less well with our solution NMR data determined for an effective internuclear distance r(N-H) = 1.023 Å and a constant angle β = 18° that the least shielded component (δ(11)) makes with the N-H bond. The variation of amide-(15)N CSA span obtained using quantum chemical calculations is found to be smaller than that obtained from solution NMR measurements, whereas the trends of the variations are found to be in close agreement. We believe that the results reported in this study will be useful in studying the structure and dynamics of membrane proteins and heme-containing proteins, and also membrane-bound protein-protein complexes such as cytochromes-b5-P450.

Interpretation of crystal effects on NMR chemical shift tensors: electron and shielding deformation densities

The relationship between the NMR observables and the supramolecular structure of any system is not straightforward. In this work we examine the influence of the crystal packing for three purine derivatives (hypoxanthine, theobromine, and 6-(2-methoxy)benzylaminopurine) on the principal components of the NMR chemical shift tensors (CSTs). We employ density functional calculations to obtain various molecular properties (the ground-state electron density, the magnitudes and orientations of the components of NMR chemical shift tensor, and the spatial distribution of the isotropic magnetic shielding) for the isolated molecules and for the molecules embedded in supramolecular clusters modeling the crystal environment and evaluate their differences. The concept has enabled us to rationalize the effect of the crystal packing on the NMR CSTs in terms of the redistribution of the ground-state electron density induced by intermolecular interactions in the solid state.

Enzyme active site interactions by Raman/FTIR, NMR, and ab initio calculations

Characterization of enzyme active site structure and interactions at high resolution is important for the understanding of the enzyme catalysis. Vibrational frequency and NMR chemical shift measurements of enzyme-bound ligands are often used for such purpose when X-ray structures are not available or when higher resolution active site structures are desired. This review is focused on how ab initio calculations may be integrated with vibrational and NMR chemical shift measurements to quantitatively determine high-resolution ligand structures (up to 0.001 Å for bond length and 0.01 Å for hydrogen bonding distance) and how interaction energies between bound ligand and its surroundings at the active site may be determined. Quantitative characterization of substrate ionic states, bond polarizations, tautomeric forms, conformational changes and its interactions with surroundings in enzyme complexes that mimic ground state or transition state can provide snapshots for visualizing the substrate structural evolution along enzyme-catalyzed reaction pathway. Our results have shown that the integration of spectroscopic studies with theoretical computation greatly enhances our ability to interpret experimental data and significantly increases the reliability of the theoretical analysis.

A computational investigation of carbon-doped beryllium monoxide nanotubes

To investigate the influence of C-doping on the electrostatic structure properties in the frame work of density functional theory (DFT), we considered beryllium monoxide nanotubes (BeONTs), consisting of 60 Be and 60 O atoms. Full geometry optimizations are performed for all structures, i.e., all atoms are allowed to relax. Afterwards, the chemical shielding (CS) tensors are calculated for Be-9, O-17 and C-13 nuclei in the C-doped forms and also pristine models of the (10, 0) zigzag and (5, 5) armchair BeONTs. Formation energies indicate that C-doping of Be atom (CBe form) could be more favorable than C-doping of O atom (CO form) in both zigzag and armchair BeONTs. Gap energies and dipole moments detected the effects of dopant in the (5, 5) armchair models; however, those parameters did not indicate any significant changes in the C-doped (10, 0) zigzag BeONT models. The results show that the CS values for the Be and O atoms-contributed to the Be-C bonds or those atoms close to the C-doped region-in the CO form of BeONTs (zigzag and armchair) are changed. The same values only for the O atoms-contributed to the O-C bonds- in the CBe form of BeONTs (zigzag and armchair) are changed.

Toward the quantum chemical calculation of nuclear magnetic resonance chemical shifts of proteins

Despite the many protein structures solved successfully by nuclear magnetic resonance (NMR) spectroscopy, quality control of NMR structures is still by far not as well established and standardized as in crystallography. Therefore, there is still the need for new, independent, and unbiased evaluation tools to identify problematic parts and in the best case also to give guidelines that how to fix them. We present here, quantum chemical calculations of NMR chemical shifts for many proteins based on our fragment-based quantum chemical method: the adjustable density matrix assembler (ADMA). These results show that (13)C chemical shifts of reasonable accuracy can be obtained that can already provide a powerful measure for the structure validation. (1)H and even more (15)N chemical shifts deviate more strongly from experiment due to the insufficient treatment of solvent effects and conformational averaging.

Aggregation of perfluoroctanoate salts studied by 19F NMR and DFT calculations: counterion complexation, poly(ethylene glycol) addition, and conformational effects

The aggregation of perfluoroctanoate salts in H(2)O is studied by (19)F NMR on solutions of LiPFO, NaPFO, and CsPFO, without and with the addition of two poly(ethylene glycol) (PEG) oligomers of molecular weight 1500 and 3400 Da, respectively, and with the addition of suitable crown ethers. The (19)F chemical shift (cs) trends are monitored, at 25 °C, in a concentration range including the critical micellar concentration (cmc) or, in the presence of PEG, the critical aggregation concentration (cac). The cac values in the samples with PEG are lower than the cmc values of the corresponding samples without PEG; moreover, the (19)F cs trends above the cac and above the polymer saturation concentration reveal and help to explain some peculiarities of the aggregation process of PEG on PFO micelles, which, in the first step, seems to occur while the surfactant concentration in water is still increasing. Also in LiPFO/H(2)O or NaPFO/H(2)O solutions containing 12-crown-4 or 15-crown-5 ethers, suitable to complex Li(+) or Na(+) ions, respectively, the cmc decreases. On the other hand, the micellization process in the presence of crown ethers does not show other peculiarities. The prevailing conformations of the PFO chain are discussed on the basis of quantum-mechanical calculations. The theoretical chemical shifts were computed at the DFT level of theory, taking into account the effects of the environment by means of the IEF-PCM method. The helical structure is the most stable one, but anti conformations are easily accessible, in both the aqueous and fluorinated environment. The comparison between computed and experimental chemical shifts indicates that anti conformations are more important in the micelles than in water and in CsPFO micelles than in LiPFO or NaPFO ones.

Intermolecular shielding contributions studied by modeling the (13)C chemical-shift tensors of organic single crystals with plane waves

In order to predict accurately the chemical shift of NMR-active nuclei in solid phase systems, magnetic shielding calculations must be capable of considering the complete lattice structure. Here we assess the accuracy of the density functional theory gauge-including projector augmented wave method, which uses pseudopotentials to approximate the nodal structure of the core electrons, to determine the magnetic properties of crystals by predicting the full chemical-shift tensors of all (13)C nuclides in 14 organic single crystals from which experimental tensors have previously been reported. Plane-wave methods use periodic boundary conditions to incorporate the lattice structure, providing a substantial improvement for modeling the chemical shifts in hydrogen-bonded systems. Principal tensor components can now be predicted to an accuracy that approaches the typical experimental uncertainty. Moreover, methods that include the full solid-phase structure enable geometry optimizations to be performed on the input structures prior to calculation of the shielding. Improvement after optimization is noted here even when neutron diffraction data are used for determining the initial structures. After geometry optimization, the isotropic shift can be predicted to within 1 ppm.