Nonadiabatic Molecular Dynamics for Thousand Atom Systems: A Tight-Binding Approach toward PYXAID

Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics

Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.

Emergence of Z-Scheme Photocatalysis for Total Water Splitting: An Improvised Route to High Efficiency

Photocatalytic water splitting to spontaneously produce H2 and O2 is a long-standing goal in solar energy conversion, presenting a significant challenge without using sacrificial electron donors or external biases. Inspired by natural photosynthesis, the design of artificial Z-scheme photocatalytic systems is at the forefront of this field. These systems achieve higher redox potential by separating photogenerated electrons and holes through a fast interlayer recombination process between valence and conduction band edges. Z-scheme photocatalysis involves using two different semiconductors with distinct bandgap energies. Here, we explore potential systems based on two-dimensional (2D) heterostructures composed of carbon, nitrogen, or similar main group elements. The advantages and disadvantages of these systems are discussed, with a focus on enhancing their efficiency through strategic design. Special emphasis is placed on the dynamics of excited charge carrier transfer and recombination processes, which are crucial for developing efficient photocatalytic systems for overall water splitting.

Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene

In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.

Interfacial Charge Transfer in Photoexcited QD-Molecule Composite of Tetrahedral CdSe Quantum Dot Coupled with Carbazole

Photoexcited charge transfer dynamics in CdSe quantum dots (QDs) coupled with carbazole were explored to model QD-molecule systems for light-harvesting applications. The absorption spectra of QDs with different sizes, i.e., Cd35Se20X30L30 (T1), Cd56Se35X42L42 (T2), and Cd84Se56X56L56 (T3) were simulated with quantum dynamical methods, which qualitatively match the reported experimental spectra. The carbazole is attached with a 3-amino group at the apex position of T1 (namely T1-3A-Cz), establishing proper electronic communication between T1 and carbazole. The spectra of T1-3A-Cz is 0.22 eV red-shifted compared to T1. A time-dependent perturbation was applied in tune with the lowest energy peak (3.63 eV) of T1-3A-Cz to investigate the charge transfer dynamics, which revealed an ultrafast charge separation within the femtosecond time scale. The electronic structure showed a favorable energy alignment between T1 and carbazole in T1-3A-Cz. The LUMO of carbazole was situated below the conduction band of the QD, while the HOMO of carbazole mixed perfectly with the top of the valence band of the QD, developing the interfacial charge transfer states. These states promoted the photoexcited electron transfer directly from the CdSe core to carbazole. A rapid and enhanced charge separation occurred with the laser field strength increasing from 0.001 to 0.005 V/Å. However, T1 connected to the other positions of carbazole did not show charge separation effectively. The photoinduced charge transfer is negligible in the case of T2-carbazole systems due to poor electronic coupling, and it is not observed in T3-carbazole systems. So, the T1-3A-Cz model acts as a perfect donor-acceptor QD-molecule nanocomposite that can harvest photon energy efficiently. Further enhancement of charge transfer can be achieved by coupling more carbazoles to the T1 QD (e.g., T1-3A-Cz2) due to the extension of hole delocalization between T1 and the carbazoles.

Prolonged Exciton Lifetime Is Achieved in Porphyrin Nanoring by Template Engineering: A Nonadiabatic Tight Binding Approach

Porphyrin nanoring has been attracting immense attention due to its light harvesting capacity and potential applications in optical, catalysis, sensor, and electronic devices. We demonstrate by nonadiabatic quantum dynamics simulations that the photovoltaic efficiency can be enhanced by template engineering. Altering the hexadentate template (T6) with two tridentate templates (2T3) within the porphyrin ring (P6) cavity accelerated the electron transfer twice and suppressed the electron-hole recombination by nearly three times. The atomistic tight-binding simulation rationalized the dynamics by different localizations of charge of the band edge states, changes in nonadiabatic coupling, alteration in quantum coherence, and involvement of diverse electron-phonon vibrational modes. Further 2T3 templates more strongly hold the P6 ring than T6, reducing the structural fluctuation. As a result, the nonadiabatic coupling becomes weaker and suppresses the carrier recombination. Current atomistic simulation presents a template engineering strategy to enhance the exciton lifetime along with ultrafast charge separation, crucial factors for photovoltaic applications.

Addressing electronic and dynamical evolution of molecules and molecular clusters: DFTB simulations of energy relaxation in polycyclic aromatic hydrocarbons

