-
Teaching Higher-Order Logic Using Isabelle
Authors:
Simon Tobias Lund,
Jørgen Villadsen
Abstract:
We present a formalization of higher-order logic in the Isabelle proof assistant, building directly on the foundational framework Isabelle/Pure and developed to be as small and readable as possible. It should therefore serve as a good introduction for someone looking into learning about higher-order logic and proof assistants, without having to study the much more complex Isabelle/HOL with heavier…
▽ More
We present a formalization of higher-order logic in the Isabelle proof assistant, building directly on the foundational framework Isabelle/Pure and developed to be as small and readable as possible. It should therefore serve as a good introduction for someone looking into learning about higher-order logic and proof assistants, without having to study the much more complex Isabelle/HOL with heavier automation. To showcase our development and approach we explain a sample proof, describe the axioms and rules of our higher-order logic, and discuss our experience with teaching the subject in a classroom setting.
△ Less
Submitted 8 April, 2024;
originally announced April 2024.
-
The Influence of Validation Data on Logical and Scientific Interpretations of Forensic Expert Opinions
Authors:
Steven P. Lund,
Hari Iyer
Abstract:
Forensic experts use specialized training and knowledge to enable other members of the judicial system to make better informed and more just decisions. Factfinders, in particular, are tasked with judging how much weight to give to experts' reports and opinions. Many references describe assessing evidential weight from the perspective of a forensic expert. Some recognize that stakeholders are each…
▽ More
Forensic experts use specialized training and knowledge to enable other members of the judicial system to make better informed and more just decisions. Factfinders, in particular, are tasked with judging how much weight to give to experts' reports and opinions. Many references describe assessing evidential weight from the perspective of a forensic expert. Some recognize that stakeholders are each responsible for evaluating their own weight of evidence. Morris (1971, 1974, 1977) provided a general framework for recipients to update their own uncertainties after learning an expert's opinion. Although this framework is normative under Bayesian axioms and several forensic scholars advocate the use of Bayesian reasoning, few resources describe its application in forensic science. This paper addresses this gap by examining how recipients can combine principles of science and Bayesian reasoning to evaluate their own likelihood ratios for expert opinions. This exercise helps clarify how an expert's role depends on whether one envisions recipients to be logical and scientific or deferential. Illustrative examples with an expert's opinion expressed as a categorical conclusion, likelihood ratio, or range of likelihood ratios, or with likelihood ratios from multiple experts, each reveal the importance and influence of validation data for logical recipients' interpretations.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
The Early Universe as an Open Quantum System: Complexity and Decoherence
Authors:
Arpan Bhattacharyya,
Suddhasattwa Brahma,
S. Shajidul Haque,
Jacob S. Lund,
Arpon Paul
Abstract:
In this work, we extend previous results, demonstrating how complexity in an open quantum system can identify decoherence between two fields, even in the presence of an accelerating background. Using the curved-space Caldeira-Leggett two-field model in de Sitter as our toy model, we discover a distinctive feature in the growth of complexity of purification, providing an alternative diagnostic for…
▽ More
In this work, we extend previous results, demonstrating how complexity in an open quantum system can identify decoherence between two fields, even in the presence of an accelerating background. Using the curved-space Caldeira-Leggett two-field model in de Sitter as our toy model, we discover a distinctive feature in the growth of complexity of purification, providing an alternative diagnostic for studying decoherence when the adiabatic perturbation is coupled to a heavy field. This paper initiates a new pathway to explore the features of quantum complexity in an accelerating background, thereby expanding our understanding of the evolution of primordial cosmological perturbations in the early universe.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
The Galactic Chemical Evolution of phosphorus observed with IGRINS
Authors:
G. Nandakumar,
N. Ryde,
M. Montelius,
B. Thorsbro,
H. Jönsson,
G. Mace,
Lund Observatory,
Department of Astronomy,
Theoretical Physics,
Lund University,
Box 43,
SE-221 00 Lund,
Sweden,
Kapteyn Astronomical Institute,
University of Groningen,
Landleven 12,
NL-9747 AD Groningen,
the Netherlands,
Department of Astronomy,
School of Science,
The University of Tokyo,
7-3-1 Hongo,
Bunkyo-ku,
Tokyo 113-0033,
Japan
, et al. (11 additional authors not shown)
Abstract:
Phosphorus (P) is considered to be one of the key elements for life, making it an important element to look for in the abundance analysis of spectra of stellar systems. Yet, there exists only a handful of spectroscopic studies to estimate the P abundances and investigate its trend across a range of metallicities. We have observed full HK band spectra at a spectral resolving power of R=45,000 with…
▽ More
Phosphorus (P) is considered to be one of the key elements for life, making it an important element to look for in the abundance analysis of spectra of stellar systems. Yet, there exists only a handful of spectroscopic studies to estimate the P abundances and investigate its trend across a range of metallicities. We have observed full HK band spectra at a spectral resolving power of R=45,000 with IGRINS instrument. Abundances are determined using SME in combination with 1D MARCS stellar atmosphere models. The investigated sample of stars have reliable stellar parameters estimated using optical FIES spectra (GILD; Jönsson et al. in prep.). In order to determine the P abundances from the 16482.92 Angstrom P line, we take special care of the CO($ν=7-4$) blend. We determine the C, N, O abundances from atomic carbon and a range of non-blended molecular lines (CO, CN, OH) which are aplenty in the H band region of K giant stars, assuring an appropriate modelling of the blending CO($ν=7-4$) line. We present [P/Fe] vs [Fe/H] trend for 38 K giant stars in the metallicity range of -1.2 dex $<$ [Fe/H] $<$ 0.4 dex. We find that our trend matches well with the compiled literature sample of prominently dwarf stars and limited number of giant stars. Our trend is found to be higher by $\sim$ 0.05 - 0.1 dex compared to the theoretical chemical evolution trend in Cescutti et al. 2012 resulting from core collapse supernova (type II) of massive stars with the P yields from Kobayashi et al. (2006) arbitrarily increased by a factor of 2.75. Thus the enhancement factor might need to be $\sim$ 0.05 - 0.1 dex higher to match our trend. We also find an empirically determined primary behaviour for phosphorus. Furthermore, the phosphorus abundance is found to be elevated by $\sim$ 0.6 - 0.9 dex in two metal poor s-enriched stars compared to the theoretical chemical evolution trend.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Snowmass'21 Accelerator Frontier Report
Authors:
S. Gourlay,
T. Raubenheimer,
V. Shiltsev,
G. Arduini,
R. Assmann,
C. Barbier,
M. Bai,
S. Belomestnykh,
S. Bermudez,
P. Bhat,
A. Faus-Golfe,
J. Galambos,
C. Geddes,
G. Hoffstaetter,
M. Hogan,
Z. Huang,
M. Lamont,
D. Li,
S. Lund,
R. Milner,
P. Musumeci,
E. Nanni,
M. Palmer,
N. Pastrone,
F. Pellemoine
, et al. (13 additional authors not shown)
Abstract:
In 2020-2022, extensive discussions and deliberations have taken place in corresponding topical working groups of the Snowmass Accelerator Frontier (AF) and in numerous joint meetings with other Frontiers, Snowmass-wide meetings, a series of Colloquium-style Agoras, cross-Frontier Forums on muon and electron-positron colliders and the collider Implementation Task Force (ITF). The outcomes of these…
▽ More
In 2020-2022, extensive discussions and deliberations have taken place in corresponding topical working groups of the Snowmass Accelerator Frontier (AF) and in numerous joint meetings with other Frontiers, Snowmass-wide meetings, a series of Colloquium-style Agoras, cross-Frontier Forums on muon and electron-positron colliders and the collider Implementation Task Force (ITF). The outcomes of these activities are summarized in this Accelerator Frontier report.
△ Less
Submitted 17 November, 2022; v1 submitted 28 September, 2022;
originally announced September 2022.
-
Summary Report of AF1 to Snowmass 2021: Beam Physics and Accelerator Education within the Accelerator Frontier
Authors:
M. Bai,
Z. Huang,
S. M. Lund
Abstract:
This report summarizes the findings of the AF1 Topical Subgroup to Snowmass 2021, which investigates beam physics and accelerator education within the accelerator frontiers (AF). The report focuses primarily on opportunities for basic accelerator and beam physics, formulated into four grand challenges. In addition, the report also studies the current status of accelerator education, outreach and d…
▽ More
This report summarizes the findings of the AF1 Topical Subgroup to Snowmass 2021, which investigates beam physics and accelerator education within the accelerator frontiers (AF). The report focuses primarily on opportunities for basic accelerator and beam physics, formulated into four grand challenges. In addition, the report also studies the current status of accelerator education, outreach and diversity issues. Physics limits of ultimate beams for future colliders are also briefly discussed.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Bayesian Reasoning and Evidence Communication
Authors:
Steven Lund,
Hari Iyer
Abstract:
Many resources for forensic scholars and practitioners, such as journal articles, guidance documents, and textbooks, address how to make a value of evidence assessment in the form of a likelihood ratio (LR) when deciding between two competing propositions. These texts often describe experts presenting their LR values to other parties in the judicial system, such as lawyers, judges, and potentially…
▽ More
Many resources for forensic scholars and practitioners, such as journal articles, guidance documents, and textbooks, address how to make a value of evidence assessment in the form of a likelihood ratio (LR) when deciding between two competing propositions. These texts often describe experts presenting their LR values to other parties in the judicial system, such as lawyers, judges, and potentially jurors, but few texts explicitly address how a recipient is expected to utilize the provided LR value. Those that do often imply, or directly suggest, a hybrid modification of Bayes' rule in which a decision maker multiplies their prior odds with another person's assessment of LR to obtain their posterior odds. In this paper, we illustrate how someone adhering to Bayesian reasoning would update their personal uncertainty in response to someone else presenting a personal LR value (or any other form of an opinion) and emphasize that the hybrid approach is a departure from Bayesian reasoning. We further consider implications of recipients adhering to Bayesian reasoning on the role and ideal content of expert's reports and testimony and address published responses to our 2017 paper (Lund and Iyer, 2017), where we previously argued that the hybrid equation is not supported by Bayesian reasoning.
△ Less
Submitted 10 May, 2022;
originally announced May 2022.
-
Strategies in Education, Outreach, and Inclusion to Enhance the US Workforce in Accelerator Science and Engineering
Authors:
M. Bai,
W. A. Barletta,
D. L. Bruhwiler,
S. Chattopadhyay,
Y. Hao,
S. Holder,
J. Holzbauer,
Z. Huang,
K. Harkay,
Y. -K. Kim,
X. Lu,
S. M. Lund,
N. Neveu,
P. Ostroumov,
J. R. Patterson,
P. Piot,
T. Satogata,
A. Seryi,
A. K. Soha,
S. Winchester
Abstract:
We summarize the community-based consensus for improvements concerning education, public outreach, and inclusion in Accelerator Science and Engineering that will enhance the workforce in the USA. The improvements identified reflect the product of discussions held within the 2021-2022 Snowmass community planning process by topical group AF1: Beam Physics and Accelerator Education within the Acceler…
▽ More
We summarize the community-based consensus for improvements concerning education, public outreach, and inclusion in Accelerator Science and Engineering that will enhance the workforce in the USA. The improvements identified reflect the product of discussions held within the 2021-2022 Snowmass community planning process by topical group AF1: Beam Physics and Accelerator Education within the Accelerator Frontier. Although the Snowmass process centers on high-energy physics, this document outlines required improvements for the entire U.S. accelerator science and engineering enterprise because education of those entering and in the field, outreach to the public, and inclusion are inextricably linked.
△ Less
Submitted 16 March, 2022;
originally announced March 2022.
-
Snowmass21 Accelerator Modeling Community White Paper
Authors:
S. Biedron,
L. Brouwer,
D. L. Bruhwiler,
N. M. Cook,
A. L. Edelen,
D. Filippetto,
C. -K. Huang,
A. Huebl,
T. Katsouleas,
N. Kuklev,
R. Lehe,
S. Lund,
C. Messe,
W. Mori,
C. -K. Ng,
D. Perez,
P. Piot,
J. Qiang,
R. Roussel,
D. Sagan,
A. Sahai,
A. Scheinker,
M. Thévenet,
F. Tsung,
J. -L. Vay
, et al. (2 additional authors not shown)
Abstract:
After a summary of relevant comments and recommendations from various reports over the last ten years, this paper examines the modeling needs in accelerator physics, from the modeling of single beams and individual accelerator elements, to the realization of virtual twins that replicate all the complexity to model a particle accelerator complex as accurately as possible. We then discuss cutting-ed…
▽ More
After a summary of relevant comments and recommendations from various reports over the last ten years, this paper examines the modeling needs in accelerator physics, from the modeling of single beams and individual accelerator elements, to the realization of virtual twins that replicate all the complexity to model a particle accelerator complex as accurately as possible. We then discuss cutting-edge and emerging computing opportunities, such as advanced algorithms, AI/ML and quantum computing, computational needs in hardware, software performance, portability and scalability, and needs for scalable I/O and in-situ analysis. Considerations of reliability, long-term sustainability, user support and training are considered next, before discussing the benefits of ecosystems with integrated workflows based on standardized input and output, and with integrated frameworks and data repositories developed as a community. Last, we highlight how the community can work more collaboratively and efficiently through the development of consortia and centers, and via collaboration with industry.
△ Less
Submitted 22 September, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Accelerator and Beam Physics: Grand Challenges and Research Opportunities
Authors:
S. Nagaitsev,
V. Shiltsev,
A. Valishev,
T. Zolkin,
J. -L. Vay,
M. Bai,
Y. Cai,
M. J. Hogan,
Z. Huang,
J. Seeman,
B. Dunham,
X. Huang,
T. Roser,
M. Minty,
J. Rosenzweig,
P. Piot,
J. Power,
J. M. Byrd,
A. Seryi,
S. Lund,
J. R. Patterson
Abstract:
Accelerator and beam physics (ABP) is the science of the motion, generation, acceleration, manipulation, prediction, observation and use of charged particle beams. It focuses on fundamental long-term accelerator and beam physics research and development. Accelerator and beam physics research has resulted in important advances in accelerator science, yet support for this research is declining. NSF…
▽ More
Accelerator and beam physics (ABP) is the science of the motion, generation, acceleration, manipulation, prediction, observation and use of charged particle beams. It focuses on fundamental long-term accelerator and beam physics research and development. Accelerator and beam physics research has resulted in important advances in accelerator science, yet support for this research is declining. NSF has terminated its program in accelerator Science and funding by DOE through GARD and Accelerator Stewardship has been steady or declining. The declining support for accelerator research will slow advances and threaten student training and work-force development in accelerator science. We propose a robust and scientifically challenging program in accelerator and beam physics, which will position the field of US High Energy Physics to be productive and competitive for decades to come.
△ Less
Submitted 13 March, 2022;
originally announced March 2022.
-
Dynamical bunching and density peaks in expanding Coulomb clouds
Authors:
Brandon Zerbe,
Xukun Xiang,
Chong-Yu Ruan,
Steve Lund,
Phil Duxbury
Abstract:
Expansion dynamics of single-species, non-neutral clouds, such as electron bunches used in ultrafast electron microscopy, show novel behavior due to high acceleration of particles in the cloud interior. This often leads to electron bunching and dynamical formation of a density shock in the outer regions of the bunch. We develop analytic fluid models to capture these effects, and the analytic predi…
▽ More
Expansion dynamics of single-species, non-neutral clouds, such as electron bunches used in ultrafast electron microscopy, show novel behavior due to high acceleration of particles in the cloud interior. This often leads to electron bunching and dynamical formation of a density shock in the outer regions of the bunch. We develop analytic fluid models to capture these effects, and the analytic predictions are validated by PIC and N-particle simulations. In the space-charge dominated regime, two and three dimensional systems with Gaussian initial densities show bunching and a strong shock response, while one dimensional systems do not; moreover these effects can be tuned using the initial particle density profile and velocity chirp.
△ Less
Submitted 27 December, 2017; v1 submitted 22 December, 2017;
originally announced December 2017.
-
Likelihood Ratio as Weight of Forensic Evidence: A Closer Look
Authors:
Steven P. Lund,
Hari K. Iyer
Abstract:
The forensic science community has increasingly sought quantitative methods for conveying the weight of evidence. Experts from many forensic laboratories summarize their findings in terms of a likelihood ratio. Several proponents of this approach have argued that Bayesian reasoning proves it to be normative. We find this likelihood ratio paradigm to be unsupported by arguments of Bayesian decision…
▽ More
The forensic science community has increasingly sought quantitative methods for conveying the weight of evidence. Experts from many forensic laboratories summarize their findings in terms of a likelihood ratio. Several proponents of this approach have argued that Bayesian reasoning proves it to be normative. We find this likelihood ratio paradigm to be unsupported by arguments of Bayesian decision theory, which applies only to personal decision making and not to the transfer of information from an expert to a separate decision maker. We further argue that decision theory does not exempt the presentation of a likelihood ratio from uncertainty characterization, which is required to assess the fitness for purpose of any transferred quantity. We propose the concept of a lattice of assumptions leading to an uncertainty pyramid as a framework for assessing the uncertainty in an evaluation of a likelihood ratio. We demonstrate the use of these concepts with illustrative examples regarding the refractive index of glass and automated comparison scores for fingerprints.
△ Less
Submitted 26 April, 2017;
originally announced April 2017.
-
Likelihood Ratio as Weight of Forensic Evidence: A Metrological Perspective
Authors:
Steven P. Lund,
Hari Iyer
Abstract:
In this article we provide a rebuttal against the possible perception that a single number, such as the Likelihood Ratio, can provide an objective, authoritative or definitive weight of evidence. We also illustrate the extent to which conclusions can vary depending on the assumptions used in the analysis, even under alternative assumptions that are judged to be consistent with available empirical…
▽ More
In this article we provide a rebuttal against the possible perception that a single number, such as the Likelihood Ratio, can provide an objective, authoritative or definitive weight of evidence. We also illustrate the extent to which conclusions can vary depending on the assumptions used in the analysis, even under alternative assumptions that are judged to be consistent with available empirical information. To facilitate these goals, we introduce the notion of a Lattice of Assumptions and an Uncertainty Pyramid illustrated in the context of a previously published example involving glass evidence. We take the position that rather than focusing on a single number summary as the weight of evidence it is the duty of the forensic expert to assist the trier of fact in forming their own interpretations from a clear understanding of the objective and demonstrably available information. We hope the presented arguments will inspire those in the forensic science community to pursue establishing their practice on a solid foundation of measurement science.
△ Less
Submitted 19 September, 2016; v1 submitted 26 August, 2016;
originally announced August 2016.
-
Fusion of Array Operations at Runtime
Authors:
Mads R. B. Kristensen,
Simon A. F. Lund,
Troels Blum,
James Avery
Abstract:
We address the problem of fusing array operations based on criteria such as shape compatibility, data reusability, and communication. We formulate the problem as a graph partition problem that is general enough to handle loop fusion, combinator fusion, and other types of subroutines.
We address the problem of fusing array operations based on criteria such as shape compatibility, data reusability, and communication. We formulate the problem as a graph partition problem that is general enough to handle loop fusion, combinator fusion, and other types of subroutines.
△ Less
Submitted 21 January, 2016; v1 submitted 20 January, 2016;
originally announced January 2016.
-
Assessing Technical Performance in Differential Gene Expression Experiments with External Spike-in RNA Control Ratio Mixtures
Authors:
Sarah A. Munro,
Steve P. Lund,
P. Scott Pine,
Hans Binder,
Djork-Arné Clevert,
Ana Conesa,
Joaquin Dopazo,
Mario Fasold,
Sepp Hochreiter,
Huixiao Hong,
Nederah Jafari,
David P. Kreil,
Paweł P. Łabaj,
Sheng Li,
Yang Liao,
Simon Lin,
Joseph Meehan,
Christopher E. Mason,
Javier Santoyo,
Robert A. Setterquist,
Leming Shi,
Wei Shi,
Gordon K. Smyth,
Nancy Stralis-Pavese,
Zhenqiang Su
, et al. (8 additional authors not shown)
Abstract:
There is a critical need for standard approaches to assess, report, and compare the technical performance of genome-scale differential gene expression experiments. We assess technical performance with a proposed "standard" dashboard of metrics derived from analysis of external spike-in RNA control ratio mixtures. These control ratio mixtures with defined abundance ratios enable assessment of diagn…
▽ More
There is a critical need for standard approaches to assess, report, and compare the technical performance of genome-scale differential gene expression experiments. We assess technical performance with a proposed "standard" dashboard of metrics derived from analysis of external spike-in RNA control ratio mixtures. These control ratio mixtures with defined abundance ratios enable assessment of diagnostic performance of differentially expressed transcript lists, limit of detection of ratio (LODR) estimates, and expression ratio variability and measurement bias. The performance metrics suite is applicable to analysis of a typical experiment, and here we also apply these metrics to evaluate technical performance among laboratories. An interlaboratory study using identical samples shared amongst 12 laboratories with three different measurement processes demonstrated generally consistent diagnostic power across 11 laboratories. Ratio measurement variability and bias were also comparable amongst laboratories for the same measurement process. Different biases were observed for measurement processes using different mRNA enrichment protocols.
△ Less
Submitted 18 June, 2014;
originally announced June 2014.
-
Estimating robustness of the tileShuffle method with repeated probes
Authors:
Sigrun Helga Lund,
Asgeir Sigurdsson,
Sigurjon Axel Gudjonsson,
Julius Gudmundsson,
Daniel Fannar Gudbjartsson,
Thorunn Rafnar,
Kari Stefansson,
Gunnar Stefansson
Abstract:
In this paper the TileShuffle method is evaluated as a search method for candidate lncRNAs at 8q24.2. The method is run on three microarrays. Microarrays which all contained the same sample and repeated copies of tiled probes. This allows the coherence of the selection method within and between microarrays to be estimated by Monte Carlo simulations on the repeated probes.
In this paper the TileShuffle method is evaluated as a search method for candidate lncRNAs at 8q24.2. The method is run on three microarrays. Microarrays which all contained the same sample and repeated copies of tiled probes. This allows the coherence of the selection method within and between microarrays to be estimated by Monte Carlo simulations on the repeated probes.
△ Less
Submitted 12 March, 2014;
originally announced March 2014.
-
Noiseless Vlasov-Poisson simulations with linearly transformed particles
Authors:
M. Campos Pinto,
E. Sonnendrücker,
A. Friedman,
D. Grote,
S. Lund
Abstract:
We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts…
▽ More
We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts to using deformation matrices within the particle shape functions; these matrices are updated via local evaluations of the forward numerical flow. Because it is necessary to periodically remap the particles on a regular grid to avoid excessively deforming their shapes, the method can be seen as a development of Denavit's Forward Semi-Lagrangian (FSL) scheme [J. Denavit, J. Comp. Physics 9, 75 (1972)]. However, it has recently been established [M. Campos Pinto, "Smooth particle methods without smoothing", arXiv:1112.1859 (2012)] that the underlying Linearly-Transformed Particle scheme converges for abstract transport problems, with no need to remap the particles; deforming the particles can thus be seen as a way to significantly lower the remapping frequency needed in the FSL schemes, and hence the associated numerical diffusion. To couple the method with electrostatic field solvers, two specific charge deposition schemes are examined, and their performance compared with that of the standard deposition method. Finally, numerical 1d1v simulations involving benchmark test cases and halo formation in an initially mismatched thermal sheet beam demonstrate some advantages of our LTPIC scheme over the classical PIC and FSL methods. Benchmarked test cases also indicate that, for numerical choices involving similar computational effort, the LTPIC method is capable of accuracy comparable to or exceeding that of state-of-the-art, high-resolution Vlasov schemes.
△ Less
Submitted 29 March, 2013; v1 submitted 21 November, 2012;
originally announced November 2012.
-
cphVB: A System for Automated Runtime Optimization and Parallelization of Vectorized Applications
Authors:
Mads Ruben Burgdorff Kristensen,
Simon Andreas Frimann Lund,
Troels Blum,
Brian Vinter
Abstract:
Modern processor architectures, in addition to having still more cores, also require still more consideration to memory-layout in order to run at full capacity. The usefulness of most languages is deprecating as their abstractions, structures or objects are hard to map onto modern processor architectures efficiently.
The work in this paper introduces a new abstract machine framework, cphVB, that…
▽ More
Modern processor architectures, in addition to having still more cores, also require still more consideration to memory-layout in order to run at full capacity. The usefulness of most languages is deprecating as their abstractions, structures or objects are hard to map onto modern processor architectures efficiently.
The work in this paper introduces a new abstract machine framework, cphVB, that enables vector oriented high-level programming languages to map onto a broad range of architectures efficiently. The idea is to close the gap between high-level languages and hardware optimized low-level implementations. By translating high-level vector operations into an intermediate vector bytecode, cphVB enables specialized vector engines to efficiently execute the vector operations.
The primary success parameters are to maintain a complete abstraction from low-level details and to provide efficient code execution across different, modern, processors. We evaluate the presented design through a setup that targets multi-core CPU architectures. We evaluate the performance of the implementation using Python implementations of well-known algorithms: a jacobi solver, a kNN search, a shallow water simulation and a synthetic stencil simulation. All demonstrate good performance.
△ Less
Submitted 25 March, 2013; v1 submitted 26 October, 2012;
originally announced October 2012.
-
The importance of distinct modeling strategies for gene and gene-specific treatment effects in hierarchical models for microarray data
Authors:
Steven P. Lund,
Dan Nettleton
Abstract:
When analyzing microarray data, hierarchical models are often used to share information across genes when estimating means and variances or identifying differential expression. Many methods utilize some form of the two-level hierarchical model structure suggested by Kendziorski et al. [Stat. Med. (2003) 22 3899-3914] in which the first level describes the distribution of latent mean expression lev…
▽ More
When analyzing microarray data, hierarchical models are often used to share information across genes when estimating means and variances or identifying differential expression. Many methods utilize some form of the two-level hierarchical model structure suggested by Kendziorski et al. [Stat. Med. (2003) 22 3899-3914] in which the first level describes the distribution of latent mean expression levels among genes and among differentially expressed treatments within a gene. The second level describes the conditional distribution, given a latent mean, of repeated observations for a single gene and treatment. Many of these models, including those used in Kendziorski's et al. [Stat. Med. (2003) 22 3899-3914] EBarrays package, assume that expression level changes due to treatment effects have the same distribution as expression level changes from gene to gene. We present empirical evidence that this assumption is often inadequate and propose three-level hierarchical models as extensions to the two-level log-normal based EBarrays models to address this inadequacy. We demonstrate that use of our three-level models dramatically changes analysis results for a variety of microarray data sets and verify the validity and improved performance of our suggested method in a series of simulation studies. We also illustrate the importance of accounting for the uncertainty of gene-specific error variance estimates when using hierarchical models to identify differentially expressed genes.
△ Less
Submitted 28 September, 2012;
originally announced September 2012.
-
Space-charge transport limits of ion beams in periodic quadrupole focusing channels
Authors:
Steven M. Lund,
Sugreev R. Chawla
Abstract:
It has been empirically observed in both experiments and particle-in-cell simulations that space-charge-dominated beams suffer strong growth in statistical phase-space area (degraded quality) and particle losses in alternating gradient quadrupole transport channels when the undepressed phase advance sigma_0 increases beyond about 85 degrees per lattice period. Although this criterion has been us…
▽ More
It has been empirically observed in both experiments and particle-in-cell simulations that space-charge-dominated beams suffer strong growth in statistical phase-space area (degraded quality) and particle losses in alternating gradient quadrupole transport channels when the undepressed phase advance sigma_0 increases beyond about 85 degrees per lattice period. Although this criterion has been used extensively in practical designs of strong focusing intense beam transport lattices, the origin of the limit has not been understood. We propose a mechanism for the transport limit resulting from classes of halo particle resonances near the core of the beam that allow near-edge particles to rapidly increase in oscillation amplitude when the space-charge intensity and the flutter of the matched beam envelope are both sufficiently large. When coupled with a diffuse beam edge and/or perturbations internal to the beam core that can drive particles outside the edge, this mechanism can result in large and rapid halo-driven increases in the statistical phase-space area of the beam, lost particles, and degraded transport. A core-particle model is applied to parametrically analyze this process. Extensive self-consistent particle in cell simulations are employed to better quantify space-charge limit and verify core-particle model predictions.
△ Less
Submitted 27 February, 2006;
originally announced February 2006.
-
Efficient computation of matched solutions of the Kapchinskij-Vladimirskij envelope equations for periodic focusing lattices
Authors:
Steven M. Lund,
Sven H. Chilton,
Edward P. Lee
Abstract:
A new iterative method is developed to numerically calculate the periodic, matched beam envelope solution of the coupled Kapchinskij-Vladimirskij (KV) equations describing the transverse evolution of a beam in a periodic, linear focusing lattice of arbitrary complexity. Implementation of the method is straightforward. It is highly convergent and can be applied to all usual parameterizations of t…
▽ More
A new iterative method is developed to numerically calculate the periodic, matched beam envelope solution of the coupled Kapchinskij-Vladimirskij (KV) equations describing the transverse evolution of a beam in a periodic, linear focusing lattice of arbitrary complexity. Implementation of the method is straightforward. It is highly convergent and can be applied to all usual parameterizations of the matched envelope solutions. The method is applicable to all classes of linear focusing lattices without skew couplings, and also applies to all physically achievable system parameters -- including where the matched beam envelope is strongly unstable. Example applications are presented for periodic solenoidal and quadrupole focusing lattices. Convergence properties are summarized over a wide range of system parameters.
△ Less
Submitted 22 February, 2006;
originally announced February 2006.
-
Artificial "spin ice" in a geometrically frustrated lattice of nanoscale ferromagnetic islands
Authors:
R. F. Wang,
C. Nisoli,
R. S. Freitas,
J. Li,
W. McConville,
B. J. Cooley,
M. S. Lund,
N. Samarth,
C. Leighton,
V. H. Crespi,
P. Schiffer
Abstract:
We report an artificial geometrically frustrated magnet based on an array of lithographically fabricated single-domain ferromagnetic islands. The islands are arranged such that the dipole interactions create a two-dimensional analogue to spin ice. Images of the magnetic moments of individual elements in this correlated system allow us to study the local accommodation of frustration. We see both…
▽ More
We report an artificial geometrically frustrated magnet based on an array of lithographically fabricated single-domain ferromagnetic islands. The islands are arranged such that the dipole interactions create a two-dimensional analogue to spin ice. Images of the magnetic moments of individual elements in this correlated system allow us to study the local accommodation of frustration. We see both ice-like short-range correlations and an absence of long-range correlations, behaviour which is strikingly similar to the lowtemperature state of spin ice. These results demonstrate that artificial frustrated magnets can provide an uncharted arena in which the physics of frustration can be directly visualized.
△ Less
Submitted 19 January, 2006;
originally announced January 2006.
-
Dipole Septum Magnet in the Fast Kicker System for Multi-Axis Advanced Radiography
Authors:
L. Wang,
S. M. Lund,
B. R. Poole
Abstract:
Here we present designs for a static septum magnet with two adjacent apertures where ideally one aperture has a uniform dipole field and the other zero field. Two designs are considered. One is a true septum magnet with a thin layer of coils and materials separating the dipole field region from the null field region. During the beam switching process, the intense electron beam will spray across…
▽ More
Here we present designs for a static septum magnet with two adjacent apertures where ideally one aperture has a uniform dipole field and the other zero field. Two designs are considered. One is a true septum magnet with a thin layer of coils and materials separating the dipole field region from the null field region. During the beam switching process, the intense electron beam will spray across this material septum leading to concerns on beam control, vacuum quality, radiation damage, etc. due to the lost particles. Therefore, another configuration without a material septum is also considered. With this configuration it is more difficult to achieve high field quality near the transition region. Shaped shims are designed to limit the degradation of beam quality (emittance growth). This approach is closely related to a previous septum magnet design with two oppositely oriented dipole field regions presented by the authors [1]. Simulations are performed to obtain the magnetic field profile in both designs. A PIC simulation is used to transport a beam slice consisting of several thousand particles through the magnet to estimate emittance growth in the magnet due to field non-uniformity.
△ Less
Submitted 2 October, 2000;
originally announced October 2000.
-
MOE11 Emittance Growth from the Thermalization of Space-Charge Nonuniformities
Authors:
Steven M. Lund,
John J. Barnard,
Edward P. Lee
Abstract:
Beams injected into a linear focusing channel typically have some degree of space-charge nonuniformity. In general, injected particle distributions with systematic charge nonuniformities are not equilibria of the focusing channel and launch a broad spectrum of collective modes. These modes can phase-mix and have nonlinear wave-wave interactions which, at high space-charge intensities, resu…
▽ More
Beams injected into a linear focusing channel typically have some degree of space-charge nonuniformity. In general, injected particle distributions with systematic charge nonuniformities are not equilibria of the focusing channel and launch a broad spectrum of collective modes. These modes can phase-mix and have nonlinear wave-wave interactions which, at high space-charge intensities, results in a relaxation to a more thermal-like distribution characterized by a uniform density profile. This thermalization can transfer self-field energy from the initial space-charge nonuniformity to the local particle temperature, thereby increasing beam phase space area (emittance growth). In this paper, we employ a simple kinetic model of a continuous focusing channel and build on previous work that applied system energy and charge conservation quantify emittance growth associated with the collective thermalization of an initial azimuthally symmetric, rms matched beam with a radial density profile that is hollowed or peaked. This emittance growth is shown to be surprisingly modest even for high beam intensities with significant radial structure in the initial density profile.
△ Less
Submitted 29 September, 2000;
originally announced September 2000.