Search

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn–Sham method. The light–matter interaction is included in the Jaynes–Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light–matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined. [ABSTRACT FROM AUTHOR]

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems. [ABSTRACT FROM AUTHOR]

The use of the mixed reference (MR) reduced density matrix, which combines reduced density matrices of the MS = +1 and -1 triplet-ground states, is proposed in the context of the collinear spin-fliptime-dependent density functional theory (SF-TDDFT) methodology. The time-dependent Kohn-Sham equation with the mixed state is solved by the use of spinor-like open-shell orbitals within the linear response formalism, which enables to generate additional configurations in the realm of TD-DFT. The resulting MR-SF-TDDFT computational scheme has several advantages before the conventional collinear SF-TDDFT. The spin-contamination of the response states of SF-TDDFT is nearly removed. This considerably simplifies the identification of the excited states, especially in the "black-box" type applications, such as the automatic geometry optimization, reaction path following, or molecular dynamics simulations. With the new methodology, the accuracy of the description of the excited states is improved as compared to the collinear SF-TDDFT. Several test examples, which include systems typified by strong non-dynamic correlation, orbital (near) degeneracy, and conical intersections, are given to illustrate the performance of the new method. [ABSTRACT FROM AUTHOR]

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the newmethod correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena. [ABSTRACT FROM AUTHOR]

The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham(REKS) method and its state interaction (SI) extension, SI-SA-REKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. In this work, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. The implementation and benchmarking of the derived formalism will be described in a subsequent article of this series. [ABSTRACT FROM AUTHOR]

The spin-restricted ensemble-referenced Kohn-Sham (REKS) method is based on an ensemble representation of the density and is capable of correctly describing the non-dynamic electron correlation stemming from (near-)degeneracy of several electronic configurations. The existing REKS methodology describes systems with two electrons in two fractionally occupied orbitals. In thiswork, the REKS methodology is extended to treat systems with four fractionally occupied orbitals accommodating four electrons and self-consistent implementation of the REKS(4,4) method with simultaneous optimization of the orbitals and their fractional occupation numbers is reported. The new method is applied to a number of molecular systems where simultaneous dissociation of several chemical bonds takes place, as well as to the singlet ground states of organic tetraradicals 2,4-didehydrometaxylylene and 1,4,6,9-spiro[4.4]nonatetrayl. [ABSTRACT FROM AUTHOR]

The analytical gradient for the two-component Normalized Elimination of the Small Component (2c-NESC) method is presented. The 2c-NESC is a Dirac-exact method that employs the exact two-component one-electron Hamiltonian and thus leads to exact Dirac spin-orbit (SO) splittings for one-electron atoms. For many-electron atoms and molecules, the effect of the two-electron SO interaction is modeled by a screened nucleus potential using effective nuclear charges as proposed by Boettger [Phys. Rev. B 62, 7809 (2000)]. The effect of spin-orbit coupling (SOC) on molecular geometries is analyzed utilizing the properties of the frontier orbitals and calculated SO couplings. It is shown that bond lengths can either be lengthened or shortened under the impact of SOC where in the first case the influence of low lying excited states with occupied antibonding orbitals plays a role and in the second case the j j-coupling between occupied antibonding and unoccupied bonding orbitals dominates. In general, the effect of SOC on bond lengths is relatively small (=5% of the scalar relativistic changes in the bond length). However, large effects are found for van der Waals complexes Hg2 and Cn2, which are due to the admixture of more bonding character to the highest occupied spinors. [ABSTRACT FROM AUTHOR]

State-averaged (SA) variants of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, SA-REKS and state-interaction (SI)-SA-REKS, implement ensemble density functional theory for variationally obtaining excitation energies of molecular systems. In this work, the currently existing version of the SA-REKS method, which included only one excited state into the ensemble averaging, is extended by adding more excited states to the averaged energy functional. A general strategy for extension of the REKS-type methods to larger ensembles of ground and excited states is outlined and implemented in extended versions of the SA-REKS and SI-SA-REKS methods. The newly developed methods are tested in the calculation of several excited states of ground-state multi-reference systems, such as dissociating hydrogen molecule, and excited states of donor-acceptor molecular systems. For hydrogen molecule, the new method correctly reproduces the distance dependence of the lowest excited state energies and describes an avoided crossing between the doubly excited and singly excited states. For bithiophene-perylenediimide stacked complex, the SI-SA-REKS method correctly describes crossing between the locally excited state and the charge transfer excited state and yields vertical excitation energies in good agreement with the ab initio wavefunction methods. [ABSTRACT FROM AUTHOR]

