The elusive excited states of bithiophene: a CASPT2 detective story

A systematic multi-reference perturbation theory investigation of the excitation energies and oscillator strengths for the lowest excited states of 2,2′-bithiophene unequivocally shows that its optical spectrum is produced by two ¹ B _u states separated from each other by approximately 1 eV. This picture is confirmed by additional calculations with alternative quantum chemical methods. Our findings are in strong contrast with the previous CASPT2 results of Rubio et al. [J Chem Phys 102:3580 (1995) and Chem Phys Chem 4:1308 (2003)], who predicted that the two lowest ¹ B _u states are quasi-degenerate. The methodological reasons responsible for the previous seemingly erroneous assignment of the optical spectrum of bithiophene are identified and explained in terms of unusually large coupling between the ¹ B _u states introduced by dynamical correlation effects. A general discussion of applicable computational techniques is offered aiming at avoiding similar problems for other molecular systems.     

Electronic excitation in anionic polymerization of butadiene: Nonempirical calculations of reaction complexes

A new mechanism of anionic polymerization of butadiene is proposed. In the elementary chemical act, the “living” polymer–monomer complex is excited into the low‐lying triplet state. This state has the character of charge (electron) and cation (Li⁺ or Na⁺) transfer from the terminal unit of the active center to the monomer molecule. In the framework of this concept, the probability of chemical bond formation is determined by spin density on radical centers of reagent molecules. Semiempirical and ab initio 6‐31G** quantum‐chemical calculations showed stable interaction between components of the complex in the ground electronic state (9–11 kcal/mol) and low energy levels of triplet excited states (<14 kcal/mol). This new approach is shown to be useful in the analysis of polymerization kinetics and the microstructure of polybutadiene depending on the cation type and the ion pair state. The mechanism of cis‐trans isomerization in the terminal unit of the living polymer consists in concerted rotation about the C_β? C_γ bond and the migration of Li between C_α and C_γ atoms. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

A Topological Study of the Ferromagnetic “No-Pair Bonding” in Maximum-Spin Lithium clusters: n+1Li n (n=2–6)

The nature of the bonding between lithium atoms, in low-spin and maximum-spin clusters, was investigated using the topological electron localization function (ELF) approach. The maximum-spin clusters are especially intriguing since their bonding is sustained without having even a single electron pair! Hence this type of bonding had been called “no-pair ferromagnetic-bonding” [Danovich, Wu, Shaik J Am Chem Soc 121:3165 (1999); Glokhovtsev, Schleyer Isr J Chem 33: 455 (1993); de Visser, Danovich, Wu, Shaik J Phys Chem A 106:4961 (2002)]. The following conclusions were reached in the study: (a) In the ground state of Li _n , covalent bonding between Li atoms is accounted by the presence of the disynaptic valence basins, which exhibit a significant degree of inter-basin delocalization. (b) Except for the ³Li₂ case, the valence basins of all maximum-spin clusters are populated by unpaired electrons. The valence basins are located off Li–Li axis (or Li–Li–Li plane), so that their spatial distribution minimizes the mutual Pauli repulsion and screens the electrostatic repulsion between the Li cores. The inter-basin delocalization is rather high, thereby indicating that the unpaired electrons are virtually delocalized over all the valence basins. (c) The ELF analysis shows that Li atoms in the low-spin clusters are bonded by “two-center two-electron” and “three-center two-electron” bonds. (d) In the maximum-spin species, bonding is sustained by “two-center one-electron” and “three-center one-electron” bonds. The latter picture is complementary to the valence bond picture [Danovich, Wu, Shaik J Am Chem Soc 121 3165 (1999); de Visser, Danovich, Wu, Shaik J Phys Chem A 106: 4961 (2002)], in which the bicentric ferromagnetic-bonding is delocalized over all the short Li–Li contacts, by the mixing of the ionic structures and other nonredundant structures into the repulsive high-spin covalent structure in which all the electrons populate the 2s atomic orbitals, i.e., the configuration. In such a manner bonding can be sustained from “purely ferromagnetic interactions” without electron pairing.Dedicated to Jean-Paul Malrieu, a friend and a poet-scientist     

The adsorptions of silver-doped small gold clusters toward carbon monoxide molecule