We present a review of the capabilities of the density functional based Tight Binding (DFTB) scheme to address the electronic relaxation and dynamical evolution of molecules and molecular clusters following energy deposition via either collision or photoabsorption. The basics and extensions of DFTB for addressing these systems and in particular their electronic states and their dynamical evolution are reviewed. Applications to PAH molecules and clusters, carbonaceous systems of major interest in astrochemical/astrophysical context, are reported. A variety of processes are examined and discussed such as collisional hydrogenation, fast collisional processes and induced electronic and charge dynamics, collision-induced fragmentation, photo-induced fragmentation, relaxation in high electronic states, electronic-to-vibrational energy conversion and statistical versus non-statistical fragmentation. This review illustrates how simulations may help to unravel different relaxation mechanisms depending on various factors such as the system size, specific electronic structure or excitation conditions, in close connection with experiments.

Controlling Charge Carrier Dynamics in Porphyrin Nanorings by Optically Active Templates

Understanding the dynamics of photogenerated charge carriers is essential for enhancing the performance of solar and optoelectronic devices. Using atomistic quantum dynamics simulations, we demonstrate that a short π-conjugated optically active template can be used to control hot carrier relaxation, charge carrier separation, and carrier recombination in light-harvesting porphyrin nanorings. Relaxation of hot holes is slowed by 60% with an optically active template compared to that with an analogous optically inactive template. Both systems exhibit subpicosecond electron transfer from the photoactive core to the templates. Notably, charge recombination is suppressed 6-fold by the optically active template. The atomistic time-domain simulations rationalize these effects by the extent of electron and hole localization, modification of the density of states, participation of distinct vibrational motions, and changes in quantum coherence. Extension of the hot carrier lifetime and reduction of charge carrier recombination, without hampering charge separation, demonstrate a strategy for enhancing efficiencies of energy materials with optically active templates.

Nonradiative Dynamics Induced by Vacancies in Wide-Gap III-Nitrides: Ab Initio Time-Domain Analysis

Insightful understanding of defect properties and prevention of defect damage are among the biggest issues in the development of photoelectronic devices based on wide-gap III-nitride semiconductors. Here, we have investigated the vacancy-induced carrier nonradiative dynamics in wide-gap III-nitrides (GaN, AlN, and Al_xGa_1-xN) by ab initio molecular dynamics and nonadiabatic (NA) quantum dynamics simulations since the considerable defect density in epitaxy samples. E-h recombination is hardly affected by V_cation, which created shallow states near the VBM. Our findings demonstrate that V_N in AlN creates defect-assisted nonradiative recombination centers and shortens the recombination time (τ) as in the Shockley-Read-Hall (SRH) model. In GaN, V_N improves the NA coupling between the CBM and the VBM. Additionally, increasing x in the Al_xGa_1-xN alloys accelerates nonradiative recombination, which may be an important issue in further improving the IQE of high Al-content Al_xGa_1-xN alloys. These findings have significant implications for the improvement of wide-gap III-nitrides-based photoelectronic devices.

Rapid Interlayer Charge Separation and Extended Carrier Lifetimes due to Spontaneous Symmetry Breaking in Organic and Mixed Organic-Inorganic Dion-Jacobson Perovskites

Promising alternatives to three-dimensional perovskites, two-dimensional (2D) layered metal halide perovskites have proven their potential in optoelectronic applications due to improved photo- and chemical stability. Nevertheless, photovoltaic devices based on 2D perovskites suffer from poor efficiency owing to unfavorable charge carrier dynamics and energy losses. Focusing on the 2D Dion-Jacobson perovskite phase that is rapidly rising in popularity, we demonstrate that doping of complementary cations into the 3-(aminomethyl)piperidinium perovskite accelerates spontaneous charge separation and slows down charge recombination, both factors improving the photovoltaic performance. Employing ab initio nonadiabatic (NA) molecular dynamics combined with time-dependent density functional theory, we demonstrate that cesium doping broadens the bandgap by 0.4 eV and breaks structural symmetry. Assisted by thermal fluctuations, the symmetry breaking helps to localize electrons and holes in different layers and activates additional vibrational modes. As a result, the charge separation is accelerated. Simultaneously, the charge carrier lifetime grows due to shortened coherence time between the ground and excited states. The established relationships between perovskite composition and charge carrier dynamics provide guidelines toward future material discovery and design of perovskite solar cells.

Non-negligible roles of charge transfer excitons in ultrafast excitation energy transfer dynamics of a double-walled carbon nanotube

