Role of Arbitrary Intensity Profile Laser Beam in Trapping of RBC for Phase-imaging
Role of Arbitrary Intensity Profile Laser Beam in Trapping of RBC for Phase-imaging
Journal of the Optical Society of Korea. 2016. Feb, 20(1): 78-87
Copyright © 2016, Optical Society of Korea
  • Received : November 03, 2015
  • Accepted : December 28, 2015
  • Published : February 25, 2016
Export by style
Cited by
About the Authors
Ranjeet Kumar
Department of Physics, HINDU College, University of Delhi, Delhi 110007, Indian
Vishal Srivastava
Department of Electrical and Instrument Engineering, Thapar University, Patiala 147004, India
Dalip Singh Mehta
Department of Physics, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
Chandra Shakher
Laser Applications and Holography Laboratory, Instrument Design Development Centre, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India

Red blood cells (RBCs) are customarily adhered to a bio-functionalised substrate to make them stationary in interferometric phase-imaging modalities. This can make them susceptible to receive alterations in innate morphology due to their own weight. Optical tweezers (OTs) often driven by Gaussian profile of a laser beam is an alternative modality to overcome contact-induced perturbation but at the same time a steeply focused laser beam might cause photo-damage. In order to address both the photo-damage and substrate adherence induced perturbations, we were motivated to stabilize the RBC in OTs by utilizing a laser beam of ‘ arbitrary intensity profile ’ generated by a source having cavity imperfections per se . Thus the immobilized RBC was investigated for phase-imaging with sinusoidal interferograms generated by a compact and robust Michelson interferometer which was designed from a cubic beam splitter having one surface coated with reflective material and another adjacent coplanar surface aligned against a mirror. Reflected interferograms from bilayers membrane of a trapped RBC were recorded and analyzed. Our phase-imaging set-up is limited to work in reflection configuration only because of the availability of an upright microscope. Due to RBC’s membrane being poorly reflective for visible wavelengths, quantitative information in the signal is weak and therefore, the quality of experimental results is limited in comparison to results obtained in transmission mode by various holographic techniques reported elsewhere.
Optical imaging is a burgeoning area whereby various non-invasive techniques are employed explicitly to retrieve quantitative morphological (volume, surface area and sphericity), biochemical (haemoglobin content and concentration), and mechanical properties (dynamic membrane fluctuations) of individual bio-specimens specially RBCs which are key to diagnose the pathophysiology of a number of diseases in hematology [1 - 4] . Automated blood cell counters based on the complete blood count have been most commonly used to analyze the primary and critical properties of RBCs in clinical hematology but limitations exist in accurate diagnosis [5] . Dynamism of RBCs and their related pathophysiology can be explored using detailed analysis at the single cell level. In this context, non-invasive and highly sensitive optical diagnostics modalities such as; optical tweezers (OTs) [6 - 8] , flow-cytometry [9] , atomic force microscopy and micro-Raman spectroscopy [10] have been employed. The intrinsic optical parameters of the human RBC, such as absorption coefficient, scattering coefficient, anisotropy factor and hence complex refractive index (RI) were determined both experimentally and theoretically by Friebel et al . [11] . The other optical techniques at the individual cell level are; interferometric quantitative phase-microscopy [12 , 13] , digital holographic phase-contrast and differential interference contrast microscopy [14 - 18] , lateral shearing digital holographic imaging [19] , Fourier-phase microscopy [20] , Hilbert phase-microscopy (HPM) [21 , 22] , diffraction phase microscopy and tomographic phase-microscopy [3 , 23 - 26] . All have enabled the extraction of sub-nanometer path-length changes, 3D RI map and millisecond scale fluctuation dynamics of RBC membrane to enhance the understanding of RBCs significantly [27] . Further Rappaz et al . [28 , 29] and Kemper et al . [30 - 33] in combination with microfluidics have demonstrated digital holographic microscopy (DHM) for quantitative phase imaging (QPI) of live-cells from a single-shot hologram followed by decoupling procedure conceptualized by Evans and others [34 , 35] . In the realm of technological improvements a simple, cost-effective and compact optical unit based on a common-path phase shifting interferometry for QPI has been reported [36 , 37] . Recently, real-time imaging and topography of RBCs from an entire complex optical wavefront (amplitude and phase) has been reported both in reflection and transmission mode using digital holographic reflectometry (DHR) [38] . But in DHR, RBCs are made stationary by adhering them directly to a glass surface that can invoke perturbations from two facets - physical contact and their own weight which might alter their innate morphology. Digital holography in conjugation with OTs (HOT) has been reported for 3D shape evaluation of convex shape objects by observing cells along different directions [39 - 42] . Mohanty et al . [43] have combined HOTs with quantitative phase-microscopy for decoupling geometric thickness from RI. It was realized by re-orienting the trapped RBCs by altering the power of beam and then retrieving phase-imaging at various orientations in transmission mode. But no clue of the consequences of radiation stress on the shape memory of anisotropic RBC on account of altering beam power is found. Most recently, Park et al . have reported a system employing optical diffraction tomography in conjugation with HOTs capable of simultaneous 3D visualization of the shapes and tracking positions of trapped microscopic samples with high lateral and axial resolution [44 , 45] . Other common-path setups: a Michelson interferometer and a compact, portable, and extremely stable τ-interferometer have been reported which have enabled the recording of amplitude and phase profiles of biological samples for quantitative phase microscopy [17 , 36 , 37 , 46 , 47] . But a reflective surface of cover glass above the sample was used to generate a separate reference wave to create interference [48] which was avoided in DHM using a Michelson interferometer [17] . Memmolo et al . [49] provided an outstanding modality for full 3D characterization of both concave and convex shape cells using OTs with fiber Y-coupler. Here, it is pertinent to caution that our phaseimaging outcomes should not be compared with those already reported by Memmolo et al . [49] . Because we want to investigate first the utility of ‘ arbitrary intensity profile ’ containing inhomogeneous transverse intensity distribution for photodamage-free trapping of RBC and then interferometric phase imaging of trapped RBC using probe interferograms generated by a lab-made miniaturized Michelson interferometer (MMI). Thus immobilized RBC was obviated from substrate adherence induced perturbations. The experiment could be carried out only in reflection configuration due to the availability of an upright microscope.
The aim of this paper is two-fold. First to establish the usability of a laser beam of ‘ arbitrary intensity profile ’ in photodamage-free trapping of delicate RBCs so that contact activated perturbations, as arising in DHR [38] during phase-microscopy, can be addressed. Secondly, phase imaging of thus trapped RBC by interferograms generated from a lab-made MMI. The novelty of this study is to prove firmly the re-use of defective lasers/low-cost laser sources emitting beam profiles other than Gaussian profile for making the specimen stationary whereby beam quality factor does not carry a significant role. In addition, developing an inexpensive OTs set-up for multi-particle trapping.
In general, the 2D intensity distribution of an interferogram can be expressed as [21] :
PPT Slide
Lager Image
where a ( x , y ) is the background intensity, b(x, y) is the amplitude modulation of the fringe pattern, Δ d ( x , y ) denotes the optical path-difference (OPD) between sample and reference arms, θ (x, y) is the tilt angle and k = 2 π /λ is the wavenumber ( λ being the wavelength of laser light). The last term in Eq. (1) is the interference sinusoidal function of (2 k Δ d ( x , y )+ θ (x, y)). To extract phase, the spectral interferogram is Fourier transformed then spatially high-pass filtered to isolate the sinusoidal term and eliminate the background components. The filtered sinusoidal term can be expressed as;
PPT Slide
Lager Image
Complex analytic signal associated with the real function of S ( x , y ) can be obtained as;
PPT Slide
Lager Image
In Eq. (3) the imaginary part of the right-hand side stands for a principal-value integral, easily identifiable as the Hilbert transform of c ( x , y ). HPM extends the concept of complex analytic signals to the spatial domain and measures quantitative phase images from single-shot spatial interferogram recording [22 , 33] . With the help of the Hilbert transform we obtained the complex analytical signal. Using Eq. (3), the wrapped phase map ϕ ( x , y ) associated with complex analytic signal S ( x , y ) is calculated as:
PPT Slide
Lager Image
However, the phase profile Δ φ ( x , y ) can also be retrieved with the following equation
PPT Slide
Lager Image
where Δ φ ( x , y ) is the phase-difference between sample ( φ s ( x , y )) and reference ( φR ( x , y )), ns ( x , y ) is the RI of sample (RBC), n 0 ( x , y ) is the RI of suspension medium which is taken constant at illumination wavelength and d is the thickness of discotic RBC (sample) at ith pixel that varies along the diameter. Henceforth a 3D topographic phase-map of the RBC can be computed from Eq. (5). From Eq. (4) both the real and imaginary parts of S ( x , y ) depend on c ( x , y ) which is related with ϕ ( x , y ) by Eq. (5) that is further entangled with optical thickness (Δ d ) of the RBC. Interdependency of phase-map extraction with optical thickness and hence computation of RI of the RBC is customarily done by using two suspension media of known RI. Procedures to decouple these entangled parameters is elucidated elsewhere [31 , 34 , 35 , 39 - 41 , 43] .
The schematic of simultaneous trapping of RBC using a laser beam of ‘ arbitrary intensity profile ’ and phase-microscopy using MMI is shown in Fig. 1 . The single-beam optical trap used in this study is calibrated using a fan-cooled diode-pumped solid state (DPSS) Nd 3+ :YVO 4 laser source (LASERGLOW from Richmond Hill, Canada) emitting an arbitrary beam profile of wavelength at 532 nm. It is imperative to mention here that cavity imperfections caused the generation of arbitrary beam per se . As far as the reproducibility of the arbitrary beam from the DPSS laser source is concerned, the spatial intensity distribution remains time-invariant. Various cavity-imperfections of this laser source such as; misalignment of cavity mirrors, age-related degradation of cavity material, non-homogeneity of active medium or lens aberrations could be inferred to be causing the generation of such beam of ‘ arbitrary intensity profile ’. To ensure technically the root cause of generation, the laser system was unbolted as shown in Fig. 2 and misalignment of mirrors was found clearly. In addition, age-related degradation of cavity material cannot be ruled out because measured power at the exit of the source was just 40 mW whereas company specifications pegged it to 1W. The role of the last two reasons mentioned above could not be ascertained but may also have participated. It is necessary to emphasize here that a broken laser was not chosen onpurposely nor motivated technical modifications were introduced within the cavity to generate such a beam profile of ‘ arbitrary intensity profile ’. But a careful monitoring of beam characteristics that was incidentally like this forced us to unbolt the laser for systematic investigation of its root cause whereby it was found that wavelength of the pumping-diode (@ 808 nm) was being converted into ideal narrow-band output via the host crystals Nd 3+ : YVO 4 . Further, utilization of ‘ arbitrary intensity profile ’ for trapping assay should not be regarded, in general, as a trendsetter henceforth because properties of other arbitrary profile beams from lasers might vary from case to case.
PPT Slide
Lager Image
Schematic of experimental optical trapping set-up in combination with custom-built Michelson interferometer for reflection phase-microscopy of RBC. BS: beam splitter; DM: dichroic mirror.
PPT Slide
Lager Image
Intra-cavity block diagram of diode pumped solid state laser source when unbolted to investigate the imperfections; M1 and M2 are mirrors.
The emitted beam was captured by a CCD camera at the exit of the source as shown in Fig. 3(a) . Then its intensity distribution along a line (with x, y co-ordinates inscribed in Fig. 3(a) ) was computed as shown in Fig. 3(b) . Top view and 3D side view are shown in Fig. 3(c , d) , respectively. It is quite evident from the intensity distribution that there is another intensity peak besides the central and two additional peaks of lesser magnitude (see Fig. 3(b) ). Such a non-basic Gaussian mode carries far-field divergence larger than the far-field divergence of a basic Gaussian mode provided beam waist radii are constant. Beam-divergence of the laser beam at far-field can be characterized quantitatively using the M 2 - parameter which has been computed for this case following the approach illustrated by Ruschin et al . [50] . The intensity profile of our beam (with blue colour) shown in Fig. 3(b) is partially similar to the rectangular and top-hat type of beam profiles as described in Ref. 50 (see Table 1 of ref. 50). M 2 - parameter for various beam profiles have been estimated numerically by Ruschin et al . [50] . Considering reference to these values, M 2 - parameter for this beam must lie between 1 and infinity (see Table 1 of Ref. 50) and is calculated to be 9.4. It is imperative to highlight that no curve fitting is implemented in intensity profile of beam shown in Fig. 3(b) during evaluation of 50% intensity drop and 1/e 2 intensity drop from peak value so that actual value of M 2 - parameter could be retrieved.
PPT Slide
Lager Image
Intensity profile of the ‘arbitrary beam profile’ captured by CCD camera at the exit of source, (b) its intensity profile along a line passing through co-ordinates (1,365) and evaluation of M2 parameter (red-line indicates 50% intensity of peak value and green-line indicates 86.5% drop of peak value), (c) top view of the intensity distribution, and (d) side view of 3D intensity profile.
This beam was focused into a sample contained in a microfluidic channel by 60X oil-immersion MO ( NA = 0.95) as shown in Fig. 4 . Trapping of one RBC from a flowing sample was realized as shown in Fig. 5(d) ~ (g) . Then an interferometric probe beam of sinusoidal intensity pattern was generated by a lab-made MMI unit as shown in Fig. 6(a) which was designed by coating one face (Y-Z plane) of a cubic beam splitter (CBS) with a thin layer (100 nm) of aluminium acting as a static mirror. The other adjacent face (X-Z plane) of the CBS was aligned against a fixed mirror mounted on an X-Y-Z micrometer stage acting as sample arm. Earlier, authors have reported patterned trapping of polystyrene spheres and linearization of floating tiny clusters of single walled carbon nanotubes using this compact and robust optical unit [51] . Probe interferograms were projected onto thus trapped RBC as shown in Fig. 6(b) and information carrying complex field reflected from RBC surface was captured at the CCD camera in a single-shot manner. The recorded interferogram is analysed by algorithms written in MATLAB and results are shown in Fig. 6(c) ~ (d) . In this experiment, two different wavelengths; one for trapping (wavelength @ 532 nm) and other for generating probe interferograms (wavelength @ 632 nm) were focused simultaneously by the same MO (60X) which has generated automatically two focal planes longitudinally separated by ~0.34 μm by virtue of wavelength difference. Individual focal planes can be calculated by using the relation 2 /( NA ) 2 and are found to be 1.76 μm and 2.1 μm, respectively, where ‘ n ’ is the RI of the object space, ‘NA’ is the numerical aperture of MO. Because of two different focal depths, the trapping beam was defocused down the propagation axis for simultaneous calibration of the probe beam using two independent beam-steering optical units.
PPT Slide
Lager Image
Schematic of prepared microfluidic channel.
PPT Slide
Lager Image
Bright-field images of trapping events: (a) slightly defocused spot of trapping beam created by 60X MO (captured by Lumenera camera), (b) a quad of trapped RBCs (see multimedia S1), (c) reshaping of one of the RBC along K-vector when the very same trapping beam of arbitrary beam profile was focused by 100X lens, (d) - (g) time sequence snapshots of trapping of RBC by 60X lens. Snapshots (a-b) are extracted from multimedia captured by Lumenera camera whereas (c-h) are still-snapshots captured by another CCD camera (SANSOVATION-SamBa EZ-130).
PPT Slide
Lager Image
Recording of images by CCD camera and their analysis for phase-map. (a) Recorded sinusoidal probe-interferogram generated by custom-built Michelson interferometer, and (b) reflected interferogram from surface of trapped RBC. (c) Image of 3D unwrapped phase map of RBC retrieved from (b), and (d) is 2D phase-profile along diameter of RBC.
Sample channel was prepared by cleaning first the glass-substrate (Corning glass) and then coating it with 10% BSA (Bovine Serum Albumin) to make it biocompatible (shown in Fig. 4 ). Channel-depth of 30 μm was maintained by sandwiching the Mylar spacer (30 μm thick) on three sides of a glass-substrate and then a cover slip (Corning glass, No. 1.5) was affixed over it. Channel-depth of 30 μm is quite sufficient for tumbling of a RBC about a tangent at its circumference, i.e., maximum depth required is twice the diameter of a RBC (~8 μm). The fresh blood sample was procured from a hospital in a sterile tube containing a preservative solution of citrate, phosphate, dextrose, and adenosine (CPDA-1) which prevents blood coagulation by chelating calcium and supplies the nutrient glucose and adenosine triphosphate while maintaining a low pH. Buffer was maintained isotonic to preserve the shape memory of RBCs and was further diluted (200 times approximately) in saline water (RI=1.348) composed of NaCl (136 mM), KCI (2.68 mM) as per requirement [52] . The diluted test sample was injected into the channel by micropipette which was finally mounted on a translational stage for investigation. The trapping beam was collimated to ~6 mm width using an assembly of lenses so that it can exactly fill the back aperture of 60X oil-immersion MO (5.9 mm diameter) to retain the ‘arbitrariness’ of the trapping beam without truncation otherwise an overfilled back-aperture would truncate the geometry of the beam.
A test sample containing freely floating RBCs was imaged by 60X MO as shown in Fig. 5(d and h) and then the very same MO was used to focus the trapping beam of arbitrary intensity profile as well as probe interferograms. The exact focus spot was slightly adjusted along the propagation axis synchronising with the probe beam which was imaged by the same MO as shown in Fig. 5(a) . In fact, this slightly defocused spot is the actual trapping plane which is adjusted within the channel-depth (30 μm). One of the floating RBCs (marked by the red circle in Fig. 5(d) ) was targeted and trapped ( Fig. 5(e ~ g) ) when the laser beam was turned on. As the beam was turned off, trapped RBCs detached and floated away as evident from Fig. 5(h) . Orientation of trapped RBCs shown in Fig. 5(b) is similar to Fig. 5(d ~ g) which is a snapshot extracted from multimedia S1 (1.21 GB movie) captured by a different recording camera Lumenera Infinity 2 (poor resolution) because the SANSOVATION-SamBa EZ-130 camera was capable to record only still-snapshots of good resolution ( Fig. 5(c ~ h) ). A quad of RBC in Fig. 5(b) remains trapped while other RBCs were observed moving within the sample (see multimedia S1). The behaviour of a RBC in optical trapping using Gaussian profile of laser beam has already been investigated extensively by several researchers [6 - 8 , 53 - 57] . Mohanty et al . [6 , 53] observed the reorientation and then rotation of RBC under optical trapping on account of shape change due to variation in the osmolality of the buffer. But Dharmadhikari et al . has observed the customary behaviour of folding (into rod shape) and unfolding of RBCs which was not observed by Mohanty et al . and concluded that it depends mainly on the power of trapping laser [7 , 8 , 57] . We have also observed the re-orientation behaviour of a RBC (see Fig. 5(c) marked by yellow dotted circle) which is consistent with the report by Mohanty et al . (see Fig. 2 of Ref. 6) when the very same trapping beam was tightly focused by 100X lens. No folding and unfolding was observed as invoked by Dharmadhikari et al. Similarly, we didn’t observe any photo-damage to an RBC trapped by the beam of ‘ arbitrary intensity profile ’ emitting at 532 nm that is further validated computationally using measured laser power at 2.25 mW on the exit of 60X objective (ED-300 OPHIR power-meter, Israel). This magnitude is very small in comparison to 60-85 mW (@ 1064 nm) recorded at the trapping plane using 100X lens by Mohanty et al . during viscous orientation torque calculation of RBC suspended in the isotonic buffer [6] . Real-time trapping mechanism of a RBC by 60X MO is depicted in Fig. (5) . One of the remarkable features of the present trapping outcome is that the symmetry plane of trapped RBC is lying nearly perpendicular to beam propagation which is in sharp-contrast to all the previously available literature where RBC reshapes to bullet shape along the propagation axis to minimize its free energy [7 , 8 , 56 , 57] . Factors affecting RBC under optical trapping at 532 nm wavelength have been addressed as below.
First, the RBC was optically trapped therefore, perturbations due to substrate-sample contact and its own weight [58 - 60] can be negated which is one of the advantages of this arbitrary beam in trapping.
Secondly, laser induced heating of RBCs under optical trapping depends primarily on the power density reaching to sample which is maximum if a Gaussian beam is steeply focused. But in this case, transversal distribution of beam intensity at the defocused focal plane is believed to create not a unique intensity gradient but enlarged according to the arbitrariness of beam-profile. In other words, total beam power is diffused spatially, thus creating multiple trap points in close proximity to each other that besides preventing RBC from radiation stress, maintained the radial plane of trapped RBC perpendicular to the laser direction in similar fashion as one-, two- and four-trap systems were used by G. Rusciano to trap, stretch and orient single RBCs [61] . To ascertain the magnitude of power reaching to the trapping plane, we computed it using numerical values of wavelength of trapping light (532 nm), magnification (60X), NA = 0.95 of focusing lens, spot-size of laser focus and bending elastic modulus (B = 2 × 10 -19 Nm ) of a normal RBC [7 , 55 , 62 - 64] and found it to be 3.65 mW. Further, laser power was experimentally measured to 2.25 mW (ED-300, OPHIR-Israel) at the exit of 60X objective. Thus both experimental and computed values of threshold power delivered to RBC under optical trapping are well within the reported value of 20 mW at 1064 nm wavelength using 100X lens by Dharmadhikari et al . [7 (p. 69)]. So the authors believe faithfully that such a low laser power reaching to RBC can hardly invoke any deleterious effects throughout a time-span (~2 min.) of the investigation following the aforementioned studies [6 - 8 , 53 - 57] .
Thirdly, laser induced heating of RBCs under optical trapping also depends on duration of exposure and wavelength. Approximately 2 minutes duration of laser exposure to RBC under trapping as evident from Fig. 5(e ~ g) can be believed to have negligible influence. Anatomically, RBCs are present in blood at a concentration of (4-6)×10 9 / mL with 33% of haemoglobin by volume which is the main component responsible for absorption at 532 nm. But water being the major component, laser induced heating is relatively mild for laser wavelength window 200-1100 nm where absorption by water is small as reported earlier by Tretheway et al . [65] . Further G. Rusciano performed Raman spectroscopy of trapped RBC from a physiological point of view at 532 nm wavelength yet no photo-damage was reported [61] . From the above two reports, we conclude firmly that trapping of RBC by 60X objective lens using a beam of ‘ arbitrary intensity profile ’ at 532 nm wavelength for a short-span of time could not generate any heating to induce photo-damage.
Finally, real-time rotation/buckling of a RBC under optical trapping is a direct signature of conversion of absorbed power. However, when the incident beam’s central axis is not aligned with the symmetry axis of RBC or RBC itself is asymmetric (crenate), free energy minimization constraint causes it to buckle/re-shape and rotates it about the axis. That is why malaria-infected/dead RBCs of asymmetric shape exhibit self-rotation under optical trapping as reported earlier [53] . It has already been proved in literature that the radiation force of the order of 1-50 piconewtons (pN) generated in Gaussian beam driven OTs is capable to deform an RBC physically without causing cell death [7 (p. 68), 55] . Scientifically, the electric field component (
PPT Slide
Lager Image
) of the laser beam tries to induce a dipole moment into the RBC according to the relation P = α × E , where α is the polarizability tensor of the RBC. But owing to the intrinsic asymmetry of the RBC, polarizability components parallel and perpendicular to RBC axis ( αII and α respectively) are significantly different from each other which causes light-induced torque generation τ = P × E . However, energy minimization constraint forces the RBC to reshape to a bullet and align along the direction of beam propagation [8 , 56] . According to the above relation, generated torque is equal to the product of the polarizability tensor and square of the electric field which in turn translates to time-average beam intensity (trapping power). But the above torque relation is not applicable here for rotationally asymmetric beam of ‘ arbitrary profile ’ because now the E-field is a complex tensor quantity fully described by Maxwell’s Stress tensor
PPT Slide
Lager Image
which will provide radiation pressure force and torque [66 - 69] . More pedagogic theoretical description characterizing every electromagnetic field of arbitrary profile and arbitrary orientation, and hence the estimation of radiation pressure force and torque on spheroidal particles using a set of respective beam shape coefficients (BSCs) under generalized Lorentz-Mie theory (GLMT) can be found elsewhere [70 - 79] . Re-shaping of RBCs to bullet shape was not observed in this case as recorded power at the exit of 60X MO was moderate (2.25 mW). But as the magnification of MO is changed from 60X to 100X keeping beam power unchanged, the fluence reaching to the trapping plane increases causing re-shaping of RBC to bullet shape which is evident from Fig. 5(c) . In short, the word ‘ arbitrary profile ’ and its electrodynamics in OTs is confined to the current laser beam only and can be generalized following aforementioned theoretical computation [66 - 79] . In-depth illustration characterising the arbitrary profile beam is itself a separate full-fledged research domain which is underway.
Now, to obtain the phase-map of each pixel, probe interferograms reflected from a trapped RBC were captured at a CCD camera and analysed by algorithms written in MATLAB. No reflection can occur from the intermediate cytoplasmic medium because it is transparent for visible wavelength of probe light (632 nm) thereby assuming no wavefront aberration from phase object (RBC) [80] . Total OPD between two co-propagating reflected interferograms: one from the top surface and other from the bottom surface of RBC (after transmission from the top surface) is twice the thickness of discotic RBC (0.9 μm - 2.5 μm) i.e., 1.8 μm - 5μm. Therefore, phase-change must appear significant beyond the 3 rd position of the decimal following Eq. (5) because probe wavelength (632 nm) is three orders less than the physical thickness variation of the RBC (1.8 μm - 5 μm.). In this case, we are considering the reflection of probe interferogram from the bilayer surface of the RBC which remains invariant untill the pathophysiological state changes. However, such a large fluctuation in OPD between two reflected interferograms (1.8 μm - 5 μm.) is well within the coherence length of trapping laser (10 μm) therefore interference will essentially take place. At the same time weak reflection (~4%) of the probe interferogram from glass-substrate cannot be ruled out. But taking into account the location of a trapped RBC lying anywhere in the middle of 30 μm channel depth, the OPD between glass-substrate and bottom surface of a trapped RBC is far greater than the coherence length of source and therefore no interference can take place. Phase-map along the diameter of an RBC is computed by avoiding the linear phase by proper alignment of the set-up as shown in Fig. 6(c) which depicts the biconcave shape of the upper surface of RBC. One dimensional phase-profile of RBC along its diameter is computed and shown in Fig. 6(d) where a variation in phase-profile values can be seen which is due to different constituents present within RBC. To validate our phase-imaging study statistically, several experimental results of simultaneous trapping and projection of the probe interferogram must be recorded. However, in this case, it could be made feasible either by trapping various RBCs from a single beam optical trap at various time lengths and then recording the respective interferograms or multiple trap centers must be realized at a time to trap various RBCs by enlarging the area of the probe interferogram. In the former case, a somewhat longer period of experiment is needed to cover various RBCs one by one for statistical analysis. But in the latter case, a collimated trapping beam of homogeneous transversal intensity must be projected on a diffractive optical element to realize multiple trap centers with equal trap stiffness which is not possible in the present case as the trapping beam contains transversal inhomogeneity in intensity distribution. Following the first case, we had recorded a few experimental observations for different RBCs of the same sample and one of the best results is shown in Fig. 5 with its analysis in Fig. 6 . We analysed these sets of observations and obtained maximum fluctuation in estimated phase of the order of ± 0.1 radian (in Fig. 6(c) ) at the wall of the RBC.
The spatial nature of the laser beam of ‘ arbitrary intensity profile ’ generated by cavity itself due to cumulative effects of both the misalignment of optics and inhomogeneity of active material (KTP) has been utilized for photodamage-free trapping of RBCs. The advantage to utilize such beams in OTs is to investigate the role of low-cost laser sources (defective) whereby beam quality factor does not play a significant role. The implicit repercussion of receiving noninevitable alterations in innate morphology of RBCs due to their own weight when adhered to glass substrate in interferometric phase-microscopy schemes, has been obviated here. Thus trapped RBC was probed with sinusoidal interferograms generated by lab-made Michelson interferometer and the phase-map was evaluated using the Hilbert transform fringe analysis technique from recorded single-shot reflected interferograms. In conclusion, a cost-effective and label-free phase-imaging modality of trapped RBC has been reported. Efficient reflective coating on a cubic beam-splitter to generate probe interferograms of good quality and assimilation of technically updated diode lasers making the system inexpensive can be regarded as an improved version of phase-imaging platform in the future.
Financial support in the form of fellowship from the ministry of human resource and development Government of India during Ph.D. at Indian Institute of Technology Delhi is greatly acknowledged.
Wang Y. L. , Discher D. E. 2008 Methods in Cell Biology - Cell Mechanics Elsevier Press
Suresh S. (2006) “Mechanical response of human red blood cells in health and disease: some structure-property-function relationships,” J. Mater. Res. 21 1871 - 1877    DOI : 10.1557/jmr.2006.0260
Park Y. K. , Silva M. D. , Popescu G. , Lykotrafitis G. , Choi W. , Feld M. S. , Suresh S. (2008) “Refractive index maps and membrane dynamics of human red blood cells parasitized by Plasmodium falciparum,” Proc. National Academy of Sciences 105 13730 - 13735    DOI : 10.1073/pnas.0806100105
Byun H. S. , Hillman T. R. , Higgins J. M. , Silva M. D. , Peng Z. , Dao M. , Dasari R. R. , Suresh S. , Park Y. K. (2012) “Optical measurement of biomechanical properties of individual erythrocytes from a sickle cell patient,” Acta Biomaterialia 8 4130 - 4138    DOI : 10.1016/j.actbio.2012.07.011
Williams K. K. 2010 Hematology McGraw-Hill Medical New York, USA
Cuche E. , Marquet P. , Depeursinge C. (2005) “Dynamics of interaction of RBC with optical tweezers,” Opt. Express 13 4745 - 4751    DOI : 10.1364/OPEX.13.004745
Ghosh A. , Sinha S. , Dharmadhikari J. A. , Roy S. , Dharmadhikari A. K. , Samuel J. , Sharma S. , Mathur D. (2006) “Euler buckling-induced folding and rotation of red blood cells in an optical trap,” Phys. Biol. 3 67 - 73    DOI : 10.1088/1478-3975/3/1/007
Dharmadhikari J. A. , Roy S. , Dharmadhikari A. K. , Sharma S. , Mathur D. (2004) “Naturally occurring, optically driven, cellular rotor,” Appl. Phys. Lett. 85 6048 - 6051    DOI : 10.1063/1.1836874
Gifford S. C. , Derganc J. , Shevkoplyas S. S. , Yoshida T. , Bitensky M. W. (2006) “A detailed study of time-dependent changes in human red blood cells: from reticulocyte maturation to erythrocyte senescence,” British Journal of Heamatology 135 395 - 404    DOI : 10.1111/j.1365-2141.2006.06279.x
Parshina E. Y. , Sarycheva A. S. , Yusipovich A. I. , Brazhe N. A. , Goodilin E. A. , Maksimov G. V. (2013) “Combined Raman and atomic force microscopy study of hemoglobin distribution inside erythrocytes and nanoparticle localization on the erythrocyte surface,” Laser Phys. Lett. 10 075607 -    DOI : 10.1088/1612-2011/10/7/075607
Friebel M. , Meinke M. (2005) “Determination of the complex refractive index of highly concentrated haemoglobin solutions using transmittance and reflectance measurements,” J. Biomed. Opt. 10 064019 -    DOI : 10.1117/1.2138027
Popescu G. 2011 Quantitative Phase Imaging of Cells and Tissues McGraw-Hill Professional
Lee K. , Kim K. , Jung J. , Heo J. , Cho S. , Lee S. , Chang G. , Jo Y. , Park H. , Park Y. (2013) “Quantitative phase imaging techniques for the study of cell pathophysiology: from principles to applications,” Sensors 13 4170 - 4191    DOI : 10.3390/s130404170
Cuche E. , Marquet P. , Depeursinge C. (1999) “Simultaneous amplitude contrast and quantitative phase-contrast microscopy by numerical reconstruction of Fresnel off-axis holograms,” Appl. Opt. 38 6694 - 7001
Carl D. , Kemper B. , Wernicke G. , von Bally G. (2004) “Parameter optimized digital holographic microscope for high-resolution living cell analysis,” Appl. Opt. 43 6536 - 6544    DOI : 10.1364/AO.43.006536
Carl C. J. , Yu L. F. , Lo C. M. , Kim M. K. (2005) “High-resolution quantitative phase-contrast microscopy by digital holography,” Opt. Express 13 8693 - 8698    DOI : 10.1364/OPEX.13.008693
Kemper B. , Vollmer A. , Rommel C. E. , Schnekenburger J. , von Bally G. (2011) “Simplified approach for quantitative digital holographic–phase contrast imaging of living cells,” J. Biomed. Opt. 16 026014:1 - 4
Marquet P. , Rappaz B. , Magistretti P. , Cuche E. , Emery Y. , Colomb T. , Depeursinge C. (2005) “Digital holographic microscopy: a non-invasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy,” Opt. Lett. 30 468 - 470    DOI : 10.1364/OL.30.000468
Singh A. S. G. , Anand A. , Leitgeb R. A. , Javidi B. (2012) “Lateral shearing digital holographic imaging of small biological specimens,” Opt. Express 20 23617 - 23622    DOI : 10.1364/OE.20.023617
Popescu G. , Deflores L. P. , Vaughan J. C. , Badizadegan K. , Iwai H. , Dasari R. R. , Feld M. S. (2004) “Fourier phase microscopy for investigation of biological structures and dynamics,” Opt. Lett. 29 2503 - 2505    DOI : 10.1364/OL.29.002503
Ikeda T. , Popescu G. , Dasari R. R. , Feld M. S. (2005) “Hilbert phase microscopy for investigating fast dynamics in transparent systems,” Opt. Lett. 10 1165 - 1167
Lue N. , Choi W. , Popescu G. , Yaqoob Z. , Badizadegan K. , Dasari R. R. , Feld M. S. (2009) “Live cell refractometry using Hilbert phase microscopy and confocal reflectance microscopy,” J. Phys. Chem. A 113 13327 - 13330    DOI : 10.1021/jp904746r
Popescu G. , Ikeda T. , Dasari R. R. , Feld M. S. (2006) “Diffraction phase microscopy for quantifying cell structure and dynamics,” Opt. Lett. 31 775 - 777    DOI : 10.1364/OL.31.000775
Choi W. , Fang-Yen C. , Badizadegan K. , Oh S. , Lue N. , Dasari R. R. , Feld M. S. (2007) “Tomographic phase microscopy,” Nature Meth. 4 717 - 719    DOI : 10.1038/nmeth1078
Park Y. K. , Best C. A. , Badizadegan K. , Dasari R. R. , Feld M. S. , Kuriabova T. , Henle M. L. , Levine A. J. , Popescu G. (2010) “Measurement of red blood cell mechanics during morphological changes,” Proc. National Academy of Sciences 107 6731 - 6736    DOI : 10.1073/pnas.0909533107
Sarmis M. , Simon B. , Debailleul M. , Colicchio B. , Georges V. , Delaunay J. J. , Haeberlé O. (2010) “High resolution reflection tomographic diffractive microscopy,” J. Mod. Opt. 57 740 - 745    DOI : 10.1080/09500341003624743
Yu X. , Hong J. , Liu C. , Kim M. K. (2014) “Review of digital holographic microscopy for three dimensional profiling and tracking,” Opt. Eng. 52 112306 -
Rappaz B. , Charrière F. , Depeursinge C. , Magistretti P. , Marquet P. (2008) “Simultaneous cell morphometry and refractive index measurement with dual-wavelength digital holographic microscopy and dye-enhanced dispersion of perfusion medium,” Opt Lett. 33 744 - 746    DOI : 10.1364/OL.33.000744
Charrière F. , Pavillon N. , Colomb T. , Depeursinge C. , Heger T. J. , Mitchell E. A. D. , Marquet P. , Rappaz B. (2006) “Living specimen tomography by digital holographic microscopy: Morphometry of testate amoeba,” Opt. Express 14 7005 - 7013    DOI : 10.1364/OE.14.007005
Kemper B. , Langehanenberg P. , Höink A. , von Bally G. , Wottowah F. , Schinkinger S. , Guck J. , Käs J. , Bredebusch I. , Schnekenburger J. , Schütze K. (2010) “Monitoring of laser micromanipulated optically trapped cells by digital holographic microscopy,” J. Biophoton. 3 425 - 431    DOI : 10.1002/jbio.201000035
Kemper B. , Bredebusch I. , Domschke W. , Kosmeier S. , Langehanenberg P. , Schnekenburger J. , von Bally G. (2007) “Integral refractive index determination of living suspension cells by multifocus digital holographic phase contrast microscopy,” J. Biomed. Opt. 12 054009 -    DOI : 10.1117/1.2798639
Kemper B. , Wibbeling J. , Kastl L. , Schnekenburger J. , Ketelhut S. (2015) “Continuous morphology and growth monitoring of different cell types in a single culture using quantitative phase microscopy,” Proc. SPIE 952902 9529 1 - 7
Kemper B. , von Bally G. (2008) “Digital holographic microscopy for live cell applications and technical inspection,” Appl. Opt. 47 A52 - A61    DOI : 10.1364/AO.47.000A52
Evans E. , Fung Y. C. (1972) “Improved measurements of the erythrocyte geometry,” Micro Vascular Research 4 335 - 347
Rappaz B. , Barbul A. , Emery Y. , Korenstein R. , Depeursinge C. , Magistretti P. J. , Marquet P. (2008) “Comparative study of human erythrocytes by digital holographic microscopy, confocal microscopy and impedance volume analyzer,” Cytometry A 73 895 - 903
Lee K. R. , Park Y. K. (2014) “Quantitative phase imaging unit,” Opt. Lett. 39 3630 - 3633    DOI : 10.1364/OL.39.003630
Kim K. H. , Yaqoob Z. , Lee K. R. , Kang J. W. , Choi Y. , Hosseini P. , So P. T. C. , Park Y. K. (2014) “Diffraction optical tomography using a quantitative phase imaging unit,” Opt. Lett. 39 6935 - 6938    DOI : 10.1364/OL.39.006935
Colomb T. , Krivec S. , Hutter H. , Akatay A. A. , Pavillon N. , Montfort F. , Cuche E. , Kühn J. , Depeursinge C. , Emery Y. (2010) “Digital holographic reflectometry,” Opt. Express 18 3719 - 3731    DOI : 10.1364/OE.18.003719
Merola F. , Miccio L. , Memmolo P. , Di Caprio G. , Galli A. , Puglisi R. , Balduzzi D. , Coppola G. , Netti P. , Ferraro P. (2013) “Digital holography as a method for 3D imaging and estimating the biovolume of motile cells,” Lab on a Chip 13 4512 - 4516    DOI : 10.1039/c3lc50515d
Panah M. D. , Zwick S. , Schaal F. , Warber M. , Javidi B. , Osten W. (2010) “3D holographic imaging and trapping for non-invasive cell identification and tracking,” J. Display Technol. 6 490 - 499    DOI : 10.1109/JDT.2010.2043499
Moradi A. R. , Ali M. K. , Panah M. D. , Anand A. , Javidi B. (2010) “Detection of calcium induced morphological changes of living cells using optical traps,” IEEE Photon. J. 2 775 - 783    DOI : 10.1109/JPHOT.2010.2068283
Habaza M. , Gilboa B. , Roichman Y. , Shaked N. T. (2015) “Tomographic phase microscopy with 180° rotation of live cells in suspension by holographic optical tweezers,” Opt. Lett. 40 1881 - 1884    DOI : 10.1364/OL.40.001881
Cardenas N. , Mohanty S. K. (2013) “Decoupling of geometric thickness and refractive index in quantitative phase microscopy,” Opt. Lett. 38 1007 - 1009    DOI : 10.1364/OL.38.001007
Kim Y. , Shim H. , Kim K. , Park H. J. , Jang S. , Park Y. K. (2014) “Profiling individual human red blood cells using common-path diffraction optical tomography,” Scientific Reports 4 6659 -    DOI : 10.1038/srep06659
Kim K. H. , Yoon J. G. , Park Y. K. (2015) “Simultaneous 3D visualization and position tracking of optically trapped particles using optical diffraction tomography,” Optica 2 343 - 346    DOI : 10.1364/OPTICA.2.000343
Jang J. , Bae C. Y. , Park J.-K. , Ye J. C. (2010) “Self-reference quantitative phase microscopy for microfluidic devices,” Opt. Lett. 35 514 -    DOI : 10.1364/OL.35.000514
Shaked N. T. (2012) “Quantitative phase microscopy of biological samples using a portable interferometer,” Opt. Lett. 37 2016 - 2018    DOI : 10.1364/OL.37.002016
Shaked N. T. , Zhu Y. , Badie N. , Bursac N. , Wax A. (2010) “Reflective interferometric chamber for quantitative phase imaging of biological sample dynamics,” J. Biomed. Opt. 15 030503 -    DOI : 10.1117/1.3420179
Memmolo P. , Miccio L. , Merola F. , Gennari O. , Netti P. A. , Ferraro P. (2014) “3D morphometry of red blood cells by digital holography,” Cytometry Part A 85 1030 - 1036    DOI : 10.1002/cyto.a.22570
Ruschin S. , Yaakobi E. , Shekel E. (2011) “Gaussian content as a laser beam quality parameter,” Appl. Opt. 50 4376 - 4381    DOI : 10.1364/AO.50.004376
Kumar R. , Shakher C. , Mehta D. S. (2010) “Compact interferometric optical tweezers for patterned trapping and manipulation of polystyrene spheres and SWCNTs,” J. Mod. Opt. 57 1157 - 1162    DOI : 10.1080/09500340.2010.500104
Mohanty S. K. , Gupta P. K. (2004) “Self-rotation of red blood cells in optical tweezers: prospects for high throughput malaria diagnosis,” Biotechnol. Lett. 26 971 - 974    DOI : 10.1023/B:BILE.0000030041.94322.71
Grover S. C. , Gauthier R. C. , Skirtach A. G. (2000) “Analysis of the behaviour of erythrocytes in an optical trapping system,” Opt. Express 7 533 - 539    DOI : 10.1364/OE.7.000533
Dao M. , Lim C. T. , Suresh S. (2003) “Mechanics of the human red blood cell deformed by optical tweezers,” J. Mech. Phys. Solids 51 2259 - 2280    DOI : 10.1016/j.jmps.2003.09.019
Khan M. , Soni H. , Sood A. K. (2009) “Optical tweezer for probing erythrocyte membrane deformability,” Appl. Phys. Lett. 95 233703 -    DOI : 10.1063/1.3272269
Dharmadhikari J. A. , Mathur D. (2004) “Using an optical trap to fold and align single red blood cells,” Current Science 86 1432 - 1437
Flugge W. 1962 Handbook of Engineering Mechanics McGraw-Hill New York, USA
Hutchinson J. W. (1967) “Imperfection sensitivity of externally pressurized spherical shells,” J. Appl. Mech. 8 49 - 55
Timoshenko S. , Gere J. M. 1961 Theory of Elastic Stability 2nd ed. McGraw-Hill New York, USA
Rusciano G. (2010) “Experimental analysis of Hb oxy-deoxy transition in single optically stretched red blood cells,” Physica Medica 26 233 - 239    DOI : 10.1016/j.ejmp.2010.02.001
Hochmuth R. M. , Skalak R. , Chien S. 1987 Handbook of Bioengineering McGraw-Hill New York, USA “Properties of red blood cells,”
Sleep J. , Wilson D. , Simmons R. , Gratzer W. (1999) “Elasticity of the red cell membrane and its relation to hemolytic disorders: an optical tweezers study,” Biophys. J. 77 3085 - 3095    DOI : 10.1016/S0006-3495(99)77139-0
Lim C. T. , Dao M. , Suresh S. , Sow C. H. , Chew K. T. (2004) “Large deformation of living cells using laser traps,” Acta Material 52 1837 - 1845    DOI : 10.1016/j.actamat.2003.12.028
Neve N. , Kohles S. S. , Winn S. R. , Tretheway D. C. (2010) “Manipulation of suspended single cells by microfluidics and optical tweezers,” Cellular and Molecular Bioengineering 3 213 - 228    DOI : 10.1007/s12195-010-0113-3
Wang X. , Wang X.-B. , Gascoyne P. R. C. (1997) “General expressions for dielectrophoretic force and electrorotational torque derived using the Maxwell stress tensor method,” Journal of Electrostatics 39 277 - 295    DOI : 10.1016/S0304-3886(97)00126-5
Barton J. P. , Alexander D. R. , Schaub S. A. (1989) “Theoretical determination of net radiation force and torque for a spherical particle illuminated by a focused laser beam,” J. Appl. Phys. 66 4594 - 4602    DOI : 10.1063/1.343813
Collett W. L. , Ventrice C. A. , Mahajan S. M. (2003) “Electromagnetic wave technique to determine radiation torque on micromachines driven by light,” Appl. Phys. Lett. 82 2730 - 2732    DOI : 10.1063/1.1567042
Paul R. , Kaler K. V. I. S. (1993) “Effects of particle shape on electromagnetic torques: A comparison of the effective-dipolemoment method with the Maxwell-stress-tensor method,” Phys. Rev. E 48 1491 - 1496    DOI : 10.1103/PhysRevE.48.1491
Gouesbet G. , Maheu B. , Gréhan G. (1988) “Light scattering from a sphere arbitrarily located in a Gaussian beam, using a Bromwich formulation,” J. Opt. Soc. Am. A 5 1427 - 1443    DOI : 10.1364/JOSAA.5.001427
Polaert H. , Gréhan G. , Gouesbet G. (1998) “Forces and torques exerted on a multilayered spherical particle by a focused Gaussian beam,” Opt. Commun. 155 169 - 179    DOI : 10.1016/S0030-4018(98)00220-X
Lock J. A. (2004) “Calculation of the radiation trapping force for laser tweezers by use of generalized Lorenz-Mie theory. I. Localized model description of an on-axis tightly focused laser beam with spherical aberration,” Appl. Opt. 43 2532 - 2544    DOI : 10.1364/AO.43.002532
Lock J. A. (2004) “Calculation of the radiation trapping force for laser tweezers by use of generalized Lorenz-Mie theory. II. On-axis trapping force,” Appl. Opt. 43 2545 - 2554    DOI : 10.1364/AO.43.002545
Rohrbach A. , Stelzer E. H. K. (2002) “Trapping forces, force constants, and potential depths for dielectric spheres in the presence of spherical aberrations,” Appl. Opt. 41 2494 - 2507    DOI : 10.1364/AO.41.002494
Ke P. C. , Gu M. (1998) “Characterization of trapping force in the presence of spherical aberration,” J. Mod. Opt. 45 2159 - 2168    DOI : 10.1080/09500349808231752
Going R. W. , Conover B. L. , Escuti M. J. (2008) “Electrostatic force and torque description of generalized spheroidal particles in optical landscapes,” Proc. SPIE 7038 703826 -
Bonessi D. , Bonin K. , Walker T. (2007) “Optical forces on particles of arbitrary shape and size,” J. Opt. A: Pure Appl. Opt. 9 S228 - S234    DOI : 10.1088/1464-4258/9/8/S16
Xu F. , Ren K. , Gouesbet G. , Cai X. , Gréhan G. (2007) “Theoretical prediction of radiation pressure force exerted on a spheroid by an arbitrarily shaped beam,” Phys. Rev. E 75 026613:1 - 14
Xu F. , Ren K. , Gouesbet G. , Gréhan G. , Cai X. (2007) “Generalized Lorenz-Mie theory for an arbitrarily oriented, located, and shaped beam scattered by a homogeneous spheroid,” J. Opt. Soc. Am. A 24 119 - 131    DOI : 10.1364/JOSAA.24.000119
Bae Y. S. , Song J. I. , Har D. , Kim D. Y. (2015) “Beam propagation analysis on thickness measurements in quantitative phase microscopy,” Opt. Rev. 22 532 - 538    DOI : 10.1007/s10043-015-0112-7