Evidence cross-validation and Bayesian inference of MAST plasma equilibria

In this paper, current profiles for plasma discharges on the Mega-Ampere Spherical Tokamak (MAST) are directly calculated from pickup coil, flux loop and Motional-Stark Effect (MSE) observations via methods based in the statistical theory of Bayesian analysis. By representing toroidal plasma current as a series of axisymmetric current beams with rectangular cross-section and inferring the current for each one of these beams, flux-surface geometry and q-profiles are subsequently calculated by elementary application of Biot-Savart's law. The use of this plasma model in the context of Bayesian analysis was pioneered by Svensson and Werner on the Joint-European Tokamak (JET) [J. Svensson and A. Werner. Current tomography for axisymmetric plasmas. Plasma Physics and Controlled Fusion, 50(8):085002, 2008]. In this framework, linear forward models are used to generate diagnostic predictions, and the probability distribution for the currents in the collection of plasma beams was subsequently calculated directly via application of Bayes' formula. In this work, we introduce a new diagnostic technique to identify and remove outlier observations associated with diagnostics falling out of calibration or suffering from an unidentified malfunction. These modifications enable good agreement between Bayesian inference of the last closed flux-surface (LCFS) with other corroborating data, such as such as that from force balance considerations using EFIT++ [L. Appel et al., Proc. 33rd EPS Conf., Rome, Italy, 2006]. In addition, this analysis also yields errors on the plasma current profile and flux-surface geometry, as well as directly predicting the Shafranov shift of the plasma core.This work was jointly funded by the Australian Government through International Science Linkages Grant No. CG130047, the Australian National University, the United Kingdom Engineering and Physical Sciences Research Council under Grant No. EP/G003955, and by the European Communities under the contract of Association between EURATOM and CCFE

Fermi-Polaron: Diagrammatic Monte Carlo for Divergent Sign-Alternating Series

Diagrammatic Monte Carlo approach is applied to a problem of a single spin-down fermion resonantly interacting with the sea of ideal spin-up fermions. On one hand, we develop a generic, sign-problem tolerant, method of exact numerical solution of polaron-type models. On the other hand, our solution is important for understanding the phase diagram and properties of the BCS-BEC crossover in the strongly imbalanced regime. This is the first, and possibly characteristic, example of how the Monte Carlo approach can be applied to a divergent sign-alternating diagrammatic series.Comment: 4 pages, 7 figure

Formalizing Size-Optimal Sorting Networks: Extracting a Certified Proof Checker

Since the proof of the four color theorem in 1976, computer-generated proofs have become a reality in mathematics and computer science. During the last decade, we have seen formal proofs using verified proof assistants being used to verify the validity of such proofs. In this paper, we describe a formalized theory of size-optimal sorting networks. From this formalization we extract a certified checker that successfully verifies computer-generated proofs of optimality on up to 8 inputs. The checker relies on an untrusted oracle to shortcut the search for witnesses on more than 1.6 million NP-complete subproblems.Comment: IMADA-preprint-c

Photoemission spectroscopy and sum rules in dilute electron-phonon systems

A family of exact sum rules for the one-polaron spectral function in the low-density limit is derived. An algorithm to calculate energy moments of arbitrary order of the spectral function is presented. Explicit expressions are given for the first two moments of a model with general electron-phonon interaction, and for the first four moments of the Holstein polaron. The sum rules are linked to experiments on momentum-resolved photoemission spectroscopy. The bare electronic dispersion and the electron-phonon coupling constant can be extracted from the first and second moments of spectrum. The sum rules could serve as constraints in analytical and numerical studies of electron-phonon models.Comment: 4 page

Model Data Fusion: developing Bayesian inversion to constrain equilibrium and mode structure

Recently, a new probabilistic "data fusion" framework based on Bayesian principles has been developed on JET and W7-AS. The Bayesian analysis framework folds in uncertainties and inter-dependencies in the diagnostic data and signal forward-models, together with prior knowledge of the state of the plasma, to yield predictions of internal magnetic structure. A feature of the framework, known as MINERVA (J. Svensson, A. Werner, Plasma Physics and Controlled Fusion 50, 085022, 2008), is the inference of magnetic flux surfaces without the use of a force balance model. We discuss results from a new project to develop Bayesian inversion tools that aim to (1) distinguish between competing equilibrium theories, which capture different physics, using the MAST spherical tokamak; and (2) test the predictions of MHD theory, particularly mode structure, using the H-1 Heliac.Comment: submitted to Journal of Plasma Fusion Research 10/11/200

