Matteo Ambrosetti

We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.

The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.

We demonstrate the pivotal role of QM density confinement effects on solvatochromic shifts. In particular, by resorting to a QM/MM approach capable of accounting for confinement effects we successfully reproduce vacuo-to-water solvatochromic shifts for dark nπ* and bright ππ* transitions of acrolein and dark nπ* transitions of pyridine and pyrimidine without the need of including explicit water molecules in the QM portion. Remarkably, our approach is also able to dissect the effects of the single forces acting on the solute-solvent couple, and allows for a rationalization of the experimental findings in terms of physico-chemical quantities.

We present the extension of the quantum/classical polarizable fluctuating charge model to the calculation of single residues of quadratic response functions, as required for the computational modeling of two-photon absorption cross-sections. By virtue of a variational formulation of the quantum/classical polarizable coupling, we are able to exploit an atomic orbital-based quasienergy formalism to derive the additional coupling terms in the response equations. Our formalism can be extended to the calculation of arbitrary order response functions and their residues. The approach has been applied to the challenging problem of one- and two-photon spectra of rhodamine 6G (R6G) in aqueous solution. Solvent effects on one- and two-photon spectra of R6G in aqueous solution have been analyzed by considering three different approaches, from a continuum (QM/PCM) to two QM/MM models (non-polarizable QM/TIP3P and polarizable QM/FQ). Both QM/TIP3P and QM/FQ simulated OPA and TPA spectra show that the inclusion of discrete water solvent molecules is essential to increase the agreement between theory and experiment. QM/FQ has been shown to give the best agreement with experiments.

The novel polarizable FQFμ force field is proposed and coupled to a QM SCF Hamiltonian. The peculiarity of the resulting QM/FQFμ approach stands in the fact the polarization effects are modeled in terms of both fluctuating charges and dipoles, which vary as a response to the external electric field/potential. Remarkably, QM/FQFμ is defined in terms of three parameters: electronegativity and chemical hardness, which are well defined in Density Functional Theory, and polarizability, which is a physical observable. Such parameters are numerically adjusted so to reproduce full QM reference electrostatic energy values. The model is challenged against test molecular systems in aqueous solution, showing remarkable accuracy and thus highlighting its potentialities for future extensive applications.

A combination of experimental and quantum mechanical investigations is applied to the study of the optical features of 4,4′-bis(2-benzoxazolyl)stilbene (BBS) dissolved in solution or in a poly(L-lactic acid) (PLA) thermoplastic matrix at different concentrations. The experimental analyses allow the characterization of BBS solutions and dispersions in terms of absorption and emission features, along with the collection of some key parameters such as fluorescence quantum yield, anisotropy and lifetime, while the computational approach gives a detailed description of the photophysical behavior of BBS in the different environments. For the 10⁻⁵ M BBS solution, the fluorescence spectra show the expected peaks at 425 and 455 nm of the non-interacting BBS molecules with a single fluorescence lifetime of 0.85 ns without revealing any aggregation phenomena, prevented by the short lifetime and fast diffusion rate of the monomer. Moreover, the calculated spectra are in excellent agreement with the experiments, thus showing the reliability of the computational approach. In time-resolved emission experiments (TRES) on more concentrated solutions (10⁻⁴ M) and on BBS crystals, the presence of an excimer is revealed by the appearance of a broad peak around 540 nm, followed by the disappearance of the two main peaks at 460 nm on a time scale of about 10 ns. The computational analysis attributes this behavior to the formation of aggregates of different geometries. The BBS dispersions in PLA reveal the presence of different BBS architectures depending on the fluorophore content. Even at low concentrations, BBS is mainly dispersed as a monomer in the matrix, spheroid aggregates of about 800–900 nm in diameter are also present and the relevant fluorescence spectra arise from the combination of monomer and aggregate contributions. At higher concentrations, BBS starts forming crystals of a peculiar helicoidal shape, with a diameter of about 2 μm, variable length up to several hundreds of μm and emission spectra similar to those of isolated BBS crystals.

We present a computational study on the spectroscopic properties of UV‐Vis absorbing dyes in water solution. We model the solvation environment by using both continuum and discrete models, with and without polarization, to establish how the physical and chemical properties of the solute–solvent interaction may affect the spectroscopic response of aqueous systems. Seven different compounds were chosen, representing different classes of organic molecules. The classical atomistic description of the solvent molecules was enriched with polarization effects treated by means of the fluctuating charges (FQ) model, propagated to the first‐order response function of the quantum‐mechanical (QM) solute to include its effects withing the modeling of the electronic excitations of the systems. Results obtained with the QM/FQ model were compared with those from continuum solvation models as well as nonpolarizable atomistic models, and then confronted with the experimental values to determine the accuracy that can be expected with each level of theory. Moreover, a thorough structural analysis using molecular dynamics simulations is provided for each system.