An external field prior for the hidden Potts model with application to cone-beam computed tomography

In images with low contrast-to-noise ratio (CNR), the information gain from the observed pixel values can be insufficient to distinguish foreground objects. A Bayesian approach to this problem is to incorporate prior information about the objects into a statistical model. A method for representing spatial prior information as an external field in a hidden Potts model is introduced. This prior distribution over the latent pixel labels is a mixture of Gaussian fields, centred on the positions of the objects at a previous point in time. It is particularly applicable in longitudinal imaging studies, where the manual segmentation of one image can be used as a prior for automatic segmentation of subsequent images. The method is demonstrated by application to cone-beam computed tomography (CT), an imaging modality that exhibits distortions in pixel values due to X-ray scatter. The external field prior results in a substantial improvement in segmentation accuracy, reducing the mean pixel misclassification rate for an electron density phantom from 87% to 6%. The method is also applied to radiotherapy patient data, demonstrating how to derive the external field prior in a clinical context.

A meshless method based on Point Interpolation Method (PIM) for the space fractional diffusion equation

This paper aims to develop a meshless approach based on the Point Interpolation Method (PIM) for numerical simulation of a space fractional diffusion equation. Two fully-discrete schemes for the one-dimensional space fractional diffusion equation are obtained by using the PIM and the strong-forms of the space diffusion equation. Numerical examples with different nodal distributions are studied to validate and investigate the accuracy and efficiency of the newly developed meshless approach.

Investigation into chemistry of new particle formation and growth in subtropical urban environment

The role of different chemical compounds, particularly organics, involved in the new particle formation (NPF) and its consequent growth are not fully understood. Therefore, this study was conducted to investigate the chemistry of aerosol particles during NPF events in an urban subtropical environment. Aerosol chemical composition was measured along with particle number size distribution (PNSD) and several other air quality parameters at five sites across an urban subtropical environment. An Aerodyne compact Time-of-Flight Aerosol Mass Spectrometer (c-TOF-AMS) and a TSI Scanning Mobility Particle Sizer (SMPS) measured aerosol chemical composition and PNSD, respectively. Five NPF events, with growth rates in the range 3.3-4.6 nm, were detected at two sites. The NPF events happened on relatively warmer days with lower humidity and higher solar radiation. Temporal percent fractions of nitrate, sulphate, ammonium and organics were modelled using the Generalised Additive Model (GAM), with a basis of penalised spline. Percent fractions of organics increased after the NPF events, while the mass fraction of ammonium and sulphate decreased. This uncovered the important role of organics in the growth of newly formed particles. Three organic markers, factors f43, f44 and f57, were calculated and the f44 vs f43 trends were compared between nucleation and non-nucleation days. f44 vs f43 followed a different pattern on nucleation days compared to non-nucleation days, whereby f43 decreased for vehicle emission generated particles, while both f44 and f43 decreased for NPF generated particles. It was found for the first time that vehicle generated and newly formed particles cluster in different locations on f44 vs f43 plot and this finding can be used as a tool for source apportionment of measured particles.

Computational micromodel for epigenetic mechanisms

Characterization of the epigenetic profile of humans since the initial breakthrough on the human genome project has strongly established the key role of histone modifications and DNA methylation. These dynamic elements interact to determine the normal level of expression or methylation status of the constituent genes in the genome. Recently, considerable evidence has been put forward to demonstrate that environmental stress implicitly alters epigenetic patterns causing imbalance that can lead to cancer initiation. This chain of consequences has motivated attempts to computationally model the influence of histone modification and DNA methylation in gene expression and investigate their intrinsic interdependency. In this paper, we explore the relation between DNA methylation and transcription and characterize in detail the histone modifications for specific DNA methylation levels using a stochastic approach.

Computational analysis of epigenetic information in human DNA sequences

Over the last few years, investigations of human epigenetic profiles have identified key elements of change to be Histone Modifications, stable and heritable DNA methylation and Chromatin remodeling. These factors determine gene expression levels and characterise conditions leading to disease. In order to extract information embedded in long DNA sequences, data mining and pattern recognition tools are widely used, but efforts have been limited to date with respect to analyzing epigenetic changes, and their role as catalysts in disease onset. Useful insight, however, can be gained by investigation of associated dinucleotide distributions. The focus of this paper is to explore specific dinucleotides frequencies across defined regions within the human genome, and to identify new patterns between epigenetic mechanisms and DNA content. Signal processing methods, including Fourier and Wavelet Transformations, are employed and principal results are reported.

Whole-brain imaging with single-cell resolution using chemical cocktails and computational analysis