Electron transfer in the ground and excited states of a model donor-acceptor (D-A) system is investigated using the single-reference and multi-reference density functional theory (DFT) methods. To analyze the results of the calculations, a simple two-site multi-reference model was derived that predicts a stepwise electron transfer in the S0 state and a wave-like dependence of the S1 electron transfer on the external stimulus. The standard single-reference Kohn-Sham (KS) DFT approach and the time-dependent DFT (TDDFT) method failed to describe the correct dependence of the S0 and S1 electron transfer on the external electric field applied along the donor-acceptor system. The multi-reference DFT approach, the spin-restricted ensemble-referenced KS (REKS) method, was able to successfully reproduce the correct behavior of the S0 and S1 electron transfer on the applied field. The REKS method was benchmarked against experimentally measured gas phase charge transfer excitations in a series of organic donor-acceptor complexes and displayed its ability to describe this type of electronic transitions with a very high accuracy, mean absolute error of 0.05 eV with the use of the standard range separated density functionals. On the basis of the calculations undertaken in this work, it is suggested that the non-adiabatic coupling between the S0 and S1 states may interfere with the electron transfer in a weakly coupled donor-acceptor system. It is also suggested that the electronic excitation of a D+-A- system may play a dual role by assisting the further electron transfer at certain magnitudes of the applied electric field and causing the backward transfer at lower electric field strengths. [ABSTRACT FROM AUTHOR]

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 À or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems. [ABSTRACT FROM AUTHOR]

Computational investigation of the longest wavelength excitations in a series of cyanines and linear n-acenes is undertaken with the use of standard spin-conserving linear response time-dependent density functional theory (TD-DFT) as well as its spin-flip variant and a ΔSCF method based on the ensemble DFT. The spin-conserving linear response TD-DFT fails to accurately reproduce the lowest excitation energy in these π-conjugated systems by strongly overestimating the excitation energies of cyanines and underestimating the excitation energies of n-acenes. The spin-flip TD-DFT is capable of correcting the underestimation of excitation energies of n-acenes by bringing in the non-dynamic electron correlation into the ground state; however, it does not fully correct for the overestimation of the excitation energies of cyanines, for which the non-dynamic correlation does not seem to play a role. The ensemble DFT method employed in this work is capable of correcting for the effect of missing non-dynamic correlation in the ground state of n-acenes and for the deficient description of differential correlation effects between the ground and excited states of cyanines and yields the excitation energies of both types of extended π-conjugated systems with the accuracy matching high-level ab initio multireference calculations. [ABSTRACT FROM AUTHOR]

Analytic second derivatives of the relativistic energy for the calculation of electric response properties are derived utilizing the normalized elimination of the small component (NESC) method. Explicit formulas are given for electric static dipole polarizabilities and infrared intensities by starting at the NESC representation of electric dipole moments. The analytic derivatives are implemented in an existing NESC program and applied to calculate dipole moments, polarizabilities, and the infrared spectra of gold- and mercury-containing molecules as well as some actinide molecules. Comparison with experiment reveals the accuracy of NESC second order electric response properties. [ABSTRACT FROM AUTHOR]

Based on the analytic derivatives formalism for the spin-free normalized elimination of the small component method, a new computational scheme for the calculation of the electric field gradient at the atomic nuclei was developed and presented. The new computational scheme was tested by the calculation of the electric field gradient at the mercury nucleus in a series of Hg-containing inorganic and organometallic compounds. The benchmark calculations demonstrate that the new formalism is capable of reproducing experimental and theoretical reference data with high accuracy. The method developed can be routinely applied to the calculation of large and very large molecules and holds considerable promise for the interpretation of the experimental data of biologically relevant compounds containing heavy elements. [ABSTRACT FROM AUTHOR]

