Next Article in Journal
Supermode Characteristics of Nested Multiple Hollow-Core Anti-Resonant Fibers
Next Article in Special Issue
Spatially Separated Heterodyne Grating Interferometer for In-Plane and Out-of-Plane Displacement Measurements
Previous Article in Journal
Two-Dimensional Elliptical Microresonator Arrays for Wide Flat Bandwidth and Boxlike Filter Response
Previous Article in Special Issue
Sensing Performance of Triple-Band Terahertz Metamaterial Absorber Based on Snowflake-Shaped Resonators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Surface Depth-Mapping of Material via the Transport-of-Intensity Equation

1
Laser and Optoelectronic Systems Department, Bauman Moscow State Technical University, 2nd Baumanskaya St. 5/1, 105005 Moscow, Russia
2
Lebedev Physical Institute, Russian Academy of Sciences, Leninskiy Prospekt 53, 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Photonics 2022, 9(11), 815; https://doi.org/10.3390/photonics9110815
Submission received: 20 September 2022 / Revised: 20 October 2022 / Accepted: 27 October 2022 / Published: 28 October 2022
(This article belongs to the Special Issue Micro-Nano Optics and High-End Measurement Instruments)

Abstract

:
We present a new approach for a surface characterization based on the TIE method combined with the SEM. Experimental verification is carried out on the example of characterization of a crater on the surface of monocrystalline silicon (111). The approach is universal and can be used for any opaque object. It improves the robustness and stability of the quantitative phase retrieval process and has two important features. Firstly, it allows one to quantitatively retrieve the phase in a region of arbitrarily chosen dimensions. Secondly, phase retrieval process does not require the choice of boundary conditions.

1. Introduction

One of the leading areas in modern laser physics research is the study of the interaction of ultrashort laser pulses with the surface of various materials, since these processes are fundamental in modern medicine [1], femtochemistry [2], creation of new nano- and microstructures for electronics [3,4], photonics [5], etc. The obvious advances in such studies are linked to the development of solid-state physics and theoretical methods for its modeling [6]. Combined, these methods can be used to determine surface structure [7,8] and chemical composition [9], study the electronic structure of solids, acquire images with atomic resolution and study processes occurring on the surface. It is generally accepted that there is no universal experimental method, and a set of mutually complementary methods should be used to solve these problems. For example, through combining methods of scanning probe microscopy, scanning electronic microscopy (SEM), energy dispersive X-ray spectroscopy and Raman spectroscopy, one can obtain complete information about the surface structure and its characteristics [10,11,12,13].
Among the abovementioned methods a special place is occupied by SEM, a powerful method of physical surface analysis allowing the establishment of composition, structure and some other properties of near-surface layers of the sample in the form of changes of intensity distribution in the observation plane. In most cases the interpretation of such images can be explained by phase contrast effects [14]. In addition to the resolution SEM, the phase contrast method was historically established as a qualitative tool to improve the image detail of the object under study [15]. At the same time, there are methods relating to objects’ phase retrieval which are based on the registration of several intensity distributions, acquired at displacement relative to the focus plane. They belong to the category of methods called inline holography [16,17,18], and one of them—based on the transport-of-intensity equation—will be discussed in this article.
The transport-of-intensity equation (TIE) is a well-established deterministic phase retrieval approach which can be used for quantitative measurements of bulk and surface characteristics of transparent media by analyzing phase distributions [19]. For instance, quantitative phase imaging of objects under study can be used to measure their refractive index [20,21], thickness, mass, birefringence and magnetic properties [22,23].
This paper proposes the use of TIE combined with SEM as an alternative to atomic force microscopy (AFM) [24,25] to address problem of surface characterization of solid opaque materials (on the example of monocrystalline silicon). The main result of this work is to provide a new approach to improve the robustness and stability of the phase retrieval process which can be used in case of uneven background and object illumination. It is shown that this approach solves the main difficulty when working with opaque objects, namely the need for a clear formulation of the boundary conditions.

