Showing 1–50 of 52 results for author: Woods, C

Iof-maint -- Modular maintenance ontology

Authors: Melinda Hodkiewicz, Caitlin Woods, Matt Selway, Markus Stumptner

Abstract: In this paper we present a publicly-available maintenance ontology (Iof-maint). Iof-maint is a modular ontology aligned with the Industrial Ontology Foundry Core (IOF Core) and contains 20 classes and 2 relations. It provides a set of maintenance-specific terms used in a wide variety of practical data-driven use cases. Iof-maint supports OWL DL reasoning, is documented, and is actively maintained… ▽ More In this paper we present a publicly-available maintenance ontology (Iof-maint). Iof-maint is a modular ontology aligned with the Industrial Ontology Foundry Core (IOF Core) and contains 20 classes and 2 relations. It provides a set of maintenance-specific terms used in a wide variety of practical data-driven use cases. Iof-maint supports OWL DL reasoning, is documented, and is actively maintained on GitHub. In this paper, we describe the evolution of the Iof-maint reference ontology based on the extraction of common concepts identified in a number of application ontologies working with industry maintenance work order, procedure and failure mode data. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Multi-objective optimisation using expected quantile improvement for decision making in disease outbreaks

Authors: Daria Semochkina, Alexander I. J. Forrester, David C Woods

Abstract: Optimization under uncertainty is important in many applications, particularly to inform policy and decision making in areas such as public health. A key source of uncertainty arises from the incorporation of environmental variables as inputs into computational models or simulators. Such variables represent uncontrollable features of the optimization problem and reliable decision making must accou… ▽ More Optimization under uncertainty is important in many applications, particularly to inform policy and decision making in areas such as public health. A key source of uncertainty arises from the incorporation of environmental variables as inputs into computational models or simulators. Such variables represent uncontrollable features of the optimization problem and reliable decision making must account for the uncertainty they propagate to the simulator outputs. Often, multiple, competing objectives are defined from these outputs such that the final optimal decision is a compromise between different goals. Here, we present emulation-based optimization methodology for such problems that extends expected quantile improvement (EQI) to address multi-objective optimization. Focusing on the practically important case of two objectives, we use a sequential design strategy to identify the Pareto front of optimal solutions. Uncertainty from the environmental variables is integrated out using Monte Carlo samples from the simulator. Interrogation of the expected output from the simulator is facilitated by use of (Gaussian process) emulators. The methodology is demonstrated on an optimization problem from public health involving the dispersion of anthrax spores across a spatial terrain. Environmental variables include meteorological features that impact the dispersion, and the methodology identifies the Pareto front even when there is considerable input uncertainty. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Interpreting Sequence-Levenshtein distance for determining error type and frequency between two embedded sequences of equal length

Authors: Robert Logan, Amy W. Wehe, Dori C. Woods, Jon Tilly, Konstantin Khrapko

Abstract: Levenshtein distance is a commonly used edit distance metric, typically applied in language processing, and to a lesser extent, in molecular biology analysis. Biological nucleic acid sequences are often embedded in longer sequences and are subject to insertion and deletion errors that introduce frameshift during sequencing. These frameshift errors are due to string context and should not be counte… ▽ More Levenshtein distance is a commonly used edit distance metric, typically applied in language processing, and to a lesser extent, in molecular biology analysis. Biological nucleic acid sequences are often embedded in longer sequences and are subject to insertion and deletion errors that introduce frameshift during sequencing. These frameshift errors are due to string context and should not be counted as true biological errors. Sequence-Levenshtein distance is a modification to Levenshtein distance that is permissive of frameshift error without additional penalty. However, in a biological context Levenshtein distance needs to accommodate both frameshift and weighted errors, which Sequence-Levenshtein distance cannot do. Errors are weighted when they are associated with a numerical cost that corresponds to their frequency of appearance. Here, we describe a modification that allows the use of Levenshtein distance and Sequence-Levenshtein distance to appropriately accommodate penalty-free frameshift between embedded sequences and correctly weight specific error types. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 10 pages, 8 figures

Bayes Linear Analysis for Statistical Modelling with Uncertain Inputs

Authors: Samuel E. Jackson, David C. Woods

Abstract: Statistical models typically capture uncertainties in our knowledge of the corresponding real-world processes, however, it is less common for this uncertainty specification to capture uncertainty surrounding the values of the inputs to the model, which are often assumed known. We develop general modelling methodology with uncertain inputs in the context of the Bayes linear paradigm, which involves… ▽ More Statistical models typically capture uncertainties in our knowledge of the corresponding real-world processes, however, it is less common for this uncertainty specification to capture uncertainty surrounding the values of the inputs to the model, which are often assumed known. We develop general modelling methodology with uncertain inputs in the context of the Bayes linear paradigm, which involves adjustment of second-order belief specifications over all quantities of interest only, without the requirement for probabilistic specifications. In particular, we propose an extension of commonly-employed second-order modelling assumptions to the case of uncertain inputs, with explicit implementation in the context of regression analysis, stochastic process modelling, and statistical emulation. We apply the methodology to a regression model for extracting aluminium by electrolysis, and emulation of the motivating epidemiological simulator chain to model the impact of an airborne infectious disease. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 37 pages, 2 figures

Semiconductor waveguide quantum well lasers

Authors: Jonathan R. C. Woods, Jon Gorecki, Roman Bek, Stephen C. Richardson, Jake Daykin, Grace Hooper, Emelia Branagan-Harris, Anne C. Tropper, James S. Wilkinson, Michael Jetter, Vasilis Apostolopoulos

Abstract: Coherent laser arrays compatible with silicon photonics are demonstrated in a waveguide geometry in epitaxially grown semiconductor membrane quantum well lasers transferred on substrates of silicon carbide and oxidised silicon; we record lasing thresholds as low as 60 mW of pump power. We study the emission of single lasers and arrays of lasers in the sub-mm range. We are able to create waveguide… ▽ More Coherent laser arrays compatible with silicon photonics are demonstrated in a waveguide geometry in epitaxially grown semiconductor membrane quantum well lasers transferred on substrates of silicon carbide and oxidised silicon; we record lasing thresholds as low as 60 mW of pump power. We study the emission of single lasers and arrays of lasers in the sub-mm range. We are able to create waveguide laser arrays with modal widths of approximately 5-10 μm separated by 10-20 μm, using real and reciprocal space imaging we study their emission characteristics and find that they maintain their mutual coherence while operating on either single or multiple longitudinal modes per lasing cavity. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 12 pages, 7 figures

Optimal designs for experiments for scalar-on-function linear models

Authors: Damianos Michaelides, Maria Adamou, David C. Woods, Antony M. Overstall

Abstract: The aim of this work is to extend the usual optimal experimental design paradigm to experiments where the settings of one or more factors are functions. Such factors are known as profile factors, or as dynamic factors. For these new experiments, a design consists of combinations of functions for each run of the experiment. After briefly introducing the class of profile factors, basis functions are… ▽ More The aim of this work is to extend the usual optimal experimental design paradigm to experiments where the settings of one or more factors are functions. Such factors are known as profile factors, or as dynamic factors. For these new experiments, a design consists of combinations of functions for each run of the experiment. After briefly introducing the class of profile factors, basis functions are described with primary focus given on the B-spline basis system, due to its computational efficiency and useful properties. Basis function expansions are applied to a functional linear model consisting of profile factors, reducing the problem to an optimisation of basis coefficients. The methodology developed comprises special cases, including combinations of profile and non-functional factors, interactions, and polynomial effects. The method is finally applied to an experimental design problem in a Biopharmaceutical study that is performed using the Ambr250 modular bioreactor. △ Less

Submitted 17 November, 2023; v1 submitted 18 October, 2021; originally announced October 2021.

Van der Waals interaction affects wrinkle formation in two-dimensional materials

Authors: Pablo Ares, Yi Bo Wang, Colin R. Woods, James Dougherty, Laura Fumagalli, Francisco Guinea, Benny Davidovitch, Kostya S. Novoselov

Abstract: Nonlinear mechanics of solids is an exciting field that encompasses both beautiful mathematics, such as the emergence of instabilities and the formation of complex patterns, as well as multiple applications. Two-dimensional crystals and van der Waals (vdW) heterostructures allow revisiting this field on the atomic level, allowing much finer control over the parameters and offering atomistic interp… ▽ More Nonlinear mechanics of solids is an exciting field that encompasses both beautiful mathematics, such as the emergence of instabilities and the formation of complex patterns, as well as multiple applications. Two-dimensional crystals and van der Waals (vdW) heterostructures allow revisiting this field on the atomic level, allowing much finer control over the parameters and offering atomistic interpretation of experimental observations. In this work, we consider the formation of instabilities consisting of radially-oriented wrinkles around mono- and few-layer "bubbles" in two-dimensional vdW heterostructures. Interestingly, the shape and wavelength of the wrinkles depend not only on the thickness of the two-dimensional crystal forming the bubble, but also on the atomistic structure of the interface between the bubble and the substrate, which can be controlled by their relative orientation. We argue that the periodic nature of these patterns emanates from an energetic balance between the resistance of the top membrane to bending, which favors large wavelength of wrinkles, and the membrane-substrate vdW attraction, which favors small wrinkle amplitude. Employing the classical "Winkler foundation" model of elasticity theory, we show that the number of radial wrinkles conveys a valuable relationship between the bending rigidity of the top membrane and the strength of the vdW interaction. Armed with this relationship, we use our data to demonstrate a nontrivial dependence of the bending rigidity on the number of layers in the top membrane, which shows two different regimes driven by slippage between the layers, and a high sensitivity of the vdW force to the alignment between the substrate and the membrane. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Journal ref: Proceedings of the National Academy of Sciences 2021, 118, e2025870118

The alignment of interstellar dust grains: thermal flipping and the Davis-Greenstein mechanism

Authors: Joseph C. Weingartner, Erald Kolasi, Cameron Woods

