Importance of Triples Contributions to NMR Spin-Spin Coupling Constants Computed at the CC3 and CCSDT Levels

On the performance of HRPA(D) for NMR spin-spin coupling constants: Smaller molecules, aromatic and fluoroaromatic compounds

In this study, the performance of the doubles-corrected higher random-phase approximation [HRPA(D)] has been investigated in calculations of nuclear magnetic resonance spin-spin coupling constants (SSCCs) for 58 molecules with the experimental values used as the reference values. HRPA(D) is an approximation to the second-order polarization propagator approximation (SOPPA) and is, therefore, computationally less expensive than SOPPA. HRPA(D) performs comparable and sometimes even better than SOPPA, and therefore, when calculating SSCCs, it should be considered as an alternative to SOPPA. Furthermore, it was investigated whether a coupled-cluster singles, doubles and perturbative triples [CCSD(T)] or Møller-Plesset second order (MP2) geometry optimization was optimal for a SOPPA and a HRPA(D) SSCC calculation for eight smaller molecules. CCSD(T) is the optimal geometry optimization for the SOPPA calculation, and MP2 was optimal for HRPA(D) SSCC calculations.

Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants

Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.

Fluorine spin-spin coupling constants of pentafluorobenzene revisited at the ab initio correlated levels

All possible spin-spin coupling constants, ¹⁹ F-¹⁹ F, ¹⁹ F-¹³ C, and ¹⁹ F-¹ H, of pentafluorobenzene were calculated at five different levels of theory, HF, DFT, SOPPA (CCSD), CCSD, and the SOPPA (CCSD)-based composite scheme with taking into account solvent, vibrational, relativistic, and correlation corrections. Most corrections were next to negligible for the long-range couplings but quite essential for the one-bond carbon-fluorine coupling constants. Hartree-Fock calculations were found to be entirely unreliable, while DFT results were comparable in accuracy with the data obtained using the wave function-based methods.

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational ¹ H and ¹³ C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation

We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.

On the Unexpected Accuracy of the M06L Functional in the Calculation of 1JFC Spin-Spin Coupling Constants

One-bond spin-spin coupling constants (SSCCs) between F and C are computed with density functional theory (DFT). Surprisingly, M06L stands out for its striking accuracy, outperforming any other investigated functional, including PBE0, otherwise considered one of the most reliable for couplings involving F. Although the computation of nuclear magnetic resonance (NMR) parameters involving F is known to be a challenging task, even with a rather small basis set as pcJ-1, M06L provides results with a MAD = 11.7 Hz, whereas the average deviation gets as much as 5 times larger for PBE0 (MAD = 60.0 Hz). In the context of SSCCs on the order of 300 Hz, this is particularly remarkable. We find that the accuracy of M06L/pcJ-1 in predicting ¹J_FC constants does not stem from a well-suited exchange or correlation part of the functional. Instead, it is believed to arise from a fortuitous cancellation of errors, as revealed by investigating the convergence of the basis set. Our findings also indicate that ¹J_FC constants are highly dependent on the amount of exact exchange included in the expression of the functional, with large fractions being critically important to achieving satisfactory results. Studying the effects of the geometry on ¹J_FC, we find that optimizing the geometry at the level of theory used to calculate SSCCs generally improves the quality of the results, although the combination of a M06-2X/aug-cc-pVTZ geometry with M06L/pcJ-1 ¹J_FC constants best reproduces the experimental data for organofluorine compounds (with the exception of fluoroalkenes).

The aug-cc-pVTZ-J basis set for the p-block fourth-row elements Ga, Ge, As, Se, and Br

The aug-cc-pVTZ-J basis set family is extended to include the fourth-row p-block elements Ga, Ge, As, Se, and Br. We use the established approach outlined by Sauer and coworkers (J. Chem. Phys. 115, 1324 [2001], J. Chem. Phys. 133, 054308 [2010], J. Chem. Theory Comput. 7, 4070 [2011], and J. Chem. Theory Comput. 7, 4077 [2011]) where the completely uncontracted aug-cc-pVTZ basis set is saturated with tight s-, p-, d-, and f-functions to form the aug-cc-pVTZ-Juc basis set for the tested elements. The saturation is carried out on the simplest hydrides possible for the tested elements GaH, GeH₄ , AsH₃ , H₂ Se, and HBr until an improvement is less than 0.01% for all s-, p-, and d-functions added. f-Functions are added to an improvement less than or equal to 1.0% due to the computational expense these functions add. The saturated aug-cc-pVTZ-Juc (26s16p12d5f) is then recontracted using the molecular orbital coefficients from self-consistent field calculations on the simple hydrides to improve computational efficiency. During contraction of the basis set, we observe that the linear hydrogen bromide molecule has a slower convergence than the other tested molecules which sets a limit on the accuracy obtained. All calculations with the contracted aug-cc-pVTZ-J [17s10p7d5f] gives results that are within 1.0% of the uncontracted results at considerable computational savings.

