A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

I. Electronic states of heavily-doped molecular crystals--naphthalene. I. Theoretical. II. Experimental. II. Lowest singlet and triplet excited states of pyrazine

PART I

The energy spectrum of heavily-doped molecular crystals was treated in the Green’s function formulation. The mixed crystal Green’s function was obtained by averaging over all possible impurity distributions. The resulting Green’s function, which takes the form of an infinite perturbation expansion, was further approximated by a closed form suitable for numerical calculations. The density-of-states functions and optical spectra for binary mixtures of normal naphthalene and deuterated naphthalene were calculated using the pure crystal density-of-state functions. The results showed that when the trap depth is large, two separate energy bands persist, but when the trap depth is small only a single band exists. Furthermore, in the former case it was found that the intensities of the outer Davydov bands are enhanced whereas the inner bands are weakened. Comparisons with previous theoretical calculations and experimental results are also made.

PART II

The energy states and optical spectra of heavily-doped mixed crystals are investigated. Studies are made for the following binary systems: (1) naphthalene-h₈ and d₈, (2) naphthalene--h₈ and αd₄, and (3) naphthalene--h₈ and βd₁, corresponding to strong, medium and weak perturbations. In addition to ordinary absorption spectra at 4˚K, band-to-band transitions at both 4˚K and 77˚K are also analyzed with emphasis on their relations to cooperative excitation and overall density-of-states functions for mixed crystals. It is found that the theoretical calculations presented in a previous paper agree generally with experiments except for cluster states observed in system (1) at lower guest concentrations. These features are discussed semi-quantitatively. As to the intermolecular interaction parameters, it is found that experimental results compare favorably with calculations based on experimental density-of-states functions but not with those based on octopole interactions or charge-transfer interactions. Previous experimental results of Sheka and the theoretical model of Broude and Rashba are also compared with present investigations.

PART III

The phosphorescence, fluorescence and absorption spectra of pyrazine-h₄ and d₄ have been obtained at 4˚K in a benzene matrix. For comparison, those of the isotopically mixed crystal pyrazine-h₄ in d₄ were also taken. All these spectra show extremely sharp and well-resolved lines and reveal detailed vibronic structure.

The analysis of the weak fluorescence spectrum resolves the long-disputed question of whether one or two transitions are involved in the near-ultraviolet absorption of pyrazine. The “mirror-image relationship” between absorption and emission shows that the lowest singlet state is an allowed transition, properly designated as ¹B_3u ← ¹A_1g. The forbidden component ¹B_2g, predicted by both “exciton” and MO theories to be below the allowed component, must lie higher. Its exact location still remains uncertain.

The phosphorescence spectrum when compared with the excitation phosphorescence spectra, indicates that the lowest triplet state is also symmetry allowed, showing a strong 0-0 band and a “mirror-image relationship” between absorption and emission. In accordance with previous work, the triplet state is designated as ³B_3u.

The vibronic structure of the phosphorescence spectrum is very complicated. Previous work on the analysis of this spectrum all concluded that a long progression of v_6a exists. Under the high resolution attainable in our work, the supposed v_6a progression proves to have a composite triplet structure, starting from the second member of the progression. Not only is the v_9a hydrogen-bending mode present as shown by the appearance of the C-D bending mode in the d₄ spectrum, but a band of 1207 cm^-1 in the pyrazine in benzene system and 1231 cm^-1 in the mixed crystal system is also observed. This band is assigned as 2v_6b and of a_1g symmetry. Its anonymously strong intensity in the phosphorescence spectrum is interpreted as due to the Fermi resonance with the 2v_6a and v_9a band.

To help resolve the present controversy over the crystal phosphorescence spectrum of pyrazine, detailed vibrational analyses of the emission spectra were made. The fluorescence spectrum has essentially the same vibronic structure as the phosphorescence spectrum.

Properties of the correlated electronic states of pyrene and hydrogenated pyrene

Approximate calculations are reported on pyrene within the PPP model Hamiltonian using a novel restricted CI scheme which employs both molecular orbital and valence bond techniques. Also reported are detailed full CI results of the PPP model on 2,7-dihydropyrene obtained using the valence bond method. Spectral studies, charge and spin density calculations in ground and excited states, and ring current calculations in the ground state of the molecules are presented. In pyrene, the calculated excitation energies are in good agreement with experiment. The closed structure pi-conjugated molecule pyrene appears to show smaller distortions from the ground state geometry compared with the open structure pi-conjugated molecule 2,7-dihydropyrene. The ground state equilibrium structure of 2,7-dihydropyrene can be viewed as two hexatriene molecules connected by a vinyl crosslink, as is evident from bond order and ring current calculations. This is consistent with the only Kekule resonant structure possible for this molecule.

Excited-state structures of molecules with nearby electronic states

Vibronic coupling among the nearby excited electronic states via the in-plane and the out-of-plane nuclear motions is examined in benzene, pyrazine, formaldehyde and thioformaldehyde. Results reveal that in benzene the structure distorts via the most active nuclear bending (planar) motion while in the other molecules the structures distort through an out-of-plane bending motion in their respective lowest excited states.

Correlated electronic states of C-60 fragments