Systems-level identification and analysis of cellular circuits in the brain will require the development of whole-brain imaging with single-cell resolution. To this end, we performed comprehensive chemical screening to develop a whole-brain clearing and imaging method, termed CUBIC (clear, unobstructed brain imaging cocktails and computational analysis). CUBIC is a simple and efficient method involving the immersion of brain samples in chemical mixtures containing aminoalcohols, which enables rapid whole-brain imaging with single-photon excitation microscopy. CUBIC is applicable to multicolor imaging of fluorescent proteins or immunostained samples in adult brains and is scalable from a primate brain to subcellular structures. We also developed a whole-brain cell-nuclear counterstaining protocol and a computational image analysis pipeline that, together with CUBIC reagents, enable the visualization and quantification of neural activities induced by environmental stimulation. CUBIC enables time-course expression profiling of whole adult brains with single-cell resolution.

An improved computational method in structural dynamics

Purpose – In structural, earthquake and aeronautical engineering and mechanical vibration, the solution of dynamic equations for a structure subjected to dynamic loading leads to a high order system of differential equations. The numerical methods are usually used for integration when either there is dealing with discrete data or there is no analytical solution for the equations. Since the numerical methods with more accuracy and stability give more accurate results in structural responses, there is a need to improve the existing methods or develop new ones. The paper aims to discuss these issues. Design/methodology/approach – In this paper, a new time integration method is proposed mathematically and numerically, which is accordingly applied to single-degree-of-freedom (SDOF) and multi-degree-of-freedom (MDOF) systems. Finally, the results are compared to the existing methods such as Newmark’s method and closed form solution. Findings – It is concluded that, in the proposed method, the data variance of each set of structural responses such as displacement, velocity, or acceleration in different time steps is less than those in Newmark’s method, and the proposed method is more accurate and stable than Newmark’s method and is capable of analyzing the structure at fewer numbers of iteration or computation cycles, hence less time-consuming. Originality/value – A new mathematical and numerical time integration method is proposed for the computation of structural responses with higher accuracy and stability, lower data variance, and fewer numbers of iterations for computational cycles.

Biophysical correlates of intrinsic and stress-induced morphological variability in lateral amygdaloid neurons: A computational study

Morphological and physiological characteristics of neurons located in the dorsolateral and two ventral subdivisions of the lateral amygdala (LA) have been compared in order to differentiate their roles in the formation and storage of fear memories (Alphs et al, SfN abs 623.1, 2003). Briefly, in these populations, significant differences are observed in input resistance, membrane time constant, firing frequency, dendritic tortuosity, numbers of primary dendrites, dendritic segments and dendritic nodes...

Numerical investigation of motion and deformation of a single red blood cell in a stenosed capillary

It is generally assumed that influence of the red blood cells (RBCs) is predominant in blood rheology. The healthy RBCs are highly deformable and can thus easily squeeze through the smallest capillaries having internal diameter less than their characteristic size. On the other hand, RBCs infected by malaria or other diseases are stiffer and so less deformable. Thus it is harder for them to flow through the smallest capillaries. Therefore, it is very important to critically and realistically investigate the mechanical behavior of both healthy and infected RBCs which is a current gap in knowledge. The motion and the steady state deformed shape of the RBCs depend on many factors, such as the geometrical parameters of the capillary through which blood flows, the membrane bending stiffness and the mean velocity of the blood flow. In this study, motion and deformation of a single two-dimensional RBC in a stenosed capillary is explored by using smoothed particle hydrodynamics (SPH) method. An elastic spring network is used to model the RBC membrane, while the RBC's inside fluid and outside fluid are treated as SPH particles. The effect of RBC's membrane stiffness (kb), inlet pressure (P) and geometrical parameters of the capillary on the motion and deformation of the RBC is studied. The deformation index, RBC's mean velocity and the cell membrane energy are analyzed when the cell passes through the stenosed capillary. The simulation results demonstrate that the kb, P and the geometrical parameters of the capillary have a significant impact on the RBCs' motion and deformation in the stenosed section.

New insights into the molecular structure of kaolinite–methanol intercalation complexes

A series of kaolinite–methanol complexes with different basal spacings were synthesized using guest displacement reactions of the intercalation precursors kaolinite–N-methyformamide (Kaol–NMF), kaolinite–urea (Kaol–U), or kaolinite–dimethylsulfoxide (Kaol–DMSO), with methanol (Me). The interaction of methanol with kaolinite was examined using X-ray diffraction (XRD), infrared spectroscopy (IR), and nuclear magnetic resonance (NMR). Kaolinite (Kaol) initially intercalated with N-methyformamide (NMF), urea (U), or dimethylsulfoxide (DMSO) before subsequent reaction with Me formed final kaolinite–methanol (Kaol–Me) complexes characterized by basal spacing ranging between 8.6 Å and 9.6 Å, depending on the pre-intercalated reagent. Based on a comparative analysis of the three Kaol–Me displacement intercalation complexes, three types of Me intercalation products were suggested to have been present in the interlayer space of Kaol: (1) molecules grafted onto a kaolinite octahedral sheet in the form of a methoxy group (Al-O-C bond); (2) mobile Me and/or water molecules kept in the interlayer space via hydrogen bonds that could be partially removed during drying; and (3) a mixture of types 1 and 2, with the methoxy group (Al-O-C bond) grafted onto the Kaol sheet and mobile Me and/or water molecules coexisted in the system after the displacement reaction by Me. Various structural models that reflected four possible complexes of Kaol–Me were constructed for use in a complimentary computational study. Results from the calculation of the methanol kaolinite interaction indicate that the hydroxyl oxygen atom of methanol plays the dominant role in the stabilization and localization of the molecule intercalated in the interlayer space, and that water existing in the intercalated Kaol layer is inevitable.

