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

Random object optical field diagnostics by using carbon nanoparticles

Open Access Open Access

Abstract

We propose a new approach of using carbon nanoparticles for correlation optical diagnostics of а complex scalar optical field created by scattering and diffraction of radiation off a rough surface. This surface is simulated and we generate a diffraction pattern of the amplitude and phase distribution in the far field. Carbon nanoparticles of a certain size and concentration are obtained by the bottom-up methods of hydrothermal synthesis of citric acid and urea followed by centrifugation. The optical properties of carbon nanoparticles, such as luminescence and absorption in the visible spectrum that essentially differs for different wavelengths, as well as particle size of about dozen nanometers, are the determining criteria for using these particles as probes for the optical speckle field. Luminescence made it possible to register the coordinate position of carbon nanoparticles in real time. The algorithm for reconstruction of the scalar optical field intensity distribution through the analysis of the nanoparticle positions is here displayed. The skeleton of the optical speckle field is analyzed by Hilbert transform to restore the phase. Special attention is paid to the restoration of the speckle field’s phase singularities.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

A complex optical scalar field, here a fully developed speckle field, generally formed by coherent light, because of radiation scattering on a rough surface, is an erratic distribution of maxima and minima of intensity. The basic element of this field is an optical vortex, i.e. a phase singularity, an element of the field, where the phase is undefined. At these points, the intensity takes on a zero value, and the phase gradient (local wave vector) circulates around this point. Optical vortices [13], as a feature of an optical field structure, are an integral part of the optical field and exist independently of any special excitation conditions. As a rule, optical singularities appear near the lens focus, because of edge diffraction effects, or due to the superposition of plane waves.

 A number of research works have been devoted to the study of optical vortices, and, in general, the approach of singular optics is used [3,4]. Singularities are considered either in a 2D plane or in 3D space. The singularity representation differs in the way they are created as a point in a plane, or a complex dark line in 3D, which determines a certain skeleton of the wavefront. The way of creating the singularities will determine their visualization method.

The use of a reference beam is the traditional (classical) way for studying the optical field singularity for identifying its location and the skeleton (base) of the optical field, which allows us to identify the object. However, such a visualization method is not always applicable; any diffraction restriction reduces the probability of determining the localization of basic elements – here the singularity points. Here, the study of remote objects, which are the targets for searching the solutions of reverse problems, is complicated when an additional reference beam is needed for obtaining and recovering of the phase information from the recorded interference pattern. The diagnostics of remote objects is not possible in this case.

In addition, the use of the interference approach to visualize optical vortices such as interference forks is a rather difficult task, since the shift of the forks with respect to the center of the vortex is conceivable, which depends on the period of the interference pattern.

As pointed out in the review [4], both generating structured light, and its application for solving the practical problems for the singularity reproduction in the interference approach, is possible with a second beam containing an optical vortex at a nodal point. This allows one to observe rather small, specifically dark regions in the object under study having optical vortices.

To avoid the necessity of applying an additional second beam, the method of phase singularity reproduction by optical correlation technique [5,6] is used here. In general, this approach enables us to describe in statistical terms a surface with roughness, determining the transmission or reflection off the surface due to scattering of the radiation on the structure inhomogeneity. The separation of the real and imaginary parts of the complex field amplitude [6] allows one to calculate the phase of the field and to reconstruct a skeleton through the spatial distribution of the lines of phase singularities. Accordingly, the points of intersection of the lines of the real and imaginary parts of the complex field amplitude determine the position of the phase singularities. However, this empirical approach reproduces both the singularity points, and the points of the intensity minima without optical vortices, which reduces the accuracy of the performed diagnostics.

In our recent works [711], we proposed to use the two-dimensional discrete Hilbert transform to restore phase information, in compliance with the sign principle, by sequentially convolving the fixed amplitude distribution of the speckle field with the selected transformation core, facilitating an approach for the reconstruction of the phase object. The mathematical processing of the registered skeleton by the Hilbert transform results in a reconstructed phase map of the entire object.

The basis of phase singularities (the skeleton of the optical field) reproduces both the internal structure of the object and its spatial shape. The change of the basis in time makes it possible to follow any changes in the object in real time. The task of reproducing the phase singularity localization is formulated, both theoretically and experimentally, and there is a need to use an approach that would significantly reduce the error in identifying the position of singularities.

To restore the amplitude and phase information of the object and to expand the existing approach for optical field restoration, we use nanoscale test particles with special properties. Selection of usable luminophores is not easy. They should be quite small but larger than organic dye molecules, in order to provide the necessary diffusional pattern. In addition, they should be strongly absorbing in yellow-green wavelengths, brightly luminescing, and weakly absorbing for He-Ne laser irradiation wavelength. Essentially, their diffusional motion should be efficiently guided by the optical energy flow, and because of that they should be dipoles. Due to these reasons, for singular optical field diagnostics and the identification of singularity points’ position, we chose fluorescent carbon nanoparticles.

Fluorescent carbon nanostructures, as a new class of nanostructured materials, are of considerable interest due to their unique and exciting properties, such as tunable luminescence, easy synthesis, chemical inertness to surface modification, low cost, low toxicity and biocompatibility [12]. Essentially, in our knowledge, being of spherical structure of the order of ten nanometers, they possess optical symmetry that was discovered in fluorescence polarization measurements [13] and described on the level of single nanoparticles, presenting them as collective quantum emitters [14]. As it was suggested, the origin of this unique behavior is the existence of Frenkel-type molecular excitons [15]. Such properties of the new carbon materials lead to potential applications for creating light-emitting diodes [16], solar cells (absorbing light energy), sensors [17], biosensors for bio-visualization in medical diagnostics [18]. Optical trapping, as a new tool for manipulating and moving nanostructures, has been successfully applied to carbon nano-objects [19]. The very mechanism of spatial displacement of such nanostructures is determined by the interaction of the laser field with instantly induced dipole moments of individual nanoparticles [14].

