1. Multi-replica biased sampling for photoswitchable π-conjugated polymers
- Author
-
Mariagrazia Fortino, Concetta Cozza, Adriana Pietropaolo, Massimiliano Bonomi, Università degli Studi 'Magna Graecia' di Catanzaro [Catanzaro, Italie] (UMG), Bioinformatique structurale - Structural Bioinformatics, Centre National de la Recherche Scientifique (CNRS)-Institut Pasteur [Paris], Università degli Studi 'Magna Graecia' di Catanzaro = University of Catanzaro (UMG), and Institut Pasteur [Paris] (IP)-Centre National de la Recherche Scientifique (CNRS)
- Subjects
Chemical Physics (physics.chem-ph) ,Physics ,010304 chemical physics ,Metadynamics ,FOS: Physical sciences ,General Physics and Astronomy ,Dihedral angle ,010402 general chemistry ,01 natural sciences ,Molecular physics ,0104 chemical sciences ,[CHIM.THEO]Chemical Sciences/Theoretical and/or physical chemistry ,symbols.namesake ,Physics - Chemical Physics ,Excited state ,0103 physical sciences ,symbols ,Density functional theory ,Emission spectrum ,Physical and Theoretical Chemistry ,Perturbation theory ,Ground state ,Hamiltonian (quantum mechanics) - Abstract
In recent years, pi-conjugated polymers are attracting considerable interest in view of their light-dependent torsional reorganization around the pi-conjugated backbone, which determines peculiar light-emitting properties. Motivated by the interest in designing conjugated polymers with tunable photoswitchable pathways, we devised a computational framework to enhance the sampling of the torsional conformational space and at the same time estimate ground to excited-state free-energy differences. This scheme is based on a combination of Hamiltonian Replica Exchange (REM), Parallel Bias metadynamics, and free-energy perturbation theory. In our scheme, each REM replica samples an intermediate unphysical state between the ground and the first two excited states, which are characterized by TD-DFT simulations at the B3LYP/6-31G* level of theory. We applied the method to a 5-mer of 9,9-dioctylfluorene and found that upon irradiation this system can undergo a dihedral inversion from 155 to -155 degrees crossing a barrier that decreases from 0.1 eV in the ground state (S0) to 0.05 eV and 0.04 eV in the first (S1) and second (S2) excited states. Furthermore, S1 and even more S2 were predicted to stabilize coplanar dihedrals, with a local free-energy minimum located at +-44 degrees. The presence of a free-energy barrier of 0.08 eV for the S1 and 0.12 eV for the S2 state can trap this conformation in a basin far from the global free-energy minimum located at 155 degrees. The simulation results were compared with the experimental emission spectrum, showing a quantitative agreement with the predictions provided by our framework., Comment: 14 pages, 6 figures. The present article has been accepted in The Journal of Chemical Physics
- Published
- 2021