Calculating Off-Site Excitations in Symmetric Donor–Acceptor Systems via Time-Dependent Density Functional Theory with Range-Separated Density Functionals

Role of Dielectric Screening in Calculating Excited States of Solvated Azobenzene: A Benchmark Study Comparing Quantum Embedding and Polarizable Continuum Model for Representing the Solvent

The low energy excited states of the conformational isomers of solvated azobenzene are calculated with several DFT methods accounting for the solute-solvent interaction implicitly with the polarizable continuum model or explicitly with subsystem DFT. For the latter, embedding potentials are calculated for 21 sampled snapshots of the solvent molecules. First, we find that accounting for the solvent implicitly or explicitly has little effect on the predicted cis-trans S₁ excitation energy gap. Second, we find that azobenzene's S₁ cis and trans energies are accurate as long as a screened range-separated hybrid exchange-correlation functional is employed. Finally, we also tested a simplified workflow whereby a single, averaged, embedding potential is used. Unfortunately, we find larger deviations against the experiment for the simplified workflow. This highlights a basic flaw in the approach, where the time scale of solvent averaging is much longer than that of the solute's electronic polarization.

Quantum dynamics of the photoinduced charge separation in a symmetric donor–acceptor–donor triad: The role of vibronic couplings, symmetry and temperature

The photoinduced charge separation in a symmetric donor–acceptor–donor (D–A–D) triad is studied quantum mechanically using a realistic diabatic vibronic coupling model. The model includes a locally excited DA*D state and two charge-transfer states D+A−D and DA−D+ and is constructed according to a procedure generally applicable to semirigid D–A–D structures and based on energies, forces, and force constants obtained by quantum chemical calculations. In this case, the electronic structure is described by time-dependent density functional theory, and the corrected linear response is used in conjunction with the polarizable continuum model to account for state-specific solvent effects. The multimode dynamics following the photoexcitation to the locally excited state are simulated by the hybrid Gaussian-multiconfigurational time-dependent Hartree method, and temperature effects are included using thermo field theory. The dynamics are connected to the transient absorption spectrum obtained in recent experiments, which is simulated and fully assigned from first principles. It is found that the charge separation is mediated by symmetry-breaking vibrations of relatively low frequency, which implies that temperature should be accounted for to obtain reliable estimates of the charge transfer rate.

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso-meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4-phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi's golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

How Well Does a Solvated Octa-acid Capsule Shield the Embedded Chromophore? A Computational Analysis Based on an Anisotropic Dielectric Continuum Model

The optical properties of chromophores embedded in a water-solvated dimer of octa-acid that forms a molecular-shaped capsule are investigated. In particular, we address the anisotropic dielectric environment that appears to blue-shift excitation energies compared to the free aqueous chromophores. Recently we reported that using an effective scalar dielectric constant ε ≈ 3 appears to reproduce the measured spectra of the embedded coumarins, suggesting that the capsule provides a significant, albeit not perfect, screening of the aqueous dielectric environment. Here, we report absorption energies using a theoretical treatment that includes continuum solvation affected by an anisotropic dielectric function reflecting the high-dielectric environment outside of the capsule and the low-dielectric region within. We report time-dependent density functional theory calculations using a range-separated functional with the Poisson boundary conditions that model the anisotropic dielectric environment. Our calculations find that the anisotropic environment due to the water-solvated hydrophobic capsule is equivalent to a homogeneous effective dielectric constant of ≈3. The calculated values also appear to reproduce measured absorption of the embedded coumarin, where we study the effect of the hydrophobic capsule on the excited state.

On the Role of the Special Pair in Photosystems as a Charge Transfer Rectifier

The special pair, a bacteriochlorophyll a (BChl) dimer found at the core of bacterial reaction centers, is known to play a key role in the functionality of photosystems as a precursor to the photosynthesis process. In this paper, we analyze the inherent affinity of the special pair to rectify the intrapair photo-induced charge transfer (CT). In particular, we show that the molecular environment affects the nuclear geometry, resulting in symmetry breaking between the two possible intrapair CT processes. To this end, we study the relationships of the intrapair CT and the molecular geometry with respect to the effective dielectric constant provided by the molecular environment. We identify the special pair structural feature that breaks the symmetry between the two molecules, leading to CT rectification. Excited state energies, oscillator strengths, and electronic coupling values are obtained via time-dependent density functional theory, employing a recently developed framework based on a screened range-separated hybrid functional within a polarizable continuum model (SRSH-PCM). We analyze the rectification capability of the special pair by calculating the CT rates using a first-principles-based Fermi's golden rule approach.

Blue M2: an intermediate melanoidin studied via conceptual DFT

