Abstract

Large-scale simulations using interatomic potentials provide deep insight into the processes occurring in solids subject to external perturbations. The atomistic description of laser-induced ultrafast nonthermal phenomena, however, constitutes a particularly difficult case and has so far not been possible on experimentally accessible length scales and timescales because of two main reasons: (i) ab initio simulations are restricted to a very small number of atoms and ultrashort times and (ii) simulations relying on electronic temperature- (Te) dependent interatomic potentials do not reach the necessary ab initio accuracy. Here we develop a self-learning method for constructing Te-dependent interatomic potentials which permit ultralarge-scale atomistic simulations of systems suddenly brought to extreme nonthermal states with density-functional theory (DFT) accuracy. The method always finds the global minimum in the parameter space. We derive a highly accurate analytical Te-dependent interatomic potential Φ(Te) for silicon that yields a remarkably good description of laser-excited and -unexcited Si bulk and Si films. Using Φ(Te) we simulate the laser excitation of Si nanoparticles and find strong damping of their breathing modes due to nonthermal melting.

Bernd Bauerhenne et al., Self-Learning Method for Construction of Analytical Interatomic Potentials to Describe Laser-Excited Materials, Phys. Rev. Lett.