Modelling solvent effects on the absorption and emission spectra of constrained cyanines with both implicit and explicit QM/EFP models

Quantum Simulation of Molecules in Solution

Quantum chemical calculations on quantum computers have been focused mostly on simulating molecules in the gas phase. Molecules in liquid solution are, however, most relevant for chemistry. Continuum solvation models represent a good compromise between computational affordability and accuracy in describing solvation effects within a quantum chemical description of solute molecules. In this work, we extend the variational quantum eigensolver to simulate solvated systems using the polarizable continuum model. To account for the state dependent solute-solvent interaction we generalize the variational quantum eigensolver algorithm to treat non-linear molecular Hamiltonians. We show that including solvation effects does not impact the algorithmic efficiency. Numerical results of noiseless simulations for molecular systems with up to 12 spin-orbitals (qubits) are presented. Furthermore, calculations performed on a simulated noisy quantum hardware (IBM Q, Mumbai) yield computed solvation free energies in fair agreement with the classical calculations.

The Impact of Retinal Configuration on the Protein-Chromophore Interactions in Bistable Jumping Spider Rhodopsin-1

Bistable rhodopsins have two stable forms that can be interconverted by light. Due to their ability to act as photoswitches, these proteins are considered as ideal candidates for applications such as optogenetics. In this work, we analyze a recently crystalized bistable rhodopsin, namely the jumping spider rhodopsin-1 (JSR1). This rhodopsin exhibits identical absorption maxima for the parent and the photoproduct form, which impedes its broad application. We performed hybrid QM/MM simulations to study three isomers of the retinal chromophore: the 9-cis, 11-cis and all-trans configurations. The main aim was to gain insight into the specific interactions of each isomer and their impact on the absorption maximum in JSR1. The absorption spectra were computed using sampled snapshots from QM/MM molecular dynamics trajectories and compared to their experimental counterparts. The chromophore-protein interactions were analyzed by visualizing the electrostatic potential of the protein and projecting it onto the chromophore. It was found that the distance between a nearby tyrosine (Y126) residue plays a larger role in the predicted absorption maximum than the primary counterion (E194). Geometric differences between the isomers were also noted, including a structural change in the polyene chain of the chromophore, as well as changes in the nearby hydrogen bonding network.

New insights into quantifying the solvatochromism of BODIPY based fluorescent probes

A simple semiempiric phenomenological approach is developed for quantifying the solvent effect on the absorption and emission properties of BODIPYs. It is based on a new rule describing the linear relationship between the difference (Stokes shift) and the sum (double Gibbs free energy of electron transfer) for absorption and emission wavenumbers derived from a combination of solvent functions of Liptay theory. This rule is correspondent to changes of dipole moments in the ground and excited states. High reliability and advantages of the developed approach in comparison with traditional methods of the analysis of the solvatochromism based on Dimroth-Reichard and Lippert-Mataga solvent scales are illustrated for selected BODIPYs exhibiting positive, negative, and near-zero solvatochromism.

The TDDFT Excitation Energies of the BODIPYs; The DFT and TDDFT Challenge Continues

The derivatives of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY) are pivotal ingredients for a large number of functional, stimuli-responsive materials and therapeutic molecules based on their photophysical properties, and there is a urgent need to understand and predict their optical traits prior to investing a large amount of resources in preparing them. Density functional theory (DFT) and time-dependent DFT (TDDFT) computations were performed to calculate the excitation energies of the lowest-energy singlet excited state of a large series of common BODIPY derivatives employing various functional aiming at the best possible combination providing the least deviations from the experimental values. Using the common "fudge" correction, a series of combinations was investigated, and a methodology is proposed offering equal or better performances than what is reported in the literature.

Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding

Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.

Detection of Chemical Weapon Agents Using Spectroscopic Probes: A Computational Study

