Computational Design of (B)Chl Models: Structural and Chemical Modifications toward Enriched Properties.

The functional units of natural photosynthetic systems control the process of converting sunlight into chemical energy. In this article, we explore a series of chemically and structurally modified bacteriochlorophyll and chlorophyll pigments through computational chemistry to evaluate their electronic spectroscopy properties. More specifically, we use multiconfigurational and time-dependent density functional theory methods, along with molecular dynamics simulations, to compute the models' energetics both in an implicit and explicit solvent environment.

Discriminating Clockwise and Counterclockwise Photoisomerization Paths in Achiral Photoswitches by Excited-State Electronic Circular Dichroism.

Despite the numerous investigations of photoisomerization reactions from both the computational and experimental points of view, even in complex environments, to date there is no direct demonstration of the direction of rotation of the retinal chromophore, initiating the vision process in several organisms, occurring upon light irradiation. In the literature, many proposals have been formulated to shed light on the details of this process, most of which are extracted from semiclassical simulations. Although high hopes are held in the development of time-resolved X-ray spectroscopy, I argue in this work that simpler but less known techniques can be used to unravel the details of this fascinating photochemical process.

A novel bioresponsive self-immolative spacer based on aza-quinone methide reactivity for the controlled release of thiols, phenols, amines, sulfonamides or amides.

A stimuli-sensitive linker is one of the indispensable components of prodrugs for cancer therapy as it covalently binds the drug and releases it upon external stimulation at the tumour site. Quinone methide elimination has been widely used as the key transformation to release drugs based on their nucleofugacity. The usual approach is to bind the drug to the linker as a carbamate and release it as a free amine after a self-immolative 1,6-elimination.

A Computational Perspective on the Reactivity of π-spacers in Self-Immolative Elimination Reactions.

The controlled release of chemicals, especially in drug delivery, is crucial, often employing "self-immolative" spacers to enhance reliability. These spacers separate the payload from the protecting group, ensuring a more controlled release. Over the years, design rules have been proposed to improve the elimination process's reaction rate by modifying spacers with electron-donating groups or reducing their aromaticity.

Formation of electron traps in semiconducting polymers via a slow triple-encounter between trap precursor particles.

Already in 2012, Blom et al. reported (Nature Materials 2012, , 882) in semiconducting polymers on a general electron-trap density of ≈3 × 10 cm, centered at an energy of ≈3.6 eV below vacuum.

Identification of Unknown Inverted Singlet-Triplet Cores by High-Throughput Virtual Screening.

Molecules where the energy of the lowest excited singlet state is found below the energy of the lowest triplet state (inverted singlet-triplet molecules) are extremely rare. It is particularly challenging to discover new ones through virtual screening because the required wavefunction-based methods are expensive and unsuitable for high-throughput calculations. Here, we devised a virtual screening approach where the molecules to be considered with advanced methods are pre-selected with increasingly more sophisticated filters that include the evaluation of the HOMO-LUMO exchange integral and approximate CASSCF calculations.

The OpenMolcas : A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Automated QM/MM Screening of Rhodopsin Variants with Enhanced Fluorescence.

We present a computational protocol for the fast and automated screening of excited-state hybrid quantum mechanics/molecular mechanics (QM/MM) models of rhodopsins to be used as fluorescent probes based on the automatic rhodopsin modeling protocol (-ARM). Such "-ARM fluorescence screening protocol" is implemented through a general Python-based driver, PyARM, that is also proposed here. The implementation and performance of the protocol are benchmarked using different sets of rhodopsin variants whose absorption and, more relevantly, emission spectra have been experimentally measured.

Retinal chromophore charge delocalization and confinement explain the extreme photophysics of Neorhodopsin.

The understanding of how the rhodopsin sequence can be modified to exactly modulate the spectroscopic properties of its retinal chromophore, is a prerequisite for the rational design of more effective optogenetic tools. One key problem is that of establishing the rules to be satisfied for achieving highly fluorescent rhodopsins with a near infrared absorption. In the present paper we use multi-configurational quantum chemistry to construct a computer model of a recently discovered natural rhodopsin, Neorhodopsin, displaying exactly such properties.

On the fluorescence enhancement of arch neuronal optogenetic reporters.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore.

Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter.

The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility.

Exciton transport in molecular organic semiconductors boosted by transient quantum delocalization.

Designing molecular materials with very large exciton diffusion lengths would remove some of the intrinsic limitations of present-day organic optoelectronic devices. Yet, the nature of excitons in these materials is still not sufficiently well understood. Here we present Frenkel exciton surface hopping, an efficient method to propagate excitons through truly nano-scale materials by solving the time-dependent Schrödinger equation coupled to nuclear motion.

Organic materials repurposing, a data set for theoretical predictions of new applications for existing compounds.

