Time-dependent Density Functional Theory Simulations Research Articles

Electronically induced defect creation at semiconductor/oxide interface revealed by time-dependent density functional theory

Carrier induced defect creation at the semiconductor-oxide interface has been known as the origin of electronic device degradation for a long time, but how exactly the interface lattice can be damaged by carriers (especially low-energy ones) remains unclear. Here we carry out real-time time-dependent density functional theory simulations on concrete $\mathrm{Si}/\mathrm{Si}{\mathrm{O}}_{2}$ interfaces to study the interaction between excited electrons and interface bonds. We show that the normal interface Si-H bonds are generally resistant to electrons due to the delocalized nature and high energy level of the Si-H antibonding states, and due to the high-energy barrier to break the Si-H bond. However, if an additional hydrogen atom exists by attaching to a nearby oxygen atom (forming a ``Si-H\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}H-O'' complex), the Si-H bond will be greatly weakened, including the reduction of energy barrier for bond breaking, and the lowering of the antibonding state energy level which favors electron injection. Together with the multiple vibrational excitation process, the corresponding Si-H bond can be broken much more easily. Thus we propose that the Si-H\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}H-O complex will be the center for defect creation and device degradation. We also explain why such a center might be relatively easy to form during the hydrogen annealing process.

Efficient enantioresolution of aromatic α-hydroxy acids with Cinchona alkaloid-based zwitterionic stationary phases and volatile polar-ionic eluents

Single enantiomers of mandelic acid (1), 3-phenyllactic acid (2), and 3-(4-hydroxyphenyl)lactic acid (3) are the subject of many fields of investigation, spanning from the pharmaceutical synthesis to that of biocompatible and biodegradable polymers, while passing from the interest towards their antimicrobial activity to their role as biomarkers of particular pathological conditions or occupational exposures to specific xenobiotics. All above mentioned issues justify the need for accurate analytical methods enabling the correct determination of the individual enantiomers. So far, all the developed liquid chromatography (LC) methods were not or hardly compatible with mass spectrometry (MS) detection. In this paper, a commercially available Cinchona-alkaloid derivative zwitterionic chiral stationary phase [that is, the CHIRALPAK® ZWIX(−)] was successfully used to optimize the enantioresolution of compounds 1–3 under polar-ionic (PI) conditions with a mobile phase consisting of an acetonitrile/methanol 95/5 (v/v) mixture with 80 mM formic acid. With the optimized conditions, enantioseparation and enantioresolution values up to 1.46 and 4.41, respectively, were obtained. In order to assess the applicability of the optimized enantioselective chromatography conditions in real-life scenarios and on MS-based systems, a proof-of-concept application was efficiently carried out by analysing dry urine spot samples spiked with 1 by means of a LC-MS system. The (S)<(R) enantiomer elution order (EEO) was established for compounds 1 and 2 by analysing a pure enantiomeric standard of known configuration. This was not possible for 3 because not commercially available. For this compound, the same EEO was identified applying a procedure based on ab initio time-dependent density-functional theory simulations coupled to electronic circular dichroism analyses. Moreover, a molecular dynamics simulation unveiled the role of the phenolic OH in compound 3 in the retention mechanism.

Synthesis, characterization, surface analysis, optical activity and solvent effects on the electronic absorptions of Schiff base-functionalized amino thiophene derivatives: Experimental and TD-DFT investigations

