Full-range parallel Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We demonstrate a full-range parallel Fourier-domain optical coherence tomography (FD-OCT) in which a tomogram free of mirror images as well as DC and autocorrelation terms is obtained in parallel. The phase and amplitude of two-dimensional spectral interferograms are accurately detected by using sinusoidal phase-modulating interferometry and a two-dimensional CCD camera, which allows for the reconstruction of two-dimensional complex spectral interferograms. By line-by-line inverse Fourier transformation of the two-dimensional complex spectral interferogram, a full-range parallel FD-OCT is realized. Tomographic images of two separated glass coverslips obtained with our method are presented as a proof-of-principle experiment.

Advanced Monte Carlo Simulation and Machine Learning for Frequency Domain Optical Coherence Tomography

Optical Coherence Tomography(OCT) is a popular, rapidly growing imaging technique with an increasing number of bio-medical applications due to its noninvasive nature. However, there are three major challenges in understanding and improving an OCT system: (1) Obtaining an OCT image is not easy. It either takes a real medical experiment or requires days of computer simulation. Without much data, it is difficult to study the physical processes underlying OCT imaging of different objects simply because there aren't many imaged objects. (2) Interpretation of an OCT image is also hard. This challenge is more profound than it appears. For instance, it would require a trained expert to tell from an OCT image of human skin whether there is a lesion or not. This is expensive in its own right, but even the expert cannot be sure about the exact size of the lesion or the width of the various skin layers. The take-away message is that analyzing an OCT image even from a high level would usually require a trained expert, and pixel-level interpretation is simply unrealistic. The reason is simple: we have OCT images but not their underlying ground-truth structure, so there is nothing to learn from. (3) The imaging depth of OCT is very limited (millimeter or sub-millimeter on human tissues). While OCT utilizes infrared light for illumination to stay noninvasive, the downside of this is that photons at such long wavelengths can only penetrate a limited depth into the tissue before getting back-scattered. To image a particular region of a tissue, photons first need to reach that region. As a result, OCT signals from deeper regions of the tissue are both weak (since few photons reached there) and distorted (due to multiple scatterings of the contributing photons). This fact alone makes OCT images very hard to interpret.

This thesis addresses the above challenges by successfully developing an advanced Monte Carlo simulation platform which is 10000 times faster than the state-of-the-art simulator in the literature, bringing down the simulation time from 360 hours to a single minute. This powerful simulation tool not only enables us to efficiently generate as many OCT images of objects with arbitrary structure and shape as we want on a common desktop computer, but it also provides us the underlying ground-truth of the simulated images at the same time because we dictate them at the beginning of the simulation. This is one of the key contributions of this thesis. What allows us to build such a powerful simulation tool includes a thorough understanding of the signal formation process, clever implementation of the importance sampling/photon splitting procedure, efficient use of a voxel-based mesh system in determining photon-mesh interception, and a parallel computation of different A-scans that consist a full OCT image, among other programming and mathematical tricks, which will be explained in detail later in the thesis.

Next we aim at the inverse problem: given an OCT image, predict/reconstruct its ground-truth structure on a pixel level. By solving this problem we would be able to interpret an OCT image completely and precisely without the help from a trained expert. It turns out that we can do much better. For simple structures we are able to reconstruct the ground-truth of an OCT image more than 98% correctly, and for more complicated structures (e.g., a multi-layered brain structure) we are looking at 93%. We achieved this through extensive uses of Machine Learning. The success of the Monte Carlo simulation already puts us in a great position by providing us with a great deal of data (effectively unlimited), in the form of (image, truth) pairs. Through a transformation of the high-dimensional response variable, we convert the learning task into a multi-output multi-class classification problem and a multi-output regression problem. We then build a hierarchy architecture of machine learning models (committee of experts) and train different parts of the architecture with specifically designed data sets. In prediction, an unseen OCT image first goes through a classification model to determine its structure (e.g., the number and the types of layers present in the image); then the image is handed to a regression model that is trained specifically for that particular structure to predict the length of the different layers and by doing so reconstruct the ground-truth of the image. We also demonstrate that ideas from Deep Learning can be useful to further improve the performance.

