Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Sub-nanometer measurement of transient structural changes in dye-doped polystyrene microspheres

Open Access Open Access

Abstract

We present an interferometric spectral-domain optical coherence tomography microscopy setup to detect structural changes using interference of light reflected from different interfaces of the sample. We induce a reproducible nanometer-scale size change in dye-doped 10-µm polystyrene microspheres by the release of Stokes shift energy of dye molecules inside the microspheres, excited by a modulated 532-nm laser. The resulting optical path length difference was measured with a sensitivity of $0.4~{{\mathrm {pm}}/{\sqrt {\mathrm {Hz}}}}$ limited by photodetection noise, and reveals elastic as well as inelastic responses, which opens up possibilities for measuring the response of cell-sized biological objects.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

The microscopic physical properties of single cells attract significant research interest, because of the ability to measure e.g. stiffness and motility of cells [1,2]. The methods to detect these physical properties are usually invasive, such as AFM to measure the stiffness [2,3] or electrode and fluorescence based [4] electrophysiology. Label-free optical techniques could provide an alternative to non-invasively recording of neuronal activities.

Quantitative measurement of the optical phase is emerging as a powerful, label-free approach to imaging cells and tissues, especially as it combines qualities found in microscopy, holography, and light scattering techniques: nanoscale sensitivity of morphology and dynamics, non-destructive imaging of completely transparent structures and quantitative signals based on intrinsic contrast [5]. By combining holography and microscopy one can perform sensitive measurements of the thickness and refractive index of biological specimens, with sub-wavelength accuracy [6]. Early on, such quantitative phase measurements were performed by single point scanning techniques, building on the advancement of optical coherence tomography (OCT) [7]. Two classes of OCT are swept-source OCT (SS-OCT), which uses a tunable swept laser [810] and spectral-domain OCT (SD-OCT) which uses a broadband light source [11,12]. While SS-OCT has some advantages in imaging applications, SD-OCT is preferable in applications with time-dependent signals as it can take an entire spectrum in a single shot. In 2001, Yang et al. developed an OCT system that could detect motion within a cell with 3.8 nm sensitivity for a slow sampling [13]. Many studies were carried out using the same methods to detect the optical and mechanical change of neurons during action potentials (APs) [1418]. In all the studies conducted to detect the cellular activities and motion, the reference beam was a mirror in a separate beam path such as a Michelson configuration or common path configurations, for instance, one side of the coverslip, which holds the sample. In these configurations besides the influence of the noise on the sensitivity of the experiment, the motion of the object relative to the coverslip exhibit crosstalk with changes within the object itself. This crosstalk problem becomes worse if the object is closer to the coverslip than the axial resolution of the SD-OCT, as in this case it is not possible to separate the desired signal from the whole object’s motion.

In this work, we aim to measure sub-nanometer transient path length changes of cell-sized objects, where we use polymer microspheres as a model system. We have constructed an OCT setup without an external reference beam, instead employing a common path configuration in which the phase is retrieved from the interference of the reflected light from the top and bottom surface of the microsphere. To show the high sensitivity of SD-OCT in a dynamical setting we use a high sampling rate of 80 kHz.

As reference objects of cellular size which can be deformed in a controllable way in real-time, we use polystyrene (PS) microspheres, with a diameter of 10 µm, doped with fluorescent dye. When the dye molecules absorb light, part of the energy (corresponding to the Stokes shift) is released as heat [19]. The OCT system does not detect fluorescence, as the fluorescent light is not coherent with the source [20], moreover, the fluorescent light is in a different wavelength range. The released heat changes the density, volume, and refractive index of the object and we measure the resulting optical path length change.

2. Methods and materials

2.1 Spectral-domain OCT setup

Light from a superluminescent diode (Superlum Inc., $\lambda _0$=840 nm, $\Delta \lambda$=160 nm) is divided by a fiber coupler, where one output is coupled into an inverted microscope (Nikon E2000 with a 20x NA=0.75 objective), which contains the sample (Fig. 1). The reflected light from the sample is coupled back into the same fiber, and part of it passes the fiber coupler to reach a spectrometer (Wasatch Photonics, Cobra S-800) with integrated line camera (E2V OCTAPLUS). The spectrum is recorded by the spectrometer with an 80-kHz line rate and a spectral resolution of 0.15 nm. An 80/20 beamsplitter inside the microscope side port allows for simultaneous sample visualization and interferometry. To expand the microspheres by heat, a fiber-coupled, 532-nm laser (Coherent Obis) is connected to the excitation port of the microscope, overlapping with the OCT beam using a Nikon G-2B cube (Long-pass Emission). This cube normally has a dichroic mirror, an excitation filter, and a transmission filter, however in this experiment we have removed the transmission filter, so the whole bandwidth of the SLD is transmitted through the cube. While this allows the fluorescence emission to reach the spectrometer, we have confirmed that its intensity is negligible in the wavelength range used for OCT. We modulated the excitation laser with a function generator, controlled by a microcontroller (Arduino), with a modulation period of 88 ms and a duty cycle of 50%. The excitation light is collimated and focused through the (underfilled) microscope objective, so that the beam waist on the sample is $\approx$4 µm and the Rayleigh length is larger than 10 µm.

 figure: Fig. 1.

