An efficient and state-of-the-art real-time time-dependent density functional theory (rt-TDDFT) method is presented, as implemented in the time-dependent ab initio package (TDAP), which aims at performing accurate simulations of the interaction between laser fields and solid-state materials. The combination of length-gauge and velocity-gauge electromagnetic field has extended the diversity of materials under consideration, ranging from low dimensional systems to periodic solids. Meanwhile, by employing a local basis presentation, systems of a large size are simulated for long electronic propagation time, with moderate computational cost while maintaining a relatively high accuracy. Non-perturbative phenomena in materials under a strong laser field and linear responses in a weak field can be simulated, either in the presence of ionic motions or not. Several quintessential works are introduced as examples for applications of this approach, including photoabsorption properties of armchair graphene nanoribbon, hole-transfer ultrafast dynamics between MoS 2 /WS 2 interlayer heterojunction, laser-induced nonthermal melting of silicon, and high harmonic generation in monolayer MoS 2 . The method demonstrates great potential for studying ultrafast electron-nuclear dynamics and nonequilibrium phenomena in a wide range of quantum systems.ultrafast dynamics and phenomena either in perturbative or non-perturbative regimes. [3][4][5][6][7][8] Therefore, it has been a unique ab initio quantum method applicable for the exploring of strong field physics beyond linear response theory, for instance, high harmonic generation [9,10] and ultrafast photoelectron emission. [11] Recently, the scope of rt-TDDFT applications has greatly extended from treating isolated atomic and molecular systems to condensed phase materials. In most previous works, numerical implementations of rt-TDDFT that aim at handling solid materials were built on real-space grids, [12,13] including some well-known program packages such as OCTOPUS [9,14] and SALMON. [15,16] Real-time TDDFT has also been implemented in plane wave codes, for example, the ELK FP-LAPW [17] and FPSID, [18] whose encouraging results have shown the effectiveness of the rt-TDDFT approaches. However, if one is interested in high energy excitation that is on the energy scale of tens to hundreds of electron volts (eV), extremely dense real-space grids and high kinetic energy of plane waves are indispensable. Meanwhile, to describe a system with N a atoms, 10 3 to 10 4 × N a real-space grids or plane waves have to be used, which makes the simulation of large-size systems impractical using computer resources available at the present stage. The above two factors will significantly increase the computational cost and in turn limit the practicability of the rt-TDDFT methods.Here we introduce a real-time ab initio approach based on local atomic basis for simulating electron-nuclear dynamics under laser excited conditions. This approach has been successfully implemented in the time-dependent ab initi...