Search

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 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]

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]

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]

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]