Structure, solvent, and relativistic effects on the NMR


Download Structure, solvent, and relativistic effects on the NMR


Preview text

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

PCCP
PAPER

View Article Online
View Journal | View Issue

Cite this: Phys. Chem. Chem. Phys., 2015, 17, 24944

Structure, solvent, and relativistic effects on the NMR chemical shifts in square-planar transitionmetal complexes: assessment of DFT approaches†
Jan V´ıcha,‡ab Jan Novotny´,‡a Michal Straka,ac Michal Repisky,d Kenneth Ruud,d Stanislav Komorovsky*d and Radek Marek*ae

Received 18th July 2015, Accepted 24th August 2015 DOI: 10.1039/c5cp04214c
www.rsc.org/pccp

The role of various factors (structure, solvent, and relativistic treatment) was evaluated for square-planar 4d and 5d transition-metal complexes. The DFT method for calculating the structures was calibrated using a cluster approach and compared to X-ray geometries, with the PBE0 functional (def2-TZVPP basis set) providing the best results, followed closely by the hybrid TPSSH and the MN12SX functionals. Calculations of the NMR chemical shifts using the two-component (2c, Zeroth-Order Regular Approximation as implemented in the ADF package) and four-component (4c, Dirac–Coulomb as implemented in the ReSpect code) relativistic approaches were performed to analyze and demonstrate the importance of solvent corrections (2c) as well as a proper treatment of relativistic effects (4c). The importance of increased exact-exchange admixture in the functional (here PBE0) for reproducing the experimental data using the current implementation of the 2c approach is partly rationalized as a compensation for the missing exchange–correlation response kernel. The kernel contribution was identified to be about 15–20% of the spin–orbit-induced NMR chemical shift, DdSO, which roughly corresponds to an increase in DdSO introduced by the artificially increased exact-exchange admixture in the functional. Finally, the role of individual effects (geometry, solvent, relativity) in the NMR chemical shift is discussed in selected complexes. Although a fully relativistic DFT approach is still awaiting the implementation of GIAOs for hybrid functionals and an implicit solvent model, it nevertheless provides reliable NMR chemical shift data at an affordable computational cost. It is expected to outperform the 2c approach, in particular for the calculation of NMR parameters in heavy-element compounds.

1. Introduction

a CEITEC – Central European Institute of Technology, Masaryk University, Kamenice 5/A4, CZ-62500 Brno, Czech Republic. E-mail: [email protected]
b Centre of Polymer Systems, University Institute, Tomas Bata University in Zlin, Trida T. Bati 5678, CZ-76001 Zlin, Czech Republic
c Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo na´m. 2, CZ-16610 Prague, Czech Republic d Centre for Theoretical and Computational Chemistry, Department of Chemistry, UiT – The Arctic University of Norway, N-9037 Tromsø, Norway. E-mail: [email protected] e Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, CZ-62500 Brno, Czech Republic † Electronic supplementary information (ESI) available. Fig. S1: role of dispersion (empirical D3 correction) in cluster-based geometry optimizations; Fig. S2: the RMSDs (pm) for the interatomic distances according to the type of central metal; Fig. S3: the RMSDs (pm) for the interatomic distances according to the type of the light spectator atom; Fig. S4: comparison of total RMSDs (in pm) for interatomic distances calculated in vacuo and in cluster; Fig. S5: the 4c mDKS NMR chemical shifts with and without empirical ‘‘CGO corrections’’. Table S1: the Dd (ppm), in vacuo; Table S2: the Dd (ppm), COSMO; Table S3: the DdXkeCrnel values (ppm). See DOI: 10.1039/c5cp04214c ‡ These two authors contributed equally.

Square-planar transition-metal complexes, molecules with ligands arranged in a square-like pattern around the central metal, are typically preferred by d8 elements such as platinum(II). Their unique chemistry is dominated by a trans effect that makes them well suited for various practical applications ranging from cancer treatment (cisplatin, oxaliplatin)1 to chemical catalysis,2 molecular switches and beyond.3 Generally, a detailed analysis of the structure and various properties is required to understand their biological effects in the human body as well as their role in catalytic processes. In addition to experimental methods (e.g., X-ray diffraction, NMR spectroscopy) theoretical methods are also frequently being employed to investigate their structure, chemistry and properties, often supplementing experimental NMR studies. The synergy between theoretical and experimental approaches allows us to combine the experimental precision with the predictive and explanatory power of theoretical methods, yielding new insight unachievable by the experimental or theoretical approaches alone.4 Relativistic effects

24944 | Phys. Chem. Chem. Phys., 2015, 17, 24944--24955

This journal is © the Owner Societies 2015

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Paper

View Article Online
PCCP

have a significant influence on the NMR chemical shifts of light ligand atoms (LA) in the vicinity of the heavy-atom center (HA). Particular attention in NMR calculations should be paid to the spin–orbit effects, part of so-called HALA effects,5 typically propagated from the heavy element to light spectator atoms via the spin–orbit/Fermi-contact (SO/FC) mechanism.5,6 The spin–orbit contribution to the NMR chemical shift, DdSO, represents a relativistic correction to the NMR chemical shift that can amount from ten to hundreds of ppm.7
Calculated relativistic effects in the vicinity of the heavy metal are generally very sensitive to (i) the character of metal–ligand bonding,8,9 requiring the use of correct and accurate structures, (ii) the inclusion of environmental effects (e.g., solvent),10 and (iii) reliable methods for treating the relativistic effects.11
(i) Density functional theory (DFT) is nowadays the most commonly used method for optimizing the structures and calculating various properties of molecules and molecular systems, mainly due to its favorable scaling with system size and the rather good accuracy that can be achieved. A large number of density functionals have been developed over the years for calculating the molecular structures and/or molecular properties. The performance of individual functionals differs dramatically for various properties and systems, and careful evaluation/calibration of the functional behavior must be performed before production DFT calculation can be used for interpreting any chemical problem.
Nice examples of such assessments are the calibration studies by Bu¨hl and coworkers,8,12 which focused on the optimization of structures of transition-metal complexes, with DFT structures optimized in vacuo and referenced to structural parameters obtained in the gas phase. However, the geometries obtained in vacuo (or in an implicit solvent) are frequently inappropriately referenced to the structures determined by X-ray diffraction, neglecting any crystal effects. This approximation can lead to substantial errors and even to incorrect conclusions.13,14
To achieve the best available precision in the geometry optimizations, either molecular clusters based on experimental X-ray structures or periodic-boundary calculations should be used for calibrating the methods for solid-state structure optimization. In this work we focus on the cluster approach to assess the performance of the DFT functionals in optimizing the molecular geometry.
Recently, we applied this approach to a very limited number of complexes (mainly octahedral),10,15 but a wider selection of compounds and new functionals are required in order to draw more general conclusions. Therefore, the cluster-based comparison of ten DFT functionals (with or without D3 dispersion corrections) for structure optimization is presented in this work for seven square-planar complexes with various ligands and central metal atoms (Pt, Pd, Au, Rh, Fig. 1). These complexes were selected primarily based on the availability of high-quality X-ray structures and the availability of complete NMR data and variability of central metals.
(ii) The solvent effects can be simulated using explicit solvent molecules or implicit solvent models. Whereas implicit solvent models in general are limited to accounting for

