Furthermore,its analytical form allows for efficient numerical implementation, making it a practical and reliable choice in constructing microscopic OMP [36]. The SHF method can be used to determine essential nuclear properties such as binding energy, nucleon density, electromagnetic moments, and mass and charge radii of certain heavy, radioactive nuclei [19; 29]. Essentially, a nucleon incident on a nucleus can either scatter elastically or induce various nuclear reactions. To model this phenomenon, the many-body nature of the interaction is approximated by a single-body (effective) potential. [PERSON]'s [37] study of neutron reactions using a purely real potential failed to reproduce physical observables consistent with experimental data. The inability of such a real potential to account for inelastic processes was addressed by [PERSON]'s compound nucleus model [38]. However, the [PERSON] model was also found to be inadequate at higher incident energies. The OMP was introduced as a solution, defined as a complex potential with both real and imaginary components--an analogy to the scattering of light in an optical medium. This concept was first implemented in a semi-classical framework by [PERSON] et al. [39]. The use of a complex potential to describe proton scattering was initially applied by [PERSON] et al. [40]. The first formal expression of the OMP in the literature used a square well potential [41]. Later, [PERSON] and [PERSON] [42] performed a complete quantum mechanical calculation using a complex potential to describe the scattering of 17 MeV protons by aluminum, successfully reproducing the then imprecise experimental data. OMP remains a widely used model in nuclear physics today. For example, [PERSON] and [PERSON] [43] calculated global OMP parameters for spherical nuclei with mass numbers in the range \(24\leq\mathrm{A}\leq 209\) and for incident energies between \(0.001\leq\mathrm{E}\leq 200\) MeV. [PERSON] and [PERSON] [44] employed the optical model (OM) in studies related to nuclear astrophysics, while [PERSON] [45] applied it in inelastic scattering calculations. In the folding model, the nn potential can be obtained by integrating the effective NN interaction over the matter distributions of the two nuclei, as detailed by [PERSON] and [PERSON] [46; 47]. In the specific case of scattering involving a single nucleon and a target nucleus, the real part of the OMP for nucleon-nucleus (Nn) scattering can be derived by integrating the effective NN interaction with the density distribution of the target nucleus [46]. [PERSON] et al. [48] used the folding model potential to successfully explain nn elastic scattering data. In this approach, the Nn potential is calculated as a single-fold integral of the target's density distribution with a suitable effective NN interaction. Accordingly, the Single Folding (SF) model is established using the standard single-fold formulation of the OMP [46]. In this study, the elastic scattering of neutrons from tungsten isotopes was investigated through a simulation-based approach comprising three main stages:  
The nucleon density distributions for the selected tungsten isotopes (\({}^{182}\)W, \({}^{183}\)W, \({}^{184}\)W, and \({}^{186}\)W) were first calculated using the Skyrme-Hartree-Fock (SHF) method with a Woods-Saxon (WS) potential framework. Based on the obtained density distributions, Nn interaction potentials were derived using the SF method and compared with values reported in the literature. These interaction potentials were then implemented as input into the TALYS 1.96 nuclear reaction simulation code to compute the angular differential cross sections for neutron elastic scattering. To evaluate the accuracy and validity of the modeling approach, the simulation results were compared graphically with both theoretical predictions from the literature and experimental data.