Showing 1–50 of 118 results for author: Armstrong, T

ALOHA 2: An Enhanced Low-Cost Hardware for Bimanual Teleoperation

Authors: ALOHA 2 Team, Jorge Aldaco, Travis Armstrong, Robert Baruch, Jeff Bingham, Sanky Chan, Kenneth Draper, Debidatta Dwibedi, Chelsea Finn, Pete Florence, Spencer Goodrich, Wayne Gramlich, Torr Hage, Alexander Herzog, Jonathan Hoech, Thinh Nguyen, Ian Storz, Baruch Tabanpour, Leila Takayama, Jonathan Tompson, Ayzaan Wahid, Ted Wahrburg, Sichun Xu, Sergey Yaroshenko, Kevin Zakka , et al. (1 additional authors not shown)

Abstract: Diverse demonstration datasets have powered significant advances in robot learning, but the dexterity and scale of such data can be limited by the hardware cost, the hardware robustness, and the ease of teleoperation. We introduce ALOHA 2, an enhanced version of ALOHA that has greater performance, ergonomics, and robustness compared to the original design. To accelerate research in large-scale bim… ▽ More Diverse demonstration datasets have powered significant advances in robot learning, but the dexterity and scale of such data can be limited by the hardware cost, the hardware robustness, and the ease of teleoperation. We introduce ALOHA 2, an enhanced version of ALOHA that has greater performance, ergonomics, and robustness compared to the original design. To accelerate research in large-scale bimanual manipulation, we open source all hardware designs of ALOHA 2 with a detailed tutorial, together with a MuJoCo model of ALOHA 2 with system identification. See the project website at aloha-2.github.io. △ Less

Submitted 7 February, 2024; originally announced May 2024.

Comments: Project website: aloha-2.github.io

Open X-Embodiment: Robotic Learning Datasets and RT-X Models

Authors: Open X-Embodiment Collaboration, Abby O'Neill, Abdul Rehman, Abhinav Gupta, Abhiram Maddukuri, Abhishek Gupta, Abhishek Padalkar, Abraham Lee, Acorn Pooley, Agrim Gupta, Ajay Mandlekar, Ajinkya Jain, Albert Tung, Alex Bewley, Alex Herzog, Alex Irpan, Alexander Khazatsky, Anant Rai, Anchit Gupta, Andrew Wang, Andrey Kolobov, Anikait Singh, Animesh Garg, Aniruddha Kembhavi, Annie Xie , et al. (267 additional authors not shown)

Abstract: Large, high-capacity models trained on diverse datasets have shown remarkable successes on efficiently tackling downstream applications. In domains from NLP to Computer Vision, this has led to a consolidation of pretrained models, with general pretrained backbones serving as a starting point for many applications. Can such a consolidation happen in robotics? Conventionally, robotic learning method… ▽ More Large, high-capacity models trained on diverse datasets have shown remarkable successes on efficiently tackling downstream applications. In domains from NLP to Computer Vision, this has led to a consolidation of pretrained models, with general pretrained backbones serving as a starting point for many applications. Can such a consolidation happen in robotics? Conventionally, robotic learning methods train a separate model for every application, every robot, and even every environment. Can we instead train generalist X-robot policy that can be adapted efficiently to new robots, tasks, and environments? In this paper, we provide datasets in standardized data formats and models to make it possible to explore this possibility in the context of robotic manipulation, alongside experimental results that provide an example of effective X-robot policies. We assemble a dataset from 22 different robots collected through a collaboration between 21 institutions, demonstrating 527 skills (160266 tasks). We show that a high-capacity model trained on this data, which we call RT-X, exhibits positive transfer and improves the capabilities of multiple robots by leveraging experience from other platforms. More details can be found on the project website https://robotics-transformer-x.github.io. △ Less

Submitted 1 June, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Project website: https://robotics-transformer-x.github.io

The impact of wettability on the co-moving velocity of two-fluid flow in porous media

Authors: Fatimah Alzubaidi, James E. McClure, Håkon Pedersen, Alex Hansen, Carl Fredrik Berg, Peyman Mostaghimi, Ryan T. Armstrong

Abstract: The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. The simulations w… ▽ More The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. The simulations were conducted using the Lattice-Boltzmann color-fluid model and evaluated the relative permeability for different wetting conditions on the same rock. The analysis of the simulation data followed the methodology proposed by Roy et al. (2022) to reconstruct a constitutive equation for the co-moving velocity. Surprisingly, it was found that the coefficients of the constitutive equation were nearly the same for all wetting conditions. Based on these results, a simple approach was proposed to reconstruct the oil phase relative permeability using only the co-moving velocity relationship and water phase relative permeability. This proposed method provides new insights into the dependency of relative permeability curves, which has implications for the history matching of production data and solving the associated inverse problem. The research findings contribute to a better understanding of the impact of wettability on fluid flow in porous media and provide a practical approach for estimating relative permeability based on the co-moving velocity relationship, which has never been shown before. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 14 pages, 6 figures

Reviving Shift Equivariance in Vision Transformers

Authors: Peijian Ding, Davit Soselia, Thomas Armstrong, Jiahao Su, Furong Huang

Abstract: Shift equivariance is a fundamental principle that governs how we perceive the world - our recognition of an object remains invariant with respect to shifts. Transformers have gained immense popularity due to their effectiveness in both language and vision tasks. While the self-attention operator in vision transformers (ViT) is permutation-equivariant and thus shift-equivariant, patch embedding, p… ▽ More Shift equivariance is a fundamental principle that governs how we perceive the world - our recognition of an object remains invariant with respect to shifts. Transformers have gained immense popularity due to their effectiveness in both language and vision tasks. While the self-attention operator in vision transformers (ViT) is permutation-equivariant and thus shift-equivariant, patch embedding, positional encoding, and subsampled attention in ViT variants can disrupt this property, resulting in inconsistent predictions even under small shift perturbations. Although there is a growing trend in incorporating the inductive bias of convolutional neural networks (CNNs) into vision transformers, it does not fully address the issue. We propose an adaptive polyphase anchoring algorithm that can be seamlessly integrated into vision transformer models to ensure shift-equivariance in patch embedding and subsampled attention modules, such as window attention and global subsampled attention. Furthermore, we utilize depth-wise convolution to encode positional information. Our algorithms enable ViT, and its variants such as Twins to achieve 100% consistency with respect to input shift, demonstrate robustness to cropping, flipping, and affine transformations, and maintain consistent predictions even when the original models lose 20 percentage points on average when shifted by just a few pixels with Twins' accuracy dropping from 80.57% to 62.40%. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 9 pages, 3 figures

Adapting to Misspecification

Authors: Timothy B. Armstrong, Patrick Kline, Liyang Sun

Abstract: Empirical research typically involves a robustness-efficiency tradeoff. A researcher seeking to estimate a scalar parameter can invoke strong assumptions to motivate a restricted estimator that is precise but may be heavily biased, or they can relax some of these assumptions to motivate a more robust, but variable, unrestricted estimator. When a bound on the bias of the restricted estimator is ava… ▽ More Empirical research typically involves a robustness-efficiency tradeoff. A researcher seeking to estimate a scalar parameter can invoke strong assumptions to motivate a restricted estimator that is precise but may be heavily biased, or they can relax some of these assumptions to motivate a more robust, but variable, unrestricted estimator. When a bound on the bias of the restricted estimator is available, it is optimal to shrink the unrestricted estimator towards the restricted estimator. For settings where a bound on the bias of the restricted estimator is unknown, we propose adaptive shrinkage estimators that minimize the percentage increase in worst case risk relative to an oracle that knows the bound. We show that adaptive estimators solve a weighted convex minimax problem and provide lookup tables facilitating their rapid computation. Revisiting five empirical studies where questions of model specification arise, we examine the advantages of adapting to -- rather than testing for -- misspecification. △ Less

Submitted 11 July, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: 69 pages, 12 figures

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

DeepAngle: Fast calculation of contact angles in tomography images using deep learning

Authors: Arash Rabbani, Chenhao Sun, Masoud Babaei, Vahid J. Niasar, Ryan T. Armstrong, Peyman Mostaghimi