Eukaryotic genome is a highly compacted nucleoprotein complex organized in a hierarchical structure based on nucleosomes. Detailed organization of this structure remains unknown. In the present work we developed algorithms for geometry modeling of the supernucleosomal chromatin structure and for computing distance distribution functions and small-angle neutron scattering (SANS) spectra of the genome-scale (∼106 nucleosomes) chromatin structure at residue resolution. Our physical nucleosome model was based on the mononucleosome crystal structure. A nucleosome was assumed to be rigid within a local coordinate system. Interface parameters between nucleosomes can be set for each nucleosome independently. Pair distance distributions were computed with Monte Carlo simulation. SANS spectra were calculated with Fourier transformation of weighted distance distribution; the concentration of heavy water in solvent and probability of H/D exchange were taken into account. Two main modes of supernucleosomal structure generation were used. In a free generation mode all interface parameters were chosen randomly, whereas nucleosome self-intersections were not allowed. The second generation mode (generation in volume) enabled spherical or cubical wall restrictions. It was shown that calculated SANS spectra for a number of our models were in general agreement with available experimental data. [ABSTRACT FROM AUTHOR]

The analytical energy gradient of the normalized elimination of the small component (NESC) method is derived for the first time and implemented for the routine calculation of NESC geometries and other first order molecular properties. Essential for the derivation is the correct calculation of the transformation matrix U relating the small component to the pseudolarge component of the wavefunction. The exact form of ∂U/∂λ is derived and its contribution to the analytical energy gradient is investigated. The influence of a finite nucleus model and that of the picture change is determined. Different ways of speeding up the calculation of the NESC gradient are tested. It is shown that first order properties can routinely be calculated in combination with Hartree-Fock, density functional theory (DFT), coupled cluster theory, or any electron correlation corrected quantum chemical method, provided the NESC Hamiltonian is determined in an efficient, but nevertheless accurate way. The general applicability of the analytical NESC gradient is demonstrated by benchmark calculations for NESC/CCSD (coupled cluster with all single and double excitation) and NESC/DFT involving up to 800 basis functions. [ABSTRACT FROM AUTHOR]

The isomer shift for the 23.87 keV M1 resonant transition in the 119Sn nucleus is calibrated with the help of ab initio calculations. The calibration constant α(119Sn) obtained from Hartree–Fock (HF) calculations (αHF(119Sn)=(0.081±0.002)a0-3 mm/s) and from second-order Mo\ller–Plesset (MP2) calculations (αMP2(119Sn)=(0.091±0.002)a0-3 mm/s) are in good agreement with the previously obtained values. The importance of a proper treatment of electron correlation effects is demonstrated on the basis of a statistical analysis of the results of the calibration. The approach used in the calibration is applied to study the 119Sn isomer shift in CaSnO3 perovskite under pressure. Comparison with the experimental results for the pressure range of 0–36 GPa shows that the current methodology is capable of describing tiny variations of isomer shift with reasonable accuracy. [ABSTRACT FROM AUTHOR]

The orbital products of occupied and virtual orbitals are employed as an expansion basis for the charge density generating the local potential in the optimized effective potential method thus avoiding the use of auxiliary basis sets. The high computational cost arising from the quadratic increase of the dimension of this product basis with system size can be greatly reduced by elimination of the linearly dependent products according to a procedure suggested by Beebe and Linderberg [Int. J. Quantum Chem. 12, 683 (1977)]. Numerical results from this approach show a very good agreement with those obtained from balancing the auxiliary basis for the expansion of the local potential with the orbital basis set. [ABSTRACT FROM AUTHOR]

The optimized effective potential (OEP) equations are solved in a matrix representation using the orbital products of occupied and virtual orbitals for the representation of both the local potential and the response function. This results in a direct relationship between the matrix elements of local and nonlocal operators for the exchange-correlation potential. The effect of the truncation of the number of such products in the case of finite orbital basis sets on the OEP orbital and total energies and on the spectrum of eigenvalues of the response function is examined. Test calculations for Ar and Ne show that rather large AO basis sets are needed to obtain an accurate representation of the response function. [ABSTRACT FROM AUTHOR]

A quantum chemical computational scheme for the calculation of isomer shift in Mössbauer spectroscopy is suggested. Within the described scheme, the isomer shift is treated as a derivative of the total electronic energy with respect to the radius of a finite nucleus. The explicit use of a finite nucleus model in the calculations enables one to incorporate straightforwardly the effects of relativity and electron correlation. The results of benchmark calculations carried out for several iron complexes as well as for a number of atoms and atomic ions are presented and compared with the available experimental and theoretical data. [ABSTRACT FROM AUTHOR]