The peculiar absorption and emission properties of carbon nanoparticles, which determines the nature of optical capture, promote the use of these particles for analysis and diagnostics of complex optical fields, including speckle fields, assuming that the particles are localized not just in the phase gradient regions, but at the singularity points of the studied optical field. The intensity minima that do not contain singularities are characterized by smaller intensity gradients, the energy flows here being much weaker than in the singularity regions. Accordingly, the number and concentration of nanoparticles will increase significantly in the minimum intensity regions with singularities.

This paper demonstrates the possibility of amplitude-phase optical field reconstruction for a random scattering object over a field skeleton visualization by observing spatially redistributed carbon nanoparticles in the points of minimum intensity with singularities due to internal optical flow action. The full intensity distribution is restored from the obtained singularity points, which is subjected to Hilbert processing. The result of this action is a phase map of the object under study. Processing the field of optical radiation scattered by an object is recorded in real time. The proposed approaches here are expected to have many applications. They may be applied in the study of spatial structure of the field of a sharply focused beam [20], where the use of interference methods is impossible; in the analysis of dynamic liquid media, including the restoration of the size distribution function of micro and nanoparticles in dynamic light scattering technologies [21], in life sciences; in the study of non-stationary processes in cells; in the study of turbulent gaseous media [22], in aerodynamic applications, as well as for the solution of problems of digital holographic interference in the analysis of non-stationary objects and scenes [23,24].

2. Study of an optical speckle field by correlation-optical method using carbon nanoparticles

In this work, to diagnose the optical scalar field (speckle field) obtained by scattering and diffraction of radiation off a rough surface, carbon nanoparticles were applied. The reason for their choice was explained above. In our case, they were prepared using a one-step hydrothermal synthesis from citric acid and urea in autoclaves with following centrifugation at t = 190°C for 2 hours, which makes it quite easy to obtain bright carbon nanoparticles with stabilized carboxy-, amino- and other groups with a diameter from 2-3 nm to hundred nanometers [2527]. In order to verify, analyze, and estimate the size of nanoparticles, atomic force microscopy measurements were performed. We use water solution of prepared carbon particles. The used technological process for the production of carbon nanoparticles ensured a high stability of the formed nanostructures; their parameters did not change during the measurements. The entire technological process of manufacturing carbon nanoparticles, approaches to assessing the size and optical characteristics, which is actually necessary for studying speckle fields, are given in our paper [27].

The evaluation of the absorption and luminescence spectra of particles obtained during our experiment sets the absorption maximum at 405 nm, the luminescence maximum at a wavelength of about 450 nm, and a weak absorption at the wavelengths of He-Ne laser radiation (633 nm), which is crucial for using these particles as a speckle field skeleton probe. Absorption rate (A) of carbon the prepared nanoparticles [27] at He-Ne laser irradiation, in the field of which the speckle field under study is formed, is minimal and is in the order of 0.2${\; }{\textrm{cm}^{ - 1}}$. The transmission coefficient (τ$\; = \textrm{exp}({ - D} )$) was evaluated at the given experimental conditions, and was in the order of 80%. Determination of these optical parameters was made by estimating the optical density (D) of the used solutions.

For these carbon nanoparticles, sized from ten to hundred nanometers, a size-dependent effect [19] is not observed, which greatly simplifies the practical implementation of the experiment to detect luminescence of carbon nanoparticles followed by particle localization in the zones of an optical singularity.

Laser radiation acting on carbon nanoparticles can cause heating of the medium surrounding the particles and of the particles themselves, which force us to discuss the temperature effects. A consequence of this effect can be the redistribution of the optical flows, which again can cause a change the direction of the particle motion and eventually their capture. It is possible to quench fluorescence at temperatures above 50°, which is accompanied by increased diffusion, phonon relaxation of the shell excitation energy, consisting of active polar groups. As a result, the emission intensity decreases [28]. In order to avoid the temperature effects including destroyed changes or possible particle oscillation due to their Brownian motion leading to the essential change of particle position, the temperature regime and the exposure time are controlled. Therefore, we use laser sources with a maximum wavelength of 405 nm and a power of 5 mW for the observation of luminescence, and a He-Ne Laser for obtaining the scattered field in which the motion of carbon particles occur. Time necessary for particle redistribution and capture is about 30 seconds, which is confirmed theoretically.

The exceptional optical properties of carbon nanoparticles, specifically their interaction with external laser radiation, allow one to use them for diagnostics of optical fields, and consequently to recover the field skeleton under study, and the information retrieval of the object.

In order to study the optical fields by correlation-optical methods, we address the following problems:

  • - simulation of an arbitrary surface with given roughness;
  • - obtaining a diffraction pattern of the simulated object with corresponding speckle field in the far zone;
  • - visualization of singular points due to the absorption of internal optical flows by carbon nanoparticles with their redistribution into the regions of singularity (luminescence of particles is registered);
  • - object phase retrieval by the use of 2-D Hilbert transform without reference wave.

Thus, to study optical field carbon nanoparticles of size about $\lambda /10$ with strong absorption in the yellow-green region (405 nm), luminescence at wavelengths of 530 nm and weak absorption at the red wavelength of speckle field formation (633 nm) are used. The particles of this size are easily controlled by the speckle field by moving along the intensity gradient [29,30]. The structure of the optical speckle field is visualized and captured by a CCD camera.

