On the geometry dependence of the nuclear magnetic resonance chemical shift of mercury in thiolate complexes: A relativistic density functional theory study

Thiolate containing mercury(II) complexes of the general formula [Hg(SR)n ]2 - n have been of great interest since the toxicity of mercury was recognized. 199Hg nuclear magnetic resonance spectroscopy (NMR) is a powerful tool for characterization of mercury complexes. In this work, the Hg shielding constants in a series of [Hg(SR)n ]2 - n complexes are therefore investigated computationally with particular emphasis on their geometry dependence. Geometry optimizations and NMR chemical shift calculations are performed at the density functional theory (DFT) level with both the zeroth-order regular approximation (ZORA) and four-component relativistic methods. The four exchange-correlation (XC) functionals PBE0, PBE, B3LYP, and BLYP are used in combination with either Dyall's Gaussian-type (GTO) or Slater-type orbitals (STOs) basis sets. Comparing ZORA and four-component calculations, one observes that the calculated shielding constants for a given molecular geometry have a constant difference of ∼ 1070 ppm. This confirms that ZORA is an acceptable relativistic method to compute NMR chemical shifts. The combinations of four-component/PBE0/v3z and ZORA/PBE0/QZ4P are applied to explore the geometry dependence of the isotropic shielding. For a given coordination number, the distance between mercury and sulfur is the key factor affecting the shielding constant, while changes in bond and dihedral angles and even different side groups have relatively little impact.

Mechanisms of Ligand Hyperfine Coupling in Transition-Metal Complexes: σ and π Transmission Pathways

Theoretical interpretation of hyperfine interactions was pioneered in the 1950s-1960s by the seminal works of McConnell, Karplus, and others for organic radicals and by Watson and Freeman for transition-metal (TM) complexes. In this work, we investigate a series of octahedral Ru(III) complexes with aromatic ligands to understand the mechanism of transmission of the spin density from the d-orbital of the metal to the s-orbitals of the ligand atoms. Spin densities and spin populations underlying ligand hyperfine couplings are analyzed in terms of π-conjugative or σ-hyperconjugative delocalization vs spin polarization based on symmetry considerations and restricted open-shell vs unrestricted wave function analysis. The transmission of spin density is shown to be most efficient in the case of symmetry-allowed π-conjugative delocalization, but when the π-conjugation is partially or fully symmetry-forbidden, it can be surpassed by σ-hyperconjugative delocalization. Despite a lower spin population of the ligand in σ-hyperconjugative transmission, the hyperfine couplings can be larger because of the direct involvement of the ligand s-orbitals in this delocalization pathway. We demonstrate a quantitative correlation between the hyperfine couplings of aromatic ligand atoms and the characteristics of the metal-ligand bond modulated by the trans substituent, a hyperfine trans effect.

Photoluminescence lifetime stability studies of β-diketonate europium complexes based phenanthroline derivatives in poly(methyl methacrylate) films

In this work, five phenanthroline derivatives substituted with different methyl groups have been selected to synthesize β-diketonate-based europium complexes to check the influence of the substitutions on the degradation effect of those complexes in poly(methyl methacrylate) (PMMA) films. The photophysical properties of Eu(III) complexes, including absorbance, excitation, and emission have been carefully investigated in solution, solid-state, and doped in PMMA film. In all these states, the complexes exhibit an impressive red emission at 614 nm with a high photoluminescence quantum yield of up to 85 %. The films have been exposed under outdoor, indoor, and dark storage stability lifetime conditions for 1200 hours. The photoluminescence measurements recorded every 400, 800, and 1200 hours demonstrated that the film containing europium complex with phenanthroline ligand substituted by a high number of methyl groups (Eu(TTA)3L5) showed good photoluminescent stability in indoor and dark conditions, and exhibited better resistance to degradation in outdoor conditions compared to other complexes. This study has proved that phenanthroline ligands could be tuned chemically leading to better stability of those types of complexes in films which can be end-used for future stable optoelectronic devices such as luminescent solar concentrators.

Dynamic and relativistic effects on Pt-Pt indirect spin-spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent-solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin-orbit Kohn-Sham (2c-ZKS) and four-component Dirac-Kohn-Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1-Pt2(NH3)4(Am)2-L2]n+ (Am = α-pyrrolidonate and pivalamidate; L = H2O, Cl-, and Br-) are studied. The computed Pt-Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N-C-O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X-) and aquahalo (L1 = H2O; L2 = X-) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Inquiring 199Hg NMR Parameters by Combining Ab Initio Molecular Dynamics and Relativistic NMR Calculations

