Calculation of electric field gradients with the exact two-component (X2C) quasi-relativistic method and its local approximations

When calculating electric field gradients (EFGs), relativistic and electron correlation effects are crucial for obtaining accurate results, and the commonly used density functional methods produce unsatisfactory results, especially for heavy elements and/or strongly correlated systems. In this work, a stand-alone program is presented, which enables calculation of EFGs from the molecular orbitals supplied by an external high accuracy quantum chemical calculation and includes relativistic effects through the exact two-component (X2C) formalism and efficient local approximations to it. Application to BiN and BiP molecules shows that a high precision can be achieved in the calculation of nuclear quadrupole coupling constants of 209Bi by combining advanced ab initio methods with the X2C approach. For seventeen iron compounds, the Mössbauer nuclear quadrupole splittings (NQS) of 57Fe calculated using a double-hybrid functional method are in very good agreement with the experimental values. It is shown that, for strongly correlated molecules, the double-hybrid functionals are much more accurate than the commonly used hybrid functionals. The computer program developed in this study furnishes a useful utility for obtaining EFGs and related nuclear properties with high accuracy.

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.

Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C₅Me₅)₂Y(μ-S)₂Mo(μ-S)₂Y(C₅Me₅)₂]^-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C₆Cl₅)₄]^- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.

Stable copernicium hexafluoride (CnF6) with an oxidation state of VI

As the heaviest group 12 element known currently, copernicium (Cn) often presents the oxidation states of I+, II+, and rarely IV+ as in its homologue mercury. In this work we systematically studied the stability of some oxides, fluorides, and oxyfluorides of Cn by two-component relativistic calculations and found that the CnF₆ molecule with an oxidation state of VI+ has an extraordinary stability. CnF₆ may decompose into CnF₄ by conquering an energy barrier of about 34 kcal mol^-1 without markedly releasing heat. Our results indicate that CnF₆ may exist under some special conditions.

Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds

We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi₁₂]^3-. [U@Bi₁₂]^3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi₁₂ torus similar to the all-metal aromatic heavy-element cluster [Th@Bi₁₂]^4-.

Theoretical Prediction and Interpretation of 237Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of ²³⁷Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the ²³⁷Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin-orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

NMR Indirect Spin-Spin Coupling Constants in a Modern Quasi-Relativistic Density Functional Framework

A quasi-relativistic implementation of NMR indirect spin-spin coupling constants is presented. The exact two-component (X2C) Hamiltonian and its diagonal local approximation to the unitary decoupling transformation (DLU) are utilized together with the (modified) screened nuclear spin-orbit approach. In a restricted kinetic balance, the finite nucleus model is available for both the scalar and vector potentials. The implementation supports density functionals up to the fourth rung of Jacob's ladder, i.e., (range-separated) hybrid and local hybrid functionals based on a seminumerical ansatz. We assess the quality of our quasi-relativistic X2C approach by comparison with "fully" relativistic four-component results for small main-group molecules and alkynyl compounds. The mean absolute error introduced by the DLU scheme is less than 0.05 × 10¹⁹ T J^-2 of the reduced coupling constant for the small main-group molecules and 0.5 Hz for the alkynyl compounds. Thus, the error is significantly smaller than finite nucleus size effects for heavy elements. The basis set convergence and the impact of different density functional approximations are further studied. We propose a simple scheme to develop segmented-contracted relativistic all-electron basis sets for NMR spin-spin couplings. Our implementation allows us to perform calculations of extended molecules with reasonable computational effort, which is illustrated for the ¹J(¹¹⁹Sn, ³¹P) coupling constant of a low-valent tin phosphinidenide complex. The corresponding results are in good agreement with the experimental findings.