镱硫属化合物的密度泛函理论研究

用密度泛函理论(DFT)研究镱硫属化合物的电子结构和性质,通过与实验比较考察了现有的几种近似密度泛函公式对镧系元素化合物的适用程度和相对论效应的影响.结果表明,用DFT计算的YbO键长对实验值的偏差约为0.002nm;但得到的键能即使在考虑梯度校正和相对论效应之后,仍比实验值高,在定域密度近似基础上引入交换梯度校正使键能计算值减小,其中PW86x使键能计算值减小稍多些,结果更接近实验值;相关梯度校正使键能计算值升高.相对论效应使键长缩短0.004～0.006nm,键能减小约0.5eV.计算结果的分析表明,Yb的5d轨道和配体的np轨道间形成σ键和π键.在所研究的分子体系中,配体原子从O到Te、Yb原子的5d轨道布后数依次减少,同键能减弱的顺序一致.相对论效应使键能减小的主要原因是在成键过程中发生了Yb的6s电子向5d轨道的转移,而相对论效应使该过程能量增加.偶极矩和电荷分布的计算表明,Yb-L键以共价性为主,相对论效应使共价性成份增加.     

Calculations of nuclear quadrupole coupling in noble gas-noble metal fluorides: interplay of relativistic and electron correlation effects

The nuclear quadrupole coupling constants (NQCCs) of noble gas and noble metal nuclei in the recently found noble gas-noble metal fluorides (NgMF, where Ng=Ar,Kr,Xe and M=Cu,Ag,Au) are obtained theoretically by high-level ab initio calculations, where both relativistic and electron correlation effects are included, and compared to experimental results. Fully relativistic four-component Dirac-Hartree-Fock (DHF) calculations are carried out at the basis set limit for electric field gradient that couples with the electric quadrupole moment of the nucleus, and uncorrelated relativistic effects are extracted by comparing DHF results to nonrelativistic (NR) HF calculations. Electron correlation effects are investigated both at fully relativistic second-order Moller-Plesset (DMP2) and at NR MP2 levels of theory, as well as at the NR coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level. The validity of the approximation where relativistic effects, on the one hand, and nonrelativistically obtained correlation effects, on the other hand, are evaluated separately and assumed to be additive, is investigated by comparison with the DMP2 results. Inclusion of relativistic effects is shown to be necessary for obtaining the correct NQCC trends as the nucleus of interest and/or its neighbors become heavier. Electron correlation treatment is needed for approaching quantitative agreement with the experimental NQCCs. The assumption of additive electron correlation and relativistic effects, corresponding to the NR correlation treatment added on top of relativistic DHF data, gives qualitatively correct noble gas NQCCs. For noble metal NQCCs, correlation treatment at the relativistic level of theory is mandatory for reaching agreement with experimental results. Current work also confirms the experimental trends of NQCCs, which have been taken as an indication of nearly covalent interaction between noble gas and noble metal in the heaviest present systems, especially in XeAuF.     

Theoretical investigation of the apparently irregular behavior of pt-pt nuclear spin-spin coupling constants

One-bond Pt-Pt nuclear spin-spin coupling constants J(Pt-Pt) for closely related dinuclear Pt complexes can differ by an order of magnitude without any obvious correlation with Pt-Pt distances. As representative examples, the spin-spin couplings of the dinuclear Pt(I) complexes [Pt(2)(CO)(6)](2+) (1) and [Pt(2)(CO)(2)Cl(4)](2-) (2) have been computationally studied with a recently developed relativistic density functional method. The experimental values are (1)J((195)Pt-(195)Pt) = 5250 Hz for 2 but 551 Hz for 1. Many other examples are known in the literature. The experimental trends are well reproduced by the computations and can be explained based on the nature of the ligands that are coordinated to the Pt-Pt fragment. The difference for J(Pt-Pt) of an order of magnitude is caused by a sensitive interplay between the influence of different ligands on the Pt-Pt bond, and relativistic effects on metal-metal and metal-ligand bonds as well as on "atomic orbital contributions" to the nuclear spin-spin coupling constants. The results can be intuitively rationalized with the help of a simple qualitative molecular orbital diagram.     

水溶液中碳酸铀酰化合物的电子结构