Abstract: DeepAngle is a machine learning-based method to determine the contact angles of different phases in the tomography images of porous materials. Measurement of angles in 3--D needs to be done within the surface perpendicular to the angle planes, and it could become inaccurate when dealing with the discretized space of the image voxels. A computationally intensive solution is to correlate and vectori… ▽ More DeepAngle is a machine learning-based method to determine the contact angles of different phases in the tomography images of porous materials. Measurement of angles in 3--D needs to be done within the surface perpendicular to the angle planes, and it could become inaccurate when dealing with the discretized space of the image voxels. A computationally intensive solution is to correlate and vectorize all surfaces using an adaptable grid, and then measure the angles within the desired planes. On the contrary, the present study provides a rapid and low-cost technique powered by deep learning to estimate the interfacial angles directly from images. DeepAngle is tested on both synthetic and realistic images against the direct measurement technique and found to improve the r-squared by 5 to 16% while lowering the computational cost 20 times. This rapid method is especially applicable for processing large tomography data and time-resolved images, which is computationally intensive. The developed code and the dataset are available at an open repository on GitHub (https://www.github.com/ArashRabbani/DeepAngle). △ Less

Submitted 28 November, 2022; originally announced November 2022.

Angular Diameters and Fundamental Parameters of Forty-Four Stars from the Navy Precision Optical Interferometer

Authors: Ellyn K. Baines, J. Thomas Armstrong, James H. Clark III, Jim Gorney, Donald J. Hutter, Anders M. Jorgensen, Casey Kyte, David Mozurkewich, Ishara Nisley, Jason Sanborn, Henrique R. Schmitt, Gerard T. van Belle

Abstract: We measured the angular diameters of 44 stars with the Navy Precision Optical Interferometer, obtaining uncertainties on the limb darkened diameter of 2% or less for all but four stars. We then used our diameters with Gaia or Hipparcos parallaxes to calculate each star's physical radius. We gathered information from the literature to determine bolometric flux and luminosity, and combined that with… ▽ More We measured the angular diameters of 44 stars with the Navy Precision Optical Interferometer, obtaining uncertainties on the limb darkened diameter of 2% or less for all but four stars. We then used our diameters with Gaia or Hipparcos parallaxes to calculate each star's physical radius. We gathered information from the literature to determine bolometric flux and luminosity, and combined that with our diameters to produce an effective temperature. Our sample consists of mostly giant stars, and spans a wide range of spectral classes from B to M. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 13 pages, 3 figures, 6 tables. arXiv admin note: substantial text overlap with arXiv:1712.08109

Journal ref: 2021AJ....162..198B

Robust Estimation and Inference in Panels with Interactive Fixed Effects

Authors: Timothy B. Armstrong, Martin Weidner, Andrei Zeleneev

Abstract: We consider estimation and inference for a regression coefficient in panels with interactive fixed effects (i.e., with a factor structure). We show that previously developed estimators and confidence intervals (CIs) might be heavily biased and size-distorted when some of the factors are weak. We propose estimators with improved rates of convergence and bias-aware CIs that are uniformly valid regar… ▽ More We consider estimation and inference for a regression coefficient in panels with interactive fixed effects (i.e., with a factor structure). We show that previously developed estimators and confidence intervals (CIs) might be heavily biased and size-distorted when some of the factors are weak. We propose estimators with improved rates of convergence and bias-aware CIs that are uniformly valid regardless of whether the factors are strong or not. Our approach applies the theory of minimax linear estimation to form a debiased estimate using a nuclear norm bound on the error of an initial estimate of the interactive fixed effects. We use the obtained estimate to construct a bias-aware CI taking into account the remaining bias due to weak factors. In Monte Carlo experiments, we find a substantial improvement over conventional approaches when factors are weak, with little cost to estimation error when factors are strong. △ Less

Submitted 4 July, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Interactive Language: Talking to Robots in Real Time

Authors: Corey Lynch, Ayzaan Wahid, Jonathan Tompson, Tianli Ding, James Betker, Robert Baruch, Travis Armstrong, Pete Florence

Abstract: We present a framework for building interactive, real-time, natural language-instructable robots in the real world, and we open source related assets (dataset, environment, benchmark, and policies). Trained with behavioral cloning on a dataset of hundreds of thousands of language-annotated trajectories, a produced policy can proficiently execute an order of magnitude more commands than previous wo… ▽ More We present a framework for building interactive, real-time, natural language-instructable robots in the real world, and we open source related assets (dataset, environment, benchmark, and policies). Trained with behavioral cloning on a dataset of hundreds of thousands of language-annotated trajectories, a produced policy can proficiently execute an order of magnitude more commands than previous works: specifically we estimate a 93.5% success rate on a set of 87,000 unique natural language strings specifying raw end-to-end visuo-linguo-motor skills in the real world. We find that the same policy is capable of being guided by a human via real-time language to address a wide range of precise long-horizon rearrangement goals, e.g. "make a smiley face out of blocks". The dataset we release comprises nearly 600,000 language-labeled trajectories, an order of magnitude larger than prior available datasets. We hope the demonstrated results and associated assets enable further advancement of helpful, capable, natural-language-interactable robots. See videos at https://interactive-language.github.io. △ Less

Submitted 12 October, 2022; originally announced October 2022.

False Discovery Rate Adjustments for Average Significance Level Controlling Tests

Authors: Timothy B. Armstrong

Abstract: Multiple testing adjustments, such as the Benjamini and Hochberg (1995) step-up procedure for controlling the false discovery rate (FDR), are typically applied to families of tests that control significance level in the classical sense: for each individual test, the probability of false rejection is no greater than the nominal level. In this paper, we consider tests that satisfy only a weaker noti… ▽ More Multiple testing adjustments, such as the Benjamini and Hochberg (1995) step-up procedure for controlling the false discovery rate (FDR), are typically applied to families of tests that control significance level in the classical sense: for each individual test, the probability of false rejection is no greater than the nominal level. In this paper, we consider tests that satisfy only a weaker notion of significance level control, in which the probability of false rejection need only be controlled on average over the hypotheses. We find that the Benjamini and Hochberg (1995) step-up procedure still controls FDR in the asymptotic regime with many weakly dependent $p$-values, and that certain adjustments for dependent $p$-values such as the Benjamini and Yekutieli (2001) procedure continue to yield FDR control in finite samples. Our results open the door to FDR controlling procedures in nonparametric and high dimensional settings where weakening the notion of inference allows for large power improvements. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Authors: Fatimah Alzubaidi, Peyman Mostaghimi, Yufu Niu, Ryan T. Armstrong, Gelareh Mohammadi, Steffen Berg, James E. McClure

Abstract: Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process/path dependent with an underlying dependency on pore structure. For applications, fuel cells to underground $CO_2$ storage, it is imperative to determine the effective phase permeability relationships where the traditional approach is based on t… ▽ More Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process/path dependent with an underlying dependency on pore structure. For applications, fuel cells to underground $CO_2$ storage, it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modelling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an Artificial Neural Network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. The developed ANN is based on a prescribed set of state variables based on physical insights that predicts the effective permeability of 4,500 unseen pore-scale geometrical states with $R^2 = 0.98$. △ Less

Submitted 22 June, 2023; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: 6 Pages, 2 Figures, and Supporting Material

MSC Class: 76S05 (Primary); 76T06 (Secondary) ACM Class: J.2

Relative permeability as a stationary process: energy fluctuations in immiscible displacement

Authors: James E. McClure, Ming Fan, Steffen Berg, Ryan T. Armstrong, Carl Fredrik Berg, Zhe Li, Thomas Ramstad

Abstract: Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) the… ▽ More Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ($\mathsf{Ca}\sim1\times10^5$), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways. △ Less

Submitted 2 July, 2022; originally announced July 2022.

Visuomotor Control in Multi-Object Scenes Using Object-Aware Representations

Authors: Negin Heravi, Ayzaan Wahid, Corey Lynch, Pete Florence, Travis Armstrong, Jonathan Tompson, Pierre Sermanet, Jeannette Bohg, Debidatta Dwibedi

Abstract: Perceptual understanding of the scene and the relationship between its different components is important for successful completion of robotic tasks. Representation learning has been shown to be a powerful technique for this, but most of the current methodologies learn task specific representations that do not necessarily transfer well to other tasks. Furthermore, representations learned by supervi… ▽ More Perceptual understanding of the scene and the relationship between its different components is important for successful completion of robotic tasks. Representation learning has been shown to be a powerful technique for this, but most of the current methodologies learn task specific representations that do not necessarily transfer well to other tasks. Furthermore, representations learned by supervised methods require large labeled datasets for each task that are expensive to collect in the real world. Using self-supervised learning to obtain representations from unlabeled data can mitigate this problem. However, current self-supervised representation learning methods are mostly object agnostic, and we demonstrate that the resulting representations are insufficient for general purpose robotics tasks as they fail to capture the complexity of scenes with many components. In this paper, we explore the effectiveness of using object-aware representation learning techniques for robotic tasks. Our self-supervised representations are learned by observing the agent freely interacting with different parts of the environment and is queried in two different settings: (i) policy learning and (ii) object location prediction. We show that our model learns control policies in a sample-efficient manner and outperforms state-of-the-art object agnostic techniques as well as methods trained on raw RGB images. Our results show a 20 percent increase in performance in low data regimes (1000 trajectories) in policy training using implicit behavioral cloning (IBC). Furthermore, our method outperforms the baselines for the task of object localization in multi-object scenes. △ Less