New derivatives of 3,3′-di-(((E)-(5-substituted-thiophen-2-yl)methylene)amino)-N-methylpropylamines (3a,b) were synthesized from the dehydration of N-methyl-diaminopropylamine and 2-thielylcarboxaldehyde derivatives; spectroscopically characterized by 1H- and 13C-NMR, IR, LC-MS, UV-Vis and elemental analysis. The nature of the electronic transitions of the SB compounds was investigated using Time-Dependent Density-Functional Theory (TD-DFT). Surface analysis and influence of solvent polarity on spectral properties were examined and established consequently. The molecular electrostatic potential (MEP) revealed that the two geometrical structures were found to be quite similar in term of electronic distributions. The presence of different electrophilic and nucleophilic sites located on the surfaces was suggested to stabilize the structures via classical H-bond and non-classical C-H…π interactions, in addition to interact with assorted solvent molecules. On the other hand, the solvatochromism of compound 3b revealed a gradual shift to the red region through the increase of the solvent polarity, recording a 12 nm of bathochromic shift. The solvation relationship between the experimental λmax and Gutmann's donicity numbers displayed a sense of positive linear behavior with a fluctuation, which was ascribed to a week interaction between molecules of solute and selected solvent. The band gap energy of compound 3a was evaluated experimentally and computationally. Using optical absorption spectra, a value of ‒ 3.801 eV was estimated following Tauc approach, while ‒ 3.720 eV was resulting from TD-DFT simulation.

Intersystem Crossing in Boron-Based Donor-Spiro-Acceptor Organic Chromophore: A Detailed Theoretical Study.

Intersystem crossing and reverse intersystem crossing (rISC) processes were investigated in a boron-based donor-spiro-acceptor organic chromophore which shows thermally activated delayed fluorescence. Due to the perpendicular arrangement between donor and acceptor moieties, the HOMO and the LUMO are spatially separated, and the compound shows charge transfer (CT) transitions. We found both S1 and T1 excited states are CT in nature (i.e., electron and hole wave functions are localized on acceptor and donor units, respectively) and T2, which is higher in energy than S1 and T1, is locally excited in nature (i.e., both electron and hole wave functions are localized on an acceptor unit). Because of the same nature of excitation (i.e., CT here), the spin-orbit coupling matrix element between S1 and T1 is very low and insignificant exciton conversion occurs from the T1 state to the S1 state (and vice versa). Our combined time-dependent density functional theory and quantum dynamics simulation shows that the rISC process from the T1 state to the S1 state can be enhanced by the presence of a nearby local excited triplet state (i.e., T2 state here). A smaller gap between the T1 and T2 states efficiently establishes the rISC route.

The critical role of hot carrier cooling in optically excited structural transitions

The hot carrier cooling occurs in most photoexcitation-induced phase transitions (PIPTs), but its role has often been neglected in many theoretical simulations as well as in proposed mechanisms. Here, by including the previously ignored hot carrier cooling in real-time time-dependent density functional theory (rt-TDDFT) simulations, we investigated the role of hot carrier cooling in PIPTs. Taking IrTe2 as an example, we reveal that the cooling of hot electrons from the higher energy levels of spatially extended states to the lower energy levels of the localized Ir–Ir dimer antibonding states strengthens remarkably the atomic driving forces and enhances atomic kinetic energy. These two factors combine to dissolute the Ir–Ir dimers on a timescale near the limit of atomic motions, thus initiating a deterministic kinetic phase transition. We further demonstrate that the subsequent cooling induces nonradiative recombination of photoexcited electrons and holes, leading to the ultrafast recovery of the Ir–Ir dimers observed experimentally. These findings provide a complete picture of the atomic dynamics in optically excited structural phase transitions.

Circularly Polarized Plasmons in Chiral Gold Nanowires via Quantum-Mechanical Design.

Time-dependent density functional theory (TDDFT) simulations are conducted on a series of chiral gold nanowires to explore whether an enhancement of circular dichroism at the plasmon resonance is possible and identify its quantum-mechanical origin. We find that in linear two-dimensional chiral nanowires the dichroic response is suppressed by destructive interference of nearly degenerate components with opposite signs, pointing to this phenomenon as a common and likely origin of the difficulty encountered so far in achieving a plasmonic CD response in experiment and suggesting nevertheless that these opposite components could be "decoupled" by using multiwall arrangements. In contrast, we predict a giant dichroic response for nanowires with three-dimensional helical coiling. We rationalize this finding via an electronic structure analysis of longitudinal and transversal plasmonic excitations and their coupling into chiral components, and we propose a simple formula for the chiral response as a function of structural parameters (nanowire length and coiling number).