It is demonstrated that the LYP correlation functional is not suited to be used for the calculation of electron spin resonance hyperfine structure (HFS) constants, nuclear magnetic resonance spin-spin coupling constants, magnetic, shieldings and other properties that require a balanced account of opposite- and equal-spin correlation, especially in the core region. In the case of the HFS constants of alkali atoms, LYP exaggerates opposite-spin correlation effects thus invoking too strong in-out correlation effects, an exaggerated spin-polarization pattern in the core shells of the atoms, and, consequently, too large HFS constants. Any correlation functional that provides a balanced account of opposite- and equal-spin correlation leads to improved HFS constants, which is proven by comparing results obtained with the LYP and the PW91 correlation functional. It is suggested that specific response properties are calculated with the PW91 rather than the LYP correlation functional. [ABSTRACT FROM AUTHOR]

The regular approximation to the normalized elimination of the small component (NESC) in the modified Dirac equation has been developed and presented in matrix form. The matrix form of the infinite-order regular approximation (IORA) expressions, obtained in [Filatov and Cremer, J. Chem. Phys. 118, 6741 (2003)] using the resolution of the identity, is the exact matrix representation and corresponds to the zeroth-order regular approximation to NESC (NESC-ZORA). Because IORA (=NESC-ZORA) is a variationally stable method, it was used as a suitable starting point for the development of the second-order regular approximation to NESC (NESC-SORA). As shown for hydrogenlike ions, NESC-SORA energies are closer to the exact Dirac energies than the energies from the fifth-order Douglas–Kroll approximation, which is much more computationally demanding than NESC-SORA. For the application of IORA (=NESC-ZORA) and NESC-SORA to many-electron systems, the number of the two-electron integrals that need to be evaluated (identical to the number of the two-electron integrals of a full Dirac–Hartree–Fock calculation) was drastically reduced by using the resolution of the identity technique. An approximation was derived, which requires only the two-electron integrals of a nonrelativistic calculation. The accuracy of this approach was demonstrated for heliumlike ions. The total energy based on the approximate integrals deviates from the energy calculated with the exact integrals by less than 5×10-9 hartree units. NESC-ZORA and NESC-SORA can easily be implemented in any nonrelativistic quantum chemical program. Their application is comparable in cost with that of nonrelativistic methods. The methods can be run with density functional theory and any wave function method. NESC-SORA has the advantage that it does not imply a picture change. [ABSTRACT FROM AUTHOR]