Submitted 12 March, 2023; v1 submitted 12 May, 2022; originally announced May 2022.

Asymptotic Efficiency Bounds for a Class of Experimental Designs

Authors: Timothy B. Armstrong

Abstract: We consider an experimental design setting in which units are assigned to treatment after being sampled sequentially from an infinite population. We derive asymptotic efficiency bounds that apply to data from any experiment that assigns treatment as a (possibly randomized) function of covariates and past outcome data, including stratification on covariates and adaptive designs. For estimating the… ▽ More We consider an experimental design setting in which units are assigned to treatment after being sampled sequentially from an infinite population. We derive asymptotic efficiency bounds that apply to data from any experiment that assigns treatment as a (possibly randomized) function of covariates and past outcome data, including stratification on covariates and adaptive designs. For estimating the average treatment effect of a binary treatment, our results show that no further first order asymptotic efficiency improvement is possible relative to an estimator that achieves the Hahn (1998) bound in an experimental design where the propensity score is chosen to minimize this bound. Our results also apply to settings with multiple treatments with possible constraints on treatment, as well as covariate based sampling of a single outcome. △ Less

Submitted 5 May, 2022; originally announced May 2022.

A comparative study of paired versus unpaired deep learning methods for physically enhancing digital rock image resolution

Authors: Yufu Niu, Samuel J. Jackson, Naif Alqahtani, Peyman Mostaghimi, Ryan T. Armstrong

Abstract: X-ray micro-computed tomography (micro-CT) has been widely leveraged to characterise pore-scale geometry in subsurface porous rock. Recent developments in super resolution (SR) methods using deep learning allow the digital enhancement of low resolution (LR) images over large spatial scales, creating SR images comparable to the high resolution (HR) ground truth. This circumvents traditional resolut… ▽ More X-ray micro-computed tomography (micro-CT) has been widely leveraged to characterise pore-scale geometry in subsurface porous rock. Recent developments in super resolution (SR) methods using deep learning allow the digital enhancement of low resolution (LR) images over large spatial scales, creating SR images comparable to the high resolution (HR) ground truth. This circumvents traditional resolution and field-of-view trade-offs. An outstanding issue is the use of paired (registered) LR and HR data, which is often required in the training step of such methods but is difficult to obtain. In this work, we rigorously compare two different state-of-the-art SR deep learning techniques, using both paired and unpaired data, with like-for-like ground truth data. The first approach requires paired images to train a convolutional neural network (CNN) while the second approach uses unpaired images to train a generative adversarial network (GAN). The two approaches are compared using a micro-CT carbonate rock sample with complicated micro-porous textures. We implemented various image based and numerical verifications and experimental validation to quantitatively evaluate the physical accuracy and sensitivities of the two methods. Our quantitative results show that unpaired GAN approach can reconstruct super-resolution images as precise as paired CNN method, with comparable training times and dataset requirement. This unlocks new applications for micro-CT image enhancement using unpaired deep learning methods; image registration is no longer needed during the data processing stage. Decoupled images from data storage platforms can be exploited more efficiently to train networks for SR digital rock applications. This opens up a new pathway for various applications of multi-scale flow simulation in heterogeneous porous media. △ Less

Submitted 16 December, 2021; originally announced December 2021.

Comments: 26 pages, 11 figures, 4 tables

Deep learning of multi-resolution X-Ray micro-CT images for multi-scale modelling

Authors: Samuel J. Jackson, Yufu Niu, Sojwal Manoorkar, Peyman Mostaghimi, Ryan T. Armstrong

Abstract: Field-of-view and resolution trade-offs in X-Ray micro-computed tomography (micro-CT) imaging limit the characterization, analysis and model development of multi-scale porous systems. To this end, we developed an applied methodology utilising deep learning to enhance low resolution images over large sample sizes and create multi-scale models capable of accurately simulating experimental fluid dyna… ▽ More Field-of-view and resolution trade-offs in X-Ray micro-computed tomography (micro-CT) imaging limit the characterization, analysis and model development of multi-scale porous systems. To this end, we developed an applied methodology utilising deep learning to enhance low resolution images over large sample sizes and create multi-scale models capable of accurately simulating experimental fluid dynamics from the pore (microns) to continuum (centimetres) scale. We develop a 3D Enhanced Deep Super Resolution (EDSR) convolutional neural network to create super resolution (SR) images from low resolution images, which alleviates common micro-CT hardware/reconstruction defects in high-resolution (HR) images. When paired with pore-network simulations and parallel computation, we can create large 3D continuum-scale models with spatially varying flow & material properties. We quantitatively validate the workflow at various scales using direct HR/SR image similarity, pore-scale material/flow simulations and continuum scale multiphase flow experiments (drainage immiscible flow pressures and 3D fluid volume fractions). The SR images and models are comparable to the HR ground truth, and generally accurate to within experimental uncertainty at the continuum scale across a range of flow rates. They are found to be significantly more accurate than their LR counterparts, especially in cases where a wide distribution of pore-sizes are encountered. The applied methodology opens up the possibility to image, model and analyse truly multi-scale heterogeneous systems that are otherwise intractable. △ Less

Submitted 15 March, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 21 pages, 9 figures

MSC Class: 76S05 ACM Class: I.4.3

Searching for TeV gamma-ray emission from SGR\,1935+2154 during its 2020 X-ray and radio bursting phase

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Anguner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlohr, B. Bi, M. Bottcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose , et al. (230 additional authors not shown)

Abstract: Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preced… ▽ More Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preceded by two gamma-ray outburst alerts by the BAT instrument aboard the Swift satellite, which triggered follow-up observations by the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. has observed SGR\,1935+2154 for 2 hr on 2020 April 28. The observations are coincident with X-ray bursts from the magnetar detected by INTEGRAL and Fermi-GBM, thus providing the first very high energy (VHE) gamma-ray observations of a magnetar in a flaring state. High-quality data acquired during these follow-up observations allow us to perform a search for short-time transients. No significant signal at energies $E>0.6$~TeV is found and upper limits on the persistent and transient emission are derived. We here present the analysis of these observations and discuss the obtained results and prospects of the H.E.S.S. follow-up program for soft gamma-ray repeaters. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Received 2021 June 13; accepted 2021 June 28; published 2021 September 29

Development of an advanced SiPM camera for the Large Size Telescope of the Cherenkov Telescope Array

Authors: M. Heller, T. Armstrong, M. Bellato, A. Bergnoli, M. Bernardos, E. Bernasconi, A. Biland, E. Charbon, D. Corti, M. Dalchenko, D. della Volpe, D. Depaoli, F. Di Pierro, G. Emery, D. Gascón, S. Gómez, R. López-Coto, M. Mariotti, L. D. M. Miranda, T. Montaruli, A. Nagai, R. Rando, T. Saito, H. Tajima, K. Ziȩtara

Abstract: Silicon photomultipliers (SiPMs) have become the baseline choice for cameras of the small-sized telescopes (SSTs) of the Cherenkov Telescope Array (CTA). On the other hand, SiPMs are relatively new to the field and covering large surfaces and operating at high data rates still are challenges to outperform photomultipliers (PMTs). The higher sensitivity in the near infra-red and longer signals comp… ▽ More Silicon photomultipliers (SiPMs) have become the baseline choice for cameras of the small-sized telescopes (SSTs) of the Cherenkov Telescope Array (CTA). On the other hand, SiPMs are relatively new to the field and covering large surfaces and operating at high data rates still are challenges to outperform photomultipliers (PMTs). The higher sensitivity in the near infra-red and longer signals compared to PMTs result in higher night sky background rate for SiPMs. However, the robustness of the SiPMs represents a unique opportunity to ensure long-term operation with low maintenance and better duty cycle than PMTs. The proposed camera for large size telescopes will feature $0.05 degree pixels, low power and fast front-end electronics and a fully digital readout. In this work, we present the status of dedicated simulations and data analysis for the performance estimation. The design features and the different strategies identified, so far, to tackle the demanding requirements and the improved performance are described. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 8 pages, Proceedings for 37th ICRC 2021

Universal description of wetting on multiscale surfaces using integral geometry

Authors: Chenhao Sun, James McClure, Steffen Berg, Peyman Mostaghimi, Ryan T. Armstrong