Consequently, the speckle field is restored by a set of spatially redistributed carbon nanoparticles in the areas of minimum intensity with singularity, thereby reproducing the structure base (skeleton) of the optical field and the phase map in real time.

2.1. Rough surface modeling [6]

Here, it is proposed to use the simulation approach of a surface with given roughness and the results of the experiment, which were tested by the authors and previously presented [31]. There a multifunctional system for the diagnostics of rough surfaces of various structure is presented. Modeling of a rough surface was carried out according to the appropriate algorithm [31]. Simulation took place in the MATHEMATICA software package. We use the MATHEMATICA program to simulate a speckle field in the present work. A surface of the field with a square of 100 by 100 µm2 is derived; the number of pixels in a given area is determined by a step length of $\lambda /3$. The choice of step length is determined by:

  • i) there should be no mutual amplification or quenching of the signal from the neighboring pixels (the wavelength is divided by an odd number),
  • ii) the obtained points should determine the pixels, in accordance with the resolution of the camera, used to record the image (1024 pixels).

During the simulation, the points with a step of about of 5 µm are considered and the height of the inhomogeneity is randomly generated within the interval 0 - hmax (hmax = 2 µm). The height of the surface inhomogeneity is obtained using spline interpolation of the generated points. Intermediate values of heights were obtained by the interpolation of periodically located random heights. At this stage of modeling, the standard MATHEMATICA function - ListInterpolation was used in the context of the 1st order linear B-spline interpolation. To obtain a complete picture of the simulated surface roughness, a Gaussian filter with a radius of 3 µm was used, with smoothing by 30 points, giving the possibility to reproduce the relief of the rough surface (Fig. 1). The choice of parameters used in the simulation is carried out through a systematic analysis of their values, followed by matching the simulation results with an experimentally recorded rough surface profile [1,6,31]. The resulting smoothed distribution of heights was the object of our detailed analysis.

 figure: Fig. 1.

Fig. 1. Two-dimensional (a) and three-dimensional (b) distributions of the heights of the object surface.

Download Full Size | PDF

Figure 1 shows the two-dimensional and three-dimensional coordinate distributions of the surface inhomogeneities obtained in the modeling process.

2.2 Obtaining a diffraction image of a simulated surface

To simulate the diffraction image of a surface with a given roughness, the diffraction integral in the Rayleigh-Sommerfeld approximation [6] is used, and in the discrete approximation, the double integral turns into a double sum:

$$U({\xi ,\zeta } )= \frac{z}{{i\lambda }}\sum\limits_{x = 1}^X {\sum\limits_{y = 1}^Y {\frac{{F({x,y} )}}{{{R^2}(x,y,z,\xi ,\zeta )}}} } \exp \{{ - ik[{R(x,y,z,\xi ,\zeta ) + 2h(x,y)} ]} \},$$
where: $F(x,y)$– is the aperture function, which determines the amplitude transmission of a rough surface; $h(x,y)$ describes the relief of the analyzed surface, and determines the phase delay of the beams diffracted on the surface $\varphi (x,y) = 2kh(x,y)$, $k = {{2\pi } / \lambda }$ is the wave number; $\lambda $ – wavelength. $R(x,y,z,\xi ,\zeta ) = \sqrt {{z^2} + {{({x - \xi } )}^2} + {{(y - \zeta )}^2}} $ is the distance between a point at the object and a point in the observation plane; z is the distance between the plane where the object is localized and the observation plane; $x,y,\xi ,\zeta $ – rectangular coordinates in the plane of the object and in the plane of observation, respectively. In our case z = 1000 µm = 10a, where а – is a linear size of the object. Here X, Y are the number of points of the analyzed object selected in x and y directions, respectively, with a step size of $\lambda /3$.

Determining the real $\textrm{Re} [{U({\xi ,\zeta } )} ]$ and imaginary ${\mathop{\rm Im}\nolimits} [{U({\xi ,\zeta } )} ]$ parts of the calculated field makes it possible to evaluate the field amplitude modulus, intensity, and phase. In the far field, a simulated diffraction pattern of 240 by 240 µm2 (400 × 400 pixels) and the corresponding phase map are shown in Fig. 2. The field is based on the radiation from a He-Ne laser (λ=633 nm).

 figure: Fig. 2.

Fig. 2. Simulated diffraction pattern (a) and calculated phase map (b). The white square shows the part of the pattern undergoing further analysis: 1 – area of 30 × 30 µm2, 2 - area of 60 × 60 µm2.

Download Full Size | PDF

The area of the diffraction pattern of 30 × 30 µm2 (60 × 60 µm2) marked by white squares in Fig. 2 was used to demonstrate the effectiveness of the given model to reconstruct the optical field intensity distribution by a coordinate-wise analysis of the tracked motion of carbon nanoparticles until they are eventually captured by regions of intensity minima having singularities. For a better assessment, the simulation results obtained for the areas of different spatial scales (Fig. 2, squares 1,2). Marked area of about 30 × 30 µm2 is later used for phase information retrieval using Hilbert analysis.

2.3. Reconstructing of the intensity distribution (computer simulation)

The next step is the search for ways to recover the intensity distribution through processing of the recorded tracks of the carbon nanoparticles’ motion in the complex optical field, which is obtained by the light scattered off the rough surface, e.g. the diffraction pattern described above (Fig. 2, selected parts)).

