“…The formation of microstructures such as grooves and spikes was then explained by Tsibidis et al based on the solution of incompressible Navier-Stokes equations, indicating the generation of hydrothermal waves and convection rolls for low and high Prandtl number fluids accordingly due to thermal gradients created by previously formed LSFL [39,40]. Nevertheless, the existing models fail to explain the formation of regular HSFL with orientations both perpendicular and parallel to laser polarization and periods down to λ/15 on metal surfaces [22,43], as well as more complex morphologies, consisting of hierarchical periodic surface structures [19,21,22,43,59,60,[67][68][69] and quasihexagonal arrangements [36, 47, 49-53, 70, 71]. For instance, the transitions between HSFL and LSFL are contradictory discussed in the literature, where both lithographic scenarios, based on electromagnetic field enhancement between regular rectangu-lar grooves and period reduction or grating splitting [57,[72][73][74][75], and hydrodynamic scenarios based on the instability [48,59,60] were proposed.…”