Search

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin–orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm–Dancoff approximation to TDDFT, spin–flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine–Ag20 complex), we emphasize that caution is needed to properly account for system–environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In order to understand the mechanism of the electrochemical reduction of CO* on Cu(111), which competitively generates two intermediates, CHO* and COH*, we performed a first principles calculation on these two electrocatalytic reactions, including the solvent effect and imposing a constant potential. The transition states of the two reactions under constant potential conditions were located by the electrochemical nudged elastic band (eNEB) method and the charge effect in the two reactions was elucidated by charge correction in the constant potential model and Bader charge analysis. It was found that the two reactions have different potential dependencies and involve different degrees of partial charge transfer. In addition, we confirmed that the COH* reaction pathway exists only in a certain potential range when using the implicit solvent model. [ABSTRACT FROM AUTHOR]

It is of great importance and worthy of efforts to give a clear structure–property relationship and microscopic mechanism of fluorescence emitters with high quantum yield. In this work, we perform a detailed computational investigation to give an explanation to the high efficiency of a fluorescence emitter XBTD‐NPh based TADF sensitized fluorescence (TSF) OLEDs, and construct a symmetry structure DSBNA‐BTD. Theoretical calculations show that XBTD‐NPh is a long‐time phosphorescent material at 77 K and TADF is attributed to the RISC of T1 to S1 state. For DSBNA‐BTD, excitons arrived at T1 state comes to a large rate of nonradiatively path to the ground state, meaning it is may not be an efficient TADF molecule. For both molecules, the fast IC between T2 and T1 state results in that the hot exciton channel T1‐Tn‐S1 makes no contribution to the TADF. [ABSTRACT FROM AUTHOR]

The nonadiabatic phenomena, which are characterized by a strong coupling between electronic and nuclear motions, are ubiquitous. The nonadiabatic effect of the studied system can be significantly affected by the surrounding environment, such as solvents, in which such nonadiabatic process takes place. It is essential to develop the theoretical models to simulate these processes while accurately modeling the solvent environment. The time-dependent density functional theory (TDDFT) is currently the most efficient approach to describe the electronic structures and dynamics of complex systems, while the polarizable continuum model (PCM) represents one of the most successful examples among continuum solvation models. Here, we formulate the first-order derivative couplings (DCs) between the ground and excited states as well as between two excited states by utilizing time-independent equation of motion formalism within the framework of both linear response and spin flip formulations of TDDFT/CPCM (the conductor-like PCM), and implement the analytical DCs into the Q-CHEM electronic structure software package. The analytic implementation is validated by the comparison of the analytical and finite-difference results, and reproducing geometric phase effect in the protonated formaldimine test case. Taking 4-(N,N-dimethylamino)benzonitrile and uracil in the gas phase and solution as an example, we demonstrate that the solvent effect is essential not only for the excitation energies of the low-lying excited-states but also for the DCs between these states. Finally, we calculate the internal conversion rate of benzophenone in a solvent with DC being used. The current implementation of analytical DCs together with the existing analytical gradient and Hessian of TDDFT/PCM excited states allows one to study the nonadiabatic effects of relatively large systems in solutions with low computational cost. [ABSTRACT FROM AUTHOR]

This work is devoted to deriving and implementing analytic second- and third-order energy derivatives with respect to the nuclear coordinates and external electric field within the framework of the hybrid quantum mechanics/molecular mechanics method with induced charges and dipoles (QM/DIM). Using these analytic energy derivatives, one can efficiently compute the harmonic vibrational frequencies, infrared (IR) and Raman scattering (RS) spectra of the molecule in the proximity of noble metal clusters/nanoparticles. The validity and accuracy of these analytic implementations are demonstrated by the comparison of results obtained by the finite-difference method and the analytic approaches and by the full QM and QM/DIM calculations. The complexes formed by pyridine and two sizes of gold clusters (Au18 and Au32) at varying intersystem distances of 3, 4, and 5 Å are used as the test systems, and Raman spectra of 4,4′-bipyridine in the proximity of Au2057 and Ag2057 metal nanoparticles (MNP) are calculated by the QM/DIM method and compared with experimental results as well. We find that the QM/DIM model can well reproduce the IR spectra obtained from full QM calculations for all the configurations, while although it properly enhances some of the vibrational modes, it artificially overestimates RS spectral intensities of several modes for the systems with very short intersystem distance. We show that this could be improved, however, by incorporating the hyperpolarizability of the gold metal cluster in the evaluation of RS intensities. Additionally, we address the potential impact of charge migration between the adsorbate and MNPs. [ABSTRACT FROM AUTHOR]

We comment on an excited‐state localization method recently proposed by Blanc et al. (J. Comput. Chem. 2023, 44, 105). Elaborate comparisons are made to demonstrate that their method is a less‐comprehensive version of the diabatization method proposed by us 2 years earlier (J. Phys. Chem. Lett. 2021, 12, 1032). [ABSTRACT FROM AUTHOR]

Natural light‐harvesting (LH) systems can divide identical dyes into unequal aggregate states, thereby achieving intelligent "allocation of labor". From a synthetic point of view, the construction of such kinds of unequal and integrated systems without the help of proteinaceous scaffolding is challenging. Here, we show that four octatetrayne‐bridged ortho‐perylene diimide (PDI) dyads (POPs) self‐assemble into a quadruple assembly (POP)4 both in solution and in the solid state. The two identical PDI units in each POP are compartmentalized into weakly coupled PDIs (P520) and closely stacked PDIs (P550) in (POP)4. The two extreme pools of PDI chromophores were unambiguously confirmed by single‐crystal X‐ray crystallography and NMR spectroscopy. To interpret the formation of the discrete quadruple assembly, we also developed a two‐step cooperative model. Quantum‐chemical calculations indicate the existence of multiple couplings within and across P520 and P550, which can satisfactorily describe the photophysical properties of the unequal quadruple assembly. This finding is expected to help advance the rational design of dye stacks to emulate functions of natural LH systems. [ABSTRACT FROM AUTHOR]

Natural light‐harvesting (LH) systems can divide identical dyes into unequal aggregate states, thereby achieving intelligent "allocation of labor". From a synthetic point of view, the construction of such kinds of unequal and integrated systems without the help of proteinaceous scaffolding is challenging. Here, we show that four octatetrayne‐bridged ortho‐perylene diimide (PDI) dyads (POPs) self‐assemble into a quadruple assembly (POP)4 both in solution and in the solid state. The two identical PDI units in each POP are compartmentalized into weakly coupled PDIs (P520) and closely stacked PDIs (P550) in (POP)4. The two extreme pools of PDI chromophores were unambiguously confirmed by single‐crystal X‐ray crystallography and NMR spectroscopy. To interpret the formation of the discrete quadruple assembly, we also developed a two‐step cooperative model. Quantum‐chemical calculations indicate the existence of multiple couplings within and across P520 and P550, which can satisfactorily describe the photophysical properties of the unequal quadruple assembly. This finding is expected to help advance the rational design of dye stacks to emulate functions of natural LH systems. [ABSTRACT FROM AUTHOR]