Search

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

We report a computational study and analysis of the optical absorption and photodecay processes in two subnanometer metal complexes deposited on an oxide support, the regular MgO(100) surface: (i) Ag3(HCO3)(C2H4)2(O) and (ii) Ag3(CO2F)(C2H4)2(O). These aggregates are chosen as derivatives of a Ag3(CO3)(C2H4)2(O) ligand/metal-cluster/support complex, previously singled out as a key intermediate in the path of ethylene partial oxidation to ethylene epoxide catalyzed by Ag3/MgO(100), and serve as model systems to investigate photochemical phenomena in ligand/metal-cluster/support complexes by subnanometer metal catalysts, an appealing field for future research. After generating optimized initial configurations and building cluster models that take properly into account the effect of the charge-separated oxide support, we use time-dependent density-functional theory (TDDFT) to determine first the photoabsorption spectra of the two aggregates and then to follow the evolution of their excited states in the optical region. We show that complexes containing such bicarbonate and fluorocarbonate groups are sensitive to optical adsorption, often leading to ligand detachment and/or cluster disaggregation, thus pointing to an "optical frailty"of these subnanometer cluster species, possibly rationalizing previous experimental observations. Additionally, we correlate the nature of the given excitations and of the corresponding photoinduced reaction products via an analysis of overlap population-density of states (OP-DOS), geometric parameters, and spatial distribution of the molecular orbitals involved in the excitation, thus providing the set of methodological tools needed to explore this novel field.

