First-principles Time-dependent Density Functional Theory Research Articles

Extending the limit of LR-TDDFT on two different approaches: Numerical algorithms and new Sunway heterogeneous supercomputer

First-principles time-dependent density functional theory (TDDFT) is a powerful tool to accurately describe the excited-state properties of molecules and solids in condensed matter physics, computational chemistry, and materials science. However, a perceived drawback in TDDFT calculations is its ultrahigh computational cost O(N5∼N6) and large memory usage O(N4) especially for plane-wave basis set, confining its applications to large systems containing thousands of atoms. Here, we present a massively parallel implementation of linear-response TDDFT (LR-TDDFT) and accelerate LR-TDDFT in two different aspects: (1) numerical algorithms on the X86 supercomputer and (2) optimizations on the heterogeneous architecture of the new Sunway supercomputer. Furthermore, we carefully design the parallel data and task distribution schemes to accommodate the physical nature of different computation steps. By utilizing these two different methods, our implementation can gain an overall speedup of 10x and 80x and efficiently scales to large systems up to 4096 and 2744 atoms within dozens of seconds.

Solid-state high harmonic spectroscopy for all-optical band structure probing of high-pressure quantum states

High pressure has triggered various novel states/properties in condensed matter, as the most representative and dramatic example being near-room-temperature superconductivity in highly pressured hydrides (~200 GPa). However, the mechanism of superconductivity is not confirmed, due to the lacking of effective approach to probe the electronic band structure under such high pressures. Here, we theoretically propose that the band structure and electron-phonon coupling (EPC) of high-pressure quantum states can be probed by solid-state high harmonic generation (sHHG). This strategy is investigated in high-pressure Im-3m H3S by the state-of-the-art first-principles time-dependent density-functional theory simulations, where the sHHG is revealed to be strongly dependent on the electronic structures and EPC. The dispersion of multiple bands near the Fermi level is effectively retrieved along different momentum directions. Our study provides unique insights into the potential all-optical route for band structure and EPC probing of high-pressure quantum states, which is expected to be helpful for the experimental exploration of high-pressure superconductivity in the future.

Improving dynamic collision frequencies: Impacts on dynamic structure factors and stopping powers in warm dense matter

Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, we systematically investigate variations in dynamic electron–ion collision frequencies ν(ω) in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in ν(ω). These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium.

Plasmonic features of a hybrid system comprising borophene and aluminium

Due to their intriguing light–matter interactions at the nanoscale, hybrid systems made up of dissimilar elements are currently receiving a lot of attention. In this study, we systematically examine the optical and plasmonic properties of a hybrid system made of borophene and Aluminium (Al) film using first-principles time-dependent density functional theory calculations. We discovered that, despite the weak charge transfer between the constituent systems, the borophene plasmon is totally quenched or damped by the proximity to the Al film for a gap distance (d) below 9 Å, and the Al plasmon dominates the optical response. When the d value is raised to 9 Å and higher, the borophene plasmon starts to reemerge. This finding is consistent with earlier finding on graphene on metallic film, which demonstrated that the graphene plasmons are quenched by the metallic film’s proximity and recovered at high d values.

Graphene Quantum Dot with Divacancy and Topological Defects: A Novel Material for Promoting Prompt and Delayed Fluorescence of Tunable Wavelengths

This work demonstrates the unique approach of introducing divacancy imperfections in topological Stone-Wales type defected graphene quantum dots for harvesting both singlet and triplet excitons, essential for fabricating fluorescent organic light-emitting diodes. Here, we first reveal that structural relaxation of these systems establishes the high-spin triplet state as the stable ground state at room temperature, thereby significantly increasing their potential in designing spintronic devices. Our extensive electron-correlated computations then demonstrate that the energetic ordering of the singlet and triplet states in these relaxed structures can trigger both prompt and delayed fluorescence of different wavelengths through various decay channels. Particularly, the position of divacancy determines the tunability range of the emission wavelengths. In addition, our results obtained from both multi-reference singles-doubles configurationinteraction (MRSDCI) and first-principles time-dependent density functional theory (TDDFT) methodologies highlight that the synergetic effects of divacancy-position,structural relaxation and spin multiplicity critically govern the nature and magnitude of shift exhibited by the most intense peak of the absorption profile, crucial for designing optoelectronic devices.

