An e ﬃ cient DFT method of predicting the one, two-and three-bond indirect spin – spin coupling constants involving a ﬂ uorine nucleus in ﬂ uoroalkanes

The values of the indirect nuclear spin – spin coupling constants for a series of aliphatic ﬂ uorocompounds have been calculated using DFT-based methods and compared with the experimental values of these parameters. The set of the molecular objects contained four ﬂ uoromethanes, ﬁ ve ﬂ uoroethanes, two ﬂ uorocyclopropanes, and eleven ﬂ uorocompounds containing either ﬁ ve-membered or six-membered rings. The e ﬀ ectiveness of three hybrid functionals, B3LYP, PBE0 and BHandH and three basis sets, 6-311++G(2d,p) (s), 6-311++G(3df,3pd) (m) and aug-pcJ-3-2006 (l) has been checked. In order to compare the results concerning various types of coupling constants and obtained by various methods, a prediction-quality criterion has been proposed. It has been found that only the BHandH functional ensures calculating the proper values of one-bond ﬂ uorine – carbon and two-bond ﬂ uorine – ﬂ uorine coupling constants. For this functional application of the s basis, the smallest of the bases tested, has already yielded acceptable good results. The DFT BHandH/s PCM method has also provided the proper values of n J (F,H) ( n ¼ 1, 2, 3) and n J (F,C) ( n ¼ 2, 3) parameters. On the other hand, the analysis of the limited number of the results concerning 1 J (C,H) coupling constants has pointed out that in this case this method is less e ﬀ ective than the DFT PBE0/l PCM method.


Introduction
Fluoroorganic chemistry has been developing for a long time.Numerous synthetic methods aimed at massive or selective uorination of various organic compounds have been elaborated.The uorine atom can replace virtually any hydrogen in a molecule and a huge number of uoroorganic compounds have been synthesized. 1,2The uorine substituent is strongly electronegative and, simultaneously, p-electron donating, as it possesses lone electron pairs which can efficiently conjugate to molecular p-electron systems.Introduction of such a substituent into a molecule usually remarkably modies its physicochemical and biological properties.Substances obtained this way play an important role in chemistry as synthons or solvents and in industry as materials of unique features, such as various uoropolymers. 3Moreover, many uoro derivatives are biologically active substances and are used in medicine as important drugs or anesthetics. 4Some uorocompounds are highly toxic and, unfortunately, have been used as chemical weapons (sarin, soman) or their components.At the same time, uorocompounds are very interesting objects for NMR investigations owing to 100% natural abundance of 19 F isotope, the nucleus of which possesses 1/2 spin and high magnetic moment.][7][8][9][10][11][12] Thus, the chemistry of uorinecontaining compounds is an important branch of chemistry and, simultaneously, NMR of such compounds is an interesting and important part of NMR spectroscopy.
It is thus obvious that it would be desirable to have theoretical methods of predicting NMR parameters, whose effectiveness would be well-checked for uorocompounds.4][15][16][17][18][19] Unfortunately, those methods are inefficient, due to practical reasons, even for the medium size molecules.For example, some difficulties already arose during application of such a method to CHF 3 molecule. 13It is true that computer techniques are continuously developing, but on the other hand, it is clear that high level theoretical methods are not a real way out.A good practical solution would be elaboration of a suitable DFT-based method.Such methods are presently standard in calculations of NMR parameters for organic compounds composed of C, H, N, Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warszawa, Poland.E-mail: agryff@ch.pw.edu.pl† Electronic supplementary information (ESI) available: Calculated and experimental values of selected NMR parameters for uoromethanes (Tables S1-S4), 1,1,2-triuoro-2-methylcyclopropane (Table S5) and for 7-uoronorbornane derivatives (Table S6).1][32][33] The problem of predicting 19 F chemical shis has recently been tested in detail on the large body of experimental data. 33,34On the other hand, the rst test calculations of the indirect spin-spin coupling constants in such compounds yielded denitely poor values of 1 J(F,C) parameter. 35It has also been noticed that the same problem concerns two-bond indirect spin-spin coupling constants between geminal uorines. 15Below we show some more examples illustrating such problems.Let us remind, however, that already 20 years ago the hope was expressed that aer construction of suitable functionals the DFT should work properly and that large divergences between theory and experiment would be eliminated also in the case of NMR parameters of uorine compounds. 35Indeed, in recent years it has been shown that values of one-bond uorine-carbon spin-spin coupling constants can be predicted quite well when using the BHandH functional, 36,37 already known for a long time.To the best of our knowledge, however, most hitherto examples showing the effectiveness of this method concern molecules in which uorine is bonded to sp 2 hybridized carbons. 11,12,18,38,39lso some interesting theoretical papers concerning calculations of various uorine-uorine coupling constants, including 2 J(F,F) constant, have appeared in recent years. 18,40Results reported in this work show that also for uoroalkanes and uorocycloalkanes good values of 1 J(F,C) spin-spin coupling constants are obtained when using BHandH functional.This approach appears to be effective also in the case of some other uorine-X (X ¼ H, C, F) coupling constants, including the difficult 2 J(F,F) case.

