Showing 1–50 of 74 results for author: Holmes, S

Effects of Strain Compensation on Electron Mobilities in InAs Quantum Wells Grown on InP(001)

Authors: C. P. Dempsey, J. T. Dong, I. Villar Rodriguez, Y. Gul, S. Chatterjee, M. Pendharkar, S. N. Holmes, M. Pepper, C. J. Palmstrøm

Abstract: InAs quantum wells (QWs) grown on InP substrates are interesting for their applications in devices with high spin-orbit coupling (SOC) and their potential role in creating topologically nontrivial hybrid heterostructures. The highest mobility QWs are limited by interfacial roughness scattering and alloy disorder scattering in the cladding and buffer layers. Increasing QW thickness has been shown t… ▽ More InAs quantum wells (QWs) grown on InP substrates are interesting for their applications in devices with high spin-orbit coupling (SOC) and their potential role in creating topologically nontrivial hybrid heterostructures. The highest mobility QWs are limited by interfacial roughness scattering and alloy disorder scattering in the cladding and buffer layers. Increasing QW thickness has been shown to reduce the effect of both of these scattering mechanisms. However, for current state-of-the-art devices with As-based cladding and barrier layers, the critical thickness is limited to $\leq7$ nm. In this report, we demonstrate the use of strain compensation techniques in the InGaAs cladding layers to extend the critical thickness well beyond this limit. We induce tensile strain in the InGaAs cladding layers by reducing the In concentration from In$_{0.81}$Ga$_{0.19}$As to In$_{0.70}$Ga$_{0.30}$As and we observe changes in both the critical thickness of the well and the maximum achievable mobility. The peak electron mobility at 2 K is $1.16\times10^6$ cm$^2/$Vs, with a carrier density of $4.2\times10^{11}$ /cm$^2$. Additionally, we study the quantum lifetime and Rashba spin splitting in the highest mobility device as these parameters are critical to determine if these structures can be used in topologically nontrivial devices. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Development of novel ionization chambers for reference dosimetry in electron FLASH radiotherapy

Authors: Kevin Liu, Shannon Holmes, Ahtesham Ullah Khan, Brian Hooten, Larry DeWerd, Emil Schüler, Sam Beddar