A cuvette with a water solution of carbon nanoparticles of size $\lambda /10$, with known properties is placed in the optical field. Randomly suspended particles move in this field because of the internal optical flows, and concentrate in regions of minimum intensity with singularities. It is then possible to reconstruct the tracks of nanoparticle motion and the intensity distribution. Luminescence reveals the motion of particles.

The measured value of the absorption rate (A = 0.2${\; }{\textrm{cm}^{ - 1}}$) at λ = 633 nm, allows one, using the Fresnel theory, to calculate the imaginary part of the refractive index (n), i.e. the absorption coefficient, as k = (A · λ) / 4π [32,33].

According to the law of energy conservation, the reflection coefficient $\rho = 1 - k - \tau $, made it possible to determine the real part of the refractive index (n) from the relation $\rho = \frac{{{{(n - 1)}^2} + {k^2}}}{{{{(n + 1)}^2} + {k^2}}},$ which is sufficient to use for strongly absorbing elements, carbon nanoparticles, and in order to simplify the calculated ratios.

Using values for n and k, the real and imaginary parts of the relative permittivity are calculated as, ${\varepsilon _r} = {n^2} - {k^2}\textrm{ and }{\varepsilon _i} = 2n\,k.$ Particle polarizability $\alpha $ (in the Gauss system) is:

$$\begin{array}{l} \textrm{Re} [\alpha ]= \alpha ^{\prime} = {r^3}\frac{{({{\varepsilon_r} - 1} )({\varepsilon _r} + 2) + \varepsilon _i^2}}{{{{({\varepsilon _r} + 2)}^2} + \varepsilon _i^2}},\\ {\mathop{\rm Im}\nolimits} [\alpha ]= \alpha ^{\prime\prime} = {r^3}\frac{{3\,{\varepsilon _i}}}{{{{({\varepsilon _r} + 2)}^2} + \varepsilon _i^2}}. \end{array}$$

We use the obtained optical parameters of the applied carbon nanoparticles for the calculation of the optical force that determines the action due to the internal energy flows.

The presented nanoparticles are described by the Rayleigh light scattering mechanism [3437], which is combined with Mie theory, giving each component of the resulting optical force acting on the particle - gradient component ${\vec{F}_{grad}}$, absorption ${\vec{F}_{abs}}\; $and scattering ${\vec{F}_{scatt}}\; $components [2730]:

$${\vec{F}_{grad}} ={-} \frac{{\alpha ^{\prime}}}{2}n\nabla {E^2},\; \; {\vec{F}_{abs}} = \frac{{n\vec{S}{C_{abs}}}}{c},\; \; {\vec{F}_{scatt}} = \frac{{n\vec{S}{C_{scatt}}}}{c},$$
where ${C_{scatt}} = \frac{{{k^4}{{|\alpha |}^2}}}{{4\pi }}$, ${C_{abs}} = k\alpha ^{\prime\prime}\; $ are the scattering and absorption cross section (in a Gaussian system of units), $k = \frac{{2\pi }}{\lambda }$ is the wavenumber.
$$\textrm{Here }{E^2} = \frac{{{{|E |}^2}}}{2},\;\textrm{a}\;\vec{S} \equiv I = \frac{c}{{8\pi }}{|E |^2}.$$

Figure 3 shows the ratio of the components of the optical force for different size of carbon nanoparticles at an observation wavelength for nanoparticles in a speckle field.

 figure: Fig. 3.

Fig. 3. Ratio of the optical force components as a function of the size of carbon nanoparticles.

Download Full Size | PDF

The distribution of the gradient force is analyzed in the x(y) directions, i.e. in the cross section of the beam. The scattering and absorbing components exert their effect in the direction of the action of the Pointing vector (longitudinal z-direction). Therefore, when analyzing the distribution and influence of these forces in the transverse plane, the action of these components does not appear. An analysis of the magnitude of the active component of the gradient force shows that the minimum regions with singularities determine a force value that can change the particles’ position in time.

Naturally, in the case of three-dimensional visualization of the complex dynamics of particle motion, all three components need to be taken into account. An analysis of the components of the optical force demonstrates that the value of the gradient component of the optical field far exceeds the absorbing component, i.e. temperature effects due to absorption of radiation in the medium will be absent.

At the same time, the value of the scattering component is such that it reveals the displacement of the nanoparticle in the longitudinal direction. The nanoparticle motion will be observed until they are captured by the intensity minimum. Preliminary modeling showed that the influence of optical energy flows on nanoparticles is best manifested for particles with a size of $\lambda /10$ which allows us to use them to restore the object’s phase information.

The equation of motion of the i-particle (i=1..N, where N is the total number of nanoparticles analyzed in the speckle field) under the action of optical force becomes:

$${m_i}\frac{{d\; {{\vec{v}}_i}}}{{dt}} = {\vec{F}_{op{t_i}}} + {\vec{F}_{s{t_i}}}, $$
where ${\vec{F}_{op{t_i}}} = {\vec{F}_{gra{d_i}}} + {\vec{F}_{ab{s_i}}} + \; {\vec{F}_{scat{t_i}}}$ – the resulting optical force, the components of which are set according to (1), (2):
  • ${\vec{F}_{s{t_i}}} = 6\pi {r_i}\eta {\vec{v}_i}$Z – Stokes force,
  • ${m_i} = \frac{4}{3}\pi {r_i}^3\rho $Z – mass of i-carbon nanoparticle, here ${r_i}$i-particle radius,
  • $\rho = 2.26\; g/\textrm{s}{\textrm{m}^3}$Z - water density [38],
  • $\eta = 8,9 \cdot {10^{ - 4}}\; Pa \cdot \textrm{sec}$Z – dynamic viscosity of the medium (water) at T=25°C [39].