Ab initio molecular dynamics (AIMD) sampling followed by relativistic density functional theory (DFT) 199Hg NMR calculations were performed for Hg organometallic complexes in water, dimethyl sulfoxide, and chloroform. The spin-orbit coupling, a relativistic effect, is a key factor for predicting δ(Hg) and 1J(Hg-C) accurately, in conjunction with a dynamic treatment of the systems. Good agreement between the theoretical and experimental results is reached by adopting implicit (based on a continuum model) and explicit (solvent molecules treated quantum mechanically) solvation models. Broader trends appearing in the experimental data available in the literature are reproduced by the calculations, and therefore, quantum chemistry is able to assist in the assignment and interpretation of 199Hg NMR data. Less pronounced trends, such as changes in the 199Hg chemical shift in different systems with the same atom types bound to Hg, are too weak to be predicted reliably by the current state-of-the-art theoretical methods based on AIMD sampling and relativistic DFT with hybrid functionals for NMR calculations.

Exact two-component theory becoming an efficient tool for NMR shieldings and shifts with spin-orbit coupling

We present a gauge-origin invariant exact two-component (X2C) approach within a modern density functional framework, supporting meta-generalized gradient approximations such as TPSS and range-separated hybrid functionals such as CAM-B3LYP. The complete exchange-correlation kernel is applied, including the direct contribution of the field-dependent basis functions and the reorthonormalization contribution from the perturbed overlap matrix. Additionally, the finite nucleus model is available for the electron-nucleus potential and the vector potential throughout. Efficiency is ensured by the diagonal local approximation to the unitary decoupling transformation in X2C as well as the (multipole-accelerated) resolution of the identity approximation for the Coulomb term (MARI-J, RI-J) and the seminumerical exchange approximation. Errors introduced by these approximations are assessed and found to be clearly negligible. The applicability of our implementation to large-scale calculations is demonstrated for a tin pincer-type system as well as low-valent tin and lead complexes. Here, the calculation of the Sn nuclear magnetic resonance shifts for the pincer-type ligand with about 2400 basis functions requires less than 1 h for hybrid density functionals. Further, the impact of spin-orbit coupling on the nucleus-independent chemical shifts and the corresponding ring currents of all-metal aromatic systems is studied.

How the Metal Ion Affects the 1H NMR Chemical Shift Values of Schiff Base Metal Complexes: Rationalization by DFT Calculations

The chemical shift (CS) values obtained by 1H NMR spectroscopy for the hydrogen atoms of a tetradentate N2O2-substituted Salphen ligand (H2L1) are differently shifted in its complexes of nickel(II), palladium(II), platinum(II), and zinc(II), all bearing the same charge on the metal ions. To rationalize the observed trends, DFT calculations have been performed in the implicit d6-DMSO solvent in terms of the electronic effects induced by the metal ion and of the nature and strength of the metal-N and metal-O bonds. Overall, the results obtained point out that, in the complexes involving group 10 elements, the CS values show the greater shift when considering the two hydrogen atoms at a shorter distance from the coordinated metal center and follow the decreasing metal charge in the order Ni > Pd > Pt. This trend suggests a more covalent character of the ligand-metal bonds with the increase of the metal atomic number. Furthermore, a slightly poorer agreement between experimental and calculated data is observed in the presence of the nickel(II) ion. Such discrepancy is explained by the formation of stacked oligomers, aimed at minimizing the repulsive interactions with the polar DMSO solvent.

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

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

Nature of NMR Shifts in Paramagnetic Octahedral Ru(III) Complexes with Axial Pyridine-Based Ligands

In recent decades, transition-metal coordination compounds have been extensively studied for their antitumor and antimetastatic activities. In this work, we synthesized a set of symmetric and asymmetric Ru(III) and Rh(III) coordination compounds of the general structure (Na⁺/K⁺/PPh₄⁺/LH⁺) [trans-M^IIIL(eq)_nL(ax)₂]^- (M = Ru^III or Rh^III; L(eq) = Cl, n = 4; L(eq) = ox, n = 2; L(ax) = 4-R-pyridine, R = CH₃, H, C₆H₅, COOH, CF₃, CN; L(ax) = DMSO-S) and systematically investigated their structure, stability, and NMR properties. ¹H and ¹³C NMR spectra measured at various temperatures were used to break down the total NMR shifts into the orbital (temperature-independent) and hyperfine (temperature-dependent) contributions. The hyperfine NMR shifts for paramagnetic Ru(III) compounds were analyzed in detail using relativistic density functional theory (DFT). The effects of (i) the 4-R substituent of pyridine, (ii) the axial trans ligand L(ax), and (iii) the equatorial ligands L(eq) on the distribution of spin density reflected in the "through-bond" (contact) and the "through-space" (pseudocontact) contributions to the hyperfine NMR shifts of the individual atoms of the pyridine ligands are rationalized. Further, we demonstrate the large effects of the solvent on the hyperfine NMR shifts and discuss our observations in the general context of the paramagnetic NMR spectroscopy of transition-metal complexes.

