An Accelerated Chow and Liu Algorithm: Fitting Tree Distributions to High Dimensional Sparse Data

Chow and Liu introduced an algorithm for fitting a multivariate distribution with a tree (i.e. a density model that assumes that there are only pairwise dependencies between variables) and that the graph of these dependencies is a spanning tree. The original algorithm is quadratic in the dimesion of the domain, and linear in the number of data points that define the target distribution $P$. This paper shows that for sparse, discrete data, fitting a tree distribution can be done in time and memory that is jointly subquadratic in the number of variables and the size of the data set. The new algorithm, called the acCL algorithm, takes advantage of the sparsity of the data to accelerate the computation of pairwise marginals and the sorting of the resulting mutual informations, achieving speed ups of up to 2-3 orders of magnitude in the experiments.

Computer Recognition of Three-Dimensional Objects in a Visual Scene

Methods are presented (1) to partition or decompose a visual scene into the bodies forming it; (2) to position these bodies in three-dimensional space, by combining two scenes that make a stereoscopic pair; (3) to find the regions or zones of a visual scene that belong to its background; (4) to carry out the isolation of objects in (1) when the input has inaccuracies. Running computer programs implement the methods, and many examples illustrate their behavior. The input is a two-dimensional line-drawing of the scene, assumed to contain three-dimensional bodies possessing flat faces (polyhedra); some of them may be partially occluded. Suggestions are made for extending the work to curved objects. Some comparisons are made with human visual perception. The main conclusion is that it is possible to separate a picture or scene into the constituent objects exclusively on the basis of monocular geometric properties (on the basis of pure form); in fact, successful methods are shown.

The mixing of bubbles in two-dimensional bidisperse foams under extensional shear

Cox, S.J. (2006) The mixing of bubbles in two-dimensional bidisperse foams under extensional shear. Journal of Non-Newtonian Fluid Mechanics . 137:39-45.

Three-dimensional bubble clusters : Shape, packing, and growth rate

Cox, S.J., and Graner, F. (2004) Three-dimensional bubble clusters: shape, packing and growth-rate. Physical review. E, Statistical, nonlinear, and soft matter physics . 69:031409.

Instabilities in two-dimensional flower and chain clusters of bubbles

M.A. Fortes et al., Instabilities in two-dimensional flower and chain clusters of bubbles, Colloids and Surfaces A: Physicochemical and Engineering Aspects Volume 309, Issues 1-3, 1 November 2007, Pages 64-70 A Collection of Papers Presented at the 6th Eufoam Conference, Potsdam, Germany, 2-6 July, 2006

Two-Dimensional?Three-Dimensional Correspondence in Mammography

R. Marti, C. Rubin, E. Denton and R. Zwiggelaar, '2D-3D correspondence in mammography', Cybernetics and Systems 35 (1), 85-105 (2004)

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.

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.