In the scalar case, taking into account the Cartesian x(y) components, the equation of motion has the form $m\frac{{d{v_{x,y}}_i}}{{dt}} = {F_{op{t_{x,y}}_i}} - 6\pi r\eta {v_{x,y}}_i$. Accordingly, the acceleration that carbon particles acquire during their motion is ${a_{x,y}}_i = \frac{{d{v_{x,y}}_i}}{{dt}}$, which allows us to determine the coordinates of the particle at a given time t. The initial moment of observation is taken as the zero value.

$${x_i}(t )= {x_i}({{t_0}} )+ {v_{{x_i}}}t + \frac{{{a_{{x_i}}}{t^2}}}{2},\;{y_i}(t )= {y_i}({{t_0}} )+ {v_{{y_i}}}t + \frac{{{a_{{y_i}}}{t^2}}}{2}.$$
A computer simulation is performed (Fig. 4), where time t - are running values for each nanoparticle, given by the observation interval, which is determined by the time of capture of the nanoparticle by the intensity minimum with singularities in the area of 30 × 30 µm2 (Fig. 5). Figure 4 demonstrates the changes in carbon nanoparticles position in time: a) the particles are randomly localized in a scalar complex optical field, b) based on the calculated coordinates we determine the track length (Δs) that the nanoparticle passes over a selected period before its capture, c) the calculated final position of the carbon nanoparticles.

 figure: Fig. 4.

Fig. 4. The particles’ position in time: the initial random position of carbon nanoparticles (a); the tracks of nanoparticles during their motion (b); the position of the particles in the speckle field after their redistribution into the areas of minimum intensity (observation time 5 second) (c).

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Speckle field (a, b, c) with gradient lines of intensity (white lines with arrows) determining the motion of carbon nanoparticles (green tracks) (c) into the areas with singularities (red point) (a,b). The size of analyzed optical field is 1.8 × 1.8 µm2 (red square pointed in Fig. 5(a), (b), (c)).

Download Full Size | PDF

The found values of, ${\alpha _{x,{y_i}}},{\nu _{x,{y_i}}}$, i.e. the solution of the equation of motion (3), where active components of the optical force (1) are taken into account, allows one to restore the distribution of the optical field intensity over the known coordinate distribution of nanoparticle motion (Figs. 4 and 5). In this case, the local extrema of the functional dependence of the intensity on the coordinate are sought for restoring both the regions of the minimum, maximum intensity, and its saddle points.

Figures 5(a) and (b) presents the distribution of intensity gradient lines (white lines with arrows) in the simulated speckle field (the size of the analyzed region is 1.8 × 1.8 µm2 (Figs. 5(a) (red square), (b), and (c)). The singularity point is marked by a red point (Figs. 5(a) and (b)). In 30 second (Fig. 5(c)) all carbon nanoparticles redistribute into the minimum area with singularities moving along the gradient lines of intensity (green tracks).

Figure 6 demonstrates the position of the particles in the intensity area of 60 × 60 µm2 at 0.375 sec. The video presents the motion of particles during the period of observation (Visualization 1) until their capture in the intensity minima. The observation time here is 5 second.

 figure: Fig. 6.

Fig. 6. Demonstration of particles’ motion in time (see Visualization 1)

Download Full Size | PDF

The resulting coordinate set of intensity values, which is generally random, and defines an unstructured data set, is converted into a regular continuous array by using a linear B-spline interpolation of the 1st order.

Further, by mathematical processing, from a continuous distribution, a discrete periodic distribution of intensity points is extracted, and using the higher order interpolation, the intensity distribution is refined (Fig. 7(b)).

 figure: Fig. 7.

Fig. 7. Original optical field (a) (30 × 30 µm2), (c) (60 × 60 µm2) and corresponding reconstructed optical field (b), (d) through analyzing coordinate distribution of the nanoparticle motion tracks.

Download Full Size | PDF

3. Phase restoring by using a Hilbert transform

The final step of our investigation is to extract the phase information from the reconstructed intensity distribution of the field (area of about 30 × 30 µm2 was analyzed). MATHEMATICA software package allows us to restore the phase by applying a Hilbert filter [40,41] to the intensity distribution of the speckle field under study. The corresponding phase map is presented in Fig. 8(a). As it is shown by the simulation results, singularity information in Hilbert transform action is lost. Here, the points of phase singularities, marked in the picture by red points, are determined by the final position of carbon nanoparticles in the optical field (Fig. 8(b)). The same points are obtained through the intersection of real and imaginary parts of the complex field amplitude (Fig. 8(c), green and blue lines). Red points in white squares determine the location of intensity minimum without singularities.

 figure: Fig. 8.

Fig. 8. Phase information: a) the restored phase map by applying of Hilbert transform, b) phase map with phase singularities (red points), identification of which is a result of the redistribution of carbon nanoparticles in the optical field; c) phase singularities obtained at the intersection of real (green lines) and imaginary (blue lines). Red points in white squares indicate the position of an intensity minimum without a singularity.

Download Full Size | PDF

The error in the results of reconstructing the phase using traditional mathematical methods of phase recalculating through the values of the real and imaginary parts of the complex field amplitude, described in 2.2, and using a Hilbert transform, proposed in MATHEMATICA, is about 15%. The estimation of this error was carried out pixel by pixel, by determining the deviation of the phase on the phase map reconstructed using the Hilbert transform from the values obtained using traditional mathematical methods.

