Algorithmic procedures for Bayesian MEG/EEG source reconstruction in SPM

AbstractThe MEG/EEG inverse problem is ill-posed, giving different source reconstructions depending on the initial assumption sets. Parametric Empirical Bayes allows one to implement most popular MEG/EEG inversion schemes (Minimum Norm, LORETA, etc.) within the same generic Bayesian framework. It also provides a cost-function in terms of the variational Free energy—an approximation to the marginal likelihood or evidence of the solution. In this manuscript, we revisit the algorithm for MEG/EEG source reconstruction with a view to providing a didactic and practical guide. The aim is to promote and help standardise the development and consolidation of other schemes within the same framework. We describe the implementation in the Statistical Parametric Mapping (SPM) software package, carefully explaining each of its stages with the help of a simple simulated data example. We focus on the Multiple Sparse Priors (MSP) model, which we compare with the well-known Minimum Norm and LORETA models, using the negative variational Free energy for model comparison. The manuscript is accompanied by Matlab scripts to allow the reader to test and explore the underlying algorithm

Sources of Differences in On-Orbit Total Solar Irradiance Measurements and Description of Proposed Laboratory Intercomparison

There is a 5 W/sq m (about 0.35 %) difference between current on-orbit Total Solar Irradiance (TSI) measurements. On 18-20 July 2005, a workshop was held at the National Institute of Standards and Technology (NIST) in Gaithersburg, Maryland that focused on understanding possible reasons for this difference, through an examination of the instrument designs, calibration approaches, and appropriate measurement equations. The instruments studied in that workshop included the Active Cavity Radiometer Irradiance Monitor III (ACRIM III) on the Active Cavity Radiometer Irradiance Monitor SATellite (ACRIMSAT), the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE), the Variability of solar IRradiance and Gravity Oscillations (VIRGO) on the Solar and Heliospheric Observatory (SOHO), and the Earth Radiation Budget Experiment (ERBE) on the Earth Radiation Budget Satellite (ERBS). Presentations for each instrument included descriptions of its design, its measurement equation and uncertainty budget, and the methods used to assess on-orbit degradation. The workshop also included a session on satellite- and ground-based instrument comparisons and a session on laboratory-based comparisons and the application of new laboratory comparison techniques. The workshop has led to investigations of the effects of diffraction and of aperture area measurements on the differences between instruments. In addition, a laboratory-based instrument comparison is proposed that uses optical power measurements (with lasers that underEll the apertures of the TSI instruments), irradiance measurements (with lasers that overfill the apertures of the TSI instrument), and a cryogenic electrical substitution radiometer as a standard for comparing the instruments. A summary of the workshop and an overview of the proposed research efforts are presented here

Non-invasive laminar inference with MEG: comparison of methods and source inversion algorithms

Magnetoencephalography (MEG) is a direct measure of neuronal current flow; its anatomical resolution is therefore not constrained by physiology but rather by data quality and the models used to explain these data. Recent simulation work has shown that it is possible to distinguish between signals arising in the deep and superficial cortical laminae given accurate knowledge of these surfaces with respect to the MEG sensors. This previous work has focused around a single inversion scheme (multiple sparse priors) and a single global parametric fit metric (free energy). In this paper we use several different source inversion algorithms and both local and global, as well as parametric and non-parametric fit metrics in order to demonstrate the robustness of the discrimination between layers. We find that only algorithms with some sparsity constraint can successfully be used to make laminar discrimination. Importantly, local t-statistics, global cross-validation and free energy all provide robust and mutually corroborating metrics of fit. We show that discrimination accuracy is affected by patch size estimates, cortical surface features, and lead field strength, which suggests several possible future improvements to this technique. This study demonstrates the possibility of determining the laminar origin of MEG sensor activity, and thus directly testing theories of human cognition that involve laminar- and frequency-specific mechanisms. This possibility can now be achieved using recent developments in high precision MEG, most notably the use of subject-specific head-casts, which allow for significant increases in data quality and therefore anatomically precise MEG recordings

Auxiliary particle theory of threshold singularities in photoemission and X-ray absorption spectra: Test of a conserving T-matrix approximation