It is worth pointing out that solving the inverse problem automatically improves the imaging depth, since previously the lower half of an OCT image (i.e., greater depth) can be hardly seen but now becomes fully resolved. Interestingly, although OCT signals consisting the lower half of the image are weak, messy, and uninterpretable to human eyes, they still carry enough information which when fed into a well-trained machine learning model spits out precisely the true structure of the object being imaged. This is just another case where Artificial Intelligence (AI) outperforms human. To the best knowledge of the author, this thesis is not only a success but also the first attempt to reconstruct an OCT image at a pixel level. To even give a try on this kind of task, it would require fully annotated OCT images and a lot of them (hundreds or even thousands). This is clearly impossible without a powerful simulation tool like the one developed in this thesis.

Primary and secondary threshold intensities of ultraviolet-laser-induced domain nucleation in nearly stoichiometric LiTaO3

The primary and secondary threshold intensities of ultraviolet-laser-induced preferential domain nucleation in nearly stoichiometric LiTaO3 is observed. The primary threshold is the minimum intensity to achieve the instantaneous preferential domain nucleation within the focus by the combined action of irradiation and electric fields. The secondary threshold is the minimum intensity to achieve the memory effect without any irradiation within the original focus. The space charge field created by the photoionization carriers is thought to be responsible for the instantaneous effect. The explanation based on the formation and transformation of extrinsic defect is presented for the memory effect. (c) 2008 American Institute of Physics.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Theory of high-power mode-locked lasers with a slow absorber

We develop an analytical theory of high-power passively mode-locked lasers with a slow absorber; the theory is valid at pulse energies well exceeding the saturation energy. We analyze the Haus modelocking master equation in the pulse-energy-domain representation, approximating the intensity profile function by a series in the vicinity of its peak value. We consider the high-power operation regime of subpicosecond blue-violet GaN mode-locked diode lasers, using the approach developed. © 2010 Springer Science+Business Media, Inc.

Mode Simulation for Midinfrared Microsquare Resonators With Sloped Sidewalls and Confined Metals

EQUILATERAL-TRIANGLE; MU-M; LASERS; MICROLASERS; MICRODISK Abstract: Mode characteristics for midinfrared microsquare resonators with sloped sidewalls and confined metal layers are investigated by finite-difference time-domain (FDTD) techniques. For a microsquare with a side length of 10 mu m, the mode quality (Q)-factors of 8329, 4772, and 2053 are obtained for TM5,7 mode at wavelength 7.1 mu m by three-dimensional FDTD simulations, as the tilting angles of the side walls are 90 degrees, 88 degrees, and 86 degrees, respectively. Furthermore, microsquare resonators laterally surrounded by SiO2 and metal layers are investigated by the two-dimensional FDTD technique for the metal layers of Au, Ti-Au, Ag-Au, and Ti-Ag-Au, respectively.

Mode Characteristics for Square Resonators With a Metal Confinement Layer

Microsquare resonators laterally confined by SiO2/Au/air multilayer structure are investigated by light ray method with reflection phase-shift of the multiple layers and two-dimensional (2-D) finite-difference time-domain (FDTD) technique. The reflectivity and phase shift of the mode light ray on the sides of the square resonator with the semiconductor/SiO2/Au/air multilayer structure are calculated for TE and TM modes by transfer matrix method. Based on the reflection phase shift and the reflectivity, the mode wavelength and factor are calculated by the resonant condition and the mirror loss, which are in agreement well with that obtained by the FDTD simulation. We find that the mode factor increases greatly with the increase of the SiO2 layer thickness, especially as d < 0.3 mu m. For the square resonator with side length 2 mu m and refractive index 3.2, anticrossing mode couplings are found for confined TE modes at wavelength about 1.6 mu m at d = 0.11 mu m, and confined TM modes at d = 0.71 mu m, respectively.

