Experimental and theoretical UV characterizations of acetylacetone and its isomers

Cryogenic matrix isolation experiments have allowed the measurement of the UV absorption spectra of the high-energy non-chelated isomers of acetylacetone, these isomers being produced by UV irradiation of the stable chelated form. Their identification has been done by coupling selective UV-induced isomerization, infrared spectroscopy, and harmonic vibrational frequency calculations using density functional theory. The relative energies of the chelated and non-chelated forms of acetylacetone in the S0 state have been obtained using density functional theory and coupled-cluster methods. For each isomer of acetylacetone, we have calculated the UV transition energies and dipole oscillator strengths using the excited-state coupled-cluster methods, including EOMCCSD (equation-of-motion coupled-cluster method with singles and doubles) and CR-EOMCCSD(T) (the completely renormalized EOMCC approach with singles, doubles, and non-iterative triples). For dipole-allowed transition energies, there is a very good agreement between experiment and theory. In particular, the CR-EOMCCSD(T) approach explains the blue shift in the electronic spectrum due to the formation of the non-chelated species after the UV irradiation of the chelated form of acetylacetone. Both experiment and CR-EOMCCSD(T) theory identify two among the seven non-chelated forms to be characterized by red-shifted UV transitions relative to the remaining five non-chelated isomers.     

Open-shell localized Hartree-Fock method based on the generalized adiabatic connection Kohn-Sham formalism for a self-consistent treatment of excited states

An effective exact-exchange Kohn-Sham approach for the treatment of excited electronic states, the generalized adiabatic connection open-shell localized Hartree-Fock (GAC-OSLHF) method is presented. The GAC-OSLHF method is based on the generalized adiabatic connection Kohn-Sham formalism and therefore capable of treating excited electronic states, which are not the energetically lowest of their symmetry. The method is self-interaction free and allows for a fully self-consistent computation of excited valence as well as Rydberg states. Results for atoms and small- and medium-size molecules are presented and compared to restricted open-shell Hartree-Fock (ROHF) and time-dependent density-functional results as well as to experimental data. While GAC-OSLHF and ROHF results are quite close to each other, the GAC-OSLHF method shows a much better convergence behavior. Moreover, the GAC-OSLHF method as a Kohn-Sham method, in contrast to the ROHF approach, represents a framework which allows also for a treatment of correlation besides an exchange by appropriate functionals. In contrast to the common time-dependent density-functional methods, the GAC-OSLHF approach is capable of treating doubly or multiply excited states and can be easily applied to molecules with an open-shell ground state. On the nodal planes of the energetically highest occupied orbital, the local multiplicative GAC-OSLHF exchange potential asymptotically approaches a different, i.e., nonzero, value than in other regions, an asymptotic behavior which is known from exact Kohn-Sham exchange potentials of ground states of molecules.     

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.     

A coupled cluster approach with a hybrid treatment of connected triple excitations: implementation and applications for open-shell systems

An implementation of the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples [denoted as CCSD(T)-h], based on the unrestricted Hartree-Fock (UHF) reference, is presented. Based on the spin-integrated formulation, we have developed a computer program to achieve the automatic derivation and implementation of the CCSD(T)-h approach. The CCSD(T)-h approach computationally scales as the seventh power of the system size, and is affordable for many medium-sized systems. The present approach has been applied to study the equilibrium geometries and harmonic vibrational frequencies in a number of open-shell diatomic molecules and bond breaking potential energy profiles in several open-shell molecules, including CH(3), NH(2), and SiH(2). For all systems under study, the overall performance of the UHF-based CCSD(T)-h approach is very close to that of the corresponding CCSDT (CC singles, doubles, and triples), and much better than that of the UHF-based CCSD(T) (CC singles, doubles, and perturbative triples).     

Spin-adapted open-shell time-dependent density functional theory. III. An even better and simpler formulation

The recently proposed spin-adapted time-dependent density functional theory (S-TD-DFT) [Z. Li and W. Liu, J. Chem. Phys. 133, 064106 (2010)] resolves the spin-contamination problem in describing singly excited states of high spin open-shell systems. It is an extension of the standard restricted open-shell Kohn-Sham-based TD-DFT which can only access those excited states due to singlet-coupled single excitations. It is also far superior over the unrestricted Kohn-Sham-based TD-DFT (U-TD-DFT) which suffers from severe spin contamination for those excited states due to triplet-coupled single excitations. Nonetheless, the accuracy of S-TD-DFT for high spin open-shell systems is still inferior to TD-DFT for well-behaved closed-shell systems. The reason can be traced back to the violation of the spin degeneracy conditions (SDC) by approximate exchange-correlation (XC) functionals. Noticing that spin-adapted random phase approximation (S-RPA) can indeed maintain the SDC by virtue of the Wigner-Eckart theorem, a hybrid ansatz combining the good of S-TD-DFT and S-RPA can immediately be envisaged. The resulting formalism, dubbed as X-TD-DFT, is free of spin contamination and can also be viewed as a S-RPA correction to the XC kernel of U-TD-DFT. Compared with S-TD-DFT, X-TD-DFT leads to much improved results for the low-lying excited states of, e.g., N(2)(+), yet with much reduced computational cost. Therefore, X-TD-DFT can be recommended for routine calculations of excited states of high spin open-shell systems.     

