Computing the1H NMR Spectrum of a Bulk Ionic Liquid from Snapshots of Car–Parrinello Molecular Dynamics Simulations

Structural Features of the [C4mim][Cl] Ionic Liquid and Its Mixtures with Water: Insight from a 1H NMR Experimental and QM/MD Study

The 1H NMR chemical shift of water exhibits non-monotonic dependence on the composition of an aqueous mixture of 1-butyl-3-methylimidazolium chloride, [C4mim][Cl], ionic liquid (IL). A clear minimum is observed for the 1H NMR chemical shift at a molar fraction of the IL of 0.34. To scrutinize the molecular mechanism behind this phenomenon, extensive classical molecular dynamics simulations of [C4mim][Cl] IL and its mixtures with water were carried out. A combined quantum mechanics/molecular mechanics approach based on the density functional theory was applied to predict the NMR chemical shifts. The proliferation of strongly hydrogen-bonded complexes between chloride anions and water molecules is found to be the reason behind the increasing 1H NMR chemical shift of water when its molar fraction in the mixture is low and decreasing. The model shows that the chemical shift of water molecules that are trapped in the IL matrix without direct hydrogen bonding to the anions is considerably smaller than the 1H NMR chemical shift predicted for the neat water. The structural features of neat IL and its mixtures with water have also been analyzed in relation to their NMR properties. The 1H NMR spectrum of neat [C4mim][Cl] was predicted and found to be in very reasonable agreement with the experimental data. Finally, the experimentally observed strong dependence of the chemical shift of the proton at position 2 in the imidazolium ring on the composition of the mixture was rationalized.

Phosphonate Chelators for Medicinal Metal Ions

A family of phosphonate-bearing chelators was synthesized to study their potential in metal-based (radio)pharmaceuticals. Three ligands (H₆phospa, H₆dipedpa, H₆eppy; structures illustrated in manuscript) were fully characterized, including X-ray crystallographic structures of H₆phospa and H₆dipedpa. NMR spectroscopy techniques were used to confirm the complexation of each ligand with selected trivalent metal ions. These methods were particularly useful in discerning structural information for Sc³⁺ and La³⁺ complexes. Solution studies were conducted to evaluate the complex stability of 15 metal complexes. As a general trend, H₆phospa was noted to form the most stable complexes, and H₆eppy associated with the least stable complexes. Moreover, In³⁺ complexes were determined to be the most stable, and complexes with La³⁺ were the least stable, across all metals. Density functional theory (DFT) was employed to calculate structures of H₆phospa and H₆dipedpa complexes with La³⁺ and Sc³⁺. A comparison of experimental ¹H NMR spectra with calculated ¹H NMR spectra using DFT-optimized structures was used as a method of structure validation. It was noted that theoretical NMR spectra were very sensitive to a number of variables, such as ligand configuration, protonation state, and the number/orientation of explicit water molecules. In general, the inclusion of an explicit second shell of water molecules qualitatively improved the agreement between theoretical and experimental NMR spectra versus a polarizable continuum solvent model alone. Formation constants were also calculated from DFT results using potential-energy optimized structures. Strong dependence of molecular free energies on explicit water molecule number, water molecule configuration, and protonation state was observed, highlighting the need for dynamic data in accurate first-principles calculations of metal-ligand stability constants.

Computational NMR Study of Ion Pairing of 1-Decyl-3-methyl-imidazolium Chloride in Molecular Solvents

The 1H NMR spectra of 10-5 mole fraction solutions of 1-decyl-3-methyl-imidazolium chloride ionic liquid in water, acetonitrile, and dichloromethane have been measured. The chemical shift of the proton at position 2 in the imidazolium ring of 1-decyl-3-methyl-imidazolium (H2) is rather different for all three samples, reflecting the shifting equilibrium between the contact pairs and free fully solvated ions. Classical molecular dynamics simulations of the 1-decyl-3-methyl-imidazolium chloride contact ion pair as well as of free ions in water, acetonitrile, and dichloromethane have been conducted, and the quantum mechanics/molecular mechanics methods have been applied to predict NMR chemical shifts for the H2 proton. The chemical shift of the H2 proton was found to be primarily modulated by hydrogen bonding with the chloride anion, while the influence of the solvents-though differing in polarity and capabilities for hydrogen bonding-is less important. By comparing experimental and computational results, we deduce that complete disruption of the ionic liquid into free ions takes place in an aqueous solution. Around 23% of contact ion pairs were found to persist in acetonitrile. Ion-pair breaking into free ions was predicted not to occur in dichloromethane.

Computational NMR Spectroscopy of Ionic Liquids: [C4C1im]Cl/Water Mixtures