Exact Two-Component TDDFT with Simple Two-Electron Picture-Change Corrections: X-ray Absorption Spectra Near L- and M-Edges of Four-Component Quality at Two-Component Cost

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

Cationic Gold(II) Complexes: Experimental and Theoretical Study

Gold(II) complexes are rare, and their application to the catalysis of chemical transformations is underexplored. The reason is their easy oxidation or reduction to more stable gold(III) or gold(I) complexes, respectively. We explored the thermodynamics of the formation of [Au^II (L)(X)]⁺ complexes (L=ligand, X=halogen) from the corresponding gold(III) precursors and investigated their stability and spectral properties in the IR and visible range in the gas phase. The results show that the best ancillary ligands L for stabilizing gaseous [Au^II (L)(X)]⁺ complexes are bidentate and tridentate ligands with nitrogen donor atoms. The electronic structure and spectral properties of the investigated gold(II) complexes were correlated with quantum chemical calculations. The results show that the molecular and electronic structure of the gold(II) complexes as well as their spectroscopic properties are very similar to those of analogous stable copper(II) complexes.

Bonding and 13 C-NMR properties of coinage metal tris(ethylene) and tris(norbornene) complexes: Evaluation of the role of relativistic effects from DFT calculations

The π-complexes of cationic coinage metal ions (Cu(I), Ag(I), Au(I)) provide useful experimental support for understanding fundamental characteristics of bonding and ¹³ C-NMR patterns of the group 11 triad. Here, we account for the role of relativistic effects on olefin-coinage metal ion interaction for cationic, homoleptic tris-ethylene, and tris-norbornene complexes, [M(η² -C₂ H₄ )₃ ]⁺ and [M(η² -C₇ H₁₀ )₃ ]⁺ (M = Cu, Ag, Au), as representative case of studies. The M-(CC) bond strength in the cationic, tris-ethylene complexes is affected sizably for Au and to a lesser extent for Ag and Cu (48.6%, 16.7%, and 4.3%, respectively), owing to the influence on the different stabilizing terms accounting for the interaction energy in the formation of coinage metal cation-π complexes. The bonding elements provided by olefin → M σ-donation and olefin ← M π-backbonding are consequently affected, leading to a lesser covalent interaction going down in the triad if the relativistic effects are ignored. Analysis of the ¹³ C-NMR tensors provides further understanding of the observed experimental values, where the degree of backbonding charge donation to π₂ *-olefin orbital is the main influence on the observed high-field shifts in comparison to the free olefin. This donation is larger for ethylene complexes and lower for norbornene counterparts. However, the bonding energy in the later complexes is slightly stabilized given by the enhancement in the electrostatic character of the interaction. Thus, the theoretical evaluation of metal-alkene bonds, and other metal-bonding situations, benefits from the incorporation of relativistic effects even in lighter counterparts, which have an increasing role going down in the group.

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.

