Thesis Haichun Configuratia Electronica Tabel
Thesis Haichun Configuratia Electronica Tabel
Thesis Haichun Configuratia Electronica Tabel
Doctoral Thesis
2014
To my father
Abstract
During the past decade, upconverting nanoparticles (UCNPs)
doped with rare earth ions have become an important class of
fluorescence contrast agents for molecular imaging, due to their
unique properties. Their property of anti-Stokes luminescence,
with both the excitation and emission wavelengths close to the
optimal for biomedical imaging, has been extensively explored in
various biomedical applications. The work described in this thesis
is mainly concerned with the investigation of other unique properties of UCNPs, including nonlinearity and saturation, to improve
fluorescence diffuse imaging and tomography. The aim of this work
was also to develop a suitable method of characterizing the powerdensity-dependent quantum yield of UCNPs, and to optimize the
excitation scheme in order to facilitate their use in deep tissues.
Upconversion emission has a nonlinear dependence on the excitation intensity. This nonlinearity was exploited to improve the
image reconstruction quality in fluorescence diffuse optical tomography, by increasing the orthogonal information density through simultaneous multibeam excitation. It is also demonstrated that the
use of nonlinear UCNPs as contrast agents in fluorescence tomography can breach the current limit on spatial resolution encountered
when using linear fluorophores. UCNPs can emit multiple emission
bands with large differences in tissue attenuation under excitation
of near-infrared light. Such multispectral information was used to
create a regularization map to guide fluorescence diffuse optical
tomography, which yielded significantly better axial resolution in
the reconstructed images than using standard Tikhonov regularization.
The quantum yield of upconversion emission increases with excitation intensity and gradually reaches a plateau. In this work, an
initial model was developed to describe the quantum yield of twophoton upconversion emission as a function of the excitation intensity, based on rate equation analysis. A balancing power density
was identified, which characterizes the excitation-intensity dependence of the quantum yield. At such a power density, the quantum
yield reaches half the maximum attainable value, occurring at very
v
Abstract
high excitation intensities when UCNPs are saturated. Pulsed excitation is proposed as a superior mode of excitation to continuous
wave excitation, due to the potential of achieving higher intrinsic
quantum yields from UNCPs without increasing the average excitation power. This will fundamentally increase the applicability of
UCNPs in deep tissues.
vi
Abstract
viii
List of Publications
This thesis is based on the following papers, which will be referred
to in the text by their Roman numerals.
I Drug quantification in turbid media by fluorescence
imaging combined with light-absorption correction
using white Monte Carlo simulation
H. Xie, H. Liu, P. Svenmarker, J. Axelsson, C. T. Xu,
S. Gr
afe, J. H. Lundeman, H. P. H. Cheng, S. Svanberg,
N. Bendsoe, P. E. Andersen, K. Svanberg,
S. Andersson-Engels.
Journal of Biomedical Optics 16(6), 066002-1 - 066002-11
(2011).
II Synthesis of NaYF4 :Yb3+ , Er3+ upconverting
nanocrystals in a capillary-based continuous
microfluidic reaction system
H. Liu, O. Jakobsson, C. T. Xu, H. Xie, T. Laurell,
S. Andersson-Engels.
Proc. of SPIE 7909, 790917-1 - 790917-6 (2011).
III Upconverting nanoparticles for pre-clinical diffuse
optical imaging, microscopy and sensing: Current
trends and future challenges
C. T. Xu, Q. Zhan, H. Liu, G. Somesfalean, J. Qian, S. He,
S. Andersson-Engels.
Laser & Photonics Reviews 7(5), 663-697 (2013).
IV Balancing power density based quantum yield
characterization of upconverting nanoparticles for
arbitrary excitation intensities
H. Liu, C. T. Xu, D. Lindgren, H. Xie, D. Thomas,
C. Gundlach, S. Andersson-Engels.
Nanoscale 5, 4770-4775 (2013).
ix
List of Publications
List of Publications
xi
Abbreviations
CT
computed tomography
CW
continuous wave
ED
electric dipole
ESA
ETU
FDOT
MC
Monte Carlo
MPE
MRI
NIR
near-infrared
QY
quantum yield
RE
rare earth
RTE
TTA
triplet-triplet annihilation
UCNP
upconverting nanoparticle
xiii
Contents
1
Introduction
1.1
Nonlinearity and quantum yield characterization of upconverting nanoparticles . . . . . . . . . . . . . . . . . . . . . . .
1.2
Excitation scheme optimization of upconverting nanoparticles
1.3
Aims and outline of this thesis . . . . . . . . . . . . . . . . . .
7
7
7
10
12
13
15
18
19
Upconverting nanoparticles
3.1
Rare-earth luminescence . . . . . . . . . . . . . . . . . . .
3.1.1
Spectroscopic properties of rare earths in solids .
3.1.2
Upconversion mechanisms . . . . . . . . . . . . .
3.1.3
Efficient upconversion systems . . . . . . . . . . .
3.2
Material engineering of upconverting nanoparticles . . . .
3.2.1
Size and morphology control . . . . . . . . . . . .
3.2.2
Luminescence enhancement . . . . . . . . . . . .
3.2.3
Excitation and emission wavelength optimization
3.2.4
Biocompatibility . . . . . . . . . . . . . . . . . . .
3.2.5
Multi-functionalization . . . . . . . . . . . . . . .
3.3
Biomedical applications . . . . . . . . . . . . . . . . . . .
3.3.1
Toxicity assessment . . . . . . . . . . . . . . . . .
3.3.2
Promising applications . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
23
23
23
26
29
31
31
32
34
35
35
36
36
38
Quantum yield characterization and excitation scheme optimization for upconverting nanoparticles
4.1
Quantum yield characterization . . . . . . . . . . . . . . . . .
4.1.1
Definition of the quantum yield . . . . . . . . . . . .
4.1.2
Power-density dependence of the quantum yield . . .
4.1.3
Two-parameter characterization . . . . . . . . . . . .
4.1.4
Experimental measurements . . . . . . . . . . . . . .
4.1.5
Extension to triplet-triplet annihilation upconversion
4.2
Excitation scheme optimization . . . . . . . . . . . . . . . . .
4.2.1
Pulsed excitation for higher quantum yield . . . . . .
4.2.2
Low light limit . . . . . . . . . . . . . . . . . . . . . .
41
41
42
43
43
45
46
48
48
49
.
.
.
.
.
.
.
.
.
.
.
.
.
2
3
4
20
Contents
4.2.3
4.2.4
5
ANSI standard . . . . . . . . . . . . . . . . . . . . .
Single-shot imaging . . . . . . . . . . . . . . . . . . .
50
51
53
53
54
55
57
58
61
61
62
62
63
64
67
Acknowledgements
73
References
77
Contents
Papers
I
105
Synthesis of NaYF4 :Yb3+ , Er3+ upconverting nanocrystals in a capillary-based continuous microfluidic reaction
system
119
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current trends and
future challenges
127
Balancing power density based quantum yield characterization of upconverting nanoparticles for arbitrary excitation intensities
165
Deep tissue optical imaging of upconverting nanoparticles enabled by exploiting higher intrinsic quantum yield
through use of millisecond single pulse excitation with
high peak power
173
183
189
199
II
III
IV
VI
VII
IX
205
Chapter
Introduction
Molecular imaging is a rapidly developing biomedical research discipline, the aim of which is the early detection of diseases via visual
representation, characterization, and quantification of physiology
and molecular signaling at the cellular and subcellular levels within
intact living organisms [1, 2]. It has now become an indispensable
tool, particularly in cancer research, clinical trials and medical
practice. Among these techniques, fluorescence imaging has become one of the most commonly used modalities, due to its ability
to visualize, for example, gene delivery [3], angiogenesis [4], and
apoptosis [5, 6], together with its simplicity and cost-effectiveness.
In fluorescence molecular imaging, contrast agents are of essential importance, as they are the final reporters of the biological
processes being probed. These can be either natural chromophores
in the body, or externally introduced fluorophores with the desired
spectroscopic properties. The latter are particularly widely employed, as they provide a large degree of freedom in tuning the fluorescence pathways in the way people expect. The most commonly
investigated and applied exogenous contrast agents include fluorescent dyes (e.g. indocyanine green and Cy5.5) [6, 7], fluorescent proteins [4, 8, 9], and semiconductor quantum dots [1017]. Despite
considerable success, the use of these fluorophores is limited in
fluorescence molecular imaging techniques by problems associated
with tissue autofluorescence, poor spatial resolution and limited
light penetration, making it difficult to image deeply located regions with high sensitivity [18]. Furthermore, photobleaching of
traditionally used organic fluorophores limits the possibility of the
repeated imaging required in longitudinal studies.
An emerging group of contrast agents - upconverting nanoparticles (UCNPs) have the potential to overcome the limitations discussed above due to their unique properties. UCNPs are nanosized inorganic crystals, doped with certain trivalent rare earth
(RE) ions that act as sensitizers and activators [19]. The most no1
1.1
Introduction
both fluorescence imaging and tomography, providing higher spatial resolution than conventional linear fluorescent contrast agents
such as dyes and quantum dots. In addition, the nonlinear fluorophores being probed in fluorescence diffuse optical tomography
(FDOT) using UCNPs as contrast agents could make different excitation beams interfere with each other, thus providing additional information useful for image reconstruction. Such additional
information is contained in the difference between the fluorescence
images collected when using the individual excitation beams and
that obtained with simultaneous, multiple excitation beams.
In contrast to other multi-photon processes, which are instantaneous, such as second-harmonic generation and two-photon fluorescence, the upconversion process is relatively slow, as it involves
real intermediate energy states. These intermediate levels may become saturated at high excitation intensities, leading to a less steep
excitation-power-density dependence of the fluorescence intensity.
Due to the nonlinearity of UCNPs, exhibiting saturation, their
quantum yield, which is defined as the ratio between the number of excitation photons and the number of emission photons,
increases with excitation intensity, finally approaching a constant
value [25]. Characterization of such a dependence on the excitation intensity of upconversion emission is important in guiding the
practical use of UCNPs in biomedical applications. In addition, it
could potentially provide standard indices for characterizing UCNPs with respect to energy conversion, whereas at present, local
references are employed in different labs.
1.2
1.3
Introduction
Chapter
2.1
2.1.1
Light-tissue interactions
Tissue scattering
(2.1)
(2.2)
Here, the coefficient a is related to the number density of the scatterers, and b to the size of the scatterer.
The scattering phase function, p(s0 , s), describes the normalized angular distribution of the power of the scattered field, fulfilling the expression:
Z
p(s0 , s)d 0 = 1.
(2.3)
4
The phase function is usually assumed to depend only on the deflection angle , and not on the azimuthal angle , yielding:
p(s0 , s) p(s0 s) = p(cos ).
8
(2.4)
10
10
p()
650 nm
800 nm
975 nm
10
10
180
90
Multiple scattering
For a collection of scatterers, each scatterer will contribute to the
final scattered light field. If the number density of scatterers, Nsca
[cm3 ], is small, the contribution of each scatterer can be regarded
as linearly independent, and thus the scattering property of the
medium can be characterized by a scattering coefficient, s [cm1 ]:
s = Nsca sca ,
0
[]
90
180
(2.5)
2
10
p(cos ) =
1
1 g2
.
4 (1 + g 2 2g cos )3/2
g=0.9
g=0.75
g=0
10
p()
10
10
10
10
180
90
0
[]
90
180
(2.6)
where a0 is a scaling factor, fRay is the fraction of Rayleigh scattering, (1 fRay ) is the fraction of Mie scattering, and bMie is the
scattering power of Mie scattering.
2.1.2
Tissue absorption
(2.9)
10
10
sum
H2O
10
Hb
HbO2
10
a (cm1)
10
lipid
10
10
10
10
300
400
500
600
700
800
Wavelengthp(nm)
900
1000
Figure 2.3. Absorption coefficients for the main absorbing chromophores in biological tissues and the corresponding synthetic absorption coefficient. The absorption coefficients were calculated using 70 %
water [53], 3 % blood volume fraction with 60 % oxygen saturation [54],
and 15 % lipids [55].
(2.12)
11
2.1.3 Luminescence
S2
IC
ISC
S1
Fluorescence
Absorption
T1
en
ph
os
Ph
sc
re
ce
S0
Ground state
2.1.3
Luminescence
After the energy of the light has been absorbed by the chromophores, one portion of the absorbed energy is generally dissipated nonradiatively within the chromophores through the
medium, finally resulting in heat. The other portion of the absorbed energy is emitted as light or luminescence. This luminescence may be named differently depending on the mechanisms involved, sometimes even for the same chromophore. For instance,
the transition between singlet states in molecules yields fluorescence, which has a lifetime typically on the order of nanoseconds.
The transition between an excited triplet state and a lower-lying
singlet state generates phosphorescence, with much longer lifetimes
than fluorescence due to the violation of the transition selection
rule. A Jablonski diagram illustrating the processes of light absorption and luminescence is shown in Figure 2.4. A chromophore
that can generate luminescence upon excitation is usually called
as a fluorophore.
Compared with light absorption and scattering, luminescence can provide structural and functional information on tissues with much higher spatial resolution, sometimes even at
the single molecule level. Hence, the luminescence from endogenous and exogenous fluorophores is used extensively in various biomedical applications based on optical techniques. Endogenous fluorophores that exist naturally in biological tissues
include reduced nicotinamide-adenine dinucleotide and reduced
nicotinamide-adenine dinucleotide phosphate [59, 60], flavin adenine dinucleotide [61], collagen and elastin [62, 63], lipopigments
[61, 64], amino acids and porphyrins. Exogenous fluorophores
12
are introduced into living organisms using biotechniques, and include, for example, fluorescent dyes [7, 65], fluorescent proteins
[4, 9, 66, 67] and semiconductor quantum dots [68, 69]. In the
present work, RE-ion-doped inorganic UCNPs were mainly utilized
as exogenous luminescent biomarkers in various biomedical applications, including fluorescence diffuse optical imaging and FDOT.
The UCNPs will be introduced in detail in Chapter 3, and their
optical properties and applications are presented in Chapters 4
and 5.
2.2
S
iii)
V
i)
ii)
iv)
i)
v)
i)
cs
ii)
iv)
iii)
qdV ,
(2.13)
| {z }
v)
(2.14)
Another quantity, L(r, s, t) [Wm2 sr1 ], called the radiance,
is often of interest, and describes the radiant power flow along s
through unit area at position r and time t. It is related to the
photon distribution function by [78]:
L = N hvc,
(2.15)
(2.16)
The RTE is analytically solvable only in very simple cases such
as one-dimensional geometries [79]. In most cases, the RTE is
solved using numerical methods, or is simplified based on certain
approximations so that simple analytical solutions may be found.
Such techniques include the spherical harmonics method [80, 81]
and the diffusion approximation [80], the simplified spherical harmonics method [8286], the discrete ordinates method [87], and
Monte Carlo (MC) simulations [88, 89]. In this thesis, the diffusion approximation and MC simulations were mainly used, and
they will be introduced in the following sections.
1 Here,
that (N
s) = N
s + N
s = N
s, as
s is fixed.
14
2.3
Diffusion theory
3
1
(r, t) +
J(r, t) s,
4
4
(2.17)
where (r, t) [Wm2 ] is the fluence rate, and J(r, t) [Wm2 ] is the
flux, fulfilling the expression:
Z
(r, t) = L(r, s, t)d,
(2.18)
4
and
J(r, t) =
L(r, s, t)sd,
(2.19)
1
= D,
3(a + 0s )
(2.23)
(2.24)
(2.25)
2 The
15
(2.26)
(2.27)
16
2/(1 R0 ) 1 + | cos c |3
,
1 | cos c |2
(2.28)
where
R0 =
(n n0 )2
,
(n + n0 )2
(2.29)
and
c = arcsin
n
0
(2.30)
z0-zb
0
z+0
n0
r0-
r0+
r
n
Analytical solutions
For some simple geometries, there are analytical solutions to the
diffusion equation. For a homogeneous infinite medium, with a
point source expressed as q0 (r, t) = P0 (r)(t), where P0 [W] is
the power of the source, the solution is [104]:
(r, t) = cP0 (
1
r2
)3/2 exp(
a ct).
4Dct
4Dct
(2.31)
P0 exp(eff r)
,
4D
r
z
Figure 2.6. A homogeneous
medium with a semi-infinite
geometry.
(2.32)
,
(2.34)
(r) =
4D
r0+
r0
q
where r0+, = r2 + (z z0+, )2 .
In an infinite slab geometry as shown in Figure 2.7, the extrapolated boundary condition must be fulfilled at both interfaces, and
the z-positions of the resulting positive and negative sources can
17
z-1z-1+
zm
(2.35)
z0+ ,
m = 0, 1, 2, 3, . . . ,
(2.36)
(2.37)
n0
z0z
+
0
r
n
z1z+1
z
(r) =
10
0
10
1
10
2
10
3
10
4
10
5
10
5
10 15 20 25 30
Depth (mm)
18
(2.38)
q
+, 2
+,
) .
= r2 + (z zm
where rm
The inclusion of mirrored sources of the first few orders is usually sufficient to provide good accuracy. Figure 2.8 presents the
fluence rate as a function of depth in a homogeneous slab, calculated using Equation (2.38) with m = 0 and 1.
2.4
Figure 2.7. A homogeneous
medium with a slab geometry.
m=+
+
P0 X exp(eff rm
) exp(eff rm
)
,
+
4D m=
rm
rm
2.4.1
Initialize'particle
(W,r,s)
(W,r',s)
a
W,
t
(2.39)
where W is the particles current weight, and t is the total interaction coefficient, equal a + s . This fraction of weight will
be deposited in the local grid element. The particle with the new
weight, i.e., W W , will undergo scattering.
Photon absorption is also involved in the probability density
function in the sampling of the step size, s, for particle motion
[111], i.e.:
p(s) = t exp(t s),
(2.40)
originating from the Beer-Lambert law for light propagation.
A pseudo-random number, , uniformly distributed over the
interval [0, 1] is generated. A non-uniformly distributed variable,
, can generally be obtained by one-to-one mapping between
and . The mapping function is defined by the equality of the
cumulative probabilities for these two variables:
F (1 ) = F (1 ).
(2.41)
For the step size obeying the probability density function stated
in Equation (2.40), it thus can be mapped by:
ln(1 )
,
s=
t
(2.42)
Absorption'E'Scattering
(W',r'',s')
Transmit?'
N
Weight'small?'
(W',r''',s'')
Y
Survive'roulette?'
Last'particle?'
Y
End
ln()
,
t
(2.43)
(2.44)
2 1,
if g=0
The simulation of reflection and transmission
When the particle hits the boundary of the medium, internal reflection and transmission occur. The reflection probability is determined by the Fresnel reflection coefficient R():
1 sin2 (i t ) tan2 (i t )
R(i ) =
+
,
(2.46)
2 sin2 (i + t ) tan2 (i + t )
where i is the incident angle, and t is the transmission angle,
given by Snells law:
ni sin i = nt sin t .
(2.47)
2.4.2
21
Chapter
Upconverting nanoparticles
Exogenous luminescent biomarkers are often introduced to enhance fluorescence imaging and tomography. UCNPs, composed
of RE ions doped into an inorganic host material, are an emerging group of biomarkers with excellent spectroscopic and physicochemical properties. In this chapter, the spectroscopic properties
of RE ions are first discussed. Then, the upconversion mechanisms,
based on the interactions between dopant ions, are introduced.
Finally, the material engineering and biomedical applications of
UCNPs are reviewed.
Rare-earth luminescence
Spectroscopic properties of rare earths in
solids
4f0
4f1
4f 5d 6s
4f 5d 6s
Pr
4f2
Nd
4f3
Pm
4f4
Eu
Gd
3.1.1
Ce
Sm
3.1
RE3+
RE
0
4f5
4f6
4f7
4f 5d 6s
4f 5d 6s
4f 5d 6s
Tb
4f 5d 6s
Dy
4f8
4f9
Ho
4f10
Er
4f11
Tm
4f12
Yb
Lu
14
4f 5d 6s
4f14 5d1 6s2
4f13
4f14
23
4f
104 cm-1
H0
HC
103 cm-1
HSO
102 cm-1
HCF
(3.1)
24
Upconverting nanoparticles
25
moment operator (Deff ), expanded in tensor operator space, is employed in order to connect 4f states, i.e. [122124]:
ESA
X | >< | V
(0)
M
/
GSA
E0 E
+V
X | >< |
(0)
M
/
E0 E
D(1) + . . .
(3.2)
Here, D(1) is the ED operator, | > are excited states of the
RE ions or ligand, and M is the model space. The potential V
is the sum of the crystal field, Coulomb interaction and spin-orbit
interaction. At first order in V , only the odd parity part of the
crystal-field potential can connect the 4f configurations of RE3+
ions. The parameters in the expansion are determined by fitting
with experimental data.
Early work by Judd and Ofelt yielded the Judd-Ofelt theory
[125, 126], which has been extensively used in areas such as laser
and phosphor design for calculating the total intensities of transitions between J -multiplets [127]. The key accomplishment of
the Judd-Ofelt theory is that the ED transition strength can be
expanded as the sum of three even-ranked tensors, with certain
assumptions, and the standard parameters are usually denoted
( = 2, 4, 6). However, the Judd-Ofelt theory should be applied
with care, as it is based on a poor approximation even for room
temperature conditions, i.e., that there is equal occupations of the
states of the ground multiplets [128].
The symmetry and strength of the crystal field play an important role in the splitting of 2S+1 LJ levels and radiative decay
rates between 4f-levels. The local crystal field around RE ions can
be tailored by ion doping, leading to altered transition intensities
[129, 130].
3.1.2
Upconversion mechanisms
When RE3+ ions are incorporated into solids, upconversion emission with anti-Stokes shift can occur in the absence of direct multiphoton processes, through the accumulation of excitation energy.
Upconversion is complicated in practice, but several basic mechanisms contribute to upconversion emission, including direct photoexcitation and energy transfer between ions. These will be discussed in the following.
Excited state absorption
The first identified upconversion mechanism is excited state absorption (ESA), following ground state absorption, as illustrated
in Figure 3.3. In 1959, Bloembergen proposed that infrared photons could be detected and counted through sequential absorption
within the levels of a given ion in a solid [131]. This can be re-
26
Upconverting nanoparticles
E2
Energy transfer upconversion (ETU) is a far more efficient upconversion process than ESA, and was first proposed by Auzel
[132, 133] and by Ovsyankin et al. [134] independently in 1966.
The process is illustrated in Figure 3.4. In a microscopic view, the
ETU process is achieved by nonradiative energy transfer between
two neighboring ions. Such energy transfer can take place not only
between excited- and ground-state ions (e.g., ETU1 in Figure 3.4),
but also between excited ions (e.g., ETU2 in Figure 3.4) [132, 133].
In addition, the participating ions are not limited to ions of the
same type [132, 134]. The energy transfer process between the
same type of ions that results in the population of an intermediate level is known as cross relaxation. In the field of upconversion
emission, the energy donor and accepter involved in the ETU process are usually termed the sensitizer and activator, respectively.
The mutual interactions between the two ions are Coulomb interactions of the van der Waals type. F
orster first treated such
interactions theoretically using quantum mechanical theory, considering the dipole-dipole interaction [135]. Later, this theory was
extended by Dexter to include higher multipole and exchange interactions [136]. Generally, the probability of energy transfer for
electric multipolar interactions can be expressed as [136, 137]:
s
1 R0
pSA =
,
(3.3)
s R
E1
GSA
ETU1
where s [s] is the actual lifetime of the excited state of the sensitizer including the contribution of multiphonon radiative decay, R
[m] is the distance between the two ions, while R0 [m] is the critical
transfer distance for which energy transfer and spontaneous deactivation of the sensitizer have equal probability, and s is a positive integer associated with the index of the multipolar interaction, with
s = 6 for dipole-dipole interactions, s = 8 for dipole-quadrupole
interactions, and s = 10 for quadrupole-quadrupole interactions.
However, this form of power law is less applicable when calculating the energy transfer probability, because the critical distance
cannot be easily obtained theoretically. Another calculation technique was developed by Kushida [138] for multipolar interactions,
applying Racahs tensorial methods. Pouradier and Auzel treated
magnetostatic and exchange interactions in a similar fashion [139].
In many cases, the ETU process is non-resonant, as illustrated
in Figure 3.4, where phonons play a critical role to account for the
energy mismatches. Actually, energy transfer with energy mismatch as high as several thousand inverse centimeters can take
27
place, with the assistance of multi-phonons [140, 141]. Unawareness of this once led to the misinterpretation of upconversion emission from a Yb3+ /Tm3+ pair as being due to cooperative sensitization [142]. The phonon-assisted energy transfer probability
depends exponentially on the energy mismatch E [143]:
pphonon (E) = pphonon (0) exp(E),
(3.4)
where pphonon (0) is the energy transfer probability when the zero
phonon lines overlap, and is a constant determined by the characteristics of the host lattice and the electron-phonon interaction.
A tremendous number of ions are assembled in a macroscopic
system, even in a very local region. In such cases, the energy transfer between ions is very complicated. Basically, three types of energy transfer can occur: sensitizer-sensitizer, activator-activator,
and sensitizer-activator transfer. The first two represent excitation migration or energy diffusion through the crystal. In order to
obtain an association with experimental data acquired from real
macroscopic samples, a statistical analysis of the energy transfer
between many ions is needed, taking into account all three types
of energy transfer. Average transfer probabilities can be calculated using certain assumptions regarding the concentrations of
the doping ions and their spatial distributions [136, 137, 144, 145].
Excitation migration could be a limiting step, or not, depending on
the sensitizer concentration, leading to different decay behavior of
the luminescence [145]. The case of high sensitizer concentration,
when the migration rate is faster than spontaneous sensitizer decay
or sensitizer-activator energy transfer, is called the fast migration
case. In such a case, the energy transfer probability is proportional to the activator concentration, but not depending on the
sensitizer concentration [145]. Otherwise, the probability is also
dependent on the sensitizer concentration. Rate equations dealing
with the populations of ions in given states can also be used for the
macroscopic system, especially in the fast migration case [19]. The
validity of this has been discussed from first principles by Grant
[146].
Excitation migration is of essential importance in upconversion
emission, as it allows energy transfer over a long distance in a
stepwise manner. This makes the pumping of activators possible,
even if they are located far away from the sensitizer ions excited
directly by photons. Such a feature opens up the possibility of
designing new efficient upconverting nanomaterials with core-shell
structure, in which the sensitizers and activators are isolated in
different layers to suppress concentration quenching1 [147149].
1 High doping levels of dopant ions lead to deleterious energy transfer
among activator ions or back energy transfer from activator ions to sensitizer
ions, resulting in quenching of upconversion emission.
28
Upconverting nanoparticles
Cooperative upconversion
Another upconversion mechanism consists of cooperative processes, including cooperative sensitization (Figure 3.5) [150, 151]
and cooperative luminescence (Figure 3.6) [152], which are about
3-5 orders of magnitude less effective than ETU [19]. The cooperative processes involve cooperative pair states, and are thus
likely to take place within clusters, and must only be considered
practically when ETU cannot take place [143]. Examples of such
cases are when the real single-ion levels allowing energy transfer do
not exist, and when the concentration is too low to allow efficient
transfer. Apart from the detection of RE ion clusters [153], very
few applications of cooperative processes exist.
The synergistic effect of all the upconversion mechanisms discussed above and the basic ground state absorption process can
lead to complex upconversion phenomena, such as the photon
avalanche effect [154158] and the looping upconversion processes
[159, 160].
3.1.3
Host materials
The host material is required to have good lattice matching with
the dopant ions. Inorganic compounds containing RE3+ , alkaline
earth and some transition metal ions (e.g., Zr4+ , Ti4+ and Mn2+ )
are thus suitable host candidates. The cations of the host materials
are generally chosen so as to be optically inert in the spectral range
of interest, in order to suppress luminescence due to the host ions.
Sc3+ , Y3+ , Gd3+ , La3+ and Lu3+ are commonly used as host ions,
but not other RE3+ ions.
As discussed in Subsection 3.1.1, the crystal field of the host
lattice has an important influence on the optical properties of a
given ion, as it determines the spectral positions and splitting of
certain optical transitions. In addition, the uneven components
of the crystal field, present when a RE ion occupies a crystallographic site without inversion symmetry, are able to relax the
parity selection rule by mixing a small amount of opposite-parity
wave functions into the 4f wave functions [162]. Thus, a host
lattice with low symmetry is usually preferable [163]. Although
29
Intensityo(a.u.)
YbEr
YbTm
YbHo
0.5
0
400
600
800 1000
Wavelengtho(nm)
30
Upconverting nanoparticles
3.2
3.2.1
Size and morphology are two major concerns when using nanoparticles as exogenous contrast agents in living organisms [187, 188].
Different sizes are required, depending on the applications. In
order to trace biological activity on the subcellular level, particles
are generally required in the range of 4-10 nm, comparable to most
31
RE3+
NaF
Mixer
180
Oil baths
Microcapillary
110
Nanocrystals
Figure 3.9. Schematic description
of the microfluidic synthesis of
upconverting nanoparticles.
3.2.2
Luminescence enhancement
(3.5)
32
Upconverting nanoparticles
shell
core
dopant
shell
4
4
3.2.3
core
F5/2
F3/2
I15/2
I13/2
I11/2
4
I9/2
Nd3+
Nd3+
Yb3+
A3+
34
version process, such as the host lattice, doping level and energy
dissipation on the nanoparticle surface, as well as the excitation
intensity. Several approaches for upconversion emission enhancement have been developed by considering these factors. Typical
methods include the use of ion doping, e.g., lanthanide and transition ions, to promote the formation of the hexagonal phase of
NaREF4 nanoparticles [182, 220, 221], and alkali ions to adjust
the local environment of dopant ions [129, 130, 222, 223]. Another
important concept for the enhancement of upconversion emission
is the adoption of a core-shell structure (Figure 3.10). A shielding
layer on the surface of the UCNPs can reduce crystal defects and
protect the optically active ions from coupling with the vibrational
modes of, for example, -OH or -NH2 groups in the solvent, thereby
reducing the nonradiative energy losses [224, 225]. The shell can
be either inert or active [172, 209, 210, 226, 227]. In addition, the
concentration quenching threshold of upconversion emission can be
also exceeded via spatial separation of dopant ions [147, 228230].
When designing new UCNPs for biomedical applications, the optical properties of tissues discussed in Chapter 2 are of major importance. Ideal UCNPs for such applications should have excitation
and emission bands located in the tissue optical window, in order to achieve maximum penetration depth. Strictly speaking, the
975 nm excited UCNPs are less optimal, because the excitation
wavelength overlaps an absorption peak of water (as shown in Figure 2.3), which constitutes a major part of most tissues. Zhan
et al. suggested the use of 915 nm excitation [231]. The problem of overheating can be partly alleviated, but this wavelength is
associated with decreased emission efficiency.
The energy transfer properties between RE3+ ions provide opportunities for tuning the luminescence pathway. A novel upconversion emission pathway was recently achieved by cascade sensitization through introducing Nd3+ ions as co-sensitizers [232],
where 800 nm light was used as excitation. Bright upconversion
emission with efficiency comparable to that under 975 nm excitation can be obtained by using a core-shell structure with Nd3+
ions in the shell and the emitting ions in the core [148, 149], as
illustrated in Figure 3.11. Compared with the most promising upconversion system for biological applications discovered previously,
i.e., Yb3+ /Tm3+ -codoped fluoride nanoparticles with excitation
and emission wavelengths of 975 nm and 800 nm, respectively
[21], the major emission bands (green and red) of such new types
of nanoparticles are less optimal. However, the optimization of
the excitation wavelength will still provide significant signal gain,
Upconverting nanoparticles
3.2.4
Biocompatibility
3.2.5
Multi-functionalization
beneficial to combine the optical imaging modality with other bioimaging modalities that have higher spatial resolution but lower
sensitivity, such as X-ray computed tomography (CT) and magnetic resonance imaging (MRI). In many cases, it is thus favorable
to have a contrast agent that can be used for multiple modalities
[259]. For UCNP-based imaging, multimodality imaging can be
enabled by modifying the chemical composition of the nanoparticles. The high atomic number of the RE elements can lead to
effective attenuation of X-rays. Thus, the use of UCNPs for Xray CT imaging has recently been proposed [260263]. Ytterbium
(Yb) and lutetium (Lu) have received most attention due to their
high atomic numbers. The UCNPs can be modified to provide
contrast in MRI by incorporating Gd3+ ions or a Gd3+ complex
into the crystal host or in a core-shell structure [264266]. Another approach for achieving MRI contrast is through the use of
iron oxide nanomaterials [255, 267270].
Doping or coating UCNPs with radioisotopes or agents can
endow them with radioactivity contrast for nuclear imaging, e.g.,
positron emission tomography [271274] and single-photon emission computed tomography [179, 275, 276]. Nuclear imaging has
a sensitivity comparable to that of fluorescence molecular imaging. However, nuclear imaging and fluorescence imaging can still
provide independent and complementary information, as the data
acquisition for these two modalities has a different time window.
3.3
Biomedical applications
3.3.1
Toxicity assessment
36
Upconverting nanoparticles
metal impurities [301], agglomeration and dispersion [302], degradation [303], and the formation of a protein corona2 [304, 305].
The toxicity of UCNPs has been less systematically investigated
than that of quantum dots and metal nanomaterials, due to their
relatively short history of existence. This topic has, however, attracted great interest.
The in vitro cytotoxicity of UCNPs with different surface modifications has been investigated in many studies, on a large number
of both human and animal cell lines [20, 238, 253, 272, 306, 307].
It was found that the cytotoxicity was very low, generally with cell
viability remaining above 80%, provided the nanoparticle concentration was kept below 100 g/ml [248].
Studies on the long-term effects of UCNPs on small animals have recently started to appear. Zhangs group investigated the biodistribution and clearance of polyethyleneiminecoated NaYF4 :Yb,Er nanoparticles in female Wistar rats following
an intravenous injection [20]. They found that the nanoparticles
were accumulated in the lungs immediately after injection, but
the amount of nanoparticles was significantly reduced in all tissues 24 h post-injection; the highest concentration being found in
the spleen. After 7 days, the nanoparticles were undetectable in
the rats. This group also studied the biocompatibility of silicacoated NaYF4 :Yb,Er nanoparticles using the same rat model, and
found accumulation in the lungs and the heart [248], and that
clearance time was significantly less than 7 days. The study by
Lis group on polyacrylic-acid-coated NaYF4 :Yb,Tm nanoparticles
in athymic nude mice showed the UCNPs to be mainly accumulated in the spleen and liver, with a clearance time longer than 7
days [308]. The same group recently studied the time-dependent
biodistribution of poly(ethylene glycol)-coated NaYF4 :Yb,Er labeled with a radioisotope (153 Sm) using single-photon emission
computed tomography imaging and -counter analysis [309], and
found that they had a long retention time in the blood, and were
partly eliminated through urinary excretion in vivo.
In addition to small animal studies, the toxicity of UCNPs
has also been assessed in other models. Yans group recently investigated the toxicity of polyethyleneimine-coated NaYF4 :Yb,Tm
nanoparticles in the roundworm Caenorhabditis elegans [310], and
found that feeding them 100 g UCNPs had no obvious toxic effect
based on observations on the growth and procreation. In a study
by Wang et al., on zebrafish embryos, they reported that concentrations of silica-coated LaF3 :Yb,Er nanoparticles under 100
g/ml had no obvious toxic effects, but concentrations above 200
g/ml led to chronic toxicity, which resulted in delayed hatching
and embryonic and larval development, and malformation [311].
2 A corona formed by adsorbed proteins surrounding a nanoparticle in
physiological environment.
37
3.3.2
Promising applications
38
Upconverting nanoparticles
In 2008, Zhou et al. demonstrated an upconversion luminescence switch using intermolecular energy transfer, employing
a diarylethene derivative and LaF3 :Yb,Ho nanoparticles loaded
in a poly(methyl methacrylate) film [332]. Since then, UCNPbased photoactivation has gained a considerable amount of interest. The group led by Branda achieved reversible ring-closing and
ring-opening reactions of dithienylethene photoswitches using NIR
light together with NaYF4 :Yb,Tm and NaYF4 :Yb,Er nanoparticles [285, 333, 334], as well as NIR-light-triggered release of a carboxylic acid from caged compounds3 [286]. Zhaos group reported NIR-triggered release of caged nitric oxide using UCNPs,
facilitating the therapeutic delivery of nitric oxide to physiological targets [335]. Zhangs group recently demonstrated the photorelease of small interfering RNA and plasmid DNA from caged
compounds in tissue phantoms with thicknesses up to 4 mm [284].
Xings group also reported the controlled photorelease of small
interfering RNA [336], and a system for the controlled uncaging
of D-luciferin and bioluminescence imaging based on D-luciferinconjugated UCNPs [337]. Other interesting progress in this field
can be found in recent publications [338342].
It is worth mentioning that the low quantum yield of UCNPs
is still a major limiting factor for applications such as these, especially when used in deep tissue, where the excitation fluence rate
is low. The use of such photoactivation systems in deep tissues
could be greatly facilitated by the adoption of pulsed excitation,
which will be discussed in detail in Chapter 4.
39
Chapter
4.1
The characterization of the QY of UCNPs, accounting for the excitation power density dependence, is of paramount importance,
both for materials development and practical applications. There
is currently no standard for characterizing and presenting the luminescence efficiency of UCNPs. Each group uses their own local
references, and compares the relative brightness of different sam41
4.1.1
Iem
,
n
Iex
(4.2)
which has the dimensions of [cm2 W1 ]n1 . However, for upconversion emission based on ETU, the n-power law for the luminescence intensity can become degraded even at moderate excitation
intensities.
In a non-saturating excitation regime, the power density dependence of the emission intensity of an n-photon upconversion
band can be expressed by:
n
Iem = k1 Iex
,
(4.3)
4.1.2
Due to the nonlinear dependence of the upconversion emission intensity on the excitation power density, the QY defined in percent
will be dependent on the power density, rather than constant [25].
In addition, the degradation of the power law caused by saturation
makes this dependence even more complex. For this reason, the
linear model for the QY characterization of two-photon upconversion emission, as demonstrated by Faulkner et al. [345], will fail
when saturation is present.
The saturation of upconversion emission is governed by the
competition between the ETU rate and the linear decay rate for the
depletion of the intermediate energy state involved in the upconversion process [343, 344]. This saturation feature has been used
to evaluate the quantum efficiency of UCNPs by many researchers
[346350], showing a small slope in the double-logarithmic plot,
and a small saturating power density indicating a high QY. However, the use of saturation to assess QY can lead to ambiguous
results unless a thorough theoretical analysis is performed.
Numerically solving the rate equations describing upconversion
processes and extracting the QY at steady state will, in principle,
provide accurate QY data for arbitrary excitation intensities, as
demonstrated in Paper V. However, this approach is extremely
demanding, and may not be feasible in practical applications, because too many parameters are required for the rate equations
and some of them are difficult to obtain. An analytical expression describing the excitation-intensity-dependent QY, with few
parameters, that can be easily obtained experimentally, will be
valuable, although the use of assumptions in achieving such an
expression would reduce the accuracy.
4.1.3
106
105
Intensity [a.u.]
104
103
102
101
100
10-1
10-2
102
103
104
105
Two-parameter characterization
Yb3+ -sensitized two-photon upconversion emission is the most efficient identified so far, and has attracted the greatest interest in
biomedical applications. A two-parameter model for QY characterization of such emission, accounting for the saturation property
of upconversion process, is proposed in Paper IV.
The mechanism of Yb3+ -sensitized two-photon upconversion
emission can be generally depicted by a quasi-three-level model for
the main upconversion activators, including Er3+ , Tm3+ and Ho3+
43
ETU2
2 (N2)
F5/2 (NYb1)
1 (N1)
ETU1
975 nm
dNYb1
NYb1
= NYb0
= 0,
dt
h
Yb1
dN1
N1
= C0 N0 NYb1 C1 N1 NYb1
= 0,
dt
1
N2
dN2
= C1 N1 NYb1
= 0.
dt
2
F7/2(NYb0)
Yb3+
Activator
0 (N0)
(4.4a)
(4.4b)
(4.4c)
= s
/b
,
1 + /b
(4.5)
with
s = C0 N0 Yb1 2 /2rad ,
and
b =
(4.6)
,
(4.7)
1 C1 Yb1 NYb0
where s is the maximum attainable QY of the system, b is the
balancing power density, reflecting the excitation intensity dependence of the QY, and 2rad is the radiative lifetime of state 2. At b ,
the ETU rate and the linear decay rate contribute equally to the
de-excitation of state 1, i.e., N1 /1 = C1 N1 NYb1 , and the slope
efficiency of the excitation-power-density dependence of the upconversion emission intensity on a log-log scale, expressed as (see
Paper IV):
k =1+
1
,
1 + 1 C1 Yb1 NYb0 /h
(4.8)
has a value of 1.5. The terms NYb0 , 1 and Yb1 , and C1 characterize the absorption of excitation light by Yb3+ ions, and the
retention times of the intermediate state of activator ions and the
excited state of Yb3+ ions, respectively, and the energy transfer
44
0.03
0.025
Quantum yield
rate between excited Yb3+ ions and activator ions. The balancing power density is thus a measure of the effectiveness of the
upconversion system; a lower value implying better performance.
Considering the fact that the upconversion emission is generated
through a stepwise process, the term s can be simply regarded
as the probability of reaching state 1 from state 0, while the term
/b
1+/b regarded as the probability of reaching state 2 from state 1.
Experimentally, s can be obtained by measuring the QY at very
high excitation intensity, and b can be extracted from routine
power dependence measurements of upconversion emission intensity.
Experimental QY data can be reasonably well fitted by Equation (4.5), as shown in Figure 4.3. The values of s and
b are estimated to be 0.91% and 3.8 W/cm2 , and 2.6% and
1.3 W/cm2 for the core (NaYF4 :Yb3+ /Tm3+ ) and core-shell
(NaYF4 :Yb3+ /Tm3+ @NaYF4 ) samples, respectively. The discrepancy could be due to the omission of certain energy transfer
processes in the rate equation, and the inaccuracy in the excitation beam size measurement involved in the measurement of the
excitation power density.
coreshell
core
0.02
0.015
0.01
0.005
0
0
0.5
1.5
x 10
4.1.4
Experimental measurements
Intensity (a.u.)
10
10
Gaussian
Uniform
10
10
4
10
10
Excitation power (W)
10
45
Quantum yield
0.08
0.06
0.04
0.02
0
0.2
0.4
Excitation power (W)
S2
S1
IC
Sn
S1
ISC
T1
T1
S0
S0
Donor
Acceptor
4.1.5
sity threshold by Monguzzi et al. [361]. This threshold is determined by a few parameters of the constituent molecules. Above
this threshold, triplet bimolecular annihilation becomes dominant
over spontaneous decay, resulting in efficient upconversion generation. Experimentally, the threshold is extracted by fitting the
power density dependence curve of upconversion emission in low
and high excitation intensity regimes; the threshold is then given
at the intersection of the fitted lines. Monguzzi et al. presented an
expression for the excitation-power-density-dependent intensity of
upconverted light [362].
As the TTA-based upconversion emission has a similar mechanism to that based on RE ions, the two-parameter QY characterization approach can be applied with some modifications. The
dynamics of a typical TTA system can be described by the following rate equations [361]:
TD
T
= (E)ISC Iex kD
TD ktr TD ,
t
(4.9a)
TA
T
= ktr TD kA
TA TT TA2 ,
(4.9b)
t
SA
1
S
= f TT TA2 kA
SA ,
(4.9c)
t
2
with the same notation as used by Monguzzi et al. [361], and where
ISC denotes the quantum efficiency of donor intersystem crossing.
Back energy transfer from acceptor to donor is not considered here.
The steady-state intensity can be expressed as:
s
#2
"
1
4(E)
k
I
f
TT
tr
ISC
ex
S
I = A kA
SA = A
,
1 + 1 +
T )2 (k + k T )
8
TT
(kA
tr
D
(4.10)
where A is the quantum efficiency of acceptor photoluminescence.
The slope efficiency on a double logarithmic scale can be described
by:
kTTA
dlogI
=1+ q
dlogIex
1+
1
4(E)TT ktr ISC Iex
T )2 (k +k T )
(kA
tr
D
(4.11)
When the contributions of spontaneous decay and bimolecular annihilation to the deactivation of the acceptor triplet states are
T
equal: kA
TA = TT TA2 , the excitation intensity threshold is given
by:
T
T 2
ktr + kD
2(kA
)
(
),
(4.12)
Ith =
(E)TT ISC
ktr
which gives a slope efficiency k = 4/3.
47
0.2
0.1
0
200
400
600
0.1
0.05
0
0
500
1000
1 A f ISC ktr
.
T )2
T
2 (kA
ktr + kD
(4.14)
sat
can
As can be seen in Equation (4.13), parameters Ith and TTA
be used as two standard indices to describe a TTA system, as they
fully characterize the QY, i.e., both the amplitude and the power
density dependence.
The QY data for PdOEP/DPA-based TTA systems reported
previously [357, 360] are well fitted by Equation (4.13), as shown
in Figure 4.7 and Figure 4.8, respectively. The parameters Ith and
sat
were estimated to be approximately 35 mW/cm2 and 27.5%,
TTA
and 5 mW/cm2 and 9% by fitting these two systems. Interestingly,
Equation (4.5) also gave a reasonably good fit to the reported QY
data, probably due to the similarity of the upconversion mechanisms in the ETU-based and TTA-based systems.
4.2
Upconversion emission can be excited using CW lasers or noncoherent lamps with excitation power densities well below 1
W/cm2 [132134], sometimes as low as 50 mW/cm2 [199]. This
is a considerable advantage compared to other types of anti-Stokes
emission such as second-harmonic generation and two-photon fluorescence, as it alleviates the need for high excitation intensity
allowing their use in deep biological tissues. However, the conventional CW excitation scheme is less optimal than pulsed excitation
from the perspective of energy conversion and possibility of suppressing the thermal effects of excitation light. Optimization of
the excitation scheme for UCNPs will be discussed in this section.
4.2.1
Due to the excitation power density dependence of the QY of upconverting materials, as shown in Figure 4.3 and Figure 4.7, higher
excitation intensity will lead to a higher QY. Considering a CW
source and a pulsed source that provide identical average power
density, a larger overall QY and thus a larger upconversion signal gain can be expected from the pulsed source, as the excitation
48
Quantum yield
10
10
10
10
10
10
10 15 20 25 30
Depth (mm)
Quantum yield
4.2.2
pulsed
CW
photons are confined to a narrow time window, resulting in an instantly higher excitation power density, as illustrated in Figure 4.9.
Upconversion processes are relatively slow with typical population times on the order of milliseconds [363, 364], because many
long-lived intermediate energy states of RE ions are involved. The
transient behavior of the upconversion system in each pulse period
must be taken into account when evaluating the performance of
pulsed excitation. The upconversion system has a lower QY at the
very beginning of each pulse than in a later phase at steady state,
since the energy states involved are far less populated. Thus, a
long rise time of upconversion emission, related to the lifetimes of
intermediate energy levels, would reduce some of the benefit of the
high excitation power density of pulsed excitation.
The applied average excitation power density is another critical
parameter determining the signal gain. If it is already located in
a saturating regime, no gain in signal will be obtained by using
pulsed excitation, relative to equivalent CW excitation.
It is worth mentioning that pulsed excitation is even better for
multi-photon upconversion emission such as blue and ultraviolet
emission from Tm3+ ions, due to their high-order dependence on
excitation power density. Such short wavelength emission has been
shown to be very promising in the photoactivation of biomolecules.
It is expected that the pulsed excitation approach will facilitate the
implementation of NIR-triggered photoactivation in deep tissues,
while the depth is limited to only a few millimeters at present
[284].
The pulsed excitation approach is described in detail and compared with CW excitation in Paper V.
Intensity
10
2
10
3
10
4
10
5
10
6
10
5
10 15 20 25 30
Depth (mm)
49
Apart from approaches to improve the efficiency of the materials, such as surface modification with optical antennas including
fluorescent dyes [205] and noble metal nanostructures [219], another solution to enhance the upconversion emission is to increase
the excitation light intensity. However, this is not possible when
using CW excitation due to the risk of tissue damage. This is
regulated by the ANSI standard as described below.
MPE ( W/cm )
1.5
1
0.5
0
500
1000
(nm)
1500
4.2.3
ANSI standard
50
Wavelength
(m)
Exposure duration, t
(s)
0.400-1.400
0.400-1.400
0.400-1.400
109 to 107
107 to 10
10 to 3104
(Jcm
MPE
)
(Wcm2 )
2C A 102
1.1CA t 0.25
0.2C A
tects against sub-threshold pulse-cumulative thermal injury. Figure 4.12 presents the MPEs for skin exposure to CW lasers with
different wavelengths. As can be seen, for a CW 975 nm laser, the
MPE is around 710 mW/cm2 , and around 320 mW/cm2 for 800
nm. It should be noted that the ANSI standard is rather conservative. Considering the optical properties of tissues at these two
wavelengths, as discussed in Chapter 2, a higher MPE should be
allowed for 800 nm radiation.
MPE ( W/cm2 )
10
10
10
4.2.4
Single-shot imaging
1.1CA t0.25
.
t
10
t (ms)
10
(4.15)
51
Chapter
5.1
In the field of tissue optics, the aim of solving the inverse problem is to extract the interior optical properties of a medium noninvasively through measurements performed on the tissue boundaries, as illustrated in Figure 5.1. The optical properties in question include the absorption, a [cm1 ], and reduced scattering
coefficients, 0s [cm1 ], the diffusion coefficient D [cm], the fluorescence lifetime, [s], and the fluorescence yield (i.e., the product
of the QY and the absorption coefficient of the fluorophore), af
53
(a)
Iexc
a , s'
5.1.1
(b)
Measurables
Iexc
54
Excitation schemes
A point-like source with a CW output is typically used to provide
the excitation energy, while point-like short-pulse laser sources
are used for time-resolved measurements [371]. Spatially modulated illumination, developed by Cuccia et al. [389, 390], has
also been exploited and has been found to be useful in alleviating
the ill-posedness of the optical tomography problem [391395]. In
the frequency domain technique, intensity-modulated laser sources
are utilized, irradiating the tissue with sinusoidally varying source
power [386, 387].
(a)
(b)
Imaging geometries
Imaging geometries can be divided into three main categories: epiillumination, transillumination, and projection geometries. In the
epi-illumination geometry (Figure 5.2(a)), with the source and the
detector placed on the same side of the tissue, the sensitivity at
shallow depths is superior to that at greater depths. The fluorescence signal is usually considerable if fluorophores are not deeply
located, as the fluorescent light will propagate only a relatively
short distance before exiting from the surfaces of the tissue. The
background from the excitation light is often relatively high, for
the same reason. Some large tissues only allow this geometry due
to their size, e.g., in fluorescence tomography of most parts of the
human body. In the transillumination geometry (Figure 5.2(b)),
the source and detector are placed on opposite sides of the tissue.
The imaging sensitivity is then higher near the source and the detector than deep in the tissue. The fluorescent light is generally
weak, as is the excitation background, due to the attenuation of
light by the tissue. In the projection geometry (Figure 5.2(c)),
either the source-detector pair or the tissue can be rotated, usually by 360 , yielding a relatively uniform sensitivity throughout
the whole tissue.
The appropriate geometry will depend on the application. In
very deep tissue, such as the human prostate, other geometries may
be necessary, e.g. interstitial measurements, using fibres inserted
into the tissue to deliver and collect the light [396].
5.1.2
(c)
In the framework of FDOT, the measurable, denoted by y, is compared with the modeled value F (x), calculated by a forward model.
The forward model is given by:
F (x) = M[(x)],
(5.1)
If
(5.2)
The solution is obtained by finding the value of x giving the minimal value of 2 , using least-squares minimization.
Due to the high scattering of biological tissues, the problem at
hand is fundamentally ill-posed, i.e., the solution is non-unique,
and a small perturbation in the data can therefore cause a significant alteration in the solution with minimal 2 . The origin of
the ill-posedness caused by high scattering can be understood by
considering the following two facts. (1) A small, deeply located
fluorophore could render the same surface fluorescence exitance as
a large, shallow one, as illustrated in Figure 5.3; and (2) a perturbation in a small volume inside the tissue causes signal changes
in many source-detector projections. In addition, this problem is
under-determined, as the number of unknowns usually exceeds the
number of measurements. Even if a large data set is obtained, for
example, by employing a charge-coupled device as the detector,
the problem remains, because only linearly independent measurements contribute to the constraint to the solution.
The ill-posedness and under-determination discussed above
make the inversion process unstable, causing the solution to fluctuate considerably. In order to stabilize and smooth the solution, a
Tikhonov regularization term is generally included in the objective
function [397]:
2 = ky F (x)k2 + kL(x x0 )k2 .
(5.3)
Here, L is a regularization matrix, and is the Tikhonov regularization parameter, giving the weight of the solution norm in the
objective function. x0 is the initial estimate of the fluorophore distribution, often set to zero. The fluorophore distribution is thus
reconstructed by minimizing the objective function. The regularization may be performed in different ways, and its meaning will
be discussed in detail in Section 5.1.4. Adopting a Euclidean norm
provides a solution when the derivative with respect to x is zero:
T
F
2
= 2
[y F (x)] + 2LT L(x x0 ) = 0.
x
x
(5.4)
F
x
T
[y F (xk+1 )] LT L(xk+1 x0 ) = 0.
(5.5)
F (xk+1 ) is expanded in Taylor series around xk and only the firstorder terms are retained, i.e.:
F (xk+1 ) = F (xk ) +
F
(xk+1 xk ),
x
(5.6)
(5.7)
where J = F
x is the Jacobian matrix, x = xk+1 xk , is the
update for the fluorophore distribution, and y = yF (xk ) is the
difference or misfit between the data and the model in the current
iteration. A Levenberg-Marquardt procedure is often applied in
the iteration, i.e., assuming x = xk x0 [398, 399], which leads
to:
T L]1 JT y.
x = [JT J + L
(5.8)
= 2.
Note
5.1.3
Explicit expression
The Jacobian introduced in Equation 5.7 gives the rate of change
of the modeled values relative to the parameters. In the calculation
of the Jacobian, the adjoint method is employed in order to reduce
computational cost [400]. In addition, the Born approximation is
usually applied, which neglects the influence of the fluorophore
distribution on the excitation field and the emission field [370].
This is valid for low fluorophore concentrations [375]. The explicit
expression of the Jacobian can thus be described by:
J=
=
(Cx xm ) = Cx m .
x
x
(5.9)
5.1.4 Regularization
The Jacobian is thus normalized in the same way. Such modifications afford advantages such as avoiding absolute photon-field
measurements and certain robustness to heterogeneities in the optical properties [402].
Singular-value analysis
Singular-value decomposition of the Jacobian is often performed
for the purpose of analyzing and optimizing an experimental setup
[403406]. This factorization yields:
J = USVT ,
(5.10)
5.1.4
Regularization
The purpose of regularization is to both improve the spatial resolution and reduce the noise in FDOT reconstruction. The difference
in the sensitivities at different depths can also be equilibrated by
manipulating the regularization. Regularization is controlled by
the regularization parameter and the regularization matrix.
The regularization parameter
The regularization parameter, , controls the weight given to minimization of the solution norm relative to minimization of the residual norm in Equation (5.3). A large favors a small solution norm
although the residual norm is large, while a small allows a large
solution norm. In addition, also controls the sensitivity of the
regularized solution to perturbations in the weight matrix and the
data [408, 409]. Regularization can be viewed as a smoothing of
the solution. Generally, a large regularization parameter results
in a smooth image with low spatial resolution, while a small regularization parameter can sharpen the solution with noise of high
spatial frequency components, as shown in Figure 5.5.
58
Target
=10-6
=10-5
=10
-1
=10
=10
10
10
10
10
10
10
Residual norm y-F(x)2
59
5.1.4 Regularization
Li,j
1,
if i = j,
(5.11)
60
within the same tissue region should have similar optical properties, and thus the sought quantity varies smoothly.
Multispectral fluorescent emission recordings provide depth information about the fluorophores embedded in tissues [379, 380,
425]. Axelsson et al. presented a method of creating a spatially
varying regularization matrix using the intensity ratio of two emission wavelengths [381]. This approach was also employed in the
present work (Paper IX) with UCNP-mediated FDOT. Both improved resolution and contrast were obtained in the reconstructed
images, compared to using standard Tikhonov regularization.
Fluorescent dyes and quantum dots are commonly used as contrast agents in FDOT. Recently, UCNPs have been introduced
as contrast agents in FDOT [282]. UCNPs have many benefits
in FDOT. Both the excitation and emission wavelengths can be
tuned to the tissue optical window, where biological tissue has relatively low absorption and scattering. This enables imaging to
a considerable depth. Upconversion emissions from UCNPs are
anti-Stokes shifted, and this makes the complete separation of upconversion emission from tissue autofluorescence possible by using
an appropriate filter system, thus achieving autofluorescence-free
imaging, as demonstrated in Paper VI. As FDOT is susceptible
to noise, the suppression of background tissue autofluorescence
leads to reconstructions with high quality [282]. In addition, the
nonlinear power dependence of upconversion is also advantageous
in FDOT, as will be discussed in detail below. Based on these
factors, high-resolution FDOT has been developed using NIRemitting Yb3+ /Tm3+ -codoped nanoparticles.
5.2.1
NaYF4:Yb3+/Tm3+@NaYF4+
Intensity+ [a.u.]
5.2
50 nm
JCPDS+281192
10
20
30
2 []
Figure 5.8. X-ray diffraction
pattern of the
NaYF4 :Yb3+ /Tm3+ @NaYF4
UCNPs used as the contrast
agents in FDOT.
The nanoparticles used in this study were synthesized primarily using the oleate route proposed by the group of Zhang [199, 228, 426]
with slight modifications, except in Paper II. The particles have a
core-shell structure, and the average size can be tuned from 20 nm
to 50 nm by varying reaction parameters. Figure 5.7 shows a typical transmission electron microscopy image. The particles appear
spherical in shape. The prepared UCNPs have a hexagonal phase,
with an X-ray diffraction pattern well consistent with the JCPDS
standard card (28-1192), as shown in Figure 5.8. Oleic acid was
used as surfactant in the synthesis. The prepared nanoparticles
could be well dispersed in commonly used nonpolar solvents, such
61
quadratic
5.2.2
Selective excitation
5.2.3
62
by acquiring a third image under excitation by two beams simultaneously at the two positions, as the emission intensity, , can be
expressed as:
(x,1 + x,2 )2 af m
= (2x,1 + 2x,2 + 2x,1 x,2 )af m ,
(5.12)
5.2.4
Multispectral regularization
5.2.5 Outlook
5.2.5
Outlook
wavelength, i.e., certain emission bands are generated under excitation at one wavelength, and the response to excitation power
density, i.e., the nonlinear excitation-power-density dependence of
upconversion emission intensity. As an optical system, UCNPs are
relatively complex, as many steps are involved in the generation
of upconversion emission, and some of them are rather slow,
characterized by the long lifetimes of the energy levels of RE ions.
The response of UCNPs to time-varying excitation light could
provide contrast or other interesting information for biomedical
imaging on an easy-to-operate time scale, i.e., sms. This could
be an interesting topic of research.
The use of UCNPs for biomedical imaging is still in its infancy,
and most work so far has been for the proof of concept. Functional
imaging using UCNPs in different biological models is becoming
and will remain a hot research topic in next stage.
65
67
68
69
70
71
Acknowledgements
First of all, I would like to express my deepest gratitude to
my supervisor, Professor Stefan Andersson-Engels, who was my
mentor throughout this endeavor. I still remember how green I
was when I started my PhD studies, but thanks to him, I now
feel much more confident in scientific research. His patience,
guidance and encouragement have been indispensable. Thank
you for encouraging me to pursue my own ideas, and for giving
me a little nudge when necessary. Thank you also for all
the illuminating discussions we have had, especially when I was
confused or frustrated. Your kindness and support will never be
forgotten.
I would also like to express my thanks to my co-supervisor,
Professor Thomas Laurell, for giving me the opportunity to dabble in the fantastic field of microfluidics, and for his wonderful
support and rewarding discussions throughout the period of these
studies.
Special thanks also to my mentor, Dr. Gabriel Somesfalean.
I am truly grateful to you for introducing me to a splendid world,
for your invaluable help in research throughout the years, and for
being a good friend.
Sincere thanks also to Professor Peter E. Andersen, for fruitful
collaboration, and for organizing the wonderful summer school at
Hven.
I would also like to express my gratitude to the Svanberg
family. Thank you all for your warm welcome when I set foot
in a foreign land. Sune, thank you for organizing the very nice
apartment. I will never forget the night when you painted the
floor. You painted my family in the bedroom, and painted yourself
out, in order to make our living environment even nicer. Thank
you also for your truly inspiring and contagious enthusiasm for
life and science. Katarina, thank you for always being friendly
and helpful, in both professional and personal matters. Emilie
73
Acknowledgements
Acknowledgements
75
References
[1].
[2].
[3].
[4].
[5].
[6].
[7].
[8].
M. Andresen, A. C. Stiel, J. F
olling, D. Wenzel, A. Sch
onle, A. Egner,
C. Eggeling, S. W. Hell and S. Jakobs. Photoswitchable fluorescent
proteins enable monochromatic multilabel imaging and dual color fluorescence nanoscopy. Nat. Biotechnol. 26, 10351040 (2008).
[9].
[10].
[11].
[12].
[13].
[14].
Y. Zhong, N. Kaji, M. Tokeshi and Y. Baba. Nanobiotechnology: quantum dots in bioimaging. Expert Rev. Proteomic. 4, 565572 (2007).
77
References
78
[15].
[16].
Y. Xing and J. H. Rao. Quantum dot bioconjugates for in vitro diagnostics & in vivo imaging. Cancer Biomark. 4, 307319 (2008).
[17].
[18].
S. C. Davis, B. W. Pogue, H. Dehghani and K. D. Paulsen. Contrastdetail analysis characterizing diffuse optical fluorescence tomography
image reconstruction. J. Biomed. Opt. 10, 050501 (2005).
[19].
[20].
[21].
[22].
[23].
J. Zhou, Z. Liu and F. Li. Upconversion nanophosphors for smallanimal imaging. Chem. Soc. Rev. 41, 13231349 (2012).
[24].
P. Svenmarker, C. T. Xu and S. Andersson-Engels. Use of nonlinear upconverting nanoparticles provides increased spatial resolution in
fluorescence diffuse imaging. Opt. Lett. 35, 27892791 (2010).
[25].
[26].
A. H. Li, Q. L
u, Z. R. Zheng, L. Sun, W. Z. Wu, W. L. Liu, Y. Q. Yang
and T. Q. L
u. Yellow-green upconversion luminescence of Dy3+ ion in
LiNbO3 crystal heavily codoped with ZnO. J. Appl. Phys. 102, 113102
(2007).
[27].
A. H. Li, Z. R. Zheng, Q. L
u, L. Sun, Y. H. Xu, W. L. Liu, W. Z. Wu,
Y. Q. Yang and T. Q. L
u. Sensitized holmium upconversion emission
in LiNbO3 triply doped with Ho3+ , Yb3+ , and Nd3+ . J. Appl. Phys.
104, 063526 (2008).
[28].
[29].
C. F. Gainer, G. S. Joshua and M. Romanowski. Toward the use of twocolor emission control in upconverting NaYF4 : Er3+ ,Yb3+ nanoparticles for biomedical imaging. Proc. SPIE 8231, 82310I (2012).
[30].
W. Zhang, F. Ding and S. Y. Chou. Large enhancement of upconversion luminescence of NaYF4 :Yb3+ /Er3+ nanocrystal by 3D plasmonic
nano-antennas. Adv. Mater. 24, OP236OP241 (2012).
[31].
[32].
References
[33].
G. Mie. Beitr
age zur optik tr
uber medien, speziell kolloidaler metall
osungen. Ann. Phys. 330, 377445 (1908).
[34].
[35].
[36].
[37].
[38].
[39].
[40].
[41].
[42].
[43].
[44].
[45].
[46].
S. J. Matcher and C. E. Cooper. Absolute quantification of deoxyhaemoglobin concentration in tissue near infrared spectroscopy. Phys.
Med. Biol. 39, 1295 (1994).
[47].
[48].
[49].
[50].
[51].
[52].
79
References
80
[53].
[54].
[55].
[56].
[57].
[58].
[59].
[60].
B. R. Masters, P. T. So and E. Gratton. Multiphoton excitation fluorescence microscopy and spectroscopy of in vivo human skin. Biophys.
J. 72, 24052412 (1997).
[61].
R. Richards-Kortum and E. Sevick-Muraca. Quantitative optical spectroscopy for tissue diagnosis. Annu. Rev. Phys. Chem. 47, 555606
(1996).
[62].
M. J. Koehler, K. K
onig, P. Elsner, R. B
uckle and M. Kaatz. In vivo
assessment of human skin aging by multiphoton laser scanning tomography. Opt. Lett. 31, 28792881 (2006).
[63].
I. Georgakoudi, B. C. Jacobson, M. G. M
uller, E. E. Sheets,
K. Badizadegan, D. L. Carr-Locke, C. P. Crum, C. W. Boone, R. R.
Dasari, J. Van Dam and M. S. Feld. NAD(P)H and collagen as in vivo
quantitative fluorescent biomarkers of epithelial precancerous changes.
Cancer Res. 62, 682687 (2002).
[64].
[65].
G. R. Cherrick, S. W. Stein, C. M. Leevy and C. S. Davidson. Indocyanine green - Observations on its physical properties, plasma decay,
and hepatic extraction. J. Clin. Invest. 39, 592600 (1960).
[66].
[67].
[68].
L. Ye, K.-T. Yong, L. Liu, I. Roy, R. Hu, J. Zhu, H. Cai, W.-C. Law,
J. Liu, K. Wang, J. Liu, Y. Liu, Y. Hu, X. Zhang, M. T. Swihart and
P. N. Prasad. A pilot study in non-human primates shows no adverse
response to intravenous injection of quantum dots. Nat. Nano. 7,
453458 (2012).
[69].
Y. Su, M. Hu, C. Fan, Y. He, Q. Li, W. Li, L.-H. Wang, P. Shen and
Q. Huang. The cytotoxicity of CdTe quantum dots and the relative
contributions from released cadmium ions and nanoparticle properties.
References
[71].
[72].
[73].
[74].
[75].
R. L. Longini and R. Zdrojkowski. A note on the theory of backscattering of light by living tissue. IEEE Transactions on Biomedical Engineering BME-15, 410 (1968).
[76].
L. Reynolds, C. Johnson and A. Ishimaru. Diffuse reflectance from a finite blood medium: applications to the modeling of fiber optic catheters.
Appl. Opt. 15, 20592067 (1976).
[77].
A. Ishimaru. Wave propagation and scattering inrandom media. Academic Press (1978).
[78].
T. Svensson. Pharmaceutical and biomedical applications of spectroscopy in the photon migration regime. PhD thesis Lund University
(2008).
[79].
S. R. Arridge and M. Schweiger. Image reconstruction in optical tomography. Phil. Trans. R. Soc. Lond. B 352, 717726 (1997).
[80].
[81].
[82].
[83].
[84].
J. B. Domnguez and Y. B
erub
e-Lauzi
ere. Diffuse light propagation in
biological media by a time-domain parabolic simplified spherical harmonics approximation with ray-divergence effects. Appl. Opt. 49,
14141429 (2010).
[85].
[86].
A. Liemert and A. Kienle. Analytical solutions of the simplified spherical harmonics equations. Opt. Lett. 35, 35073509 (2010).
[87].
[88].
[89].
81
References
[90].
E. Alerstam. Optical spectroscopy of turbid media: time-domain measurements and accelerated Monte Carlo modelling. PhD thesis Lund
University (2011).
[91].
E. Alerstam. Anisotropic diffusive transport: a rigorous theory for connecting microscopic scattering and macroscopic transport properties.
arXiv:1311.3603v3 (2013).
[92].
M. Schweiger, S. R. Arridge, M. Hiraoka and D. T. Delpy. The finite element method for the propagation of light in scattering media:
boundary and source conditions. Med. Phys. 22, 17791792 (1995).
[93].
[94].
L. Mart-L
opez, J. C. Hebden and J. B. Domnguez. Estimates of
minimum pulse width and maximum modulation frequency for diffusion
optical tomography. Opt. Lasers Eng. 44, 11721184 (2006).
[95].
[96].
[97].
T. Durduran, A. G. Yodh, B. Chance and D. A. Boas. Does the photondiffusion coefficient depend on absorption? J. Opt. Soc. Am. A 14,
33583365 (1997).
[98].
[99].
82
References
[108]. F. Martelli, D. Contini, A. Taddeucci and G. Zaccanti. Photon migration through a turbid slab described by a model based on diffusion
approximation. II. Comparison with Monte Carlo results. Appl. Opt.
36, 46004612 (1997).
[109]. E. Alerstam, S. Andersson-Engels and T. Svensson. White Monte Carlo
for time-resolved photon migration. J. Biomed. Opt. 13, 041304 (2008).
[110]. E. Alerstam, T. Svensson and S. Andersson-Engels. Parallel computing
with graphics processing units for high-speed Monte Carlo simulation
of photon migration. J. Biomed. Opt. 13, 060504 (2008).
[111]. S. A. Prahl, M. Keijzer, S. L. Jacques and A. J. Welch. A Monte Carlo
model of light propagation in tissue. IS 5, 102111 (1989).
[112]. J. Swartling, A. Pifferi, A. M. K. Enejder and S. Andersson-Engels.
Accelerated Monte Carlo models to simulate fluorescence spectra from
layered tissues. J. Opt. Soc. Am. A 20, 714727 (2003).
[113]. X. Teng, Y. Zhu, W. Wei, S. Wang, J. Huang, R. Naccache, W. Hu,
A. I. Y. Tok, Y. Han, Q. Zhang, Q. Fan, W. Huang, J. A. Capobianco
and L. Huang. Lanthanide-doped Nax ScF3+x nanocrystals: Crystal
structure evolution and multicolor tuning. J. Am. Chem. Soc. 134,
83408343 (2012).
[114]. Q. M. Huang, J. C. Yu, E. Ma and K. M. Lin. Synthesis and characterization of highly efficient near-infrared upconversion Sc3+ /Er3+ /Yb3+
tridoped NaYF4 . J. Phys. Chem. C 114, 47194724 (2010).
[115]. D. Hartree. The Calculation of atomic Structures. John Wiley & Sons
Ltd., New York (1957).
[116]. J. C. Slater. Quantum theory of atomic structure. McGraw-Hill, New
York (1960).
[117]. X. Xu and M. Su, editors. Luminescence and luminescent materials
(in Chinese). Chemical Industry Press (China) (2004).
[118]. G. H. Dieke and H. M. Crosswhite. The spectra of the doubly and triply
Ionized rare earths. Appl. Opt. 2, 675686 (1963).
[119]. W. T. Carnall, G. L. Goodman, K. Rajnak and R. S. Rana. A systematic analysis of the spectra of the lanthanides doped Into single-crystal
LaF3 . J. Chem. Phys. 90, 34433457 (1989).
[120]. G. Liu. Spectroscopic properties of rare earths in optical materials
chapter Electronic energy level structure. Springer (2005).
[121]. B. Henderson and G. F. Imbusch. Optical Spectroscopy for Inorganic
Solids. Oxford Univ Press Inc (1989).
[122]. G. W. Burdick, S. M. Crooks and M. F. Reid. Ambiguities in
the parametrization of 4f N 4f N electric-dipole transition intensities.
Phys. Rev. B 59, R7789R7792 (1999).
[123]. D. J. Newman and G. Balasubramanian. Parametrization of rare-earth
ion transition intensities. J. Phys. C 8, 3744 (1975).
[124]. M. F. Reid and F. S. Richardson. Lanthanide 4f4f electric dipole
intensity theory. J. Phys. Chem. 88, 35793586 (1984).
[125]. B. R. Judd. Optical absorption intensities of rare-earth ions. Phys.
Rev. 127, 750761 (1962).
[126]. G. S. Ofelt. Intensities of crystal spectra of rare-earth ions. J. Chem.
Phys. 37, 511520 (1962).
[127]. M. P. Hehlen, M. G. Brik and K. W. Kramer. 50th anniversary of the
Judd-Ofelt theory: An experimentalists view of the formalism and its
application. J. Lumin. 136, 221239 (2013).
[128]. M. F. Reid. Spectroscopic properties of rare earths in optical materials
chapter Transition Intensities. Springer (2005).
83
References
84
References
[149]. Y.-F. Wang, G.-Y. Liu, L.-D. Sun, J.-W. Xiao, J.-C. Zhou and C.-H.
Yan. Nd3+ -sensitized upconversion nanophosphors: efficient in vivo
bioimaging probes with minimized heating effect. ACS Nano 7, 7200
7206 (2013).
[150]. F. W. Ostermayer and L. G. Van Uitert. Cooperative energy transfer
from Yb 3+ to Tb 3+ in YF3 . Phys. Rev. B 1, 42084212 (1970).
[151]. H. Liang, G. Chen, L. Li, Y. Liu, F. Qin and Z. Zhang. Upconversion
luminescence in Yb3+ /Tb3+ -codoped monodisperse NaYF4 nanocrystals. Opt. Commun. 282, 3028 3031 (2009).
[152]. E. Nakazawa and S. Shionoya. Cooperative luminescence in YbPO4 .
Phys. Rev. Lett. 25, 17101712 (1970).
[153]. F. Auzel, D. Meichenin, F. Pell
e and P. Goldner. Cooperative luminescence as a defining process for RE-ions clustering in glasses and
crystals. Opt. Mater. 4, 3541 (1994).
[154]. J. S. Chivian, W. E. Case and D. D. Eden. The photon avalanche:
A new phenomenon in Pr3+ -based infrared quantum counters. Appl.
Phys. Lett. 35, 124125 (1979).
[155]. M. F. Joubert, S. Guy and B. Jacquier. Model of the photon-avalanche
effect. Phys. Rev. B 48, 1003110037 (1993).
[156]. M. F. Joubert. Photon avalanche upconversion in rare earth laser materials. Opt. Mater. 11, 181203 (1999).
[157]. R. Scheps. Upconversion laser processes. Prog. Quantum. Electron. 20,
271 358 (1996).
[158]. Ph. Goldner and F. Pell
e. Photon avalanche fluorescence and lasers.
Opt. Mater. 5, 239249 (1996).
[159]. S. Sivakumar, F. C. J. M van Veggel and P. S. May. Near-Infrared
(NIR) to red and green up-conversion emission from silica sol-gel thin
films made with La0.45 Yb0.50 Er0.05 F3 nanoparticles, hetero-loopingenhanced energy transfer (Hetero-LEET): A new up-conversion process. J. Am. Chem. Soc. 129, 620625 (2007).
[160]. G. Chen, H. Liang, H. Liu, G. Somesfalean and Z. Zhang. Anomalous power dependence of upconversion emissions in Gd2 O3 : Er3+
nanocrystals under diode laser excitation of 970 nm. J. Appl. Phys.
105, 114315 (2009).
[161]. E. Downing, L. Hesselink, J. Ralston and R. Macfarlane. A three-color,
solid-state, three-dimensional display. Science 273, 11851189 (1996).
[162]. G. Blasse and B. C. Grabmaier. Luminescent materials. SpringerVerlag Berlin Heidelberg (1994).
[163]. F. Wang and X. Liu. Recent advances in the chemistry of lanthanidedoped upconversion nanocrystals. Chem. Soc. Rev. 38, 976989 (2009).
[164]. J. M. F. Vandijk and M. F. H. Schuurmans. On the nonradiative and
radiative decay-rates and a modified exponential energy-gap law for 4f4f transitions in rare-earth ions. J. Chem. Phys. 78, 53175323 (1983).
[165]. K. W. Kr
amer, D. Biner, G. Frei, H. U. G
udel, M. P. Hehlen and S. R.
L
uthi. Hexagonal sodium yttrium fluoride based green and blue emitting
upconversion phosphors. Chem. Mater. 16, 12441251 (2004).
[166]. S. Heer, K. K
ompe, H. U. G
udel and M. Haase. Highly efficient multicolour upconversion emission in transparent colloids of lanthanidedoped NaYF4 nanocrystals. Adv. Mater. 16, 21022105 (2004).
[167]. D. K. Chatterjee, M. K. Gnanasammandhan and Y. Zhang. Small upconverting fluorescent nanoparticles for biomedical applications. Small
6, 27812795 (2010).
[168]. A. Gnach and A. Bednarkiewicz. Lanthanide-doped up-converting
nanoparticles: Merits and challenges. Nano Today 7, 532563 (2012).
85
References
[169]. G. Wang, Q. Peng and Y. Li. Upconversion luminescence of monodisperse CaF2 :Yb3+ /Er3+ nanocrystals. J. Am. Chem. Soc. 131, 14200
14201 (2009).
[170]. M. Ito, C. Goutaudier, Y. Guyot, K. Lebbou, T. Fukuda and G. Boulon.
Crystal growth, Yb3+ spectroscopy, concentration quenching analysis
and potentiality of laser emission in Ca1X YbX F2+X . J. Phys.: Condens. Matter 16, 15011521 (2004).
[171]. N.-N. Dong, M. Pedroni, F. Piccinelli, G. Conti, A. Sbarbati, J. E.
Ramrez-Hern
andez, L. M. Maestro, M. C. Iglesias-de la Cruz, F. SanzRodriguez, A. Juarranz, F. Chen, F. Vetrone, J. A. Capobianco, J. G.
Sol
e, M. Bettinelli, D. Jaque and A. Speghini. NIR-to-NIR two-photon
excited CaF2 :Tm3+ ,Yb3+ nanoparticles: multifunctional nanoprobes
for highly penetrating fluorescence bio-imaging. ACS Nano 5, 8665
8671 (2011).
[172]. G. Chen, J. Shen, T. Y. Ohulchanskyy, N. J. Patel, A. Kutikov,
Z. Li, J. Song, R. K. Pandey, H.
Agren, P. N. Prasad and G. Han.
(-NaYbF4 :Tm3+ )/CaF2 core/shell nanoparticles with efficient nearInfrared to near-infrared upconversion for high-contrast deep tissue
bioimaging. ACS Nano 6, 82808287 (2012).
[173]. J. Shen, G. Chen, T. Y. Ohulchanskyy, S. J. Kesseli, S. Buchholz,
Z. Li, P. N. Prasad and G. Han. Tunable near infrared to ultraviolet
upconversion luminescence enhancement in (-NaYbF4 :Tm3+ )/CaF2
core/shell nanoparticles for in situ real-time recorded biocompatible
photoactivation. Small 9, 32133217 (2013).
[174]. F. Shi, J. Wang, X. Zhai, D. Zhao and W. Qin. Facile synthesis of
-NaLuF4 : Yb/Tm hexagonal nanoplates with intense ultraviolet upconversion luminescence. CrystEngComm 13, 37823787 (2011).
[175]. Q. Liu, Y. Sun, T. Yang, W. Feng, C. Li and F. Li. Sub-10 nm hexagonal lanthanide-doped NaLuF4 upconversion nanocrystals for sensitive
bioimaging in vivo. J. Am. Chem. Soc. 133, 1712217125 (2011).
[176]. J. Zhou, X. Zhu, M. Chen, Y. Sun and F. Li. Water-stable NaLuF4 based upconversion nanophosphors with long-term validity for multimodal lymphatic imaging. Biomaterials 33, 62016210 (2012).
[177]. T. S. Yang, Y. Sun, Q. Liu, W. Feng, P. Y. Yang and F. Y. Li. Cubic sub-20 nm NaLuF4 -based upconversion nanophosphors for highcontrast bioimaging in different animal species. Biomaterials 33, 3733
3742 (2012).
[178]. S. Zeng, J. Xiao, Q. Yang and J. Hao.
Bi-functional
NaLuF4 :Gd3+ /Yb3+ /Tm3+ nanocrystals: structure controlled synthesis, near-infrared upconversion emission and tunable magnetic properties. J. Mater. Chem. 22, 98709874 (2012).
[179]. Y. Yang, Y. Sun, T. Cao, J. Peng, Y. Liu, Y. Wu, W. Feng, Y. Zhang
and F. Li. Hydrothermal synthesis of NaLuF4 : 153 Sm,Yb,Tm nanoparticles and their application in dual-modality upconversion luminescence
and SPECT bioimaging. Biomaterials 34, 774 783 (2013).
[180]. N. Niu, P. Yang, F. He, X. Zhang, S. Gai, C. Li and J. Lin.
Tunable multicolor and bright white emission of one-dimensional
NaLuF4 :Yb3+ ,Ln3+ (Ln = Er, Tm, Ho, Er/Tm, Tm/Ho) microstructures. J. Mater. Chem. 22, 1088910899 (2012).
[181]. Y. Li, J. Zhang, Y. Luo, X. Zhang, Z. Hao and X. Wang. Color control
and white light generation of upconversion luminescence by operating
dopant concentrations and pump densities in Yb3+ , Er3+ and Tm3+
tri-doped Lu2 O3 nanocrystals. J. Mater. Chem. 21, 28952900 (2011).
[182]. F. Wang, Y. Han, C. S. Lim, Y. Lu, J. Wang, J. Xu, H. Chen, C. Zhang,
M. Hong and X. Liu. Simultaneous phase and size control of upconversion nanocrystals through lanthanide doping. Nature 463, 10611065
86
References
(2010).
[183]. F. Wang and X. Liu. Upconversion multicolor fine-tuning: Visible to
near-infrared emission from lanthanide-doped NaYF4 nanoparticles. J.
Am. Chem. Soc. 130, 56425643 (2008).
[184]. G. Chen, Y. Liu, Y. Zhang, G. Somesfalean, Z. Zhang, Q. Sun and
F. Wang. Bright white upconversion luminescence in rare-earth-iondoped Y2 O3 nanocrystals. Appl. Phys. Lett. 91, 133103 (2007).
[185]. G. Chen, H. Liu, G. Somesfalean, H. Liang and Z. Zhang. Upconversion emission tuning from green to red in Yb3+ /Ho3+ -codoped NaYF4
nanocrystals by tridoping with Ce3+ ions. Nanotechnology 20, 385704
(2009).
[186]. N. Rakov, G. S. Maciel, M. L. Sundheimer, L. de S. Menezes, A. S. L.
Gomes, Y. Messaddeq, F. C. Cassanjes, G. Poirier and S. J. L. Ribeiro.
Blue upconversion enhancement by a factor of 200 in Tm3+ -doped tellurite glass by codoping with Nd3+ ions. J. Appl. Phys. 92, 63376339
(2002).
[187]. T. S. Hauck, A. A. Ghazani and W. C. W. Chan. Assessing the effect of
surface chemistry on gold nanorod uptake, toxicity, and gene expression
in mammalian cells. Small 4, 153159 (2008).
[188]. F. Zhao, Y. Zhao, Y. Liu, X. Chang, C. Chen and Y. Zhao. Cellular uptake, intracellular trafficking, and cytotoxicity of nanomaterials.
Small 7, 13221337 (2011).
[189]. H. Kobayashi, M. Ogawa, R. Alford, P. L. Choyke and Y. Urano. New
strategies for fluorescent probe design in medical diagnostic imaging.
Chem. Rev. 110, 26202640 (2010).
[190]. J. A. Champion and S. Mitragotri. Role of target geometry in phagocytosis. Proc. Natl. Acad. Sci. U. S. A. 103, 49304934 (2006).
[191]. A. Verma and F. Stellacci. Effect of surface properties on nanoparticlecell interactions. Small 6, 1221 (2010).
[192]. Y. Qiu, Y. Liu, L. Wang, L. Xu, R. Bai, Y. Ji, X. Wu, Y. Zhao, Y. Li
and C. Chen. Surface chemistry and aspect ratio mediated cellular
uptake of Au nanorods. Biomaterials 31, 76067619 (2010).
[193]. J.-C. Boyer and F. C. J. M. van Veggel. Absolute quantum yield measurements of colloidal NaYF4 : Er3+ , Yb3+ upconverting nanoparticles. Nanoscale 2, 14171419 (2010).
[194]. G. Yi, H. Lu, S. Zhao, Y. Ge, W. Yang, D. Chen and L.-H. Guo.
Synthesis, characterization, and biological application of size-controlled
nanocrystalline NaYF4 :Yb,Er infrared-to-visible up-conversion phosphors. Nano Lett. 4, 21912196 (2004).
[195]. X. Wang, J. Zhuang, Q. Peng and Y. D. Li. A general strategy for
nanocrystal synthesis. Nature 437, 121124 (2005).
[196]. J.-C. Boyer, L. A. Cuccia and J. A. Capobianco. Synthesis of colloidal upconverting NaYF4 : Er3+ /Yb3+ and Tm3+ /Yb3+ monodisperse nanocrystals. Nano Lett. 7, 847852 (2007).
[197]. H.-X. Mai, Y.-W. Zhang, R. Si, Z.-G. Yan, L.-D. Sun, L.-P. You and
C.-H. Yan. High-quality sodium rare-earth fluoride nanocrystals: Controlled synthesis and optical properties. J. Am. Chem. Soc. 128, 6426
6436 (2006).
[198]. G. S. Yi and G. M. Chow. Synthesis of hexagonal-phase NaYF4 :Yb,Er
and NaYF4 :Yb,Tm nanocrystals with efficient up-conversion fluorescence. Adv. Funct. Mater. 16, 23242329 (2006).
[199]. Z. Li and Y. Zhang. An efficient and user-friendly method for the
synthesis of hexagonal-phase NaYF4 :Yb, Er/Tm nanocrystals with
controllable shape and upconversion fluorescence. Nanotechnology 19,
345606 (2008).
87
References
88
References
89
References
90
References
91
References
92
References
93
References
[290]. X. Huang, S. Han, W. Huang and X. Liu. Enhancing solar cell efficiency: the search for luminescent materials as spectral converters.
Chem. Soc. Rev. 42, 173201 (2013).
[291]. G. Chen, J. Seo, C. Yang and P. N. Prasad. Nanochemistry and nanomaterials for photovoltaics. Chem. Soc. Rev. 42, 83048338 (2013).
[292]. W. Qin, D. Zhang, D. Zhao, L. Wang and K. Zheng. Near-infrared photocatalysis based on YF3 : Yb3+ ,Tm3+ /TiO2 core/shell nanoparticles.
Chem. Commun. 46, 23042306 (2010).
[293]. M. Zhang, Y. Lin, T. J. Mullen, W.-F. Lin, L.-D. Sun, C.-H. Yan, T. E.
Patten, D. Wang and G.-Y. Liu. Improving Hematites solar water splitting efficiency by incorporating rare-earth upconversion nanomaterials.
J. Phys. Chem. Lett. 3, 31883192 (2012).
[294]. Y. Zhang and Z. Hong. Synthesis of lanthanide-doped NaYF4 @TiO2
core-shell composites with highly crystalline and tunable TiO2 shells
under mild conditions and their upconversion-based photocatalysis.
Nanoscale 5, 89308933 (2013).
[295]. D.-X. Xu, Z.-W. Lian, M.-L. Fu, B. Yuan, J.-W. Shi and H.-J. Cui. Advanced near-infrared-driven photocatalyst: Fabrication, characterization, and photocatalytic performance of -NaYF2 :Yb3+ ,Tm3+ @TiO2
core@shell microcrystals. Appl. Catalysis B: Environ. 142143, 377
386 (2013).
[296]. J. Park, D.-H. Lim, H.-J. Lim, T. Kwon, J.-S. Choi, S. Jeong, I.-H. Choi
and J. Cheon. Size dependent macrophage responses and toxicological
effects of Ag nanoparticles. Chem. Commun. 47, 43824384 (2011).
[297]. B. D. Chithrani, A. A. Ghazani and W. C. W. Chan. Determining the
size and shape dependence of gold nanoparticle uptake into mammalian
cells. Nano Lett. 6, 662668 (2006).
[298]. C. M. Goodman, C. D. McCusker, T. Yilmaz and V. M. Rotello. Toxicity of gold nanoparticles functionalized with cationic and anionic side
chains. Bioconjugate Chem. 15, 897900 (2004).
[299]. B. Ballou, B. C. Lagerholm, L. A. Ernst, M. P. Bruchez and A. S.
Waggoner. Noninvasive imaging of quantum dots in mice. Bioconjugate
Chem. 15, 7986 (2004).
[300]. S. Harper, C. Usenko, J. E. Hutchison, B. L. S. Maddux and R. L.
Tanguay. In vivo biodistribution and toxicity depends on nanomaterial
composition, size, surface functionalisation and route of exposure. J.
Exp. Nanosci. 3, 195206 (2008).
[301]. Y. Liu, Y. Zhao, B. Sun and C. Chen. Understanding the toxicity of
carbon nanotubes. Acc. Chem. Res. 46, 702713 (2013).
[302]. S. Liu, L. Wei, L. Hao, N. Fang, M. W. Chang, R. Xu, Y. Yang and
Y. Chen. Sharper and faster nano darts kill more bacteria: A study
of antibacterial activity of individually dispersed pristine single-walled
carbon nanotube. ACS Nano 3, 38913902 (2009).
[303]. V. E. Kagan, N. V. Konduru, W. Feng, B. L. Allen, J. Conroy,
Y. Volkov, I. I. Vlasova, N. A. Belikova, N. Yanamala, A. Kapralov,
Y. Y. Tyurina, J. Shi, E. R. Kisin, A. R. Murray, J. Franks, D. Stolz,
P. Gou, J. Klein-Seetharaman, B. Fadeel, A. Star and A. A. Shvedova.
Carbon nanotubes degraded by neutrophil myeloperoxidase induce less
pulmonary inflammation. Nat. Nano. 5, 354359 (2010).
[304]. C. D. Walkey and W. C. W. Chan. Understanding and controlling the
interaction of nanomaterials with proteins in a physiological environment. Chem. Soc. Rev. 41, 27802799 (2012).
[305]. Y. Xu, X. Lin and C. Chen. Key factors influencing the Toxicity of
nanomaterials. Chin. Sci. Bull. 58, 24662478 (2013).
[306]. J. Shan, J. Chen, J. Meng, J. Collins, W. Soboyejo, J. S. Friedberg and
94
References
95
References
96
References
97
References
98
References
T. Krucker and V. Ntziachristos. Hybrid FMT-CT imaging of amyloidbeta plaques in a murine Alzheimers disease model. Neuroimage 44,
13041311 (2009).
[370]. A. Corlu, R. Choe, T. Durduran, M. A. Rosen, M. Schweiger, S. R. Arridge, M. D. Schnall and A. G. Yodh. Three-dimensional in vivo fluorescence diffuse optical tomography of breast cancer in humans. Opt.
Express 15, 66966716 (2007).
[371]. M. J. Niedre, R. H. de Kleine, E. Aikawa, D. G. Kirsch, R. Weissleder
and V. Ntziachristos. Early photon tomography allows fluorescence
detection of lung carcinomas and disease progression in mice in vivo.
Proc. Natl. Acad. Sci. U. S. A. 105, 1912619131 (2008).
[372]. F. Leblond, H. Dehghani, D. Kepshire and B. W. Pogue. Early-photon
fluorescence tomography: spatial resolution improvements and noise
stability considerations. J. Opt. Soc. Am. A 26, 14441457 (2009).
[373]. P. I. H. Bastiaens and A. Squire. Fluorescence lifetime imaging microscopy: spatial resolution of biochemical processes in the cell. Trends
Cell Biol. 9, 4852 (1999).
[374]. L. Zhao, K. Abe, M. Barroso and X. Intes. Active wide-field illumination for high-throughput fluorescence lifetime imaging. Opt. Lett. 38,
39763979 (2013).
[375]. M. A. OLeary, D. A. Boas, X. D. Li, B. Chance and A. G. Yodh.
Fluorescence lifetime imaging in turbid media. Opt. Lett. 21, 158160
(1996).
[376]. J. McGinty, V. Y. Soloviev, K. B. Tahir, R. Laine, D. W. Stuckey,
J. V. Hajnal, A. Sardini, P. M. W. French and S. R. Arridge. Threedimensional imaging of F
orster resonance energy transfer in heterogeneous turbid media by tomographic fluorescent lifetime imaging. Opt.
Lett. 34, 27722774 (2009).
[377]. V. Gaind, S. Kularatne, P. S. Low and K. J. Webb. Deep-tissue imaging
of intramolecular fluorescence resonance energy-transfer parameters.
Opt. Lett. 35, 13141316 (2010).
[378]. K. Abe, L. Zhao, A. Periasamy, X. Intes and M. Barroso. Non-invasive
in vivo imaging of near infrared-labeled transferrin in breast cancer
cells and tumors using fluorescence lifetime FRET. PLOS ONE 8,
e80269 (2013).
[379]. J. Swartling, J. Svensson, D. Bengtsson, K. Terike and S. AnderssonEngels. Fluorescence spectra provide information on the depth of fluorescent lesions in tissue. Appl. Opt. 44, 19341941 (2005).
[380]. J. Svensson and S. Andersson-Engels. Modeling of spectral changes for
depth localization of fluorescent inclusion. Opt. Express 13, 42634274
(2005).
[381]. J. Axelsson, J. Svensson and S. Andersson-Engels. Spatially varying
regularization based on spectrally resolved fluorescence emission in fluorescence molecular tomography. Opt. Express 15, 1357413584 (2007).
[382]. A. D. Zacharopoulos, P. Svenmarker, J. Axelsson, M. Schweiger, S. R.
Arridge and S. Andersson-Engels. A matrix-free algorithm for multiple wavelength fluorescence tomography. Opt. Express 17, 30423051
(2009).
[383]. A. D. Klose. Hyperspectral excitation-resolved fluorescence tomography
of quantum dots. Opt. Lett. 34, 24772479 (2009).
[384]. A. J. Chaudhari, S. Ahn, R. Levenson, R. D. Badawi, S. R. Cherry and
R. M. Leahy. Excitation spectroscopy in multispectral optical fluorescence tomography: methodology, feasibility and computer simulation
studies. Phys. Med. Biol. 54, 46874704 (2009).
[385]. A. D. Klose and T. P
oschinger. Excitation-resolved fluorescence tomog-
99
References
100
References
[403]. J. P. Culver, V. Ntziachristos, M. J. Holboke and A. G. Yodh. Optimization of optode arrangements for diffuse optical tomography: A
singular-value analysis. Opt. Lett. 26, 701703 (2001).
[404]. H. Xu, H. Dehghani, B. W. Pogue, R. Springett, K. D. Paulsen and J. F.
Dunn. Near-infrared imaging in the small animal brain: optimization
of fiber positions. J. Biomed. Opt. 8, 102110 (2003).
[405]. T. Lasser and V. Ntziachristos. Optimization of 360 projection fluorescence molecular tomography. Med. Image Anal. 11, 389399 (2007).
[406]. F. Leblond, K. M. Tichauer and B. W. Pogue. Singular value decomposition metrics show limitations of detector design in diffuse fluorescence
tomography. Biomed. Opt. Express 1, 15141531 (2010).
[407]. E. E. Graves, J. P. Culver, J. Ripoll, R. Weissleder and V. Ntziachristos.
Singular-value analysis and optimization of experimental parameters in
fluorescence molecular tomography. J. Opt. Soc. Am. A 21, 231241
(2004).
[408]. P. C. Hansen. REGULARIZATION TOOLS: A Matlab package for
analysis and solution of discrete ill-posed problems. Numer. Algorithms
6, 135 (1994).
101
References
102
Papers
Paper I
Drug quantification in turbid media by fluorescence
imaging combined with light-absorption correction
using white Monte Carlo simulation
H. Xie, H. Liu, P. Svenmarker, J. Axelsson, C. T. Xu, S. Gr
afe,
J. H. Lundeman, H. P. H. Cheng, S. Svanberg, N. Bendsoe,
P. E. Andersen, K. Svanberg, S. Andersson-Engels.
Journal of Biomedical Optics 16(6), 066002-1 - 066002-11 (2011).
Paper I
Lund, Sweden
Abstract. Accurate quantification of photosensitizers is in many cases a critical issue in photodynamic therapy.
As a noninvasive and sensitive tool, fluorescence imaging has attracted particular interest for quantification in
pre-clinical research. However, due to the absorption of excitation and emission light by turbid media, such
as biological tissue, the detected fluorescence signal does not have a simple and unique dependence on the
fluorophore concentration for different tissues, but depends in a complex way on other parameters as well. For
this reason, little has been done on drug quantification in vivo by the fluorescence imaging technique. In this paper
we present a novel approach to compensate for the light absorption in homogeneous turbid media both for the
excitation and emission light, utilizing time-resolved fluorescence white Monte Carlo simulations combined with
the BeerLambert law. This method shows that the corrected fluorescence intensity is almost proportional to the
absolute fluorophore concentration. The results on controllable tissue phantoms and murine tissues are presented
and show good correlations between the evaluated fluorescence intensities after the light-absorption correction
and absolute fluorophore concentrations. These results suggest that the technique potentially provides the means
to quantify the fluorophore concentration from fluorescence images. C 2011 Society of Photo-Optical Instrumentation Engineers
(SPIE). [DOI: 10.1117/1.3585675]
Keywords: fluorescence; imaging; biomedical optics; tissues; optical properties; absorption; Monte Carlo; photon migration.
Paper 11047R received Feb. 2, 2011; revised manuscript received Apr. 6, 2011; accepted for publication Apr. 8, 2011; published online
Jun. 1, 2011.
1 Introduction
Photodynamic therapy (PDT) has been clinically accepted to
treat certain types of malignant tumors as well as some nonmalignant diseases.1, 2 In PDT a photosensitizer (PS) is administrated either systemically or topically. It is activated by irradiating appropriate light to the sensitized tumor. As the PS
absorbs light, the gained energy can be transferred to nearby
oxygen molecules, leading to the formation of highly reactive
oxygen radicals and thereafter tissue damage. PDT is a nonthermal photochemical reaction, which requires the presence of a
photosensitizing agent (i.e., PS), oxygen and light, simultaneously. In PDT, quantification of the PS in a noninvasive way is
in many cases a critical issue, since light dosimetry, irradiation
parameters, and therapeutic outcome depends significantly on
the PS quantities distributed in the region of interest (ROI).
Various techniques have been proposed to quantify the PS
concentration. In ex vivo animal experiments, high performance
liquid chromatography (HPLC) of excised tissues is conventionally used as a gold standard for quantitative analysis of the PS
concentration as well as for its pharmacokinetic behavior.3 Optical techniques can offer alternatives and be used in vivo. Among
Address all correspondence to: Haiyan Xie, Lund University, Department of
Physics, P.O. Box 118, 22100 Lund, Sweden. Tel: 0046-46-222 3119; Fax: 004646-222 4250; E-mail: [email protected].
these techniques, absorption spectroscopy47 provides a noninvasive tool for PS concentration studies. However, it suffers from
a relatively poor detection sensitivity, limiting its applicability.4
Most photosensitizing agents are, however, fluorescent, providing another possibility for measuring its concentration. Fluorescence has already been extensively used for tumor localization
and to assess treatment progression during diagnostic screening
or image-guided surgery to improve clinical decision-making
and the therapeutic outcome (see e.g., Refs. 811).
Fluorescence spectroscopy (either in the point-monitoring or
imaging mode) can also be used as a tool for PS concentration
measurements.4, 1215 No tissue excision would be required in
contrast to HPLC. Thus, it may constitute a tool for minimally
invasive quantification studies providing in vivo capabilities.
In these types of measurements, it is, however, a challenge to
reduce the influence of the attenuation of the probing light.
The signal depends not only on the concentration of the fluorophore but also on the optical properties, the detection geometry, and the tissue autofluorescence [i.e., the fluorescence
from endogenous tissue fluorophores such as collagen or nicotinamide adenine dinucleotide (NADH)16 ], etc. Ultimately, these
dependencies tend to decrease the correlation coefficient between the fluorescence signal and the true PS concentration.
C 2011 SPIE
1083-3668/2011/16(6)/066002/11/$25.00
066002-1
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
107
Drug quantification in turbid media by fluorescence imaging combined with light-absorption correction
using white Monte Carlo simulation
2.2 Photosensitizer
R
In the animal measurement, Fospeg
(Biolitec AG, Jena,
Germany), a liposomal formulation of the active ingredient mTHPC, or temoporfin, was employed as the PS.26, 27 Liposomes
are designed as carrier and delivery systems with the aim of
improving the tumor accumulation behavior of the PS during
PDT. The PS was diluted in 50 l of 5% glucose. All compounds
were stored at 4 C in the darkness. The extinction coefficients
and fluorescence emission spectrum of m-THPC dissolved in
ethanol were measured using a conventional spectrofluorometer
(HORIBA Jobin Yvon GmbH, Unterhaching, Germany).
066002-2
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
108
(1)
Paper I
(2)
optical fiber. The distal end of the fiber was positioned to obtain
either a transillumination or an epi-fluorescence geometry.
For the phantom experiments, Rhodamine 6G was excited by
an Nd:YAG laser (Viasho VA-I-N-532200mW, Beijing, China)
at 532 nm with a spot size of 3 mm in diameter from the bottom
of the container. The fluorescence was spectrally filtered using
a liquid crystal tunable filter (LCTF VIS 20-35, Varispec, CRI,
Inc., Woburn, California). Images of transmitted fluorescence
were acquired with the LCTF set to 600 nm, using a chargecoupled device (CCD) camera (Andor iXon DU-897, Belfast,
Ireland) with a standard camera lens (50 mm focal length and
f/1.8, Nikon, Tokyo, Japan). To suppress the transmitted laser
light from reaching the camera, a long-pass, cut-off color-glass
filter (OG-550, Schott Inc., W. Germany) was fixed between the
sample and the LCTF.
For the animal measurement, a continuous-wave laser source
at 405 nm developed at DTU Fotonik was used for excitation of the PS. It was based on a frequency doubled tapered
diode amplifier placed in an external grating cavity with Littrow feedback, described in Ref. 32. The entire laser system
was built on a breadboard and placed on a mobile cart. The
output power was 130 mW, out of which 70% was coupled
into and delivered through an optical fiber mounted above the
target. The spot size on the tissue was approximately 4 cm in
diameter. The fluorescence images were captured by the CCD
camera with the LCTF set to 652 nm (corresponding to mTHPC fluorescence peak) and 525 nm (tissue autofluorescence),
respectively.
The fluorescence images were acquired in a dimmed room to
avoid any artifacts from background light. Background images
without the excitation light were also acquired using the same
filter wavelength.
066002-3
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
109
Drug quantification in turbid media by fluorescence imaging combined with light-absorption correction
using white Monte Carlo simulation
(a)
(b)
Fig. 2 Schematic diagram of the geometry used for fluorescence white Monte Carlo simulations. (a) The medium is divided into volume elements
using small grids along the r and z axes. Similarly, time is divided into intervals with a size of dt. (b) View of the coordinate system used to calculate
2 + 2r r cos )1/2 .
the convolution of excitation and emission light for a slab of thickness dz at z, where d = (r 2 + r D
D
2.7.1
Image ratio
(3)
where Fimage (652 nm) and Fimage (525 nm) denote the mean value
of the fluorescence intensities within the ROI at the two wavelengths for each animal tissue sample.
2.7.2
R
dz
2
dr
d AW xa = 0, xs , r, z, t x
m
m
, (4)
eff (r, , z)E W m
a = 0, s , d, z, t
where A W (ax = 0, sx , r, z, t x ) is the probability per unit
volume and unit time for an excitation photon to be absorbed
at a fluorophore position (r , , z), after a time delay t x from
m
the injection point; E W (am = 0, m
s , d, z, t ) is the probability
066002-4
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
110
Paper I
Grid size
No. photons
eff
dz (m)
dr (m)
dt (ps)
2106
0.25
500
500
10
61
120
400
400
107
0.25
40
40
125
200
150
150
Phantoms
Murine tissues
(5)
m
FA (rD , t x, t m ) = FW
,
(rD , t x, t m ) exp xa vt x exp m
a vt
nz
nr
nt 1
nt 2
3 Results
3.1 Time-Resolved Fluorescence
Distribution from the FWMC Simulations
The simulated temporal distribution of the fluorescence signal
F(r D , t x , t m , z D ) is plotted in Fig. 3. This corresponds to the
detected signal at a surface of the simulated phantom at a radial
distance r D = 5.0 mm, where z D indicates if the detector is at
the bottom (the left column, z D = 3.1 cm, corresponding to the
transillumination geometry) or the top surface (the right column,
(6)
where v = c/n is the light speed in the medium with a refractive
index n. Integrating over the time, we will get a signal corresponding to the recorded fluorescence signal from a pixel at a
radius r D , or a Cartesian position (x, y) in the image, where
r D = (x 2 + y 2 )1/2 . We denote this Fimage (x, y):
FA (rD , t x , t m )dt x dt m .
(7)
Fimage (x, y) = Fimage (rD ) =
tx tm
(9)
Fimage
FW
(rD , t x , t m )dt x dt m ds eff Conc,
t x t m ROI
(10)
Journal of Biomedical Optics
066002-5
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
111
Drug quantification in turbid media by fluorescence imaging combined with light-absorption correction
using white Monte Carlo simulation
m
a [cm1]
25
20
15
10
5
ax [cm1]
11
(b) 10 x 10
Fluorescence intensities [a.u.]
11
10
10
10
10
a=0.012 cm
10
a=0.14 cm
10
2
1
0
0.025
0.05
0.075
0.1
Rhodamine 6G concentration [M]
0.5
400
500
600
800
0
900
the same ink concentration were the same at both 532 and
600 nm, as the India ink absorption is very flat over this wavelength range. For different ink concentrations, the results were
0.37, 0.26, 0.14, and 0.012 cm1 , respectively. The reduced
scattering coefficients were measured and extrapolated to be
10.10 cm1 at 532 nm and 8.30 cm1 at 600 nm.
The simulated dependence of the correction factor on different absorption coefficients at the excitation and emission
wavelength is shown in Fig. 4(a). For the four sets of phantoms
with different absorptions, the correction factor was = 5303,
949, 87, and 2, respectively. The large variation in illustrates
that the fluorescence signal from the surface is heavily dependent
on the optical properties of the absorbing medium. Figure 4(b)
shows the correlation between the fluorescence signal and the
fluorescent dye (Rhodamine 6G) concentrations before and after
the light-absorption correction, respectively. The markers without circles represent the raw signals directly from the images,
while the markers with circles representing the corrected fluorescence intensities, which were achieved by multiplying the
measured fluorescence from the images by the corresponding
simulated . The correlation between the fluorescence signal
and the fluorophore concentration is dramatically improved from
3000:1 to 1.3:1.
066002-6
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
112
700
Wavelength [nm]
=0.26 cm1
a=0.37 cm
0
300
12
10
log() [a.u.]
(a)
Paper I
Fig. 6 Fluorescence images of some animal organs ex vivo captured at different emission wavelengths with an exposure time of 11 s. The color bars
indicate fluorescence intensity after background subtraction. (a) The top row shows the images taken with the LCTF set to 652 nm (corresponding
to the maximal drug fluorescence), while (b) the bottom row shows the images at 525 nm (tissue autofluorescence).
(a)
(b)
5000
(c)
2000
1000
0.4
0.8
1.2
8
Corrected Fluo @652nm [a.u.]
3000
Fluo@652nm/Fluo@525nm [a.u.]
Muscle
Liver
Kidney
Lung
4000
x 10
400
300
200
100
2
0
0
0
1.6
0.5
HPLC [ng/mg]
1.5
0.4
HPLC [ng/mg]
0.8
1.2
1.6
HPLC [ng/mg]
Fig. 7 Scatter plots showing (a) the fluorescence signals captured at 652 nm versus the HPLC values for individual organs, where a linear fit of the
data points for each type of organs (solid line) is also shown; (b) the spectral ratio of the fluorescence intensities at 652 nm to that at 525 nm; and
(c) the fluorescence signals after the light-absorption correction versus the HPLC values.
6000
5000
250
x 10
(c) Fluorescence*
200
4000
Slope k [a.u.]
4
150
3000
3
100
2000
2
50
1000
Muscle
Liver
Kidney
Lung
Muscle
Liver
Kidney
Lung
Muscle
Liver
Kidney
Lung
Fig. 8 Comparison of the linear fit slopes between the fluorescence signals and the HPLC values. The methods are the same as stated in Fig. 7.
066002-7
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
113
Drug quantification in turbid media by fluorescence imaging combined with light-absorption correction
using white Monte Carlo simulation
6.2
0.4
88.7
23.0
0.8
26.5
1.8
57.5
34.9
0.8
4.9
0.3
240.9
117.4
0.8
13.5
0.6
142.0
110.0
0.8
the fluorescence imaging data and the HPLC values, mainly due
to the varying optical properties of the tissue under investigation.
As shown in the middle columns, the slope from the spectral
ratio varies from 240 for muscle to 60 for liver (4.4:1). After the
absorbed light is compensated for, using the optical properties
of a different type of murine organs listed in Table 2, the result
on the right shows the slope varying from 53,500 for liver to
21,400 for lung (2.5:1).
4 Discussion
The simulated time-resolved fluorescence signal provides us
with a clear picture of how the excitation and emitted photons
migrate in the medium. With the information of the photomigration time, one can simulate the amount of light which is absorbed
within the turbid medium, using the proposed absorption correction method. As demonstrated with tissue-like phantoms in
Fig. 4(b), the correlation between the corrected (intrinsic) fluorescence signal and the dye concentration has been dramatically
improved by approximately 2000-fold in comparison to the uncorrected signal, when absorption of the phantoms was varied
over a wide range (approximately by a factor of 30). However,
the corrected correlation is not a perfect single line (about 30%
error in the correlation) for different sets of optical properties.
This results mainly from about a 10% relative error in the optical properties measured by the TOF spectroscopy system.34
Therefore, we tested the sensitivity of due to the error in the
measured optical properties by altering them to a 10% variation, while the other parameters remained the same in the simulations. How is affected for the set of phantoms with the highest
absorption and scattering coefficients (ax = am = 0.37 cm1 ,
1
m
1
x
s = 10.10 cm , s = 8.30 cm ) is shown in Tables 3
and 4, respectively. The relative change in is within the range
of 18% to +36%. For the other phantoms, a smaller variation
in was obtained for lower values of the optical properties.
Provided more accurate optical properties, one could expect
a better correlation between the fluorescence and fluorophore
concentrations. In addition, photobleaching of Rhodamine 6G
(which also could provide an error) is not relevant due to the
low light fluence and short acquisition time used.
R
For the animal study, Fospeg
content ex vivo in murine
organs was quantified. This drug contains the active ingredient
m-THPC, and therefore is interesting as a PDT sensitizer. The
observed signal level directly from the fluorescence imaging
measurement is affected by many factors. It generally leads to
a poor correlation between the signal level and drug concentration, as shown in Figs. 7(a) and 8(a). The fitted linear slope of the
raw fluorescence signals at 652 nm varies drastically between
different organs. The large variations in the raw signals indicate
that absorption caused by tissue should be taken into account to
allow quantification of the PS in the tissue samples. Tissue absorption at the interesting wavelengths is mostly dependent on
the amount of hemoglobin (both oxy- and deoxy-hemoglobin)
in the tissue. The high absorbance in lung, kidney, and liver is
caused by high blood contents, attenuating both the excitation
and emission light. This effect explains the lower slope of the
correlation curves for these organs as compared to those for
muscle and skin. The detailed selectivity and biodistribution of
R
Fospeg
following systemic administration were studied separately in Ref. 31.
To compensate for the tissue absorbed light, the optical properties for different murine organs need to be known. Values from
literatures vary a lot and seem not fully consistent between tissue types. The large variation in literature values results from
the fact that measurements were performed on tissues from different species, in vitro or in vivo, and at a great variety of sample
preparation techniques and other experimental conditions. In the
preparation period of the tissue samples in this study, bleeding
during organ removal caused partial loss of the blood resulting
in a variation of absorption coefficient even for the same type of
tissues. As mentioned above, we do not have the precise values
of the tissue optical properties. Instead, the tabulated values from
Ref. 23 were used. They could only be interpreted as reasonable
but rough estimates of the true mouse tissue optical properties.
This is one of the main reasons why the correlation is not as
good as that for the phantom. The sensitivity of the correction
factor to tissue absorption coefficients are also examined for the
liver with highest absorption (ax = 26.5 cm1 , am = 1.8 cm1 ,
1
m
1
x
s = 11.5 cm , s = 7.0 cm ). In Table 5, the absorption
coefficients have been altered by 30%, resulting in a variation
in from 34% to +41%. Apparently this variation is less
sensitive than that of the phantoms due to the much shorter migration pathlengths the photons spend in the epi-fluorescence
geometry. Furthermore, if the tissue optical properties could be
measured individually for each tissue investigated, simultaneously with the fluorescence imaging measurements, the correlation could be further improved.
The absorption correction method proposed in this study
is based on assumptions that the fluorescent drug is homogeneously distributed within these tissues, and that all tissues
have homogeneous optical properties. Unfortunately very little
is known about distribution of this new drug following system-
Table 3 Sensitivity of the simulated correction factor to absorption coefficients of the phantoms.
Variation in
(ax , am )
(ax 10%, am )
(ax + 10%, am )
(ax , am 10%)
(ax , am + 10%)
23%
+ 26%
27%
+ 32%
066002-8
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
114
Paper I
Variation in
(sx , sm )
22%
+ 25%
18%
+ 19%
Table 5 Sensitivity of the simulated correction factor to absorption coefficients of murine livers.
Variation in
(ax , am )
(ax 30%, am )
(ax + 30%, am )
(ax , am 30%)
(ax , am + 30%)
34%
+ 41%
10%
+ 9%
066002-9
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
115
Drug quantification in turbid media by fluorescence imaging combined with light-absorption correction
using white Monte Carlo simulation
5 Conclusion
In conclusion, this paper has shown that the fluorescence imaging technique can be used as a noninvasive and sensitive tool to
quantify the fluorescent markers in homogeneous turbid media,
using the novel light-absorption correction approach combining
the fluorescence imaging and FWMC simulations. The results
on both the well-controlled tissue-like phantoms and ex vivo animal tissues have shown that this method provides an acceptable
quantification of fluorescent molecule markers in media with
known geometry and optical properties at both the excitation and
emission wavelengths. An improved linear correlation with the
true concentrations is obtained independent of the tissue optical
properties, since this method efficiently compensates for light
attenuation and thus more directly relates to the intrinsic fluorescence signal levels from the fluorophores. This approach offers
the advantages of minimizing the dependence on the tissue optical properties, a very low concentration detection limit, and wide
spectral bands. The sensitivity of the results to the medium optical properties variation are presented and discussed, in order to
point toward the possible future improvement of this technique.
In future work, it is highly desirable to combine the sensitivity of the fluorescence imaging technique with the tissue
optical properties measurements to constitute an even better and
more reliable fluorophore concentration estimate. The absorption correction and image ratio methods could also be combined
to compensate for tissue autofluorescence.
Acknowledgments
The authors are grateful to Biolitec AG for the realization of
animal studies, HPLC analysis of tissue samples, and providing photosensitizer formulations (a drug they have a commercial interest in), and to Erik Alerstam for his assistance with
the time-of-flight measurements. This work was funded by the
EUBrighter project (FP6-IST-035266) (all groups), and a
Linnaeus grant for the Lund Laser Center.
1. A. Johansson and S. Andersson-Engels, Chapter 8 in Laser Imaging
and Manipulation in Cell Biology, pp. 167199, Wiley-VCH, Verlag
(2010).
2. K. Svanberg, N. Bendsoe, J. Axelsson, S. Andersson-Engels, and
S. Svanberg, Photodynamic therapy: superficial and interstitial illumination, J. Biomed. Opt. 15(4), 041502 (2010).
3. H.-P. Lassalle, D. Dumas, S. Grafe, M.-A. DHallewin, F. Guillemin,
and L. Bezdetnaya, Correlation between in vivo pharmacokinetics,
intratumoral distribution and photodynamic efficiency of liposomal
mTHPC, J. Controlled Release 134(2), 118124 (2009).
4. A. Johansson, J. Svensson, N. Bendsoe, K. Svanberg, E. Alexandratou,
M. Kyriazi, D. Yova, S. Grafe, T. Trebst, and S. Andersson-Engels,
Fluorescence and absorption assessment of a lipid mTHPC formulation
following topical application in a non-melanotic skin tumor model,
J. Biomed. Opt. 12(3), 034026 (2007).
5. J. R. Mourant, I. J. Bigio, D. A. Jack, T. M. Johnson, and H. D. Miller,
Measuring absorption coefficients in small volumes of highly scattering media: source-detector separations for which path lengths do
not depend on scattering properties, Appl. Opt. 36(22), 56555661
(1997).
6. J. R. Mourant, T. M. Johnson, G. Los, and I. J. Bigio, Non-invasive
measurement of chemotherapy drug concentrations in tissue: prelimi-
066002-10
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
116
Paper I
066002-11
Downloaded from SPIE Digital Library on 30 Nov 2011 to 130.235.39.169. Terms of Use: http://spiedl.org/terms
117
Paper II
Synthesis of NaYF4 :Yb3+ , Er3+ upconverting
nanocrystals in a capillary-based continuous
microfluidic reaction system
H. Liu, O. Jakobsson, C. T. Xu, H. Xie, T. Laurell,
S. Andersson-Engels.
Proc. of SPIE 7909, 790917-1 - 790917-6 (2011).
Paper II
2 Department
ABSTRACT
We report for the rst time continuous ow synthesis of NaYF4 : Yb3+ , Er3+ upconverting nanocrystals in a
capillary-based microuidic reaction system. Two sequential temperature steps were employed for heating with
an initial high temperature (180 C) step to burst nuclei and a subsequent low temperature (110 C) step to
promote growth of the nanocrystals in order to obtain high-performance nanocrystals. The prepared nanocrystals
exhibit green and red emissions under excitation of 974 nm diode laser. Our research opens the door for the
synthesis of upconverting nanocrystals in microuidic systems.
Keywords: upconverting nanocrystals, continuous ow, microuidic
1. INTRODUCTION
Upconverting nanocrystals doped with lanthanide ions have drawn much attention due to their potential as optical imaging probes in biomedical applications. These contrast agents have several advantages over conventional
uorescent biomarkers (e.g., organic dyes and semiconductor nanocrystals), such as dramatically reduced autouorescence from cells or tissues, better light penetration depth, no photo-damage to living organisms, and higher
spatial resolution in the captured images.13 Upconverting nanocrystals have come to be widely used in such
diverse elds as, luminescence microscopy,4 reectance and transillumination imaging of biological tissues,57
and uorescence diuse optical tomography8, 9during the last decade.
Among upconverting nanomaterials, lanthanide-ions-doped NaYF4 nanocrystals (NaYF4 : Yb3+ , Er3+ /Tm3+ )
have been demonstrated to be the most ecient to date. Many eorts are being made in order to develop ecient
and less complex methods for synthesizing high quality and monodispersed NaYF4 : Yb3+ , Er3+ /Tm3+ nanocrystals for various biomedical applications. Currently, NaYF4 nanocrystals are mainly synthesized in batch-control
modes in small volumes, for example, using the hydrothermal method,10 solvothermal method,11 microwave
synthesis,12 and pyrolysis of triuoroacetates.13 Although achievements have been made, batch reactors often
suer from the disability to rapidly establish homogenous reaction conditions, e.g., during temperature changes.
This problem is commonly addressed by minimizing reactor volumes, thus leading to a negative impact on the
production rate. It is also dicult to implement fast screening and optimize the synthesis conditions in batch
modes.14
Microuidic reaction systems oer a solution to these challenges and have an increasingly important role in
the synthesis of nanocrystals as highly controlled thermal and stoichiometric microenvironments can be obtained
in the synthesis process.15 Particularly the synthesis of quantum dots was early targeted, and the strategy of
using microuidic reactors was rst presented by Edel J.B. et al, at Imperial College London in 2002.16 In
this paper, we report continuous ow based synthesis of NaYF4 : Yb3+ , Er3+ nanocrystals in a capillary-based
microuidic reaction system. In order to obtain high-performance nanocrystals, a typical synthesis scheme for
colloidal nanocrystals was utilized, which separated the nucleation and growth phase of the nanocrystals by
employing two temperature steps.17, 18
Further author information: (Send correspondence to Haichun Liu)
Haichun Liu: E-mail: [email protected], Telephone: +46 46 222 7471
Colloidal Quantum Dots/Nanocrystals for Biomedical Applications VI, edited by Wolfgang J. Parak, Kenji Yamamoto,
Marek Osinski, Proc. of SPIE Vol. 7909, 790917 2011 SPIE CCC code: 1605-7422/11/$18 doi: 10.1117/12.874348
121
2. EXPERIMENTAL
2.1 Chemicals
All chemicals were purchased from Sigma-Aldrich, Ltd., and used as received without further purication. Ln2 O3
(Ln = Y, Yb, Er) (99.99%), NaF (99.9%), and HNO3 (70%) were used as precursor materials, while polyethylene
glycol (PEG, 99%) with an average molecular weight of 200 g/mol and ethylene glycol (EG, 99%) were used as
the solvents. Oleic acid (OA, 90%) was used as the surfactant.
C oil bath, and subsequently heated in a 110 C oil bath. Retention times of the reaction mixture in the two
oil baths were determined by the ow rates and capillary lengths. The ow rates of Ln3+ and NaF solutions
were 50 l/min and 200 l/min, respectively, while the capillary lengths in the rst and second oil bath were 8
cm and 100 cm, respectively, resulting in retention times of 10 s in the rst oil bath and 120 s in the second oil
bath. The product were collected from the outlet of the capillary in a test tube. A schematic representation of
the reaction system is shown in Fig. 1. In the synthesis of sample B, the PEG solution of lanthanide ions with
the same concentration and molar ratio of lanthanide ions was used to replace corresponding EG solution, and
OA was further added in with a volume ratio of OA/PEG 1:1. The nanocrystals were separated by diluting the
obtained suspension using acetone followed by centrifugation, and washed using ethanol and deionized water for
several times. The nanocrystals were nally redispersed in ethanol for further characterization.
2.3 Characterization
The size and morphology of the nanocrystals were observed using a JEOL 3000SFF transmission electron microscope (TEM). Samples were prepared by placing a drop of ethanol dispersed nanocrystals on the surface of a
copper specimen grid coated with a perforated carbon lm. The uorescence spectra were recorded on a USB
6500 Ocean Optics spectrometer equipped with a 974 nm diode laser as the excitation source with nanocrystals
redispersed in ethanol.
122
Paper II
50 nm
50 nm
(a)
(b)
or ellipsoidal precursors synthesized in the microcapillary, as reported in the synthesis of worm and chain-like
metal nanoparticles in Ref.19. One possible reason is that the growth of nanocrystals were not properly quenched
due to the lack of protection of surfactants. In view of this, OA was added in the subsequent experiment, the
product of which was shown in Fig. 2 (b). As can be seen, the synthesized nanocrystals (sample B) show more
regular shapes and have much smaller size (less than 10 nm) as compared to sample A. However, they tend to
aggregate when dispersed in ethanol, probably due to the adsorption of OA on the surface of the nanocrystals.
These results reveals that OA can be employed to modify the shape and size of upconverting nanocrystals in
the microuidic reaction. The inuence of OA on the microuidic synthesis of upconverting nanocrystals is
currently being thoroughly investigated, and the postprocessing procedures are being improved in order to wash
and redisperse the nanocrystals.
The upconversion uorescence spectra of sample A and B under excitation of 974 nm are presented in Fig. 3
(a) and (b), respectively. As seen, both sample A and B have dominant red emissions, while sample A has a
much larger red/green emission ratio (10:1), which can be explained by the dierence in the size and morphology
of the nanocrystals obtained in the two dierent batches. The emission bands can easily be assigned to intra-4f
electronic transitions of the Er3+ ions. The green emissions between 510 and 530 nm and between 530 and 570
nm were assigned to the 2 H11/2 4 I15/2 and 4 S3/2 4 I15/2 transitions, respectively. The red emission between
635 and 700 nm originated from the 4 F9/2 4 I15/2 transition.20
In order to determine the number of photons responsible for the upconversion mechanism, the dependence of
the upconversion emissions of sample A on the excitation power was investigated. As presented in Fig. 4, both the
green and red upconversion emission intensities demonstrated second order power dependence at low excitation
densities, indicating a two-photon upconversion mechanism. The power dependencies of Er3+ upconversion
emissions became linear at high excitation densities due to saturationof the upconversion processes.21 The
upconversion emissions of sample B have similar power dependences, not shown here.
The upconversion excitation pathways in the Yb3+ /Er3+ system are well-known and are shown in Fig. 5.
The Er3+ ions are rstly excited to 4 I11/2 via an energy transfer from Yb3+ ions in 2 F5/2 state. A second 974
nm photon, or energy transfer from Yb3+ ions can then pump Er3+ ions to the 4 F7/2 level. Subsequently the
Er3+ ions undergo a multi-phonon assisted relaxation to the 2 H11/2 and 4 S3/2 levels and the green emissions
occur via transitions of 2 H11/2 /4 S3/2 4 I15/2 . Alternatively, the ions can further relax and populate the 4 F9/2
leading to the red emission via transition of 4 F9/2 4 I15/2 .20
123
6000
4500
5500
4000
5000
3500
4500
3000
15/2
F I
9/2
4
15/2
15/2
3/2
500
S I
1000
1500
500
2000
11/2
1000
4I
1500
S3/2 I15/2
15/2
2000
2500
2500
11/2
F9/24I15/2
3000
Intensity [a.u.]
3500
Intensity
[a.u.]
4000
0
500
550
600
650
Wavelength [nm]
700
750
500
550
600
650
Wavelength [nm]
(a)
700
750
(b)
Fig. 3. Measured upconvertsion emission in NaYF4 : Yb3+ , Er3+ nanocrystals under excitation of 974 nm diode laser.
13
4
F9/24I15/2
12
S3/24I15/2
11
10
9
8
7
6
5
5.5
6.5
7
ln(Laser Power (mW))
7.5
Fig. 4. Power dependence of the upconversion emissions of sample A when excited at 974 nm.
25
4
F7/2
2
H11/2
4
S3/2
20
15
2
10
F9/2
I9/2
I11/2
F5/2
I13/2
I15/2
5
2
Yb3+
F7/2
Er3+
Fig. 5. The energy level diagrams of the Er3+ and Yb3+ dopant ions and the upconversion mechanisms for the green
and red emissions under excitation of 974 nm diode laser. The solid, dotted, and curly arrows represent emission, energy
transfer, and multiphonon relaxation processes, respectively.
124
Paper II
4. CONCLUSIONS
The NaYF4 : Yb3+ , Er3+ nanocrystals were successfully synthesized in a capillary-based continuous microuidic
reaction system. The nanocrystals prepared without the presence of OA shows rod-like shapes with a large
length distribution but a relative narrow diameter distribution. The shape and size of the nanocrystals can
be dramatically modied by addition of surfactant OA. The prepared nanocrystals can emit green and red
upconversion emissions under excitation of 974 nm, and both the green and red emissions originate from twophoton processes. On-going research targets the further improvements of the upconverting nanocrystals synthesis
process.
ACKNOWLEDGMENTS
This work was supported by a Swedish Research Council grant (VR 2007-4214) and a Linnaeus grant for the
Lund Laser Centre.
REFERENCES
[1] Auzel, F., Upconversion and anti-stokes processes with f and d ions in solids, Chem. Rev. 104(1), 139173
(2004).
[2] Xu, C. T., Svensson, N., Axelsson, J., Svenmarker, P., Somesfalean, G., Chen, G., Liang, H., Liu, H.,
Zhang, Z., and Andersson-Engels, S., Autouorescence insensitive imaging using upconverting nanocrystals
in scattering media, Appl. Phys. Lett. 93(17), 171103 (2008).
[3] Svenmarker, P., Xu, C. T., and Andersson-Engels, S., Use of nonlinear upconverting nanoparticles provides
increased spatial resolution in uorescence diuse imaging, Opt. Lett. 35(16), 27892791 (2010).
[4] Yu, M., Li, F., and Chen, Z., Laser scanning up-conversion luminescence microscopy for imaging cells
labeled with rare-earth nanophosphors, Anal. Chem. 81(3), 930935 (2009).
[5] Nyk, M., Kumar, R., Ohulchanskyy, T. Y., Bergey, E. J., and Prasad, P. N., High contrast in vitro and in
vivo photoluminescence bioimaging using near infrared to near infrared up-conversion in Tm3+ and Yb3+
doped uoride nanophosphors, Nano Lett. 8(11), 38343838 (2008).
[6] Chatterjee, D., Rufaihah, A., and Zhang, Y., Upconversion uorescence imaging of cells and small animals
using lanthanide doped nanocrystals, Biomaterials 29(7), 937943 (2008).
[7] Vinegoni, C., Razansky, D., Hilderbrand, S. A., Shao, F., Ntziachristos, V., and Weissleder, R., Transillumination uorescence imaging in mice using biocompatible upconverting nanoparticles, Opt. Lett. 34(17),
25662568 (2009).
[8] Xu, C. T., Axelsson, J., and Andersson-Engels, S., Fluorescence diuse optical tomography using upconverting nanoparticles, Appl. Phys. Lett. 94(25), 251107 (2009).
[9] Liu, H., Xu, C. T., and Andersson-Engels, S., Multibeam uorescence diuse optical tomography using
upconverting nanoparticles, Opt. Lett. 35(5), 718720 (2010).
[10] Sun, Y., Chen, Y., Tian, L., Yu, Y., Kong, X., Zhao, J., and Zhang, H., Controlled synthesis and morphology dependent upconversion luminescence of NaYF4 : Yb, Er nanocrystals, Nanotechnology 18(27), 275609
(2007).
[11] Wang, L. and Li, Y., Controlled synthesis and luminescence of lanthanide doped NaYF4 nanocrystals,
Chem. Mater. 19(4), 727734 (2007).
[12] Wang, H.-Q. and Nann, T., Monodisperse upconverting nanocrystals by microwave-assisted synthesis,
ACS Nano 3(11), 38043808 (2009).
[13] Boyer, J.-C., Cuccia, L. A., and Capobianco, J. A., Synthesis of colloidal upconverting NaYF4 : Er3+ /Yb3+
and Tm3+ /Yb3+ monodisperse nanocrystals, Nano Lett. 7(3), 847852 (2007).
[14] Song, Y., Hormes, J., and Kumar, C. S. S. R., Microuidic synthesis of nanomaterials, Small 4(6), 698711
(2008).
[15] Marre, S. and Jensen, K. F., Synthesis of micro and nanostructures in microuidic systems, Chem. Soc.
Rev. 39(3), 11831202 (2010).
[16] Edel, J. B., Fortt, R., deMello, J. C., and deMello, A. J., Microuidic routes to the controlled production
of nanoparticles, Chem. Commun. (10), 11361137 (2002).
125
[17] Yin, Y. and Alivisatos, A. P., Colloidal nanocrystal synthesis and the organic-inorganic interface, Nature 437(7059), 664670 (2005).
[18] Zhu, X., Zhang, Q., Li, Y., and Wang, H., Redispersible and water-soluble LaF3 : Ce, Tb nanocrystals via
a microuidic reactor with temperature steps, J. Mater. Chem. 18(42), 50605062 (2008).
[19] Song, Y., Sun, Q., Zhang, T., Jin, P., and Han, L., Synthesis of worm and chain-like nanoparticles by a
microuidic reactor process, J. Nanopart. Res. 12(7), 26892697 (2010).
[20] Suyver, J. F., Grimm, J., van Veen, M. K., Biner, D., Kramer, K. W., and G
udel, H. U., Upconversion
spectroscopy and properties of NaYF4 doped with Er3+ , Tm3+ and/or Yb3+ , J. Lumin. 117(1), 112
(2006).
[21] Suyver, J. F., Aebischer, A., Garca-Revilla, S., Gerner, P., and G
udel, H. U., Anomalous power dependence
of sensitized upconversion luminescence, Phys. Rev. B 71(12), 125123 (2005).
126
Paper III
Upconverting nanoparticles for pre-clinical diffuse
optical imaging, microscopy and sensing: Current
trends and future challenges
C. T. Xu, Q. Zhan, H. Liu, G. Somesfalean, J. Qian, S. He,
S. Andersson-Engels.
Laser & Photonics Reviews 7(5), 663-697 (2013).
Paper III
REVIEW
ARTICLE
1. Introduction
Upconverting nanoparticles (UCNPs) constitute a novel
type of contrast agent with highly interesting and unique
properties for luminescence bioimaging. The aim of this review paper is to illustrate the great potential of this emerging field. The phenomenon of upconverting luminescence
has been studied for decades, with the earliest experimental
works presented in, e.g., Refs. [1, 2]. The idea originated
from Bloembergen in 1959 [3], who proposed that infrared
(IR) light could be detected by sequential stepwise absorption of an ion in a solid material. An extensive review of
the early work was given by Auzel in 2004 [4]. Most efforts were focused on rare-earth (RE) ions doped into a
lattice host. Such ions are ideal because they exhibit very
long lifetimes for the intermediate states in the upconver-
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
129
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
130
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
665
Figure 1 (online color at: www.lpr-journal.org) Energy level structure and proposed UC mechanisms of the Yb3+ , Er3+ / Tm3+ /Ho3+
co-doped UCNPs under excitation at 975 nm.
2.1.2. Activators
www.lpr-journal.org
located within the window of optical transparencyfor biological tissues, in which both light absorption and scattering
are significantly reduced. This feature makes Tm3+ -doped
UCNPs particularly interesting for imaging of deeply located bio-targets. The Tm3+ ions have another two generally less efficient UC emission bands at around 479 and
648 nm, generated by the transitions 1 G4 3 H6 and
1
G4 3 F4 , respectively [82]. The blue emission band
is less suitable for bioimaging due to much higher light
absorption and scattering in biological tissues at this wavelength. The intensity ratio between different UC bands is
host and doping concentration dependent, even for the same
activator.
For other RE ions, Wang et al. [83] has made pioneering
work to exploit their UC emission ability. They suggested
a general approach for realizing efficient UC emissions
through gadolinium sublattice-mediated energy migration,
by incorporating a set of RE ions into separated layers
at precisely defined concentrations. In this way, they have
demonstrated efficient UC emissions from a wide range of
activators, such as Tb3+ , Eu3+ , Dy3+ and Sm3+ . Generally, the doping concentration of activators is relatively low
(usually < 2 mol%) in order to minimize cross-relaxation
energy losses.
2.1.3. Sensitizers
In single-doped nanoparticles, the efficiencies of UC emissions are relatively low. This is due to the difficulty in
finding an equilibrium point for minimizing the quenching
effect by reducing the RE-ion concentration, and maximizing the absorption of pump energy by increasing its concentration. To enhance the UC luminescence efficiency, a
sensitizer with a sufficiently large absorption cross-section
in the NIR region is usually incorporated together with the
activator. Certainly, it is required that efficient energy transfer between the sensitizer and activator can occur. These
types of sensitizers can be called direct sensitizers. For
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
131
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
In a typical thermal decomposition procedure, metal trifluoroacetates are thermally decomposed to corresponding
metal fluorides. Zhang et al. [55] first reported the synthesis of single-crystalline and monodisperse LaF3 triangular
nanoplates via the thermal decomposition of lanthanum trifluoroacetates (La(CF3 COO)3 ) in a mixture of oleic acid
(OA) and octadecene (ODE). The approach was later developed as a general route to synthesize high-quality REF3
and NaREF4 nanoparticles [61, 62, 67, 80, 98] (Fig. 2(a-b)).
In this approach, ODE with a high boiling temperature (315 C) is used to provide a high-temperature environment, while OA having a good coordinate capability,
is acting as the coordinating solvent which aids in capping the surface of UCNPs to prevent agglomeration. Be-
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
132
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
667
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
133
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
Figure 3 (online color at: www.lpr-journal.org) Photoluminescence spectra emitted by UCNPs doped with different RE-ions.
The NaYbF4 :Yb3+ /Er3+ /Ho3+ /Tm3+ nanoparticles are suspended
in a chloroform solution and being irradiated by a 915-nm diode
laser with a power density of 500 mW/cm2 . (Reprinted with permission from Ref. [139]. Copyright 2011, American Chemical
Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
134
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
669
N j Nl C ET
ji,lk
i j,kl
NR
Ni Wi,i1
Ni AiED
j
n
I Iex
.
MD
+ Ai j
Ni Nk CiET
j,kl
(1)
i j,kl
MD
where AiED
j and Ai j are Einstein coefficients for electric
dipole (ED) and magnetic dipole (MD) radiative transitions
NR
from manifold i to j; Wi,i1
is the nonradiative multiphonon
relaxation (MPR) rate constant from manifold i to i 1;
ET
Ci j,kl is the microscopic energy transfer parameter for the
transfer of energy via the donor transition i j and the
acceptor kl transition. In this model, the interactions
among more than two ions (such as the CS process) are not
considered. The intensity of any given UC emission peak
is proportional to the product of the population density of
the emitting state and the microscopic rate constants for the
radiative transition. Obviously, these rate constants play
critical roles in depicting the UC emissions. Once they are
known, the characteristics of UC emissions, including the
QYs and spectral purities, can be obtained by calculation
directly, as shown in a recent work of Chan et al. [142]. This
model can also be used to determine critical energy transfer
transitions involved in UC process [142]. However, the use
of this model is demanding mainly due to the difficulty
in determining the rate constants, although in theory the
ED, MD transition rates and energy transfer rate CiET
j,kl can
be calculated using the Judd-Ofelt theory [39, 40], while
the quantum mechanical magnetic dipole operator [144]
and the nonradiative MPR can be treated with a modified
energy gap law [145]. Thus, this model is usually simplified
to only include the major transitions and MPR processes
identified by previous studies.
The investigation on the time dependent behavior of
UC emissions is a good way to verify the validity of the
proposed UC pathways, and the rate constants could be
extracted by fitting the measured time-dependent emission
intensity with time-dependent rate equations [146148].
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
135
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
Figure 5 (online color at: www.lpr-journal.org) Quantum yield of the 800 nm emission band in core-shell
NaYF4 :Yb3+ /Tm3+ @NaYF4 nanoparticles and core nanoparticles NaYF4 :Yb3+ /Tm3+ . The quantum yield increases linearly
with the excitation intensity until a saturation point, from which
the quantum yield approaches a constant value. Solid lines show
simulated corresponding QYs for highly efficient two-photon
dyes under identical experimental conditions. (Reprinted with
permission from Ref. [151]. Copyright 2012, American Chemical
Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
136
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
671
A shield layer on the surface of the UCNPs can reduce crystal defects and protect the active optical ions
from the coupling with vibrational modes in the solvent,
thereby reducing the non-radiative energy losses [42].
The shell could be either inert (i.e., with no optical active ions) or active (i.e., containing active ions, mainly
Yb3+ ). For example, increased QYs were observed after growing of an undoped NaYF4 shell over Er3+ - or
Tm3+ -doped core UCNPs [50, 151]. In another report,
a 50-fold higher QY was observed for sub-10 nm UCNPs, and the emission from 9-nm core-shell UCNPs was
larger than that of comparable 37-nm cores when normalized to the absorbance at 980 nm [126]. For the use of
active shells, examples include comparison of the pure
BaGdF5 :Yb3+ /Er3+ core with the active-shell coated counterparts, where it was shown that the luminescence intensity could be enhanced by several hundreds of times [158].
In another report, a significant increase in the UC intensity was measured in NaGdF4 :Er3+ , Yb3+ @NaGdF4 :Yb3+
mol 20% active-core/active-shell UCNPs compared to
NaGdF4 :Er3+ ,Yb3+ @NaGdF4 active-core/inert-shell and
NaGdF4 :Er3+ /Yb3+ core-only nanoparticles [159]. In addi-
www.lpr-journal.org
tion, by spatially separating dopant ions using a layer (laminated) structure, the effect of the concentration quenching
can be suppressed, and an enhancement of the UC emission
can be expected [83, 159, 160].
3.4.3. Surface plasmonics coupling
It is well known that the unique surface plasmon properties of metallic structures can be exploited to enhance
the fluorescence from adjacent fluorophores (organic dyes
and inorganic QDs) [161,162]. Similarly, surface plasmons
with strong local field can also be used to enhance the
efficiency of UC emission [163,164]. Zhang et al. and Sudheendra et al. have successfully attached gold nanoparticles onto UCNPs to modulate the UC emission [165, 166].
A specifically designed plasmonic gold surface coupling
with 980-nm radiation was shown to clearly enhance NIRto-visible upconversion luminescence from the nanocrystalline layer [167]. It is very important to study the surface
plasmon enhancement mechannism at single nanoparticle
level in order to provide an optimal design of hybrid UCNPs
and metallic nanoparticles. Schietinger et al. coupled single NaYF4 :Yb3+ /Er3+ nanocrystals with gold nanospheres
(30 and 60 nm in diameter) to enhance UC emission in a
combined optical and atomic force microscopy (AFM) system, gaining an overall enhancement factor of 3.8 [168].
Time-resolved measurements revealed that both the excitation and the emission processes are influenced by the
coupling to plasmon resonance in the gold nanospheres.
Several studies have shown that the separation critically
determines whether enhancement or quenching eventually
dominates [161, 169]. The enhancement of UC emission
is highly spectral dependent. A more-photon-involved UC
process resulted in a larger enhancement factor under the
same excitation power density [170].
The enhancement is known to result from the modification of the radiative and nonradiative decay rates and
the enhancement of the excitation intensity by the localized surface plasmonic resonance of metallic nanostructures [163, 168, 171, 172]. However, some details remain
unknown, e.g., is the change of decay rates or the increase
of excitation intensity contributing the most to the enhancement of the intrinsic QY of the nanocrystals, and by how
much the ultimate (maximum achievable) intrinsic QY will
be increased due to the change of the decay rates (the enhancement of the excitation intensity will not change the
ultimate intrinsic QY). Further experimental studies on the
QY enhancement are needed. Regarding the QY of the
system, although it always tends to increase due to the enhanced local field, the enhancement is, however, related
with the intrinsic power-density dependent QY of the crystal. Thus the understanding of the influence of the surface
plasmonic coupling on the intrinsic QY could address questions related to the amount of gain that can be expected.
Although the UC emission enhancement depends on
the fourth or higher power of the local field, solid theoretical analysis revealed much weaker gains than that of
Raman enhancement [163]. Furthermore, chemical sample
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
137
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
4. Surface modification
Similar to other nanosized biomarkers, UCNPs also need to
be water dispersible, highly stable, biocompatible, sensitive
and biotargeting for successful bio-applicability. However,
most of the commonly used UCNPs were synthesized using
oil-phase methods which render UCNPs with no intrinsic
aqueous dispersibility and functional organic groups on the
surface. Recently, several single-step synthesis methods,
mainly including polyol process [173,174], one-pot synthesis assisted by hydrophilic ligands [175,176], hydrothermal
microemulsion synthesis [177] and ionic liquid-based synthesis [94, 178], have been developed to directly synthesize
hydrophilic UCNPs. However, it is still very demanding to
obtain monodisperse and hydrophilic UCNPs with small
size via single-step methods, although they simplify the
reaction procedure and reduce the post-processing. For hydrophobic UCNP synthesis strategies, it is easier to control
the size, morphology, phase and crystallinity of the UCNPs,
and the procedure of post-processing is straight forward
and controllable. Furthermore, the directly synthesized hydrophilic UCNPs without any treatment also show toxicity
to cells, tissues or whole organisms. Therefore, UCNPs
synthesized from hydrophobic methods are still the most
commonly used in biological applications. In the following, more emphasis will be placed on how to re-process
the as-synthesized hydrophobic UCNPs into a hydrophilic
state with high stability, high biocompatibility and bioconjugation platform.
Ligand exchange is an effective method for modificating the surface of nanoparticles. When the hydrophobic ligands on the as-prepared nanoparticles are replaced
by some hydrophilic ligands, the nanoparticles will become water dispersible. Due to the surfactant OA ad-
Figure 6 (online color at: www.lpr-journal.org) FTIR spectra of (a) OA-capped NaYbF4 :Yb3+ /Er3+ sample, (b) MSANaYbF4 :Yb3+ /Er3+ and (c) PAH-MSA- NaYbF4 :Yb3+ /Er3+ .The
alternation of FTIR spectrum peaks (corresponding to the specific chemical bonds) revealed the process of surface modification.(Reprinted with permission from Ref. [139]. Copyright 2011,
American Chemical Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
138
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
673
4.1.3. Silanization
www.lpr-journal.org
The described surface modification methods have different applicabilities due to their distinctive features, and thus
the method selection should rely on the specific application of interest. As previously discussed, a large number
of successful bioapplications of UCNPs were substantially
demonstrated with the assist of these surface modification
methods. It is reasonable to believe that these methods are
essential and effective, and they hardly have considerable
side effect on the optical properties of UCNPs. However, it
is still of significance to investigate the impact of these processes on the luminescence intensity in order to improve or
maintain their brightness. In very recent years, there have
been an increasing number of studies revealing the impact
of surface modification on the optical proterties of UCNPs,
although they are not covered in most review papers. According to some reports, the luminescence intensity, to some
extent, was actually decreased after treatment with certain
modification processes [113, 115, 196, 197]. As a case of
point, decreased luminescence intensity can be observed
for PAA-coated UCNPs due to the interaction between the
surface of UCNPs and PAA [196, 197]. Besides, in the
silanization method the silica shell can slightly scatter both
excitation and emission light and thus weaken the luminescence to some extent [115]. The type of solvent is also
an important factor affecting the brightness. UCNPs dispersed in aqueous solvent were reported to have a decreased
brightness when compared to the same UCNPs dispersed
in organic solvent [113]. The reason is that water molecules
have high energy vibrational modes, which probably results
in an increased nonradiative relaxation of the excited states
and thus an overall quenching of the luminescence [113].
These problems could be overcome partially by a protective
shell of NaYF4 , which has low phonon energy and can act
as an isolation layer in order to greatly weaken the negative interaction between UCNPs and surface ligands/local
environment [113, 196]. The other strategies of enhancing
QY previously discussed in section 3 can also be exerted to
compensate the possible brightness decrease caused by surface modification and aqueous solvent. The key is to keep
or even improve the unique optical properties of UCNPs
while performing surface modification.
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
139
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
Despite all of the advantages of UCNPs as compared to conventional fluorophores, in order for them to be an attractive
choice for biological imaging studies, it is obviously of utmost importance that the toxicity and potential health issues
are thoroughly investigated. Even though UCNPs have not
been available for very long, the importance of this topic
has already attracted a great amount of interest. Most of
the studies have focused on in vitro cytotoxicity, however,
results concerning long-term effects within small animals
have recently started to appear. In this section, a few of
the most important results concerning the cytotoxicity of
UCNPs will be highlighted.
Following the increasing interest of UCNPs, several
studies of in-vitro cytotoxicity of UCNPs have been conducted. Already in 2008, Chatterjee et al. [203] used murine
stem cells to assess the cell biocompatibility of PEI-coated
UCNPs, while Shan et al. [204] studied the toxicity of silica
coated UCNPs in human osteosarcoma (HOS) cells, using
the methylthiazol tetrazolium (MTT) assay to evaluate the
cytotoxicity. Both of these early studies found the cytotoxicity to be very low. To date, the cytotoxicity of a large
number of both human and animal cell lines have been
studied, including HeLa cells [177], human glioblastoma
(U87MG) cells [199], human nasopharyngeal epidermal
carcinoma (KB) cells [205], and human hepatic (L02) cells
[206]. An example of the cell viability of HeLa, KB and
L02 cells are shown in Fig. 7. Currently, no severe adverse
effects have been found that can be directly related to the
UCNPs, indicating them to be of high biocompatibility.
An important aspect that in-vitro cell studies cannot provide an immediate answer to is the biodistribution and clearance of UCNPs following an injection into an animal. The
long-term biodistribution of intravenously injected UCNPs
has been reported by Xiong et al. [205]. In this study, the
animals were followed for 115 days and the results showed
that the UCNPs mainly accumulate within the spleen and
liver. Furthermore, the UCNPs were found to have a clearing time longer than 7 days, in contrast to a previous study
which showed a more rapid clearance speed [186]. This
shows that sample preparation is of high importance in order to determine the pharmacokinetics of the UCNPs within
an organism. However, perhaps more importantly, for these
Figure 7 (online color at: www.lpr-journal.org) In vitro cell viability of HeLa, KB and L02 cells incubated with UCNP-OA-CDT
at different concentration for 5 (a) and 24 h (b), respectively.
(Reprinted with permission from Ref. [206]. Copyright 2008, Elsevier Ltd.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
140
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
675
6. Applications of UCNPs
In recent years, UCNPs have attracted remarkable attention in the biophotonics area due to their merits of autofluorescence free, large anti-Stokes shifts, sharp emission
bandwidths, high resistance to photobleaching, nonblinking
emission behavior, deep detection ability and high spatial
resolution. As shown in Fig. 8, UCNPs have widely been
employed in in vitro cell microscopy, in vivo animal diffuse imaging, luminescence diffuse optical tomograpy, in
vivo multimodal animal imaging (MRI/PET), highly sensitive bioassays (luminescence resonance energy transfer
(LRET)), in vitro temperature sensing and photodynamic
therapy (PDT). In the following, all the above biological
applications will be covered and discusssed in detail.
6.1. Bioassays
Bioassay techniques are of fundamental importance in bioanalytical chemistry and biological sciences. They can of-
www.lpr-journal.org
Figure 8 (online color at: www.lpr-journal.org) Schematic summary of biological applications using UCNPs.
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
141
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
where I1 and I2 are the integrated intensities of the emissions from the higher state and the lower state, respectively;
C is a constant which depends on the degeneracy, spontaneous emission rate, and photon energies of the emitting
states in the host materials; k is the Boltzmann constant,
and T is the absolute temperature. Up to date, several emission bands of different activators have been used in temperature sensing, including the green emission bands of Er3+
[219225] and Ho3+ [226], red emission bands of Er3+
[227], and blue emission bands of Tm3+ [228]. In addition, the ratio between two well-separated emission bands
has also been used as temperature-sensitive measurables
[218,229]. A theoretical model is, however, needed to evaluate such information.
Recently, optical thermometry based on UCNPs has
been used in cell models. Vetrone et al. [217] reported
the use of green UC emissions from NaYF4 :Yb3+ , Er3+
nanoparticles for temperature sensing in HeLa cervical cancer cells (Fig. 10(a)). In this study, the excitation light at
920 nm with excitation intensity well below 0.5 kW/cm2
was used in order to avoid any pump-induced heating.
Fischer et al. [218] reported the use of the same kind of
nanoparticles for temperature sensing in human embryo
kidney cells (Fig. 10(b)).
It is worth pointing out that the ratiometric optical thermometry using UCNPs, although reliable, is primarily applicable for superficial imaging. When the UCNPs are embedded in biological tissue, the luminescence light will be
spectrally distorted by absorption and scattering of the tissue. The intensity ratio of two emission bands can thus not
directly be related to a temperature. In addition, considering different power-density dependencies of various UC
emission bands, even in superficial imaging, optical thermometry based on the use of two well-separated UC bands
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
142
(3)
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
677
Figure 10 (online color at: www.lpr-journal.org) Temperature sensing using NaYF4 :Yb3+ , Er3+ nanoparticles in cell models. (a) Top:
Optical transmission images of an individual HeLa cell at three inner temperatures. Cell death is observed at 45 C. Bottom: Temperature
of the HeLa cell, determined by the fluorescence intensity ratio of the NaYF4 :Yb3+ , Er3+ UCNPs at 525 and 545 nm, as a function
of the applied voltage. (b) Temperature-dependent images of human embryo kidney (HEK) cells transfected with NaYF4 :Yb3+ , Er3+
UCNPs showing sub-micrometer resolution. (c) Calibration plot for the temperature-sensitive nanoparticles inside HEK cells (Adapted
with permission from ref. [217] and [218]. Copyright 2010, 2011, Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, and American
Chemical Society.)
Due to the multiphoton excitation process involved, UCNPbased microscopy can yield high resolution images under
CW excitation. In most cases, a Gaussian beam of light
is used to excite the emissive samples and this Gaussian
intensity profile (unsaturated level) of the excitation beam
can be expressed as
Iex (r ) = I0 exp
2r 2
02
(4)
www.lpr-journal.org
Figure 11 (online color at: www.lpr-journal.org) Calculated photoluminescence intensity profiles upon excitation by a Gaussian
beam for three cases: one-photon (conventional fluorophores),
two-photon (red and green emission of UCNP) and three-photon
luminescence (blue emission of UCNP).
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
143
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
144
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
679
Figure 13 (online color at: www.lpr-journal.org) (a) MSA-UCNPincubated non-specific cell imaging; (b)(d) Charge-dependent
cellular imaging (PEI-UCNPs, PVP-UCNPs, PAA-UCNPs); photoluminescence images (upper row) and bright field images (bottom row).(Reprinted with permission from Ref. [82]. Copyright
2008; [236], Copyright 2009, American Chemical Society.)
cell membrane, as shown in Fig. 14. In their report, the assynthesized OA-capped UCNPs were first rendered aqueous dispersible through MSA encapsulation. Then negatively charged MSA-UCNPs were further polymer-coated
through physical adsorption by positively charged PAH,
which allows for antibody protein bioconjugation. Other
representative work have been reported by Zhang et al.
[177, 187], where FA-modified UCNPs were introduced
to specifically bind to the folate receptor overexpressed by
cancer cells for targeted imaging. High-affinity polypeptide
have in recent years been shown to be effective agents for
probing biological systems with high specificity. Utilized as
peptides, RGD-conjugated UCNPs have been successfully
used for cell targeted imaging [199,237]. Yu et al. [202] conjugated PEI-coated NaYF4 : Yb3+ , Er3+ /Ce3+ with recombinant chlorotoxin a typical peptide neurotoxin that could
bind with high specificity to many types of cancer cells
and incubated the modified UCNPs with C6 glioma cells for
targeted imaging. Utilized as target agents, UCNPs could
also be nanocarriers and indicators for some molecules
and drugs. In a recent report, Zhang et al. [238, 239]
used UCNPs for the intracellular investigation of smallinterference (siRNA) in living cells. Their results have
shown that UCNPs are capable of delivering and tracking
siRNA.
Figure 14 (online color at: www.lpr-journal.org) In vitro cancer cell imaging using 915-nm laser-excited UCNPs: images of HeLa cells
separately incubated with (a) blank, (b) anti-CEA8-NaYbF4 :Yb3+ /Ho3+ , (c) anti-CEA8-NaYbF4 :Yb3+ /Er3+ , (d) NaYbF4 :Yb3+ /Er3+ , and
(e) anti-CEA8-NaYbF4 :Yb3+ /Er3+ in the presence of 10-fold excess unlabeled antiCEA8. Bright field images (upper-row), photoluminescence images (middle-row) and superimposed images(bottom-row). (Reprinted with permission from Ref. [139]. Copyright 2011,
American Chemical Society.)
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
145
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
excitation light around 920 nm (where the optical absorption coefficient of water is much lower and the excitation
can still occur effectively) to replace 980-nm light in order
to avoid cell damage [139].
Figure 15 (online color at: www.lpr-journal.org) Calculated results of cell-in-celldish imaging model: excitation laser induced
temperature elevation ( C) after (a, b) 30 s and (c, d) 5 min
of irradiation with (a, c) a 915 nm laser and (b, d) a 980 nm
laser(power 300mW).(Adapted with permission from Ref. [139].
Copyright 2011, American Chemical Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
146
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
681
luminescence diffuse optical imaging employing an exogenous contrast agent, the ever-present tissue autofluorescence will certainly deteriorate the signal-to-background
ratio and limit the sensitivity of the system, and can at the
same time cause severe artifacts in the reconstructed tomographs [261].
In order to obtain accurate representations of the fluorophore distribution as well as to increase the sensitivities, significant efforts have been invested to develop
methods to overcome the tissue autofluorescence. Suggested approaches include subtraction-methods which
model the tissue autofluorescence [262], time-domain separation using fluorophores with suitable lifetimes [263],
spectral unmixing and multispectral methods [264, 265],
transillumination and normalized approaches [266], and
large Stokes-shifting markers, such as QDs [265, 267].
Although these methods can reduce the effects of unwanted tissue autofluorescence, they are usually associated with complex instrumentation and heavy computational needs. The use of UCNPs can, however, completely
eliminate the tissue autofluorescence since the emission
from endogenous fluorophores is Stokes shifted, as shown
in Fig. 16 [19, 38, 199, 236]. As illustrated in Fig. 17,
Figure 17 (online color at: www.lpr-journal.org) LDOT reconstruction of two cylindrical luminescent targets in a tissue phantom. (a) Reconstruction using UCNPs. (b) Reconstruction using
rhodamine 6G. It is demonstrated that the presence of even a
weak autofluorescent background can cause severe artifacts in
the tomographic reconstructions due to the ill-conditioned mathematical formulation. (Reprinted with permission from Ref. [269].
Copyright 2009, American Institute of Physics.)
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
147
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
Figure 18 (online color at: www.lpr-journal.org) Experimental results demonstrating that due to the two-photon-dependent emission,
the use of UCNPs can significantly increase the ultimate spatial resolution obtainable in LDOT. The experiments were carried out in
tissue phantoms, and two glass tubes filled with either DY-781 dye or UCNPs were placed with different center-to-center separation to
evaluate the spatial resolution obtained. (Reprinted with permission from Ref. [151]. Copyright 2012, American Chemical Society.)
Figure 19 (online color at: www.lpr-journal.org) Simulated sensitivity profiles for three different source positions and a fixed
detector position. The sharp gradients associated with UCNPs
(quadratic luminescence markers), enable LDOT reconstructions
with superior spatial resolutions as compared with conventional
fluorophores (linear luminescence markers). (Reprinted with permission from Ref. [151]. Copyright 2012, American Chemical Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
148
(5)
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
683
Figure 20 (online color at: www.lpr-journal.org) Experimental results demonstrating that a dual-beam excitation scheme can be
used for UCNPs to extract additional spatial information facilitating FDOT reconstructions of higher qualities. (a) Single-beam
setup. (b) Dual-beam setup. (Reprinted with permission from Ref.
[271]. (Copyright 2010, Optical Society of America).
www.lpr-journal.org
Figure 22 (online color at: www.lpr-journal.org) In vivo temperature ( C) monitoring for (a) a 915-nm laser irradiated mouse
and (b) a 980-nm laser irradiated mouse and (c) the corresponding temperature line profiles (curve A and B in (c) corresponding to line A in (a) and line B in (b), respectively). The experiment was performed under a stable room temperature of 18
C. (Reprinted with permission from Ref. [139]. Copyright 2011,
American Chemical Society.)
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
149
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
150
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
685
using histology since the staining agents used typically operate within the UVblue range of the spectrum. As known,
light at these wavelengths has extremely high attenuation
and thus in-vivo imaging is difficult except for very superficial targets. In addition, these dyes are highly sensitive to
photobleaching, which sets a limit on the time span for longitudinal studies in contrast to the stable UCNPs. Idris et al.
[112] first demonstrated the use of UCNPs for both in-vitro
and in-vivo tracking of stem cells over a time span of more
than one week using confocal microscopy. On the macroscopic scale, Cheng et al. [236] have shown that human
cancer cells can be tracked after injection into an animal
and the further development of the tumor itself can also be
monitored (Fig. 25).
www.lpr-journal.org
6.5.1. Dosimetry
For deep-tissue PDT with the currently available photosensitizers, it is of utmost importance to understand the lighttissue interaction in order to provide a useful dosimetry
[288290]. Various kinds of photosensitizers are available,
and we do not intend to review the current state-of-theart photosensitizers. Instead, using simple calculations, we
would like to highlight the importance of understanding
light-tissue interaction in order to estimate the outcome of
a treatment.
For simplicity, we will assume an infinite homogeneous
medium and a conservative light-dose threshold of 1 J/cm2
to reach a satisfactory treatment. The distribution of light
for a source of power P0 for steady state diffusion is given
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
151
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
by [291]
P0 eff r
e
(r ) =
,
4 Dr
(6)
em
,
ex
(7)
where denotes the QY, and the subscripts ex and em denote the excitation light and the emission light, respectively.
Furthermore, we can estimate the optical properties at the
excitation wavelength of UCNPs at 975 nm [245,292] to be
a (975nm) = 0.7 cm1 and s (975nm) = 5 cm1 . Using
these values we can then again calculate the treatment time
needed to reach a threshold dose for UCNP-mediated PDT.
Figure 26(a) shows the calculated fluence rates of the
excitation light at 660 nm and 975 nm, with the inset showing the corresponding QYs for each excitation depth for the
UCNPs. The strength of the source was chosen to be 500
mW and the QY of the UCNPs was obtained by a linear
fit (in a log-log scale) of the low-intensity values from previously published work [151]. The time needed to reach a
threshold dose for direct excitation of the photosensitizer
and UCNP-mediated PDT is shown in Fig. 26(b). A few
conclusions can be immediately drawn from these results.
Firstly, for prostate tissue, the attenuation factor of light at
660 nm and 975 nm is similar. This certainly depends on
the tissue type, with prostate tissue being relatively rich of
hemoglobin and water. However, even for tissue types with
less amount of water and hemoglobin, the difference in the
fluence rates between the two excitation wavelengths will
still be relatively negligible for depths less than 1 cm. Perhaps the most remarkable result is the treatment time that
will be needed if UCNPs are used to excite the photosensi-
tizer drugs. From Fig. 26, it can be seen that the treatment
time while using UCNPs will be at least 4 orders of magnitude longer. This treatment time can be shortened by either
increasing the excitation fluence rates (for instance by employing a pulsed excitation source) or by increasing the QY
of the UCNPs. However, it is quite clear that it is unrealistic to increase either of these factors by several orders of
magnitudes directly, especially since the excitation power
is already chosen to be very high. It may be possible to use a
pulsed light source to keep the average power down, while
using the peak power to excite the UCNPs. Pulsed light
will allow higher fluence rates in the tissue without causing
tissue heating. Ideally one would like to work with fluence
rates well above the saturation of the UCNPs, meaning that
the QY ideally would be in a region where it is not power dependent. For tissue regions with such high fluence rates, the
PDT efficiency will not decrease with depth as quickly, and
would have a depth profile more similar to the case of direct excitation of a photosensitizer. However, this approach
most probably cannot fully close the gap between direct
excitation and UCNP-mediated excitation, as the highest
reported QY (obtained under very high intensities) for UCNPs are still on the order of a few percents [50, 151]. Thus,
the treatment time while using UCNPs will still be longer
than direct excitation for most practical cases.
The efficiency with depth depends on two things, attenuation of the excitation light (in favor for long excitation
wavelengths) and the fact that the excitation is a two-photon
process for UCNPs and a linear process for conventional
PDT drugs (in favor for the linear excitation process). Our
simulations clearly indicate that the non-linear excitation
process hampers the excitation efficiency more than the
gain in better light penetration. This is a very similar observation as for molecular two-photon PDT sensitizers. The
calculations above are performed under extremely simplified conditions. Clearly, they are not meant to serve as an
accurate representation of the reality, but rather to highlight
the trends and needs in order to make UCNP-mediated PDT
feasible. The quality and efficiency of UCNPs have been
increasing in a very rapid fashion over the last years. Even
though the results above may indicate that it is difficult to
perform PDT with UCNPs directly, with the ever-increasing
knowledge of the UCNPs, and the use of properly pulsed
light sources, UCNP-mediated PDT may still be feasible
within the foreseeable future.
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
152
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
687
Figure 26 (online color at: www.lpr-journal.org) (a) Calculated fluence rates within human prostate tissue under excitation by 660 nm
and 975 nm, respectively, based on the assumptions given in the text. The inset shows the corresponding QY for the UCNPs (the QY
for the photosensitizer at 660 nm is constant) as a function of depth. (b) The treatment time needed to reach the light threshold dose
of 1 J/cm2 .
interest of enabling optical molecular imaging for multimodal contrast agents [297,298]. The high sensitivity, compactness and speed of optical imaging systems have made
luminescence molecular imaging extremely attractive as a
compliment to the more established modalities. However,
luminescent contrast agents that are optically bright and
stable are difficult to find. Thus the luminescent probes on
the multimodal contrast agents tend to degrade and bleach,
making it difficult to perform truly longitudinal studies in
parallel with other modalities. UCNPs, on the other hand,
are extremely stable and do not in general suffer from
bleaching and aging effects. In addition, since UCNPs are
based upon crystal hosts, multimodality imaging can, in
some cases, be enabled by modifying the crystal host.
Diffuse optical imaging is often associated with poor
spatial resolution, although the use of UCNPs can significantly improve this characteristic [151]. The spatial resolution of CT and MRI are in general higher, however, the
information obtained is quite different. While CT and MRI
can provide detailed anatomical information, for molecular
imaging, luminescent optical contrast agents usually have
a much higher sensitivity, comparable with PET. However,
while PET can provide detailed pre-operative data, they
are not practical for acquiring intra-operative data due to
both complexity and the relatively short half lives of the radioactive tracers. In this section, the current state-of-the-art
multimodality imaging agents which are based on UCNPs
will be discussed.
UCNPs have been demonstrated for both T1 -and T2 weighted MRI. For T1 imaging, the UCNP crystal can be
modified to provide contrast in an MR system directly. For
MRI, gadolinium (Gd) has been an extensively used re-
www.lpr-journal.org
Nuclear imaing is extensively used within the field of biology and medicine. PET imaging, for example, can achieve
sensitivities down to picomolar range. Recently, the use
of multimodal UCNPs for PET imaging was demonstrated
[206, 301], thereby the commonly used 18 F isotope was
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
153
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
chosen as the PET agent and very promising in vivo experiments were conducted. In addition, Zhou et al. used
a gadolinium-based crystal host, thus simultaneously enabling tri-modal imaging (further discussed below). Such
multimodal agents may prove to be very useful in future
clinical development. However, due to the limited half-lives
of the radioactive tracers, it is important to develop reaction
processes with high yield. Another aspect that is worth to
consider is information orthogonality. It is obvious that the
largest gain from a certain contrast agent can be obtained
if each modality provides independent information as compared with the other modalities. For the case of PET and UC
emission, although the information space may overlap, the
time window for acquiring PET data is much more limited
as compared with UC emission, however, PET may be able
to provide more detailed data.
make these particles compatible with currently used labeling and imaging technology and to be employed in important biomedical studies. This will require reproducible and
commercially available nanoparticles in user-friendly kits
and in large quantities. Also corresponding imaging systems must be developed for such applications. Apart from
these obvious technical aspects, a few research challenges
remain to be solved in order to fully explore the potential
of these very promising particles as contrast agents. These
challenges as well as ideas about how these could be addressed will be discussed below.
One challenge despite the strong and rapid development of the particles is their still relatively low QYs.
We believe that this can be partly overcome by utilizing
pulsed excitation with pulse lengths that match the relatively long lifetimes of UCNPs, as indicated in our discussion about QY in section 3. A similar case is the use of
femtosecond pulsed excitation in two-photon fluorescence
microscopy. Since it is known that the QY is power density
dependent as clearly illustrated in Fig. 5 consequently,
a high QY could be achieved by increasing the fluence rate
of the excitation light. In order to limit the heating of the
tissue to acceptable levels, strictly speaking, the only possibility is to utilize pulsed excitation with an acceptable
average power. A difference from the two-photon fluorescence microscopy is that the excitation of UCNPs relies
on long-lived intermediate states, which mean that pulses
shorter than a few milliseconds do not improve the QY
substantially. By utilizing pulsed excitation we foresee that
improved sensitivity of deep tissue targets could be accomplished without much side-effects in terms of heating.
This is a simple improvement which could be realized with
off-the-shelf diode lasers, thereby drastically increasing the
applicability of UCNPs for diffuse imaging.
For PDT applications, the potential of using UCNPs
would also be drastically improved by pulsed excitation.
The presented UCNP-PDT approaches seem to be severely
limited by the relatively low efficiency in generating cytotoxic agents in the PDT process. The low efficiency in
the excitation process using UCNP is here a limiting factor.
The possibility to produce a sufficient PDT effect could
probably be drastically improved by utilizing millisecond
pulsed excitation. We suggest to further study this approach
and to investigate its potential for improved efficiency. It
seems more feasible to reach the threshold dose with improved QY of the production of cytotoxic agents in the PDT
process via pulsed excitation.
The commonly used UCNPs employ Yb3+ ions as sensitizers, which have a strong absorption band at 975 nm.
However, this absorption band coincides with that of water which causes light attenuation and heating. To improve
the penetration depth and reduce the heating of tissue in
imaging applications, it is possible to shift the excitation
wavelength to the 920-nm absorption band of Yb3+ . Although this absorption band is significantly weaker than
the 975-nm band, excitation at this wavelength has been
demonstrated to be very feasible and motivated for in
vitro cell imaging, and imaging of very deeply embedded
targets within tissues. The conditions determining which
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
154
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
689
Appendix
List of abbreviations
AFM
atomic force microscopy
CR
cross relaxation
CS
cooperative sensitization
CT
computed tomography
CTAB cetyltrimethylammonium bromide
CTX
chlorotoxin
DOI
diffuse optical imaging
www.lpr-journal.org
ED
EDC
EDTA
ESA
ETU
FA
FDOT
FRET
FTIR
GFP
GSA
IR
LDOT
LRET
LSS
MD
MPR
MRI
MTT
NADH
NHS
OA
ODE
OM
PA
PAA
PAH
PBS
PCI
PDT
PEG
PEI
PET
QD
QY
RE
RET
RGD
SNR
SPECT
TEOS
TOP
TOPO
TSC
UC
UCL
UCNP
UV
electric dipole
ethyl dimethylaminopropyl carbodiimide
ethylenediamine tetraacetic acid
excited state absorption
energy transfer upconversion
folic acid
fluorescence diffuse optical tomography
Forster resonance energy transfer
fourier transform infrared
green fluorescent protein
ground state absorption
infrared
luminescence diffuse optical tomography
luminescence resonance energy transfer
liquid-solid-solution
magnetic dipole
multiphonon relaxation
magnetic resonance imaging
methylthiazol tetrazolium
nicotinamide adenine dinucleotide
N-hydroxysulfosuccinimide
oleic acid
octadecene
oleylamine
photon avalance
polyacrylic acid
polyallylamine hydrochloride
phosphate buffer solution
photochemical internalization
photodynamic therapy
polyethylene glycol
polyethylenimine
positron emission tomography
quantum dot
quantum yield
rare-earth
resonance energy transfer
arginineglycineaspartic acid
signal-to-noise ratio
single-photon emission computed tomography
tetraethyl orthosilicate
trioctyphosphine
trioctylphosphine oxide
trisodium citrate
upconversion
upconversion luminescence
upconverting nanoparticles
ultraviolet
Acknowledgements. The authors would like to thank any collaborator in this field for nice and fruitful collaboration, including Pontus Svenmarker, Johan Axelsson, Dan Wang, Fuhong
Cai, Ola Jakobsson, Thomas Laurell, Maria Messing, Reine Wallenberg, Bjorn
Dumlupinar and the dedi Thomasson, Gokhan
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
155
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
References
Gabriel Somesfalean received the Ph.D.
degree in Engineering Physics from Lund
University, Sweden, in 2005. Thereafter
he worked with optical environmental
monitoring and biophotonics at Harbin Institute of Technology, Lund University and
Zhejiang University, where he became an
Associate Professor in 2007. Currently,
he is the Vice Director of the Joint Research Center of Photonics between the Royal Institute of
Technology, Zhejiang University and Lund University. His current research interests include applied spectroscopy towards
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
156
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
[7]
[8]
[9]
[10]
[11]
[12]
[13]
[14]
[15]
[16]
[17]
[18]
[19]
[20]
[21]
[22]
[23]
[24]
[25]
[26]
[27]
[28]
691
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
157
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
158
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
693
[105] F. Vetrone, V. Mahalingam, and J. A. Capobianco, Chemistry of Materials 21(9), 18471851 (2009).
[106] X. Wang, J. Zhuang, Q. Peng, and Y. D. Li, Nature
437(7055), 121124 (2005).
[107] J. H. Zeng, Z. H. Li, J. Su, L. Wang, R. Yan, and Y. Li,
Nanotechnology 17(14), 3549 (2006).
[108] L. Wang and Y. Li, Nano Letters 6(8), 16451649
(2006).
[109] F. Zhang, J. Li, J. Shan, L. Xu, and D. Zhao, Chemistry A
European Journal 15(41), 1101011019 (2009).
[110] Z. Li and Y. Zhang, Nanotechnology 19(34), 345606
(2008).
[111] H. S. Qian and Y. Zhang, Langmuir 24(21), 1212312125
(2008).
[112] N. M. Idris, Z. Q. Li, L. Ye, E. K. W. Sim, R. Mahendran,
P. C. L. Ho, and Y. Zhang, Biomaterials 30(28), 51045113
(2009).
[113] J. C. Boyer, C. J. Carling, B. D. Gates, and N. R. Branda, J
Am Chem Soc 132(44), 1576615772 (2010).
[114] F. Wang, Y. Han, C. S. Lim, Y. Lu, J. Wang, J. Xu, H.
Chen, C. Zhang, M. Hong, and X. Liu, Nature 463(7284),
10611065 (2010).
[115] Z. Li, L. Wang, Z. Wang, X. Liu, and Y. Xiong, The Journal
of Physical Chemistry C 115(8), 32913296 (2011).
[116] H. Q. Wang and T. Nann, ACS Nano 3(11), 38043808
(2009).
[117] C. Chen, L. D. Sun, Z. X. Li, L. L. Li, J. Zhang, Y. W. Zhang,
and C. H. Yan, Langmuir 26(11), 87978803 (2010).
[118] C. Mi, Z. Tian, C. Cao, Z. Wang, C. Mao, and S. Xu,
Langmuir 27(23), 1463214637 (2011).
[119] D. Wang, L. Ren, X. Zhou, X. Z. Wang, J. Zhou,
Y. Han, and N. Kang, Nanotechnology 23(22), 225705
(2012).
[120] J. B. Edel, R. Fortt, J. C. deMello, and A. J. deMello, Chem.
Commun. (10), 11361137 (2002).
[121] E. M. Chan, A. P. Alivisatos, and R. A. Mathies, Journal
of the American Chemical Society. 127(40), 1385413861
(2005).
[122] X. Z. Lin, A. D. Terepka, and H. Yang, Nano letters 4(11),
22272232 (2004).
[123] J. Wagner and J. M. Kohler, Nano letters 5(4), 685691
(2005).
[124] Y. Song, J. Hormes, and C. S. S. R. Kumar, Small 4(6),
698711 (2008).
[125] H. C. Liu, O. Jakobsson, C. T. Xu, H. Y. Xie, T.
Laurell, and S. Andersson-Engels, Colloidal Quantum
Dots/nanocrystals For Biomedical Applications Vi 7909,
790917 (2011).
[126] A. D. Ostrowski, E. M. Chan, D. J. Gargas, E. M. Katz, G.
Han, P. J. Schuck, D. J. Milliron, and B. E. Cohen, ACS
Nano 6(3), 26862692 (2012).
[127] G. S. Yi and G. M. Chow, J. Mater. Chem. 15(41), 4460
4464 (2005).
[128] J. Shan and Y. Ju, Nanotechnology 20(27), 275603 (2009).
[129] P. Ghosh and A. Patra, The Journal of Physical Chemistry
C 112(9), 32233231 (2008).
[130] Z. Wang, F. Tao, L. Yao, W. Cai, and X. Li, Journal of
Crystal Growth 290(1), 296300 (2006).
[131] X. Yu, M. Li, M. Xie, L. Chen, Y. Li, and Q. Wang, Nano
Research 3, 5160 (2010), 10.1007/s12274-010-1008-2.
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
159
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
160
www.lpr-journal.org
Paper III
REVIEW
ARTICLE
Laser Photonics Rev. 7, No. 5 (2013)
695
www.lpr-journal.org
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
161
Upconverting nanoparticles for pre-clinical diffuse optical imaging, microscopy and sensing: Current
trends and future challenges
[254]
[255]
[256]
[257]
[258]
[259]
[260]
[261]
[262]
[263]
[264]
[265]
[266]
[267]
[268]
[269]
[270]
[271]
[272]
[273]
[274]
[275]
[276]
[277]
[278]
[279]
[280]
[281]
[282]
[283]
[284]
[285]
[286]
[287]
[288]
[289]
[290]
[291]
[292]
[293]
[294]
[295]
[296]
[297]
[298]
C 2013 The Authors. Laser Photonics Rev. published by Wiley-VCH Verlag GmbH & Co. KGaA Weinheim.
162
www.lpr-journal.org
Paper III
163
Paper IV
Balancing power density based quantum yield
characterization of upconverting nanoparticles for
arbitrary excitation intensities
H. Liu, C. T. Xu, D. Lindgren, H. Xie, D. Thomas, C. Gundlach,
S. Andersson-Engels.
Nanoscale 5, 4770-4775 (2013).
Paper IV
Nanoscale
View Article Online
PAPER
process. Here we show that the QYs of Yb3+ sensitized two-photon upconversion emissions can be well
characterized by the balancing power density, at which the ETU rate and the linear decay rate have
equal contributions, and its corresponding QY. The results in this paper provide a method to fully
DOI: 10.1039/c3nr00469d
describe the QY of upconverting nanoparticles for arbitrary excitation power densities, and is a fast and
www.rsc.org/nanoscale
simple approach for assessing the applicability of UCNPs from the perspective of energy conversion.
Introduction
Division of Solid State Physics, Department of Physics, Lund University, P. O. Box 118,
S-221 00 Lund, Sweden
Electronic supplementary
10.1039/c3nr00469d
information
(ESI)
available.
See
DOI:
been reported that UCNPs could have QYs on the order of a few
percent at high excitation intensities where the UCNPs are
saturated,14 while the QYs could decline by many orders of
magnitude when they are used at low excitation intensities.14
Obviously, such power density dependent QYs need to be
properly evaluated in order to assess the applicability of UCNPs
in specic biomedical areas. In spite of the interest, this crucial
issue has not been addressed in any satisfactory manner,
neither theoretically nor experimentally.22 To date, the reports
on the QYs of UCNPs are surprisingly scarce in the literature,14,2326 and even in the few publications available, the QY
data are usually provided at a specic excitation intensity,
ignoring their power density dependency.2326 Although full QY
information can be obtained by extensive measurements at all
excitation intensities, obviously this approach is not ideal
because of the accompanying burden of such measurements. In
addition, large errors would be also introduced in the
measurements at high excitation intensities due to saturation
eects of the optical equipment, such as the attenuator and the
power meter typically needed for such measurements. This will
ruin the accuracy of the QY data. Hence, a better understanding
of the power density dependency of the QY for a particular
design of UCNPs and thus characterizing this in a convenient
way is highly desirable, and will be of major importance for the
future development and applications of UCNPs in general. This
is addressed in this paper.
The power density dependence of the QY is governed by the
competition between the two major relaxation mechanisms
involved at the intermediate energy state in the UC process, i.e.,
167
Balancing power density based quantum yield characterization of upconverting nanoparticles for arbitrary
excitation intensities
View Article Online
Paper
Experimental
2.1
3.1
168
Nanoscale
Fig. 1 Schematic energy level diagrams of the Yb3+ and activator ions and the
proposed UC mechanism following laser diode excitation at 975 nm. The variables
used in the text for the population densities of dierent levels are indicated within
the parentheses.
Paper IV
View Article Online
Nanoscale
Paper
(1a)
dN1
N1
C0 N0 NYb1 C1 N1 NYb1
0;
dt
s1
(1b)
dN2
N2
C1 N1 NYb1
0;
dt
s2
(1c)
NYb1 sYb1sNYb0r,
(2)
dlog I
1
:
1
dlog r
1 s1 $C1 sYb1 sNYb0 r
(4)
1
;
s1 $C1 sYb1 sNYb0
(5)
where the ETU rate and linear decay rate equally contribute to
1
the depopulation of state 1, i.e., C1 sYb1 sNYb0 rb , the power
s1
density dependence curve has a slope eciency of 1.5.
Based on eqn (3), the QY, h, of the two-photon UC emission
at any power density can be dened by
hh
C0 C1 sYb1 2 s2 =srad
N0 sNYb0 r
I
2
:
1
sNYb0 rhn
C1 sYb1 sNYb0 r
s1
(6)
(7)
(8)
hs
:
2
(9)
169
Balancing power density based quantum yield characterization of upconverting nanoparticles for arbitrary
excitation intensities
View Article Online
Paper
Nanoscale
Fig. 3 (a) The upconversion and (b) the infrared PL spectra of the core NaYF4:Yb3+,Tm3+ nanoparticles. The insets in (a) and (b) present the zoomed-in
spectra in order to visualize the emission peaks around 700 nm and 1200 nm,
respectively. Both spectra were recorded at a power density of 125 W cm2 under
excitation of a CW 975 nm laser diode.
Fig. 2 TEM images of (a) the core NaYF4:Yb3+,Tm3+ nanoparticles and (b) the
coreshell NaYF4:Yb3+,Tm3+@NaYF4 nanoparticles. (c) XRD patterns of the
synthesized core and coreshell UCNPs.
170
Paper IV
View Article Online
Nanoscale
Fig. 4 The power density dependencies of the NIR UC emission band at 800 nm
of the core NaYF4:Yb3+,Tm3+ nanoparticles (red diamonds) and the coreshell
NaYF4:Yb3+,Tm3+@NaYF4 nanoparticles (blue circles). The black solid lines
represent the tangents of the power density dependence curves.
Paper
quantitative QY measurements can be dramatically reduced.
Especially, measurements under harsh pump conditions (in a
saturating range) can be avoided because the balancing power
density is signicantly lower than the saturation power density.
Noticing that the QY starts to decline dramatically when the
excitation intensity decreases to below the balancing point
according to eqn (8), a low balancing power density implies that
a considerable QY can be achieved under mild pump conditions. In this sense, the determination of the balancing power
density can be used as a fast and simple approach to evaluate
the applicability of UCNPs in applications where low excitation
intensities are required, such as deep tissue imaging in biological applications. The merit by doing so is that no absolute
measurements on luminescence intensities need to be performed, since the balancing power density depends on the
trend of the luminescence intensity change instead of absolute
intensities.
4
longer lifetimes s1 and sYb1, caused by the protection of the
shielding layer epitaxially grown on the core nanoparticles.38,43
The QYs of the synthesized UCNPs were measured in a power
density range of 0.02720 W cm2 on a spectrouorometerbased setup reported in our previous work,14 as shown in Fig. 5.
At the balancing power density, the core and coreshell nanoparticles have QYs of approximately 0.45% and 1.2%, respectively. The ttings of the QY data with eqn (8) were subsequently
implemented with the parameter rb locked to the experimentally obtained values. As is seen, the QY data can be well tted
both for the core and coreshell nanoparticles, with the tted
maximum attainable QY of 0.91% for the core and 2.6% for
coreshell nanoparticles, which can be well estimated by the
twofold QY at the balancing power density.
One main advantage of the proposed approach for QY
characterization by providing (rb, 2hb) is that the number of
Conclusions
Acknowledgements
M. E. Messing and L. R. Wallenberg are gratefully acknowledged
for the help with the TEM measurements. J. Larsson is
acknowledged for the help with the XRD measurements.
D. Hessman is acknowledged for the help with infrared PL
measurements. We thank S. Fredriksson, F. Olsson, A. Gisselsson, G. Dumlupinar and X. Wu for the help with the
synthesis of the nanoparticles. This work was supported by a
grant from the Swedish Research Council (grant no. 621-20114265) and a Linneaus grant to the Lund Laser Centre.
References
Fig. 5 The quantum yields of the NIR UC emission band at 800 nm of the core
NaYF4:Yb3+,Tm3+ nanoparticles (red diamonds) and the coreshell NaYF4:Yb3+,Tm3+@NaYF4 nanoparticles (blue circles) at various excitation power
densities. The black solid lines stand for the tted data.
171
Balancing power density based quantum yield characterization of upconverting nanoparticles for arbitrary
excitation intensities
View Article Online
Paper
172
Nanoscale
24 Q. Liu, Y. Sun, T. Yang, W. Feng, C. Li and F. Li, J. Am. Chem.
Soc., 2011, 133, 1712217125.
25 A. D. Ostrowski, E. M. Chan, D. J. Gargas, E. M. Katz, G. Han,
P. J. Schuck, D. J. Milliron and B. E. Cohen, ACS Nano, 2012,
6, 26862692.
26 G. Chen, J. Shen, T. Y. Ohulchanskyy, N. J. Patel, A. Kutikov,
Paper V
Deep tissue optical imaging of upconverting
nanoparticles enabled by exploiting higher intrinsic
quantum yield through use of millisecond single
pulse excitation with high peak power
H. Liu, C. T. Xu, G. Dumlupinar, O. B. Jensen, P. E. Andersen,
S. Andersson-Engels.
Nanoscale 5, 10034-10040 (2013).
Paper V
Nanoscale
View Article Online
PAPER
depths and shorter data acquisition times compared with continuous wave excitation, while
simultaneously keeping the possible thermal side-eects of the excitation light moderate. These key
results facilitate means to break through the general shallow depth limit of upconverting-nanoparticle-
DOI: 10.1039/c3nr01917a
based uorescence techniques, necessary for a range of biomedical applications, including diuse
www.rsc.org/nanoscale
optical imaging, photodynamic therapy and remote activation of biomolecules in deep tissues.
Introduction
Department of Physics, Lund University, P.O. Box 118, S-221 00 Lund, Sweden.
E-mail: [email protected]; Fax: +46 46 222 4250; Tel: +46 46 222 7471
information
(ESI)
University
available.
of
Denmark,
See
DOI:
175
Deep tissue optical imaging of upconverting nanoparticles enabled by exploiting higher intrinsic quantum
yield through use of millisecond single pulse excitation with high peak power
View Article Online
Paper
increase enables us to implement single-shot imaging of
UCNPs, shortening data acquisition time by orders of magnitude while simultaneously improving imaging depth as
compared to CW excitation. These results have the potential to
fundamentally broaden the applicability of UCNPs in deep
tissue regions relying on diuse light excitation, breaking the
shallow-depth limitation in UCNP-based imaging.
Experimental
2.3
2.1
Synthesis of UCNPs
Coreshell NaYF4:Yb3+,Tm3+@NaYF4 nanoparticles, synthesized through a recently reported protocol,24 are used as a
representative of UCNPs in this work.
All the chemicals were purchased from Sigma-Aldrich and
used without further purication. The core nanoparticles
NaYF4:Yb3+,Tm3+ were rst synthesized using the protocol
reported in ref. 25. In a typical synthesis, anhydrous powders of
YCl3 (0.75 mmol), YbCl3 (0.25 mmol) and TmCl3 (0.003 mmol)
were dissolved in 6 mL oleic acid and 17 mL octadecene in a
250 mL ask at 160 C for 30 min. Aer the clear solution cooled
down to room temperature, 10 mL of a methanol solution containing 4 mmol NH4F and 2.5 mmol NaOH was added, and the
mixture was stirred for 30 min at 50 C. The methanol was
removed from the system by slowly heating it, and the resulting
solution was heated to 300 C for 1.5 h under argon atmosphere.
Aer the mixture cooled to room temperature, the nanoparticles
were precipitated with ethanol and washed with an ethanol
water mixture several times, and then redispersed in hexane to
form a nanoparticle suspension. The coreshell nanoparticles
were subsequently produced by slightly modifying the above
procedure through incorporation of the prepared core nanoparticles as the seeds in the synthesis.24 1 mmol YCl3 was solely
used to provide rare-earth ions for the shielding layer. Other steps
were kept the same as the synthesis of core nanoparticles.
Nanoscale
176
Fig. 1 Schematic energy level diagrams of Yb3+ and Tm3+ ions and the proposed
UC mechanism following the excitation at 975 nm. The variables used in the text
for the population densities of dierent levels are indicated within the
parentheses.
Paper V
View Article Online
Nanoscale
Paper
dNYb1 sr
NYb1
NYb0 C0 N0 C1 N1 C2 N2 NYb1
;
hn
dt
sYb1
(1a)
dN10
C0 N0 NYb1 b01 N10 ;
dt
(1b)
dN1
N1
b01 N10 C1 N1 NYb1
;
dt
s1
(1c)
dN20
C1 N1 NYb1 b02 N20 ;
dt
(1d)
dN2
N2
b02 N20 C2 N2 NYb1
;
dt
s2
(1e)
Table 1
Fig. 2 (a) Simulated power density dependence of the QY of the NIR UC emission in the steady state under CW excitation. (b) The temporally cumulative QYs
under CW excitation and under pulsed excitation in the rst pulse period. The
pulsed excitation had a duty cycle of 4% and a repetition rate of 2 Hz. Both the
CW and pulsed excitation approaches provided an average power density of 1 W
cm2. (c) The temporally cumulative QYs under CW excitation and under pulsed
excitation in multiple periods. The pulsed excitation had a xed duty cycle of 4%
and various repetition rates. All the excitation approaches provided the same
average power density of 0.1 W cm2.
s (cm )
N0 (cm3)
NYb0 (cm3)
sYb1 (ms)
s1 (ms)
s2 (ms)
C0 (cm3 s1)
C1 (cm3 s1)
C2 (cm3 s1)
b10 (s1)
b20 (s1)
1.69 1020a
1.25 1019b
1.52 1021b
1.32c
7.43d
1.49c
1.6 1018e
6.2 1016e
1.6 1018e
1.7 104d
1 105d
a
e
From Jiang et al.26 b Calculated (see the ESI for the calculation of ion concentrations). c From measurement (see Fig. S2). d From Ivanova et al.27
Estimated from Braud et al.28 and Ivanova et al.27 (see the ESI for the selection of the ETU rates).
177
Deep tissue optical imaging of upconverting nanoparticles enabled by exploiting higher intrinsic quantum
yield through use of millisecond single pulse excitation with high peak power
View Article Online
Paper
178
Nanoscale
Paper V
View Article Online
Nanoscale
Paper
Fig. 5 The dependency of the blue (at 474 nm) and red (at 644 nm) UC signal
gains by the pulsed excitation on the average power density. The pulsed excitation was set to have a pulse width of 20 ms and a repetition rate of 2 Hz.
Fig. 4 (a) The NIR UC signal gain by the pulsed excitation with dierent pulse
widths. The data were measured with an average excitation power density of
0.12 W cm2. Inset: the response of the NIR UC emission to a square-wave excitation with IRF denoting the instrument response function. (b) The dependence
of the NIR UC signal gain by the pulsed excitation on the average power density.
The pulsed excitation was set to have a 20 ms pulse width and 2 Hz repetition
rate. Inset: the average power density dependence of the NIR UC emission
intensity under CW and pulsed excitations.
179
Deep tissue optical imaging of upconverting nanoparticles enabled by exploiting higher intrinsic quantum
yield through use of millisecond single pulse excitation with high peak power
View Article Online
Paper
Nanoscale
Fig. 6 The NIR UC emission images taken for a luminescent inclusion located at a
depth of 10 mm under the (a) CW and (b) pulsed excitation, detected in a transillumination geometry. The average power density was 1.2 W cm2, and the
exposure time was 10 s for both excitation approaches. The NIR UC emission
images taken for a luminescent inclusion located at a depth of 13 mm under the
(c) CW and (d) 50 ms single pulse excitation, detected in an epi-illumination
geometry. Both excitation sources provided the maximum permissible power
density allowed by the ANSI standard for exposure to human skin. The exposure
time was 10 s for the CW excitation and 1 s for the single pulse excitation.
180
Conclusions
Acknowledgements
S. Fredriksson, F. Olsson, A. Gisselsson, and P. Kjellman are
gratefully acknowledged for the help with the synthesis of the
UCNPs and ICP-OES measurement. M. E. Messing and L. R.
Wallenberg are acknowledged for the help with the TEM and
XEDS measurements. B. Thomasson is acknowledged for the
help in the single-shot imaging experiment. S. Johansson and A.
Shaharin are acknowledged for the help with the pTOFS
measurements. We thank D. Kroon and D. Guenot for the
assistance in the laser power measurement, and J. Axelsson and
P. Svenmarker for their helpful discussions. This work was supported by a grant from the Swedish Research Council (grant no.
621-2011-4265) and a Linneaus grant to the Lund Laser Centre.
References
1 X. Wang, J. Zhuang, Q. Peng and Y. D. Li, Nature, 2005, 437,
121124.
2 H.-X. Mai, Y.-W. Zhang, R. Si, Z.-G. Yan, L.-D. Sun, L.-P. You
and C.-H. Yan, J. Am. Chem. Soc., 2006, 128, 64266436.
3 J.-C. Boyer, L. A. Cuccia and J. A. Capobianco, Nano Lett.,
2007, 7, 847852.
4 E. M. Chan, C. Xu, A. W. Mao, G. Han, J. S. Owen, B. E. Cohen
and D. J. Milliron, Nano Lett., 2010, 10, 18741885.
5 F. Wang, Y. Han, C. S. Lim, Y. Lu, J. Wang, J. Xu, H. Chen,
C. Zhang, M. Hong and X. Liu, Nature, 2010, 463, 10611065.
6 F. Zhang, G. B. Braun, A. Pallaoro, Y. Zhang, Y. Shi, D. Cui,
M. Moskovits, D. Zhao and G. D. Stucky, Nano Lett., 2012,
12, 6167.
Paper V
View Article Online
Nanoscale
7 L. Wang, R. Yan, Z. Huo, L. Wang, J. Zeng, J. Bao, X. Wang,
Q. Peng and Y. Li, Angew. Chem., Int. Ed., 2005, 44, 60546057.
8 S. F. Lim, R. Riehn, W. S. Ryu, N. Khanarian, C. K. Tung,
D. Tank and R. H. Austin, Nano Lett., 2006, 6, 169174.
9 M. Nyk, R. Kumar, T. Y. Ohulchanskyy, E. J. Bergey and
P. N. Prasad, Nano Lett., 2008, 8, 38343838.
10 H. P
akkil
a, M. Ylih
arsil
a, S. Lahtinen, L. Hattara,
N. Salminen, R. Arppe, M. Lastusaari, P. Saviranta and
T. Soukka, Anal. Chem., 2012, 84, 86288634.
11 M. Yu, F. Li, Z. Chen, H. Hu, C. Zhan, H. Yang and C. Huang,
Anal. Chem., 2009, 81, 930935.
12 C. T. Xu, Q. Zhan, H. Liu, G. Somesfalean, J. Qian, S. He and
S. Andersson-Engels, Laser Photonics Rev., 2012, DOI:
10.1002/lpor201200052.
13 J. Zhou, Z. Liu and F. Li, Chem. Soc. Rev., 2012, 41, 13231349.
14 N. M. Idris, M. K. Gnanasammandhan, J. Zhang, P. C. Ho,
R. Mahendran and Y. Zhang, Nat. Med., 2012, 18, 15801585.
15 F. C. J. M. van Veggel, C. Dong, N. J. J. Johnson and
J. Pichaandi, Nanoscale, 2012, 4, 73097321.
16 Laser Institute of America, Americal National Standard for
Safe Use of Lasers; in ANSI Z136.1-2000, Laser Institute of
America, Orlando, Florida, 2000.
17 C. T. Xu, N. Svensson, J. Axelsson, P. Svenmarker,
G. Somesfalean, G. Chen, H. Liang, H. Liu, Z. Zhang and
S. Andersson-Engels, Appl. Phys. Lett., 2008, 93, 171103.
18 Q. Zhan, S. He, J. Qian, H. Cheng and F. Cai, Theranostics,
2013, 3, 303316.
19 H. Liu, C. T. Xu and S. Andersson-Engels, Opt. Lett., 2010, 35,
718720.
20 P. Svenmarker, C. T. Xu and S. Andersson-Engels, Opt. Lett.,
2010, 35, 27892791.
21 F. Auzel, Chem. Rev., 2004, 104, 139173.
22 L. M. Maestro, E. M. Rodriguez, F. Vetrone, R. Naccache,
H. L. Ramirez, D. Jaque, J. A. Capobianco and J. G. Sol
e,
Opt. Express, 2010, 18, 2354423553.
Paper
23 C. F. Gainer, G. S. Joshua, C. R. De Silva and M. Romanowski,
J. Mater. Chem., 2011, 21, 1853018533.
24 H.-S. Qian and Y. Zhang, Langmuir, 2008, 24, 1212312125.
25 Z. Li and Y. Zhang, Nanotechnology, 2008, 19, 345606.
26 X. P. Jiang, Z. M. Yang, T. Liu and S. H. Xu, J. Appl. Phys.,
2009, 105, 103113.
27 S. E. Ivanova, A. M. Tkachuk, A. Mirzaeva and F. Pelle, Opt.
Spectrosc., 2008, 105, 228241.
28 A. Braud, S. Girard, J. L. Doualan, M. Thuau, R. Moncorg'e
and A. M. Tkachuk, Phys. Rev. B: Condens. Matter Mater.
Phys., 2000, 61, 52805292.
181
Paper VI
Autofluorescence insensitive imaging using
upconverting nanoparticles in scattering media
C. T. Xu, N. Svensson, J. Axelsson, P. Svenmarker, G. Somesfalean,
G. Chen, H. Liang, H. Liu, Z. Zhang, S. Andersson-Engels.
Applied Physics Letters 93(17), 171103-1 - 171103-3 (2008).
Paper VI
Department of Physics, Lund University, P.O. Box 118, S-221 00 Lund, Sweden
Department of Physics, Harbin Institute of Technology, P.O. Box 3025, Harbin 150080,
Peoples Republic of China
Received 27 August 2008; accepted 1 October 2008; published online 27 October 2008
Autofluorescence is a nuisance in the field of fluorescence imaging and tomography of exogenous
molecular markers in tissue, degrading the quality of the collected data. In this letter, we report
autofluorescence insensitive imaging using highly efficient upconverting nanocrystals
NaYF4 : Yb3+ / Tm3+ in a tissue phantom illuminated with near-infrared radiation of 85 mW/ cm2.
It was found that imaging with such nanocrystals leads to an exceptionally high contrast compared
to traditional downconverting fluorophores due to the absence of autofluorescence. Upconverting
nanocrystals may be envisaged as important biological markers for tissue imaging purposes. 2008
American Institute of Physics. DOI: 10.1063/1.3005588
In recent years, the interest for fluorescence diffuse optical imaging and tomography has grown tremendously.13
Much of the work has been focused on developing inexpensive and compact systems for macroscopic imaging of fluorophores embedded in small animals. The systems could, for
example, be used to monitor the effects from drugs on cancer
tumors over a period of time, without sacrificing the animals.
Extensive research has been performed both on the practical
and the theoretical side in this area. Currently, mainly traditional dyes are used as fluorophores. These dyes emit Stoke
shifted light upon excitation, and can be very effective with
quantum efficiencies close to unity. However, since tissue
itself autofluoresces due to several endogenous fluorophores
mainly in the skin, a background fluorescence from the bulk
tissue will always exist.4,5 In the presence of tissue autofluorescence, different spectral unmixing algorithms can also be
used to extract the signal. These algorithms utilize the spectral characteristics of the fluorophores and the
autofluorescence.6 However, the measurement procedure can
be quite complex since one needs to acquire the emission at
multiple wavelengths. To suppress the effects of autofluorescence, several approaches have been suggested. The two
most promising ones are based on quantum dots and upconverting markers.
Quantum dots have been used as fluorophores in a great
number of studies.79 Their large Stoke shift and narrow
emission spectra allow the emission to be detected in a spectral band with low tissue autofluorescence. Quantum dots are
very bright fluorophores due to their high absorption, mainly
in the UV region. Their emission wavelength is dependent on
the core size and can be selected over a wide range while the
same excitation wavelength can be used.7 The main drawback at this stage for these fluorophores is their toxicity,
together with the fact that penetration of the UV light is low
in tissue. Typical quantum dots are based on CdSe due to the
well established fabrication technology.10 Recent studies
have shown that these quantum dots are potentially harmful
a
0003-6951/2008/9317/171103/3/$23.00
93, 171103-1
Downloaded 27 Oct 2008 to 130.235.188.55. Redistribution subject to AIP license or copyright; see http://apl.aip.org/apl/copyright.jsp
185
171103-2
Xu et al.
1600
1400
1200
1000
800
1800
10
Slope = 2.0
10
10
10
1
600
10
Intensity (mW/cm)
400
200
0
500
600
700
Wavelength (nm)
800
FIG. 2. Color online Schematic of the imaging setup. Light from a fibercoupled laser is scanned in an array on the surface of the phantom. An air
cooled CCD camera is used to capture an image for every scanned position.
FIG. 3. Images comparing the DY-781 dye and the nanoparticles with and
without autofluorescence, along with plots showing the sums in the vertical
directions. The white dots have been added artificially and represent the
positions used for the excitation light. The left column shows the results
using DY-781, and the right column shows the results using upconverting
nanoparticles. a and b are taken without any added autofluorophores. c
and d are taken with a background autofluorophore concentration of
40 nM.
the surface of the tissue phantom. The second tube was filled
with a solution of ordinary fluorophores DY-781, Dyomics
GmbH dissolved in ethanol with a concentration of 1 M.
The fluorophores were excited with a 780 nm laser diode
with a power of 40 mW on the surface of the phantom. The
concentration of the nanoparticles is reasonably consistent
with studies performed using quantum dots in vivo. In those
studies, approximately 1 nmol of CdSe quantum dots were
injected into a mouse 18 g, giving an approximate concentration of 0.01 wt % if distributed homogeneously.8,9
With functionalized quantum dots, selective accumulation in
tumors can be achieved. Taking this into consideration and
the fact that the nanoparticles have molar mass of the same
order of magnitude as the quantum dots, the concentration of
1 wt % used in this study seems acceptable.
The lasers were raster scanned in a 4.4 4.4 cm2 array
consisting of 121 positions, and an image was acquired for
every laser position. In order to minimize the effects of random bright pixels, a median filter with a mask of 3
3 pixels was applied to all images. The individual images
were then summed, and a representation of the photon distribution on the surface was obtained. Even if this does not
accurately reflect the fluorophore distribution, it enables detection of fluorescent inclusions. To mimic tissue autofluorescence, a small amount of DY-781 was also added into the
phantom.
Figure 3 shows images taken with and without autofluorescence. As expected, the autofluorescence effectively hides
the signal from the inclusion with the ordinary fluorophores
Fig. 3c, while no background appears on the images using
the upconversion scheme Fig. 3d. It is worth to notice
that using traditional fluorophores, even without any artificial
Downloaded 27 Oct 2008 to 130.235.188.55. Redistribution subject to AIP license or copyright; see http://apl.aip.org/apl/copyright.jsp
186
Paper VI
171103-3
Xu et al.
Downloaded 27 Oct 2008 to 130.235.188.55. Redistribution subject to AIP license or copyright; see http://apl.aip.org/apl/copyright.jsp
187
Paper VII
High-resolution fluorescence diffuse optical
tomography developed with nonlinear upconverting
nanoparticles
C. T. Xu, P. Svenmarker, H. Liu, X. Wu, M. E. Messing,
L. R. Wallenberg, S. Andersson-Engels.
ACS Nano 6(6), 4788-4795 (2012).
Paper VII
Can T. Xu,,* Pontus Svenmarker, Haichun Liu, Xia Wu, Maria E. Messing, L. Reine Wallenberg, and
Stefan Andersson-Engels
ARTICLE
Department of Physics, Lund University, Box 118, S-221 00 Lund, Sweden, Division of Solid State Physics, Lund University, Box 118, S-221 00 Lund, Sweden, and
Polymer and Materials Chemistry/nCHREM, Lund University, Box 124, S-221 00 Lund, Sweden
ABSTRACT
Fluorescence diuse optical tomography (FDOT) is an emerging biomedical imaging technique that
can be used to localize and quantify deeply situated uorescent molecules within tissues. However, the
potential of this technique is currently limited by its poor spatial resolution. In this work, we
demonstrate that the current resolution limit of FDOT can be breached by exploiting the nonlinear
power-dependent optical emission property of upconverting nanoparticles doped with rare-earth
elements. The rare-earth-doped coreshell nanoparticles, NaYF4:Yb3/Tm3@NaYF4 of hexagonal
phase, are synthesized through a stoichiometric method, and optical characterization shows that the
upconverting emission of the nanoparticles in tissues depends quadratically on the power of excitation.
In addition, quantum-yield measurements of the emission from the synthesized nanoparticles are
performed over a large range of excitation intensities, for both core and coreshell particles. The
measurements show that the quantum yield of the 800 nm emission band of coreshell upconverting
nanoparticles is 3.5% under an excitation intensity of 78 W/cm2. The FDOT reconstruction experiments
are carried out in a controlled environment using liquid tissue phantoms. The experiments show that
the spatial resolution of the FDOT reconstruction images can be signicantly improved by the use
of the synthesized upconverting nanoparticles and break the current spatial resolution limits of FDOT
images obtained from using conventional linear uorophores as contrast agents.
KEYWORDS: upconversion nanoparticles . bioimaging . quantum yield .
resolution . diuse imaging
from the use of linear uorophores in diusive tissues. In this work, we report on the
exploration of rare-earth-doped coreshell
nanoparticles with upconverting optical
properties as contrast agents for FDOT and
demonstrate that the spatial resolution of
FDOT can be signicantly improved by exploiting their nonlinear power-dependent
optical properties. The experiments were
VOL. 6
NO. 6
* Address correspondence to
[email protected].
Received for review December 19, 2011
and accepted May 7, 2012.
Published online May 08, 2012
10.1021/nn3015807
C 2012 American Chemical Society
47884795
2012
4788
www.acsnano.org
191
ARTICLE
NO. 6
47884795
4789
2012
www.acsnano.org
192
Paper VII
VOL. 6
NO. 6
47884795
ARTICLE
Figure 2. Emission spectrum (blue solid line) for the synthesized NaYF4:Yb3/Tm3@NaYF4 nanoparticles of hexagonal
phase under excitation at 975 nm and attenuation spectrum
(green solid line) of typical murine muscle tissue. The
excitation intensity was 1 W/cm2, leading to a strong
emission band at 800 nm. The weak emission bands around
470 and 650 nm are less suitable for use within tissues, since
the attenuation at those wavelengths is signicantly higher,
as shown. The inset shows the power dependence of the
emission band at 800 nm, demonstrating that the emission
originates from a two-photon process.
4790
2012
www.acsnano.org
193
ARTICLE
Figure 3. Quantum yield of the 800 nm emission band for the synthesized hexagonal coreshell NaYF4:Yb3/Tm3@NaYF4
nanoparticles and hexagonal core NaYF4:Yb3/Tm3 nanoparticles. It can be seen that the quantum yield increases linearly
(slope of 1) with the excitation intensity until a saturation point, from which the quantum yield approaches a constant value,
which is the expected behavior for a linear uorophore. Solid lines show simulated corresponding quantum yields for highly
ecient two-photon dyes under identical experimental conditions, which in contrast to upconverting nanoparticles require
two photons to be simultaneously absorbed. Compared with upconverting nanoparticles, even two-photon dyes with cross
sections on the order of 100 kGM will be 108 times less bright.
the linear uorophores had a reduced scattering coecient of 0 s(785 nm) = 10.1 cm1 and an absorption
coecient of a(785 nm) = 0.51 cm1, while the phantom
for the upconverting nanoparticles had a reduced scattering coecient of 0 s(975 nm) = 10.1 cm1 and an
absorption coecient of a(975 nm) = 0.52 cm1, determined by a time-of-ight spectroscopy (TOFS) system.54
The experimental setup is schematically illustrated in
Figure 4, where two capillary tubes with inner diameters
of 2.0 mm, lled with either the DY-781 uorophores or
the synthesized upconverting nanoparticles, were used
to simulate uorescent targets. By adjusting the centerto-center distance of the two uorescent tubes and
performing one tomographic reconstruction for each
separation distance, the obtainable spatial resolution in
the reconstruction images could be evaluated. It is
important to note that autouorescence can cause severe artifacts in the reconstructions for linear Stokesshifting uorophores.5 Thus, a phantom material with
low autouorescence was selected, and any remaining
autouorescence eects or other background signals
were removed by subtracting from each image its corresponding image obtained from a scan without any
uorescent targets.
Cross-sectional slices from the tomographic reconstructions for dierent center-to-center distances of the uorescent tubes are shown in Figure 5. The reconstructions
for the linearly power-dependent DY-781 uorophores
XU ET AL.
are shown in Figure 5a to e, with decreasing tube separations from 10 mm in Figure 5a to 6 mm in Figure 5e.
The corresponding reconstructions for the quadratically
power-dependent NaYF4:Yb3/Tm3@NaYF4 upconverting nanoparticles are shown in Figure 5f to m, with
VOL. 6
NO. 6
47884795
4791
2012
www.acsnano.org
194
Paper VII
ARTICLE
Figure 5. Cross-sectional slices (two-dimensional plots) of the FDOT reconstructions with the linearly power-dependent
DY-781 uorophores and the quadratically power-dependent NaYF4:Yb3/Tm3@NaYF4 upconverting nanoparticles as
contrast agents and their corresponding intensity proles (line plots). The true depth was z = 7 mm. The separation distance
between the uorescent tubes was varied from 10 to 6 mm (step sizes of 1 mm) for the case of the linear uorophores, shown
in (a)(e). The use of quadratic upconverting nanoparticles clearly leads to reconstructions with higher spatial resolutions
and qualities. Thus, the separation distance between the uorescent tubes was varied from 10 to 3 mm (stepsizes of 1 mm),
shown in (f)(m). Using the linear uorophores, already at a separation distance of 8 mm (c), the two uorescent tubes can no
longer be separated. However, in the images obtained with the use of the quadratic upconverting nanoparticles, besides the
signicantly higher qualities of the reconstructions, the two uorescent tubes can still be separated at a separation distance of
5 mm (k).
NO. 6
47884795
4792
2012
www.acsnano.org
195
ARTICLE
Figure 6. Simulated sensitivity proles of linearly powerdependent uorophores and quadratically power-dependent
upconverting nanoparticles for three source positions and a
xed detector position. The circles represent the experimental
case of two uorescent tubes with a center-to-center separation of 5 mm. (ac) Sensitivity proles of linear uorophores.
(df) Sensitivity proles of the upconverting nanoparticles.
From the gure, it is clear that the quadratic power dependence of the emission results in sensitivity proles that are
much more conned and sharply dened as compared with
linear uorophores. Thus, the two uorescent targets can
eectively be selectively excited if upconverting nanoparticles
are employed as contrast agents.
EXPERIMENTAL METHODS
Synthesis and Characterization of the Nanoparticles. In the synthesis of the core NaYF4:Yb3/Tm3 nanoparticles, stoichiometric
amounts of YCl3 (0.75 mmol), YbCl3 (0.25 mmol), and TmCl3
(0.003 mmol) were mixed with 6 mL of oleic acid (OA) and 17 mL
of octadecene (ODE) in a 250 mL flask. The mixture was heated
to 160 C for 30 min, forming a clear solution. Then, 4 mmol of
NH4F (0.1482 g) and 2.5 mmol of NaOH (0.1 g) dissolved in 10 mL
of methanol was aspirated into a syringe and added into the
solution by a syringe pump using a flow rate of 200 L/min. The
obtained solution was first slowly heated to evaporate the
methanol and degassed at 100 C for 10 min, followed by a
second step of heating to 300 C for 1.5 h under an argon
atmosphere. After the solution had cooled to room temperature, the nanoparticles were precipitated and washed with
100 mL of ethanol/water mixture (1:1 v/v) three times and
centrifuged at 5100 rpm for 10 min before being collected
and redispersed in 6 mL of n-hexane. In order to reduce the
nonradiative losses caused by surface effects, the nanoparticles
XU ET AL.
VOL. 6
NO. 6
47884795
4793
2012
www.acsnano.org
196
Paper VII
NP
As FNP
s
ANP Fs
(1)
(2)
(3)
XU ET AL.
VOL. 6
NO. 6
47884795
ARTICLE
standard uorophores calibrated using the integrating-spherebased Hamamatsu C9920 quantum yield measurement system.
The resulting quantum yield of the upconverting nanoparticles
was thus given by
4794
2012
www.acsnano.org
197
XU ET AL.
VOL. 6
NO. 6
47884795
4795
2012
www.acsnano.org
198
ARTICLE
18. Mazhar, A.; Cuccia, D. J.; Gioux, S.; Durkin, A. J.; Frangioni,
J. V.; Tromberg, B. J. Structured Illumination Enhances
Resolution and Contrast in Thick Tissue Fluorescence
Imaging. J. Biomed. Opt. 2010, 15, 010506.
19. Menyuk, N.; Pierce, J.; Dwight, K. NaYF4:Yb,Er - Ecient
Upconversion Phosphor. Appl. Phys. Lett. 1972, 21, 159161.
20. Haase, M.; Schafer, H. Upconverting Nanoparticles. Angew.
Chem., Int. Ed. 2011, 50, 58085829.
21. Mader, H. S.; Kele, P.; Saleh, S. M.; Wolfbeis, O. S. Upconverting Luminescent Nanoparticles for Use in Bioconjugation
and Bioimaging. Curr. Opin. Chem. Biol. 2010, 14, 582596.
22. Ong, L. C.; Gnanasammandhan, M. K.; Nagarajan, S.; Zhang,
Y. Upconversion: Road to El Dorado of the Fluorescence
World. Luminescence 2010, 25, 290293.
23. Wang, F.; Banerjee, D.; Liu, Y. S.; Chen, X. Y.; Liu, X. G.
Upconversion Nanoparticles in Biological Labeling, Imaging, and Therapy. Analyst (Cambridge, U. K.) 2010, 135,
18391854.
24. Liu, H. C.; Xu, C. T.; Andersson-Engels, S. Multibeam
Fluorescence Diuse Optical Tomography Using Upconverting Nanoparticles. Opt. Lett. 2010, 35, 718720.
25. Svenmarker, P.; Xu, C. T.; Andersson-Engels, S. Use of
Nonlinear Upconverting Nanoparticles Provides Increased
Spatial Resolution in Fluorescence Diuse Imaging. Opt.
Lett. 2010, 35, 27892791.
26. Nyk, M.; Kumar, R.; Ohulchanskyy, T. Y.; Bergey, E. J.; Prasad,
P. N. High Contrast in Vitro and in Vivo Photoluminescence
Bioimaging Using Near Infrared to Near Infrared UpConversion in Tm3 and Yb3 Doped Fluoride Nanophosphors. Nano Lett 2008, 8, 38343838.
27. Bouzigues, C.; Gacoin, T.; Alexandrou, A. Biological Applications of Rare-Earth Based Nanoparticles. ACS Nano
2011, 5, 84888505.
28. Wang, F.; Deng, R. R.; Wang, J.; Wang, Q. X.; Han, Y.; Zhu,
H. M.; Chen, X. Y.; Liu, X. G. Tuning Upconversion through
Energy Migration in Core-Shell Nanoparticles. Nat. Mater.
2011, 10, 968973.
29. Cao, T. Y.; Yang, Y.; Gao, Y. A.; Zhou, J.; Li, Z. Q.; Li, F. Y. HighQuality Water-Soluble and Surface-Functionalized Upconversion Nanocrystals as Luminescent Probes for Bioimaging. Biomaterials 2011, 32, 29592968.
30. Bunzli, J. C. G. Lanthanide Luminescence for Biomedical
Analyses and Imaging. Chem. Rev. (Washington, DC, U. S.)
2010, 110, 27292755.
31. Auzel, F. Upconversion and Anti-Stokes Processes with
f and d Ions in Solids. Chem. Rev. (Washington, DC, U. S.)
2004, 104, 139173.
32. Yi, G.; Lu, H.; Zhao, S.; Yue, G.; Yang, W.; Chen, D.; Guo, L.
Synthesis, Characterization, and Biological Application of
Size-Controlled Nanocrystalline NaYF4:Yb,Er Infrared-toVisible Up-Conversion Phosphors. Nano Lett. 2004, 4,
21912196.
33. Lim, S. F.; Riehn, R.; Tung, C. K.; Ryu, W. S.; Zhuo, R.; Dalland,
J.; Austin, R. H. Upconverting Nanophosphors for Bioimaging. Nanotechnology 2009, 20, 405701.
34. Nam, S. H.; Bae, Y. M.; Il Park, Y.; Kim, J. H.; Kim, H. M.; Choi,
J. S.; Lee, K. T.; Hyeon, T.; Suh, Y. D. Long-Term Real-Time
Tracking of Lanthanide Ion Doped Upconverting Nanoparticles in Living Cells. Angew. Chem., Int. Ed. 2011, 50,
60936097.
35. Zhou, J.; Yu, M. X.; Sun, Y.; Zhang, X. Z.; Zhu, X. J.; Wu, Z. H.;
Wu, D. M.; Li, F. Y. Fluorine-18-Labeled Gd3/Yb3/Er3 CoDoped NaYF4 Nanophosphors for Multimodality PET/MR/
UCL Imaging. Biomaterials 2011, 32, 11481156.
36. Liu, Q.; Chen, M.; Sun, Y.; Chen, G. Y.; Yang, T. S.; Gao, Y.;
Zhang, X. Z.; Li, F. Y. Multifunctional Rare-Earth SelfAssembled Nanosystem for Tri-Modal Upconversion Luminescence/Fluorescence/Positron Emission Tomography Imaging. Biomaterials 2011, 32, 82438253.
37. Wang, J.; Wang, F.; Wang, C.; Liu, Z.; Liu, X. G. Single-Band
Upconversion Emission in Lanthanide-Doped KMnF3
Nanocrystals. Angew. Chem., Int. Ed. 2011, 50, 10369
10372.
38. Li, Z. Q.; Zhang, Y. An Ecient and User-Friendly Method
for the Synthesis of Hexagonal-Phase NaYF4:Yb, Er/Tm
Paper VIII
Multibeam fluorescence diffuse optical tomography
using upconverting nanoparticles
H. Liu, C. T. Xu, S. Andersson-Engels.
Optics Letters 35(5), 718-720 (2010).
Paper VIII
718
in FDOT [13]. These nanoparticles can emit antiStokes shifted light when excited at 980 nm [14],
which enables the signal to be detected in an
autofluorescence-free environment [15]. This leads to
a significant reduction of artifacts in the reconstructions. In addition, owing to the quadratic power dependence of the nanoparticles, the reconstructions
are more sharply defined compared with the reconstructions of a linear fluorophore [13].
In this Letter, we present an approach to exploit
the quadratic power dependence of upconverting
nanoparticles to gain additional information by utilizing two beams simultaneously for excitation in
FDOT. The effect of the images taken with dual-beam
excitation (named type-D images) on the reconstructions of the nanoparticle number density distribution, n, is demonstrated. In addition, comparisons of
reconstructed results between the linear rhodamine
6G and the quadratic upconverting nanoparticles are
made.
The excitation and emission fields can be modeled
by two coupled diffusion equations [13]. For quadratic fluorophores, the fluorescence signal detected
at a fixed detector position under excitation of the kth
beam, k, can be described by the forward model
N
k = Uf*rd,rinriUersk,ri2Vi ,
i=1
201
719
10
= 2 Uf*rd,rinriUersk,riUersj,riVi
i=1
+ k + j ,
with = 1 for linear fluorophores and = 2 for quadratic fluorophores. Calculations were performed using the NIRFAST package implementing the finiteelement method [16]. W was factorized according to
W = UV* ,
202
i=1
10
10
10
Singlebeam, = 1
Singlebeam, = 2
Singleanddualbeam, = 1
15
Singleanddualbeam, = 2
10
20
10
20
40
60
80
Singularvalue index
100
120
140
Paper VIII
720
well reconstructed at the center of the fluorescent lesion (represented by the circles) for both reconstructions, the reconstructed fluorescent lesion is more
confined for the case of using both type-S and type-D
images. This result confirms that the images of type
D indeed contribute to the inverse problem and lead
to better reconstructions for the quadratic upconverting nanoparticles. The corresponding reconstructions
for the linear rhodamine 6G were also carried out,
whose cross-sectional slices are presented in Fig. 4.
Compared with the results for the nanoparticles, the
reconstructions for rhodamine 6G do not benefit from
adding the type-D images, which is in agreement
with the theory. The true depth of the fluorescent lesion is also poorly reconstructed.
In summary, based on previous work regarding the
employment of upconverting nanoparticles in FDOT,
we propose and demonstrate an additional unique
advantage of the nonlinear power dependence of upconverting nanoparticles. This advantage enables the
possibility to obtain additional information for the inverse problem by using images taken with two or
more excitation beams simultaneously. We found that
this resulted in improved reconstructions. The same
advantage could not be found when using linear fluorophores, e.g., rhodamine 6G.
This work was supported by a Swedish Research
Council grant (VR 2007-4214) and a Linnaeus grant
for the Lund Laser Centre. The authors thank Prof.
Zhiguo Zhang and his group from the Harbin Institute of Technology and Dr. Gabriel Somesfalean for
their collaboration work.
Fig. 4. (Color online) Cross-sectional slices of the reconstructed relative rhodamine 6G distribution at x = 17 mm.
(a) Reconstruction with only type-S images. (b) Reconstruction with both type-S and type-D images.
References
1. A. Hielscher, Curr. Opin. Biotechnol. 16, 79 (2005).
2. A. Corlu, R. Choe, T. Durduran, M. A. Rosen, M.
Schweiger, S. R. Arridge, M. D. Schnall, and A. G.
Yodh, Opt. Express 15, 6696 (2007).
3. V. Ntziachristos, C.-H. Tung, C. Bremer, and R.
Weissleder, Nat. Med. 8, 757 (2002).
4. V. Ntziachristos, E. A. Schellenberger, J. Ripoll, D.
Yessayan, E. Graves, A. Bogdanov, L. Josephson, and
R. Weissleder, Proc. Natl. Acad. Sci. USA 101, 12294
(2004).
5. S. R. Arridge, Inverse Probl. 15, R41 (1999).
6. A. Soubret and V. Ntziachristos, Phys. Med. Biol. 51,
3983 (2006).
7. D. C. Comsa, T. J. Farrell, and M. S. Patterson, Phys.
Med. Biol. 53, 5797 (2008).
8. J. R. Mansfield, K. W. Gossage, C. C. Hoyt, and R. M.
Levenson, J. Biomed. Opt. 10, 041207 (2005).
9. G. Y. Panasyuk, Z.-M. Wang, J. C. Schotland, and V. A.
Markel, Opt. Lett. 33, 1744 (2008).
10. S. C. Davis, B. W. Pogue, R. Springett, C. Leussler, P.
Mazurkewitz, S. B. Tuttle, S. L. Gibbs-Strauss, S. S.
Jiang, H. Dehghani, and K. D. Paulsen, Rev. Sci.
Instrum. 79, 064302 (2008).
11. S. F. Lim, R. Riehn, W. S. Ryu, N. Khanarian, C.-K.
Tung, D. Tank, and R. H. Austin, Nano Lett. 6, 169
(2006).
12. C. Vinegoni, D. Razansky, S. A. Hilderbrand, F. Shao,
V. Ntziachristos, and R. Weissleder, Opt. Lett. 34, 2566
(2009).
13. C. T. Xu, J. Axelsson, and S. Andersson-Engels, Appl.
Phys. Lett. 94, 251107 (2009).
14. H. J. Liang, G. Y. Chen, L. Li, Y. Liu, F. Qin, and Z. G.
Zhang, Opt. Commun. 282, 3028 (2009).
15. C. T. Xu, N. Svensson, J. Axelsson, P. Svenmarker, G.
Somesfalean, G. Y. Chen, H. J. Liang, H. C. Liu, Z. G.
Zhang, and S. Andersson-Engels, Appl. Phys. Lett. 93,
171103 (2008).
16. H. Dehghani, M. E. Eames, P. K. Yalavarthy, S. C.
Davis, S. Srinivasan, C. M. Carpenter, B. W. Pogue,
and K. D. Paulsen, Commun. Numer. Methods Eng. 25,
711 (2009).
17. J. P. Culver, V. Ntziachristos, M. J. Holboke, and A. G.
Yodh, Opt. Lett. 26, 701 (2001).
18. E. Alerstam, S. Andersson-Engels, and T. Svensson, J.
Biomed. Opt. 13, 041304 (2008).
Fig. 3. (Color online) Cross-sectional slices of the reconstructed relative nanoparticle distribution at x = 17 mm. (a)
Reconstruction with only type-S images. (b) Reconstruction
with both type-S and type-D images.
203
Paper IX
Multispectral guided fluorescence diffuse optical
tomography using upconverting nanoparticles
P. Svenmarker, C. T. Xu, H. Liu, X. Wu, S. Andersson-Engels.
(2013) Accepted for publication by Applied Physics Letters.
Paper IX
Multispectral Guided Fluorescence Diffuse Optical Tomography using Upconverting Nanoparticles
Pontus Svenmarker,1, 2, 3, a) Can T. Xu,1 Haichun Liu,1 Xia Wu,1 and Stefan Andersson-Engels1
1)
Department of Physics, Lund University, P.O. Box 118, SE-221 00 Lund,
Sweden
2)
Department of Physics, Ume
a University, SE-901 87 Ume
a Sweden
3)
Centre for Mircobial Research (UCMR), Ume
a University, SE-901 87 Ume
a Sweden
(Dated: 7 October 2013)
We report improved image detectability for fluorescence diffuse optical tomography using upconverting
nanoparticles doped with rare-earth elements. Core-shell NaYF4 : Yb3+ /Er3+ @NaYF4 upconverting nanoparticels were synthesized through a stoichiometric method. The Yb3+ /Er3+ sensitizer-activator pair yielded
two anti-Stokes shifted fluorescence emission bands at 540 nm and 660 nm, here used to a priori estimate the
fluorescence source depth with sub-millimeter precision. A spatially varying regularization incorporated the
a priori fluorescence source depth estimation into the tomography reconstruction scheme. Numerically it is
fast, simple and direct to implement. Tissue phantom experiments showed both an improved resolution and
contrast in the reconstructed images as compared to not using any a priori information.
been proposed as an alternative.13 They offer a wider
excitation band, which may extend over a few hundreds
of nanometers, while being accompanied by large Stokes
shifted narrow band emission. Quantum dots, however,
typically suffers from limited imaging depth, since their
excitation band in the ultra-violet to the blue part of the
spectrum coincides with strong tissue attenuation. Upconverting nanoparticles, on the other hand, enable deep
tissue imaging due to a combination of NIR excitation
and large anti-Stokes shifted NIR and visible emission
bands.18
In this letter, we report improved image detectability
for fDOT using core-shell NaYF4 : Yb3+ /Er3+ @NaYF4
upconverting nanoparticles.
By the use of the
Yb3+ /Er3+ sensitizer-activator pair, two anti-Stokes
shifted emission bands (540 nm and 660 nm) upon 980
nm excitation, were available for fluorescence source
depth estimation. An advantageous approach of incorporating the a priori fluorescence source depth information
into the fDOT image reconstruction is through a spatially
varying regularization.15 To evaluate the performance
gained, a comparison between a multispectral guided spatially varying regularization and a non-spatially varying
Tikhonov regularization were made in a well controlled
tissue phantom experiment.
The method used for reconstructing the fluorescence distribution follows the procedure undertaken by
Zacharopoulos et al..12 In brief, it involves modeling the
photon migration as a set of diffusion equations expressed
as
a) Electronic
(1)
where Dx/f is the diffusion coefficient, a is the absorption coefficient at the excitation/fluorescence wavelength,
respectively; q(rq ) is the excitation source at position rq
and is the power-dependence exponent ( = 2 for the
upconverting nanoparticles used in this work ). A Robin
mail: [email protected]
207
z (mm)
(a)
17
4 mm
z=0
(b)
(c)
17 mm
(d)
10
0
-10
10
-10
10
-10
10
y (mm)
q,m (1 ) fm (r, 1 )
q,m (2 ) f (r, 2 ) < .
m
kJ(r) bk + kL(r)k
(3)
208
(4)
Paper IX
to be measured and the conjugate view (symmetrically
mirror the phantom at z = 0 mm, see Fig. 1a) could be
obtained by duplicating the data.
The first row in Fig. 3 displays the reconstructed images when using the multispectral regularization to guide
the reconstruction process. Each image shows an x-z
slice at y = 0 mm together with a cross-sectional plot.
For comparison, reconstructed images using Tikhonov
regularization are presented in the second row. Both
the resolution and the contrast were higher when using the multispectral regularization compared to the use
of Tikhonov regularization, leading to an increased detectability for the multispectral regularization guided reconstructions. The rods spaced 4 mm apart could be
resolved by the multispectral regularization guided reconstructions, while the Tikhonov regularization guided
reconstructions could only resolve the rods at an 8 mm
center-center distance. A comparison of the image contrast, defined as IC = (Imax Imin )/Imax , for the case
when the rods were placed at an 8 mm center-center distance gave 0.46 compared to 0.24, for the multispectral
regularization and the Tikhonov regularization, respectively.
For creation of a successful multispectral regularization, accurate and a large difference in optical properties
between the two wavelengths detected are crucial. In
a simplistic analysis, assuming an infinite homogeneous
geometry, the fluorescence ratio depth dependence can
be written as; log(1 /2 ) = log(D2 /D1 ) r(1eff 2eff ).
From this expression, it is clear that the difference in
effective attenuation coefficients between wavelengths 1
and 2 sets the depth sensitivity. A large difference, typically achieved by utilizing the sharp absorption profile
of haemoglobin, gives a high sensitivity and hence facilitates precise depth localization. An error in the deduced
optical properties will naturally lead to an inaccurate estimated depth.17 Any autofluorescence present will further misguide the depth localization.
In conclusion,
we reported improved image
detectability
for
fDOT
using
core-shell
NaYF4 : Yb3+ /Er3+ @NaYF4 upconverting nanoparticles, while guiding the reconstruction procedure with
a priori fluorescence source depth estimations, here
implemented through a spatially varying regularization.
Numerically it was fast, simple and direct to implement.
The use of NaYF4 : Yb3+ /Er3+ @NaYF4 upconverting
nanoparticles as contrast agent, with two emission bands
experiencing drastically different tissue attenuation,
yielded a multispectral regularization with precise depth
localization. Careful characterization of the tissue attenuation by time-of-flight and steady-state spectroscopy
resulted in an accurate depth localization. Autofluorescence insensitive measurements achieved through the use
of upconverting nanoparticles ensured contamination
free data, which otherwise may misguide the source
depth estimation and also the fDOT reconstruction.
The authors would like to acknowledge Erik Alerstam
for assisting with the TOFS measurements.
(a)
20
15
10
5
100rnm
0
10
(b)
s (cm-1)
eff (cm-1)
100
diameterr(nm)
1000
1
15
10
5
500 600 700 800 900
0.5
1
0
500
600
700
800
wavelengthr(nm)
900
fluorescencerintensityr(-)
numberr(d)
25
FIG.
2.
Characterization
of
the
core-shell
NaYF4 : Yb3+ /Er3+ @NaYF4 nanoparticles and optical
properties of the tissue phantom. (a) Nanoparticle size
distribution measured with DLS (diameter average equals
28 nm) together with an inset containing a corresponding
TEM image with a median particle diameter of 22 nm. (b)
A Fluorescence emission spectrum of the nanoparticles upon
975 nm excitation is given by the dotted line. The solid
line indicates the effective attenuation of the liquid phantom
measured with steady-state spectroscopy and the red crosses
give the same quantity measured by TOFS. The inset show
the reduced scattering coefficient obtained by TOFS together
with an extrapolation.
209
4 mm
17
5 mm
6 mm
7 mm
8 mm
IC = 0.10
IC = 0.19
IC = 0.37
IC = 0.37
IC = 0.46
IC = 0
IC = 0
IC = 0
IC = 0
IC = 0.24
z (mm)
10
0
17
10
0
-10
10
-10
10
-10
0
x (mm)
10
-10
10
-10
10
FIG. 3. Reconstructed fluorescence diffuse optical tomography images for different center-center distances of the fluorescent
contrast agent. All images show a x-z slice at y=0 mm together with a cross-section plot along the z-direction. The first row
displays multispectral regularization guided reconstructions and the second row Tikhonov regularization guided reconstructions.
REFERENCES
1
210
Paper IX
sub-surface imaging. Physics in Medicine and Biology,
56(21):6823, 2011.
J. Svensson and S. Andersson-Engels. Modeling of
spectral changes for depth localization of fluorescent inclusion. Optics Express, 13(11):42634274, May 2005.
18
Qiuqiang Zhan, Jun Qian, Huijuan Liang, Gabriel
Somesfalean, Dan Wang, Sailing He, Zhiguo Zhang,
and Stefan Andersson-Engels. Using 915 nm laser excited Tm3+ /Er3+ /Ho3+ -doped NaYbF4 upconversion
17
211