Calculations
Looking for an efficient and practicable theoretical method of calculating the values of nuclear indirect spin-spin coupling constants involving uorine nuclei for aliphatic uorocompounds, we have performed several test DFT calculations for the series of uoromethanes and uoroethanes.On the basis of the obtained results some tentative recommendations concerning the choice of theoretical methods suitable for such calculations could be formulated.The effectiveness of the selected methods has been then veried by applying them in calculations of NMR parameters for several larger objects.Taking into account some literature hints and also our earlier results we have decided to test the effectiveness of three hybrid functionals: B3LYP, 41 PBE0 (BPE1BPE in Gaussian 03) 42 and BHandH. 36,37The two rst functionals are very popular in NMR calculations 20,24,25,43,44 and the third was shown to yield the best values of one-bond uorine-carbon spin-spin coupling constants for uoroaromatics. 11,12,18,39In our calculations we have used three basis sets: 6-311++G(2d,p), 6-311++G(3df,3pd) and aug-pcJ-3-2006 abbreviated from here on: s, m and l, respectively.The rst two basis sets are of moderate size and more or less standard for NMR calculations, 20,24,25 whereas the third basis has been dedicated especially to calculations of the spin-spin coupling constants. 45It would be the most comfortable situation if a universal and relatively small basis set appears effective for NMR calculations, although some Jdedicated basis can also be considered.The problem with good bases of this type is that they are usually much larger than s basis, sufficient in most magnetic shielding calculations.Nevertheless, if the quality of the prediction were substantially better thanks to such larger basis (or perhaps appropriate locally-dense bases), the method would indeed be useful.That is why we performed some calculations with l basis, too.
It is a common knowledge that the values of NMR parameters obtained in calculation are highly dependent on the molecular geometry introduced in the input data.Most frequently, the calculations is performed for the theoretical minimum-energy molecular structures and the molecular geometry optimization is the rst step of the whole calculation. 20,24,25In the primary tests performed in this work the same functional and basis set have been used during both calculation steps (molecular geometry optimization and NMR parameters computation).For uoromethanes the calculations have been performed using all the combinations of three mentioned functionals and three basis sets.
7][48][49][50][51][52][53][54][55][56][57][58][59] In most of these sources only the absolute values of the indirect spin-spin coupling constants have been determined.Throughout this work the signs of these constants have been assumed to be the same as those of calculated values.6][57] NMR spectra for other compounds of interest were measured in various solvents (from CCl 4 to methanol-d 4 ).In such a situation, in order to compensate, at least partially, for medium effects, most of our calculations have been performed using polarizable continuum model of the solvent (PCM) 60 with the set of atomic radii collected in the Gaussian soware under the keyword "radii ¼ UAKS". 36ll quantum chemical calculations were performed using GAUSSIAN 03 soware, 61 which was also the source of functionals as well as s and m basis sets. 36,37The basis set l was taken from ref. 45.