In this work, I have analyzed the structure of binary mixtures of 1-butyl-3-methylimidazolium chloride ionic liquid, [C₄C₁im]Cl, and water, using computational NMR spectroscopy. The structure of the complex fluid phase, where the ionic and hydrophobic nature of ionic liquids is further complicated by the addition of water, is first generated by classical Molecular Dynamics (MD) and then validated by calculating the NMR properties with DFT at the ONIOM(B3LYP/cc-pVTZ//B3LYP/3-21G) on clusters extracted during the MD trajectories. Three ionic liquid/water mixtures have been considered with the [C₄C₁im]Cl mole fraction of 1.00, 0.50, and 0.01, that is the pure ionic liquid [C₄C₁im]Cl, the equimolar [C₄C₁im]Cl/water mixture, and a diluted solution of [C₄C₁im]Cl in water. A good agreement is obtained with published experimental data that, at the same time, validates the structural features obtained from the MD and the force field used, and provides an example of the power of NMR spectroscopy applied to complex fluid phases.

Bioionic Liquid Conjugation as Universal Approach To Engineer Hemostatic Bioadhesives

Adhesion to wet and dynamic surfaces is vital for many biomedical applications. However, the development of effective tissue adhesives has been challenged by the required combination of properties, which includes mechanical similarity to the native tissue, high adhesion to wet surfaces, hemostatic properties, biodegradability, high biocompatibility, and ease of use. In this study, we report a novel bioinspired design with bioionic liquid (BIL) conjugated polymers to engineer multifunctional highly sticky, biodegradable, biocompatible, and hemostatic adhesives. Choline-based BIL is a structural precursor of the phospholipid bilayer in the cell membrane. We show that the conjugation of choline molecules to naturally derived polymers (i.e., gelatin) and synthetic polymers (i.e., polyethylene glycol) significantly increases their adhesive strength and hemostatic properties. Synthetic or natural polymers and BILs were mixed at room temperature and cross-linked via visible light photopolymerization to make hydrogels with tunable mechanical, physical, adhesive, and hemostatic properties. The hydrogel adhesive exhibits a close to 50% decrease in the total blood volume loss in tail cut and liver laceration rat animal models compared to the control. This technology platform for adhesives is expected to have further reaching application vistas from tissue repair to wound dressings and the attachment of flexible electronics.

Simulating structure and dynamics in small droplets of 1-ethyl-3-methylimidazolium acetate

To investigate the structure and dynamics of small ionic liquid droplets in gas phase, we performed a DFT-based ab initio molecular dynamics study of several 1-ethyl-3-methylimidazolium acetate clusters in vacuum as well as a bulk phase simulation. We introduce an unbiased criterion for average droplet diameter and density. By extrapolation of the droplet densities, we predict the experimental bulk phase density with a deviation of only a few percent. The hydrogen bond geometry between cations and anions is very similar in droplets and bulk, but the hydrogen bond dynamics is significantly slower in the droplets, becoming slower with increasing system size, with hydrogen bond lifetimes up to 2000 ps. From a normal mode analysis of the trajectories, we identify the modes of the ring proton C-H stretching, which are strongly affected by hydrogen bonding. From analyzing these, we find that the hydrogen bond becomes weaker with increasing system size. The cations possess an increased concentration inside the clusters, whereas the anions show an excess concentration on the outside. Almost all anions point towards the droplet center with their carboxylic groups. Ring stacking is found to be a very important structural motif in the droplets (as in the bulk), but side chain interactions are only of minor importance. By using Voronoi tessellation, we define the exposed droplet surface and find that it consists mainly of hydrogen atoms from the cation's and anion's methyl and ethyl groups. Polar atoms are rarely found on the surface, such that the droplets appear completely hydrophobic on the outside.

Aggregation behavior of dihexadecylviologen bistriflimide ionic liquid crystal in different solvents: influence of polarity and concentration

Solutions of 1,1'-dihexadecyl-4,4'-bipyridinium di[bis(trifluoromethanesulfonyl)imide] salt, also known as dihexadecylviologen bistriflimide, in deuterated acetonitrile (ACN), dichloromethane (DCM) and chloroform (CDCl3), respectively, were investigated by the combination of 1H and DOSY NMR spectroscopy, DFT calculations and MD simulation to understand the influence of solvent polarity and solute concentration (10-5-10-1 M) on its aggregation behavior. We found that the polar solvent acetonitrile (ACN) does not favor ion aggregation and cluster formation. In the whole range of concentrations investigated, the system appears to be dominated by neutral ion pairs composed of one cation and two anions, possibly in fast equilibrium (on the NMR time scale) with small or slightly larger aggregates. The diffusion coefficient of the cationic species is only weakly affected by concentration. In contrast, the low-polar solvents of chloroform (CDCl3) and dichloromethane (DCM) strongly favor cluster formation above a certain concentration and the viologen diffusion coefficient in CDCl3 is much smaller and more strongly dependent on concentration than that in ACN. The information obtained from the MD simulations suggests that the aggregates have a structure similar to the isotropic liquid phase of the viologen-based ionic liquids and ionic liquid crystals. The lifetimes of such large clusters appear to be relatively long, beyond the time scale of tens of nanoseconds. Moreover, the results from the aromatic proton NMR resonances provide some insights on the dielectric constants inside the viologen aggregates.

Prediction of 1H NMR chemical shifts for clusters of imidazolium-based ionic liquids

Nuclear magnetic resonance (NMR) has been widely used to elucidate the bulk structure of ionic liquids. In this work, we calculated ¹H NMR chemical shifts of 1-ethyl-3-methylimidazolium (C₂mim⁺) ionic liquids combined with various anions such as chloride (Cl), tetrafluoroborate (BF₄), hexafluorophosphate (PF₆), acetate (OAc), trifluoroacetate (TFA), and dicyanamide (DCA). The previously established level of theory, HF/6-311G+(3df,2p), was used for the accurate prediction of NMR chemical shifts both in gas phase and in solvents with varying dielectric constant such as CHCl₃ and ethanol. The following factors affecting the predicted proton chemical shifts were considered. Firstly, ionic clusters consisting of 2, 8 and 16 ion pairs were optimized to model interionic interactions present in the bulk of ionic liquids. In larger clusters the distribution of the calculated chemical shifts of individual protons in the C₂mim⁺ cation was examined with respect to the position of the cation in the cluster. We further confirmed that electronic properties of ionic liquids such as magnetic shielding had local nature, thus allowing us to accurately predict proton NMR chemical shifts of ionic liquids from relatively small-sized clusters. Secondly, solvent effects in single ion pairs as well as larger ionic clusters were accounted through a Conductor-like Polarisable Continuum Model (CPCM). Solvent effects generated through a dielectric constant of either chloroform or ethanol were found to be important in single ion pairs due to improved description of interionic distances. With increasing cluster size the difference between gas-phase and CPCM optimized structures became minimal, thus resulting in similar values for calculated ¹H NMR chemical shifts. We also established that the model size that produced the best results for imidazolium ionic liquids strongly depended on the anion type. Strongly coordinating anions such as chloride and acetate require calculations of clusters consisting of at least 8 ion pairs, whereas weakly coordinating anions produce excellent accuracy for single ion pairs optimized in the presence of solvent. The polarity of the solvent was found to play a minor role.

Ion pairing in 1-butyl-3-methylpyridinium halide ionic liquids studied using NMR and DFT calculations

Experimental and calculated ¹H, ¹³C and ¹⁵N NMR data of bulk 1-butyl-3-methylimidazolium halides ionic liquids provide key insights on their ion pairing.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Bio ionic liquids and water mixtures: a structural study

In this study we have explored, by means of ab initio molecular dynamics, a subset of three different water/cho⁺–phe⁻ mixtures.

Ionic Liquid Crystals: Versatile Materials

This Review covers the recent developments (2005-2015) in the design, synthesis, characterization, and application of thermotropic ionic liquid crystals. It was designed to give a comprehensive overview of the "state-of-the-art" in the field. The discussion is focused on low molar mass and dendrimeric thermotropic ionic mesogens, as well as selected metal-containing compounds (metallomesogens), but some references to polymeric and/or lyotropic ionic liquid crystals and particularly to ionic liquids will also be provided. Although zwitterionic and mesoionic mesogens are also treated to some extent, emphasis will be directed toward liquid-crystalline materials consisting of organic cations and organic/inorganic anions that are not covalently bound but interact via electrostatic and other noncovalent interactions.

Temperature Dependence of NMR Parameters Calculated from Path Integral Molecular Dynamics Simulations

The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with conventional DFT molecular dynamics or with 1D scans through the potential energy surface. Predicted NMR observables for the model systems were in excellent agreement with experimental data.

Quantitative prediction of physical properties of imidazolium based room temperature ionic liquids through determination of condensed phase site charges: a refined force field

Quantitative prediction of physical properties of room temperature ionic liquids through nonpolarizable force field based molecular dynamics simulations is a challenging task. The challenge lies in the fact that mean ion charges in the condensed phase can be less than unity due to polarization and charge transfer effects whose magnitude cannot be fully captured through quantum chemical calculations conducted in the gas phase. The present work employed the density-derived electrostatic and chemical (DDEC/c3) charge partitioning method to calculate site charges of ions using electronic charge densities obtained from periodic density functional theory (DFT) calculations of their crystalline phases. The total ion charges obtained thus range between -0.6e for chloride and -0.8e for the PF6 ion. The mean value of the ion charges obtained from DFT calculations of an ionic liquid closely matches that obtained from the corresponding crystal thus confirming the suitability of using crystal site charges in simulations of liquids. These partial charges were deployed within the well-established force field developed by Lopes et al., and consequently, parameters of its nonbonded and torsional interactions were refined to ensure that they reproduced quantum potential energy scans for ion pairs in the gas phase. The refined force field was employed in simulations of seven ionic liquids with six different anions. Nearly quantitative agreement with experimental measurements was obtained for the density, surface tension, enthalpy of vaporization, and ion diffusion coefficients.

Physical, structural, and dehydrogenation properties of ammonia borane in ionic liquids

Hydrogen desorption profiles of AB–ILs with H₂ yield.

A simple AIMD approach to derive atomic charges for condensed phase simulation of ionic liquids

The atomic charges for two ionic liquids (ILs), 1-n-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]) and 1-ethyl-3-methylimidazolium hexafluorophosphate ([EMIM][PF6]), were derived from periodic crystal phase calculations with density functional theory (DFT) and plane wave basis sets (denoted as "AIMD-c charge"). For both ILs, the total charge was found to be ±0.8 e for the cation and anion, respectively, due to the charge transfer between ions and polarization caused by the environment. These atomic charges were used in a force field developed within the general Amber force field framework. Using this force field, static, dynamic, and thermodynamic properties were computed for the two ILs using molecular dynamics simulation. The results were compared against results obtained using the same Amber force field but four different sets of partial charges, denoted as full charge, scaled charge, AIMD-l charge, and AIMD-b charge, respectively. The full charge was derived from quantum chemistry calculation of isolated ions in a vacuum and resulted in a total charge of unity on each ion. The scaled charge was obtained by uniformly scaling the full charge by 0.8. AIMD-l and AIMD-b charges were derived from liquid phase ab initio molecular dynamics simulations. The scaled charges have the same total charge on the ions as the AIMD-c charge but different distributions. It was found that simulation results not only depend on the total charge of each ion, but they are also sensitive to the charge distribution within an ion, especially for dynamic and thermodynamic properties. Overall, for the two ILs under study, the AIMD-c charge was found to predict experimental results better than the other four sets of charges, indicating that fitting charges from crystal phase DFT calculations, instead of extensive sampling of the liquid phase configurations, is a simple and reliable way to derive atomic charges for condensed phase ionic liquid simulations.