An all-electron scalar relativistic calculation on Au _n AgCO (n = 1–12) clusters has been performed using density functional theory with the generalized gradient approximation at PW91 level. The introduction of impurity silver weakens the adsorption, and, however, promotes the reactivity enhancement of CO molecule. The CO molecule is relatively more favorable to be adsorbed by the odd-numbered Au _n Ag clusters with closed-shell electronic structure. The values of chemical hardness indicate that the Au _n AgCO cluster is less stable than the corresponding Au _n+1CO cluster chemically. This picture of the influence of impurity silver on the adsorption behavior of Au _n Ag (n = 1–12) clusters toward CO molecule is consistent with previous experimental work (Haeck et al. in J Phys Chem A 115:2103, 2011), in which the cluster’s reaction probability toward CO molecule is reduced upon substitution of gold atoms for silver and the clusters with closed electronic shell are the most reactive toward CO molecule.     

Electronic,spectra, and spin orbit interaction for FrAr van der Waals system

The potential energy curves and spectroscopic constants of the ground and many excited states of the FrAr van der Waals system have been determined using a one‐electron pseudopotential approach. The Fr⁺ core and the electron–Ar interactions are replaced by effective potentials. The Fr⁺Ar core–core interaction is incorporated using the accurate CCSD(T) potential of Hickling et al. (Phys. Chem. Chem. Phys. 2004, 6, 4233). This approach reduces the number of active electrons of the FrAr van der Waals system to only one valence electron, which permits the use of very large basis sets for the Fr and Ar atoms. Using this technique, the potential energy curves of the ground and many excited states are calculated at the self consistent field (SCF) level. In addition, the spin–orbit interaction is also considered using the semiempirical scheme for the states dissociating into Fr (7p) and Fr (8p). The FrAr system is not studied previously and its potential interactions, spectroscopic constants and dipole functions are presented here for the first time. Furthermore, we have predicted the X²Σ⁺–A²Π_1/2, X²Σ⁺–AΠ_3/2, X²Σ⁺–B²Σ_1/2⁺, X²Σ⁺–3²Π_1/2, X²Σ⁺–3²Π_3/2, and X²Σ⁺–5²Σ_1/2⁺ absorption spectra. © 2012 Wiley Periodicals, Inc.     

Quantum optimal control strategies for photoisomerization via electronically excited states

Optimal control of the photoisomerization of Li₂Na from the stable acute to the near-degenerate obtuse triangular configuration is simulated by means of representative wave packet dynamics on two ab initio potential energy surfaces for the electronic ground (X²A′) and excited (4²A′) states. Product state specifity is achieved by means of new iteration methods [W. Zhu, J. Botina and H. Rabitz, J. Chem. Phys. 108 (1998) 1953] which incorporate feedback from the control field. An additional restriction yields a smooth switch on and off behaviour of the optimized pulses. A windowed Fourier transform decomposes the optimized laser field into efficient pump–dump pulses     

A localized electrons detector for atomic and molecular systems

The local value of the single-particle momentum provides a direct three-dimensional representation of bonding interactions in molecules. It is given exclusively in terms of the electron density and its gradient, and therefore is an ideal localized electrons detector (LED). The results introduced here extend to molecular systems our study of the single-particle local momentum in atomic systems (Bohórquez and Boyd in J Chem Phys 129:024110, 2008; Chem Phys Lett 480:127, 2009). LED is able to clearly identify covalent and hydrogen bonding interactions by depicting distinctive regions around the bond critical points, emerging as a complementary tool in conventional atoms in molecules studies. The local variable we introduce here is an intuitively interpretable 3D electron-pairs locator in atoms and molecules that can be computed either from theoretical or experimentally derived electron densities.     

Photophysics of Schiff Bases: Theoretical Study of Salicylidene Methylamine

The proton‐transfer reaction in a model aromatic Schiff base, salicylidene methylamine (SMA), in the ground and in the lowest electronically‐excited singlet states, is theoretically analyzed with the aid of second‐order approximate coupled‐cluster model CC2, time‐dependent density functional theory (TD‐DFT) using the Becke, three‐parameter Lee–Yang–Parr (B3LYP) functional, and complete active space perturbation theory CASPT2 electronic structure methods. Computed vertical‐absorption spectra for the stable ground‐state isomers of SMA fully confirm the photochromism of SMA. The potential‐energy profiles of the ground and the lowest excited singlet state are calculated and four photophysically relevant isomeric forms of SMA; α, β, γ, and δ are discussed. The calculations indicate two S₁/S₀ conical intersections which provide non‐adiabatic gates for a radiationless decay to the ground state. The photophysical scheme which emerges from the theoretical study is related to recent experimental results obtained for SMA and its derivatives in the low‐temperature argon matrices (J. Grzegorzek, A. Filarowski, Z. Mielke, Phys. Chem. Chem. Phys. 2011 , 13, 16596–16605). Our results suggest that aromatic Schiff bases are potential candidates for optically driven molecular switches.     

Dispersion coefficients for Li+–H and Li+–He systems with coulomb and screened coulomb potentials