Multi-configurational short-range density functional theory can describe spin-spin coupling constants of transition metal complexes

The multi-configurational short-range (sr) density functional theory has been extended to the calculation of indirect spin-spin coupling constants (SSCCs) for nuclear magnetic resonance spectroscopy. The performance of the new method is compared to Kohn-Sham density functional theory and the ab initio complete active space self-consistent field for a selected set of molecules with good reference values. Two density functionals have been considered, the local density approximation srLDA and srPBE from the GGA class of functionals. All srDFT calculations are of Hartree-Fock-type HF-srDFT or complete active space-type CAS-srDFT. In all cases, the calculated SSCC values are of the same quality for srLDA and srPBE functionals, suggesting that one should use the computationally cost-effective srLDA functionals in applications. For all the calculated SSCCs in organic compounds, the best choice is HF-srDFT; the more expensive CAS-srDFT does not provide better values for these single-reference molecules. Fluorine is a challenge; in particular, the FF, FC, and FO couplings have much higher statistical errors than the rest. For SSCCs involving fluorine and a metal atom CAS-srDFT with singlet, generalized Tamm-Dancoff approximation is needed to get good SSCC values although the reference ground state is not a multi-reference case. For VF₆ ^-1, all other considered models fail blatantly.

An efficient method for generating property-energy consistent basis sets. New pecJ-n (n = 1, 2) basis sets for high-quality calculations of indirect nuclear spin-spin coupling constants involving 1H, 13C, 15N, and 19F nuclei

This paper presents a new method of generating property-energy consistent (PEC) basis sets that can be applied to any arbitrary molecular property. The PEC method generates a basis set that is optimized for the molecular property under interest, providing the least possible total molecular energy. The main algorithm of the PEC approach involves Monte Carlo simulations to generate random exponents in the predetermined range. In this work, the PEC method is introduced in the example of generation of new pecJ-n (n = 1, 2) basis sets suited for high-quality correlated calculations of indirect nuclear spin-spin coupling constants involving the most popular NMR-active nuclei: 1H, 13C, 15N, and 19F.

Benchmarking Correlated Methods for Static and Dynamic Polarizabilities: The T145 Data Set Evaluated with RPA, RPA(D), HRPA, HRPA(D), SOPPA, SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD

Due to the importance of predicting static and dynamic polarizabilities, the performance of various correlated linear response methods including random phase approximation (RPA), RPA(D), higher-order random phase approximation (HRPA), HRPA(D), second-order polarization propagator approximation (SOPPA), SOPPA(CC2), SOPPA(CCSD), CC2, and CCSD has been evaluated against CCSD(T) (static case) and CCSD (dynamic cases) for the T145 set of 145 organic molecules. The benchmark reveals that the HRPA(D) method has the best performance for both static and dynamic polarizabilities apart from CCSD. RPA(D) ranks second for the dynamic cases and third for the static case. Using coupled-cluster amplitudes in SOPPA(CCSD) and SOPPA(CC2), the SOPPA results are significantly improved. The HRPA method has the largest deviations from the reference values for both cases. In general, according to the performance and computational cost of the methods, the HRPA(D) and RPA(D) methods are proposed for calculations of static and dynamic polarizabilities of this and similar sets of molecules.

Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties

A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.

Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives

This review is written amid a marked progress in the calculation of NMR parameters of carbohydrates substantiated by a vast amount of experimental data coming from several laboratories worldwide. By no means are we trying to cover in the present compilation a huge amount of all available data. The main idea of the present review was only to outline general trends and perspectives in this dynamically developing area on the background of a marked progress in theoretical and computational NMR. Presented material is arranged in three basic sections: (1)-a brief theoretical introduction; (2)-applications and perspectives in computational NMR of monosaccharides; and (3)-calculation of NMR chemical shifts and spin-spin coupling constants of di- and polysaccharides.

Karplus relationships of the 2 JHNα and 3 JΗΝβ couplings in organic azides

Density functional theory (DFT) and second-order polarization propagator approximation (SOPPA) computations in model organic azides revealed a Karplus-like dependence not only of the vicinal ³ J_H-C-Nα-Nβ coupling but also of the geminal ² J_H-C-Nα one, with the H-C-N_α N_β dihedral angle. Karplus equations were derived from the DFT computations on the isopropylazide model system. In light of these stablished relationships, natural abundance ¹ H-¹⁵ N couplings obtained for the azide group of the zidovudine antiviral helped to probe its conformation around the C-N_α bond as being of the synclinal type.