应用相对论密度泛函理论系统研究了水溶液中非水合化和水合化碳酸铀酰化合物Cn/m(其中n和m分别为结构中碳酸配体和水配体的个数)的结构.溶剂效应采用类导体屏蔽模型(COSMO),并采用零级规整近似(ZORA)方法考虑标量相对论效应和旋-轨耦合相对论效应.电子跃迁采用包含旋-轨耦合相对论效应的含时密度泛函理论并在相关交换势中采用轨道势能统计平均(SAOP)做近似计算.结果表明碳酸配体对配合物结构和电子跃迁有很大的影响.C3/0配合物的稳定性可归于5f轨道参与了高占据轨道的成键作用.增加碳酸盐配体导致最大波长的蓝移,并在近可见光区域出现高强度的吸收.     

Quantum relativistic investigation about the coordination and bonding effects of different ligands on uranyl complexes

The coordination and bonding effects of equatorial ligands such as fluoride (F⁻), chloride (Cl⁻), cyanide (CN⁻), isocyanide (NC⁻), and carbonate (CO₃⁻²) on uranyl dication (UO₂²⁺) has been studied using relativistic density functional theory. The ZORA Hamiltonian was applied for the inclusion of relativistic effects taking into account all the electrons for the optimization and the explicit inclusion of spin–orbit coupling effects. Geometry optimizations including the counterions and frequencies analysis were carried out with PW91 and PBE functional. Solvents effects were considered by using the conductor like screening model (COSMO) for water and acetonitrile. The Time-Dependent Density Functional Theory (TDDFT) was used to calculate the excitation energies with GGA SAOP functional and the electronic transitions were analyzed using double group irreducible representations. The theoretical results are in a good agreement with experimental IR, Raman and EXAFS spectra and previous theoretical results. New information about the effect of different (donor and acceptors) ligands on the bonding of uranyl ion and on the electronic transitions involved in these complexes is provided with a possible impact on the understanding of the uranyl coordination chemistry.     

MX~4(M=Ti,Zr,Hf;X=Cl,Br)的电子结构研究

本文利用单电子非相对论Hartree-Fock-Slater和完全相对论Dirac-Fock-Slater两种离散变分局域密度泛函方法(DV-Xα), 对MX~4(M=Ti,Zr,Hf;X=Cl,Br)的电子基态和相应于低能UV光谱的激发态进行了计算, 结果与实验符合得较好。用Mulliken布居分析方法研究了分子的共价性质, 发现除HfBr~4外,相对论效应对金属与配体之间的键级影响很小。     

Comparative studies of quasi-relativistic density functional methods for the description of lanthanide and actinide complexes

We present a comparative Density Functional Theory (DFT) study based on two different implementations of relativistic effects within the Kohn-Sham (KS) approach, to describe the metal-ligand interaction in I(3)M-L complexes (L = NH(3), NCCH(3), CO and M = La, Nd, U). In the first model, the scalar corrections were included by a quasi-relativistic approach (QR) via the so-called ZORA or Pauli Hamiltonians, while in the second, these effects are taken into account in a quasi-Relativistic Effective Core Potential (RECP). These relativistic approaches were used in conjunction with various gradient corrected (GGA) or hybrid (SCH) functionals. The structural parameters obtained from geometry optimizations have been compared to experimental structural trends, and rationalized by a KS orbital analysis. Both approaches provide similar results for mainly ionic metal-ligand bonds (e.g., for the sigma-donor ligand L = NH(3)). For the pi-acceptor ligands (NCCH(3), CO), the QR approach is in agreement with experimental trends and consistent with the presence of a backbonding interaction between U(III) and the neutral ligand, which does not exist in the lanthanide homologues. The GGA/RECP methods also reproduce this phenomenon, while the SCH/RECP scheme fails to describe this interaction. The role of the RECP, of its size, and of additional polarization functions has also been examined. Finally, the failure of the SCH/RECP approach was interpreted as a consequence of a bad estimation of frontier orbital energy levels in the uranium and ligand species.     

Organolanthanide mediated catalytic cycles: a computational perspective

In this perspective the contribution of recent theoretical studies to our understanding of lanthanide (Ln) catalysis is explored. In general, the results of computational studies have proven consistent with available experimental evidence. Considerable success has been obtained in elucidating the mechanisms for C-H bond activation (sigma-bond metathesis in particular) and the addition of C-X bonds across an unsaturated functionality (and the hydroamination of alkenes in particular). Ln catalysts are computationally challenging because relativistic effects are important, and large ligands are required to restrict high coordination numbers, in addition to limiting facile redistribution processes. Thus, key technical issues relating to the computational investigation of organolanthanide complexes are discussed. Increasing computational resources have seen studies expand from the optimisation of simple molecules to the study of catalytic cycles where the Ln is coordinated by larger and more complex ligands. The ability of theoretical studies to complement experimental developments by supplying a deeper understanding of the mechanistic process is reviewed with emphasis on the elucidation of transition state structures, intermediates, spectator ligand coordination, and negative entropy steps. Recent computational investigations of the catalytic cycle for Ln mediated hydroamination are a focus, as these have provided substantial and detailed rationalisations for the regio- and stereo-selectivity of inter- and intra-molecular hydroamination. Examination of transition state geometries and electronic structure appears to offer insights that could be used to facilitate the rational design of new Ln-based catalysts.     

