More Like this

1. Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory

   https://doi.org/10.1021/acs.jctc.4c00737  

   Dickinson, Joseph A; Hammes-Schiffer, Sharon (September 2024, Journal of Chemical Theory and Computation)

   Proton transfer and hydrogen tunneling play key roles in many processes of chemical and biological importance. The generalized nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method was developed in order to capture hydrogen tunneling effects in systems involving the transfer and tunneling of one or more protons. The generalized NEO-MSDFT method treats the transferring protons quantum mechanically on the same level as the electrons and obtains the delocalized vibronic states associated with hydrogen tunneling by mixing localized NEO-DFT states in a nonorthogonal configuration interaction scheme. Herein, we present the derivation and implementation of analytical gradients for the generalized NEO-MSDFT vibronic state energies and the nonadiabatic coupling vectors between these vibronic states. We use this methodology to perform adiabatic and nonadiabatic dynamics simulations of the double proton transfer reactions in the formic acid dimer and the heterodimer of formamidine and formic acid. The generalized NEO-MSDFT method is shown to capture the strongly coupled synchronous or asynchronous tunneling of the two protons in these processes. Inclusion of vibronically nonadiabatic effects is found to significantly impact the double proton transfer dynamics. This work lays the foundation for a variety of nonadiabatic dynamics simulations of multiple proton transfer systems, such as proton relays and hydrogen-bonding networks. 

   more » « less
2. Organic single crystals of charge-transfer complexes: model systems for the study of donor/acceptor interactions

   https://doi.org/10.1039/D1MH01214B  

   Goetz, Katelyn P.; Iqbal, Hamna F.; Bittle, Emily G.; Hacker, Christina A.; Pookpanratana, Sujitra; Jurchescu, Oana D. (January 2022, Materials Horizons)

   The charge-transfer (CT) state arising as a hybrid electronic state at the interface between charge donor and charge acceptor molecular units is important to a wide variety of physical processes in organic semiconductor devices. The exact nature of this state depends heavily on the nature and co-facial overlap between the donor and acceptor; however, altering this overlap is usually accompanied by extensive confounding variations in properties due to extrinsic factors, such as microstructure. As a consequence, establishing reliable relationships between donor/acceptor molecular structures, their molecular overlap, degree of charge transfer and physical properties, is challenging. Herein, we examine the electronic structure of a polymorphic system based on the donor dibenzotetrathiafulvalene (DBTTF) and the acceptor 7,7,8,8-tetracyanoquinodimethane (TCNQ) in the form of high-quality single crystals varying in the donor–acceptor overlap. Using angle-resolved photoemission spectroscopy, we resolve the highest occupied molecular orbital states of the CT crystals. Analysis based on field-effect transistors allows us to probe the sub-gap states impacting hole and electron transport. Our results expand the understanding on the impact of donor and acceptor interactions on electronic structure and charge transport. 

   more » « less
3. Vibronic structure of photosynthetic pigments probed by polarized two-dimensional electronic spectroscopy and ab initio calculations

   https://doi.org/10.1039/C9SC02329A  

   Song, Yin; Schubert, Alexander; Maret, Elizabeth; Burdick, Ryan K.; Dunietz, Barry D.; Geva, Eitan; Ogilvie, Jennifer P. (September 2019, Chemical Science)

   Bacteriochlorophyll a (Bchl a) and chlorophyll a (Chl a) play important roles as light absorbers in photosynthetic antennae and participate in the initial charge-separation steps in photosynthetic reaction centers. Despite decades of study, questions remain about the interplay of electronic and vibrational states within the Q-band and its effect on the photoexcited dynamics. Here we report results of polarized two-dimensional electronic spectroscopic measurements, performed on penta-coordinated Bchl a and Chl a and their interpretation based on state-of-the-art time-dependent density functional theory calculations and vibrational mode analysis for spectral shapes. We find that the Q-band of Bchl a is comprised of two independent bands, that are assigned following the Gouterman model to Q x and Q y states with orthogonal transition dipole moments. However, we measure the angle to be ∼75°, a finding that is confirmed by ab initio calculations. The internal conversion rate constant from Q x to Q y is found to be 11 ps −1 . Unlike Bchl a, the Q-band of Chl a contains three distinct peaks with different polarizations. Ab initio calculations trace these features back to a spectral overlap between two electronic transitions and their vibrational replicas. The smaller energy gap and the mixing of vibronic states result in faster internal conversion rate constants of 38–50 ps −1 . We analyze the spectra of penta-coordinated Bchl a and Chl a to highlight the interplay between low-lying vibronic states and their relationship to photoinduced relaxation. Our findings shed new light on the photoexcited dynamics in photosynthetic systems where these chromophores are primary pigments. 

   more » « less