An unexpected success
For a long time a simple HF molecule have been considered to be a troublesome case.In a fundamental work in which Tozer et al. 44 tested performance of various exchange-correlation functionals in predicting values of indirect nuclear spin-spin coupling constants, the authors decided to "omit the challenging HF molecule from the analysis".The reason was that as opposite to the other molecules tested, in the case of HF, the calculated 1 J(F,H) values were not satisfying, independently of the DFT method used.Indeed, the calculated value was at best 458.0 Hz, while the experimental one was 538.0 Hz. 62 Prompted by the encouraging results in predicting 1 J(F,C) in aromatic uorocompounds by DFT methods when using BHandH functional we have optimized the HF molecular geometry using BHandH/s method and then calculated 1 J(F,H) coupling constant for this molecule using the same functional and s and l bases.The results have been 469.2Hz and 543.3 Hz, respectively.Such an excellent result is probably an effect of error cancellation.Even if it is so, it will be worth knowing how frequently such error cancellations occur when calculating n J(F,X) parameters.

The simplest objects
Looking for an effective DFT method of calculating nuclear spin-spin coupling constants for aliphatic uorocompounds we have rst performed several test calculations of NMR parameters for the series of uoromethanes and four uoroethanes (CH 2 FCH 3 , CHF 2 CH 3 , CF 3 CH 3 and CF 3 CH 2 F).These objects were selected for two obvious reasons: the simplicity of their molecular structures and the accessibility of the precise experimental NMR data for them.
Fluoromethanes are the simplest uoroalkanes and possess rigid structures.Let us note, too, that the molecules of the chosen uoroethanes contain either methyl or triuoromethyl group.Owing to the symmetry of these groups, the three staggered conformers of these compounds have identical equilibrium molecular geometry and so are equally populated at any temperature.Moreover, at temperatures normally used in liquid state NMR measurements the internal rotation about C-C bonds in these molecules is rapid and averages out the magnetic environments of three protons (or uorines) as well as appropriate NMR parameters.This feature simplies the structures of NMR spectra and makes the comparison of the experimentally obtained and calculated values of NMR parameters much easier.
6][57] For tetrauoromethane 31 and selected uoroethanes 8,31,47,48 only the solution data are available.Fortunately, for molecules whose structures (including conformation) are solvent independent, the spin-spin coupling constants are also poorly solvent dependent.Nevertheless, as it was mentioned, in order to compensate for these small medium effects the polarizable continuum solvent model (PCM) 60 has been used in our calculations.Still, when comparing such theoretical data with experimental ones, the inherent limitations of the PCM method have to be kept in mind.
The values of the isotropic indirect spin-spin coupling constants calculated for uoromethanes are collected in Table 1.These values have been obtained with the aid of 9 theoretical methods in which the same functional/basis pair has been used at both calculation steps.The table also contains the values of these parameters determined experimentally as well as some high-level ab initio theoretical data.Since it is believed that the gas phase molecular geometries are usually very well reproduced already when using the s basis, the smallest one we have used, we additionally performed some test NMR calculations using such geometries independently of the basis used in the second step of computation.Indeed, the differences between the results obtained this way and by the more time consuming methods mentioned above were usually small, although not always negligible.
Similar data for the selected uoroethanes are collected in Table 2.The table contains the data concerning only these J(F,X) [or 1 J(C,H)] coupling constants the experimental values of which are available in the cited literature sources.Some additional calculated and experimental values of NMR parameters for View Article Online uoromethanes, uoroethanes discussed in this work can be found in the ESI (Tables S1-S4 †).

Prediction-quality criterion
Before we formulate a nal recommendation concerning theoretical calculations of NMR parameters, being of interest in this study, it seems desirable to assume a criterion allowing the quality of a theoretical prediction to be expressed in a more or less objective way.This seemingly trivial question is not such, actually, when the substantial amount of data obtained by various methods and results concerning various types of coupling constants are to be compared.It is obvious that, say, 5 Hz difference between theoretical and experimental value denotes denitely poor prediction in the case of a vicinal coupling constant and an excellent prediction of one-bond 19 F-13 C coupling constant.Again, the relative error is not an adequate criterion for small coupling constants.Moreover, it has to be taken into account that the available experimental values are of very different precision.The values of uorinecarbon coupling constants, for example, happen to be read out directly from the spectrum and reported as spectral line separations (without line-shape analysis). 47,50,53The error of such data can be of the order of the linewidth or higher.Taking all these arguments into account we have assumed for coupling constants expressed in Hz the prediction-quality indicator as: The good quality prediction will denote from here on C # 1, and poor prediction C $ 3.
Eventually, it is obvious that the size of the basis set determines the computer time needed to complete the calculation.Frequently, just the size of the basis set required for getting reliable results decides whether the proposed method is useful for molecules being of interest to chemists, or not.For instance, calculation of NMR parameters for triuoromethane (without molecular geometry optimization) demanded 4 minutes, 15 minutes and almost 15 hours of the nominal computer time, when using B3LYP functional and the three basis sets mentioned above.It is worth keeping in mind these numbers when comparing the effectiveness of particular calculation methods.