We calculate the exponents of the threshold singularities in the photoemission spectrum of a deep core hole and its X-ray absorption spectrum in the framework of a systematic many-body theory of slave bosons and pseudofermions (for the empty and occupied core level). In this representation, photoemission and X-ray absorption can be understood on the same footing; no distinction between orthogonality catastrophe and excitonic effects is necessary. We apply the conserving slave particle T-matrix approximation (CTMA), recently developed to describe both Fermi and non-Fermi liquid behavior systems with strong local correlations, to the X-ray problem as a test case. The numerical results for both photoemission and X-ray absorption are found to be in agreement with the exact infrared powerlaw behavior in the weak as well as in the strong coupling regions. We point out a close relation of the CTMA with the parquet equation approach of Nozi{\`e}res et al.Comment: 10 pages, 9 figures, published versio

Vortex states in superconducting rings

The superconducting state of a thin superconducting disk with a hole is studied within the non-linear Ginzburg-Landau theory in which the demagnetization effect is accurately taken into account. We find that the flux through the hole is not quantized, the superconducting state is stabilized with increasing size of the hole for fixed radius of the disk, and a transition to a multi-vortex state is found if the disk is sufficiently large. Breaking the circular summetry through a non central location of the hole in the disk enhances the multi-vortex state.Comment: 11 pages, 23 figures (postscript). To appear in Physical Review B, Vol. 61 (2000

Dark mammoth trunks in the merging galaxy NGC 1316 and a mechanism of cosmic double helices

NGC 1316 is a giant, elliptical galaxy containing a complex network of dark, dust features. The morphology of these features has been examined in some detail using a Hubble Space Telescope, Advanced Camera for Surveys image. It is found that most of the features are constituted of long filaments. There also exist a great number of dark structures protruding inwards from the filaments. Many of these structures are strikingly similar to elephant trunks in H II regions in the Milky Way Galaxy, although much larger. The structures, termed mammoth trunks, generally are filamentary and often have shapes resembling the letters V or Y. In some of the mammoth trunks the stem of the Y can be resolved into two or more filaments, many of which showing signs of being intertwined. A model of the mammoth trunks, related to a recent theory of elephant trunks, is proposed. Based on magnetized filaments, the model is capable of giving an account of the various shapes of the mammoth trunks observed, including the twined structures.Comment: Accepted for publication in Astrophysics & Space Scienc

Confirmation of a pi_1^0 Exotic Meson in the \eta \pi^0 System

The exclusive reaction $\pi^- p \to \eta \pi^0 n$, $\eta \to \pi^+ \pi^- \pi^0$ at 18 GeV$/c$ has been studied with a partial wave analysis on a sample of 23~492 $\eta \pi^0 n$ events from BNL experiment E852. A mass-dependent fit is consistent with a resonant hypothesis for the $P_+$ wave, thus providing evidence for a neutral exotic meson with $J^{PC} = 1^{-+}$, a mass of $1257 \pm 20 \pm 25$ MeV$/c^2$, and a width of $354 \pm 64 \pm 60$ MeV$/c^2$. New interpretations of the meson exotics in neutral $\eta \pi^0$ system observed in E852 and Crystal Barrel experiments are discussed.Comment: p3, rewording the paragraph (at the bottom) about the phase variations. p4, rewording paragrath "The second method ..." . p4, at the bottom of paragrath "The third method ..." added consistent with the results of methods 1 and 2

The Nucleon-Nucleon Interaction in the Chromo-Dielectric Soliton Model: Dynamics

The present work is an extension of a previous study of the nucleon-nucleon interaction based on the chromo-dielectric soliton model. The former approach was static, leading to an adiabatic potential. Here we perform a dynamical study in the framework of the Generator Coordinate Method. In practice, we derive an approximate Hill-Wheeler differential equation and obtain a local nucleon-nucleon potential as a function of a mean generator coordinate. This coordinate is related to an effective separation distance between the two nucleons by a Fujiwara transformation. This latter relationship is especially useful in studying the quark substructure of light nuclei. We investigate the explicit contribution of the one-gluon exchange part of the six-quark Hamiltonian to the nucleon-nucleon potential, and we find that the dynamics are responsible for a significant part of the short-range N-N repulsion.Comment: 16 pages (REVTEX), 6 figures (uuencoded Postscript) optionally included using epsfig.st

Towers and fibered products of model categories

Given a left Quillen presheaf of localized model structures, we study the homotopy limit model structure on the associated category of sections. We focus specifically on towers and fibered products of model categories. As applications we consider Postnikov towers of model categories, chromatic towers of spectra and Bousfield arithmetic squares of spectra. For spectral model categories, we show that the homotopy fiber of a stable left Bousfield localization is a stable right Bousfield localization