Fig. 1 Structures of selected square-planar complexes. The transitionmetal centers are shown in blue and the ligand spectator NMR atoms (LA) are shown in green.
electrostatic effects of the solvent environment (continuum) on the NMR parameters of the solute,16 an explicit solvent model accounts also for specific weak interactions (hydrogen bonding, stacking),17 which can significantly alter the NMR chemical shifts of atoms involved in these interactions. This applies in particular to the hydrogen atoms and the easily polarizable heavy element.18 The application of an explicit solvent model is, however, beyond the scope of everyday calculations due to its complexity. For a proper determination of the NMR chemical shifts, classical molecular dynamic (MD) or QM/MD simulations must be performed to determine the positions of the solvent molecules relative to the solute and to calculate NMR chemical shifts averaged over individual snapshots, which can be very time and resource consuming, although great advances are being made for NMR parameters of molecules containing light elements.19,20 Furthermore, computationally demanding QM/MD simulations are limited by the size of the model that can be evaluated in a reasonable period of time, whereas specialized force-field parameters for MD simulations must be developed for each structurally different molecule or non-standard solvent.21 Therefore, the explicit solvent model is used predominantly in detailed studies of individual systems, whereas the implicit solvent model is

This journal is © the Owner Societies 2015

Phys. Chem. Chem. Phys., 2015, 17, 24944--24955 | 24945

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

PCCP

View Article Online
Paper

adopted more generally, being usually sufficient for calculations of many properties, including NMR chemical shifts of HALAinfluenced light atoms as we demonstrate in this work.
(iii) The relativistic effects can be treated in several ways at various levels of theory. A majority of the current computational codes operate with an approximate treatment of relativity. One of the most widely used approaches is the Zeroth-Order Regular Approximation. The ZORA Hamiltonian can be defined as:22

H^ZORA ¼ T^ZORA þ V^ ¼ r Á p^ c2 r Á p^ þ V^

(1)

2mc2 À V^

The relatively simple insertion of an electronic mean-field or the full Kohn–Sham potential into Tˆ ZORA allows for the inclu-
sion of SO contributions, yielding the two-component (2c) spin–orbit ZORA (SO-ZORA) approach.23,24
The SO-ZORA approach provides a good description of
the valence-shell orbitals, often comparable with that of fourcomponent methods.23 Because the NMR chemical shifts of the light atoms (e.g., 1H, 13C, 15N) are significantly influenced by
the type of chemical bonding, determined by the valence
orbitals, the performance of the SO-ZORA approach for NMR chemical shift calculations of light atoms is good.4,7,25,26 However, it has been demonstrated10,11,15,27 that increasing the
admixture of exact exchange in the functional is required to
obtain correct NMR chemical shifts of atoms strongly affected
by HALA effects. Recently, the role of the missing self-consistent
first-order response of the DFT exchange–correlation (XC) potential
to the external magnetic-field perturbation, the so-called XC
response kernel ( fxc), has been highlighted in NMR chemical shift calculations using the SO-ZORA approach.28 The missing terms arising from fxc were found to be quite sizable for hydrogen and mercury NMR chemical shifts. However, to the
best of our knowledge, these effects have not been evaluated for light, non-hydrogen NMR atoms (13C, 15N).
For a more complete inclusion of relativistic effects, the Dirac Hamiltonian must be considered:29

Hˆ D = caÁpˆ + bmc2 + Vˆ

(2)

In the Dirac formalism, operators and wave functions are represented as four-component objects. This will ultimately lead to more computationally demanding methods compared

to approaches utilizing more approximate Hamiltonians. Until recently, such methods have been considered applicable only to relatively small molecular systems (up to 10–20 atoms), which greatly limited their general applicability to solve reallife chemical problems. However, modern four-component methods have now reached a level of maturity that make them very useful tools for modeling and understanding chemically interesting systems.30,31
Therefore, in addition to the two-component (2c) SO-ZORA approach, we employed four-component (4c) calculations in the Dirac–Coulomb framework. The fully relativistic NMR calculations (ReSpect program) use restricted magnetically balanced basis sets for the small-component wave function.32 To further improve basis-set convergence, gauge-including atomic orbitals (GIAO) are employed for the pure DFT functionals.33 In the case of the hybrid functionals, only the common gauge origin (CGO) methodology is implemented in the ReSpect program, and thus special attention must be paid to errors arising from basis-set incompleteness in these calculations. In contrast, the implementation of the state-of-the-art non-collinear DFT kernel for GGA (generalized gradient approximation) functionals is used.34
Here we evaluate the influence of structure, solvent, spin– orbit contribution, exact-exchange admixture, role of fxc, and CGO approximation for hybrid functionals on relativistic NMR chemical shift calculations (2c and 4c) of light atoms in the vicinity of a heavy element for a series of square-planar transitional-metal complexes.
2. Methods
2.1 Selection of complexes and experimental data
A set of square-planar complexes with high-quality X-ray structures, complete NMR chemical shift data, and a diversity of the central metal atom was selected. The experimental data for the selected complexes are summarized in Table 1. In addition to ‘‘regular’’ square-planar platinum and palladium complexes with aromatic and aliphatic ligands (Pt and Pd), also complexes with ideal square-planar geometries distorted by the shape of the ligands, namely Au1 and Rh (see Fig. 1) are investigated. In Au1, the 2-benzoylpyridine ligand is bent out of the plane by 301, whereas for Rh, where the square-planar shape is enforced