An increase of the probability of the retrieval of amplitude zeros (phase singularities) is possible by the approach proposed in [4246]. Here, the zeros of $U({\xi ,\eta } )$ are transformed into branch points of $\log [{U({\xi ,\eta } )} ]$ by applying the function $\log [{U({\xi ,\eta } )} ]= \log |{U({\xi ,\eta } )} |+ i\varphi ({\xi ,\eta } ).$ Such an approximation allows one to use the two-dimensional logarithmic Hilbert transform that involves extracting the zero information in the complex lower half-plane, in which a real variable is changed into a complex one. At this stage of data processing, information about zeros is also lost. Further, applying of exponential filters $\exp [{ - 2\pi \,{y_c}\,\xi } ]$, $\exp [{ - 2\pi \,{y_c}\,\eta } ]$, ${y_c} \gt 0$ (${y_c}$– interpolation constant) in the object plane (${\xi ,\eta }$) makes it possible to leave the zero information in the complex lower half-plane. This revises the obtained Hilbert phase on each line parallel to the ${x_1},{x_2}$ axis in the Fourier transform plane, as described in detail in [42]. In this case, the accuracy of reproducing the phase map with singularities of the optical field increases to 92%.

As a conclusion, it can be noted that the approach used to reconstruct the amplitude and phase information of an optical field based on stationary scattering object opens up additional possibilities both for non-stationary scattering object diagnostics and macro-analysis of the shape of the remote object for solving the aerodynamics and space problems as well as for micro-diagnostics of the image in non-stationary processes in cells in the life sciences.

4. Conclusion

A new approach of using carbon nanoparticle to study the optical field obtained because of diffraction and scattering on a stationary object in real time is demonstrated. Сarbon nanoparticles of size about $\lambda /10$, characterized by strong absorption and luminescence at yellow-green wavelengths and weak absorption at He-Ne wavelength were used for the first time for probing a speckle field. Visualization of singularity points by these nanoparticles, and their motion in the optical field due to the action of internal optical flows, makes it possible to reconstruct the intensity distribution of an optical field in real time. Applying Hilbert transform for phase retrieval of the optical field, describing the random light scattering object, provides the information about the phase map of the optical field. The error in reconstructing of the phase map using both the Hilbert transform and traditional approaches for estimating the components of the complex amplitude does not exceed 15%. Obtained results open up new possibilities in the diagnostics of non-stationary objects in liquid and turbulent media, as well as for study the non-stationary processes in cells in the life sciences.

Acknowledgments

The authors thank Prof. Alexander P. Demchenko for helpful assistance and discussions

Disclosures

The authors declare no conflicts of interest.

References

1. Introduction to Singular Correlation Optics, Oleg V. Angelsky, ed. (SPIE Press Book, 2019), 252 p.

2. M. S. Soskin and M. V. Vasnetsov, Singular Optics: Progress in Optics, E. Wolf, ed., (1st Edition – Elsevier, 2001), pp. 219–276

3. I. I. Mokhun, Introduction to Linear Singular Optics: Optical Correlation: Techniques and Applications, O. Angelsky, ed. (SPIE Press, Bellingham, 2007), pp. 1–131.

4. H. Rubinsztein-Dunlop, A. Forbes, M. V. Berry, M. R. Dennis, D. L. Andrews, M. Mansuripur, C. Denz, C. Alpmann, P. Banzer, T. Bauer, E. Karimi, L. Marrucci, M. Padgett, M. Ritsch-Marte, N. M. Litchinitser, N. P. Bigelow, C. Rosales-Guzmáán, A. Belmonte, J. P. Torres, T. W. Neely, M. Baker, R. Gordon, A. B. Stilgoe, J. Romero, A. G. White, R. Fickler, A. E. Willner, G. Xie, B. McMorran, and A. M. Weiner, “Roadmap on structured light,” J. Opt. 19(1), 013001 (2017). [CrossRef]  

5. O. V. Angelsky, M. P. Gorsky, S. G. Hanson, V. P. Lukin, I. I. Mokhun, P. V. Polyanskii, and P. A. Ryabiy, “Optical correlation algorithm for reconstructing phase skeleton of complex optical fields for solving the phase problem,” Opt. Express 22(5), 6186–6193 (2014). [CrossRef]  

6. O. V. Angelsky, P. P. Maksimyak, and S. G. Hanson. “The Use of Optical-Correlation Techniques for Characterizing Scattering Object and Media.” Bellingham: SPIE Press PM71 1999.

7. C. Y. Zenkova, M. P. Gorsky, and P. A. Ryabiy, “Different approaches to phase restoration of distant complex optical fields,” Opt. Appl. 45(2), 139–150 (2015). [CrossRef]  

8. C. Y. Zenkova, M. P. Gorsky, and P. A. Ryabiy, “Phase retrieval of speckle fields based on 2D Hilbert transform,” Opt. Mem. Neural Networks 24(4), 303–308 (2015). [CrossRef]  

9. C. Y. Zenkova, M. P. Gorsky, and P. A. Ryabiy, “Methods of restoring spatial phase distribution of complex optical fields in the approximation of singular optics,” Romanian Rep. in Phys. 67(4), 1401–1411 (2015).

10. C. Y. Zenkova, M. P. Gorsky, P. A. Ryabiy, and A. O. Angelskaya, “Additional approaches to solving the phase problem in optics,” Appl. Opt. 55(12), B78–B84 (2016). [CrossRef]  

11. C. Y. Zenkova, M. P. Gorsky, and P. A. Ryabiy, “Pseudo-phase mapping of speckle fields using 2D Hilbert transformation,” Opt. Appl. 46(1), 153–162 (2016). [CrossRef]  

