Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-Time Simulation of Ultrafast Electronic Dynamics of Nanoscale Systems Involving an Organic Molecule and a Nanoparticle Dimer

Understanding and predicting the behavior of nanomaterials composed of plasmons interacting with quantum emitters at ultrafast timescales is crucial for the better manipulation of light at the nanoscale and advancing technologies like ultrafast communication and computing. Here we perform a simulation of the "real-time" electronic dynamics of a coupled molecule-metal nanoparticle dimer interacting with an ultrashort resonant laser pulse by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency-dependent fluctuating charge (TD-ωFQ) model, an atomistic electromagnetic (AEM) model for the dynamic plasmonic response of nanoparticles. It is shown that the induced dipoles evolve from an exponential decay pattern to a beat pattern with an increase in coupling strength, which is altered by changing the molecular orientation relative to the dimer axis. It is further shown that in the strong coupling regime, both the excited molecule and the plasmon relax rapidly due to the molecule-plasmon interaction, and the efficient coherent energy exchange between the interacting molecule and plasmon modes occurs on a femtosecond (fs) timescale. This work provides guidance on manipulating light-matter interaction and studying molecular plasmonics at extremely fast timescales.

Time-dependent Kohn-Sham electron dynamics coupled with nonequilibrium plasmonic response via atomistic electromagnetic model

Computational modeling of plasmon-mediated molecular photophysical and photochemical behaviors can help us better understand and tune the bound molecular properties and reactivity and make better decisions to design and control nanostructures. However, computational investigations of coupled plasmon-molecule systems are challenging due to the lack of accurate and efficient protocols to simulate these systems. Here, we present a hybrid scheme by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency dependent fluctuating charge (TD-ωFQ) model. At first, we transform ωFQ in the frequency-domain, an atomistic electromagnetic model for the plasmonic response of plasmonic metal nanoparticles (PMNPs), into the time-domain and derive its equation-of-motion formulation. The TD-ωFQ introduces the nonequilibrium plasmonic response of PMNPs and atomistic interactions to the electronic excitation of the quantum mechanical (QM) region. Then, we combine TD-ωFQ with RT-TDDFT. The derived RT-TDDFT/TD-ωFQ scheme allows us to effectively simulate the plasmon-mediated "real-time" electronic dynamics and even the coupled electron-nuclear dynamics by combining them with the nuclear dynamics approaches. As a first application of the RT-TDDFT/TD-ωFQ method, we study the nonradiative decay rate and plasmon-enhanced absorption spectra of two small molecules in the proximity of sodium MNPs. Thanks to the atomistic nature of the ωFQ model, the edge effect of MNP on absorption enhancement has also been investigated and unveiled.

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge-Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6-0.7 bohr in OPAW-TDDFT rather than the 0.4-0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe-Salpeter equation.

Tracking electron motion within and outside of Floquet bands from attosecond pulse trains in time-resolved ARPES

Floquet engineering has recently emerged as a technique for controlling material properties with light. Floquet phases can be probed with time- and angle-resolved photoelectron spectroscopy (Tr-ARPES), providing direct access to the laser-dressed electronic bands. Applications of Tr-ARPES to date focused on observing the Floquet-Bloch bands themselves, and their build-up and dephasing on sub-laser-cycle timescales. However, momentum and energy resolved sub-laser-cycle dynamics between Floquet bands have not been analyzed. Given that Floquet theory strictly applies in time-periodic conditions, the notion of resolving sub-laser-cycle dynamics between Floquet states seems contradictory-it requires probe pulse durations below a laser cycle that inherently cannot discern the time-periodic nature of the light-matter system. Here we propose to employ attosecond pulse train probes with the same temporal periodicity as the Floquet-dressing pump pulse, allowing both attosecond sub-laser-cycle resolution and a proper projection of Tr-ARPES spectra on the Floquet-Bloch bands. We formulate and employ this approach inab-initiocalculations in light-driven graphene. Our calculations predict significant sub-laser-cycle dynamics occurring within the Floquet phase with the majority of electrons moving within and in-between Floquet bands, and a small portion residing and moving outside of them in what we denote as 'non-Floquet' bands. We establish that non-Floquet bands arise from the pump laser envelope that induces non-adiabatic electronic excitations during the pulse turn-on and turn-off. By performing calculations in systems with poly-chromatic pumps we also show that Floquet states are not formed on a sub-laser-cycle level. This work indicates that the Floquet-Bloch states are generally not a complete basis set for sub-laser-cycle dynamics in steady-state phases of matter.

Optical properties and exciton transfer between N-heterocyclic carbene iridium(III) complexes for blue light-emitting diode applications from first principles

N-heterocyclic carbene (NHC) iridium(III) complexes are considered as promising candidates for blue emitters in organic light-emitting diodes. They can play the roles of the emitter as well as of electron and hole transporters in the same emission layer. We investigate optical transitions in such complexes with account of geometry and electronic structure changes upon excitation or charging and exciton transfer between the complexes from first principles. It is shown that excitation of NHC iridium complexes is accompanied by a large reorganization energy ∼0.7 eV and a significant loss in the oscillator strength, which should lead to low exciton diffusion. Calculations with account of spin-orbit coupling reveal a small singlet-triplet splitting ∼0.1 eV, whereas the oscillator strength for triplet excitations is found to be an order of magnitude smaller than for the singlet ones. The contributions of the Förster and Dexter mechanisms are analyzed via the explicit integration of transition densities. It is shown that for typical distances between emitter complexes in the emission layer, the contribution of the Dexter mechanism should be negligible compared to the Förster mechanism. At the same time, the ideal dipole approximation, although giving the correct order of the exciton coupling, fails to reproduce the result taking into account spatial distribution of the transition density. For charged NHC complexes, we find a number of optical transitions close to the emission peak of the blue emitter with high exciton transfer rates that can be responsible for exciton-polaron quenching. The nature of these transitions is analyzed.

