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.

Investigations of secondary organic aerosol formation in laboratory chambers

Secondary organic aerosol (SOA) is produced in the atmosphere by oxidation of volatile organic compounds. Laboratory chambers are used understand the formation mechanisms and evolution of SOA formed under controlled conditions. This thesis presents studies of SOA formed from anthropogenic and biogenic precursors and discusses the effects of chamber walls on suspended vapors and particles.

During a chamber experiment, suspended vapors and particles can interact with the chamber walls. Particle wall loss is relatively well-understood, but vapor wall losses have received little study. Vapor wall loss of 2,3-epoxy-1,4-butanediol (BEPOX) and glyoxal was identified, quantified, and found to depend on chamber age and relative humidity.

Particles reside in the atmosphere for a week or more and can evolve chemically during that time period, a process termed aging. Simulating aging in laboratory chambers has proven to be challenging. A protocol was developed to extend the duration of a chamber experiment to 36 h of oxidation and was used to evaluate aging of SOA produced from m-xylene. Total SOA mass concentration increased and then decreased with increasing photooxidation suggesting a transition from functionalization to fragmentation chemistry driven by photochemical processes. SOA oxidation, measured as the bulk particle elemental oxygen-to-carbon ratio and fraction of organic mass at m/z 44, increased continuously starting after 5 h of photooxidation.

The physical state and chemical composition of an organic aerosol affect the mixing of aerosol components and its interactions with condensing species. A laboratory chamber protocol was developed to evaluate the mixing of SOA produced sequentially from two different sources by heating the chamber to induce particle evaporation. Using this protocol, SOA produced from toluene was found to be less volatile than that produced from a-pinene. When the two types of SOA were formed sequentially, the evaporation behavior most closely represented that of SOA from the second parent hydrocarbon, suggesting that the structure of the mixed SOA particles resembles a core of SOA from the first precursor coated by a layer of SOA from the second precursor, indicative of limiting mixing.

Silicide formation and the interaction of metals with polycrystalline Si

The main factors affecting solid-phase Si-metal interactions are reported in this work. The influence of the orientation of the Si substrates and the presence of impurities in metal films and at the Si-metal interface on the formation of nickel and chromium silicides have been demonstrated. We have observed that the formation and kinetic rate of growth of nickel silicides is strongly dependent on the orientation and crystallinity of the Si substrates; a fact which, up to date, has never been seriously investigated in silicide formation. Impurity contaminations in the Cr film and at the Si-Cr interface are the most dominant influencing factors in the formation and kinetic rate of growth of CrSi₂. The potentiality and use of silicides as a diffusion barrier in metallization on silicon devices were also investigated.

Two phases, Ni₂Si and NiSi, form simultaneously in two distinct sublayers in the reaction of Ni with amorphous Si, while only the former phase was observed on other substrates. On (111) oriented Si substrates the growth rate is about 2 to 3 times less than that on <100> or polycrystalline Si. Transmission electron micrographs establish-·that silicide layers grown on different substrates have different microcrystalline structures. The concept of grain-boundary diffusion is speculated to be an important factor in silicide formation.

The composition and kinetic rate of CrSi₂ formation are not influenced by the underlying Si substrate. While the orientation of the Si substrate does not affect the formation of CrSi₂ , the purity of the Cr film and the state of Si-Cr interface become the predominant factors in the reaction process. With an interposed layer of Pd₂Si between the Cr film and the Si substrate, CrSi₂ starts to form at a much lower temperature (400°C) relative to the Si-Cr system. However, the growth rate of CrSi₂ is observed to be independent of the thickness of the Pd2Si layer. For both Si-Cr and Si-Pd₂Si-Cr samples, the growth rate is linear with time with an activation energy of 1.7 ± 0.1 ev.

A tracer technique using radioactive ³¹Si (T_1/2 = 2.26 h) was used to study the formation of CrSi₂ on Pd₂Si. It is established from this experiment that the growth of CrSi₂ takes place partly by transport of Si directly from the Si substrate and partly by breaking Pd₂Si bonds, making free Si atoms available for the growth process.

