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.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

A comparative relativistic DFT and ab initio study on the structure and thermodynamics of the oxofluorides of uranium(IV), (V) and (VI)

All the possible uranium(VI, V, IV) oxides, fluorides and oxofluorides were studied theoretically by using density functional theory (DFT) in the generalised gradient approximation (GGA), and three different relativistic methods (all-electron scalar four component Dyall RESC method (AE), relativistic small-core ECPs, and zeroth order regular approximation ZORA). In order to test different correlation methods, for the two former relativistic methods hybrid DFT, and, for the AE method, MP2 molecular orbital calculations were performed as well. Single-point AE-CCSD(T) energies were calculated on MP2 geometries as well. Energies of the uranium(VI) and (V) oxofluorides dissociation, uranium(VI) fluoride hydrolysis and oxofluoride disproportionation were calculated and compared against the available experimental thermochemical data. AE-CCSD(T) energies were the closest to the experiment. For GGA DFT methods, all the relativistic methods used yield similar results. For thermochemistry, the best quantitative agreement with the experimental and CCSD(T) values for both U=O and U-F bond strengths was obtained with hybrid DFT methods, provided that a reliable basis set was used. Both the GGA DFT and MP2 MO methods show overbinding of these bonds; moreover, this overbinding was found to be not uniform but strongly dependent on the coordination environment of the uranium atom in each case. U=O vibrational frequencies given by hybrid DFT, however, are systematically overestimated, and are better reproduced by GGA DFT; MP2 values usually fall in-between. Reaction enthalpies, U=O frequencies and complex geometries given by the PBE, MPBE, BPBE, BLYP and OLYP GGA functionals are quite similar, with OLYP performing slightly better than the others but still not as good as hybrid DFT. The geometries of the molecules are found to be influenced by the following factors: the inverse transinfluence (ITI) of the oxygen ligand and, for U(V), and U(IV), the Jahn-Teller distortion.     

Focal point analysis of the singlet-triplet energy gap of octacene and larger acenes

A benchmark theoretical study of the electronic ground state and of the vertical and adiabatic singlet-triplet (ST) excitation energies of n-acenes (C(4n+2)H(2n+4)) ranging from octacene (n = 8) to undecacene (n = 11) is presented. The T1 diagnostics of coupled cluster theory and further energy-based criteria demonstrate that all investigated systems exhibit predominantly a (1)A(g) singlet closed-shell electronic ground state. Singlet-triplet (S(0)-T(1)) energy gaps can therefore be very accurately determined by applying the principle of a focal point analysis (FPA) onto the results of a series of single-point and symmetry-restricted calculations employing correlation consistent cc-pVXZ basis sets (X = D, T, Q, 5) and single-reference methods [HF, MP2, MP3, MP4SDQ, CCSD, and CCSD(T)] of improving quality. According to our best estimates, which amount to a dual extrapolation of energy differences to the level of coupled cluster theory including single, double, and perturbative estimates of connected triple excitations [CCSD(T)] in the limit of an asymptotically complete basis set (cc-pV∞Z), the S(0)-T(1) vertical (adiabatic) excitation energies of these compounds amount to 13.40 (8.21), 10.72 (6.05), 8.05 (3.67), and 7.10 (2.58) kcal/mol, respectively. In line with the absence of Peierls distortions (bond length alternations), extrapolations of results obtained at this level for benzene (n = 1) and all studied n-acenes so far (n = 2-11) indicate a vanishing S(0)-T(1) energy gap, in the limit of an infinitely large polyacene, within an uncertainty of 1.5 kcal/mol (0.06 eV). Lacking experimental values for the S(0)-T(1) energy gaps of n-acenes larger than hexacene, comparison is made with recent optical and electrochemical determinations of the HOMO-LUMO band gap. Further issues such as scalar relativistic, core correlation, and diagonal Born-Oppenheimer corrections (DBOCs) are tentatively examined.     

