Sebastian Fernandez-Alberti

Nonadiabatic Excited-State Molecular Dynamics Modeling of Photoinduced Dynamics in Conjugated Molecules

Nonadiabatic dynamics generally defines the entire evolution of electronic excitations in optically active molecular materials. It is commonly associated with a number of fundamental and complex processes such as intraband relaxation, energy transfer, and light harvesting influenced by the spatial evolution of excitations and transformation of photoexcitation energy into electrical energy via charge separation (e.g., charge injection at interfaces). To treat ultrafast excited-state dynamics and exciton/charge transport we have developed a nonadiabatic excited-state molecular dynamics (NA-ESMD) framework incorporating quantum transitions. Our calculations rely on the use of the Collective Electronic Oscillator (CEO) package accounting for many-body effects and actual potential energy surfaces of the excited states combined with Tullys fewest switches algorithm for surface hopping for probing nonadiabatic processes. This method is applied to model the photoinduced dynamics of distyrylbenzene (a small oligomer of polyphenylene vinylene, PPV). Our analysis shows intricate details of photoinduced vibronic relaxation and identifies specific slow and fast nuclear motions that are strongly coupled to the electronic degrees of freedom, namely, torsion and bond length alternation, respectively. Nonadiabatic relaxation of the highly excited mA(g) state is predicted to occur on a femtosecond time scale at room temperature and on a picosecond time scale at low temperature.

Evolutionary conservation of protein vibrational dynamics

The aim of the present work is to study the evolutionary divergence of vibrational protein dynamics. To this end, we used the Gaussian Network Model to perform a systematic analysis of normal mode conservation on a large dataset of proteins classified into homologous sets of family pairs and superfamily pairs. We found that the lowest most collective normal modes are the most conserved ones. More precisely, there is, on average, a linear correlation between normal mode conservation and mode collectivity. These results imply that the previously observed conservation of backbone flexibility (B-factor) profiles is due to the conservation of the most collective modes, which contribute the most to such profiles. We discuss the possible roles of normal mode robustness and natural selection in the determination of the observed behavior. Finally, we draw some practical implications for dynamics-based protein alignment and classification and discuss possible caveats of the present approach.

Self-Trapping of Excitons, Violation of Condon Approximation, and Efficient Fluorescence in Conjugated Cycloparaphenylenes

Cycloparaphenylenes, the simplest structural unit of armchair carbon nanotubes, have unique optoelectronic properties counterintuitive in the class of conjugated organic materials. Our time-dependent density functional theory study and excited state dynamics simulations of cycloparaphenylene chromophores provide a simple and conceptually appealing physical picture explaining experimentally observed trends in optical properties in this family of molecules. Fully delocalized degenerate second and third excitonic states define linear absorption spectra. Self-trapping of the lowest excitonic state due to electron-phonon coupling leads to the formation of spatially localized excitation in large cycloparaphenylenes within 100 fs. This invalidates the commonly used Condon approximation and breaks optical selection rules, making these materials superior fluorophores. This process does not occur in the small molecules, which remain inefficient emitters. A complex interplay of symmetry, π-conjugation, conformational distortion and bending strain controls all photophysics of cycloparaphenylenes.

Evolutionary Conservation of Protein Backbone Flexibility

Internal protein dynamics is essential for biological function. During evolution, protein divergence is functionally constrained: properties more relevant for function vary more slowly than less important properties. Thus, if protein dynamics is relevant for function, it should be evolutionary conserved. In contrast with the well-studied evolution of protein structure, the evolutionary divergence of protein dynamics has not been addressed systematically before, apart from a few case studies. X-Ray diffraction analysis gives information not only on protein structure but also on B-factors, which characterize the flexibility that results from protein dynamics. Here we study the evolutionary divergence of protein backbone dynamics by comparing the Cα flexibility (B-factor) profiles for a large dataset of homologous proteins classified into families and superfamilies. We show that Cα flexibility profiles diverge slowly, so that they are conserved at family and superfamily levels, even for pairs of proteins with nonsignificant sequence similarity. We also analyze and discuss the correlations among the divergences of flexibility, sequence, and structure.