An implementation of Deflate in Coq

The widely-used compression format "Deflate" is defined in RFC 1951 and is based on prefix-free codings and backreferences. There are unclear points about the way these codings are specified, and several sources for confusion in the standard. We tried to fix this problem by giving a rigorous mathematical specification, which we formalized in Coq. We produced a verified implementation in Coq which achieves competitive performance on inputs of several megabytes. In this paper we present the several parts of our implementation: a fully verified implementation of canonical prefix-free codings, which can be used in other compression formats as well, and an elegant formalism for specifying sophisticated formats, which we used to implement both a compression and decompression algorithm in Coq which we formally prove inverse to each other -- the first time this has been achieved to our knowledge. The compatibility to other Deflate implementations can be shown empirically. We furthermore discuss some of the difficulties, specifically regarding memory and runtime requirements, and our approaches to overcome them

Binary pattern tile set synthesis is NP-hard

In the field of algorithmic self-assembly, a long-standing unproven conjecture has been that of the NP-hardness of binary pattern tile set synthesis (2-PATS). The $k$-PATS problem is that of designing a tile assembly system with the smallest number of tile types which will self-assemble an input pattern of $k$ colors. Of both theoretical and practical significance, $k$-PATS has been studied in a series of papers which have shown $k$-PATS to be NP-hard for $k = 60$, $k = 29$, and then $k = 11$. In this paper, we close the fundamental conjecture that 2-PATS is NP-hard, concluding this line of study. While most of our proof relies on standard mathematical proof techniques, one crucial lemma makes use of a computer-assisted proof, which is a relatively novel but increasingly utilized paradigm for deriving proofs for complex mathematical problems. This tool is especially powerful for attacking combinatorial problems, as exemplified by the proof of the four color theorem by Appel and Haken (simplified later by Robertson, Sanders, Seymour, and Thomas) or the recent important advance on the Erd\H{o}s discrepancy problem by Konev and Lisitsa using computer programs. We utilize a massively parallel algorithm and thus turn an otherwise intractable portion of our proof into a program which requires approximately a year of computation time, bringing the use of computer-assisted proofs to a new scale. We fully detail the algorithm employed by our code, and make the code freely available online

Hyperfine interaction and electronic spin fluctuation study on Sr2−x_{2-x}Lax_xFeCoO6_6 (x = 0, 1, 2) by high-resolution back-scattering neutron spectroscopy

The study of hyperfine interaction by high-resolution inelastic neutron scattering is not very well known compared to the other competing techniques viz. NMR, M\"ossbauer, PACS etc. Also the study is limited mostly to magnetically ordered systems. Here we report such study on Sr$_{2-x}$La$_x$FeCoO$_6$ (x = 0, 1, 2) of which first (Sr$_2$FeCoO$_6$ with x = 0) has a canonical spin spin glass, the second (SrLaFeCoO$_6$ with x = 1) has a so-called magnetic glass and the third (La$_2$FeCoO$_6$ with x = 2) has a magnetically ordered ground state. Our present study revealed clear inelastic signal for SrLaFeCoO$_6$, possibly also inelastic signal for Sr$_2$FeCoO$_6$ below the spin freezing temperatures $T_{sf}$ but no inelastic signal at all for for the magnetically ordered La$_2$FeCoO$_6$ in the neutron scattering spectra. The broadened inelastic signals observed suggest hyperfine field distribution in the two disordered magnetic glassy systems and no signal for the third compound suggests no or very small hyperfine field at the Co nucleus due to Co electronic moment. For the two magnetic glassy system apart from the hyperfine signal due only to Co, we also observed electronic spin fluctuations probably from both Fe and Co electronic moments. \end{abstract

Single Event Effects in the Pixel readout chip for BTeV

In future experiments the readout electronics for pixel detectors is required to be resistant to a very high radiation level. In this paper we report on irradiation tests performed on several preFPIX2 prototype pixel readout chips for the BTeV experiment exposed to a 200 MeV proton beam. The prototype chips have been implemented in commercial 0.25 um CMOS processes following radiation tolerant design rules. The results show that this ASIC design tolerates a large total radiation dose, and that radiation induced Single Event Effects occur at a manageable level.Comment: 15 pages, 6 Postscript figure