## 1 Introduction
[19] discussed morphea treatment using pulsed dye lasers (585 and 595 nm), fractional lasers (CO\({}_{2}\) and Er:YAG), Excimer lasers (308 nm), and Nd:YAG (1064 nm) lasers. [PERSON] et al. [20] provided an updated overview of laser treatments for skin scars and burn wounds, covering indications, modes of operation, effectiveness, and side effects. [PERSON] [21] outlined guidelines for fractional laser use in individuals with darker skin tones, suggesting it as a preferred treatment for various dermatological conditions (Fitzpatrick skin phototypes IV to VI). Finally, [PERSON] and [PERSON] [22] introduced ablative, non-ablative, and fractional photothermolysis procedures for skin rejuvenation. [PERSON] et al. [23] employed a finite element analysis (FEA) using the COMSOL package (Version 5.4) to evaluate the impact of transitional photodynamic degradation on wine stain laser treatment. The study involved the utilization of laser pulses with a frequency of up to three per millisecond to determine the depth of thermal penetration within the tissue layer through the application of diffusion approximation theory. [PERSON] et al. [24] conducted a simulation to analyze the temperature changes and phases in materials subjected to a KrF-pulsed laser beam on a silicon surface. The pulse duration was 27 nanoseconds, and the impacts were evaluated at varying time intervals, investigating the impact of changing laser parameters. [PERSON] et al. [25] and their associates performed a study that involved modeling the impact of laser irradiation on a metal surface. They utilized an Nd:YAG laser pulse with a duration of 180 ns, which resulted in the melting of the metal surface. The process of solidification was then monitored after the laser was no longer in use. [PERSON] et al. [26] emphasized the influence of thermomechanical dynamics associated with laser therapy by utilizing a finite element model. The study utilized a three-layered model of the skin that was exposed to a focused Nd:YAG infrared laser beam. The results gained provided an understanding of how laser parameters, such as repetition rate, laser power, and pulse width, impact the ablation process. [PERSON] and [PERSON] [27] used bioheat transfer models to simulate the thermal behavior in the tumor and estimate the power dissipation requirement for the specified conditions during magnetic hyperthermia. Thermal damage was also evaluated based on the Arrhenius equation to build a modified thermal damage model with regeneration terms, and the equivalent thermal dose equation. [PERSON] et al. [28] introduced a model of magnetic nanoparticle-assisted thermal therapy for minimally invasive, non-surgical cancer treatment. Pennes bio-heat and Arrhenius equations were implemented in COMSOL to model heat transfer and evaluate thermal damage in healthy and tumorous models. [PERSON] et al. [29] proposed a mathematical heat transfer model to predict the control temperature profile at the targeted position. This technique used a dual phase-lag (DPL) model of heat transfer in multilayer tissues with a modified Gaussian distribution heat source subjected to the most generalized boundary condition and interface at the adjacent layers. [PERSON] et al. [30] suggested a thermomechanical model to simulate laser heating in shrinking tissue, exploring the impact of wavelength, laser irradiation intensity, and irradiation beam area on the process. Higher laser irradiation intensity resulted in increased temperature elevation and tissue shrinkage, with notable differences observed across varying wavelengths. Skin tumors represent a significant health concern globally, necessitating efficient and minimally invasive treatment modalities. However, existing approaches often pose challenges due to their limited precision and potential damage to surrounding healthy tissue. By delving into the optimization of laser parameters and assessing the resultant thermal effects, our study seeks to address these critical gaps in current methodologies [31].