Theoretical and experimental study of dual-fiber laser ablation for prostate cancer
Theoretical and experimental study of dual- fiber laser ablation for prostate cancer
Xing Wu 0 1
Kangwei ZhangID 1
Yini Chen 0 1
Ren Wang 0 1
Lei Chen 0 1
Aili Zhang 1
Bing Hu 0 1
0 Department of Ultrasound in Medicine, Shanghai Jiao Tong University Affiliated 6th People's Hospital , Shanghai , China , 2 Shanghai Institute of Ultrasound in Medicine , Shanghai , China , 3 School of Biomedical Engineering, Shanghai Jiao Tong University , Shanghai , China
1 Editor: Xiaoming He, University of Maryland , UNITED STATES
Single-fiber laser treatment of the prostate has been widely accepted in the clinic due to its minimal invasiveness and high controllability. However, for large tumors, multiple insertions of the laser probe would be needed to achieve full coverage of the tumor, increasing the complexity of the treatment and occasionally resulting in the incomplete killing of tumor cells due to a mismatch between the planned insertion location and the actual probe insertion location. Treatment with a dual-fiber laser results in greater lesion coverage following a single insertion of the probe, with the lesion coverage being even greater than the sum of the coverage of two sequential insertion of a single-fiber laser probe, potentially reducing treatment time and clinical complications. Both theoretical and experimental analyses have been performed to evaluate the proposed dual-fiber laser treatment. A finite element model was established to simulate the treatment process. The simulation results indicated that there is a clear difference between the ablation coverage created using dual-fiber laser ablation and that created using the superposition of sequential single-fiber laser ablation. In addition, the coverage is dependent on the spacing distance between the two fibers. Both ex vivo and in vivo canine prostate tissues were treated by dual-fiber laser ablation, with lesions analyzed by magnetic resonance imaging (MRI), ultrasound imaging, and pathology. The results demonstrate that dual-fiber laser ablation can markedly increase the range of the ablation zone when compared with single-fiber modes. The safety and feasibility of dual-fiber laser treatment has been confirmed, and a treatment plan using dual-fiber laser ablation has also been proposed.
Data Availability Statement: All relevant data are
within the paper.
Funding: This work was supported by the National
Natural Science Foundation of China
(NSFC51476101, 51521005), the Shanghai
Science and Technology Commission of Shanghai
Municipality (16441903100), the Ministry of
Science and Technology of China
(2016YFC0106200), and a clinical support
department capacity building project
(SHDC22015001). The funders had no role in
Laser ablation, which allows for concentrated energy pulses and precise control over the
thermal field, has been widely accepted in clinical treatments of prostate cancer [
]. Studies on
laser ablation in prostate cancer have included clinical trials for low-grade human prostate
cancer, experimental investigations in the canine prostate, and theoretical models in the rat
]. Among these studies, a single-fiber laser has been used to create ablation lesions
with clear local boundaries around a target by moving the fiber to multiple positions under the
study design, data collection and analysis, decision
to publish, or preparation of the manuscript.
guidance of surgical imaging techniques. Such procedures induce large-scale tissue damage by
creating multiple lesions. However, because the ablation process creates many bubbles,
physicians must wait approximately 15 minutes before re-positioning the laser probe to produce a
second lesion, which significantly extends the duration of surgery. Most importantly, cancer
cells may survive in the margins between adjacent lesions.
By comparison, multiple laser fibers can be used (e.g., a dual-fiber laser source), and due to
the interactions of the two fibers, the resulting lesion volume will be much larger than the pure
superposition of two single-fiber ablations. When planning for serial ablative lesion using
single-fiber probe treatments, the second lesion must be adjusted on the basis of the previous
lesion which is difficult to image using current techniques. In contrast, during dual-fiber laser
treatment, the second laser probe is inserted relative to the first laser probe, which can be
clearly imaged. Thus, for targeted prostate cancer tissues larger than a single-fiber ablative
lesion, the use of dual fibers may greatly decrease the duration of surgery and the risk of
Due to the complex physiological environment in which the prostate is located, careful
pretreatment planning is crucial to the final outcome of prostate cancer treatment. Several
theoretical models have been proposed to simulate laser ablation for better clinical planning. S.C.
Jiang has reported a simulation of laser thermotherapy and discussed the dependence of
thermal damage on laser wavelength, power, duration, and other factors [
]. Paola Saccomandi
and others have proposed the use of Gaussian distributions and the Beer-Lambert law to
describe the distribution of laser power [
]. F. M. Di Matteo et al. has confirmed, through
several experiments, that numerical simulation results based on the same model agree well
with the experimental results , suggesting that numerical simulations can be used to guide
laser ablation treatments. However, the aforementioned studies do not take the synergistic
effect of multiple optical fibers into account, and the outcome of dual-fiber laser ablations
Therefore, to realize a more efficient and minimally invasive treatment of large prostate
cancer, the dual-fiber laser is proposed in this study. Both finite element simulations and
canine prostate ablation experiments are performed to investigate the feasibility of the
proposed treatment. A model has been established to predict the treatment outcome and to
examine the effects of the distances between the two fibers of the laser. Both ex vivo and in vivo
canine prostate dual-fiber ablation experiments are performed, and the lesion volumes are
studied using MRI, ultrasound imaging, and histological analysis.
2. Materials and methods
2.1 Theoretical model
Heat transfer in prostate tissue is described by the well-known Pennes bioheat transfer
¼ r ðkrTðx; y; z; tÞÞ þ Qblood þ Qmet þ Qlaser
where ρ, c, and k are the density (kg/m3), specific heat capacity (W/m3 K), and thermal
conductivity (W/m K) of prostate tissue, respectively, and T(x,y,z,t) is the tissue temperature.
The patterns in Eq (1) are defined below:
(1) Qblood (W/m3) is the heat transfer resulting from blood perfusion, which can be
Qblood ¼ obrbcb ðT TbÞ
The blood perfusion rate of prostate tissue (ωb) was set to be 0.0056 1/s [
] for in vivo tissue
and 0 1/s for ex vivo tissue.
(2) Qmet (W/m3) is the metabolically generated heat. The transient heat flux produced by
the laser is extremely large and the heat flux resulting from metabolism is very small.
Therefore, the metabolic heating rate was neglected, i.e., Qmet = 0.
(3) Qe (W/m3) is the energy change due to water evaporation and can be expressed by the
where hfg) (J/kg) is the latent heat for water evaporation, and ρw denotes the water content of
the tissue (kg/m3). Assuming that all water, including water vapor, remained in the prostate
tissue, ρw can be expressed as follows [
>8>>>>>< 0:778 0@1 eT 3:41206A1 T 103
>>>> 0:0289 T3
Qe can be expressed as the following form:
8:924 T2 þ 919:6T
Qlaser ¼ a Iðx; yÞ e az
Iðx; yÞ ¼ I0 e ðx2þy2Þ=2s2
where I0 = P/(2 π σ2). P(W) is the laser power in continuous output mode.
In Eq (8), the Beer-Lambert law describes the absorption of laser power along the z
direction; α denoted the absorption coefficient (m−1), which is dependent on the laser wavelength
and tissue properties. As the scattering coefficient in tissues was large, it could not be ignored;
Therefore, Eq (1) can be re-written in the following form:
¼ r ðkrTðx; y; z; tÞÞ þ Qblood þ Qlaser
(4) Qlaser (W/m3) denotes the laser power absorbed by tissues and can be described by the
Beer-Lambert law [
where I(x,y)(W/m2) denotes the power distribution of the laser along the cross section of the
incident direction and can be described by a two-dimensional Gaussian distribution. σ = rf/3
was used to ensure that 99% of the laser power was concentrated on the cross section of the
optical fiber core (rf is the radius of the optical fiber):
thus, an effective attenuation coefficient αeff was used to replace α in Eq (8):
aeff ¼ 3aða þ asð1 gÞÞ ð10Þ
where g is the anisotropic coefficient and αs is the scattering coefficient (m−1). The variables α,
αs and g are temperature dependent, due to tissue denaturation, and values corresponding to
the temperature range were applied.
Eq (1) can be solved using finite element analysis software; the parameters in the model
were summarized in Table 1. To investigate the properties of the heating process in dual-fiber
laser ablation, two optical fibers were implanted into the tissue in parallel. Four different
spacing distances between the optical fibers (3, 6, 9, and 12 mm) were evaluated. The initial
temperature of the in vivo tissue was 37˚C, the initial temperature of the ex vivo tissue was 20˚C,
and the external boundary condition was set at a constant temperature (20˚C). In this study, a
laser power of 5 W was applied with duration of 100 s, which was commonly used in clinic.
The finite element numerical method was applied, and all the simulations were carried out in
COMSOL Multiphysics 5.2.
2.2 Experimental evaluation
This research was conducted in accordance with the guidelines of the National Institutes of
Health of China for the care and use of laboratory animals. All of the experimental protocols
were approved by the local animal ethics committee affiliated with Shanghai Sixth People’s
2.2.1 Experimental methods for ex vivo canine prostate tissue. Three 5-year-old beagle
dogs were euthanized by injections of potassium cyanide (KCN) under general anesthesia, and
the prostate was excised. Each individual ex vivo canine prostate was embedded in 2% agar
and placed in a custom-made acrylic box. One side of the box served as a guide-groove
substrate, which was used to position the optical fibers. The guide-groove substrate consisted of
two plates, on each of which 10×10 mesh holes with diameters of approximately 1.5 mm were
evenly distributed, with 3-mm spacing between holes (Fig 1A and 1B).
Under ultrasound guidance, two optical fibers were inserted into one lateral lobe of the ex
vivo prostate; the two laser fibers were parallel, the insertion depth was 12 mm, and the
distance between the top and bottom fibers was 6 mm (Fig 1C). The laser therapeutic system used
was an Echolaser X4 emitting a Nd: YAG laser set at a wavelength of 1064±10 nm. The length
and diameter of the laser fiber used were 15 cm and 300 μm, respectively, and the external
diameter of the equipped guided percutaneous needle was 18 G. The power and energy output
of each laser fiber were 5 W and 500 J, respectively. The Mylab Twice ultrasound
interventional system (5 MHz ultrasound) was used to monitor the ablation process (Fig 1D). After
30 (T<65˚C); 40 (T>65˚C)
8000 (T < 65˚C); 18000 (T > 65˚C)
4 / 13
Fig 1. Ex vivo prostate placement and optical fibers arrayed in ex vivo and in vivo experiments. (a) The prostate was sealed in a box filled with agar. (b) One
side of the box consisted of a guide-groove substrate. (c) The positions of the two optical fibers were monitored using ultrasound in ex vivo experiments. (d) The
ablation process revealed a large number of strong-echo light spots. (e) The positions of the two optical fibers were monitored using ultrasound in in vivo
experiments. (f) The cross section of the short axis revealed the biopsy needles as two strong point echoes (white arrow).
ablation, the prostate tissues were dissected along the puncture needle. The section was
considered to be the largest ablation plane.
2.2.2 Experimental methods for in vivo canine prostate tissue. In this experiment, four
5-year-old beagle dogs were used. The following experimental procedures were performed on
each dog after anesthesia by an intravenous injection of sodium pentobarbital at 30 mg/kg.
Among the four dogs, three were ablated by double laser fibers (A-D), and one by single fiber
ablation (E). Under the guidance of a TRT23-type biplane transrectal ultrasound probe, the
fibers were inserted into the right lobe of the prostate after transrectal biopsy; the two fibers
were parallel with a spacing distance of 7 mm (Fig 1E and 1F). The power of a single optical
fiber was 5 W, and the energy output was 500 J. The vital signs of the canines were closely
monitored throughout the ablation process.
After treatment, an enhanced dynamic MRI scan was performed after an intravenous bolus
injection of the contrast agent, gadolinium-diethylenetriamine Penta-acetic acid (Gd–DTPA),
at 0.1 mmol/kg. A SIGNA-type superconductive MRI from GE was used; the magnetic field
strength was 1.5 T, and the coil used was a surface coil. After the MRI scan, the dogs rested for
30 min and then underwent transrectal prostatic contrast-enhanced ultrasound sonography
(CEUS). A TRT23-type biplane transrectal ultrasound probe (Esaote, Italy) was used; the
5 / 13
ultrasound used was a Mylab Twice (Esaote, Italy). An intravenous bolus injection was used to
inject the contrast agent, Sonovue (Bracco, Italy).
The dogs were put in general anesthesia during the whole process (about one and a half
hours), and no sign of stress was observed. After experiments, all dogs were immediately
euthanized by the injection of potassium cyanide (250mg) under general anesthesia, and the
prostate was excised for the observation of immediate effects by ablation. The right lobe of the
prostate was longitudinally sliced along the needle direction; one slice was obtained every 3
3.1 Simulation results
The simulated temperature distribution with either a single fiber or dual fibers is shown in Fig
2. A critical temperature (i.e. 60˚C) was used to define the boundary of tissue necrosis [
Fig 2A and 2B illustrated the temperature profiles for two serial single-fiber laser treatments in
a central cross-sectional plane parallel and perpendicular to the laser fiber. A single-fiber laser
lesion exhibited a drop-like shape, consistent with the findings of other studies [
Fig 2. Temperature distributions after heating for 100 s with single- and dual-fiber lasers. The isothermal
temperature is 60˚C; the red region in the center has a temperature greater than 350˚C. (a, c) Temperature
distributions in the maximum ablation plane parallel to the plane of the fibers; (b, d) Temperature distributions of the
maximum ablation plane perpendicular to the plane of the fibers; (a, b) The ablation region from two serial single-fiber
ablations with a spacing distance of 9 mm, as indicated by the dotted line. (c, d) The ablation region of dual-fiber laser
ablation with a spacing distance of 9 mm.
6 / 13
Fig 3. Temperature distributions after heating for 100 s with dual-fiber lasers with different spacing distances. (a,
b) Spacing distance of 3 mm; (c, d) Spacing distance of 12 mm.
serial fibers were spaced 9 mm apart. The total superposed lesion is illustrated by a dotted line
in Fig 2A and 2B. As seen in Fig 2A, the majority of the area between the two fibers was not
within the ablated range. In contrast, in a lesion created using a dual-fiber laser with the same
spacing of 9 mm, a much larger ablated region was created between the two fibers, particularly
in the area closer to the distal end of the laser, as shown in Fig 2C and 2D, and the lesion
exhibited a kettle-like shape. The tissue volume with a temperature higher than 60˚C defined as
ablated region can be directly calculated by the “Volume Integration” function of COMSOL.
The total ablated region produced by two serial single-fiber treatments was approximately
3.8 × 10−6 m3, while that of the dual-fiber treatment was approximately 4.85 × 10−6 m3 when
the fibers were spaced 6 mm apart.
When fibers were spaced 3 mm apart (Fig 3A and 3B), the shape of the ablated region
roughly resembled that of the ablation produced by a single-fiber laser, albeit with a larger size
(the area of the dual-fiber ablated region was approximately 50% larger). As the spacing
distance increasing, the ablated region between the two optical fibers gradually separated. When
the spacing was 12 mm, as shown in Fig 3C and 3D, a large area between the two fibers was
not ablated, and two independent teardrop shapes were formed.
To investigate the effect of the spacing distance between fibers on laser heating, spacing
distances from 3 to 12 mm with increments of 1 mm were studied. Fig 4 shows the dependence of
both the ablated volume (Fig 4A) and the maximum ablation range (Fig 4B) along the x, y, and
z directions (shown in Fig 3C and 3D) on the spacing distance. As shown in Fig 4A, the
7 / 13
Fig 4. Dependence of both ablation volume and ablation range on the spacing between fibers. (a) The dependence of ablation volume on the spacing between the
two fibers in dual-fiber ablation. The dotted line indicates the sum of the volumes of two single-fiber ablations, regarded as distance infinity. (b) The dependence of the
maximum ablation range on the spacing between the two fibers along the x, y, and z directions.
ablation volume increased dramatically with the spacing distance when the two fibers were
relatively close. The rate of increase gradually slowed when the two fibers were separated further,
and the ablated volume reached its maximum. The dotted line represented the sum of the
volumes of two single-fiber ablations (regarded as distance infinity), which was 4.8 × 10−6 m3. We
found that the ablated volume resulting from a dual-fiber ablation was significantly greater
than that of two single-fiber ablations. The maximum volume ablated by dual-fiber treatment
was approximately 5.3 × 10−6 m3, which was 10.4% higher than that produced by the sum of
two single-fiber treatments. When the two fibers were infinitely far from one another, the total
ablation volume of the dual-fiber ablation would be the superposition of the two single-fiber
ablations. Therefore, After the dual-fiber ablation volume reached its maximum value, it
would gradually decline and eventually approach the dotted line. From Fig 4B, it can be
observed that the incremental increases in the x direction were proportional to the spacing
distance between the two fibers, with the exception of cases when the two fibers were very close.
The x direction range for dual-fiber ablation was slightly larger than that for treatment with
two serial single fibers at a distance of 3 mm, due to the overlapping thermal fields. The
ablation ranges along the y and z directions decreased with an increase in the spacing distance.
Compared with two serial-fiber treatments, the ablated ranges in these two directions were
much larger for dual-fiber ablation. For example, at a distance of 7 mm, the ablation ranges
increased by 33% in the y and z directions when dual-fiber treatment was used. With increased
spacing distances this effect was more obvious, and dual-fiber ablation significantly increased
the ablation range in the y and z directions.
3.2 Ex vivo experiments
To study the feasibility of using dual-fiber lasers for treating prostate cancer, ex vivo dual-fiber
laser ablation experiments with a spacing of 6 mm were performed, and the results are shown
in Fig 5. Normal prostate tissue was pink; in contrast, the ablated regions were pale white due
8 / 13
Fig 5. Comparison of the ablation zone in ex vivo experiment. (a) The shape of simulated ablation zone; (b) The
shape of experimental zone.
to protein denaturation. Fig 5 compared the simulated ablative regions with the ex vivo
ablations performed on canine prostate tissue. We found that the shapes of the cross sections of
the maximum ablations in simulations and experiments resembled each other.
The transverse diameter (x) and anteroposterior diameter (z) of the ablation regions were
measured and found to be 17±0.5 mm and 19±0.7 mm, respectively. The simulated values of
the transverse diameter (x) and anteroposterior diameter (z) were found to be 16.30 mm and
18.11 mm, respectively, which were close to the experimental values. The maximum ablation
area was defined by ImageJ. The maximum ablation areas in ex vivo tissue was 238 mm2 for
the simulations and 259±67 mm2 for the experimental results.
3.3 In vivo experiments
Contrast-enhanced MRI and CEUS images of in vivo canine prostate tissues after dual-fiber
laser ablation treatment (spacing distance 7 mm) are shown in Fig 6. As the tissues in the
prostate ablation area were inactivated and lacked a blood supply, they were observed as filling
defects, while the surrounding normal prostate tissues exhibited hyper-enhancement.
Therefore, there was a distinct contrast between the two types of tissues. The three maximum radial
lines of the ablation zone were measured via MRI and CEUS imaging, shown in Table 2.
In excised pathological prostate tissues after treatment, four regions were generally
observed: a cavitation region containing the fibers, a carbonization region between the fibers, a
pale region of inactivated tissue, and the edema cincture interfacing with normal tissues (Fig
7A). As the edema cincture was clearly visible, the size of the ablation zone could also be
accurately measured. The histopathological results revealed that the tissues within the ablation
zone were all coagulation necrotic, and no area of incomplete ablation was observed between
the two fibers (Fig 7C and 7D), which was consistent with the predicted results. Fig 7B showed
a pathological section after single-fiber laser ablation, where the teardrop shape agreed well
with the simulated results shown in Fig 2A. ImageJ was used to calculate the area of the lesions
in the x-z plane in Fig 7A and 7B. The results showed that the lesion area resulting from
dualfiber treatment was 207.12 ± 10.13 mm2, while the simulated value was 239.07 mm2. The
ablated area produced by single ablations was 104.10 mm2, while the simulated value was
125.66 mm .
9 / 13
Fig 6. Comparison of the ablation zone in in vivo experiments. (a, b) Contrast-enhanced MRI imaging of the
prostate after dual-fiber laser ablation. The ablation zone shows no filling of contrast agents (red arrow) and a weak
signal; longitudinal (a) and transverse (b) sections are shown respectively; (c, d) CEUS imaging of the prostate after
dual-fiber laser ablation. The ablation zone shows no filling of contrast agents (red arrow) and a weak signal;
longitudinal (c) and transverse (d) sections are shown respectively.
The three radial lines of the ablation obtained in the simulation were similar to the actual
measurements determined by MRI, CEUS, and pathological examination (Table 2), thereby
validating the model. When using two serial-fiber treatments, the simulated y range was 0.99
cm, and z range was 1.25 cm, which were much smaller than those using dual-fiber treatments.
The maximum lesion area was 175.87 ± 8.70 mm2 in the transverse plane (x-y) and 163.81 ±
7.65 mm2 in the coronal plane (y-z) for in vivo experiments. The simulated values in the
transverse and coronal planes were 191.99 mm2 and 164.27 mm2, respectively. The experimental
area in the coronal plane was very similar to the simulated area, while the results for the
transverse plane differed by 9.1%.
x × z × y (cm)
1.82 × 1.59 × 1.35
1.70 × 1.45 × 1.12
1.95 × 1.70 × 1.40
x × z × y (cm)
1.82 × 1.72 × 1.61
1.78 × 1.65 × 1.55
2.01 × 1.75 × 1.53
x × z × y (cm)
1.71 × 1.45 × 1.20
1.65 × 1.32 × 1.13
1.80 × 1.58 × 1.25
10 / 13
Fig 7. Comparison of in vivo experimental and simulated results. (a) A pathological section after dual-fiber laser
ablation. Four regions were generally observed: a cavitation region containing the fibers, a carbonization region
between the fibers, a pale region of inactivated tissue, and the edema cincture interfacing with normal tissues (b) A
pathological section after single fiber laser ablation. (c) H&E staining of normal zone. (d) H&E staining of ablated
Both the ex vivo and in vivo experiments validated the model and confirmed that dual-fiber
treatment can widen the treatment range in all three dimensions. The extensions of the
ablation range in the y direction and z direction are significant. However, the ablation range in the
x direction does not change significantly and is determined only by the distance between the
two fibers. Dual-fiber treatment results in a larger ablation range than that resulting from two
serial single ablations, the difference in which is primarily concentrated in the central region
between the two fibers, which is shown by an increase in the y and z ranges. This difference is
because the central region is at a relatively high temperature due to the heat transfer between
the two fibers in the dual-fiber treatment, but the range at the edge region depends on heat
transfer from the unilateral fiber, which is the same as that in a single-fiber ablation. Thus,
compared to single-fiber ablation, dual-fiber ablation can achieve a greater ablation range in
less time, especially in the areas between the two fibers.
For prostate ablations, it is important to protect the urethra and surrounding vessels. The
curves in Fig 4 can be used in treatment planning, particularly for determining the spacing
distance when using dual-fiber lasers, to prevent damage to normal tissues and organs. When the
11 / 13
three dimensions of the tumor are known, an appropriate dual-fiber spacing can be chosen by
referring to these curves. When the targeted tissue exceeds the maximum ablation range,
treatment for more lesions may be considered.
This study also shows that both contrast-enhanced MRI and CEUS are effective methods to
evaluate clinical treatment outcomes and can also be used for re-examination/follow-up.
Ultrasound can guide the puncture needle into the prostate tissue for real-time positioning
and dynamic imaging. MRI has a high soft tissue resolution and can provide multi-directional,
multi-sequence imaging, resulting in a high image resolution when used in cases of prostate
disease. CEUS and contrast MRI are capable of distinguishing lesions from normal tissue by
differences in blood supply. Moreover, the two imaging modalities do not use radiation, can
be operated repeatedly, and accurately measure the lesion range. Thus, MRI and CEUS are
suitable choices for efficient evaluations and long-term follow-up.
In this study, a laser power of 5 W was chosen based on reports of clinical applications and
in vivo experiments; laser ablation powers of 4 W to 15 W have been widely used [
However, the optimal laser power for ablation remains debatable, and various laser powers
should be characterized in future studies. Based on the excellent agreement between the
simulation and the experiments, even higher energies may be used, and more fibers may be
introduced to simulate and study the ablation of larger areas. For follow-up studies of multi-focal
prostate cancer, additional fibers can be used for ablation under the guidance of simulation
results when dual-fiber ablation cannot cover all foci.
The results of the study presented herein indicate that dual-fiber laser ablation can remarkably
increase the range of the ablation zone compared with single-fiber modes. This study has
found the shape of the ablations created using dual-fiber laser ablation is quite different to the
superposition of ablations using sequential single-fiber laser ablation. In addition, the shape of
the ablation is dependent on the spacing distance between the two fibers. The results from the
theoretical model and the experimental results were found to be in good agreement. The
three-dimensional ablation ranges depended on the spacing distance between multiple fibers
were also obtained which can be used for future treatment planning. The performed in vivo
canine prostate experiments performed confirm the safety and feasibility of dual-fiber laser
treatment for enlarging the ablation range, particularly in the y and z directions.
This work was supported by the National Natural Science Foundation of China (NSFC5147
6101, 51521005), the Shanghai Science and Technology Commission of Shanghai Municipality
(16441903100), the Ministry of Science and Technology of China (2016YFC0106200), and a
clinical support department capacity building project (SHDC22015001). We would like to
thank Wei Liu for his assistance with the experiments and American Journal Experts (AJE) for
helping with language editing.
Conceptualization: Xing Wu, Aili Zhang, Bing Hu.
Funding acquisition: Bing Hu.
Investigation: Xing Wu, Kangwei Zhang, Yini Chen, Ren Wang, Lei Chen.
Methodology: Xing Wu, Kangwei Zhang, Yini Chen, Ren Wang, Lei Chen.
12 / 13
Validation: Aili Zhang, Bing Hu.
Writing – original draft: Xing Wu, Kangwei Zhang.
Writing – review & editing: Aili Zhang, Bing Hu.
1. Lee T , Mendhiratta N , Sperling D , Lepor H . Focal laser ablation for localized prostate cancer: principles, clinical trials, and our initial experience . Rev Urol . 2014 ; 16 ( 2 ): 55 . PMID: 25009445
2. Rosa MRD , Trachtenberg J , Chopra R , Haider MA . Early experience in MRI-guided therapies of prostate cancer: HIFU, laser and photodynamic treatment . Cancer Image . 2011 ; 11 ( 1A ): S3 - S8 .
Wenger H , Yousuf A , Oto A , Eggener S. Laser ablation as focal therapy for prostate cancer . Curr Opin Urol . 2014 ; 24 ( 3 ): 236 . https://doi.org/10.1097/MOU.0000000000000044 PMID: 24625427
4. Ahmed HU , Hindley RG , Dickinson L , Freeman A , Kirkham AP , Sahu M , et al. Focal therapy for localized unifocal and multifocal prostate cancer: a prospective development study . Lancet Oncol . 2012 ; 13 ( 6 ): 622 - 632 . https://doi.org/10.1016/S1470-2045( 12 ) 70121 - 3 PMID: 22512844
5. Marqa MF , Colin P , Nevoux P , Mordon SR , Betrouni N. Focal laser ablation of prostate cancer: numerical simulation of temperature and damage distribution . Biomed Eng Online . 2011 ; 10 ( 1 ): 45 - 45 .
6. Raz O , Haider MA , Davidson SR , Lindner U , Hlasny E , Weersink R , et al. Real-time magnetic resonance imaging-guided focal laser therapy in patients with low-risk prostate cancer . Eur Urol . 2010 ; 58 ( 1 ): 173 - 177 . https://doi.org/10.1016/j.eururo. 2010 . 03 .006 PMID: 20334965
7. Oto A , Sethi I , Karczmar G , McNichols R , Ivancevic MK , Stadler WM , et al. MR imaging-guided focal laser ablation for prostate cancer: phase I trial . Radiology . 2013 ; 267 ( 3 ): 932 - 940 . https://doi.org/10. 1148/radiol.13121652 PMID: 23440319
8. Jiang SC , Zhang XX . Dynamic modeling of photothermal interactions for laser-induced interstitial thermotherapy: parameter sensitivity analysis . Lasers Med Sci . 2005 ; 20 ( 3-4 ): 122 - 131 . https://doi.org/10. 1007/s10103-005 -0359-5 PMID: 16328097
9. Saccomandi P , Schena E , Di MF , Pandolfi M , Martino M , Rea R , et al. Laser interstitial thermotherapy for pancreatic tumor ablation: theoretical model and experimental validation . International Conference of the IEEE Engineering in Medicine & Biology Society, Embc Conf Proc IEEE Eng Med Biol Soc . 2011 : 5585 .
10. Saccomandi P , Schena E , Caponero MA , Di MF , Martino M , Pandolfi M , et al. Theoretical analysis and experimental evaluation of laser-induced interstitial thermotherapy in ex vivo porcine pancreas . IEEE Trans Biomed Eng . 2012 ; 59 ( 10 ): 2958 - 2964 . https://doi.org/10.1109/TBME. 2012 .2210895 PMID: 22929361
11. Di MF , Martino M , Rea R , Pandolfi M , Panzera F , Stigliano E , et al. US-guided application of Nd:YAG laser in porcine pancreatic tissue: an ex vivo study and numerical simulation . Gastrointest Endosc . 2013 ; 78 ( 5 ): 750 - 755 . https://doi.org/10.1016/j.gie. 2013 . 04 .178 PMID: 23680175
12. Duck FA . Physical Properties Of Tissue: A Comprehensive Reference Book . Ta Egerev Sv Lyamshev Lm Naugol'nykh Ka Pashin Ae . 1990 .
13. Mohammed Y , Verhey JF . A finite element method model to simulate laser interstitial thermo therapy in anatomical inhomogeneous regions . Biomed Eng Online . 2005 ; 4 ( 1 ): 1 - 16 .
14. Berjano EJ . Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future . BioMed Eng Online . 2006 ; 5 ( 1 ): 24 - 24 .
15. Diller KR , Pearce JA . Issues in modeling thermal alterations in tissues . Ann NY Acad Sci . 2010 ; 888 ( 1 ): 153 - 164 .
16. Mckenzie AL . Physics of thermal processes in laser-tissue interaction . Physics Med Biol . 1990 ; 35 ( 9 ): 1175 .
17. Panescu D , Whayne JG , Fleischman SD , Mirotznik MS , Swanson DK , Webster JG . Three-dimensional finite element analysis of current density and temperature distributions during radio-frequency ablation . IEEE Trans Biomed Eng . 1995 ; 42 ( 9 ): 879 - 890 . https://doi.org/10.1109/10.412649 PMID: 7558062
18. Mohammed Y , Verhey JF . A finite element method model to simulate laser interstitial thermo therapy in anatomical inhomogeneous regions . Biomed Eng Online . 2005 ; 4 ( 1 ): 1 - 16 .