Organophosphorus compounds are organic compounds widely employed in agriculture as well as in chemical weapons. The use in agriculture is due to their insecticidal properties. However, in chemical warfare, the use of organophosphorus is associated with acetylcholinesterase inhibition, which promotes the cholinergic syndromes. In this line, the fast detection of this class of compound is crucial for the determination of environmental exposure. This improved detection will naturally allow for more prompt courses of treatment depending on the contaminant findings. In this perspective, the dipyrrinone oxime (1) was employed for the detection of organophosphorus compounds that are employed as nerve agents, such as cyclosarin, sarin, soman, diethyl chlorophosphate, diisopropylfluorophosphate, 2-(dimethylamino)ethyl N,N-dimethylphosphoramidofluoridate, O-ethyl-S-[2-(diethylamino)ethyl]methylphosphonothioate, O-ethyl-S-[2(diisopropylamino)ethyl] methylphosphonothioate, and O,O-diethyl-S-[2-(diethylamino)ethyl] phosphorothioate, through fluorescent emission. The thermodynamics and kinetic parameters as well as spectroscopic properties of the complexes formed for 1 and all organophosphorus compounds previously cited were investigated by means of theoretical calculations. From our findings, only the diethyl chlorophosphate, 2-(dimethylamino)ethyl N,N-dimethylphosphoramidofluoridate, and O,O-diethyl-S-[2-(diethylamino)ethyl] phosphorothioate emitted fluorescence in the hexane, toluene, chloroform, dichloromethane, methanol, acetonitrile, water, and dimethyl sulfoxide solvents. The study of the absorption wavelength with the most polar solvent showed higher values compared to apolar solvents. In the same solvent, for instance, soman in hexane showed the lowest absorption wavelength value, 324.5 nm, and DCP the highest value, 330.8 nm. This behavior was observed in other tested solvents. The thermodynamic parameters indicate negative Gibbs free energy (ΔG) values for the O-ethyl-S-[2(diisopropylamino)ethyl] methylphosphonothioate with 1 reaction. On the other hand, the sarin and cyclosarin revealed the lowest Gibbs free energy (ΔG‡) values, being kinetically favorable and presenting more reactivity.

Ab Initio Prediction of Fluorescence Lifetimes Involving Solvent Environments by Means of COSMO and Vibrational Broadening

The fluorescence lifetime is a key property of fluorophores that can be utilized for microenvironment probing, analyte sensing, and multiplexing as well as barcoding applications. For the rational design of lifetime probes and barcodes, theoretical methods have been developed to enable the ab initio prediction of this parameter, which depends strongly on interactions with solvent molecules and other chemical species in the emitteŕs immediate environment. In this work, we investigate how a conductor-like screening model (COSMO) can account for variations in fluorescence lifetimes that are caused by such fluorophore-solvent interactions. Therefore, we calculate vibrationally broadened fluorescence spectra using the nuclear ensemble method to obtain distorted molecular geometries to sample the electronic transitions with time-dependent density functional theory (TDDFT). The influence of the solvent on fluorescence lifetimes is accounted for with COSMO. For example, for 4-hydroxythiazole fluorophore containing different heteroatoms and acidic and basic moieties in aprotic and protic solvents of varying polarity, this approach was compared to experimentally determined lifetimes in the same solvents. Our results demonstrate a good correlation between theoretically predicted and experimentally measured fluorescence lifetimes except for the polar solvents ethanol and acetonitrile that can specifically interact with the heteroatoms and the carboxylic acid of the thiazole derivative.

Exchange-repulsion energy in QM/EFP

The effective fragment potential (EFP) is a quantum mechanics (QM)-based model designed to accurately describe intermolecular interactions. Hybrid QM/EFP calculations combine quantum mechanical methods with an EFP embedding to study complex systems in which many-body effects are relevant. As in EFP-only calculations, non-bonded interactions between the QM region and EFP fragments are computed as a sum of electrostatic, polarization, dispersion, and exchange-repulsion energies. The exchange-repulsion term is a computational bottleneck of the EFP calculations. Here, we present a general procedure for computing the QM/EFP exchange-repulsion interactions based on one-electron contributions to the QM Hamiltonian, by using Gaussian functions to represent localized molecular orbitals of the effective fragments. The accuracy of the exchange-repulsion and total QM/EFP interaction energies is evaluated on a diverse set of dimers, including complexes from the S22 dataset of non-covalent interactions. In most cases, the QM/EFP energies are at least as accurate as corresponding EFP energies. A simple and computationally efficient form of the introduced QM/EFP exchange-repulsion term will facilitate further developments and applications of QM/EFP methods.

Relation of molecular structure to Franck-Condon bands in the visible-light absorption spectra of symmetric cationic cyanine dyes

A Franck-Condon (FC) model is used to study the solution-phase absorbance spectra of a series of seven symmetric cyanine dyes having between 22 and 77 atoms. Electronic transition energies were obtained from routine visible-light absorbance and fluorescence emission spectra. Harmonic normal modes were computed using density functional theory (DFT) and a polarizable continuum solvent model (PCM), with frequencies corrected using measured mid-infrared spectra. The model predicts the relative energies of the two major vibronic bands to within 5% and 11%, respectively, and also reproduces structure-specific differences in vibronic band shapes. The bands themselves result from excitation of two distinct subsets of normal modes, one with frequencies between 150 and 625cm(-1), and the other between 850 and 1480cm(-1). Vibronic transitions excite symmetric in-plane bending of the polymethine chain, in-plane bends of the polymethine and aromatic C-H bonds, torsions and deformations of N-alkyl substituents, and in the case of the indocyanines, in-plane deformations of the indole rings. For two dyes, the model predicts vibronic coupling into symmetry-breaking torsions associated with trans-cis photoisomerization.