Simulation workflows to predict the circular dichroism and circularly polarized luminescence of chiral materials

Chiral materials are attracting considerable interest in various fields in view of their unique properties and optical activity. Indeed, the peculiar features of chiral materials to absorb and emit circularly polarized light enable their use in an extensive range of applications. Motivated by the interest in boosting the development of chiral materials characterized by enhanced chiroptical properties such as circular dichroism (CD) and circular polarized luminescence (CPL), we herein illustrate in this tutorial how theoretical simulations can be used for the predictions and interpretations of chiroptical data and for the identification of chiral geometries. We are focusing on computational frameworks that can be used to investigate the theoretical aspects of chiral materials' photophysical and conformational characteristics. We will then illustrate ab initio methods based on density functional theory (DFT) and its time-dependent extension (TD-DFT) to simulate CD and CPL signals, and we will exemplify a variety of enhanced sampling techniques useful for an adequate sampling of the configurational space for chiral systems.

Biocompatible and Photostable Photoacoustic Contrast Agents as Nanoparticles Based on Bodipy Scaffold and Polylactide Polymers: Synthesis, Formulation, and In Vivo Evaluation

We have designed a new Bodipy scaffold for efficient in vivo photoacoustic (PA) imaging of nanoparticles commonly used as drug nanovectors. The new dye has an optimized absorption band in the near-infrared window in biological tissue and a low fluorescence quantum yield that leads to a good photoacoustic generation efficiency. After Bodipy-initiated ring-opening polymerization of lactide, the polylactide-Bodipy was formulated into PEGylated nanoparticles (NPs) by mixing with PLA-PEG at different concentrations. Formulated NPs around 100 nm exhibit excellent PA properties: an absorption band at 760 nm and a molar absorption coefficient in between that of molecular PA absorbers and gold NPs. Highly improved photostability compared to cyanine-labeled PLA NPs as well as innocuity in cultured macrophages were demonstrated. After intravenous injection in healthy animals, NPs were easily detected using a commercial PA imaging system and spectral unmixing, opening the way to their use as theranostic agents.

The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy

The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.

Supramolecular BODIPY based dimers: synthesis, computational and spectroscopic studies

The synthetic procedures for the preparation of supramolecular BODIPY dimers decorated with complementary patterns able to induce the formation of a triple hydrogen bond through mutual interactions are here reported. The BODIPY and styryl-equipped BODIPY species have been suitably functionalized in meso position with 2,6-diacetamido-4-pyridyl and 1-butyl-6-uracyl moieties. Dimers and monomers have been subjected to computational and photophysical investigations in solvent media. Various peculiarities concerning the effects of the interaction geometry on the stability of the H-bonded systems have also been investigated. The combination of modelling and experimental data provides a paradigm for improving and refining the BODIPY synthetic pathway to have chromophoric architectures with a programmable supramolecular identity. Furthermore, the possibility of assembling dimers of different dyes through H-bonds could be appealing for a systematic investigation of the principal factors affecting the dynamics of the energy migration and possibly driving coherent transfer mechanisms. Our work highlights how the chemical versatility of these dyes can be exploited to design new BODIPY-based supramolecular architectures.

Ultrafast Fluctuations in PM6 Domains of Binary and Ternary Organic Photovoltaic Thin Films Probed with Two-Dimensional White-Light Spectroscopy

We present two-dimensional white-light spectroscopy (2DWL) measurements of binary and ternary bulk heterojunctions of the polymer donor PM6 mixed with state-of-the-art nonfullerene acceptors Y6 or IT4F. The ternary film has a shorter lifetime and faster spectral diffusion than either of the binary films. 2D line shape analysis of the PM6 ground state bleach with a Kubo model determines that all three films have similar amplitudes of fluctuations (Δ = 0.29 fs^-1) in their transition frequencies, but different relaxation times (ranging from 102 to 24 fs). The ternary film exhibits faster dynamics than either of the binary films. The short lifetime of the ternary blend is consistent with increased photoexcitation transfer and the fast frequency fluctuations are consistent with structural dynamics of aliphatic side chains. These results suggest that the femtosecond fluctuations of PM6 are impacted by the choice of the acceptor molecules. We hypothesize that those dynamics are either indicative, or perhaps the initial source, of structural dynamics that ultimately contribute to solar cell operation.

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Q_x and Q_y states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Q_x-Q_y crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

The effect of hydrogen bonds on the ultrafast relaxation dynamics of a BODIPY dimer

The influence of hydrogen bonds (H-bonds) in the structure, dynamics, and functionality of biological and artificial complex systems is the subject of intense investigation. In this broad context, particular attention has recently been focused on the ultrafast H-bond dependent dynamical properties in the electronic excited state because of their potentially dramatic consequences on the mechanism, dynamics, and efficiency of photochemical reactions and photophysical processes of crucial importance for life and technology. Excited-state H-bond dynamics generally occur on ultrafast time scales of hundreds of femtoseconds or less, making the characterization of associated mechanisms particularly challenging with conventional time-resolved techniques. Here, 2D electronic spectroscopy is exploited to shed light on this still largely unexplored dynamic mechanism. An H-bonded molecular dimer prepared by self-assembly of two boron-dipyrromethene dyes has been specifically designed and synthesized for this aim. The obtained results confirm that upon formation of H-bonds and the dimer, a new ultrafast relaxation channel is activated in the ultrafast dynamics, mediated by the vibrational motions of the hydrogen donor and acceptor groups. This relaxation channel also involves, beyond intra-molecular relaxations, an inter-molecular transfer process. This is particularly significant considering the long distance between the centers of mass of the two molecules. These findings suggest that the design of H-bonded structures is a particularly powerful tool to drive the ultrafast dynamics in complex materials.

Triangulenium dyes: the comprehensive photo-absorption and emission story of a versatile family of chromophores

The triangulenium dyes constitute a family of versatile chromophores whose impressive photo-absorption and emission properties are currently highlighted in numerous novel experimental applications. In this investigation, we provide a comprehensive TDDFT characterization of their spectroscopic properties elucidating the origin of their large and complex absorption and emission vibronic spectra spread over the (whole) visible region. More precisely, by benchmarking the performance of 10 commonly-used exchange-correlation density functionals belonging to different classes of approximation, we develop and validate a computational protocol allowing the accurate modeling of both the position and optical line-shape of their vibrationally-resolved absorption and emission band structures. We find that semilocal approximations provide the best estimate of the structure of the vibronic spectra, however they spuriously and strongly underestimate their position. We finally show that global-hybrid density functionals mixing between 20 and 30% of exact-like exchange are an excellent compromise to get a satisfactory estimate of both of these properties.

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.

Sonochemical synthesis of microscale Zn(ii)-MOF with dual Lewis basic sites for fluorescent turn-on detection of Al3+ and methanol with low detection limits

A microscale Zn-MOF containing dual Lewis basic sites for the efficient bifunctional fluorescent detection of Al³⁺ and methanol.