Table 1 Acronym, formula, R-factor, Cambridge structural database code (CSD code), NMR chemical shift (ppm) for selected atoms bonded to heavy metal (highlighted in Fig. 1), NMR solvent used in the study, and corresponding references for selected square-planar complexes

Acronym
Pt1
Pt2 Pd1 Pd2 Au1
Au2 Rh

Formula
trans(S,N)-[Pt(2-ppy*)(DMSO)Cl]
[Pt(oxalato)(1R,2R-cyclohexanediamine)] [Pd(2,20-biquinoline)Cl2] [Pd(4,4-di-tert-butylbipyridine)Cl2] [Au(2-benzoylpyridine)Cl2]
[Au(2-phenylpyridine)Cl3] [Rh(dipyrrolylphoshinoxylene) CO]

R-factor (%)
1.8
1.9 2.7 4.8 2.9
2.2 3.3

CSD code
JISPAD01
CUHKEV YASPAK MOYWIG01 PUKYAV
YIDMAA SOWDUE

a All 15N NMR chemical shifts are reported relative to liquid NH3. b This work.

NMR chemical shifta
15N 220.6 13C 140.2 15N 9.1b 15N 224.0 15N 211.1 15N 215.2 13C 136.7 15N 227.7 31P 119.8 13C 172.2

NMR solvent
DMSO-d6
DMSO-d6 CDCl3 DMSO-d6 DMSO-d6
CDCl3 THF-d8

Ref.
35 and 36
37 38 and 39 40 and 41 42 and 43
44 and 45 46

24946 | Phys. Chem. Chem. Phys., 2015, 17, 24944--24955

This journal is © the Owner Societies 2015

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Paper

View Article Online
PCCP

by the so-called PCP ‘‘pincer’’ ligand, the bond angles around the central metal are about 121 distorted from the ideal 901.
2.2 Geometry optimization – DFT functional and basis set
Molecular cluster models were prepared for each complex individually, with a single central molecule (color-coded in Fig. 2) completely surrounded by crystal-packing molecules (orangecoded in Fig. 2) in a way that all close-contact interactions between the central molecule and the packing ones were treated in the final cluster. As a result, the central molecule was completely surrounded in a crystal-like environment with a packing sphere up to 12 Å, depending on the size of the central and packing molecules. Subsequently, the coordinates of the packing molecules were fixed, whereas the geometry of the central molecule was optimized using different DFT functionals in the field created by packing ones. This approach is assumed to recover the majority of the crystal lattice effects10,15,47–50 allowing for a direct comparison of calculated parameters to experimental X-ray data. An example of the molecular cluster for Au1 is shown in Fig. 2.
The geometries of the central molecules in the clusters were optimized using the def2-TZVPP basis set51 for all atoms, with corresponding relativistic effective core potentials (def2ECPs)52 for the metal center (ECP substituting 60 electrons for Pt and Au and 28 electrons for Rh and Pd). The packing molecules with fixed coordinates were treated using the def2SVP basis set51 for all atoms, with corresponding ECP for the metal center.
Each cluster was optimized using ten selected density functionals. The LYP-based53 functionals BLYP,54 B3LYP,55 and BHLYP56 were selected to map the effect of exact-exchange admixture on the structural parameters. Due to the known unbalanced performance of the B3LYP functional in calculations involving transition metals,57 the CAM-B3LYP,58 the B3LYP functional with corrected long-range exchange, was included in the test set. From the GGA family, PBE59 and BP8654,60 are

