Enlarging the depth of focus by filtering in the phase-space domain

The ambiguity function was employed as a merit function to design an optical system with a high depth of focus. The ambiguity function with the desired enlarged-depth-of-focus characteristics was obtained by using a properly designed joint filter to modify the ambiguity function of the original pupil in the phase-space domain. From the viewpoint of the filter theory, we roughly propose that the constraints of the spatial filters that are used to enlarge the focal depth must be satisfied. These constraints coincide with those that appeared in the previous literature on this topic. Following our design procedure, several sets of apodizers were synthesized, and their performances in the defocused imagery were compared with each other and with other previous designs. (c) 2005 Optical Society of America.

Electrochromism accompanying ferroelectric domain inversion in congruent RuO2 : LiNbO3 crystal

Electrochromic phenomena accompanying the ferroelectric domain inversion in congruent RuO2-doped z-cut LiNbO3 crystals at room temperature are observed in experiments. During the electric poling process, the electrochromism accompanies the ferroelectric domain inversion simultaneously in the same poled area. The electrochromism is completely reversible when the domain is inverted from the reverse direction. The influences of electric field and annealing conditions on domain inversion and electrochromism are also discussed. We propose the reasonable assumption that charge redistribution within the crystal structure caused by domain inversion is the source for electrochemically oxidation and reduction of Ru ion to produce the electrochromic effect. (c) 2005 Optical Society of America.

Reversible electrochromic effect accompanying domain-inversion in LiNbO3 : Ru : Fe crystals

A reversible electrochromic effect during the electrical poling process in LiNbO3:Ru:Fe crystals at room temperature is observed. In electrode area, both electrochromism and domain-inversion occur mutually and electrochromism is reversible during back-switch poling, which are experimentally verified, and a microstructure model to explain the mechanism is proposed. In addition, different from the undoped LiNbO3 crystals, the breakdown field (> 25.0 kV/mm) is higher than the coercive (21.0 kV/mm) measured in LiNbO3:Ru:Fe, which proves a possible new technique to realize domain-inversion by constant electric field rather than pulsed one. (c) 2005 American Institute of Physics.

Symmetric representation of an integral domain over a subdomain

Let F(θ) be a separable extension of degree n of a field F. Let Δ and D be integral domains with quotient fields F(θ) and F respectively. Assume that Δ ᴝ D. A mapping φ of Δ into the n x n D matrices is called a Δ/D rep if (i) it is a ring isomorphism and (ii) it maps d onto dI_n whenever d ϵ D. If the matrices are also symmetric, φ is a Δ/D symrep.

Every Δ/D rep can be extended uniquely to an F(θ)/F rep. This extension is completely determined by the image of θ. Two Δ/D reps are called equivalent if the images of θ differ by a D unimodular similarity. There is a one-to-one correspondence between classes of Δ/D reps and classes of Δ ideals having an n element basis over D.

The condition that a given Δ/D rep class contain a Δ/D symrep can be phrased in various ways. Using these formulations it is possible to (i) bound the number of symreps in a given class, (ii) count the number of symreps if F is finite, (iii) establish the existence of an F(θ)/F symrep when n is odd, F is an algebraic number field, and F(θ) is totally real if F is formally real (for n = 3 see Sapiro, “Characteristic polynomials of symmetric matrices” Sibirsk. Mat. Ž. 3 (1962) pp. 280-291), and (iv) study the case D = Z, the integers (see Taussky, “On matrix classes corresponding to an ideal and its inverse” Illinois J. Math. 1 (1957) pp. 108-113 and Faddeev, “On the characteristic equations of rational symmetric matrices” Dokl. Akad. Nauk SSSR 58 (1947) pp. 753-754).

The case D = Z and n = 2 is studied in detail. Let Δ’ be an integral domain also having quotient field F(θ) and such that Δ’ ᴝ Δ. Let φ be a Δ/Z symrep. A method is given for finding a Δ’/Z symrep ʘ such that the Δ’ ideal class corresponding to the class of ʘ is an extension to Δ’ of the Δ ideal class corresponding to the class of φ. The problem of finding all Δ/Z symreps equivalent to a given one is studied.

Quantitative measurement of domain inversion in RuO2 : LiNbO3 crystal by digital holographic interferometry

We quantitatively study the domain inversion in a RuO2:LiNbO3 crystal wafer by the digital holographic interferometry. The crystal wafer is placed into one arm of a Mach-Zehnder-type interferometer to record a series of holograms. Making use of the angular spectrum backward propagation algorithm, we reconstruct the optical wave field in the crystal plane. The extracted phase difference from the reconstructed optical wave field is a well linear function of the applied external voltage. We deduce that the linear electro-optic coefficient of the detected RuO2:LiNbO3 crystal sample is 9.1x10(-12) m/V. An unexpected phase contrast at the antiparallel domain wall is observed and the influence of the applied external voltage on it is studied in detail. Also the built-in internal field is quantitatively measured as 0.72 kV/mm. (c) 2006 American Institute of Physics.