Abstract: The aim of this study was to optimize the design and performance of parallel plate ion chambers for use in ultra-high dose rate (UHDR) dosimetry applications, and evaluate their potential as reference class chambers for calibration purposes. Three chambers were designed and produced: the A11-VAR (0.2-1.0 mm electrode gap, 20 mm diameter collector), the A11-TPP (0.3 mm electrode gap, 20 mm diameter… ▽ More The aim of this study was to optimize the design and performance of parallel plate ion chambers for use in ultra-high dose rate (UHDR) dosimetry applications, and evaluate their potential as reference class chambers for calibration purposes. Three chambers were designed and produced: the A11-VAR (0.2-1.0 mm electrode gap, 20 mm diameter collector), the A11-TPP (0.3 mm electrode gap, 20 mm diameter collector), and the A30 (0.3 mm electrode gap, 5.4 mm diameter collector).The chambers underwent full characterization using an UHDR 9 MeV electron beam with individually varied beam parameters of pulse repetition frequency (PRF, 10-120Hz), pulse width (PW, 0.5-4us), and pulse amplitude (0.01-9 Gy/pulse). The response of the ion chambers was evaluated as a function of the dose per pulse (DPP), PRF, PW, dose rate, electric field strength, and electrode gap. The chamber response was found to be dependent on DPP and PW, whose dependencies were mitigated with larger electric field strengths and smaller electrode spacing. At a constant electric field strength, we measured a larger charge collection efficiency (CCE) as a function of DPP for ion chambers with a smaller electrode gap in the A11-VAR. For ion chambers with identical electrode gap (A11-TPP and A30), higher electric field strengths were found to yield better CCE at higher DPP. A PW dependence was observed at low electric field strengths (500 V/mm) for DPP values ranging from 1-5 Gy at PWs ranging from 0.5-4 μs, but at electric field strengths of 1000 V/mm and higher, these effects become negligible. This study confirmed that the charge collection efficiency of ion chambers depends strongly on the electrode spacing and the electric field strength, and also on the DPP and the PW of the UHDR beam. The new finding of this study is that the PW dependence becomes negligible with reduced electrode spacing and increased electric field. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 29 pages, 9 figures

Generative Models: An Interdisciplinary Perspective

Authors: Kris Sankaran, Susan P. Holmes

Abstract: By linking conceptual theories with observed data, generative models can support reasoning in complex situations. They have come to play a central role both within and beyond statistics, providing the basis for power analysis in molecular biology, theory building in particle physics, and resource allocation in epidemiology, for example. We introduce the probabilistic and computational concepts und… ▽ More By linking conceptual theories with observed data, generative models can support reasoning in complex situations. They have come to play a central role both within and beyond statistics, providing the basis for power analysis in molecular biology, theory building in particle physics, and resource allocation in epidemiology, for example. We introduce the probabilistic and computational concepts underlying modern generative models and then analyze how they can be used to inform experimental design, iterative model refinement, goodness-of-fit evaluation, and agent-based simulation. We emphasize a modular view of generative mechanisms and discuss how they can be flexibly recombined in new problem contexts. We provide practical illustrations throughout, and code for reproducing all examples is available at https://github.com/krisrs1128/generative_review. Finally, we observe how research in generative models is currently split across several islands of activity, and we highlight opportunities lying at disciplinary intersections. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Statistical Modeling for Practical Pooled Testing During the COVID-19 Pandemic

Authors: Saskia Comess, Hannah Wang, Susan Holmes, Claire Donnat

Abstract: Pooled testing offers an efficient solution to the unprecedented testing demands of the COVID-19 pandemic, although with potentially lower sensitivity and increased costs to implementation in some settings. Assessments of this trade-off typically assume pooled specimens are independent and identically distributed. Yet, in the context of COVID-19, these assumptions are often violated: testing done… ▽ More Pooled testing offers an efficient solution to the unprecedented testing demands of the COVID-19 pandemic, although with potentially lower sensitivity and increased costs to implementation in some settings. Assessments of this trade-off typically assume pooled specimens are independent and identically distributed. Yet, in the context of COVID-19, these assumptions are often violated: testing done on networks (housemates, spouses, co-workers) captures correlated individuals, while infection risk varies substantially across time, place and individuals. Neglecting dependencies and heterogeneity may bias established optimality grids and induce a sub-optimal implementation of the procedure. As a lesson learned from this pandemic, this paper highlights the necessity of integrating field sampling information with statistical modeling to efficiently optimize pooled testing. Using real data, we show that (a) greater gains can be achieved at low logistical cost by exploiting natural correlations (non-independence) between samples -- allowing improvements in sensitivity and efficiency of up to 30% and 90% respectively; and (b) these gains are robust despite substantial heterogeneity across pools (non-identical). Our modeling results complement and extend the observations of Barak et al (2021) who report an empirical sensitivity well beyond expectations. Finally, we provide an interactive tool for selecting an optimal pool size using contextual information △ Less

Submitted 12 July, 2021; originally announced July 2021.

A Statistical Perspective on the Challenges in Molecular Microbial Biology

Authors: Pratheepa Jeganathan, Susan P. Holmes

Abstract: High throughput sequencing (HTS)-based technology enables identifying and quantifying non-culturable microbial organisms in all environments. Microbial sequences have enhanced our understanding of the human microbiome, the soil and plant environment, and the marine environment. All molecular microbial data pose statistical challenges due to contamination sequences from reagents, batch effects, une… ▽ More High throughput sequencing (HTS)-based technology enables identifying and quantifying non-culturable microbial organisms in all environments. Microbial sequences have enhanced our understanding of the human microbiome, the soil and plant environment, and the marine environment. All molecular microbial data pose statistical challenges due to contamination sequences from reagents, batch effects, unequal sampling, and undetected taxa. Technical biases and heteroscedasticity have the strongest effects, but different strains across subjects and environments also make direct differential abundance testing unwieldy. We provide an introduction to a few statistical tools that can overcome some of these difficulties and demonstrate those tools on an example. We show how standard statistical methods, such as simple hierarchical mixture and topic models, can facilitate inferences on latent microbial communities. We also review some nonparametric Bayesian approaches that combine visualization and uncertainty quantification. The intersection of molecular microbial biology and statistics is an exciting new venue. Finally, we list some of the important open problems that would benefit from more careful statistical method development. △ Less

Submitted 6 March, 2021; originally announced March 2021.

Comments: To appear in the Journal of Agricultural, Biological and Environmental Statistics

Non-planar geometrical effects on the magnetoelectrical signal in a three-dimensional nanomagnetic circuit

Authors: Fanfan Meng, Claire Donnelly, Claas Abert, Luka Skoric, Stuart Holmes, Zhuocong Xiao, Jung-Wei Liao, Peter J. Newton, Crispin H. W. Barnes, Dédalo Sanz-Hernández, Aurelio Hierro-Rodriguez, Dieter Suess, Russell P. Cowburn, Amalio Fernández-Pacheco

Abstract: Expanding nanomagnetism and spintronics into three dimensions (3D) offers great opportunities for both fundamental and technological studies. However, probing the influence of complex 3D geometries on magnetoelectrical phenomena poses important experimental and theoretical challenges. In this work, we investigate the magnetoelectrical signals of a ferromagnetic 3D nanodevice integrated into a micr… ▽ More Expanding nanomagnetism and spintronics into three dimensions (3D) offers great opportunities for both fundamental and technological studies. However, probing the influence of complex 3D geometries on magnetoelectrical phenomena poses important experimental and theoretical challenges. In this work, we investigate the magnetoelectrical signals of a ferromagnetic 3D nanodevice integrated into a microelectronic circuit using direct-write nanofabrication. Due to the 3D vectorial nature of both electrical current and magnetisation, a complex superposition of several magnetoelectrical effects takes place. By performing electrical measurements under the application of 3D magnetic fields, in combination with macrospin simulations and finite element modelling, we disentangle the superimposed effects, finding how a 3D geometry leads to unusual angular dependences of well-known magnetotransport effects such as the anomalous Hall effect. Crucially, our analysis also reveals a strong role of the noncollinear demagnetising fields intrinsic to 3D nanostructures, which results in an angular dependent magnon magnetoresistance contributing strongly to the total magnetoelectrical signal. These findings are key to the understanding of 3D spintronic systems and underpin further fundamental and device-based studies. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Modeling the Heterogeneity in COVID-19's Reproductive Number and its Impact on Predictive Scenarios

Authors: Claire Donnat, Susan Holmes

Abstract: The correct evaluation of the reproductive number $R$ for COVID-19 -- which characterizes the average number of secondary cases generated by each typical primary case -- is central in the quantification of the potential scope of the pandemic and the selection of an appropriate course of action. In most models, $R$ is modeled as a universal constant for the virus across outbreak clusters and indivi… ▽ More The correct evaluation of the reproductive number $R$ for COVID-19 -- which characterizes the average number of secondary cases generated by each typical primary case -- is central in the quantification of the potential scope of the pandemic and the selection of an appropriate course of action. In most models, $R$ is modeled as a universal constant for the virus across outbreak clusters and individuals -- effectively averaging out the inherent variability of the transmission process due to varying individual contact rates, population densities, demographics, or temporal factors amongst many. Yet, due to the exponential nature of epidemic growth, the error due to this simplification can be rapidly amplified and lead to inaccurate predictions and/or risk evaluation. From the statistical modeling perspective, the magnitude of the impact of this averaging remains an open question: how can this intrinsic variability be percolated into epidemic models, and how can its impact on uncertainty quantification and predictive scenarios be better quantified? In this paper, we propose to study this question through a Bayesian perspective, creating a bridge between the agent-based and compartmental approaches commonly used in the literature. After deriving a Bayesian model that captures at scale the heterogeneity of a population and environmental conditions, we simulate the spread of the epidemic as well as the impact of different social distancing strategies, and highlight the strong impact of this added variability on the reported results. We base our discussion on both synthetic experiments -- thereby quantifying of the reliability and the magnitude of the effects -- and real COVID-19 data. △ Less

Submitted 14 April, 2021; v1 submitted 10 April, 2020; originally announced April 2020.

Geomstats: A Python Package for Riemannian Geometry in Machine Learning

Authors: Nina Miolane, Alice Le Brigant, Johan Mathe, Benjamin Hou, Nicolas Guigui, Yann Thanwerdas, Stefan Heyder, Olivier Peltre, Niklas Koep, Hadi Zaatiti, Hatem Hajri, Yann Cabanes, Thomas Gerald, Paul Chauchat, Christian Shewmake, Bernhard Kainz, Claire Donnat, Susan Holmes, Xavier Pennec

Abstract: We introduce Geomstats, an open-source Python toolbox for computations and statistics on nonlinear manifolds, such as hyperbolic spaces, spaces of symmetric positive definite matrices, Lie groups of transformations, and many more. We provide object-oriented and extensively unit-tested implementations. Among others, manifolds come equipped with families of Riemannian metrics, with associated expone… ▽ More We introduce Geomstats, an open-source Python toolbox for computations and statistics on nonlinear manifolds, such as hyperbolic spaces, spaces of symmetric positive definite matrices, Lie groups of transformations, and many more. We provide object-oriented and extensively unit-tested implementations. Among others, manifolds come equipped with families of Riemannian metrics, with associated exponential and logarithmic maps, geodesics and parallel transport. Statistics and learning algorithms provide methods for estimation, clustering and dimension reduction on manifolds. All associated operations are vectorized for batch computation and provide support for different execution backends, namely NumPy, PyTorch and TensorFlow, enabling GPU acceleration. This paper presents the package, compares it with related libraries and provides relevant code examples. We show that Geomstats provides reliable building blocks to foster research in differential geometry and statistics, and to democratize the use of Riemannian geometry in machine learning applications. The source code is freely available under the MIT license at \url{geomstats.ai}. △ Less

Submitted 7 April, 2020; originally announced April 2020.

AutoCogniSys: IoT Assisted Context-Aware Automatic Cognitive Health Assessment

Authors: Mohammad Arif Ul Alam, Nirmalya Roy, Sarah Holmes, Aryya Gangopadhyay, Elizabeth Galik

Abstract: Cognitive impairment has become epidemic in older adult population. The recent advent of tiny wearable and ambient devices, a.k.a Internet of Things (IoT) provides ample platforms for continuous functional and cognitive health assessment of older adults. In this paper, we design, implement and evaluate AutoCogniSys, a context-aware automated cognitive health assessment system, combining the sensin… ▽ More Cognitive impairment has become epidemic in older adult population. The recent advent of tiny wearable and ambient devices, a.k.a Internet of Things (IoT) provides ample platforms for continuous functional and cognitive health assessment of older adults. In this paper, we design, implement and evaluate AutoCogniSys, a context-aware automated cognitive health assessment system, combining the sensing powers of wearable physiological (Electrodermal Activity, Photoplethysmography) and physical (Accelerometer, Object) sensors in conjunction with ambient sensors. We design appropriate signal processing and machine learning techniques, and develop an automatic cognitive health assessment system in a natural older adults living environment. We validate our approaches using two datasets: (i) a naturalistic sensor data streams related to Activities of Daily Living and mental arousal of 22 older adults recruited in a retirement community center, individually living in their own apartments using a customized inexpensive IoT system (IRB #HP-00064387) and (ii) a publicly available dataset for emotion detection. The performance of AutoCogniSys attests max. 93\% of accuracy in assessing cognitive health of older adults. △ Less

Submitted 16 March, 2020; originally announced March 2020.

MSC Class: 62H20 ACM Class: I.2.7

Arcsine laws for random walks generated from random permutations with applications to genomics

Authors: Xiao Fang, Han Liang Gan, Susan Holmes, Haiyan Huang, Erol Peköz, Adrian Röllin, Wenpin Tang

Abstract: A classical result for the simple symmetric random walk with $2n$ steps is that the number of steps above the origin, the time of the last visit to the origin, and the time of the maximum height all have exactly the same distribution and converge when scaled to the arcsine law. Motivated by applications in genomics, we study the distributions of these statistics for the non-Markovian random walk g… ▽ More A classical result for the simple symmetric random walk with $2n$ steps is that the number of steps above the origin, the time of the last visit to the origin, and the time of the maximum height all have exactly the same distribution and converge when scaled to the arcsine law. Motivated by applications in genomics, we study the distributions of these statistics for the non-Markovian random walk generated from the ascents and descents of a uniform random permutation and a Mallows($q$) permutation and show that they have the same asymptotic distributions as for the simple random walk. We also give an unexpected conjecture, along with numerical evidence and a partial proof in special cases, for the result that the number of steps above the origin by step $2n$ for the uniform permutation generated walk has exactly the same discrete arcsine distribution as for the simple random walk, even though the other statistics for these walks have very different laws. We also give explicit error bounds to the limit theorems using Stein's method for the arcsine distribution, as well as functional central limit theorems and a strong embedding of the Mallows$(q)$ permutation which is of independent interest. △ Less

Submitted 23 January, 2020; originally announced January 2020.

MSC Class: 60C05; 60J65; 05A05

Learning Weighted Submanifolds with Variational Autoencoders and Riemannian Variational Autoencoders

Authors: Nina Miolane, Susan Holmes

Abstract: Manifold-valued data naturally arises in medical imaging. In cognitive neuroscience, for instance, brain connectomes base the analysis of coactivation patterns between different brain regions on the analysis of the correlations of their functional Magnetic Resonance Imaging (fMRI) time series - an object thus constrained by construction to belong to the manifold of symmetric positive definite matr… ▽ More Manifold-valued data naturally arises in medical imaging. In cognitive neuroscience, for instance, brain connectomes base the analysis of coactivation patterns between different brain regions on the analysis of the correlations of their functional Magnetic Resonance Imaging (fMRI) time series - an object thus constrained by construction to belong to the manifold of symmetric positive definite matrices. One of the challenges that naturally arises consists of finding a lower-dimensional subspace for representing such manifold-valued data. Traditional techniques, like principal component analysis, are ill-adapted to tackle non-Euclidean spaces and may fail to achieve a lower-dimensional representation of the data - thus potentially pointing to the absence of lower-dimensional representation of the data. However, these techniques are restricted in that: (i) they do not leverage the assumption that the connectomes belong on a pre-specified manifold, therefore discarding information; (ii) they can only fit a linear subspace to the data. In this paper, we are interested in variants to learn potentially highly curved submanifolds of manifold-valued data. Motivated by the brain connectomes example, we investigate a latent variable generative model, which has the added benefit of providing us with uncertainty estimates - a crucial quantity in the medical applications we are considering. While latent variable models have been proposed to learn linear and nonlinear spaces for Euclidean data, or geodesic subspaces for manifold data, no intrinsic latent variable model exists to learn nongeodesic subspaces for manifold data. This paper fills this gap and formulates a Riemannian variational autoencoder with an intrinsic generative model of manifold-valued data. We evaluate its performances on synthetic and real datasets by introducing the formalism of weighted Riemannian submanifolds. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Estimation of Orientation and Camera Parameters from Cryo-Electron Microscopy Images with Variational Autoencoders and Generative Adversarial Networks

Authors: Nina Miolane, Frédéric Poitevin, Yee-Ting Li, Susan Holmes

Abstract: Cryo-electron microscopy (cryo-EM) is capable of producing reconstructed 3D images of biomolecules at near-atomic resolution. As such, it represents one of the most promising imaging techniques in structural biology. However, raw cryo-EM images are only highly corrupted - noisy and band-pass filtered - 2D projections of the target 3D biomolecules. Reconstructing the 3D molecular shape starts with… ▽ More Cryo-electron microscopy (cryo-EM) is capable of producing reconstructed 3D images of biomolecules at near-atomic resolution. As such, it represents one of the most promising imaging techniques in structural biology. However, raw cryo-EM images are only highly corrupted - noisy and band-pass filtered - 2D projections of the target 3D biomolecules. Reconstructing the 3D molecular shape starts with the removal of image outliers, the estimation of the orientation of the biomolecule that has produced the given 2D image, and the estimation of camera parameters to correct for intensity defects. Current techniques performing these tasks are often computationally expensive, while the dataset sizes keep growing. There is a need for next-generation algorithms that preserve accuracy while improving speed and scalability. In this paper, we combine variational autoencoders (VAEs) and generative adversarial networks (GANs) to learn a low-dimensional latent representation of cryo-EM images. We perform an exploratory analysis of the obtained latent space, that is shown to have a structure of "orbits", in the sense of Lie group theory, consistent with the acquisition procedure of cryo-EM images. This analysis leads us to design an estimation method for orientation and camera parameters of single-particle cryo-EM images, together with an outliers detection procedure. As such, it opens the door to geometric approaches for unsupervised estimations of orientations and camera parameters, making possible fast cryo-EM biomolecule reconstruction. △ Less

Submitted 23 May, 2021; v1 submitted 19 November, 2019; originally announced November 2019.

Constrained Bayesian ICA for Brain Connectome Inference

Authors: Claire Donnat, Leonardo Tozzi, Susan Holmes

Abstract: Brain connectomics is a developing field in neurosciences which strives to understand cognitive processes and psychiatric diseases through the analysis of interactions between brain regions. However, in the high-dimensional, low-sample, and noisy regimes that typically characterize fMRI data, the recovery of such interactions remains an ongoing challenge: how can we discover patterns of co-activit… ▽ More Brain connectomics is a developing field in neurosciences which strives to understand cognitive processes and psychiatric diseases through the analysis of interactions between brain regions. However, in the high-dimensional, low-sample, and noisy regimes that typically characterize fMRI data, the recovery of such interactions remains an ongoing challenge: how can we discover patterns of co-activity between brain regions that could then be associated to cognitive processes or psychiatric disorders? In this paper, we investigate a constrained Bayesian ICA approach which, in comparison to current methods, simultaneously allows (a) the flexible integration of multiple sources of information (fMRI, DWI, anatomical, etc.), (b) an automatic and parameter-free selection of the appropriate sparsity level and number of connected submodules and (c) the provision of estimates on the uncertainty of the recovered interactions. Our experiments, both on synthetic and real-life data, validate the flexibility of our method and highlight the benefits of integrating anatomical information for connectome inference. △ Less

Submitted 13 November, 2019; originally announced November 2019.

Convex Hierarchical Clustering for Graph-Structured Data

Authors: Claire Donnat, Susan Holmes

Abstract: Convex clustering is a recent stable alternative to hierarchical clustering. It formulates the recovery of progressively coalescing clusters as a regularized convex problem. While convex clustering was originally designed for handling Euclidean distances between data points, in a growing number of applications, the data is directly characterized by a similarity matrix or weighted graph. In this pa… ▽ More Convex clustering is a recent stable alternative to hierarchical clustering. It formulates the recovery of progressively coalescing clusters as a regularized convex problem. While convex clustering was originally designed for handling Euclidean distances between data points, in a growing number of applications, the data is directly characterized by a similarity matrix or weighted graph. In this paper, we extend the robust hierarchical clustering approach to these broader classes of similarities. Having defined an appropriate convex objective, the crux of this adaptation lies in our ability to provide: (a) an efficient recovery of the regularization path and (b) an empirical demonstration of the use of our method. We address the first challenge through a proximal dual algorithm, for which we characterize both the theoretical efficiency as well as the empirical performance on a set of experiments. Finally, we highlight the potential of our method by showing its application to several real-life datasets, thus providing a natural extension to the current scope of applications of convex clustering. △ Less

Submitted 10 December, 2019; v1 submitted 8 November, 2019; originally announced November 2019.

Uncertainty Quantification in Multivariate Mixed Models for Mass Cytometry Data

Authors: Christof Seiler, Lisa M. Kronstad, Laura J. Simpson, Mathieu Le Gars, Elena Vendrame, Catherine A. Blish, Susan Holmes

Abstract: Mass cytometry technology enables the simultaneous measurement of over 40 proteins on single cells. This has helped immunologists to increase their understanding of heterogeneity, complexity, and lineage relationships of white blood cells. Current statistical methods often collapse the rich single-cell data into summary statistics before proceeding with downstream analysis, discarding the informat… ▽ More Mass cytometry technology enables the simultaneous measurement of over 40 proteins on single cells. This has helped immunologists to increase their understanding of heterogeneity, complexity, and lineage relationships of white blood cells. Current statistical methods often collapse the rich single-cell data into summary statistics before proceeding with downstream analysis, discarding the information in these multivariate datasets. In this article, our aim is to exhibit the use of statistical analyses on the raw, uncompressed data thus improving replicability, and exposing multivariate patterns and their associated uncertainty profiles. We show that multivariate generative models are a valid alternative to univariate hypothesis testing. We propose two models: a multivariate Poisson log-normal mixed model and a logistic linear mixed model. We show that these models are complementary and that either model can account for different confounders. We use Hamiltonian Monte Carlo to provide Bayesian uncertainty quantification. Our models applied to a recent pregnancy study successfully reproduce key findings while quantifying increased overall protein-to-protein correlations between first and third trimester. △ Less

Submitted 19 March, 2019; originally announced March 2019.

Zero-Magnetic Field Fractional Quantum States

Authors: S. Kumar, M. Pepper, S. N. Holmes, H. Montagu, Y. Gul, D. A. Ritchie, I. Farrer

Abstract: Since the discovery of the Fractional Quantum Hall Effect in 1982 there has been considerable theoretical discussion on the possibility of fractional quantization of conductance in the absence of Landau levels formed by a quantizing magnetic field. Although various situations have been theoretically envisaged, particularly lattice models in which band flattening resembles Landau levels, the predic… ▽ More Since the discovery of the Fractional Quantum Hall Effect in 1982 there has been considerable theoretical discussion on the possibility of fractional quantization of conductance in the absence of Landau levels formed by a quantizing magnetic field. Although various situations have been theoretically envisaged, particularly lattice models in which band flattening resembles Landau levels, the predicted fractions have never been observed. In this Letter, we show that odd and even denominator fractions can be observed, and manipulated, in the absence of a quantizing magnetic field, when a low-density electron system in a GaAs based one-dimensional quantum wire is allowed to relax in the second dimension. It is suggested that such a relaxation results in formation of a zig-zag array of electrons with ring paths which establish a cyclic current and a resultant lowering of energy. The behavior has been observed for both symmetric and asymmetric confinement but increasing the asymmetry of the confinement potential, to result in a flattening of confinement, enhances the appearance of new fractional states. We find that an in-plane magnetic field induces new even denominator fractions possibly indicative of electron pairing. The new quantum states described here have implications both for the physics of low dimensional electron systems and also for quantum technologies. This work will enable further development of structures which are designed to electrostatically manipulate the electrons for the formation of particular configurations. In turn, this could result in a designer tailoring of fractional states to amplify particular properties of importance in future quantum computation. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 4 figures, 16 pages including supplementary material

Journal ref: Phys. Rev. Lett. 122, 086803 (2019)

The Block Bootstrap Method for Longitudinal Microbiome Data

Authors: Pratheepa Jeganathan, Benjamin J. Callahan, Diana M. Proctor, David A. Relman, Susan P. Holmes

Abstract: Microbial ecology serves as a foundation for a wide range of scientific and biomedical studies. Rapidly-evolving high-throughput sequencing technology enables the comprehensive search for microbial biomarkers using longitudinal experiments. Such experiments consist of repeated biological observations from each subject over time and are essential in accounting for the high between-subject and withi… ▽ More Microbial ecology serves as a foundation for a wide range of scientific and biomedical studies. Rapidly-evolving high-throughput sequencing technology enables the comprehensive search for microbial biomarkers using longitudinal experiments. Such experiments consist of repeated biological observations from each subject over time and are essential in accounting for the high between-subject and within-subject variability. Unfortunately, many of the statistical tests based on parametric models rely on correctly specifying temporal dependence structure which is unavailable in most microbiome data. In this paper, we propose an extension of the nonparametric bootstrap method that enables inference on these types longitudinal data. The proposed moving block bootstrap (MBB) method accounts for within-subject dependency by using overlapping blocks of repeated observations within each subject to draw valid inferences based on approximately pivotal statistics. Our simulation studies show an increase in power compared to merge-by-subject (MBS) strategies. We also show that compared to tests that presume independent samples (PIS), our proposed method reduces false microbial biomarker discovery rates. In this paper, we illustrated the MBB method using three different pregnancy data and an oral microbiome data. We provide an open-source R package https://github.com/PratheepaJ/bootLong to make our method accessible and the study in this paper reproducible. △ Less

Submitted 29 November, 2018; v1 submitted 6 September, 2018; originally announced September 2018.

Self-Organised Fractional Quantisation in a Hole Quantum Wire

Authors: Y. Gul, S. N. Holmes, M. Myronov, S. Kumar, M. Pepper

Abstract: We have investigated hole transport in one-dimensional quantum wires in strained germanium two-dimensional layers. The ballistic conductance characteristics show the regular quantised plateaux in units of n2e2/h, where n is an integer. Additionally, new quantised levels are formed which correspond to values of n = 1/4 reducing to 1/8 in the presence of a strong parallel magnetic field which lifts… ▽ More We have investigated hole transport in one-dimensional quantum wires in strained germanium two-dimensional layers. The ballistic conductance characteristics show the regular quantised plateaux in units of n2e2/h, where n is an integer. Additionally, new quantised levels are formed which correspond to values of n = 1/4 reducing to 1/8 in the presence of a strong parallel magnetic field which lifts the spin degeneracy but does not quantise the wavefunction. A further plateau is observed corresponding to n = 1/32 which does not change in the presence of a parallel magnetic field. These values indicate that the system is behaving as if charge was fractionalised with values e/2 and e/4, possible mechanisms are discussed. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Comments: 12 pages, 4 figures

Journal ref: J. Phys.: Condens. Matter, 30,09LT01 (2018)

Tracking network dynamics: a survey of distances and similarity metrics

Authors: Claire Donnat, Susan Holmes

Abstract: From longitudinal biomedical studies to social networks, graphs have emerged as a powerful framework for describing evolving interactions between agents in complex systems. In such studies, after pre-processing, the data can be represented by a set of graphs, each representing a system's state at different points in time. The analysis of the system's dynamics depends on the selection of the approp… ▽ More From longitudinal biomedical studies to social networks, graphs have emerged as a powerful framework for describing evolving interactions between agents in complex systems. In such studies, after pre-processing, the data can be represented by a set of graphs, each representing a system's state at different points in time. The analysis of the system's dynamics depends on the selection of the appropriate analytical tools. After characterizing similarities between states, a critical step lies in the choice of a distance between graphs capable of reflecting such similarities. While the literature offers a number of distances that one could a priori choose from, their properties have been little investigated and no guidelines regarding the choice of such a distance have yet been provided. In particular, most graph distances consider that the nodes are exchangeable and do not take into account node identities. Accounting for the alignment of the graphs enables us to enhance these distances' sensitivity to perturbations in the network and detect important changes in graph dynamics. Thus the selection of an adequate metric is a decisive --yet delicate--practical matter. In the spirit of Goldenberg, Zheng and Fienberg's seminal 2009 review, the purpose of this article is to provide an overview of commonly-used graph distances and an explicit characterization of the structural changes that they are best able to capture. We use as a guiding thread to our discussion the application of these distances to the analysis of both a longitudinal microbiome dataset and a brain fMRI study. We show examples of using permutation tests to detect the effect of covariates on the graphs' variability. Synthetic examples provide intuition as to the qualities and drawbacks of the different distances. Above all, we provide some guidance for choosing one distance over another in certain types of applications. △ Less

Submitted 9 March, 2018; v1 submitted 22 January, 2018; originally announced January 2018.

Inference of Dynamic Regimes in the Microbiome

Authors: Kris Sankaran, Susan P. Holmes

Abstract: Many studies have been performed to characterize the dynamics and stability of the microbiome across a range of environmental contexts [Costello et al., 2012, Faust et al., 2015]. For example, it is often of interest to identify time intervals within which certain subsets of taxa have an interesting pattern of behavior. Viewed abstractly, these problems often have a flavor not just of time series… ▽ More Many studies have been performed to characterize the dynamics and stability of the microbiome across a range of environmental contexts [Costello et al., 2012, Faust et al., 2015]. For example, it is often of interest to identify time intervals within which certain subsets of taxa have an interesting pattern of behavior. Viewed abstractly, these problems often have a flavor not just of time series modeling but also of regime detection, a problem with a rich history across a variety of applications, including speech recognition [Fox et al., 2011], finance [Lee, 2009], EEG analysis [Camilleri et al., 2014], and geophysics [Weatherley and Mora, 2002]. In spite of the parallels, regime detection methods are rarely used in microbiome analysis, most likely due to the fact that references for these methods are scattered across several literatures, descriptions are inaccessible outside limited research communities, and implementations are difficult to come across. We distill the core ideas of different regime detection methods, provide example applications, and share reproducible code, making these techniques more accessible to microbiome researchers. We re-analyze data of Dethlefsen and Relman [2011], a study of the effects of antibiotics on the microbiome, using Classification and Regression Trees (CART) [Breiman et al., 1984], Hidden Markov Models (HMMs) [Rabiner and Juang, 1986], Bayesian nonparametric HMMs [Teh and Jordan, 2010, Fox et al., 2008], mixtures of Gaussian Processes (GPs) [Rasmussen and Ghahramani, 2002], switching dynamical systems [Linderman et al., 2016], and multiple changepoint detection [Fan and Mackey, 2015]. Along the way, we summarize each method, their relevance to the microbiome, and tradeoffs associated with using them. Ultimately, our goal is to describe types of temporal or regime switching structure that can be incorporated into studies of microbiome dynamics. △ Less

Submitted 30 November, 2017; originally announced December 2017.

Latent Variable Modeling for the Microbiome

Authors: Kris Sankaran, Susan P. Holmes

Abstract: The human microbiome is a complex ecological system, and describing its structure and function under different environmental conditions is important from both basic scientific and medical perspectives. Viewed through a biostatistical lens, many microbiome analysis goals can be formulated as latent variable modeling problems. However, although probabilistic latent variable models are a cornerstone… ▽ More The human microbiome is a complex ecological system, and describing its structure and function under different environmental conditions is important from both basic scientific and medical perspectives. Viewed through a biostatistical lens, many microbiome analysis goals can be formulated as latent variable modeling problems. However, although probabilistic latent variable models are a cornerstone of modern unsupervised learning, they are rarely applied in the context of microbiome data analysis, in spite of the evolutionary, temporal, and count structure that could be directly incorporated through such models. We explore the application of probabilistic latent variable models to microbiome data, with a focus on Latent Dirichlet Allocation, Nonnegative Matrix Factorization, and Dynamic Unigram models. To develop guidelines for when different methods are appropriate, we perform a simulation study. We further illustrate and compare these techniques using the data of [10], a study on the effects of antibiotics on bacterial community composition. Code and data for all simulations and case studies are available publicly. △ Less

Submitted 15 November, 2017; v1 submitted 15 June, 2017; originally announced June 2017.

Comments: 31 pages, 16 figures

Template shape estimation: correcting an asymptotic bias

Authors: Nina Miolane, Susan Holmes, Xavier Pennec

Abstract: We use tools from geometric statistics to analyze the usual estimation procedure of a template shape. This applies to shapes from landmarks, curves, surfaces, images etc. We demonstrate the asymptotic bias of the template shape estimation using the stratified geometry of the shape space. We give a Taylor expansion of the bias with respect to a parameter $σ$ describing the measurement error on the… ▽ More We use tools from geometric statistics to analyze the usual estimation procedure of a template shape. This applies to shapes from landmarks, curves, surfaces, images etc. We demonstrate the asymptotic bias of the template shape estimation using the stratified geometry of the shape space. We give a Taylor expansion of the bias with respect to a parameter $σ$ describing the measurement error on the data. We propose two bootstrap procedures that quantify the bias and correct it, if needed. They are applicable for any type of shape data. We give a rule of thumb to provide intuition on whether the bias has to be corrected. This exhibits the parameters that control the bias' magnitude. We illustrate our results on simulated and real shape data. △ Less

Submitted 2 February, 2017; v1 submitted 6 September, 2016; originally announced October 2016.

WASP-92b, WASP-93b and WASP-118b: Three new transiting close-in giant planets

Authors: K. L. Hay, A. Collier-Cameron, A. P. Doyle, G. Hébrard, I. Skillen, D. R. Anderson, S. C. C. Barros, D. J. A. Brown, F. Bouchy, R. Busuttil, P. Delorme, L. Delrez, O. Demangeon, R. F. Díaz, M. Gillon, E. Gonzàlez, C. Hellier, S. Holmes, J. F. Jarvis, E. Jehin, Y. C. Joshi, U. Kolb, M. Lendl, P. F. L. Maxted, J. McCormac , et al. (14 additional authors not shown)

Abstract: We present the discovery of three new transiting giant planets, first detected with the WASP telescopes, and establish their planetary nature with follow up spectroscopy and ground-based photometric lightcurves. WASP-92 is an F7 star, with a moderately inflated planet orbiting with a period of 2.17 days, which has $R_p = 1.461 \pm 0.077 R_{\rm J}$ and $M_p = 0.805 \pm 0.068 M_{\rm J}$. WASP-93b or… ▽ More We present the discovery of three new transiting giant planets, first detected with the WASP telescopes, and establish their planetary nature with follow up spectroscopy and ground-based photometric lightcurves. WASP-92 is an F7 star, with a moderately inflated planet orbiting with a period of 2.17 days, which has $R_p = 1.461 \pm 0.077 R_{\rm J}$ and $M_p = 0.805 \pm 0.068 M_{\rm J}$. WASP-93b orbits its F4 host star every 2.73 days and has $R_p = 1.597 \pm 0.077 R_{\rm J}$ and $M_p = 1.47 \pm 0.029 M_{\rm J}$. WASP-118b also has a hot host star (F6) and is moderately inflated, where $R_p = 1.440 \pm 0.036 R_{\rm J}$ and $M_p = 0.513 \pm 0.041 M_{\rm J}$ and the planet has an orbital period of 4.05 days. They are bright targets (V = 13.18, 10.97 and 11.07 respectively) ideal for further characterisation work, particularly WASP-118b, which is being observed by K2 as part of campaign 8. WASP-93b is expected to be tidally migrating outwards, which is divergent from the tidal behaviour of the majority of hot Jupiters discovered. △ Less

Submitted 4 July, 2016; originally announced July 2016.

Comments: 14 pages, 10 figures; submitted to MNRAS

Bayesian Nonparametric Ordination for the Analysis of Microbial Communities

Authors: Boyu Ren, Sergio Bacallado, Stefano Favaro, Susan Holmes, Lorenzo Trippa

Abstract: Human microbiome studies use sequencing technologies to measure the abundance of bacterial species or Operational Taxonomic Units (OTUs) in samples of biological material. Typically the data are organized in contingency tables with OTU counts across heterogeneous biological samples. In the microbial ecology community, ordination methods are frequently used to investigate latent factors or clusters… ▽ More Human microbiome studies use sequencing technologies to measure the abundance of bacterial species or Operational Taxonomic Units (OTUs) in samples of biological material. Typically the data are organized in contingency tables with OTU counts across heterogeneous biological samples. In the microbial ecology community, ordination methods are frequently used to investigate latent factors or clusters that capture and describe variations of OTU counts across biological samples. It remains important to evaluate how uncertainty in estimates of each biological sample's microbial distribution propagates to ordination analyses, including visualization of clusters and projections of biological samples on low dimensional spaces. We propose a Bayesian analysis for dependent distributions to endow frequently used ordinations with estimates of uncertainty. A Bayesian nonparametric prior for dependent normalized random measures is constructed, which is marginally equivalent to the normalized generalized Gamma process, a well-known prior for nonparametric analyses. In our prior the dependence and similarity between microbial distributions is represented by latent factors that concentrate in a low dimensional space. We use a shrinkage prior to tune the dimensionality of the latent factors. The resulting posterior samples of model parameters can be used to evaluate uncertainty in analyses routinely applied in microbiome studies. Specifically, by combining them with multivariate data analysis techniques we can visualize credible regions in ecological ordination plots. The characteristics of the proposed model are illustrated through a simulation study and applications in two microbiome datasets. △ Less

Submitted 20 January, 2017; v1 submitted 19 January, 2016; originally announced January 2016.

An 800-MeV superconducting LINAC to support megawatt proton operations at Fermilab

Authors: Paul Derwent, Steve Holmes, Valeri Lebedev

Abstract: Active discussion on the high energy physics priorities in the US carried out since summer of 2013 resulted in changes in Fermilab plans for future development of the existing accelerator complex. In particular, the scope of Project X was reduced to the support of the Long Base Neutrino Facility (LBNF) at the project first stage. The name of the facility was changed to the PIP-II (Proton Improveme… ▽ More Active discussion on the high energy physics priorities in the US carried out since summer of 2013 resulted in changes in Fermilab plans for future development of the existing accelerator complex. In particular, the scope of Project X was reduced to the support of the Long Base Neutrino Facility (LBNF) at the project first stage. The name of the facility was changed to the PIP-II (Proton Improvement Plan). This new facility is a logical extension of the existing Proton Improvement Plan aimed at doubling average power of the Fermilab's Booster and Main Injector (MI). Its design and required R&D are closely related to the Project X. The paper discusses the goals of this new facility and changes to the Project X linac introduced to support the goals. △ Less

Submitted 5 February, 2015; originally announced February 2015.

Comments: 4 pp

Report number: FERMILAB-CONF-14-305-AD

A Staged Muon Accelerator Facility For Neutrino and Collider Physics

Authors: Jean-Pierre Delahaye, Charles Ankenbrandt, Stephen Brice, Alan David Bross, Dmitri Denisov, Estia Eichten, Stephen Holmes, Ronald Lipton, David Neuffer, Mark Alan Palmer, S. Alex Bogacz, Patrick Huber, Daniel M. Kaplan, Pavel Snopok, Harold G. Kirk, Robert B. Palmer, Robert D. Ryne

Abstract: Muon-based facilities offer unique potential to provide capabilities at both the Intensity Frontier with Neutrino Factories and the Energy Frontier with Muon Colliders. They rely on a novel technology with challenging parameters, for which the feasibility is currently being evaluated by the Muon Accelerator Program (MAP). A realistic scenario for a complementary series of staged facilities with in… ▽ More Muon-based facilities offer unique potential to provide capabilities at both the Intensity Frontier with Neutrino Factories and the Energy Frontier with Muon Colliders. They rely on a novel technology with challenging parameters, for which the feasibility is currently being evaluated by the Muon Accelerator Program (MAP). A realistic scenario for a complementary series of staged facilities with increasing complexity and significant physics potential at each stage has been developed. It takes advantage of and leverages the capabilities already planned for Fermilab, especially the strategy for long-term improvement of the accelerator complex being initiated with the Proton Improvement Plan (PIP-II) and the Long Baseline Neutrino Facility (LBNF). Each stage is designed to provide an R&D platform to validate the technologies required for subsequent stages. The rationale and sequence of the staging process and the critical issues to be addressed at each stage, are presented. △ Less

Submitted 5 February, 2015; originally announced February 2015.

Comments: 5 pp

Report number: FERMILAB-CONF-14-211-APC-E-T

Status and Opportunities at Project X: A Multi-MW Facility for Intensity Frontier Research

Authors: S. D. Holmes, M. Kaducak, R. Kephart, I. Kourbanis, V. Lebedev, S. Mishra, S. Nagaitsev, N. Solyak, R. Tschirhart

Abstract: Project X is a multi-megawatt proton facility being developed to support a world-leading program in Intensity Frontier physics at Fermilab. The facility will support programs in elementary particle and nuclear physics, with the potential for broader applications in materials and energy research. Project X is in the development stage with a R&D program focused on front end and superconducting RF ac… ▽ More Project X is a multi-megawatt proton facility being developed to support a world-leading program in Intensity Frontier physics at Fermilab. The facility will support programs in elementary particle and nuclear physics, with the potential for broader applications in materials and energy research. Project X is in the development stage with a R&D program focused on front end and superconducting RF acceleration technologies, and with design concepts for a staged implementation. This paper will review the status of the Project X conceptual development and the associated R&D programs. △ Less

Submitted 19 September, 2014; originally announced September 2014.

Comments: 3 p. Presented at North American Particle Accelerator Conference (PAC 2013) Conference

Report number: FERMILAB-CONF-13-364-AD-APC-CD-TD

The Project-x Injector Experiment: A Novel High Performance Front-end For A Future High Power Proton Facility At Fermilab

Authors: S. Nagaitsev, S. Holmes, D. Johnson, M. Kaducak, R. Kephart, V. Lebedev, S. Mishra, A. Shemyakin, N. Solyak, R. Stanek, V. Yakovlev, P. Ostroumov, D. Li, P. Singh, M. Pande, S. Malhotra

Abstract: A multi-MW proton facility, Project X, has been proposed and is currently under development at Fermilab. We are carrying out a program of research and development aimed at integrated systems testing of critical components comprising the front end of Project X. This program, known as the Project X Injector Experiment (PXIE), is being undertaken as a key component of the larger Project X R&D program… ▽ More A multi-MW proton facility, Project X, has been proposed and is currently under development at Fermilab. We are carrying out a program of research and development aimed at integrated systems testing of critical components comprising the front end of Project X. This program, known as the Project X Injector Experiment (PXIE), is being undertaken as a key component of the larger Project X R&D program. The successful completion of this program will validate the concept for the Project X front end, thereby minimizing a primary technical risk element within Project X. PXIE is currently under construction at Fermilab and will be completed over the period FY12-17. PXIE will include an H- ion source, a CW 2.1-MeV RFQ and two superconductive RF (SRF) cryomodules providing up to 25 MeV energy gain at an average beam current of 1 mA (upgradable to 2 mA). Successful systems testing will also demonstrate the viability of novel front end technologies that are expected find applications beyond Project X. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: 4 p. Presented at North American Particle Accelerator Conference (PAC13) Conference

Report number: FERMILAB-CONF-13-415-AD-APC-TD

Curvature and Concentration of Hamiltonian Monte Carlo in High Dimensions

Authors: Susan Holmes, Simon Rubinstein-Salzedo, Christof Seiler

Abstract: In this article, we analyze Hamiltonian Monte Carlo (HMC) by placing it in the setting of Riemannian geometry using the Jacobi metric, so that each step corresponds to a geodesic on a suitable Riemannian manifold. We then combine the notion of curvature of a Markov chain due to Joulin and Ollivier with the classical sectional curvature from Riemannian geometry to derive error bounds for HMC in imp… ▽ More In this article, we analyze Hamiltonian Monte Carlo (HMC) by placing it in the setting of Riemannian geometry using the Jacobi metric, so that each step corresponds to a geodesic on a suitable Riemannian manifold. We then combine the notion of curvature of a Markov chain due to Joulin and Ollivier with the classical sectional curvature from Riemannian geometry to derive error bounds for HMC in important cases, where we have positive curvature. These cases include several classical distributions such as multivariate Gaussians, and also distributions arising in the study of Bayesian image registration. The theoretical development suggests the sectional curvature as a new diagnostic tool for convergence for certain Markov chains. △ Less

Submitted 19 May, 2015; v1 submitted 3 July, 2014; originally announced July 2014.

Comments: Comments welcome

Long-range electronic reconstruction to a $d_{xz,yz}$-dominated Fermi surface below the LaAlO$_3$/SrTiO$_3$ interface

Authors: A. P. Petrović, A. Paré, T. R. Paudel, K. Lee, S. Holmes, C. H. W. Barnes, A. David, T. Wu, E. Y. Tsymbal, C. Panagopoulos

Abstract: Low dimensionality, broken symmetry and easily-modulated carrier concentrations provoke novel electronic phase emergence at oxide interfaces. However, the spatial extent of such reconstructions - i.e. the interfacial "depth" - remains unclear. Examining LaAlO$_3$/SrTiO$_3$ heterostructures at previously unexplored carrier densities $n_{2D}\geq6.9\times10^{14}$ cm$^{-2}$, we observe a Shubnikov-de… ▽ More Low dimensionality, broken symmetry and easily-modulated carrier concentrations provoke novel electronic phase emergence at oxide interfaces. However, the spatial extent of such reconstructions - i.e. the interfacial "depth" - remains unclear. Examining LaAlO$_3$/SrTiO$_3$ heterostructures at previously unexplored carrier densities $n_{2D}\geq6.9\times10^{14}$ cm$^{-2}$, we observe a Shubnikov-de Haas effect for small in-plane fields, characteristic of an anisotropic 3D Fermi surface with preferential $d_{xz,yz}$ orbital occupancy extending over at least 100~nm perpendicular to the interface. Quantum oscillations from the 3D Fermi surface of bulk doped SrTiO$_3$ emerge simultaneously at higher $n_{2D}$. We distinguish three areas in doped perovskite heterostructures: narrow ($<20$ nm) 2D interfaces housing superconductivity and/or other emergent phases, electronically isotropic regions far ($>120$ nm) from the interface and new intermediate zones where interfacial proximity renormalises the electronic structure relative to the bulk. △ Less

Submitted 20 June, 2014; originally announced June 2014.

Comments: Supplementary material available at Scientific Reports website

Journal ref: Scientific Reports 4, 5338 (2014)

Application of CellDesigner to the Selection of Anticancer Drug Targets: Test Case using P53

Authors: Raul Isea, Johan Hoebeke, Rafael Mayo, Fernando Alvarez, David S. Holmes

Abstract: Cancer is a disease involving many genes, consequently it has been difficult to design anticancer drugs that are efficacious over a broad range of cancers. The robustness of cellular responses to gene knockout and the need to reduce undesirable side effects also contribute to the problem of effective anti-cancer drug design. To promote the successful selection of drug targets, each potential targe… ▽ More Cancer is a disease involving many genes, consequently it has been difficult to design anticancer drugs that are efficacious over a broad range of cancers. The robustness of cellular responses to gene knockout and the need to reduce undesirable side effects also contribute to the problem of effective anti-cancer drug design. To promote the successful selection of drug targets, each potential target should be subjected to a systems biology scrutiny to locate effective and specific targets while minimizing undesirable side effects. The gene p53 is considered a good candidate for such a target because it has been implicated in 50% of all cancers and is considered to encode a hub protein that is highly connected to other proteins. Using P53 as a test case, this paper explores the capacity of the systems biology tool, CellDesigner, to aid in the selection of anticancer drug targets and to serve as a teaching resource for human resource development. △ Less

Submitted 11 November, 2013; originally announced November 2013.

Comments: 5 pages, 2 figures

Journal ref: In book: Lineas estrategicas de las Tecnologias de la Informacion y las Comunicaciones en Salud. Editorial Cyted (2013), pp: 75-82. ISBN: 978-84-15413-18-9

The Vortex Signature of Discrete Ferromagnetic Dipoles at the LaAlO$_3$/SrTiO$_3$ Interface

Authors: A. P. Petrović, A. Paré, T. R. Paudel, K. Lee, S. Holmes, C. H. W. Barnes, A. David, T. Wu, E. Y. Tsymbal, C. Panagopoulos

Abstract: A hysteretic in-plane magnetoresistance develops below the superconducting transition of LaAlO$_3$/SrTiO$_3$ interfaces for $\left|H_{/\!/}\right|<$ 0.15 T, independently of the carrier density or oxygen annealing. We show that this hysteresis arises from vortex depinning within a thin superconducting layer, in which the vortices are created by discrete ferromagnetic dipoles located solely above t… ▽ More A hysteretic in-plane magnetoresistance develops below the superconducting transition of LaAlO$_3$/SrTiO$_3$ interfaces for $\left|H_{/\!/}\right|<$ 0.15 T, independently of the carrier density or oxygen annealing. We show that this hysteresis arises from vortex depinning within a thin superconducting layer, in which the vortices are created by discrete ferromagnetic dipoles located solely above the layer. We find no evidence for finite-momentum pairing or bulk magnetism and hence conclude that ferromagnetism is strictly confined to the interface, where it competes with superconductivity. △ Less

Submitted 10 November, 2013; originally announced November 2013.

Comments: 10 pages, including Supplementary Material

Journal ref: New J. Phys. 16 103012 (2014)

Discussion of "Geodesic Monte Carlo on Embedded Manifolds"

Authors: Simon Byrne, Mark Girolami, Persi Diaconis, Christof Seiler, Susan Holmes, Ian L. Dryden, John T. Kent, Marcelo Pereyra, Babak Shahbaba, Shiwei Lan, Jeffrey Streets, Daniel Simpson

Abstract: Contributed discussion and rejoinder to "Geodesic Monte Carlo on Embedded Manifolds" (arXiv:1301.6064) Contributed discussion and rejoinder to "Geodesic Monte Carlo on Embedded Manifolds" (arXiv:1301.6064) △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: Discussion of arXiv:1301.6064. To appear in the Scandinavian Journal of Statistics. 18 pages

Journal ref: Scandinavian Journal of Statistics (2014) 41 (1), pages 1-21

Waste Not, Want Not: Why Rarefying Microbiome Data is Inadmissible

Authors: Paul J. McMurdie, Susan Holmes

Abstract: The interpretation of count data originating from the current generation of DNA sequencing platforms requires special attention. In particular, the per-sample library sizes often vary by orders of magnitude from the same sequencing run, and the counts are overdispersed relative to a simple Poisson model These challenges can be addressed using an appropriate mixture model that simultaneously accoun… ▽ More The interpretation of count data originating from the current generation of DNA sequencing platforms requires special attention. In particular, the per-sample library sizes often vary by orders of magnitude from the same sequencing run, and the counts are overdispersed relative to a simple Poisson model These challenges can be addressed using an appropriate mixture model that simultaneously accounts for library size differences and biological variability. This approach is already well-characterized and implemented for RNA-Seq data in R packages such as edgeR and DESeq. We use statistical theory, extensive simulations, and empirical data to show that variance stabilizing normalization using a mixture model like the negative binomial is appropriate for microbiome count data. In simulations detecting differential abundance, normalization procedures based on a Gamma-Poisson mixture model provided systematic improvement in performance over crude proportions or rarefied counts -- both of which led to a high rate of false positives. In simulations evaluating clustering accuracy, we found that the rarefying procedure discarded samples that were nevertheless accurately clustered by alternative methods, and that the choice of minimum library size threshold was critical in some settings, but with an optimum that is unknown in practice. Techniques that use variance stabilizing transformations by modeling microbiome count data with a mixture distribution, such as those implemented in edgeR and DESeq, substantially improved upon techniques that attempt to normalize by rarefying or crude proportions. Based on these results and well-established statistical theory, we advocate that investigators avoid rarefying altogether. We have provided microbiome-specific extensions to these tools in the R package, phyloseq. △ Less

Submitted 12 December, 2013; v1 submitted 1 October, 2013; originally announced October 2013.

Comments: 22 pages, 5 figures, 2 supplementary sections

MSC Class: 62P10

Measures of dependence between random vectors and tests of independence. Literature review

Authors: Julie Josse, Susan Holmes

Abstract: Simple correlation coefficients between two variables have been generalized to measure association between two matrices in many ways. Coefficients such as the RV coefficient, the distance covariance (dCov) coefficient and kernel based coefficients have been adopted by different research communities. Scientists use these coefficients to test whether two random vectors are linked. If they are, it is… ▽ More Simple correlation coefficients between two variables have been generalized to measure association between two matrices in many ways. Coefficients such as the RV coefficient, the distance covariance (dCov) coefficient and kernel based coefficients have been adopted by different research communities. Scientists use these coefficients to test whether two random vectors are linked. If they are, it is important to uncover what patterns exist in these associations. We discuss the topic of measures of dependence between random vectors and tests of independence and show links between different approaches. We document some of the interesting rediscoveries and lack of interconnection between bodies of literature. After providing definitions of the coefficients and associated tests, we present the recent improvements that enhance their statistical properties and ease of interpretation. We summarize multi-table approaches and provide scenarii where the indices can provide useful summaries of heterogeneous multi-block data. We illustrate these different strategies on several examples of real data and suggest directions for future research. △ Less

Submitted 17 August, 2014; v1 submitted 28 July, 2013; originally announced July 2013.

Comments: Incorporated new section on actual examples of data analyses

MSC Class: 62H20

Discovery of WASP-65b and WASP-75b: Two Hot Jupiters Without Highly Inflated Radii

Authors: Y. Gómez Maqueo Chew, F. Faedi, D. Pollacco, D. J. A. Brown, A. P. Doyle, A. Collier Cameron, M. Gillon, M. Lendl, B. Smalley, A. H. M. J. Triaud, R. G. West, P. J. Wheatley, R. Busuttil, C. Liebig, D. R. Anderson, D. J. Armstrong, S. C. C. Barros, J. Bento, J. Bochinski, V. Burwitz, L. Delrez, B. Enoch, A. Fumel, C. A. Haswell, G. Hébrard , et al. (16 additional authors not shown)

Abstract: We report the discovery of two transiting hot Jupiters, WASP-65b (M_pl = 1.55 +/- 0.16 M_J; R_pl = 1.11 +/- 0.06 R_J), and WASP-75b (M_pl = 1.07 +/- 0.05 M_J; R_pl = 1.27 +/- 0.05 R_J). They orbit their host star every 2.311, and 2.484 days, respectively. The planet host WASP-65 is a G6 star (T_eff = 5600 K, [Fe/H] = -0.07 +/- 0.07, age > 8 Gyr); WASP-75 is an F9 star (T_eff = 6100 K, [Fe/H] = 0.0… ▽ More We report the discovery of two transiting hot Jupiters, WASP-65b (M_pl = 1.55 +/- 0.16 M_J; R_pl = 1.11 +/- 0.06 R_J), and WASP-75b (M_pl = 1.07 +/- 0.05 M_J; R_pl = 1.27 +/- 0.05 R_J). They orbit their host star every 2.311, and 2.484 days, respectively. The planet host WASP-65 is a G6 star (T_eff = 5600 K, [Fe/H] = -0.07 +/- 0.07, age > 8 Gyr); WASP-75 is an F9 star (T_eff = 6100 K, [Fe/H] = 0.07 +/- 0.09, age of 3 Gyr). WASP-65b is one of the densest known exoplanets in the mass range 0.1 and 2.0 M_J (rho_pl = 1.13 +/- 0.08 rho_J), a mass range where a large fraction of planets are found to be inflated with respect to theoretical planet models. WASP-65b is one of only a handful of planets with masses of around 1.5 M_J, a mass regime surprisingly underrepresented among the currently known hot Jupiters. The radius of Jupiter-mass WASP-75b is slightly inflated (< 10%) as compared to theoretical planet models with no core, and has a density similar to that of Saturn (rho_pl = 0.52 +/- 0.06 rho_J). △ Less

Submitted 1 October, 2013; v1 submitted 24 July, 2013; originally announced July 2013.

Comments: 11 pages, 6 figures. Version 2 only minor changes; accepted for publication to A&A

Project X: Accelerator Reference Design

Authors: S. D. Holmes, R. Alber, B. Chase, K. Gollwitzer, D. Johnson, M. Kaducak, A. Klebaner, I. Kourbanis, V. Lebedev, A. Leveling, D. Li, S. Nagaitsev, P. Ostroumov, R. Pasquinelli, J. Patrick, L. Prost, V. Scarpine, A. Shemyakin, N. Solyak, J. Steimel, V. Yakovlev, R. Zwaska

Abstract: Part 1 of "Project X: Accelerator Reference Design, Physics Opportunities, Broader Impacts". Part 1 contains the volume Preface and a description of the conceptual design for a high-intensity proton accelerator facility being developed to support a world-leading program of Intensity Frontier physics over the next two decades at Fermilab. Subjects covered include performance goals, the accelerator… ▽ More Part 1 of "Project X: Accelerator Reference Design, Physics Opportunities, Broader Impacts". Part 1 contains the volume Preface and a description of the conceptual design for a high-intensity proton accelerator facility being developed to support a world-leading program of Intensity Frontier physics over the next two decades at Fermilab. Subjects covered include performance goals, the accelerator physics design, and the technological basis for such a facility. Part 2 is available as arXiv:1306.5009 [hep-ex] and Part 3 is available as arXiv:1306.5024 [physics.acc-ph]. △ Less

Submitted 15 July, 2013; v1 submitted 20 June, 2013; originally announced June 2013.

Comments: Edited by S. D. Holmes 212 pp. with many figures; prepared in part for the DPF Community Summer Study; v2 corrects typos and conforms with the version being printed

Report number: FERMILAB-TM-2557; ANL/PHY-13/2; BNL-101116-2013-BC/81834; JLAB-ACP-13-1725; LBNL-6334E; PNNL-22523; UASLP-IF-13-001; SLAC-R-1029

Project X: A Flexible High Power Proton Facility

Authors: Stephen Holmes, Sergei Nagaitsev, Robert Tschirhart

Abstract: Project X is a high intensity proton facility that will support a world-leading Intensity Frontier research program over the next several decades at Fermilab. When compared to other facilities in the planning stages elsewhere in the world Project X is completely unique in its ability to deliver, simultaneously, up to 6 MW of site-wide beam power to multiple experiments, at multiple energies, and w… ▽ More Project X is a high intensity proton facility that will support a world-leading Intensity Frontier research program over the next several decades at Fermilab. When compared to other facilities in the planning stages elsewhere in the world Project X is completely unique in its ability to deliver, simultaneously, up to 6 MW of site-wide beam power to multiple experiments, at multiple energies, and with flexible beam formats. Project X will support a wide range of experiments based on neutrinos, muons, kaons, nucleons, and nuclei. In addition, Project X will lay the foundation for the long-term development of a Neutrino Factory and/or Muon Collider. A complete concept for Project X has been developed and is documented in the Project X Reference Design Report. △ Less

Submitted 20 June, 2013; v1 submitted 10 May, 2013; originally announced May 2013.

Comments: Whitepaper submitted to the Snowmass 2013 Summer Study. Five pages, one figure

The Legacy of the Tevatron in the Area of Accelerator Science

Authors: Stephen D. Holmes, Vladimir D. Shiltsev

Abstract: For more than 25 years the Tevatron was the highest energy accelerator in the world, providing the first access to particle collisions beyond 1 TeV and achieving an ultimate performance a factor of four hundred beyond the original design goals. This article reviews the many formidable challenges that were overcome, and the knowledge gained, in building, operating, and improving the Tevatron over i… ▽ More For more than 25 years the Tevatron was the highest energy accelerator in the world, providing the first access to particle collisions beyond 1 TeV and achieving an ultimate performance a factor of four hundred beyond the original design goals. This article reviews the many formidable challenges that were overcome, and the knowledge gained, in building, operating, and improving the Tevatron over its lifetime. These challenges included: the first operations of an accelerator based on superconducting magnets, production of antiprotons in sufficient numbers to support a useable luminosity, management of beam-beam, intrabeam, and other collective effects, novel manipulations of the beam longitudinal phase space, and development and application of a wide variety of innovative technologies. These achievements established the legacy of the Tevatron as the progenitor of all subsequently constructed high energy hadron colliders. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: Submitted to Annual Reviews of Nuclear and Particle Science

Report number: FERMILAB-PUB-13-003-APC

PXIE: Project X Injector Experiment

Authors: P. N. Ostroumov, S. D. Holmes, R. D. Kephart, J. S. Kerby, V. A. Lebedev, C. S. Mishra, S. Nagaitsev, A. V. Shemyakin, N. Solyak, R. P. Stanek, D. Li

Abstract: A multi-MW proton facility, Project X has been proposed and is currently under development at Fermilab. As part of this development program, we are constructing a prototype of the front end of the Project X linac at Fermilab. The construction and successful operations of this facility will validate the concept for the Project X front end, thereby minimizing the primary technical risk element withi… ▽ More A multi-MW proton facility, Project X has been proposed and is currently under development at Fermilab. As part of this development program, we are constructing a prototype of the front end of the Project X linac at Fermilab. The construction and successful operations of this facility will validate the concept for the Project X front end, thereby minimizing the primary technical risk element within the Project. The Project X Injector Experiment (PXIE) can be constructed over the period FY12-16 and will include an H- ion source, a CW 2.1-MeV RFQ and two SC cryomodules providing up to 30 MeV energy gain at an average beam current of 1 mA. Successful operations of the facility will demonstrate the viability of novel front end technologies that will find applications beyond Project X in the longer term. △ Less

Submitted 30 January, 2013; originally announced January 2013.

Comments: 3 pp. 3rd International Particle Accelerator Conference (IPAC 2012) 20-25 May 2012, New Orleans, Louisiana

Report number: FERMILAB-CONF-12-172-AD

Near-UV Absorption, Chromospheric Activity, and Star-Planet Interactions in the WASP-12 system

Authors: C. A. Haswell, L. Fossati, T. Ayres, K. France, C. S. Froning, S. Holmes, U. C. Kolb, R. Busuttil, R. A. Street, L. Hebb, A. Collier Cameron, B. Enoch, V. Burwitz, J. Rodriguez, R. G. West, D. Pollacco, P. J. Wheatley, A. Carter

Abstract: We observed the extreme close-in hot Jupiter system WASP-12 with HST. Near-UV transits up to three times deeper than the optical transit of WASP-12b reveal extensive diffuse gas, extending well beyond the Roche lobe. The distribution of absorbing gas varies between visits. The deepest NUV transits are at wavelength ranges with strong photospheric absorption, implying the absorbing gas may have tem… ▽ More We observed the extreme close-in hot Jupiter system WASP-12 with HST. Near-UV transits up to three times deeper than the optical transit of WASP-12b reveal extensive diffuse gas, extending well beyond the Roche lobe. The distribution of absorbing gas varies between visits. The deepest NUV transits are at wavelength ranges with strong photospheric absorption, implying the absorbing gas may have temperature and composition similar to the stellar photosphere. Our spectra reveal significantly enhanced absorption (greater than 3 σbelow the median) at ~200 wavelengths on each of two HST visits; 65 of these wavelengths are consistent between the two visits, using a strict criterion for velocity matching which excludes matches with velocity shifts exceeding ~20 km/s. Excess transit depths are robustly detected throughout the inner wings of the MgII resonance lines independently on both HST visits. We detected absorption in FeII 2586A, the heaviest species yet detected in an exoplanet transit. The MgII line cores have zero flux, emission cores exhibited by every other observed star of similar age and spectral type are conspicuously absent. WASP-12 probably produces normal MgII profiles, but the inner portions of these strong resonance lines are likely affected by extrinsic absorption. The required Mg+ column is an order of magnitude greater than expected from the ISM, though we cannot completely dismiss that possibility. A more plausible source of absorption is gas lost by WASP-12b. We show that planetary mass loss can produce the required column. Our Visit 2 NUV light curves show evidence for a stellar flare. We show that some of the possible transit detections in resonance lines of rare elements may be due instead to non-resonant transitions in common species. We present optical observations and update the transit ephemeris. △ Less

Submitted 25 January, 2013; v1 submitted 9 January, 2013; originally announced January 2013.

Comments: 79 pages, 23 figures, published in ApJ "replaced" version is identical except for corrected punctuation in the arXiv abstract

Journal ref: 2012, The Astrophysical Journal, 760, 79

Solitary Matter Waves in Combined Linear and Nonlinear Potentials: Detection, Stability, and Dynamics

Authors: Scott Holmes, Mason A. Porter, Peter Krüger, Panayotis G. Kevrekidis

Abstract: We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and consider, in particular, how to reveal it by exploiting the… ▽ More We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and consider, in particular, how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov-de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step. △ Less

Submitted 24 September, 2013; v1 submitted 8 January, 2013; originally announced January 2013.

Comments: 7 pages (published version), 4 figures

Journal ref: Phys. Rev. A, Vol. 88, article 033627, 2013

WASP-54b, WASP-56b and WASP-57b: Three new sub-Jupiter mass planets from SuperWASP

Authors: F. Faedi, D. Pollacco, S. C. C. Barros, D. Brown, A. Collier Cameron, A. P. Doyle, R. Enoch, M. Gillon, Y. Gomez Maqueo Chew, G. Hebrard, M. Lendl, C. Liebig, B. Smalley, A. H. M. J. Triaud, R. G. West, P. J. Wheatley, K. A. Alsubai, D. R. Anderson, D. J. Armstrong, J. Bento, J. Bochinski, F. Bouchy, R. Busuttil, L. Fossati, A. Fumel , et al. (16 additional authors not shown)

Abstract: We present three newly discovered sub-Jupiter mass planets from the SuperWASP survey: WASP-54b is a heavily bloated planet of mass 0.636$^{+0.025}_{-0.024}$ \mj and radius 1.653$^{+0.090}_{-0.083}$ \rj. It orbits a F9 star, evolving off the main sequence, every 3.69 days. Our MCMC fit of the system yields a slightly eccentric orbit ($e=0.067^{+0.033}_{-0.025}$) for WASP-54b. We investigated furthe… ▽ More We present three newly discovered sub-Jupiter mass planets from the SuperWASP survey: WASP-54b is a heavily bloated planet of mass 0.636$^{+0.025}_{-0.024}$ \mj and radius 1.653$^{+0.090}_{-0.083}$ \rj. It orbits a F9 star, evolving off the main sequence, every 3.69 days. Our MCMC fit of the system yields a slightly eccentric orbit ($e=0.067^{+0.033}_{-0.025}$) for WASP-54b. We investigated further the veracity of our detection of the eccentric orbit for WASP-54b, and we find that it could be real. However, given the brightness of WASP-54 V=10.42 magnitudes, we encourage observations of a secondary eclipse to draw robust conclusions on both the orbital eccentricity and the thermal structure of the planet. WASP-56b and WASP-57b have masses of 0.571$^{+0.034}_{-0.035}$ \mj and $0.672^{+0.049}_{-0.046}$ \mj, respectively; and radii of $1.092^{+0.035}_{-0.033}$ \rj for WASP-56b and $0.916^{+0.017}_{-0.014}$ \rj for WASP-57b. They orbit main sequence stars of spectral type G6 every 4.67 and 2.84 days, respectively. WASP-56b and WASP-57b show no radius anomaly and a high density possibly implying a large core of heavy elements; possibly as high as $\sim$50 M$_{\oplus}$ in the case of WASP-57b. However, the composition of the deep interior of exoplanets remain still undetermined. Thus, more exoplanet discoveries such as the ones presented in this paper, are needed to understand and constrain giant planets' physical properties. △ Less

Submitted 16 January, 2013; v1 submitted 8 October, 2012; originally announced October 2012.

Comments: 15pages, 19 figures, Submitted to A&A 8-10-2012 Accepted for publication 23-12-2012

New Measurements of the Transverse Beam Asymmetry for Elastic Electron Scattering from Selected Nuclei

Authors: The HAPPEX, PREX Collaborations, :, S. Abrahamyan, A. Acha, A. Afanasev, Z. Ahmed, H. Albataineh, K. Aniol, D. S. Armstrong, W. Armstrong, J. Arrington, T. Averett, B. Babineau, S. L. Bailey, J. Barber, A. Barbieri, A. Beck, V. Bellini, R. Beminiwattha, H. Benaoum, J. Benesch, F. Benmokhtar, P. Bertin, T. Bielarski , et al. (173 additional authors not shown)

Abstract: We have measured the beam-normal single-spin asymmetry $A_n$ in the elastic scattering of 1-3 GeV transversely polarized electrons from $^1$H and for the first time from $^4$He, $^{12}$C, and $^{208}$Pb. For $^1$H, $^4$He and $^{12}$C, the measurements are in agreement with calculations that relate $A_n$ to the imaginary part of the two-photon exchange amplitude including inelastic intermediate st… ▽ More We have measured the beam-normal single-spin asymmetry $A_n$ in the elastic scattering of 1-3 GeV transversely polarized electrons from $^1$H and for the first time from $^4$He, $^{12}$C, and $^{208}$Pb. For $^1$H, $^4$He and $^{12}$C, the measurements are in agreement with calculations that relate $A_n$ to the imaginary part of the two-photon exchange amplitude including inelastic intermediate states. Surprisingly, the $^{208}$Pb result is significantly smaller than the corresponding prediction using the same formalism. These results suggest that a systematic set of new $A_n$ measurements might emerge as a new and sensitive probe of the structure of heavy nuclei. △ Less

Submitted 12 October, 2012; v1 submitted 30 August, 2012; originally announced August 2012.

Comments: 5 pages, 3 figures, accepted by PRL. v3: fixed one author name and affiliation, otherwise no change

Journal ref: Phys. Rev. Lett. 109, (2012) 192501

Sampling From A Manifold

Authors: Persi Diaconis, Susan Holmes, Mehrdad Shahshahani

Abstract: We develop algorithms for sampling from a probability distribution on a submanifold embedded in Rn. Applications are given to the evaluation of algorithms in 'Topological Statistics'; to goodness of fit tests in exponential families and to Neyman's smooth test. This article is partially expository, giving an introduction to the tools of geometric measure theory. We develop algorithms for sampling from a probability distribution on a submanifold embedded in Rn. Applications are given to the evaluation of algorithms in 'Topological Statistics'; to goodness of fit tests in exponential families and to Neyman's smooth test. This article is partially expository, giving an introduction to the tools of geometric measure theory. △ Less

Submitted 4 July, 2012; v1 submitted 28 June, 2012; originally announced June 2012.

Project X functional requirements specification

Authors: The Project X Collaboration, S. D. Holmes, S. D. Henderson, R. Kephart, J. Kerby, I. Kourbanis, V. Lebedev, S. Mishra, S. Nagaitsev, N. Solyak, R. Tschirhart

Abstract: Project X is a multi-megawatt proton facility being developed to support a world-leading program in Intensity Frontier physics at Fermilab. The facility is designed to support programs in elementary particle and nuclear physics, with possible applications to nuclear energy research. A Functional Requirements Specification has been developed in order to establish performance criteria for the Projec… ▽ More Project X is a multi-megawatt proton facility being developed to support a world-leading program in Intensity Frontier physics at Fermilab. The facility is designed to support programs in elementary particle and nuclear physics, with possible applications to nuclear energy research. A Functional Requirements Specification has been developed in order to establish performance criteria for the Project X complex in support of these multiple missions, and to assure that the facility is designed with sufficient upgrade capability to provide U.S. leadership for many decades to come. This paper will briefly review the previously described Functional Requirements, and then discuss their recent evolution. △ Less

Submitted 14 May, 2012; originally announced May 2012.

Comments: 3 pp

Report number: FERMILAB-CONF-12-079-APC

The duality diagram in data analysis: Examples of modern applications

Authors: Omar De la Cruz, Susan Holmes

Abstract: Today's data-heavy research environment requires the integration of different sources of information into structured data sets that can not be analyzed as simple matrices. We introduce an old technique, known in the European data analyses circles as the Duality Diagram Approach, put to new uses through the use of a variety of metrics and ways of combining different diagrams together. This issue of… ▽ More Today's data-heavy research environment requires the integration of different sources of information into structured data sets that can not be analyzed as simple matrices. We introduce an old technique, known in the European data analyses circles as the Duality Diagram Approach, put to new uses through the use of a variety of metrics and ways of combining different diagrams together. This issue of the Annals of Applied Statistics contains contemporary examples of how this approach provides solutions to hard problems in data integration. We present here the genesis of the technique and how it can be seen as a precursor of the modern kernel based approaches. △ Less

Submitted 24 February, 2012; originally announced February 2012.

Comments: Published in at http://dx.doi.org/10.1214/10-AOAS408 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS408

Journal ref: Annals of Applied Statistics 2011, Vol. 5, No. 4, 2266-2277

Project X functional requirements specification

Authors: S. D. Holmes, S. D. Henderson, R. Kephart, J. Kerby, S. Mishra, S. Nagaitsev, R. Tschirhart

Abstract: Project X is a multi-megawatt proton facility being developed to support intensity frontier research in elementary particle physics, with possible applications to nuclear physics and nuclear energy research, at Fermilab. A Functional Requirements Specification has been developed in order to establish performance criteria for the Project X complex in support of these multiple missions. This paper w… ▽ More Project X is a multi-megawatt proton facility being developed to support intensity frontier research in elementary particle physics, with possible applications to nuclear physics and nuclear energy research, at Fermilab. A Functional Requirements Specification has been developed in order to establish performance criteria for the Project X complex in support of these multiple missions. This paper will describe the Functional Requirements for the Project X facility and the rationale for these requirements. △ Less

Submitted 17 February, 2012; originally announced February 2012.

Comments: 3 pp. Particle Accelerator, 24th Conference (PAC'11) 28 Mar - 1 Apr 2011: New York, USA

Report number: FERMILAB-CONF-11-050-APC

Muon Collider

Authors: Stephen D. Holmes, Vladimir D. Shiltsev

Abstract: Both e+e- and μ+μ- colliders have been proposed as possible candidates for a lepton collider to complement and extend the reach of the Large Hadron Collider (LHC) at CERN. The physics program that could be pursued by a new lepton collider (e+e- or μ+μ-) with sufficient luminosity would include understanding the mechanism behind mass generation and electroweak symmetry breaking; searching for, and… ▽ More Both e+e- and μ+μ- colliders have been proposed as possible candidates for a lepton collider to complement and extend the reach of the Large Hadron Collider (LHC) at CERN. The physics program that could be pursued by a new lepton collider (e+e- or μ+μ-) with sufficient luminosity would include understanding the mechanism behind mass generation and electroweak symmetry breaking; searching for, and possibly discovering, supersymmetric particles; and hunting for signs of extra spacetime dimensions and quantum gravity. However, the appropriate energy reach for such a collider is currently unknown, and will only be determined following initial physics results at the LHC. It is entirely possible that such results will indicate that a lepton collider with a collision energy well in excess of 1 TeV will be required to illuminate the physics uncovered at LHC. Such a requirement would require consideration of muons as the lepton of choice for such a collider. △ Less

Submitted 31 July, 2017; v1 submitted 16 February, 2012; originally announced February 2012.

Comments: v.2., 6 pp. To appear in the 2nd edition of the book Elementary Particles, Landolt-Boernstein Series published by Springer. arXiv admin note: text overlap with arXiv:physics/9901022 by other author

Report number: FERMILAB-PUB-10-539-APC

Measurement of the Neutron Radius of 208Pb Through Parity-Violation in Electron Scattering

Authors: S. Abrahamyan, Z. Ahmed, H. Albataineh, K. Aniol, D. S. Armstrong, W. Armstrong, T. Averett, B. Babineau, A. Barbieri, V. Bellini, R. Beminiwattha, J. Benesch, F. Benmokhtar, T. Bielarski, W. Boeglin, A. Camsonne, M. Canan, P. Carter, G. D. Cates, C. Chen, J. -P. Chen, O. Hen, F. Cusanno, M. M. Dalton, R. De Leo , et al. (110 additional authors not shown)

Abstract: We report the first measurement of the parity-violating asymmetry A_PV in the elastic scattering of polarized electrons from 208Pb. A_PV is sensitive to the radius of the neutron distribution (Rn). The result A_PV = 0.656 \pm 0.060 (stat) \pm 0.014 (syst) ppm corresponds to a difference between the radii of the neutron and proton distributions Rn - Rp = 0.33 +0.16 -0.18 fm and provides the first e… ▽ More We report the first measurement of the parity-violating asymmetry A_PV in the elastic scattering of polarized electrons from 208Pb. A_PV is sensitive to the radius of the neutron distribution (Rn). The result A_PV = 0.656 \pm 0.060 (stat) \pm 0.014 (syst) ppm corresponds to a difference between the radii of the neutron and proton distributions Rn - Rp = 0.33 +0.16 -0.18 fm and provides the first electroweak observation of the neutron skin which is expected in a heavy, neutron-rich nucleus. △ Less

Submitted 13 January, 2012; v1 submitted 12 January, 2012; originally announced January 2012.

Comments: 6 pages, 1 figure

Report number: JLAB-PHY-12-1480

Journal ref: Phys. Rev. Lett. 108, 112502 (2012)