present because they are known to perform reasonably well for transition-metal complexes.12 The TPSSH61,62 functional is an example of a successful meta-GGA hybrid functional, and it was shown to be a very good choice for optimizing transitionmetal complexes.8,12 The PBE0,63,64 the hybrid version of the ‘‘parameter-free’’ PBE functional, has demonstrated superb performance for geometry optimization of heavy transitionmetal compounds.10,12,15 Advanced long-range corrected hybrid functional with empirical dispersion corrections, oB97XD,65,66 was recently tested with good results in structure optimizations of transition-metal complexes.67 Furthermore, a recent addition to the family of ‘‘Minnesota functionals’’, the so-called screened-exchange density functional, MN12SX,68 which was reported to offer very good performance in calculating various properties,68 is included. The effect of the D3 dispersion correction69 on the geometry was tested in several cases, see Fig. S1 (ESI†). The observed effect was found negligible for the PBE0 functional (it even deteriorates the results by 0.2 pm for PBE and BP86), and the D3 correction is thus not further considered in this work.
A convergence of the basis set was estimated in the series def2-SVP - def2-QZVPP using the PBE0 functional, see Fig. 3. The def2-TZVPP basis provides results almost identical (difference o1 pm) to those of much larger (and computationally demanding) QZ bases, confirming our previous results.10
Based on the performance of individual functionals (vide infra), the basis-set effects (Fig. 3), computational costs, and the marginal effects of dispersion correction, the PBE0/def2-TZVPP/ ECP approach was used to optimize molecular geometries in the production calculations. For all the in-solution optimizations discussed below, the COSMO (COnductor-like Screening MOdel)16 solvent model was used, see Section 2.3.
2.3 NMR chemical shift calculations – solvent effects and method selection
Geometry. To calculate in-solution NMR chemical shifts, individual single molecules were optimized using the PBE0 with the def2-TZVPP basis set for all atoms and corresponding ECP for the transition metal (the approach calibrated in Section 2.2

Fig. 2 Molecular cluster of Au1 with central (in color) and packing (in orange) molecules.

Fig. 3 Basis-set effects on the M–L distance in the optimization of compound Au1 using the PBE0 functional with def2-ECP.

This journal is © the Owner Societies 2015

Phys. Chem. Chem. Phys., 2015, 17, 24944--24955 | 24947

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

PCCP

View Article Online
Paper

with respect to the experimental X-ray reference geometries). In addition, the implicit COSMO solvent model16 was adopted throughout the following geometry optimizations. The standard COSMO parameters as implemented in the ADF program (DMSO – the relative dielectric constant e = 46.7, the radius of the solvent probe rsolv = 3.04 Å; CHCl3 – e = 4.8, rsolv = 3.17 Å; THF – e = 7.6, rsolv = 3.18 Å; benzene – e = 2.3, rsolv = 3.28 Å) were used.
Two-component (2c) relativistic NMR calculations. Twocomponent (2c) relativistic NMR calculations were performed using the SO-ZORA approach (as implemented in the ADF package)24,70,71 and PBE, PBE0, or PBE0 with an exact-exchange admixture of 40% (PBE-40). The choice of the functional is justified by our previous studies10,11 and the dependence on the basis set was tested in this study (Fig. 4). As the difference between the NMR chemical shifts (13C and 15N) for Au1 calculated using TZP and QZ4P basis sets is B1 ppm (Fig. 4), the standard TZP basis set from the ADF library was used for all atoms in the production calculations. All calculations were performed in vacuo as well as using the COSMO solvent model in order to estimate the effect of the solvent on the NMR chemical shifts. The gauge-origin dependence was handled using gaugeincluding atomic orbitals (GIAOs).72
Four-component (4c) relativistic NMR calculations. Fourcomponent (4c) relativistic NMR calculations were performed using the Dirac–Coulomb Hamiltonian (as implemented in the ReSpect code)32,33,73 using PBE, PBE0, and PBE-40 functionals, with the uncontracted pc-2 basis set74,75 for the light atoms and the uncontracted Dyall’s valence triple-zeta (dyall-vtz) basis set76–78 for the central metal. All shielding calculations using the ReSpect code were performed in vacuo; the solvent corrections were estimated at the SO-ZORA level as (sCOSMO À svacuo) and added to the NMR chemical shifts calculated at the fourcomponent mDKS level. The PBE calculations were performed using GIAOs to treat the gauge-origin dependence whereas for the PBE0 and PBE-40 functionals, the Common Gauge Origin (CGO) approximation, with the GO placed on the metal center, was used. Because the CGO and GIAO calculations are directly

compared, we used empirical CGO - GIAO correction to account for the basis set incompleteness error in the CGO calculations. This correction is obtained for each atom as a difference between shielding calculated using GIAO and CGO approximations at the PBE level (sPBE(CGO) À sPBE(GIAO)) and was added to the shielding constant calculated at the PBE0 or PBE-40 level. Inclusion of the ‘‘CGO correction’’ significantly improves the agreement with experiment (see Section 3.2 and ESI†).
Effect of the kernel of the exchange–correlation functional ( fxc) on the NMR chemical shifts. Non-vanishing two-electron contributions ( f2e) to the NMR shielding response equations consist of the exchange–correlation ( fxc) and exact-exchange ( fex) kernels ( f2e = fxc + fex). In the absence of spin–orbit (SO) effects, there is no contribution from the kernel of the exchange–correlation functional ( fxc) to the NMR shielding calculations. Thus in the framework of pure DFT, NMR shielding constants can be calculated in non-relativistic or scalarrelativistic theories using simple sum-over-states expressions (note that in the Hartree–Fock or hybrid DFT framework, the kernel arising from the exact exchange ( fex) cannot be neglected at any relativistic or non-relativistic level). However, in the presence of the SO coupling, the contribution from the exchange–correlation kernel ( fxc) is non-vanishing and can give significant contributions to the NMR shielding constants, especially for the light atom (LA) directly bonded to the heavy-atom center. For example, in the case of hydrogen iodide, the exchange–correlation kernel contributes more than 3 ppm to the total 1H NMR shielding constant.32 Because this effect is missing from the calculations of standard reference compounds (e.g., tetramethylsilane for 13C, ammonia for 15N) it fully manifests itself in the relativistic NMR chemical shift calculations.
To determine the importance of fxc, 4c calculations using the ReSpect code were performed also without fxc, and the corresponding correction to the NMR chemical shift, DdXkeCrnel, was obtained as a difference between the NMR chemical shift calculated with and without kernel contribution at a given level of theory. With DdXkeCrnel being unaffected by the gauge origin (GIAO vs. CGO), (tested using GIAO/CGO calculations at the 4c-PBE level, data not shown) we assume its additivity for the 2c SO-ZORA results (vide infra).
NMR chemical shift referencing. To reduce systematic errors in calculating the NMR chemical shifts, the following secondary references were used. For 15N, N-methylpyridinium in DMSO (d = 203.5 ppm relative to liquid ammonia);79,80 for 13C, benzene in benzene (d = 127.83 ppm relative to TMS);4 and for 31P, tert-butylphosphine in CHCl3 (d = 63.0 ppm relative to 85% H3PO4).81 The NMR chemical shifts, di, were obtained using the equation:

di = sref À si + dref

Fig. 4 Basis set effects on the 13C and 15N NMR chemical shifts in compound Au1 calculated at the SO-ZORA level using the PBE-40 functional (PBE0/def2-TZVPP/def2-ECP geometry).

where sref is the calculated NMR shielding of the corresponding nucleus in the reference compound, si is the calculated NMR shielding of the investigated nucleus, and dref is the experimental NMR chemical shift of the secondary reference relative
to the primary standard.

24948 | Phys. Chem. Chem. Phys., 2015, 17, 24944--24955