Valuation of the results for uoromethanes and uoroethanes
Application of the prediction-quality indicator (eqn (1)) to the data of Tables 1 and 2 immediately shows that in most cases the values of spin-spin coupling constants involving uorine nuclei can be successfully predicted using the BHandH/s method (the last entries in Tables 1-3).It is noteworthy that for CH 3 F, CH 2 F 2 and CHF 3 these DFT results are as good as or better than those obtained by high level ab initio calculations (Table 1). 13Especially in the case of 1 J(F,C) coupling constants the prevalence of BHandH/s calculation method over other methods tested is well visible (Table 3).For this method, only in the case of uoromethane, a remarkable 23.5 Hz difference between calculated and experimental results is observed.Nevertheless, this deviation is still smaller than any one obtained by methods using other tested functionals, and still represents only 14% of the   a In this case the data for 1,1-diuoro-2-chlorocyclopropane (Table 5) have also been included.
experimental parameter value (Q ¼ 2.7).The predictions of all other 1 J(F,C) values by the BHandH/s method are good, whereas methods using two other functionals have yielded remarkably poorer results.We will show that similar situation occurs for 2 J(F,F) coupling constants, although at the moment we have had at our disposal only two ab initio results, 13 which could be confronted with the results of our calculations.The results concerning other types of J(F,X) constants in uoromethanes and uoroethanes are also encouraging (Table 3), however, some types of coupling constants are poorly represented in the data set of Tables 1 and 2 and so conclusions concerning effectiveness of the BHandH/s method in those cases will be veried below.

1,2-Diuoroethane
It has been mentioned that non-rigidity of molecules of uoroethanes possessing methyl or triuoromethyl group simplies a comparison of the results of NMR measurements and theoretical calculations.The molecule of 1,2-diuoroethane, however, represents quite a different case.The antiperiplanar (trans) minimum-energy conformer of this compound differs from two enantiomeric synclinal (gauche) conformers.The populations of trans and gauche conformers are different and are temperature and solvent dependent.The molecule of 1,2-diuoroethane has been investigated many times in the past both experimentally and theoretically. 8,63-654][65] We have concluded that it could be interesting whether the values of spin-spin coupling constants calculated by BHandH/s method conrm this nding.Table 4 contains the results of our calculations performed for both types of conformers of this molecule.More precisely, the values calculated for antiperiplanar conformer and mean values of those calculated for two synclinal conformers are given in the table .First of all, one can notice that three of four vicinal spinspin coupling constants are strongly conformation dependent.
Comparing them with the appropriate experimental data clearly shows that the values of coupling constants calculated for gauche conformations are close to the experimental ones.Our results thus remain in full agreement with the earlier ndings concerning equilibrium conformation of this molecule.
It is perhaps noteworthy that the achieved agreement between calculated and experimental values of the indirect spinspin coupling constants for all uoroethanes is actually better than expected.As usual, all calculations have been done for equilibrium geometries of minimum-energy conformations.It is obvious, however, that the structure of all uoroethanes is not rigid, the barriers to internal rotations are low and even in the ground vibrational state molecules librate around equilibrium geometry.Moreover, at room temperature probably the higher vibrational states of these modes are populated.There is no guarantee that averaging of the spin-spin coupling constants by librations yields the same values as those for the equilibrium geometry.Neglect of these internal motions during NMR parameters calculation is probably an oversimplication. 66A much more sophisticated treatment of this problem, such as that developed by Schaefer, 66 would probably be more adequate.