Herein, we employed a developed linear response time dependent density functional theory-based nonadiabatic dynamics simulation method that explicitly takes into account the excitonic effects to investigate photoinduced excitation energy transfer dynamics of a double-walled carbon nanotube (CNT) model with different excitation energies. The E₁₁ excitation of the outer CNT will generate a local excitation (LE) |out*〉 exciton due to its low energy, which does not induce any charge separation. In contrast, the E₁₁ excitation of the inner CNT can generate four kinds of excitons with the LE exciton |in*〉 dominates. In the 500-fs dynamics simulation, the LE exciton |in*〉 and charge transfer (CT) excitons |out^-in⁺〉 and |out⁺in^-〉 are all gradually converted to the |out*〉 exciton, corresponding to a photoinduced excitation energy transfer, which is consistent with experimental studies. Finally, when the excitation energy is close to the E₂₂ state of the outer CNT (∼1.05 eV), a mixed population of different excitons, with the |out*〉 exciton dominated, is generated. Then, photoinduced energy transfer from the outer to inner CNTs occurs in the first 50 fs, which is followed by an inner to outer excitation energy transfer that is completed in 400 fs. The present work not only sheds important light on the mechanistic details of wavelength-dependent excitation energy transfer of a double-walled CNT model but also demonstrates the roles and importance of CT excitons in photoinduced excitation energy transfer. It also emphasized that explicitly including the excitonic effects in electronic structure calculations and nonadiabatic dynamics simulations is significant for correct understanding/rational design of optoelectronic properties of periodically extended systems.

Passivation of Hematite by a Semiconducting Overlayer Reduces Charge Recombination: An Insight from Nonadiabatic Molecular Dynamics

Hematite (α-Fe₂O₃) is a promising photoanode material for photoelectrochemical water splitting. Surface-passivating layers are effective in improving water oxidation kinetics; however, the passivation mechanism is not fully understood due to the complexity of interfacial reactions. Focusing on the Fe-terminated Fe₂O₃ (0001) surface that exhibits surface states in the band gap, we perform ab initio quantum dynamics simulations to study the effect of an α-Ga₂O₃ overlayer on charge recombination. The overlayer eliminates surface states and suppresses charge recombination 4-fold. This explains in part the observed cathodic shift in the onset potential for water oxidation. The increased charge carrier lifetime is an outcome of two factors, energy gap and electron-vibrational coupling, with a positive contribution from the former but a negative contribution from the latter. This work presents an advance in the atomistic time-domain understanding of the influence of surface passivation on charge recombination dynamics and provides guidance for designing novel α-Fe₂O₃ photoanodes.

Fluctuations at Metal Halide Perovskite Grain Boundaries Create Transient Trap States: Machine Learning Assisted Ab Initio Analysis

All-inorganic perovskites are promising candidates for solar energy and optoelectronic applications, despite their polycrystalline nature with a large density of grain boundaries (GBs) due to facile solution-processed fabrication. GBs exhibit complex atomistic structures undergoing slow rearrangements. By studying evolution of the Σ5(210) CsPbBr₃ GB on a nanosecond time scale, comparable to charge carrier lifetimes, we demonstrate that GB deformations appear every ∼100 ps and increase significantly the probability of deep charge traps. However, the deep traps form only transiently for a few hundred femtoseconds. In contrast, shallow traps appear continuously at the GB. Shallow traps are localized in the GB layer, while deep traps are in a sublayer, which is still distorted from the pristine structure and can be jammed in unfavorable conformations. The GB electronic properties correlate with bond angles, with notable exception of the Br-Br distance, which provides a signature of halide migration along GBs. The transient nature of trap states and localization of electrons and holes at different parts of GBs indicate that charge carrier lifetimes should be long. At the same time, charge mobility can be reduced. The complex, multiscale evolution of geometric and electronic structures of GBs rationalize the contradictory statements made in the literature regarding both benign and detrimental roles of GBs in perovskite performance and provide new atomistic insights into perovskite properties.

How Good Is the Vibronic Hamiltonian Repetition Approach for Long-Time Nonadiabatic Molecular Dynamics?

Multiple applied studies of slow nonadiabatic processes in nanoscale and condensed matter systems have adopted the "repetition" approximation in which long trajectories for such simulations are obtained by concatenating shorter trajectories, directly available from ab initio calculations, many times. Here, we comprehensively assess this approximation using model Hamiltonians with parameters covering a wide range of regimes. We find that state transition time scales may strongly depend on the length of the repeated data, although the convergence is not monotonic and may be slow. The repetition approach may under- or overestimate the time scales by a factor of ≤7-8, does not directly depend on the dispersion of energy gap and nonadiabatic coupling (NAC) frequencies, but may depend on the magnitude of the NACs. We suggest that the repetition-based nonadiabatic dynamics may be inaccurate in simulations with very small NACs, where intrinsic transition times are on the order of ≥100 ps.

Symmetrical Linkage in Porphyrin Nanoring Suppressed the Electron-Hole Recombination Demonstrated by Nonadiabatic Molecular Dynamics

