Konstantin
Tamarov
*ab,
Maxim
Gongalsky
ac,
Liubov
Osminkina
ae,
Yuanhui
Huang
cd,
Murad
Omar
cd,
Valery
Yakunin
a,
Vasilis
Ntziachristos
cd,
Daniel
Razansky
cd and
Victor
Timoshenko
ae
aFaculty of Physics, M. V. Lomonosov Moscow State University, 119991 Moscow, Russia. E-mail: konstantin.tamarov@uef.fi
bDepartment of Applied Physics, University of Eastern Finland, 70211 Kuopio, Finland
cInstitute of Biological and Medical Imaging, Helmholtz Zentrum München, 85764 Neuherberg, Germany
dChair of Biological Imaging, Technische Universität München, 80333 Munich, Germany
eNational Research Nuclear University “MEPhI”, International Laboratory “Bionanophotonics”, 115409 Moscow, Russia
First published on 24th March 2017
The development of suitable contrast agents can significantly enhance the efficiency of modern imaging and treatment techniques, such as thermoacoustic (TA) tomography and radio-frequency (RF) hyperthermia of cancer. Here, we examine the heating of aqueous suspensions of silicon (Si) and gold (Au) nanoparticles (NPs) under RF irradiation in the MHz frequency range. The heating rate of aqueous suspensions of Si NPs exhibited non-monotonic dependency on the electrical conductivity of the suspension. The experimental results were explained by the mathematical model considering oscillating solvated ions as the main source of Joule heating. These ions could be the product of the dissolution of Si NPs or organic coating of Au NPs. Thus, the ions governed the conductivity of the suspensions, which in turn governs both the heating rate and the near-field RF TA response. The model predicted the contrast in different tissues taking into account both Joule heating and dielectric losses.
Gold nanoparticles (Au NPs) have attracted attention as possible RF-absorbing agents in the last decade,8–10 and they were successfully used for non-invasive treatment of cancer.11 However, the underlying mechanisms of RF heating of Au NPs remain controversial: some groups have reported the absence of heating for uncoated Au NPs in deionized water.12,13 A number of explanations were proposed to describe the Au NP hyperthermia effect, including Joule heating,14 classical Mie's theory,15 quantum mechanisms,15 electrophoresis16 and polarization.17 Also, Si-based NPs have been observed to provide strong temperature increases of suspensions and then suppression of tumor growth in mice models under RF irradiation at 27.12 MHz.18 There, the heating effect was qualitatively explained by taking into account the ionic currents on the surface of NPs.
In the present work, we investigated the mechanisms responsible for the heating of aqueous suspensions of solid-state NPs, i.e. porous silicon (PSi) NPs, ablated Si (ASi) NPs, Si nanowires (NWs), and Au NPs as a result of their interaction with RF fields of various frequencies. Furthermore, Si NPs and NWs were assessed as possible TA contrast agents. Finally, we proposed a model considering collective oscillations of solvated ions (e.g. protons) in a viscous liquid medium as a key mechanism of heating. In order to verify the model, we additionally carried out experiments with electrolyte solutions of different pH values and conductivities. Considering the conductivity of the electrolytes allowed the calculation of the contrast ratio between the medium containing NPs and the reference medium, i.e. deionized water and biological tissues. The obtained results are discussed from the viewpoint of application as RF ablation sensitizers and TA contrast agents.
Si NWs were fabricated by metal (silver)-assisted chemical etching (MACE) of Si wafers of specific resistivity ρ = 1–10 mΩ cm and 10 Ω cm.22 Firstly, to form silver NPs on the surface of the Si wafers, they were immersed into a 0.03 M aqueous solution of silver nitrate (AgNO3) and 5 M of HF in a volume ratio of 1:1. Secondly, Ag-coated Si wafers were placed into the solution of 5 M HF and 30% H2O2 in the volume ratio 10:1 for 1 hour. Finally, the layer of Si NWs grown on the substrate was rinsed several times and dried in the air at room temperature. Si NWs were then detached from the substrate by sonication in an ultrasound bath (Elmasonic S10H, Elma Schmidbauer, Germany) in distilled water for 3 hours. Note that MACE of highly doped (HD) wafers (ρ = 1–10 mΩ cm) led to the formation of porous Si NWs, whereas etching of lowly doped (LD) wafers (ρ = 1–10 Ω cm) resulted in non-porous Si NWs.
To prepare NPs by laser ablation, a Si or Au target was placed at the bottom of a glass vessel filled with 20 mL of deionized water. Radiation from a Yb:KGW femtosecond laser (Amplitude Systems, Pessac, France, 1025 nm, 480 fs, 500 mJ, 1–5 kHz) was focused with the help of a 750 mm lens onto the target surface to provide the ablation of the target.23–25 The target was moved at a scanning velocity of 0.35 mm s−1 in the focusing plane to obtain identical surface conditions during the laser ablation, while the thickness of the water layer above the target was about 1 cm.
The RF heating experiments at 27.12 MHz and 40.68 MHz were performed by applying standard medical apparatus for physiotherapy with maximum powers of 66 W and 60 W, respectively. Each apparatus consists of a power source and a pair of electrodes of diameter 40 mm with a distance of 40 mm between them. The heating systems for frequencies of 1.2 MHz and 4.1 MHz consisted of a generator for the corresponding frequency, an amplifier, and two electrodes. For 1.2 MHz the same electrodes as described previously were used, while for 4.1 MHz two electrodes of diameter 78 mm with a separation of 35 mm were used. The estimated field strengths at 1.2 MHz, 4.1 MHz, 27.12 MHz and 40.68 MHz were equal to 20 V cm−1, 90 V cm−1, 90 V cm−1 and 30 V cm−1, respectively, which were similar to the values measured using an antenna connected to an oscilloscope. The temperature of the aqueous suspensions of Si NPs and NWs and saline solutions was measured using an infra-red thermometer (A&D DT–633, A&D, Japan) after switching off the power sources. To investigate the heating process, 8 ml glass cuvettes were filled with the solution and placed in the center, so there always was an air gap between the electrodes and the glass cuvette. The electrical conductivity of the solutions was measured by using an HP 4192A LF impedance analyzer.
A hand-made near-field RF TA (NRT) imaging setup consisted of a custom-made high-voltage impulse generator providing 50 ns (FWHM) impulses carrying 50 mJ of energy. The emission spectrum is presented in the ESI† (Fig. S1). The induced TA response was detected using a PZT ultrasound transducer with a central frequency of 3.5 MHz, which was attached to automatic translational stages for 2D scanning. The suspensions were put into a rigid low-density polyethylene tube with matched acoustic impedance. The content of the tube could be replaced by two syringes connected to its holes. The coil for inducing the thermoacoustic signals, the transducer and the tube were all submerged in deionized water.
Fig. 2 (a) The dependency of heating rates (S(σ)) on conductivity of the electrolyte (black squares), conductivity of the aqueous suspension of Si (blue circles) and Au (green triangles) NPs, and fitting of the electrolyte heating rates according to the model (eqn (8), red line) at 4.1 MHz. (b) The same at 27.12 MHz. |
Similar results were obtained by the near-field RF TA tomography (NRT) experiments. The set of aqueous suspensions was irradiated by RF pulses in the window of 3–12 MHz and thermally induced acoustic waves were detected. Basically, the method provided 2D or 3D-images of an object, governed by its TA response. A typical image of a straw filled with a suspension of NPs is shown in the inset of Fig. 3. Such images were used to obtain average integrated TA responses for each sample. Fig. 3 shows the conductivity dependency of Si NP suspensions and saline solutions. All the points lie perfectly on the red line, which can be considered as the calibration curve. This led to non-linear, but monotonic dependency of the TA response on the sample conductivity. The RF irradiation spectrum and horizontal TA profiles of the tubes are presented in the ESI† (Fig. S1).
Fig. 3 Thermoacoustic response of NP suspensions (blue squares) and saline solutions (black circles) versus their electrolytic conductivity. The red curve shows the interpolation. |
Nevertheless, conductivity remained the major factor influencing TA intensity and, consequently, the efficiency of heating. This also confirmed the idea that both RF heating and TA waves were governed mostly by ionic oscillations in solutions rather than electrons in solid NPs. The latter are not supposed to contribute to RF heating, because their plasmonic frequency is too high, therefore even low-conducting NPs perfectly reflect electromagnetic waves instead of absorbing them.
In the experimental results, it was supposed that the solvated ions were the primary source of temperature increase in the electrolytes exposed to the radiation at RF frequencies below 1 GHz. The oscillating ions converted the energy of the electromagnetic field into heat via viscous friction with surrounding liquid (Fig. 4). The behavior of a single ion can be described by Newton's second law of motion:
mẍ = FE − Fvis, | (1) |
• Positive and negative ions are uniformly arranged in the volume, have charges +q and −q, respectively, and are represented as two areas equal in magnitude and opposite in charge density (blue and brown areas in Fig. 4). The electrolyte is neutral in the absence of external forces. For convenience of description, these uniformly charged areas are considered as undeformable regions, which can move relative to each other. When the external electromagnetic field is applied, the charged areas are displaced by an offset from their equilibrium position in the opposite directions, i.e. displacements of positively and negatively charged areas are x+ and −x−, respectively.
• The movement of the areas instantaneously achieves a quasi-stationary regime because of a strong viscous friction, i.e. the acceleration is neglected and the averaged Brownian random force is equal to zero. In this case, one can use the concept of mobility, μ, which can be expressed as: μ = ẋ/Etotal. μ+ and μ− values are used for positive and negative ions, respectively.
• The relative dielectric permittivity of the electrolytes, ε, equals the relative permittivity of water and does not depend on RF frequency in the range of 1–100 MHz and concentration of solvated ions. More careful calculations can be potentially performed by using the experimental data summarized in ref. 26.
• The skin-effect in the electrolyte is ignored, because its depth is estimated to be in the range of 0.6 m and 160 m in the investigated range of conductivities and RF frequencies. The estimation was done according to the equation for the skin depth: δ = (2/σω0μ0)1/2, where σ is the electrical conductivity of the electrolyte, ω0 is the frequency of RF field and μ0 = 4π × 10−7 H m−1.
The total electric field strength, Etotal, consists of the external RF field, Eext and the internal field of the electrolyte polarization, Eint. For calculation of Eint one can use the formula for two oppositely charged planes Eint = −ΣE/εε0, where ΣE is the surface charge density, ε is the relative permittivity of the electrolytes and ε0 = 8.85 × 10−12 F m−1. ΣE can be considered as ΣE = qn(x+ − x−) and the expression for Etotal will be:
(2) |
Eqn (2) can be rewritten in the following way:
(3) |
Eqn (2) shows, that ẋ+/μ+ = −ẋ−/μ−, and subsequently x+/μ+ = −x−/μ−. Therefore, relative displacement between both areas is the following:
(4) |
Taking into account that the electrical conductivity related to both positive and negative ions σ = qnμ+ + qnμ− and eqn (4), eqn (3) can be written as:
(5) |
By solving eqn (5) the following expression can be obtained:
(6) |
The heating rate produced by viscous friction of ions is obviously proportional to PVI(t) = qnEext(t)(ẋ+(t) − ẋ−(t)). And the average absorbed power density will be:
(7) |
Or in terms of conductivities:
(8) |
Note, that despite there being no direct dependency of PVI(σ) on viscosity, there is indirect dependency, because σ = qnμ and μ decreases with increase in viscosity.
The heating rate is calculated from the heat equation taking into account the negligible temperature exchange between the electrolyte and ambient air:
(9) |
Table 1 demonstrates that the conductivities σ0 corresponding to maximal heating obtained from approximation of the experimental data are in good agreement with the model. However, at 1.2 MHz the shape of the calculated peak differs from the experimental one due to slight deviation from the quasi-stationary regime.
Frequency (MHz) | σ 0 from experiment (S m−1) | σ 0 from model (S m−1) |
---|---|---|
1.2 | 0.0053 ± 0.0023 | 0.0051 |
4.1 | 0.0176 ± 0.0021 | 0.0173 |
27.12 | 0.121 ± 0.011 | 0.115 |
40.68 | 0.183 ± 0.045 | 0.173 |
The developed model shows a linear dependency of the maximum heating rate on frequency: Smax(ω0) = E02 × εε0 × ω0/4cρ. Fitting of the experimental data allowed acquiring values ωmax and heating rate Smax. Fig. 5 shows the frequency dependence of maximum heating normalized by the corresponding field strength. The modeling and the experimental data seem to agree reasonably well with each other confirming the model in the RF region (1–40 MHz).
Fig. 5 Dependence of the maximum heating rate on the external electromagnetic field frequency, Smax(ω0). |
(10) |
(11) |
It is known that the medical devices for physiotherapy are designed to produce the most efficient heating on biological tissues because the conductivities of different parts of the human body lie in the RF range of the corresponding frequencies.30 In order to find the most suitable therapeutic frequency, the contrast η between the heating rates of different types of tissues and suspensions of Si NPs was calculated:
(12) |
Fig. 6 shows the simulated frequency dependence of the contrast for four types of tissues: muscle, bone, fat and blood; the conductivities of the tissues were taken from ref. 30. On the one hand, the contrast is quite low in the range of 2.5–60 MHz for bone, fat and other tissues, which have conductivities in the 0.03–0.06 S m−1 range. On the other hand, the contrast becomes low at frequencies higher than 70 MHz for high conductivity tissues like muscle or blood. There is no measurable contrast between heating of DI water, suspensions of Si NPs and the electrolytes at 2.5 GHz according to the measurements using MW heating equipment.
Fig. 6 Contrast values of the heating rate between a suspension of Si NPs and different tissues: muscle (black line), bone (red line), fat (green line) and blood (blue line) versus RF irradiation frequency, η(ω0). The contrast was calculated by using eqn (12). |
Therefore, it might be concluded that to obtain perfect contrast close to 1 and sufficient hyperthermia effect in the local area, the selected frequency should be about 1 MHz or less. However, some difficulties arise at low frequencies. In particular, at lower frequencies, the conductivity value, where maximum heating is observed, is also low, which leads to a small temperature rise. Additionally, the shielding of the electromagnetic field by surrounding high conductivity tissues may occur, and it can reduce the heating rate of low conductivity materials. Consequently, the power of RF radiation must be increased in order to provide the same temperatures as at high frequencies. This gives rise to another problem, because the creation of such intense electromagnetic fields requires powerful sources and necessary protection to ensure the safety of patients. Thus, we suppose that the most convenient frequency is 70 MHz as it gives the contrast of 25% between Si NPs and a tissue and this value does not depend on the type of biological tissue.
Some early publications attributed the temperature rise of Au NPs suspensions to the absorption of RF radiation by NPs themselves via Joule-like or inductive mechanisms.14,31 But although the high conductivity of gold can yield high heating rates when a current passes through it, in the external electromagnetic field free electrons in metals polarize easily, screen it, and prevent further heating. The negligible heating of Au NPs dispersed in deionized water was observed in our experiments and supports this assumption (Fig. S7, ESI†).
In the second approach, the absorption cross-section of metal NPs suspensions is calculated according to Mie scattering theory.15,32 The authors reported that the negligible absorption of RF radiation by Au NPs is insufficient to heat the entire suspension, whereas the heating is more significant when the NPs are dispersed in a conductive medium; and this conclusion agrees with the developed model.
Thirdly, electrophoretic motion can be considered as a heat source because, typically, the commercially available Au NPs and Si NPs have significant net negative charge at neutral pH values.33,34 Under RF irradiation, charged particles can follow alternating electromagnetic fields. The imaginary component corresponding to power loss during electrophoresis was theoretically considered by Sassaroli32 and experimentally investigated by Corr.16 Although in the latter work Au NPs with diameters less than 10 nm had observable heat production at 13.56 MHz, there is still no full understanding of the influence of the citrate coating. The electrophoresis seem to be more pronounced at lower frequencies, but the calculations suggest that this effect is not sufficient enough to explain the heating rates observed in our experiments.
Fourthly, the presence of a NP in the electrolyte can influence its polarization in the alternating RF field.17 There, the authors demonstrated that in the presence of the charged particle, the local conductivity near the particle is different from the conductivity of the bulk electrolyte due to the formation of an electrical double layer.35 However, although the particle did alter the heating, the high temperature rise was proposed at low conductivity values, which does not correspond to our data.
In the present work, it is argued if solid-state NPs are the main absorbers of RF radiation. Indeed, in the experiments, where significant heating was observed, citrate-coated Au NPs were used.8–11,14,16 On the contrary, several authors reported negligible heating rates for suspensions of non-functionalized Au NPs suspended in pure deionized water.12,13,36 Thus, the citrate coating dissolves and produces ions, which start to play a major role in the heating of suspensions.16 This conclusion is also supported by the work of Ogunlade and Beard,37 who tested gadolinium, iron oxide NPs and single-wall carbon nanotubes as potential TA CAs and demonstrated, that the NP coating was the main source of the contrast.
In order to exclude any influence of surface chemistry, which is usually fairly complicated for particles obtained by wet chemical approaches, ultrapure Au NPs were prepared by laser ablation of a gold target in deionized water.38 Hence, the suspension of Au NPs was obtained in low conductive medium, therefore the effect of intrinsic Au NPs heating would be evident, if it existing. However, the green point in Fig. 2 and the comparison with deionized water (Fig. S7, ESI†) shows that the heating of the Au NP suspension can be completely described by the proposed model. Thus, the remarkable heating reported in the literature is governed mostly by contamination of the solvent by radicals which moved from surface coating to the solution. Similar results were observed for the TA response of Au NPs, which was also negligible.
In the case of Si NPs, they tend to dissolve in aqueous media forming considerable amounts of silicic acid.39 It leads to subsequent conductivity and heating rate increase of the surrounding liquid. To confirm the impact of the dissolution, the suspensions of Si NPs were centrifuged at 20000g and the supernatant was filtered through a 0.2 μm filter. The heating rates of the supernatant (Fig. S8, ESI†) were similar to that of the suspensions shown in Fig. 1. When the sedimented NPs were resuspended in fresh deionized water, heating similar to water was observed (Fig. S9, ESI†).
Despite the fact that the heating rate mostly depends on the ionic conductivity, Si NPs can be potentially used as a source of local conductivity increase and thus act as contrast agents for both therapy and diagnostics. The heating rate of Si NPs in low conductivity tissues could be much higher than the heating rate of the tissue itself, if appropriate frequency tuning is applied. Porous Si NPs can also be considered as compartments for liquid filling their pores. In this case, diffusion inside Si NPs is changed and intraporous liquid may maintain its conductivity for a long time. This may be also beneficial for the delivery of kosmotropes and triggered drug release techniques, because the surface of Si NPs can be coated by thermosensitive polymers, which change their conformation and open pore gates after reaching a certain temperature.40,41 Note, that the concentration dependence diffusion of liquids inside the pores42 should be taken into account for a correct description of RF heating of porous NPs.
Another promising application of solid-state NPs is NRT.4,7 Indeed, different materials have been tested for that purpose: graphene-based NPs,43,44 iron oxide NPs,45 and a Gd-based medical contrast agent.46 Graphene-based NPs demonstrated at best two-fold enhancement, while Gd NPs and iron oxide amplified the TA signal 3-fold and 6-fold, respectively, due to magnetic heating. However, in these studies MW radiation of 3 GHz was used, which penetrates around 0.7 cm into muscle tissue.45 Since Si NPs also demonstrated efficient TA response, they can be considered as a flexible framework for contrast enhancement in the MHz frequency region. The advantage of Si NPs over inorganic NPs mentioned above is their dissolution, which provides significantly better local heating around themselves. That in turn governs a stronger TA response of Si NPs. Moreover, the porous structure of Si NPs provides additional opportunities for contrast enhancement. TA contrast can be estimated, ηRA, as follows:4
(13) |
Footnote |
† Electronic supplementary information (ESI) available: Includes the spectrum of TA emission, structural characterization of Si NPs, heating rates of NaCl solutions at 27.12 MHz, heating rates of electrolytes and suspensions of Si NPs at 1.2 MHz and 40.68 MHz, and heating rates of aqueous suspensions of Si NPs before and after centrifugation. See DOI: 10.1039/c7cp00728k |
This journal is © the Owner Societies 2017 |