This journal is © the Owner Societies 2015

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Paper

View Article Online
PCCP

3. Results and discussion
3.1 Geometry optimization – estimating the structural effects on the NMR chemical shift
The molecular clusters were built based on X-ray diffraction geometries (for details, see Methods Section) and the central molecule of each cluster was optimized in the field created by fixed packing molecules. The performance of individual functionals was evaluated using bond distances for non-hydrogen atoms (comparison with the experimental values for 161 bonds) with special focus on 28 bonds between the central metal and the ligand atoms. Root-mean-square deviations (RMSD) between calculated and experimental values for all 161 non-hydrogen bonds (RMSD161) and for 28 bonds around the metal center (RMSD28) from all molecules of our set are summarized in Fig. 5.
Poor results are obtained by using the BLYP functional, which greatly overestimates the M–L bonds (RMSD28 = 5.2 pm) and thus should not be used for structure optimizations of transition-metal complexes.10 Increasing the exact-exchange admixture in the LYPbased functionals improves the correlation with experimental data for M–L bonds [BLYP (0%, 5.2 pm) - B3LYP (20%, 3.3 pm) BHLYP (50%, 2.4 pm)] notably. When all bond lengths are considered, both B3LYP and BHLYP are performing with the same level of accuracy (RMSD161 B 2 pm), which is due to the fact that the impact of non-metal ligand bonds (mostly organics, in which B3LYP is known to perform excellently) compensates for the inaccuracy in describing the M–L bonds. The best of the LYP-based functionals, BHLYP, produced moderately accurate results (RMSD28 = 2.4 pm, RMSD161 = 1.9 pm) compared to the rest of the functionals. However, almost identical performance was observed for PBE and BP86, at notably lower computational costs. Both of these GGA functionals can be recommended for geometry optimization of transition-metal complexes in cases where the use of hybrid GGAs is inconvenient, e.g. for very large molecules or molecular clusters.
The oB97XD functional produces good overall results (RMSD161 = 1.7 pm), whereas the description of the M–L bonding is comparable to BHLYP, BP86, and PBE.

The importance of corrected long-range exchange in the CAM-B3LYP functional to optimize the geometry of the transitionmetal complexes is obvious from a comparison of the B3LYP and CAM-B3LYP results. The CAM version improves the results of B3LYP by 40% for M–L bonds (RMSD28 is 2.0 pm vs. 3.3 pm for B3LYP) and by about 20% for all non-hydrogen bonds (RMSD161 is 1.7 pm vs. 2.1), see Fig. 5.
The best results for the tested set were produced by three functionals – MN12SX (RMSD28 = 1.7 pm and RMSD161 = 1.6 pm), TPSSH (with RMSD28 = RMSD161 = 1.6 pm), and PBE0, which performs excellently in describing the M–L bonds (RMSD28 = 1.2 pm, RMSD161 = 1.5 pm). For a more detailed analysis of the results according to the atom type (Fig. S2 and S3, ESI†) and for a comparison between the cluster approach and the in vacuo results (Fig. S4), see ESI.†
To estimate the magnitude of the structural effects on the NMR chemical shifts (NMR CS) separately from all the other factors, we calculated the NMR CS of metal-bonded light atoms (highlighted in green in Fig. 1) using non-relativistic and SO-ZORA approaches at the PBE-40/TZP level in vacuo for BLYP, B3LYP, PBE0, and MN12SX geometries (def2-TZVPP/def2-ECP basis set) optimized in vacuo (see Table 2). The difference between the NMR CS calculated using BLYP and MN12SX functionals (Ddgeom, see Table 2) ranges between 5–14 ppm at the non-relativistic level, and 2–22 ppm at the 2c SO-ZORA level. Also, the deficiency of the popular B3LYP functional in describing the M–L bond length can result in a 19 ppm difference in the relativistic 13C NMR chemical shift (e.g., 13C in Au1 calculated at the 2c SO-ZORA level). We conclude that the inaccurate geometry can be reflected in the 13C and 15N NMR CS deviations, amounting easily to as much as 10–20 ppm. These deviations are more pronounced for the NMR chemical shifts calculated at the relativistic level (SO-ZORA) compared to those calculated using the non-relativistic approach (Table 2). This relativistic effect is particularly evident for the 13C NMR CS in Au1 due to the more covalent character of the M–C bond compared to that of the M–N bond (Table 2).11,82 As expected, the PBE0, TPSSH, and MN12SX functionals, producing the most reliable geometries in the test set, give very similar calculated NMR CS.
In this work, the PBE0/def2-TZVPP/def2-ECP approach was selected for all further geometry optimizations in Sections 3.2–3.5.

Fig. 5 The averaged RMSDs (in pm) for interatomic distances relative to the experimental X-ray data, calculated using various density functionals. The M–L bonds (28) are in red, and all non-hydrogen bonds (161) are in green.

3.2 Relativistic calculations of NMR chemical shifts (NMR CS) – role of exact-exchange (EE) admixture
All the NMR chemical shifts were calculated for single molecules optimized using the PBE0 functional (for details, see Methods Section and Section 3.1), however, only those that require relativistic treatment, namely atoms significantly affected by the HALA effects from the metal atom (see Fig. 1 and Introduction part), are reported here for clarity. It is known that the calculated NMR CS are significantly modulated by the amount of exact-exchange admixture in the hybrid GGA functional.11,27 To analyze the role of exact-exchange admixture, we performed relativistic calculations of the light-atom NMR CS using 2c SO-ZORA and 4c mDKS with the PBE, PBE0,

This journal is © the Owner Societies 2015

Phys. Chem. Chem. Phys., 2015, 17, 24944--24955 | 24949

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

PCCP

View Article Online
Paper

Table 2 NMR chemical shifts calculated at non-relativistic and relativistic (SO-ZORA) levels of theory using PBE-40/TZP in vacuo for compounds Pt1 and Au1 optimized using BLYP, B3LYP, PBE0, or MN12SX functionals (def2-TZVPP/def2-ECP). NMR chemical shifts in ppm (for shift referencing, see Methods Section)

