1
|
Ahn DH, Nakajima T, Hirao K, Song JW. Long-range Corrected Density Functional Theory Including a Two-Gaussian Hartree-Fock Operator for High Accuracy Core-excitation Energy Calculations of Both the Second- and Third-Row Atoms (LC2gau-core-BOP). J Chem Theory Comput 2024. [PMID: 39106473 DOI: 10.1021/acs.jctc.4c00651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/09/2024]
Abstract
In the previous work, LCgau-core-BOP, which includes the short-range interelectronic Gaussian attenuating Hartree-Fock (HF) exchange to the long-range HF exchange, showed high accuracy core-excitation energies from 1s orbitals of the second-row atoms (1s → π*, 1s → σ*, 1s → n*, and 1s → Rydberg), but underestimates the core-excitation energies from 1s orbitals of the third-row atoms. To improve this, we added one more Gaussian attenuating HF exchange to LCgau-core-BOP. We named it LC2gau-core-BOP, which achieves a mean absolute error (MAE) of 0.6 and 0.3 eV for core excitation energies of the second- and third-row atoms of the tested small molecules, respectively. We found that the inclusion of the short-range interelectronic HF exchange at a distance ranging from 0.2 to 0.6 a.u. contributes to the increase of performances on 1s orbital energy calculations of the second-row atoms, while the inclusion of more short-range interelectronic HF exchange at a distance ranging from 0 to 0.2 a.u. does to the increase of performance on 1s orbital energy calculations of the third-row atoms. It is notable that all of these improvements were accomplished using flexible Gaussian attenuating HF exchange inclusion. LC2gau-core-BOP shows deviations of less than 0.8 eV from experimental values for all of the core-excitation energies of the tested medium-size molecules consisting of thymine, oxazole, glycine, and dibenzothiophene sulfone. Moreover, by optimizing one parameter of the OP correlation functional, LC2gau-core-BOP provides atomization energies over the G3 test set with an accuracy comparable to that of B3LYP.
Collapse
Affiliation(s)
- Dae-Hwan Ahn
- Department of Chemistry Education, Daegu University, Gyeongsan-si 113-8656, Korea
| | | | - Kimihiko Hirao
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8501, Japan
| | - Jong-Won Song
- Department of Chemistry Education, Daegu University, Gyeongsan-si 113-8656, Korea
| |
Collapse
|
2
|
Park W, Lashkaripour A, Komarov K, Lee S, Huix-Rotllant M, Choi CH. Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT. J Chem Theory Comput 2024; 20:5679-5694. [PMID: 38902891 DOI: 10.1021/acs.jctc.4c00640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2024]
Abstract
Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Alireza Lashkaripour
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
- Department of Chemistry, University of Zürich, Zürich 8057, Switzerland
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
3
|
Fransson T, Pettersson LGM. TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional. J Chem Phys 2024; 160:234105. [PMID: 38884399 DOI: 10.1063/5.0209719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 06/03/2024] [Indexed: 06/18/2024] Open
Abstract
We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| |
Collapse
|
4
|
Gaba NP, de Moura CEV, Majumder R, Sokolov AY. Simulating transient X-ray photoelectron spectra of Fe(CO) 5 and its photodissociation products with multireference algebraic diagrammatic construction theory. Phys Chem Chem Phys 2024; 26:15927-15938. [PMID: 38805029 DOI: 10.1039/d4cp00801d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Accurate simulations of transient X-ray photoelectron spectra (XPS) provide unique opportunities to bridge the gap between theory and experiment in understanding the photoactivated dynamics in molecules and materials. However, simulating X-ray photoelectron spectra along a photochemical reaction pathway is challenging as it requires accurate description of electronic structure incorporating core-hole screening, orbital relaxation, electron correlation, and spin-orbit coupling in excited states or at nonequilibrium ground-state geometries. In this work, we employ the recently developed multireference algebraic diagrammatic construction theory (MR-ADC) to investigate the core-ionized states and X-ray photoelectron spectra of Fe(CO)5 and its photodissociation products (Fe(CO)4, Fe(CO)3) following excitation with 266 nm light. The simulated transient Fe 3p and CO 3σ XPS spectra incorporating spin-orbit coupling and high-order electron correlation effects are shown to be in a good agreement with the experimental measurements by Leitner et al. [J. Chem. Phys., 2018, 149, 044307]. Our calculations suggest that core-hole screening, spin-orbit coupling, and ligand-field splitting effects are similarly important in reproducing the experimentally observed chemical shifts in transient Fe 3p XPS spectra of iron carbonyl complexes. Our results also demonstrate that the MR-ADC methods can be very useful in interpreting the transient XPS spectra of transition metal compounds.
Collapse
Affiliation(s)
- Nicholas P Gaba
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA.
| | - Carlos E V de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA.
| | - Rajat Majumder
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA.
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA.
| |
Collapse
|
5
|
Safari AA, Anderson RJ, Manni GL. Toward a Stochastic Complete Active Space Second-Order Perturbation Theory. J Phys Chem A 2024; 128:281-291. [PMID: 38154124 PMCID: PMC10788896 DOI: 10.1021/acs.jpca.3c05109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 11/14/2023] [Accepted: 11/14/2023] [Indexed: 12/30/2023]
Abstract
In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.
Collapse
Affiliation(s)
- Arta A. Safari
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| | | | - Giovanni Li Manni
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Segatta F, Aranda D, Aleotti F, Montorsi F, Mukamel S, Garavelli M, Santoro F, Nenov A. Time-Resolved X-ray Absorption Spectroscopy: An MCTDH Quantum Dynamics Protocol. J Chem Theory Comput 2024; 20:307-322. [PMID: 38101807 PMCID: PMC10782456 DOI: 10.1021/acs.jctc.3c00953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 11/04/2023] [Accepted: 11/09/2023] [Indexed: 12/17/2023]
Abstract
Expressions for linear and nonlinear spectroscopy simulation in the X-ray window in which the time evolution of a photoexcited molecular system is treated via quantum dynamics are derived. By leveraging on the peculiar properties of core-excited/ionized states, first- and third-order response functions are recast in the limit of time-scale separation between the extremely short core-state lifetime and the (comparably longer) electronic-state transfer and nuclear vibrational motion. This work is a natural extension of Segatta et al. (J. Chem. Theory Comput. 2023, 19, 2075-2091), in which some of the present authors coupled MCTDH quantum dynamics to spectroscopy simulation at different levels of sophistication. Full quantum dynamics and approximate expressions are compared by simulating X-ray transient absorption spectroscopy at the carbon K-edge in the pyrene molecule.
Collapse
Affiliation(s)
- Francesco Segatta
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Daniel Aranda
- ICMol, Universidad de Valencia, c/Catedrático José
Beltrán,
2, 46980 Paterna, Spain
- Istituto
di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via
Moruzzi 1, I-56124 Pisa, Italy
| | - Flavia Aleotti
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Francesco Montorsi
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Shaul Mukamel
- Department
of Chemistry and Department of Physics and Astronomy, University of California, Irvine, 92697 California, United States
| | - Marco Garavelli
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Fabrizio Santoro
- Istituto
di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via
Moruzzi 1, I-56124 Pisa, Italy
| | - Artur Nenov
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| |
Collapse
|
7
|
Scott M, Delcey MG. Complex Linear Response Functions for a Multiconfigurational Self-Consistent Field Wave Function in a High Performance Computing Environment. J Chem Theory Comput 2023; 19:5924-5937. [PMID: 37596971 PMCID: PMC10500980 DOI: 10.1021/acs.jctc.3c00317] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Indexed: 08/21/2023]
Abstract
We present novel developments for the highly efficient evaluation of complex linear response functions of a multiconfigurational self-consistent field (MCSCF) wave function as implemented in MultiPsi. Specifically, expressions for the direct evaluation of linear response properties at given frequencies using the complex polarization propagator (CPP) approach have been implemented, within both the Tamm-Dancoff approximation (TDA) and the random phase approximation (RPA). Purely real algebra with symmetric and antisymmetric trial vectors in a shared subspace is used wherein the linear response equations are solved. Two bottlenecks of large scale MC-CPP calculations, namely, the memory footprint and computational time, are addressed. The former is addressed by limiting the size of the subspace of trial vectors by using singular value decomposition (SVD) on either orbital or CI subspaces. The latter is addressed using an efficient parallel implementation as well as the strategy of dynamically adding linear response equations at near-convergence to neighboring roots. Furthermore, a novel methodology for decomposing MC-CPP spectra in terms of intuitive orbital excitations in an approximate fashion is presented. The performance of the code is illustrated with several numerical examples, including the X-ray spectrum of a molecule with nearly one hundred atoms. Additionally, for X-ray spectroscopy, the effect of including or excluding the core orbital in the active space on small covalent metal complexes is discussed.
Collapse
Affiliation(s)
- Mikael Scott
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Mickael G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
- Division
of Theoretical Chemistry, Department of Chemistry, Lund University, SE-221 00 Lund, Sweden
| |
Collapse
|
8
|
Kaczun T, Dempwolff AL, Huang X, Gelin MF, Domcke W, Dreuw A. Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation. J Phys Chem Lett 2023:5648-5656. [PMID: 37310800 DOI: 10.1021/acs.jpclett.3c01018] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.
Collapse
Affiliation(s)
- Tobias Kaczun
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| | - Xiang Huang
- Department of Chemistry, Technical University of Munich, Garching D-85747, Germany
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching D-85747, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| |
Collapse
|
9
|
Carlini L, Montorsi F, Wu Y, Bolognesi P, Borrego-Varillas R, Casavola AR, Castrovilli MC, Chiarinelli J, Mocci D, Vismarra F, Lucchini M, Nisoli M, Mukamel S, Garavelli M, Richter R, Nenov A, Avaldi L. Electron and ion spectroscopy of azobenzene in the valence and core shells. J Chem Phys 2023; 158:054201. [PMID: 36754795 DOI: 10.1063/5.0133824] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Azobenzene is a prototype and a building block of a class of molecules of extreme technological interest as molecular photo-switches. We present a joint experimental and theoretical study of its response to irradiation with light across the UV to x-ray spectrum. The study of valence and inner shell photo-ionization and excitation processes combined with measurement of valence photoelectron-photoion coincidence and mass spectra across the core thresholds provides a detailed insight into the site- and state-selected photo-induced processes. Photo-ionization and excitation measurements are interpreted via the multi-configurational restricted active space self-consistent field method corrected by second order perturbation theory. Using static modeling, we demonstrate that the carbon and nitrogen K edges of azobenzene are suitable candidates for exploring its photoinduced dynamics thanks to the transient signals appearing in background-free regions of the NEXAFS and XPS.
Collapse
Affiliation(s)
- L Carlini
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - F Montorsi
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - Y Wu
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - P Bolognesi
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - R Borrego-Varillas
- CNR-Istituto di Fotonica e Nanotecnologie, CNR-IFN, Piazza Leonardo da Vinci 32, Milano, Italy
| | - A R Casavola
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - M C Castrovilli
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - J Chiarinelli
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| | - D Mocci
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - F Vismarra
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - M Lucchini
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - M Nisoli
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - S Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| | - M Garavelli
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - R Richter
- Elettra Sincrotrone Trieste, Area Science Park, Basovizza, Italy
| | - A Nenov
- Dipartimento di Chimica Industriale, Università Degli Studi di Bologna, Bologna, Italy
| | - L Avaldi
- CNR-Istituto di Struttura Della Materia, CNR-ISM, Area Della Ricerca di Roma 1, Monterotondo, Italy
| |
Collapse
|
10
|
Park W, Alías-Rodríguez M, Cho D, Lee S, Huix-Rotllant M, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory for Accurate X-ray Absorption Spectroscopy. J Chem Theory Comput 2022; 18:6240-6250. [PMID: 36166346 DOI: 10.1021/acs.jctc.2c00746] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
It is demonstrated that the challenging core-hole particle (CHP) orbital relaxation for core electron spectra can be readily achieved by the mixed-reference spin-flip (MRSF)-time-dependent density functional theory (TDDFT). With the additional scalar relativistic effects on K-edge excitation energies of 24 second- and 17 third-row molecules, the particular ΔCHP-MRSF(R) exhibited near perfect predictions with RMSE ∼0.5 eV, featuring a median value of 0.3 and an interquartile range of 0.4. Overall, the CHP effect is 2-4 times stronger than relativistic ones, contributing more than 20 eV in the cases of sulfur and chlorine third-row atoms. Such high precision allows to explain the splitting and spectral shapes of O, N, and C atom K-edges in the ground state of thymine with atom as well as orbital specific accuracy. The same protocol with a double hole particle relaxation also produced remarkably accurate K-edge spectra of core to valence hole excitation energies from the first (nO8π*) and second (ππ*) excited states of thymine, confirming the assignment of 1s → n excitation for the experimentally observed 526.4 eV peak. Regarding both accuracy and practicality, therefore, MRSF-TDDFT provides a promising protocol for core electron spectra of both ground and excited electronic states alike.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Marc Alías-Rodríguez
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Daeheum Cho
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
11
|
Cheng X, Wei M, Tian G, Luo Y, Hua W. Vibrationally-Resolved X-ray Photoelectron Spectra of Six Polycyclic Aromatic Hydrocarbons from First-Principles Simulations. J Phys Chem A 2022; 126:5582-5593. [PMID: 35959595 DOI: 10.1021/acs.jpca.2c04426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Vibrationally resolved C 1s X-ray photoelectron spectra (XPS) of a series of six polycyclic aromatic hydrocarbons (PAHs; phenanthrene, coronene, naphthalene, anthracene, tetracene, and pentacene) were computed by combining the full core hole density functional theory and the Franck-Condon simulations with the inclusion of the Duschinsky rotation effect. Simulated spectra of phenanthrene, coronene, and naphthalene agree well with experiments both in core binding energies (BEs) and profiles, which validate the accuracy of our predictions for the rest molecules with no high-resolution experiments. We found that three types of carbons i (inner C), p (peripheral C bonded to three C atoms), and h (peripheral C bonded to an H atom) show decreasing BEs. In linear PAHs (the latter four), h-type carbons further split into h1 or h2 (on inner or edge benzene ring) subtypes with chemical shifts of ca. 0.2-0.4 eV. All major Franck-Condon-active modes are characterized to be in-plane vibrations: low-frequency (<800 cm-1) C-C ring deformation modes play an essential role in determining the peak asymmetries; and for each h-type carbon a high-frequency (ca. 3600 cm-1) C*-H stretching mode is responsible for the high-energy tail. We found that core ionization leads to reduction of all C*-C and C*-H bond lengths and ring deformation with a definite direction. Based on theoretical spectra of four linear PAHs, we found asymptotic relations and anticipated possible spectral features for even larger linear PAHs. Our calculations provide accurate reference spectra for XPS characterizations of PAHs, which are useful in understanding the vibronic coupling effects in this family.
Collapse
Affiliation(s)
- Xiao Cheng
- Hefei National Laboratory for Physical Science at the Microscale, University of Science and Technology of China, 230026 Hefei, China.,Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, Royal Institute of Technology, SE-106 91, Stockholm, Sweden
| | - Minrui Wei
- MIIT Key Laboratory of Semiconductor Microstructure and Quantum Sensing, Department of Applied Physics, School of Science, Nanjing University of Science and Technology, 210094 Nanjing, China
| | - Guangjun Tian
- Key Laboratory for Microstructural Material Physics of Hebei Province, School of Science, Yanshan University, Qinhuangdao 066004, China
| | - Yi Luo
- Hefei National Laboratory for Physical Science at the Microscale, University of Science and Technology of China, 230026 Hefei, China
| | - Weijie Hua
- MIIT Key Laboratory of Semiconductor Microstructure and Quantum Sensing, Department of Applied Physics, School of Science, Nanjing University of Science and Technology, 210094 Nanjing, China
| |
Collapse
|