In this computational study, ten density functionals, viz. CAM-B3LYP, LC-ω PBE, M11, M11L, MN12L, MN12SX, N12, N12SX, ω B97X, and ω B97XD, related to the Def2TZVP basis sets, are assessed together with the SMD solvation model for calculation of the molecular properties and structure of blue-M2 intermediate melanoidin pigment. All the chemical reactivity descriptors for the system are calculated via conceptual density functional theory (DFT). The active sites suitable for nucleophilic, electrophilic, and radical attacks are selected by linking them with the Fukui function indices, electrophilic Parr functions, and condensed dual descriptors Δf(r), respectively. The prediction of the maximum absorption wavelength is considerably accurate relative to its experimental value. The study reveals that the MN12SX and N12SX density functionals are the most appropriate density functionals for predicting the chemical reactivity of the molecule under study.

Conceptual DFT Study of the Local Chemical Reactivity of the Colored BISARG Melanoidin and Its Protonated Derivative

This computational study assessed eight fixed RSH (range-separated hybrid) density functionals that include CAM-B3LYP, LC-ωPBE, M11, MN12SX, N12SX, ωB97, ωB97X, and ωB97XD related to the Def2TZVP basis sets together with the SMD solvation model in the calculation the molecular structure and reactivity properties of the BISARG intermediate melanoidin pigment (5-(2-(E)-(Z)-5-[(2-furyl)methylidene]-3-(4-acetylamino-4-carboxybutyl)-2-imino-1,3-dihydroimidazol-4-ylideneamino(E)-4-[(2-furyl)methylidene]-5-oxo-1H-imidazol-1-yl)-2-acetylaminovaleric acid) and its protonated derivative, BISARG(p). The chemical reactivity descriptors for the systems were calculated via the Conceptual Density Functional Theory. The choice of active sites applicable to nucleophilic, electrophilic as well as radical attacks were made by linking them with Fukui functions indices, electrophilic and nucleophilic Parr functions, and the condensed Dual Descriptor Δf(r). The study found the MN12SX and N12SX density functionals to be the most appropriate in predicting the chemical reactivity of the molecular systems under study starting from the knowledge of the HOMO, LUMO, and HOMO-LUMO gap energies.

Molecular Reactivity and Absorption Properties of Melanoidin Blue-G1 through Conceptual DFT

This computational study presents the assessment of eleven density functionals that include CAM-B3LYP, LC-wPBE, M11, M11L, MN12L, MN12SX, N12, N12SX, wB97, wB97X and wB97XD related to the Def2TZVP basis sets together with the Solvation Model Density (SMD) solvation model in calculating the molecular properties and structure of the Blue-G1 intermediate melanoidin pigment. The chemical reactivity descriptors for the system are calculated via the conceptual Density Functional Theory (DFT). The choice of the active sites related to the nucleophilic, electrophilic, as well as radical attacks is made by linking them with the Fukui function indices, the electrophilic Parr functions and the condensed dual descriptor Δ f ( r ) . The prediction of the maximum absorption wavelength tends to be considerably accurate relative to its experimental value. The study found the MN12SX and N12SX density functionals to be the most appropriate density functionals in predicting the chemical reactivity of the studied molecule.

Enhancing charge mobilities in organic semiconductors by selective fluorination: a design approach based on a quantum mechanical perspective

Selective fluorination of organic semiconducting molecules is proposed as a means to achieving enhanced hole mobility. Naphthalene is examined here as a root molecular system with fluorination performed at various sites. Our quantum chemical calculations show that selective fluorination can enhance attractive intermolecular interactions while reducing charge trapping. Those observations suggest a design principle whereby fluorination is utilized for achieving high charge mobilities in the crystalline form. The utility of this design principle is demonstrated through an application to perylene, which is an important building block of organic semiconducting materials. We also show that a quantum mechanical perspective of nuclear degrees of freedom is crucial for a reliable description of charge transport.

Simultaneous Determination of Structures, Vibrations, and Frontier Orbital Energies from a Self-Consistent Range-Separated Hybrid Functional

A self-consistent optimally tuned range-separated hybrid density functional (scOT-RSH) approach is developed. It can simultaneously predict accurate geometries, vibrational modes, and frontier orbital energies. This is achieved by optimizing the range-separation parameter, γ, to both satisfy the ionization energy theorem and minimize interatomic forces. We benchmark our approach against an established hybrid functional, B3LYP, using the G2 test set. scOT-RSH greatly improves the accuracy of occupied frontier orbital energies, with a mean absolute error (MAE) of only 0.2 eV relative to experimental ionization energies compared to 2.96 eV with B3LYP. Geometries do not change significantly compared to those obtained from B3LYP, with a bond length MAE of 0.012 Å compared to 0.008 Å for B3LYP, and a 6.5% MAE for zero-point energies, slightly larger than that of B3LYP (3.1%). scOT-RSH represents a new paradigm in which accurate geometries and ionization energies can be predicted simultaneously from a single functional approach.

Pair-Wise and Many-Body Dispersive Interactions Coupled to an Optimally Tuned Range-Separated Hybrid Functional

We propose a nonempirical, pair-wise or many-body dispersion-corrected, optimally tuned range-separated hybrid functional. This functional retains the advantages of the optimal-tuning approach in the prediction of the electronic structure. At the same time, it gains accuracy in the prediction of binding energies for dispersively bound systems, as demonstrated on the S22 and S66 benchmark sets of weakly bound dimers.