A simple modification of the zeroth-order regular approximation (ZORA) in relativistic theory is suggested to suppress its erroneous gauge dependence to a high level of approximation. The method, coined gauge-independent ZORA (ZORA-GI), can be easily installed in any existing nonrelativistic quantum chemical package by programming simple one-electron matrix elements for the quasirelativistic Hamiltonian. Results of benchmark calculations obtained with ZORA-GI at the Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) level for dihalogens X2 (X=F,Cl,Br,I,At) are in good agreement with the results of four-component relativistic calculations (HF level) and experimental data (MP2 level). ZORA-GI calculations based on MP2 or coupled-cluster theory with single and double perturbations and a perturbative inclusion of triple excitations [CCSD(T)] lead to accurate atomization energies and molecular geometries for the tetroxides of group VIII elements. With ZORA-GI/CCSD(T), an improved estimate for the atomization energy of hassium (Z=108) tetroxide is obtained. © 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The infinite-order regular approximation (IORA) and IORA with modified metric (IORAmm) is used to develop an algorithm for calculating relativistically corrected isotropic hyperfine structure (HFS) constants. The new method is applied to the calculation of alkali atoms Li–Fr, coinage metal atoms Cu, Ag, and Au, the Hg+ radical ion, and the mercury containing radicals HgH, HgCH3, HgCN, and HgF. By stepwise improvement of the level of theory from Hartree–Fock to second-order Møller–Plesset theory and to quadratic configuration interaction theory with single and double excitations, isotropic HFS constants of high accuracy were obtained for atoms and for molecular radicals. The importance of relativistic corrections is demonstrated. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A new method for calculating the indirect nuclear spin–spin coupling constant within the regular approximation to the exact relativistic Hamiltonian is presented. The method is completely analytic in the sense that it does not employ numeric integration for the evaluation of relativistic corrections to the molecular Hamiltonian. It can be applied at the level of conventional wave function theory or density functional theory. In the latter case, both pure and hybrid density functionals can be used for the calculation of the quasirelativistic spin–spin coupling constants. The new method is used in connection with the infinite-order regular approximation with modified metric (IORAmm) to calculate the spin–spin coupling constants for molecules containing heavy elements. The importance of including exact exchange into the density functional calculations is demonstrated. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The exact relativistic Hamiltonian for electronic states is expanded in terms of energy-independent linear operators within the regular approximation. An effective relativistic Hamiltonian has been obtained, which yields in lowest order directly the infinite-order regular approximation (IORA) rather than the zeroth-order regular approximation method. Further perturbational expansion of the exact relativistic electronic energy utilizing the effective Hamiltonian leads to new methods based on ordinary (IORAn) or double [IORAn(2)] perturbation theory (n: order of expansion), which provide improved energies in atomic calculations. Energies calculated with IORA4 and IORA3(2) are accurate up to c[sup -20]. Furthermore, IORA is improved by using the IORA wave function to calculate the Rayleigh quotient, which, if minimized, leads to the exact relativistic energy. The outstanding performance of this new IORA method coined scaled IORA is documented in atomic and molecular calculations. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Analytic expressions for the derivatives of the total molecular energy with respect to external electric field are derived within the regular approximation to the full four-component relativistic Hamiltonian and presented in matrix form suitable for implementation in standard quantum-chemical codes. Results of benchmark calculations using the infinite-order regular approximation with modified metric method are presented and discussed. The static electric dipole polarizabilities of group VIII metal tetroxides MO[SUB4] for M = Ru, Os, Hs (Z5108) are studied with the help of second-order Møller-Plesset perturbation theory using the infinite-order regular approximation with modified metric Hamiltonian. The polarizabilities obtained vary in the sequence RuO[SUB4] > OsO[SUB4] > HsO[SUB4], which is different from those obtained in other studies. However, it is in line with calculated [SUP1]T[SUB2] ← [SUP1]A[SUB1] excitation energies of the group VIII tetroxides, which provide a measure for the magnitude of their polarizabilities. [ABSTRACT FROM AUTHOR]

A new method for relativistically corrected nuclear magnetic resonance (NMR) chemical shifts is developed by combining the individual gauge for the localized orbital approach for density functional theory with the normalized elimination of a small component using an effective potential. The new method is used for the calculation of the NMR chemical shifts of [SUP95]Mo and [SUP183]W in various molybdenum and tungsten compounds. It is shown that quasirelativistic corrections lead to an average improvement of calculated NMR chemical shift values by 300 and 120 ppm in the case of [SUP95]Mo and [SUP183]W, respectively, which is mainly due to improvements in the paramagnetic contributions. The relationship between electronic structure of a molecule and the relativistic paramagnetic corrections is discussed. Relativistic effects for the diamagnetic part of the magnetic shielding caused by a relativistic contraction of the s,p orbitals in the core region concern only the shielding values, however, have little consequence for the shift values because of the large independence from electronic structure and a cancellation of these effects in the shift values. It is shown that the relativistic corrections can be improved by level shift operators and a B3LYP hybrid functional, for which Hartree-Fock exchange is reduced to 15%. [ABSTRACT FROM AUTHOR]

Analytic expressions are derived for the evaluation of derivatives of the total molecular energy with respect to external parameters (nuclear coordinates, external electric fields, etc.) within the relativistic regular approximation. The presented formalism employs the spectral resolution of the identity avoiding, however, the explicit use of an auxiliary basis set in the calculation of the matrix elements of the regular relativistic Hamiltonian. The final formulas for the total energy and energy derivatives are presented in matrix form suitable for implementation into standard quantum chemical packages. Results of benchmark calculations for gold containing diatomic molecules and for xenone hexafluoride performed at the Hartree-Fock and various correlation corrected levels of theory are presented and discussed. [ABSTRACT FROM AUTHOR]

Investigates various atomic and molecular multiplet states using a proposed spin-restricted open-shell Kohn-Sham (ROKS) method. Range of multiplets considered; ROKS equations; States of atomic species; Angular momentum states of linear molecules; Diradicaloid states of perpendicular ethylene and of square cyclobutadiene.

