6,607 research outputs found
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 -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 colors. Of both theoretical and practical significance, -PATS
has been studied in a series of papers which have shown -PATS to be NP-hard
for , , and then . 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 SrLaFeCoO (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
SrLaFeCoO (x = 0, 1, 2) of which first (SrFeCoO with x
= 0) has a canonical spin spin glass, the second (SrLaFeCoO with x = 1) has
a so-called magnetic glass and the third (LaFeCoO with x = 2) has a
magnetically ordered ground state. Our present study revealed clear inelastic
signal for SrLaFeCoO, possibly also inelastic signal for SrFeCoO
below the spin freezing temperatures but no inelastic signal at all
for for the magnetically ordered LaFeCoO 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
- …