2. Materials and Methods

The propagation of an electromagnetic wave with a wavelength λ 0 , a complex field E ( x , y , z ) along the z axis in an isotropic medium is described by a paraxial differential equation
i · E z + 1 2 k E + k E = 0
where E ( x , y , z ) = I ( x , y , z ) · e x p [ i φ ( x , y , z ) ] , I ( x , y , z ) —intensity of the wave, and φ ( x , y , z ) —phase of the wave; —two-dimensional Laplacian, k—wavenumber, E / z —longitudinal derivative.
After a series of transitions and simplifications in Equation (1) according to [26,27,28], one can obtain the following equation:
[ I φ ] = k · I z
where ∇ = {∂/∂x; ∂/∂y}—Nabla operator for the two-dimensional space; I / z —longitudinal derivative.
The resulting equation is nothing but the real part of the Helmholtz equation and is called the transport-of-intensity equation. It determines the one-to-one relationship between intensity I , its longitudinal derivative I / z and the phase φ of the electromagnetic wave.
There are many methods and approaches used to solve Equation (2). For instance, in [27] an approach based on the Helmholtz decomposition theorem, which allows the representation of the vector field as a scalar and vector potential, was proposed [29]. Since quantitative phase imaging via TIE belongs to the deterministic phase retrieval approach [30], it is necessary to restrict it only to cases where the field is free from phase singularities [31]. Under such constraints, the vector potential can be neglected and the TIE solution can be obtained using only the phase scalar potential without any generality restriction. Then the original Equation (2) is simplified to the standard Poisson equation written with respect to the phase.
φ ( x , y , z ) = k 2 [ · ( 1 I ( x , y ) 2 I ( x , y ) z ) ]
One of the most popular methods for solving Equation (3) is a method based on expressing differential operators through the fast Fourier transform (FFT) [32]. In [33] it is shown that this approach bypasses the serious mathematical limitations in solving TIE which arise when using finite elements methods. Using the properties of the Fourier transform, the differential operators in Equation (3) can be written as [34]:
{ f ( x , y ) } = F 1 { [ j k x + j k y ] · F { f ( x , y ) } }
2 { f ( x , y ) } = F 1 { [ 1 ( j k x ) 2 + ( j k y ) 2 ] · F { f ( x , y ) } }
where F и F−1—operators of direct and inverse Fourier transforms, respectively; kx,y = 2πνx,y—frequency coefficients, νx,y—frequency grids.
In [35] it was proved that this transition preserves the correctness and uniqueness of the solution to an additive constant. It should be noted that this approach requires strict adherence to boundary conditions. Violation of boundary conditions leads to distortion of boundaries and overlapping of spectra in the whole selected area. Due to the periodic nature of the Fourier transform, periodic boundary conditions (PBC) must be satisfied, where the function value at the boundaries must cyclically repeat [36], before applying relations (4) and (5).
In general, the choice of the boundary conditions when solving TIE is a very difficult task, because the effect of the boundary conditions is independent of the chosen solution method. Due to the periodic nature of the FFT, the boundary conditions are implicit. However, this implication can be eliminated with an extended FFT method that uses a symmetrization rule [37]. This allows the use of special cases of Neumann boundary conditions (NBC, Equation (6)), Dirichlet boundary conditions (DBC, Equation (7)), PBC or any other mixed boundary conditions [38]. The correctness of using NBC, DBC and PBC when solving TIE using the extended FFT method has been proven by numerical analysis, which investigates the effect of boundary conditions on the accuracy of phase retrieval [36]. In this case, it is assumed that the intensity distribution turns to zero outside the selected area [35].
I φ n | Ω = g
where g is a smooth function on the boundary Ω ; φ n is the outward normal derivative.
φ | Ω = g
where g is a smooth function on the boundary Ω .
In all cases, the derivation of the boundary conditions follows from the law of energy conservation, which in a more general case is replaced by the law of intensity conservation which is assumed to be true for TIE [26,35]. According to [35,39], the law of intensity conservation will be satisfied only if the dissipation of energy flux through the selected circuit is negligible. In practice, however, this condition is met by increasing the ratio between the size of the selected area and the size of the object itself. This minimizes the change in overall illumination in the image plane in the selected area.