The connection between the exact quasirelativistic approach developed in the title reference [W. Kutzelnigg and W. Liu, J. Chem. Phys. 123, 241102 (2005)] and the method of elimination of the small component in matrix form developed previously by Dyall is explicitly worked out. An equation that links Hermitian and non-Hermitian formulations of the exact quasirelativistic theory is derived. Besides establishing a kinship between the existing formulations, the proposed equation can be employed for the derivation of new formulations of the exact quasirelativistic theory. [ABSTRACT FROM AUTHOR]

A new algorithm for the two-component Normalized Elimination of the Small Component (2cNESC) method is presented and tested in the calculation of spin-orbit (SO) splittings for a series of heavy atoms and their molecules. The 2cNESC is a Dirac-exact method that employs the exact two-component one-electron Hamiltonian and thus leads to exact Dirac SO splittings for one-electron atoms. For many-electron atoms and molecules, the effect of the two-electron SO interaction is modeled by a screened nucleus potential using effective nuclear charges as proposed by Boettger [Phys. Rev. B 62, 7809 (2000)]. The use of the screened nucleus potential for the two-electron SO interaction leads to accurate spinor energy splittings, for which the deviations from the accurate Dirac Fock-Coulomb values are on the average far below the deviations observed for other effective one-electron SO operators. For hydrogen halides HX (X = F, Cl, Br, I, At, and Uus) and mercury dihalides HgX2 (X = F, Cl, Br, I) trends in spinor energies and SO splittings as obtained with the 2cNESC method are analyzed and discussed on the basis of coupling schemes and the electronegativity of X. [ABSTRACT FROM AUTHOR]

The dependence of the nuclear quadrupole coupling constants (NQCC) on the interaction between electrons and a nucleus of finite size is theoretically analyzed. A deviation of the ratio of the NQCCs obtained from two different isotopomers of a molecule from the ratio of the corresponding bare nuclear electric quadrupole moments, known as quadrupole anomaly, is interpreted in terms of the logarithmic derivatives of the electric field gradient at the nuclear site with respect to the nuclear charge radius. Quantum chemical calculations based on a Dirac-exact relativistic methodology suggest that the effect of the changing size of the Au nucleus in different isotopomers can be observed for Au-containing molecules, for which the predicted quadrupole anomaly reaches values of the order of 0.1%. This is experimentally detectable and provides an insight into the charge distribution of non-spherical nuclei. [ABSTRACT FROM AUTHOR]

The use of the mixed reference (MR) reduced density matrix, which combines reduced density matrices of the M S = +1 and -1 triplet-ground states, is proposed in the context of the collinear spin-flip-time-dependent density functional theory (SF-TDDFT) methodology. The time-dependent Kohn-Sham equation with the mixed state is solved by the use of spinor-like open-shell orbitals within the linear response formalism, which enables to generate additional configurations in the realm of TD-DFT. The resulting MR-SF-TDDFT computational scheme has several advantages before the conventional collinear SF-TDDFT. The spin-contamination of the response states of SF-TDDFT is nearly removed. This considerably simplifies the identification of the excited states, especially in the "black-box" type applications, such as the automatic geometry optimization, reaction path following, or molecular dynamics simulations. With the new methodology, the accuracy of the description of the excited states is improved as compared to the collinear SF-TDDFT. Several test examples, which include systems typified by strong non-dynamic correlation, orbital (near) degeneracy, and conical intersections, are given to illustrate the performance of the new method.

The isomer shift for the 23.87 keV M1 resonant transition in the (119)Sn nucleus is calibrated with the help of ab initio calculations. The calibration constant alpha((119)Sn) obtained from Hartree-Fock (HF) calculations (alpha(HF)((119)Sn) = (0.081+/-0.002)a(0)(-3) mm/s) and from second-order Moller-Plesset (MP2) calculations (alpha(MP2)((119)Sn) = (0.091+/-0.002)a(0)(-3) mm/s) are in good agreement with the previously obtained values. The importance of a proper treatment of electron correlation effects is demonstrated on the basis of a statistical analysis of the results of the calibration. The approach used in the calibration is applied to study the (119)Sn isomer shift in CaSnO(3) perovskite under pressure. Comparison with the experimental results for the pressure range of 0-36 GPa shows that the current methodology is capable of describing tiny variations of isomer shift with reasonable accuracy.