Macromolecular porphyrin nanorings are receiving significant attention because of their excellent optoelectronic properties. However, their efficiencies as potential solar materials are significantly affected by nonradiative charge recombination. To understand the recombination mechanism by alternating structural parameters and using tight-binding nonadiabatic molecular dynamics, we demonstrate that charge recombination depends strongly on the mode of the linker in the porphyrin nanoring. The nanoring having all-butadiyne-linkage (pristine-P8) inhibits carrier relaxation. In contrast, a partially fused nanoring (fused-P8) expedites the rate of quantum transition. An extension of the lifetime by a factor of 4 is due to the larger optical gap in pristine-P8 that reduces the NA coupling by decreasing the overlap between band edge states. Additionally, an intense phonon peak in the low-frequency region and rapid coherence loss within the electronic subsystem favors prolonging the carrier lifetime. This study provides an atomistic realization for the design of macromolecular porphyrin nanorings for the potential use in photovoltaic materials.

Detailed Balance and Independent Electron Surface-Hopping Method: The Importance of Decoherence and Correct Calculation of Diabatic Populations

We benchmark and improve the independent electron surface-hopping (IESH) method developed by J. C. Tully's group for nonadiabatic simulations near metal surfaces. We have incorporated decoherence within the IESH method as well as implemented a scheme for the accurate calculation of diabatic populations. We benchmark the original IESH method with the above inclusions for a model system to calculate rate constants and long-time populations. The original IESH method fails to capture the detailed balance for some of the parameters, which is corrected with the inclusion of decoherence and accurate calculation of diabatic populations. Total rate constants are well captured both within the original IESH method as well as within our modified IESH.

Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids

In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.

Simulations of molecular photodynamics in long timescales

Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Ag-Bi Charge Redistribution Creates Deep Traps in Defective Cs2AgBiBr6: Machine Learning Analysis of Density Functional Theory

Lead-free double perovskites hold promise for stable and environmentally benign solar cells; however, they exhibit low efficiencies because defects act as charge recombination centers. Identifying trap-assisted loss mechanisms and developing defect passivation strategies constitute an urgent goal. Applying unsupervised machine learning to density functional theory and nonadiabatic molecular dynamics, we demonstrate that negatively charged Br vacancies in Cs₂AgBiBr₆ create deep hole traps through charge redistribution between the adjacent Ag and Bi atoms. Vacancy electrons are first accepted by Bi and then shared with Ag, as the trap transforms from shallow to deep. Subsequent charge losses are promoted by Ag and Bi motions perpendicular to rather than along the Ag-Bi axis, as can be expected. In contrast, charge recombination in pristine Cs₂AgBiBr₆ correlates most with displacements of Cs atoms and Br-Br-Br angles. Doping with In to replace Ag at the vacancy maintains the electrons at Bi and keeps the trap shallow.

Challenges of modeling nanostructured materials for photocatalytic water splitting

Understanding the water splitting mechanism in photocatalysis is a rewarding goal as it will allow producing clean fuel for a sustainable life in the future. However, identifying the photocatalytic mechanisms by modeling photoactive nanoparticles requires sophisticated computational techniques based on multiscale modeling. In this review, we will survey the strengths and drawbacks of currently available theoretical methods at different length and accuracy scales. Understanding the surface-active site through Density Functional Theory (DFT) using new, more accurate exchange-correlation functionals plays a key role for surface engineering. Larger scale dynamics of the catalyst/electrolyte interface can be treated with Molecular Dynamics albeit there is a need for more generalizations of force fields. Monte Carlo and Continuum Modeling techniques are so far not the prominent path for modeling water splitting but interest is growing due to the lower computational cost and the feasibility to compare the modeling outcome directly to experimental data. The future challenges in modeling complex nano-photocatalysts involve combining different methods in a hierarchical way so that resources are spent wisely at each length scale, as well as accounting for excited states chemistry that is important for photocatalysis, a path that will bring devices closer to the theoretical limit of photocatalytic efficiency.

All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning

Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.

Scalable Ehrenfest Molecular Dynamics Exploiting the Locality of Density-Functional Tight-Binding Hamiltonian

To explore the science behind excited-state dynamics in high-complexity chemical systems, a scalable nonadiabatic molecular dynamics (MD) technique is indispensable. In this study, by treating the electronic degrees of freedom at the density-functional tight-binding level, we developed and implemented a reduced scaling and multinode-parallelizable Ehrenfest MD method. To achieve this goal, we introduced a concept called patchwork approximation (PA), where the effective Hamiltonian for real-time propagation of the electronic density matrix is partitioned into a set of local parts. Numerical results for giant icosahedral fullerenes, which comprise up to 6000 atoms, suggest that the scaling of the present PA-based method is less than quadratic, which yields a significant advantage over the conventional cubic scaling method in terms of computational time. The acceleration by the parallelization on multiple nodes was also assessed. Furthermore, the electronic and structural dynamics resulting from the perturbation by the external electric field were accurately reproduced with the PA, even when the electronic excitation was spatially delocalized.