Investigation of mode radiation loss for microdisk resonators with pedestals by FDTD technique

Mode radiation loss for microdisk resonators with pedestals is investigated by three-dimensional (3D) finite-difference time-domain (FDTD) technique. For the microdisk with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, on a pedestal with a refractive index of 3.17, the mode quality (Q) factor of the whispering-gallery mode (WGM) quasi-TE7,1 first increases with the increase of the radius of the pedestal, and then quickly decreases as the radius is larger than 0.75 mu m. The mode radiation loss is mainly the vertical radiation loss induced by the mode coupling between the WGM and vertical radiation mode in the pedestal, instead of the scattering loss around the perimeter of the round pedestal. The WG M can keep the high Q factor when the mode coupling is forbidden.

Mode characteristics of metallically coated square microcavity connected with an output waveguide

Mode characteristics of a square microcavity with an output waveguide on the middle of one side, laterally confined by an insulating layer SiO2 and a p-electrode metal Au, are investigated by two-dimensional finite-difference time-domain technique. The mode quality (Q) factors versus the width of the output waveguide are calculated for Fabry-Peacuterot type and whispering-gallery type modes in the square cavity. Mode coupling between the confined modes in the square cavity and the guided modes in the output waveguide determines the mode Q factors, which is greatly influenced by the symmetry behaviors of the modes. Fabry-Peacuterot type modes can also have high Q factors due to the high reflectivity of the Au layer for the vertical incident mode light rays. For the square cavity with side length 4 mu m and refractive index 3.2, the mode Q factors of the Fabry-Peacuterot type modes can reach 10(4) at the mode wavelength of 1.5 mu m as the output waveguide width is 0.4 mu m.

Investigation on multiple-port microcylinder lasers based on coupled modes

Microcylinder resonators with multiple ports connected to waveguides are investigated by 2D finite-difference time-domain (FDTD) simulation for realizing microlasers with multiple outputs. For a 10 mu m radius microcylinder with a refractive index of 3.2 and three 2 mu m wide waveguides, confined mode at the wavelength of 1542.3 nm can have a mode Q factor of 6.7 x 10(4) and an output coupling efficiency of 0.76. AlGaInAs/InP microcylinder lasers with a radius of 10 mu m and a 2 mu m wide output waveguide are fabricated by planar processing techniques. Continuous-wave electrically injected operation is realized with a threshold current of 4 mA at room temperature, and the jumps of output power are observed accompanying a lasing mode transformation.

Mode coupling and vertical radiation loss for whispering-gallery modes in 3-D microcavities

The characteristics of whispering-gallery modes (WGMs) in 3-D cylindrical, square, and triangular microcavities with vertical optical confinement of semiconductors are numerically investigated by the finite-difference time-domain (FDTD) technique. For a microcylinder with a vertical refractive index 3.17/3.4/3.17 and a center layer thickness 0.2 mu m, Q-factors of transverse electric (TE) WGMs around wavelength 1550 nm are smaller than 10(3), as the radius R < 4 mu m and reach the orders of 10(4) and 10(6) as R = 5 and 6 mu m, respectively. However, the Q-factor of transverse magnetic (TM) WGMs at wavelength 1.659 mu m reaches 7.5 x 10(5) as R = 1 mu m. The mode coupling between the WGMs and vertical radiation modes in the cladding layer results in vertical radiation loss for the WGMs. In the microcylinder, the mode wavelength of TM WGM is larger than the cutoff wavelength of the vertical radiation mode with the same mode numbers, so TM WGMs cannot couple with the vertical radiation mode and have high Q-factor. In contrast, TE WGMs can couple with the corresponding vertical radiation mode in the 3-D microcylinder as R < 5 mu m. However, the mode wavelength of the TE WGM approaches (is larger than) the cutoff wavelength of the corresponding radiation modes at R = 5 mu m (6 mu m), so TE WGMs have high Q-factors in such microcylinders too. The results show that a critical lateral size is required for obtaining high, Q-factor TE WGMs in the 3-D microcylinder. For 3-D square and triangular microcavities, we also find that the Q-factor of TM WGM is larger than that of TE WGM.

