RESUMEN
Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.
Asunto(s)
Simulación de Dinámica Molecular , Péptidos , Fluorescencia , Espectrometría de FluorescenciaRESUMEN
Photoelectron angular distribution (PAD) in the laboratory frame for randomly oriented molecules is typically described by a single anisotropy parameter, the so-called asymmetry parameter. However, especially from a theoretical perspective, it is more natural to consider molecular photoionization by using a molecular frame. The molecular frame PADs (MFPADs) may be used to extract information about the electronic structure of the system studied. In the last decade, significant experimental efforts have been directed to MFPAD measurements. MFPADs are highly characterizing signatures of the final ionic states. In particular, they are very sensitive to the nature of the final state, which is embodied in the corresponding Dyson orbital. In our previous work on acetylacetone, a prototype system for studying intra-molecular hydrogen bond interactions, we followed the dynamics of the excited states involved in the photoexcitation-deexcitation process of this molecule. It remains to be explored the possibility of discriminating between different excited states through the MFPAD profiles. The calculation of MFPADs to differentiate excited states can pave the way to the possibility of a clear discrimination for all the cases where the recognition of excited states is otherwise intricate.
RESUMEN
Two algorithms for calculating overlaps between CIS (or TDDFT) type excited state wave functions are presented, one based on an expansion of overlap determinants into level 2 minors (OL2M) and the other based on an expansion of the wave functions into natural transition orbitals (ONTO). Both algorithms are significantly faster than previously available algorithms, with the ONTO algorithm reducing the cost of a single overlap element calculation by a factor of the square of the number of occupied orbitals in the system. The algorithm exhibits orders of magnitude faster calculations for large systems and significantly increases the size of systems for which TDDFT based nonadiabatic dynamics simulations can be performed. The OL2M algorithm is substantially slower for a single overlap matrix element but becomes preferred when overlaps between large numbers of states are required. Additionally, we test the accuracy of approximate overlaps calculated using truncated wave functions and show that truncation can lead to large errors in the overlaps. Lastly, we provide examples of applications for wave function overlaps outside the context of nonadiabatic dynamics.
RESUMEN
A set of mass-weighted internal coordinates was derived and applied to the double proton transfer reaction in the formic acid dimer (FAD). The coordinate set was obtained starting from the Hirschfelder "mobile" by an optimization procedure consisting of a sequence of kinematic rotations. In FAD, the optimization procedure leads to three coordinates that do change significantly along the reaction path. These coordinates span the reaction space, whereas the remaining modes are treated in a harmonic approximation. The effect that the dimer dissociative motion has on the ground and excited vibrational states dynamics was explored. In the frequency region corresponding to the symmetric OH-stretch vibration four doublets have been identified with splittings of 2.76, 0.07, 0.60, and 4.03 cm(-1).