Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study

We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.

New tetralactam hosts for squaraine dyes

The photophysical properties of a deep-red fluorescent squaraine dye can be improved by encapsulating it within a tetralactam macrocycle. Three new tetralactams are described with different substituents (methyl, methoxy, methylenedioxy) on the macrocycle aromatic sidewalls. The capability of each tetralactam to encapsulate a squaraine dye in chloroform solution was determined experimentally using absorption, fluorescence, and NMR spectroscopy. Two of the tetralactams were found to thread a squaraine dye with association constants on the order of 106 M-1, while a third macrocycle exhibited no squaraine affinity. An X-ray crystal structure of the third tetralactam showed that the substituents sterically blocked squaraine association. Of the two tetralactams that encapsulate a squaraine, one induces an increase in squaraine fluorescence quantum yield, while the other quenches the squaraine fluorescence. The results suggest that these new squaraine binding systems will be useful for biological imaging and diagnostics applications.

Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules

Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.

A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches

A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.

High Affinity Macrocycle Threading by a Near-Infrared Croconaine Dye with Flanking Polymer Chains

Croconaine dyes have narrow and intense absorption bands at ∼800 nm, very weak fluorescence, and high photostabilities, which combine to make them very attractive chromophores for absorption-based imaging or laser heating technologies. The physical supramolecular properties of croconaine dyes have rarely been investigated, especially in water. This study focuses on a molecular threading process that encapsulates a croconaine dye inside a tetralactam macrocycle in organic or aqueous solvent. Macrocycle association and rate constant data are reported for a series of croconaine structures with different substituents attached to the ends of the dye. The association constants were highest in water (Ka ∼ 10(9) M(-1)), and the threading rate constants (kon) increased in the solvent order H2O > MeOH > CHCl3. Systematic variation of croconaine substituents located just outside the croconaine/macrocycle complexation interface hardly changed Ka but had a strong influence on kon. A croconaine dye with N-propyl groups at each end of the structure exhibited a desirable mixture of macrocycle threading properties; that is, there was rapid and quantitative croconaine/macrocycle complexation at relatively high concentrations in water, and no dissociation of the preassembled complex when it was diluted into a solution of fetal bovine serum, even after laser-induced photothermal heating of the solution. The combination of favorable near-infrared absorption properties and tunable mechanical stability makes threaded croconaine/macrocycle complexes very attractive as molecular probes or as supramolecular composites for various applications in absorption-based imaging or photothermal therapy.

0-0 Energies Using Hybrid Schemes: Benchmarks of TD-DFT, CIS(D), ADC(2), CC2, and BSE/GW formalisms for 80 Real-Life Compounds

The 0-0 energies of 80 medium and large molecules have been computed with a large panel of theoretical formalisms. We have used an approach computationally tractable for large molecules, that is, the structural and vibrational parameters are obtained with TD-DFT, the solvent effects are accounted for with the PCM model, whereas the total and transition energies have been determined with TD-DFT and with five wave function approaches accounting for contributions from double excitations, namely, CIS(D), ADC(2), CC2, SCS-CC2, and SOS-CC2, as well as Green's function based BSE/GW approach. Atomic basis sets including diffuse functions have been systematically applied, and several variations of the PCM have been evaluated. Using solvent corrections obtained with corrected linear-response approach, we found that three schemes, namely, ADC(2), CC2, and BSE/GW allow one to reach a mean absolute deviation smaller than 0.15 eV compared to the measurements, the two former yielding slightly better correlation with experiments than the latter. CIS(D), SCS-CC2, and SOS-CC2 provide significantly larger deviations, though the latter approach delivers highly consistent transition energies. In addition, we show that (i) ADC(2) and CC2 values are extremely close to each other but for systems absorbing at low energies; (ii) the linear-response PCM scheme tends to overestimate solvation effects; and that (iii) the average impact of nonequilibrium correction on 0-0 energies is negligible.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

Accurate simulation of geometry, singlet-singlet and triplet-singlet excitation of cyclometalated iridium(III) complex

