Finite element analysis of silver nanorods, spheres, ellipsoids and core-shell structures for hyperthermia treatment of cancer

  Citation: Daud, M.U.; Abbas, G.; Afzaal, M.; Naz, M.Y.; Fatima, N.G.; Ghuffar, A.; Irfan, M.; Mahnashi, M.H.; Legutko, S.; Petr˚u, J.; et al. Finite Element Analysis of Silver Nanorods, Spheres, Ellipsoids and Core–Shell Structures for Hyperthermia Treatment of Cancer. Materials 2022,15, 1786. https:// doi.org/10.3390/ma15051786 Academic Editors: Weronika Kruszelnicka and Andrzej Tomporowski Received: 21 December 2021 Accepted: 21 February 2022 Published: 26 February 2022 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Copyright: © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). materials Article Finite Element Analysis of Silver Nanorods, Spheres, Ellipsoids and Core–Shell Structures for Hyperthermia Treatment of Cancer Muhammad Usama Daud 1, Ghulam Abbas 1,*, Muhammad Afzaal 1,*, Muhammad Yasin Naz 2, Nazma Goher Fatima 1, Abdul Ghuffar 1, Muhammad Irfan 3, Mater H. Mahnashi 4, Stanislaw Legutko 5, Jana Petr˚u 6, JiˇríKratochvíl6and Usama Muhammad Niazi 7 1Department of Physics, Faisalabad Campus, Riphah International University, Faisalabad 44000, Pakistan; [email protected] (M.U.D.); [email protected] (N.G.F.); [email protected] (A.G.) 2Department of Physics, University of Agriculture, Faisalabad 38040, Pakistan; [email protected] 3Electrical Engineering Department, College of Engineering, Najran University, Najran 61441, Saudi Arabia; [email protected] 4Department of Pharmaceutical Chemistry, College of Pharmacy, Najran University, Najran 11001, Saudi Arabia; [email protected] 5Faculty of Mechanical Engineering, Poznan University of Technology, 60-965 Poznan, Poland; stanislaw[email protected] 6Faculty of Mechanical Engineering, VSB—Technical University of Ostrava, Poruba, 708 00 Ostrava, Czech Republic; [email protected] (J.P.); [email protected] (J.K.) 7Department of Mechanical Engineering Technology, National Skills University, Islamabad 44000, Pakistan; [email protected] *Correspondence: [email protected] (G.A.); [email protected] (M.A.) Abstract: The finite element analysis technique was used to investigate the suitability of silver nanorods, spheres, ellipsoids and core–shell structures for the hyperthermia treatment of cancer. The temperature of the silver nanostructures was raised from 42 to 46 ◦ C, in order to kill the cancerous cells. The time taken by the nanostructures to attain this temperature, with external source heating, was also estimated. The heat transfer module in COMSOL Multiphysics was used for the finite element analysis of hyperthermia, based on silver nanostructures. The thermal response of different shapes of silver nanostructures was evaluated by placing them inside the spherical domain of the tumor tissue. The proposed geometries were heated at different time intervals. Optimization of the geometries was performed to achieve the best treatment temperature. It was observed that silver nanorods quickly attain the desired temperature, as compared to other shapes. The silver nanorods achieved the highest temperature of 44.3 ◦ C among all the analyzed geometries. Moreover, the central volume, used to identify the thermal response, was the maximum for the silver nano-ellipsoids. Thermal equilibrium in the treatment region was attained after 0.5 µs of heating, which made these structures suitable for hyperthermia treatment. Keywords: COMSOL Multiphysics; hyperthermia; surface coating; finite element analyses; silver nanostructures 1. Introduction Cancer is a multifactorial illness produced by a complex combination of hereditary and environmental variables [ 1 – 3 ]. Treating cancer is the most challenging health issue in the 21st century [ 4 – 7 ]. The cells of cancer are found to aggressively invade other bodily regions. These cells either form a tumor together, or can disseminate to the blood stream or lymph system [ 8 , 9 ]. These cells can metastasize to other organs and develop new tumors in places far away from the initial illness site. The ability of cancer cells to metastasize to other areas of the body is determined by a variety of factors, including blood flow and the type of cancer Materials 2022,15, 1786. https://doi.org/10.3390/ma15051786 https://www.mdpi.com/journal/materials Materials 2022,15, 1786 2 of 13 cells, and also the initial location of the cancer [ 8 – 10 ]. The substantial advancements in cancer research have resulted in a better understanding of cancer at the genetic, molecular, and cellular levels, allowing for novel therapeutic targets and procedures [ 11 ]. Cancer treatments include chemotherapy, surgery and radiation therapy [ 12 , 13 ]. Another way of treating cancer is by using hyperthermia, which involves heating the tumor region without damaging the normal cells [ 14 ]. Using hyperthermia, the temperature at the tumor location is raised to a certain level to kill the cancer cells [ 15 ]. The targeted tumor sites must attain a temperature in the range of 42–46 ◦ C for hyperthermia to be therapeutically effective in cancer therapy [ 16 ]. Protein denaturation occurs when cells are exposed to this temperature range, resulting in a high fraction of co-aggregated denatured proteins [ 15 ]. Furthermore, a high temperature influences the cellular structure function and alters intracellular processes, ultimately leading to cancer mortality [16,17], as illustrated in Figure 1. Materials 2022, 15, x FOR PEER REVIEW 2 of 14 tumors in places far away from the initial illness site. The ability of cancer cells to metastasize to other areas of the body is determined by a variety of factors, including blood flow and the type of cancer cells, and also the initial location of the cancer [8–10]. The substantial advancements in cancer research have resulted in a better understanding of cancer at the genetic, molecular, and cellular levels, allowing for novel therapeutic targets and procedures [11]. Cancer treatments include chemotherapy, surgery and radiation therapy [12,13]. Another way of treating cancer is by using hyperthermia, which involves heating the tumor region without damaging the normal cells [14]. Using hyperthermia, the temperature at the tumor location is raised to a certain level to kill the cancer cells [15]. The targeted tumor sites must attain a temperature in the range of 42–46 °C for hyperthermia to be therapeutically effective in cancer therapy [16]. Protein denaturation occurs when cells are exposed to this temperature range, resulting in a high fraction of co-aggregated denatured proteins [15]. Furthermore, a high temperature influences the cellular structure function and alters intracellular processes, ultimately leading to cancer mortality [16,17], as illustrated in Figure 1. Figure 1. Mechanism of hyperthermia therapy of cancer cells. Figure 2 shows the following two examples of hyperthermia-mediated cancer treatment: photothermal therapy (PTT) and magnetic hyperthermia therapy (MHT) [18–21]. Insufficient blood flow, nourishment and oxygen supply are established inside the blood vessels in the tumor environment, due to the rapid temperature change in malignant cells; however, tumors are more resistant to temperature fluctuations [18,19]. Hyperthermia has been used in cancer therapy to improve the therapeutic effectiveness. Whole-body hyperthermia, local hyperthermia and regional hyperthermia are all possible through hyperthermia therapy [7]. Whole-body hyperthermia involves heating through an external heat source, such as microwaves or radio frequencies, which may have unhealthy side effects because of non-selective heating through non-selective process [21]. Regional hyperthermia heats a large area of cells, such as a body cavity, a limb or an organ. The regional perfusion technique, or the continuous hyperthermic peritoneal perfusion (CHPP) strategy, can be used to apply regional hyperthermia. Drawing blood from the patient’s body, heating it, and then pumping it back into the damaged organ constitutes the regional perfusion approach. The CHPP technique is used to treat peritoneal cancers, including primary peritoneal mesothelioma and stomach cancer [22]. If present in small areas, the cancer cells are killed by introducing heat carriers (Fe, Co, Ni, Ag, Au, etc.) into the body. This is referred to as local body hyperthermia, where nanoparticles are used as heat susceptors [23]. Figure 1. Mechanism of hyperthermia therapy of cancer cells. Figure 2shows the following two examples of hyperthermia-mediated cancer treatment: photothermal therapy (PTT) and magnetic hyperthermia therapy (MHT) [ 18 – 21 ]. Insufficient blood flow, nourishment and oxygen supply are established inside the blood vessels in the tumor environment, due to the rapid temperature change in malignant cells; however, tumors are more resistant to temperature fluctuations [ 18 , 19 ]. Hyperthermia has been used in cancer therapy to improve the therapeutic effectiveness. Whole-body hyperthermia, local hyperthermia and regional hyperthermia are all possible through hyperthermia therapy [ 7 ]. Whole-body hyperthermia involves heating through an external heat source, such as microwaves or radio frequencies, which may have unhealthy side effects because of non-selective heating through non-selective process [ 21 ]. Regional hyperthermia heats a large area of cells, such as a body cavity, a limb or an organ. The regional perfusion technique, or the continuous hyperthermic peritoneal perfusion (CHPP) strategy, can be used to apply regional hyperthermia. Drawing blood from the patient’s body, heating it, and then pumping it back into the damaged organ constitutes the regional perfusion approach. The CHPP technique is used to treat peritoneal cancers, including primary peritoneal mesothelioma and stomach cancer [ 22 ]. If present in small areas, the cancer cells are killed by introducing heat carriers (Fe, Co, Ni, Ag, Au, etc.) into the body. This is referred to as local body hyperthermia, where nanoparticles are used as heat susceptors [23]. Materials 2022,15, 1786 3 of 13 Materials 2022, 15, x FOR PEER REVIEW 3 of 14 Figure 2. Schematic of hyperthermia utilizing nanoparticles: (a) photothermal treatment with laser irradiation to the cancer site, and (b) magnetic hyperthermia therapy with magnetic field exposure [21]. The nanoparticles (NPs) can be heated through high-intensity focused ultrasounds, magnetic hyperthermia, microwave/radio frequencies and plasmonic photo-thermal therapy. The heating method is selected by considering the given conditions. Since all hyperthermia treatment techniques have drawbacks, there is a need to develop new cancer treatment approaches that are more effective and less harmful to the healthy cells. When a magnetic field, in alternating mode, is applied to the NPs in magnetic hyperthermia, electromagnetic energy is transformed into heat. The heating of NPs is possible due to three mechanisms, namely, frictional heating in an anisotropic magnetic particle, eddy currents that have high electrical conductivity [24], and hysteresis [25]. The effects of eddy currents are generally negligible, due to the small size of NPs. The heating process depends upon the shape, nature, size and thermal characteristics of the tissue, and the frequency and magnitude of the magnetic field. In designing NPs for practical use, the major issue is the biocompatibilization and functionalization of the surface coating materials, along with the selection of a suitable core. Noble metals, long-chain organic ligands, inorganic polymers and organic polymers are some examples of surface coatings. The importance of such coatings is that these are used to anchor the functional groups, for example, biomarkers, peptides and antibodies. These coatings are also used to prevent clustering of NPs, due to interactions between particles, which eventually provides stability to the colloidal solutions prepared with NPs. Furthermore, these coatings enhance the biocompatibility of NPs for preventing the leakage of toxic ions from the magnetic core into the biological system [26,27]. Due to their strong magnetic characteristics, outstanding biocompatibility and low cost, magnetic iron oxide nanoparticles are the most extensively employed magnetic nanomaterials. Iron oxide nanoparticles are being researched for their applications in a variety of fields, including biomedical, environmental science, sensing, electronic devices and energy storage [28,29]. Several researcher groups have reported the potential uses of iron oxide nanoparticles in hyperthermia. Table 1 provides a comparison of past studies on different nanostructured materials for the treatment of cancer. Figure 2. Schematic of hyperthermia utilizing nanoparticles: ( a ) photothermal treatment with laser irradiation to the cancer site, and ( b ) magnetic hyperthermia therapy with magnetic field exposure [21]. The nanoparticles (NPs) can be heated through high-intensity focused ultrasounds, magnetic hyperthermia, microwave/radio frequencies and plasmonic photo-thermal therapy. The heating method is selected by considering the given conditions. Since all hyperthermia treatment techniques have drawbacks, there is a need to develop new cancer treatment approaches that are more effective and less harmful to the healthy cells. When a magnetic field, in alternating mode, is applied to the NPs in magnetic hyperthermia, electromagnetic energy is transformed into heat. The heating of NPs is possible due to three mechanisms, namely, frictional heating in an anisotropic magnetic particle, eddy currents that have high electrical conductivity [ 24 ], and hysteresis [ 25 ]. The effects of eddy currents are generally negligible, due to the small size of NPs. The heating process depends upon the shape, nature, size and thermal characteristics of the tissue, and the frequency and magnitude of the magnetic field. In designing NPs for practical use, the major issue is the biocompatibilization and functionalization of the surface coating materials, along with the selection of a suitable core. Noble metals, long-chain organic ligands, inorganic polymers and organic polymers are some examples of surface coatings. The importance of such coatings is that these are used to anchor the functional groups, for example, biomarkers, peptides and antibodies. These coatings are also used to prevent clustering of NPs, due to interactions between particles, which eventually provides stability to the colloidal solutions prepared with NPs. Furthermore, these coatings enhance the biocompatibility of NPs for preventing the leakage of toxic ions from the magnetic core into the biological system [26,27]. Due to their strong magnetic characteristics, outstanding biocompatibility and low cost, magnetic iron oxide nanoparticles are the most extensively employed magnetic nanomaterials. Iron oxide nanoparticles are being researched for their applications in a variety of fields, including biomedical, environmental science, sensing, electronic devices and energy storage [ 28 , 29 ]. Several researcher groups have reported the potential uses of iron oxide nanoparticles in hyperthermia. Table 1provides a comparison of past studies on different nanostructured materials for the treatment of cancer. Materials 2022,15, 1786 4 of 13 Table 1. A comparison of the past studies utilizing nanostructured materials for the treatment of cancer. Study Type Particle Type and Size Coating Results and Drawbacks Reference Magnetic Hyperthermia (Experimental) 3–4 nm of Fe2O3Dextran Post-treatment tumor development is slowed. Induces inadvertent MF infiltration. Jordan et al. (1997) [30] Photo-thermal hyperthermia (theoretical) Silica nano shell of 20 nm radius Au coated Suitable to obtain thermal regime. There is inconsideration of the role of blood perfusion rate and metabolic heat. Dombrovsky et al. (2011) [31] Theoretical model Fe3O4of 0.9 mm The appropriate dose of nanoparticles for hyperthermia. The blood vessels close to the tumor reduce the temperature achieved in the tissue, whereas there are no large blood vessels around the tumor. Pavel M. et al. (2009) [32] Theoretical + experimental Fe3O4of 10 nm The highest achievable temperature depends on surface-to-volume ratio. It is possible to describe patient-specific models. Henrich F., Rahn H. and Odenbachs. (2015) [33] Experimental Photothermal hyperthermia Lipos AuNPs of 5–8 nm Gold coated Hybrid NPs (biodegradable) for treatment. Rengan et al. (2015) [34] Theoretical Fe3O4of 18 nm They found that steady stable temperature achieved after 200 s was at the center of the tumor. The temperature curve declines with increasing distance from the centre at different exposure times. Wu, L., Cheng, J., Liu, W., & Chen, X. (2015) [35] Silver nanoparticles (AgNPs) stand out among the metallic nanomaterials because of their uses in health-care products, textiles, consumer products, medical devices and biosensing, due to their unique physical and chemical characteristics. AgNPs show good thermal, optical and electrical properties and activity against fungus, bacteria, and even viruses [ 36 – 38 ]. Recently, AgNPs have piqued the interest of researchers in nanomedicine, since multiple studies have shown that these NPs can generate antitumoral effects in in vitro and in vivo tumor models, potentially benefiting a variety of oncotherapy modalities and diagnostic tools [ 25 – 28 ]. In addition to antibacterial properties, AgNPs have unique cytotoxic effects against mammalian cells, making silver-based nanoparticles potentially useful in tumor treatment. The therapeutic efficacy of AgNPs is based on their distinct way of inducing cell death in mammalian cells. Despite the physical and chemical features, such as size, shape and heterogeneity of the capping material, their mechanism of action to promote cancer cell death is quite deterministic [ 39 ]. AgNPs are assembled in endosomes after being taken up by endocytosis-related processes, and the organelles are subsequently guided to undergo lysosomal fusion. The lysosomal acidic environment causes an increase in the release of silver ions from AgNPs, which then unbalances cellular homeostasis and leads to apoptotic cell death, depending on the biological aspect of the targeted cell [40]. Crystalline and/or amorphous core–shell structures have distinct properties, along with their medical applications. There is gradual deposition, according to the Stöber method, of small gold colloids onto the surface of cores [ 41 ]. The gold nanoparticles then grow in number and consolidate by forming an isolated island, to create an imperfect, uneven coating, which eventually forms a continuous complete shell that surrounds the core. Wu et al. [ 35 ] conducted a theoretical study on magnetite nanoparticles to analyze the temperature distribution in the tumor. They concluded that the constant temperature was obtained after 200 s, which is a very long time. The motivation behind our work was to achieve the desired temperature in a shorter time period, by using AgNPs. A comparative study on the heat generation, using a single silver nanoparticle, nanosphere, nanorod and Materials 2022,15, 1786 5 of 13 nano-ellipsoid, was conducted in the reported work. The volume of these shapes is set as the following: Vsphere ∼ =Vrod ∼ =Vellipsoid i . e ., Vsphere =4 3πr3= 33510.32 nm3 , Vellipsoid = 4 3πabc = 33401.41 nm3 , Vrod =4 3πr2(r+h)= 33324.96 nm3 . The thermal effect of heat propagation in the tumor cell, as well as the spatial–temporal distributions of temperature during the treatment, are described in this model. Furthermore, the effect of different coating materials, such as gold (Au) and polymer (PEG), is analyzed. The thickness of the shell is discovered to be an essential component that influences the thermal response through the thermal characteristics of the materials utilized in the treatment system. Finally, varying quantities of NPs are attached to the core surface to simulate an incomplete coating surface. The heat transfer module in COMSOL Multiphysics is used to simulate the heating process of the nanostructures, using finite element simulations [42]. 2. Methods A 0.5 µm spherical domain of tissue enclosed the nanoparticle. COMSOL Multiphysics was used to identify the spatial and temporal distribution of temperature in this domain. The major goal of this research was to examine the thermal response of living tissue when the heating source (silver nanostructure) was of various forms, i.e., nanosphere, nanorod, nano-ellipsoid and complex core–shell structure. Due to differences in the structure and function of tissues, balancing thermal energy in different biological tissues is a difficult process. Energy balancing is affected by the relative relevance of a heat transfer technique, the relevant time scale of the deposited energy, and changes in the boundary and initial conditions [ 43 ]. In order to simulate the fundamental features of the thermal state of the organism (or its components), as well as the impacts of the boundary and initial conditions, simplifying assumptions are frequently required. The application of the energy conservation law to a control volume is typically the first step in the design of such models. Qgain =Qstorage +Qloss +W(1) where Qgain is the heat gained by the tissue, Qstorage is the heat stored in the tissue, Qloss is heat loss by conduction and W is work conducted by the tissue. Heat loss by conduction or heat exchange with flowing fluids, as well as the work conducted by the tissue, balance out the heat storage. Heat gain and heat energy storage due to heat generated by unit tissue segment q(r,t) can be represented as an integral over the control volume, as follows: Qgain =Zq(r,t)dV (2) Qstorage =Zρc∂ ∂tT(r,t)dV (3) When no inertia is present, the heat conduction process is often represented using the Fourier law for bio-heat transfer problems. Qconduction =−kt ∆Lt At(T1−T2),T1>T2(4) where kt is the thermal conductivity of tissue, with dimensions [W/m*K]. It is based on the biological materials’ microscopic structures. The differential form of the conduction component of heat flux is as follows: qcond.=−kt∂ ∂xTt(x,t)(5) The convection term is expressed as follows: qconv.=hTw−Tf luid(6) Materials 2022,15, 1786 6 of 13 where his known as the coefficient of heat transfer, and it depends upon fluid velocity. The heat transfer due to fluid flow should contribute to blood perfusion and flow distribution in the biological tissue, due to the diversity of living structures. The flow is proportional to the difference in artery (Tart or Ta) and venous (Tven) blood temperature. qb=ωbρbcb(Tart −Tven)(7) It is possible that blood flows extremely slowly in capillary beds to achieve perfect thermal equilibrium with the tissues. Therefore, the heat flux related to perfusion (Equation (7)) can be approximated as follows: qb=ωbρbcb(Tart −Tt)(8) Under these conditions, the transfer of energy by the stream of blood throughout the entire volume is as follows: Qb=Zωbρbcb[Tart(r,t)−Tt(r,t)]dV (9) Even though it is commonly observed in model formulations, the result of Equation (9) is not unconditional; in larger arteries, it would be invalid, owing to intense blood mixing. By analyzing the various heat transmission modes, it is possible to create a thermal energy balance by ignoring the work conducted by the tissue across an arbitrary volume element, as follows: Zρtct∂ ∂tTt(r,t)dV =Z−kt∇Tt(r,t)dV +Zωbρbcb[Tart(r,t)−Tt(r,t)]dV +ZQm(r,t)dV (10) We obtain, in a one-dimensional scenario with a homogenously distributed source of metabolic heat, Qm. ρtct∂ ∂tTt(x,t)=kt∂2 ∂x2Tt(x,t)+ωbρbcb[Tart(r,t)−Tt(r,t)] +Qm(11) This is the most common form of the heat transfer equation for living organisms, also known as the bio-heat transfer equation. Through scaling and dimensional analysis, we seek to examine its components and assumptions, which may alter depending on the conditions imposed by the modelled item, as well as the contribution of the transport mechanism involved. The temperature distribution in the tumor cell can be modeled using the Fourier heat equation, as follows: ρCp∂T ∂t+ρCpu.∇T+∇.q=Q+Qbio (12) where, q=−k∇T(13) Qbio =ρbCp,bwb(Tb−T)+Qmet (14) which results in the following: δtzρCp∂T ∂t+∇.(−k∇T)=ρbCbωb(Tb−T)+Qmet +Qext (15) where Cp is tissue-specific heat capacity at constant pressure, ρ is tissue density, k is cell thermal conductivity, ρb is blood density, which is 1000 kg/m3 , Cb is specific heat of blood, which is 4180 J/(kg ∗K) , ωb is blood perfusion rate of value 0.0064 1 /s , Tb is arterial blood temperature, which is approximately equal to core body temperature, i.e., 37 °C [ 44 ], Tis the local temperature, and Q= 10 16 W/m3 is heat dissipated by the nanoparticles in the volume of the cell [ 45 ], Qmet = 5790 W/m3 for cancerous cells [ 46 – 50 ] and Qext is the heat Materials 2022,15, 1786 7 of 13 generated by loss powers. Initially, the temperature of the tissue was taken as normal for the human body (Ti=37 °C). Ti(r, 0)=T0i,∂Ti(r,0) ∂t=0, while i=1, 2. The geometry, which consisted of the tissue and nanoparticle, was discretized on all domains by free tetrahedral elements, as shown in Figure 3. There were 60,246 meshing domain elements after meshing the geometry, while the elemental size was chosen as finer. The following boundary conditions were applied to complete the procedure: 1. The tumor cell receives the heat flux from the particle in its entirety, i.e., continuity. 2. The temperature of the outer surface of the tissue is maintained at T=T0=37 °C. Materials 2022, 15, x FOR PEER REVIEW 7 of 14 where 𝐶 is tissue-specific heat capacity at constant pressure, 𝜌 is tissue density, 𝑘 is cell thermal conductivity, 𝜌 is blood density, which is 1000 kg/m, 𝐶 is specific heat of blood, which is 4180 𝐽/(kg ∗ 𝐾), 𝜔 is blood perfusion rate of value 0.0064 1/𝑠, 𝑇 is arterial blood temperature, which is approximately equal to core body temperature, i.e., 37 ℃ [44], T is the local temperature, and 𝑄 =10 (W/m) is heat dissipated by the nanoparticles in the volume of the cell [45], 𝑄 =5790 W/m for cancerous cells [46–50] and 𝑄 is the heat generated by loss powers. Initially, the temperature of the tissue was taken as normal for the human body (𝑇=37 ℃). 𝑇(𝑟,0)=𝑇,(,)  =0,while 𝑖=1,2. The geometry, which consisted of the tissue and nanoparticle, was discretized on all domains by free tetrahedral elements, as shown in Figure 3. There were 60,246 meshing domain elements after meshing the geometry, while the elemental size was chosen as finer. The following boundary conditions were applied to complete the procedure: 1. The tumor cell receives the heat flux from the particle in its entirety, i.e., continuity. 2. The temperature of the outer surface of the tissue is maintained at 𝑇=𝑇=37 ℃. Figure 3. The discretized geometry of tissue and nanoparticle. Silver nanostructures of different shapes (sphere, rod and ellipsoid) were analyzed, as shown in Figure 4. Firstly, the spherical form of silver, with a radius of 20 nm, was analyzed. The volume of the particle was kept the same as the other shapes by computing the dimensions for the rod and ellipsoid. For the nanorod, the length of the cylinder was taken as 𝐿 =73 nm, radius for the hemispherical caps 𝑅 =𝑅 =11 nm, and the dimensions for the ellipsoid were taken as 12 −15 − 44.3 nm. To study the coating effect, the thickness of the Au and PEG polymer coating was set as 5, 10, 20, 30 and 40 nm. The core–shell structure, with a radius of 20 nm, was chosen for the purpose of analyzing the thermal evolution of the proposed structure. Three different coating surfaces (one spherical and two ellipsoidal) were also simulated on the silver nanoparticle core to compare the thermal responses of the ellipsoidal and spherical surface coatings. A silver nanoparticle, which provided a core with a radius of 20 nm, was coated with gold. The simulation growth of gold nanoparticles on the core of AgNP can be observed in Figure 5. The thermal properties of various materials used in this work are given in Table 2. Figure 3. The discretized geometry of tissue and nanoparticle. Silver nanostructures of different shapes (sphere, rod and ellipsoid) were analyzed, as shown in Figure 4. Firstly, the spherical form of silver, with a radius of 20 nm, was analyzed. The volume of the particle was kept the same as the other shapes by computing the dimensions for the rod and ellipsoid. For the nanorod, the length of the cylinder was taken as Lcyl = 73 nm , radius for the hemispherical caps Rcyl =Rcap = 11 nm , and the dimensions for the ellipsoid were taken as 12 − 15 − 44.3 nm . To study the coating effect, the thickness of the Au and PEG polymer coating was set as 5, 10, 20, 30 and 40 nm. The core–shell structure, with a radius of 20 nm, was chosen for the purpose of analyzing the thermal evolution of the proposed structure. Three different coating surfaces (one spherical and two ellipsoidal) were also simulated on the silver nanoparticle core to compare the thermal responses of the ellipsoidal and spherical surface coatings. A silver nanoparticle, which provided a core with a radius of 20 nm, was coated with gold. The simulation growth of gold nanoparticles on the core of AgNP can be observed in Figure 5. The thermal properties of various materials used in this work are given in Table 2. Materials 2022, 15, x FOR PEER REVIEW 8 of 14 Figure 4. Different shapes used in simulations: (a) nano-ellipsoid, (b) nanorod and (c) nanosphere. Figure 5. The core–shell structure simulated in COMSOL Multiphysics: (a) naked core with radius of 20 nm, with AgNPs with 4 nm radius attached to the core surface, (b) core (0 NP) and (c–e) 10, 40 and 70 nanoparticles attached to core, (f) complete shell (100 NPs). Table 2. Thermal response of source materials used in this work. Thermal Conductivity [W/m*k] Mass Density [kg/m3] Specific Heat Capacity [J/kg*K] Tissue [46] 0.512 1000 3800 Tumor [47] 71 21,500 132 Gold [48] 317 19,300 129 Silver [49] 429 10,500 235 Polymer [48] 0.2 1000 1000 3. Results and Discussion To investigate the potential use of silver nanostructures, with different shapes, in hyperthermia, a nanorod, nano-ellipsoid and nanosphere were placed in a tumor with a spherical shape and a 500 nm radius. Their heating effect was simulated to kill the tumor cells. The first simulation showed the maximum temperature of all geometries in the tissue, when heated using an external source. The heating of the nanoparticles caused the temperature and thermal equilibrium of the tumoral cell to change over time, as observed in Figure 6. The temperatures attained by using various simulated forms differed notably. The variation in temperatures can be explained by the difference in the proportionate surface of each form, given that the thermal response is determined from the heat generation multiplied by the volume and surface (𝑄=10 (W/m). The maximum temperature, attained using the nanosphere, was 43 ℃. Similarly, the nano-ellipsoid and nanorod resulted in maximum attainable temperatures of 43.1 ℃ and 44.3 ℃, respectively. These temperature values were taken from the center of the particles, which then propagated into the surrounding medium, as shown in Figure 7. The thermal field distribution of an ellipse is larger than that of a sphere or a rod. Figure 4. Different shapes used in simulations: (a) nano-ellipsoid, (b) nanorod and (c) nano-sphere. Materials 2022,15, 1786 8 of 13 Materials 2022, 15, x FOR PEER REVIEW 8 of 14 Figure 4. Different shapes used in simulations: (a) nano-ellipsoid, (b) nanorod and (c) nanosphere. Figure 5. The core–shell structure simulated in COMSOL Multiphysics: (a) naked core with radius of 20 nm, with AgNPs with 4 nm radius attached to the core surface, (b) core (0 NP) and (c–e) 10, 40 and 70 nanoparticles attached to core, (f) complete shell (100 NPs). Table 2. Thermal response of source materials used in this work. Thermal Conductivity [W/m*k] Mass Density [kg/m3] Specific Heat Capacity [J/kg*K] Tissue [46] 0.512 1000 3800 Tumor [47] 71 21,500 132 Gold [48] 317 19,300 129 Silver [49] 429 10,500 235 Polymer [48] 0.2 1000 1000 3. Results and Discussion To investigate the potential use of silver nanostructures, with different shapes, in hyperthermia, a nanorod, nano-ellipsoid and nanosphere were placed in a tumor with a spherical shape and a 500 nm radius. Their heating effect was simulated to kill the tumor cells. The first simulation showed the maximum temperature of all geometries in the tissue, when heated using an external source. The heating of the nanoparticles caused the temperature and thermal equilibrium of the tumoral cell to change over time, as observed in Figure 6. The temperatures attained by using various simulated forms differed notably. The variation in temperatures can be explained by the difference in the proportionate surface of each form, given that the thermal response is determined from the heat generation multiplied by the volume and surface (𝑄=10 (W/m). The maximum temperature, attained using the nanosphere, was 43 ℃. Similarly, the nano-ellipsoid and nanorod resulted in maximum attainable temperatures of 43.1 ℃ and 44.3 ℃, respectively. These temperature values were taken from the center of the particles, which then propagated into the surrounding medium, as shown in Figure 7. The thermal field distribution of an ellipse is larger than that of a sphere or a rod. Figure 5. The core–shell structure simulated in COMSOL Multiphysics: ( a ) naked core with radius of 20 nm, with AgNPs with 4 nm radius attached to the core surface, ( b ) core (0 NP) and ( c – e ) 10, 40 and 70 nanoparticles attached to core, (f) complete shell (100 NPs). Table 2. Thermal response of source materials used in this work. Thermal Conductivity [W/m*k] Mass Density [kg/m3] Specific Heat Capacity [J/kg*K] Tissue [46] 0.512 1000 3800 Tumor [47] 71 21,500 132 Gold [48] 317 19,300 129 Silver [49] 429 10,500 235 Polymer [48] 0.2 1000 1000 3. Results and Discussion To investigate the potential use of silver nanostructures, with different shapes, in hyperthermia, a nanorod, nano-ellipsoid and nanosphere were placed in a tumor with a spherical shape and a 500 nm radius. Their heating effect was simulated to kill the tumor cells. The first simulation showed the maximum temperature of all geometries in the tissue, when heated using an external source. The heating of the nanoparticles caused the temperature and thermal equilibrium of the tumoral cell to change over time, as observed in Figure 6. The temperatures attained by using various simulated forms differed notably. The variation in temperatures can be explained by the difference in the proportionate surface of each form, given that the thermal response is determined from the heat generation multiplied by the volume and surface ( Q= 10 8W/m2 . The maximum temperature, attained using the nanosphere, was 43 °C . Similarly, the nano-ellipsoid and nanorod resulted in maximum attainable temperatures of 43.1 °C and 44.3 °C , respectively. These temperature values were taken from the center of the particles, which then propagated into the surrounding medium, as shown in Figure 7. The thermal field distribution of an ellipse is larger than that of a sphere or a rod. Materials 2022, 15, x FOR PEER REVIEW 9 of 14 Figure 6. 2D spatial temperature distributions in tumoral cell using (a) nano-ellipsoid, (b) nanorod and (c) nanosphere. The temperature distributions were obtained after 0.5 µm of heating. A functional surface is provided by the NPs coated with various materials in hyperthermia applications. The influence of the coating surface on thermal dissipation in the surrounding medium was investigated using simulations. The silver magnetic core was used as a heat source, with a polymeror gold-like shell. The thermal conductivity coefficients for the two simulated shell materials resulted in opposing heat dissipation effects in relation to the shell thickness. The thermal conductivity for gold decreased, while it increased for the polymer. Figure 7. The time-dependent and radial–temporal distribution for three different shapes: (a) temperature evolution at the center of tumor (x = 0) and (b) the radial distribution after 3 µs of heating process. The hyperthermia process can be related to the functionality of core–shell structures. To achieve the maximum temperature, a sphere, with a radius of 20 nm, was chosen for the assigned nanoparticles. In this simulation, the core of AgNP was covered by the shell of gold or PEG polymer, and then the temperature induced by the core was analyzed. Different values of shell thickness were considered, and the effect of the temperature induced by the core of the nanoparticle is given in Figure 8. The temperature falls with an increase in the thickness of the Au shell, and rises with an increase in the thickness of the PEG polymer shell. As the thermal conductivity of gold is high, it rapidly transfers the heat to its surroundings. While the polymer, which has low thermal conductivity, preserves more heat inside the particle, which results in a rise in temperature. The effect and influence of the heat directly depend on the conductivity. These findings clearly show the possibility of controlling the temperature by changing the thickness of the shell. In Figure 6. 2D spatial temperature distributions in tumoral cell using ( a ) nano-ellipsoid, ( b ) nanorod and (c) nanosphere. The temperature distributions were obtained after 0.5 µm of heating. Materials 2022,15, 1786 9 of 13 A functional surface is provided by the NPs coated with various materials in hyperthermia applications. The influence of the coating surface on thermal dissipation in the surrounding medium was investigated using simulations. The silver magnetic core was used as a heat source, with a polymeror gold-like shell. The thermal conductivity coefficients for the two simulated shell materials resulted in opposing heat dissipation effects in relation to the shell thickness. The thermal conductivity for gold decreased, while it increased for the polymer. Materials 2022, 15, x FOR PEER REVIEW 9 of 14 Figure 6. 2D spatial temperature distributions in tumoral cell using (a) nano-ellipsoid, (b) nanorod and (c) nanosphere. The temperature distributions were obtained after 0.5 µm of heating. A functional surface is provided by the NPs coated with various materials in hyperthermia applications. The influence of the coating surface on thermal dissipation in the surrounding medium was investigated using simulations. The silver magnetic core was used as a heat source, with a polymeror gold-like shell. The thermal conductivity coefficients for the two simulated shell materials resulted in opposing heat dissipation effects in relation to the shell thickness. The thermal conductivity for gold decreased, while it increased for the polymer. Figure 7. The time-dependent and radial–temporal distribution for three different shapes: (a) temperature evolution at the center of tumor (x = 0) and (b) the radial distribution after 3 µs of heating process. The hyperthermia process can be related to the functionality of core–shell structures. To achieve the maximum temperature, a sphere, with a radius of 20 nm, was chosen for the assigned nanoparticles. In this simulation, the core of AgNP was covered by the shell of gold or PEG polymer, and then the temperature induced by the core was analyzed. Different values of shell thickness were considered, and the effect of the temperature induced by the core of the nanoparticle is given in Figure 8. The temperature falls with an increase in the thickness of the Au shell, and rises with an increase in the thickness of the PEG polymer shell. As the thermal conductivity of gold is high, it rapidly transfers the heat to its surroundings. While the polymer, which has low thermal conductivity, preserves more heat inside the particle, which results in a rise in temperature. The effect and influence of the heat directly depend on the conductivity. These findings clearly show the possibility of controlling the temperature by changing the thickness of the shell. In Figure 7. The time-dependent and radial–temporal distribution for three different shapes: ( a ) temperature evolution at the center of tumor (x = 0) and ( b ) the radial distribution after 3 µ s of heating process. The hyperthermia process can be related to the functionality of core–shell structures. To achieve the maximum temperature, a sphere, with a radius of 20 nm, was chosen for the assigned nanoparticles. In this simulation, the core of AgNP was covered by the shell of gold or PEG polymer, and then the temperature induced by the core was analyzed. Different values of shell thickness were considered, and the effect of the temperature induced by the core of the nanoparticle is given in Figure 8. The temperature falls with an increase in the thickness of the Au shell, and rises with an increase in the thickness of the PEG polymer shell. As the thermal conductivity of gold is high, it rapidly transfers the heat to its surroundings. While the polymer, which has low thermal conductivity, preserves more heat inside the particle, which results in a rise in temperature. The effect and influence of the heat directly depend on the conductivity. These findings clearly show the possibility of controlling the temperature by changing the thickness of the shell. In addition, the material and shell thickness determine the required temperature for hyperthermia, which is dependent on the specific location of use in the human body. Materials 2022, 15, x FOR PEER REVIEW 10 of 14 addition, the material and shell thickness determine the required temperature for hyperthermia, which is dependent on the specific location of use in the human body. Figure 8. (a) The core–shell structure, (b) maximum temperature obtained by coating of gold and PEG polymer shell of thickness 5, 10, 20, 30 and 40 nm. The shape of the coating surface plays an important role in the functionality of nanostructures. To study the effect of the coating surface, two forms of coating surfaces, namely, spherical (30 nm radius) and two ellipsoids ( 𝑒𝑙𝑙𝑖𝑝𝑠𝑜𝑖𝑑:25−25− 43.2 nm,𝑒𝑙𝑙𝑖𝑝𝑠𝑜𝑖𝑑:22− 25 − 49 nm), were considered. From Figure 9, it is clear that there is no major difference in the temperature of the coating surfaces, which shows that anisotropy of coating surfaces is not important for the hyperthermia process, since gold nanoparticles grow and form an isolated island, to create an incomplete irregular coating, which is then transformed into a complete shell to cover the core. For this purpose, small AuNPs, with a 4 nm radius, were attached to the core surface of AgNPs, with a radius of 20 nm, which was then embedded into tissue with a 0.5 µm radius. The volume coverage ratio of AuNPs, compared to the volume of the full shell, was determined to describe the temperature profile for incompletely covered nanoparticles with varying amounts of nanoparticles attached, as shown in Figure 10. It is observed that the maximum temperature at the center of the nanoparticle was 42.3 ℃ for the naked core, while there was a gradual decrease in temperature with an increase in the amount of surface coating. The minimum temperature at the center of the nanoparticle was 39.9 ℃ when a complete shell was formed. Figure 9. Radial temperature distributions of different coated shapes. Figure 8. ( a ) The core–shell structure, ( b ) maximum temperature obtained by coating of gold and PEG polymer shell of thickness 5, 10, 20, 30 and 40 nm. The shape of the coating surface plays an important role in the functionality of nanostructures. To study the effect of the coating surface, two forms of coating surfaces, namely, spherical (30 nm radius) and two ellipsoids ( ellipsoid1 : 25 − 25 − 43.2 nm , ellipsoid2 :