Time-dependent Density Functional Theory Approach Research Articles

Dipole instability in molecules irradiated by XUV pulses

We study the response of small covalent molecules to XUV laser pulses. The theoretical description relies on a real-time and real-space Time-Dependent Density Functional Theory (TDDFT) approach at the level of the local density approximation complemented by an efficient self-interaction correction. We observe the development of a dipole instability well after the laser pulse has died out. We find that this instability mechanism is robust with respect to ionic motion, to a wide variety of laser characteristics and to the inclusion of incoherent correlations at the level of a relaxation time ansatz. To rule out any potential numerical effects, we use two independent computational implementations of the TDDFT approach. A comparison of the various laser parameters together with the widely used model approach consisting in an instantaneous hole excitation shows the generic character of this instability in terms of the level depletion of a deep lying electron state. An experimental verification of the phenomenon is proposed in terms of a time-resolved measurement of the photoelectron spectrum.

Role of Dipolar Organic Cations on Light-triggered Charge Transfer at TiO2 /CH3 NH3 PbI3 Interfaces.

The TiO2 /MAPbI3 (MA=CH3 NH3 ) interfaces have manifested correlation with current-voltage hysteresis in perovskite solar cells (PSCs) under light illumination conditions, but the relations between the photo-induced charge transfer and the collective polarization response of the dipolar MA cations are largely unexplored. In this work, we adopt density functional theory (DFT) and time-dependent DFT approach to study the light-triggered charge transfer across the TiO2 /MAPbI3 interfaces with MAI- and PbI-exposed terminations. It is found that regardless of the surface exposure of the MAPbI3 , the photo-induced charge transfer varies when going from the ground-state geometries to the excited-state configurations. Besides, thanks to the electrostatic interactions between the ends of MA cations and the photogenerated electrons, the photo-induced charge transfer across the interfaces is enhanced (weakened) by the negatively (positively) charged CH3 (NH3 ) moieties of the MA species. Resultantly, the positively charged iodine vacancies at the TiO2 /MAPbI3 interfaces tend to inhibit the charge transfer induced by light. Combining with the energy level alignment which is significantly modulated by the orientation of the MA species at the interfaces, the dipolar MA cations might be a double-edge sword for the hysteresis in PSCs with the TiO2 /MAPbI3 interfaces.

Ultrafast Photocontrolled Rotation in a Molecular Motor Investigated by Machine Learning-Based Nonadiabatic Dynamics Simulations.

The thermal helix inversion (THI) of the overcrowded alkene-based molecular motors determines the speed of the unidirectional rotation due to the high reaction barrier in the ground state, in comparison with the ultrafast photoreaction process. Recently, a phosphine-based motor has achieved all-photochemical rotation experimentally, promising to be controlled without a thermal step. However, the mechanism of this photochemical reaction has not yet been fully revealed. The comprehensive computational studies on photoisomerization still resort to nonadiabatic molecular dynamics (NAMD) simulations based on electronic structure calculations, which remains a high computational cost for large systems such as molecular motors. Machine learning (ML) has become an accelerating tool in NAMD simulations recently, where excited-state potential energy surfaces (PESs) are constructed analytically with high accuracy, providing an efficient approach for simulations in photochemistry. Herein the reaction pathway is explored by a spin-flip time-dependent density functional theory (SF-TDDFT) approach in combination with ML-based NAMD simulations. According to our computational simulations, we notice that one of the key factors of fulfilling all-photochemical rotation in the phosphine-based motor is that the excitation energies of four isomers are similar. Additionally, a shortcut photoinduced transformation between unstable isomers replaces the THI step, which shares the conical intersection (CI) with photoisomerization. In this study, we provide a practical approach to speed up the NAMD simulations in photochemical reactions for a large system that could be extended to other complex systems.

Resonance Raman spectra and excited state properties of methyl viologen and its radical cation from time-dependent density functional theory.