Abstract: Interstellar dust grains are non-spherical and, in some environments, partially aligned along the direction of the interstellar magnetic field. Numerous alignment theories have been proposed, all of which examine the grain rotational dynamics. In 1999, Lazarian & Draine introduced the important concept of thermal flipping, in which internal relaxation processes induce the grain body to flip while… ▽ More Interstellar dust grains are non-spherical and, in some environments, partially aligned along the direction of the interstellar magnetic field. Numerous alignment theories have been proposed, all of which examine the grain rotational dynamics. In 1999, Lazarian & Draine introduced the important concept of thermal flipping, in which internal relaxation processes induce the grain body to flip while its angular momentum remains fixed. Through detailed numerical simulations, we study the role of thermal flipping on the grain dynamics during periods of relatively slow rotation, known as `crossovers', for the special case of a spheroidal grain with a non-uniform mass distribution. Lazarian & Draine proposed that rapid flipping during a crossover would lead to `thermal trapping', in which a systematic torque, fixed relative to the grain body, would time average to zero, delaying spin-up to larger rotational speeds. We find that the time-averaged systematic torque is not zero during the crossover and that thermal trapping is not prevalent. As an application, we examine whether the classic Davis-Greenstein alignment mechanism is viable, for grains residing in the cold neutral medium and lacking superparamagnetic inclusions. We find that Davis-Greenstein alignment is not hindered by thermal trapping, but argue that it is, nevertheless, too inefficient to yield the alignment of large grains responsible for optical and infrared starlight polarization. Davis-Greenstein alignment of small grains could potentially contribute to the observed ultraviolet polarization. The theoretical and computational tools developed here can also be applied to analyses of alignment via radiative torques and rotational disruption of grains. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Comments: 21 pages, 13 figures, accepted for publication in MNRAS

Optical Magnetic Multipolar Resonances in Large Dynamic Metamolecules

Authors: Omar Ibrahim, Sunghee Lee, Sung Wook Kim, Seung Beom Pyun, Connor Woods, Eun Chul Cho, So-Jung Park, Zahra Fakhraai

Abstract: Dynamic metamolecules (DMMs) are composed of a dielectric core made of hydrogel surrounded by randomly-packed plasmonic beads that can display magnetic resonances when excited by light at optical frequencies. Their optical properties can be controlled by controlling their core diameter through temperature variations. We have recently shown that DMMs display strong optical magnetism, including magn… ▽ More Dynamic metamolecules (DMMs) are composed of a dielectric core made of hydrogel surrounded by randomly-packed plasmonic beads that can display magnetic resonances when excited by light at optical frequencies. Their optical properties can be controlled by controlling their core diameter through temperature variations. We have recently shown that DMMs display strong optical magnetism, including magnetic dipole and magnetic quadrupole resonances, offering significant potential for novel applications. Here, we use a T-matrix approach to characterize the magnetic multipole resonance modes of model metamolecules and explore their presence in experimental data. We show that high-order multipole resonances become prominent as the bead size and the overall structure sizes are increased, and when the the inter-bead gap is decreased. In this limit, mode mixing among high-order magnetic multipole modes also become significant, particularly in the directional scattering spectra. We discuss trends in magnetic scattering observed in both experiments and simulations, and provide suggestions for experimental design and verification of high-order optical magnetic resonances in the forward or backward scattering spectra. In addition, angular scattering of higher-order magnetic modes can display Fano-like interference patterns that should be experimentally detectable. △ Less

Submitted 24 May, 2021; v1 submitted 7 March, 2021; originally announced March 2021.

Comments: Revised version of paper: added references, rearranged content and improved discussions, added new figures, and improved existing figures

Statistical methods used to combine the effective reproduction number, R(t), and other related measures of COVID-19 in the UK

Authors: T. Maishman, S. Schaap, D. S. Silk, S. J. Nevitt, D. C. Woods, V. E. Bowman

Abstract: In the recent COVID-19 pandemic, a wide range of epidemiological modelling approaches have been used to predict the effective reproduction number, R(t), and other COVID-19 related measures such as the daily rate of exponential growth, r(t). These candidate models use different modelling approaches or differing assumptions about spatial or age mixing, and some capture genuine uncertainty in scienti… ▽ More In the recent COVID-19 pandemic, a wide range of epidemiological modelling approaches have been used to predict the effective reproduction number, R(t), and other COVID-19 related measures such as the daily rate of exponential growth, r(t). These candidate models use different modelling approaches or differing assumptions about spatial or age mixing, and some capture genuine uncertainty in scientific understanding of disease dynamics. Combining estimates using appropriate statistical methodology from multiple candidate models is important to better understand the variation of these outcome measures to help inform decision making. In this paper, we combine these estimates for specific UK nations and regions using random effects meta analyses techniques, utilising the restricted maximum likelihood (REML) method to estimate the heterogeneity variance parameter, and two approaches to calculate the confidence interval for the combined estimate: the standard Wald-type intervals; and the Knapp and Hartung (KNHA) method. As estimates in this setting are derived using model predictions, each with varying degrees of uncertainty, equal weighting is favoured over the more standard inverse-variance weighting in order avoid potential up-weighting of models providing estimates with lower levels of uncertainty that are not fully accounting for inherent uncertainties. Utilising these meta-analysis techniques has allowed for statistically robust combined estimates to be calculated for key COVID-19 outcome measures. This in turn allows timely and informed decision making based on all of the available information. △ Less

Submitted 12 August, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

MSC Class: 62P10

Charge-polarized interfacial superlattices in marginally twisted hexagonal boron nitride

Authors: C. R. Woods, P. Ares, H. Nevison-Andrews, M. J. Holwill, R. Fabregas, F. Guinea, A. K. Geim, K. S. Novoselov, N. R. Walet, L. Fumagalli

Abstract: When two-dimensional crystals are brought into close proximity, their interaction results in strong reconstruction of electronic spectrum and local crystal structure. Such reconstruction strongly depends on the twist angle between the two crystals and has received growing attention due to new interesting electronic and optical properties that arise in graphene and transitional metal dichalcogenide… ▽ More When two-dimensional crystals are brought into close proximity, their interaction results in strong reconstruction of electronic spectrum and local crystal structure. Such reconstruction strongly depends on the twist angle between the two crystals and has received growing attention due to new interesting electronic and optical properties that arise in graphene and transitional metal dichalcogenides. Similarly, novel and potentially useful properties are expected to appear in insulating crystals. Here we study two insulating crystals of hexagonal boron nitride (hBN) stacked at a small twist angle. Using electrostatic force microscopy, we observe ferroelectric-like domains arranged in triangular superlattices with a large surface potential that is independent on the size and orientation of the domains as well as the thickness of the twisted hBN crystals. The observation is attributed to interfacial elastic deformations that result in domains with a large density of out-of-plane polarized dipoles formed by pairs of boron and nitrogen atoms belonging to the opposite interfacial surfaces. This effectively creates a bilayer-thick ferroelectric with oppositely polarized (BN and NB) dipoles in neighbouring domains, in agreement with our modelling. The demonstrated electrostatic domains and their superlattices offer many new possibilities in designing novel van der Waals heterostructures. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Adaptive multi-channel event segmentation and feature extraction for monitoring health outcomes

Authors: Xichen She, Yaya Zhai, Ricardo Henao, Christopher W. Woods, Christopher Chiu, Geoffrey S. Ginsburg, Peter X. K. Song, Alfred O. Hero

Abstract: $\textbf{Objective}$: To develop a multi-channel device event segmentation and feature extraction algorithm that is robust to changes in data distribution. $\textbf{Methods}… ▽ More $\textbf{Objective}$: To develop a multi-channel device event segmentation and feature extraction algorithm that is robust to changes in data distribution. $\textbf{Methods}$: We introduce an adaptive transfer learning algorithm to classify and segment events from non-stationary multi-channel temporal data. Using a multivariate hidden Markov model (HMM) and Fisher's linear discriminant analysis (FLDA) the algorithm adaptively adjusts to shifts in distribution over time. The proposed algorithm is unsupervised and learns to label events without requiring $\textit{a priori}$ information about true event states. The procedure is illustrated on experimental data collected from a cohort in a human viral challenge (HVC) study, where certain subjects have disrupted wake and sleep patterns after exposure to a H1N1 influenza pathogen. $\textbf{Results}$: Simulations establish that the proposed adaptive algorithm significantly outperforms other event classification methods. When applied to early time points in the HVC data the algorithm extracts sleep/wake features that are predictive of both infection and infection onset time. $\textbf{Conclusion}$: The proposed transfer learning event segmentation method is robust to temporal shifts in data distribution and can be used to produce highly discriminative event-labeled features for health monitoring. $\textbf{Significance}$: Our integrated multisensor signal processing and transfer learning method is applicable to many ambulatory monitoring applications. △ Less

Submitted 19 November, 2020; v1 submitted 20 August, 2020; originally announced August 2020.

Journal ref: IEEE Transactions on Biomedical Engineering, Nov. 17 2020

Holographic reconstruction of the interlayer distance of bilayer two-dimensional crystal samples from their convergent beam electron diffraction patterns

Authors: Tatiana Latychevskaia, Yichao Zou, Colin Robert Woods, Yi Bo Wang, Matthew Holwill, Eric Prestat, Sarah J. Haigh, Kostya S. Novoselov

Abstract: The convergent beam electron diffraction (CBED) patterns of twisted bilayer samples exhibit interference patterns in their CBED spots. Such interference patterns can be treated as off-axis holograms and the phase of the scattered waves, meaning the interlayer distance can be reconstructed. A detailed protocol of the reconstruction procedure is provided in this study. In addition, we derive an exac… ▽ More The convergent beam electron diffraction (CBED) patterns of twisted bilayer samples exhibit interference patterns in their CBED spots. Such interference patterns can be treated as off-axis holograms and the phase of the scattered waves, meaning the interlayer distance can be reconstructed. A detailed protocol of the reconstruction procedure is provided in this study. In addition, we derive an exact formula for reconstructing the interlayer distance from the recovered phase distribution, which takes into account the different chemical compositions of the individual monolayers. It is shown that one interference fringe in a CBED spot is sufficient to reconstruct the distance between the layers, which can be practical for imaging samples with a relatively small twist angle or when probing small sample regions. The quality of the reconstructed interlayer distance is studied as a function of the twist angle. At smaller twist angles, the reconstructed interlayer distance distribution is more precise and artefact free. At larger twist angles, artefacts due to the moiré structure appear in the reconstruction. A method for the reconstruction of the average interlayer distance is presented. As for resolution, the interlayer distance can be reconstructed by the holographic approach at an accuracy of 0.5 A, which is a few hundred times better than the intrinsic z-resolution of diffraction limited resolution, as expressed through the spread of the measured k-values. Moreover, we show that holographic CBED imaging can detect variations as small as 0.1 A in the interlayer distance, though the quantitative reconstruction of such variations suffers from large errors. △ Less