Effect of spin-orbit coupling on the high harmonics from the topological Dirac semimetal Na3Bi

In this work, we performed extensive first-principles simulations of high-harmonic generation in the topological Diract semimetal Na3Bi using a first-principles time-dependent density functional theory framework, focusing on the effect of spin-orbit coupling (SOC) on the harmonic response. We also derived an analytical model describing the microscopic mechanism of strong-field dynamics in presence of spin-orbit coupling, starting from a locally U(1) × SU(2) gauge-invariant Hamiltonian. Our results reveal that SOC: (i) affects the strong-field excitation of carriers to the conduction bands by modifying the bandstructure of Na3Bi, (ii) makes each spin channel reacts differently to the driven laser by modifying the electron velocity (iii) changes the emission timing of the emitted harmonics. Moreover, we show that the SOC affects the harmonic emission by directly coupling the charge current to the spin currents, paving the way to the high-harmonic spectroscopy of spin currents in solids.

First-principles method for nonlinear light propagation at oblique incidence.

We have developed a computational method to describe the nonlinear light propagation of an intense and ultrashort pulse at oblique incidence on a flat surface. In the method, coupled equations of macroscopic light propagation and microscopic electron dynamics are simultaneously solved using a multiscale modeling. The microscopic electronic motion is described by first-principles time-dependent density functional theory. The macroscopic Maxwell equations that describe oblique light propagation are transformed into one-dimensional wave equations. As an illustration of the method, light propagation at oblique incidence on a silicon thin film is presented.

Light-Induced Charge Transfer from Transition-Metal-Doped Aluminum Clusters to Carbon Dioxide.

Charge transfer between molecules and catalysts plays a critical role in determining the efficiency and yield of photochemical catalytic processes. In this paper, we study light-induced electron transfer between transition-metal-doped aluminum clusters and CO2 molecules using first-principles time-dependent density-functional theory. Specifically, we carry out calculations for a range of dopants (Zr, Mn, Fe, Ru, Co, Ni, and Cu) and find that the resulting systems fall into two categories: Cu- and Fe-doped clusters exhibit no ground-state charge transfer, weak CO2 adsorption, and light-induced electron transfer into the CO2. In all other systems, we observe ground-state electron transfer into the CO2 resulting in strong adsorption and predominantly light-induced electron back-transfer from the CO2 into the cluster. These findings pave the way toward a rational design of atomically precise aluminum photocatalysts.

First-principles study of ultrafast and nonlinear optical properties of graphite thin films

We theoretically investigate ultrafast and nonlinear optical properties of graphite thin films based on first-principles time-dependent density functional theory. We first calculate electron dynamics in a unit cell of graphite under a strong pulsed electric field and explore the transient optical properties of graphite. The optical response of graphite shows a sudden change from the conducting to the insulating phase at a certain intensity range of the applied electric field. It also appears to be a saturable absorption (SA) in the energy transfer from the electric field to electrons. We next investigate a light propagation in graphite thin films by solving the coupled dynamics of the electrons and the electromagnetic fields simultaneously. It is observed that the SA manifests in the propagation with small attenuation in the spatial region where the electric field amplitude is about $4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}2}\phantom{\rule{4.pt}{0ex}}\text{to}\phantom{\rule{4.pt}{0ex}}7\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}2}\phantom{\rule{0.16em}{0ex}}\mathrm{V}/\AA{}$.

Excited-State Charge Distribution of a Donor-π-Acceptor Zn Porphyrin Probed by N K-Edge Transient Absorption Spectroscopy.