Vibrational averaging of the chemical shift in crystalline α-glycine

Averaging of the chemical shift over the molecular motion improves the simulated data and provides additional information about the temperature dependence and system dynamics. However, crystal modeling is difficult due to the limited precision of the plane-wave density functional theory (DFT) methods and approximate vibrational schemes. On the glycine example, we investigate how the averaging can be achieved within the periodic boundary conditions at the DFT level. The nuclear motion is modeled with the vibrational configuration interaction, with other simplified quantum anharmonic schemes, and the classical Born-Oppenheimer molecular dynamics (BOMD). The results confirm a large vibrational contribution to the isotropic shielding values. Both the first and second derivatives of the shielding were found important for the quantum averaging. The first derivatives influence the shielding mostly due to the anharmonic character of the CH and NH stretching modes, whereas second derivatives produce most vibrational corrections associated with the lower-frequency vibrational modes. Temperature excitations of the lowest-frequency vibrational states and the expansion of the crystal cell both determine the temperature dependence of nuclear magnetic resonance parameters. The vibrational quantum approach as well as classical BOMD schemes provided temperature dependencies of the chemical shifts that are consistent with the previous experimental data.

Dynamical effects in ab initio NMR calculations: classical force fields fitted to quantum forces

NMR chemical shifts for an L-alanine molecular crystal are calculated using ab initio plane wave density functional theory. Dynamical effects including anharmonicity may be included by averaging chemical shifts over an ensemble of structural configurations generated using molecular dynamics (MD). The time scales required mean that ab initio MD is prohibitively expensive. Yet the sensitivity of chemical shifts to structural details requires that the methodologies for performing MD and calculating NMR shifts be consistent. This work resolves these previously competing requirements by fitting classical force fields to reproduce ab initio forces. This methodology is first validated by reproducing the averaged chemical shifts found using ab initio molecular dynamics. Study of a supercell of L-alanine demonstrates that finite size effects can be significant when accounting for dynamics.

