We investigate the spectroscopy and photoinduced electron dynamics within the conduction band of reduced rutile TiO 2 (110) surface by multiphoton photoemission (mPP) spectroscopy with wavelength tunable ultrafast (~20 fs) laser pulse excitation. Tuning the mPP photon excitation energy between 2.9 and 4.6 eV reveals a nearly degenerate pair of new unoccupied states located at 2.73±0.05 and 2.85±0.05 eV above the Fermi level, which can be analyzed through the polarization and sample azimuthal orientation dependence of the mPP spectra. Based on the calculated electronic structure and optical transition moments, as well as related spectroscopic evidence, we assign these resonances to transitions between Ti 3d-bands of nominally t 2g and e g symmetry, which are split by crystal-field. The initial states for the optical transition are the reduced Ti 3+ states of t 2g symmetry populated by formation oxygen vacancy defects, which exist within the band gap of TiO 2 . Furthermore, we studied the electron dynamics within the conduction band of TiO 2 by three-dimensional (3D) time-resolved pump-probe interferometric mPP measurements. The spectroscopic and time-resolved studies reveal competition between 2PP and 3PP processes where the t 2g -e g transitions in the 2PP process saturate, and are overtaken by the 3PP process initiated by the band gap excitation from the valence band of TiO 2 .

Within density-functional theory, perturbation theory(PT) is the state-of-the-art formalism for assessing the response to homogeneous electric fields and the associated material properties, e.g., polarizabilities, dielectric constants, and Raman intensities. Here, we derive a real-space formulation of PT and present an implementation within the all-electron, numeric atom-centered orbitals electronic structure code FHI-aims that allows for massively parallel calculations. As demonstrated by extensive validation, we achieve a rapid computation of accurate response properties of molecules and solids. As an application showcase, we present harmonic and anharmonic Raman spectra, the latter obtained by combining hundreds of thousands of PT calculations with ab initio molecular dynamics. By using the PBE exchange-correlation functional with many-body van der Waals corrections, we obtain spectra in good agreement with experiment especially with respect to lineshapes for the isolated paracetamol molecule and two polymorphs of the paracetamol crystal.

HONPAS is an ab initio electronic structure program for linear scaling or O(N) first-principles calculations of large and complex systems using standard norm-conserving pseudopotentials, numerical atomic orbitals (NAOs) basis sets, and periodic boundary conditions. HONPAS is developed in the framework of the SIESTA methodology and focuses on the development and implementation of efficient O(N) algorithms for ab initio electronic structure calculations. The Heyd-Scuseria-Ernzerhof (HSE) screened hybrid density functional has been implemented using a NAO2GTO scheme to evaluate the electron repulsion integrals (ERIs) with NAOs. ERI screening techniques allow the HSE functional calculations to be very efficient and scale linearly. The density matrix purification algorithms have been implemented, and the PSUTC2 and SUTC2 methods have been developed to deal with spin unrestricted systems with or without predetermined spin multiplicity, respectively. After the self-consistent field (SCF) process, additional O(N) post-SCF calculations for frontier molecular orbitals and maximally localized Wannier functions are also developed and implemented. Finally, an O(N) method based on the density matrix perturbation theory has been proposed and implemented to treat electric field in solids. This article provides an overall introduction to capabilities of HONPAS and implementation details of different O(N) algorithms.

A real-space formalism for density-functional perturbation theory (DFPT) is derived and applied for the computation of harmonic vibrational properties in molecules and solids. The practical implementation using numeric atom-centered orbitals as basis functions is demonstrated exemplarily for the all-electron Fritz Haber Institute ab initio molecular simulations (FHI-aims) package. The convergence of the calculations with respect to numerical parameters is carefully investigated and a systematic comparison with finite-difference approaches is performed both for finite (molecules) and extended (periodic) systems. Finally, the scaling tests and scalability tests on massively parallel computer systems demonstrate the computational efficiency