Zinc porphyrin solar cell dyes with donor-π-acceptor architectures combine light absorber (π), electron-donor, and electron-acceptor moieties inside a single molecule with atomic precision. The donor-π-acceptor design promotes the separation of charge carriers following optical excitation. Here, we probe the excited-state electronic structure within such molecules by combining time-resolved X-ray absorption spectroscopy at the N K-edge with first-principles time-dependent density functional theory (TD-DFT) calculations. Customized Zn porphyrins with strong-donor triphenylamine groups or weak-donor tri-tert-butylbenzene groups were synthesized. Energetically well-separated N K-edge absorption features simultaneously probe the excited-state electronic structure from the perspectives of the macrocycle and triphenylamine N atoms. New absorption transitions between the macrocycle N atoms and the excited-state HOMO vacancy are observed, and the triphenylamine associated absorption feature blue-shifts, consistent with partial oxidation of the donor groups in the excited state.

Carrier Dynamics and Transfer across the CdS/MoS2 Interface upon Optical Excitation.

Carrier dynamics across the interface of heterostructures have important technological, photovoltaic, and catalytic implications. Using first-principles time-dependent density functional theory, we have systematically investigated the charge transfer of excited carriers from CdS to MoS2 and found that two interdependent mechanisms are responsible for the transfer, one slow and one fast. While the slower process may be attributed to typical electron-phonon coupling, the interfacial dipole resulting from this transfer enables a fast secondary process involving a level crossing of the excited carrier state in CdS with receiving states in MoS2. An analysis based on the interfacial binding energy reveals that the Cd-terminated (001) interface is by far the most energetically favorable, which in addition to its calculated fast resonant electron transfer suggests it is a good candidate to explain the experimentally observed charge transfer between CdS and MoS2.

Pariser–Parr–Pople Model Based Configuration-Interaction Study of Linear Optical Absorption in Lower-Symmetry Polycyclic Aromatic Hydrocarbon Molecules

The electronic and optical properties of various polycyclic aromatic hydrocarbons (PAHs) with lower symmetry, namely, benzo[ghi]perylene (C22H12), benzo[a]coronene (C28H14), naphtho[2,3a]coronene (C32H16), anthra[2,3a]coronene (C36H18), and naphtho[8,1,2-abc]coronene (C30H14), were investigated. We performed electron-correlated calculations using screened and standard parameters in the π-electron Pariser–Parr–Pople (PPP) Hamiltonian, and the correlation effects were included, both for ground and excited states, using multireference singles–doubles configuration-interaction (MRSDCI) methodology. The PPP model Hamiltonian includes long-range Coulomb interactions, which increase the accuracy of our calculations. The results of our calculations predict that, with the increasing sizes of the coronene derivatives, optical spectra are red-shifted, and the optical gaps decrease. In each spectrum, the first peak representing the optical gap is of moderate intensity, while the more intense peaks appear at higher energies. Our computed spectra are in good agreement with the available experimental data. For the purpose of comparison, we also performed first-principles time-dependent density-functional theory (TDDFT) calculations of the optical gaps of these molecules using Gaussian basis functions and found that they yielded values lower than our configuration-interaction (CI) results.

Microscopic Theory of Plasmons in Substrate-Supported Borophene.

We compute the dielectric properties of freestanding and metal-supported borophene from first-principles time-dependent density functional theory. We find that both the low- and high-energy plasmons of borophene are fully quenched by the presence of a metallic substrate at borophene-metal distances smaller than ≃9 Å. Based on these findings, we derive an electrodynamic model of the interacting, momentum-dependent polarizability for a two-dimensional metal on a model metallic substrate, which quantitatively captures the evolution of the dielectric properties of borophene as a function of metal-borophene distance. Applying this model to a series of metallic substrates, we show that maximizing the plasmon energy detuning between borophene and substrate is the key material descriptor for plasmonic performance.

Photoabsorption Spectra of Aqueous Oxyluciferin Anions Elucidated by Explicit Quantum Solvent.