Non-relativistic

BLYP

B3LYP

PBE0

MN12SX

Ddgeoma

SO-ZORA BLYP

B3LYP

PBE0

Pt1

13C

182.8

179.1

177.0

177.0

5.8

141.6

141.1

143.2

15N

270.0

261.6

257.1

256.6

13.4

234.8

229.3

228.8

Au1

13C

176.6

172.1

168.9

168.6

8.0

147.0

144.6

125.6

15N

250.4

242.3

237.2

236.8

13.6

244.2

236.0

230.5

a Ddgeom calculated as d(BLYP) – d(MN12SX); note that d(BLYP) – d(PBE0) would result in almost identical numbers.

MN12SX
143.4 228.2
125.3 230.2

Ddgeoma
À1.8 6.6
21.7 14.0

and PBE-40 functionals. Results (in vacuo), plotted as deviations from the experimental values (Dd), are shown in Fig. 6 and more details are given in Tables S1 and S2 (ESI†); for the influence of the CGO correction in the 4c mDKS PBE0 and PBE-40 calculations, see Fig. S5 (ESI†).
Generally, increasing the EE admixture in the functional increases the nuclear magnetic shielding of the light atoms bound to the metal center in platinum and gold compounds considerably.10,82 However, this effect is marginal for palladium and even reversed for 31P in the rhodium complex. The increased EE admixture improves the results for some cases already in vacuo, particularly for 13C NMR CS, where the solvent effects obviously play a less significant role compared to 15N NMR CS.80,83–85 For instance, the PBE functional results in Dd(13C) = 21.4 ppm for Pt1-C, whereas PBE-40 gives Dd(13C) = 1.5 ppm (see Fig. 6). Conversely, the improvement is small or negligible for 15N NMR CS, where solvent effects are known to play a significant role and the effects of the EE admixture are less important (see the palladium and gold complexes in Fig. 6).
3.3 Solvent correction to the NMR CS
Due to the sensitivity of the NMR CS (particularly 15N)85 to the solvent, it is necessary to include bulk environmental effects in

the NMR calculations.86 In this study, we analyzed the effect of the implicit solvent model on the NMR CS of LA (13C and 15N). Solvent corrections to the NMR CS (Ddsolv), calculated at the SO-ZORA level using the PBE, PBE0, and PBE-40 functionals with
the COSMO model for the solvent, are summarized in Table 3. The implicit solvent improved the RMSDSPOBE-ZORA from
21.5 ppm in vacuo to 14 ppm with COSMO (for RMSDm PBDEKS, 18.6 ppm in vacuo and 11.4 ppm including Ddsolv). In parallel, the RMSDSPOBE-Z0ORA was reduced from 16.5 ppm to 4.4 ppm (RMSDm PBDEK0S from 17.5 ppm to 5.1) and the RMSDSPOBE-Z-4O0RA from 14.6 ppm to 4.5 ppm (RMSDm PBDEK-4S0 15.1 ppm to 5.2 ppm).
Note that the magnitude of Ddsolv is increasing considerably with increasing exact-exchange admixture in the functional; the
more exact exchange in the functional, the more important the
solvent correction is, which is partially reflected in the decrease
of the RSMD for the individual approaches. Whereas PBE results
for both the 2c SO-ZORA and 4c mDKS approaches still deviate
significantly from the experimental values, which is partially due to the differences in the description of DdSO by the individual functionals (see Table 4), PBE0 and PBE-40 results are consider-
ably closer to the experimental data. For the differences between
calculated (PBE, PBE0, or PBE-40 functional using implicit
solvent) and experimental NMR CS, see Fig. 7.

Fig. 6 Total NMR chemical shifts calculated using 2c SO-ZORA and 4c mDKS methods at various levels of theory in vacuo (for details, see Methods section). The calculated values are reported as deviations (Dd) from the experimental NMR chemical shifts.

24950 | Phys. Chem. Chem. Phys., 2015, 17, 24944--24955

This journal is © the Owner Societies 2015

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Paper

View Article Online
PCCP

Table 3 The Ddsolv calculated at PBE, PBE0, and PBE-40 levels of theory using the 2c SO-ZORA method

Ddsolv
Atom
Pt1-N Pt2-N Pd1-N Pd2-N Au1-N Au2-N Pt1-C Au1-C Rh-C Rh-P

PBE
À5.0 À1.3 À4.6 À6.0 À12.3 À14.5 À4.1 À2.5 À0.2 +2.7

PBE0
À5.7 À4.4 À14.6 À16.8 À15.7 À17.9 À4.4 À3.4 À0.3 +2.2

PBE-40
À6.1 À6.3 À20.4 À22.2 À17.5 À20.1 À4.5 À4.1 À0.2 +2.0

Although at the SO-ZORA relativistic level with solvent
correction, the PBE0 and PBE-40 calculations result in almost identical averaged RMSD values (RMSDSPOBE-Z0ORA = 4.4 ppm and RMSDSPOBE-Z-4O0RA = 4.5 ppm), the PBE-40 provides better results in 7 out of 10 systems, which is in agreement with previous studies.10,11 The better average RMSD achieved by the PBE0
functional can be ascribed to quite poor results of the PBE-40
approach for Pt2-N, NH2 group in oxaliplatin, which is probably

involved in protonation equilibria affecting its experimental 15N NMR CS, and overestimated 31P NMR chemical shift for Rh-P. When the two abovementioned problematic Pt2-N and Rh-P cases are excluded, RMSDSPOBE-Z-4O0RA decreases to 2.4 ppm, whereas RMSDSPOBE-Z0ORA = 4.7 ppm.
Similar results were achieved by the 4c mDKS approach, where slightly better RMSD was obtained by using PBE0 (RMSDm PBDEK0S = 5.1 ppm) than by PBE-40 (RMSDm PBDEK-4S0 = 5.2 ppm), however, one should keep in mind that this was achieved by applying several external corrections (solvent, CGO).
Given the good results obtained for PBE0 and PBE-40 with an implicit solvent correction, the use of explicit solvent models would not bring any considerable improvement at reasonable costs for routine NMR applications; therefore, we use the COSMO implicit solvent model for calculating the NMR CS in subsequent sections.
3.4 Spin–orbit contribution to the NMR CS
It is well known that most of the HA-LA effect is due to spin– orbit coupling.86–89 Hence it is vitally important to observe and analyze the difference between calculations with and without the SO contribution (DdSO), i.e., the SO in ZORA as well as the

