Ab initio configuration-interaction study of the ground and low-lying electronic states of NiI

For the first time, we have studied the potential-energy curves, spectroscopic terms, vibrational levels, and the spectroscopic constants of the ground and low-lying excited states of NiI by employing the complete active space self-consistent-field method with relativistic effective core potentials followed by multireference configuration-interaction calculations. We have identified six low-lying electronic states of NiI with doublet spin multiplicities, including three states of Delta symmetry and three states of Pi symmetry of the molecule within 15 000 cm(-1). The lowest (2)Delta state is identified as the ground state of NiI, and the lowest (2)Pi state is found at 2174.56 cm(-1) above it. These results fully support the previous conclusion of the observed spectra although our computational energy separation of the two states is obviously larger than that of the experimental values. The present calculations show that the low-lying excited states [13.9] (2)Pi and [14.6] (2)Delta are 3 (2)Pi and 3 (2)Delta electronic states of NiI, respectively. Our computed spectroscopic terms, vibrational levels, and spectroscopic constants for them are in good agreement with the experimental data available at present. In the present work we have not only suggested assignments for the observed states but also computed more electronic states that are yet to be observed experimentally. (c) 2005 American Institute of Physics.

Quantum electronic transport in a configuration interaction basis

An overview of a many-body approach to calculation of electronic transport in molecular systems is given. The physics required to describe electronic transport through a molecule at the many-body level, without relying on commonly made assumptions such as the Landauer formalism or linear response theory, is discussed. Physically, our method relies on the incorporation of scattering boundary conditions into a many-body wavefunction and application of the maximum entropy principle to the transport region. Mathematically, this simple physical model translates into a constrained nonlinear optimization problem. A strategy for solving the constrained optimization problem is given. (C) 2004 Wiley Periodicals, Inc.

Tools for analysing configuration interaction wavefunctions