Fig. 1. Setup of the SD-OCT system. (a) The light source and detection part that includes a superluminescent light source (SLD), a fiber coupler (FC) (two practical options 90:10 or 50:50), a collimator (CL1), and spectrometer (SP). (b) Microscopy section where the light is directed to the sample by a collimated lens (CL2), BS1 (P), microscope objective, and arrives at the microsphere. (c) excitation arm with a 532-nm laser, with collimation lens (CL3). (d) Schematic of the sample (e) Measured interferogram formed by the top and bottom reflections from the microsphere and (f) its intensity Fourier transform

Download Full Size | PDF

The sample is a single polystyrene microsphere (FluoSpheres, ThermoFisher) doped with orange fluorescent dye. The excitation peak is at 540 nm, with the emission peak at 560 nm. The microspheres were suspended in 1% low melt, bio-compatible agarose gel, which is commonly used tissue imaging and cell biology [2123]. It keeps the microspheres fixed in space and reduces their thermal motion motion during heating.

The desired contribution to our signal is formed by the reflections from the bottom and top surface of the microsphere and their interference term (Fig. 1). The interferogram due to two reflectors, here arbitrarily identified as reference and sample, is given by [24,25],

$$I(k)= (\rho/e) \ s(k) \delta k \Delta t \sqrt{R_r R_s}\cos(i2nk [z_0 +\delta z (t)])$$
where $z_0 + \delta z$ is the distance between the reference and sample reflectors, $\rho$ is the detector responsivity, $e$ is electron charge, $\delta k$ is is the spectrometer spectral resolution, $\Delta t$ is exposure time, and n is the refractive index between the reference and sample reflectors. $z_0$ denotes the initial position of the sample reflector, which is determined within the source coherence length, while $\delta z(t)$ accounts for time-dependent optical path length deviations. $s(k)$ is the source spectral power density function; $R_r$ and $R_s$ are the reference and sample reflectivity, respectively. Fourier transformation of $I(k)$ yields an A-scan $\mathfrak {I}(nz)$ that has peak values at $2nz_0$,
$$\mathfrak{I}(\pm2n z_0)= (\rho/2e) S \Delta t E(2n z_0) \sqrt{R_r R_s} \exp{[{\pm} i2 k_c n \delta z(t)]},$$
where $S$ is total power, $E(2n z_0)$ is the unity-amplitude coherence envelope function and $k_c$ is the center wave number of the source. Since the phase of $\mathfrak {I}(\pm 2n z_0)$ is a linear function of $\delta z$, time-varying changes in the distance between the sample and reference reflectors can be measured with reference to the distance at time $t_0$,
$$\delta \varphi(t)=\angle \mathfrak{I}(\pm2n z_0,t) - \angle \mathfrak{I}(\pm2n z_0,t_0)=2k_c n \delta z,$$

Equation (3) shows the time varying phase of the interferograms where $\angle$ denotes angle operator. Since we use the interfaces on the top and bottom of the microsphere itself rather than using an external reference surface such as coverslip, the measured OPL, optical path difference, directly relates to the microsphere diameter and is insensitive to its position relative to the coverslip. The signal processing is exactly analogous to Ref. [24].

3. Results

3.1 Transient structure analysis of microspheres

In Fig. 2 we show the measured changes in OPL while illuminating a single microsphere with the modulated excitation laser. We observe that every pulse of the excitation laser induces a temporary reduction of the OPL, even though a thermal expansion of the microsphere is expected. The change of the OPL is proportional to the excitation laser power. To explain the negative sign of the change we invoke the Clausius-Mossotti equation which relates the refractive index to density of the object, leading to a decrease of the refractive index with temperature [26].

 figure: Fig. 2.

Fig. 2. Optical path length of two microspheres exposed to various power, as example of two different microspheres behavior (a) Microsphere showing a low amount of inelastic creep. (b) Microsphere showing significant creep. Note the different scale of the y-axis. Raw data is shown in lighter color behind the smoothed data (moving average 50 data points)

