Nonparametric Bayesian Density Modeling with Gaussian Processes

We present the Gaussian process density sampler (GPDS), an exchangeable generative model for use in nonparametric Bayesian density estimation. Samples drawn from the GPDS are consistent with exact, independent samples from a distribution defined by a density that is a transformation of a function drawn from a Gaussian process prior. Our formulation allows us to infer an unknown density from data using Markov chain Monte Carlo, which gives samples from the posterior distribution over density functions and from the predictive distribution on data space. We describe two such MCMC methods. Both methods also allow inference of the hyperparameters of the Gaussian process.

Imaging the femoral cortex: thickness, density and mass from clinical CT.

There is growing evidence that focal thinning of cortical bone in the proximal femur may predispose a hip to fracture. Detecting such defects in clinical CT is challenging, since cortices may be significantly thinner than the imaging system's point spread function. We recently proposed a model-fitting technique to measure sub-millimetre cortices, an ill-posed problem which was regularized by assuming a specific, fixed value for the cortical density. In this paper, we develop the work further by proposing and evaluating a more rigorous method for estimating the constant cortical density, and extend the paradigm to encompass the mapping of cortical mass (mineral mg/cm(2)) in addition to thickness. Density, thickness and mass estimates are evaluated on sixteen cadaveric femurs, with high resolution measurements from a micro-CT scanner providing the gold standard. The results demonstrate robust, accurate measurement of peak cortical density and cortical mass. Cortical thickness errors are confined to regions of thin cortex and are bounded by the extent to which the local density deviates from the peak, averaging 20% for 0.5mm cortex.

Imaging the femoral cortex: Thickness, density and mass from clinical CT

There is growing evidence that focal thinning of cortical bone in the proximal femur may predispose a hip to fracture. Detecting such defects in clinical CT is challenging, since cortices may be significantly thinner than the imaging system's point spread function. We recently proposed a model-fitting technique to measure sub-millimetre cortices, an ill-posed problem which was regularized by assuming a specific, fixed value for the cortical density. In this paper, we develop the work further by proposing and evaluating a more rigorous method for estimating the constant cortical density, and extend the paradigm to encompass the mapping of cortical mass (mineral mg/cm 2) in addition to thickness. Density, thickness and mass estimates are evaluated on sixteen cadaveric femurs, with high resolution measurements from a micro-CT scanner providing the gold standard. The results demonstrate robust, accurate measurement of peak cortical density and cortical mass. Cortical thickness errors are confined to regions of thin cortex and are bounded by the extent to which the local density deviates from the peak, averaging 20% for 0.5mm cortex. © 2012 Elsevier B.V.

Density stratified environments: The double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the 'double-tank' method proposed by Oster (Sci Am 213:70-76, 1965). We refer to Oster's method as the 'forced-drain' approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the 'free-drain' approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the 'free-drain' extension of Oster's method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established. © 2008 Springer-Verlag.

Random function priors for exchangeable arrays with applications to graphs and relational data

A fundamental problem in the analysis of structured relational data like graphs, networks, databases, and matrices is to extract a summary of the common structure underlying relations between individual entities. Relational data are typically encoded in the form of arrays; invariance to the ordering of rows and columns corresponds to exchangeable arrays. Results in probability theory due to Aldous, Hoover and Kallenberg show that exchangeable arrays can be represented in terms of a random measurable function which constitutes the natural model parameter in a Bayesian model. We obtain a flexible yet simple Bayesian nonparametric model by placing a Gaussian process prior on the parameter function. Efficient inference utilises elliptical slice sampling combined with a random sparse approximation to the Gaussian process. We demonstrate applications of the model to network data and clarify its relation to models in the literature, several of which emerge as special cases.

Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band-gap edge: An investigation into the role of the density of photon states near resonance

In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e., at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called "mirrorless"cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in two-dimensional photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first-order differential equations, solved using the Green's-function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data. © 2013 American Physical Society.

Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band-gap edge: an investigation into the role of the density of photon states near resonance.

In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e., at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called "mirrorless"cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in two-dimensional photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first-order differential equations, solved using the Green's-function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data.