Abstract: Hypothesis Emerging energy-related technologies deal with multiscale hierarchical structures, intricate surface morphology, non-axisymmetric interfaces, and complex contact lines where wetting is difficult to quantify with classical methods. We hypothesis that a universal description of wetting on multiscale surfaces can be developed by using integral geometry coupled to thermodynamic laws. The… ▽ More Hypothesis Emerging energy-related technologies deal with multiscale hierarchical structures, intricate surface morphology, non-axisymmetric interfaces, and complex contact lines where wetting is difficult to quantify with classical methods. We hypothesis that a universal description of wetting on multiscale surfaces can be developed by using integral geometry coupled to thermodynamic laws. The proposed approach separates the different hierarchy levels of physical description from the thermodynamic description, allowing for a universal description of wetting on multiscale surfaces. Theory and Simulations The theoretical framework is presented followed by application to limiting cases of various wetting states. The wetting states include those considered in the Wenzel, Cassie-Baxter and wicking state models. The wetting behaviour of multiscale surfaces is then explored by conducting direct simulations of a fluid droplet on a structurally rough surface and a chemically heterogeneous surface. Findings The underlying origin of the classical wetting models is shown to be rooted within the proposed theoretical framework. In addition, integral geometry provides a topological-based wetting metric that is not contingent on any type of wetting state. Being of geometrical origin the wetting metric is applicable to describe any type of wetting phenomena on complex surfaces. The proposed framework is applicable to any complex fluid topology and multiscale surface. △ Less

Submitted 21 August, 2021; originally announced August 2021.

LMC N132D: A mature supernova remnant with a power-law gamma-ray spectrum extending beyond 8 TeV

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (212 additional authors not shown)

Abstract: We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for mod… ▽ More We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. We unambiguously detect N132D at very high energies (VHE) with a significance of 5.7 sigma. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position. [Abridged] △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 655, A7 (2021)

Monte Carlo Simulations and Validation of NectarCAM, a Medium Sized Telescope Camera for CTA

Authors: Thomas Armstrong, Heide Costantini, Jean-François Glicenstein, Jean-Philippe Lenain, Ullrich Schwanke, Thomas Tavernier

Abstract: The upcoming Cherenkov Telescope Array (CTA) ground-based gamma-ray observatory will open up our view of the very high energy Universe, offering an improvement in sensitivity of 5-10 times that of previous experiments. NectarCAM is one of the proposed cameras for the Medium-Sized Telescopes (MST) which have been designed to cover the core energy range of CTA, from 100 GeV to 10 TeV. The final came… ▽ More The upcoming Cherenkov Telescope Array (CTA) ground-based gamma-ray observatory will open up our view of the very high energy Universe, offering an improvement in sensitivity of 5-10 times that of previous experiments. NectarCAM is one of the proposed cameras for the Medium-Sized Telescopes (MST) which have been designed to cover the core energy range of CTA, from 100 GeV to 10 TeV. The final camera will be capable of GHz sampling and provide a field of view of 8 degrees with its 265 modules of 7 photomultiplier each (for a total of 1855 pixels). In order to validate the performance of NectarCAM, a partially-equipped prototype has been constructed consisting of only the inner 61-modules. It has so far undergone testing at the integration test-bench facility in CEA Paris-Saclay (France) and on a prototype of the MST structure in Adlershof (Germany). To characterize the performance of the prototype, Monte Carlo simulations were conducted using a detailed model of the 61 module camera in the CORSIKA/sim_telarray framework. This contribution provides an overview of this work including the comparison of trigger and readout performance on test-bench data and trigger and image parameterization performance during on-sky measurements. △ Less

Submitted 1 August, 2021; originally announced August 2021.

Comments: 8 pages, 5 figures

Toward Interlanguage Parallel Scripting for Distributed-Memory Scientific Computing

Authors: Justin M. Wozniak, Timothy G. Armstrong, Ketan C. Maheshwari, Daniel S. Katz, Michael Wilde, Ian T. Foster

Abstract: Scripting languages such as Python and R have been widely adopted as tools for the productive development of scientific software because of the power and expressiveness of the languages and available libraries. However, deploying scripted applications on large-scale parallel computer systems such as the IBM Blue Gene/Q or Cray XE6 is a challenge because of issues including operating system limitat… ▽ More Scripting languages such as Python and R have been widely adopted as tools for the productive development of scientific software because of the power and expressiveness of the languages and available libraries. However, deploying scripted applications on large-scale parallel computer systems such as the IBM Blue Gene/Q or Cray XE6 is a challenge because of issues including operating system limitations, interoperability challenges, parallel filesystem overheads due to the small file system accesses common in scripted approaches, and other issues. We present here a new approach to these problems in which the Swift scripting system is used to integrate high-level scripts written in Python, R, and Tcl, with native code developed in C, C++, and Fortran, by linking Swift to the library interfaces to the script interpreters. In this approach, Swift handles data management, movement, and marshaling among distributed-memory processes without direct user manipulation of low-level communication libraries such as MPI. We present a technique to efficiently launch scripted applications on large-scale supercomputers using a hierarchical programming model. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Comments: 2015 IEEE International Conference on Cluster Computing

TeV emission of Galactic plane sources with HAWC and H.E.S.S

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, P. Brun , et al. (299 additional authors not shown)

Abstract: The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their… ▽ More The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques. △ Less

Submitted 8 September, 2021; v1 submitted 3 July, 2021; originally announced July 2021.

Evidence of 100 TeV $γ$-ray emission from HESS J1702-420: A new PeVatron candidate

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun, P. Brun, M. Bryan , et al. (211 additional authors not shown)

Abstract: The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and p… ▽ More The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and processed using improved analysis techniques. The analysis configuration was optimized to enhance the collection area at the highest energies. We applied a three-dimensional (3D) likelihood analysis to model the source region and adjust non thermal radiative spectral models to the $γ$-ray data. We also analyzed archival data from the Fermi Large Area Telescope (LAT) to constrain the source spectrum at $γ$-ray energies >10 GeV. We report the detection of a new source component called HESS J1702-420A, that was separated from the bulk of TeV emission at a $5.4σ$ confidence level. The power law $γ$-ray spectrum of HESS J1702-420A extends with an index of $Γ=1.53\pm0.19_\text{stat}\pm0.20_\text{sys}$ and without curvature up to the energy band 64-113 TeV, in which it was detected by H.E.S.S. at a $4.0σ$ confidence level. This brings evidence for the source emission up to $100\,\text{TeV}$, which makes HESS J1702-420A a compelling candidate site for the presence of extremely high energy cosmic rays. Remarkably, in a hadronic scenario, the cut-off energy of the proton distribution powering HESS J1702-420A is found to be higher than 0.5 PeV at a 95% confidence level. HESS J1702-420A becomes therefore one of the most solid PeVatron candidates detected so far in H.E.S.S. data, altough a leptonic origin of its emission could not be ruled out either. △ Less

Submitted 14 June, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the 2. Astrophysical processes section of Astronomy & Astrophysics

Search for dark matter annihilation signals from unidentified Fermi-LAT objects with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun , et al. (205 additional authors not shown)

Abstract: Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for… ▽ More Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for very-high-energy (E $\gtrsim$ 100 GeV) gamma-ray emissions using H.E.S.S. observations towards four selected UFOs. Since no significant very-high-energy gamma-ray emission is detected in any dataset of the four observed UFOs nor in the combined UFO dataset, strong constraints are derived on the product of the velocity-weighted annihilation cross section $\langle σv \rangle$ by the $J$-factor for the dark matter models. The 95% C.L. observed upper limits derived from combined H.E.S.S. observations reach $\langle σv \rangle J$ values of 3.7$\times$10$^{-5}$ and 8.1$\times$10$^{-6}$ GeV$^2$cm$^{-2}$s$^{-1}$ in the $W^+W^-$ and $τ^+τ^-$ channels, respectively, for a 1 TeV dark matter mass. Focusing on thermal WIMPs, the H.E.S.S. constraints restrict the $J$-factors to lie in the range 6.1$\times$10$^{19}$ - 2.0$\times$10$^{21}$ GeV$^2$cm$^{-5}$, and the masses to lie between 0.2 and 6 TeV in the $W^+W^-$ channel. For the $τ^+τ^-$ channel, the $J$-factors lie in the range 7.0$\times$10$^{19}$ - 7.1$\times$10$^{20}$ GeV$^2$cm$^{-5}$ and the masses lie between 0.2 and 0.5 TeV. Assuming model-dependent predictions from cosmological N-body simulations on the $J$-factor distribution for Milky Way-sized galaxies, the dark matter models with masses greater than 0.3 TeV for the UFO emissions can be ruled out at high confidence level. △ Less

Submitted 15 June, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 11 pages, 7 figures, matches accepted version in The Astrophysical Journal

Journal ref: Astrophys.J., 918, 17 (2021)

