QUANTITATIVE THREE DIMENSIONAL ELASTICITY IMAGING

Neoplastic tissue is typically highly vascularized, contains abnormal concentrations of extracellular proteins (e.g. collagen, proteoglycans) and has a high interstitial fluid pres- sure compared to most normal tissues. These changes result in an overall stiffening typical of most solid tumors. Elasticity Imaging (EI) is a technique which uses imaging systems to measure relative tissue deformation and thus noninvasively infer its mechanical stiffness. Stiffness is recovered from measured deformation by using an appropriate mathematical model and solving an inverse problem. The integration of EI with existing imaging modal- ities can improve their diagnostic and research capabilities. The aim of this work is to develop and evaluate techniques to image and quantify the mechanical properties of soft tissues in three dimensions (3D). To that end, this thesis presents and validates a method by which three dimensional ultrasound images can be used to image and quantify the shear modulus distribution of tissue mimicking phantoms. This work is presented to motivate and justify the use of this elasticity imaging technique in a clinical breast cancer screening study. The imaging methodologies discussed are intended to improve the specificity of mammography practices in general. During the development of these techniques, several issues concerning the accuracy and uniqueness of the result were elucidated. Two new algorithms for 3D EI are designed and characterized in this thesis. The first provides three dimensional motion estimates from ultrasound images of the deforming ma- terial. The novel features include finite element interpolation of the displacement field, inclusion of prior information and the ability to enforce physical constraints. The roles of regularization, mesh resolution and an incompressibility constraint on the accuracy of the measured deformation is quantified. The estimated signal to noise ratio of the measured displacement fields are approximately 1800, 21 and 41 for the axial, lateral and eleva- tional components, respectively. The second algorithm recovers the shear elastic modulus distribution of the deforming material by efficiently solving the three dimensional inverse problem as an optimization problem. This method utilizes finite element interpolations, the adjoint method to evaluate the gradient and a quasi-Newton BFGS method for optimiza- tion. Its novel features include the use of the adjoint method and TVD regularization with piece-wise constant interpolation. A source of non-uniqueness in this inverse problem is identified theoretically, demonstrated computationally, explained physically and overcome practically. Both algorithms were test on ultrasound data of independently characterized tissue mimicking phantoms. The recovered elastic modulus was in all cases within 35% of the reference elastic contrast. Finally, the preliminary application of these techniques to tomosynthesis images showed the feasiblity of imaging an elastic inclusion.

Correlations in low dimensional quantum systems

In this thesis I present the work done during my PhD in the area of low dimensional quantum gases. The chapters of this thesis are self contained and represent individual projects which have been peer reviewed and accepted for publication in respected international journals. Various systems are considered, the first of which is a two particle model which possesses an exact analytical solution. I investigate the non-classical correlations that exist between the particles as a function of the tunable properties of the system. In the second work I consider the coherences and out of equilibrium dynamics of a one-dimensional Tonks-Girardeau gas. I show how the coherence of the gas can be inferred from various properties of the reduced state and how this may be observed in experiments. I then present a model which can be used to probe a one-dimensional Fermi gas by performing a measurement on an impurity which interacts with the gas. I show how this system can be used to observe the so-called orthogonality catastrophe using modern interferometry techniques. In the next chapter I present a simple scheme to create superposition states of particles with special emphasis on the NOON state. I explore the effect of inter-particle interactions in the process and then characterise the usefulness of these states for interferometry. Finally I present my contribution to a project on long distance entanglement generation in ion chains. I show how carefully tuning the environment can create decoherence-free subspaces which allows one to create and preserve entanglement.

Fabrication, characterisation and electroanalysis at 1-dimensional nanostructures

Integrated nanowire electrodes that permit direct, sensitive and rapid electrochemical based detection of chemical and biological species are a powerful emerging class of sensor devices. As critical dimensions of the electrodes enter the nanoscale, radial analyte diffusion profiles to the electrode dominate with a corresponding enhancement in mass transport, steady-state sigmoidal voltammograms, low depletion of target molecules and faster analysis. To optimise these sensors it is necessary to fully understand the factors that influence performance limits including: electrode geometry, electrode dimensions, electrode separation distances (within nanowire arrays) and diffusional mass transport. Therefore, in this thesis, theoretical simulations of analyte diffusion occurring at a variety of electrode designs were undertaken using Comsol Multiphysics®. Sensor devices were fabricated and corresponding experiments were performed to challenge simulation results. Two approaches for the fabrication and integration of metal nanowire electrodes are presented: Template Electrodeposition and Electron-Beam Lithography. These approaches allow for the fabrication of nanowires which may be subsequently integrated at silicon chip substrates to form fully functional electrochemical devices. Simulated and experimental results were found to be in excellent agreement validating the simulation model. The electrochemical characteristics exhibited by nanowire electrodes fabricated by electronbeam lithography were directly compared against electrochemical performance of a commercial ultra-microdisc electrode. Steady-state cyclic voltammograms in ferrocenemonocarboxylic acid at single ultra-microdisc electrodes were observed at low to medium scan rates (≤ 500 mV.s-1). At nanowires, steady-state responses were observed at ultra-high scan rates (up to 50,000 mV.s-1), thus allowing for much faster analysis (20 ms). Approaches for elucidating faradaic signal without the requirement for background subtraction were also developed. Furthermore, diffusional process occurring at arrays with increasing inter-electrode distance and increasing number of nanowires were explored. Diffusion profiles existing at nanowire arrays were simulated with Comsol Multiphysics®. A range of scan rates were modelled, and experiments were undertaken at 5,000 mV.s-1 since this allows rapid data capture required for, e.g., biomedical, environmental and pharmaceutical diagnostic applications.