Nonadiabatic excited-state molecular dynamics: Numerical tests of convergence and parameters

Nonadiabatic molecular dynamics simulations, involving multiple Born-Oppenheimer potential energy surfaces, often require a large number of independent trajectories in order to achieve the desired convergence of the results, and simulation relies on different parameters that should be tested and compared. In addition to influencing the speed of the simulation, the chosen parameters combined with the frequently reduced number of trajectories can sometimes lead to unanticipated changes in the accuracy of the simulated dynamics. We have previously developed a nonadiabatic excited state molecular dynamics methodology employing Tullys fewest switches surface hopping algorithm. In this study, we seek to investigate the impact of the number of trajectories and the various parameters on the simulation of the photoinduced dynamics of distyrylbenzene (a small oligomer of polyphenylene vinylene) within our developed framework. Various user-defined parameters are analyzed: classical and quantum integration time steps, the value of the friction coefficient for Langevin dynamics, and the initial seed used for stochastic thermostat and hopping algorithms. Common approximations such as reduced number of nonadiabatic coupling terms and the classical path approximation are also investigated. Our analysis shows that, at least for the considered molecular system, a minimum of ~400 independent trajectories should be calculated in order to achieve statistical averaging necessary for convergence of the calculated relaxation timescales.

Nonadiabatic excited-state molecular dynamics: Treatment of electronic decoherence

Within the fewest switches surface hopping (FSSH) formulation, a swarm of independent trajectories is propagated and the equations of motion for the quantum coefficients are evolved coherently along each independent nuclear trajectory. That is, the phase factors, or quantum amplitudes, are retained. At a region of strong coupling, a trajectory can branch into multiple wavepackets. Directly following a hop, the two wavepackets remain in a region of nonadiabatic coupling and continue exchanging population. After these wavepackets have sufficiently separated in phase space, they should begin to evolve independently from one another, the process known as decoherence. Decoherence is not accounted for in the standard surface hopping algorithm and leads to internal inconsistency. FSSH is designed to ensure that at any time, the fraction of classical trajectories evolving on each quantum state is equal to the average quantum probability for that state. However, in many systems this internal consistency requirement is violated. Treating decoherence is an inherent problem that can be addressed by implementing some form of decoherence correction to the standard FSSH algorithm. In this study, we have implemented two forms of the instantaneous decoherence procedure where coefficients are reinitialized following hops. We also test the energy-based decoherence correction (EDC) scheme proposed by Granucci et al. and a related version where the form of the decoherence time is taken from Truhlars Coherent Switching with Decay of Mixing method. The sensitivity of the EDC results to changes in parameters is also evaluated. The application of these computationally inexpensive ad hoc methods is demonstrated in the simulation of nonradiative relaxation in two conjugated oligomer systems, specifically poly-phenylene vinylene and poly-phenylene ethynylene. We find that methods that have been used successfully for treating small systems do not necessarily translate to large polyatomic systems and their success depends on the particular system under study.

Instantaneous normal modes, resonances, and decay channels in the vibrational relaxation of the amide I mode of N-methylacetamide-D in liquid deuterated water

A nonequilibrium molecular dynamics (MD) study of the vibrational relaxation of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution is carried out using instantaneous normal modes (INMs). The identification of the INMs as they evolve over time, which is necessary to analyze the energy fluxes, is made by using a novel algorithm which allows us to assign unequivocally each INM to an individual equilibrium normal mode (ENM) or to a group of ENMs during the MD simulations. The time evolution of the energy stored in each INM is monitored and the occurrence of resonances during the relaxation process is then investigated. The decay of the amide I mode, initially excited with one vibrational quantum, is confirmed to fit well to a biexponential function, implying that the relaxation process involves at least two mechanisms with different rate constants. By freezing the internal motions of the solvent, it is shown that the intermolecular vibration-vibration channel to the bending modes of the solvent is closed. The INM analysis reveals then the existence of a major and faster decay channel, which corresponds to an intramolecular vibrational redistribution process and a minor, and slower, decay channel which involves the participation of the librational motions of the solvent. The faster relaxation pathway can be rationalized in turn using a sequential kinetic mechanism of the type P-->M+L-->L, where P (parent) is the initially excited amide I mode, and M (medium) and L (low) are specific midrange and lower-frequency NMAD vibrational modes, respectively.