Experimental photoabsorption spectra of three possible isomers (phenolate-keto, phenolate-enol, and phenol-enolate) of oxyluciferin anions in aqueous solution were reproduced by first-principles time-dependent density functional theory simulations in which the entire system including the oxyluciferin anion and 64 water molecules were modeled by full quantum mechanics (full QM), unlike the conventional hybrid method, where the surrounding water molecules are modeled by molecular mechanics (MM) or a continuum solvent model. The full QM photoabsorption spectra were calculated from 1000 structures that had been obtained using the first-principles Born-Oppenheimer molecular dynamics simulations, which included the van der Waals correction, to take into account the effect of dynamical fluctuations of the hydration structure. The full QM calculation with CAM-B3LYP functional, which is the most elaborate one and is apparently the most consistent with experiment, is compared to others obtained with different levels of the functional and the solvent model. The amount of charge leakage from the oxyluciferin anions to the aqueous solution is found to differ significantly between the ground and excited states and is strongly dependent on the simulation method. The conventional solvent models do not take this into account, but the QM/MM can do it appropriately when including more than 10 water molecules into the QM region.

Heavy ion ranges from first-principles electron dynamics

The effects of incident energetic particles, and the modification of materials under irradiation, are governed by the mechanisms of energy losses of ions in matter. The complex processes affecting projectiles spanning many orders of magnitude in energy depend on both ion and electron interactions. Developing multi-scale modeling methods that correctly capture the relevant processes is crucial for predicting radiation effects in diverse conditions. In this work, we obtain channeling ion ranges for tungsten, a prototypical heavy ion, by explicitly simulating ion trajectories with a method that takes into account both the nuclear and the electronic stopping power. The electronic stopping power of self-ion irradiated tungsten is obtained from first-principles time-dependent density functional theory (TDDFT). Although the TDDFT calculations predict a lower stopping power than SRIM by a factor of three, our result shows very good agreement in a direct comparison with ion range experiments. These results demonstrate the validity of the TDDFT method for determining electronic energy losses of heavy projectiles, and in turn its viability for the study of radiation damage.

Nonlinear polarization evolution using time-dependent density functional theory.

We propose a theoretical and computational approach to investigate temporal behavior of a nonlinear polarization in a perturbative regime induced by an intense and ultrashort pulsed electric field. First-principles time-dependent density functional theory is employed to describe the electron dynamics. Temporal evolution of third-order nonlinear polarization is extracted from a few calculations of electron dynamics induced by pulsed electric fields with the same time profile but different amplitudes. We discuss characteristic features of the nonlinear polarization evolution as well as an extraction of nonlinear susceptibilities and time delays by fitting the polarization. We also carry out a decomposition of temporal and spatial changes in the electron density in power series with respect to the field amplitude. It helps to get insight into the origin of the nonlinear polarization in atomic scale.

Photochemical Etching of Carbonyl Groups from a Carbon Matrix: The (001) Diamond Surface.

The surface of diamond is reported to undergo nonablative photochemical etching when exposed to ultraviolet (UV) radiation which allows controlled single and partial layer removal of lattice layers. Oxygen termination of surface dangling bonds is known to be crucial for the etching process; however, the exact mechanism of carbon ejection remains unclear. We investigate the interaction of UV laser pulses with oxygen-terminated diamond surfaces using atomic-scale surface characterization combined with first-principles time-dependent density functional theory calculations. We present evidence for laser-induced desorption (LID) from carbonyl functional groups at the diamond {001} surface. The doubly bonded carbonyl group is photoexcited into a triply bonded CO-like state, including scission of the underlying C─C bonds. The carbon removal process in LID is atom by atom; therefore, this mechanism provides a novel "top-down" approach for creating nanostructures on the surface of diamond and other carbon-containing semiconductors.

A first-principles simulation method for ultra-fast nano-optics

We develop a computational approach for ultrafast nano-optics based on first-principles time-dependent density functional theory. Solving Maxwell equations for light propagation and time-dependent Kohn-Sham equation for electron dynamics simultaneously, intense and ultrashort laser pulse interaction with a dielectric nano-structure is described taking full account of nonlinear effects. As an illustrative example, irradiation of a pulsed light on silicon nano-sphere system is presented.

Time-dependent density functional theory for interaction of ultrashort light pulse with thin materials