Modeling Non-adiabatic Dynamics in Nanoscale and Condensed Matter Systems

Rapid, far-from-equilibrium processes involving excitation of electronic, vibrational, spin, photon, topological, and other degrees of freedom form the basis of modern technologies, including electronics and optoelectronics, solar energy harvesting and conversion to electrical and chemical energy, quantum information processing, spin- and valleytronics, chemical detection, and medical therapies. Such processes are studied experimentally with various time-resolved spectroscopies that allow scientists to track system's evolution on ultrafast time scales and at close to atomistic level of detail. The availability of various forms of lasing has made such measurements easily accessible to many experimental groups worldwide, to study atoms and small molecules, nanoscale and condensed matter systems, proteins, cells, and mesoscopic materials. The experimental work necessitates parallel theoretical efforts needed to interpret the experiments and to provide insights that cannot be gained through measurements due to experimental limitations.Non-adiabatic (NA) molecular dynamics (MD) allows one to study processes at the atomistic level and in the time domain most directly mimicking the time-resolved experiments. Atomistic modeling takes full advantage of chemical intuition and principles that guide design and fabrication of molecules and materials. It provides atomistic origins of quasi-particles, such as holes, excitons, trions, plasmons, phonons, polarons, polaritons, spin-waves, momentum-resolved and topological states, electrically and magnetically polarized structures, and other abstract concepts. An atomistic description enables one to study realistic aspects of materials, which necessarily contain defects, dopants, surfaces, interfaces, passivating ligands, and solvent layers. Often, such realistic features govern material properties and are hard to account for phenomenologically. NA-MD requires few approximations and assumptions. It does not need to assume that atomic motions are harmonic, that electrons are Drude oscillators, that coupling between different degrees of freedom is weak, that dynamics is Markovian or has short memory, or that evolution occurs by exponential kinetics of transitions between few states. The classical or semiclassical treatment of atomic motions constitutes the main approximation of NA-MD and is used because atoms are 3-5 orders of magnitude heavier than electrons. NA-MD is limited by system size, typically hundreds or thousands of atoms, and time scale, picoseconds to nanoseconds. The quality of NA-MD simulations depends on the electronic structure method used to obtain excited state energies and NA couplings.NA-MD has been largely popularized and advanced in the chemistry community that focuses on molecules. Modeling far-from-equilibrium dynamics in nanoscale and condensed matter systems often has to account for other types of physics. At the same time, condensed phase NA-MD allows for approximations that may not work in molecules. Focusing on the recent NA-MD developments aimed at studying excited state processes in nanoscale and condensed phases, this Account considers how the phenomena important on the nanoscale can be incorporated into NA-MD and what approximations can be made to increase its efficiency with complex systems and processes.

Accelerating molecular property calculations with semiempirical preconditioning

Computing ab initio molecular linear response properties, e.g., electronic excitation energies and transition dipole moments, requires the solution of large eigenvalue problems or large systems of equations. These large eigenvalue problems or large systems of equations are commonly solved iteratively using Krylov space algorithms, such as the Davidson algorithm for eigenvalue problems. A critical ingredient in Krylov space algorithms is the preconditioner, which is used to generate optimal update vectors in each iteration. We propose to use semiempirical approximations as preconditioners to accelerate the calculation of ab initio properties. The crucial advantage to improving the preconditioner is that the converged result is unchanged, so there is no trade-off between accuracy and speedup. We demonstrate our approach by accelerating the calculation of electronic excitation energies and electric polarizabilities from linear response time-dependent density functional theory using the simplified time-dependent density functional theory semiempirical model. For excitation energies, the semiempirical preconditioner reduces the number of iterations on average by 37% and up to 70%. The semiempirical preconditioner reduces the number of iterations for computing the polarizability by 15% on average and up to 33%. Moreover, we show that the preconditioner can be further improved by tuning the empirical parameters that define the semiempirical model, leading to an additional reduction in the number of iterations by about 20%. Our approach bridges the gap between semiempirical models and ab initio methods and charts a path toward combining the speed of semiempirical models with the accuracy of ab initio methods.

Common Defects Accelerate Charge Carrier Recombination in CsSnI3 without Creating Mid-Gap States