Shishiodoshi unidirectional energy transfer mechanism in phenylene ethynylene dendrimers.

Non-adiabatic excited-state molecular dynamics is used to study the ultrafast intramolecular energy transfer between two-, three-, and four-ring linear polyphenylene ethynylene chromophore units linked through meta-substitutions. Twenty excited-state electronic energies, with their corresponding gradients and nonadiabatic coupling vectors were included in the simulations. The initial laser excitation creates an exciton delocalized between the different absorbing two-ring linear PPE units. Thereafter, we observe an ultrafast directional change in the spatial localization of the transient electronic transition density. The analysis of the intramolecular flux of the transition density shows a sequential through-bond two-ring→three-ring→four-ring transfer as well as an effective through-space direct two-to-four ring transfer. The vibrational excitations of C≡C stretching motions change according to that. Finally, a mechanism of unidirectional energy transfer is presented based on the variation of the energy gaps between consecutive electronic excited states in response to the intramolecular flux of the transition density. The mechanism resembles a Shishiodoshi Japanese bamboo water fountain where, once the electronic population has been transferred to the state directly below in energy, the two states decouple thereby preventing energy transfer in the opposite direction.

Analysis of state-specific vibrations coupled to the unidirectional energy transfer in conjugated dendrimers.

The nonadiabatic excited-state molecular dynamics (NA-ESMD) method and excited-state instantaneous normal modes (ES-INMs) analyses have been applied to describe the state-specific vibrations that participate in the unidirectional energy transfer between the coupled chromophores in a branched dendrimeric molecule. Our molecule is composed of two-, three-, and four-ring linear poly(phenyleneethynylene) (PPE) units linked through meta-substitutions. After an initial laser excitation, an ultrafast sequential S(3) → S(2) → S(1) electronic energy transfer from the shortest to longest segment takes place. During each S(n) → S(n-1) (n = 3, 2) transition, ES-INM(S(n)) and ES-INM(S(n-1)) analyses have been performed on S(n) and S(n-1) states, respectively. Our results reveal a unique vibrational mode localized on the S(n) state that significantly matches with the corresponding nonadiabatic coupling vector d(n,(n-1)). This mode also corresponds to the highest frequency ES-INM(S(n)) and it is seen mainly during the electronic transitions. Furthermore, its absence as a unique ES-INM(S(n-1)) reveals that state-specific vibrations play the main role in the efficiency of the unidirectional S(n) → S(n-1) electronic and vibrational energy funneling in light-harvesting dendrimers.

Protein Conformational Diversity Modulates Sequence Divergence

It is well established that the conservation of protein structure during evolution constrains sequence divergence. The conservation of certain physicochemical environments to preserve protein folds and then the biological function originates a site-specific structurally constrained substitution pattern. However, protein native structure is not unique. It is known that the native state is better described by an ensemble of conformers in a dynamic equilibrium. In this work, we studied the influence of conformational diversity in sequence divergence and protein evolution. For this purpose, we derived a set of 900 proteins with different degrees of conformational diversity from the PCDB database, a conformer database. With the aid of a structurally constrained protein evolutionary model, we explored the influence of the different conformations on sequence divergence. We found that the presence of conformational diversity strongly modulates the substitution pattern. Although the conformers share several of the structurally constrained sites, 30% of them are conformer specific. Also, we found that in 76% of the proteins studied, a single conformer outperforms the others in the prediction of sequence divergence. It is interesting to note that this conformer is usually the one that binds ligands participating in the biological function of the protein. The existence of a conformer-specific site-substitution pattern indicates that conformational diversity could play a central role in modulating protein evolution. Furthermore, our findings suggest that new evolutionary models and bioinformatics tools should be developed taking into account this substitution bias.