Table 4 The DdSO calculated at SO-ZORA and mDKS levels of theory using PBE, PBE0, and PBE-40 functionals in vacuo

Atom

DdSO (ZORA)

PBE

PBE0

PBE-40

DDdSOa

DDdtotalb

DdSO (mDKS)

PBE

PBE0

Pt1-N Pt2-N Pd1-N Pd2-N Au1-N Au2-N Pt1-C Au1-C Rh-C Rh-P

À9.3 À24.5
À5.1 À7.1 À2.7
6.0 À14.0
À8.0 À3.1 À17.0

À15.2 À27.5
À9.3 À11.5
À3.5
5.6 À23.2 À12.1
À4.3 À22.0

À19.1 À30.3 À12.4 À15.1
À4.5
5.0 À29.1 À16.2
À5.9 À24.9

À9.8 À5.8 À7.3 À8.0 À1.8 À1.0 À15.1 À8.2 À2.8 À7.9

À14.4 À27.2
+2.4 À2.7 À4.6 À3.6 À19.9 À14.9 À5.5
+14.9

À10.7 À24.7
À5.6 À8.6
0.6
9.9 À18.3
À9.0 À2.5 À5.1

À15.7 À27.5
À9.5 À13.8
À1.5
7.3 À28.1 À14.3
À4.5 À10.3

a Calculated as DDdSO = DdSPOBE-40 À DdSPOBE. b Calculated as DDdtotal = DdtPoBtEal-40 À DdtPoBtEal.

PBE-40
À19.0 À29.6 À12.6 À17.6
À3.0 6.3
À29.5 À18.4
À5.8 À12.7

DDdSOa
À8.3 À4.0 À7.0 À9.0 À3.6 À3.6 À11.2 À9.4 À3.3 À7.6

DDdtotalb
À9.0 À23.3
+5.7 À1.0 À1.9 À3.9 À14.1 À16.8 À4.8 +15.3

Fig. 7 NMR chemical shifts calculated using 2c SO-ZORA and 4c mDKS methods at various levels of theory with solvent correction (Ddsolv). Results are reported as standard deviation from the experimental value.

This journal is © the Owner Societies 2015

Phys. Chem. Chem. Phys., 2015, 17, 24944--24955 | 24951

PCCP

View Article Online
Paper

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Fig. 8 The NMR chemical shifts calculated at the 2c level using PBE0 (red) and PBE-40 (green) functionals (a) without DdXkeCrnel contribution (dark colors) and (b) with the addition of DdXkeCrnel values (light colors).

difference between the full mDKS and that without SO. Analogously to Ddsolv, the shielding character of DdSO is increasing considerably with increasing EE admixture, see Table 4.
The increase in the shielding character of DdSO is an important factor for explaining the differences between the values calculated at the PBE, PBE0, and PBE-40 levels (vide infra). Comparing DdSO calculated at these three levels of theory with total NMR chemical shift differences (in vacuo, Table S1, ESI†) indicates that DDdSO (Table 4) induced by altering the exact-exchange admixture (40 - 0%) can represent a significant portion of the DDdtotal.
The DdSO values in Table 4 calculated using the ZORA and mDKS methods are of a comparable magnitude. Generally, the DdSO values vary between +5 ppm for Au2 (deshielding) and À30 ppm for Pt2 (strong shielding contribution). It should be noted that the electronic factors responsible for the SO deshielding (Au2-N) or small DdSO (Au1-N) in the gold complexes compared to their platinum counterparts were interpreted very recently by our group as a result of the involvement of 6p orbitals in Au–N bonding82 caused by the electrostatic potential of Au3+ and large scalar-relativistic effects90 found in the gold complexes.
3.5 The effect of response exchange–correlation kernel
The kernel-corrected SO-ZORA approach for NMR CS calculations is currently missing from the available code (ADF 2014). Therefore, we are using the 4c mDKS approach at the PBE, PBE0, and PBE-40 levels to estimate an approximate influence of terms arising from the missing XC response kernel ( fxc) to the NMR chemical shift (see Table S3 in ESI†), denoted as DdXkeCrnel. The DdXkeCrnel was found to range between +0.3 and À5.2 ppm, with an average value of about À2.1 ppm for PBE and between +2.7 and À5.3 ppm (À1.8 ppm in average) for PBE40. Similarly to DdSO,82 the DdXkeCrnel values cause deshielding in Au complexes, but they have shielding character for the Pt, Pd, and Rh compounds (see Table S3, ESI†). On average, the exchange–correlation kernel corresponds to 15–20% of DdSO found for a given spectator atom.

Table 5 Difference between DdSO calculated using the PBE-40 and PBE0 functionals at the 2c SO-ZORA level (DDdSO). DdXkeCrnel for the PBE0 level
calculated using the 4c mDKS approach

Atom

DDdSOa

DdXkeCrnelb

Pt1-N Pt2-N Pd1-N Pd2-N Au1-N Au2-N Pt1-C Au1-C Rh-C Rh-P

À3.9 À2.8 À3.1 À3.6 À1.0 À0.6 À5.9 À4.1 À1.6 À2.9

À3.0 À3.2 À0.8 À1.6
0.9
2.7 À5.6 À3.6 À1.3 À2.4

a Calculated at the SO-ZORA level as DDdSO = DdSPOBE0 À DdSPOBE-40. b DdXkeCrnel estimated at the mDKS level by using the PBE0 functional.