Download Full Size | PDF

Since the OPL is the product of the geometric length of the optical path and the refractive index,

$$\delta \mathrm{OPL}=2n\delta z(T) + 2z \delta n(T).$$

Using the Clausius-Mossotti relation we find,

$$\frac{d \delta \mathrm{OPL}}{dT}=2 [- \frac{(n^2-1)(n^2+2)}{2n}+n]\frac{d\delta z}{dT}$$
where the first term in brackets represents the effect of change in the refractive index and the second one the thermal expansion of the radius of the microsphere on the OPL. In the case of a polystyrene sample with n=$1.57$, the term in brackets evaluates to $-0.53$, as the effect of change in refractive index is dominant.

In addition to the consistent step-wise changes attributed to heating, we also observe a slow shift of the OPL baseline which is different for each microsphere we observed. This slow change is not due to a drift in position, as this would cause a loss of fringe visibility which we do not observe. We tentatively attribute the slow drift to inelastic deformation [27,28] of the sphere which is heated inhomogeneously by our laser. This slow creep deformation depends on the molecular properties of the polymer and the stress/strain history that polymer has been through [27,29]. In Fig. 2 we show the response of two microspheres from the same batch. Some microspheres show predominantly elastic behavior as shown in Fig. 2(a), while in Fig. 2(b) we show an example of a microsphere that undergoes a significant amount of permanent path length change. This strong difference in the inelastic behavior could be attributed to a possible different annealing history of the microspheres. Of the 36 microspheres we measured, 15 showed behavior qualitatively similar to Fig. 2(a) while 21 were similar to Fig. 2(b). The amplitude of the fast elastic response is equal in all observed microspheres.

To rule out contributions of absorption elsewhere, we repeated the experiments on non-doped PS microspheres at the same size and kept all the parameters constant. From the response curves in Fig. 3 (which are representative of several different microspheres) we see that the dye-doped microspheres show transient changes that follow the excitation laser profile, while samples without dye showed neither ageing nor any elastic response.

 figure: Fig. 3.

Fig. 3. Comparison of OPL change in dye-doped and non-absorbing microspheres experiencing the same laser power of 3.12 mW.

Download Full Size | PDF

3.2 Sensitivity

We define sensitivity as the change in the $\delta \mathrm {OPL}$ that corresponds to a signal-to-noise ratio (SNR) of 1. The SD-OCT setup samples the spectra at 80 kHz. To characterize the sensitivity, we show the measured amplitude spectral density of the OPL in the absence of laser excitation (Fig. 4). The spectral sensitivity averaged over the relevant interval of 0.5-40 kHz is $0.4~{{\mathrm {pm}}/{\sqrt {\mathrm {Hz}}}}$. As a result, at an 80-kHz sampling rate and with 2.7 mW IR power on the sample, the RMS noise of the OPL is 25 pm, while with a reduced power of 0.7 mW we reach 40 pm. The peak at 148 Hz is caused by the acoustic noise of ventilation in the building.

 figure: Fig. 4.

Fig. 4. Measured noise spectrum of our SD-OCT setup, obtained by Fourier transforming the measured OPL of a non-doped sphere, without excitation light.

Download Full Size | PDF

To verify photodetection-noise limited operation we show the sensitivity as a function of integration time in a semilogarithmic plot in Fig. 5 which has a slope of $-0.57$ close to the value of $-0.5$ characteristic of shot noise [24].

 figure: Fig. 5.

Fig. 5. Semilogarithmic plot of sensitivity of the $\delta \mathrm {OPL}$ versus integration time $\Delta t$. The power of SLD on the sample is 0.7 mW.

Download Full Size | PDF

4. Conclusion and discussion

We demonstrate a sensitive SD-OCT to measure sub-nanometer transient configurational alterations. We use a superluminescent light source with broad spectral bandwidth combined with a high-speed and high-resolution spectrometer. The system can detect changes in diameter of PS microspheres with a photodetection-noise limited sensitivity of $0.4~{{\mathrm {pm}}/{\sqrt {\mathrm {Hz}}}}$, which reveals the effects of thermal expansion and annealing in dye-doped microspheres caused by mild laser heating. This experiment paves the road to detecting sub-nm dimensional and sub-millisecond changes in cells and similarly sized objects, decoupled from their center of mass motion.

Funding

Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Neurophotonics program 6NEPH03, Vici 68047618).

Acknowledgments

The authors would like to thank Dave van den Heuvel, Aron Opheij, Dante Killian, and Paul Jurrius for technical support, and Marcel van der Heijden, Arnout Imhof, Dries van Oosten and Tayebeh Saghaei for discussions.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. S. Kwon, W. Yang, D. Moon, and K. S. Kim, “Comparison of cancer cell elasticity by cell type,” J. Cancer 11(18), 5403–5412 (2020). [CrossRef]  

2. K. Iida, R. Sakai, S. Yokoyama, N. Kobayashi, S. Togo, H. Y. Yoshikawa, A. Rawangkan, K. Namiki, and M. Suganuma, “Cell softening in malignant progression of human lung cancer cells by activation of receptor tyrosine kinase axl,” Sci. Rep. 7(1), 17770 (2017). [CrossRef]  

3. P. K. V. Babu, C. Rianna, U. Mirastschijski, and M. Radmacher, “Nano-mechanical mapping of interdependent cell and ecm mechanics by afm force spectroscopy,” Sci. Rep. 9(1), 1518 (2019). [CrossRef]  

4. J. M. Kralj, A. D. Douglass, D. R. Hochbaum, D. Maclaurin, and A. E. Cohen, “Optical recording of action potentials in mammalian neurons using a microbial rhodopsin,” Nat. Methods 9(1), 90–95 (2012). [CrossRef]  

5. Y. Park, C. Depeursinge, and G. Popescu, “Quantitative phase imaging in biomedicine,” Nat. Photonics 12(10), 578–589 (2018). [CrossRef]  

6. P. Hariharan, Optical Interferometry (Elsevier, 2003), 2nd ed.

7. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991). [CrossRef]  

8. W. Wang, S. Liu, Z. Qiu, M. He, L. Wang, Y. Li, and W. Huang, “Choroidal thickness in diabetes and diabetic retinopathy: a swept source oct study,” Invest. Ophthalmol. Vis. Sci. 61(4), 29 (2020). [CrossRef]  

9. S. Yu, B. E. Frueh, D. Steinmair, A. Ebneter, S. Wolf, M. S. Zinkernagel, and M. R. Munk, “Cataract significantly influences quantitative measurements on swept-source optical coherence tomography angiography imaging,” PLoS One 13(10), e0204501 (2018). [CrossRef]  

10. F. Zheng, Q. Zhang, Y. Shi, J. F. Russell, E. H. Motulsky, J. T. Banta, Z. Chu, H. Zhou, N. A. Patel, L. de Sisternes, M. K. Durbin, W. Feuer, G. Gregori, R. Wang, and P. J. Rosenfeld, “Age-dependent changes in the macular choriocapillaris of normal eyes imaged with swept-source optical coherence tomography angiography,” Am. J. Ophthalmol. 200, 110–122 (2019). [CrossRef]  

11. V. D. Nguyen, B. I. Akca, K. Wörhoff, R. M. De Ridder, M. Pollnau, T. G. van Leeuwen, and J. Kalkman, “Spectral domain optical coherence tomography imaging with an integrated optics spectrometer,” Opt. Lett. 36(7), 1293–1295 (2011). [CrossRef]  

12. J. van der Horst, A. K. Trull, and J. Kalkman, “Deep-tissue label-free quantitative optical tomography,” Optica 7(12), 1682–1689 (2020). [CrossRef]  

13. C. Yang, A. Wax, M. S. Hahn, K. Badizadegan, R. R. Dasari, and M. S. Feld, “Phase-referenced interferometer with subwavelength and subhertz sensitivity applied to the study of cell membrane dynamics,” Opt. Lett. 26(16), 1271–1273 (2001). [CrossRef]  

14. T. Akkin, C. Joo, and J. F. De Boer, “Depth-resolved measurement of transient structural changes during action potential propagation,” Biophys. J. 93(4), 1347–1353 (2007). [CrossRef]  

15. S. Batabyal, S. Satpathy, L. Bui, Y.-T. Kim, S. Mohanty, R. Bachoo, and D. P. Davé, “Label-free optical detection of action potential in mammalian neurons,” Biomed. Opt. Express 8(8), 3700–3713 (2017). [CrossRef]  

16. C. Fang-Yen, S. Oh, Y. Park, W. Choi, S. Song, H. S. Seung, R. R. Dasari, and M. S. Feld, “Imaging voltage-dependent cell motions with heterodyne mach-zehnder phase microscopy,” Opt. Lett. 32(11), 1572–1574 (2007). [CrossRef]  

17. Y. Zhang, B. Shen, T. Wu, J. Zhao, J. C. Jing, P. Wang, K. Sasaki-Capela, W. G. Dunphy, D. Garrett, K. Maslov, W. Wang, and L. V. Wang, “Ultrafast and hypersensitive phase imaging of propagating internodal current flows in myelinated axons and electromagnetic pulses in dielectrics,” Nat. Commun. 13(1), 5247 (2022). [CrossRef]  