Further objects
In order to test further the effectiveness of the selected theoretical method we have included into the analyzed data set the spin-spin coupling constants for two rather non-standard uorinated alkanes, namely, 1,1-diuoro-2-chlorocyclopropane and 1,1,2-triuoro-2-methylcyclopropane.A large number of spin-spin coupling constants for these compounds have been reported by Brey. 50,53We have found out that also in this case the values of most of the reported coupling constants have been reproduced well when the BHandH/s method is used (Tables 5  and S5 †).
The selected theoretical methods have also been applied to calculate NMR parameters for axial and equatorial conformers of uorocyclohexane investigated at low temperatures and ve uorinated norbornane derivatives.These seven objects have been interesting from our point of view due to their rigid structures.In this case, the reliability of our results could be veried by confronting them with the experimental data reported by Abraham et al. 47,59 The results of our calculations of the indirect 19 F- 13 C spin-spin coupling constants for these molecules and the reported experimental values of these parameters are collected in Tables 6 and S6.† Again, the agreement between calculated and experimental values, in general, is good.
Our collection of the experimental coupling constant data has contained until now only three two-bond coupling constants involving geminal uorines, those in cyclopropanes and 2,2-diuoronorbornane.In order to get more data for experiment/theory comparison we have performed calculations for a few other molecules possessing -CF 2fragment: 5,5-diuoro-1,3-dioxane, whose spectrum was recorded at very low temperature, 54 1,1-diuoro-3-methylcyclohexane which is expected to occur in solutions mostly as a conformer with the methyl substituent occupying the equatorial position 47 2,2-diuoro-cis-3-methylcyclopentanol 51 and 3,3-diuoro-5-methylcyclopentene. 52Three of these compounds probably occur in solution in measurement conditions as equilibrium conformer mixtures.Calculations show, however, that the 2 J(F,F) parameter is poorly conformation-dependent.This observation is in accord with the nding of Tormena et al. 67 who have shown that in CF 2 group the geminal F-F coupling constant depends mainly on the angle between C-F bonds involved.