Submitted 9 August, 2020; originally announced August 2020.

Uncertainty quantification for epidemiological forecasts of COVID-19 through combinations of model predictions

Authors: D. S. Silk, V. E. Bowman, D. Semochkina, U. Dalrymple, D. C. Woods

Abstract: Scientific advice to the UK government throughout the COVID-19 pandemic has been informed by ensembles of epidemiological models provided by members of the Scientific Pandemic Influenza group on Modelling (SPI-M). Among other applications, the model ensembles have been used to forecast daily incidence, deaths and hospitalizations. The models differ in approach (e.g. deterministic or agent-based) a… ▽ More Scientific advice to the UK government throughout the COVID-19 pandemic has been informed by ensembles of epidemiological models provided by members of the Scientific Pandemic Influenza group on Modelling (SPI-M). Among other applications, the model ensembles have been used to forecast daily incidence, deaths and hospitalizations. The models differ in approach (e.g. deterministic or agent-based) and in assumptions made about the disease and population. These differences capture genuine uncertainty in the understanding of disease dynamics and in the choice of simplifying assumptions underpinning the model. Although analyses of multi-model ensembles can be logistically challenging when time-frames are short, accounting for structural uncertainty can improve accuracy and reduce the risk of over-confidence in predictions. In this study, we compare the performance of various ensemble methods to combine short-term (14 day) COVID-19 forecasts within the context of the pandemic response. We address practical issues around the availability of model predictions and make some initial proposals to address the short-comings of standard methods in this challenging situation. △ Less

Submitted 12 August, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

MSC Class: 62P10

Computer-aided whole-cell design: taking a holistic approach by integrating synthetic with systems biology

Authors: Lucia Marucci, Matteo Barberis, Jonathan Karr, Oliver Ray, Paul R. Race, Miguel de Souza Andrade, Claire Grierson, Stefan Andreas Hoffmann, Sophie Landon, Elibio Rech, Joshua Rees-Garbutt, Richard Seabrook, William Shaw, Christopher Woods

Abstract: Computer-aided design for synthetic biology promises to accelerate the rational and robust engineering of biological systems; it requires both detailed and quantitative mathematical and experimental models of the processes to (re)design, and software and tools for genetic engineering and DNA assembly. Ultimately, the increased precision in the design phase will have a dramatic impact on the produc… ▽ More Computer-aided design for synthetic biology promises to accelerate the rational and robust engineering of biological systems; it requires both detailed and quantitative mathematical and experimental models of the processes to (re)design, and software and tools for genetic engineering and DNA assembly. Ultimately, the increased precision in the design phase will have a dramatic impact on the production of designer cells and organisms with bespoke functions and increased modularity. Computer-aided design strategies require quantitative representations of cells, able to capture multiscale processes and link genotypes to phenotypes. Here, we present a perspective on how whole-cell, multiscale models could transform design-build-test-learn cycles in synthetic biology. We show how these models could significantly aid in the design and learn phases while reducing experimental testing by presenting case studies spanning from genome minimization to cell-free systems, and we discuss several challenges for the realization of our vision. The possibility to describe and build in silico whole-cells offers an opportunity to develop increasingly automatized, precise and accessible computer-aided design tools and strategies throughout novel interdisciplinary collaborations. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Nonmyopic and pseudo-nonmyopic approaches to optimal sequential design in the presence of covariates

Authors: Mia S. Tackney, David C. Woods, Ilya Shpitser

Abstract: In sequential experiments, subjects become available for the study over a period of time, and covariates are often measured at the time of arrival. We consider the setting where the sample size is fixed but covariate values are unknown until subjects enrol. Given a model for the outcome, a sequential optimal design approach can be used to allocate treatments to minimize the variance of the treatme… ▽ More In sequential experiments, subjects become available for the study over a period of time, and covariates are often measured at the time of arrival. We consider the setting where the sample size is fixed but covariate values are unknown until subjects enrol. Given a model for the outcome, a sequential optimal design approach can be used to allocate treatments to minimize the variance of the treatment effect. We extend existing optimal design methodology so it can be used within a nonmyopic framework, where treatment allocation for the current subject depends not only on the treatments and covariates of the subjects already enrolled in the study, but also the impact of possible future treatment assignments. The nonmyopic approach is computationally expensive as it requires recursive formulae. We propose a pseudo-nonmyopic approach which has a similar aim to the nonmyopic approach, but does not involve recursion and instead relies on simulations of future possible decisions. Our simulation studies show that the myopic approach is the most efficient for the logistic model case with a single binary covariate and binary treatment. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: 22 pages, 5 figures

Convergent beam electron diffraction of multilayer van der Waals structures

Authors: Tatiana Latychevskaia, Colin Robert Woods, Yi Bo Wang, Matthew Holwill, Eric Prestat, Sarah J. Haigh, Kostya S. Novoselov

Abstract: Convergent beam electron diffraction is routinely applied for studying deformation and local strain in thick crystals by matching the crystal structure to the observed intensity distributions. Recently, it has been demonstrated that CBED can be applied for imaging two-dimensional (2D) crystals where a direct reconstruction is possible and three-dimensional crystal deformations at a nanometre resol… ▽ More Convergent beam electron diffraction is routinely applied for studying deformation and local strain in thick crystals by matching the crystal structure to the observed intensity distributions. Recently, it has been demonstrated that CBED can be applied for imaging two-dimensional (2D) crystals where a direct reconstruction is possible and three-dimensional crystal deformations at a nanometre resolution can be retrieved. Here, we demonstrate that second-order effects allow for further information to be obtained regarding stacking arrangements between the crystals. Such effects are especially pronounced in samples consisting of multiple layers of 2D crystals. We show, using simulations and experiments, that twisted multilayer samples exhibit extra modulations of interference fringes in CBED patterns, i. e., a CBED moiré. A simple and robust method for the evaluation of the composition and the number of layers from a single-shot CBED pattern is demonstrated. △ Less

Submitted 28 March, 2020; originally announced March 2020.

Journal ref: Ultramicroscopy 212, 112976 (2020)

Nonlinear optical study of commensurability effects in graphene-hBN heterostructures

Authors: E. A. Stepanov, S. V. Semin, C. R. Woods, M. Vandelli, A. V. Kimel, K. S. Novoselov, M. I. Katsnelson

Abstract: Second-order nonlinear optical response allows to detect different properties of the system associated with the inversion symmetry breaking. Here, we use a second harmonic generation effect to investigate the alignment of a graphene/hexagonal Boron Nitride heterostructure. To achieve that, we activate a commensurate-incommensurate phase transition by a thermal annealing of the sample. We find that… ▽ More Second-order nonlinear optical response allows to detect different properties of the system associated with the inversion symmetry breaking. Here, we use a second harmonic generation effect to investigate the alignment of a graphene/hexagonal Boron Nitride heterostructure. To achieve that, we activate a commensurate-incommensurate phase transition by a thermal annealing of the sample. We find that this structural change in the system can be directly observed through a strong modification of a nonlinear optical signal. This result reveals the potential of a second harmonic generation technique for probing structural properties of layered systems. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Journal ref: ACS Applied Materials & Interfaces 12 (24), 27758-27764 (2020)

Composite super-moiré lattices in double aligned graphene heterostructures

Authors: Zihao Wang, Yi Bo Wang, J. Yin, E. Tóvári, Y. Yang, L. Lin, M. Holwill, J. Birkbeck, D. J. Perello, Shuigang Xu, J. Zultak, R. V. Gorbachev, A. V. Kretinin, T. Taniguchi, K. Watanabe, S. V. Morozov, M. Anđelković, S. P. Milovanović, L. Covaci, F. M. Peeters, A. Mishchenko, A. K. Geim, K. S. Novoselov, Vladimir I. Fal'ko, Angelika Knothe , et al. (1 additional authors not shown)

Abstract: When two-dimensional atomic crystals are brought into close proximity to form a van der Waals heterostructure, neighbouring crystals can start influencing each others electronic properties. Of particular interest is the situation when the periodicity of the two crystals closely match and a moiré pattern forms, which results in specific electron scattering, reconstruction of electronic and excitoni… ▽ More When two-dimensional atomic crystals are brought into close proximity to form a van der Waals heterostructure, neighbouring crystals can start influencing each others electronic properties. Of particular interest is the situation when the periodicity of the two crystals closely match and a moiré pattern forms, which results in specific electron scattering, reconstruction of electronic and excitonic spectra, crystal reconstruction, and many other effects. Thus, formation of moiré patterns is a viable tool of controlling the electronic properties of 2D materials. At the same time, the difference in the interatomic distances for the two crystals combined, determines the range in which the electronic spectrum is reconstructed, and thus is a barrier to the low energy regime. Here we present a way which allows spectrum reconstruction at all energies. By using graphene which is aligned simultaneously to two hexagonal boron nitride layers, one can make electrons scatter in the differential moiré pattern, which can have arbitrarily small wavevector and, thus results in spectrum reconstruction at arbitrarily low energies. We demonstrate that the strength of such a potential relies crucially on the atomic reconstruction of graphene within the differential moiré super-cell. Such structures offer further opportunity in tuning the electronic spectra of two-dimensional materials. △ Less

Submitted 27 December, 2019; originally announced December 2019.

Comments: 35 pages, 15 figures

Journal ref: Science Advances 20 Dec 2019: Vol. 5, no. 12, eaay8897

Piezoelectricity in monolayer hexagonal boron nitride

Authors: P. Ares, T. Cea, M. Holwill, Y. B. Wang, R. Roldan, F. Guinea, D. V. Andreeva, L. Fumagalli, K. S. Novoselov, C. R. Woods