Synthesis and applications of ionic liquids derived from natural sugars

Aiming to develop environmentally compatible chemical syntheses, the replacement of traditional organic solvents with ionic liquids (ILs) has attracted considerable attention. ILs are special molten salts with melting points below 100 degrees C that are typically constituted of organic cations (imidazolium, pyridinium, sulfonium, phosphonium, etc.) and inorganic anions. Due to their ionic nature, they are endowed with high chemical and thermal stability, good solvent properties, and non-measurable vapor pressure. Although the recovery of unaltered ILs and recycling partly compensate their rather high cost, it is important to develop new synthetic approaches to less expensive and environmentally sustainable ILs based on renewable raw materials. In fact, most of these alternative solvents are still prepared starting from fossil feedstocks. Until now, only a limited number of ILs have been prepared from renewable sources. Surprisingly, the most available and inexpensive raw material, i.e., carbohydrates, has been hardly exploited in the synthesis of ILs. In 2003 imidazolium-based ILs were prepared from o-fructose and used as solvents in Mizoroki-Heck and Diels-Alder reactions. Later on, the first chiral ILs derived from sugars were prepared from methyl D-glucopyranoside. In the same year, a family of new chiral ILs, obtained from commercial isosorbide (dianhydro-D-glucitol), was described. A closely related approach was followed by other researchers to synthesize mono- and bis-ammonium ILs from isomannide (dianhydro-D-mannitol). Finally, a few ILs bearing a pentofuranose unit as the chiral moiety were prepared using sugar phosphates as glycosyl donors and 1-methylimidazole as the acceptor.

