A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

Studies of the trigonal field parameters and g factors for magnetic semiconductor NaCrS2

The trigonal field parameters v and v' of magnetic semiconductor NaCrS2 are calculated from the superposition model, and its g factors g parallel and g perpendicular are calculated from the high-order perturbation formulas of 3d3 ions in trigonal symmetry obtained by the one and two spin-orbit coupling parameter models. These calculations are based on the structural data of NaCrS2 crystal. The calculated results suggest that the superposition model is effective to the analyses of the low-symmetry field parameters of 3dn ions in crystals and that the two spin-orbit coupling parameter model (where the contributions from both the spin-orbit coupling parameter of 3dn ion and that of ligand are considered) is preferable to the one spin-orbit coupling parameter model (in which the contribution from only the spin-orbit coupling parameter of 3dn ion is considered) in the explanations of g factors in the cases that ligands have large spin-orbit coupling parameter in 3dn clusters.     

Determination of spin-orbit coupling contributions in the framework of density functional theory

We present a noniterative method to calculate spin-orbit coupling by means of a theoretical approach that provides the use of the full Breit-Pauli operator. This method was applied to compute one and two-electron spin-orbit coupling contributions between singlet and triplet, and doublet and doublet states, respectively. These states have been represented by monodeterminantal wave functions and optimized using the PW91 gradient-corrected exchange-correlation functional and the hybrid B3LYP one. They have been supplied by the conventional density functional theory packages, and thus coupled by our spin-orbit coupling code. Different size basis sets have been employed and the obtained results have been compared with the corresponding ones provided by some of the already existing methods and with the experimental data. They have been found to be in good quantitative agreement.     

Structural and spectroscopic study of the excited electronic states of silver dihalides by quantum chemical methods

The structural and electronic properties of the excited electronic states of AgX(2) (X = F, Cl, Br, and I), have been calculated, taking electron correlation and spin-orbit coupling into account and employing improved relativistic-effective-core potentials for silver and the halogen atoms. The relative ordering of the excited states of these molecules has been discussed via molecular-orbital arguments. The spin-orbit splittings of three degenerate electronic states ((2)Pi(g), (2)Pi(u), and (2)Delta(g)) have been calculated and the spin-orbit induced inter-state (Sigma - Pi) coupling has been discussed. The composition of the spin-orbit eigenstates is analyzed in terms of scalar-relativistic electronic states. Finally, a theoretical prediction of the photodetachment bands of the title molecules has been accomplished.     

Ab initio calculation of anisotropic magnetic properties of complexes. I. Unique definition of pseudospin Hamiltonians and their derivation

A methodology for the rigorous nonperturbative derivation of magnetic pseudospin Hamiltonians of mononuclear complexes and fragments based on ab initio calculations of their electronic structure is described. It is supposed that the spin-orbit coupling and other relativistic effects are already taken fully into account at the stage of quantum chemistry calculations of complexes. The methodology is based on the establishment of the correspondence between the ab initio wave functions of the chosen manifold of multielectronic states and the pseudospin eigenfunctions, which allows to define the pseudospin Hamiltonians in the unique way. Working expressions are derived for the pseudospin Zeeman and zero-field splitting Hamiltonian corresponding to arbitrary pseudospins. The proposed calculation methodology, already implemented in the SINGLE_ANISO module of the MOLCAS-7.6 quantum chemistry package, is applied for a first-principles evaluation of pseudospin Hamiltonians of several complexes exhibiting weak, moderate, and very strong spin-orbit coupling effects.     

Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules

Terms arising from the relativistic spin-orbit effect on both hyperfine and Zeeman interactions are introduced to density-functional theory calculation of nuclear magnetic shielding in paramagnetic molecules. The theory is a generalization of the former nonrelativistic formulation for doublet systems and is consistent to O(alpha4), the fourth power of the fine structure constant, for the spin-orbit terms. The new temperature-dependent terms arise from the deviation of the electronic g tensor from the free-electron g value as well as spin-orbit corrections to hyperfine coupling tensor A, the latter introduced in the present work. In particular, the new contributions include a redefined isotropic pseudocontact contribution that consists of effects due to both the g tensor and spin-orbit corrections to hyperfine coupling. The implementation of the spin-orbit terms makes use of all-electron atomic mean-field operators and/or spin-orbit pseudopotentials. Sample results are given for group-9 metallocenes and a nitroxide radical. The new O(alpha4) corrections are found significant for the metallocene systems while they obtain small values for the nitroxide radical. For the isotropic shifts, none of the three beyond-leading-order hyperfine contributions are negligible.     