Search for dark matter annihilation in the dwarf irregular galaxy WLM with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (211 additional authors not shown)

Abstract: We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dar… ▽ More We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dark matter dominated objects with well measured kinematics and small uncertainties on their dark matter distribution profiles. In 2018, the H.E.S.S. five-telescope array observed the dwarf irregular galaxy WLM for 18 hours. We present the first analysis based on data obtained from an imaging atmospheric Cherenkov telescope for this subclass of dwarf galaxy. As we do not observe any significant excess in the direction of WLM, we interpret the result in terms of constraints on the velocity-weighted cross section for dark matter pair annihilation as a function of the dark matter particle mass for various continuum channels as well as the prompt gamma-gamma emission. For the $τ^+τ^-$ channel the limits reach a $\langle σv \rangle$ value of about $4\times 10^{-22}$ cm3s-1 for a dark matter particle mass of 1 TeV. For the prompt gamma-gamma channel, the upper limit reaches a $\langle σv \rangle$ value of about $5 \times10^{-24}$ cm3s-1 for a mass of 370 GeV. These limits represent an improvement of up to a factor 200 with respect to previous results for the dwarf irregular galaxies for TeV dark matter search. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. D 103, 102002 (2021)

Deep learning with photosensor timing information as a background rejection method for the Cherenkov Telescope Array

Authors: Samuel Spencer, Thomas Armstrong, Jason Watson, Salvatore Mangano, Yves Renier, Garret Cotter

Abstract: New deep learning techniques present promising new analysis methods for Imaging Atmospheric Cherenkov Telescopes (IACTs) such as the upcoming Cherenkov Telescope Array (CTA). In particular, the use of Convolutional Neural Networks (CNNs) could provide a direct event classification method that uses the entire information contained within the Cherenkov shower image, bypassing the need to Hillas para… ▽ More New deep learning techniques present promising new analysis methods for Imaging Atmospheric Cherenkov Telescopes (IACTs) such as the upcoming Cherenkov Telescope Array (CTA). In particular, the use of Convolutional Neural Networks (CNNs) could provide a direct event classification method that uses the entire information contained within the Cherenkov shower image, bypassing the need to Hillas parameterise the image and allowing fast processing of the data. Existing work in this field has utilised images of the integrated charge from IACT camera photomultipliers, however the majority of current and upcoming generation IACT cameras have the capacity to read out the entire photosensor waveform following a trigger. As the arrival times of Cherenkov photons from Extensive Air Showers (EAS) at the camera plane are dependent upon the altitude of their emission and the impact distance from the telescope, these waveforms contain information potentially useful for IACT event classification. In this test-of-concept simulation study, we investigate the potential for using these camera pixel waveforms with new deep learning techniques as a background rejection method, against both proton and electron induced EAS. We find that a means of utilising their information is to create a set of seven additional 2-dimensional pixel maps of waveform parameters, to be fed into the machine learning algorithm along with the integrated charge image. Whilst we ultimately find that the only classification power against electrons is based upon event direction, methods based upon timing information appear to out-perform similar charge based methods for gamma/hadron separation. We also review existing methods of event classifications using a combination of deep learning and timing information in other astroparticle physics experiments. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: Full paper accepted in Astroparticle Physics. 39 Pages with 11 Figures. Minimal code to reproduce results in the paper available at: https://github.com/STSpencer/wavelearn_release. Some early results previously presented at ICRC2019 (doi:10.22323/1.358.0798)

Report number: 102579

Journal ref: Astroparticle Physics 129C (2021) 102579

Bias-Aware Inference in Regularized Regression Models

Authors: Timothy B. Armstrong, Michal Kolesár, Soonwoo Kwon

Abstract: We consider inference on a scalar regression coefficient under a constraint on the magnitude of the control coefficients. A class of estimators based on a regularized propensity score regression is shown to exactly solve a tradeoff between worst-case bias and variance. We derive confidence intervals (CIs) based on these estimators that are bias-aware: they account for the possible bias of the esti… ▽ More We consider inference on a scalar regression coefficient under a constraint on the magnitude of the control coefficients. A class of estimators based on a regularized propensity score regression is shown to exactly solve a tradeoff between worst-case bias and variance. We derive confidence intervals (CIs) based on these estimators that are bias-aware: they account for the possible bias of the estimator. Under homoskedastic Gaussian errors, these estimators and CIs are near-optimal in finite samples for MSE and CI length. We also provide conditions for asymptotic validity of the CI with unknown and possibly heteroskedastic error distribution, and derive novel optimal rates of convergence under high-dimensional asymptotics that allow the number of regressors to increase more quickly than the number of observations. Extensive simulations and an empirical application illustrate the performance of our methods. △ Less

Submitted 10 August, 2023; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: 49 pages, including all appendices

Observation of a sudden cessation of a very-high-energy gamma-ray flare in PKS 1510-089 with H.E.S.S. and MAGIC in May 2016

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Arm, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, M. de Bony de Lavergne, J. Bregeon , et al. (409 additional authors not shown)

Abstract: The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in… ▽ More The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in this source. While a common variability timescale of $1.5\,$hr is found, there is a significant deviation near the end of the flare with a timescale of $\sim 20\,$min marking the cessation of the event. The peak flux is nearly two orders of magnitude above the low-level emission. For the first time, curvature is detected in the VHE $γ$-ray spectrum of PKS 1510-089, which is fully explained through absorption by the extragalactic background light. Optical R-band observations with ATOM reveal a counterpart of the $γ$-ray flare, even though the detailed flux evolution differs from the VHE ightcurve. Interestingly, a steep flux decrease is observed at the same time as the cessation of the VHE flare. In the high energy (HE, $E>100\,$MeV) $γ$-ray band only a moderate flux increase is observed with Fermi-LAT, while the HE $γ$-ray spectrum significantly hardens up to a photon index of 1.6. A search for broad-line region (BLR) absorption features in the $γ$-ray spectrum indicates that the emission region is located outside of the BLR. Radio VLBI observations reveal a fast moving knot interacting with a standing jet feature around the time of the flare. As the standing feature is located $\sim 50\,$pc from the black hole, the emission region of the flare may have been located at a significant distance from the black hole. If this correlation is indeed true, VHE $γ$ rays have been produced far down the jet where turbulent plasma crosses a standing shock. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 25 pages, 17 figures. Accepted for publication by Astronomy & Astrophysics. Corresponding authors: M. Zacharias, J. Sitarek, D. Sanchez, T. Terzic

Journal ref: A&A 648, A23 (2021)

Capillary fluctuations and energy dynamics for flow in porous media

Authors: James E. McClure, Steffen Berg, Ryan T. Armstrong

Abstract: Capillary energy barriers have important consequences for immiscible fluid flow in porous media. We derive time-and-space averaging theory to account for non-equilibrium behavior and understand the role of athermal capillary fluctuations in the context of their relationship to larger scale phenomenological equations. The formulation resolves several key challenges associated with two-fluid flow in… ▽ More Capillary energy barriers have important consequences for immiscible fluid flow in porous media. We derive time-and-space averaging theory to account for non-equilibrium behavior and understand the role of athermal capillary fluctuations in the context of their relationship to larger scale phenomenological equations. The formulation resolves several key challenges associated with two-fluid flow in porous media: (1) geometric and thermodynamic quantities are constructed as smooth functions of time based on time-and space averages; (2) averaged thermodynamics are developed for films; (3) multi-scale fluctuation terms are identified, which account for transient behaviours of interfaces and films that occur due to pore-scale events; (4) geometric constraints are derived and imposed on the averaged thermodynamics; (5) a new constitutive model is proposed for capillary pressure dynamics that includes contributions from films; and (6) a time-and-space criterion for representative elementary volume (REV) is established based on capillary fluctuations. Capillary fluctuations are assessed quantitatively based on pore-scale simulations and experimental core-flooding data. △ Less

Submitted 19 May, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Thermodynamics of fluctuations based on time-and-space averages

Authors: James E. McClure, Steffen Berg, Ryan T. Armstrong

Abstract: We develop non-equilibrium theory by using averages in time and space as a generalized way to upscale thermodynamics in non-ergodic systems. The approach offers a classical perspective on the energy dynamics in fluctuating systems. The rate of entropy production is shown to be explicitly scale dependent when considered in this context. We show that while any stationary process can be represented a… ▽ More We develop non-equilibrium theory by using averages in time and space as a generalized way to upscale thermodynamics in non-ergodic systems. The approach offers a classical perspective on the energy dynamics in fluctuating systems. The rate of entropy production is shown to be explicitly scale dependent when considered in this context. We show that while any stationary process can be represented as having zero entropy production, second law constraints due to the Clausius theorem are preserved due to the fact that heat and work are related based on conservation of energy. As a demonstration we consider the energy dynamics for the Carnot cycle and for Maxwell's demon. We then consider non-stationary processes, applying time-and-space averages to characterize non-ergodic effects in heterogeneous systems where energy barriers such as compositional gradients are present. We show that the derived theory can be used to understand the origins of anomalous diffusion phenomena in systems where Fick's law applies at small length scales but not at large length scales. We then characterize fluctuations in capillary-dominated systems, which are non-stationary due to the irreversibility of cooperative events. △ Less