The configuration interaction (CI) approach to quantum chemical calculations is a well-established means of calculating accurately the solution to the Schrodinger equation for many-electron systems. It represents the many-body electron wavefunction as a sum of spin-projected Slater determinants of orthogonal one-body spin-orbitals. The CI wavefunction becomes the exact solution of the Schrodinger equation as the length of the expansion becomes infinite, however, it is a difficult quantity to visualise and analyse for many-electron problems. We describe a method for efficiently calculating the spin-averaged one- and two-body reduced density matrices rho(psi)((r) over bar; (r) over bar' ) and Gamma(psi)((r) over bar (1), (r) over bar (2); (r) over bar'(1), (r) over bar'(2)) of an arbitrary CI wavefunction Psi. These low-dimensional functions are helpful tools for analysing many-body wavefunctions; we illustrate this for the case of the electron-electron cusp. From rho and Gamma one can calculate the matrix elements of any one- or two-body spin-free operator (O) over cap. For example, if (O) over cap is an applied electric field, this field can be included into the CI Hamiltonian and polarisation or gating effects may be studied for finite electron systems. (C) 2003 Elsevier B.V. All rights reserved.

Configuration-interaction effects in and intensity dependence of F- multiphoton detachment

Two- and three-photon detachment rates have been obtained for F- using several expansions in the R-matrix Floquet approach. These rates are compared with other theoretical and experimental results. The use of Hartree-Fock wavefunctions for the ground state of F with addition of continuum electrons does not lead to agreement with experiment for two- and three-photon detachment. By adding correlation terms, agreement with experiment and other theoretical results is improved considerably, demonstrating the importance of electron correlation effects. However, convergence with respect to the wavefunction expansion cannot be established, we also study the intensity dependence of multiphoton detachment rates for F- at the Nd-YAG frequency. Due to the ponderomotive shift the three-photon detachment channel closes at an intensity of 8.5 x 10(11) W cm(-2) and the influence of this channel closure on the multiphoton detachment peaks is illustrated by determining the heights of the excess-photon peaks obtained using a Gaussian laser pulse.

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Calculation of converged rovibrational energies and partition function for methane using vibrational-rotational configuration interaction

The rovibration partition function of CH4 was calculated in the temperature range of 100-1000 K using well-converged energy levels that were calculated by vibrational-rotational configuration interaction using the Watson Hamiltonian for total angular momenta J=0-50 and the MULTIMODE computer program. The configuration state functions are products of ground-state occupied and virtual modals obtained using the vibrational self-consistent field method. The Gilbert and Jordan potential energy surface was used for the calculations. The resulting partition function was used to test the harmonic oscillator approximation and the separable-rotation approximation. The harmonic oscillator, rigid-rotator approximation is in error by a factor of 2.3 at 300 K, but we also propose a separable-rotation approximation that is accurate within 2% from 100 to 1000 K. (C) 2004 American Institute of Physics.

Analytical potential energy function and spectroscopic parameters for the ground and excited states of NaH

set of energies at different internuclear distances for the ground electronic state and two excited electronic states of NaH molecule have been calculated using valence internally contracted multireference configuration interaction(MRCI) including Davidson correction and three basis sets. Then, a potential energy curve (PEC) for each state was determined by extrapolating MRCI energies to the complete basis sets limit. Based on the PECs, accurate vibrational energy levels and rotational constants were determined. The computational PECs are were fitted to analytical potential energy functions using the Murrell-Sorbie potential function. Then, accurate spectroscopic parameters were calculated. Compared with experimental results, values obtained with the basis set extrapolation yield a potential energy curve that gives accurate vibrational energy levels, rotational constants and spectroscopic parameters for the NaH molecule. (C) 2004 Elsevier B.V. All rights reserved.

Lowest energy potentials of Li2 + F accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Ab initio potentials of F + Li-2 accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li-2 + F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li-2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Structure and interactions of ultracold Yb ions and Rb atoms

In order to study ultracold charge-transfer processes in hybrid atom-ion traps, we have mapped out the potential-energy curves and molecular parameters for several low-lying states of the Rb, Yb+ system. We employ both a multireference configuration interaction and a full configuration interaction (FCI) approach. Turning points, crossing points, potential minima, and spectroscopic molecular constants are obtained for the lowest five molecular states. Long-range parameters, including the dispersion coefficients, are estimated from our ab initio data. The separated-atom ionization potentials and atomic polarizability of the ytterbium atom (ad=128.4 atomic units) are in good agreement with experiment and previous calculations. We present some dynamical calculations for (adiabatic) scattering lengths for the two lowest (Yb, Rb+) channels that were carried out in our work. However, we find that the pseudopotential approximation is rather limited in validity and only applies to nK temperatures. The adiabatic scattering lengths for both the triplet and singlet channels indicate that both are large and negative in the FCI approximation.

Ab initio prediction of the infrared-absorption spectrum of the C2Cl radical

The three lowest (1(2)A('), 2(2)A('), and 1(2)A(')) potential-energy surfaces of the C2Cl radical, correlating at linear geometries with (2)Sigma(+) and (2)Pi states, have been studied ab initio using a large basis set and multireference configuration-interaction techniques. The electronic ground state is confirmed to be bent with a very low barrier to linearity, due to the strong nonadiabatic electronic interactions taking place in this system. The rovibronic energy levels of the (CCCl)-C-12-C-12-Cl-35 isotopomer and the absolute absorption intensities at a temperature of 5 K have been calculated, to an upper limit of 2000 cm(-1), using diabatic potential-energy and dipole moment surfaces and a recently developed variational method. The resulting vibronic states arise from a strong mixture of all the three electronic components and their assignments are intrinsically ambiguous. (c) 2005 American Institute of Physics.

Electronic structure and spectra of a new molecular species: SI. A theoretical contribution

A high level theoretical approach is used to characterize for the first time a manifold of doublet and quartet A + S and Omega states correlating with the first two dissociation channels of an as yet experimentally unknown molecular species, SI, sulfur monoidide. A set of spectroscopic constants is determined, including vibrationally averaged spin-orbit coupling constants, vibrationally averaged dipole moments, and dissociation energies. The transition dipole moment function for the spin-forbidden transition a (4)Sigma -X (2)Pi, and the associated radiative lifetimes were also evaluated. Two possibilities to detect transitions experimentally and to derive spectroscopic constants are suggested. (C) 2011 Elsevier B. V. All rights reserved.