1
|
Zapata Escobar AD, Maldonado AF, Aucar GA. The LRESC-Loc Model to Analyze Magnetic Shieldings with Localized Molecular Orbitals. J Phys Chem A 2022; 126:9519-9534. [PMID: 36512732 DOI: 10.1021/acs.jpca.2c05604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The leading electronic mechanisms of relativistic effects in the NMR magnetic shieldings of heavy-atom (HA) containing molecules are well described by the linear response with elimination of small components model (LRESC). We show here first results from a new version of the LRESC model written in terms of localized molecular orbitals (LMOs) which is coined as LRESC-Loc. Those LMOs resemble "chemist's orbitals", representing lone-pairs, atomic cores, and bonds. The whole set of relativistic effects are expressed in terms of non-ligand-dependent and ligand-dependent contributions. We show the electronic origin of trends and behavior of different mechanisms in molecular systems which contain heavy elements that belong to any of the IB to VIIA groups of the periodic table. The SO mechanism has a well-defined dependence with the LPs (LPσ and LPπ) when the HAs have them, but the non-SO mechanisms mostly depend on other LMOs. In addition we propose here that the SO mechanism can be used to characterize interactions involving LPs and the non-SO mechanisms to characterize covalent and close-shell interactions. All our main results are in accord with previous findings, though we are now able to analyze them in a different manner.
Collapse
Affiliation(s)
- Andy D Zapata Escobar
- Physics Department, Natural and Exact Science Faculty, Northeastern University of Argentina, Corrientes, W3404AAS, Argentina.,Institute for Modeling and Innovative Technology, IMIT (CONICET-UNNE), Corrientes, W3404AAS, Argentina
| | - Alejandro F Maldonado
- Institute for Modeling and Innovative Technology, IMIT (CONICET-UNNE), Corrientes, W3404AAS, Argentina
| | - Gustavo A Aucar
- Physics Department, Natural and Exact Science Faculty, Northeastern University of Argentina, Corrientes, W3404AAS, Argentina.,Institute for Modeling and Innovative Technology, IMIT (CONICET-UNNE), Corrientes, W3404AAS, Argentina
| |
Collapse
|
2
|
Melo JI, Maldonado AF, Aucar GA. Performance of the LRESC Model on top of DFT Functionals for Relativistic NMR Shielding Calculations. J Chem Inf Model 2020; 60:722-730. [PMID: 31877038 DOI: 10.1021/acs.jcim.9b00912] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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 SnX4 (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 (SnI4).
Collapse
Affiliation(s)
- Juan I Melo
- Departamento de Física, Facultad de Ciencias Exactas y Naturales , Universidad de Buenos Aires and IFIBA CONICET , Buenos Aires 1428 , Argentina
| | - Alejandro F Maldonado
- Physics Department, Natural and Exact Science Faculty, Northeastern University of Argentina and Institute of Modelling and Innovation on Technology , IMIT CONICET-UNNE , Corrientes W3404AAS , Argentina
| | - Gustavo A Aucar
- Physics Department, Natural and Exact Science Faculty, Northeastern University of Argentina and Institute of Modelling and Innovation on Technology , IMIT CONICET-UNNE , Corrientes W3404AAS , Argentina
| |
Collapse
|
3
|
Melo JI, Maldonado AF, Aucar GA. Relativistic effects on nuclear magnetic shieldings of CHnX4−n and CHXYZ (X, Y, Z = H, F, Cl, Br, I). J Chem Phys 2012; 137:214319. [DOI: 10.1063/1.4768470] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
4
|
Arcisauskaite V, Melo JI, Hemmingsen L, Sauer SPA. Nuclear magnetic resonance shielding constants and chemical shifts in linear 199Hg compounds: A comparison of three relativistic computational methods. J Chem Phys 2011; 135:044306. [DOI: 10.1063/1.3608153] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Melo JI, Maldonado A, Aucar GA. Relativistic effects on the shielding of SnH2XY and PbH2XY (X, Y = F, Cl, Br and I) heavy atom–containing molecules. Theor Chem Acc 2011. [DOI: 10.1007/s00214-010-0886-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
6
|
Roukala J, Maldonado AF, Vaara J, Aucar GA, Lantto P. Relativistic effects on group-12 metal nuclear shieldings. Phys Chem Chem Phys 2011; 13:21016-25. [DOI: 10.1039/c1cp22043h] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Iliaš M, Saue T, Enevoldsen T, Jensen HJA. Gauge origin independent calculations of nuclear magnetic shieldings in relativistic four-component theory. J Chem Phys 2009; 131:124119. [DOI: 10.1063/1.3240198] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
|
8
|
Zaccari D, Melo JI, Ruiz de Azúa MC, Giribet CG. Relativistic two-component geometric approximation of the electron-positron contribution to magnetic properties in terms of Breit-Pauli spinors. J Chem Phys 2009; 130:084102. [PMID: 19256592 DOI: 10.1063/1.3063639] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
An alternative approach for the calculation of the electron-positron (e-p) contribution to magnetic properties based on two-component Breit-Pauli spinors is presented. In it, the elimination of the small component scheme is applied to the inverse propagator matrix of e-p pairs. The effect of the positronic manifold is expressed as an operator acting on Breit-Pauli spinors. The operator form thus obtained sums up the relativistic correction as a geometric series and as a result a totally different behavior in the vicinity of a nucleus is obtained as compared to the one of the linear response approximation. This feature has deep influence in numerical values of the e-p contribution to the nuclear magnetic shielding of heavy atoms. Numerical calculations carried out for Kr, Xe, and I show that with this approach, the e-p contributions to this property are in good agreement with those of four-component methods.
Collapse
Affiliation(s)
- Daniel Zaccari
- Departamento de Fisica, Facultad de Ciencias Exactas, Fisicoquimicas y Naturales, Universidad Nacional de Rio Cuarto, Ruta 36, km 601, 5800 Rio Cuarto, Cordoba, Argentina
| | | | | | | |
Collapse
|
9
|
Casabianca LB, de Dios AC. Ab initiocalculations of NMR chemical shifts. J Chem Phys 2008; 128:052201. [DOI: 10.1063/1.2816784] [Citation(s) in RCA: 100] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|
10
|
Vaara J. Theory and computation of nuclear magnetic resonance parameters. Phys Chem Chem Phys 2007; 9:5399-418. [PMID: 17925967 DOI: 10.1039/b706135h] [Citation(s) in RCA: 192] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The art of quantum chemical electronic structure calculation has over the last 15 years reached a point where systematic computational studies of magnetic response properties have become a routine procedure for molecular systems. One of their most prominent areas of application are the spectral parameters of nuclear magnetic resonance (NMR) spectroscopy, due to the immense importance of this experimental method in many scientific disciplines. This article attempts to give an overview on the theory and state-of-the-art of the practical computations in the field, in terms of the size of systems that can be treated, the accuracy that can be expected, and the various factors that would influence the agreement of even the most accurate imaginable electronic structure calculation with experiment. These factors include relativistic effects, thermal effects, as well as solvation/environmental influences, where my group has been active. The dependence of the NMR spectra on external magnetic and optical fields is also briefly touched on.
Collapse
Affiliation(s)
- Juha Vaara
- Laboratory of Physical Chemistry, Department of Chemistry, FIN-00014, University of Helsinki, Finland.
| |
Collapse
|
11
|
Lantto P, Romero RH, Gómez SS, Aucar GA, Vaara J. Relativistic heavy-atom effects on heavy-atom nuclear shieldings. J Chem Phys 2006; 125:184113. [PMID: 17115744 DOI: 10.1063/1.2378737] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The principal relativistic heavy-atom effects on the nuclear magnetic resonance (NMR) shielding tensor of the heavy atom itself (HAHA effects) are calculated using ab initio methods at the level of the Breit-Pauli Hamiltonian. This is the first systematic study of the main HAHA effects on nuclear shielding and chemical shift by perturbational relativistic approach. The dependence of the HAHA effects on the chemical environment of the heavy atom is investigated for the closed-shell X(2+), X(4+), XH(2), and XH(3) (-) (X=Si-Pb) as well as X(3+), XH(3), and XF(3) (X=P-Bi) systems. Fully relativistic Dirac-Hartree-Fock calculations are carried out for comparison. It is necessary in the Breit-Pauli approach to include the second-order magnetic-field-dependent spin-orbit (SO) shielding contribution as it is the larger SO term in XH(3) (-), XH(3), and XF(3), and is equally large in XH(2) as the conventional, third-order field-independent spin-orbit contribution. Considering the chemical shift, the third-order SO mechanism contributes two-thirds of the difference of approximately 1500 ppm between BiH(3) and BiF(3). The second-order SO mechanism and the numerically largest relativistic effect, which arises from the cross-term contribution of the Fermi contact hyperfine interaction and the relativistically modified spin-Zeeman interaction (FC/SZ-KE), are isotropic and practically independent of electron correlation effects as well as the chemical environment of the heavy atom. The third-order SO terms depend on these factors and contribute both to heavy-atom shielding anisotropy and NMR chemical shifts. While a qualitative picture of heavy-atom chemical shifts is already obtained at the nonrelativistic level of theory, reliable shifts may be expected after including the third-order SO contributions only, especially when calculations are carried out at correlated level. The FC/SZ-KE contribution to shielding is almost completely produced in the s orbitals of the heavy atom, with values diminishing with the principal quantum number. The relative contributions converge to universal fractions for the core and subvalence ns shells. The valence shell contribution is negligible, which explains the HAHA characteristics of the FC/SZ-KE term. Although the nonrelativistic theory gives correct chemical shift trends in present systems, the third-order SO-I terms are necessary for more reliable predictions. All of the presently considered relativistic corrections provide significant HAHA contributions to absolute shielding in heavy atoms.
Collapse
Affiliation(s)
- Perttu Lantto
- NMR Research Group, Department of Physical Sciences, University of Oulu, P.O. Box 3000, FIN-90014 Oulu, Finland.
| | | | | | | | | |
Collapse
|
12
|
Roura PG, Melo JI, Ruiz de Azúa MC, Giribet CG. Mean field linear response within the elimination of the small component formalism to evaluate relativistic effects on magnetic properties. J Chem Phys 2006; 125:64107. [PMID: 16942273 DOI: 10.1063/1.2244572] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The linear response within the elimination of the small component formalism is aimed at obtaining the leading order relativistic corrections to magnetic molecular properties in the context of the elimination of the small component approximation. In the present work we extend the method in order to include two-body effects in the form of a mean field one-body operator. To this end we consider the four-component Dirac-Hartree-Fock operator as the starting point in the evaluation of the second order relativistic expression of magnetic properties. The approach thus obtained is the fully consistent leading order approximation of the random phase approximation four-component formalism. The mean field effect on the relativistic corrections to both the diamagnetic and paramagnetic terms of magnetic properties taking into account both the Coulomb and Breit two-body interactions is considered.
Collapse
Affiliation(s)
- P G Roura
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pabellon I, 1428 Buenos Aires, Argentina
| | | | | | | |
Collapse
|
13
|
|
14
|
Manninen P, Ruud K, Lantto P, Vaara J. Leading-order relativistic effects on nuclear magnetic resonance shielding tensors. J Chem Phys 2005; 122:114107. [PMID: 15836201 DOI: 10.1063/1.1861872] [Citation(s) in RCA: 105] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.
Collapse
Affiliation(s)
- Pekka Manninen
- Laboratory of Physical Chemistry, Department of Chemistry, University of Helsinki, Finland.
| | | | | | | |
Collapse
|
15
|
Gómez SS, Melo JI, Romero RH, Aucar GA, de Azúa MR. Relativistic corrections to the diamagnetic term of the nuclear magnetic shielding: Analysis of contributions from localized orbitals. J Chem Phys 2005; 122:064103. [PMID: 15740363 DOI: 10.1063/1.1845391] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have calculated the relativistic corrections to the diamagnetic term of the nuclear magnetic shielding constants for a series of molecules containing heavy atoms. An analysis of the contributions from localized orbitals is performed. We establish quantitatively the relative importance of inner core and valence shell molecular orbitals in each correcting term. Contributions from the latter are much less important than those from the former. The calculated values of the correction sigma(L-PSO), first derived within the linear response elimination of small component formalism, show a power-law dependence on the nuclear charge approximately Z(3.5), in contrast with the approximately Z(3.1) behavior of the mass-velocity external-field correction to the paramagnetic term previously reported.
Collapse
Affiliation(s)
- Sergio S Gómez
- Physics Department, FCENA, Northeastern University, Avda. Libertad 5500, W 3404 AAS Corrientes, Argentina
| | | | | | | | | |
Collapse
|
16
|
Melo JI, Ruiz de Azua MC, Giribet CG, Aucar GA, Provasi PF. Relativistic effects on nuclear magnetic shielding constants in HX and CH3X (X=Br,I) based on the linear response within the elimination of small component approach. J Chem Phys 2004; 121:6798-808. [PMID: 15473737 DOI: 10.1063/1.1787495] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Numerical calculations of relativistic effects on nuclear magnetic shielding constants sigma corresponding to all one-body operators obtained within a formalism developed in previous work were carried out. In this formalism, the elimination of small component scheme is applied to evaluate all quantities entering a four-component RSPT(2) expression of magnetic molecular properties. HX and CH3X (X=Br,I) were taken as model compounds. Calculations were carried out at the Hartree-Fock level for first-order quantities, and at the random-phase approximation (RPA) level for second- and third-order ones. It was found that values of sigma(X) are largely affected by several relativistic corrections not previously considered in the bibliography. sigma Values of the H nucleus are in close agreement with four-component RPA ones. Overall relativistic effects on the shift of sigma(X) from HX to CH3X are smaller than the nonrelativistic shifts.
Collapse
Affiliation(s)
- J I Melo
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Cdad. Universitaria, Pab. 1, 1428 Buenos Aires, Argentina
| | | | | | | | | |
Collapse
|