A density functional study of uranyl monocarboxylates

We studied uranium(VI) monocarboxylate complexes by a relativistic density functional method using simple carboxylic acids as ligands, i.e. [UO2(OOCR)]+ (R = H, CH3, CH2CH3). These complexes exist in aqueous solution and, for R = CH3 and CH2CH3, may also be considered as models of uranyl complexated by humic substances. We investigated mono- and bidentate coordination modes. Short-range solvent effects were accounted for explicitly via aqua ligands of the first hydration shell and long-range electrostatic interactions were described via a polarizable continuum model. The calculated results for the uranyl U=O bond, the bond to aqua ligands, and the averaged uranium distances to equatorial oxygen atoms, U-Oeq, agreed quite well with EXAFS-derived interatomic distances. However, the uranyl-carboxylate bond was calculated to be notably shorter than the experimentally determined value. Experimental differences between mono- and bidentate coordination, obtained mainly from crystal structures, were qualitatively reproduced for the U-C distance but not for the average bond length, U-Oeq. We discuss these discrepancies between calculated and experimental results in some detail and suggest changes in the coordination number rather than variations of the coordination geometry as the main source of the experimentally observed variation of the U-Oeq distance.     

On the unusual weak intramolecular C...C interactions in Ru3(CO)12: a case of bond path artifacts introduced by the multipole model?

In a recent publication in this journal, an experimental charge density analysis on the triruthenium cluster Ru(3)(CO)(12) showed unusual C...C bond paths linking the axial carbonyl ligands [Gervasio, G.; Marabello, D.; Bianchi, R.; Forni, A. J. Phys. Chem. A 2010, 114, 9368, hereafter GMBF]. These were also observed in one theoretical DFT calculation, and are associated with very low values of ρ(r(b)) and ?(2)ρ(r(b)). Our independent experimental charge density analysis on Ru(3)(CO)(12) is entirely consistent with GMBF and confirms the presence of these apparent weak interactions in the multipole model density. However, we conclusively demonstrate that these unusual C...C bond paths between the axial carbonyl ligands are in fact artifacts arising from the Hansen-Coppens multipole model, which is used to analyze the experimental data. Numerous relativistic and nonrelativistic gas-phase DFT calculations, using very extensive basis sets and with corrections for dispersion effects, uniformly fail to reproduce these intramolecular features in the QTAIM topology of the electron density. Moreover, multipole fitting of theoretical static structure factors computed from these quantum electron densities results in the reappearance of the C...C bond paths between the axial carbonyl ligands in the derived molecular graphs. On the other hand, using the experimental structure factors to generate "experimental" X-ray constrained DFT wave functions once again yields molecular graphs which do not show these secondary C...C bond paths. The evidence therefore strongly implicates the multipole model as the source of these spurious features and in turn suggests that great caution should be applied in the interpretation of bond paths where the values of ρ(r(b)) and ?(2)ρ(r(b)) are very low.     

Simulations of 129Xe NMR chemical shift of atomic xenon dissolved in liquid benzene

The isotropic ¹²⁹Xe NMR chemical shift of atomic Xe dissolved in liquid benzene was simulated by combining classical molecular dynamics and quantum chemical calculations of ¹²⁹Xe nuclear magnetic shielding. Snapshots from the molecular dynamics trajectory of xenon atom in a periodic box of benzene molecules were used for the quantum chemical calculations of isotropic ¹²⁹Xe chemical shift using nonrelativistic density functional theory as well as relativistic Breit?CPauli perturbation corrections. Thus, the correlation and relativistic effects as well as the temperature and dynamics effects could be included in the calculations. Theoretical results are in a very good agreement with the experimental data. The most of the experimentally observed isotropic ¹²⁹Xe shift was recovered in the nonrelativistic dynamical region, while the relativistic effects explain of about 8% of the total ¹²⁹Xe chemical shift.     

Solvent effects on heavy atom nuclear spin-spin coupling constants: a theoretical study of Hg-C and Pt-P couplings.