Theoretical calculation of the photodetachment spectra of XAuY- (X, Y = Cl, Br, and I)

The photodetachment spectra of the title molecules have been calculated, taking electron correlation and spin-orbit coupling into account and employing improved relativistic effective core potentials for gold and the halogen atoms. The calculated spectra have been compared with existing experimental spectra. The spin-orbit splitting of several degenerate electronic states has been calculated. The composition of the spin-orbit eigenstates are analyzed in terms of scalar relativistic electronic states. A comparison of the relative position of peaks in the calculated photodetachment spectra of the title molecules has been made.     

Calculation of the vibronic structure of the photodetachment spectra of CCCl- and CCBr-

The vibronic structure of the closely spaced and strongly coupled X 2Sigma+ and A 2Pi states in the photodetachment spectra of CCCl- and CCBr- has been calculated by considering Sigma-Pi vibronic coupling together with spin-orbit coupling. The stretching modes are treated within the so-called linear-vibronic-coupling model. The vibronic and spin-orbit parameters have been determined by accurate ab initio electronic-structure calculations. While the nonrelativistic vibronic-coupling parameters are of approximately equal strength in CCCl and CCBr, the vibronic-coupling parameters of spin-orbit origin are found to be larger in the latter. The calculated photodetachment spectra of both systems are shown to exhibit a complicated vibronic structure due to strong Sigma-Pi vibronic coupling. The spectral envelopes of the calculated photodetachment spectra exhibit a double-hump reminiscent of strongly coupled Exe Jahn-Teller systems.     

Reinvestigation of the electronic spectroscopy of the Au-Ar complex

The Au-Ar complex is reinvestigated employing resonance-enhanced multiphoton ionization spectroscopy. Spectra are reported, corresponding to the atomic transition Au(6p<--6s). This electronic excitation yields (2)Pi and (2)Sigma(+) states of Au-Ar, which interact under the influence of spin-orbit coupling. The spectra are consistent with strong sigma-pi mixing induced by the large spin-orbit coupling of Au, leading to strong interaction of the two Omega=12 states, which arise from the Ar((1)S(0))+Au((2)P(12,32)) asymptotes, and the consequent formation of a "shelf" on the outer wall of the lowest Omega=12 state. In addition, high-level ab initio calculations are reported on the ground electronic state, X (2)Sigma(+), including extrapolation to the basis set limit.     

Extensive theoretical study on the low-lying electronic states of silicon monofluoride cation including spin-orbit coupling

Ab initio calculations on the low-lying electronic states of SiF+ are performed using the internally contracted multireference configuration interaction method with the Davidson correction and entirely uncontracted aug-cc-pV5Z basis set. The effects of spin-orbit coupling are accounted for by the state interaction approach with the full Breit-Pauli Hamiltonian. The entire 23 Omega states generated from the 12 valence Lambda-S states, which correlate with the first dissociation channel are studied for the first time. Good agreement is found between the calculated results and the available experimental data. The spin-orbit coupling effects on the potential energy curves and spectroscopic properties are studied. Various curve crossings are revealed, which could lead to the predissociation of the a3Pi, A1Pi, and (2)3Sigma+ states and the predissociation pathways are analyzed based upon the calculated spin-orbit matrix elements. The calculated ionization potentials of the ground-state SiF to a few states of SiF+ are in good agreement with the available experimental measurements. Moreover, the transition dipole moments of the dipole-allowed transitions and the transition properties for the A3Pi0+ -X1Sigma+ 0+ and B3Pi1-X1Sigma+ 0+ transitions are predicted, including the Franck-Condon factors and the radiative lifetimes.     

Spin-orbit coupling in O2(upsilon)+O2 collisions: I. Electronic structure calculations on dimer states involving the X 3Sigmag-, a 1Deltag, and b 1Sigmag+ states of O2