Submitted 21 May, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. E 104, 035106 (2021)

Transporter Networks: Rearranging the Visual World for Robotic Manipulation

Authors: Andy Zeng, Pete Florence, Jonathan Tompson, Stefan Welker, Jonathan Chien, Maria Attarian, Travis Armstrong, Ivan Krasin, Dan Duong, Ayzaan Wahid, Vikas Sindhwani, Johnny Lee

Abstract: Robotic manipulation can be formulated as inducing a sequence of spatial displacements: where the space being moved can encompass an object, part of an object, or end effector. In this work, we propose the Transporter Network, a simple model architecture that rearranges deep features to infer spatial displacements from visual input - which can parameterize robot actions. It makes no assumptions of… ▽ More Robotic manipulation can be formulated as inducing a sequence of spatial displacements: where the space being moved can encompass an object, part of an object, or end effector. In this work, we propose the Transporter Network, a simple model architecture that rearranges deep features to infer spatial displacements from visual input - which can parameterize robot actions. It makes no assumptions of objectness (e.g. canonical poses, models, or keypoints), it exploits spatial symmetries, and is orders of magnitude more sample efficient than our benchmarked alternatives in learning vision-based manipulation tasks: from stacking a pyramid of blocks, to assembling kits with unseen objects; from manipulating deformable ropes, to pushing piles of small objects with closed-loop feedback. Our method can represent complex multi-modal policy distributions and generalizes to multi-step sequential tasks, as well as 6DoF pick-and-place. Experiments on 10 simulated tasks show that it learns faster and generalizes better than a variety of end-to-end baselines, including policies that use ground-truth object poses. We validate our methods with hardware in the real world. Experiment videos and code are available at https://transporternets.github.io △ Less

Submitted 5 January, 2022; v1 submitted 27 October, 2020; originally announced October 2020.

Comments: Project webpage: https://transporternets.github.io Summary video: https://youtu.be/8afHfReCfPo?t=12214

An extreme particle accelerator in the Galactic plane: HESS J1826$-$130

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, P. Bordas, M. Breuhaus , et al. (215 additional authors not shown)

Abstract: The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $γ$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $γ$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The… ▽ More The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $γ$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $γ$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The source spectrum is best fit with either a power-law function with a spectral index $Γ$ = 1.78 $\pm$ 0.10$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ and an exponential cut-off at 15.2$^{+5.5}_{-3.2}$ TeV, or a broken power-law with $Γ_{1}$ = 1.96 $\pm$ 0.06$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$, $Γ_{2}$ = 3.59 $\pm$ 0.69$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ for energies below and above $E_{\rm{br}}$ = 11.2 $\pm$ 2.7 TeV, respectively. The VHE flux from HESS J1826$-$130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula (PWN), HESS J1825$-$137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826$-$130 VHE emission related to PSR J1826$-$1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826$-$130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to $\gtrsim$200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants (SNRs), molecular clouds, and counterparts detected in radio, X-rays, and TeV energies. △ Less

Submitted 25 October, 2020; originally announced October 2020.

Comments: 9 Pages, 5 Figures

Journal ref: A&A 644, A112 (2020)

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

The human behavioural immune system is a product of cultural evolution

Authors: Edwin S. Dalmaijer, Thomas Armstrong

Abstract: Disgust is a basic emotion that serves to avoid contaminants, and is central to the behavioural immune system. While disgust-motivated avoidance occurs in bonobos and chimpanzees, humans show uniquely high levels of contamination sensitivity. Current theory postulates that human disgust is primarily a genetic adaptation, side-lining social transmission and denying parental modelling. Here, we test… ▽ More Disgust is a basic emotion that serves to avoid contaminants, and is central to the behavioural immune system. While disgust-motivated avoidance occurs in bonobos and chimpanzees, humans show uniquely high levels of contamination sensitivity. Current theory postulates that human disgust is primarily a genetic adaptation, side-lining social transmission and denying parental modelling. Here, we test whether this strong view is warranted by simulating 100000 years of cultural and biological evolution in various scenarios. We modelled disgust as a trait that governed the extent to which individuals forewent potentially contaminated nutrition. This lowered their risk of gastrointestinal illness, at the cost of increasing starvation risk and birth-interval. As environmental levels and costs of contamination rose, so did relative fitness of high-disgust individuals. Evolutionary shifts in disgust as polygenic trait occurred as a consequence of germ-cell mutations, but were most prominent in populations with high initial genetic variance. Crucially, cultural transmission between generations operated at a higher rate, even if parental modelling was eliminated. This study serves not only as evidence of cultural evolution shaping the behavioural immune system, but also as an illustration of emerging theories that paint affective and cognitive mechanisms as socially transmitted rather than biologically determined functions. △ Less

Submitted 14 February, 2022; v1 submitted 30 August, 2020; originally announced August 2020.

Comments: 18 pages, 2 figures

Search for dark matter signals towards a selection of recently-detected DES dwarf galaxy satellites of the Milky Way with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdallah, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, M. Breuhaus , et al. (206 additional authors not shown)

Abstract: Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-base… ▽ More Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-based observations with the imaging atmospheric Cherenkov telescope array H.E.S.S. We present a search for very-high-energy ($E\gtrsim100$ GeV) gamma-ray emission using H.E.S.S. observations carried out recently towards Reticulum II, Tucana II, Tucana III, Tucana IV and Grus II satellites. No significant very-high-energy gamma-ray excess is found from the observations on any individual object nor in the combined analysis of all the datasets. Using the most recent modeling of the dark matter distribution in the dwarf galaxy halo, we compute for the first time on DES satellites individual and combined constraints from Cherenkov telescope observations on the annihilation cross section of dark matter particles in the form of Weakly Interacting Massive Particles. The combined 95% C.L. observed upper limits reach $\langle σv \rangle \simeq 1 \times 10^{-23}$ cm$^3$s$^{-1}$ in the $W^+W^-$ channel and $4 \times 10^{-26}$ cm$^3$s$^{-1}$ in the $γγ$ channels for a dark matter mass of 1.5 TeV. The H.E.S.S. constraints well complement the results from Fermi-LAT, HAWC, MAGIC and VERITAS and are currently the most stringent in the $γγ$ channels in the multi-GeV/multi-TeV mass range. △ Less

Submitted 4 September, 2020; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: 20 pages, 12 figures

Journal ref: Phys. Rev. D 102, 062001 (2020)

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

Characterization of wetting using topological principles

Authors: Chenhao Sun, James E. McClure, Peyman Mostaghimi, Anna L. Herring, Douglas E. Meisenheimer, Dorthe Wildenschild, Steffen Berg, Ryan T. Armstrong

Abstract: Hypothesis Understanding wetting behavior is of great importance for natural systems and technological applications. The traditional concept of contact angle, a purely geometrical measure related to curvature, is often used for characterizing the wetting state of a system. It can be determined from Young's equation by applying equilibrium thermodynamics. However, whether contact angle is a represe… ▽ More Hypothesis Understanding wetting behavior is of great importance for natural systems and technological applications. The traditional concept of contact angle, a purely geometrical measure related to curvature, is often used for characterizing the wetting state of a system. It can be determined from Young's equation by applying equilibrium thermodynamics. However, whether contact angle is a representative measure of wetting for systems with significant complexity is unclear. Herein, we hypothesize that topological principles based on the Gauss-Bonnet theorem could yield a robust measure to characterize wetting. Theory and Experiments We introduce a macroscopic contact angle based on the deficit curvature of the fluid interfaces that are imposed by contacts with other immiscible phases. We perform sessile droplet simulations followed by multiphase experiments for porous sintered glass and Bentheimer sandstone to assess the sensitivity and robustness of the topological approach and compare the results to other traditional approaches. Findings We show that the presented topological principle is consistent with thermodynamics under the simplest conditions through a variational analysis. Furthermore, we elucidate that at sufficiently high image resolution the proposed topological approach and local contact angle measurements are comparable. While at lower resolutions, the proposed approach provides more accurate results being robust to resolution-based effects. Overall, the presented concepts open new pathways to characterize the wetting state of complex systems and theoretical developments to study multiphase systems. △ Less

Submitted 6 May, 2020; originally announced May 2020.