We have studied electronic states of various fragments of C-60 within the Pariser-Parr-Pople (PPP) model and have obtained structural, magnetic and spectral properties of these molecules. The fragments studied include corannulene, fluoranthene and pyracylene. Pyracylene is studied using the exact valence bond (VB) approach while fluoranthene and corannulene are studied using a novel restricted CI technique which employs molecular orbitals for constructing the VB functions. The electronic excitations, bond order and ring currents are calculated for these systems. From these studies, the wide range of absorptions in C-60 can be viewed as those localized on pyracylene units or on the corannulene/fluoranthene units. The bond orders and ring currents show the hexagons to be similar to benzene rings. The pentagon-hexagon bonds are also found to be longer than the hexagon-hexagon bonds.

Optical study of electronic states in GaAsN

GaAs1-xNx alloys with small N composition (x<1%) and GaAsN/GaAs quantum wells (QWs) were studied by continuous wave photoluminescence (PL), pulse wave excitaiton PL and time-resolved PL. In the PL spectra an extra transition located at the higher energy side of the commonly reported N-related emissions was observed. By measuring the PL dependence on temperature and excitation power, the new PL peak was identified as a transition of alloy band edge-related recombination in GaAsN and delocalized transition in QWs. The PL dynamics further confirms its intrinsic nature of band edge states rather than N-related bound states.

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. II. Coulomb interaction effects in single conjugated polymer chains

Conjugated polymers have attracted considerable attention in the last few decades due to their potential for optoelectronic applications. A key step that needs optimisation is charge carrier separation following photoexcitation. To understand better the dynamics of the exciton prior to charge separation, we have performed simulations of the formation and dynamics of localised excitations in single conjugated polymer strands. We use a nonadiabatic molecular dynamics method which allows for the coupled evolution of the nuclear degrees of freedom and of multiconfigurational electronic wavefunctions. We show the relaxation of electron-hole pairs to form excitons and oppositely charged polaron pairs and discuss the modifications to the relaxation process predicted by the inclusion of the Coulomb interaction between the carriers. The issue of charge photogeneration in conjugated polymers in dilute solution is also addressed. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600404]

Electronic Structure of the Ground and Low-Lying Electronic States of MoB and MoB(+)