In the current contribution, we present a critical study of the theoretical protocol used for the determination of the electronic spectra properties of luminescent cyclometalated iridium(III) complex, [Ir(III)(ppy)₂H₂dcbpy]⁺ (where, ppy = 2-phenylpyridine, H₂dcbpy = 2,2'-bipyridine-4,4'-dicarboxylic acid), considered as a representative example of the various problems related to the prediction of electronic spectra of transition metal complex. The choice of the exchange-correlation functional is crucial for the validity of the conclusions that would be drawn from the numerical results. The influence of the exchange-correlation on geometry parameter and absorption/emission band, the role of solvent effects on time-dependent density function theory (TD-DFT) calculations, as well as the importance of the chosen proper procedure to optimize triplet excited geometry, have been thus examined in detail. From the obtained results, some general conclusions and guidelines are presented: i) PBE0 functional is the most accurate in prediction of ground state geometry; ii) the well-established B3LYP, B3P86, PBE0, and X3LYP have similar accuracy in calculation of absorption spectrum; and iii) the hybrid approach TD-DFT//CIS gives out excellent agreement in the evaluation of triplet excitation energy.

The calculations of excited-state properties with Time-Dependent Density Functional Theory

In this tutorial review, we show how Time-Dependent Density Functional Theory (TD-DFT) has become a popular tool for computing the signatures of electronically excited states, and more specifically, the properties directly related to the optical (absorption and emission) spectra of molecules. We discuss the properties that can be obtained with widely available programs as well as how to account for the environmental effects (solvent and surfaces) and present recent applications in these fields. We next expose the transformation of the TD-DFT results into chemically intuitive parameters (colours as well as charge-transfer distances). Eventually, the non-specialised reader will find a series of advices and warnings necessary to perform her/his first TD-DFT calculations.

Theoretical modeling of large molecular systems. Advances in the local self consistent field method for mixed quantum mechanics/molecular mechanics calculations

Molecular mechanics methods can efficiently compute the macroscopic properties of a large molecular system but cannot represent the electronic changes that occur during a chemical reaction or an electronic transition. Quantum mechanical methods can accurately simulate these processes, but they require considerably greater computational resources. Because electronic changes typically occur in a limited part of the system, such as the solute in a molecular solution or the substrate within the active site of enzymatic reactions, researchers can limit the quantum computation to this part of the system. Researchers take into account the influence of the surroundings by embedding this quantum computation into a calculation of the whole system described at the molecular mechanical level, a strategy known as the mixed quantum mechanics/molecular mechanics (QM/MM) approach. The accuracy of this embedding varies according to the types of interactions included, whether they are purely mechanical or classically electrostatic. This embedding can also introduce the induced polarization of the surroundings. The difficulty in QM/MM calculations comes from the splitting of the system into two parts, which requires severing the chemical bonds that link the quantum mechanical subsystem to the classical subsystem. Typically, researchers replace the quantoclassical atoms, those at the boundary between the subsystems, with a monovalent link atom. For example, researchers might add a hydrogen atom when a C-C bond is cut. This Account describes another approach, the Local Self Consistent Field (LSCF), which was developed in our laboratory. LSCF links the quantum mechanical portion of the molecule to the classical portion using a strictly localized bond orbital extracted from a small model molecule for each bond. In this scenario, the quantoclassical atom has an apparent nuclear charge of +1. To achieve correct bond lengths and force constants, we must take into account the inner shell of the atom: for an sp(3) carbon atom, we consider the two core 1s electrons and treat that carbon as an atom with three electrons. This results in an LSCF+3 model. Similarly, a nitrogen atom with a lone pair of electrons available for conjugation is treated as an atom with five electrons (LSCF+5). This approach is particularly well suited to splitting peptide bonds and other bonds that include carbon or nitrogen atoms. To embed the induced polarization within the calculation, researchers must use a polarizable force field. However, because the parameters of the usual force fields include an average of the induction effects, researchers typically can obtain satisfactory results without explicitly introducing the polarization. When considering electronic transitions, researchers must take into account the changes in the electronic polarization. One approach is to simulate the electronic cloud of the surroundings by a continuum whose dielectric constant is equal to the square of the refractive index. This Electronic Response of the Surroundings (ERS) methodology allows researchers to model the changes in induced polarization easily. We illustrate this approach by modeling the electronic absorption of tryptophan in human serum albumin (HSA).