Abstract: Two-dimensional (2D) hexagonal boron nitride (hBN) is a wide-bandgap van der Waals crystal with a unique combination of properties, including exceptional strength, large oxidation resistance at high temperatures and optical functionalities. Furthermore, in recent years hBN crystals have become the material of choice for encapsulating other 2D crystals in a variety of technological applications, fr… ▽ More Two-dimensional (2D) hexagonal boron nitride (hBN) is a wide-bandgap van der Waals crystal with a unique combination of properties, including exceptional strength, large oxidation resistance at high temperatures and optical functionalities. Furthermore, in recent years hBN crystals have become the material of choice for encapsulating other 2D crystals in a variety of technological applications, from optoelectronic and tunnelling devices to composites. Monolayer hBN, which has no center of symmetry, has been predicted to exhibit piezoelectric properties, yet experimental evidence is lacking. Here, by using electrostatic force microscopy, we observed this effect as a strain-induced change in the local electric field around bubbles and creases, in agreement with theoretical calculations. No piezoelectricity was found in bilayer and bulk hBN, where the centre of symmetry is restored. These results add piezoelectricity to the known properties of monolayer hBN, which makes it a desirable candidate for novel electromechanical and stretchable optoelectronic devices, and pave a way to control the local electric field and carrier concentration in van der Waals heterostructures via strain. The experimental approach used here also shows a way to investigate the piezoelectric properties of other materials on the nanoscale by using electrostatic scanning probe techniques. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Stacking Order in Graphite Films Controlled by Van der Waals Technology

Authors: Yaping Yang, Yi-Chao Zou, Colin R. Woods, Yanmeng Shi, Jun Yin, Shuigang Xu, Servet Ozdemir, Takashi Taniguchi, Kenji Watanabe, Andre K. Geim, Kostya S. Novoselov, Sarah J. Haigh, Artem Mishchenko

Abstract: In graphite crystals, layers of graphene reside in three equivalent, but distinct, stacking positions typically referred to as A, B, and C projections. The order in which the layers are stacked defines the electronic structure of the crystal, providing an exciting degree of freedom which can be exploited for designing graphitic materials with unusual properties including predicted high-temperature… ▽ More In graphite crystals, layers of graphene reside in three equivalent, but distinct, stacking positions typically referred to as A, B, and C projections. The order in which the layers are stacked defines the electronic structure of the crystal, providing an exciting degree of freedom which can be exploited for designing graphitic materials with unusual properties including predicted high-temperature superconductivity and ferromagnetism. However, the lack of control of the stacking sequence limits most research to the stable ABA form of graphite. Here we demonstrate a strategy to control the stacking order using van der Waals technology. To this end, we first visualise the distribution of stacking domains in graphite films and then perform directional encapsulation of ABC-rich graphite crystallites with hexagonal boron nitride (hBN). We found that hBN-encapsulation which is introduced parallel to the graphite zigzag edges preserves ABC stacking, while encapsulation along the armchair edges transforms the stacking to ABA. The technique presented here should facilitate new research on the important properties of ABC graphite. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Bayes Linear Emulation of Simulator Networks

Authors: Samuel E. Jackson, David C. Woods

Abstract: Computationally expensive simulators, implementing mathematical models in computer codes, are commonly approximated using statistical emulators. We develop and assess novel emulation methods for systems best modelled via a chain, series or network of simulators. Using a Bayes linear framework, we link statistical emulators of the component simulators to explicitly account for the simulator input u… ▽ More Computationally expensive simulators, implementing mathematical models in computer codes, are commonly approximated using statistical emulators. We develop and assess novel emulation methods for systems best modelled via a chain, series or network of simulators. Using a Bayes linear framework, we link statistical emulators of the component simulators to explicitly account for the simulator input uncertainty induced by links between models in arbitrarily large networks. We demonstrate the advantages of these methods compared to use of a single emulator of the composite simulator network for a variety of examples, including the motivating epidemiological simulator chain to model the impact of an airborne infectious disease. △ Less

Submitted 25 August, 2021; v1 submitted 17 October, 2019; originally announced October 2019.

Comments: 25 pages, 7 figures

An unsupervised transfer learning algorithm for sleep monitoring

Authors: Xichen She, Yaya Zhai, Ricardo Henao, Christopher W. Woods, Geoffrey S. Ginsburg, Peter X. K. Song, Alfred O. Hero

Abstract: Objective: To develop multisensor-wearable-device sleep monitoring algorithms that are robust to health disruptions affecting sleep patterns. Methods: We develop an unsupervised transfer learning algorithm based on a multivariate hidden Markov model and Fisher's linear discriminant analysis, adaptively adjusting to sleep pattern shift by training on dynamics of sleep/wake states. The proposed algo… ▽ More Objective: To develop multisensor-wearable-device sleep monitoring algorithms that are robust to health disruptions affecting sleep patterns. Methods: We develop an unsupervised transfer learning algorithm based on a multivariate hidden Markov model and Fisher's linear discriminant analysis, adaptively adjusting to sleep pattern shift by training on dynamics of sleep/wake states. The proposed algorithm operates, without requiring a priori information about true sleep/wake states, by establishing an initial training set with hidden Markov model and leveraging a taper window mechanism to learn the sleep pattern in an incremental fashion. Our domain-adaptation algorithm is applied to a dataset collected in a human viral challenge study to identify sleep/wake periods of both uninfected and infected participants. Results: The algorithm successfully detects sleep/wake sessions in subjects whose sleep patterns are disrupted by respiratory infection (H3N2 flu virus). Pre-symptomatic features based on the detected periods are found to be strongly predictive of both infection status (AUC = 0.844) and infection onset time (AUC = 0.885), indicating the effectiveness and usefulness of the algorithm. Conclusion: Our method can effectively detect sleep/wake states in the presence of sleep pattern shift. Significance: Utilizing integrated multisensor signal processing and adaptive training schemes, our algorithm is able to capture key sleep patterns in ambulatory monitoring, leading to better automated sleep assessment and prediction. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Copula-based robust optimal block designs

Authors: W. G. Mueller, A. Rappold, D. C. Woods

Abstract: Blocking is often used to reduce known variability in designed experiments by collecting together homogeneous experimental units. A common modelling assumption for such experiments is that responses from units within a block are dependent. Accounting for such dependencies in both the design of the experiment and the modelling of the resulting data when the response is not normally distributed can… ▽ More Blocking is often used to reduce known variability in designed experiments by collecting together homogeneous experimental units. A common modelling assumption for such experiments is that responses from units within a block are dependent. Accounting for such dependencies in both the design of the experiment and the modelling of the resulting data when the response is not normally distributed can be challenging, particularly in terms of the computation required to find an optimal design. The application of copulas and marginal modelling provides a computationally efficient approach for estimating population-average treatment effects. Motivated by an experiment from materials testing, we develop and demonstrate designs with blocks of size two using copula models. Such designs are also important in applications ranging from microarray experiments to experiments on human eyes or limbs with naturally occurring blocks of size two. We present methodology for design selection, make comparisons to existing approaches in the literature and assess the robustness of the designs to modelling assumptions. △ Less

Submitted 6 November, 2018; originally announced November 2018.

MSC Class: 62K05; 62K10

Convergent and divergent beam electron holography and reconstruction of adsorbates on free-standing two-dimensional crystals

Authors: Tatiana Latychevskaia, Colin Robert Woods, Yi Bo Wang, Matthew Holwill, Eric Prestat, Sarah J. Haigh, Kostya S. Novoselov

Abstract: Van der Waals heterostructures have been lately intensively studied because they offer a large variety of properties that can be controlled by selecting 2D materials and their sequence in the stack. The exact arrangement of the layers as well as the exact arrangement of the atoms within the layers, both are important for the properties of the resulting device. Recently it has been demonstrated tha… ▽ More Van der Waals heterostructures have been lately intensively studied because they offer a large variety of properties that can be controlled by selecting 2D materials and their sequence in the stack. The exact arrangement of the layers as well as the exact arrangement of the atoms within the layers, both are important for the properties of the resulting device. Recently it has been demonstrated that convergent beam electron diffraction (CBED) allows quantitative three-dimensional mapping of atomic positions in three-dimensional materials from a single CBED pattern. In this study we investigate CBED in more detail by simulating and performing various CBED regimes, with convergent and divergent wavefronts, on a somewhat simplified system: a 2D monolayer crystal. In CBED, each CBED spot is in fact an in-line hologram of the sample, where in-line holography is known to exhibit high intensity contrast in detection of weak phase objects that are not detectable in conventional in-focus imaging mode. Adsorbates exhibit strong intensity contrast in zero and higher order CBED spots, whereas lattice deformation such as strain or rippling cause noticeable intensity contrast only in the first and higher order CBED spots. The individual CBED spots can be reconstructed as typical in-line holograms, and the resolution of 2.13 A can be in principle achieved in the reconstructions. We provide simulated and experimental examples of CBED of a 2D monolayer crystal. The simulations show that individual CBED spots can be treated as in-line holograms and sample distributions such as adsorbates, can be reconstructed. Individual atoms can be reconstructed from a single CBED pattern provided the later exhibits high-order CBED spots. Examples of reconstructions obtained from experimental CBED patterns, at a resolution of 2.7 A, are shown. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Journal ref: Frontiers of Physics 14(1), 13606 (2019)

Convergent beam electron holography for analysis of van der Waals heterostructures

Authors: Tatiana Latychevskaia, Colin Robert Woods, Yi Bo Wang, Matthew Holwill, Eric Prestat, Sarah J. Haigh, Kostya S. Novoselov