3. Experimental Demonstration

An experimental demonstration of the method was carried out by the characterization of an ablation crater obtained on the surface of single-crystal silicon with a crystallographic orientation (111). This crater was obtained by tight focusing of laser radiation with wavelength λ = 515 nm and pulse duration τ = 300 fs through a micro-objective with numerical aperture NA = 0.65 into a spot with radius R1/e ≈ 0.6 μm. Ablation was performed in a single-pulse mode, with the energy per pulse being E = 450 nJ.
Pre-characterization of the ablation crater was performed using Certus Standard V AFM («Nano Scan Technology» Ltd., Dolgoprudnyy, Russia), which has lateral spatial resolution <0.1 nm and vertical resolution <0.01 nm. AFM scanning was performed in the semi-contact mode using a NSG10 cantilever. Figure 1 shows a normalized 2D image of the crater under study with a depth profile. Image size is 6.05 × 6.05 μm and the spatial resolution is 108 nm.
To demonstrate the applicability of TIE for three-dimensional measurement of the surface depth profile with atomic resolution, SEM images (Tescan Vega, Brno, Czech Republic) were used. A series of images was made with 3 keV accelerating voltage and 30 pA beam current. According to [40], this mode provides high spatial resolution and minimum depth of focus, and is optimal for this particular sample.
The series of images was formed by tuning the configuration of the electron column, namely by changing the working distance (WD) of the system (the distance between the lower pole tip of the lens and the focusing plane of the electron beam). The nominal WD value was ~11.9 mm. The focusing plane was displaced relative to the initial WD position by 100 μm with 5 μm steps, with a displacement accuracy of ~1 nm. The choice of the displacement range and steps was based on two limiting factors: to prevent high spatial frequency attenuation due to spatial coherence; and to ensure that the change in the intensity distribution was sufficient to calculate the longitudinal intensity derivative. Some of the obtained SEM-images are shown on Figure 2a–c.

4. Results and Discussion