The computation of indirect nuclear spin-spin coupling constants, based on the relativistic two-component zeroth order regular approximate Hamiltonian, has been recently implemented by us into the Amsterdam Density Functional program. Applications of the code for the calculation of one-bond metal-ligand couplings of coordinatively unsaturated compounds containing (195)Pt and (199)Hg, including spin-orbit coupling or coordination effects by solvent molecules, show that relativistic density functional calculations are able to reproduce the experimental findings with good accuracy for the systems under investigation. Spin-orbit effects are rather small for these cases, while coordination of the heavy atoms by solvent molecules has a great impact on the calculated couplings. Experimental trends for different solvents are reproduced. An orbital-based analysis of the solvent effect is presented. The scalar relativistic increase of the coupling constants is of the same order of magnitude as the nonrelativistically obtained values, making a relativistic treatment essential for obtaining quantitatively correct results. Solvent effects can be of similar importance.     

103Rh NMR chemical shifts in organometallic complexes: a combined experimental and density functional study

Experimental 103Rh NMR chemical shifts of mono- and binuclear rhodium(I) complexes containing s- or as-hydroindacenide and indacenediide bridging ligands with different ancillary ligands (1,5-cyclooctadiene, ethylene, carbonyl) are presented. A protocol, based on density functional theory calculations, was established to determine 103Rh NMR shielding constants in order to rationalise the effects of electronic and structural variations on the spectroscopic signal, and to gain insight into the efficiency of this computational method when applied to organometallic systems. Scalar and spin-orbit relativistic effects based on the ZORA (zeroth order regular approximation) level have been taken into account and discussed. A good agreement was found for model compounds over a wide range of chemical shifts of rhodium (approximately 10,000 ppm). This allowed us to discuss the experimental and calculated delta(103Rh) in larger complexes and to relate it to their electronic structure.     

Relativistic effects in gas-phase ion chemistry: an experimentalist's view

Gas-phase experiments provide information which, in conjunction with results from electronic structure calculations, help to unravel the critical role relativistic effects play in many areas of transition-metal chemistry. Examples include the thermochemical data of gold halides in different oxidation states, the fascinating structural properties of gold(I) complexes, the dramatic effects of ligands on the ionization energy of gold, or the binding in cationic metal-carbene complexes. Furthermore, in the context of methane functionalization, special emphasis is paid to the chemistry of cationic metal-carbene complexes, and at uncovering the mechanistic details of important carbon-heteroatom coupling reactions. It is the interplay of conducting experiments of "isolated" molecules under well-defined conditions with reliable electronic structure calculations that has considerably improved our understanding of the role relativistic effects play in the context of transition-metal chemistry, catalysis, and beyond.     

Ionization potentials and electron affinities of Cu,Ag, and Au: Electron correlation and relativistic effects

The electron correlation and relativistic effects on ionization potentials and electron affinities of Cu, Ag, and Au are investigated in the framework of the coupled cluster method and different 1-component approximations to the relativistic Dirac-Coulomb Hamiltonian. The first-order perturbation approach based on the massvelocity and Darwin terms is found to be sufficiently accurate for Cu and Ag while it fails for Au. The spin-averaged Douglas-Kroll no-pair method gives excellent results for the studied atomic properties. The ionization potentials obtained within this method and the coupled cluster scheme for the electron correlation effects are 7.733(7.735) eV for Cu, 7.461(7.575) eV for Ag, and 9.123(9.225) eV for Au (experimental values given in parentheses). The calculated (experimental) electron affinity results for Cu, Ag, and Au are 1.236(1.226), 1.254(1.303), and 2.229(2.309) eV, respectively. There is a marked relativistic effect on both the ionization potential and electron affinity of Ag which sharply increases for Au while Cu exhibits only a little relativistic character. A similar pattern of relativistic effects is also observed for electric dipole polarizabilities of the coinage metal atoms and their ions. The coupled cluster dipole polarizabilities of the coinage metal atoms calculated in this article in the Douglas-Kroll no-pair formalism (Cu: 46.50 au; Ag: 52.46 au; Au: 36.06 au) are compared with our earlier data for their singly positive and singly negative ions. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 557–565, 1997     

Theoretical predictions of nuclear magnetic resonance parameters in a novel organo-xenon species: chemical shifts and nuclear quadrupole couplings in HXeCCH