Exploiting powerful computational aspects and highly correlated exponential wave functions for two‐electron atoms, we have investigated the effects of screened Coulomb interaction on the hexadecapole polarizability of Li⁺(1¹S), and the dispersion coefficients C₆, C₈, C₁₀, and C₁₂ for interaction of Li⁺ with H and He atoms in their ground states. The dispersion coefficients and hexadecapole polarizability for different screening parameters ranging from 0 to 1.0 a are reported. In the unscreened case, the hexadecapole polarizability of Li⁺, and the dispersion C₁₂ coefficients for Li⁺–H and Li⁺–He system are reported for the first time in the literature. The C₆, C₈, and C₁₀ coefficients for the unscreened cases are comparable with the reported results. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Applications of n-electron wave functions built up with two kinds of geminals

The ground states of atoms and molecules Li^?, Be, LiH, LiH, and Li₂ have been calculated using the n-electron wave functions built up with two kinds of geminals. As a comparison, the above systems have been calculated with the Hartree–Fock self-consistent field and the multi-configuration self-consistent field method as well. The results show that the wave functions in this work are capable of describing the electron correlations, and both kinds of geminals can be taken as a starting point in building up n-electron ground states.     

Density functional study of AgnPd and AgnPdH clusters

Small Ag_nPd (n = 5) clusters and their hydrides Ag_nPdH (n = 5) have been studied by density functional theory calculations. For bare clusters, the structures in which the Pd atom has a maximum number of neighboring Ag atoms tend to be energetically favorable. Hydrogen prefers binding to Ag? Pd bridge site of Ag_nPd clusters except for Ag₅Pd. The binding energy has a strong odd–even oscillation. The electron transfers are from Ag atoms to Pd in bare clusters and are from metal clusters to H in cluster hydrides. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical study of laser light generation and color image formation: FA1:Cs+ and FA2:Li+ centers at the low coordination (100) and (110) surfaces of AgCl and AgBr

The F_A1:Cs⁺ and F_A2:Li⁺ color centers at the low coordination (100) and (110) surfaces of AgCl and AgBr play important roles in laser light generation and color image formation. Double‐well potentials at these surfaces are investigated by using ab initio calculations. Quantum clusters were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces, and ions that are the nearest neighbors to the F_A ? defect site are allowed to relax to equilibrium. The calculated Stokes shifts suggest that laser light generation is sensitive to the simultaneous effects of the vibrational coupling mode, the impurity cation, the coordination number of the surface ion, the lattice anion, and the choice of the basis set centered on the anion vacancy. An attempt has been made to explain these effects in terms of Madelung potential, electron affinity, and optical–optical conversion efficiency. All relaxed excited states of the defect‐containing surfaces are deep below the lower edges of the conduction bands of the ground‐state defect‐free surfaces, suggesting that the F_A(I):Cs⁺ and F_A(II):Li⁺ centers are suitable laser defects. The dependence of orientational destruction, recording sensitivity, and exciton (energy) transfer on the empty cation; the coordination number of the surface ion; and the lattice anion is clarified. The Glasner–Tompkins empirical rule was generalized to include the impurity cation and the coordination number of the surface ion. As far as color image formation is concerned, the supersensitizer was found to increase the sensitizing capabilities of two primary dyes in the excited states by increasing the relative yield of quantum efficiency. The (110) surfaces of AgBr and AgCl were more sensitive than the corresponding (100) surfaces, and AgBr thin film was found to be more sensitive than that of AgCl. On the basis of quasi‐Fermi levels, the difference in the sensitizing capabilities between the examined dyes in the excited states is determined. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical investigation of the adsorption of oxygen on small bimetallic LimCun (m,n ≤ 4) clusters

A theoretical study of the adsorption of molecular oxygen on small bimetallic Li_mCu_n (m, n ≤ 4) clusters was carried out using density functional methods, and it was compared with the adsorption of O₂ on copper (Cu_n, n ≤ 8) clusters. The study of O₂‐Li_mCu_n system is important to understand the promotion effects of the alkali atoms on the copper surface participating in the catalytic processes. Adsorption energies ranging from 7.9 to 51 kcal/mol were found, which represented values over 30% to those calculated for the adsorption of O₂ on copper clusters in a previous study. Thus, the reactivity of molecular oxygen on bimetallic clusters is more favorable with high tendency being in favor of the dissociation of the O₂ molecule. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Evidence for long-lived excited states of [CnH2]2+ carbodications