The AFM data has a non-normalized dynamic range and so before analysis it was normalized to a minimum (that is, the point with the greatest depth) in order to simplify the analysis. Figure 3 shows cross-sections of the depth profile along the 0X and 0Y axes through the zero-depth point. The horizontal dashed line indicates the silicon surface level, while vertical lines show the abstract intersection of the crater’s inner surface with silicon surface level. The resulting crater depth is 432 ± 1 nm.
Quantitative phase imaging requires a priori information about the wavelength of the radiation. For SEM the theoretical value of the wavelength λ 0 will be equal to the de Broglie wavelength of the electron passing the accelerating potential U
λ 0 h 2 m e e U
where h—Planks constant, m e —mass of the electron, e—elementary charge, U —acceleration voltage.
The boundary conditions and the size of the area under study are the parameters when solving a TIE. Since the boundary conditions of the problem to be solved are indefinite, we consider them in the framework of the Dirichlet-Neumann boundary problem. In turn, the choice of the size of the area under study is derived from the law of energy conservation, which requires a negligibly small value of energy flux dissipation through the selected contour. Since such a calculation is not possible in this case, we will empirically determine the optimal ratio between the area under study and the object size.
Figure 4a–c show a normalized 2D image with a depth profile for the case when the size of the area under study corresponds to the size of the object for three different cases: without imposing any boundary conditions; boundary conditions according to Neumann; and boundary conditions according to Dirichlet, respectively. Figure 4d–f show the results for the case where the size of the area under study was increased by a factor of 2 for each of the axis.
In particular, in Figure 4a–c the object is positioned so that it is almost in contact with the image boundaries, resulting in serious boundary artifacts that affect the accuracy of phase retrieval and the quantification of crater depth. Figure 4d–f show that the influence of the boundary conditions on the measuring of the crater depth and the contribution of heterogeneity (non-uniformity of the background value and object illumination) decreases as the area under study increases, at least visually. However, it will be numerically shown below that an increase in the area under study increases the error for some cases of boundary conditions while reconstructing the depth maps.
The accuracy of phase retrieval using TIE depends not only on the appropriately chosen boundary conditions, as demonstrated earlier, but also on the quasi-uniformity of the intensity distribution. It is worth noting that the background intensity values on the SEM-images have a normal Gaussian intensity distribution with a mean value of 45 and dispersion of 22.5. In practice, this condition is unfeasible, so it was ensured by methods of digital image processing. The use of such methods does not contradict the conditions of intensity singularity (presence of small and/or zero intensity values). Figure 5 shows 2D crater images obtained after SEM-images were processed digitally. The boundary conditions and TIE calculation are similar to those in Figure 4.
It can be seen from Figure 5 that images with a uniform background resulted in a more stable TIE solution. As a result, the error of the depth profile reconstruction decreased. This finding is confirmed by the calculation of the root-mean-square error (RMSE) between depth maps obtained by AFM h A F M ( x , y ) and SEM h S E M ( x , y ) with TIE, according to
R M S E =   [ h S E M ( x , y ) h A F M ( x , y ) ] 2 d x d y   h A F M ( x , y ) 2 d x d y
The RMSE was calculated for all cases of boundary conditions and several values of the size of the area under study, which ranged from 4.875 to 15.6 μm. Its calculation was performed only within the area where the ablation crater was located, regardless of the size of the area under study used for the TIE. The resulting dependences of the RMSE on the size of the studied region used for TIE for different boundary conditions are shown in Figure 6a–c.
The results of the calculations show that before the digital processing of SEM images, the best result was achieved using Dirichlet boundary conditions. Figure 6c shows that the error value decreases as the size of the study area increases, while for no boundary (Figure 6a) and Neumann boundary (Figure 6b) conditions the error monotonically increases. If we speculate on the area size which could net the minimum RSME an extrapolation of the data from Figure 6c is needed. Assuming the dependence is linear the optimal size of the region under study for the Dirichlet boundary conditions would be ~50 μm. If we are considering the results obtained after digital processing, again Dirichlet boundary conditions give the best results. In this case, the error of the depth map measurement remains practically unchanged in the selected range of the areas under study (except for some points).
As mentioned earlier, AFM has lateral resolution <0.1 nm (X and Y directions), and longitudinal resolution <0.01 nm (Z direction). By comparison, the SEM has high lateral resolution (~2 nm [40]) and low longitudinal resolution [41]. Therefore, combining SEM-imaging with the TIE phase retrieval method proposed in this paper, we can obtain the depth map of objects with much higher longitudinal resolution. As seen in Figure 6a–c, the cumulative RSME of depth maps reconstruction was no more than 6% for all cases considered, and shows that the proposed method allows the carrying out of surface characterization compatible with that of AFM.

5. Conclusions

We presented a new approach for quantitative surface characterization based on applying TIE to SEM-images, and demonstrated its feasibility using the example of a crater on the surface of monocrystalline silicon. The proposed approach is universal and can be used for any opaque object, can increase the robustness and stability of the quantitative phase retrieval process and has two important features. Firstly, it allows one to retrieve the quantitative phase in an area with arbitrarily chosen dimensions. Secondly, no choice of boundary conditions is required for phase retrieval. We proved that if the intensity is positive in a given region of the object and converts to a constant value there is no violation of the intensity singularity condition. The solution of the TIE in such a case is singular to the additive constant in the absence of any boundary conditions. We have demonstrated that by applying this TIE solution to SEM-images it is possible to produce high-resolution surface characterization compatible with that of AFM. The cumulative error in depth map reconstruction was no more than 6% for all cases considered, but the best results were obtained by solving this problem using Dirichlet boundary conditions.