Electron cascades and secondary electron emission in graphene under energetic ion irradiation

Highly energetic ions traversing a two-dimensional material such as graphene produce strong electronic excitations. Electrons excited to energy states above the work function can give rise to secondary electron emission, reducing the amount of energy that remains the graphene after the ion impact. Electrons can either be emitted (kinetic energy transfer) or captured by the passing ion (potential energy transfer). To elucidate this behavior that is absent in three-dimensional materials, we simulate the electron dynamics in graphene during the first femtoseconds after ion impact. We employ two conceptually different computational methods: a Monte Carlo (MC) based one, where electrons are treated as classical particles, and time-dependent density functional theory (TDDFT), where electrons are described quantum-mechanically. We observe that the linear dependence of electron emission on deposited energy, emerging from MC simulations, becomes sublinear and closer to the TDDFT values when the electrostatic interactions of emitted electrons with graphene are taken into account via complementary particle-in-cell simulations. Our TDDFT simulations show that the probability for electron capture decreases rapidly with increasing ion velocity, whereas secondary electron emission dominates in the high velocity regime. We estimate that these processes reduce the amount of energy deposited in the graphene layer by 15\,\% to 65\,\%, depending on the ion and its velocity. This finding clearly shows that electron emission must be taken into consideration when modelling damage production in two-dimensional materials under ion irradiation.

Femtosecond responses of hydrated DNA irradiated by ionizing rays focus on the sugar-phosphate part

In this article, we investigate the mechanisms of DNA ionization upon irradiation by 0.5 meV alpha particles. We focus on the sugar-phosphate group and its hydration shell. In radiation chemistry, the term quasi-direct effect refers the physical and chemical responses taking place after irradiation of solvent molecules pertaining to the solvation shells of solutes. The molecular mechanisms accounting for the quasi-direct effect are actually largely elusive, especially for those prevailing in the early timescales (< 10–12 s). We report Real-Time Time-Dependent Auxiliary Density Functional Theory simulations carried out within the framework of hybrid QM/MM scheme (Quantum Mechanics/Molecular Mechanics) with polarizable and non-polarizable embedding. Ten water molecules from the solvation shell of DNA backbone are independently irradiated. We find that during the first femtoseconds after irradiation, the holes formed on the irradiated water remain at their sites of formation. Electrostatic induction within the environment does not significantly impact charge migrations. We address the hypothesis that charge migration driven by electron correlation is responsible for an ultrafast H2O+ to DNA charge transfers, which would account for a quasi-direct effect. We find that pure charge migration at fixed nuclear positions is not responsible for the quasi-direct effect when considering sugar-phosphate solvation shells.

First-Principles Demonstration of Nonadiabatic Thouless Pumping of Electrons in a Molecular System.

We demonstrate nonadiabatic Thouless pumping of electrons in trans-polyacetylene in the framework of Floquet engineering using first-principles theory. We identify the regimes in which the quantized pump is operative with respect to the driving electric field for a time-dependent Hamiltonian. By employing the time-dependent maximally localized Wannier functions in real-time time-dependent density functional theory simulation, we connect the winding number, a topological invariant, to a molecular-level understanding of the quantized pumping. While the pumping dynamics constitutes the opposing movement of the Wannier functions that represent both double and single bonds, the resulting current is unidirectional due to the greater number of double-bond electrons. Using a gauge-invariant formulation called dynamical transition orbitals, an alternative viewpoint on the nonequilibrium dynamics is obtained in terms of the particle-hole excitation. A single time-dependent transition orbital is found to be largely responsible for the observed quantized pumping. In this representation, the pumping dynamics manifests itself in the dynamics of this single orbital as it undergoes changes from its π bonding orbital character at equilibrium to acquiring resonance and antibonding character in the driving cycle. The work demonstrates the Floquet engineering of the nonadiabatic topological state in an extended molecular system, paving the way for experimental realization of the new quantum material phase.