First-principles modeling of dye-sensitized solar cells: challenges and perspectives

Since dye-sensitized solar cells (DSSCs) appeared as a promising inexpensive alternative to the traditional silicon-based solar cells, DSSCs have attracted a considerable amount of experimental and theoretical interest. In contrast with silicon-based solar cells, DSSCs use different components for the light-harvesting and transport functions, which allow researchers to fine-tune each material and, under ideal conditions, to optimize their overall performance in assembled devices. Because of the variety of elementary components present in these cells and their multiple possible combinations, this task presents experimental challenges. The photoconversion efficiencies obtained up to this point are still low, despite the significant experimental efforts spent in their optimization. The development of a low-cost and efficient computational protocol that could qualitatively (or even quantitatively) identify the promising semiconductors, dyes, and electrolytes, as well as their assembly, could save substantial experimental time and resources. In this Account, we describe our computational approach that allows us to understand and predict the different elementary mechanisms involved in DSSC working principles. We use this computational framework to propose an in silico route for the ab initio design of these materials. Our approach relies on a unique density functional theory (DFT) based model, which allows for an accurate and balanced treatment of electronic and spectroscopic properties in different phases (such as gas, solution, or interfaces) and avoids or minimizes spurious computational effects. Using this tool, we reproduced and predicted the properties of the isolated components of the DSSC assemblies. We accessed the microscopic measurable characteristics of the cells such as the short circuit current (J(sc)) or the open circuit voltage (V(oc)), which define the overall photoconversion efficiency of the cell. The absence of empirical or material-related parameters in our approach should allow for its wide application to the optimization of existing devices or the design of new ones.

Excited-state calculations with TD-DFT: from benchmarks to simulations in complex environments

In this perspective, we present an overview of recent progress on Time-Dependent Density Functional Theory (TD-DFT) with a specific focus on its accuracy and on models able to take into account environmental effects, including complex media. To this end, we first summarise recent benchmarks and define an average TD-DFT accuracy in reproducing excitation energies when a conventional approach is used. Next, coupling of TD-DFT with models able to account for different kinds of interactions between a central chromophore and nearby chemical objects (solvent, organic cage, metal as well as semi-conducting surface) is investigated. Examples of application to excitation properties are presented, allowing to briefly describe several recent computational strategies. In addition, an extension of TD-DFT to describe a phenomenon involving interacting chromophores, e.g. the electronic energy transfer (EET), is presented to illustrate that this methodology can be applied to processes beyond the vertical excitation. This perspective therefore aims to provide to non-specialists a flavour of recent trends in the field of simulations of excited states in "realistic" situations.

Absorption and emission spectra of fluorescent silica nanoparticles from TD-DFT/MM/PCM calculations

A multi-scale computational protocol, which combines Quantum Mechanics and Molecular Mechanics (QM/MM) calculations with the polarisable continuum model (PCM), has been used to study the tetramethylrhodamine isothiocyanate (TRITC) fluorophore, embedded in three different environments, namely in water, on an amorphous silica surface and covalently encapsulated in a silica nanoparticle (C dot). Absorption and emission spectra have been simulated by using TD-B3LYP/PCM calculations, performed on the TRITC ground and excited state geometries, optimized at the QM/MM level. The results are in good agreement with experimental data confirming the caging effect played by the silica shell on the mobility of the TRITC molecule when covalently encapsulated in silica nanoparticles. This could result in a decrease of the nonradiative decay rate and thus an increase of the quantum yield of the molecule.

Theoretical study of new ruthenium-based dyes for dye-sensitized solar cells

