ClassSTRONG: Classical simulations of strong field processes

A set of Mathematica functions is presented to model classically two of the most important processes in strong field physics, namely high-order harmonic generation (HHG) and above-threshold ionization (ATI). Our approach is based on the numerical solution of the Newton–Lorentz equation of an electron moving on an electric field and takes advantage of the symbolic languages features and graphical power of Mathematica. Like in the Strong Field Approximation (SFA), the effects of atomic potential on the motion of electron in the laser field are neglected. The SFA was proven to be an essential tool in strong field physics in the sense that it is able to predict with great precision the harmonic (in the HHG) and energy (in the ATI) limits. We have extended substantially the conventional classical simulations, where the electric field is only dependent on time, including spatial nonhomogeneous fields and spatial and temporal synthesized fields. Spatial nonhomogeneous fields appear when metal nanosystems interact with strong and short laser pulses and temporal synthesized fields are routinely generated in attosecond laboratories around the world. Temporal and spatial synthesized fields have received special attention nowadays because they would allow to exceed considerably the conventional harmonic and electron energy frontiers. Classical simulations are an invaluable tool to explore exhaustively the parameters domain at a cheap computational cost, before massive quantum mechanical calculations, absolutely indispensable for the detailed analysis, are performed.