The near-edge X-ray absorption fine structure (NEXAFS and X-ray photoelectron (XP) spectra of gas-phase 2,8-bis-(diphenylphosphoryl)dibenzo[b,d]thiophene (PPT) and triphenylphosphine oxide (TPPO) have been measured at the S and P L-II,L-III-edge regions. The time-dependent density functional theory (TDDFT) based on the relativistic two-component zeroth-order regular approximation approach has been used to provide an assignment of the experimental spectra, giving the contribution of the spin-orbit splitting and of the molecular-field splitting to the sulfur and phosphor binding energies. Computed XP and NEXAFS spectra agree well with the experimental measurements. In going from dibenzothiophene and TPPO to PPT, the nature of the most intense S 2p and P 2p NEXAFS features are preserved; this trend suggests that the electronic and geometric behaviors of the S and P atoms in the two building block moieties are conserved in the more complex system of PPT. This work enables us to shed some light onto the structure of the P-O bond, a still highly debated topic in the chemical literature. Since the S 2p and P 2p NEXAFS intensities provide specific information on the higher-lying localized sigma*(C-S) and sigma*(P-O) virtual MOs, we have concluded that P 3d AOs are not involved in the formation of the P-O bond. Moreover, the results support the mechanism of negative hyperconjugation, by showing that transitions toward sigma*(P-O) states occur at lower energies with respect to those toward it pi*(P-O) states.

The photoionization dynamics of OsO4 and RuO4, chosen as model systems of small-size mononuclear heavy-metal complexes, has been theoretically studied by the time-dependent density functional theory (TDDFT). Accurate experimental measurements of photoionization dynamics as a benchmarking test for the theory are reported for the photoelectron asymmetry parameters of outer valence ionizations of OsO4, measured in the 17-90 eV photon energy range. The theoretical results are in good agreement with the available experimental data. The observed dynamical behavior of partial cross sections and asymmetry parameters has been related to both the coupling to the continuum of discrete excited states, giving strong modulations in the photon energy dependency, and the atomic composition of the initial ionized states, which determines the rate of decay of ionization probability for increasing excitation energies. Overall, an extensive analysis of the photoionization dynamics for valence and core orbitals is presented, showing good agreement with all the available experimental data. This provides confidence for the validity of the TDDFT approach in describing photoionization of heavy transition element compounds, with the perspective of being used for larger systems. Further experimental work is suggested for RuO4 to gather evidence of the sensitivity of the theoretical method to the nature of the metal atom.

The valence and core-level photoelectron spectra of gaseous indole, 2,3-dihydro-7-azaindole, and 3-formylindole have been investigated using VUV and soft X-ray radiation supported by both an ab initio electron propagator and density functional theory calculations. Three methods were used to calculate the outer valence band photoemission spectra: outer valence Green function, partial third order, and renormalized partial third order. While all gave an acceptable description of the valence spectra, the last method yielded very accurate agreement, especially for indole and 3-formylindole. The carbon, nitrogen, and oxygen is core-level spectra of these heterocycles were measured and assigned. The double ionization appearance potential for indole has been determined to be 21.8 +/- 0.2 eV by C is and N is Auger photoelectron spectroscopy. Theoretical analysis identifies the doubly ionized states as a band consisting of two overlapping singlet states and one triplet state with dominant configurations corresponding to holes in the two uppermost molecular orbitals. One of the singlet states and the triplet state can be described as consisting largely of a single configuration, but other doubly ionized states are heavily mixed by configuration interactions. This work provides full assignment of the relative binding energies of the core level features and an analysis of the electronic structure of substituted indoles in comparison with the parent indole.

We performed a theoretical investigation on the influence of electronic correlation effects on the B1s NEXAFS spectrum of boronic acid derivatives, namely, boric acid [B(OH)(3)], phenyl boronic acid (PBA), and 1,4-phenyl diboronic acid (PDBA), employing different computational schemes of increasing complexity, ranging from the purely one-electron scheme based on the transition potential method of density functional theory (DFT-TP), time-dependent DFT (TDDFT), and multiconfigurational self-consistent field (MCSCF). We also report experimental measurements of the B1s NEXAFS spectra of the aforementioned molecules together with the high-resolution C1s NEXAFS spectrum of PBA. We demonstrate that due to the shallow B1s core energy levels compared to C, O, and N, the inclusion of static correlation effects, which can be incorporated by using multireference approaches to excited states, assumes a decisive role in reconciling experiment and theory on B1s core-electron excitation energies and oscillator strengths to valence states. This claim is corroborated by the good agreement that we find between the DFT-TP calculated C1s NEXAFS spectrum and that experimentally measured for PBA and by the failure of both DFT-TP and TDDFT approaches with a selection of xc functionals kernels to properly describe the B1s NEXAFS spectrum of PBA and PDBA, at variance with the good agreement with the experiment that is found by employing the MCSCF wave function approach.

Understanding the adsorption properties of organic molecules on graphene-based substrates is important for such applications as air and water filters. Pristine graphene is often the model substrate used in the theoretical investigations of this problem. While useful, pristine single-layer graphene is however an idealized model. In this work, we assess the effect of the presence of point defects (single vacancy, divacancy, and the Stone-Wales defect) in single-layer and bilayer graphene on the energetics of adsorption of benzene and toluene. Our calculations benchmark three different dispersion-corrected DFT schemes, namely PBE-D2, vdW-DF1, and vdW-DF2-C09. Whereas the presence of the single vacancy and the double vacancy does not appear to alter the adsorption energies of the aromatic molecules by an appreciable amount, the Stone-Wales defect and the addition of a second graphene layer stabilizes their interaction with the substrate by several tens of meV.

The near-edge x-ray-absorption fine-structure (NEXAFS) and X-ray photoelectron spectroscopy (XPS) spectra of benzo[b]thiophene (BBT) and dibenzothiophene (DBT) in the gas phase have been measured at the carbon Kedge and sulfur L-II,L-III-edge regions. The assignment of the spectral features has been provided by theoretical calculations based on density functional theory (DFT) and its time-dependent generalization (TDDFT) in the linear response regime. Observed trends in computed C 1s and S 2p ionization potentials (IPs) have been rationalized in terms of both the inductive effects due to the presence of S and the increased pi-electrons delocalization arising from the benzoannulation process. The analysis of the NEXAFS carbon K-edge and sulfur L-II,L-III-edge regions regions provided information on both low-lying delocalized virtual pi orbitals, and higher-lying localized sigma*(C-S) states. The evolution of the NEXAFS carbon K-edge spectral features along the series thiophene (T) and derivatives, BBT and DBT, is informative of a stabilizing effect due to increased aromaticity. This effect is however more pronounced in going from T to BBT compared to the introduction of a second annulated phenyl ring in DBT. The nature of the most intense sulfur L-II,L-III-edge NEXAFS spectral features is instead conserved along the series reflecting thus the localized nature of the virtual states involved in the S 2p core-excitation process.

The Near-Edge X-ray-Absorption Fine-Structure (NEXAFS) and X-ray Photoemission Spectroscopy (XPS) of short-chain oligothiophenes (thiophene, 2,2'-bithiophene, and 2,2':5',2 ''-terthiophene) in the gas phase have been measured in the sulfur L-2,L-3-edge region. The assignment of the spectral features is based on the relativistic two-component zeroth-order regular approximation time dependent density functional theory approach. The calculations allow us to estimate both the contribution of the spin-orbit splitting and of the molecular-field splitting to the sulfur binding energies and give results in good agreement with the experimental measurements. The deconvolution of the calculated S2p NEXAFS spectra into the two manifolds of excited states converging to the L-III and L-II edges facilitates the attribution of the spectral structures. The main S2p NEXAFS features are preserved along the series both as concerns the energy positions and the nature of the transitions. This behaviour suggests that the electronic and geometrical environment of the sulfur atom in the three oligomers is relatively unaffected by the increasing chain length. This trend is also observed in the XPS spectra. The relatively simple structure of S2p NEXAFS spectra along the series reflects the localized nature of the virtual states involved in the core excitation process. Published by AIP Publishing.

A general implementation of the lowest nonvanishing order perturbation theory for the calculation of molecular multiphoton ionization cross sections is proposed in the framework of density functional theory. Bound and scattering wave functions are expanded in a multicentric basis set and advantage is taken of the full molecular point group symmetry, thus enabling the application of the formalism to medium-size molecules. Multiphoton ionization cross sections and angular asymmetry parameters have been calculated for the two- and four-photon ionization of the H-2(+) molecule, for linear and circular light polarizations. Both fixed and random orientations of the target molecule have been considered. To demonstrate the efficiency of the proposed methodology, the two-photon cross section and angular asymmetry parameters for the HOMO and HOMO-1 orbital ionization of benzene are also presented. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4754820]