The role of CrSi₂ in Pd-Al metallization on Si was studied. It is established that a thin CrSi₂ layer can be used as a diffusion barrier to prevent Al from interacting with Pd₂Si in the Pd-Al metallization on Si.

As a generalization of what has been observed for polycrystalline-Si-Al interaction, the reactions between polycrystalline Si (poly Si) and other metals were studied. The metals investigated include Ni, Cr, Pd, Ag and Au. For Ni, Cr and Pd, annealing results in silicide formation, at temperatures similar to those observed on single crystal Si substrates. For Al, Ag and Au, which form simple eutectics with Si annealing results in erosion of the poly Si layer and growth of Si crystallites in the metal films.

Backscattering spectrometry with 2.0 and 2.3 MeV ⁴He ions was the main analytical tool used in all our investigations. Other experimental techniques include the Read camera glancing angle x-ray diffraction, scanning electron, optical and transmission electron microscopy. Details of these analytical techniques are given in Chapter II.

Part I: The abortive infection of bacteriophage øx174 at low temperatures. Part II: The early stages in the process of infection by bacteriophage øx174

Part I

The infection of E. coli by ΦX174 at 15°C is abortive; the cells are killed by the infection but neither mature phage nor SS (single-stranded) DNA are synthesized. Parental RF (replicative form) is formed and subsequently replicated at 15°C. The RF made at 15°C shows normal infectivity and full competence to act as precursor to progeny SS DNA after an increase in temperature to 37°C. The investigations suggest that all of the proteins required for SS DNA synthesis and phage maturation are present in the abortive infection at 15°C.

Three possible causes are suggested for the abortive infection at 15°C: (a) A virus-coded protein whose role is essential to the infection is made at 15°C and assumes its native conformation, but its rate of activity is too low at this temperature to sustain the infection process. (b) Virus maturation may involve the formation of a DNA-protein complex and conformational changes which have an energy threshold infrequently reached at 15°C. (c) A host-coded protein present in uninfected cells, and whose activity is essential to the infection at all temperatures, but not to the host at 15°C, is inactive at 15°C. An hypothesis of this type is offered which proposes that the temperature-limiting factor in SS DNA synthesis in vivo may reflect a temperature-dependent property of the host DNA polymerase.

Part II

Three distinct stages are demonstrated in the process whereby ΦX174 invades its host: (1) Attachment: The phage attach to the cell in a manner that does not irreversibly alter the phage particle and which exhibits "single-hit" kinetics. The total charge on the phage particle is demonstrated to be important in determining the rate at which stable attachment is effected. The proteins specified by ΦX cistrons II, III and VII play roles, which may be indirect, in the attachment reaction. (2) Eclipse: 'The attached phage undergo a conformational change. Some of the altered phage particles spontaneously detach from the cell (in a non-infective form) while the remainder are more tightly bound to the cell. The altered phage particles detached (spontaneously or chemically) from such complexes have at least 40% of their DNA extruded from the phage coat. It is proposed that this particle is, or derives from, a direct intermediate in the penetration of the viral DNA.

The kinetics for the eclipse of attached phage particles are first-order with respect to phage concentration and biphasic; about 85% of the phage eclipse at one rate (k = 0.86 min^-1) and the remainder do so at a distinctly lesser rate (k = 0.21 min^-1).

The eclipse event is very temperature-dependent and has the relatively high Arrhenius activation energy of 36.6 kcal/mole, indicating the cooperative nature of the process. The temperature threshold for eclipse is 17 to 18°C.

At present no specific ΦX cistron is identified as affecting the eclipse process. (3) DNA penetration: A fraction of the attached, eclipsed phage particles corresponding in number to the plaque-forming units complete DNA penetration. The penetrated DNA is found in the cell as RF, and the empty phage protein coat remains firmly attached to the exterior of the cell. This step is inhibited by prior irradiation of the phage with relatively high doses of UV light and is insensitive to the presence of KCN and NaN₃. Temporally excluded superinfecting phages do not achieve DNA penetration.

Both eclipsed phage particles and empty phage protein coats may be dissociated from infected cells; some of their properties are described.