Two relevant, recently reported, ruthenium-based complexes to be used as sensitizers in Grätzel photovoltaic cells are theoretically studied. The UV/vis absorption spectra have been computed within the time-dependent density functional theory formalism. The obtained excitation energies are compared with the experimental results, and the nature of the transition is analyzed in terms of the electronic density. A preliminary study on the performance of different functionals against the equation of motion coupled cluster is performed on a smaller model system.

New Class of Hydroxy-Substituted Squaraine Rotaxane

A templated macrocyclization reaction was used to permanently encapsulate a highly fluorescent hydroxy-substituted squaraine dye inside a tetralactam macrocycle. The free squaraine dye is quite rigid due to internal hydrogen bonding and its photophysical properties hardly change upon encapsulation. A combination of X-ray and NMR data show that the surrounding tetralactam macrocycle adopts an unusually rigid chair conformation and does not undergo rapid pirouetting. Because of its large size and conformational rigidity, the macrocycle creates anisotropic NMR shielding zones that extend over the N,N-dibutylamino groups at each end of the squaraine thread. This shielding anisotropy allows hindered aryl-N rotation to be observed by NMR and provides direct experimental evidence that quinoid-like resonance structures are major contributors to the bis(N,N-dialkylaminophenyl)squaraine resonance hybrid.

Modeling an electronic conductor based on natural peptide sequences

In this work we used quantum mechanical calculations, molecular dynamics simulations, and QM/MM methods to examine the formation of a pi-stacking ladder and the existence of charge transfer processes in a tubular nanostructure constructed using protein building blocks. Initially, the conformational properties of beta-3-thienylalanine, a synthetic amino acid with an aromatic side group, were studied using quantum mechanical calculations. Next, this amino acid was used to substitute the natural residues at selected positions of the nanotube. Molecular dynamics simulations showed that a rational design of the targeted replacements enhances the stability of the tubular nanostructure and forms a pi-stacking ladder in the inner core of the tube. Finally, QM/MM calculations on the oxidized nanotube evidenced the delocalization of the pi-electron deficiency between the thienyl side groups of the different synthetic residues. The overall results reflected that the system under study is a potential candidate to be used as a nanowire.

Discovery and early development of squaraine rotaxanes

The chemical and photophysical properties of a fluorescent squaraine dye are greatly enhanced when it is mechanically encapsulated inside a tetralactam macrocycle. This feature article describes the synthesis, structure, and photophysical performance of first-generation squaraine rotaxanes, and shows how they can be used as fluorescent imaging probes and chemosensors.

Squaraine rotaxanes with boat conformation macrocycles

Mechanical encapsulation of fluorescent, deep-red bis(anilino)squaraine dyes inside Leigh-type tetralactam macrocycles produces interlocked squaraine rotaxanes. The surrounding macrocycles are flexible and undergo rapid exchange of chair and boat conformations in solution. A series of X-ray crystal structures show how the rotaxane co-conformational exchange process involves simultaneous lateral oscillation of the macrocycle about the center of the encapsulated squaraine thread. Rotaxane macrocycles with 1,4-phenylene sidewalls and 2,6-pyridine dicarboxamide bridging units are more likely to adopt boat conformations in the solid state than analogous squaraine rotaxane systems with isophthalamide-containing macrocycles. A truncated squaraine dye, with a secondary amine attached directly to the central C(4)O(2) core, is less electrophilic than the extended bis(anilino)squaraine analogue, but it is still susceptible to chemical and photochemical bleaching. Its stability is greatly enhanced when it is encapsulated as an interlocked squaraine rotaxane. An X-ray crystal structure of this truncated squaraine rotaxane shows the macrocycle in a boat conformation, and NMR studies indicate that the boat is maintained in solution. Encapsulation as a rotaxane increases the dye's brightness by a factor of 6. The encapsulation process appears to constrain the dye and reduce deformation of the chromophore from planarity. This study shows how mechanical encapsulation as a rotaxane can be used as a rational design parameter to fine-tune the chemical and photochemical properties of squaraine dyes.