Correlated ab initio calculations of one-bond 31 P77 Se and 31 P125 Te spin-spin coupling constants in a series of PSe and PTe systems accounting for relativistic effects (part 2)

Synthetic chalcogen-phosphorus chemistry permanently makes new challenges to computational Nuclear Magnetic Resonance (NMR) spectroscopy, which has proven to be a powerful tool of structural analysis of chalcogen-phosphorus compounds. This paper reports on the calculations of one-bond ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te NMR spin-spin coupling constants (SSCCs) in the series of phosphine selenides and tellurides. The applicability of the combined computational approach to the one-bond ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te SSCCs, incorporating the composite nonrelativistic scheme, built of high-accuracy correlated SOPPA (CC2) and Coupled Cluster Single and Double (CCSD) methods and the Density Functional Theory (DFT) relativistic corrections (four-component level), was examined against the experiment and another scheme based on the four-component relativistic DFT method. A special J-oriented basis set (acv3z-J) for selenium and tellurium atoms, developed previously by the authors, was used throughout the NMR calculations in this work at the first time. The proposed computational methodologies (combined and 'pure') provided a reasonable accuracy for ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te SSCCs against experimental data, characterizing by the mean absolute percentage errors of about 4% and 1%, and 12% and 8% for selenium and tellurium species, respectively. The present study reports typical relativistic corrections to ⁷⁷ Se³¹ P and ¹²⁵ Te³¹ P SSCCs, calculated within the four-component DFT formalism for a broad series of tertiary phosphine selenides and tellurides with different substituents at phosphorus.

NMR parameters of FNNF as a test for coupled-cluster methods: CCSDT shielding and CC3 spin-spin coupling

NMR shielding and spin-spin coupling constants of cis and trans isomers of FNNF have been determined to near-quantitative accuracy from ab initio calculations. The FNNF system, containing multiple N-F bonds and fluorine atoms, provides a severe test of computational methods. Coupled-cluster methods were used with large basis sets and complete basis set (CBS) extrapolations of the equilibrium geometry results, with vibrational and relativistic corrections. Shielding constants were calculated with basis sets as large as aug-cc-pCV7Z, together with coupled-cluster expansions up to CCSDT, at the all-electron CCSD(T)/aug-cc-pCVQZ optimized geometries. Spin-spin coupling constants have been determined with specialized versions of the correlation consistent basis sets ccJ-pVXZ, further augmented with diffuse functions. All-electron coupled-cluster methods up to CC3 were applied in these calculations. The results of this work highlight the application of state-of-the-art theoretical techniques, and provide the most accurate NMR properties of FNNF to date, which can serve to guide and supplement NMR experimentation.

Calculated coupling constants 1 J(X-Y) and 1 K(X-Y), and fundamental relationships among the reduced coupling constants for molecules Hm X-YHn , with X, Y ═ 1 H, 7 Li, 9 Be, 11 B, 13 C, 15 N, 17 O, 19 F, 31 P, 33 S, and 35 Cl

Equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) calculations have been performed to determine coupling constants ¹ J(X-Y) for 65 molecules H_m X-YH_n , with X,Y ═ ¹ H, ⁷ Li, ⁹ Be, ¹¹ B, ¹³ C, ¹⁵ N, ¹⁷ O, ¹⁹ F, ³¹ P, ³³ S, and ³⁵ Cl. The computed ¹ J(X-Y) values are in good agreement with available experimental data. The reduced coupling constants ¹ K(X-Y) have been derived from ¹ J(X-Y) by removing the dependence on the magnetogyric ratios of X and Y. Patterns are found for the reduced coupling constants on a ¹ K(X-Y) surface that are related to the positions of X and Y in the periodic table.

How large a 3 JHCCH coupling can be? Prediction of giant vicinal couplings in quinonoid systems

Second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) B3LYP computations revealed the presence of giant ³ J_HCCH couplings, up to more than 60 Hz, in conjugated quinonoid systems, very far above the typical limit of ~15 Hz. Strong hyperconjugative interactions of the C-H σ orbitals with the quinonoid 8-e^- π system, which allows effective propagation of the spin polarization, seem to be responsible for those extraordinarily large couplings. Computation on several model systems showed the additive nature of the contributions from the available coupling paths to the total coupling.

Computational 1 H NMR: Part 1. Theoretical background

This is the first one of the three closely interrelated reviews to be published in Magnetic Resonance in Chemistry dealing with accordingly theoretical background, chemical applications, and biochemical studies of and by means of computational ¹ H NMR. Presented in the first part of the review is a general outline of the modern theoretical methods and accuracy factors of computational ¹ H NMR involving locally dense basis set schemes, solvent effects, vibrational corrections, and relativistic effects performed at the density functional theory and/or nonempirical levels. This review is dedicated to Prof. Stephan Sauer in view of his invaluable contribution to the field of computational nuclear magnetic resonance.

