Search

This article shows a very satisfactory performance of the DFT Koopmans theorem analogue was demonstrated with respect to the energy intervals between the electronic levels of 4Me-BODIPY derivatives.

Vibronic interactions in the pyridine radical cation ground state, 2A1, and its lowest excited states, 2A2 and 2B1, are studied theoretically. These states originate from the ionization out of the highest occupied orbitals of pyridine, 7a1 (nσ), 1a2 (π), and 2b1 (π), respectively, and give rise to the lowest two photoelectron maxima. According to our previous high-level ab initio calculations [Trofimov et al., J. Chem. Phys. 146, 244307 (2017)], the 2A2 (π−1) excited state is very close in energy to the 2A1 (nσ−1) ground state, which suggests that these states could be vibronically coupled. Our present calculations confirm that this is indeed the case. Moreover, the next higher excited state, 2B1 (π−1), is also involved in the vibronic interaction with the 2A1 (nσ−1) and 2A2 (π−1) states. The three-state vibronic coupling problem was treated within the framework of a linear vibronic coupling model employing parameters derived from the ionization energies of pyridine computed using the linear response coupled-cluster method accounting for single, double, and triple excitations (CC3). The potential energy surfaces of the 2A1 and 2A2 states intersect in the vicinity of the adiabatic minimum of the 2A2 state, while the surfaces of the 2A2 and 2B1 states intersect near the 2B1 state minimum. The spectrum computed using the multi-configuration time-dependent Hartree (MCTDH) method accounting for 24 normal modes is in good qualitative agreement with the experimental spectrum of pyridine obtained using high-resolution He I photoelectron spectroscopy and allows for some assignment of the observed features.

The third-order algebraic-diagrammatic construction method for studies of electron detachment processes within the electron propagator framework [IP-ADC(3)] was extended to treat the properties of molecular states with a detached electron using the intermediate state representation (ISR) formalism. The second-order ISR(2) equations for the one-particle (transition) density matrix have been derived and implemented as an extension of the IP-(U)ADC(3) method available in the Q-CHEM program. As a first systematic test of the present IP-(U)ADC(3)/ISR(2) method, the dipole moments of various electronic states of closed- and open-shell molecules have been computed and compared to full configuration interaction (FCI) results. The present study employing FCI benchmarks also provides the first rigorous estimates for the accuracy of electron detachment energies obtained using the IP-ADC(3) method.

Cell and gene damage caused by ionizing radiation has been studied for many years. It is accepted that DNA lesions (single- and double-strand breaks, for example) are induced by secondary species such as radicals, ions and the abundant low-energy secondary electrons generated by the primary radiation. Particularly harmful are dense ionization clusters of several ionization processes within a volume typical for the biomolecular system. Here we report the observation of a damage mechanism in the form of a non-local autoionizing process called intermolecular Coulombic decay (ICD). It directly involves DNA constituents or other organic molecules in an aqueous environment. The products are two energetic ions and three reactive secondary electrons that can cause further damage in their vicinity. Hydrogen-bonded complexes that consist of one tetrahydrofuran (THF) molecule—a surrogate of deoxyribose in the DNA backbone—and one water molecule are used as a model system. After electron impact ionization of the water molecule in the inner-valence shell the vacancy is filled by an outer-valence electron. The released energy is transferred across the hydrogen bridge and leads to ionization of the neighbouring THF molecule. This energy transfer from water to THF is faster than the otherwise occurring intermolecular proton transfer. The signature of the ICD reaction is identified in triple-coincidence measurements of both ions and one of the final state electrons. These results could improve the understanding of radiation damage in biological tissue. The authors study intermolecular Coulomb decay that occurs in a sample of THF and water in a reaction microscope employing triple-coincidence measurements of two ions and one electron. They find that ICD is a previously unconsidered effect between water and other organic molecules that are hydrogen-bonded, with ICD outpacing proton transfer.

Valence photoelectron spectra and photoelectron angular distributions of trans-dichloroethene have been measured with vibrational resolution at photon energies between 19 eV and 90 eV. Calculations of photoelectron anisotropy parameters, β, and harmonic vibrational modes help provide initial insight into the molecular structure. The photon energy range encompasses the expected position of the atomic Cl 3p Cooper minimum. A corresponding dip observed here in the anisotropy of certain photoelectron bands permits the identification and characterization of those molecular orbitals that retain a localized atomic Cl character. The adiabatic approximation holds for the X 2Au state photoelectron band, but vibronic coupling was inferred within the A–B–C and the D–E states by noting various failures of the Franck–Condon model, including vibrationally dependent β-parameters. This is further explored using the linear vibronic coupling model with interaction parameters obtained from ab initio calculations. The A/B photoelectron band is appreciably affected by vibronic coupling, owing to the low-lying conical intersection of the A 2Ag and B 2Bu states. The C 2Bg band is also affected, but to a lesser extent. The adiabatic minima of the D 2Au and E 2Ag states are almost degenerate, and the vibronic interaction between these states is considerable. The potential energy surface of the D 2Au state is predicted to have a double-minimum shape with respect to the au deformations of the molecular structure. The irregular vibrational structure of the resulting single photoelectron band reflects the non-adiabatic nuclear dynamics occurring on the two coupled potential energy surfaces above the energy of their conical intersection.