Comments: 11 pages, 9 figures, 1 table

ML-LBM: Machine Learning Aided Flow Simulation in Porous Media

Authors: Ying Da Wang, Traiwit Chung, Ryan T. Armstrong, Peyman Mostaghimi

Abstract: Simulation of fluid flow in porous media has many applications, from the micro-scale (cell membranes, filters, rocks) to macro-scale (groundwater, hydrocarbon reservoirs, and geothermal) and beyond. Direct simulation of flow in porous media requires significant computational resources to solve within reasonable timeframes. An integrated method combining predictions of fluid flow (fast, limited acc… ▽ More Simulation of fluid flow in porous media has many applications, from the micro-scale (cell membranes, filters, rocks) to macro-scale (groundwater, hydrocarbon reservoirs, and geothermal) and beyond. Direct simulation of flow in porous media requires significant computational resources to solve within reasonable timeframes. An integrated method combining predictions of fluid flow (fast, limited accuracy) with direct flow simulation (slow, high accuracy) is outlined. In the tortuous flow paths of porous media, Deep Learning techniques based on Convolutional Neural Networks (CNNs) are shown to give an accurate estimate of the steady state velocity fields (in all axes), and by extension, the macro-scale permeability. This estimate can be used as-is, or as initial conditions in direct simulation to reach a fully accurate result in a fraction of the compute time. A Gated U-Net Convolutional Neural Network is trained on a datasets of 2D and 3D porous media generated by correlated fields, with their steady state velocity fields calculated from direct LBM simulation. Sensitivity analysis indicates that network accuracy is dependent on (1) the tortuosity of the domain, (2) the size of convolution filters, (3) the use of distance maps as input, (4) the use of mass conservation loss functions. Permeability estimation from these predicted fields reaches over 90\% accuracy for 80\% of cases. It is further shown that these velocity fields are error prone when used for solute transport simulation. Using the predicted velocity fields as initial conditions is shown to accelerate direct flow simulation to physically true steady state conditions an order of magnitude less compute time. Using Deep Learning predictions (or potentially any other approximation method) to accelerate flow simulation to steady state in complex pore structures shows promise as a technique push the boundaries fluid flow modelling. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 23 pages, 16 figures

Probing the magnetic field in the GW170817 outflow using H.E.S.S. observations

Authors: H. E. S. S. Collaboration, :, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa-Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (209 additional authors not shown)

Abstract: The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $γ$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which laste… ▽ More The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $γ$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios. △ Less

Submitted 18 May, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 7 pages, 2 figures, 1 table

Journal ref: The Astrophysical Journal Letters, 894, L16, 2020

Robust Empirical Bayes Confidence Intervals

Authors: Timothy B. Armstrong, Michal Kolesár, Mikkel Plagborg-Møller

Abstract: We construct robust empirical Bayes confidence intervals (EBCIs) in a normal means problem. The intervals are centered at the usual linear empirical Bayes estimator, but use a critical value accounting for shrinkage. Parametric EBCIs that assume a normal distribution for the means (Morris, 1983b) may substantially undercover when this assumption is violated. In contrast, our EBCIs control coverage… ▽ More We construct robust empirical Bayes confidence intervals (EBCIs) in a normal means problem. The intervals are centered at the usual linear empirical Bayes estimator, but use a critical value accounting for shrinkage. Parametric EBCIs that assume a normal distribution for the means (Morris, 1983b) may substantially undercover when this assumption is violated. In contrast, our EBCIs control coverage regardless of the means distribution, while remaining close in length to the parametric EBCIs when the means are indeed Gaussian. If the means are treated as fixed, our EBCIs have an average coverage guarantee: the coverage probability is at least $1 - α$ on average across the $n$ EBCIs for each of the means. Our empirical application considers the effects of U.S. neighborhoods on intergenerational mobility. △ Less

Submitted 14 May, 2022; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: 45 pages plus a 25-page supplemental appendix

Journal ref: Econometrica, Volume 90, Issue 6, November 2021, pages 2567-2602

Very high energy $γ$-ray emission from two blazars of unknown redshift and upper limits on their distance

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus , et al. (204 additional authors not shown)

Abstract: We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $γ$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum t… ▽ More We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $γ$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $γ$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$. △ Less

Submitted 20 April, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS. 13 pages, 4 figures

Big Data Staging with MPI-IO for Interactive X-ray Science

Authors: Justin M. Wozniak, Hemant Sharma, Timothy G. Armstrong, Michael Wilde, Jonathan D. Almer, Ian Foster

Abstract: New techniques in X-ray scattering science experiments produce large data sets that can require millions of high-performance processing hours per week of computation for analysis. In such applications, data is typically moved from X-ray detectors to a large parallel file system shared by all nodes of a petascale supercomputer and then is read repeatedly as different science application tasks proce… ▽ More New techniques in X-ray scattering science experiments produce large data sets that can require millions of high-performance processing hours per week of computation for analysis. In such applications, data is typically moved from X-ray detectors to a large parallel file system shared by all nodes of a petascale supercomputer and then is read repeatedly as different science application tasks proceed. However, this straightforward implementation causes significant contention in the file system. We propose an alternative approach in which data is instead staged into and cached in compute node memory for extended periods, during which time various processing tasks may efficiently access it. We describe here such a big data staging framework, based on MPI-IO and the Swift parallel scripting language. We discuss a range of large-scale data management issues involved in X-ray scattering science and measure the performance benefits of the new staging framework for high-energy diffraction microscopy, an important emerging application in data-intensive X-ray scattering. We show that our framework accelerates scientific processing turnaround from three months to under 10 minutes, and that our I/O technique reduces input overheads by a factor of 5 on 8K Blue Gene/Q nodes. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 10 pages, submitted to International Symposium on Big Data Computing (BDC 2014)

Physical Accuracy of Deep Neural Networks for 2D and 3D Multi-Mineral Segmentation of Rock micro-CT Images

Authors: Ying Da Wang, Mehdi Shabaninejad, Ryan T. Armstrong, Peyman Mostaghimi

Abstract: Segmentation of 3D micro-Computed Tomographic uCT) images of rock samples is essential for further Digital Rock Physics (DRP) analysis, however, conventional methods such as thresholding, watershed segmentation, and converging active contours are susceptible to user-bias. Deep Convolutional Neural Networks (CNNs) have produced accurate pixelwise semantic segmentation results with natural images an… ▽ More Segmentation of 3D micro-Computed Tomographic uCT) images of rock samples is essential for further Digital Rock Physics (DRP) analysis, however, conventional methods such as thresholding, watershed segmentation, and converging active contours are susceptible to user-bias. Deep Convolutional Neural Networks (CNNs) have produced accurate pixelwise semantic segmentation results with natural images and $μ$CT rock images, however, physical accuracy is not well documented. The performance of 4 CNN architectures is tested for 2D and 3D cases in 10 configurations. Manually segmented uCT images of Mt. Simon Sandstone are treated as ground truth and used as training and validation data, with a high voxelwise accuracy (over 99%) achieved. Downstream analysis is then used to validate physical accuracy. The topology of each segmented phase is calculated, and the absolute permeability and multiphase flow is modelled with direct simulation in single and mixed wetting cases. These physical measures of connectivity, and flow characteristics show high variance and uncertainty, with models that achieve 95\%+ in voxelwise accuracy possessing permeabilities and connectivities orders of magnitude off. A new network architecture is also introduced as a hybrid fusion of U-net and ResNet, combining short and long skip connections in a Network-in-Network configuration. The 3D implementation outperforms all other tested models in voxelwise and physical accuracy measures. The network architecture and the volume fraction in the dataset (and associated weighting), are factors that not only influence the accuracy trade-off in the voxelwise case, but is especially important in training a physically accurate model for segmentation. △ Less

Submitted 15 February, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: 16 figures

Detection of very-high-energy γ-ray emission from the colliding wind binary η Car with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun , et al. (210 additional authors not shown)

Abstract: Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) γ-ray emission from η Car around the last per… ▽ More Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) γ-ray emission from η Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around η Car was observed with H.E.S.S. between orbital phase p = 0.78 - 1.10, with a closer sampling at p {\approx} 0.95 and p {\approx} 1.10 (assuming a period of 2023 days). Optimised hardware settings as well as adjustments to the data reduction, reconstruction, and signal selection were needed to suppress and take into account the strong, extended, and inhomogeneous night sky background (NSB) in the η Car field of view. Tailored run-wise Monte-Carlo simulations (RWS) were required to accurately treat the additional noise from NSB photons in the instrument response functions. Results. H.E.S.S. detected VHE γ-ray emission from the direction of η Car shortly before and after the minimum in the X-ray light-curve close to periastron. Using the point spread function provided by RWS, the reconstructed signal is point-like and the spectrum is best described by a power law. The overall flux and spectral index in VHE γ rays agree within statistical and systematic errors before and after periastron. The γ-ray spectrum extends up to at least ~400 GeV. This implies a maximum magnetic field in a leptonic scenario in the emission region of 0.5 Gauss. No indication for phase-locked flux variations is detected in the H.E.S.S. data. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 9 pages, 4 figures, 3 tables, in press with A&A

