Research Article  Open Access
A New 3D Simulation Method for the Construction of Optical Phase Contrast Images of Gold Nanoparticle Clusters in Biological Cells
Abstract
A new 3D simulation method based on the finitedifference time domain (FDTD) approach in combination with Fourier optic techniques is applied to the modeling of optical phase contrast microscope (OPCM) imaging of gold nanoparticles (NPs) in singe biological cells. We consider a realistic size 3D cell model at optical immersion conditions, that is, when the refractive index values of the cytoplasm and of the extracellular medium are equal. For the first time, an FDTDbased OPCM model is applied to visualize the presence of a cluster of gold NPs in the cytoplasm at both resonant and nonresonant conditions. The results demonstrate the capability to model OPCM image enhancement by optically controlling the resonant properties of the NPs. Our research study extends the applicability of the FDTD modeling approach into a new biomedical optics research area.
1. Introduction
One of the growing areas in biomedical optics and photonics [1, 2] is the application of optical software simulation and modeling tools to provide deeper understanding of newly developed optical diagnostics and imaging techniques. There are a number of numerical modeling approaches that can be used for the modeling of the light scattering from biological cells [3] including the finitedifference time domain (FDTD) method [4–9]. Recently, we demonstrated the applicability of the 2D FDTD approach to numerically study the changes in the forward scattered light phase and intensity distributions due to different cell membrane thicknesses and different configurations of gold NPs in the cytoplasm of small size cells using the optical immersion technique (OIT) [7, 8]. The OIT is associated with the socalled “optical clearing” (OC) effect consisting in the increased light transmission through cells or tissues due to the matching of the refractive indices of some of their morphological components to that of the extracellular or interstitial media [10–12]. In this article, we present the first results of a research study applying a 3D FDTDbased optical phase contrast microscope (OPCM) imaging approach to the modeling of a cluster of gold nanoparticles (NPs) in realistic size single cells. We demonstrate a new simple way of using the FDTD approach for the construction of OPCM images by moving its modeling power from the calculation of phase distribution curves [7] to realistic cell images. We validate our approach by modeling a gold nanoparticle cluster in the cell cytoplasm to provide a valuable design methodology for further studies in optical nanobioimaging and nanotherapeutics.
2. The FDTD Approach
The FDTD technique is a numerical solution of Maxwell’s equations [13]. In most light scattering simulation experiments, one usually uses the socalled total field/scattered field (TFSF) formulation [8, 13] to excite the magnetic and electric fields and simulate a linearly polarized plane wave propagating in a finite region containing the simulation equivalent of the object under investigation [13]. A schematic representation of the 2D FDTD computational domain is shown in Figure 1(a). TFSF sources are used to separate the FDTD computation domain into two distinct regions. The first one is internal to the computational domain and contains the total near fields and , that is, the sum of the incident field and the scattered near field. The second one occupies the rest of the computational domain and contains only the scattered near fields and . The final FDTD simulation results require additional postprocessing computational procedures to transform the real values of the calculated near fields into their complex far field counterparts. The far fields include the transverse distribution of the forward scattered light from the cell which is of high relevance for the modeling of optical phase contrast imaging. At the edges, the entire FDTD computational domain is truncated by the socalled perfectly matched layer (PML) boundary conditions [13].
(a)
(b)
The 3D case of the TFSF formulation described above follows a similar logic. The 3D simulation results provided here are based on a modified 3D TFSF formulation that could be more appropriately called total field/reflected field (TFRF). The 3D “TFRF” formulation uses a TFSF region which contains the biological cell and extends beyond the limits of the simulation domain (Figure 1(b)). The extension of the transverse dimension of the input field beyond the limits of the computational domain through the PML boundaries would lead to distortions of its ideal plane wave shape and eventually distort the simulation results. To avoid these distortions, one must use Bloch periodic boundary conditions (Figure 1(b)) in the lateral x and ydirections which are perpendicular to the direction of propagationz [13]. Bloch boundary conditions are periodic boundary conditions which take into account the phase effects due to the tilting of the input plane waves incoming at periodic structures, that is, what we are actually modeling is a periodic row of biological cells. The near scattered fields, however, are calculated in the transverse planes located in the close proximity to the cell where the coupling effect due to waves scattered from adjacent cells is negligible. This effect can be further minimized or completely removed by controlling the lateral dimension of computational domain by using a large enough period of the periodic cell structure. The larger is this period, the smaller is the coupling effect. In the 3D TFRF formulation, the location in the computational domain corresponding to the forward scattered light is positioned within the total field region (Figure 1(b)). The OPCM simulation model requires the explicit availability of the forward scattered transverse distribution of the fields. The phase of the scattered field accumulated by a plane wave propagating through a biological cell will be used in the FDTD model of the OPCM that will be described in Section 3.
3. Optical Phase Contrast Microscope Model
3.1. FDTD OPCM Principle
Phase contrast microscopy is utilized to produce highcontrast images of transparent specimens such as microorganisms, thin tissue slices, living cells, and subcellular components such as nuclei and organelles. It translates smallphase variations into corresponding changes in amplitude visualized as differences in image contrast. A standard phase contrast microscope design is shown in Figure 2(a), where an image with a strong contrast ratio is created by coherently interfering a reference (R) with a diffracted beam (D) from the specimen (http://www.microscopyu.com/articles/phasecontrast/phasemicroscopy.html). Relative to the reference beam, the diffracted beam has lower amplitude and is retarded in phase by approximately through interaction with the specimen.
(a)
(b)
The main concept underlying the design of the phase contrast microscope is the spatial separation of the R beam from D wave front emerging from the specimen. In addition, the amplitude of the R beam light must be reduced and the phase should be advanced or retarded by another in order to maximize the differences in the intensity between the specimen and the background in the image plane. The mechanism for generating relative phase retardation has two steps: (i) the D beam is being retarded in phase by a quarter wavelength (i.e., ) at the specimen, and (ii) the R beam is advanced (or retarded) in phase by a phase plate positioned in or very near the objective rear focal plane (Figure 2(a)). This twostep process is enabled by a specially designed annular diaphragm—the condenser annulus—which is placed in the condenser front focal plane, is matched in diameter, and optically conjugates to the phase plate residing in the objective rear focal plane. The resulting image, where the total phase difference is translated by interference at the image plane into an amplitude variation, can have a high contrast ratio, particularly if both beams have the same amplitude.
3.2. Excitation and Simulation of the FDTD near Fields
Figure 2(a) illustrates the part of the microscope that will become the subject of FDTD modeling combined with Fourier optics. Figure 2(b) provides a visual representation illustrating the major steps in the FDTD OPCM model. The phase contrast microscope uses incoherent annular illumination that could be approximately modeled by adding up the results of eight different simulations using ideal input plane waves incident at a given polar angle (30 degrees), an azimuthal angle (0, 90, 180, or 270 degrees), and a specific light polarization (parallel or perpendicular to the plane of the graph). Every single FDTD simulation provides the near field components in a transverse monitoring plane located right behind the cell (Figure 1(b)).
3.3. FarField Transformation
The far field transformations (Figure 3) use the FDTD calculated near fields right behind the cell and return the three complex components of the electromagnetic fields at 1 meter distance (i.e., long enough) from the location of the near fields, that is, in the far field [13, 14]: and , where refer to a spherical coordinate system and the variables and are the x and y components of the unit direction vector u. The unit direction vector is related to the angular variables by with . The inplane wave vectors for each plane wave are given by and , where . To scale the far field to a different distance, we can recognize that in 3D, the electric field intensity decreases like : , where is the field at 1 m distance resulting from the far field transformation and R is measured in meters. The specific details of the far field transformation can be found elsewhere [13].
It is important to note that the neartofar field projection must take into account the Bloch periodic boundary conditions in the lateral dimension and is calculated only at angles that correspond to diffracted orders of the periodic structure as defined by the Bragg conditions. This is done by calculating the direction cosines of all the diffracted orders that meet the Bragg condition, and interpolating the previously farfield distributions onto those specific directions. The neartofar field projection, therefore, provides the electric field amplitude and phase corresponding to each diffracted order. The zeroth order, that is, the light travelling through the scattering object without any deviation in angle, is the reference beam; and the phase contrast microscope is designed to provide a phase delay to this order.
3.4. Optical Magnification
Now that we have the field components , we can use their amplitudes and phases to determine all the properties of the polarization as well as do Fourier optics for both the scattered and reference beams. We assume an ideal optical lens system that could be characterized by a given magnification factor. This simple model could be easily extended to include the numerical equivalent of the two lenses (Figure 2) together with an additional model to take into account any aberrational effects. The magnification was implemented by modifying the angle of light propagation, that is, by multiplying the direction cosines and by the inverse value of the desired magnification factor : and . In any other circumstances, the modification of the direction cosines would lead to complications because of the vectorial nature of the field. In our case, however, working in spherical coordinates ( and ) leads to the advantage that the vectorial components do not change when and are modified because they are part of a local coordinate system that is tied to the values of and . The magnification factor was applied to the far fields before the interference of the diffracted (D) and reference (R) beams at the image plane.
3.5. Numerical Aperture
It was also possible to apply the effect of a numerical aperture which clips any light that has too steep angle and would not be collected by the lens system. This means that all the light with is being clipped. The effect of the aperture is defined by the inequality applied to the corrected aperture angles and : where and the “sqrt” labels a square root mathematical operation. The magnified field components will then have the following form.
Diffracted Beam:
Reference Beam:
where are are the far field components of the reference
beam.
3.6. OPCM Image Construction
The fields given above are then used to calculate back the Fourier inverse transform of the far field transformed fields leading to the distribution of the scattered and the references beams in the image plane.
Diffracted Beam:
Reference Beam:
where the summation is over all angles.
The OPCM images at the image plane are calculated by adding up the scattered and the reference beam at any desired phase offset : The coefficient and the phase correspond to the ability of the OPCM to adjust the relative amplitudes and the phase difference between the two beams. In all simulation results presented in this article, and is used as parameter.
4. Cell Structure, Gold Nanoparticle Model and Simulation Results
4.1. Cell Structure
The cell is modeled as a dielectric sphere with a realistic radius . It is critically important to pay special attention to the membrane model since assigning the membrane the precisely desired thickness value leads to inadequate modeling of optical immersion experiments. This is due to the fact that the FDTD technique averages out the refractive index values at the interface of two different materials effectively reducing the membrane thickness value and eventually, due to the staircase approximation error, destroying the continuity of the membrane. This is the reason for us to select a membrane thickness nm which corresponds to effective thickness of approximately 10 nm. The cell nucleus is also spherical with a radius centered at a position which is 2.0 shifted from the cell center in a direction perpendicular to the direction of light propagation. The refractive index of the cytoplasm is , of the nucleus , of the membrane , and of the extracellular material (no refractive index matching) or 1.36 (refractive index matching).
4.2. Optical Properties of Gold Nanoparticles
Gold NPs exhibit theability to resonantly scatter visible and near infrared light due to the excitation of surface plasmonresonances (SPRs). The scattering ability is extremely sensitive to their size, shape, andaggregation state offering a great potential for optical cellular imaging and detection labeling studies [15–19]. Our FDTD approach uses the dispersion model for gold derived from the experimental data provided by Johnson and Christy [7, 20] where the total, complexvalued permittivity is given as Each of the three contributions to the permittivity arises from a different material model. The first term represents the contribution due to the basic, background permittivity. The second and third terms represent Lorentz and plasma contributions: where all material constants are summarized in Table 1.

We have modeled two different cases—resonant and nonresonant. The ability to model both these cases together with the effect of optical clearing effect provides the opportunity to numerically study the possibility for imaging of the uptake of clusters of NPs—a scenario which, to the best of our knowledge, was not studied before. We used the FDTD technique to calculate the scattering and absorption crosssections over a 400–900 nm wavelength range for a single 50 nm diameter gold NP immersed in a material having the properties of the cytoplasm () and resolution nm. The scattering crosssection is defined as where is the total scattered power and is the intensity of the incident light in (W/m^{2}). It was calculated by applying the total field/scattered field FDTD formulation [13] and by creating 12 field power monitoring planes around the nanoparticle in the form of a box: 6 in the total field region and 6 in the scattered field region. The total scattered power is calculated by summing up the power flowing outward through 6 scattered field power monitors located in the scattered field region.
The absorption crosssection is similarly defined as where is the total power absorbed by the particle. The power absorbed by the particle is calculated by calculating the net power flowing inward through the 6 total field power monitors located in the total field region. The extinction crosssection is the sum of the absorption and scattering crosssections Figure 4 shows that the extinction crosssection has a maximum of 3.89 at around 543.0 nm corresponding to one of the radiation wavelengths of HeNe lasers. In this article, we also present results for nm (a Krypton laser wavelength) which corresponds to the nonresonant case (extinction crosssection value 0.322, ~12 times smaller than 3.89). The FDTD results are compared with the theoretical curve calculated by Mie theory. The main reason for the slight discrepancy between the theoretical and FDTD results for the extinction crosssection is the finite mesh size.
To consistency of the results could be improved by reducing the mesh size. In the FDTD simulations including both the cell and the nanoparticles, the numerical resolution of the nanoparticles was hardcoded to nm to make sure that their numerically manifested optical resonant properties will be the same as the ones shown in Figure 4. However, Figure 5 visualizes the positioning of a cluster of 42 nanoparticles in the cytoplasm that was used to produce the FDTD simulation results presented in Section 5. It was randomly generated without controlling the possibility for spatial overlaps or contacts between the nanoparticles. The closer visual inspection showed that there were no such overlaps and contacts.
4.3. OPCM Image Simulation Results
The simulation results presented here are based on the 3D FDTD OPCM model. The simulations use nonuniform meshing where the number of mesh points in space is automatically calculated to ensure a higher number of mesh points in materials with higher values of the refractive index [14]. The realistic cell dimensions (including both cell radius and membrane) require a very fine numerical resolution making the simulations computationally intensive. The locally refined resolution of the nanoparticles leads to additional requirements for the CPU time and memory (~120 Gbs RAM) requiring high performance computing resources. The time step used during the simulation was defined by means of a factor of 0.99 in the Courant stability limit: . The cell center is located in the middle () of the computational domain with dimensions . The nucleus’ center is located at . The cluster of gold nanoparticles is located at (Figure 5).
The refractive index matching (RIM) between the cytoplasm and the extracellular medium leads to the optical clearing of the image and a significantly better OPCM image of the nucleus and of the membrane—the only sources of phase contrast at optical immersion conditions. Based on the fact that RIM enhances significantly the imaging of the cell components, we have used the FDTD OPCM model to create the OPCM images of the same cell at optical immersion conditions including the cluster of 42 gold NPs (Figure 5, see also Figure 1(b)) and for different phase offsets between the reference beam and the scattered beam (assuming ). The insertion of the nanoparticles increases significantly the CPU and memory requirements. The dimensions of the computational domain are . The OPCM images in Figure 6 show that the gold NP cluster can be clearly identified at both the nonresonant ( nm) and resonant ( nm) wavelengths.
(a)
(b)
(c)
(d)
The graphs in Figure 7 compare the geometrical crosssections of the three OPCM images shown at the bottom of Figure 6 at showing the relevant half of the images (right) where the gold NP cluster is located. At resonance, the optical contrast of the gold NP peak is ~2.24 times larger than the one at no resonance and 24.79 times larger than the background optical contrast corresponding to the case when there are no nanoparticles. The enhanced imaging of the gold NP cluster at resonant conditions is of no surprise. It, however, needs to be further studied as a function of particular phase delay between the reference beam and the scattered beam. Figure 8 visualizes the optical contrast due to the gold nanoparticle cluster as a function of the phase offsets between the reference (R) and the diffracted (D) beams of the optical phase microscope. It shows that the enhancement of the optical contrast due to the nanoparticle resonance changes significantly from a minimum of 0.0 () to a maximum of 3.60 ().
(a)
(b)
The correct modeling of the gold NPs required making sure that the FDTD resolution is the same as in the calculation of their optical scattering/absorption properties, that is, nm. This was done by using the builtin features of our FDTD simulation software [14] enabling the possibility for using both nonuniform meshing in the entire computational domain and local hard coding of the resolution of the NPs. The results presented in Figures 6–8 clearly demonstrate the capability of the proposed FDTD OPCM approach to model the visual effect of gold NPs at nonresonant and resonant conditions. A future extension of this research will be to study the capability of the model to provide valuable insights for the application of gold NPs in optical nanotherapeutics.
5. Conclusions
In this article, we provided simulation results on the application of the FDTD approach to the modeling of the OPCM imaging of gold NPs in singe biological cells. We have first reproduced the effect of optical immersion on the OPCM images of a realistic size cell containing a cytoplasm, a nucleus, and a membrane. The model was then extended to include the presence of a cluster of gold NPs in the cytoplasm at optical immersion conditions as well as the enhancing imaging effect of the optical resonance of the nanoparticles. The results do not allow analyzing the scaling of the NP imaging effect as a function of the number of the NPs in the cluster. However, the validation of the model provides a basis for future studies of OPCM nanobioimaging including the effects of NP cluster size, NP size and number, as well as average dimension between the NPs. To the best of our knowledge, this is the first time that FDTD simulation results were postprocessed in a way providing realistic OPCM images of biological cells including a cluster of gold NPs. The results demonstrate the potential of the FDTD approach to model advanced optical nanobioimaging instruments.
Acknowledgments
S. Tanev, J. Pond, and P. Paddon acknowledge the use of the computing resources of WestGrid (Western Canada Research Grid)—a $50million project to operate gridenabled high performance computing and collaboration infrastructure at institutions across Canada. V. V. Tuchin was supported by grants of Federal Agency of Education of RF no. 1.4.06, RNP.2.1.1.4473, RFBR no. 060216740, and Photonics4Life of FP7. The authors are thankful to Dr. Vladimir P. Zharov for the initial discussion on the relevance of the problem of nanoparticle clustering in single cells.
References
 P. N. Prasad, Introduction to Biophotonics, chapter 7, John Wiley & Sons, New York, NY, USA, 2003.
 V. V. Tuchin, Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis, vol. 166 of Press Monograph Series, SPIE Press, Bellingham, Wash, USA, 2nd edition, 2007.
 F. M. Kahnert, “Numerical methods in electromagnetic scattering theory,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 7980, pp. 775–824, 2003. View at: Publisher Site  Google Scholar
 R. Drezek, A. Dunn, and R. RichardsKortum, “A pulsed finitedifference timedomain (FDTD) method for calculating light scattering from biological cells over broad wavelength ranges,” Optics Express, vol. 6, no. 7, pp. 147–157, 2000. View at: Google Scholar
 T. Tanifuji and M. Hijikata, “Finite difference time domain (FDTD) analysis of optical pulse responses in biological tissues for spectroscopic diffused optical tomography,” IEEE Transactions on Medical Imaging, vol. 21, no. 2, pp. 181–184, 2002. View at: Publisher Site  Google Scholar
 R. Drezek, M. Guillaud, T. Collier et al., “Light scattering from cervical cells throughout neoplastic progression: influence of nuclear morphology, DNA content, and chromatin texture,” Journal of Biomedical Optics, vol. 8, no. 1, pp. 7–16, 2003. View at: Publisher Site  Google Scholar
 S. Tanev, V. V. Tuchin, and P. Paddon, “Light scattering effects of gold nanoparticles in cells: FDTD modeling,” Laser Physics Letters, vol. 3, no. 12, pp. 594–598, 2006. View at: Publisher Site  Google Scholar
 S. Tanev, V. V. Tuchin, and P. Paddon, “Cell membrane and gold nanoparticles effects on optical immersion experiments with noncancerous and cancerous cells: finitedifference timedomain modelling,” Journal of Biomedical Optics, vol. 11, no. 6, Article ID 064037, 6 pages, 2006. View at: Publisher Site  Google Scholar
 X. Li, A. Taflove, and V. Backman, “Recent progress in exact and reducedorder modeling of lightscattering properties of complex structures,” IEEE Journal on Selected Topics in Quantum Electronics, vol. 11, no. 4, pp. 759–765, 2005. View at: Publisher Site  Google Scholar
 R. Barer, K. F. A. Ross, and S. Tkaczyk, “Refractometry of living cells,” Nature, vol. 171, no. 4356, pp. 720–724, 1953. View at: Publisher Site  Google Scholar
 B. A. Fikhman, “Microbiological refractometry,” Medicine, pp. 81–94, 1967 (Russian). View at: Google Scholar
 V. V. Tuchin, Optical Clearing of Tissues and Blood, vol. 154 of Press Monograph Series, SPIE Press, Bellingham, Wash, USA, 2005.
 A. Taflove and S. Hagness, Computational Electrodynamics: The FiniteDifference Time Domain Method, Artech House, Boston, Mass, USA, 3rd edition, 2005.
 The simulations were performed by the FDTD ${\text{Solutions}}^{\text{TM}}$ software developed by Lumerical Solutions Inc., Vancouver, BC, Canada, http://www.lumerical.com.
 K. Sokolov, M. Follen, J. Aaron et al., “Realtime vital optical imaging of precancer using antiepidermal growth factor receptor antibodies conjugated to gold nanoparticles,” Cancer Research, vol. 63, no. 9, pp. 1999–2004, 2003. View at: Google Scholar
 I. H. ElSayed, X. Huang, and M. A. ElSayed, “Surface plasmon resonance scattering and absorption of antiEGFR antibody conjugated gold nanoparticles in cancer diagnostics: applications in oral cancer,” Nano Letters, vol. 5, no. 5, pp. 829–834, 2005. View at: Publisher Site  Google Scholar
 N. G. Khlebtsov, A. G. Melnikov, L. A. Dykman, and V. A. Bogatyrev, “Optical properties and biomedical applications of nanostructures based on gold and silver bioconjugates,” in Photopolarimetry in Remote Sensing, G. Videen, Ya. S. Yatskiv, and M. I. Mishchenko, Eds., vol. 161 of NATO Science Series II: Mathematics, Physics and Chemistry, pp. 265–308, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2004. View at: Publisher Site  Google Scholar
 V. P. Zharov, J.W. Kim, D. T. Curiel, and M. Everts, “Selfassembling nanoclusters in living systems: application for integrated photothermal nanodiagnostics and nanotherapy,” Nanomedicine: Nanotechnology, Biology and Medicine, vol. 1, no. 4, pp. 326–345, 2005. View at: Publisher Site  Google Scholar
 V. P. Zharov, K. E. Mercer, E. N. Galitovskaya, and M. S. Smeltzer, “Photothermal nanotherapeutics and nanodiagnostics for selective killing of bacteria targeted with gold nanoparticles,” Biophysical Journal, vol. 90, no. 2, pp. 619–627, 2006. View at: Publisher Site  Google Scholar
 P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Physical Review B, vol. 6, no. 12, pp. 4370–4379, 1972. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2008 Stoyan Tanev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.