Tracking Charge Migration with Frequency-Matched Strobo-Spectroscopy

We present frequency-matched strobo-spectroscopy (FMSS) of charge migration (CM) in bromobutadiyne, simulated with time-dependent density functional theory. CM + FMSS is a pump-probe scheme that uses a frequency-matched high harmonic generation (HHG)-driving laser as an independent probe step, following the creation of a localized hole on the bromine atom that induces CM dynamics. We show that the delay-dependent harmonic yield tracks the phase of the CM dynamics through its sensitivity to the amount of electron density on the bromine end of the molecule. FMSS takes advantage of the intrinsic attosecond time resolution of the HHG process in which different harmonics are emitted at different times and thus probe different locations of the electron hole. Finally, we show that the CM-induced modulation of the HHG signal is dominated by the recombination step of the HHG process, with a negligible contribution from the ionization step.

Mapping Light-Dressed Floquet Bands by Highly Nonlinear Optical Excitations and Valley Polarization

Ultrafast nonlinear optical phenomena in solids have been attracting a great deal of interest as novel methodologies for the femtosecond spectroscopy of electron dynamics and control of the properties of materials. Here, we theoretically investigate strong-field nonlinear optical transitions in a prototypical two-dimensional material, hBN, and show that the k-resolved conduction band charge occupation patterns induced by an elliptically polarized laser can be understood in a multiphoton resonant picture, but, remarkably, only if using the Floquet light-dressed states instead of the undressed matter states. Our work demonstrates that Floquet dressing affects ultrafast charge dynamics and photoexcitation even from a single pump pulse and establishes a direct measurable signature for band dressing in nonlinear optical processes in solids, opening new paths for ultrafast spectroscopy and valley manipulation.

Velocity-Gauge Real-Time Time-Dependent Density Functional Tight-Binding for Large-Scale Condensed Matter Systems

