Laser Pulse Duration Optimization for Photothermal Therapy with Gold Nanostars ()
1. Introduction
Gold nanoparticles (GNP) have attracted a large amount of attention in cancer research in recent decades [1] [2] [3] [4] . When stimulated with light radiation at the proper wavelengths, the conduction electrons on the GNP surface display unique behavior of collective oscillation, which is known as the Localized Surface Plasmon Resonance (LSPR) [3] [5] [6] [7] . At the LSPR wavelengths, GNPs interact strongly with the incident light, and absorb and scatter much more efficiently than other nanoparticles of similar sizes [6] [8] [9] [10] . Another important feature of GNPs is the strong correlation between its physical shape and its optical properties. Slight differences in size, shape, geometry and the surrounding environments of GNP shift its optical spectrum dramatically. Controlling the shape and size of the GNP becomes an important handle for researchers to manipulate the optical responses. These two important features make GNP an ideal candidate for photothermal imaging [11] [12] [13] [14] and therapy [15] [16] [17] [18] [19] .
In biological application settings, light therapy can be a particularly powerful tool, because it is effective, accurate, cost-effective, and induces little side effects. However, not all light interacts with the biological tissues in the same manner. Ultraviolet and visible lights can deliver a large amount of energy, but they do not penetrate deep into the tissue. For infrared light beyond 1300 nm, it is strongly absorbed by water, which in turn restricts the light penetration depth as well. There is a range of wavelength, between 800 and 1100 nm, where biological tissues do not absorb much light, which is known as the “biological window”. An ideal therapeutic or imaging agent would be designed to interact strongly with the light in this window, because it would enable the light-nanoparticle interaction to take place deeper inside the tissue and induce fewer side effects.
Because GNPs have enhanced optical response to specific wavelength which can be manipulated through changing its physical shape and geometry, several different GNP structures have been demonstrated to be promising candidates for photo-therapy and imaging in the biological window. Gold nanorods [12] [15] [20] , gold nanoshells [21] [22] and gold nanostars [16] [17] [23] [24] are examples of GNP with LSPR in the near-infrared range. Among these different gold structures, gold nanostars stand out as an ideal candidate. It can be easily synthesized, batch produced and purified [25] [26] . It also possesses different curvatures on different surfaces at the tip, on the arm and at the center of the star [27] . These unique different curvatures offer different functionalities for drug delivery purpose. Photothermal therapy with gold nanostars has been performed on animals, and promising results have been shown [16] . However, there are a few technical difficulties concerning the current medical applications of gold nanostars. One of the unique features about gold nanostars is their sharp corners. It has been well-demonstrated that there is a strong plasmon resonance field enhancement effect at the tip of gold nanostructures [28] [29] . When this strong tip enhancement effect happens in a photothermal therapy setting, it gives rise to the temperature heterogeneity on a single nanostar. In other words, the temperature of the tips of the gold nanostars could potentially be higher than other parts of the gold nanostar. This effect is generally a beneficial one because it gives rise to localized heating effect in a confined area. In the case of photothermal therapy of cancer, gold nanostars are delivered to tumor sites and the significantly over-heated tips of gold nanostars will disrupt the micro-environment of tumor tissues, while the normal tissues with no gold nanostars stay relatively unperturbed.
In order to achieve the maximum amount of temperature difference around gold nanostars when they are in tissue, pulse lasers are usually used because they deliver a large amount of power within a short period of time and let the normal tissue rest for an extended period of time [17] [19] . However, the energy in the laser pulses can be so high that it elevates the temperature of the GNP to above its melting point, which would cause it to melt and undergo morphology changes. Similar effects have been observed and heavily studied on gold nanorods. El-Sayed and coworkers studied the melting of gold nanorods and found out the surface defects are strongly correlated with the melting phenomenon [30] . They further determined the melting energy of gold nanorods [31] . In the case of gold nanostars, the tips are likely to reach the melting point at a lower energy density from the laser than the rest of the nanostar body. The morphology change is especially important for gold nanostars, because the sharp tips of the nanostars result in a large aspect ratio, which pins the LSPR of the gold nanostars in the near infrared range. In order to avoid this shape change, Continuous Wave (CW) lasers have been proposed for photothermal therapies. It is far safer for normal tissues and it is far less likely to melt down the gold nanostars. However, it means that we are not close to exploiting the full potential of the photothermal therapeutic effect of gold nanostars. There is a delicate balance that we need to strike here to achieve the optimal result: the laser pulses should be strong enough to significantly heat up the photothermal agents, while keeping the normal tissues and the shape of the nanostars relative unperturbed. In this manuscript, we aim to explore the optimal conditions for striking this delicate balance with our numerical calculations. We first simulate the interaction between gold nanostars and pulse laser light with different pulse duration and pulse energy levels. Then we calculate the heat transfer process when the absorbed optical power of the gold nanostars is dissipated into its surrounding environment. Temperatures of the tips of the gold nanostars and the average temperature of the structure can be obtained and compared to the melting point of gold. After a series of comparisons, we can find out the optimal pulse duration and energy level for gold nanostars in photothermal therapy settings.
2. Material and Methods
In this manuscript, we use Finite Element Analysis with Comsol Multiphysics 5.3 to simulate the gold nanostar interaction with light and the heat transfer process that follows. First of all, the shape of the gold nanostar was constructed with a center sphere fused together with six cones pointing at two opposite directions of all three axes of the Cartesian coordinates. Gold nanostars with different sizes have been used for photothermal therapy. Vo-Dinh and coworkers studied gold nanostars with 30 nm diameter and 60 nm diameter, and their therapeutic effects. Results showed that 30 nm gold nanostars are more effective [16] . Particles with similar sizes were used in a different study and were proven again to be robust therapeutic agents [17] . Therefore we model the gold nanostars in our study to be similar to the ones reported in literature. The diameter of the center sphere is 20 nm, the diameters of the cones are 12 nm at the axes, and the height of the cones is 30 nm. This shape matches what was reported in literature for gold nanostars used in photothermal therapy. The refractive index of gold reported by Johnson and Christy [7] , and the surrounding environment is set to be water with refractive index of 1.33 + 0*i. The light absorption of water between 500 nm to 1300 nm is ignored.
First, we simulate the interaction between gold nanostar and the light in a range of wavelengths. A gold nanostar is placed at origin of the geometry and a concentric sphere encloses the gold nanostar, which represents the homogeneous aqueous environment that the nanoparticle is suspended in. A layer of Perfectly Matched Layer is wrapped around the sphere. The Electromagnetic Waves, Frequency Domain (emw) physics is selected to simulate the light-matter interaction process. The strength of the electric field is specified at different levels to simulate the different illumination intensities that the nanostar is exposed to under different laser intensities. The polarization of the light is along the z axis, direction of traveling is along the x axis. A scattering boundary is selected to be the outmost surface of the Perfectly Matched Layer. The wavelengths for calculation are from 500 nm to 1400 nm, which is wider than the “optical window” that we are interested in.
Then we simulate the heat transfer process between the gold nanostar and the water environment. In order to simulate the pulsed heating process, we first calculate the absorbed power from gold nanostar under peak radiance of pulsed laser radiation. The absorbed optical power of the gold nanostars is assumed to be the heating power in the heat transfer process. Then in the “Definition” section of the Comsol Multiphysics software, we define an analytical function, which denotes a Gaussian distribution of power over the period of pulse duration, with the highest value of intensity as unity. We then select the Heat Transfer in Fluids (ht) physics and define a heat source as the gold nanostar. The power of the heat source is defined to be the absorbed optical power from the light-particle calculation multiplied by the pulse function that we defined earlier. The heat equation in the heat transfer process is:
(1)
(2)
where
is the fluid density;
is the fluid heat capacity at given pressure; k is the fluid thermal conductivity; u is the fluid velocity field; Q is the heat source. An important advantage of Comsol Multiphysics pertaining to our study is that it is able to couple the absorbed optical power from the Electromagnetic Waves (emw) calculation directly in to the Heat Transfer in Fluid (ht) study in every finite element in three dimensions. In other words, with Comsol Multiphysics, we are able to readily map out the different heating power from different parts of the gold nanostars and their resulting different temperatures. The gold nanostars are modeled as pulsed heating sources with a Gaussian pulse shape and pulse durations of 100 femtoseconds, 100 picoseconds, and 100 nanoseconds respectively. It reaches the maximum heating power, as the absorbed power from the Electromagnetic Waves study (emw∙Qh), at the center of the pulse. The repetition rate of the laser is set to be 10 Hz, which is common for lasers for phototherapies.
It is worth noting that this treatment is a little different from the realistic physical process. In reality, when the nanostar interacts with the pulse laser, the nanostar interacts with the light under different intensities in the pulsing period of the laser, which in turn translates to different heating powers. However, since the absorbed power from nanostar is proportional to the illumination power, and the simulation of light-matter interaction over time is much more difficult because it involves calculations of both frequency-domain and time-domain at the same time, we simply simulate the light-matter interaction under peak radiance and obtained the peak heating power and treated the gold nanostar as a pulsed heating source in the time domain in the heat transfer study. Our treatment does not make any approximation and is effectively equivalent to the physical process that we are interested. This treatment saves a large amount of computational power without making any compromises.
3. Results and Discussion
In Figure 1(b), we can clearly see that this gold nanostar structure has a strong LSPR peak in the near infrared range as we desired. The absorption peak is around 830 nm. In Figure 2, the electric field distribution on and near the gold nanostar are shown for (a) femtosecond, (b) picosecond and (c)nanosecond pulses. It is clearly shown that the electric field is much stronger at the tips of the gold nanostars than other parts no matter what the pulse duration is. This observation agrees with what has been reported before, known as the tip enhancement effect [16] [23] [28] . The degree of heterogeneity is larger for (a) femtosecond laser, compared to (b) picosecond and (c) nanosecond, because the peak radiance of the incident light is far more intense in the case of femtosecond laser.
As a result of the heterogeneity of electric field in space, the temperature profile under the peak radiance intensity is also different at the tips. However, as shown in Figures 2(d)-(f), temperatures at the nanostar tips reaches a significantly higher level under femtosecond laser radiation, while this heterogeneity in
Figure 1. (a) Geometry and meshing of the simulated model; (b) Simulated absorption spectrum of gold nanostar in water.
Figure 2. (a)-(c) are the distribution of electric field in space for femtosecond, picosecond and nanosecond pulses respectively; (d)-(f) are the temperature distribution in space at the peak temperature of gold nanostar for femtosecond, picosecond and nanosecond pulses respectively.
temperature is not visible in the cases of picosecond laser and nanosecond lasers. When paring the electric field distribution graphs (a)-(c), together with the corresponding temperature profile graphs of (d)-(f), it is clear that the distribution heterogeneity in electric field translates to temperature heterogeneity in case of femtosecond laser radiation, but not in the cases of picosecond and nanosecond lasers. This interesting observation can be explained by the heat transfer process.
As we know, the temperature of an object is determined by the amount of heat that this object takes in, the amount of heat that the object gives out and the heat capacity of the object. When the gold nanostar interacts with the electromagnetic field, the electric field is going to be much more intense at the tip of the gold nanostar, because the electric field on the surface of gold nanostar is determined by the collective oscillation of the conduction electrons on the surface, also known as the Surface Plasmon Enhancement (SPR) effect. The surface to volume ratio at the tip of the nanostar is far larger than anywhere else on the body. Therefore the electrons are far more crowded at the tips, which resulted in tip enhancement effect. However, because the surface to volume ratio is high at the tips, it is much easier for the tip to cool down through convective heat transfer with a given amount of heat input. Therefore although it is true that the electric field is always going to be much more intense at the tips, which will result in much more intense heat generation, no matter what pulse duration is, cooling is also much quicker at the tip as well. This competing effect manifests itself clearly in Figure 3.
Tables 1-3 show the highest temperatures that the tips and the body of the gold nanostar reach under different peak radiance and pulse durations. In Table 1, under 100 femtosecond pulse radiation, the difference of the temperature rises between the tip and the body of the gold nanostar is very large. For example, under 663.1 TW/m2 peak radiance, the tip temperature rise can reach as high as 46321 K, while the average temperature of the nanostar body only rises by 1814 K. In Table 2, under 100 picosecond pulse radiation, such difference is still clearly present, but is far less dramatic. In Table 3, under 100 nanosecond pulse
Figure 3. (a)-(c) are exemplary temperature decay curves of gold nanostar in water after heating by femtosecond, picosecond and nanosecond laser pulses.
Table 1. Temperature rise of the tip and the nanostar body under different peak radiance with 100 femtosecond pulse duration. The number marked in yellow are above the temperature rise need to melt gold.
Table 2. Temperature rise of the tip and the nanostar body under different peak radiance with 100 picosecond pulse duration. The number marked in yellow are above the temperature rise need to melt gold.
Table 3. Temperature rise of the tip and the nanostar body under different peak radiance with 100 nanosecond pulse duration. The number marked in yellow are above the temperature rise need to melt gold.
radiation, the temperature rises of the tip of the gold nanostar and the average body almost the same. This comparison is also shown in Figure 3 with the time decay curves.
Figure 4 displays the results from Tables 1-3 and compare them to the melting temperature of gold. In Figure 4(a), it is clearly shown that femtosecond pulses will melt the tips much earlier than it melt the body of the nanostars. Figure 4(b) shows that with picosecond pulses, the tips will still melt significantly earlier than the nanostar body, but it is possible to heat the gold nanostar body up to a relatively high temperature before melting its tips. Figure 4(c) shows that nanosecond pulses do not induce much temperature heterogeneity
Figure 4. (a)-(c) are the peak temperatures of the tip of the gold nanostar versus the average temperature of the entire gold nanostar at different laser peak radiance levels for femtosecond, picosecond and nanosecond pulses. The melting point of gold is specified in the graphs for comparison.
on the gold nanostar. It resembles the behavior that one would expect from a Continuous Wave laser. Therefore, there is no need to widen the pulse duration further down.
Based on these tables and graphs, we can clearly observe that under femtosecond laser pulses, the temperature difference between the tips of the gold nanostars and the gold nanostar body is very significant. Such difference is also present under picosecond pulse radiation, but is very minimal under nanosecond laser pulses. This important finding indicates that femtosecond lasers, which are very popular in laser ablation applications, might not be the best option for photothermal therapy with gold nanostars. Picosecond lasers, on the other hand, can generate enough localized heating effect without easily melting the tips of the gold nanostars. Nanosecond lasers behave very similarly to Continuous Wave laser, which does not induce too much temperature heterogeneity on gold nanostars.
4. Conclusion
In this manuscript, the temperature change of gold nanostars under pulse laser with different durations is studied. The temperature heterogeneity caused by the tip enhancement effect of gold nanostar is carefully investigated with different peak radiance and pulse durations. Based on the results, femtosecond laser heats up the tips of the gold nanostars to a much higher (~20×) temperature rise compared to the gold nanostar body. Picosecond laser creates similar temperature heterogeneity, but the difference is not nearly as dramatic (~2×). Nanosecond laser does not induce much temperature heterogeneity. Therefore, femtosecond laser, which is one of the most popular choices for phototherapy, might not be suitable for photothermal therapy with gold nanostars. Because it is likely to melt the tips of the nanostar, which in turn changes the absorption spectrum entirely, before the body of the nanostar reaches a reasonably high temperature. Picosecond second lasers, on the other hand, induce a manageable temperature heterogeneity which can be an ideal option for this application. Nanosecond lasers behave very similarly to CW lasers in this case, which might not be a good candidate because of the elevated costs compared to CW lasers. In the future, experiments will be conducted to validate this simulation. This study offers important insights into selecting the optimal laser for photothermal therapy, which can potentially be a highly specific and minimally invasive tool for cancer therapy.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.