A model electronic Hamiltonian to study low-lying electronic states of [Fe(bpy)3]2+ in aqueous solution

A simple model electronic Hamiltonian to describe the potential energy surfaces of several low-lying d-d states of the [Fe(bpy)(3)](2+) complex is developed for use in molecular dynamics (MD) simulation studies. On the basis of a method proposed previously for first-row transition metal ions in aqueous solution, the model Hamiltonian is constructed using density functional theory calculations for the lowest singlet and quintet states. MD simulations are then carried out for the two spin states in aqueous solution in order to examine the performance of the model Hamiltonian. The simulation results indicate that the present model electronic Hamiltonian reasonably describes the potential energy surfaces of the two spin states of the aqueous [Fe(bpy)(3)](2+) system, while retaining sufficient simplicity for application in simulation studies on excited state dynamics.     

Theoretical study of low-lying electronic states of Mn2 using a newly developed relativistic model core potential

We investigated the electronic structure of low-lying electronic states of Mn₂ using a newly developed relativistic model core potential (spdsMCP). Calculations were performed at complete active space self-consistent field (CASSCF) and second-order multiconfiguration quasidegenerate perturbation theory (MCQDPT2) levels. The MCQDPT2 calculations reveal that the 1Σg⁺ state is the ground state. Calculated spectroscopic constants are very similar to the results of recent all-electron calculations and experimental values, indicating that the spdsMCP works well for Mn₂, which requires a highly correlated calculation. The wave functions of low-lying states are also analyzed at the CASSCF level.     

Electronic structure of CeF from frozen-core four-component relativistic multiconfigurational quasidegenerate perturbation theory

We have investigated the ground state and the two lowest excited states of the CeF molecule using four-component relativistic multiconfigurational quasidegenerate perturbation theory calculations, assuming the reduced frozen-core approximation. The ground state is found to be (4f(1))(5d(1))(6s(1)), with Omega = 3.5, where Omega is the total electronic angular momentum around the molecular axis. The lowest excited state with Omega = 4.5 is calculated to be 0.104 eV above the ground state and corresponds to the state experimentally found at 0.087 eV. The second lowest excited state is experimentally found at 0.186 eV above the ground state, with Omega = 3.5 based on ligand field theory calculations. The corresponding state having Omega = 3.5 is calculated to be 0.314 eV above the ground state. Around this state, we also have the state with Omega = 4.5. The spectroscopic constants R(e), omega(e), and nu(1-0) calculated for the ground and first excited states are in almost perfect agreement with the experimental values. The characteristics of the CeF ground state are discussed, making comparison with the LaF(+) and LaF molecules. We denote the d- and f-like polarization functions as d(*) and f(*). The chemical bond of CeF is constructed via {Ce(3.6+)(5p(6)d(*0.3)f(*0.1))F(0.6-)(2p(5.6))}(3+) formation, which causes the three valence electrons to be localized at Ce(3.6+).     

Perturbation method to calculate the interaction potentials and electronic excitation spectra of atoms in He nanodroplets

A method is proposed for the calculation of potential energy curves and related electronic excitation spectra of dopant atoms captured in/on He nanodroplets and is applied to alkali metal atoms. The method requires knowledge of the droplet density distribution at equilibrium (here calculated within a bosonic-He density functional approach) and of a set of valence electron orbitals of the bare dopant atom (here calculated by numeric solution of the Schr?dinger equation in a suitably parametrized model potential). The electron-helium interaction is added as a perturbation, and potential energy curves are obtained by numeric diagonalization of the resulting Hamiltonian as a function of an effective coordinate z(A) (here the distance between the dopant atom and center of mass of the droplet, resulting in a pseudodiatomic potential). Excitation spectra are calculated for Na in the companion paper as the Franck-Condon factors between the v = 0 vibrational state in the ground electronic state and excited states of the pseudodiatomic molecule. They agree well with available experimental data, even for highly excited states where a more traditional approach fails.     

Low-lying excited states of C120 and C151: a multireference perturbation theory study