We present a new velocity-gauge real-time, time-dependent density functional tight-binding (VG-rtTDDFTB) implementation in the open-source DFTB+ software package (https://dftbplus.org) for probing electronic excitations in large, condensed matter systems. Our VG-rtTDDFTB approach enables real-time electron dynamics simulations of large, periodic, condensed matter systems containing thousands of atoms with a favorable computational scaling as a function of system size. We provide computational details and benchmark calculations to demonstrate its accuracy and computational parallelizability on a variety of large material systems. As a representative example, we calculate laser-induced electron dynamics in a 512-atom amorphous silicon supercell to highlight the large periodic systems that can be examined with our implementation. Taken together, our VG-rtTDDFTB approach enables new electron dynamics simulations of complex systems that require large periodic supercells, such as crystal defects, complex surfaces, nanowires, and amorphous materials.

Phonon-enhanced nonlinearities in hexagonal boron nitride

Polar crystals can be driven into collective oscillations by optical fields tuned to precise resonance frequencies. As the amplitude of the excited phonon modes increases, novel processes scaling non-linearly with the applied fields begin to contribute to the dynamics of the atomic system. Here we show two such optical nonlinearities that are induced and enhanced by the strong phonon resonance in the van der Waals crystal hexagonal boron nitride (hBN). We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies beyond the third.

Defect-Mediated Atomic Layer Etching Processes on Cl-Si(100): An Atomistic Insight

Defects play a significant role in atomic layer etching (ALE) processes; however, a fundamental understanding at the atomic level is still lacking. To bridge this knowledge gap, this study investigated the role of point defects in the laser-induced ALE of Cl-Si(100) using density functional theory (DFT) and real-time time-dependent DFT calculations. In the calculations, both the pristine surface and the defective surface were considered for comparative analysis. The key finding is the enhanced desorption of SiCl molecules, facilitated by point defects under laser pulse irradiation. The presence of point defects was found to effectively reduce both the desorption energy barrier and the laser intensity threshold required for desorption. Additionally, extra defective levels within the band gap were observed through the density-of-state diagram. Based on these findings, a defect-mediated etching regime was proposed to elucidate the layer-by-layer etching process. This study provides atomistic insight into understanding the role of defects in laser-induced ALE processes. The presence of point defects can enhance the etching selectivity between the topmost layer and the underlying layers, thereby contributing to highly efficient and damage-free etching processes through the defect-mediated etching mechanism.

Non-Linear Optics at Twist Interfaces in h-BN/SiC Heterostructures

Understanding the emergent electronic structure in twisted atomically thin layers has led to the exciting field of twistronics. However, practical applications of such systems are challenging since the specific angular correlations between the layers must be precisely controlled and the layers have to be single crystalline with uniform atomic ordering. Here, an alternative, simple, and scalable approach is suggested, where nanocrystallinetwo-dimensional (2D) film on 3D substrates yields twisted-interface-dependent properties. Ultrawide-bandgap hexagonal boron nitride (h-BN) thin films are directly grown on high in-plane lattice mismatched wide-bandgap silicon carbide (4H-SiC) substrates to explore the twist-dependent structure-property correlations. Concurrently, nanocrystalline h-BN thin film shows strong non-linear second-harmonic generation and ultra-low cross-plane thermal conductivity at room temperature, which are attributed to the twisted domain edges between van der Waals stacked nanocrystals with random in-plane orientations. First-principles calculations based on time-dependent density functional theory manifest strong even-order optical nonlinearity in twisted h-BN layers. This work unveils that directly deposited 2D nanocrystalline thin film on 3D substrates could provide easily accessible twist-interfaces, therefore enabling a simple and scalable approach to utilize the 2D-twistronics integrated in 3D material devices for next-generation nanotechnology.

Circularly polarized attosecond light generation from OCS molecules irradiated by the combination of linear polarized infrared and orthogonal terahertz fields

Researching ultrafast dynamics and creating coherent light sources will both benefit significantly from the establishment of polarization control in high-order harmonic generation (HHG). By employing the time-dependent density functional theory method, we investigate HHG of carbonyl sulfide molecules using a combination of a linear polarized infrared (IR) laser and a weaker orthogonal Terahertz (THz) field. Our findings show that by adjusting the amplitude of the THz field, the movement scale of electrons in the THz direction can be tuned, thereby one can control the harmonic intensity in the IR laser direction. This method allows for the creation of near-circularly polarized attosecond pulses. Furthermore, the ellipticity of the attosecond pulse may be changed by modifying the carrier-envelope phase of the IR laser pulse.

Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction

We revisit the use of Meta-Generalized Gradient Approximations (mGGAs) in time-dependent density functional theory, reviewing conceptual questions and solving the generalized Kohn-Sham equations by real-time propagation. After discussing the technical aspects of using mGGAs in combination with pseudopotentials and comparing real-space and basis set results, we focus on investigating the importance of the current-density based gauge invariance correction. For the two modern mGGAs that we investigate in this work, TASK and r2SCAN, we observe that for some systems, the current density correction leads to negligible changes, but for others, it changes excitation energies by up to 40% and more than 0.8 eV. In the cases that we study, the agreement with the reference data is improved by the current density correction.

Ultrafast Spin Dynamics and Photoinduced Insulator-to-Metal Transition in α-RuCl3

Laser-induced ultrafast demagnetization is a phenomenon of utmost interest and attracts significant attention because it enables potential applications in ultrafast optoelectronics and spintronics. As a spin-orbit coupling assisted magnetic insulator, α-RuCl3 provides an attractive platform to explore the physics of electronic correlations and unconventional magnetism. Using time-dependent density functional theory, we explore the ultrafast laser-induced dynamics of the electronic and magnetic structures in α-RuCl3. Our study unveils that laser pulses can introduce ultrafast demagnetizations, accompanied by an out-of-equilibrium insulator-to-metal transition in a few tens of femtoseconds. The spin response significantly depends on the laser wavelength and polarization on account of the electron correlations, band renormalizations, and charge redistributions. These findings provide physical insights into the coupling between the electronic and magnetic degrees of freedom in α-RuCl3 and shed light on suppressing the long-range magnetic orders and reaching a proximate spin liquid phase for two-dimensional magnets on an ultrafast time scale.

Minimum structure of high-order harmonic spectrum from molecular multi-orbital effects involving inner-shell orbitals

The spectral features of high-order harmonic spectra can provide rich information for probing the structure and dynamics of molecules in intense laser fields. We theoretically study the high harmonic spectrum with the laser polarization direction perpendicular to the N2O molecule and find a minimum structure in the plateau region of the harmonic spectrum. Through analyzing the time-dependent survival probability of different electronic orbitals and the time-dependent wave packet evolution, it is found that this minimum position is caused by the harmonic interference of HOMO a, HOMO-1, and HOMO-3 a orbitals. Moreover, this interference minimum is discovered over a wide frequency range of 0.087 a.u. to 0.093 a.u., as well as a range of driving laser intensities with peak amplitudes between 0.056 a.u. and 0.059 a.u.. This study sheds light on the multi-electron effects and ultrafast dynamics of inner-shell electrons in intense laser pulses, which are crucial for understanding and controlling chemical reactions in molecules.

Band Nonlinearity-Enabled Manipulation of Dirac Nodes, Weyl Cones, and Valleytronics with Intense Linearly Polarized Light

We study low-frequency linearly polarized laser-dressing in materials with valley (graphene and hexagonal-Boron-Nitride) and topological (Dirac- and Weyl-semimetals) properties. In Dirac-like linearly dispersing bands, the laser substantially moves the Dirac nodes away from their original position, and the movement direction can be fully controlled by rotating the laser polarization. We prove that this effect originates from band nonlinearities away from the Dirac nodes. We further demonstrate that this physical mechanism is widely applicable and can move the positions of the valley minima in hexagonal materials to tune valley selectivity, split and move Weyl cones in higher-order Weyl semimetals, and merge Dirac nodes in three-dimensional Dirac semimetals. The model results are validated with ab initio calculations. Our results directly affect efforts for exploring light-dressed electronic structure, suggesting that one can benefit from band nonlinearity for tailoring material properties, and highlight the importance of the full band structure in nonlinear optical phenomena in solids.

Ultrafast control over chiral sum-frequency generation

We introduce an ultrafast all-optical approach for efficient chiral recognition that relies on the interference between two low-order nonlinear processes that are ubiquitous in nonlinear optics: sum-frequency generation and third-harmonic generation. In contrast to traditional sum-frequency generation, our approach encodes the medium's handedness in the intensity of the emitted harmonic signal, rather than in its phase, and it enables full control over the enantiosensitive response. We show how, by sculpting the sub-optical-cycle oscillations of the driving laser field, we can force one molecular enantiomer to emit bright light while its mirror twin remains dark, thus reaching the ultimate efficiency limit of chiral sensitivity via low-order nonlinear light-matter interactions. Our work paves the way for ultrafast and highly efficient imaging and control of the chiral electronic clouds of chiral molecules using lasers with moderate intensities, in all states of matter: from gases to liquids to solids, with molecular specificity and on ultrafast time scales.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Ab initio simulations of ultrashort laser pulse interaction with Cl-Si(100): implications for atomic layer etching

Due to the remarkable resistance of SiCl against photo-induced desorption, achieving atomic layer etching (ALE) of Cl-Si(100) through a laser-based method has remained a formidable challenge. In this study, we investigate the interaction between ultrashort laser pulses and the Cl-Si(100) surface via ab initio simulations that combine real-time time-dependent density functional theory and molecular dynamics. Our results demonstrate the direct desorption of the stubborn SiCl layer through the application of appropriate femtosecond laser pulses. Notably, the desorption process is enhanced by pulses with higher laser intensity, shorter wavelength, and longer pulse duration. There is a threshold intensity beyond which the SiCl can be directly desorbed under laser pulses with a wavelength of 488 nm and a pulse duration of 40 ℏ eV-1 (26.3 fs). Analysis of electron localization function reveals a critical bond breaking length of 2.98 Å between Si-Si, connecting SiCl to the bulk material. The time evolution of bond lengths and forces reveals that the desorption of SiCl is primarily driven by repulsive forces generated within the Si-Si bond. Furthermore, electron density difference analysis and Keldysh factor calculations indicate that these repulsive forces arise from multiphoton ionization. This study provides crucial atomic-level insights into the interaction of ultrashort laser pulses with Cl-Si(100), thereby propelling the advancement of laser-induced atomic layer etching techniques.

First-Principles Simulations of Tip Enhanced Raman Scattering Reveal Active Role of Substrate on High-Resolution Images

Tip-enhanced Raman scattering (TERS) has emerged as a powerful tool to obtain subnanometer spatial resolution fingerprints of atomic motion. Theoretical calculations that can simulate the Raman scattering process and provide an unambiguous interpretation of TERS images often rely on crude approximations of the local electric field. In this work, we present a novel and first-principles-based method to compute TERS images by combining Time Dependent Density Functional Theory (TD-DFT) and Density Functional Perturbation Theory (DFPT) to calculate Raman cross sections with realistic local fields. We present TERS results on free-standing benzene and C60 molecules, and on the TCNE molecule adsorbed on Ag(100). We demonstrate that chemical effects on chemisorbed molecules, often ignored in TERS simulations of larger systems, dramatically change the TERS images. This observation calls for the inclusion of chemical effects for predictive theory-experiment comparisons and an understanding of molecular motion at the nanoscale.

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn-Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix-vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh-Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Plasmon-Mediated Hydrogen Dissociation with Symmetry Tunability

The atomic-scale mechanism of plasmon-mediated H₂ dissociation on gold nanoclusters is investigated using time-dependent density functional theory. The position relationship between the nanocluster and H₂ has a strong influence on the reaction rate. When the hydrogen molecule is located in the interstitial center of the plasmonic dimer, the hot spot here has a great field enhancement, which can promote dissociation effectively. The change in the molecular position results in symmetry breaking, and the molecular dissociation is inhibited. For the asymmetric structure, direct charge transfer from the gold cluster to the antibonding state of the hydrogen molecule by plasmon decay makes a prominent contribution to the reaction. The results provide deep insights into the influence of structural symmetry on plasmon-assisted photocatalysis in the quantum regime.

Ab initio real-time quantum dynamics of charge carriers in momentum space

Application of the non-adiabatic molecular dynamics (NAMD) approach is limited to studying carrier dynamics in the momentum space, as a supercell is required to sample the phonon excitation and electron-phonon (e-ph) interaction at different momenta in a molecular dynamics simulation. Here we develop an ab initio approach for the real-time charge carrier quantum dynamics in the momentum space (NAMD_k) by directly introducing e-ph coupling into the Hamiltonian based on the harmonic approximation. The NAMD_k approach maintains the zero-point energy and includes memory effects of carrier dynamics. The application of NAMD_k to the hot carrier dynamics in graphene reveals the phonon-specific relaxation mechanism. An energy threshold of 0.2 eV-defined by two optical phonon modes-separates the hot electron relaxation into fast and slow regions with lifetimes of pico- and nanoseconds, respectively. The NAMD_k approach provides an effective tool to understand real-time carrier dynamics in the momentum space for different materials.

OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources

We have proposed, for the first time, an OpenCL implementation for the all-electron density-functional perturbation theory (DFPT) calculations in FHI-aims, which can effectively compute all its time-consuming simulation stages, i.e., the real-space integration of the response density, the Poisson solver for the calculation of the electrostatic potential, and the response Hamiltonian matrix, by utilizing various heterogeneous accelerators. Furthermore, to fully exploit the massively parallel computing capabilities, we have performed a series of general-purpose graphics processing unit (GPGPU)-targeted optimizations that significantly improved the execution efficiency by reducing register requirements, branch divergence, and memory transactions. Evaluations on the Sugon supercomputer have shown that notable speedups can be achieved across various materials.

Generation of the isolated highly elliptically polarized attosecond pulse using the polarization gating technique: TDDFT approach

This paper theoretically investigates the generation of isolated elliptically polarized attosecond pulses with a tunable ellipticity from the interaction of Cl₂ molecule and a polarization-gating laser pulse. A three-dimensional calculation based on the time-dependent density functional theory is done. Two different methods are proposed for generating elliptically polarized single attosecond pulses. The first method is based on applying a single-color polarization gating laser and controlling the orientation angle of the Cl₂ molecule with respect to the polarization direction of the laser at the gate window. An attosecond pulse with an ellipticity of 0.66 and a pulse duration of 275 as is achieved by tuning the molecule orientation angle to 40° in this method and superposing harmonics around the harmonic cutoff. The second method is based on irradiating an aligned Cl₂ molecule with a two-color polarization gating laser. The ellipticity of the attosecond pulses obtained by this method can be controlled by adjusting the intensity ratio of the two colors. Employing an optimized intensity ratio and superposing harmonics around the harmonic cutoff would lead to the generation of an isolated, highly elliptically polarized attosecond pulse with an ellipticity of 0.92 and a pulse duration of 648 as.

First-principles study of electronic and optical properties in 1-dimensional oligomeric derivatives of telomestatin

Real-space self-interaction corrected (time-dependent) density functional theory has been used to investigate the ground-state electronic structure and optical absorption profiles of a series of linear oligomers inspired by the natural product telomestatin. Length-dependent development of plasmonic excitations in the UV region is seen in the neutral species which is augmented by polaron-type absorption with tunable wavelengths in the IR when the chains are doped with an additional electron/hole. Combined with a lack of absorption in the visible region this suggests these oligomers as good candidates for applications such as transparent antennae in dye-sensitised solar energy collection materials. Due to strong longitudinal polarisation in their absorption spectra, these compounds are also indicated for use in nano-structured devices displaying orientation-sensitive optical responses.

Gaussian-Approximated Poisson Preconditioner for Iterative Diagonalization in Real-Space Density Functional Theory

Various real-space methods optimized on massive parallel computers have been developed for efficient large-scale density functional theory (DFT) calculations of materials and biomolecules. The iterative diagonalization of the Hamiltonian matrix is a computational bottleneck in real-space DFT calculations. Despite the development of various iterative eigensolvers, the absence of efficient real-space preconditioners has hindered their overall efficiency. An efficient preconditioner must satisfy two conditions: appropriate acceleration of the convergence of the iterative process and inexpensive computation. This study proposed a Gaussian-approximated Poisson preconditioner (GAPP) that satisfied both conditions and was suitable for real-space methods. A low computational cost was realized through the Gaussian approximation of a Poisson Green's function. Fast convergence was achieved through the proper determination of Gaussian coefficients to fit the Coulomb energies. The performance of GAPP was evaluated for several molecular and extended systems, and it showed the highest efficiency among the existing preconditioners adopted in real-space codes.

Energy Transfer to Molecular Adsorbates by Transient Hot Electron Spillover

Hot electron (HE) photocatalysis is one of the most intriguing fields of nanoscience, with a clear potential for technological impact. Despite much effort, the mechanisms of HE photocatalysis are not fully understood. Here we investigate a mechanism based on transient electron spillover on a molecule and subsequent energy release into vibrational modes. We use state-of-the-art real-time Time Dependent Density Functional Theory (rt-TDDFT), simulating the dynamics of a HE moving within linear chains of Ag or Au atoms, on which CO, N₂, or H₂O are adsorbed. We estimate the energy a HE can release into adsorbate vibrational modes and show that certain modes are selectively activated. The energy transfer strongly depends on the adsorbate, the metal, and the HE energy. Considering a cumulative effect from multiple HEs, we estimate this mechanism can transfer tenths of an eV to molecular vibrations and could play an important role in HE photocatalysis.

Time-Dependent Density-Functional Theory for Determining the Electron-Capture Cross Section for Protons Impacting on Atoms and Molecules

The use of the Time-Dependent Density-Functional Theory (TDDFT) has increased in the atomic collision field. Calculating the electron-capture cross section (ECCS) for protons is an important question in hadrontherapy and plasma physics, among other areas. In previous studies, it was shown that the approach based on the Local Density Approximation (LDA) fails in the 1-50 keV region, requiring the use of the Optimized Effective Potential (OEP) method. In this work, the ECCS values for 1-50 keV protons impacting on isolated hydrogen, carbon, nitrogen, oxygen, and nitrogenous atoms were determined using the TDDFT. It is shown that adding the Self Interaction Correction to the LDA (LDA-Sic) allows obtaining results close to those provided by the OEP and experiments, with the advantage that the LDA-Sic consumes less computational time. In addition, it was demonstrated that it is imperative to include the spin correction for the specific helium and oxygen cases, in order to get good results for the ECCS using the TDDFT. Theoretical results obtained in this work show excellent agreement with experimental values.

Real-time dynamic simulation of laser-induced N2 dissociation on two-dimensional graphene sheets

Due to its relatively high inertness, nitrogen dissociation at ambient temperature and pressure has always been a challenging task. Plasmon driven photocatalysis has proved to be an effective method. Owing to their unique physical, chemical, and electronic properties, two-dimensional planar materials have become the most promising candidates to replace noble metal catalytic nitrogen reduction. In this study, real-time dynamics of N₂ dissociation on graphene sheets under femtosecond laser irradiation was studied by using time-dependent density functional theory. We confirm that electrons generated by plasmon excitation of graphene transfer to the N₂ molecular antibonding orbital and activate the N-N bond. The threshold of laser intensity of N₂ dissociation can be effectively reduced by mixing CO molecules. This work provides basic insights for understanding the plasmon induced N₂ activation process at the atomic scale and proves that graphene can be used as one of the candidate materials for N₂ reduction photocatalysts with excellent performance.

Cavity-Free Continuum Solvation: Implementation and Parametrization in a Multiwavelet Framework

We present a multiwavelet-based implementation of a quantum/classical polarizable continuum model. The solvent model uses a diffuse solute-solvent boundary and a position-dependent permittivity, lifting the sharp-boundary assumption underlying many existing continuum solvation models. We are able to include both surface and volume polarization effects in the quantum/classical coupling, with guaranteed precision, due to the adaptive refinement strategies of our multiwavelet implementation. The model can account for complex solvent environments and does not need a posteriori corrections for volume polarization effects. We validate our results against a sharp-boundary continuum model and find a very good correlation of the polarization energies computed for the Minnesota solvation database.

Dynamic Precision Approach for Accelerating Large-Scale Eigenvalue Solvers in Electronic Structure Calculations on Graphics Processing Units

Single precision (SP) arithmetic can be greatly accelerated as compared to double precision (DP) arithmetic on graphics processing units (GPUs). However, the use of SP in the whole process of electronic structure calculations is inappropriate for the required accuracy. We propose a 3-fold dynamic precision approach for accelerated calculations but still with the accuracy of DP. Here, SP, DP, and mixed precision are dynamically switched during an iterative diagonalization process. We applied this approach to the locally optimal block preconditioned conjugate gradient method to accelerate a large-scale eigenvalue solver for the Kohn-Sham equation. We determined a proper threshold for switching each precision scheme by examining the convergence pattern on the eigenvalue solver only with the kinetic energy operator of the Kohn-Sham Hamiltonian. As a result, we achieved up to 8.53× and 6.60× speedups for band structure and self-consistent field calculations, respectively, for test systems under various boundary conditions on NVIDIA GPUs.

The contribution of inner electron excitation to the electronic stopping power of palladium for protons

The electronic stopping power of palladium (Pd) for protons is investigated based on time-dependent density functional theory combined with Ehrenfest molecular dynamics simulations. The electronic stopping power of Pd with explicitly considering inner electrons for protons is calculated and the excitation mechanism for the inner electrons of Pd is revealed. The velocity proportionality of the low-energy stopping power of Pd is reproduced. Our study verified that the inner electron excitation contributes significantly to the electronic stopping power of Pd in the high energy range, which is strongly dependent on the impact parameter. The electronic stopping power obtained from the off-channeling geometry is in quantitative agreement with the experimental data in a wide velocity range, and the discrepancy around the stopping maximum is further reduced by considering the relativistic correction on the binding energy of inner electrons. The velocity dependence of the mean steady-state charge of protons is quantified, and the results showed that the participation of 4p-electrons reduces the mean steady-state charge of protons, and consequently decreases the electronic stopping power of Pd in the low energy range.

Solid state-like high harmonic generation from cluster molecules with rotational periodicities

High harmonic generation (HHG) from solid-state crystals in strong laser fields has been understood by the band structure of the solids, which is based on the periodic boundary condition (PBC) due to translational invariance. For the systems with PBC due to rotational invariance, an analogous Bloch theorem can be applied. Considering a ring-type cluster of cyclo[18]carbon as an example, we develop a quasi-band model and predict the solid state-like HHG in this system. Under the irradiation of linearly polarized laser field, cyclo[18]carbon exhibits solid state-like HHG originated from intraband oscillations and interband transitions, which, in turn, is promising to optically detect the symmetry and geometry of molecular or material structures. Our results based on the Liouville-von Neumann equations are well reproduced by the time-dependent density functional theory calculations and are foundational in providing a connection linking the HHG physics of gases and solids.

Fast Exchange with Gaussian Basis Set Using Robust Pseudospectral Method

In this article, we present an algorithm to efficiently evaluate the exchange matrix in periodic systems when a Gaussian basis set with pseudopotentials is used. The usual algorithm for evaluating exchange matrix scales cubically with the system size because one has to perform O(N²) fast Fourier transform (FFT). Here, we introduce an algorithm that retains the cubic scaling but reduces the prefactor significantly by eliminating the need to do FFTs during each exchange build. This is accomplished by representing the products of Gaussian basis function using a linear combination of an auxiliary basis the number of which scales linearly with the size of the system. We store the potential due to these auxiliary functions in memory, which allows us to obtain the exchange matrix without the need to do FFT, albeit at the cost of additional memory requirement. Although the basic idea of using auxiliary functions is not new, our algorithm is cheaper due to a combination of three ingredients: (a) we use a robust pseudospectral method that allows us to use a relatively small number of auxiliary basis to obtain high accuracy; (b) we use occ-RI exchange, which eliminates the need to construct the full exchange matrix; and (c) we use the (interpolative separable density fitting) ISDF algorithm to construct these auxiliary basis sets that are used in the robust pseudospectral method. The resulting algorithm is accurate, and we note that the error in the final energy decreases exponentially rapidly with the number of auxiliary functions.

Analyzing Excitation-Energy Transfer Based on the Time-Dependent Density Functional Theory in Real Time

Excitation-energy transfer is a key step in processes such as photosynthesis that convert light into other forms of energy. Time-dependent density functional theory (DFT) in real time is ideal for the first-principles simulation of such processes due to its computational efficiency. We here demonstrate how real-time DFT can be used for analyzing excitation-energy transfer from first-principles. We discuss several measures of energy transfer that are based solely on the time-dependent density, are well founded in the DFT framework, allow for intuitive understanding and visualization, and reproduce important limiting cases of an analytical model. We demonstrate their usefulness in calculations for model systems, both with static nuclei and in the context of DFT-based Ehrenfest dynamics.

Simple Exchange-Correlation Energy Functionals for Strongly Coupled Light-Matter Systems Based on the Fluctuation-Dissipation Theorem

Recent experimental advances in strongly coupled light-matter systems have sparked the development of general ab initio methods capable of describing interacting light-matter systems from first principles. One of these methods, quantum-electrodynamical density-functional theory (QEDFT), promises computationally efficient calculations for large correlated light-matter systems with the quality of the calculation depending on the underlying approximation for the exchange-correlation functional. So far no true density-functional approximation has been introduced limiting the efficient application of the theory. In this Letter, we introduce the first gradient-based density functional for the QEDFT exchange-correlation energy derived from the adiabatic-connection fluctuation-dissipation theorem. We benchmark this simple-to-implement approximation on small systems in optical cavities and demonstrate its relatively low computational costs for fullerene molecules up to C_{180} coupled to 400 000 photon modes in a dissipative optical cavity. This Letter now makes first principle calculations of much larger systems possible within the QEDFT framework effectively combining quantum optics with large-scale electronic structure theory.

Neural network-based pseudopotential: development of a transferable local pseudopotential

Transferable local pseudopotentials (LPPs) are essential for fast quantum simulations of materials. However, various types of LPPs suffer from low transferability, especially since they do not consider the norm-conserving condition. Here we propose a novel approach based on a deep neural network to produce transferable LPPs. We introduced a generalized Kerker method expressed with the deep neural network to represent the norm-conserving pseudo-wavefunctions. Its unique feature is that all necessary conditions of pseudopotentials can be explicitly considered in terms of a loss function. Then, it can be minimized using the back-propagation technique just with single point all-electron atom data. To assess the transferability and accuracy of the neural network-based LPPs (NNLPs), we carried out density functional theory calculations for the s- and p-block elements of the second to the fourth periods. The NNLPs outperformed other types of LPPs in both atomic and bulk calculations for most elements. In particular, they showed good transferability by predicting various properties of bulk systems including binary alloys with higher accuracy than LPPs tailored to bulk data.

Application of two-component neural network for exchange-correlation functional interpolation

Density functional theory (DFT) is one of the primary approaches to solving the many-body Schrodinger equation. The essential part of the DFT theory is the exchange-correlation (XC) functional, which can not be obtained in analytical form. Accordingly, the accuracy improvement of the DFT is mainly based on the development of XC functional approximations. Commonly, they are built upon analytic solutions in low- and high-density limits and result from quantum Monte Carlo or post-Hartree-Fock numerical calculations. However, there is no universal functional form to incorporate these data into XC functional. Instead, various parameterizations use heuristic rules to build a specific XC functional. The neural network (NN) approach to interpolate the data from higher precision theories can give a unified path to parametrize an XC functional. Moreover, data from many existing quantum chemical databases could provide the XC functional with improved accuracy. We develop NN XC functional, which gives exchange potential and energy density without direct derivatives of exchange-correlation energy density. Proposed NN architecture consists of two parts NN-E and NN-V, which could be trained in separate ways, adding new flexibility to XC functional. We also show that the developed NN XC functional converges in the self-consistent cycle and gives reasonable energies when applied to atoms, molecules, and crystals.

Probing phonon dynamics with multidimensional high harmonic carrier-envelope-phase spectroscopy

High harmonic generation (HHG) has recently been established as a powerful method for probing ultrafast electron dynamics in solids. However, it remains unknown if HHG can be similarly applied for probing lattice distortions such as phonons. Specifically, it is unclear if the extreme nonlinearity of HHG can contribute to enhanced temporal resolution or sensitivity for probing lattice dynamics (compared to other, perturbative, methods). Here, we theoretically explore HHG in solids with active phonons. We present a pump-probe and multidimensional spectroscopy approach that relies on carrier-envelope-phase-sensitivity, in which HHG is highly sensitive for phonon dynamics. Strikingly, the predicted temporal resolution is ∼1 femtosecond, much below the probe pulse duration, owing to the subcycle nature of the approach.

We explore pump-probe high harmonic generation (HHG) from monolayer hexagonal-boron-nitride, where a terahertz pump excites coherent optical phonons that are subsequently probed by an intense infrared pulse that drives HHG. We find, through state-of-the-art ab initio calculations, that the structure of the emission spectrum is attenuated by the presence of coherent phonons and no longer comprises discrete harmonic orders, but rather a continuous emission in the plateau region. The HHG yield strongly oscillates as a function of the pump-probe delay, corresponding to ultrafast changes in the lattice such as specific bond compression or stretching dynamics. We further show that in the regime where the excited phonon period and the pulse duration are of the same order of magnitude, the HHG process becomes sensitive to the carrier-envelope phase (CEP) of the driving field, even though the pulse duration is so long that no such sensitivity is observed in the absence of coherent phonons. The degree of CEP sensitivity versus pump-probe delay is shown to be a highly selective measure for instantaneous structural changes in the lattice, providing an approach for ultrafast multidimensional HHG spectroscopy. Remarkably, the obtained temporal resolution for phonon dynamics is ∼1 femtosecond, which is much shorter than the probe pulse duration because of the inherent subcycle contrast mechanism. Our work paves the way toward routes of probing phonons and ultrafast material structural changes with subcycle temporal resolution and provides a mechanism for controlling the HHG spectrum.

Ab Initio Cluster Approach for High Harmonic Generation in Liquids

High harmonic generation (HHG) takes place in all phases of matter. In gaseous atomic and molecular media, it has been extensively studied and is very well understood. In solids, research is ongoing, but a consensus is forming for the dominant microscopic HHG mechanisms. In liquids, on the other hand, no established theory yet exists, and approaches developed for gases and solids are generally inapplicable, hindering our current understanding. We develop here a powerful and reliable ab initio cluster-based approach for describing the nonlinear interactions between isotropic bulk liquids and intense laser pulses. The scheme is based on time-dependent density functional theory and utilizes several approximations that make it feasible yet accurate in realistic systems. We demonstrate our approach with HHG calculations in water, ammonia, and methane liquids and compare the characteristic response of polar and nonpolar liquids. We identify unique features in the HHG spectra of liquid methane that could be utilized for ultrafast spectroscopy of its chemical and physical properties, including a structural minimum at 15-17 eV that is associated solely with the liquid phase. Our results pave the way to accessible calculations of HHG in liquids and illustrate the unique nonlinear nature of liquid systems.

New opportunities for ultrafast and highly enantio-sensitive imaging of chiral nuclear dynamics enabled by synthetic chiral light

Synthetic chiral light [D. Ayuso et al., Nat. Photon., 2019, 13, 866-871] has opened up new opportunities for ultrafast and highly efficient imaging and control of chiral matter. Here we show that the giant enantio-sensitivity enabled by such light could be exploited to probe chiral nuclear rearrangements during chemical reactions in a highly enantio-sensitive manner. Using a state-of-the-art implementation of real-time time-dependent density functional theory, we explore how the nonlinear response of the prototypical chiral molecule H₂O₂ changes as a function of its dihedral angle, which defines its handedness. The macroscopic intensity emitted from randomly oriented molecules at even harmonic frequencies (of the fundamental) depends strongly on this nuclear coordinate. Because of the ultrafast nature of such nonlinear interactions, the direct mapping between the dissymmetry factor and the nuclear geometry provides a way to probe chiral nuclear dynamics at their natural time scales. Our work paves the way for ultrafast and highly efficient imaging of enantio-sensitive dynamics in more complex chiral systems, including biologically relevant molecules.

Nanovectorization of Ivermectin to avoid overdose of drugs

Ivermectin is an antiparasitic drug that results in the death of the targeted parasites using several mechanical actions. While very well supported, it can induce in rare cases, adverse effects including coma and respiratory failure in case of overdose. This problem should be solved especially in an emergency situation. For instance, the first pandemic of the 21th century was officially declared in early 2020, and while several vaccines around the worlds have been used, an effective treatment against this new strain of coronavirus, better known as SARS-CoV-2, should also be considered, especially given the massive appearance of variants. From all the tested therapies, Ivermectin showed a potential reduction of the viral portability, but sparked significant debate around the dose needed to achieve these positive results. To answer this general question, we propose, using simulations, to show that the nanovectorization of Ivermectin on BN oxide nanosheets can increase the transfer of the drug to its target and thus decrease the quantity of drug necessary to cope with the disease. This first application could help science to develop such nanocargo to avoid adverse effects.Communicated by Ramaswamy H. Sarma.

Machine learning the derivative discontinuity of density-functional theory

Machine learning is a powerful tool to design accurate, highly non-local, exchange-correlation functionals for density functional theory. So far, most of those machine learned functionals are trained for systems with an integer number of particles. As such, they are unable to reproduce some crucial and fundamental aspects, such as the explicit dependency of the functionals on the particle number or the infamous derivative discontinuity at integer particle numbers. Here we propose a solution to these problems by training a neural network as the universal functional of density-functional theory that (a) depends explicitly on the number of particles with a piece-wise linearity between the integer numbers and (b) reproduces the derivative discontinuity of the exchange-correlation energy. This is achieved by using an ensemble formalism, a training set containing fractional densities, and an explicitly discontinuous formulation.

Detecting multiple chiral centers in chiral molecules with high harmonic generation

Characterizing chiral is highly important for applications in the pharmaceutical industry, as well as in the study of dynamical chemical and biological systems. However, this task has remained challenging, especially due to the ongoing increasing complexity and size of the molecular structure of drugs and active compounds. In particular, large molecules with many active chiral centers are today ubiquitous, but remain difficult to structurally analyze due to their high number of stereoisomers. Here we theoretically explore the sensitivity of high harmonic generation (HHG) to the chiral of molecules with a varying number of active chiral centers. We find that HHG driven by bi-chromatic non-collinear lasers is a sensitive probe for the stereo-configuration of a chiral molecule. We first show through calculations (from benchmark chiral molecules with up to three chiral centers) that the HHG spectrum is imprinted with information about the handedness of each chiral center in the driven molecule. Next, we show that using both classical- and deep-learning-based reconstruction algorithms, the composition of an unknown mixture of stereoisomers can be reconstructed with high fidelity by a single-shot HHG measurement. Our work illustrates how the combination of non-linear optics and machine learning might open routes for ultra-sensitive sensing in chiral systems.

Indirect to Direct Charge Transfer Transition in Plasmon-Enabled CO2 Photoreduction

Understanding hot carrier dynamics between plasmonic nanomaterials and its adsorbate is of great importance for plasmon-enhanced photoelectronic processes such as photocatalysis, optical sensing and spectroscopic analysis. However, it is often challenging to identify specific dominant mechanisms for a given process because of the complex pathways and ultrafast interactive dynamics of the photoelectrons. Here, using CO₂ reduction as an example, the underlying mechanisms of plasmon-driven catalysis at the single-molecule level using time-dependent density functional theory calculations is clearly probed. The CO₂ molecule adsorbed on two typical nanoclusters, Ag₂₀ and Ag₁₄₇ , is photoreduced by optically excited plasmon, accompanied by the excitation of asymmetric stretching and bending modes of CO₂ . A nonlinear relationship has been identified between laser intensity and reaction rate, demonstrating a synergic interplay and transition from indirect hot-electron transfer to direct charge transfer, enacted by strong localized surface plasmons. These findings offer new insights for CO₂ photoreduction and for the design of effective pathways toward highly efficient plasmon-mediated photocatalysis.