The importance of vibrational-to-electronic (V-E) energy transfer mediated by spin-orbit coupling in the collisional removal of O2(X 3Sigmag-,upsilon>or=26) by O2 has been reported in a recent communication [F. Dayou, J. Campos-Martinez, M. I. Hernandez, and R. Hernandez-Lamoneda, J. Chem. Phys. 120, 10355 (2004)]. The present work provides details on the electronic properties of the dimer (O2)2 relevant to the self-relaxation of O2(X 3Sigmag-,upsilon>0) where V-E energy transfer involving the O2(a 1Deltag) and O2(b 1Sigmag+) states is incorporated. Two-dimensional electronic structure calculations based on highly correlated ab initio methods have been carried out for the potential-energy and spin-orbit coupling surfaces associated with the ground singlet and two low-lying excited triplet states of the dimer dissociating into O2(X 3Sigmag-)+O2(X 3Sigmag-), O2(a 1Deltag)+O2(X 3Sigmag-), and O2(b 1Sigmag+)+O2(X 3Sigmag-). The resulting interaction potentials for the two excited triplet states display very similar features along the intermolecular separation, whereas differences arise with the ground singlet state for which the spin-exchange interaction produces a shorter equilibrium distance and higher binding energy. The vibrational dependence is qualitatively similar for the three studied interaction potentials. The spin-orbit coupling between the ground and second excited states is already nonzero in the O2+O2 dissociation limit and keeps its asymptotic value up to relatively short intermolecular separations, where the coupling increases for intramolecular distances close to the equilibrium of the isolated diatom. On the other hand, state mixing between the two excited triplet states leads to a noticeable collision-induced spin-orbit coupling between the ground and first excited states. The results are discussed in terms of specific features of the dimer electronic structure (including a simple four-electron model) and compared with existing theoretical and experimental data. This work gives theoretical insight into the origin of electronic energy-transfer mechanisms in O2+O2 collisions.     