Lead-free metal halide perovskites are environmentally friendly and have favorable electro-optical properties; however, their efficiencies are significantly below the theoretical limit. Using ab initio nonadiabatic molecular dynamics, we show that common intrinsic defects accelerate nonradiative charge recombination in CsSnI₃ without creating midgap traps. This is in contrast to Pb-based perovskites, in which many defects have little influence on and even prolong carrier lifetimes. Sn-related defects, such as Sn vacancies and replacement of Sn with Cs are most detrimental, since Sn removal breaks the largest number of bonds and strongly perturbs the Sn-I lattice that supports the carriers. The defects increase the nonadiabatic electron-vibrational coupling and interact strongly with free carrier states. Point defects associated with I atoms are less detrimental, and therefore, CsSnI₃ synthesis should be performed in Sn rich conditions. The study provides an atomistic rationalization of why lead-free CsSnI₃ exhibits lower photovoltaic efficiency compared to some lead-based perovskites.

On application of deep learning to simplified quantum-classical dynamics in electronically excited states

Deep learning (DL) is applied to simulate non-adiabatic molecular dynamics of phenanthrene, using the time-dependent density functional based tight binding (TD-DFTB) approach for excited states combined with mixed quantum–classical propagation. Reference calculations rely on Tully’s fewest-switches surface hopping (FSSH) algorithm coupled to TD-DFTB, which provides electronic relaxation dynamics in fair agreement with various available experimental results. Aiming at describing the coupled electron-nuclei dynamics in large molecular systems, we then examine the combination of DL for excited-state potential energy surfaces (PESs) with a simplified trajectory surface hopping propagation based on the Belyaev–Lebedev (BL) scheme. We start to assess the accuracy of the TD-DFTB approach upon comparison of the optical spectrum with experimental and higher-level theoretical results. Using the recently developed SchNetPack (Schütt et al 2019 J. Chem. Theory Comput. 15 448–55) for DL applications, we train several models and evaluate their performance in predicting excited-state energies and forces. Then, the main focus is given to the analysis of the electronic population of low-lying excited states computed with the aforementioned methods. We determine the relaxation timescales and compare them with experimental data. Our results show that DL demonstrates its ability to describe the excited-state PESs. When coupled to the simplified BL scheme considered in this study, it provides reliable description of the electronic relaxation in phenanthrene as compared with either the experimental data or the higher-level FSSH/TD-DFTB theoretical results. Furthermore, the DL performance allows high-throughput analysis at a negligible cost.

Ab initio symmetric quasi-classical approach to investigate molecular Tully models

We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.

Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.

Common Defects Accelerate Charge Separation and Reduce Recombination in CNT/Molecule Composites: Atomistic Quantum Dynamics

Carbon nanotubes (CNTs) are appealing candidates for solar and optoelectronic applications. Traditionally used as electron sinks, CNTs can also perform as electron donors, as exemplified by coupling with perylenediimide (PDI). To achieve high efficiencies, electron transfer (ET) should be fast, while subsequent charge recombination should be slow. Typically, defects are considered detrimental to material performance because they accelerate charge and energy losses. We demonstrate that, surprisingly, common CNT defects improve rather than deteriorate the performance. CNTs and other low dimensional materials accommodate moderate defects without creating deep traps. At the same time, charge redistribution caused by CNT defects creates an additional electrostatic potential that increases the CNT work function and lowers CNT energy levels relative to those of the acceptor species. Hence, the energy gap for the ET is decreased, while the gap for the charge recombination is increased. The effect is particularly important because charge acceptors tend to bind near defects due to enhanced chemical interactions. The time-domain simulation of the excited-state dynamics provides an atomistic picture of the observed phenomenon and characterizes in detail the electronic states, vibrational motions, inelastic and elastic electron-phonon interactions, and time scales of the charge separation and recombination processes. The findings should apply generally to low-dimensional materials, because they dissipate defect strain better than bulk semiconductors. Our calculations reveal that CNT performance is robust to common defects and that moderate defects are essential rather than detrimental for CNT application in energy, electronics, and related fields.

A Simple Solution to Trivial Crossings: A Stochastic State Tracking Approach

We present a new state tracking algorithm based on a stochastic state reassignment that reflects the quantum mechanical interpretation of the state time-overlaps. We assess the new method with a range of model Hamiltonians and demonstrate that it yields the results generally consistent with the deterministic min-cost algorithm. However, the stochastic state tracking algorithm reduces magnitudes of the state population fluctuations as the quantum system evolves toward its equilibrium. The new algorithm facilitates the thermalization of quantum state populations and suppresses the population revivals and oscillations near the equilibrium in many-state systems. The new stochastic algorithm has a favorable computational scaling, is easy to implement due to its conceptual transparency, and treats various types of state identity changes (trivial or avoided crossings and any intermediate cases) on equal footing.

Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings

Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.

Photoinduced Dynamics of Charge Carriers in Metal Halide Perovskites from an Atomistic Perspective