Excited states of two 7-aminocoumarin derivatives, coumarin 120 (7-amino-4-methylcoumarin) and coumarin 151 (7-amino-4-trifluoromethylcoumarin), were investigated using generalized multiconfigurational quasidegenerate perturbation theory (GMC-QDPT), multiconfigurational quasidegenerate perturbation theory (MC-QDPT) and time-dependent density functional theory (TDDFT) with the B3LYP and CAM-B3LYP functionals. The absorption and fluorescence spectra of C120 and C151 were calculated. We elucidated the characters of the low-lying states of C120 and C151. The absorption spectra calculated with GMC-QDPT and TDDFT B3LYP agreed well with the experimental data, while for the fluorescence spectra, the TDDFT calculations overestimated the fluorescence spectra compared to GMC-QDPT calculations. Utilizing active spaces with large numbers of electrons and orbitals for reference functions, GMC-QDPT showed a better performance than MC-QDPT with a complete active space self-consistent field (CASSCF) reference of active space with smaller number of electrons and orbitals. In our gas phase calculation, we found that the optimized structures for the first excited states have a planar amino group with a CN single bond, while the amino group is pyramidal in the ground state.     

Intramolecular charge-transfer state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution: RISM-SCF study

Intramolecular charge-transfer (ICT) state formation of 4-(N,N-dimethylamino)benzonitrile in acetonitrile solution is studied by the reference interaction site model self-consistent field (RISM-SCF) method. Geometry optimizations are performed for each electronic state in solution with the complete-active-space SCF wave functions. Dynamic electron correlation effects are taken into account by using the multiconfigurational quasidegenerate perturbation theory. Two-dimensional free energy surfaces are constructed as the function of the twisting and wagging angles of the dimethylamino group for the ground and locally excited (LE) states. The calculated absorption and fluorescence energies are in good agreement with experiments. The validity of the twisted ICT (TICT) model is confirmed in explaining the dual fluorescence, and the possibility of the planar ICT model is ruled out. To examine the mechanism of the TICT state formation, a "crossing" seam between the LE and charge-transfer (CT) state surfaces is determined. The inversion of two electronic states occurs at a relatively small twisting angle. The effect of solvent reorganization is also examined. It is concluded that the intramolecular twisting coordinate is more important than the solvent fluctuation for the TICT state formation, because the energy difference between the two states is minimally dependent on the solvent configuration.     

An Anderson impurity model for efficient sampling of adiabatic potential energy surfaces of transition metal complexes