12. Y-P. Sun, Carbon Dots, (Springer, 2020), 203 p.

13. M. O. Dekaliuk, O. Viagin, Y. V. Malyukin, and A. P. Demchenko, “Fluorescent carbon nanomaterials:“quantum dots” or nanoclusters?” Phys. Chem. Chem. Phys. 16(30), 16075–16084 (2014). [CrossRef]  

14. S. Ghosh, A. M. Chizhik, N. Karedla Mariia, O. Dekaliuk, I. Gregor, H. Schuhmann, M. Seibt, K. Bodensiek, I. A. T. Schaap, O. Schulz, A. P. Demchenko, J. Enderlein, and A. I. Chizhik, “Photoluminescence of carbon nanodots: Dipole emission centers and electron–phonon coupling,” Nano Lett. 14(10), 5656–5661 (2014). [CrossRef]  

15. A. P. Demchenko, “Excitons in Carbonic Nanostructures,” C 5(4), 71 (2019). [CrossRef]  

16. E. Liu, D. Li, X. Zhou, G. Zhou, H. Xiao, D. Zhou, P. Tian, R. Guo, and S. Qu, “Highly emissive carbon dots in solid state and their applications in light-emitting devices and visible light communication,” ACS Sustainable Chem. Eng. 7(10), 9301–9308 (2019). [CrossRef]  

17. J. B. Essner and G. A. Baker, “The emerging roles of carbon dots in solar photovoltaics: a critical review,” Environ. Sci.: Nano 4(6), 1216–1263 (2017). [CrossRef]  

18. K. Hola, Y. Zhang, Y. Wang, E. P. Giannelis, R. Zboril, and A. L. Rogach, “Carbon dots—Emerging light emitters for bioimaging, cancer therapy and optoelectronics,” Nano Today 9(5), 590–603 (2014). [CrossRef]  

19. Y. Sun, P. Wang, Z. Lu, F. Yang, M. Meziani, G. LeCroy, Y. Liu, and H. Qian, “Host-guest carbon dots for enhanced optical properties and beyond,” Sci. Rep. 5(1), 12354 (2015). [CrossRef]  

20. T. G. Brown, M. A. Alonso, A. Vella, M. J. Theisen, S. T. Head, S. R. Gillmer, and J. D. Ellis, “Focused beam scatterometry for deep subwavelength metrology,” Proc. SPIE 8949, 89490Y (2014). [CrossRef]  

21. J. Stetefeld, S. A. McKenna, and T. R. Patel, “Dynamic light scattering: a practical guide and applications in biomedical sciences,” Biophys. Rev. 8(4), 409–427 (2016). [CrossRef]  

22. V. P. Lukin and B. V. Fortes, “Phase-correction of turbulent distortions of an optical wave propagating under conditions of strong intensity fluctuations,” Appl. Opt. 41(27), 5616 (2002). [CrossRef]  

23. T. Latychevskaia, “Iterative phase retrieval for digital holography: tutorial,” J. Opt. Soc. Am. A 36(12), D31–D40 (2019). [CrossRef]  

24. Y. Tanaka, S. Tani, and S. Murata, “Phase retrieval method for digital holography with two cameras in particle measurement,” Opt. Express 24(22), 25233 (2016). [CrossRef]  

25. X. Liu, H.-B. Li, L. Shi, X. Meng, Y. Wang, X. Chen, H. Xu, W. Zhang, X. Fang, and T. Ding, “Structure and photoluminescence evolution of nanodots during pyrolysis of citric acid: from molecular nanoclusters to carbogenic nanoparticles,” J. Phys. Chem. C 5(39), 10302–10312 (2017). [CrossRef]  

26. J. Schneider, C. J. Reckmeier, Y. Xiong, M. von Seckendorff, A. S. Susha, P. Kasák, and A. L. Rogach, “Molecular Fluorescence in Citric Acid-Based Carbon Dots,” J. Phys. Chem. C 121(3), 2014–2022 (2017). [CrossRef]  

27. P. P. Maksimyak, C. Y. Zenkova, and V. M. Tkachuk, “Carbon Nanoparticles. Production, properties, perspectives of use,” Phys. Chem. Solid State 21(1), 13–18 (2020). [CrossRef]  

28. C. J. Reckmeier, J. Schneider, A. S. Susha, and A. L. Rogach, “Luminescent colloidal carbon dots: optical properties and effects of doping,” Opt. Express 24(2), A312–A340 (2016). [CrossRef]  

29. V. Angelsky, M. P. Gorsky, P. P. Maksimyak, A. P. Maksimyak, S. G. Hanson, and C. Y. Zenkova, “Investigation of optical currents in coherent and partially coherent vector fields,” Opt. Express 19(2), 660–672 (2011). [CrossRef]  

30. C. Y. Zenkova, M. P. Gorsky, P. P. Maksimyak, and A. P. Maksimyak, “Optical currents in vector fields,” Appl. Opt. 50(8), 1105–1112 (2011). [CrossRef]  

31. O. V. Angelsky, P. P. Maksimyak, V. V. Ryukhtin, and S. G. Hanson, “New feasibilities for characterizing rough surfaces by optical-correlation techniques,” Appl. Opt. 40(31), 5693–5707 (2001). [CrossRef]  

32. M. Fox, Optical properties of solids (Oxford University Press Inc, New York, 2001), pp. 1–8

33. E. M. Gullikson, “13. Optical Properties of Materials,” Exp. Methods Phys. Sci. 31, 257–270 (1998). [CrossRef]  