Fiber-optic interferometric two-dimensional scattering-measurement system.

We present a fiber-optic interferometric system for measuring depth-resolved scattering in two angular dimensions using Fourier-domain low-coherence interferometry. The system is a unique hybrid of the Michelson and Sagnac interferometer topologies. The collection arm of the interferometer is scanned in two dimensions to detect angular scattering from the sample, which can then be analyzed to determine the structure of the scatterers. A key feature of the system is the full control of polarization of both the illumination and the collection fields, allowing for polarization-sensitive detection, which is essential for two-dimensional angular measurements. System performance is demonstrated using a double-layer microsphere phantom. Experimental data from samples with different sizes and acquired with different polarizations show excellent agreement with Mie theory, producing structural measurements with subwavelength accuracy.

3D refraction correction and extraction of clinical parameters from spectral domain optical coherence tomography of the cornea.

Capable of three-dimensional imaging of the cornea with micrometer-scale resolution, spectral domain-optical coherence tomography (SDOCT) offers potential advantages over Placido ring and Scheimpflug photography based systems for accurate extraction of quantitative keratometric parameters. In this work, an SDOCT scanning protocol and motion correction algorithm were implemented to minimize the effects of patient motion during data acquisition. Procedures are described for correction of image data artifacts resulting from 3D refraction of SDOCT light in the cornea and from non-idealities of the scanning system geometry performed as a pre-requisite for accurate parameter extraction. Zernike polynomial 3D reconstruction and a recursive half searching algorithm (RHSA) were implemented to extract clinical keratometric parameters including anterior and posterior radii of curvature, central cornea optical power, central corneal thickness, and thickness maps of the cornea. Accuracy and repeatability of the extracted parameters obtained using a commercial 859nm SDOCT retinal imaging system with a corneal adapter were assessed using a rigid gas permeable (RGP) contact lens as a phantom target. Extraction of these parameters was performed in vivo in 3 patients and compared to commercial Placido topography and Scheimpflug photography systems. The repeatability of SDOCT central corneal power measured in vivo was 0.18 Diopters, and the difference observed between the systems averaged 0.1 Diopters between SDOCT and Scheimpflug photography, and 0.6 Diopters between SDOCT and Placido topography.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.

Parity-violating vertices for spin-3 gauge fields

The problem of constructing consistent parity-violating interactions for spin-3 gauge fields is considered in Minkowski space. Under the assumptions of locality, Poincaré invariance, and parity noninvariance, we classify all the nontrivial perturbative deformations of the Abelian gauge algebra. In space-time dimensions n=3 and n=5, deformations of the free theory are obtained which make the gauge algebra non-Abelian and give rise to nontrivial cubic vertices in the Lagrangian, at first order in the deformation parameter g. At second order in g, consistency conditions are obtained which the five-dimensional vertex obeys, but which rule out the n=3 candidate. Moreover, in the five-dimensional first-order deformation case, the gauge transformations are modified by a new term which involves the second de Wit-Freedman connection in a simple and suggestive way. © 2006 The American Physical Society.

Experimental investigation of thermal structures in regular three-dimensional falling films

Interfacial waves on the surface of a falling liquid film are known to modify heat and mass transfer. Under non-isothermal conditions, the wave topology is strongly influenced by the presence of thermocapillary (Marangoni) forces at the interface which leads to a destabilization of the film flow and potentially to critical film thinning. In this context, the present study investigates the evolution of the surface topology and the evolution of the surface temperature for the case of regularly excited solitary-type waves on a falling liquid film under the influence of a wall-side heat flux. Combining film thickness (chromatic confocal imaging) and surface temperature information (infrared thermography), interactions between hydrodynamics and thermocapillary forces are revealed. These include the formation of rivulets, film thinning and wave number doubling in spanwise direction. Distinct thermal structures on the films’ surface can be associated to characteristics of the surface topology.