The third-order non-Dyson algebraic-diagrammatic construction approach to the electron propagator [IP-ADC(3)] is extended using the intermediate state representation (ISR) formalism, allowing the wave functions and properties of molecular states with detached electron to be studied. The second-order ISR equations [ISR(2)] for the one-particle (transition) density matrix have been derived and implemented in the Q-CHEM program. The approach is completely general and enables evaluation of arbitrary one-particle operators and interpretation of electron detachment processes in terms of density-based quantities. The IP-ADC(3)/ISR(2) equations were implemented for Ŝz-adapted intermediate states, allowing open-shell molecules to be studied using unrestricted Hartree-Fock references. As a first test for computations of ground state properties, dipole moments of various closed- and open-shell molecules have been computed by means of electron detachment from the corresponding anions. The results are in good agreement with experimental data. The potential of IP-ADC(3)/ISR(2) for the interpretation of photoelectron spectra is demonstrated for the galvinoxyl free radical.

International audience; The influence of vibronic coupling on the outer valence ionic states of cis-dichloroethene has been investigated by recording photoelectron spectra over the excitation range 19–90 eV using plane polarized synchrotron radiation, for two polarization orientations. The photoelectron anisotropy parameters and electronic state branching ratios derived from these spectra have been compared to theoretical predictions obtained with the continuum multiple scattering approach. This comparison shows that the photoionization dynamics of the à 2B2, B̃ 2A1, C̃ 2A2, and D̃ 2B1 states, all of which are formed through the ejection of an electron from a nominally chlorine lone-pair orbital, exhibit distinct evidence of the Cooper minimum associated with the halogen atom. While retaining a high degree of atomic character, these orbital ionizations nevertheless display clear distinctions. Simulations, assuming the validity of the Born-Oppenheimer and the Franck-Condon approximations, of the X̃ 2B1, à 2B2, and D̃ 2B1 state photoelectron bands have allowed some of the vibrational structure observed in the experimental spectra to be assigned. The simulations provide a very satisfactory interpretation for the X̃ 2B1 state band but appear less successful for the à 2B2 and D̃ 2B1 states, with irregularities appearing in both. The B̃ 2A1 and C̃ 2A2 state photoelectron bands exhibit very diffuse and erratic profiles that cannot be reproduced at this level. Photoelectron anisotropy parameters, β, have been evaluated as a function of binding energy across the studied photon energy range. There is a clear step change in the β values of the à 2B2 band at the onset of the perturbed peak intensities, with β evidently adopting the value of the B̃ 2A1 band β. The D̃ 2B1 band β values also display an unexpected vibrational level dependence, contradicting Franck-Condon expectations. These various behaviors are inferred to be a consequence of vibronic coupling in this system

The valence shell electronic structure and photoionisation dynamics of ortho-dichlorobenzene, ortho-bromochlorobenzene and trichlorobenzene have been investigated both experimentally and theoretica ...

The valence shell ionization spectrum of pyridine was studied using the third-order algebraic-diagrammatic construction approximation scheme for the one-particle Green’s function and the outer-valence Green’s function method. The results were used to interpret angle resolved photoelectron spectra recorded with synchrotron radiation in the photon energy range of 17–120 eV. The lowest four states of the pyridine radical cation, namely, 2A2 (1a 2 −1 1a2−1 ), 2A1(7a 1 −1 7a1−1), 2B1(2b 1 −1 2b1−1), and 2B2(5b 2 −1 5b2−1), were studied in detail using various high-level electronic structure calculation methods. The vertical ionization energies were established using the equation-of-motion coupled-cluster approach with single, double, and triple excitations (EOM-IP-CCSDT) and the complete basis set extrapolation technique. Further interpretation of the electronic structure results was accomplished using Dyson orbitals, electron density difference plots, and a second-order perturbation theory treatment for the relaxation energy. Strong orbital relaxation and electron correlation effects were shown to accompany ionization of the 7a1 orbital, which formally represents the nonbonding σ-type nitrogen lone-pair (nσ) orbital. The theoretical work establishes the important roles of the π-system (π-π* excitations) in the screening of the nσ-hole and of the relaxation of the molecular orbitals in the formation of the 7a1(nσ)−1 state. Equilibrium geometric parameters were computed using the MP2 (second-order Moller-Plesset perturbation theory) and CCSD methods, and the harmonic vibrational frequencies were obtained at the MP2 level of theory for the lowest three cation states. The results were used to estimate the adiabatic 0-0 ionization energies, which were then compared to the available experimental and theoretical data. Photoelectron anisotropy parameters and photoionization partial cross sections, derived from the experimental spectra, were compared to predictions obtained with the continuum multiple scattering approach.