We present a comprehensive theoretical description for an irradiation of an ultrashort light pulse normally on thin materials based on first-principles time-dependent density functional theory. As the most elaborate scheme, we develop a microscopic description solving Maxwell equations for light electromagnetic fields and the time-dependent Kohn-Sham equation for electron dynamics simultaneously in the time domain using a common spatial grid. We call it the microscopic Maxwell-TDDFT scheme. We test this scheme for silicon thin films of various thickness, from a few atomic layers to a few tens of nm. We show that the microscopic Maxwell-TDDFT scheme provides a satisfactory description incorporating the electronic structure of thin films in the first-principles level, multiple reflections of the electromagnetic fields at the surfaces, and nonlinear light-matter interaction when the incident light pulse is strong. However, the calculation becomes expensive as the thickness increases. We then consider two limiting cases of extremely thin and sufficiently thick films and develop approximate schemes. For the extremely thin case including two-dimensional atomic-layered materials, a two-dimensional macroscopic electromagnetism is developed: a two-dimensional susceptibility is introduced for a weak field, while time evolution equation is derived for an intense field. For sufficiently thick films, the microscopic Maxwell-TDDFT scheme is expected to coincide with a description utilizing ordinary macroscopic electromagnetism. We numerically confirm it comparing the calculated results: For a weak field, a comparison is made with a description using the bulk dielectric susceptibility. For a strong field, a comparison is made with a multiscale Maxwell-TDDFT scheme which the authors' group developed previously.

SALMON: Scalable Ab-initio Light–Matter simulator for Optics and Nanoscience

SALMON (Scalable Ab-initio Light–Mattersimulator for Optics and Nanoscience, http://salmon-tddft.jp) is a software package for the simulation of electron dynamics and optical properties of molecules, nanostructures, and crystalline solids based on first-principles time-dependent density functional theory. The core part of the software is the real-time, real-space calculation of the electron dynamics induced in molecules and solids by an external electric field solving the time-dependent Kohn–Sham equation. Using a weak instantaneous perturbing field, linear response properties such as polarizabilities and photoabsorptions in isolated systems and dielectric functions in periodic systems are determined. Using an optical laser pulse, the ultrafast electronic response that may be highly nonlinear in the field strength is investigated in time domain. The propagation of the laser pulse in bulk solids and thin films can also be included in the simulation via coupling the electron dynamics in many microscopic unit cells using Maxwell’s equations describing the time evolution of the electromagnetic fields. The code is efficiently parallelized so that it may describe the electron dynamics in large systems including up to a few thousand atoms. The present paper provides an overview of the capabilities of the software package showing several sample calculations. Program summaryProgram Title: SALMON: Scalable Ab-initio Light–Matter simulator for Optics and NanoscienceProgram Files doi:http://dx.doi.org/10.17632/8pm5znxtsb.1Licensing provisions: Apache-2.0Programming language: Fortran 2003Nature of problem: Electron dynamics in molecules, nanostructures, and crystalline solids induced by an external electric field is calculated based on first-principles time-dependent density functional theory. Using a weak impulsive field, linear optical properties such as polarizabilities, photoabsorptions, and dielectric functions are extracted. Using an optical laser pulse, the ultrafast electronic response that may be highly nonlinear with respect to the exciting field strength is described as well. The propagation of the laser pulse in bulk solids and thin films is considered by coupling the electron dynamics in many microscopic unit cells using Maxwell’s equations describing the time evolution of the electromagnetic field.Solution method: Electron dynamics is calculated by solving the time-dependent Kohn–Sham equation in real time and real space. For this, the electronic orbitals are discretized on a uniform Cartesian grid in three dimensions. Norm-conserving pseudopotentials are used to account for the interactions between the valence electrons and the ionic cores. Grid spacings in real space and time, typically 0.02 nm and 1 as respectively, determine the spatial and temporal resolutions of the simulation results. In most calculations, the ground state is first calculated by solving the static Kohn–Sham equation, in order to prepare the initial conditions. The orbitals are evolved in time with an explicit integration algorithm such as a truncated Taylor expansion of the evolution operator, together with a predictor–corrector step when necessary. For the propagation of the laser pulse in a bulk solid, Maxwell’s equations are solved using a finite-difference scheme. By this, the electric field of the laser pulse and the electron dynamics in many microscopic unit cells of the crystalline solid are coupled in a multiscale framework.