34. M. Dienerowitz, M. Mazilu, P. J. Reece, T. F. Krauss, and K. Dholakia, “Optical vortex trap for resonant confinement of metal nanoparticles,” Opt. Express 16(7), 4991–4999 (2008). [CrossRef]  

35. C. Y. Zenkova, M. P. Gorsky, I. V. Soltys, and P. O. Angelsky, “On the possibilities of using inhomogeneity in light energy distribution for estimating the degree of coherence of superposing waves,” Appl. Opt. 51(10), C38–C43 (2012). [CrossRef]  

36. O. V. Angelskii, O. G. Ushenko, D. N. Burkovets, O. D. Arkhelyuk, and Y. A. Ushenko, “Polarization-correlation studies of multifractal structures in biotissues and diagnostics of their pathologic changes,” Laser phys.-Lawrence 10(5), 1136–1142 (2000).

37. C. Zenkova, I. Soltys, and P. Angelsky, “The use of motion peculiarities of particles of the Rayleigh light scattering mechanism for defining the coherence properties of optical fields,” Opt. Appl. 43(2), 297–312 (2013). [CrossRef]  

38. C. Y. Zenkova, “Interconnection of polarization properties and coherence of optical fields,” Appl. Opt. 53(10), B43–B52 (2014). [CrossRef]  

39. https://www.americanelements.com/carbon-nanoparticles-7440-44-0

40. http://hyperphysics.phy-astr.gsu.edu/hbase/Tables/viscosity.html

41. O. Sheehan, “Computing the Hilbert transform and its inverse,” Math. Comp. 80(275), 1745–1767 (2011). [CrossRef]  

42. J. K. Costain and C. Coruh, “Basic Theory in Reflection Seismology: with MATHEMATICA Notebooks and Examples” on CD-ROM, Elsevier, 2004, 576 p.

43. N. Nakajima and T. Asakura, “Two dimensional phase retrieval using the logarithmic Hilbert transform and the estimation technique of zero information,” J. Phys. D: Appl. Phys. 19(3), 319–331 (1986). [CrossRef]  

44. R. E. Burge, M. A. Fiddy, A. H. Greenaway, and G. Ross, “The phase problem,” Proc. R. Soc. London. 350, 191–212 (1976).

45. G. Ross, M. A. Fiddy, M. Nieto-Vesperinas, and M. W. L. Wheeler, “The Phase Problem in Scattering Phenomenon: The Zeros of Entire Functions and their Significance,” Proc. R. Soc. Lond. A 360(1700), 25–45 (1978). [CrossRef]  

46. Inverse Source Problems in Optics, H.P. Baltes, ed. (Publisher Springer-Verlag, Berlin Heidelberg, 1978), 204p.

Supplementary Material (1)

NameDescription
Visualization 1       An added video has been uploaded. This was specifically asked for by one of the reviewers.

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 (8)

Fig. 1.
Fig. 1. Two-dimensional (a) and three-dimensional (b) distributions of the heights of the object surface.
Fig. 2.
Fig. 2. Simulated diffraction pattern (a) and calculated phase map (b). The white square shows the part of the pattern undergoing further analysis: 1 – area of 30 × 30 µm2, 2 - area of 60 × 60 µm2.
Fig. 3.
Fig. 3. Ratio of the optical force components as a function of the size of carbon nanoparticles.
Fig. 4.
Fig. 4. The particles’ position in time: the initial random position of carbon nanoparticles (a); the tracks of nanoparticles during their motion (b); the position of the particles in the speckle field after their redistribution into the areas of minimum intensity (observation time 5 second) (c).
Fig. 5.
Fig. 5. Speckle field (a, b, c) with gradient lines of intensity (white lines with arrows) determining the motion of carbon nanoparticles (green tracks) (c) into the areas with singularities (red point) (a,b). The size of analyzed optical field is 1.8 × 1.8 µm2 (red square pointed in Fig. 5(a), (b), (c)).
Fig. 6.
Fig. 6. Demonstration of particles’ motion in time (see Visualization 1)
Fig. 7.
Fig. 7. Original optical field (a) (30 × 30 µm2), (c) (60 × 60 µm2) and corresponding reconstructed optical field (b), (d) through analyzing coordinate distribution of the nanoparticle motion tracks.
Fig. 8.
Fig. 8. Phase information: a) the restored phase map by applying of Hilbert transform, b) phase map with phase singularities (red points), identification of which is a result of the redistribution of carbon nanoparticles in the optical field; c) phase singularities obtained at the intersection of real (green lines) and imaginary (blue lines). Red points in white squares indicate the position of an intensity minimum without a singularity.

Equations (6)

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

U ( ξ , ζ ) = z i λ x = 1 X y = 1 Y F ( x , y ) R 2 ( x , y , z , ξ , ζ ) exp { i k [ R ( x , y , z , ξ , ζ ) + 2 h ( x , y ) ] } ,
Re [ α ] = α = r 3 ( ε r 1 ) ( ε r + 2 ) + ε i 2 ( ε r + 2 ) 2 + ε i 2 , Im [ α ] = α = r 3 3 ε i ( ε r + 2 ) 2 + ε i 2 .
F g r a d = α 2 n E 2 , F a b s = n S C a b s c , F s c a t t = n S C s c a t t c ,
Here  E 2 = | E | 2 2 , a S I = c 8 π | E | 2 .
m i d v i d t = F o p t i + F s t i ,
x i ( t ) = x i ( t 0 ) + v x i t + a x i t 2 2 , y i ( t ) = y i ( t 0 ) + v y i t + a y i t 2 2 .
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.