Evaluation of local and global atrophy measurement techniques with simulated Alzheimer's disease data

The main goal of this work was to evaluate several well-known methods which provide global (BSI and SIENA) or local (Jacobian integration) estimates of atrophy in brain structures using Magnetic Resonance images. For that purpose, we have generated realistic simulated Alzheimer's disease images in which volume changes are modelled with a Finite Element thermoelastic model, which mimic the patterns of change obtained from a cohort of 19 real controls and 27 probable Alzheimer's disease patients. SIENA and BSI results correlate very well with gold standard data (BSI mean absolute error <0.29%; SIENA <0.44%). Jacobian integration was guided by both fluid and FFD-based registration techniques and resulting deformation fields and associated Jacobians were compared, region by region, with gold standard ones. The FFD registration technique provided more satisfactory results than the fluid one. Mean absolute error differences between volume changes given by the FFD-based technique and the gold standard were: sulcal CSF <2.49%; lateral ventricles 2.25%; brain <0.36%; hippocampi <0.42%

Longitudinal Voxel-based morphometry with unified segmentation: evaluation on simulated Alzheimer’s disease

The goal of this work is to evaluate Voxel-Based Morphometry and three longitudinally-tailored methods of VBM.We use a cohort of simulated images produced by deforming original scans using a Finite Element Method, guided to emulate Alzheimer-like changes. The simulated images provide quite realistic data with a known pattern of spatial atrophy, with which VBM’s findings can be meaningfully compared. We believe this is the first evaluation of VBM for which anatomically-plausible ‘gold-standard’ results are available. The three longitudinal VBM methods have been implemented within the unified segmentation framework of SPM5; one of the techniques is a newly developed procedure, which shows promising potential

Phenomenological model of diffuse global and regional atrophy using finite-element methods

The main goal of this work is the generation of ground-truth data for the validation of atrophy measurement techniques, commonly used in the study of neurodegenerative diseases such as dementia. Several techniques have been used to measure atrophy in cross-sectional and longitudinal studies, but it is extremely difficult to compare their performance since they have been applied to different patient populations. Furthermore, assessment of performance based on phantom measurements or simple scaled images overestimates these techniques' ability to capture the complexity of neurodegeneration of the human brain. We propose a method for atrophy simulation in structural magnetic resonance (MR) images based on finite-element methods. The method produces cohorts of brain images with known change that is physically and clinically plausible, providing data for objective evaluation of atrophy measurement techniques. Atrophy is simulated in different tissue compartments or in different neuroanatomical structures with a phenomenological model. This model of diffuse global and regional atrophy is based on volumetric measurements such as the brain or the hippocampus, from patients with known disease and guided by clinical knowledge of the relative pathological involvement of regions and tissues. The consequent biomechanical readjustment of structures is modelled using conventional physics-based techniques based on biomechanical tissue properties and simulating plausible tissue deformations with finite-element methods. A thermoelastic model of tissue deformation is employed, controlling the rate of progression of atrophy by means of a set of thermal coefficients, each one corresponding to a different type of tissue. Tissue characterization is performed by means of the meshing of a labelled brain atlas, creating a reference volumetric mesh that will be introduced to a finite-element solver to create the simulated deformations. Preliminary work on the simulation of acquisition artefa- - cts is also presented. Cross-sectional and

High Pressure Thermoelasticity of Body-centered Cubic Tantalum

We have investigated the thermoelasticity of body-centered cubic (bcc) tantalum from first principles by using the linearized augmented plane wave (LAPW) and mixed--basis pseudopotential methods for pressures up to 400 GPa and temperatures up to 10000 K. Electronic excitation contributions to the free energy were included from the band structures, and phonon contributions were included using the particle-in-a-cell (PIC) model. The computed elastic constants agree well with available ultrasonic and diamond anvil cell data at low pressures, and shock data at high pressures. The shear modulus $c_{44}$ and the anisotropy change behavior with increasing pressure around 150 GPa because of an electronic topological transition. We find that the main contribution of temperature to the elastic constants is from the thermal expansivity. The PIC model in conjunction with fast self-consistent techniques is shown to be a tractable approach to studying thermoelasticity.Comment: To be appear in Physical Review

The size of the proton - closing in on the radius puzzle

