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. ; ; ; ; (, 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
  3. ; ; ; (, 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
  4. ; (, The Journal of Chemical Physics)
    PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule’s electron density. 
    more » « less
  5. ; ; ; ; ; (, The Journal of Physical Chemistry B)
    We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically di!erent depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/ molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture di!erent from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email