Perovskite solar cells have attracted intense attention over the past decade because of their low cost, abundant raw materials, and rapidly growing power conversion efficiency (PCE). However, nonradiative charge carrier losses still constitute a major factor limiting the PCE to well below the Shockley-Queisser limit. This Perspective summarizes recent atomistic quantum dynamics studies on the photoinduced excited-state processes in metal halide perovskites (MHPs), including both hybrid organic-inorganic and all-inorganic MHPs and three- and two-dimensional MHPs. The simulations, performed using a combination of time-domain ab initio density functional theory and nonadiabatic molecular dynamics, allow emphasis on various intrinsic and extrinsic features, such as components, structure, dimensionality and interface engineering, control and exposure to various environmental factors, defects, surfaces, and their passivation. The detailed atomistic simulations advance our understanding of electron-vibrational dynamics in MHPs and provide valuable guidelines for enhancing the performance of perovskite solar cells.

Structural rigidity accelerates quantum decoherence and extends carrier lifetime in porphyrin nanoballs: a time domain atomistic simulation

Nonradiative electron-hole (e-h) recombination is the primary source of energy loss in photovoltaic cells and inevitably, it competes with the charge transfer process, leading to poor device performance. Therefore, much attention has to be paid for delaying such processes; increasing the excitonic lifetime may be a solution for this. Using the real-time, density functional tight-binding theory (DFTB) combined with nonadiabatic molecular dynamics (NAMD) simulations, we demonstrate the exciton relaxation phenomena of different metal-centered porphyrin nanoballs, which are supposed to be very important for the light-harvesting process. It has been revealed that the carrier recombination rate gradually decreases with the increase in the molecular stiffness by introducing metal-coordinating templating agents into the nanoball. Our simulation demonstrates that the lower atomic fluctuations lead to poorer electron-phonon nonadiabatic coupling in association with weak phonon modes and these as a whole are responsible for shorter quantum coherence and hence delayed recombination events. Our analysis is in good agreement with the recent experimental observation. By replacing the Zn metal center with a heavier Cd atom, a similar trend is observed; however, the rate slows down abruptly. The present simulation study provides the fundamental mechanism in detail behind the undesired energy loss during exciton recombination and suggests a rational design of impressive nanosystems for future device fabrication.

Theoretical investigation of the electronic relaxation in highly excited chrysene and tetracene: The effect of armchair vs zigzag edge

Non-adiabatic molecular dynamics of neutral chrysene and tetracene molecules is investigated using Tully's fewest switches surface hopping algorithm coupled to the time-dependent density functional based tight-binding (TD-DFTB) method for electronic structure calculations. We first assess the performance of two DFTB parameter sets based on the computed TD-DFTB absorption spectra. The main focus is given to the analysis of the electronic relaxation from the brightest excited state following absorption of a UV photon. We determine the dynamical relaxation times and discuss the underlying mechanisms. Our results show that the electronic population of the brightest excited singlet state in armchair-edge chrysene decays an order-of-magnitude faster than the one in zigzag-edge tetracene. This is correlated with a qualitatively similar difference of energy gaps between the brightest state and the state lying just below in energy, which is also consistent with our previous study on polyacenes.

Hot Electron Cooling in Silicon Nanoclusters via Landau-Zener Nonadiabatic Molecular Dynamics: Size Dependence and Role of Surface Termination

We develop a new express methodology for modeling excited-state dynamics occurring in dense manifolds of electronic states in atomistic systems. The approach leverages a modified Landau-Zener formula, the neglect of a back-reaction approximation, and the highly efficient density functional tight-binding method. We study the hot electron dynamics in a series of H- and F-terminated silicon nanocrystals (NCs) containing up to several hundred atoms. We explain the slower electron cooling dynamics in F-terminated NCs by the larger energy gaps between the adjacent electronic states in these systems as well as their slower fluctuations. We conclude that both the mass and chemical identity of the surface termination groups equally influence the electron dynamics, on average. However, the mass effect becomes dominant for higher-energy excitations. We find that the electron decay dynamics in F-terminated NCs has a greater sensitivity to the mass of the surface ligands than do the H-terminated NCs and explain this observation by the details of the electron-phonon coupling in the systems. We find that in the H-terminated NCs, electronic transitions in the cooling process occur predominantly between the surface states, whereas in F-terminated Si NCs, both surface and NC core states are coupled to the nuclear vibrations. We find that electron energy relaxation is accelerated in larger NCs and attribute this effect to the higher densities of states and smaller energy gaps in these systems.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Coherent Exciton Dynamics in Ensembles of Size-Dispersed CdSe Quantum Dot Dimers Probed via Ultrafast Spectroscopy: A Quantum Computational Study

