Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.
more »
« less
This content will become publicly available on February 5, 2026
AIMD‐Based Protocols for Modeling Exciplex Fluorescence Spectra and Inter‐System Crossing in Photocatalytic Chromophores
ABSTRACT This study introduces a computational protocol for modeling the emission spectra of exciplexes using excited‐state ab initio molecular dynamics (AIMD) simulations. The protocol is applied to a model exciplex formed by oligo‐p‐phenylenes (OPPs) and triethylamine (TEA), which is of interest in the context of photocatalytic reduction of . AIMD facilitates efficient sampling of the conformational space of OPP3 and OPP4 exciplexes with TEA, offering a dynamic alternative to previously employed static methods. The AIMD‐based protocol successfully reproduces experimental emission spectra for OPP‐TEA exciplexes, agreeing with previous computational and experimental findings. The results show that AIMD simulations provide an efficient means of sampling the conformational space of these exciplexes, requiring less user input and, in some instances, fewer computational resources than multiple excited‐state optimizations initiated from user‐specified initial structures. The study also evaluates the yield of intersystem crossing (ISC) using AIMD and Landau‐Zener probability. The results suggest that ISC is a minor decay channel for OPP3 and OPP4. This work provides new insights into the structural flexibility and emission characteristics of OPP‐TEA photoredox catalyst systems, potentially contributing to improved design strategies for organic chromophores in reduction applications.
more »
« less
- Award ID(s):
- 2154482
- PAR ID:
- 10596567
- Publisher / Repository:
- J Comput Chem
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 46
- Issue:
- 4
- ISSN:
- 0192-8651
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This computational study characterises charge-transfer-to-solvent (CTTS) states of aqueous thiocyanate anion using equation-of-motion coupled-cluster methods combined with electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. Equilibrium sampling was carried out using classical molecular dynamics (MD) with standard force-fields and QM/MM ab initio molecular dynamics (AIMD) using density functional theory. The two calculations yield significantly different local structure around solvated SCN− . Because of the diffuse character of CTTS states, they are very sensitive to the local structure of solvent around the solute and its dynamic fluctuations. Owing to this sensitivity, the spectra computed using MD and AIMD based snapshots differ considerably. This sensitivity suggests that the spectroscopy exploiting CTTS transitions can provide an experimental handle for assessing the quality of force-fields and density functionals. By combining CTTS-based spectroscopies with reliable theoretical modeling, detailed microscopic information of the solvent structure can be obtained. We present a robust computational protocol for modeling spectra of solvated anions and emphasise the use of an ab initio characterization of individual electronic transitions as CTTS or local excitations.more » « less
-
Shea, Joan-Emma (Ed.)We aim to automatize the identification of collective variables to simplify and speed up enhanced sampling simulations of conformational dynamics in biomolecules. We focus on anharmonic low-frequency vibrations that exhibit fluctuations on timescales faster than conformational transitions but describe a path of least resistance towards structural change. A key challenge is that harmonic approximations are ill-suited to characterize these vibrations, which are observed at far-infrared frequencies and are easily excited by thermal collisions at room temperature. Here, we approached this problem with a frequency-selective anharmonic (FRESEAN) mode analysis that does not rely on harmonic approximations and successfully isolates anharmonic low-frequency vibrations from short molecular dynamics simulation trajectories. We applied FRESEAN mode analysis to simulations of alanine dipeptide, a common test system for enhanced sampling simulation protocols, and compare the performance of isolated low-frequency vibrations to conventional user-defined collective variables (here backbone dihedral angles) in enhanced sampling simulations. The comparison shows that enhanced sampling along anharmonic low-frequency vibrations not only reproduces known conformational dynamics but can even further improve sampling of slow transitions compared to user-defined collective variables. Notably, free energy surfaces spanned by low-frequency anharmonic vibrational modes exhibit lower barriers associated with conformational transitions relative to representations in backbone dihedral space. We thus conclude that anharmonic low-frequency vibrations provide a promising path for highly effective and fully automated enhanced sampling simulations of conformational dynamics in biomolecules.more » « less
-
Abstract ortho‐Phosphinophenol (oPP) is an unusual example of an air‐stable primary phosphine and a valuable precursor to a variety of useful organophosphorus compounds. The presence of PH2and OH functionalities offers the possibility of intermolecular and intramolecular P⋅⋅⋅HO hydrogen bonding (HB). The close proximity of these two groups also offers the opportunity for intramolecular PH2⋅⋅⋅HO dihydrogen bonding (DHB). This work provides experimental and computational evidence for these various types of interactions. In the solid state,oPPis associated by significant intermolecular P⋅⋅⋅HO hydrogen bonds as revealed by a single crystal X‐ray structural determination. Multinuclear NMR and IR spectroscopic studies, coupled with DFT computational studies, suggest thatoPPadopts multiple conformations in solution whose nature varies with the identity of the solvent. In the gas phase or non‐polar solvents (such as cyclohexane) an equilibrium between four conformations ofoPPis proposed. Interestingly,in silico, a conformational isomer having bifurcated intramolecular PH2⋅⋅⋅HODHB(PP4) is found to be more stable than a conformational isomer having intramolecular P⋅⋅⋅HOHB(PP1). In polar solvents (S), NMR studies indicate intermolecular OH⋅⋅⋅S HBplays a dominant role in modulating31P NMR chemical shifts over a 17 ppm range.more » « less
-
The photophysical properties of 2,4-dithiouracil (2,4-DTU) in the gas phase are studied by time-resolved photoelectron spectroscopy (TRPES) with three different excitation wavelengths in direct extension of previous work on uracil (U), 2-thiouracil (2-TU) and 4-thiouracil (4-TU). Non-radiative deactivation in the canonical nucleobases like uracil mainly occurs via internal conversion (IC) along singlet excited states, although intersystem crossing (ISC) to a long-lived triplet state was confirmed to play a minor role. In thionated uracils, ISC to the triplet state becomes ultrafast and highly efficient with a quantum yield near unity; however, the lifetime of the triplet state is strongly dependent on the position of the sulfur atom. In 2-TU, ISC back to the ground state occurs within a few hundred picoseconds, whereas the population remains trapped in the lowest triplet state in the case of 4-TU. Upon doubling the degree of thionation, ISC remains highly efficient and dominates the photophysics of 2,4-DTU. However, several low-lying excited states contribute to competing IC and ISC pathways and a complex deactivation mechanism, which is evaluated here based on TRPES measurements and discussed in the context of the singly thionated uracils.more » « less
An official website of the United States government