SCF?CI and SCF open-shell studies of the base components of the nucleic acids

The π-electron structure of adenine, guanine, cytosine, thymine and uracil in their ground, ionized, singlet and triplet excited states are investigated by means of the SCF ? CI and SCF open-shell methods. The calculations for singlets fit the maxima of the absorption bands well. The energy difference between the first and the second singlet states of adenine is found to be very small. The open-shell method leads to the same relative ionization potential as does the SCF (with the integrals empirically corrected). The calculated energies of the triplet states almost coincide in the SCF open-shell and the SCF ? CI approximation. The calculated transition energies to the first triplet state of the pyrimidines are higher than in the case of the purines. The value of the singlet–triplet separation energy of purines is in agreement with experimental data.     

Spin-Adapted restricted Hartree–Fock reference coupled-cluster theory for open-shell systems: Noniterative triples for noncanonical orbitals

The coupled clusters singles and doubles (CCSD ) method for calculations of open-shell systems with the single restricted Hartree–Fock (ROHF ) reference determinant is extended by the noniterative triples to give CCSD(T) . Our approach profits from the fact that (a) single- and double-excitation amplitudes are spin-adapted, which directly leads to a computationally less demanding algorithm than are nonadapted procedures and produces the spin-adapted CCSD wave function and (b) triple excitations calculated from converged spin-adapted (SA ) CCSD amplitudes are also obtained more effectively. Altogether, computer demands of our SA CCSD(T) approach, applicable to high-spin open-shell cases which are well represented by a single-determinant reference is comparable to that for closed-shell systems. Our approach is not based on semicanonical orbitals, applied by Bartlett's group. However, we compare some other possible choices of ROHF orbitals to this “standard.” Numerical results for a series of atoms and molecules demonstrate little sensitivity to this selection. © 1995 John Wiley & Sons, Inc.     

Two-hole-one-particle configuration interaction approach for the ionization of open-shell molecules: Application to NO2

A special configuration interaction approach for the ionic states of open-shell molecules with a spatially nondegenerate ground state is presented. This approach takes into account all genuine one-hole and two-hole-one-particle excitations with respect to the ground state and allows for the determination of ionic energies and relative intensities for the whole energy range. A numerical application for NO₂ is reported providing in particular theoretical access to the inner-valence region where the interpretation of experimental data hitherto has been open to question.     

Analytic UHF-CCSD(T) second derivatives: implementation and application to the calculation of the vibration-rotation interaction constants of NCO and NCS

An implementation of analytic open-shell UHF-CCSD(T) second derivatives is presented. To demonstrate applicability and test the accuracy of the UHF-CCSD(T) approach for the determination of spectroscopical parameters, vibration-rotation interaction constants were calculated for the ground (1²Π) and first electronically excited (1²Σ) states of the NCO and NCS radicals. In addition, harmonic vibrational frequencies for both states, the Renner-Teller parameter for the ground state, as well as the 1²Π→1²Σ excitation energy are reported. While the computed values are in good agreement with reliable experimental information for NCO, most of the data presented for NCS are predictions of quantities not well known from experiment. Received: 4 May 1998 / Accepted: 27 July 1998 / Published online: 12 October 1998     

Extension of the renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed hamiltonian to open-shell systems: a benchmark study

The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; W?och, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; W?och, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the F2+ ion and singlet-triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree--Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H --> C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to show that the dependence of the ROHF-based CR-CC(2,3) energies on the method of canonicalization of the ROHF orbitals is, for all practical purposes, negligible.     

Generalization of Dewar's half-electron method for calculating energies of open-shell electronic states

Dewar's “half-electron” model for calculating electronic energies of certain open-shell doublet and triplet states is extended so as to be applicable to the lowest-energy open-shell state of any given symmetry and multiplicity.     

Open-shell singlet character of stable derivatives of nonacene, hexacene and teranthene