Abstract: Van der Waals heterostructures, which explore the synergetic properties of two-dimensional (2D) materials when assembled into three-dimensional stacks, have already brought to life a number of exciting new phenomena and novel electronic devices. Still, the interaction between the layers in such assembly, possible surface reconstruction, intrinsic and extrinsic defects are very difficult to charact… ▽ More Van der Waals heterostructures, which explore the synergetic properties of two-dimensional (2D) materials when assembled into three-dimensional stacks, have already brought to life a number of exciting new phenomena and novel electronic devices. Still, the interaction between the layers in such assembly, possible surface reconstruction, intrinsic and extrinsic defects are very difficult to characterise by any method, because of the single-atomic nature of the crystals involved. Here we present a convergent beam electron holographic technique which allows imaging of the stacking order in such heterostructures. Based on the interference of electron waves scattered on different crystals in the stack, this approach allows one to reconstruct the relative rotation, stretching, out-of-plane corrugation of the layers with atomic precision. Being holographic in nature, our approach allows extraction of quantitative information about the three-dimensional structure of the typical defects from a single image covering thousands of square nanometres. Furthermore, qualitative information about the defects in the stack can be extracted from the convergent diffraction patterns even without reconstruction - simply by comparing the patterns in different diffraction spots. We expect that convergent beam electron holography will be widely used to study the properties of van der Waals heterostructures. △ Less

Submitted 5 July, 2018; originally announced July 2018.

Journal ref: Proceedings of the National Academy of Sciences, Volume 115, Issue 29, pages 7473-7478, 2018

Blip: JIT and Footloose On The Edge

Authors: Andy Edmonds, Chris Woods, Ana Juan Ferrer, Juan Francisco Ribera, Thomas Micheal Bohnert

Abstract: Edge environments offer a number of advantages for software developers including the ability to create services which can offer lower latency, better privacy, and reduced operational costs than traditional cloud hosted services. However large technical challenges exist, which prevent developers from utilising the Edge; complexities related to the heterogeneous nature of the Edge environment, issue… ▽ More Edge environments offer a number of advantages for software developers including the ability to create services which can offer lower latency, better privacy, and reduced operational costs than traditional cloud hosted services. However large technical challenges exist, which prevent developers from utilising the Edge; complexities related to the heterogeneous nature of the Edge environment, issues with orchestration and application management and lastly, the inherent issues in creating decentralised distributed applications which operate at a large geographic scale. In this conceptual and architectural paper we envision a solution, Blip, which offers an easy to use programming and operational environment which addresses the these issues. It aims to remove the technical barriers which will inhibit the wider adoption Edge application development. This paper validates the Blip concept by demonstrating how it will deliver on the advantages of the Edge for a familiar scenario. △ Less

Submitted 31 May, 2018; originally announced June 2018.

Summary of effect aliasing structure (SEAS): new descriptive statistics for factorial and supersaturated designs

Authors: Frederick Kin Hing Phoa, Yi-Hua Liao, David C. Woods, Shah-Kae Chou

Abstract: In the assessment and selection of supersaturated designs, the aliasing structure of interaction effects is usually ignored by traditional criteria such as $E(s^2)$-optimality. We introduce the Summary of Effect Aliasing Structure (SEAS) for assessing the aliasing structure of supersaturated designs, and other non-regular fractional factorial designs, that takes account of interaction terms and pr… ▽ More In the assessment and selection of supersaturated designs, the aliasing structure of interaction effects is usually ignored by traditional criteria such as $E(s^2)$-optimality. We introduce the Summary of Effect Aliasing Structure (SEAS) for assessing the aliasing structure of supersaturated designs, and other non-regular fractional factorial designs, that takes account of interaction terms and provides more detail than usual summaries such as (generalized) resolution and wordlength patterns. The new summary consists of three criteria, abbreviated as MAP: (1) Maximum dependency aliasing pattern; (2) Average square aliasing pattern; and (3) Pairwise dependency ratio. These criteria provided insight when traditional criteria fail to differentiate between designs. We theoretically study the relationship between the MAP criteria and traditional quantities, and demonstrate the use of SEAS for comparing some example supersaturated designs, including designs suggested in the literature. We also propose a variant of SEAS to measure the aliasing structure for individual columns of a design, and use it to choose assignments of factors to columns for an $E(s^2)$-optimal design. △ Less

Submitted 15 January, 2021; v1 submitted 30 November, 2017; originally announced November 2017.

MSC Class: 62K15

Large tunable valley splitting in edge-free graphene quantum dots on boron nitride

Authors: Nils M. Freitag, Tobias Reisch, Larisa A. Chizhova, Peter Nemes-Incze, Christian Holl, Colin R. Woods, Roman V. Gorbachev, Yang Cao, Andre K. Geim, Kostya S. Novoselov, Joachim Burgdörfer, Florian Libisch, Markus Morgenstern

Abstract: Coherent manipulation of binary degrees of freedom is at the heart of modern quantum technologies. Graphene offers two binary degrees: the electron spin and the valley. Efficient spin control has been demonstrated in many solid state systems, while exploitation of the valley has only recently been started, yet without control on the single electron level. Here, we show that van-der Waals stacking… ▽ More Coherent manipulation of binary degrees of freedom is at the heart of modern quantum technologies. Graphene offers two binary degrees: the electron spin and the valley. Efficient spin control has been demonstrated in many solid state systems, while exploitation of the valley has only recently been started, yet without control on the single electron level. Here, we show that van-der Waals stacking of graphene onto hexagonal boron nitride offers a natural platform for valley control. We use a graphene quantum dot induced by the tip of a scanning tunneling microscope and demonstrate valley splitting that is tunable from -5 to +10 meV (including valley inversion) by sub-10-nm displacements of the quantum dot position. This boosts the range of controlled valley splitting by about one order of magnitude. The tunable inversion of spin and valley states should enable coherent superposition of these degrees of freedom as a first step towards graphene-based qubits. △ Less

Submitted 20 March, 2018; v1 submitted 30 August, 2017; originally announced August 2017.

Journal ref: Nature Nanotechnology (2018)

A Novel Potential Field Controller for Use on Aerial Robots

Authors: Alexander C. Woods, Hung M. La

Abstract: Unmanned Aerial Vehicles (UAV), commonly known as drones, have many potential uses in real world applications. Drones require advanced planning and navigation algorithms to enable them to safely move through and interact with the world around them. This paper presents an extended potential field controller (ePFC) which enables an aerial robot, or drone, to safely track a dynamic target location wh… ▽ More Unmanned Aerial Vehicles (UAV), commonly known as drones, have many potential uses in real world applications. Drones require advanced planning and navigation algorithms to enable them to safely move through and interact with the world around them. This paper presents an extended potential field controller (ePFC) which enables an aerial robot, or drone, to safely track a dynamic target location while simultaneously avoiding any obstacles in its path. The ePFC outperforms a traditional potential field controller (PFC) with smoother tracking paths and shorter settling times. The proposed ePFC's stability is evaluated by Lyapunov approach, and its performance is simulated in a Matlab environment. Finally, the controller is implemented on an experimental platform in a laboratory environment which demonstrates the effectiveness of the controller. △ Less

Submitted 15 April, 2017; originally announced April 2017.

Phase segregation facilitates exfoliation of franckeite crystals to a single unit cell thickness

Authors: Matěj Velický, Peter S. Toth, Alexander M. Rakowski, Aidan P. Rooney, Aleksey Kozikov, Artem Mishchenko, Colin R. Woods, Thanasis Georgiou, Sarah J. Haigh, Kostya S. Novoselov, Robert A. W. Dryfe

Abstract: Weak interlayer van der Waals interactions in bulk crystals facilitate their mechanical exfoliation to monolayer and few-layer two-dimensional (2D) materials, which exhibit striking physical phenomena absent in their bulk form. Here we study a 2D form of a mineral franckeite and show that phase segregation into discrete layers at the sub-nanometre scale facilitates its layered structure and basal… ▽ More Weak interlayer van der Waals interactions in bulk crystals facilitate their mechanical exfoliation to monolayer and few-layer two-dimensional (2D) materials, which exhibit striking physical phenomena absent in their bulk form. Here we study a 2D form of a mineral franckeite and show that phase segregation into discrete layers at the sub-nanometre scale facilitates its layered structure and basal cleavage. This behaviour is likely to be common in a wider family of complex crystals and could be exploited for a single-step synthesis of van der Waals heterostructures, as an alternative to stacking of 2D materials. Mechanical exfoliation allowed us to produce crystals down to a single unit cell thickness and rationalise its basal cleavage by atomic-resolution scanning transmission electron microscopy (STEM). We demonstrate p-type electrical conductivity and remarkable electrochemical properties in exfoliated crystals, which shows promise for energy storage applications. △ Less

Submitted 8 December, 2016; originally announced December 2016.

Comments: Main text: 19 pages, 5 figures; Supplementary Information: 19 pages, 10 figures

Journal ref: Nature Communications 8 (2017) 14410

Stacking transition in bilayer graphene caused by thermally activated rotation

Authors: Mengjian Zhu, Davit Ghazaryan, Seok-Kyun Son, Colin R. Woods, Abhishek Misra, Lin He, Takashi Taniguchi, Kenji Watanabe, Kostya S. Novoselov, Yang Cao, Artem Mishchenko

Abstract: Crystallographic alignment between two-dimensional crystals in van der Waals heterostructures brought a number of profound physical phenomena, including observation of Hofstadter butterfly and topological currents, and promising novel applications, such as resonant tunnelling transistors. Here, by probing the electronic density of states in graphene using graphene-hexagonal boron nitride tunnellin… ▽ More Crystallographic alignment between two-dimensional crystals in van der Waals heterostructures brought a number of profound physical phenomena, including observation of Hofstadter butterfly and topological currents, and promising novel applications, such as resonant tunnelling transistors. Here, by probing the electronic density of states in graphene using graphene-hexagonal boron nitride tunnelling transistors, we demonstrate a structural transition of bilayer graphene from incommensurate twisted stacking state into a commensurate AB stacking due to a macroscopic graphene self-rotation. This structural transition is accompanied by a topological transition in the reciprocal space and by pseudospin texturing. The stacking transition is driven by van der Waals interaction energy of the two graphene layers and is thermally activated by unpinning the microscopic chemical adsorbents which are then removed by the self-cleaning of graphene. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 13 pages, 6 figures

Journal ref: 2D Mater. 4 (2017) 011013

Tuning the pseudospin polarization of graphene by a pseudo-magnetic field