Calculated spectroscopic and electric properties of the alkali metal-ammonia complexes from Kn-NH3 to Frn-NH3 (n=0,+1)

The newly developed Stuttgart small-core scalar relativistic pseudopotentials for the alkali metals are used to study spectroscopic and electric properties of the heavier alkali metal-ammonia complexes from K(n)-NH(3) to Fr(n)-NH(3) (n=0,+1) at the second-order Moller-Plesset (MP2) and coupled cluster [CCSD(T)] levels of theory. Equilibrium geometries and dissociation energies computed at the MP2 level are in reasonable agreement with their CCSD(T) counterparts, whereas for the dipole polarizabilities MP2 is not performing well overestimating significantly electron correlation effects. The bond distances increase monotonically with increasing mass of the metal atom as relativistic effects are small in these systems. However, the dipole polarizabilities are more sensitive to such effects and we find a decrease in this property from Cs-NH(3) to Fr-NH(3). Combination of CCSD(T) harmonic frequencies and MP2 anharmonic corrections obtained from a perturbative vibrational treatment leads to fundamental frequencies in good agreement with experimental results obtained by Suzer and Andrews [J. Am. Chem. Soc. 109, 300 (1986)]. We also present the results of variational calculations with a three-dimensional vibrational Hamiltonian, making use of CCSD(T) potential energy and electric dipole moment surfaces. Complexation of NH(3) to the metal causes a strong infrared intensification of the symmetric NH(3) stretching mode in the neutral complexes, which is absent in the charged species.     

Evaluation of DFT methods for computing the interaction energies of homomolecular and heteromolecular dimers of monosubstituted benzene

We present density functional theory (DFT) interaction energies for the sandwich and T‐shaped conformers of substituted benzene dimers. The DFT functionals studied include TPSS, HCTH407, B3LYP, and X3LYP. We also include Hartree–Fock (HF) and second‐order Møller–Plesset perturbation theory calculations (MP2), as well as calculations using a new functional, P3LYP, which includes PBE and HF exchange and LYP correlation. Although DFT methods do not explicitly account for the dispersion interactions important in the benzene–dimer interactions, we find that our new method, P3LYP, as well as HCTH407 and TPSS, match MP2 and CCSD(T) calculations much better than the hybrid methods B3LYP and X3LYP methods do. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

Studies of neutral and ionic CuAr and CuKr van der Waals complexes

Ab initio calculations of the interaction potential between Cu (or Cu(+)) and Ar (or Kr) have been carried out. A range of theoretical methods, including Hartree-Fock (HF), Moeller-Plesset perturbation methods to second order (MP2), and single and double excitation coupled cluster methods, with the perturbational effect of triple excitations (CCSD(T)), were employed with relativistic pseudopotential basis sets. The effects of bond functions and diffuse polarization (f, g, h) functions were tested on the calculation of the weak intermolecular interactions. Potential energy curves were obtained for all four complexes by MP2 and CCSD(T) methods. For CuKr, even with the largest basis set used, a binding energy that is only 37% of the measured value was obtained. Possible reasons for the disagreement are discussed.     

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

Geometries,vibrational frequencies,and electron affinities of X2Cl (X=C,Si,Ge) clusters

Ab initio quantum chemical calculations have been performed on X₂Cl^? and X₂Cl (X = C, Si, Ge) clusters. The geometrical structures, vibrational frequencies, electronic properties and dissociation energies are investigated at the Hartree–Fock (HF), Møller–Plesset second‐ and fourth‐order (MP2, MP4), CCSD(T) level with the 6‐311+G(d) basis set. The X₂Cl (X = C, Si, Ge) and X₂Cl^? (X = Si, Ge) take a bent shape obtained at the ground state, while C₂Cl^? has a linear structure. The impact on internal electron transfer between the X₂Cl and the corresponding anional clusters is studied. The three different types of electron affinities (EAs) at the CCSD(T) are reported. The most reliable adiabatic electronic affinities, obtained at the CCSD(T)/cc‐pvqz level of theory, are predicted to be 3.30, 2.62, and 1.98 eV for C₂Cl, Si₂Cl, and Ge₂Cl, respectively. The calculated EAs of C₂Cl and Ge₂Cl are in good agreement with theoretical results reported. The correlation effects and basis sets effects on the geometrical structures and dissociation energies are discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete

The conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero-point vibrational energies are computed for all pertinent molecular systems using SCF theory, second-order perturbation theory, and density functional theory and employing the correlation consistent basis sets cc-pVDZ, cc-pVTZ, and cc-pVQZ. Single-point fourth-order perturbation theory, CCSD, and CCSD(T) calculations employing the cc-pVTZ and the cc-pVQZ basis sets are computed using the MP2/cc-pVTZ and MP2/cc-pVQZ optimized geometries, respectfully, to ascertain the contribution of higher order correlation. Three DFT functionals, B3LYP, wB97XD, and M06-2X, are employed to determine whether they can yield results similar to those obtained at the CCSD(T) level.     

Relativistic calculation of nuclear magnetic shielding using normalized elimination of the small component

The normalized elimination of the small component (NESC) theory, recently proposed by Filatov and Cremer, is extended to include magnetic interactions and applied to the calculation of the nuclear magnetic shielding in HX (X=F, Cl, Br, I) systems. The NESC calculations are performed at the levels of the zeroth-order regular approximation (ZORA) and the second-order regular approximation (SORA). The calculations show that the NESC-ZORA results are very close to the NESC-SORA results, except for the shielding of the I nucleus. Both the NESC-ZORA and NESC-SORA calculations yield very similar results to the previously reported values obtained using the relativistic infinite-order two-component coupled Hartree-Fock method. The difference between NESC-ZORA and NESC-SORA results is significant for the shieldings of iodine.     

A coupled-cluster study of the structure and vibrational spectra of pyrazole and imidazole

Harmonic force fields were calculated at the corresponding optimized geometries for pyrazole and imidazole at the HF, B3LYP, MP2, CCSD and CCSD(T) levels using the 6-31G* basis set and at the HF and B3LYP levels using the cc-pVTZ basis set. The agreement between the calculated and experimental geometries by the CCSD and CCSD(T) methods was generally similar to that obtained with the B3LYP and MP2 methods. The force fields were scaled using one-scale-factor (1SF), 3SF and 7SF scaling schemes. The scale factors were varied with respect to the experimental frequencies. Using 7SF scaling, the root-mean-square (RMS) deviation of the calculated frequencies from the experimental frequencies by the HF, B3LYP, MP2, CCSD and CCSD(T) methods and the 6-31G* basis set was 16, 7, 13, 11 and 11 cm(-1), respectively. This shows that the B3LYP method is preferred for force field calculations over the perturbative MP2, CCSD and CCSD(T) methods. Using 1SF scaling, the CCSD(T) scale factor was 0.931, the highest among the five methods used but close to that obtained with the B3LYP method and the cc-pVTZ basis set with lower RMS deviation.     

Analytic calculation of first-order molecular properties at the explicitly correlated second-order Moller-Plesset level: basis-set limits for the molecular quadrupole moments of BH and HF

The analytic calculation of first-order properties has been implemented in the DALTON program at the level of explicitly correlated second-order Moller-Plesset perturbation theory (MP2-R12). The implementation has been accomplished for MP2-R12 theory based on standard approximations A, A', and B, using an auxiliary basis for the resolution-of-the-identity approximation, with and without a frozen core. MP2-R12 first-order molecular properties have been calculated analytically for a few small test molecules. For BH and HF, the MP2-R12 results were supplemented with explicitly correlated coupled-cluster calculations (but at this level from numerical derivatives) including vibrational and relativistic corrections.     

Thermodynamic properties of the C5, C6, and C8 n-alkanes from ab initio electronic structure theory