Author Contributions

Conceptualization, M.K. and S.K.; methodology, I.P. and Y.G.; software, I.P. and N.S.; validation, M.K. and Y.G.; formal analysis, G.K. and S.K.; investigation, N.S. and G.K.; resources, S.K.; data curation, M.K. and S.K.; writing—original draft preparation, N.S. and G.K.; writing—review and editing, N.S., G.K. and M.K.; visualization, N.S.; supervision, G.K.; project administration, M.K. and S.K.; funding acquisition, S.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research was financially supported by the Russian Science Foundation (Project No. 21-79-30063).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are available from the corresponding author upon reasonable request.

Acknowledgments

We thank Iliya Gritsenko for his technical support with the software.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vogel, A.; Venugopalan, V. Mechanisms of Pulsed Laser Ablation of Biological Tissues. Chem. Rev. 2003, 103, 577–644. [Google Scholar] [CrossRef] [Green Version]
  2. Assion, A.; Baumert, T.; Bergt, M.; Brixner, T.; Kiefer, B.; Seyfried, V.; Strehleand, M.; Gerber, G. Control of Chemical Reactions by Feedback-Optimized Phase-Shaped Femtosecond Laser Pulses. Science 1998, 282, 919–922. [Google Scholar] [CrossRef]
  3. Yangdong, W.; Haibo, Y.; Yuzhao, Z.; Ye, Q.; Peiwen, L.; Xiaoduo, W.; Boliang, J.; Lianqing, L.; Wen Jung, L. Recent Advances in Femtosecond Laser Fabrication: From Structures to Applications. IEEE Open J. Nanotechnol. 2021, 2, 161–177. [Google Scholar]
  4. Gattass, R.; Mazur, E. Femtosecond laser micromachining in transparent materials. Nat. Photonics 2008, 2, 219–225. [Google Scholar] [CrossRef]
  5. Ponkratova, E.; Ageev, E.; Komissarenko, F.; Koromyslov, S.; Kudryashov, D.; Mukhin, I.; Veiko, V.; Kuchmizhak, A.; Zuev, D. Femtosecond Laser Fabrication of Hybrid Metal-Dielectric Structures with Nonlinear Photoluminescence. Photonics 2021, 8, 121. [Google Scholar] [CrossRef]
  6. Patterson, J.; Bailey, B. Solid-State Physics Introduction to the Theory, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  7. Wedler, H. Information on Surface Structure by Low Energy Electron Diffraction. Vakuum 1995, 7, 107–114. [Google Scholar] [CrossRef]
  8. Galiy, P.V.; Losovy, J.; Ya, B.; Nenchuk, T.; Yarovets, I.R. Low-Energy-Electron-Diffraction Structural Studies of (100) Cleavage Surfaces of In4Se3 Layered Crystals. Ukr. J. Phys. 2014, 59, 612–621. [Google Scholar] [CrossRef]
  9. Chen, L.X.; Lee, P.L.; Gosztola, D.; Svec, W.A.; Montano, P.A.; Wasielewski, M.R. Time-Resolved X-ray Absorption Determination of Structural Changes following Photoinduced Electron Transfer within Bis-porphyrin Heme Protein Models. J. Phys. Chem. B 1999, 103, 3270–3274. [Google Scholar] [CrossRef]
  10. Gawelda, W.; Saes, M.; Kaiser, M.; Tarnovsky, A.; Johnson, S.L.; Grolimund, D.; Abela, R.; Chergui, M.; Bressler, C. Structural dynamics and electronic structure changes probed with lasers and X-rays. In Femtochemistry and Femtobiology, Proceedings of the Ultrafast Events in Molecular Science VIth International Conference on Femtochemistry Maison de la Chimie, Paris, France, 6–10 July 2003; Martin, M.M., Hynes, J.T., Eds.; Elsevier Science: Amsterdam, The Netherlands, 2004. [Google Scholar]
  11. Petrone, L.; Ragg, N.L.C.; Girvan, L.; McQuillan, A.J. Scanning Electron Microscopy and Energy Dispersive X-Ray Microanalysis of Perna canaliculus Mussel Larvae Adhesive Secretion. J. Adhes. 2009, 85, 78–96. [Google Scholar] [CrossRef]
  12. Sharma, S.; Green, M.; Hyer, R.; Dark, C.; Black, T.; Chourasia, A.; Mishra, K. Growth of diamond films and characterization by Raman, scanning electron microscopy, and X-ray photoelectron spectroscopy. J. Mater. Res. 1990, 5, 2424–2432. [Google Scholar] [CrossRef]
  13. Sharma, S.C.; Dark, C.A.; Hyer, R.C.; Green, M.; Black, T.D. Deposition of diamond films at low pressures and their characterization by positron annihilation, Raman, scanning electron microscopy, and X-ray photoelectron spectroscopy. Appl. Phys. Lett. 1990, 56, 1781. [Google Scholar] [CrossRef]
  14. Hecht, B.; Sick, B.; Wild, U.P. Scanning near-field optical microscopy with aperture probes: Fundamentals and applications. J. Chem. Phys. 2000, 112, 7761. [Google Scholar] [CrossRef]
  15. Maurer, C.; Jesacher, A.; Bernet, S.; Ritsch-Marte, M. Phase contrast microscopy with full numerical aperture illumination. Opt. Express 2008, 16, 19821–19829. [Google Scholar] [CrossRef] [Green Version]
  16. Latychevskaia, T.; Formanek, P.; Koch, C.T.; Lubk, A. Off-axis and inline electron holography: Experimental comparison. Ultramicroscopy 2010, 110, 472–482. [Google Scholar] [CrossRef] [Green Version]
  17. Zhou, W.-J.; Guan, X.; Liu, F.; Yu, Y.; Zhang, H.; Poon, T.-C.; Banerjee, P.P. Phase retrieval based on transport of intensity and digital holography. Appl. Opt. 2018, 57, A229–A234. [Google Scholar] [CrossRef]
  18. Zhou, W.; Liu, S.; Wang, C.; Zhang, H.; Yu, Y.; Poon, T.-C. Elimination of Quadratic Phase Aberration in Digital Holographic Microscopy by Using Transport of Intensity. Front. Photonics 2022, 3, 848453. [Google Scholar] [CrossRef]
  19. Zuo, C.; Li, J.; Sun, J.; Fan, Y.; Zhang, J.; Lu, L.; Zhang, R.; Wang, B.; Huang, L.; Chen, Q. Transport of intensity equation: A tutorial. Opt. Lasers Eng. 2020, 135, 106187. [Google Scholar] [CrossRef]
  20. Gritsenko, I.V.; Kovalev, M.S.; Stsepuro, N.G.; Gulina, Y.S.; Krasin, G.K.; Gonchukov, S.A.; Kudryashov, S.I. The optical refractometry using transport-of-intensity equation. Laser Phys. Lett. 2022, 19, 076201. [Google Scholar] [CrossRef]
  21. Kovalev, M.; Gritsenko, I.; Stsepuro, N.; Nosov, P.; Krasin, G.; Kudryashov, S. Reconstructing the Spatial Parameters of a Laser Beam Using the Transport-of-Intensity Equation. Sensors 2022, 22, 1765. [Google Scholar] [CrossRef]
  22. McCray, A.R.C.; Cote, T.; Li, Y.; Petford-Long, A.K.; Phatak, C. Understanding Complex Magnetic Spin Textures with Simulation-Assisted Lorentz Transmission Electron Microscopy. Phys. Rev. Applied 2012, 15, 044025. [Google Scholar] [CrossRef]
  23. Zhidong, H.; Dongsheng, S.; Weiwei, W.; Ning, W.; Binghui, G.; Shouguo, W.; Mingliang, T.; Haifeng, D. Visualizing Emergent Magnetic Flux of Antiskyrmions in Mn1.4PtSn Magnet. Adv. Funct. Mater. 2022, 32, 2112661. [Google Scholar]
  24. Cao, F.; Donnarumma, F.; Murray, K.K. Wavelength dependent atomic force microscope tip-enhanced laser ablation. Appl. Surf. Sci. 2018, 447, 437–441. [Google Scholar] [CrossRef]
  25. Smirnov, N.A.; Kudryashov, S.I.; Rudenko, A.A.; Zayarny, D.A.; Ionin, A.A. Pulsewidth and ambient medium effects during ultrashort-pulse laser ablation of silicon in air and water. Appl. Surf. Sci. 2021, 5621, 150243. [Google Scholar] [CrossRef]
  26. Teague, M.R. Irradiance moments: Their propagation and use for unique retrieval of phase. J. Opt. Soc. Am. 1982, 72, 1199–1209. [Google Scholar] [CrossRef]
  27. Teague, M.R. Deterministic phase retrieval: A green’s function solution. J. Opt. Soc. Am. 1983, 73, 1434–1441. [Google Scholar] [CrossRef]
  28. Teague, M.R. Image formation in terms of the transport equation. J. Opt. Soc. Am. A 1985, 2, 2019–2026. [Google Scholar] [CrossRef]
  29. Philip, M.M.; Feshbach, H. Methods of Theoretical Physics, Part I.; McGraw-Hill: New York, USA, 1953. [Google Scholar]
  30. Semichaevsky, A.; Testorf, M. Phase-space interpretation of deterministic phase retrieval. JOSA A 2004, 21, 2173–2179. [Google Scholar] [CrossRef]
  31. Gureyev, T.E.; Roberts, A.; Nugent, K.A. Partially coherent fields, the transport of intensity equation and phase uniqueness. JOSA A 1995, 12, 1942–1946. [Google Scholar] [CrossRef]
  32. Paganin, D.; Nugent, K.A. Noninterferometric Phase Imaging with Partially Coherent Light. Phys. Rev. Lett. 1998, 80, 2586–2589. [Google Scholar] [CrossRef]
  33. De Graef, M. Lorentz microscopy: Theoretical basis and image simulations. In Magnetic Imaging and Its Applications to Materials, 1st ed.; De Graef, M., Zhu, Y., Eds.; Academic Press: Cambridge, MA, USA, 2001; Volume 38, pp. 27–67. [Google Scholar]
  34. Paganin, D. Coherent X-ray optics. Oxford University Press: Oxford, UK, 2006. [Google Scholar]
  35. Gureyev, T.E.; Nugent, K.A. Phase retrieval with the transport-of-intensity equation. II. Orthogonal series solution for nonuniform illumination. J. Opt. Soc. Am. A 1996, 13, 1670–1682. [Google Scholar] [CrossRef]
  36. Martinez-Carranza, J.; Falaggis, K.; Kozacki, T.; Kujawinska, M. Effect of imposed boundary conditions on the accuracy of transport of intensity equation based. Model. Asp. Opt. Metrol. IV 2013, 8789, 195–208. [Google Scholar]
  37. Volkov, V.; Zhu, Y.; De Graef, M. A new symmetrized solution for the phase retrieval using the Transport of Intensity Equation. Micron 2002, 33, 411–416. [Google Scholar] [CrossRef]
  38. Zuo, C.; Chen, Q.; Li, H.; Qu, W.; Asundi, A. Boundary-artifact-free phase retrieval with the transport of intensity equation II: Applications to microlens characterization. Opt. Express 2014, 22, 18310–18324. [Google Scholar] [CrossRef] [PubMed]
  39. Piskunov, N. Differential and Integral Calculus; Noordhoff: Groningen, The Netherlands, 1965. [Google Scholar]
  40. Available online: https://www.tescan.com (accessed on 20 September 2022).
  41. Palchoudhury, S.; Baalousha, M.; Lead, J.R. Methods for Measuring Concentration (Mass, Surface Area and Number) of Nanomaterials. In Frontiers of Nanoscience; Baalousha, M., Lead, J.R., Eds.; Elsevier: Amsterdam, The Netherlands, 2015; Volume 8, pp. 153–181. [Google Scholar]