RPA(D) and HRPA(D): Two new models for calculations of NMR indirect nuclear spin-spin coupling constants

In this article, the RPA(D) and HRPA(D) models for the calculation of linear response functions are presented. The performance of the new RPA(D) and HRPA(D) models is compared to the performance of the established RPA, HRPA, and SOPPA models in calculations of indirect nuclear spin-spin coupling constants using the CCSD model as a reference. The doubles correction offers a significant improvement on both the RPA and HRPA models; however, the improvement is more dramatic in the case of the RPA model. For all coupling types investigated in this study, the results obtained using the HRPA(D) model are comparable in accuracy to those given by the SOPPA model, while requiring between 30% and 90% of the calculation time needed for SOPPA. The RPA(D) model, while of slightly lower accuracy compared to the CCSD model than HRPA(D), offered calculation times of only approximately 25% of those required for SOPPA for all the investigated molecules. © 2018 Wiley Periodicals, Inc.

Computational protocols for calculating 13C NMR chemical shifts

The most recent results dealing with the computation of ¹³C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.

Exploring Through-Space Spin-Spin Couplings for Quantum Information Processing: Facing the Challenge of Coherence Time and Control Quantum States

Nuclear magnetic resonance (NMR) is a powerful tool for studying quantum information processing (QIP). Recently quantum technologies have been proposed to overcome the challenges in large-scale NMR QIP. Furthermore, computational chemistry can promote its improvement. Nuclear spins-¹/₂ are natural qubits and have been used in most NMR quantum computation experiments. However, molecules that enable many qubits NMR QIP implementations should meet some requirements regarding their spectroscopic properties. Exceptionally large through-space (TS) P-P spin-spin coupling constants (SSCC or J) observed in 1,8-diphosphanaphthalenes (PPN) and in naphtho[1,8- cd]-1,2-dithiole phenylphosphines (NTP) were proposed and investigated to provide more accurate control within large-scale NMR QIP. Spectroscopic properties of PPN and NTP derivatives were explored by theoretical strategies using locally dense basis sets (LDBS). ³¹P chemical shifts (δ) calculated at the B3LYP/aug-cc-pVTZ-J level and TS P-P SSCCs at the PBE1PBE/pcJ-2 (LDBS-1) level are very close to the experimental data for the PPN molecule. Differently, for the NTP dimer, PBE1PBE/pcJ-2 (LDBS-2) predicts more accurate ³¹P δ, whereas PBE1PBE/Def2-TZVP (LDBS-1) forecasts more accurate TS P-P SSCCs. From our results, PPN_o-F, PPN_o-ethyl, and PPN_o-NH₂ were the best candidates for NMR QIP, in which the large TS SSCCS could face the need of long-time quantum gates implementations. Therefore, it could overcome natural limitations concerning the development of large-scale NMR.

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.

Value of NMR Parameters and DFT Calculations for Quantum Information Processing Utilizing Phosphorus Heterocycles

Quantum computing is the field of science that uses quantum-mechanical phenomena, such as superposition and entanglement, to perform operations on data. The fundamental information unit used in quantum computing is the quantum bit or qubit. It is well-known that quantum computers could theoretically be able to solve problems much more quickly than any classical computers. Currently, the first and still the most successful implementations of quantum information processing (QIP) have been based on nuclear spins in liquids. However, molecules that enable many qubits NMR QIP implementations should meet some conditions: have large chemical shifts and be appropriately dispersed for qubit addressability, appreciable spin-spin coupling between any pair of spins, and a long relaxation time. In this line, benzyldene-2,3-dihydro-1H-[1,3]diphosphole (BDF) derivatives have been theoretically tested for maximizing large chemical shifts, spin-spin coupling, and minimizing the hyperfine coupling constant. Thus, the structures were optimized at the B3LYP/6-311G(d,p) level and showed a significant similarity with the experimental geometrical parameters. The NMR spectroscopic parameters (δ and J) were calculated with six different DFT functionals. The τ-HCTH/6-31G(2d) level is in better agreement with the experimental data of ³¹P and ¹³C chemical shifts, while PCM-B3LYP/cc-pVDZ level shows a decrease on deviation between calculated and experimental values for P-P and P-C SSCC. The surface response technique was employed to rationalize how the hyperfine constant varies with the chemical shifts and coupling constants values. From our findings, BDF-NO₂ was the best candidate for NMR quantum computations (NMR-QC) among the studied series.