The electronic ground states of the recently synthesized stable nonacene derivatives (J. Am. Chem. Soc. 2010, 132, 1261) are open-shell singlets with a polyradical nature instead of closed-shell singlets as originally assumed, according to the unrestricted broken spin-symmetry density functional theory (UBS-DFT) computations (at B3LYP/6-31G*). It is the bulky protecting groups, not the transfer from the open-shell singlet to closed-shell singlet ground state, that stabilizes these longest characterized acenes. Similar analyses also confirmed the open-shell singlet character of the hexacene and teranthene derivatives.     

Excited state geometries within time-dependent and restricted open-shell density functional theories

Singlet excited state geometries of a set of medium sized molecules with different characteristic lowest excitations are studied. Geometry optimizations of excited states are performed with two closely related restricted open-shell Kohn–Sham methods and within linear response to time-dependent density functional theory. The results are compared to wave-function based methods. Excitation energies (vertical and adiabatic) calculated from the open-shell methods show systematic errors depending on the type of excitation. However, for all states accessible by the restricted methods a good agreement for the geometries with time-dependent density functional theory and wave-function based methods is found. An analysis of the energy with respect to the mixing angle for the singly occupied orbitals reveals that some states (mostly [n→π^*]) are stable when symmetry constraints are relaxed and others (mostly [π→π^*]) are instable. This has major implications on the applicability of the restricted open-shell methods in molecular dynamics simulations.     

Implementation of a general multireference configuration interaction procedure with analytic gradients in a semiempirical context using the graphical unitary group approach

The graphical unitary group approach has been applied in an efficient implementation of a general multireference configuration interaction (MRCI) method for use with small active molecular orbital spaces in a semiempirical framework. Gradients can be computed analytically for molecular orbitals from a closed-shell or a half-electron open-shell Hartree-Fock calculation. CPU times for single point energy and gradient calculations are reported. The code allows MRCI geometry optimizations of large molecules, as illustrated for the singlet ground state and the four lowest triplet states of fullerene C(76).     

Multireference coupled-cluster study of the symmetry breaking in the C2B radical

The potential energy surfaces (PESs) for both the ground and the excited electronic states of the C(2)B radical are investigated using various multireference (MR) coupled-cluster (CC) approaches. In the ground state case we employ the reduced MR (RMR) CC approach with singles (S) and doubles (D), the RMR CCSD method, as well as its RMR CCSD(T) version corrected for secondary triples, relying on various model spaces and basis sets. The reliability of this approach is also tested against the benchmark full configuration interaction results obtained for a small Dunning-Hay (DH) basis set. The results imply a clear preference for a cyclic structure which, however, breaks the C(2v) symmetry. This symmetry breaking manifests itself strongly at the level of the independent particle model, as represented by the restricted open-shell Hartree-Fock approximation, but the tendency toward symmetry breaking diminishes with the increasing size of the basis set employed as well as with the enhanced account of the correlation effects. It is likely to disappear in the complete basis set limit. The general model space CCSD method is then used to compute vertical excitation energies for a number of excited states as well as the cuts of the PES as the boron atom moves around the C(2) fragment. These results also explain why no symmetry breaking is found when relying on a spin contaminated unrestricted Hartree-Fock reference, as in the UMP2 method.     

Open-shell singlet character of cyclacenes and short zigzag nanotubes

The electronic ground states of [n]cyclacenes, as well as short-zigzag nanotubes, computed at unrestricted broken spin-symmetry density functional theory (UBS-DFT), were found to be open-shell singlets, rather than triplets. Computations for [6]cyclacene at complete active space self-consistent field (CASSCF) and multireference perturbation theory (MRMP2) levels support this conclusion. Along with strain, the radical character of the open-shell singlet with antiferromagnetically coupled electron spins may contribute to the difficulties in synthesizing [n]cyclacenes.     

Delayed Doublet Emission in a Cerium(III) Complex

Doublet emission from open-shell molecules has demonstrated its research and application value in recent years. However, understandings of the photoluminescence mechanism of open-shell molecules are far less than that of closed-shell molecules, leading to challenges in molecular design of efficient doublet emission systems. Here we report a cerium(III) 4-(9H-carbozol-9-yl)phenyl-tris(pyrazolyl)borate complex Ce(CzPhTp)₃ with a new luminescence mechanism of delayed doublet emission, which also represents the first example with metal-centered delayed photoluminescence. The energy gap between the doublet and triplet excited states of Ce(CzPhTp)₃ is reduced by the management of the inner and outer coordination spheres, thereby promoting efficient energy transfer between the two excited states and activating the delayed emission. The photoluminescence mechanism discovered may provide a new way for the design of efficient doublet emission and bring insights into rational molecular design and energy level regulation in open-shell molecules.