A multi-scale atomistic-continuum modelling of crack propagation in a two-dimensional macroscopic plate

A novel multi-scale seamless model of brittle-crack propagation is proposed and applied to the simulation of fracture growth in a two-dimensional Ag plate with macroscopic dimensions. The model represents the crack propagation at the macroscopic scale as the drift-diffusion motion of the crack tip alone. The diffusive motion is associated with the crack-tip coordinates in the position space, and reflects the oscillations observed in the crack velocity following its critical value. The model couples the crack dynamics at the macroscales and nanoscales via an intermediate mesoscale continuum. The finite-element method is employed to make the transition from the macroscale to the nanoscale by computing the continuum-based displacements of the atoms at the boundary of an atomic lattice embedded within the plate and surrounding the tip. Molecular dynamics (MD) simulation then drives the crack tip forward, producing the tip critical velocity and its diffusion constant. These are then used in the Ito stochastic calculus to make the reverse transition from the nanoscale back to the macroscale. The MD-level modelling is based on the use of a many-body potential. The model successfully reproduces the crack-velocity oscillations, roughening transitions of the crack surfaces, as well as the macroscopic crack trajectory. The implications for a 3-D modelling are discussed.

Details of an integrated approach to three-dimensional dynamic fluid structure interaction

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. Numerical modelling of dynamic fluid-structure interaction (DFSI) involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge and until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. A single, finite volume unstructured mesh (FV-UM) spatial discretisation method has been employed on a single mesh for the entire domain. The Navier Stokes equations for fluid flow are solved using a SIMPLE type procedure and the Newmark b algorithm is employed for solving the dynamic equilibrium equations for linear elastic solid mechanics and mesh movement is achieved using a spring based mesh procedure for dynamic mesh movement. In the paper we describe a number of additional computation issues for the efficient and accurate modelling of three-dimensional, dynamic fluid-structure interaction problems.

Central nervous system drugs: Low temperature X-ray structures of (I) the base 5-(2,3-dichlorophenyl)-2,4-diamino-6-fluoromethyl-pyrimidine (4030W92) and (II) 5-(2,6-dichlorophenyl)-1-H-2,4-diamino-6-methyl-pyrimidine methanesulphonic acid salt (227C89)

The X-ray crystal structures of (I), the base 4030W92, 5-(2,3-dichlorophenyl)-2,4-diamino-6-fluoromethyl-pyrimidine, C11H9Cl2FN4, and (II) 227C89, the methanesulphonic acid salt of 5-(2,6-dichlorophenyl)-1-H-2,4-diamino-6-methyl-pyrimidine, C11H11Cl2N4 center dot CH3O3S, have been carried out at low temperature. A detailed comparison of the two structures is given. Structure (I) is non-centrosymmetric, crystallizing in space group P2(1) with unit cell a = 10.821(3), b = 8.290(3), c = 13.819(4) angstrom, beta = 105.980(6)degrees, V = 1191.8(6) angstrom(3), Z = 4 (two molecules per asymmetric unit) and density (calculated) = 1.600 mg/m(3). Structure (II) crystallizes in the triclinic space group P (1) over bar with unit cell a = 7.686(2), b = 8.233(2), c = 12.234(2) angstrom, alpha = 78.379(4), beta = 87.195(4), gamma = 86.811(4)degrees, V = 756.6(2) angstrom(3), Z = 2, density (calculated) = 1.603 mg/m(3). Final R indices [I > 2sigma(I)] are R1 = 0.0572, wR2 = 0.1003 for (I) and R1 = 0.0558, wR2 = 0.0982 for (II). R indices (all data) are R1 = 0.0983, wR2 = 0.1116 for (I) and R1 = 0.1009, wR2 = 0.1117 for (II). 5- Phenyl-2,4 diaminopyrimidine and 6-phenyl-1,2,4 triazine derivatives, which include lamotrigine (3,5-diamino-6-(2,3-dichlorophenyl)-1,2,4-triazine), have been investigated for some time for their effects on the central nervous system. The three dimensional structures reported here form part of a newly developed data base for the detailed investigation of members of this structural series and their biological activities.

Low temperature X-ray structures of two crystalline forms (I) and (II) of the pyrimidine derivative and sodium channel blocker BW1003C87: 5-(2,3,5-trichlorophenyl)-2,4-diaminopyrimidine

