skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges
Title: Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges
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–Ag 20 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.  more » « less
Award ID(s):
2102071
PAR ID:
10337062
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
156
Issue:
21
ISSN:
0021-9606
Page Range / eLocation ID:
210901
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Molecules)
    FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin’s isoalloxazine ring in the diagnostic 1450–1750 cm−1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution. 
    more » « less
  2. (, International Journal of Quantum Chemistry)
    Abstract The environment may significantly affect molecular properties. Thus, it is desirable to account explicitly for these effects on the wave function and its derivatives, especially when the latter are evaluated with accurate methods, such as those belonging to coupled cluster (CC) theory. In this tutorial review, we discuss how to combine CC methods with the polarizable continuum model of solvation (PCM). We describe useful approximations that include the solvent response to the correlation and excited state equations while maintaining the computational cost comparable to in vacuo calculations. Although applied to PCM, the theoretical framework presented in this review is general and can be used with any polarizable embedding model. Representative applications of the CC‐PCM method to ground and excited state properties of solvated molecules are presented, and comparisons with experiment, and between the full and approximate schemes are discussed. 
    more » « less
  3. ; ; (, WIREs Computational Molecular Science)
    Abstract Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed‐phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long‐range effects, to polarization and advanced embedding potentials. This article is categorized under:Electronic Structure Theory > Combined QM/MM MethodsMolecular and Statistical Mechanics > Molecular InteractionsSoftware > Simulation Methods 
    more » « less
  4. ; ; ; ; (, Journal of Computational Chemistry)
    Abstract 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. 
    more » « less
  5. ; ; ; (, Journal of physical chemistry C)
    The robustness of nickelocene’s (NiCp2, Cp = cyclopentadienyl) magnetic anisotropy and addressability of its spin states make this molecular magnet attractive as a spin sensor. However, microscopic understanding of its magnetic anisotropy is still lacking, especially when NiCp2 is deposited on a surface to make quantum sensing devices. Quantum chemical calculations of such molecule/solid-state systems are limited to density functional theory (DFT) or DFT+U (Hubbard correction to DFT). We investigate the magnetic behavior of NiCp2 using the spin-flip variant of the equation-of-motion coupled-cluster (EOM-SF-CC) method and use the EOM-SF-CC results to benchmark SF-TD-DFT. Our first-principle calculations agree well with experimentally derived magnetic anisotropy and susceptibility values. The calculations show that magnetic anisotropy in NiCp2 originates from a large spin–orbit coupling (SOC) between the triplet ground state and the third singlet state, whereas the coupling with lower singlet excited states is negligible. We also considered a set of six ring-substituted NiCp2 derivatives and a model system of the NiCp2/MgO(001) adsorption complex, for which we used SF-TD-DFT method. To gain insight into the electronic structure of these systems, we analyze spinless transition density matrices and their natural transition orbitals (NTOs). The NTO analysis of SOCs explains how spin states and magnetic properties are retained upon modification of the NiCp2 coordination environment and upon its adsorption on a surface. Such resilience of the NiCp2 magnetic behavior supports using NiCp2 as a spin-probe molecule by functionalization of the tip of a scanning tunneling microscope. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email