The heats of formation for the n-alkanes C(n)H(n+2) for n = 5, 6, and 8 have been calculated using ab initio molecular orbital theory. Coupled-cluster calculations with perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation-consistent basis sets were used, up through the augmented quadruple zeta, to extrapolate to the complete basis set limit. Geometries were optimized at the B3LYP/TZVP and MP2/aug-cc-pVTZ levels. The MP2 geometries were used in the CCSD(T) calculations. Frequencies were determined at the density functional level (B3LYP/TZVP), and scaled zero point energies were calculated from the B3LYP frequencies. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. The core/valence corrections are not small, (approximately 1.1 kcal/mol per carbon unit) and cannot be neglected for chemical accuracy. The calculated deltaH(298)f values are -35.0, -40.2, and -50.2 kcal/mol for C5H12, C6H14, and C8H18, respectively, in excellent agreement with the respective experimental values of -35.11 +/- 0.19, -39.89 +/- 0.19, and -49.90 +/- 0.31 kcal/mol. Isodesmic reaction energies are presented for some simple reactions involving C8H18 and are shown not to be strongly method dependent.     

A combined experimental and quantum chemistry study of selenium chemical shift tensors

A comprehensive investigation of selenium chemical shift tensors is presented. Experimentally determined chemical shift tensors were obtained from solid-state 77Se NMR spectra for several organic, organometallic, or inorganic selenium-containing compounds. The first reported indirect spin-spin coupling between selenium and chlorine is observed for Ph(2)SeCl(2) where 1J(77Se,35Cl)iso is 110 Hz. Selenium magnetic shielding tensors were calculated for all of the molecules investigated using zeroth-order regular approximation density functional theory, ZORA DFT. The computations provide the orientations of the chemical shift tensors, as well as a test of the theory for calculating the magnetic shielding interaction for heavier elements. The ZORA DFT calculations were performed with nonrelativistic, scalar relativistic, and scalar with spin-orbit relativistic levels of theory. Relativistic contributions to the magnetic shielding tensor were found to be significant for (NH4)2WSe4 and of less importance for organoselenium, organophosphine selenide, and inorganic selenium compounds containing lighter elements.     

Interaction of dihydrogen with small and light molecules

Second-order M?ller-Plesset (MP2) calculations (using the approximate resolution of the identity, RI-MP2), explicitly correlated MP2 (MP2-R12) calculations, and coupled-cluster calculations including all single and double excitations with a perturbative estimate of triple excitations [CCSD(T)] are performed to study the interaction of molecular hydrogen with the small molecules HF, H2O, NH3, and LiOH. Different adsorption positions are studied. In the cases of H2O and NH3, the most favorable configuration places H2 in an end-on fashion on the O or N atom, respectively. In the cases of HF and LiOH, the H2 molecule takes a side-on position on the H atom of HF or the Li atom. With respect to MP2 calculations in a triple-zeta basis, both the enlargement of the basis set and the extension of the correlation treatment (CCSD(T) vs MP2) increase the interaction energy. The basis set limit CCSD(T) estimates of the interaction energy of H2 with the HF, H2O, NH3, and LiOH molecules amount to 4.40, 2.67, 3.02, and 10.74 kJ mol-1, respectively. The interaction energy for the simultaneous interaction of H2 with two LiOH molecules does not significantly exceed the value obtained for the interaction with a single LiOH molecule. Furthermore, the interaction energies (by MP2) of H2 with glycine, the glycine dimer, and imidazolium chloride amount to 2.78, 5.00, and 6.30 kJ mol-1, respectively.     

Calculation of electric-field gradients based on higher-order generalized Douglas-Kroll transformations