We present a data set of 48182 organic semiconductors, constituted of molecules that were prepared with a documented synthetic pathway and are stable in solid state. We based our search on the Cambridge Structural Database, from which we selected semiconductors with a computational funnel procedure. For each entry we provide a set of electronic properties relevant for organic materials research, and the electronic wavefunction for further calculations and/or analyses.

Bright Frenkel Excitons in Molecular Crystals: A Survey.

We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models.

Computational Identification of Novel Families of Nonfullerene Acceptors by Modification of Known Compounds.

We considered a database of tens of thousands of known organic semiconductors and identified those compounds with computed electronic properties (orbital energies, excited state energies, and oscillator strengths) that would make them suitable as nonfullerene electron acceptors in organic solar cells. The range of parameters for the desirable acceptors is determined from a set of experimentally characterized high-efficiency nonfullerene acceptors. This search leads to ∼30 lead compounds never considered before for organic photovoltaic applications.

Shortwave infrared-absorbing squaraine dyes for all-organic optical upconversion devices.

Shortwave infrared (SWIR) optical sensing and imaging are essential to an increasing number of next-generation applications in communications, process control or medical imaging. An all-organic SWIR upconversion device (OUC) consists of an organic SWIR sensitive photodetector (PD) and an organic light-emitting diode (OLED), connected in series. OUCs directly convert SWIR to visible photons, which potentially provides a low-cost alternative to the current inorganic compound-based SWIR imaging technology.

Amplification of the chiroptical response of UV-transparent amines and alcohols by N-phthalimide derivatization enabling absolute configuration determination through ECD computational analysis.

The stereoselective transformation of chiral UV-transparent amines and alcohols to phthalimides has proved to be a simple and efficient method to enhance the chiroptical response of these substrates allowing their reliable absolute configuration determination by computational analysis of ECD spectra. Such a transformation also leads to a significant reduction in the molecular conformational flexibility thus simplifying the conformational analysis required by the computational treatment. The method described herein thus allows the absolute configuration assignment to these challenging substrates to be much easier and reliable.

On the arrangement of chromophores in light harvesting complexes: chance versus design.

We used a homogeneous computational approach to derive the excitonic Hamiltonian for five light harvesting complexes containing only one type of chromophore and compare them in terms of statistical descriptors. We then studied the approximate exciton dynamics for the five complexes introducing a measure, the (averaged and time-dependent) inverse participation ratio, that enables the comparison between very diverse complexes on the same ground. We find that the global dynamics are very similar across the set of systems despite the variety of geometric structures of the complexes.

How fine-tuned for energy transfer is the environmental noise produced by proteins around biological chromophores?

We investigate the role of the local protein environment on the energy transfer processes in biological molecules, excluding from the analysis the effect of intra-chromophore nuclear motions, and focussing on the exciton-phonon coupling. We studied three different proteins (FMO and two variants of the WSCP protein) with different biological functions but similar chromophores, to understand whether a classification of chromophores based on the details of the environment would be possible, and whether specific environments enhance or suppress the coupling between exciton and protein dynamics. Our results show that despite the different biological role, there is no significant difference in the influence of the environment on the properties of the chromophores.

How and How Much Molecular Conformation Affects Electronic Circular Dichroism: The Case of 1,1-Diarylcarbinols.

Chiroptical spectra such as electronic circular dichroism (ECD) are said to be much more sensitive to conformation than their non-chiroptical counterparts, however, it is difficult to demonstrate such a common notion in a clear-cut way. We run DFT and TDDFT calculations on two closely related 1,1-diarylmethanols which show mirror-image ECD spectra for the same absolute configuration. We demonstrate that the main reason for the different chiroptical response of the two compounds lies in different conformational ensembles, caused by a single hydrogen-to-methyl substitution.

The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects.

Comparison between chiroptical spectra and theoretical predictions is the method of choice for the assignment of the absolute configuration of chiral compounds in solution. Here we report the case of an apparently simple biarylcarbinol, whose electronic circular dichroism (ECD) in the L region exhibits a peculiar alternation of negative and positive bands. Adopting Density Functional Theory, and describing solvent effects with implicit methods, we found three stable conformers in ethanol, each of them with two close lying states corresponding to similar local L excitations on the two phenyls.

Chromophore-Dependent Intramolecular Exciton-Vibrational Coupling in the FMO Complex: Quantification and Importance for Exciton Dynamics.

In this paper, we adopt an approach suitable for monitoring the time evolution of the intramolecular contribution to the spectral density of a set of identical chromophores embedded in their respective environments. We apply the proposed method to the Fenna-Matthews-Olson (FMO) complex, with the objective to quantify the differences among site-dependent spectral densities and the impact of such differences on the exciton dynamics of the system. Our approach takes advantage of the vertical gradient approximation to reduce the computational demands of the normal modes analysis.