Prevalence of the adiabatic exchange-correlation potential approximation in time-dependent density functional theory

On the Prediction of Spectroscopic Fingerprints of Co2+ Complexes Relevant for the ZIF Nucleation Process

The nucleation process of zeolitic imidazolate frameworks (ZIFs) is to date not completely understood. Recently, it has been found that, during the formation of Co-ZIF-67, after mixing imidazole-type ligands with octahedral precursors containing oxygen-coordinated ligands, a metal-organic pool with a diversity of transition metal complexes (TMCs) is formed showing fingerprints of octahedral and tetrahedral Co2+ complexes with both types of ligands [Filez, M. Cell Rep. Phys. Sci. 2021, 2, 100680]. In order to further unravel this process, we aim to characterize the d-d transitions of the TMCs and focus on their number, intensity, and position, which change during the process and can thus serve as a fingerprint for the formed species. It was previously shown that the number of ligands and symmetry has a detrimental influence on the ground state properties of Co2+ TMCs. Herein, we investigate how far excited state properties of TMCs relevant during nanoporous formation processes can be predicted by time-dependent density functional theory (TDDFT) and ligand field density functional theory (LFDFT). As TMCs are known to be challenging systems with possibly degenerate ground states and double excitations, we first investigate the performance of both techniques on first-row octahedral aqua-complexes. With this knowledge, we then focus on tetrahedral Co2+ complexes with aqua and imidazole-type ligands in order to investigate in how far we can propose a spectroscopic fingerprint that allows us to follow the Co2+ complexes during the formation of Co-ZIF-67. The results of TDDFT and LFDFT are qualitatively in agreement and provide complementary information. We found that various features can be used to distinguish between the species. However, as LFDFT is not suited for TMCs possessing the extended imidazole-type ligands and double and spin-flip states are not included in TDDFT, both techniques need to be complemented with more advanced methods to obtain complete insight into the d-d excitations of TMCs with imidazole ligands. Therefore, we particularly explored ab initio ligand field theory, which is capable of describing double excitations and is, in contrast to LFDFT, suitable for TMCs with extended ligands.

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree-Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

High-frequency limit of spectroscopy

We consider an arbitrary quantum mechanical system, initially in its ground-state, exposed to a time-dependent electromagnetic pulse with a carrier frequency $\omega_0$ and a slowly varying envelope of finite duration. By working out a solution to the time-dependent Schr\"odinger equation in the high-$\omega_0$ limit, we find that, to the leading order in $\omega_0^{-1}$, a perfect self-cancellation of the system's linear response occurs as the pulse switches off. Surprisingly, the system's observables are, nonetheless, describable in terms of a combination of its linear density response function and nonlinear functions of the electric field. An analysis of jellium slab and jellium sphere models reveals a very high surface sensitivity of the considered setup, producing a richer excitation spectrum than accessible within the conventional linear response regime. On this basis, we propose a new spectroscopic technique, which we provisionally name the Nonlinear High-Frequency Pulsed Spectroscopy (NLHFPS). Combining the advantages of the extraordinary surface sensitivity, the absence of constraints by the traditional dipole selection rules, and the clarity of theoretical interpretation utilizing the linear response time-dependent density functional theory, NLHFPS has a potential to evolve into a powerful characterization method for nanoscience and nanotechnology.

Real-time time-dependent density functional theory using density fitting and the continuous fast multipole method

An implementation of real-time time-dependent density functional theory (RT-TDDFT) within the TURBOMOLE program package is reported using Gaussian-type orbitals as basis functions, second and fourth order Magnus propagator, and the self-consistent field as well as the predictor-corrector time integration schemes. The Coulomb contribution to the Kohn-Sham matrix is calculated combining density fitting approximation and the continuous fast multipole method. Performance of the implementation is benchmarked for molecular systems with different sizes and dimensionalities. For linear alkane chains, the wall time for density matrix time propagation step is comparable to the Kohn-Sham (KS) matrix construction. However, for larger two- and three-dimensional molecules, with up to about 5,000 basis functions, the computational effort of RT-TDDFT calculations is dominated by the KS matrix evaluation. In addition, the maximum time step is evaluated using a set of small molecules of different polarities. The photoabsorption spectra of several molecular systems calculated using RT-TDDFT are compared to those obtained using linear response time-dependent density functional theory and coupled cluster methods.

Collision dynamics of proton with formaldehyde: fragmentation and ionization

Using time-dependent density functional theory, applied to the valence electrons and coupled non-adiabatically to molecular dynamics of the ions, we study the ionization and fragmentation of formaldehyde in collision with a proton. Four different impact energies: 35 eV, 85 eV, 135 eV, and 300 eV are chosen in order to study the energy effect in the low energy region, and ten different incident orientations at 85 eV are considered for investigating the steric effect. Fragmentation ratios, single, double, and total electron ionization cross sections are calculated. For large impact parameters, these results are close to zero irrespective of the incident orientations due to a weak projectile-target interaction. For small impact parameters, the results strongly depend on the collision energy and orientation. We also give the kinetic energy releases and scattering angles of protons, as well as the cross section of different ion fragments and the corresponding reaction channels.

Spatial distribution of electron densities during optical excitation of C60

We calculate the spatial distribution of the oscillating electron density for the low energy optical excitations in C(60), using ab initio time-dependent density functional theory. We contrast this with the unscreened free response of the Kohn-Sham eigenfunctions, in which the charge transfer is between hemispheres and the simple density oscillation patterns are very similar for electrons outside and inside the C(60) radius. In the screened case, the oscillation patterns involve charge transfer on the scale of interatomic distances. The exterior and interior densities play different roles in the oscillator strength and screening. Almost all of the dipole oscillation for optical absorption strength is from exterior electrons, a consequence of the superposition of π and σ electron oscillations. The lowest transition with appreciable strength is dominated by screening from π electrons.

Time-dependent density functional theory for open quantum systems with unitary propagation

We extend the Runge-Gross theorem for a very general class of open quantum systems under weak assumptions about the nature of the bath and its coupling to the system. We show that for Kohn-Sham (KS) time-dependent density functional theory, it is possible to rigorously include the effects of the environment within a bath functional in the KS potential. A Markovian bath functional inspired by the theory of nonlinear Schrödinger equations is suggested, which can be readily implemented in currently existing real-time codes. Finally, calculations on a helium model system are presented.