Modeling of Short-Pulse Laser Interactions with Monolithic and Porous Silicon Targets with an Atomistic–Continuum Approach

• Published in: Nanomaterials (Basel, Switzerland) Vol. 13; no. 20; p. 2809
• Main Authors: Grigoryeva, Maria S., Kutlubulatova, Irina A., Lukashenko, Stanislav Yu, Fronya, Anastasia A., Ivanov, Dmitry S., Kanavin, Andrey P., Timoshenko, Victor Yu, Zavestovskaya, Irina N.
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01.10.2023; MDPI
• Subjects: Ablation; atomistic–continuum numerical method; Conduction heating; Conductive heat transfer; crystal structure; Electronic equipment; Equilibrium; Identification and classification; Laser ablation; laser ablation of semiconductors; Laser beams; Lasers; Melting; Methods; Modelling; Molecular dynamics; Nanoparticles; Optical properties; Phase transitions; porosity; Porous materials; Porous silicon; Properties; Residual stress; Semiconductor materials; Short pulses; Silicon; Temperature; Thresholds; Russia
• ISSN: 2079-4991; 2079-4991
• DOI: 10.3390/nano13202809

The acquisition of reliable knowledge about the mechanism of short laser pulse interactions with semiconductor materials is an important step for high-tech technologies towards the development of new electronic devices, the functionalization of material surfaces with predesigned optical properties, and the manufacturing of nanorobots (such as nanoparticles) for bio-medical applications. The laser-induced nanostructuring of semiconductors, however, is a complex phenomenon with several interplaying processes occurring on a wide spatial and temporal scale. In this work, we apply the atomistic–continuum approach for modeling the interaction of an fs-laser pulse with a semiconductor target, using monolithic crystalline silicon (c-Si) and porous silicon (Si). This model addresses the kinetics of non-equilibrium laser-induced phase transitions with atomic resolution via molecular dynamics, whereas the effect of the laser-generated free carriers (electron–hole pairs) is accounted for via the dynamics of their density and temperature. The combined model was applied to study the microscopic mechanism of phase transitions during the laser-induced melting and ablation of monolithic crystalline (c-Si) and porous Si targets in a vacuum. The melting thresholds for the monolithic and porous targets were found to be 0.32 J/cm2 and 0.29 J/cm2, respectively. The limited heat conduction mechanism and the absence of internal stress accumulation were found to be involved in the processes responsible for the lowering of the melting threshold in the porous target. The results of this modeling were validated by comparing the melting thresholds obtained in the simulations to the experimental values. A difference in the mechanisms of ablation of the c-Si and porous Si targets was considered. Based on the simulation results, a prediction regarding the mechanism of the laser-assisted production of Si nanoparticles with the desired properties is drawn.