-
Pick-and-place transfer of arbitrary-metal electrodes for van der Waals device fabrication
Authors:
Kaijian Xing,
Daniel McEwen,
Weiyao Zhao,
Abdulhakim Bake,
David Cortie,
Jingying Liu,
Thi-Hai-Yen Vu,
James Hone,
Alastair Stacey,
Mark T. Edmonds,
Kenji Watanabe,
Takashi Taniguchi,
Qingdong Ou,
Dong-Chen Qi,
Michael S. Fuhrer
Abstract:
Van der Waals electrode integration is a promising strategy to create near-perfect interfaces between metals and two-dimensional materials, with advantages such as eliminating Fermi-level pinning and reducing contact resistance. However, the lack of a simple, generalizable pick-and-place transfer technology has greatly hampered the wide use of this technique. We demonstrate the pick-and-place tran…
▽ More
Van der Waals electrode integration is a promising strategy to create near-perfect interfaces between metals and two-dimensional materials, with advantages such as eliminating Fermi-level pinning and reducing contact resistance. However, the lack of a simple, generalizable pick-and-place transfer technology has greatly hampered the wide use of this technique. We demonstrate the pick-and-place transfer of pre-fabricated electrodes from reusable polished hydrogenated diamond substrates without the use of any surface treatments or sacrificial layers. The technique enables transfer of large-scale arbitrary metal electrodes, as demonstrated by successful transfer of eight different elemental metals with work functions ranging from 4.22 to 5.65 eV. The mechanical transfer of metal electrodes from diamond onto van der Waals materials creates atomically smooth interfaces with no interstitial impurities or disorder, as observed with cross-sectional high-resolution transmission electron microscopy and energy-dispersive X-ray spectroscopy. As a demonstration of its device application, we use the diamond-transfer technique to create metal contacts to monolayer transition metal dichalcogenide semiconductors with high-work-function Pd, low-work-function Ti, and semi metal Bi to create n- and p-type field-effect transistors with low Schottky barrier heights. We also extend this technology to other applications such as ambipolar transistor and optoelectronics, paving the way for new device architectures and high-performance devices.
△ Less
Submitted 21 May, 2024;
originally announced May 2024.
-
Differentiable and accelerated wavelet transforms on the sphere and ball
Authors:
Matthew A. Price,
Alicja Polanska,
Jessica Whitney,
Jason D. McEwen
Abstract:
Directional wavelet dictionaries are hierarchical representations which efficiently capture and segment information across scale, location and orientation. Such representations demonstrate a particular affinity to physical signals, which often exhibit highly anisotropic, localised multiscale structure. Many physically important signals are observed over spherical domains, such as the celestial sky…
▽ More
Directional wavelet dictionaries are hierarchical representations which efficiently capture and segment information across scale, location and orientation. Such representations demonstrate a particular affinity to physical signals, which often exhibit highly anisotropic, localised multiscale structure. Many physically important signals are observed over spherical domains, such as the celestial sky in cosmology. Leveraging recent advances in computational harmonic analysis, we design new highly distributable and automatically differentiable directional wavelet transforms on the $2$-dimensional sphere $\mathbb{S}^2$ and $3$-dimensional ball $\mathbb{B}^3 = \mathbb{R}^+ \times \mathbb{S}^2$ (the space formed by augmenting the sphere with the radial half-line). We observe up to a $300$-fold and $21800$-fold acceleration for signals on the sphere and ball, respectively, compared to existing software, whilst maintaining 64-bit machine precision. Not only do these algorithms dramatically accelerate existing spherical wavelet transforms, the gradient information afforded by automatic differentiation unlocks many data-driven analysis techniques previously not possible for these spaces. We publicly release both S2WAV and S2BALL, open-sourced JAX libraries for our transforms that are automatically differentiable and readily deployable both on and over clusters of hardware accelerators (e.g. GPUs & TPUs).
△ Less
Submitted 14 March, 2024; v1 submitted 2 February, 2024;
originally announced February 2024.
-
Differentiable and accelerated spherical harmonic and Wigner transforms
Authors:
Matthew A. Price,
Jason D. McEwen
Abstract:
Many areas of science and engineering encounter data defined on spherical manifolds. Modelling and analysis of spherical data often necessitates spherical harmonic transforms, at high degrees, and increasingly requires efficient computation of gradients for machine learning or other differentiable programming tasks. We develop novel algorithmic structures for accelerated and differentiable computa…
▽ More
Many areas of science and engineering encounter data defined on spherical manifolds. Modelling and analysis of spherical data often necessitates spherical harmonic transforms, at high degrees, and increasingly requires efficient computation of gradients for machine learning or other differentiable programming tasks. We develop novel algorithmic structures for accelerated and differentiable computation of generalised Fourier transforms on the sphere $\mathbb{S}^2$ and rotation group $\text{SO}(3)$, i.e. spherical harmonic and Wigner transforms, respectively. We present a recursive algorithm for the calculation of Wigner $d$-functions that is both stable to high harmonic degrees and extremely parallelisable. By tightly coupling this with separable spherical transforms, we obtain algorithms that exhibit an extremely parallelisable structure that is well-suited for the high throughput computing of modern hardware accelerators (e.g. GPUs). We also develop a hybrid automatic and manual differentiation approach so that gradients can be computed efficiently. Our algorithms are implemented within the JAX differentiable programming framework in the S2FFT software code. Numerous samplings of the sphere are supported, including equiangular and HEALPix sampling. Computational errors are at the order of machine precision for spherical samplings that admit a sampling theorem. When benchmarked against alternative C codes we observe up to a 400-fold acceleration. Furthermore, when distributing over multiple GPUs we achieve very close to optimal linear scaling with increasing number of GPUs due to the highly parallelised and balanced nature of our algorithms. Provided access to sufficiently many GPUs our transforms thus exhibit an unprecedented effective linear time complexity.
△ Less
Submitted 20 May, 2024; v1 submitted 24 November, 2023;
originally announced November 2023.
-
Bayesian model comparison for simulation-based inference
Authors:
A. Spurio Mancini,
M. M. Docherty,
M. A. Price,
J. D. McEwen
Abstract:
Comparison of appropriate models to describe observational data is a fundamental task of science. The Bayesian model evidence, or marginal likelihood, is a computationally challenging, yet crucial, quantity to estimate to perform Bayesian model comparison. We introduce a methodology to compute the Bayesian model evidence in simulation-based inference (SBI) scenarios (also often called likelihood-f…
▽ More
Comparison of appropriate models to describe observational data is a fundamental task of science. The Bayesian model evidence, or marginal likelihood, is a computationally challenging, yet crucial, quantity to estimate to perform Bayesian model comparison. We introduce a methodology to compute the Bayesian model evidence in simulation-based inference (SBI) scenarios (also often called likelihood-free inference). In particular, we leverage the recently proposed learnt harmonic mean estimator and exploit the fact that it is decoupled from the method used to generate posterior samples, i.e. it requires posterior samples only, which may be generated by any approach. This flexibility, which is lacking in many alternative methods for computing the model evidence, allows us to develop SBI model comparison techniques for the three main neural density estimation approaches, including neural posterior estimation (NPE), neural likelihood estimation (NLE), and neural ratio estimation (NRE). We demonstrate and validate our SBI evidence calculation techniques on a range of inference problems, including a gravitational wave example. Moreover, we further validate the accuracy of the learnt harmonic mean estimator, implemented in the HARMONIC software, in likelihood-based settings. These results highlight the potential of HARMONIC as a sampler-agnostic method to estimate the model evidence in both likelihood-based and simulation-based scenarios.
△ Less
Submitted 8 November, 2023; v1 submitted 8 July, 2022;
originally announced July 2022.
-
Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
Authors:
Augustin Marignier,
Jason D. McEwen,
Ana M. G. Ferreira,
Thomas D. Kitching
Abstract:
Inverse problems defined on the sphere arise in many fields, including seismology and cosmology where problems are defined on the globe and the cosmic sphere. These are generally high-dimensional and computationally very complex and, as a result, sampling the posterior of spherical inverse problems is a challenging task. In this work, we describe a framework that leverages a proximal Markov chain…
▽ More
Inverse problems defined on the sphere arise in many fields, including seismology and cosmology where problems are defined on the globe and the cosmic sphere. These are generally high-dimensional and computationally very complex and, as a result, sampling the posterior of spherical inverse problems is a challenging task. In this work, we describe a framework that leverages a proximal Markov chain Monte Carlo (MCMC) algorithm to efficiently sample the high-dimensional space of spherical inverse problems with a sparsity-promoting wavelet prior. We detail the modifications needed for the algorithm to be applied to spherical problems, and give special consideration to the crucial forward modelling step which contains spherical harmonic transforms that are computationally expensive. By sampling the posterior, our framework allows for full and flexible uncertainty quantification, something which is not possible with other methods based on, for example, convex optimisation. We demonstrate our framework in practice on full-sky cosmological mass-mapping and on a common problem in global seismic tomography. We find that our approach is potentially useful at moderate resolutions, such as those of interest in seismology. Our framework is generally limited by resolution requirements, such as those required for astrophysical applications, due to the poor scaling of the complexity of spherical harmonic transforms with resolution. A new Python package, pxmcmc, containing the proximal MCMC sampler, measurement operators, wavelet transforms and sparse priors is made publicly available.
△ Less
Submitted 18 November, 2022; v1 submitted 14 July, 2021;
originally announced July 2021.
-
Multiscale Optimal Filtering on the Sphere
Authors:
Adeem Aslam,
Zubair Khalid,
Jason D. McEwen
Abstract:
We present a framework for the optimal filtering of spherical signals contaminated by realizations of an additive, zero-mean, uncorrelated and anisotropic noise process on the sphere. Filtering is performed in the wavelet domain given by the scale-discretized wavelet transform on the sphere. The proposed filter is optimal in the sense that it minimizes the mean square error between the filtered wa…
▽ More
We present a framework for the optimal filtering of spherical signals contaminated by realizations of an additive, zero-mean, uncorrelated and anisotropic noise process on the sphere. Filtering is performed in the wavelet domain given by the scale-discretized wavelet transform on the sphere. The proposed filter is optimal in the sense that it minimizes the mean square error between the filtered wavelet representation and wavelet representation of the noise-free signal. We also present a simplified formulation of the filter for the case when azimuthally symmetric wavelet functions are used. We demonstrate the use of the proposed optimal filter for denoising of an Earth topography map in the presence of additive, zero-mean, uncorrelated and white Gaussian noise, and show that the proposed filter performs better than the hard thresholding method and weighted spherical harmonic~(weighted-SPHARM) signal estimation framework.
△ Less
Submitted 6 February, 2021; v1 submitted 15 October, 2020;
originally announced October 2020.
-
Harmonic analysis of spherical sampling in diffusion MRI
Authors:
A. Daducci,
J. D. McEwen,
D. Van De Ville,
J. -Ph. Thiran,
Y. Wiaux
Abstract:
In the last decade diffusion MRI has become a powerful tool to non-invasively study white-matter integrity in the brain. Recently many research groups have focused their attention on multi-shell spherical acquisitions with the aim of effectively mapping the diffusion signal with a lower number of q-space samples, hence enabling a crucial reduction of acquisition time. One of the quantities commonl…
▽ More
In the last decade diffusion MRI has become a powerful tool to non-invasively study white-matter integrity in the brain. Recently many research groups have focused their attention on multi-shell spherical acquisitions with the aim of effectively mapping the diffusion signal with a lower number of q-space samples, hence enabling a crucial reduction of acquisition time. One of the quantities commonly studied in this context is the so-called orientation distribution function (ODF). In this setting, the spherical harmonic (SH) transform has gained a great deal of popularity thanks to its ability to perform convolution operations efficiently and accurately, such as the Funk-Radon transform notably required for ODF computation from q-space data. However, if the q-space signal is described with an unsuitable angular resolution at any b-value probed, aliasing (or interpolation) artifacts are unavoidably created. So far this aspect has been tackled empirically and, to our knowledge, no study has addressed this problem in a quantitative approach. The aim of the present work is to study more theoretically the efficiency of multi-shell spherical sampling in diffusion MRI, in order to gain understanding in HYDI-like approaches, possibly paving the way to further optimization strategies.
△ Less
Submitted 1 June, 2011;
originally announced June 2011.