Insights into the value of statistical models, solvent, and relativistic effects for investigating Re complexes of 2-(4'-aminophenyl)benzothiazole: a potential spectroscopic probe

Cancer affects a major part of the worldwide population, and, to minimize deaths, the diagnosis in the early stages of the disease is fundamental. Thus, to improve diagnosis and treatment new potential spectroscopic probes are crucial. Benzothiazole derivates present antitumor properties and are highly selective and interact strongly with the enzyme phosphoinositide 3-kinase (PI3K), which was associated with cell proliferation and breast cancer cells. In this paper, the rhenium shielding tensors (¹⁸⁷Re(σ)) and hydrogen and carbon chemical shifts (¹H(δ) and ¹³C(δ)) of the Re(CO)₃(NNO) complex conjugated with 2-(4'-aminophenyl)benzothiazole (ReABT) were evaluated. A statistical HCA model was used to analyze the best DFT protocol to compute σ and δ values and to evaluate the relativistic effects, both in the basis set and Hamiltonian as well as the functionals M06L or PBE0. The best protocol was applied to obtain ¹⁸⁷Re(σ) of the ReABT complex in different environments (gas phase, solution, and in the active site of the PI3K enzyme). The results point out that ¹⁸⁷Re(σ) values of the ReABT complex change significantly when the complex is docked in the PI3K enzyme.

DFT Calculations of 31P NMR Chemical Shifts in Palladium Complexes

In this study, comparative analysis of calculated (GIAO method, DFT level) and experimental ³¹P NMR shifts for a wide range of model palladium complexes showed that, on the whole, the theory reproduces the experimental data well. The exceptions are the complexes with the P=O phosphorus, for which there is a systematic underestimation of shielding, the value of which depends on the flexibility of the basis sets, especially at the geometry optimization stage. The use of triple-ζ quality basis sets and additional polarization functions at this stage reduces the underestimation of shielding for such phosphorus atoms. To summarize, in practice, for the rapid assessment of ³¹P NMR shifts, with the exception of the P=O type, a simple PBE0/{6-311G(2d,2p); Pd(SDD)}//PBE0/{6-31+G(d); Pd(SDD)} approximation is quite acceptable (RMSE = 8.9 ppm). Optimal, from the point of view of “price–quality” ratio, is the PBE0/{6-311G(2d,2p); Pd(SDD)}//PBE0/{6-311+G(2d); Pd(SDD)} (RMSE = 8.0 ppm) and the PBE0/{def2-TZVP; Pd(SDD)}//PBE0/{6-311+G(2d); Pd(SDD)} (RMSE = 6.9 ppm) approaches. In all cases, a linear scaling procedure is necessary to minimize systematic errors.

Which DFT factors influence the accuracy of 1H, 13C and 195Pt NMR chemical shift predictions in organopolymetallic square-planar complexes? New scaling parameters for homo- and hetero-multimetallic compounds and their direct applications

Because of their chemical heterogeneity, stereochemical complexity and the presence of heavy atoms involving orbitals with high quantum number L, organopolymetallic complexes require considerable focus during their NMR spectral interpretation.

Accurate X-ray Absorption Spectra near L- and M-Edges from Relativistic Four-Component Damped Response Time-Dependent Density Functional Theory

The simulation of X-ray absorption spectra requires both scalar and spin–orbit (SO) relativistic effects to be taken into account, particularly near L- and M-edges where the SO splitting of core p and d orbitals dominates. Four-component Dirac–Coulomb Hamiltonian-based linear damped response time-dependent density functional theory (4c-DR-TDDFT) calculates spectra directly for a selected frequency region while including the relativistic effects variationally, making the method well suited for X-ray applications. In this work, we show that accurate X-ray absorption spectra near L_2,3- and M_4,5-edges of closed-shell transition metal and actinide compounds with different central atoms, ligands, and oxidation states can be obtained by means of 4c-DR-TDDFT. While the main absorption lines do not change noticeably with the basis set and geometry, the exchange–correlation functional has a strong influence with hybrid functionals performing the best. The energy shift compared to the experiment is shown to depend linearly on the amount of Hartee–Fock exchange with the optimal value being 60% for spectral regions above 1000 eV, providing relative errors below 0.2% and 2% for edge energies and SO splittings, respectively. Finally, the methodology calibrated in this work is used to reproduce the experimental L_2,3-edge X-ray absorption spectra of [RuCl₂(DMSO)₂(Im)₂] and [WCl₄(PMePh₂)₂], and resolve the broad bands into separated lines, allowing an interpretation based on ligand field theory and double point groups. These results support 4c-DR-TDDFT as a reliable method for calculating and analyzing X-ray absorption spectra of chemically interesting systems, advance the accuracy of state-of-the art relativistic DFT approaches, and provide a reference for benchmarking more approximate techniques.

The paper demonstrates that relativistic four-component TDDFT theory can reproduce and analyze experimental X-ray absorption spectra near L_2,3- and M_4,5-edges of transition metal and actinide compounds with different central atoms, ligands, and oxidation states. With variational inclusion of scalar and spin−orbit relativistic effects and hybrid functionals with an optimized amount of Hartee−Fock exchange (60%), it achieves relative errors below 0.2% and 2% for edge energies and spin−orbit (SO) splittings.

Fe@B6H6 aggregates: from simple building blocks to graphene analogue

We suggest the possibility to build graphene analogue with the planar hexacoordinate wheel-type Fe@B₆H₆ cluster as the building block through studying theoretically the geometry, stability, and electron structure of its dimer and trimer as well as the dimerization of the two trimers. Employing the dehydrogenation route to polymerization, we can obtain the hexagonal boron sheet that are partly and uniformly filled by Fe atoms in the center of the holes, achieving uniform chemical doping and a very large hexagonal-hole density. Thus, we may offer a novel cluster-assembled material for experimental chemists to construct graphene analogue.

Crystal and Substituent Effects on Paramagnetic NMR Shifts in Transition-Metal Complexes

Nuclear magnetic resonance (NMR) spectroscopy of paramagnetic molecules provides detailed information about their molecular and electron-spin structure. The paramagnetic NMR spectrum is a very rich source of information about the hyperfine interaction between the atomic nuclei and the unpaired electron density. The Fermi-contact contribution to ligand hyperfine NMR shifts is particularly informative about the nature of the metal–ligand bonding and the structural arrangements of the ligands coordinated to the metal center. In this account, we provide a detailed experimental and theoretical NMR study of compounds of Cr(III) and Cu(II) coordinated with substituted acetylacetonate (acac) ligands in the solid state. For the first time, we report the experimental observation of extremely paramagnetically deshielded ¹³C NMR resonances for these compounds in the range of 900–1200 ppm. We demonstrate an excellent agreement between the experimental NMR shifts and those calculated using relativistic density-functional theory. Crystal packing is shown to significantly influence the NMR shifts in the solid state, as demonstrated by theoretical calculations of various supramolecular clusters. The resonances are assigned to individual atoms in octahedral Cr(acac)₃ and square-planar Cu(acac)₂ compounds and interpreted by different electron configurations and magnetizations at the central metal atoms resulting in different spin delocalizations and polarizations of the ligand atoms. Further, effects of substituents on the ¹³C NMR resonance of the ipso carbon atom reaching almost 700 ppm for Cr(acac)₃ compounds are interpreted based on the analysis of Fermi-contact hyperfine contributions.

The ligand NMR shifts in paramagnetic acetylacetonato Cr(III) and Cu(II) complexes have been predicted and measured in the solid state and interpreted by relativistic DFT calculations. The effects of the metal atom, ligand, and crystal packing on the spin delocalization and polarization reflected in the Fermi-contact contribution to the hyperfine interaction are rationalized.

Phosphonate Chelators for Medicinal Metal Ions

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

Oxidized polysaccharides for anticancer-drug delivery: What is the role of structure?

Study provides an in-depth analysis of the structure-function relationship of polysaccharide anticancer drug carriers and points out benefits and potential drawbacks of differences in polysaccharide glycosidic bonding, branching and drug binding mode of the carriers. Cellulose, dextrin, dextran and hyaluronic acid have been regioselectively oxidized to respective dicarboxylated derivatives, allowing them to directly conjugate cisplatin, while preserving their major structural features intact. The structure of source polysaccharide has crucial impact on conjugation effectiveness, carrier capacity, drug release rates, in vitro cytotoxicity and cellular uptake. For example, while branched structure of dextrin-based carrier partially counter the undesirable initial burst release, it also attenuates the cellular uptake and the cytotoxicity of carried drug. Linear polysaccharides containing β-(1→4) glycosidic bonds and oxidized at C2 and C3 (cellulose and hyaluronate) have the best overall combination of structural features for improved drug delivery applications including potentiation of the cisplatin efficacy towards malignances.

1H and 195Pt NMR prediction for inclusion compounds formed by cisplatin and oxidized carbon nanostructures

Prediction of NMR chemical shifts can assist experimentalists in the characterization of drug delivery systems based on carbon nanocomposites. Chemical shifts are strongly correlated to the nucleus position and its chemical neighborhood. Therefore, to predict structures and NMR properties of complex chemical models, choosing a more consistent theoretical level capable of providing more realistic results and moderate computational demand is a major challenge. In this work, we predicted the NMR spectra of inclusion compounds formed by cisplatin (cDDP) and an oxidized carbon nanotube (CNTox) and nanocone (CNCox) considered by specialists as potential drug delivery systems. The ¹⁹⁵Pt NMR chemical shifts calculated at the DFT level with the new relativistic NMR-DKH basis set were −2314 ppm and −2192 ppm for cDDP@CNTox and cDDP@CNCox complexes, respectively, which are both high-field shifted relative to the free cDDP (−2110 ppm). ¹H NMR chemical shifts are also sensitive to the inclusion process. The H (NH₃) signals are found on average at +4.3 (cDDP), −5.1 (cDDP@CNTox) and +6.6 ppm (cDDP@CNCox). Interestingly, despite the similar inclusion modes in CNTox and CNCox cavities, the ¹H NMR shifts were in opposite directions. A possible reason might be the higher stability of cDDP@CNTox (ΔE_F = −19.9 kcal mol⁻¹) than that of cDDP@CNCox (ΔE_F = −5.7 kcal mol⁻¹), which suggests a short guest–host contact in the former and consequently, a more efficient shielding of hydrogen atoms due to the electron-rich carbon structure. These results may be helpful as comparison data in the NMR spectra assignment in solution and the inclusion compounds' structural elucidation.

¹⁹⁵Pt NMR chemical shifts of the cisplatin molecule in the inclusion complex formed by oxidized carbon nanotubes.

High denticity oxinate-linear-backbone chelating ligand for diagnostic radiometal ions [111In]In3+ and [89Zr]Zr4+

A potentially decadentate oxinate-containing ligand was synthesized and assessed through solution thermodynamics studies, concentration dependent radiolabeling and serum stability assays with [^nat/111In]In³⁺ and [^nat/89Zr]Zr⁴⁺.

Comprehensive Benchmark Study on the Calculation of 29Si NMR Chemical Shifts

A comprehensive and diverse benchmark set for the calculation of ²⁹Si NMR chemical shifts is presented. The SiS146 set includes 100 silicon containing compounds with 146 experimentally determined reference ²⁹Si NMR chemical shifts measured in nine different solvents in a range from -400 to +828 ppm. Silicon atoms bound to main group elements as well as transition metals with coordination numbers of 2-6 in various bonding patterns including multiple bonds and coordinative and aromatic bonding are represented. The performance of various common and specialized density functional approximations including (meta-)GGA, hybrid, and double-hybrid functionals in combination with different AO basis sets and for differently optimized geometries is evaluated. The role of scalar-relativistic effects is further investigated by inclusion of the zeroth order regular approximation (ZORA) method into the calculations. GGA density functional approximations (DFAs) are found to outperform hybrid DFAs with B97-D3 performing best with an MAD of 7.2 ppm for the subset including only light atoms (Z < 18), while TPSSh is the best tested hybrid functional with an MAD of 10.3 ppm. For ²⁹Si cores in the vicinity of heavier atoms, the application of ZORA proved indispensable. Inclusion of spin-orbit effects into the ²⁹Si NMR chemical shift calculation decreases the mean absolute deviations by up to 74% compared to calculations applying effective core potentials.

First-Principles Calculation of 1H NMR Chemical Shifts of Complex Metal Polyhydrides: The Essential Inclusion of Relativity and Dynamics

1H NMR spectroscopy has become an important technique for the characterization of transition-metal hydride complexes, whose metal-bound hydrides are often difficult to locate by X-ray diffraction. In this regard, the accurate prediction of 1H NMR chemical shifts provides a useful, but challenging, strategy to help in the interpretation of the experimental spectra. In this work, we establish a density-functional-theory protocol that includes relativistic, solvent, and dynamic effects at a high level of theory, allowing us to report an accurate and reliable interpretation of 1H NMR hydride chemical shifts of iridium polyhydride complexes. In particular, we have studied in detail the hydride chemical shifts of the [Ir6(IMe)8(CO)2H14]2+ complex in order to validate previous assignments. The computed 1H NMR chemical shifts are strongly dependent on the relativistic treatment, the choice of the DFT exchange-correlation functional, and the conformational dynamics. By combining a fully relativistic four-component electronic-structure treatment with ab initio molecular dynamics, we were able to reliably model both the terminal and bridging hydride chemical shifts and to show that two NMR hydride signals were inversely assigned in the experiment.

The accuracy challenge of the DFT-based molecular assignment of 13C MAS NMR characterization of surface intermediates in zeolite catalysis

The influence of the model and method choice on the DFT predicted ¹³C NMR chemical shifts of zeolite surface methoxide species has been systematically analyzed. Twelve ¹³C NMR chemical shift calculation protocols on full periodic and hybrid periodic-cluster DFT calculations with varied structural relaxation procedures are examined. The primary assessment of the accuracy of the computational protocols has been carried out for the Si-O(CH₃)-Al surface methoxide species in ZSM-5 zeolite with well-defined experimental NMR parameters (chemical shift, δ(¹³C) value) as a reference. Different configurations of these surface intermediates and their location inside the ZSM-5 pores are considered explicitly. The predicted δ value deviates by up to ±0.8 ppm from the experimental value of 59 ppm due to the varied confinement of the methoxide species at different zeolite sites (model accuracy). The choice of the exchange-correlation functional (method accuracy) introduces ±1.5 ppm uncertainty in the computed chemical shifts. The accuracy of the predicted ¹³C NMR chemical shifts for the computational assignment of spectral characteristics of zeolite intermediates has been further analyzed by considering the potential intermediate species formed upon methane activation by Cu/ZSM-5 zeolite. The presence of Cu species in the vicinity of surface methoxide increases the prediction uncertainty to ±2.5 ppm. The full geometry relaxation of the local environment of an active site at an appropriate level of theory is critical to ensure a good agreement between the experimental and computed NMR data. Chemical shifts (δ) calculated via full geometry relaxation of a cluster model of a relevant portion of the zeolite lattice site are in the best agreement with the experimental values. Our analysis indicates that the full geometry optimization of a cluster model at the PBE0-D3/6-311G(d,p) level of theory followed by GIAO/PBE0-D3/aug-cc-pVDZ calculations is the most suitable approach for the calculation of ¹³C chemical shifts of zeolite surface intermediates.

Towards Accurate Predictions of Proton NMR Spectroscopic Parameters in Molecular Solids

The factors contributing to the accuracy of quantum-chemical calculations for the prediction of proton NMR chemical shifts in molecular solids are systematically investigated. Proton chemical shifts of six solid amino acids with hydrogen atoms in various bonding environments (CH, CH₂ , CH₃ , OH, SH and NH₃ ) were determined experimentally using ultra-fast magic-angle spinning and proton-detected 2D NMR experiments. The standard DFT method commonly used for the calculations of NMR parameters of solids is shown to provide chemical shifts that deviate from experiment by up to 1.5 ppm. The effects of the computational level (hybrid DFT functional, coupled-cluster calculation, inclusion of relativistic spin-orbit coupling) are thoroughly discussed. The effect of molecular dynamics and nuclear quantum effects are investigated using path-integral molecular dynamics (PIMD) simulations. It is demonstrated that the accuracy of the calculated proton chemical shifts is significantly better when these effects are included in the calculations.

Performance of the LRESC Model on top of DFT Functionals for Relativistic NMR Shielding Calculations

The linear response within the elimination of the small component model (LRESC) is an insightful and computationally efficient method for including relativistic effects on molecular properties like the nuclear magnetic shielding constants, spin-rotation constant, g-tensor, and electric field gradient of heavy atom containing molecules with atoms belonging up to the sixth row of the periodic table. One of its main advantages is its capacity to analyze the electronic origin of the different relativistic correcting terms. Until now, it was always applied on top of Hartree-Fock ground-state wave functions (LRESC/HF) to calculate and analyze NMR shieldings. In this work, we show the performance of the LRESC formalism on top of some density functional theory (DFT) functionals to compute tin shielding constants in SnX₄ (X = H, F, Cl, Br, I) molecular systems. We analyze the performance of each LRESC/DFT scheme on reproducing the electronic mechanisms of the shieldings, taking as a benchmark the results of relativistic calculations at the RPA level of approach (4c/RPA). As in previous works, we divide the LRESC relativistic correcting terms into two groups: core-dependent and ligand-dependent contributions. It is shown here that core-dependent corrections are well-reproduced for the selected DFT functionals, but some differences arise in the ligand-dependent ones. We focus on the performance of different functionals, including the same electron correlation part but containing different amounts of HF exchange. The best results are obtained for the BHandHLYP functional (50% of HF exchange) and the worst for BLYP (0%). When the percentage of HF exchange increases, ligand-dependent contributions are better described, and the final LRESC/DFT results are closer to those obtained with LRESC/HF and 4c/RPA methods. The spin-orbit correction to the shielding constant is one of the main ligand-dependent contributions (there are two more) with total value depending on the amount of HF exchange included in the functional. When the amount of HF exchange decreases, the spin-orbit contribution becomes larger, overestimating the shielding constant even when nonrelativisitc DFT values are much smaller than the nonrelativistic HF ones, as it happens for the heaviest molecular system studied here (SnI₄).

NMR absolute shielding scales and nuclear magnetic dipole moments of transition metal nuclei

This work reports new, accurate nuclear magnetic dipole moments for transition metal nuclei where the long-standing systematic error due to obsolete diamagnetic correction has been eliminated by ab initio calculations of NMR shielding constants.

Computational 1 H NMR: Part 2. Chemical applications

This is the second one of three closely interrelated reviews dealing with computation of ¹ H NMR chemical shifts and ¹ H-¹ H spin-spin coupling constants prepared for Magnetic Resonance in Chemistry. Presented in this review are some basic notes and illustrative examples of how modern computational ¹ H NMR could be used for structural elucidation and stereoelectronic studies of the medium-sized organic molecules involving saturated, unsaturated, aromatic, and heteroaromatic compounds together with their functional derivatives and coordination complexes to get deeper insight into their stereochemical structure and stereodynamic behavior.

Quantum chemical calculations of 31P NMR chemical shifts of P-donor ligands in platinum(II) complexes

This work aims to find the most suitable method that is practically applicable for the calculation of ³¹P NMR chemical shifts of Pt(II) complexes. The influence of various all-electron and ECP basis sets, DFT functionals, and solvent effects on the optimized geometry was tested. A variety of combinations of DFT functionals BP86, B3LYP, PBE0, TPSSh, CAM-B3LYP, and ωB97XD with all-electron basis sets 6-31G, 6-31G(d), 6-31G(d,p), 6-311G(d,p), and TZVP and ECP basis sets SDD, LanL2DZ, and CEP-31G were used. Chemical shielding constants were then calculated using BP86, PBE0, and B3LYP functionals in combination with the TZ2P basis. The magnitude of spin-orbit interactions was also evaluated.

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.

Relativistic Effects on a Metal-Metal Bond: Osmium Corrole Dimers

A series of metal-metal bonded osmium corrole dimers, {Os[T pXPC]}₂, were synthesized in reasonably good yields (35-46%) via the interaction of the corresponding free-base meso-tris( p-X-phenyl)corroles (H₃[T pXPC], X = CF₃, H, CH₃, and OCH₃), Os₃(CO)₁₂, and potassium carbonate in 1,2,4-trichlorobenzene under an inert atmosphere at 180 °C over several hours. The complexes are only the second class of Os corroles reported to date (the first being Os^VIN corroles) and also the second class of metal-metal bonded metallocorrole dimers (the other being Ru corrole dimers). Comparison of the X-ray structures, redox potentials, and optical spectra of analogous Ru and Os corrole dimers, along with scalar-relativistic DFT calculations, has provided an experimentally calibrated account of relativistic effects in these complexes. Three of the Os corrole dimers (X = CF₃, H, and OCH₃) were analyzed with single-crystal X-ray diffraction analysis, revealing inversion-related corrole rings with eclipsed Os-N bonds and Os-Os distances of ∼2.24 Å that are ∼0.06 Å longer than the Ru-Ru distances in the analogous Ru corrole dimers. Interestingly, a comparison of scalar-relativistic and nonrelativistic DFT calculations indicates that this difference in metal-metal bond distance does not, in fact, reflect a differential relativistic effect. For a given corrole ligand, the Ru and Os corrole dimers exhibit nearly identical oxidation potentials but dramatically different reduction potentials, with the Os values ∼0.5 V lower relative to Ru, suggesting that whereas oxidation occurs in a ligand-centered manner, reduction is substantially metal-centered, which indeed was confirmed by scalar-relativistic calculations. The calculations further indicate that approximately a third of the ∼0.5 V difference in reduction potentials can be ascribed to relativity. The somewhat muted value of this relativistic effect appears to be related to the finding that reduction of an Os corrole dimer is not exclusively metal-based but that a significant amount of spin density is delocalized over to the corrole ligand; in contrast, reduction of an Ru corrole dimer occurs exclusively on the Ru-Ru linkage. For isoelectronic complexes, the Ru and Os corrole dimers exhibit substantially different UV-vis spectra. A key difference is a strong near-UV feature of the Os series, which in energy terms is blue-shifted by ∼0.55 V relative to the analogous feature of the Ru series. TDDFT calculations suggest that this difference may be related to higher-energy Os(5d)-based LUMOs in the Os case relative to analogous MOs for Ru.

Mechanisms and Beyond: Elucidation of Fluxional Dynamics by Exchange NMR Spectroscopy

Detailed mechanistic information is crucial to our understanding of reaction pathways and selectivity. Dynamic exchange NMR techniques, in particular 2D exchange spectroscopy (EXSY) and its modifications, provide indispensable intricate information on the mechanisms of organic and inorganic reactions and other phenomena, for example, the dynamics of interfacial processes. In this Review, key results from exchange NMR studies of small molecules over the last few decades are systemised and discussed. After a brief introduction to the theory, the key types of dynamic processes are identified and fundamental examples given of intra- and intermolecular reactions, which, in turn, could involve, or not, bond-making and bond-breaking events. Following that logic, internal molecular rotation, intramolecular stereomutation and molecular recognition will first be considered because they do not typically involve bond breaking. Then, rearrangements, substitution-type reactions, cyclisations, additions and other processes affecting chemical bonds will be discussed. Finally, interfacial molecular dynamics and unexpected combinations of different types of fluxional processes will also be highlighted. How exchange NMR spectroscopy helps to identify conformational changes, coordination and molecular recognition processes as well as quantify reaction energy barriers and extract detailed mechanistic information by using reaction rate theory in conjunction with computational techniques will be shown.