The energetics, structures, stabilities and reactivities of[C_nH₂]²⁺ ions have been investigated using computational methods and experimental mass spectrometric techniques. Spontaneous decompositions of [C_nH₂]²⁺ into [C_nH]⁺ + H⁺ products, observed for ions with odd-n values, have been explained by invoking the formation of excited triplet states. Even-n [C_nH]⁺ ions possess triplet ground states with low-lying excited states, whereas odd-n ions have triplet states with energies several eV above ground singlet states. Radiationless transitions of vibrationally excited long-lived triplet state ions into singlet state continua are suggested as possible mechanisms for spontaneous deprotonation processes of odd-n [C_nH₂]²⁺ ions. Evidence for these long-lived excited states has been obtained in bimolecular single electron transfer reactions.     

Dynamics of Electronically Excited Metal Atoms (Zn and Cd) in Rare Gas Matrices: Simulation of Multiple-band Emission using Molecular Dynamics with Quantum Transitions

Molecular dynamics with quantum transitions approach is employed to simulate the spectroscopic characteristics of the ¹P₁↔¹S₀ transitions in atomic zinc and cadmium in order to gain insight into the excited state behavior of these atoms isolated in solid rare gases neon, argon, and krypton. The absorption and emission spectra are simulated. Non-radiative processes play a fundamental role in the transfer of population among the three electronic states initially accessed in absorption. Three distinct relaxation pathways were identified. Two of these are related to the dynamical modes described in previous works [McCaffrey and Kerins, J. Chem. Phys. 106 , 7885 (1997); Kerins and McCaffrey, J. Chem. Phys. 109 , 3131 (1998)] in which the system evolves to form a square planar configuration around the metal atom. The third distinct pathway involves motion on a hexagonal close packed plane. The temperature dependence of complex formation was also determined for the three relaxation pathways.     

Theoretical studies on structures and electronic spectra of linear carbon chains C2nH+ (n = 1−5)

The density functional theory (DFT) and the complete active space self‐consistent‐field (CASSCF) method have been used for full geometry optimization of carbon chains C_2nH⁺ (n = 1–5) in their ground states and selected excited states, respectively. Calculations show that C_2nH⁺ (n = 1–5) have stable linear structures with the ground state of X³Π for C₂H⁺ or X³Σ^? for other species. The excited‐state properties of C_2nH⁺ have been investigated by the multiconfigurational second‐order perturbation theory (CASPT2), and predicted vertical excitation energies show good agreement with the available experimental values. On the basis of our calculations, the unsolved observed bands in previous experiments have been interpreted. CASSCF/CASPT2 calculations also have been used to explore the vertical emission energy of selected low‐lying states in C_2nH⁺ (n = 1–5). Present results indicate that the predicted vertical excitation and emission energies of C_2nH⁺ have similar size dependences, and they gradually decrease as the chain size increases. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Pseudorotational dynamics of H3+ and Li3+

The origin of the pseudoprecession phenomenon is investigated through a computational study of the time evolution of H₃⁺ and Li₃⁺ by electron nuclear dynamics theory. In particular, the pseudorotation of both molecules is shown to induce a spatial rotation, which in turn leads to Coriolis coupling of the two orthogonal nuclear shape deformation modes. This effect is rooted in an anisotropy of the molecular ground state potential energy surface that is caused by the interaction between the D_3h ground state and a twofold degenerate first excited state. Computations are performed for a variety of vibrational energies. In addition, the impact of the anharmonicity of the ground state potential surface on the shape deformation modes and the coupling between them is discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Study of electron impact excitation of atoms through lasers

We discuss three different experiments for studying electron-excitation of atoms where lasers have been used in combination. These are stepwise electron–photon excitation, superelastic electron scattering from laser excited atoms and excitation of atoms using spin polarized electrons produced by lasers. We present distorted wave calculations and compare our results with the recently reported such experimental measurements. In particular, the results for the alignment and orientation of the excited n ²P states of K ($n=4$) and Rb ($n=5$) atoms and the spin parameters for the lowest excited ¹P₁ and ³P_0,1,2 states of argon by polarized electrons are presented and discussed.     

Threshold photoionization behaviour of (Li2O)Lin clusters produced by a laser vaporization source

We report on the production of small and medium size lithium and lithium oxide clusters by a laser vaporization cluster source. The isotopomeric distribution of natural lithium allowed to identify Li_kO clusters as the most abundant components in the mass spectrum. Photoionization efficiency curves of Li_kO clusters with photon energies from 3.4 to 4.7 eV were measured for 8 ≤ k ≤ 27. Using linear extrapolation of the increase in photoionization efficiency with photon energy, ionization potentials were extracted. With the chemical bond of the O^2- anion to two Li atoms, leaving n = k-2 valence electrons in the (Li₂O)Li_n clusters, clear shell closure effects are present at n = 8 and n = 20.