Analysis of mode characteristics for equilateral-polygonal resonators with a center hole

Mode characteristics for equilateral triangles, squares, and hexagonal resonators with a center hole are numerically simulated by the finite-different time domain (FDTD) technique. The center hole does not break the symmetry behavior of the original resonators and can result in modification of the mode field patterns and mode Q factors. In an equilateral triangle resonator the center hole can suppress the symmetry state of degenerate states with the merit of single mode operation. In a square resonator, the Q factor can be enhanced for some modes with a suitable size of the hole. For a hexagonal resonator with a side length of 1 mu m and a refractive index of 3.2, the mode Q factors first gradually decrease with the increase of the hole diameter for modes at a wavelength of about 1500 nm, then the modes transform to that of a microdisk with a jump of the mode wavelength as the hole diameter approaches 0.7 mu m. Finally, the mode Q factors greatly enhance as the hole diameter reaches about 1 mu m. The results indicate that the center hole can greatly modify mode characteristics, especially that of the mode Q factor. (C) 2009 Optical Society of America

Mode Analysis for Equilateral-Triangle-Resonator Microlasers with Metal Confinement Layers

Mode characteristics are analyzed for electrically injected equilateral-triangle-resonator (ETR) semiconductor microlasers, which are laterally confined by insulating barrier SiO2 and electrode metals Ti-Au. For the ETR without metal layers, the totally confined mode field patterns are derived based on the reflection phase shifts, and the Q-factors are calculated from the far-field emission of the analytical near field distribution, which are agreement very well with the numerical results of the finite-difference time-domain (FDTD) simulation. The polarization dependence reflections for light rays incident on semiconductor-SiO2 -Ti-Au multi-layer structures are accounted in considering the confinement of TE and TM modes in the ETR with the metal layers. The reflectivity will greatly reduce with a Ti layer between SiO2 and Au for light rays with incident angle less than 30 especially for the TE mode, even the thickness of the Ti layer is only 10 nm. If the ETR is laterally confined by SiO2-Au layers without the Ti layer, the Fabry-Perot type modes with an incident angle of zero on one side of the ETR can also have high Q-factor. The FDTD simulation for the ETR confined by metal layers verifies the above analysis based on multi-layer reflections. The output spectra with mode intervals of whispering-gallery modes and Fabry-Perot type modes are observed from different ETR lasers with side length of 10 m, respectively.

Analysis of mode characteristics for equilateral triangle semiconductor microlasers with imperfect boundaries

The influence of imperfect boundaries on the mode quality factor is investigated for equilateral-triangle-resonator (ETR) semiconductor microlasers by the finite difference time domain technique and the Pade approximation with Baker's algorithm. For 2-D ETR with a refractive index of 3.2 and side length of 5 mum, the confined modes can still have a quality factor of about 1000 as small triangles with side length of 1 mum are cut from the vertices of the ETR. For a deformed 5 mum ETR with round vertices and curve sides, the simulated mode quality factors are comparable to the measured results.

Mode quality factor based on far-field emission for square resonators

The quality factors of modes in square resonators are calculated based on the far-field emission of the analytical field distribution. The obtained quality factors are in reasonable agreement with those calculated by the finite-difference time-domain (FDTD) technique and Pade approximation method. The emission power in the square diagonal directions for whispering-gallery-like modes in square resonators is zero due to the interference cancellation caused by the odd field distributions relative to the diagonal mirror planes, so they have larger quality factors than the modes with even field distribution.