Journal ref: A&A 635, A167 (2020)

Linking continuum-scale state of wetting to pore-scale contact angles in porous media

Authors: Chenhao Sun, James E. McClure, Peyman Mostaghimi, Anna L. Herring, Mehdi Shabaninejad, Steffen Berg, Ryan T. Armstrong

Abstract: Wetting phenomena play a key role in flows through porous media. Relative permeability and capillary pressure-saturation functions show a high sensitivity to wettability, which has different definitions at the continuum- and pore-scale. At the continuum-scale, the state of wetting is defined as Amott-Harvey or USBM (United States Bureau of Mines) by capillary pressure drainage and imbibition cycle… ▽ More Wetting phenomena play a key role in flows through porous media. Relative permeability and capillary pressure-saturation functions show a high sensitivity to wettability, which has different definitions at the continuum- and pore-scale. At the continuum-scale, the state of wetting is defined as Amott-Harvey or USBM (United States Bureau of Mines) by capillary pressure drainage and imbibition cycles. At the pore-scale, the concept of contact angle is used, which until recently was not experimentally possible to determine within an opaque porous medium. Recent progress on measurements of pore-scale contact angles by X-ray computed micro-tomography has therefore attracted significant attention in various research communities. In this work, the Gauss-Bonnet theorem is applied to provide a direct link between capillary pressure saturation $P_c(S_w)$ data and measured distributions of pore-scale contact angles. We propose that the wetting state of a porous medium can be described in terms of geometrical arguments that constrain the morphological state of immiscible fluids. The constraint describes the range of possible contact angles and interfacial curvatures that can exist for a given system. We present measurements in a tested sandstone for which the USBM index, $P_c(S_w)$, and pore-scale contact angles are measured. Additional studies are also performed using two-phase Lattice Boltzmann simulations to test a wider range of wetting conditions. We show that mean pore-scale contact angle measurements can be predicted from petrophysical data within a few differences. This provides a general framework on how continuum-scale data can be used to describe the geometrical state of fluids within porous media. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 8 pages, 6 figures

Interferometric Fringe Visibility Null as a Function of Spatial Frequency: a Probe of Stellar Atmospheres

Authors: J. T. Armstrong, A. M. Jorgensen, D. Mozurkewich, H. R. Neilson, E. K. Baines, H. R. Schmitt, G. T. van Belle

Abstract: We introduce an observational tool based on visibility nulls in optical spectro-interferometry fringe data to probe the structure of stellar atmospheres. In a preliminary demonstration, we use both Navy Precision Optical Interferometer (NPOI) data and stellar atmosphere models to show that this tool can be used, for example, to investigate limb darkening. Using bootstrapping with either multiple… ▽ More We introduce an observational tool based on visibility nulls in optical spectro-interferometry fringe data to probe the structure of stellar atmospheres. In a preliminary demonstration, we use both Navy Precision Optical Interferometer (NPOI) data and stellar atmosphere models to show that this tool can be used, for example, to investigate limb darkening. Using bootstrapping with either multiple linked baselines or multiple wavelengths in optical and infrared spectro-interferometric observations of stars makes it possible to measure the spatial frequency $u_0$ at which the real part of the fringe visibility ${\rm Re}(V)$ vanishes. That spatial frequency is determined by $u_0 = B_\perp/λ_0$, where $B_\perp$ is the projected baseline length, and $λ_0$ is the wavelength at which the null is observed. Since $B_\perp$ changes with the Earth's rotation, $λ_0$ also changes. If $u_0$ is constant with wavelength, $λ_0$ varies in direct proportion to $B_\perp$. Any departure from that proportionality indicates that the brightness distribution across the stellar disk varies with wavelength via variations in limb darkening, in the angular size of the disk, or both. In this paper, we introduce the use of variations of $u_0$ with $λ$ as a means of probing the structure of stellar atmospheres. Using the equivalent uniform disk diameter $θ_{\rm UD, 0}(λ_0)$, given by $θ_{\rm UD, 0} = 1.22/u_0(λ_0)$, as a convenient and intuitive parameterization of $u_0(λ_0)$, we demonstrate this concept by using model atmospheres to calculate the brightness distribution for $ν$ Ophiuchi and predict $θ_{\rm UD, 0}(λ_0)$, and then comparing the predictions to coherently averaged data from observations taken with the NPOI. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Boosting Resolution and Recovering Texture of micro-CT Images with Deep Learning

Authors: Ying Da Wang, Ryan T. Armstrong, Peyman Mostaghimi

Abstract: Digital Rock Imaging is constrained by detector hardware, and a trade-off between the image field of view (FOV) and the image resolution must be made. This can be compensated for with super resolution (SR) techniques that take a wide FOV, low resolution (LR) image, and super resolve a high resolution (HR), high FOV image. The Enhanced Deep Super Resolution Generative Adversarial Network (EDSRGAN)… ▽ More Digital Rock Imaging is constrained by detector hardware, and a trade-off between the image field of view (FOV) and the image resolution must be made. This can be compensated for with super resolution (SR) techniques that take a wide FOV, low resolution (LR) image, and super resolve a high resolution (HR), high FOV image. The Enhanced Deep Super Resolution Generative Adversarial Network (EDSRGAN) is trained on the Deep Learning Digital Rock Super Resolution Dataset, a diverse compilation 12000 of raw and processed uCT images. The network shows comparable performance of 50% to 70% reduction in relative error over bicubic interpolation. GAN performance in recovering texture shows superior visual similarity compared to SRCNN and other methods. Difference maps indicate that the SRCNN section of the SRGAN network recovers large scale edge (grain boundaries) features while the GAN network regenerates perceptually indistinguishable high frequency texture. Network performance is generalised with augmentation, showing high adaptability to noise and blur. HR images are fed into the network, generating HR-SR images to extrapolate network performance to sub-resolution features present in the HR images themselves. Results show that under-resolution features such as dissolved minerals and thin fractures are regenerated despite the network operating outside of trained specifications. Comparison with Scanning Electron Microscope images shows details are consistent with the underlying geometry of the sample. Recovery of textures benefits the characterisation of digital rocks with a high proportion of under-resolution micro-porous features, such as carbonate and coal samples. Images that are normally constrained by the mineralogy of the rock (coal), by fast transient imaging (waterflooding), or by the energy of the source (microporosity), can be super resolved accurately for further analysis downstream. △ Less

Submitted 26 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: \keywords{Digital Rock Imaging \and Super Resolution \and Convolutional Neural Networks \and Generative Adversarial Networks}

Prospects for the Use of Photosensor Timing Information with Machine Learning Techniques in Background Rejection

Authors: Samuel Spencer, Thomas Armstrong, Jason Watson, Garret Cotter

Abstract: Recent developments in machine learning (ML) techniques present a promising new analysis method for high-speed imaging in astroparticle physics experiments, for example with imaging atmospheric Cherenkov telescopes (IACTs). In particular, the use of timing information with new machine learning techniques provides a novel method for event classification. Previous work in this field has utilised ima… ▽ More Recent developments in machine learning (ML) techniques present a promising new analysis method for high-speed imaging in astroparticle physics experiments, for example with imaging atmospheric Cherenkov telescopes (IACTs). In particular, the use of timing information with new machine learning techniques provides a novel method for event classification. Previous work in this field has utilised images of the integrated charge from IACT camera photomultipliers, but the majority of current and upcoming IACT cameras have the capacity to read out the entire photosensor waveform following a trigger. As the arrival times of Cherenkov photons from extensive air showers (EAS) at the camera plane are dependent upon the altitude of their emission, these waveforms contain information useful for IACT event classification. In this work, we investigate the potential for using these waveforms with ML techniques, and find that a highly effective means of utilising their information is to create a set of seven additional two dimensional histograms of waveform parameters to be fed into the machine learning algorithm along with the integrated charge image. This appears to be superior to using only these new ML techniques with the waveform integrated charge alone. We also examine these timing-based ML techniques in the context of other experiments. △ Less

Submitted 10 July, 2019; originally announced July 2019.

Comments: 8 pages, 3 figures, Submitted to Proceedings of the 36th International Cosmic Ray Conference (ICRC 2019), Madison, WI, U.S.A

Journal ref: PoS(ICRC2019)798