Time-dependent density functional theory (TDDFT) was applied to gain insights into the electronic and vibrational spectroscopic properties of an important electron transport mediator, methyl viologen (MV2+ ). An organic dication, MV2+ has numerous applications in electrochemistry that include energy conversion and storage, environmental remediation, and chemical sensing and electrosynthesis. MV2+ is easily reduced by a single electron transfer to form a radical cation species (MV•+ ), which has an intense UV-visible absorption near 600 nm. The redox properties of the MV2+ /MV•+ couple and light-sensitivity of MV•+ have made the system appealing for photo-electrochemical energy conversion (e.g., solar hydrogen generation from water) and the study of photo-induced charge transfer processes through electronic absorption and resonance Raman spectroscopic measurements. The reported work applies leading TDDFT approaches to investigate the electronic and vibrational spectroscopic properties of MV2+ and MV•+ . Using a conventional hybrid exchange functional (B3-LYP) and a long-range corrected hybrid exchange functional (ωB97X-D3), including with a conductor-like polarizable continuum model to account for solvation, the electronic absorption and resonance Raman spectra predicted are in good agreement with experiment. Also analyzed are the charge transfer character and natural transition orbitals derived from the TDDFT vertical excitations calculated. The findings and models developed further the understanding of the electronic properties of viologens and related organic redox mediators important in renewable energy applications and serve as a reference for guiding the interpretation of electronic absorption and Raman spectra of the ions.

Turn on the luminescence channel by the suppression of conical intersection process for the sulfur-containing hydrogen-bonded molecules

The sulfur-containing hydrogen bonds exhibit notable potential in gene and protein engineering. As an important research direction of sulfur-containing hydrogen bonding dynamics, the excited-state intramolecular proton transfer (ESIPT) reaction of the thiol proton has garnered significant interest. Recently, Chou et al. synthesized a series of sulfur-containing H-bonded molecules, inclusive of 3-thiolflavone (3 TF), its trifluoromethyl derivative (3FTF), and its diethylamino derivative (3NTF), thereby inaugurating a new chapter in the ESIPT reaction of the thiol proton. Nevertheless, many open issues cannot be answered solely by experimental means and call for a detailed theoretical investigation, especially the relationship between the molecular excited-state properties and different luminescence properties of the three molecules. Herein, we validate that the molecules undergo the ESIPT reaction of the thiol proton and identify the time scale through the time-dependent density functional theory (TDDFT) calculation. Moreover, employing the spin-flip TDDFT approach, the analysis of the non-radiative transition process reveals that difference in Gibbs free energy of activation (ΔG⧧) for the transition from Frank-Condon point to minimum energy conical intersection (MECI) point primarily accounts for the different luminescence properties of the three molecules. Besides, we also reveal the reason for the abnormal blue-shifted fluorescence of the 3NTF in the aggregated state by the ONIOM (TDDFT: UFF) strategy. These findings are favorable to understanding the luminescence mechanism of sulfur-containing H-bonded molecules and provide a reference for the synthesis and design of the luminescent materials based on the molecules with thiol proton transfer properties.

Modified optoelectronic parameters by end-group engineering of A-D-A type non-fullerene-based small symmetric acceptors constituting IBDT core for high-performance photovoltaics

End-chain exploration of Non-Fullerene Acceptors is an efficient approach to amplify the photovoltaic properties of organic solar cells (OSCs). This green energy technology has attained tremendous consideration from scientists in the last few decades. Current investigations on OSCs have presented remarkable optoelectronic applications of A-D-A-type (A = acceptor and D = donor unit) molecules comprising seven newly designed molecules (INIC1–INIC7) due to their reduced band gaps, red-shifting of optical profiles towards the visible region, Light Harvesting Efficiency (LHE), improved dipole moment, ionization potential (IP), and electron affinity (EA). These parameters are analyzed using the DFT (density functional theory) and TD-DFT (time-dependent-density functional theory) approaches with varying functionals. Numerous specifications like orbital analysis including frontier molecular and natural transition orbital (FMO and NTO) evaluation, as well as density of states (DOS), transition density matrix (TDM), molecular electrostatic potential (MEP) investigation, and various other parameters have demonstrated that the refinement in these factors is attributed to the presence of an IBDT core with seven distinct NFFRAs (Non-Fullerene Fused Ring Acceptors). Amongst the designed molecules, several have exhibited superior features to the modeled molecule. The charge carrier mobilities of different proposed molecules have yielded magnificent results compared to the reference molecule, which is ultimately scrutinized by TDM graphs. The Voc (open-circuit voltage) of INIC5 capitulates preferable impact, which could lead to substitution of all-other molecules in the field of OSCs. Consequently, the impressive dipole moment cooperates in establishing the non-covalent interactions with chloroform leading to exceptional maneuverability. Nonetheless, various derived molecules have exhibited astonishing results concerning different examined parameters.

Studying excited-state-specific perturbation theory on the Thiel set.

We explore the performance of a recently introduced N5-scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set's inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles normoffers a relatively low-cost way to test for doubly excited character without the need to define an active space.