TD-DFT Simulation and Experimental Studies of a Mirrorless Lasing of Poly[(9,9-dioctylfluorenyl-2,7-diyl)-co-(1,4-diphenylene-vinylene-2-methoxy-5-{2-ethylhexyloxy}-benzene)

In this work, we investigate the TD-DFT simulation, optical, and mirrorless laser properties of conjugated polymer (CP) Poly[(9,9-dioctylfluorenyl-2,7-diyl)-co-(1,4-diphenylene-vinylene-2-methoxy-5-{2-ethylhexyloxy}-benzene)], also known as (PFO-co-PPV-MEHB) or ADS125GE. TD-DFT calculations were performed for three monomer units with truncated tails using time-dependent density functional theory (TD-DFT) calculations. The calculations showed a highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) structure and a very high oscillator strength of 6.434 for the singlet-singlet transition at 374.43 nm. Experimentally, the absorption and fluorescence spectra were examined at various concentrations in verity of solvents, such as benzene, toluene, and hexane. The experimental results obtained in hexane were comparable with theoretical UV-VIS spectra calculated under vacuum. Amplified spontaneous emission (ASE) spectra peaked at approximately 509 nm for CO PFO-co-PPV-MEHB in solution and were obtained at suitable concentrations and pump energies. Additionally, the photochemical stability of this CP and coumarin (C510) were compared. Time-resolved spectroscopy (TRS) studies with a sub-nanosecond resolution were performed for the CO under various pump energies. These results showed the excited state dynamics and single-pass optical gain of CO PFO-co-PPV-MEHB.

Manipulating magnons via ultrafast magnetization modulation

We demonstrate how fundamental properties of magnons may be manipulated using femtosecond laser pulses by performing ab initio real-time time-dependent density functional theory simulations. To illustrate this, we show how the spin-wave dynamics in ${\mathrm{Fe}}_{50}{\mathrm{Ni}}_{50}$ can be manipulated in three different ways by tailoring the applied laser pulse to excite optically induced intersublattice transfer (OISTR): (1) element-selective destruction of magnon modes depending on the laser intensity, (2) delay-dependent freezing of the magnon mode into a transient noncollinear state (where the delay is in the pulse peak with respect to the start of simulations), and (3) OISTR-driven renormalization of the optical magnon frequency. Harnessing such processes would significantly speed up magnonic devices.

Engineering the electronic and optical properties of 2D porphyrin-paddlewheel metal-organic frameworks

Metal-organic frameworks (MOFs) are promising photocatalytic materials due to their high surface area and tuneability of their electronic structure. We discuss here how to engineer the band structures and optical properties of a family of two-dimensional porphyrin-based MOFs, consisting of M-tetrakis(4-carboxyphenyl)porphyrin structures (M-TCPP, where M = Zn or Co) and metal (Co, Ni, Cu or Zn) paddlewheel clusters, with the aim of optimising their photocatalytic behaviour in solar fuel synthesis reactions (water-splitting and/or CO2 reduction). Based on density functional theory (DFT) and time-dependent DFT simulations with a hybrid functional, we studied three types of composition/structural modifications: (a) varying the metal centre at the paddlewheel or at the porphyrin centre to modify the band alignment; (b) partially reducing the porphyrin unit to chlorin, which leads to stronger absorption of visible light; and (c) substituting the benzene bridging between the porphyrin and paddlewheel, by ethyne or butadiyne bridges, with the aim of modifying the linker to metal charge transfer behaviour. Our work offers new insights on how to improve the photocatalytic behaviour of porphyrin- and paddlewheel-based MOFs.

PyRad: A software shell for simulating radiolysis with Qball package.

