The ultrafast electronic and structural dynamics invoked by photoexcitation, i.e., “dynamic exciton” phenomena, not only are important in the context of practical applications such as solar cells, but also raise many questions from the viewpoint of fundamental science. Experimental measurement, computational simulations, and theoretical interpretation will be the three pillars for deciphering the dynamic exciton phenomena.From the viewpoint of computational simulations, molecular dynamics (MD) techniques combined with quantum chemical calculations, i.e., quantum molecular dynamics (QMD), has been the popular tool to simulate the dynamic exciton phenomena. The quantum chemical calculations, which are typically conducted on the basis of the density-functional theory (DFT), require the large computational resources and time. These have been the limiting factors for the accessible spatial and time scales by the QMD simulations.To extend the coverage of simulations to more complex, large-scale systems, we have developed efficient excited-state QMD methods that can include nonadiabatic effects. Our method combines the density-functional tight binding (DFTB) method, which is an approximate DFT, and the surface hopping method, which is a theoretical framework to incorporate the nonadiabatic effects into the QMD simulations. The method was further improved to be suitable for condensed-phase simulations explicitly including the environment, i.e., solvent, by using a “divide-and-conquer” style quantum chemical calculation technique. These theoretical framework enables us to simulate the coupled electronic–structural dynamics in excited states of systems consisting of 102–103 atoms[1,2,3].In addition, using the developed method, we conducted the real-time simulations of the ultrafast processes invoked by photoexcitation of lead iodide perovskites, which are known as the key materials for perovskite solar cells. The dissociation of the exciton into the positive and negative charge carriers was observed. Moreover, the hot carrier cooling, where the charge carriers dissipate excess energy via the electron–phonon coupling and relax to the band edges, was also tracked. Finally, the direct evidence of the polaron formation, where the structural deformation is induced by the presence of charge carriers, was observed. These results highlight the importance of the coupling between electronic and structural degrees of freedom.In the talk, recent improvements in the methodology and future perspectives will also be presented[5].