Ab initio study of the X2Π and A2Σ+ states of OH. I. Potential curves and properties

Optimizing Molecular Geometries in Strong Magnetic Fields

An efficient implementation of geometrical derivatives at the Hartree-Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie-Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Quadrupole, octopole, and hexadecapole electric moments of Sigma, Pi, Delta, and Phi electronic states: cylindrically asymmetric charge density distributions in linear molecules with nonzero electronic angular momentum

The number of independent components, n, of traceless electric 2(l)-multipole moments is determined for C(infinity v) molecules in Sigma(+/-), Pi, Delta, and Phi electronic states (Lambda=0,1,2,3). Each 2(l) pole is defined by a rank-l irreducible tensor with (2l+1) components P(m)((l)) proportional to the solid spherical harmonic r(l)Y(m)(l)(theta,phi). Here we focus our attention on 2(l) poles with l=2,3,4 (quadrupole Theta, octopole Omega, and hexadecapole Phi). An important conclusion of this study is that n can be 1 or 2 depending on both the multipole rank l and state quantum number Lambda. For Sigma(+/-)(Lambda=0) states, all 2(l) poles have one independent parameter (n=1). For spatially degenerate states--Pi, Delta, and Phi (Lambda=1,2,3)--the general rule reads n=1 for l<2/Lambda/ (when the 2(l)-pole rank lies below 2/Lambda/ but n=2 for higher 2(l) poles with l>or=2/Lambda/. The second nonzero term is the off-diagonal matrix element [formula: see text]. Thus, a Pi(Lambda=1) state has one dipole (mu(z)) but two independent 2(l) poles for l>or=2--starting with the quadrupole [Theta(zz),(Theta(xx)-Theta(yy))]. A Delta(Lambda=2) state has n=1 for 2((1,2,3)) poles (mu(z),Theta(zz),Omega(zzz)) but n=2 for higher 2((l>or=4)) poles--from the hexadecapole Phi up. For Phi(Lambda=3) states, it holds that n=1 for 2(1) to 2(5) poles but n=2 for all 2((l>or=6)) poles. In short, what is usually stated in the literature--that n=1 for all possible 2(l) poles of linear molecules--only applies to Sigma(+/-) states. For degenerate states with n=2, all Cartesian 2(l)-pole components (l>or=2/Lambda/) can be expressed as linear combinations of two irreducible multipoles, P(m=0)((l)) and P/m/=2 Lambda)((l)) [parallel (z axis) and anisotropy (xy plane)]. Our predictions are exemplified by the Theta, Omega, and Phi moments calculated for Lambda=0-3 states of selected diatomics (in parentheses): X (2)Sigma(+)(CN), X (2)Pi(NO), a (3)Pi(u)(C(2)), X (2)Delta(NiH), X (3)Delta(TiO), X (3)Phi(CoF), and X (4)Phi(TiF). States of Pi symmetry are most affected by the deviation from axial symmetry.

Calculation of static molecular properties in the framework of the unitary group based coupled cluster approach

The recently developed and implemented state selective, fully spin-adapted coupled cluster (CC) method that employs a single, yet effectively multiconfigurational, spin-free reference and the formalism of the unitary group approach (UGA) to the many-electron correlation problem, has been employed to calculate static electric properties of various open-shell (OS) systems using the finite field (FF) technique. Starting with the lithium atom, the method was applied at the first-order interacting space single and double excitation level (CCSD(is)) to several first- and second-row hydrides having OS ground state, namely to the CH, NH, OH, SiH, PH, and SH radicals. In the case of NH we also considered three OS excited states. In all cases the dipole moment and polarizability were determined using a high quality basis set and compared with the experiment, whenever available, as well as with various configuration interaction results and other theoretical results that are available from the literature. The agreement of our CCSD(is) values with experiment is very satisfactory except for the 3Σ− ground state of the NH radical, where the experimentally determined dipole moment is too small. No experimental data are available for the corresponding polarizabilities. It is also shown that the FF technique is not suitable for calculations of higher order static properties, such as the hyperpolarizability β and γ tensors. For this reason we formulate the linear response version of our UGA-based CCSD approach and discuss the aspects of its future implementation. Key words: static molecular properties, dipole moments, polarizabilities, free radicals, unitary group based coupled cluster method, linear response theory, finite field technique.