Authors: Alexander Georgi, Peter Nemes-Incze, Ramon Carrillo-Bastos, Daiara Faria, Silvia Viola Kusminskiy, Dawei Zhai, Martin Schneider, Dinesh Subramaniam, Torge Mashoff, Nils M. Freitag, Marcus Liebmann, Marco Pratzer, Ludger Wirtz, Colin R. Woods, Roman V. Gorbachev, Yang Cao, Kostya S. Novoselov, Nancy Sandler, Markus Morgenstern

Abstract: One of the intriguing characteristics of honeycomb lattices is the appearance of a pseudo-magnetic field as a result of mechanical deformation. In the case of graphene, the Landau quantization resulting from this pseudo-magnetic field has been measured using scanning tunneling microscopy. Here we show that a signature of the pseudo-magnetic field is a local sublattice symmetry breaking observable… ▽ More One of the intriguing characteristics of honeycomb lattices is the appearance of a pseudo-magnetic field as a result of mechanical deformation. In the case of graphene, the Landau quantization resulting from this pseudo-magnetic field has been measured using scanning tunneling microscopy. Here we show that a signature of the pseudo-magnetic field is a local sublattice symmetry breaking observable as a redistribution of the local density of states. This can be interpreted as a polarization of graphene's pseudospin due to a strain induced pseudo-magnetic field, in analogy to the alignment of a real spin in a magnetic field. We reveal this sublattice symmetry breaking by tunably straining graphene using the tip of a scanning tunneling microscope. The tip locally lifts the graphene membrane from a SiO$_2$ support, as visible by an increased slope of the $I(z)$ curves. The amount of lifting is consistent with molecular dynamics calculations, which reveal a deformed graphene area under the tip in the shape of a Gaussian. The pseudo-magnetic field induced by the deformation becomes visible as a sublattice symmetry breaking which scales with the lifting height of the strained deformation and therefore with the pseudo-magnetic field strength. Its magnitude is quantitatively reproduced by analytic and tight-binding models, revealing fields of 1000 T. These results might be the starting point for an effective THz valley filter, as a basic element of valleytronics. △ Less

Submitted 6 March, 2017; v1 submitted 18 November, 2016; originally announced November 2016.

Comments: Revised manuscript: streamlined the abstract and introduction, added methods to supplement, Nano Letters, 2017

Electrostatically confined monolayer graphene quantum dots with orbital and valley splittings

Authors: Nils M. Freitag, Larisa A. Chizhova, Peter Nemes-Incze, Colin R. Woods, Roman V. Gorbachev, Yang Cao, Andre K. Geim, Kostya S. Novoselov, Joachim Burgdörfer, Florian Libisch, Markus Morgenstern

Abstract: The electrostatic confinement of massless charge carriers is hampered by Klein tunneling. Circumventing this problem in graphene mainly relies on carving out nanostructures or applying electric displacement fields to open a band gap in bilayer graphene. So far, these approaches suffer from edge disorder or insufficiently controlled localization of electrons. Here we realize an alternative strategy… ▽ More The electrostatic confinement of massless charge carriers is hampered by Klein tunneling. Circumventing this problem in graphene mainly relies on carving out nanostructures or applying electric displacement fields to open a band gap in bilayer graphene. So far, these approaches suffer from edge disorder or insufficiently controlled localization of electrons. Here we realize an alternative strategy in monolayer graphene, by combining a homogeneous magnetic field and electrostatic confinement. Using the tip of a scanning tunneling microscope, we induce a confining potential in the Landau gaps of bulk graphene without the need for physical edges. Gating the localized states towards the Fermi energy leads to regular charging sequences with more than 40 Coulomb peaks exhibiting typical addition energies of 7-20 meV. Orbital splittings of 4-10 meV and a valley splitting of about 3 meV for the first orbital state can be deduced. These experimental observations are quantitatively reproduced by tight binding calculations, which include the interactions of the graphene with the aligned hexagonal boron nitride substrate. The demonstrated confinement approach appears suitable to create quantum dots with well-defined wave function properties beyond the reach of traditional techniques. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Bayesian design of experiments for generalised linear models and dimensional analysis with industrial and scientific application

Authors: David C. Woods, Antony M. Overstall, Maria Adamou, Timothy W. Waite

Abstract: The design of an experiment can be always be considered at least implicitly Bayesian, with prior knowledge used informally to aid decisions such as the variables to be studied and the choice of a plausible relationship between the explanatory variables and measured responses. Bayesian methods allow uncertainty in these decisions to be incorporated into design selection through prior distributions… ▽ More The design of an experiment can be always be considered at least implicitly Bayesian, with prior knowledge used informally to aid decisions such as the variables to be studied and the choice of a plausible relationship between the explanatory variables and measured responses. Bayesian methods allow uncertainty in these decisions to be incorporated into design selection through prior distributions that encapsulate information available from scientific knowledge or previous experimentation. Further, a design may be explicitly tailored to the aim of the experiment through a decision-theoretic approach using an appropriate loss function. We review the area of decision-theoretic Bayesian design, with particular emphasis on recent advances in computational methods. For many problems arising in industry and science, experiments result in a discrete response that is well described by a member of the class of generalised linear models. We describe how Gaussian process emulation, commonly used in computer experiments, can play an important role in facilitating Bayesian design for realistic problems. A main focus is the combination of Gaussian process regression to approximate the expected loss with cyclic descent (coordinate exchange) optimisation algorithms to allow optimal designs to be found for previously infeasible problems. We also present the first optimal design results for statistical models formed from dimensional analysis, a methodology widely employed in the engineering and physical sciences to produce parsimonious and interpretable models. Using the famous paper helicopter experiment, we show the potential for the combination of Bayesian design, generalised linear models and dimensional analysis to produce small but informative experiments. △ Less

Submitted 26 October, 2016; v1 submitted 19 June, 2016; originally announced June 2016.

MSC Class: 62K05; 62K20; 62J12; 62C10

Journal ref: Quality Engineering, 29, 91-103, 2017

Macroscopic self-reorientation of interacting two-dimensional crystals

Authors: C. R. Woods, F. Withers, M. J. Zhu, Y. Cao, G. Yu, A. Kozikov, M. Ben Shalom, S. V. Morozov, M. M. van Wijk, A. Fasolino, M. I. Katsnelson, K. Watanabe, T. Taniguchi, A. K. Geim, A. Mishchenko, K. S. Novoselov

Abstract: Microelectromechanical systems, which can be moved or rotated with nanometre precision, already find applications in such fields like radio-frequency electronics, micro-attenuators, sensors and many others. Especially interesting are those which allow fine control over the motion on atomic scale due to self-alignment mechanisms and forces acting on the atomic level. Such machines can produce well-… ▽ More Microelectromechanical systems, which can be moved or rotated with nanometre precision, already find applications in such fields like radio-frequency electronics, micro-attenuators, sensors and many others. Especially interesting are those which allow fine control over the motion on atomic scale due to self-alignment mechanisms and forces acting on the atomic level. Such machines can produce well-controlled movements as a reaction to small changes of the external parameters. Here we demonstrate that, for the system of graphene on hexagonal boron nitride, the interplay between the van der Waals and elastic energies results in graphene mechanically self-rotating towards the hexagonal boron nitride crystallographic directions. Such rotation is macroscopic (for graphene flakes of tens of micrometres the tangential movement can be on hundreds of nanometres) and can be used for reproducible manufacturing of aligned van der Waals heterostructures. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 32 pages, 17 figures, includes supplementary information, published in Nature Communications (2016)

Journal ref: Nature Communications 7, Article number:10800, 2016

Model selection via Bayesian information capacity designs for generalised linear models

Authors: David C. Woods, James M. McGree, Susan M. Lewis

Abstract: The first investigation is made of designs for screening experiments where the response variable is approximated by a generalised linear model. A Bayesian information capacity criterion is defined for the selection of designs that are robust to the form of the linear predictor. For binomial data and logistic regression, the effectiveness of these designs for screening is assessed through simulatio… ▽ More The first investigation is made of designs for screening experiments where the response variable is approximated by a generalised linear model. A Bayesian information capacity criterion is defined for the selection of designs that are robust to the form of the linear predictor. For binomial data and logistic regression, the effectiveness of these designs for screening is assessed through simulation studies using all-subsets regression and model selection via maximum penalised likelihood and a generalised information criterion. For Poisson data and log-linear regression, similar assessments are made using maximum likelihood and the Akaike information criterion for minimally-supported designs that are constructed analytically. The results show that effective screening, that is, high power with moderate type I error rate and false discovery rate, can be achieved through suitable choices for the number of design support points and experiment size. Logistic regression is shown to present a more challenging problem than log-linear regression. Some areas for future work are also indicated. △ Less

Submitted 26 October, 2016; v1 submitted 29 January, 2016; originally announced January 2016.

MSC Class: 62K05; 62K20; 62J12

Emulation of multivariate simulators using thin-plate splines with application to atmospheric dispersion

Authors: Veronica E. Bowman, David C. Woods

Abstract: It is often desirable to build a statistical emulator of a complex computer simulator in order to perform analysis which would otherwise be computationally infeasible. We propose methodology to model multivariate output from a computer simulator taking into account output structure in the responses. The utility of this approach is demonstrated by applying it to a chemical and biological hazard pre… ▽ More It is often desirable to build a statistical emulator of a complex computer simulator in order to perform analysis which would otherwise be computationally infeasible. We propose methodology to model multivariate output from a computer simulator taking into account output structure in the responses. The utility of this approach is demonstrated by applying it to a chemical and biological hazard prediction model. Predicting the hazard area that results from an accidental or deliberate chemical or biological release is imperative in civil and military planning and also in emergency response. The hazard area resulting from such a release is highly structured in space and we therefore propose the use of a thin-plate spline to capture the spatial structure and fit a Gaussian process emulator to the coefficients of the resultant basis functions. We compare and contrast four different techniques for emulating multivariate output: dimension-reduction using (i) a fully Bayesian approach with a principal component basis, (ii) a fully Bayesian approach with a thin-plate spline basis, assuming that the basis coefficients are independent, and (iii) a "plug-in" Bayesian approach with a thin-plate spline basis and a separable covariance structure; and (iv) a functional data modeling approach using a tensor-product (separable) Gaussian process. We develop methodology for the two thin-plate spline emulators and demonstrate that these emulators significantly outperform the principal component emulator. Further, the separable thin-plate spline emulator, which accounts for the dependence between basis coefficients, provides substantially more realistic quantification of uncertainty, and is also computationally more tractable, allowing fast emulation. For high resolution output data, it also offers substantial predictive and computational advantages over the tensor-product Gaussian process emulator. △ Less