18. S. Oh, C. Fang-Yen, W. Choi, Z. Yaqoob, D. Fu, Y. Park, R. R. Dassari, and M. S. Feld, “Label-free imaging of membrane potential using membrane electromotility,” Biophys. J. 103(1), 11–18 (2012). [CrossRef]  

19. S. Kasap, Optoelectronics & photonics: principles & practices: international edition (Pearson, 2013).

20. V. Ntziachristos, “Going deeper than microscopy: the optical imaging frontier in biology,” Nat. Methods 7(8), 603–614 (2010). [CrossRef]  

21. L. Godinho, “Imaging zebrafish development,” Cold Spring Harb Protoc. 2011(7), pdb.prot5647 (2011). [CrossRef]  

22. J. He, D. Mo, J. Chen, and L. Luo, “Combined whole-mount fluorescence in situ hybridization and antibody staining in zebrafish embryos and larvae,” Nat. Protoc. 15(10), 3361–3379 (2020). [CrossRef]  

23. H. Hu, “Cell-surface heparan sulfate is involved in the repulsive guidance activities of slit2 protein,” Nat. Neurosci. 4(7), 695–701 (2001). [CrossRef]  

24. M. A. Choma, A. K. Ellerbee, C. Yang, T. L. Creazzo, and J. A. Izatt, “Spectral-domain phase microscopy,” Opt. Lett. 30(10), 1162–1164 (2005). [CrossRef]  

25. J. W. Goodman, Statistical optics (John Wiley & Sons, 2015).

26. H. Sasabe and S. Saito, “Effects of temperature and pressure on the dielectric constant in non-polar polymers,” Polym. J. 3(6), 749–755 (1972). [CrossRef]  

27. H. F. Brinson and L. C. Brinson, Polymer engineering science and viscoelasticity (Springer, 2008), pp. 99–157, 2nd ed.

28. E. Riande, R. Diaz-Calleja, M. Prolongo, R. Masegosa, and C. Salom, “Flow behavior of polymer melts and solutions,” in Polymer Viscoelasticity: Stress and Strain in Practice, (Marcel Dekker, New York, 2000).

29. Y.-H. Lin, Polymer viscoelasticity: basics, molecular theories, experiments, and simulations (World Scientific, 2011).

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Setup of the SD-OCT system. (a) The light source and detection part that includes a superluminescent light source (SLD), a fiber coupler (FC) (two practical options 90:10 or 50:50), a collimator (CL1), and spectrometer (SP). (b) Microscopy section where the light is directed to the sample by a collimated lens (CL2), BS1 (P), microscope objective, and arrives at the microsphere. (c) excitation arm with a 532-nm laser, with collimation lens (CL3). (d) Schematic of the sample (e) Measured interferogram formed by the top and bottom reflections from the microsphere and (f) its intensity Fourier transform
Fig. 2.
Fig. 2. Optical path length of two microspheres exposed to various power, as example of two different microspheres behavior (a) Microsphere showing a low amount of inelastic creep. (b) Microsphere showing significant creep. Note the different scale of the y-axis. Raw data is shown in lighter color behind the smoothed data (moving average 50 data points)
Fig. 3.
Fig. 3. Comparison of OPL change in dye-doped and non-absorbing microspheres experiencing the same laser power of 3.12 mW.
Fig. 4.
Fig. 4. Measured noise spectrum of our SD-OCT setup, obtained by Fourier transforming the measured OPL of a non-doped sphere, without excitation light.
Fig. 5.
Fig. 5. Semilogarithmic plot of sensitivity of the $\delta \mathrm {OPL}$ versus integration time $\Delta t$. The power of SLD on the sample is 0.7 mW.

Equations (5)

Equations on this page are rendered with MathJax. Learn more.

I ( k ) = ( ρ / e )   s ( k ) δ k Δ t R r R s cos ( i 2 n k [ z 0 + δ z ( t ) ] )
I ( ± 2 n z 0 ) = ( ρ / 2 e ) S Δ t E ( 2 n z 0 ) R r R s exp [ ± i 2 k c n δ z ( t ) ] ,
δ φ ( t ) = I ( ± 2 n z 0 , t ) I ( ± 2 n z 0 , t 0 ) = 2 k c n δ z ,
δ O P L = 2 n δ z ( T ) + 2 z δ n ( T ) .
d δ O P L d T = 2 [ ( n 2 1 ) ( n 2 + 2 ) 2 n + n ] d δ z d T
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.