Complex Frequency-Dependent Polarizability through the π → π* Excitation Energy of Azobenzene Molecules by a Combined Charge-Transfer and Point-Dipole Interaction Model

A new framework for frequency-dependent polarizable force fields

A frequency-dependent extension of the polarizable force field "Atom-Condensed Kohn-Sham density functional theory approximated to the second-order" (ACKS2) [Verstraelen et al., J. Chem. Phys. 141, 194114 (2014)] is proposed, referred to as ACKS2ω. The method enables theoretical predictions of dynamical response properties of finite systems after partitioning of the frequency-dependent molecular response function. Parameters in this model are computed simply as expectation values of an electronic wavefunction, and the hardness matrix is entirely reused from ACKS2 as an adiabatic approximation is used. A numerical validation shows that accurate models can already be obtained with atomic monopoles and dipoles. Absorption spectra of 42 organic and inorganic molecular monomers are evaluated using ACKS2ω, and our results agree well with the time-dependent DFT calculations. Also for the calculation of C₆ dispersion coefficients, ACKS2ω closely reproduces its TDDFT reference. When parameters for ACKS2ω are derived from a PBE/aug-cc-pVDZ ground state, it reproduces experimental values for 903 organic and inorganic intermolecular pairs with an MAPE of 3.84%. Our results confirm that ACKS2ω offers a solid connection between the quantum-mechanical description of frequency-dependent response and computationally efficient force-field models.

Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug

Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.

Evaluation of Molecular Polarizability and of Intensity Carrying Modes Contributions in Circular Dichroism Spectroscopies

We re-examine the theory of electronic and vibrational circular dichroism spectroscopy in terms of the formalism of frequency-dependent molecular polarizabilities. We show the link between Fermi’s gold rule in circular dichroism and the trace of the complex electric dipole–magnetic dipole polarizability. We introduce the C++ code polar to compute the molecular polarizability complex tensors from quantum chemistry outputs, thus simulating straightforwardly UV-visible absorption (UV-Vis)/electronic circular dichroism (ECD) spectra, and infrared (IR)/vibrational circular dichroism (VCD) spectra. We validate the theory and the code by referring to literature data of a large group of chiral molecules, showing the remarkable accuracy of density functional theory (DFT) methods. We anticipate the application of this methodology to the interpretation of vibrational spectra in various measurement conditions, even in presence of metal surfaces with plasmonic properties. Our theoretical developments aim, in the long run, at embedding the quantum-mechanical details of the chiroptical spectroscopic response of a molecule into the simulation of the electromagnetic field distribution at the surface of plasmonic devices. Such simulations are also instrumental to the interpretation of the experimental spectra measured from devices designed to enhance chiroptical interactions by the surface plasmon resonance of metal nanostructures.

Germanium Fluoride Nanocages as Optically Transparent n-Type Materials and Their Endohedral Metallofullerene Derivatives

Carbon- and silicon-based n-type materials tend to suffer from instability of the corresponding radical anions. With DFT calculations, we explore a promising route to overcome such challenges with molecular nanocages which utilize the heavier element Ge. The addition of fluorine substituents creates large electron affinities in the range 2.5-5.5 eV and HOMO-LUMO gaps between 1.6 and 3.2 eV. The LUMOs envelop the surfaces of these structures, suggesting extensive delocalization of injected electrons, analogous to fullerene acceptors. Moreover, these Ge _nF _n inorganic cages are found to be transparent in the UV-visible region as probed with their excited states. Their capacitance, linear polarizabilities, and dielectric constants are computed and found to be on the same order of magnitude as saturated oligomers and some extended π-organics (azobenzenes). Furthermore, we explore fullerene-type endohedral isomers, i.e., cages with internal substituents or guest atoms, and find them to be more stable than the parent exohedral isomers by up to -206.45 kcal mol^-1. We also consider the addition of Li, He, Cs, and Bi, to probe the utility of the exo/ endo cages as host-guest systems. The endohedral He/Li@F₈@Ge₆₀F₅₂ cages are significantly more stable than their parent exohedral isomers He/Li@Ge₆₀F₅₂ by -182.46 and -49.22 kcal mol^-1, respectively. The energy of formation of endohedral He@F₈@Ge₆₀F₅₂ is exothermic by -10.4 kcal mol^-1, while Cs and Bi guests are too large to be accommodated but are stable in the exohedral parent cages. Conceivable applications of these materials include n-type semiconductors and transparent electrodes, with potential for novel energy storage modalities.

Local electric field factors by a combined charge-transfer and point–dipole interaction model

A model for the local electric field as a linear response to a frequency-dependent external electric field is presented based on a combined charge-transfer and point–dipole interaction force-field model.