The X-ray crystal structures of two crystalline forms of 5-(2,3,5-trichlorophenyl)-2,4-diaminopyrimidine, C10H7Cl3N4 (code name BW1003C87) (I) and (II), have been carried out at liquid nitrogen temperature. A detailed comparison of the two structures is given. Both are centrosymmetric, with structure (I) in the triclinic space group P (1) over bar unit cell a = 6.4870(10), b = 9.216(2), c = 12.016(2) angstrom, alpha = 75.78(3)degrees, beta = 89.95(3)degrees, gamma = 83.45(3)degrees, V = 691.5(2) angstrom(3), Z = 2 and density (calculated) = 1.544 Mg/m(3); and (II) in the monoclinic space group P2(1)/c, unit cell a = 12.000(2), b = 7.518(2), c = 13.450(3) angstrom, beta = 97.87(3)degrees, V = 1202.0(5) angstrom(3), Z = 4, Density (calculated) = 1.600 Mg/m(3). Structure (I) includes a solvated CH3OH in the lattice. Final R indices [I > 2sigma(I)] are R1 = 0.0427, wR2 = 0.1075 for (I) and R1 = 0.0487, wR2 = 0.1222 for (II). R indices (all data) are R1 = 0.0470, wR2 = 0.1118 for (I) and R1 = 0.0623, wR2 = 0.1299 for (II). 5-Phenyl-2,4 diaminopyrimidine and 6-phenyl-1,2,4 triazine derivatives, which include lamotrigine (3,5-diamino-6-(2,3-dichlorophenyl)-1,2,4-triazine), have been investigated for some time for their effects on the central nervous system. Both lamotrigine and 5-(2,3,5-trichlorophenyl)-2,4-diaminopyrimidine (code name BW1003C87), the subject of the present study, are anticonvulsant as well as neuroprotective in models of brain ischaemia and in a model of white matter ischaemia. BW1003C87 is a sodium channel blocker which also reduces the release of the neurotransmitter glutamate. The three dimensional structures reported here form part of a newly developed data base for the detailed investigation of members of this drug family and their biological activities.

Exploring the Red Sea seasonal ecosystem functioning using a three-dimensional biophysical model

The Red Sea exhibits complex hydrodynamic and biogeochemical dynamics, which vary both in time and space. These dynamics have been explored through the development and application of a 3-D ecosystem model. The simulation system comprises two off-line coupled submodels: the MIT General Circulation Model (MITgcm) and the European Regional Seas Ecosystem Model (ERSEM), both adapted for the Red Sea. The results from an annual simulation under climatological forcing are presented. Simulation results are in good agreement with satellite and in situ data illustrating the role of the physical processes in determining the evolution and variability of the Red Sea ecosystem. The model was able to reproduce the main features of the Red Sea ecosystem functioning, including the exchange with the Gulf of Aden, which is a major driving mechanism for the whole Red Sea ecosystem and the winter overturning taking place in the north. Some model limitations, mainly related to the dynamics of the extended reef system located in the southern part of the Red Sea, which is not currently represented in the model, still need to be addressed.

Globcurrent: Sentinel-3 synergy in action

The ESA Data User Element (DUE) funded GlobCurrent project (http://www.globcurrent.org) aims to: (i) advance the quantitative estimation of ocean surface currents from satellite sensor synergy; and (ii) demonstrate impact in user-led scientific, operational and commercial applications that, in turn, will improve and strengthen the uptake of satellite measurements. Today, a synergetic approach for quantitative analysis can build on high-resolution imaging radar and spectrometer data, infrared radiometer data and radar altimeter measurements. It will further integrate Sentinel-3 in combination with Sentinel-1 SAR data. From existing and past missions, it is often demonstrated that sharp gradients in the sea surface temperature (SST) field and the ocean surface chlorophyll-a distribution are spatially correlated with the sea surface roughness anomaly fields at small spatial scales, in the sub-mesocale (1-10 km) to the mesoscale (30-80 km). At the larger mesoscale range (>50 km), information derived from radar altimeters often depict the presence of coherent structures and eddies. The variability often appears largest in regions where the intense surface current regimes (>100 - 200 km) are found. These 2-dimensional structures manifested in the satellite observations represent evidence of the upper ocean (~100-200 m) dynamics. Whereas the quasi geostrophic assumption is valid for the upper ocean dynamics at the larger scale (>100 km), possible triggering mechanisms for the expressions at the mesoscale-to-submesoscale may include spiraling tracers of inertial motion and the interaction of the wind-driven Ekman layer with the quasi-geostrophic current field. This latter, in turn, produces bands of downwelling (convergence) and upwelling (divergence) near fronts. A regular utilization of the sensor synergy approach with the combination of Sentinel-3 and Sentinel-1 will provide a highly valuable data set for further research and development to better relate the 2-dimensional surface expressions and the upper ocean dynamics.