Bayesian computational methods for spatial analysis of images

This thesis introduces a new way of using prior information in a spatial model and develops scalable algorithms for fitting this model to large imaging datasets. These methods are employed for image-guided radiation therapy and satellite based classification of land use and water quality. This study has utilized a pre-computation step to achieve a hundredfold improvement in the elapsed runtime for model fitting. This makes it much more feasible to apply these models to real-world problems, and enables full Bayesian inference for images with a million or more pixels.

Modelling CO2 adsorption and separation on experimentally-realized B40 fullerene

Searching for efficient solid sorbents for CO2 adsorption and separation is important for developing emergent carbon reduction and natural gas purification technology. This work, for the first time, has investigated the adsorption of CO2 on newly experimentally realized cage-like B40 fullerene (Zhai et al., 2014) based on density functional theory calculations. We find that the adsorption of CO2 on B40 fullerene involves a relatively large energy barrier (1.21 eV), however this can be greatly decreased to 0.35 eV by introducing an extra electron. A practical way to realize negatively charged B40 fullerene is then proposed by encapsulating a Li atom into the B40 fullerene (Li@B40). Li@B40 is found to be highly stable and can significantly enhance both the thermodynamics and kinetics of CO2 adsorption, while the adsorptions of N2, CH4 and H2 on the Li@B40 fullerene remain weak in comparison. Since B40 fullerene has been successfully synthesized in a most recent experiment, our results highlight a new promising material for CO2 capture and separation for future experimental validation.

Stochastic analysis of a radial-inflow turbine in the presence of parametric uncertainties

This paper presents an uncertainty quantification study of the performance analysis of the high pressure ratio single stage radial-inflow turbine used in the Sundstrand Power Systems T-100 Multi-purpose Small Power Unit. A deterministic 3D volume-averaged Computational Fluid Dynamics (CFD) solver is coupled with a non-statistical generalized Polynomial Chaos (gPC) representation based on a pseudo-spectral projection method. One of the advantages of this approach is that it does not require any modification of the CFD code for the propagation of random disturbances in the aerodynamic and geometric fields. The stochastic results highlight the importance of the blade thickness and trailing edge tip radius on the total-to-static efficiency of the turbine compared to the angular velocity and trailing edge tip length. From a theoretical point of view, the use of the gPC representation on an arbitrary grid also allows the investigation of the sensitivity of the blade thickness profiles on the turbine efficiency. The gPC approach is also applied to coupled random parameters. The results show that the most influential coupled random variables are trailing edge tip radius coupled with the angular velocity.

Sodium and lithium storage properties of spray-dried molybdenum disulfide-graphene hierarchical microspheres

Developing nano/micro-structures which can effectively upgrade the intriguing properties of electrode materials for energy storage devices is always a key research topic. Ultrathin nanosheets were proved to be one of the potential nanostructures due to their high specific surface area, good active contact areas and porous channels. Herein, we report a unique hierarchical micro-spherical morphology of well-stacked and completely miscible molybdenum disulfide (MoS2) nanosheets and graphene sheets, were successfully synthesized via a simple and industrial scale spray-drying technique to take the advantages of both MoS2 and graphene in terms of their high practical capacity values and high electronic conductivity, respectively. Computational studies were performed to understand the interfacial behaviour of MoS2 and graphene, which proves high stability of the composite with high interfacial binding energy (−2.02 eV) among them. Further, the lithium and sodium storage properties have been tested and reveal excellent cyclic stability over 250 and 500 cycles, respectively, with the highest initial capacity values of 1300 mAh g−1 and 640 mAh g−1 at 0.1 A g−1.

Influence of pipe length and flow rate on nano-particle deposition in laminar circular pipe flows

The Lagrangian particle tracking provides an effective method for simulating the deposition of nano- particles as well as micro-particles as it accounts for the particle inertia effect as well as the Brownian excitation. However, using the Lagrangian approach for simulating ultrafine particles has been limited due to computational cost and numerical difficulties. The aim of this paper is to study the deposition of nano-particles in cylindrical tubes under laminar condition using the Lagrangian particle tracking method. The commercial Fluent software is used to simulate the fluid flow in the pipes and to study the deposition and dispersion of nano-particles. Different particle diameters as well as different pipe lengths and flow rates are examined. The results show good agreement between the calculated deposition efficiency and different analytic correlations in the literature. Furthermore, for the nano-particles with higher diameters and when the effect of inertia has a higher importance, the calculated deposition efficiency by the Lagrangian method is less than the analytic correlations based on Eulerian method due to statistical error or the inertia effect.