Submitted 26 October, 2016; v1 submitted 23 December, 2015; originally announced December 2015.

MSC Class: 62-07; 62G08; 62P12

Journal ref: SIAM/ASA Journal of Uncertainty Quantification, 4, 1323-1344, 2016

WSe2 light-emitting tunneling transistors with enhanced brightness at room temperature

Authors: F. Withers, O. Del Pozo-Zamudio, S. Schwarz, S. Dufferwiel, P. M. Walker, T. Godde, A. P. Rooney, A. Gholinia, C. R. Woods, P. Blake, S. J. Haigh, K. Watanabe, T. Taniguchi, I. L. Aleiner, A. K. Geim, V. I. Falko, A. I. Tartakovskii, K. S. Novoselov

Abstract: Monolayers of molybdenum and tungsten dichalcogenides are direct bandgap semiconductors, which makes them promising for opto-electronic applications. In particular, van der Waals heterostructures consisting of monolayers of MoS2 sandwiched between atomically thin hexagonal boron nitride (hBN) and graphene electrodes allows one to obtain light emitting quantum wells (LEQWs) with low-temperature ext… ▽ More Monolayers of molybdenum and tungsten dichalcogenides are direct bandgap semiconductors, which makes them promising for opto-electronic applications. In particular, van der Waals heterostructures consisting of monolayers of MoS2 sandwiched between atomically thin hexagonal boron nitride (hBN) and graphene electrodes allows one to obtain light emitting quantum wells (LEQWs) with low-temperature external quantum efficiency (EQE) of 1%. However, the EQE of MoS2 and MoSe2-based LEQWs shows behavior common for many other materials: it decreases fast from cryogenic conditions to room temperature, undermining their practical applications. Here we compare MoSe2 and WSe2 LEQWs. We show that the EQE of WSe2 devices grows with temperature, with room temperature EQE reaching 5%, which is 250x more than the previous best performance of MoS2 and MoSe2 quantum wells in ambient conditions. We attribute such a different temperature dependences to the inverted sign of spin-orbit splitting of conduction band states in tungsten and molybdenum dichalcogenides, which makes the lowest-energy exciton in WSe2 dark. △ Less

Submitted 19 November, 2015; originally announced November 2015.

Designs for Generalized Linear Models

Authors: Anthony C. Atkinson, David C. Woods

Abstract: This paper reviews the design of experiments for generalised linear models, including optimal design, Bayesian design and designs for models with random effects. This paper reviews the design of experiments for generalised linear models, including optimal design, Bayesian design and designs for models with random effects. △ Less

Submitted 18 October, 2015; originally announced October 2015.

MSC Class: 62K05

Journal ref: Handbook of Design and Analysis of Experiments (editors A. Dean, M. Morris, J. Stufken, D. Bingham), chapter 13, 2015, Chapman and Hall/CRC Press

Design of Experiments for Screening

Authors: David C. Woods, Susan M. Lewis

Abstract: The aim of this paper is to review methods of designing screening experiments, ranging from designs originally developed for physical experiments to those especially tailored to experiments on numerical models. The strengths and weaknesses of the various designs for screening variables in numerical models are discussed. First, classes of factorial designs for experiments to estimate main effects a… ▽ More The aim of this paper is to review methods of designing screening experiments, ranging from designs originally developed for physical experiments to those especially tailored to experiments on numerical models. The strengths and weaknesses of the various designs for screening variables in numerical models are discussed. First, classes of factorial designs for experiments to estimate main effects and interactions through a linear statistical model are described, specifically regular and nonregular fractional factorial designs, supersaturated designs and systematic fractional replicate designs. Generic issues of aliasing, bias and cancellation of factorial effects are discussed. Second, group screening experiments are considered including factorial group screening and sequential bifurcation. Third, random sampling plans are discussed including Latin hypercube sampling and sampling plans to estimate elementary effects. Fourth, a variety of modelling methods commonly employed with screening designs are briefly described. Finally, a novel study demonstrates six screening methods on two frequently-used exemplars, and their performances are compared. △ Less

Submitted 18 October, 2015; originally announced October 2015.

MSC Class: 62K15

Nonlocal Response and Anamorphosis: the Case of Few-Layer Black Phosphorus

Authors: Artem Mishchenko, Yang Cao, Geliang Yu, Colin R. Woods, Roman V. Gorbachev, Kostya S. Novoselov, Andre K. Geim, Leonid S. Levitov

Abstract: Few-layer black phosphorus was recently rediscovered as a narrow-bandgap atomically thin semiconductor and has already attracted unprecedented attention due to its interesting properties. One feature of this material that sets it apart from other atomically thin crystals is its structural in-plane anisotropy which manifests in strongly anisotropic transport characteristics. However, traditional an… ▽ More Few-layer black phosphorus was recently rediscovered as a narrow-bandgap atomically thin semiconductor and has already attracted unprecedented attention due to its interesting properties. One feature of this material that sets it apart from other atomically thin crystals is its structural in-plane anisotropy which manifests in strongly anisotropic transport characteristics. However, traditional angle-resolved conductance measurements present a challenge for nanoscale systems such as black phosphorus, calling for new approaches in precision studies of transport anisotropy. Here we show that the nonlocal response, being exponentially sensitive to the anisotropy value, provides a powerful tool for determining the anisotropy. This is established by combining measurements of the orientation-dependent nonlocal resistance response with the analysis based on the anamorphosis relations. We demonstrate that the nonlocal response can differ by orders of magnitude for different crystallographic directions even when the anisotropy is at most order-one, allowing us to extract accurate anisotropy values. △ Less

Submitted 22 September, 2015; originally announced September 2015.

Journal ref: Nano Lett. 15 (10), pp 6991-6995 (2015)

Quantum oscillations of the critical current and high-field superconducting proximity in ballistic graphene

Authors: M. Ben Shalom, M. J. Zhu, V. I. Fal'ko, A. Mishchenko, A. V. Kretinin, K. S. Novoselov, C. R. Woods, K. Watanabe, T. Taniguchi, A. K. Geim, J. R. Prance

Abstract: Graphene-based Josephson junctions provide a novel platform for studying the proximity effect due to graphene's unique electronic spectrum and the possibility to tune junction properties by gate voltage. Here we describe graphene junctions with a mean free path of several micrometres, low contact resistance and large supercurrents. Such devices exhibit pronounced Fabry-Pérot oscillations not only… ▽ More Graphene-based Josephson junctions provide a novel platform for studying the proximity effect due to graphene's unique electronic spectrum and the possibility to tune junction properties by gate voltage. Here we describe graphene junctions with a mean free path of several micrometres, low contact resistance and large supercurrents. Such devices exhibit pronounced Fabry-Pérot oscillations not only in the normal-state resistance but also in the critical current. The proximity effect is mostly suppressed in magnetic fields below 10mT, showing the conventional Fraunhofer pattern. Unexpectedly, some proximity survives even in fields higher than 1 T. Superconducting states randomly appear and disappear as a function of field and carrier concentration, and each of them exhibits a supercurrent carrying capacity close to the universal quantum limit. We attribute the high-field Josephson effect to mesoscopic Andreev states that persist near graphene edges. Our work reveals new proximity regimes that can be controlled by quantum confinement and cyclotron motion. △ Less

Submitted 16 December, 2015; v1 submitted 13 April, 2015; originally announced April 2015.

Journal ref: Nature Phys 12, 318-322 (2016)

The Wigner-Witmer diatomic eigenfunction

Authors: James O. Hornkohl, Alexander C. Woods, Christian G. Parigger

Abstract: Born and Oppenheimer reported an approximate separation of molecular eigenfunctions into electronic, vibrational, and rotational parts, but at the end of their paper showed that the two angles describing rotation of the nuclei in a diatomic molecule are exactly separable. A year later in a two-part work devoted strictly to diatomic molecules, Wigner and Witmer published (1) an exact diatomic eigen… ▽ More Born and Oppenheimer reported an approximate separation of molecular eigenfunctions into electronic, vibrational, and rotational parts, but at the end of their paper showed that the two angles describing rotation of the nuclei in a diatomic molecule are exactly separable. A year later in a two-part work devoted strictly to diatomic molecules, Wigner and Witmer published (1) an exact diatomic eigenfunction and (2) the rules correlating the electronic state of a diatomic molecule to the orbital and spin momenta of the separated atoms. The second part of the Wigner-Witmer paper became famous for its correlation rules, but, oddly, the exact eigenfunction from which their rules were obtained received hardly any attention. Using three fundamental symmetries, we give a derivation of the Wigner-Witmer diatomic eigenfunction. Applications of our derivations are fundamental to predicting accurate diatomic molecular spectra that we compare with recorded spectra for diagnostic purposes, such as measurements of molecular spectra following generation of laser-induced plasma. △ Less

Submitted 6 February, 2015; originally announced April 2015.

Designs for generalized linear models with random block effects via information matrix approximations

Authors: Timothy W. Waite, David C. Woods

Abstract: The selection of optimal designs for generalized linear mixed models is complicated by the fact that the Fisher information matrix, on which most optimality criteria depend, is computationally expensive to evaluate. Our focus is on the design of experiments for likelihood estimation of parameters in the conditional model. We provide two novel approximations that substantially reduce the computatio… ▽ More The selection of optimal designs for generalized linear mixed models is complicated by the fact that the Fisher information matrix, on which most optimality criteria depend, is computationally expensive to evaluate. Our focus is on the design of experiments for likelihood estimation of parameters in the conditional model. We provide two novel approximations that substantially reduce the computational cost of evaluating the information matrix by complete enumeration of response outcomes, or Monte Carlo approximations thereof: (i) an asymptotic approximation which is accurate when there is strong dependence between observations in the same block; (ii) an approximation via Kriging interpolators. For logistic random intercept models, we show how interpolation can be especially effective for finding pseudo-Bayesian designs that incorporate uncertainty in the values of the model parameters. The new results are used to provide the first evaluation of the efficiency, for estimating conditional models, of optimal designs from closed-form approximations to the information matrix derived from marginal models. It is found that correcting for the marginal attenuation of parameters in binary-response models yields much improved designs, typically with very high efficiencies. However, in some experiments exhibiting strong dependence, designs for marginal models may still be inefficient for conditional modelling. Our asymptotic results provide some theoretical insights into why such inefficiencies occur. △ Less