We present a model intended for rapid sampling of ground and excited state potential energy surfaces for first-row transition metal active sites. The method is computationally inexpensive and is suited for dynamics simulations where (1) adiabatic states are required "on-the-fly" and (2) the primary source of the electronic coupling between the diabatic states is the perturbative spin-orbit interaction among the 3d electrons. The model Hamiltonian we develop is a variant of the Anderson impurity model and achieves efficiency through a physically motivated basis set reduction based on the large value of the d-d Coulomb interaction U(d) and a Lanczos matrix diagonalization routine to solve for eigenvalues. The model parameters are constrained by fits to the partial density of states obtained from ab initio density functional theory calculations. For a particular application of our model we focus on electron transfer occurring between cobalt ions solvated by ammonium, incorporating configuration interaction between multiplet states for both metal ions. We demonstrate the capability of the method to efficiently calculate adiabatic potential energy surfaces and the electronic coupling factor we have calculated compares well to previous calculations and experiment. (     

CASSCF and CASPT2 studies on the structures, transition energies, and dipole moments of ground and excited states for azulene

The electronic structure of azulene molecule has been studied. We have obtained the optimized structures of ground and singlet excited states by using the complete active space self-consistent-field (CASSCF) method, and calculated vertical and 0-0 transition energies between the ground and excited states with second-order M?ller-Plesset perturbation theory (CASPT2). The CASPT2 calculations indicate that the bond-equalized C(2v) structure is more stable than the bond-alternating C(s) structure in the ground state. For a physical understanding of electronic structure change from C(2v) to C(s), we have performed the CASSCF calculations of Duschinsky matrix describing mixing of the b(2) vibrational mode between the ground (1A(1)) and the first excited (1B(2)) states based on the Kekule-crossing model. The CASPT2 0-0 transition energies are in fairly good agreement with experimental results within 0.1-0.3 eV. The CASSCF oscillator strengths between the ground and excited states are calculated and compared with experimental data. Furthermore, we have calculated the CASPT2 dipole moments of ground and excited states, which show good agreement with experimental values.     

Nonadiabatic coupling and diabatic electronic population dynamics on 11A2 and 11B1 states of ozone molecule

In this work, we examine nonadiabatic population dynamics for 1¹B₁ and 1¹A₂ states of ozone molecule (O₃). In O₃, two lowest singlet excited states, ¹A₂ and ¹B₁, can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O₃, bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O₃, respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1¹B₁ and 1¹A₂ excited electronic states for the first 500 fs after photoexcitation.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.     

Electronic structure of LaF+ and LaF from frozen-core four-component relativistic multiconfigurational quasidegenerate perturbation theory

The electronic structure of the molecules LaF+ and LaF was studied using frozen-core four-component multiconfigurational quasidegenerate perturbation theory. To obtain proper excitation energies for LaF+, it was essential to include electronic correlations between the outermost valence electrons (4f, 5d, and 6s) and ionic core electrons composed of (4s, 4p, 4d, 5s, and 5p). The lowest-lying 16 excited states were examined for LaF+, and the lowest 30 states were examined for LaF. The excitation energies calculated for LaF+ agree with the available experimental values, as well as with values from ligand field theory. Errors are within 0.4 eV; for example, the highest observed state 2Pi is 3.77 eV above the ground state, and the present value is 4.09 eV. For LaF, agreement between the experimental and theoretical state assignments and between the experimental and calculated excitation energies was generally good, except for the electron configurations of certain states. Errors are within 0.4 eV except for a single anomaly; for example, the highest observed excited-state discussed in this work is 2.80 eV above the ground state, and the present value is 2.42 eV. We discuss the characteristics of the bonding in LaF+ and LaF.     

SPOCK.CI: a multireference spin-orbit configuration interaction method for large molecules

We present SPOCK.CI, a selecting direct multireference spin-orbit configuration interaction (MRSOCI) program based on configuration state functions. It constitutes an extension of the spin-free density functional theory/multireference configuration interaction (DFT/MRCI) code by Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] and includes spin-orbit interaction on the same footing with electron correlation. Key features of SPOCK.CI are a fast determination of coupling coefficients between configuration state functions, the use of a nonempirical effective one-electron spin-orbit atomic mean-field Hamiltonian, the application of a resolution-of-the-identity approximation to computationally expensive spin-free four-index integrals, and the use of an efficient multiroot Davidson diagonalization scheme for the complex Hamiltonian matrix. SPOCK.CI can be run either in ab initio mode or as semiempirical procedure combined with density functional theory (DFT/MRSOCI). The application of these techniques and approximations makes it possible to compute spin-dependent properties of large molecules in ground and electronically excited states efficiently and with high confidence. Second-order properties such as phosphorescence rates are known to converge very slowly when evaluated perturbationally by sum-over-state approaches. We have investigated the performance of SPOCK.CI on these properties in three case studies on 4H-pyran-4-thione, dithiosuccinimide, and free-base porphin. In particular, we have studied the dependence of the computed phosphorescence lifetimes on various technical parameters of the MRSOCI wave function such as the size of the configuration space, selection of single excitations, diagonalization thresholds, etc. The results are compared to the outcome of extensive quasidegenerate perturbation theory (QDPT) calculations as well as experiment. In all three cases, the MRSOCI approach is found to be superior to the QDPT expansion and yields results in very good agreement with experimental findings. For molecules up to the size of free-base porphin, MRSOCI calculations can easily be run on a single-processor personal computer. Total CPU times for the evaluation of the electronic excitation spectrum and the phosphorescence lifetime of this molecule are below 40 h.     

Evaluating First-Order Molecular Properties of Delocalized Ionic or Excited States in Molecular Aggregates by Renormalized Excitonic Method

In the past few years, the renormalized excitonic model (REM) approach was developed as an efficient low-scaling ab initio excited state method, which assumes the low-lying excited states of the whole system are a linear combination of various single monomer excitations and utilizes the effective Hamiltonian theory to derive their couplings. In this work, we further extend the REM calculations for the evaluations of first-order molecular properties (e.g. charge population and transition dipole moment) of delocalized ionic or excited states in molecular aggregates, through generalizing the effective Hamiltonian theory to effective operator representation. Results from the test calculations for four different kinds of one dimensional (1D) molecular aggregates (ammonia, formaldehyde, ethylene and pyrrole) indicate that our new scheme can efficiently describe not only the energies but also wavefunction properties of the low-lying delocalized electronic states in large systems.     