Verication of the selected method
Finally, all the data collected during this work have been subjected to the prediction-effectiveness test using the proposed criterion (eqn (1)).The results of this analysis are summarized in Table 7 and graphically presented in Fig. 1 and S1.† They conrm that there are at least two cases in which theoretical predictions of the values of spin-spin coupling constants involving uorine nucleus by DFT methods commonly used for calculating NMR parameters yield results biased by unacceptably large errors.These troublesome parameters are 1 J(F,C) and 2 J(F,F). 15,17,35The rst entry of Table 7 shows that application of the DFT BHandH/s method does solve this problem in the case of one-bond 19 F- 13 C coupling constants involving the aliphaticcarbon-bound uorines (Q av ¼ 0.9).Let us remind that in recent years it has been shown that this method worked well also for uorines bonded to aromatic carbons. 11,12,18,38,39Our calculations additionally show, which is even more novel, that this DFT method solves a similar problem occurring in the case of a twobond coupling involving geminal uorines.We have found out that, without exception, the recommended DFT BHandH/s method reproduces experimental results well (Q av ¼ 0.6), much better than the method which uses the B3LYP functional (Q av ¼ 6.2).Since, overall, quite a lot of experimental and theoretical data have been compared we may conclude that the effectiveness of the DFT BHandH/s method in predicting 1 J(F,C) and 2 J(F,F) values is well-veried.Moreover, it has been noticed that in every case the differences between theoretical parameter À3.9 À3.2 À3.7 À1.9 3 J(F-1,H-3) 10.7 11.6 12.2 11.0 3 J(F-2,H-1) 12.0 13.5 14.2 13.4 3 J(F-1,H-1) 1.1 2.9 2.5 4.6 3 J(F-2,H-2) Table 6 Comparison of the experimentally determined fluorine-carbon spin-spin coupling constants (from ref. 47) for two conformers of fluorocyclohexane, 1,1-difluoro-3-methylcyclohexane and 2,2-difluoronorbornane with the values of these parameters calculated using BHandH/s (PCM) method.Superscripts: eqequatorial, axaxial, enendo and exexo inform about the stereochemical position occupied by a given fluorine substituent values gained for these parameters when using methods involving B3LYP and BHandH functionals come almost exclusively from different calculated contributions of the Fermicontact coupling mechanism to the total constant (Table 8).This nding remains in full agreement with the results of the work by García de la Vega and San Fabián 40 who have recently treated that problem in more detail.The values of the collected two-bond uorine-hydrogen coupling constants have been predicted very well by both B3LYP/s and BHandH/s methods.Their effectiveness expressed by Q av values has been 0.5 and 0.4, respectively.Surprisingly, the PBE0/l method has been much less effective in this case.
Next, as it concerns two-bond 19 F- 13 C couplings, we have found out that, without exception, the recommended BHandH/ s method has reproduced experimental results much better than B3LYP/s method.On the other hand, the much more costly PBE0/l method provides somewhat better values of these constants, at least for uoroethanes and uorocyclopropanes, although this improvement is practically meaningless (Tables 2, 5 and S5 †).One can notice, however, that for uorocyclopropanes the calculated values of 2 J(F,C) coupling constants are systematically lower than the experimental ones and that the divergences between experimental and calculated values of these constants are larger than in remaining cases.The averaged prediction-quality indicator (Q av ) for the recommended DFT method calculated separately for uorocyclopropanes and for other objects are 2.4 and 0.4, respectively (Table 7).Probably, this divergence is a consequence of the structural specicity of cyclopropanes, although the precision of experimental data may play some role in this case as well.
Indeed, a limited precision of some spin-spin coupling data would be a serious deciency from our point of view.At the same time, it seems that in ref. 47   individual 19 F- 13 C coupling constants from such spectra is impossible.
The of line overlapping could also affect the threebond coupling values, especially 3 J(F,H) ones reported for uorocyclopropanes, where protons are parts of complex spinsystems.In the case of 3 J(F,C) several coupling constant values have been of the order of the typical linewidths in the 13 C NMR spectra.Nevertheless, examination of the results in Tables 5, S5, 6, S6 and 7 clearly shows that all trends visible in experimental data are well reproduced by the theoretical values obtained by BHandH/s method.Out of 15 collected 3 J(F,H) values only two have been poorly predicted and the averaged prediction-quality indicator Q av was 1.5 (Table 7).Similarly, only one out of 15 3 J(F,C) values has been poorly predicted and the Q av indicator amounted to 1.4 (Table 7).
We have additionally noticed that the method suitable for predictions of J(F,X) (X ¼ H, C, F) coupling constants yields poor results in the case of 1 J(C,H) constants.It seems that in this case, the largest of three basis sets used in this study is required at the nal calculation step (Tables 1 and 3).When this basis is used, all three functionals tested work sufficiently well and the result is only poorly dependent on the basis used during geometry optimization.It is somewhat surprising that the results provided by B3LYP/s method (Q av ¼ 0.9) are only slightly worse than, say, PBE0/l (Q av ¼ 0.3).These observations, however, are a by-product of this work and probably need further verication.Quite a different approach to predicting 1 J(C,H) values has been proposed recently by Fabián et al. 68

Conclusions
A reliable theoretical prediction of the values of the indirect spin-spin coupling constants involving uorine nucleus is an important problem for the structural chemistry of uoroorganic compounds.In this work we have performed theoretical calculations of these NMR parameters for a set of uoromethanes, uoroethanes and several uorocycloalkanes using various DFT methods.In order to analyze the large amount of the obtained data we have adopted a prediction-quality criterion (eqn (1)) allowing the effectiveness of various theoretical methods to be compared.We have found out that the proposed criterion is a very convenient tool, which allows the effectiveness of various theoretical methods to be evaluated in a straightforward and, to some degree, objective manner.Using this criterion, we have concluded that the values of the isotropic indirect spin-spin coupling constants involving a 19 F nucleus are predicted well when using DFT BHandH/s PCM method during molecular geometry optimization and NMR calculations.To be more precise, the use of the BHandH functional is indispensable only during the calculations of 1 J(F,C) and 2 J(F,F) parameters.In the case of other coupling constants involving uorine, the results gained when using BHandH/s method are also better than those achieved with the methods using B3LYP or PBE0, although some of the results of two latter functionals are of comparable quality.Application of larger basis sets like m, l or other specialized sets as m 0 ¼ aug-cc-pVTZ-J, 69 suggested by the referee, seems to be unjustied, as it enhances computation time and simultaneously does not guarantee achievement of results.We believe that the recommended BHandH/s method yields the results which are sufficiently accurate and reliable, so that they can be useful in solving various chemical problems by combining NMR measurements and DFT calculations.At the same time, it is quite probable that the recommended method is not yet the optimal one.It seems, however, that further improvement of the experiment/theory agreement demands, rst of all, inclusion of the effects of molecular vibrations, molecular collisions and a more elaborate description of impact of other medium effects, into the theoretical model.Finally, two problems, which are to some degree connected with the subject of this study, are to be mentioned.We have noticed that, for uorinated carbons, independently of the functional used, a large basis set is necessary to obtain 1 J(C,H) values being close to the experimental ones.Secondly, some introductory tests seem to point out that the BHandH functional is less effective than B3LYP or PBE0 predicting carbon and uorine chemical shis, but this conclusion has yet to be tested.