Ridge-shape phase distribution adjacent to 180 degrees domain wall in congruent LiNbO3 crystal

The digital holographic interferometry is used in the dynamic and static measurements of phase variation induced by domain inversion. For the first time, to the authors' knowledge, they observe the existence of ridge-shape phase distribution adjacent to 180 degrees domain wall in congruent LiNbO3 crystal. During the domain wall motion, the phase variations are not uniform but have obvious relaxations. In the static measurement, the ridge elevation can vary linearly with the uniform electric field. The reasonable assumptions are proposed to explain these effects. (c) 2006 American Institute of Physics.

Time domain synthesis of electro-optical birefringent system for optical waveform generation

In this paper is described a novel technique for producing an electro-optical intensity synthesizer which can generate different periodic time domain waveforms through only sine or cosine wave applied-voltages. The synthesizer presented here consists of a series of stages between two polarizers, with each stage consisting of an electro-optic element and a compensator. Every electro-optical element has the same applied-voltage function but different azimuth angles and ratios between the longitudinal and transverse lengths. The main principle is the synthesis of an electro-optic effect and a polarization interference effect in the time domain. This technique is based on an expanded Fourier positive-direction searching algorithm, which can not only simplify the calculation process but also produces many choices of structural parameters for different waveforms generation. A three-stage synthesis of an electro-optical birefringent system for continuous square waveform is undertaken to prove the principle.

Dynamic phase-mapping of domain nucleation on MgO : LiNbO3 crystal by digital holographic interferometry

The quantitative phase-mapping of the domain nucleation in MgO:LiNbO3 crystals is presented by using the digital holographic interferometry. An unexpected peak phase at the beginning of the domain nucleation is observed and it is lowered as the spreading of the domain nucleus. The existence of the nucleus changes the moving speed of the domain wall by pinning it for 3s. Such in-situ quantitative analysis of the domain nucleation process is a key to optimizing domain structure fabrication.

Red laser-induced domain inversion in MgO-doped lithium niobate crystals

A laser beam at wavelength 647 nm is focused on a sample of 5 mol% MgO-doped lithium niobate crystal for domain inversion by a conventional external electric field. In this case, a reduction of 36% in the electric field required for domain nucleation (nucleation field) is observed. To the best of our knowledge, it is the longest wavelength reported for laser-induced domain inversion. This extends the spectrum of laser inducing, and the experimental results are helpful to understand the nucleation dynamics under laser illumination. The dependence of nucleation fields on intensities of laser beams is analysed in experiments.

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.

Digital holographic interferometry and the angular spectrum method applied to measuring the domain inversion in ferroelectric crystal

The application of digital holographic interferometry on the quantitative measurement of the domain inversion in a RuO2: LiNbO3 crystal wafer is presented. The recorded holograms are reconstructed by the angular spectrum method. From the reconstructed phase distribution we can clearly observe the boundary between the inverted and un-inverted domain regions. Comparisons with the results reconstructed by use of the Fresnel transform method are given. Factors that influence the measurement include the spectrum filter size and the spectrum movement are discussed. The spectrum filter size has an effect on the measurement of the details. Although the spectrum movement affects every single reconstructed image, it has no influence on the final measurement.

Wavelength dependence of light-induced domain nucleation in MgO-doped congruent LiNbO3 crystal

Within the wavelength range from 351 to 799 nm, the different reductions of nucleation field induced by the focused continuous laser irradiation are achieved in the 5 mol % MgO-doped congruent LiNbO3 crystals. The reduction proportion increases exponentially with decreasing irradiation wavelength and decreases exponentially with increasing irradiation wavelength. At one given wavelength, the reduction proportion increases exponentially with increasing irradiation intensity. An assumption is proposed that the reduction of nucleation field is directly related to the defect structure of crystal lattice generated by the complex coaction of incident irradiation field and external electric field. (c) 2007 American Institute of Physics.

Dynamic full-range Fourier-domain optical coherence tomography using sinusoidal phase-modulating interferometry

We propose a technique for dynamic full-range Fourier-domain optical coherence tomography by using sinusoidal phase-modulating interferometry, where both the full-range structural information and depth-resolved dynamic information are obtained. A novel frequency-domain filtering algorithm is proposed to reconstruct a time-dependent complex spectral interferogram from the sinusoidally phase-modulated interferogram detected with a high-rate CCD camera. By taking the amplitude and phase of the inverse Fourier transform of the complex spectral interferogram, a time-dependent full-range cross-sectional image and depth-resolved displacement are obtained. Displacement of a sinusoidally vibrating glass cover slip behind a fixed glass cover slip is measured with subwavelength sensitivity to demonstrate the depth-resolved dynamic imaging capability of our system. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

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.