Dynamical rate equation model for femtosecond laser-induced breakdown in dielectrics

Abstract

Experimental and theoretical studies of laser-induced breakdown in dielectrics provide conflicting conclusions about the possibility to trigger ionization avalanche on the subpicosecond time scale and the relative importance of carrier-impact ionization over field ionization. On the one hand, current models based on a single ionization-rate equation do not account for the gradual heating of the charge carriers, which, for short laser pulses, might not be sufficient to start an avalanche. On the other hand, kinetic models based on microscopic collision probabilities have led to variable outcomes that do not necessarily match experimental observations as a whole. In this paper, we present a rate-equation model that accounts for the avalanche process phenomenologically by using an auxiliary differential equation to track the gradual heating of the charge carriers and define the collisional impact rate dynamically. The computational simplicity of this dynamical rate-equation model offers the flexibility to extract effective values from experimental data. This is demonstrated by matching the experimental scaling trends for the laser-induced damage threshold of several dielectric materials for pulse durations ranging from a few fs to a few ps. Through numerical analysis, we show that the proposed model gives results comparable to those obtained with multiple rate equations and identify potential advantages for the development of large-scale, three-dimensional electromagnetic methods for the modeling of laser-induced breakdown in transparent media.