Table 3
Comparison of the effectiveness of three DFT methods in predicting values of various types of spin-spin coupling constants for fluoromethanes and fluoroethanes.[Q avthe averaged value of the prediction-quality indicator (eqn (1)), nthe number of experimental values of a given type of coupling constant, n w , n m and n pnumber of values predicted well (Q # 1), moderately well (1 < Q < 3) and poorly (Q $ 3), respectively]

a
As reported in ref. 50, adopting signs being in accord with calculated values.b Molecular geometry optimization calculated by PBE0/m.

, 50 and 53 the values of the 19 F
-13 C coupling constants were read out directly from the line positions in the spectra, without lineshape analysis.When several magnetically nonequivalent uorine nuclei are present in a molecule the carbon signals in its13 C proton-decoupled NMR spectrum are split by several (in general different) carbon-uorine coupling constants.Such signals have a form of multiplet of several lines.Due to a limited resolution of standard13 C NMR spectra, for certain combinations of coupling constant values, some lines of such multiplets can overlap one another.For example, the C-6 signal of 2,2-diuoronorbornane (X part of ABX spin system) in a13 C NMR spectrum of standard resolution probably has a form of an apparent triplet.Similarly, C-7 signal is probably a doublet of broadened lines.Without performing lineshape analysis the determination of the accurate values of

Table 1
The experimental and theoretically predicted values of the indirect spin-spin coupling constants [Hz] for fluoromethanes.The molecular geometry optimization and NMR parameter calculations were performed using the same DFT method a Ref. 13. b Ref. 57. c Ref. 56. d Ref. 46. e The averaged value of the prediction-quality indicator of J(F,X) parameters for a given theoretical method.This journal is © The Royal Society of Chemistry 2016 RSC Adv., 2016, 6, 82783-82792 | 82785 Paper RSC Advances Open Access Article.Published on 26 August 2016.Downloaded on 10/22/2023 5:14:48 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Table 2
The experimental and theoretically predicted values of the selected indirect spin-spin coupling constants (in Hz) for a few fluoroethanes a In this case the molecular geometry was optimized using DFT PBE0/m PCM(solvent) method.bRef.47.c Ref.49.dRef.8. e The averaged value of the prediction-quality indicator for various J(F,X) parameters, including the data for uoromethanes (Table1

Table 4
Comparison of the theoretically predicted values of spinspin coupling constants for two conformers of 1,2-difluoroethane in CH 2 Cl 2 and C 6 H 12 with the experimental values of these parameters 3 J(F,H) 2 J(H,H) 3 J(H,H) 3 J(H,H) 3 J(F,F)a The values averaged by a virtual exchange between two synclinal (enantiomeric) optimum energy conformers.b Ref. 8.

Table 7
Comparison of the effectiveness of B3LYP/s and BHandH/s DFT methods in predicting values of various types of spin-spin coupling constants involving 19 F nucleus. [Q avthe averaged value of the prediction-quality indicator (eqn (1)), nthe number of experimental values of a given type of coupling constant, n w , n m and n pnumber of values predicted well (Q # 1), moderately well (1 < Q < 3) and poorly (Q $ 3), respectively] a Coupling constants for uoroalkanes other than uorocyclopropanes.b Coupling constants for uorocyclopropanes.

Table 8
Selected one-bond fluorine-carbon and two bond fluorine-fluorine spin-spin coupling constants as well as contributions of various mechanisms to these couplings calculated by DFT methods, compared with the measured values of spin-spin constants Compound Method Fermi contact Spin dipolar Paramagnetic SO Diamagnetic SO 1 J(F,C) DFT 1 J(F,C) exp Ref.