4. Analytical gradients for nuclear–electronic orbital multistate density functional theory: Geometry optimizations and reaction paths

   https://doi.org/10.1063/5.0085344  

   Yu, Qi; Schneider, Patrick E.; Hammes-Schiffer, Sharon (March 2022, The Journal of Chemical Physics)

   Hydrogen tunneling plays a critical role in many biologically and chemically important processes. The nuclear–electronic orbital multistate density functional theory (NEO-MSDFT) method was developed to describe hydrogen transfer systems. In this approach, the transferring proton is treated quantum mechanically on the same level as the electrons within multicomponent DFT, and a nonorthogonal configuration interaction scheme is used to produce delocalized vibronic states from localized vibronic states. The NEO-MSDFT method has been shown to provide accurate hydrogen tunneling splittings for fixed molecular systems. Herein, the NEO-MSDFT analytical gradients for both ground and excited vibronic states are derived and implemented. The analytical gradients and semi-numerical Hessians are used to optimize and characterize equilibrium and transition state geometries and to generate minimum energy paths (MEPs), for proton transfer in the deprotonated acetylene dimer and malonaldehyde. The barriers along the resulting MEPs are lower when the transferring proton is quantized because the NEO-MSDFT method inherently includes the zero-point energy of the transferring proton. Analysis of the proton densities along the MEPs illustrates that the proton density can exhibit symmetric or asymmetric bilobal character associated with symmetric or slightly asymmetric double-well potential energy surfaces and hydrogen tunneling. Analysis of the contributions to the intrinsic reaction coordinate reveals that changes in the C–O bond lengths drive proton transfer in malonaldehyde. This work provides the foundation for future reaction path studies and direct nonadiabatic dynamics simulations of a wide range of hydrogen transfer reactions. 

   more » « less
5. Nodeless vibrational amplitudes and quantum nonadiabatic dynamics in the nested funnel for a pseudo Jahn-Teller molecule or homodimer

   https://doi.org/10.1063/1.5009762  

   Peters, William_K; Tiwari, Vivek; Jonas, David_M (November 2017, The Journal of Chemical Physics)

   The nonadiabatic states and dynamics are investigated for a linear vibronic coupling Hamiltonian with a static electronic splitting and weak off-diagonal Jahn-Teller coupling through a single vibration with a vibrational-electronic resonance. With a transformation of the electronic basis, this Hamiltonian is also applicable to the anti-correlated vibration in a symmetric homodimer with marginally strong constant off-diagonal coupling, where the non-adiabatic states and dynamics model electronic excitation energy transfer or self-exchange electron transfer. For parameters modeling a free-base naphthalocyanine, the nonadiabatic couplings are deeply quantum mechanical and depend on wavepacket width; scalar couplings are as important as the derivative couplings that are usually interpreted to depend on vibrational velocity in semiclassical curve crossing or surface hopping theories. A colored visualization scheme that fully characterizes the non-adiabatic states using the exact factorization is developed. The nonadiabatic states in this nested funnel have nodeless vibrational factors with strongly avoided zeroes in their vibrational probability densities. Vibronic dynamics are visualized through the vibrational coordinate dependent density of the time-dependent dipole moment in free induction decay. Vibrational motion is amplified by the nonadiabatic couplings, with asymmetric and anisotropic motions that depend upon the excitation polarization in the molecular frame and can be reversed by a change in polarization. This generates a vibrational quantum beat anisotropy in excess of 2/5. The amplitude of vibrational motion can be larger than that on the uncoupled potentials, and the electronic population transfer is maximized within one vibrational period. Most of these dynamics are missed by the adiabatic approximation, and some electronic and vibrational motions are completely suppressed by the Condon approximation of a coordinate-independent transition dipole between adiabatic states. For all initial conditions investigated, the initial nonadiabatic electronic motion is driven towards the lower adiabatic state, and criteria for this directed motion are discussed. 

   more » « less