The assessment of the radiolytic stability of media is an important task in the fields of nuclear power engineering and radiochemistry. Such studies must be carried out in special laboratory conditions with the use of sources of ionizing radiation, which may increase personal doses of the staff. In addition, difficulties arise in studying the products of irradiated media. While it is impossible to abandon experiments to obtain reliable results in this area, computational methods of quantum chemistry can reduce the number of experiments and help understand the mechanisms of the reactions that occur during radiolysis. Here we would like to present a software shell of the Qb@ll program performing time-dependent density functional theory simulations of the radiolysis process.

Efficient Photoexcited Charge Separation at the Interface of a Novel 0D/2D Heterojunction: A Time-Dependent Ultrafast Dynamic Study.

To achieve efficient conversion and avoid loss of solar energy, ultrafast charge separation and slow electron-hole recombination are desired. Combining time-dependent density functional theory (TD-DFT) with nonadiabatic molecular dynamics, Au9(PH3)8/MoS2, as a prototype for zero-dimensional/two-dimensional (0D/2D) heterojunction, has been demonstrated to present excellent light absorption capacity and effective charge separation characteristics. In the heterojunction, photoexcitation of the Au9(PH3)8 nanocluster drives an ultrafast electron transfer from Au9(PH3)8 to MoS2 within 20 fs, whereas photoexcitation of the MoS2 nanosheet leads to hole transfer from MoS2 to Au9(PH3)8 within 680 fs. The strong nonadiabatic coupling and prominent density overlap are responsible for the faster electron separation relative to hole separation. In competition with the charge separation, electron-hole recombination requires 205 ns, ensuring an effective carrier separation. Our atomistic TD-DFT simulation provides valuable insights into the photocarrier dynamics at the Au9(PH3)8/MoS2 interface, which would stimulate the exploration of 0D/2D hybrid materials for photovoltaic and optoelectronic devices.

Correction: Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations.

Correction for 'Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations' by Junjie Yang et al., Phys. Chem. Chem. Phys., 2020, 22, 26838-26851, DOI: .

Synthesis, Structure Determination and Catalytic Activity of a Novel Ruthenium(II) [RuCl(dppb)(44bipy)(4-pic)]PF6 Complex

This work reports the synthesis, structure and catalytic activity of a novel ruthenium(II) complex, [RuCl(dppb)(44bipy)(4-pic)]PF6 (where dppb = 1,4-bis(diphenylphosphine)butane; 44bipy = 4,4’-dimethyl-2,2’-dipyridyl; 4-pic = 4-picoline). The molecular structure and catalytic activity were studied by Fourier transform infrared (FTIR), UV-Vis and nuclear magnetic resonance (NMR) spectroscopies, cyclic voltammetry, and X-ray crystallography, while the electronic structure was investigated by density-functional theory (DFT) and time dependent DFT (TD-DFT) methods. Electrochemical studies showed the substitution of the chlorido ligand from the precursor by the 4-pic ligand, exhibiting the RuII/RuIII process at 1.21 V. The structure of the compound was optimized using DFT simulations and showed data similar to the X-ray structure. The UV-Vis absorption spectrum showed a good agreement with TD-DFT simulations. The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies were determined at the Becke, 3-parameter, Lee-Yang-Parr (B3LYP) level. The study of the catalytic activity in the transfer hydrogenation of ketones by the 1H NMR showed efficient transfer hydrogenation reaction at 60 ºC, employing acetophenone as substrate and resulting in a high conversion. The formation of two ruthenium-hydride species was observed.

Microscopic force driving the photoinduced ultrafast phase transition: Time-dependent density functional theory simulations of IrTe2