The effect of central transition metals and electron-donating substituent on the performances of dye/TiO2 interface for dye-sensitized solar cells applications

Density functional theory (DFT) and time-dependent density functional theory (TD-DFT) approaches were applied to explore the effect of central transition metals and the dye/TiO2 interface on dye-sensitized solar cell (DSSC) performance and supply a promising way to estimate and screen possible candidates for DSSC applications. The interaction properties, bonding characteristics, sensitized mechanisms, charge transfer, frontier molecular orbitals, energy gap, partial densities of states (PDOS), non-covalent interactions (NCI), and electronic absorption spectra were examined and analyzed to provide the photovoltaic characteristics of Sc, Cu, and Ti tetrasulfonic acid phthalocyanine sensitizer@TiO2 interface in both gas phase and polar solvent as acetonitrile.The interfacial of TiPc-(SO3H)4 with TiO2, which facilitates the driving force for the electron injection of photosensitizers (ΔGinj), increases the charge separation, open-circuit voltage (Voc), and high light-harvesting efficiency (LHE) values. However, minimize the charge recombination, lifetime of the excited state (τ), regeneration driving force (ΔGreg). Our results reveal that TiPc-(SO3H)4@TiO2 is superior to those of ScPc-(SO3H)4, CuPc-(SO3H)4, TiPc-(SO3H)4 and Pc-(SO3H)4/TiO2, indicating the novel TiPc-(SO3H)4@TiO2 interface could be promising candidates for DSSC photovoltaic devices performance. In addition, the static mean polarizability and first hyperpolarizability of all six dyes elucidated that the TiPc-(SO3H)4@TiO2 interface can be regarded as a potential performer in non-linear optical (NLO) properties. These theoretical identifications may provide novel perspectives and instructions for future experimental researchers to promote the synthesis and application of TiPc-(SO3H)4@TiO2 interfaces to improve the photo-to-current conversion efficiency.

Linear-response time-dependent density functional theory approach to warm dense matter with adiabatic exchange-correlation kernels

We present a new methodology for the linear-response time-dependent density functional theory (LR-TDDFT) calculation of the dynamic density response function of warm dense matter in an adiabatic approximation that can be used with any available exchange-correlation (XC) functional across Jacob's Ladder and across temperature regimes. The main novelty of the presented approach is that it can go beyond the adiabatic local density approximation (ALDA) and generalized LDA (AGGA) while preserving the self-consistence between the Kohn-Sham (KS) response function and adiabatic XC kernel for extended systems. The key ingredient for the presented method is the combination of the adiabatic XC kernel from the direct perturbation approach with the macroscopic dynamic KS response from the standard LR-TDDFT method using KS orbitals. We demonstrate the application of the method for the example of warm dense hydrogen, for which we perform a detailed analysis of the KS density response function, the RPA result, the total density response function and of the adiabatic XC kernel. The analysis is performed using LDA, GGA, and meta-GGA level approximations for the XC effects. The presented method is directly applicable to disordered systems such as liquid metals, warm dense matter, and dense plasmas.

Theoretical Study on the Formation and Decomposition Mechanisms of Coelenterazine Dioxetanone.

Bioluminescence has been drawing broad attention due to its high signal-to-noise ratio and high bioluminescence quantum yields, which has been widely applied in the fields of biomedicine, bioanalysis, and so on. Among numerous bioluminescent substrates, coelenterazine is famous for its wide distribution. However, the oxygenation reaction mechanism of coelenterazine is far from being completely understood. In this paper, the formation and decomposition mechanisms of coelenterazine dioxetanone were investigated via density functional theory (DFT) and time-dependent (TD) DFT approaches. The results showed that the oxygenation reaction first occurred along the triplet-state potential energy surface (PES), after the intersystem crossing (ISC), second jumped to the diradical-state PES, and ultimately formed coelenterazine dioxetanone. For the decomposition mechanism of dioxetanone, the computational results showed that the chemiexcitation of neutral dioxetanone was more efficient than that of other dioxetanone species. Moreover, the diradical properties and the degree of ionic character are modified by the counter ions.

Iron(II) Complexes with Porphyrin and Tetrabenzoporphyrin: CASSCF/MCQDPT2 Study of the Electronic Structures and UV-Vis Spectra by sTD-DFT.