Interdot coherent excitonic dynamics in nanometric colloidal CdSe quantum dots (QD) dimers lead to interdot charge migration and energy transfer. We show by electronic quantum dynamical simulations that the interdot coherent response to ultrashort fs laser pulses can be characterized by pump-probe transient absorption spectroscopy in spite of the inevitable inherent size dispersion of colloidal QDs. The latter, leading to a broadening of the excitonic bands, induce accidental resonances that actually increase the efficiency of the interdot coupling. The optical electronic response is computed by solving the time-dependent Schrodinger equation including the interaction with the oscillating electric field of the pulses for an ensemble of dimers that differ by their size. The excitonic Hamiltonian of each dimer is parameterized by the QD size and interdot distance, using an effective mass approximation. Local and charge transfer excitons are included in the dimer basis set. By tailoring the QD size, the excitonic bands can be tuned to overlap and thus favor interdot coupling. Computed pump-probe transient absorption maps averaged over the ensemble show that the coherence of excitons in QD dimers that lead to interdot charge migration can survive size disorder and could be observed in fs pump-probe, four-wave mixing, or covariance spectroscopy.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

EMPIRE: a highly parallel semiempirical molecular orbital program: 3: Born-Oppenheimer molecular dynamics

Direct NDDO-based Born-Oppenheimer molecular dynamics (MD) have been implemented in the semiempirical molecular orbital program EMPIRE. Fully quantum mechanical MD simulations on unprecedented time and length scales are possible, since the calculation of self-consistent wavefunctions and gradients is performed in a massively parallel manner. MD simulations can be performed in the NVE and NVT ensembles, using either deterministic (Berendsen) or stochastic (Langevin) thermostats. Furthermore, dynamics for condensed-phase systems can be performed under periodic boundary conditions. We show three exemplary applications: the dynamics of molecular reorganization upon ionization, long timescale dynamics of an endohedral fullerene, and calculation of the vibrational spectrum of a nanoparticle consisting of more than eight hundred atoms.

Strain Controls Charge Carrier Lifetimes in Monolayer WSe2: Ab Initio Time Domain Analysis

Mono- and few-layer transition metal dichalcogenides (TMDs) are among the most appealing candidates for electronic and optoelectronic devices. During synthesis, TMDs actively interact with substrates, which induce notable strain and influence significantly charge carriers in TMDs. By performing time-domain ab initio simulations on monolayer WSe₂, we demonstrate that charge carrier lifetimes vary by a factor of 3 within a typical 1% strain range, the bandgap changes by 0.2 eV, and electron-phonon interactions vary by 60%. Fortuitously, the most common tensile strain extends the lifetimes. The changes arise because of modifications in interatomic interactions. Further, compared to the optimized structure, at ambient temperature the bandgap drops by 0.1 eV and fluctuates by 0.1 eV. WSe₂ obeys linear response within 1% strain; however, further strain leads to nonlinear qualitative changes in WSe₂ electronic properties. The conduction band is affected more strongly than the valence band. Charges couple to phonons within a 100-400 cm^-1 frequency range, with the strongest coupling to in-plane and out-of-plane modes at 250 cm^-1. The reported findings agree with the available experiments and should be generic to other 2D materials. The strain effects need to be considered during TMD synthesis and provide means to control and tune TMD properties for 2D device applications.

How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm

Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.

Electron-Phonon Scattering Is Much Weaker in Carbon Nanotubes than in Graphene Nanoribbons

Carbon nanotubes (CNTs) and graphene nanoribbons (GNRs) are lower-dimensional derivatives of graphene. Similar to graphene, they exhibit high charge mobilities; however, in contrast to graphene, they are semiconducting and thus are suitable for electronics, optics, solar energy devices, and other applications. Charge carrier mobilities, energies, and lifetimes are governed by scattering with phonons, and we demonstrate, using ab initio nonadiabatic molecular dynamics, that charge-phonon scattering is much stronger in GNRs. Focusing on a GNR and a CNT of similar size and electronic properties, we show that the difference arises because of the significantly higher stiffness of the CNT. The GNR undergoes large-scale undulating motions at ambient conditions. Such thermal geometry distortions localize wave functions, accelerate both elastic and inelastic charge-phonon scattering, and increase the rates of energy and carrier losses. Even though, formally, both CNTs and GNRs are quantum confined derivatives of graphene, charge-phonon scattering differs significantly between them. Showing good agreement with time-resolved photoconductivity and photoluminescence measurements, the study demonstrates that GNRs are quite similar to molecules, such as conjugated polymers, while CNTs exhibit extended features attributed to bulk materials. The state-of-the-art simulations alter the traditional view of graphene nanostructures and demonstrate that the performance can be tuned not only by size and composition but also by stiffness and response to thermal excitation.