In this paper, the calculation of electric-field-like properties based on higher-order Douglas-Kroll-Hess (DKH) transformations is discussed. The electric-field gradient calculated within the Hartree-Fock self-consistent field framework is used as a representative property. The properties are expressed as an analytic first derivative of the four-component Dirac energy and the nth-order DKH energy, respectively. The differences between a "forward" transformation of the relativistic energy or the "back transformation" of the wave function is discussed in some detail. Detailed test calculations were carried out on the electric-field gradient at the halogen nucleus in the series HX (X=F,Cl,Br,I,At) for which extensive reference data are available. The DKH method is shown to reproduce (spin-free) four-component Dirac-Fock results to an accuracy of better than 99% which is significantly closer than previous DKH studies. The calculations of both the Hamiltonian and the property operator are shown to be essentially converged after the second-order transformation, even for elements as heavy as At. In addition, we have obtained results within the density-functional framework using the DKHZ and zeroth-order regular approximation (ZORA) methods. The latter results included picture-change effects at the scalar relativistic variant of the ZORA-4 level and were shown to be in quantitative agreement with earlier results obtained by van Lenthe and Baerends. The picture-change effects are somewhat smaller for the ZORA method compared to DKH. For heavier elements significant differences in the field gradients predicted by the two methods were found. Based on comparison with four-component Dirac-Kohn-Sham calculations, the DKH results are more accurate. Compared to the spin-free Dirac-Kohn-Sham reference values, the ZORA-4 formalism did not improve the results of the ZORA calculations.     

The Nature of the Binding of Au, Ag, and Pd to Benzene, Coronene, and Graphene: From Benchmark CCSD(T) Calculations to Plane-Wave DFT Calculations

The adsorption of Ag, Au, and Pd atoms on benzene, coronene, and graphene has been studied using post Hartree-Fock wave function theory (CCSD(T), MP2) and density functional theory (M06-2X, DFT-D3, PBE, vdW-DF) methods. The CCSD(T) benchmark binding energies for benzene-M (M = Pd, Au, Ag) complexes are 19.7, 4.2, and 2.3 kcal/mol, respectively. We found that the nature of binding of the three metals is different: While silver binds predominantly through dispersion interactions, the binding of palladium has a covalent character, and the binding of gold involves a subtle combination of charge transfer and dispersion interactions as well as relativistic effects. We demonstrate that the CCSD(T) benchmark binding energies for benzene-M complexes can be reproduced in plane-wave density functional theory calculations by including a fraction of the exact exchange and a nonempirical van der Waals correction (EE+vdW). Applying the EE+vdW method, we obtained binding energies for the graphene-M (M = Pd, Au, Ag) complexes of 17.4, 5.6, and 4.3 kcal/mol, respectively. The trends in binding energies found for the benzene-M complexes correspond to those in coronene and graphene complexes. DFT methods that use empirical corrections to account for the effects of vdW interactions significantly overestimate binding energies in some of the studied systems.     

Time-dependent quasirelativistic density-functional theory based on the zeroth-order regular approximation

A time-dependent quasirelativistic density-functional theory for excitation energies of systems containing heavy elements is developed, which is based on the zeroth-order regular approximation (ZORA) for the relativistic Hamiltonian and a noncollinear form for the adiabatic exchange-correlation kernel. To avoid the gauge dependence of the ZORA Hamiltonian a model atomic potential, instead of the full molecular potential, is used to construct the ZORA kinetic operator in ground-state calculations. As such, the ZORA kinetic operator no longer responds to changes in the density in response calculations. In addition, it is shown that, for closed-shell ground states, time-reversal symmetry can be employed to simplify the eigenvalue equation into an approximate form that is similar to that of time-dependent nonrelativistic density-functional theory. This is achieved by invoking an independent-particle approximation for the induced density matrix. The resulting theory is applied to investigate the global potential-energy curves of low-lying LambdaS- and omega omega-coupled electronic states of the AuH molecule. The derived spectroscopic parameters, including the adiabatic and vertical excitation energies, equilibrium bond lengths, harmonic and anharmonic vibrational constants, fundamental frequencies, and dissociation energies, are in good agreement with those of time-dependent four-component relativistic density-functional theory and ab initio multireference second-order perturbation theory. Nonetheless, this two-component relativistic version of time-dependent density-functional theory is only moderately advantageous over the four-component one as far as computational efforts are concerned.