The geometry and electronic structures of iron(II) complexes with porphyrin (FeP) and tetrabenzoporphyrin (FeTBP) in ground and low-lying excited electronic states are determined by DFT (PBE0/def2-TZVP) calculations and the complete active space self-consistent field (CASSCF) method, followed by the multiconfigurational quasi-degenerate second-order perturbation theory (MCQDPT2) approach to determine the dynamic electron correlation. The minima on the potential energy surfaces (PESs) of the ground (3A2g) and low-lying, high-spin (5A1g) electronic states correspond to the planar structures of FeP and FeTBP with D4h symmetry. According to the results of the MCQDPT2 calculations, the wave functions of the 3A2g and 5A1g electronic states are single determinant. The electronic absorption (UV-Vis) spectra of FeP and FeTBP are simulated within the framework of the simplified time-dependent density functional theory (sTDDFT) approach with the use of the long-range corrected CAM-B3LYP function. The most intensive bands of the UV-Vis spectra of FeP and FeTBP occur in the Soret near-UV region of 370-390 nm.

Quantum computational investigation into structural, spectroscopic, topological and electronic properties of L-histidinium-L-tartrate hemihydrate: Nonlinear optical organic single crystal

For the first time, a systematic investigation of optimization in the geometrical, vibrational, natural bonding orbital (NBO), electronic, linear and nonlinear optical properties, and Hirshfeld surface analysis for the L-histidinium-l-tartrate hemihydrate (HT) crystal is reported by employing the density functional theory (DFT). The geometrical parameters and vibrational frequencies obtained from the B3LYP/6–311++G(d,p) level of theory are in good agreement with the experimental values. The presence of strong hydrogen bonding interactions in the molecule causes an intense absorption peak in the infrared spectrum below 2000 cm−1. Quantum Theory of Atoms in Molecules (QTAIM) has been used to evaluate the topology of the electron density of a particular molecule to identify the critical points of the system using Multiwfn 3.8. These studies included ELF, LOL, and RDG studies. A time-dependent DFT approach is employed to obtain the excitation energies, oscillator strengths, and UV–Vis spectra for different solvents, such as methanol, ethanol, and water. The NBO analysis of the chosen compound, HT, is performed in terms of atom hybridization and electronic structure. The HOMO-LUMO energies and other associated electronic parameters are also computed. The nucleophilic sites are identified from MEP and Fukui functions analysis. The electrostatic potential and total density of states spectra of HT are discussed in detail. The theoretically obtained polarizability and first order hyperpolarizability values confirm that the grown material HT has NLO efficiency 15.771 times that of urea, and is proposed to be an exceptional nonlinear optical material. In addition, Hirshfeld surface analysis is performed to determine the inter-and intramolecular interactions in the title compound.

Electronic Born-Oppenheimer approximation in nuclear-electronic orbital dynamics.

Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born-Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.

Molecular Structure, Spectral Analysis, Molecular Docking and Physicochemical Studies of 3-Bromo-2-hydroxypyridine Monomer and Dimer as Bromodomain Inhibitors.

In this paper, both methods (DFT and HF) were used in a theoretical investigation of 3-bromo-2-Hydroxypyridine (3-Br-2HyP) molecules where the molecular structures of the title compound have been optimized. Molecular electrostatic potential (MEP) was computed using the B3LYP/6-311++G(d,p) level of theory. The time-dependent density functional theory (TD-DFT) approach was used to simulate the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) on the one hand to achieve the frontier orbital gap and on the other hand to calculate the UV-visible spectrum of the compound in gas phase and for different solvents. In addition, electronic localization function and Fukui functions were carried out. Intermolecular interactions were discussed by the topological AIM (atoms in molecules) approach. The thermodynamic functions have been reported with the help of spectroscopic data using statistical methods revealing the correlations between these functions and temperature. To describe the non-covalent interactions, the reduced density gradient (RDG) analysis is performed. To study the biological activity of the compound of the molecule, molecular docking studies were executed on the active sites of BRD2 inhibitors and to explore the hydrogen bond interaction, minimum binding energies with targeted receptors such as PDB ID: 5IBN, 3U5K, 6CD5 were calculated.

New Polyketide and Sesquiterpenoid Derivatives from the Magellan Seamount-Derived Fungus Penicillium rubens AS-130.