Author(s): Liu, WH; Luo, JW; Li, SS; Wang, LW | Abstract: Photoinduced phase transitions can have complex and intriguing behaviors more than material ground-state dynamics. Understanding the underlying mechanism can help us to design new ways to manipulate the materials. A variety of mechanisms has been proposed to explain the photoinduced phase transitions of IrTe2, but a consensus has yet to be reached. Here, we study the photo-induced phase transitions of IrTe2 by performing the real-time time-dependent density functional theory (rt-TDDFT) simulations in combination with the occupation constrained DFT method. We reveal that the microscopic driving force for the photo-induced phase transitions arises from the tendency to lower the energy levels occupied by the excited carriers, through the increase or decrease of the associated atomic pair distances, depending on whether the newly occupied states are antibonding or bonding states, respectively. The geometric constraints between different bonds represented by the Poisson ratio can bring together different tendencies from different atomic pairs, thus forming a complex intriguing dynamic picture depending on the intensity of the excitation. We also find that phonons don't play a primary role, but can assist the phase transition. These findings imply that one can control the structural phase transitions by selectively exciting photocarriers into designated atomic states using appropriate photon sources.

Fluorescence Quantum Yields in Complex Environments from QM-MM TDDFT Simulations: The Case of Indole in Different Solvents.

Fluorescence is commonly exploited to probe microscopic properties. An important example is tryptophan in protein environments, where variations in fluorescence quantum yield, and in absorption and emission maxima, are used as indicators of changes in the environment. Modeling the fluorescence quantum yield requires the determination of both radiative and nonradiative decay constants, both on the potential energy surface of the excited fluorophore. Furthermore, the inclusion of complex environments implies their accurate representation as well as extensive configurational sampling. In this work, we present and test various methodologies based on time-dependent density functional theory (TDDFT) and quantum mechanics/molecular mechanics (QM/MM) dynamics that take all of these requirements into account to provide a quantitative prediction of the effect of the environment on the fluorescence quantum yield of indole, a tryptophan fluorophore. This investigation paves the way for applications to the realistic spectroscopic characterization of the local protein environment of tryptophan from computer simulations.

Ultrafast Optical Modulation of Harmonic Generation in Two-Dimensional Materials.

The modulation of optical harmonic generation in two-dimensional (2D) materials is of paramount importance in nanophotonic and nano-optoelectronic devices for their applications in optical switching and communication. However, an effective route with ultrafast modulation speed, ultrahigh modulation depth, and broad operation wavelength range is awaiting a full exploration. Here, we report that an optical pump can dynamically modulate the third harmonic generation (THG) of a graphene monolayer with a relative modulation depth above 90% at a time scale of 2.5 ps for a broad frequency ranging from near-infrared to ultraviolet. Our observation, together with the real-time, time-dependent density functional theory (TDDFT) simulations, reveals that this modulation process stems from nonlinear dynamics of the photoexcited carriers in graphene. The superior performance of the nonlinear all-optical modulator based on 2D materials paves the way for its potential applications including nanolasers and optical communication circuits.

Twisted Intramolecular Charge Transfer State of a "Push-Pull" Emitter.

The excited state Raman spectra of 4-dicyanomethylene-2-methyl-6-(p-dimethylaminostyryl)-4H-pyran (DCM) in the locally-excited (LE) and the intramolecular charge transfer (ICT) states have been separately measured by time-resolved stimulated Raman spectroscopy. In a polar dimethylsulfoxide solution, the ultrafast ICT of DCM with a time constant of 1.0 ps was observed in addition to the vibrational relaxation in the ICT state of 4–7 ps. On the other hand, the energy of the ICT state of DCM becomes higher than that of the LE state in a less polar chloroform solution, where the initially-photoexcited ICT state with the LE state shows the ultrafast internal conversion to the LE state with a time constant of 300 fs. The excited-state Raman spectra of the LE and ICT state of DCM showed several major vibrational modes of DCM in the LE and ICT conformer states coexisting in the excited state. Comparing to the time-dependent density functional theory simulations and the experimental results of similar push-pull type molecules, a twisted geometry of the dimethylamino group is suggested for the structure of DCM in the S1/ICT state.