Customising excitation properties of polycyclic aromatic hydrocarbons by rational positional heteroatom doping: the peri-xanthenoxanthene (PXX) case.

In this paper we tackle the challenge of gaining control of the photophysical properties of PAHs through a site-specific N-doping within the structural aromatic framework. By developing a simple predictive tool that identifies C(sp2)-positions that if substituted with a heteroatom would tailor the changes in the absorption and emission spectral envelopes, we predict optimal substitutional patterns for the model peri-xanthenoxanthene (PXX) PAH. Specifically, TDDFT calculations of the electron density difference between the S1 excited state and S0 ground state of PXX allowed us to identify the subtleties in the role of sites i.e., electron donating or withdrawing character on excitation. The replacement of two C(sp2)-atoms with two N-atoms, in either electron donating or withdrawing positions, shifts the electronic transitions either to low or high energy, respectively. This consequently shifts the PXX absorption spectral envelop bathochromically or hypsochromically, as demonstrated by steady-state absorption spectroscopic measurements. Within the series of synthesised N-doped PXX, we tune the optical band gap within an interval of ∼0.4 eV, in full agreement with the theoretical predictions. Relatedly, measurements show the more blueshifted the absorption/emission energies, the greater the fluorescence quantum yield value (from ∼45% to ∼75%). On the other hand, electrochemical investigations suggested that the N-pattern has a limited influence on the redox properties. Lastly, depending on the N-pattern, different supramolecular organisations could be obtained at the solid-state, with the 1,7-pattern PXX molecule forming multi-layered, graphene-like, supramolecular sheets through a combination of weak H-bonding and π-π stacking interactions. Supramolecular striped patterned sheets could also be formed with the 3,9- and 4,10-congeners when co-crystallized with a halogen-bond donor molecule.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue.

The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S1 and dark S2 surfaces near the S1 minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S1/S2 state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue.

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Laurdan and Di-4-ANEPPDHQ Influence the Properties of Lipid Membranes: A Classical Molecular Dynamics and Fluorescence Study.

Environmentally sensitive (ES) dyes have been used for many decades to study the lipid order of cell membranes, as different lipid phases play a crucial role in a wide variety of cell processes. Yet, the understanding of how ES dyes behave, interact, and affect membranes at the atomistic scale is lacking, partially due to the lack of molecular dynamics (MD) models of these dyes. Here, we present ground- and excited-state MD models of commonly used ES dyes, Laurdan and di-4-ANEPPDHQ, and use MD simulations to study the behavior of these dyes in a disordered and an ordered membrane. We also investigate the effect that these two dyes have on the hydration and lipid order of the membranes, where we see a significant effect on the hydration of lipids proximal to the dyes. These findings are combined with experimental fluorescence experiments of ordered and disordered vesicles and live HeLa cells stained by the aforementioned dyes, where the generalized polarization (GP) values were measured at different concentrations of the dyes. We observe a small but significant decrease of GP at higher Laurdan concentrations in vesicles, while the same effect is not observed in cell membranes. The opposite effect is observed with di-4-ANEPPDHQ where no significant change in GP is seen for vesicles but a very substantial and significant decrease is seen in cell membranes. Together, our results show the profound effect that ES dyes have on membranes, and the presented MD models will be important for further understanding of these effects.