We analyze the recent electron-proton scattering data from Mainz using a dispersive framework that respects the constraints from analyticity and unitarity on the nucleon structure. We also perform a continued fraction analysis of these data. We find a small electric proton charge radius, r_E^p = 0.84_{-0.01}^{+0.01} fm, consistent with the recent determination from muonic hydrogen measurements and earlier dispersive analyses. We also extract the proton magnetic radius, r_M^p = 0.86_{-0.03}^{+0.02} fm, consistent with earlier determinations based on dispersion relations.Comment: 4 pages, 2 figures, fit improved, small modifications, section on continued fractions modified, conclusions on the proton charge radius unchanged, version accepted for publication in European Physical Journal

A mean-field kinetic lattice gas model of electrochemical cells

We develop Electrochemical Mean-Field Kinetic Equations (EMFKE) to simulate electrochemical cells. We start from a microscopic lattice-gas model with charged particles, and build mean-field kinetic equations following the lines of earlier work for neutral particles. We include the Poisson equation to account for the influence of the electric field on ion migration, and oxido-reduction processes on the electrode surfaces to allow for growth and dissolution. We confirm the viability of our approach by simulating (i) the electrochemical equilibrium at flat electrodes, which displays the correct charged double-layer, (ii) the growth kinetics of one-dimensional electrochemical cells during growth and dissolution, and (iii) electrochemical dendrites in two dimensions.Comment: 14 pages twocolumn, 17 figure

Calibration of Super-Kamiokande Using an Electron Linac

In order to calibrate the Super-Kamiokande experiment for solar neutrino measurements, a linear accelerator (LINAC) for electrons was installed at the detector. LINAC data were taken at various positions in the detector volume, tracking the detector response in the variables relevant to solar neutrino analysis. In particular, the absolute energy scale is now known with less than 1 percent uncertainty.Comment: 24 pages, 16 figures, Submitted to NIM

Search for Neutral Q-balls in Super-Kamiokande II

A search for Q-balls induced groups of successive contained events has been carried out in Super-Kamiokande II with 541.7 days of live time. Neutral Q-balls would emit pions when colliding with nuclei, generating a signal of successive contained pion events along a track. No candidate for successive contained event groups has been found in Super-Kamiokande II, so upper limits on the possible flux of such Q-balls have been obtained.Comment: 5 pages, 5 figures, Submitted to Phys. Lett.

Hypogene Calcitization: Evaporite Diagenesis in the Western Delaware Basin

Evaporite calcitization within the Castile Formation of the Delaware Basin is more widespread and diverse than originally recognized. Coupled field and GIS studies have identified more than 1000 individual occurrences of calcitization within the Castile Formation outcrop area, which includes both calcitized masses (limestone buttes) and laterally extensive calcitized horizons (limestone sheets). Both limestone buttes and sheets commonly contain a central brecciated zone that we attribute to hypogene dissolution. Lithologic fabric of calcitized zones ranges from little alteration of original varved laminae to fabrics showing extensive laminae distortion as well as extensive vuggy and open cavernous porosity. Calcitization is most abundant in the western portion of the Castile outcrop region where surface denudation has been greatest. Calcitization often forms linear trends, indicating fluid migration along fractures, but also occurs as dense clusters indicating focused, ascending, hydrocarbon-rich fluids. Native sulfur, secondary tabular gypsum (i.e. selenite) and hypogene caves are commonly associated with clusters of calcitization. This assemblage suggests that calcium sulfate diagenesis within the Castile Formation is dominated by hypogene speleogemesis

On the selection of AGN neutrino source candidates for a source stacking analysis with neutrino telescopes

The sensitivity of a search for sources of TeV neutrinos can be improved by grouping potential sources together into generic classes in a procedure that is known as source stacking. In this paper, we define catalogs of Active Galactic Nuclei (AGN) and use them to perform a source stacking analysis. The grouping of AGN into classes is done in two steps: first, AGN classes are defined, then, sources to be stacked are selected assuming that a potential neutrino flux is linearly correlated with the photon luminosity in a certain energy band (radio, IR, optical, keV, GeV, TeV). Lacking any secure detailed knowledge on neutrino production in AGN, this correlation is motivated by hadronic AGN models, as briefly reviewed in this paper. The source stacking search for neutrinos from generic AGN classes is illustrated using the data collected by the AMANDA-II high energy neutrino detector during the year 2000. No significant excess for any of the suggested groups was found.Comment: 43 pages, 12 figures, accepted by Astroparticle Physic