A Fock space coupled cluster study on the electronic structure of the UO(2), UO(2) (+), U(4+), and U(5+) species

The ground and excited states of the UO(2) molecule have been studied using a Dirac-Coulomb intermediate Hamiltonian Fock-space coupled cluster approach (DC-IHFSCC). This method is unique in describing dynamic and nondynamic correlation energies at relatively low computational cost. Spin-orbit coupling effects have been fully included by utilizing the four-component Dirac-Coulomb Hamiltonian from the outset. Complementary calculations on the ionized systems UO(2) (+) and UO(2) (2+) as well as on the ions U(4+) and U(5+) were performed to assess the accuracy of this method. The latter calculations improve upon previously published theoretical work. Our calculations confirm the assignment of the ground state of the UO(2) molecule as a (3)Phi(2u) state that arises from the 5f(1)7s(1) configuration. The first state from the 5f(2) configuration is found above 10,000 cm(-1), whereas the first state from the 5f(1)6d(1) configuration is found at 5,047 cm(-1).     

A study of the ground state of manganese dimer using quasidegenerate perturbation theory

We study the electronic structure of the ground state of the manganese dimer using the state-averaged complete active space self-consistent field method, followed by second-order quasidegenerate perturbation theory. Overall potential energy curves are calculated for the 1Sigmag+, 11Sigmau+, and 11Piu states, which are candidates for the ground state. Of these states, the 1Sigmag+ state has the lowest energy and we therefore identify it as the ground state. We find values of 3.29 A, 0.14 eV, and 53.46 cm(-1) for the bond length, dissociation energy, and vibrational frequency, in good agreement with the observed values of 3.4 A, 0.1 eV, and 68.1 cm(-1) in rare-gas matrices. These values show that the manganese dimer is a van der Waals molecule with antiferromagnetic coupling.     

Interpreting ultrafast molecular fragmentation dynamics with ab initio electronic structure calculations

High-level ab initio electronic structure calculations are used to interpret the fragmentation dynamics of CHBr(2)COCF(3), following excitation with an intense ultrafast laser pulse. The potential energy surfaces of the ground and excited cationic states along the dissociative C-CF(3) bond have been calculated using multireference second order perturbation theory methods. The calculations confirm the existence of a charge transfer resonance during the evolution of a dissociative wave packet on the ground state potential energy surface of the molecular cation and yield a detailed picture of the dissociation dynamics observed in earlier work. Comparisons of the ionic spectrum for two similar molecules support a general picture in which molecules are influenced by dynamic resonances in the cation during dissociation.     

The Ni + O2 reaction: a combined IR matrix isolation and theoretical study of the formation and structure of NiO2

The reaction of Ni atoms with molecular oxygen has been reinvestigated experimentally in neon matrices and theoretically at the DFT PW91PW91/6311G(3df) level. Experimental results show that i) the nature of the ground electronic state of the superoxide metastable product is the same in neon and argon matrices, ii) two different photochemical pathways exist for the conversion of the superoxide to the dioxide ground state (involving 1.6 or 4 eV photons) and iii) an important matrix effect exists in the Ni + O(2)--> Ni(O(2)) or ONiO branching ratios. Theoretical results confirm that the electronic ground state of the metastable superoxide corresponds to the singlet state, in agreement with former CCSD(T) calculations, but in contradiction with other recent works. Our results show that the ground electronic state of the dioxide is (1)Sigma(+)(g) with the lowest triplet and quintet states at slightly higher energy, consistent with the observation of weak vibronic transitions in the near infrared. The potential energy profiles are modelled for the ground state and nine electronic excited states and a pathway for the Ni(triplet) + O(2)(triplet) --> Ni(O(2)) or ONiO (singlet) reaction is proposed, as well as for the Ni(O(2)) --> ONiO photochemical reaction, accounting for the experimental observations.