To illustrate the role of the missing fxc terms, appropriate DdXkeCrnel values were added to the 2c SO-ZORA NMR chemical shifts calculated at the PBE0 and PBE-40 levels, see Fig. 8.
Upon addition of DdXkeCrnel, the majority of the PBE0 results were improved, with the RMSDSPOBE-Z0ORA reduced to 3.4 ppm, whereas the PBE-40 results are uniformly deteriorated (RMSDSPOBE-Z-4O0RA = 5.8 ppm). The most important fact is revealed, however, when the differences in DdSO caused by the increase of the EE admixture from 25% to 40% are compared with the DdXkeCrnel calculated at the PBE0 level, see Table 5. Our observations suggest that the absence of the response fxc term in the current SO-ZORA implementation (ADF2014) results in DdXkeCrnel deviations that are almost perfectly counterbalanced by the artificial DdSO rise introduced by increasing the EE admixture in the functional (25 - 40% in PBE0), DDdSO in Table 5.
3.6 A summary: importance of the individual effects on the total NMR CS
To demonstrate the importance of the individual effects on the total NMR CS, we performed a comparative stepwise analysis. This includes the evaluation of the role of geometry, solvent,

24952 | Phys. Chem. Chem. Phys., 2015, 17, 24944--24955

This journal is © the Owner Societies 2015

Paper

View Article Online
PCCP

Open Access Article. Published on 25 August 2015. Downloaded on 7/8/2022 4:58:02 PM. This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence.

Fig. 9 The deviations of calculated 13C (blue) and 15N (green) NMR CS from experimental values for compound Au1. The geometry of Au1 was optimized using PBE0/def2-TZVPP/def2-ECP/COSMODMSO with an exception of A where the BLYP/def2-TZVPP/def2-ECP/COSMODMSO approach was used. The NMR chemical shifts were calculated in ADF at various levels of theory: (A) PBE/TZP, (B) PBE/TZP, (C) PBE/TZP/COSMODMSO, (D) PBE/TZP/ZORA/ COSMODMSO, (E) PBE/TZP/SO-ZORA/COSMODMSO, and (F) PBE-40/TZP/SO-ZORA/COSMODMSO. The differences between neighboring columns represent the effects on the NMR chemical shifts of geometry (GEOM), solvent (SOLV), scalar relativity (SR), spin–orbit coupling (SO), and 40% admixture of exact exchange in the PBE functional (EE).

scalar relativity, spin–orbit coupling, and admixture of exact exchange in the functional as discussed in detail in Sections 3.1–3.5. The case study of 13C and 15N NMR chemical shifts for compound Au1 is shown in Fig. 9 (for additional data, see ESI†).
Although individual factors play different roles in 13C and 15N NMR CS of different compounds, the EE admixture seems to be crucial and one of the most important factors for reproducing the experimental NMR values. However, as demonstrated here, when the correct geometry and the treatment of solvent effects is combined with the proper description of relativity (including terms arising from fxc), standard 25% exact-exchange admixture in the PBE0 functional (instead of somewhat artificial 40%) seems to be sufficient for the NMR CS calculations in squareplanar transition-metal complexes.
4. Conclusions
In this work, we evaluated the effects of molecular geometry, implicit solvent, role of the functional, and the treatment of relativistic effects on the NMR chemical shifts calculated for light atoms in the vicinity of a heavy-metal center in selected square-planar complexes. The selection of the DFT functional used for optimizing the geometry was performed by a cluster approach. The best geometries were obtained for the PBE0 functional, closely followed by TPSSH and MN12SX. The effect of geometry on the NMR chemical shifts (deviation between best and worst structure) was found to be rather large, 14 ppm for 15N and 22 ppm for 13C.
The role of an implicit solvent was found to be of considerable importance for the NMR chemical shifts in the square-planar complexes, mainly for 15N, with an average 13 ppm contribution to the total NMR chemical shifts (estimated at the SO-ZORA level using the PBE0 functional).
The most important relativistic NMR CS contribution is the DdSO, with its size varying between +5 ppm and À30 ppm (for both 2c SO-ZORA and 4c mDKS), which was also demonstrated

to be highly dependent on the amount of exact-exchange admixture in the functional.
The importance of increased exact-exchange admixture to 40% in the PBE0 functional for the 2c SO-ZORA approach was confirmed for the light spectator atoms (13C, 15N, and 31P), and rationalized as a compensation for the missing DdXkeCrnel contribution. Therefore, with the current implementation of the SO-ZORA approach used in this work, the PBE-40 functional slightly outperforms PBE0. However, accounting for the effect of the response exchange– correlation kernel, the standard PBE0 functional provides somewhat better agreement with the experiment. Undoubtedly, the effect of the XC kernel should not be neglected for systems with a large HALA effect, such as actinides.7
Finally, the relative importance of individual effects (geometry, solvent, relativistic approximation, DFT functional) to the total NMR CS is demonstrated by stepwise calculations. Currently, the use of hybrid GGA functionals (with the DFT functionals used) is mandatory for reproducing the experimental NMR data in heavy-element compounds.
In the fully relativistic mDKS approach using restricted magnetically balanced basis sets, an implementation of GIAOs for hybrid functionals as well as an implicit solvent model is still missing. However, using simple empirical corrections estimated at the 2c SO-ZORA level provides reliable NMR CS data at affordable computational costs. Undoubtedly, the fourcomponent mDKS method would outperform 2c approaches in calculating heavy-element NMR parameters.91 Our efforts in this direction will be reported elsewhere.
Acknowledgements
This work was supported by the Czech Science Foundation (15-09381S to RM and 14-03564S to MS) and carried out at CEITEC – the Central European Institute of Technology with research infrastructure supported by the project CZ.1.05/1.1.00/ 02.0068 financed by the European Regional Development Fund.

This journal is © the Owner Societies 2015

Phys. Chem. Chem. Phys., 2015, 17, 24944--24955 | 24953

Preparing to load PDF file. please wait...

0 of 0
100%
Structure, solvent, and relativistic effects on the NMR