Submitted 14 December, 2014; originally announced December 2014.

MSC Class: 62K05; 62K10

Journal ref: Biometrika, 2015, 102, 677-693

Twist-controlled resonant tunnelling in graphene-boron nitride-graphene heterostructures

Authors: A. Mishchenko, J. S. Tu, Y. Cao, R. V. Gorbachev, J. R. Wallbank, M. T. Greenaway, V. E. Morozov, S. V. Morozov, M. J. Zhu, S. L. Wong, F. Withers, C. R. Woods, Y. -J. Kim, K. Watanabe, T. Taniguchi, E. E. Vdovin, O. Makarovsky, T. M. Fromhold, V. I. Falko, A. K. Geim, L. Eaves, K. S. Novoselov

Abstract: Recent developments in the technology of van der Waals heterostructures made from two-dimensional atomic crystals have already led to the observation of new physical phenomena, such as the metal-insulator transition and Coulomb drag, and to the realisation of functional devices, such as tunnel diodes, tunnel transistors and photovoltaic sensors. An unprecedented degree of control of the electronic… ▽ More Recent developments in the technology of van der Waals heterostructures made from two-dimensional atomic crystals have already led to the observation of new physical phenomena, such as the metal-insulator transition and Coulomb drag, and to the realisation of functional devices, such as tunnel diodes, tunnel transistors and photovoltaic sensors. An unprecedented degree of control of the electronic properties is available not only by means of the selection of materials in the stack but also through the additional fine-tuning achievable by adjusting the built-in strain and relative orientation of the component layers. Here we demonstrate how careful alignment of the crystallographic orientation of two graphene electrodes, separated by a layer of hexagonal boron nitride (hBN) in a transistor device, can achieve resonant tunnelling with conservation of electron energy, momentum and, potentially, chirality. We show how the resonance peak and negative differential conductance in the device characteristics induces a tuneable radio-frequency oscillatory current which has potential for future high frequency technology. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: For Supplementary Information see http://www.nature.com/nnano/journal/vaop/ncurrent/extref/nnano.2014.187-s1.pdf

Journal ref: Nature Nanotechnology 9, 808-813 (2014)

Heterostructures produced from nanosheet-based inks

Authors: F. Withers, H. Yang, L. Britnell, A. P Rooney, E. Lewis, A. Felten, C. R. Woods, V. Sanchez Romaguera, T. Georgiou, A. Eckmann, Y. J. Kim, S. G. Yeates, S. J. Haigh, A. K. Geim, K. S. Novoselov, C. Casiraghi

Abstract: The new paradigm of heterostructures based on two-dimensional (2D) atomic crystals has already led to the observation of exciting physical phenomena and creation of novel devices. The possibility of combining layers of different 2D materials in one stack allows unprecedented control over the electronic and optical properties of the resulting material. Still, the current method of mechanical transf… ▽ More The new paradigm of heterostructures based on two-dimensional (2D) atomic crystals has already led to the observation of exciting physical phenomena and creation of novel devices. The possibility of combining layers of different 2D materials in one stack allows unprecedented control over the electronic and optical properties of the resulting material. Still, the current method of mechanical transfer of individual 2D crystals, though allowing exceptional control over the quality of such structures and interfaces, is not scalable. Here we show that such heterostructures can be assembled from chemically exfoliated 2D crystals, allowing for low-cost and scalable methods to be used in the device fabrication. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Journal ref: Nano Letters, 14, 3987 (2014)

Sub-diffractional, volume-confined polaritons in a natural hyperbolic material: hexagonal boron nitride

Authors: Joshua D. Caldwell, Andrey Kretinin, Yiguo Chen, Vincenzo Giannini, Michael M. Fogler, Yan Francescato, Chase T. Ellis, Joseph G. Tischler, Colin R. Woods, Alexander J. Giles, Minghui Hong, Kenji Watanabe, Takashi Taniguchi, Stefan A. Maier, Kostya S. Novoselov

Abstract: Strongly anisotropic media where the principal components of the dielectric tensor have opposite signs are called hyperbolic. Such materials length exhibit unique nanophotonic properties enabled by the highly directional propagation of slow-light modes localized at deeply sub-diffractional scales. While artificial hyperbolic metamaterials have been demonstrated, they suffer from high plasmonic los… ▽ More Strongly anisotropic media where the principal components of the dielectric tensor have opposite signs are called hyperbolic. Such materials length exhibit unique nanophotonic properties enabled by the highly directional propagation of slow-light modes localized at deeply sub-diffractional scales. While artificial hyperbolic metamaterials have been demonstrated, they suffer from high plasmonic losses and require complex nanofabrication, which in turn induces the size-dependent limitations on optical confinement. The low-loss, mid-infrared, natural hyperbolic material, hexagonal boron nitride is an attractive alternative. We observe four series of multiple (up to seven) 'hyperbolic polariton' modes in boron nitride nanocones in two spectral bands. The resonant modes obey the predicted aspect ratio dependence and exhibit record-high quality factors (Q up to 283) in the strong confinement regime (lambda/86 in the smallest structures). These observations assert hexagonal boron nitride as a promising platform for studying novel regimes of light-matter interactions and nanophotonic device engineering. △ Less

Submitted 18 August, 2014; v1 submitted 2 April, 2014; originally announced April 2014.

Journal ref: Nature Communications 5, 5221 (2014)

Electronic Properties of Graphene Encapsulated with Different Two-Dimensional Atomic Crystals

Authors: A. V. Kretinin, Y. Cao, J. S. Tu, G. L. Yu, R. Jalil, K. S. Novoselov, S. J. Haigh, A. Gholinia, A. Mishchenko, M. Lozada, T. Georgiou, C. R. Woods, F. Withers, P. Blake, G. Eda, A. Wirsig, C. Hucho, K. Watanabe, T. Taniguchi, A. K. Geim, R. V. Gorbachev

Abstract: Hexagonal boron nitride is the only substrate that has so far allowed graphene devices exhibiting micron-scale ballistic transport. Can other atomically flat crystals be used as substrates for making quality graphene heterostructures? Here we report on our search for alternative substrates. The devices fabricated by encapsulating graphene with molybdenum or tungsten disulphides and hBN are found t… ▽ More Hexagonal boron nitride is the only substrate that has so far allowed graphene devices exhibiting micron-scale ballistic transport. Can other atomically flat crystals be used as substrates for making quality graphene heterostructures? Here we report on our search for alternative substrates. The devices fabricated by encapsulating graphene with molybdenum or tungsten disulphides and hBN are found to exhibit consistently high carrier mobilities of about 60,000 cm$^{2}$V$^{-1}$s$^{-1}$. In contrast, encapsulation with atomically flat layered oxides such as mica, bismuth strontium calcium copper oxide and vanadium pentoxide results in exceptionally low quality of graphene devices with mobilities of ~ 1,000 cm$^{2}$ V$^{-1}$s$^{-1}$. We attribute the difference mainly to self-cleansing that takes place at interfaces between graphene, hBN and transition metal dichalcogenides. Surface contamination assembles into large pockets allowing the rest of the interface to become atomically clean. The cleansing process does not occur for graphene on atomically flat oxide substrates. △ Less

Submitted 24 May, 2014; v1 submitted 20 March, 2014; originally announced March 2014.

Comments: 19 pages, 11 figures, 1 table including Supporting Information

Journal ref: Nano Letters 14, 3270 - 3276 (2014)

Commensurate-incommensurate transition for graphene on hexagonal boron nitride

Authors: C. R. Woods, L. Britnell, A. Eckmann, R. S. Ma, J. C. Lu, H. M. Guo, X. Lin, G. L. Yu, Y. Cao, R. V. Gorbachev, A. V. Kretinin, J. Park, L. A. Ponomarenko, M. I. Katsnelson, Yu. N. Gornostyrev, K. Watanabe, T. Taniguchi, C. Casiraghi, H. J. Gao, A. K. Geim, K. S. Novoselov

Abstract: When a crystal is subjected to a periodic potential, under certain circumstances (such as when the period of the potential is close to the crystal periodicity; the potential is strong enough, etc.) it might adjust itself to follow the periodicity of the potential, resulting in a, so called, commensurate state. Such commensurate-incommensurate transitions are ubiquitous phenomena in many areas of c… ▽ More When a crystal is subjected to a periodic potential, under certain circumstances (such as when the period of the potential is close to the crystal periodicity; the potential is strong enough, etc.) it might adjust itself to follow the periodicity of the potential, resulting in a, so called, commensurate state. Such commensurate-incommensurate transitions are ubiquitous phenomena in many areas of condensed matter physics: from magnetism and dislocations in crystals, to vortices in superconductors, and atomic layers adsorbed on a crystalline surface. Of particular interest might be the properties of topological defects between the two commensurate phases: solitons, domain walls, and dislocation walls. Here we report a commensurate-incommensurate transition for graphene on top of hexagonal boron nitride (hBN). Depending on the rotational angle between the two hexagonal lattices, graphene can either stretch to adjust to a slightly different hBN periodicity (the commensurate state found for small rotational angles) or exhibit little adjustment (the incommensurate state). In the commensurate state, areas with matching lattice constants are separated by domain walls that accumulate the resulting strain. Such soliton-like objects present significant fundamental interest, and their presence might explain recent observations when the electronic, optical, Raman and other properties of graphene-hBN heterostructures have been notably altered. △ Less

Submitted 20 March, 2014; v1 submitted 12 January, 2014; originally announced January 2014.

Journal ref: Nature Physics 10(6), 451 - 456 (2014)