The C 1s ionization spectrum of CH3I has been studied both experimentally and theoretically. Synchrotron radiation has been employed to record polarization dependent photoelectron spectra at a photon energy of 614 eV. These spectra encompass the main-line due to the C 1s single-hole state and the peaks associated with the shake-up satellites. Vertical ionization energies and relative photoelectron intensities have been computed using the fourth-order algebraic-diagrammatic construction approximation scheme for the one-particle Green's function and the 6-311++G** basis set. The theoretical spectrum derived from these calculations agrees qualitatively with the experimental results, thereby allowing the principal spectral features to be assigned. According to our calculations, two 2A1 shake-up states of the C 1s-1 σCI → σCI * type with singlet and triplet intermediate coupling of the electron spins (S' = 0, 1) play an important role in the spectrum and contribute significantly to the overall intensity. Both of these states are expected to have dissociative diabatic potential energy surfaces with respect to the C-I separation. Whereas the upper of these states perturbs the manifold of Rydberg states, the lower state forms a band which is characterized by a strongly increased width. Our results indicate that the lowest shake-up peak with significant spectral intensity is due to the pair (S' = 0, 1) of 2E (C 1s-1 I 5p → σCI *) states. We predict that these 2E states acquire photoelectron intensity due to spin-orbit interaction. Such interactions play an important role here due to the involvement of the I 5p orbitals.

A combined experimental and theoretical investigation has been performed to study the valence shell electronic structure and photoionisation dynamics of meta-dichlorobenzene and meta-bromochloroben ...

The valence shell photoelectron spectra of uracil, 1-methyluracil and 6-methyluracil have been studied experimentally and theoretically. Synchrotron radiation has been used to record spectra at photon energies of 40 and 80 eV. Photoelectron angular distributions have been determined and these provide an experimental means of distinguishing between σ- and π-type orbitals. Vertical ionization energies and spectral intensities have been evaluated using the many-body Green’s function approach, thereby enabling theoretical photoelectron spectra to be derived. The calculated spectra display a satisfactory agreement with the experimental data and this has allowed most of the photoelectron bands to be assigned. Two of the outer-valence vertical ionization energies are similar to one another and the vibrational progressions associated with these transitions overlap strongly. Vibronic interaction between these states, induced through the excitation of out-of-plane vibrational modes, may lead to nonadiabatic effects. Preliminary theoretical investigation of this interaction has been performed.

Valence shell photoelectron spectra of the 5-halouracils have been recorded using synchrotron radiation at photon energies of 40 and 80 eV. Photoelectron angular distributions have been determined and these provide an experimental means of deducing the extent to which the molecular orbitals possess predominantly σ, π or halogen np character. In chloro-, bromo- and iodouracil two peaks, associated with the essentially halogen atom lone-pair orbitals, have been observed. The separation between these peaks, and the value of their photoelectron asymmetry parameters at an excitation energy of 40 eV, are compared with results for the corresponding orbitals in the halobenzenes and halothiophenes. Vertical ionisation energies and spectral intensities have been evaluated using the many-body Green’s function approach, thereby enabling theoretical spectra to be derived. Assignments have been proposed for most of the structure observed in the valence shell photoelectron spectra based upon the available experimental and theoretical evidence.

Angle resolved photoelectron spectra of the X̃(2)B1, Ã(2)A2, B̃(2)B2, and C̃(2)B1 states of bromobenzene have been recorded over the excitation range 20.5-94 eV using linearly polarized synchrotron radiation. The photoelectron anisotropy parameters and electronic branching ratios derived from these spectra have been compared to theoretical predictions obtained with the continuum multiple scattering approach. This comparison shows that ionization from the 8b2 orbital and, to a lesser extent, the 4b1 orbital is influenced by the Cooper minimum associated with the bromine atom. The 8b2 and 4b1 orbitals are nominally bromine lone-pairs, but the latter orbital interacts strongly with the π-orbitals in the benzene ring and this leads to a reduced atomic character. Simulations of the X̃(2)B1, B̃(2)B2, and C̃(2)B1 state photoelectron bands have enabled most of the vibrational structures appearing in the experimental spectra to be assigned. Many of the photoelectron peaks exhibit an asymmetric shape with a tail towards low binding energy. This asymmetry has been examined in the simulations of the vibrationally unexcited peak, due mainly to the adiabatic transition, in the X̃(2)B1 state photoelectron band. The simulations show that the asymmetric profile arises from hot-band transitions. The inclusion of transitions originating from thermally populated levels results in a satisfactory agreement between the experimental and simulated peak shapes.

The valence-shell ionization spectrum of bromobenzene, as a representative halogen substituted aromatic, was studied using the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator. This method, also referred to as IP-ADC(3), was implemented as a part of the Q-Chem program and enables large-scale calculations of the ionization spectra, where the computational effort scales as n(5) with respect to the number of molecular orbitals n. The IP-ADC(3) scheme is ideally suited for investigating low-lying ionization transitions, so fresh insight could be gained into the cationic state manifold of bromobenzene. In particular, the present IP-ADC(3) calculations with the cc-pVTZ basis reveal a whole class of low-lying low-intensity two-hole-one-particle (2h-1p) doublet and quartet states, which are relevant to various photoionization processes. The good qualitative agreement between the theoretical spectral profile for the valence-shell ionization transitions generated with the smaller cc-pVDZ basis set and the experimental photoelectron spectrum measured at a photon energy of 80 eV on the PLEIADES beamline at the Soleil synchrotron radiation source allowed all the main features to be assigned. Some theoretical aspects of the ionization energy calculations concerning the use of various approximation schemes and basis sets are discussed.

A polarization propagator method, referred to as algebraic-diagrammatic construction (ADC), is extended to the treatment of static and dynamic response properties of molecules. The recent intermediate state representation (ISR) concept of the ADC theory, giving direct access to excited states wave functions and properties, allows us to derive simple closed-form expressions for linear and higher response functions. The use of the band-Lanczos algorithm is proposed to evaluate efficiently the resolvent type ADC expressions. The performance of the method is tested in computations of static and dynamic polarizabilities of several small molecules at the second-order (ADC(2)) level of the theory. The ADC(2) results are compared with those of full configuration interaction (FCI), coupled cluster (CC), and SOPPA (second-order polarization propagator approximation) treatments.

A new relaxation pathway of optically allowed valence-excitations in molecules bound to metal atoms was investigated. We show that the σnonπ* excitation of HCN in MgnHCN (n = 1,2) clusters decays through the emission of an electron from the metallic component. The calculated lifetimes are of the order of 140 fs for MgHCN and 25 fs for Mg2HCN. These short lifetimes suggest that the relaxation mechanism discussed here should dominate for molecules adsorbed on metal clusters or surfaces. We find that at large distances this decay is dominated by energy transfer.

The photoelectron spectrum of hydrogen selenide has been recorded using synchrotron radiation in the photon energy range 15–110 eV and the inner valence region has been studied in detail for the first time. Green’s function methods have been employed to evaluate the ionisation energies and spectral intensities of all valence states and the results have facilitated an interpretation of the experimental spectra. Two outer valence satellites have been observed at binding energies of ∼17 and ∼18.5 eV. The continuum multiple scattering (CMS-Xα) approach has been used to calculate photoelectron asymmetry parameters and branching ratios for the outer valence orbitals and these have been compared with experimental data. The Se 4p Cooper minimum strongly affects the photoionisation dynamics of the outermost 4b 1 orbital and its influence is evident in both the theoretically predicted and the measured asymmetry parameter for the X ˜ 2 B 1 state. The experimental data for the A ˜ 2 A 1 and B ˜ 2 B 2 state asymmetry parameters also display energy dependent variations in the vicinity of the Cooper minimum. However, for these two states the agreement between the predicted and observed asymmetry parameters is less satisfactory. The present results for H 2 Se have been compared with similar data for the closely related H 2 S molecule.

The complete valence shell photoelectron spectrum of tetrafluoromethane has been recorded using synchrotron radiation and the observed structure has been interpreted using ionisation energies and relative spectral intensities computed using the third-order algebraic-diagrammatic-construction (ADC(3)) scheme for the one-particle Green’s function and the outer valence Green’s function (OVGF) method. The ADC(3) calculations were performed using both the original variant based on the Dyson equation and the recently proposed non-Dyson scheme. The theoretical predictions for the single-hole ionic states due to outer valence shell ionisation agree satisfactorily with the experimental results. Ionisation from the inner valence 2t2 and 3a1 orbitals is strongly influenced by reorganisation effects and the intensity is spread amongst numerous satellites. Highly resolved spectra have been measured for the C ∼ 2 T 2 and the D ∼ 2 A 1 states, and the vibrational structure has been associated primarily with excitation of the ν 1 + ( a 1 ) mode. However, the analyses also provide evidence for excitation of the Jahn-Teller active ν 2 + ( e ) mode in the C ∼ 2 T 2 state, and the ν 3 + ( t 2 ) mode in the D ∼ 2 A 1 state. Excitation of this latter mode can be explained in terms of vibronic coupling between the C ∼ 2 T 2 and the D ∼ 2 A 1 states. Photoelectron angular distributions and branching ratios have been determined and demonstrate that shape resonances affect the outer valence shell photoionisation dynamics. The experimental results are compared with previous theoretical predictions but a consistent interpretation has not been obtained. The major difficulty concerns the uncertainty in the locations of valence shell transitions into the 5a1 and the 5t2 virtual orbitals.

International audience; The excitation spectra and molecular dynamics of furan associated with its low-lying excited singlet states 1A2(3s), 1B2(V), 1A1(V'), and 1B1(3p) are investigated using an ab initio quantum-dynamical approach. The ab initio results of our previous work [J. Chem. Phys. 119, 737 (2003)] on the potential energy surfaces (PES) of these states indicate that they are vibronically coupled with each other and subject to conical intersections. This should give rise to complex nonadiabatic nuclear dynamics. In the present work the dynamical problem is treated using adequate vibronic coupling models accounting for up to four coupled PES and thirteen vibrational degrees of freedom. The calculations were performed using the multiconfiguration time-dependent Hartree method for wave-packet propagation. It is found that in the low-energy region the nuclear dynamics of furan is governed mainly by vibronic coupling of the 1A2(3s) and 1B2(V) states, involving also the 1A1(V') state. These interactions are responsible for the ultrafast internal conversion from the 1B2(V) state, characterized by a transfer of the electronic population to the 1A2(3s) state on a time scale of approximately 25 fs. The calculated photoabsorption spectrum of furan is in good qualitative agreement with experimental data. Some assignments of the measured spectrum are proposed.

The multi-mode and multi-state vibronic interactions in the heterocyclic molecules furan, pyrrole, thiophene and their radical cations are investigated theoretically, employing a linear vibronic coupling scheme. The underlying system parameters are determined from large-scale ab initio computations. Previous time-independent dynamical calculations on the radical cations are extended by wave-packet propagations (using the MCTDH method) confirming the strong nonadiabatic coupling effects. For the singlet excited states of furan and thiophene quantum dynamical calculations are presented which go beyond the two-state approximation frequently applied in the literature. The characteristic spectral structures are well reproduced, especially in the case of furan. The implications of these results on the photochemical reaction dynamics of these species are discussed.

2 (V), 1 A 1 (V8), respectively, at the C 2v ground-state molecular configuration# have been studied using the equation-of-motion coupled-cluster singles and doubles method ~EOM-CCSD!. Full geometry optimizations with subsequent computation of harmonic vibrational frequencies have been performed in order to locate and characterize stationary points on the potential energy surfaces ~PES!. The latter optimization work was enabled by the availability of analytic energy gradient techniques for the EOM-CCSD approach. A major new finding is that both the 1 B2(V) and 1 A1(V8) valence states are unstable with respect to non-totally symmetric distortions at the C2v configuration. The symmetry breaking in the 1 B2(V) state involves an in-plane coordinate of b2 symmetry. The relaxation process begins on the S2 adiabatic PES and, after passing through a conical intersection of the S2 and S1 PES, continues on the S1 surface, taking the system finally to the adiabatic minimum ofS1 ( 1 A2 state!. The 1 A1(V8) valence state is found to be unstable with respect to the out-of-plane bending coordinates of b1 and a2 symmetry. The resulting relaxed molecular structures have Cs and C2 symmetry, respectively. The present findings are analyzed in terms of a linear vibronic coupling model and spectroscopic implications are discussed. © 2003 American Institute of Physics. @DOI: 10.1063/1.1578051#

A recently developed consistent third-order propagator method for the treatment of electronic excitation in molecules is tested in first applications. The method referred to as third-order algebraic-diagrammatic construction [ADC(3)] extends the existing second-order approximation and aims at a more accurate computation of excitation energies and transition moments than afforded at the second-order level. For a stringent test of the method we compare the ADC(3) energies for over 40 singlet and triplet vertical transitions in H2O, HF, N2, and Ne with the results of recent full configuration interaction (FCI) and coupled cluster (CC) computations. The ADC(3) results reflect a substantial and uniform improvement with respect to the second-order description. The mean absolute deviation of the single excitation energies from the FCI results is below 0.2 eV. Although this does not equal the accuracy of the third-order CC3 model, the ADC(3) method, scaling as N6 with the number of orbitals, may be viewed as a good compromise between accuracy and computational cost.

The recently introduced synchrotron radiation-based Fourier transform spectroscopy has been employed to study the excited electronic states of thiophene. A highly resolved photoabsorption spectrum has been measured between ∼5 and 12.5 eV, providing a wealth of new data. High-level ab initio computations have been performed using the second-order algebraic-diagrammatic construction (ADC(2)) polarization propagator approach, and the equation-of-motion coupled-cluster (EOM-CC) method at the CCSD and CC3 levels, to guide the assignment of the spectrum. The adiabatic energy corrections have been evaluated, thereby extending the theoretical study beyond the vertical excitation picture and leading to a significantly improved understanding of the spectrum. The low-lying π→π* and π→σ* transitions result in prominent broad absorption bands. Two strong Rydberg series converging onto the X(~)(2)A2 state limit have been assigned to the 1a2→npb1(1)B2 and the 1a2→nda2(1)A1 transitions. A second, and much weaker, d-type series has been assigned to the 1a2→ndb1(1)B2 transitions. Excitation into some of the Rydberg states belonging to the two strong series gives rise to vibrational structure, most of which has been interpreted in terms of excitations of the totally symmetric ν4 and ν8 modes. One Rydberg series, assigned to the 3b1→nsa1(1)B1 transitions, has been identified converging onto the Ã(2)B1 state limit, and at higher energies Rydberg states converging onto the B(~)(2)A1 state limit could be identified. The present spectra reveal highly irregular vibrational structure in certain low energy absorption bands, and thus provide a new source of information for the rapidly developing studies of excited state non-adiabatic dynamics and photochemistry.

The valence shell photoelectron spectra of thiophene, 2-chlorothiophene (2-Cl-Th) and 3-chlorothiophene (3-Cl-Th) have been investigated theoretically and experimentally. The third-order algebraic-diagrammatic construction approximation scheme for the one-particle Green’s function has been employed to evaluate the vertical ionisation energies and spectral intensities. The ground state geometrical parameters of the three molecules have been optimised at the level of the second order Moller–Plesset perturbation theory, and standard cc-pvDZ basis sets have been used throughout. The results for the outer valence region of thiophene agree well with available experimental and theoretical data. Very satisfactory agreements have been obtained between the theoretical predictions for the photoelectron spectra of 2-Cl-Th and 3-Cl-Th and the corresponding experimental data recorded in the present study. Assignments have been proposed for the major spectral structures. In all three molecules the breakdown of the molecular orbital ionisation picture associated with the π1 molecular orbital is discussed. Synchrotron radiation has been used to record the complete valence shell photoelectron spectra of 2-Cl-Th and 3-Cl-Th. Photoelectron angular distributions and branching ratios have been determined in the photon energy range 13–115 eV, and those for the chlorine and sulphur 3p orbitals show effects that can be attributed to Cooper minima. A high resolution photoelectron spectrum of the outer valence shell of 2-Cl-Th has been recorded using HeI radiation, and vibrational structure has been observed and analysed in the X 2 A ′′ , A 2 A ′′ , B 2 A ′ , C 2 A ′′ and I 2 A ′ state bands.

The group of three lowest singlet C 1s-excited states of formaldehyde H2CO is studied theoretically. The equilibrium geometries are determined at the restricted open-shell Hartree–Fock (ROHF) level and refined total energies are obtained using the multireference configuration interaction (MRCI) approach. In agreement with an earlier prediction [Chem. Phys. 122, 9 (1988)] the second lowest singlet state, 1B2, is characterized by a doubly excited, “two particle–two hole” (2p–2h), configuration C 1s,n→π*2. Our calculations predict that H2CO in the 1B2(2p–2h) state has a stable pyramidal equilibrium structure with a barrier to inversion of 0.28 eV, the valence angle being close to 107°. The calculated length of the CO bond is 1.390 A. The 1B2(2p–2h) state is shown to be also bound with respect to all possible dissociation and rearrangement processes. The lowest predicted dissociation energy for the 1B2 state (H2CO*→H2+CO* reaction) is 0.29 eV (6.69 kcal/mol). The rationalization of the great stability of the ...

The valence electronic structure and momentum-space electron density distributions of n-butane have been studied by means of high-resolution (e,2e) electron momentum spectroscopy based on noncoplanar symmetric kinematics. Ionization spectra for the range of binding energies 6 to 32 eV and momenta described by azimuthal angles φ=0°, 2°, 4°, 6°, 8°, and 10° have been recorded and compared to the results of one-particle Green’s function calculations, performed using the third-order algebraic–diagrammatic construction [ADC(3)] approximation and series of basis sets of improving quality. Experimental electron momentum profiles have been determined from a set of 11 measurements and compared to theoretical results. It has been shown that despite the complex structure of the spectral bands and the conformational versatility of n-butane, the experimental electron momentum distributions are accurately described by the momentum-space form of orbital densities obtained from Becke three-parameter Lee–Yang–Parr (B3LYP)...

A propagator method referred to as third-order algebraic–diagrammatic construction [ADC(3)] for the direct computation of electronic excitation energies and transition moments is presented. This approach is based on a specific reformulation of the diagrammatic perturbation expansion for the polarization propagator, and extends the existing second-order [ADC(2)] scheme to the next level of perturbation theory. The computational scheme combines diagonalization of a Hermitian secular matrix and perturbation theory for the matrix elements. The characteristic properties of the method are compact configuration spaces, regular perturbation expansions, and size-consistent results. The configuration space is spanned by singly and doubly excited states, while the perturbation expansions in the secular matrix extend through third order in the p-h block, second order in the p-h/2p-2h coupling block, and first order in the 2p-2h block. While the simpler ADC(2) method, representing a counterpart to the MP2 (second-orde...

An efficient third-order propagator method to compute ionization potentials and electron affinities of atoms and molecules is presented. The development is based on the algebraic diagrammatic construction (ADC) representing a specific reformulation of the diagrammatic perturbation series of the electron propagator G(ω). In contrast with previous approximation schemes, relying on the Dyson equation and approximations for the self-energy part, the ADC procedure here is applied directly to the (N∓1)-electron parts G−(ω) and G+(ω), respectively, of the electron propagator. This leads to decoupled secular equations for the ionization energies ((N−1)-electron part) and electron affinities ((N+1)-electron part), respectively. In comparison with the Dyson-type approach, there is a substantial reduction of the secular matrix dimension opposed by a small additional expense in computing some second- and third-order contributions to the secular matrix elements. The relationship of the non-Dyson ADC(3) method to coupl...

The 2A2 and 2B1 states formed in the ionization of the outermost π orbitals in furan, pyrrole and thiophene are shown to interact vibronically via nontotally symmetric b2 vibrational modes. The interaction is strongest in pyrrole and thiophene, where the conical intersection between the two adiabatic surfaces occurs near the minimum of the upper (2B1) state. The resulting nonadiabatic effects manifest themselves in the 2B1 bands by a lack of resolved structure in case of pyrrole and thiophene, and by a line broadening in case of furan. The spectra are investigated using a linear vibronic coupling model. All totally symmetric a1 (tuning) modes and nontotally symmetric b2 (coupling) modes describing the ring motion are taken into account. The parameters of the model are obtained with the aid of ab initio calculations. The ground state optimized geometries and vibrational frequencies are computed at the level of the second-order Mo/ller–Plesset perturbation theory, while the dependence of the ionization ener...

The electronic spectrum of furan is investigated theoretically beyond the previous vertical-electronic description. A polarization propagator method referred to as second-order algebraic-diagrammatic construction (ADC(2)) has been used in the electronic structure calculations. The vibrational excitation accompanying the electronic transitions is described with the aid of a linear electron-vibrational coupling model. The spectral information thereby obtained permits extensive comparison with experiment. The average accuracy of the present method, estimated by comparing adiabatic transition energies, is better than 0.4 eV. Only for the lowest π-π ∗ valance transition, V′( 1 A 1 ) and V′( 1 B 2 ), and for the Rydberg excitations agree The results for the other π-π ∗ valence transitions, V( 1 B 2 ), and for the Rydberg excitations agree well with findings of previous experimental and theoretical work. A (multistate) vibronic coupling effect involving the V′( 1 A 1 ) and V( 1 B 2 ) valence transitions and the 3s( 1 A 2 Rydberg excitation is suggested as the reason for the highly diffuse character of the 5.7–6.7 eV photoabsorption band.

International audience; The potential energy surfaces (PESs) of the two lowest excited singlet states of furan [correlating with the Rydberg A2(3s) and valence B2(V) states at the C2v ground-state molecular configuration] have been studied in some detail with regard to the photoinduced ring-opening reaction. The surfaces have been characterized in terms of their stationary points and points of minimum energy conical intersections along the ring-opening pathway. The optimization of the geometrical parameters has been performed with the equation of motion coupled cluster singles and doubles method. The ab initio PESs have been modeled by energy grids and Taylor series. The resulting 11-dimensional PESs reproduce the ab initio results to a good accuracy and can be used in dynamical calculations.

The core level electron excitation and ionization spectra of thymine and adenine have been investigated by photoabsorption and photoemission spectroscopy, and the results interpreted by means of ab initio calculations using the second-order algebraic-diagrammatic construction (ADC(2)) method for the polarization propagator and the fourth-order ADC method (ADC(4)) for the one-particle Green's function. The photoabsorption spectra are dominated by transitions from core levels to unoccupied pi states, but also show clear structures due to Rydberg transitions. The calculated spectra are in good agreement with the experimental results, and many of the observed structures are assigned. (C) 2007 Elsevier B.V. All rights reserved.

An earlier proposed propagator method for the treatment of molecular ionization is tested in first applications. The method referred to as the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator represents a computationally promising alternative to the existing Dyson ADC(3) method. The advantage of the nD-ADC(3) scheme is that the (N+/-1)-electronic parts of the one-particle Green's function are decoupled from each other and the corresponding equations can be solved separately. For a test of the method the nD-ADC(3) results for the vertical ionization transitions in C(2)H(4), CO, CS, F(2), H(2)CO, H(2)O, HF, N(2), and Ne are compared with available experimental and theoretical data including results of full configuration interaction (FCI) and coupled cluster computations. The mean error of the nD-ADC(3) ionization energies relative to the experimental and FCI results is about 0.2 eV. The nD-ADC(3) method, scaling as n(5) with the number of orbitals, requires the solution of a relatively simple Hermitian eigenvalue problem. The method renders access to ground-state properties such as dipole moments. Moreover, also one-electron properties of (N+/-1) electron states can now be studied as a consequence of a specific intermediate-state representation (ISR) formulation of the nD-ADC approach. Corresponding second-order ISR equations are presented.

Propagator methods provide a direct approach to energies and transition moments for (generalized) electronic excitations from the ground state, but they do not usually allow one to determine excited state wave functions and properties. Using a specific intermediate state representation (ISR) concept, we here show how this restriction can be overcome in the case of the algebraic-diagrammatic construction (ADC) propagator approach. In the ISR reformulation of the theory the basic ADC secular matrix is written as a representation of the Hamiltonian (or the shifted Hamiltonian) in terms of explicitly constructable states, referred to as intermediate (or ADC) states. Similar intermediate state representations can be derived for operators other than the Hamiltonian. Together with the ADC eigenvectors, the intermediate states give rise to an explicit formulation of the excited wave functions and allow one to calculate physical properties of excited states as well as transition moments for transitions between different excited states. As for the ground-state excitation energies and transition moments, the ADC excited state properties are size consistent so that the theory is suitable for applications to large systems. The established hierarchy of higher-order [ADC(n)] approximations, corresponding to systematic truncations of the IS configuration space and the perturbation-theoretical expansions of the ISR matrix elements, can readily be extended to the excited state properties. Explicit ISR matrix elements for arbitrary one-particle operators have been derived and coded at the second-order [ADC(2)] level of theory. As a first computational test of the method we have carried out ADC(2) calculations for singlet and triplet excited state dipole moments in H(2)O and HF, where comparison to full CI results can be made. The potential of the ADC(2) method is further demonstrated in an exploratory study of the excitation energies and dipole moments of the low-lying excited states of paranitroaniline. We find that four triplet states, T1-T4, and two singlet states, S1 and S2, lie (vertically) below the prominent charge transfer (CT) excitation, S3. The dipole moment of the S3 state (17.0D) is distinctly larger than that of the corresponding T3 triplet state (11.7D).

The one-particle Green’s function theory in its various implementations is a well-established many-body approach for the calculation of electron ionization and attachment energies in atoms and molecules. In order to describe not only scalar-relativistic effects but also spin–orbit splitting on an equal footing an embedding of this theory in the four-component framework was carried out and fully relativistic ionization energies of the noble gas atoms Ne through Xe were calculated using the second-order algebraic diagrammatic construction [ADC(2)] approximation scheme. Comparison with nonrelativistic ADC(2) results and experimental data was made.

A new polarization propagator approach to indirect nuclear spin-spin coupling constantans is formulated within the framework of the algebraic-diagrammatic construction (ADC) approximation and implemented at the level of the strict second-order approximation scheme, ADC(2). The ADC approach possesses transparent computational procedure operating with Hermitian matrix quantities defined with respect to physical excitations. It is size-consistent and easily extendable to higher orders via the hierarchy of available ADC approximation schemes. The ADC(2) method is tested in the first applications to HF, N(2), CO, H(2)O, HCN, NH(3), CH(4), C(2)H(2), PH(3), SiH(4), CH(3)F, and C(2)H(4). The calculated indirect nuclear spin-spin coupling constants are in good agreement with the experimental data and results of the second-order polarization propagator approximation method. The computational effort of the ADC(2) scheme scales as n(5) with respect to the number of molecular orbitals n, which makes this method promising for applications to larger molecules.

An earlier proposed approach to molecular response functions based on the intermediate state representation (ISR) of polarization propagator and algebraic-diagrammatic construction (ADC) approximations is for the first time employed for calculations of nonlinear response properties. The two-photon absorption (TPA) spectra are considered. The hierarchy of the first- and second-order ADC∕ISR computational schemes, ADC(1), ADC(2), ADC(2)-x, and ADC(3/2), is tested in applications to H(2)O, HF, and C(2)H(4) (ethylene). The calculated TPA spectra are compared with the results of coupled cluster (CC) models and time-dependent density-functional theory (TDDFT) calculations, using the results of the CC3 model as benchmarks. As a more realistic example, the TPA spectrum of C(8)H(10) (octatetraene) is calculated using the ADC(2)-x and ADC(2) methods. The results are compared with the results of TDDFT method and earlier calculations, as well as to the available experimental data. A prominent feature of octatetraene and other polyene molecules is the existence of low-lying excited states with increased double excitation character. We demonstrate that the two-photon absorption involving such states can be adequately studied using the ADC(2)-x scheme, explicitly accounting for interaction of doubly excited configurations. Observed peaks in the experimental TPA spectrum of octatetraene are assigned based on our calculations.

The valence shell photoelectron spectra of 2-chloropyridine and 3-chloropyridine have been studied both experimentally and theoretically. Synchrotron radiation has been employed to record angle resolved photoelectron spectra in the photon energy range 20–100 eV, and these have enabled anisotropy parameters and branching ratios to be derived. The experimental results have been compared with theoretical predictions obtained using the continuum multiple scattering Xα approach. This comparison shows that the anisotropy parameter associated with the nominally chlorine lone-pair orbital lying in the molecular plane is strongly affected by the atomic Cooper minimum. In contrast, the photoionization dynamics of the second lone-pair orbital, orientated perpendicular to the molecular plane, seem relatively unaffected by this atomic phenomenon. The outer valence ionization has been studied theoretically using the third-order algebraic-diagrammatic construction (ADC(3)) approximation scheme for the one-particle Green’s function, the outer valence Green’s function method, and the equation-of-motion (EOM) coupled cluster (CC) theory at the level of the EOM-IP-CCSD and EOM-EE-CC3 models. The convergence of the results to the complete basis set limit has been investigated. The ADC(3) method has been employed to compute the complete valence shell ionization spectra of 2-chloropyridine and 3-chloropyridine. The relaxation mechanism for ionization of the nitrogen σ-type lone-pair orbital (σN LP) has been found to be different to that for the corresponding chlorine lone-pair (σCl LP). For the σN LP orbital, π-π* excitations play the main role in the screening of the lone-pair hole. In contrast, excitations localized at the chlorine site involving the chlorine πCl LP lone-pair and the Cl 4p Rydberg orbital are the most important for the σCl LP orbital. The calculated photoelectron spectra have allowed assignments to be proposed for most of the structure observed in the experimental spectra. The theoretical work also highlights the formation of satellite states, due to the breakdown of the single particle model of ionization, in the inner valence region.