Calculation of NMR chemical shifts in organic solids: accounting for motional effects

NMR chemical shifts were calculated from first principles for well defined crystalline organic solids. These density functional theory calculations were carried out within the plane-wave pseudopotential framework, in which truly extended systems are implicitly considered. The influence of motional effects was assessed by averaging over vibrational modes or over snapshots taken from ab initio molecular dynamics simulations. It is observed that the zero-point correction to chemical shifts can be significant, and that thermal effects are particularly noticeable for shielding anisotropies and for a temperature-dependent chemical shift. This study provides insight into the development of highly accurate first principles calculations of chemical shifts in solids, highlighting the role of motional effects on well defined systems.

Predicting the NMR spectra of nucleotides by DFT calculations: cyclic uridine monophosphate

We present an experimental and quantum chemical NMR study of the mononucleotide cyclic uridine monophosphate in water. Spectral parameters ((1)H and (13)C chemical shifts and (1)H--(1)H, (13)C--(1)H, (31)P--(13)C and (31)P--(1)H spin-spin coupling constants) have been carefully obtained experimentally and calculated using DFT methods including the solvent effect and the conformational flexibility of the solute. This study confirms that the (1)H and (13)C spectra of polar, flexible molecules in aqueous solution can be predicted with a high level of accuracy, comparable to that obtained for less complex systems.

Refined potential model for atomistic simulations of ionic liquid [bmim][PF6]

Refined parameters of an atomistic interaction potential model for the room temperature ionic liquid 1-n-butyl,3-methylimidazolium hexafluorophosphate are presented. Classical molecular dynamics simulations have been carried out to validate this fully flexible all-atom model. It predicts the density of the liquid at different temperatures between 300 and 500 K within 1.4% of the experimental value. Intermolecular radial distribution functions and the spatial distribution functions obtained from the new model are in close agreement with ab initio simulations. The calculated diffusion coefficients of ions and the surface tension of the liquid agree well with experiment.

Prediction of the 1H and 13C NMR Spectra of α-d-Glucose in Water by DFT Methods and MD Simulations

We have applied computational protocols based on DFT and molecular dynamics simulations to the prediction of the alkyl 1H and 13C chemical shifts of alpha-d-glucose in water. Computed data have been compared with accurate experimental chemical shifts obtained in our laboratory. 13C chemical shifts do not show a marked solvent effect. In contrast, the results for 1H chemical shifts provided by structures optimized in the gas phase are only fair and point out that it is necessary to take into account both the flexibility of the glucose structure and the strong effect exerted by solvent water thereupon. Thus, molecular dynamics simulations were carried out to model both the internal geometry as well as the influence of solvent molecules on the conformational distribution of the solute. Snapshots from the simulation were used as input to DFT NMR calculations with varying degrees of sophistication. The most important factor that affects the accuracy of computed 1H chemical shifts is the solute geometry; the effect of the solvent on the shielding constants can be reasonably accounted for by self-consistent reaction field models without the need of explicitly including solvent molecules in the NMR property calculation.