We calibrate the methodology for the calculation of nuclear magnetic resonance (NMR) properties in novel organo-xenon compounds. The available state-of-the-art quantum-chemical approaches are combined and applied to the HXeCCH molecule as the model system. The studied properties are (129)Xe, (1)H, and (13)C chemical shifts and shielding anisotropies, as well as (131)Xe and (2)H nuclear quadrupole coupling constants. The aim is to obtain, as accurately as currently possible, converged results with respect to the basis set, electron correlation, and relativistic effects, including the coupling of relativity and correlation. This is done, on one hand, by nonrelativistic correlated ab initio calculations up to the CCSD(T) level and, on the other hand, for chemical shifts and shielding anisotropies by the leading-order relativistic Breit-Pauli perturbation theory (BPPT) with correlated ab initio and density-functional theory (DFT) reference states. BPPT at the uncorrelated Hartree-Fock level as well as the corresponding fully relativistic Dirac-Hartree-Fock method are found to be inapplicable due to a dramatic overestimation of relativistic effects, implying the influence of triplet instability in this multiply bonded system. In contrast, the fully relativistic second-order Moller-Plesset perturbation theory method can be applied for the quadrupole coupling, which is a ground-state electric property. The performance of DFT with various exchange-correlation functionals is found to be inadequate for the nonrelativistic shifts and shielding anisotropies as compared to the CCSD(T) results. The relativistic BPPT corrections to these quantities can, however, be reasonably predicted by DFT, due to the improved triplet excitation spectrum as compared to the Hartree-Fock method, as well as error cancellation within the five main BPPT contributions. We establish three computationally feasible models with characteristic error margins for future calculations of larger organo-xenon compounds to guide forthcoming experimental NMR efforts. The predicted (129)Xe chemical shift in HXeCCH is in a novel range for this nucleus, between weakly bonded or solvated atomic xenon and xenon in the hitherto characterized molecules.     

NMR shielding constants in hydrogen molecule isotopomers

We analyze the NMR shielding constants in three isotopomers of the hydrogen molecule: H₂, HD and D₂. The results obtained within the Born?COppenheimer approximation using the coupled-cluster singles-and-doubles model are very close to the previous theoretical values. In particular, the isotope shifts computed using significantly larger basis sets agree with the earlier literature results, confirming the disagreement of these calculations with the available experimental data. To examine the accuracy of the computed isotope shifts, we analyze in addition the relativistic corrections and estimate the role of the adiabatic and nonadiabatic effects. The relativistic corrections appear to be negligible; on the other hand, the changes in the shielding constants due to the adiabatic and nonadiabatic effects may account for the discrepancies between the computed and experimental isotope shifts.     

Electronic spectrum of UO2(2+) and [UO2Cl4]2- calculated with time-dependent density functional theory

The electronic spectra of UO(2) (2+) and [UO(2)Cl(4)](2-) are calculated with a recently proposed relativistic time-dependent density functional theory method based on the two-component zeroth-order regular approximation for the inclusion of spin-orbit coupling and a noncollinear exchange-correlation functional. All excitations out of the bonding sigma(u) (+) orbital into the nonbonding delta(u) or phi(u) orbitals for UO(2) (2+) and the corresponding excitations for [UO(2)Cl(4)](2-) are considered. Scalar relativistic vertical excitation energies are compared to values from previous calculations with the CASPT2 method. Two-component adiabatic excitation energies, U-O equilibrium distances, and symmetric stretching frequencies are compared to CASPT2 and combined configuration-interaction and spin-orbit coupling results, as well as to experimental data. The composition of the excited states in terms of the spin-orbit free states is analyzed. The results point to a significant effect of the chlorine ligands on the electronic spectrum, thereby confirming the CASPT2 results: The excitation energies are shifted and a different luminescent state is found.     

DFT studies of uranyl acetate,carbonate, and malonate,complexes in solution

The aim of this work is to demonstrate that theoretical chemistry can be used as a complementary tool in determining geometric parameters of a number of uranyl complexes in solution, which are not observable by experimental methods. In addition, we propose plausible structures with partial geometric data from experimental results. A gradient corrected DFT methodology with relativistic effects is used employing a COSMO solvation model. The theoretical calculations show good agreement with experimental X-ray and EXAFS data for the triacetato-dioxo-uranium(VI) and tricarbonato-dioxo-uranium(VI) complexes and are used to assign possible geometries for dicalcium-tricarbonato-dioxo-uranium(VI) and malonato-dioxo-uranium(VI) complexes. The results of this exercise indicate that carbonate bonding in these complexes is mainly bidentate and that hydroxo bridging plays a critical role in the stabilization of the polynuclear uranyl complexes.