Multiconfigurational second-order perturbation theory (CASSCF//CASPT2) and quadruple-zeta ANO-RCC basis sets were employed to investigate the ground and low-lying electronic states of MoB and MoB(+). Spectroscopic constants, potential energy curves, wavefunctions, Mulliken population analyses, and ionization energies are given. The ground state of MoB is of X(6)Pi symmetry (R(e) = 1.968 angstrom, omega(e) = 664 cm(-1), and mu = 2.7 D), giving rise to a Omega = 7/2 ground state after including spin-orbit coupling. For MoB(+), the ground state is computed to be of X(7)Sigma(+) symmetry (R(e) = 2.224 angstrom, omega(e) = 141 cm(-1), and mu = 1.2 D), with an adiabatic ionization energy of 7.19 eV and a vertical one of 7.53 eV. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3362-3370, 2011

The low-lying electronic states of BeAs: a first principles characterization

The electronic structure and spectroscopic properties of a manifold of states of a new molecular species, BeAs, have been investigated theoretically at the complete active space self-consistent field/multireference single and double excitations configuration interaction (CASSCF/MRSDCI) approach, using the aug-cc-pV5Z-PP basis set for arsenic, which includes a relativistic effective core potential, and the cc-pV5Z set for beryllium. Potential energy curves of five quartet and eight doublet (I > + S) states correlating with the five lowest-lying dissociation limit are constructed. The effect of spin-orbit coupling is also included in the description of the ground state, and of the doublet states correlating with the second dissociation channel. Dipole moment functions and vibrationally averaged dipole moments are also evaluated. The similarities and differences between BeAs, BeP, and BeN are analyzed. Spin-orbit effects are small for the ground state close to the equilibrium distance, but avoided crossings between Omega = 1/2 states, and between Omega = 3/2 states changes significantly the I > + S curves for the lowest-lying doublets.

Excited Electronic States, Transition Probabilities, and Radiative Lifetimes of CAs: A Theoretical Contribution and Challenge to Experimentalists

High-level CASSCF/MRCI calculations with a quintuple-zeta quality basis set are reported by characterizing for the first time a manifold of electronic states of the CAs radical yet to be investigated experimentally. Along with the potential energy curves and the associated spectroscopic constants, the dipole moment functions for selected electronic states as well as the transition dipole moment functions for the most relevant electronic transitions are also presented. Estimates of radiative transition probabilities and lifetimes complement this investigation, which also assesses the effect of spin-orbit interaction on the A (2)Pi state. Whenever pertinent, comparisons of similarities and differences with the isovalent CN and CP radicals are made.

Electronic Structure and Chemical Bonding in the Lowest Electronic States of TcN

Multiconfiguration second-order perturbation theory, with the inclusion of relativistic effects and spin-orbit Coupling, was employed to investigate the nature of the ground and low-lying Lambda-S and Omega states of the TcN molecule. Spectroscopic constants, effective bond order, and potential energy curves for 13 low-lying Lambda-S states and 5 Omega states are given, The computed ground state of TcN is of Omega = 3 symmetry (R(e) = 1.605 angstrom and omega(e) = 1085 cm(-1)), originating mainly from the (3)Delta Lambda-S ground state. This result is contrasted with the nature of the ground state for other VIIB transtion-metal mononitrides, including X(3)Sigma(-) symmetry for MnN and Omega = 0(+) symmetry for ReN, derived also from a X(3)Sigma(-) state.

The low-lying electronic states of BeP: a reliable and accurate quantum mechanical prediction

A very high level of theoretical treatment (complete active space self-consistent field CASSCF/MRCI/aug-cc-pV5Z) was used to characterize the spectroscopic properties of a manifold of quartet and doublet states of the species BeP, as yet experimentally unknown. Potential energy curves for 11 electronic states were obtained, as well as the associated vibrational energy levels, and a whole set of spectroscopic constants. Dipole moment functions and vibrationally averaged dipole moments were also evaluated. Similarities and differences between BeN and BeP were analysed along with the isovalent SiB species. The molecule BeP has a X (4)Sigma(-) ground state, with an equilibrium bond distance of 2.073 angstrom, and a harmonic frequency of 516.2 cm(-1); it is followed closely by the states (2)Pi (R(e) = 2.081 angstrom, omega(e) = 639.6 cm(-1)) and (2)Sigma(-) (R(e) = 2.074 angstrom, omega(e) = 536.5 cm(-1)), at 502 and 1976 cm(-1), respectively. The other quartets investigated, A (4)Pi (R(e) = 1.991 angstrom, omega(e) = 555.3 cm(-1)) and B (4)Sigma(-) (R(e) = 2.758 angstrom, omega(e) = 292.2 cm(-1)) lie at 13 291 and 24 394 cm(-1), respectively. The remaining doublets ((2)Delta, (2)Sigma(+)(2) and (2)Pi(3)) all fall below 28 000 cm(-1). Avoided crossings between the (2)Sigma(+) states and between the (2)Pi states add an extra complexity to this manifold of states.

Electronic Structure and Chemical Bonding in the Ground and Low-Lying Electronic States of Ta(2)

The electronic structure and chemical bonding of the ground and low-lying Lambda - S and Omega states of Ta(2) were investigated at the multiconfiguration second-order perturbation theory (CASSCF//CASPT2) level. The ground state of Ta(2) is computed to be a X(3)Sigma(-)(g) state (R(e) = 2.120 angstrom, omega(e) = 323 cm(-1), and D(e) = 4.65 eV), with two low-lying singlet states close to it (a(1) Sigma(+)(g) : T(e) = 409 cm(-1), R(e) = 2.131 angstrom, and omega(e) = 313 cm(-1); b(1) Gamma(g): T(e) = 1, 038 cm(-1), R(e) = 2.127 angstrom, and omega(e) = 316 cm(-1)). These electronic states are derived from the same electronic configuration: vertical bar 13 sigma(2)(g)14 sigma(2)(g)7 delta(2)(g)13 pi(4)(u)>. The effective bond order of the X(3) Sigma(-)(g) state is 4.52, which indicates that the Ta atoms are bound by a quintuple chemical bond. The a(1) Sigma(+)(g) state interacts strongly with the X(3)Sigma(-)(g) g ground state by a second-order spin-orbit interaction, giving rise to the (1)0(g)(+) (ground state) (dominated by the X(3)Sigma(-)(g) Lambda - S ground state) and (9)0(g)(+) (dominated by the a(1) Sigma(+)(g) Lambda - S state) Omega states. These results are in line with those reported for the group 5B homonuclear transition metal diatomics. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1306-1315, 2011

Low-lying singlet and triplet electronic states of RhB

The low-lying X-1 Sigma(+), a(3)Delta, A(1)Delta, b(3)Sigma(+), B-1 Pi, c(3)Pi, C-1 Phi, D-1 Sigma(+), E-1 Pi, d(1)Phi, and e(3)Pi electronic states of RhB have been investigated at the ab initio level, using the multistate multiconfigurational second-order perturbation (MS-CASPT2) theory, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the eleven electronic states included in this work, only three (the X-1 Sigma(+), D-1 Sigma(+), and E-1 Pi states) have been investigated experimentally. Potential energy curves, spectroscopic constants, dipole moments, binding energies, and chemical bonding aspects are presented for all electronic states.

A CASSCF/MRCI Study of the Low-Lying Electronic States of the BeS Molecule

Accurate potential energy curves, dipole moment functions, dissociation energies, and molecular constants for several low-lying singlet and triplet electronic states of BeS were investigated using the CASSCF/MRCI methodology, and the cc-pV5Z basis set for beryllium, and the aug-cc-pV(5+d) Z set for sulfur. Besides presenting improved results for the three lowest lying states, this study presents the first theoretical characterization of another set of nine excited states so far unknown experimentally. Our results are sufficiently accurate to reliably guide the experimental search and characterization of these states, and also to confirm the experimental assignment of the B(1)Sigma(+) excited state. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1694-1700, 2011