Rubenpolyketone A (1), a polyketide featuring a new carbon skeleton having cyclohexenone condensed with a methyl octenone chain and a new linear sesquiterpenoid, chermesiterpenoid D (2), together with seven known secondary metabolites (3-9) were isolated and identified from the Magellan Seamount-derived fungus Penicillium rubens AS-130. Their structures were determined based on detailed analysis of NMR and mass spectroscopic data and the absolute configurations of these two new compounds were elucidated by the combination of quantum mechanical (QM)-NMR and time-dependent density functional (TDDFT) ECD calculation approaches. Chermesiterpenoids B (3) and C (4) showed potent inhibitory activities against the aquatic pathogen Vibrio anguillarum with MIC values of 0.5 and 1 μg/mL, respectively, while chermesin F (6) exhibited activity against Escherichia coli with MIC value of 1 μg/mL.

Newly Synthesized Aminothiazole Based Disazo Dyes and Their Theoretical Calculations

In this study, ten novel disazo dyes based on aminothiazoles (4a-e, 5a-e) were synthesized in the first section. By using their 1H-NMR, FT-IR, and UV-Vis spectra, these synthesized dyes’ structures were verified. In the second section, the density functional theory (DFT) was used to calculate the molecular geometry, vibrational frequency and molecular electrostatic potential (MEP) diagram in the ground state using the B3LYP level and the 6-311 G(d,p) basis set. UV-vis data were obtained using the Gaussian 09 TD-DFT Cam-B3LYP 6-311 G ++ (d,p) basis set. To explore the chemical shift values, the gauge independent atomic orbital (GIAO) approach was used. Time-dependent DFT (TD-DFT) approach was used to determine the HOMO–LUMO calculations. As a consequence of our research, it was found that the experimental values and the computational spectroscopic values of the synthesized disazo compounds were consistent.

Emergence of Wigner oscillations in a model of real time cooling process: a time-dependent density-functional theory approach

Friedel and Wigner oscillations are well known phenomena occurring in quantum systems. Specifically, in a system composed by N confined fermions, the former are characterized by the presence of peaks in the density distributions, whereas the last by N peaks. Here, we consider N = 2 fermions harmonically confined in one-dimensional quantum dots. It is known that the transition from the Friedel to the Wigner oscillations is induced by the increment of interaction between the fermions. The increment of temperature, on the other hand, acts on eliminating the oscillations. In this context, by employing a time-dependent density-functional theory formalism, we here obtain the emergence of Wigner oscillations in a model which simulates a real time cooling process.

ESIPT-induced spin-orbit coupling enhancement leads to tautomer fluorescence quenching of the 10-HHBF molecule.

We present novel insights into the interplay between excited state intramolecular proton transfer (ESIPT) and spin-orbit coupling (SOC) in the 10-hydroxy-11H-benzo[b]fluoren-11-one (10-HHBF) molecule, utilizing the time-dependent density functional theory approach and femtosecond transient absorption spectroscopy. Our discoveries entail a reassessment of the luminescence mechanism for 10-HHBF, characterizing it as an ESIPT fluorophore. Additionally, we demonstrate that the molecule undergoes intersystem crossing (ISC) following proton transfer, which quenches the fluorescence of the proton-transferred state, thus resulting in the absence of dual emission and a limited spectral range of fluorescence. Furthermore, our investigation reveals that 10-HHBF displays an SOC enhancement feature induced by ESIPT, which facilitates the ISC process. This trait serves as a barrier to the application of 10-HHBF in single-molecule white light emitters (SMWLEs). Our findings underscore the notable influence of the ESIPT-induced spin-orbit interaction enhancement on luminescent properties, which necessitates consideration in the design of SMWLEs.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

We assess most recent TD-DFT methods on excimers and highlight shortcomings of current strategies to treat dispersion in excited states.

XPS and quantum chemical analysis of 4Me-BODIPY derivatives.

The results of a X-ray photoelectron spectroscopy (XPS) and steady-state absorption spectroscopy study of the electronic structure, and cationic and excited states of a series of 1,3,5,7-tetramethyl-substituted BODIPYs (4Me,2R-BODIPYs) are presented. The experimental data were interpreted using high-level ab initio quantum chemical computations, including the algebraic diagrammatic construction method for the polarization propagator of the second order (ADC(2)), the outer-valence Green's function (OVGF) method, the density functional (DFT) approach, and the time-dependent DFT (TD-DFT) approach. Substitution effects on the XPS and absorption spectra were determined for 2,6-positions of 4Me,2R-BODIPY pyrrole nuclei (R = H, Br, Bu, benzyl). A very satisfactory performance of the DFT Koopmans theorem analogue was demonstrated with respect to the energy intervals between the electronic levels of 4Me,2R-BODIPY above 13 eV (BHHLYP functional) and the values of the HOMO-LUMO energy gap (ωB97X functional).