Constructing spin-adiabatic states for the modeling of spin-crossing reactions. I. A shared-orbital implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a “transition state” structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions (predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄), the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (ie, a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.     

Theoretical studies of nuclear magnetic resonance parameters for the proton-exchange pathways in porphyrin and porphycene

The nuclear magnetic resonance (NMR) parameters in porphyrin and porphycene have been calculated to investigate their changes during the process of proton exchange, using density-functional theory (DFT) for both the spin-spin coupling constants and the shielding constants. In addition, in calculations on the smaller 1,3-bis(arylimino)isoindoline molecule, we have tested the performance of our computational approach against experimental data. The calculated nuclear spin-spin coupling constants and shielding constants have been analyzed as functions of the progress of the proton transfer between two nitrogen atoms. The one-bond couplings between proton and nitrogen, dominated by the Fermi-contact term, decay steeply as the internuclear distance increases. The small changes in the intramolecular J(HH) coupling between two inner protons are mainly determined by the sum of relatively large spin-orbit terms. The isotropic shielding constant shows a strong deshielding of the nitrogen nuclei as the proton migrates away. Both the isotropic shielding of the exchanged protons and the shielding anisotropy exhibit a minimum close to the transition states.     

Jahn-Teller and related effects in the silver trimer. II: vibrational analysis of the A 2E"-X 2E' electronic transition

The laser-excited, jet-cooled A 2E"-X 2E' electronic spectrum of the silver trimer yields detailed information about its A- and X-state vibronic structure. Following extensive parameter fitting, the absorption and emission spectra are simulated and the bands are assigned. The Jahn-Teller analysis includes both linear and quadratic coupling terms, considered simultaneously with spin-orbit coupling. The spin-orbit splitting is shown to be largely quenched in both the A and X electronic states. The Jahn-Teller analysis of the A and X vibronic structures reveals the distortion of their corresponding potential energy surfaces.     

A theoretical study of the excited states of AmO2n+, n=1,2,3

The ground and excited states of the AmO(2) (+), AmO(2) (2+), and AmO(2) (3+) ions have been studied using the four-component configuration interaction singles doubles, spin-orbit complete active space self-consistent field, and spin-orbit complete active space-order perturbation theory methods. The roles of scalar relativistic effects and spin-orbit coupling are analyzed; results with different methods are carefully compared by a precise analysis of the wave functions. A molecular spinor diagram is used in relation to the four-component calculations while a ligand field model is used for the two-step method. States with the same number of electrons in the four nonbonding orbitals are in very good agreement with the two methods while ligand field and charge transfer states do not have the same excitation energies.     

Renner-Teller and spin-orbit vibronic coupling effects in linear triatomic molecules with a half-filled pi shell

The vibronic and spin-orbit-induced interactions among the (3)Sigma(-), (1)Delta, and (1)Sigma(+) electronic states arising from a half-filled pi orbital of a linear triatomic molecule are considered, employing the microscopic (Breit-Pauli) spin-orbit coupling operator. The 6 x 6 Hamiltonian matrix is derived in a diabatic spin-orbital electronic basis set, including terms up to fourth order in the expansion of the molecular Hamiltonian in the bending normal coordinate about the linear geometry. The symmetry properties of the Hamiltonian are analyzed. Aside from the nonrelativistic fourth-order Renner-Teller vibronic coupling within the (1)Delta state and the second-order nonrelativistic vibronic coupling between the (1)Sigma(+) and (1)Delta states, there exist zeroth-order, first-order, as well as third-order vibronic coupling terms of spin-orbit origin. The latter are absent when the phenomenological expression for the spin-orbit coupling operator is used instead of the microscopic form. The effects of the nonrelativistic and spin-orbit-induced vibronic coupling mechanisms on the (3)Sigma(-), (1)Delta, and (1)Sigma(+) adiabatic potential energy surfaces as well as on the spin-vibronic energy levels are discussed for selected parameter values.     

Spin-orbit vibronic coupling in 3Pi states of linear triatomic molecules

The Renner-Teller vibronic-coupling problem of a 3Pi electronic state of a linear molecule is analyzed with the inclusion of the spin-orbit coupling of the 3Pi electronic state, employing the microscopic (Breit-Pauli) spin-orbit coupling operator for the two unpaired electrons. The 6x6 Hamiltonian matrix in a diabatic spin-electronic basis is obtained by an expansion of the molecular Hamiltonian in powers of the bending amplitude. The symmetry properties of the Hamiltonian with respect to the time-reversal operator and the relativistic vibronic angular momentum operator are analyzed. It is shown that there exists a linear vibronic-coupling term of spin-orbit origin, which has not been considered so far in the Renner-Teller theory of 3Pi electronic states. While two of the six adiabatic electronic wave functions do not exhibit a geometric phase, the other four carry nontrivial topological phases which depend on the radius of the integration contour. The spectroscopic effects of the linear spin-orbit vibronic-coupling mechanism have been analyzed by numerical calculations of the vibronic spectrum for selected model examples.     

Investigating the M*He exciplexes, M={Li,Na,K,Rb,Cs,Fr}: density functional approach

Potential curves for the ground and the first lowest excited states of the MHe (where M={Li,Na,K,Rb,Cs,Fr}) exciplexes are calculated using the density functional theory (DFT) formalism. Relativistic calculations are carried out with and without spin-orbit (SO) coupling effect, using a zeroth order regular approximation (ZORA) approach. The depth De and position Re parameters of the potential curves for the case without spin-orbit effect are presented and compared with other works. Potential curves for Li-He, Cs-He, and Fr-He without spin-orbit effect and Cs-He with the spin-orbit effect are shown. A bond analysis is presented too, since the central interest of the study of the exciplexes is the nature and the existence of the bonding states that are the origin of emission spectra observed experimentally for these systems.     

DFT studies on the magnetic exchange across the cyanide bridge

Exchange coupling across the cyanide bridge in a series of novel cyanometalate complexes with CuII-NC-MIII (M = Cr and low-spin Mn, Fe) fragments has been studied using the broken-symmetry DFT approach and an empirical model, which allows us to relate the exchange coupling constant with sigma-, pi-, and pi*-type spin densities of the CN- bridging ligand. Ferromagnetic exchange is found to be dominated by pi-delocalization via the CN- pi pathway, whereas spin polarization with participation of sigma orbitals (in examples, where the dz2 orbital of MIII is empty) and pi* orbitals of CN- yields negative spin occupations in these orbitals, and reduces the CuII-MIII exchange coupling constant. When the dz2 orbital of MIII is singly occupied, an additional positive spin density appears in the sigma(CN) orbital and leads to an increase of the ferromagnetic Cu-NC-M exchange constant. For low-spin [MIII(CN)6]3- complexes, the dz2 orbital occupancy results in high-spin metastable excited states, and this offers interesting aspects for applications in the area of molecular photomagnetism. The DFT values of the exchange coupling parameters resulting from different occupations of the t2g orbitals of low-spin (t2g5) FeIII are used to discuss the effect of spin-orbit coupling on the isotropic and anisotropic exchange coupling in linear Cu-NC-Fe pairs.