Figure 1. 2D AFM crater image with a depth profile.
Figure 1. 2D AFM crater image with a depth profile.
Photonics 09 00815 g001
Figure 2. SEM-images of the crater under study: (a) in the focusing plane; (b) displaced by 50 μm; (c) displaced by 100 μm.
Figure 2. SEM-images of the crater under study: (a) in the focusing plane; (b) displaced by 50 μm; (c) displaced by 100 μm.
Photonics 09 00815 g002
Figure 3. Cross-sections of the AFM depth profile along the 0X and 0Y axes.
Figure 3. Cross-sections of the AFM depth profile along the 0X and 0Y axes.
Photonics 09 00815 g003
Figure 4. 2D crater images retrieved using TIE for different conditions: (a,d) no boundary conditions; (b,e) Neumann boundary conditions; (c,f) Dirichlet boundary conditions.
Figure 4. 2D crater images retrieved using TIE for different conditions: (a,d) no boundary conditions; (b,e) Neumann boundary conditions; (c,f) Dirichlet boundary conditions.
Photonics 09 00815 g004aPhotonics 09 00815 g004b
Figure 5. 2D crater images retrieved using TIE from digitally processed SEM-images: (a,d) no boundary conditions; (b,e) Neumann boundary conditions; (c,f) Dirichlet boundary conditions.
Figure 5. 2D crater images retrieved using TIE from digitally processed SEM-images: (a,d) no boundary conditions; (b,e) Neumann boundary conditions; (c,f) Dirichlet boundary conditions.
Photonics 09 00815 g005aPhotonics 09 00815 g005b
Figure 6. RMSE between depth maps obtained by AFM and SEM for: (a) no boundary conditions; (b) Neumann boundary conditions; (c) Dirichlet boundary conditions.
Figure 6. RMSE between depth maps obtained by AFM and SEM for: (a) no boundary conditions; (b) Neumann boundary conditions; (c) Dirichlet boundary conditions.
Photonics 09 00815 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stsepuro, N.; Kovalev, M.; Krasin, G.; Podlesnykh, I.; Gulina, Y.; Kudryashov, S. Surface Depth-Mapping of Material via the Transport-of-Intensity Equation. Photonics 2022, 9, 815. https://doi.org/10.3390/photonics9110815

AMA Style

Stsepuro N, Kovalev M, Krasin G, Podlesnykh I, Gulina Y, Kudryashov S. Surface Depth-Mapping of Material via the Transport-of-Intensity Equation. Photonics. 2022; 9(11):815. https://doi.org/10.3390/photonics9110815

Chicago/Turabian Style

Stsepuro, Nikita, Michael Kovalev, George Krasin, Ivan Podlesnykh, Yulia Gulina, and Sergey Kudryashov. 2022. "Surface Depth-Mapping of Material via the Transport-of-Intensity Equation" Photonics 9, no. 11: 815. https://doi.org/10.3390/photonics9110815

APA Style

Stsepuro, N., Kovalev, M., Krasin, G., Podlesnykh, I., Gulina, Y., & Kudryashov, S. (2022). Surface Depth-Mapping of Material via the Transport-of-Intensity Equation. Photonics, 9(11), 815. https://doi.org/10.3390/photonics9110815

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop