Simulation of Xray Absorption Near Edge Spectroscopy XANES
Simulation of X-ray Absorption Near Edge Spectroscopy (XANES) of Molecules Luke Campbell Shaul Mukamel Daniel Healion Rajan Pandey
Motivation • X-ray Absorption Near Edge Spectroscopy (XANES) is an attractive tool for measuring local changes in electronic structure due to geometry and charge distribution of transient species. • Recent advances in ultrashort (femtosecond to attosecond) x-ray pulses enable real time probing of optically induced electron motions and chemical processes. • Time resolved XANES measures changes in geometry and charge distribution during and after the excitation. • Theory can provide a guide for the design and interpretation of these measurements.
Basic Physics of X-ray Absorption • X-ray absorption probes the unoccupied dipole allowed one electron density of states of a molecule in the vicinity of the absorbing atom. Ø µ(ω): absorption coefficient, intensity for depth x. Ø σ(ω): absorption cross section. Ø : initial state with energy Ei. Ø : final state with energy Ef ; only transitions to unoccupied states are allowed. Ø : dipole operator (core size much smaller than x-ray wavelength). • Localized core → only local DOS contributes.
Methodology Sum Over States Method (SOS): Ø Many-electron ground states (with and without core holes) are calculated using standard quantum chemistry codes. within density functional theory or Hartree-Fock approximation, (Z+1 approximation, where Z is the nuclear charge). Ø Electronically excited states are calculated using time dependent density functional theory (TDDFT) or time dependent Hartree-Fock (TDHF) theory. Ø Computationally expensive, requires explicit calculation of excited states.
Transition Potential Method: Ø Uses a reference system with partially filled orbitals (incorporated in the Sto. Be Demon code). Ø Represents systems with different numbers of core holes by different occupation numbers of a single set of reference orbitals. Ø Computationally less expensive than SOS. Ø Works well for core level spectroscopies of small molecules.
Simulation of x-ray absorption near edge spectra (XANES) of molecules • Start with the Deep Core Hamiltonian • Neglect valence-core exchange Valence Core • Electron-electron interaction • One-electron valence terms • Core hole potential → use Z+1 approximation, core hole approximated as point charge → equivalent to nuclear charge increased by 1. Interaction
Fermi’s Golden Rule gives the absorption cross section: Dipole operator in ν direction Dipole matrix element → Initial wavefunction with energy Ei. → Final wavefunction with energy Ef. → Electron annihilation (creation) operator for orbital l.
Core-valence separation • Deep core Hamiltonian → separate eigenvalue problem for valence and core electrons → can represent as product space In the Z+1 approximation: → Initial valence wavefunction. → Fully occupied core wavefunction. → Final valence wavefunction with core hole potential present. → Core wavefunction with orbital g unoccupied.
• Effective valence Hamiltonians • Core filled (initial state) valence Hamiltonian: • Valence Hamiltonian with core hole in orbital g: • The absorption spectrum:
First principles computation of ground and excited state XANES Of chemical species Ø Use quantum chemistry code (Gaussian 03) to find electronic structure of ground and excited states. Ø Find energies and intensities of transitions from a given initial ground or excited state to possible final excited states. Ø Basis set: Selection based on kind of chemical species in a molecule Ø Level of theory: Becke 3 -parameter density functional with Lee-Yang-Parr correlation, Hartree-Fock approximation. Ø Code: GAUSSIAN-03 Ø Geometry: from x-ray crystallography data (complex molecules). Ø Ground state: • singlet spin • 5 -15 singlet and/or triplet excited states with TDDFT or TDHF Ø Core excited state: • Z+1 approximation • doublet spin • 50 or more excited states with TDDFT/TDHF
[Ru(bpy)3]2+ Experimental XANES L 3 -Edge • 1 e. V valence shift of main peak (B → B') after photoexcitation to 3 MLCT state. • Appearance of new peak A' after photoexcitation.
[Ru(bpy)3]2+ SOS Simulated XANES L 3 -Edge B 3 LYP/3 -21 G • Ground state XANES (solid line) shows peak B. • MLCT XANES (dottes) shows peak B' blue shifted by 1 e. V and appearance of peak A'. Luke Campbell and Shaul Mukamel, J. Chem. Phys. 121, 12323 (2004).
Excited State Effects on X-ray Absorption Charge transfer to or from the absorbing atom can alter the energies and intensities of transitions to the bound states. Examples: • Removing an electron makes the atom more positively charged, so more energy is needed to excite the core electron to orbitals farther from atom. Absorption peaks shift position • When electrons are taken out of previously filled orbitals, new core → valence transitions are possible. • When electrons are put into previously empty orbitals, peaks can disappear.
Single and Double Excitations Neglecting changes in orbitals due to core excitation: • From any initial optically excited state, the final XANES state (a) can be reproduced with two excitations from the lowest core excited state (b). l (b) (a) • From some initial states, such as the ground state or HOMO to LUMO excitations, the final XANES state can be represented by one excitation from the lowest core excited state (b). Transition (1) gives ground state XANES (a), transition (2) gives HOMO to LUMO excitation XANES (c). (a) l l (1) (2) (b) (c)
XANES spectra of water (O K-edge) 1. 90 e. V HF/6 -311++G** 2 b 2 4 a 1 Energy XANES Ionization potential H H X-ray photon O Absorption H H O Water monomer
Peak splitting between the lowest transitions corresponding to 1 a 1 → 4 a 1 and 1 a 1 → 2 b 2 1. 90 e. V 1. 92 e. V 2. 04 e. V 1. 83 e. V Sum Over States SOS (solid line) gives a good agreement with the experiment. Plots and numbers reproduced (except solid curve - SOS) from Ref: M. Cavalleri et al. J. Chem. Phys. Vol. 121, 10074 (2004)
Methyl Alcohol O K-Edge SOS Transition Potential
XANES of Benzonitrile (N K-edge) Method/Basis TDDFT (B 3 LYP)/D 95** Gives good agreement for the intensity ratio. However, peak splitting is not exact. TDHF/D 95** Gives good agreement in the peak splitting. However, the intensity ratio is different than experiment. Ref: S. Carniato et al. Phys. Rev A 58, 022511 (2005).
X-Ray Fluorescence Hamiltonians in the Z+1 approximation:
Fluorescence Spectrum of Water Molecule Excitation at O K-edge 1 b 1 3 a 1 Method/Basis SOS (HF)/D 95 V+* 1 b 2 Ref: J. -H. Guo et al. Phys. Rev. Lett. , Vol 89, 137402 (2002). HF/Sadlej using Dalton program
Methyl Alcohol HF/Sadlej
Fluorescence Spectra of Methyl Alcohol
Theoretical Challenges of Femtosecond X-Ray Simulations Time Resolved Geometry Changes Ø Immediately after electronic excitation, the molecule will begin to relax to a new equilibrium structure. This can involve: • photodissociation • changes in conformation • vibrations Ø Fast codes for excited state dynamics. Ø Codes for computing current profiles within molecules. Ø Simulate quantum molecular dynamics to find forces on atoms in excited state. Ø Use mixed quantum/classical molecular dynamics for solvent. Ø Study of X-ray fluorescence and four wave mixing when the molecule is initially in the optically excited state.
- Slides: 23