title: Testing for Typicality with Respect to an Ensemble of Learned Distributions

abstract: Methods of performing anomaly detection on high-dimensional data sets are needed, since algorithms which are trained on data are only expected to perform well on data that is similar to the training data. There are theoretical results on the ability to detect if a population of data is likely to come from a known base distribution, which is known as the goodness-of-fit problem. One-sample approaches to this problem offer significant computational advantages for online testing, but require knowing a model of the base distribution. The ability to correctly reject anomalous data in this setting hinges on the accuracy of the model of the base distribution. For high dimensional data, learning an accurate-enough model of the base distribution such that anomaly detection works reliably is very challenging, as many researchers have noted in recent years. Existing methods for the one-sample goodness-of-fit problem do not account for the fact that a model of the base distribution is learned. To address that gap, we offer a theoretically motivated approach to account for the density learning procedure. In particular, we propose training an ensemble of density models, considering data to be anomalous if the data is anomalous with respect to any member of the ensemble. We provide a theoretical justification for this approach, proving first that a test on typicality is a valid approach to the goodness-of-fit problem, and then proving that for a correctly constructed ensemble of models, the intersection of typical sets of the models lies in the interior of the typical set of the base distribution. We present our method in the context of an example on synthetic data in which the effects we consider can easily be seen.


title: Robust multi-stage model-based design of optimal experiments for nonlinear estimation

abstract: We study approaches to robust model-based design of experiments in the context of maximum-likelihood estimation. These approaches provide robustification of model-based methodologies for the design of optimal experiments by accounting for the effect of the parametric uncertainty. We study the problem of robust optimal design of experiments in the framework of nonlinear least-squares parameter estimation using linearized confidence regions. We investigate several well-known robustification frameworks in this respect and propose a novel methodology based on multi-stage robust optimization. The proposed methodology aims at problems, where the experiments are designed sequentially with a possibility of re-estimation in-between the experiments. The multi-stage formalism aids in identifying experiments that are better conducted in the early phase of experimentation, where parameter knowledge is poor. We demonstrate the findings and effectiveness of the proposed methodology using four case studies of varying complexity.

categories: stat.ME cs.LG cs.SY eess.SY math.OC stat.ML


title: Clustering of Big Data with Mixed Features

abstract: Clustering large, mixed data is a central problem in data mining. Many approaches adopt the idea of k-means, and hence are sensitive to initialisation, detect only spherical clusters, and require a priori the unknown number of clusters. We here develop a new clustering algorithm for large data of mixed type, aiming at improving the applicability and efficiency of the peak-finding technique. The improvements are threefold: (1) the new algorithm is applicable to mixed data; (2) the algorithm is capable of detecting outliers and clusters of relatively lower density values; (3) the algorithm is competent at deciding the correct number of clusters. The computational complexity of the algorithm is greatly reduced by applying a fast k-nearest neighbors method and by scaling down to component sets. We present experimental results to verify that our algorithm works well in practice. Keywords: Clustering; Big Data; Mixed Attribute; Density Peaks; Nearest-Neighbor Graph; Conductance.


title: A Gateway to Astronomical Image Processing: Vera C. RubinObservatory LSST Science Pipelines on AWS

abstract: The Legacy Survey of Space and Time, operated by the Vera C. Rubin Observatory, is a 10-year astronomical survey due to start operations in 2022 that will image half the sky every three nights. LSST will produce ~20TB of raw data per night which will be calibrated and analyzed in almost real time. Given the volume of LSST data, the traditional subset-download-process paradigm of data reprocessing faces significant challenges. We describe here, the first steps towards a gateway for astronomical science that would enable astronomers to analyze images and catalogs at scale. In this first step we focus on executing the Rubin LSST Science Pipelines, a collection of image and catalog processing algorithms, on Amazon Web Services (AWS). We describe our initial impressions on the performance, scalability and cost of deploying such a system in the cloud.


title: Bayesian inference for transportation origin-destination matrices: the Poisson-inverse Gaussian and other Poisson mixtures

abstract: In this paper we present Poisson mixture approaches for origin-destination (OD) modeling in transportation analysis. We introduce covariate-based models which incorporate different transport modeling phases and also allow for direct probabilistic inference on link traffic based on Bayesian predictions. Emphasis is placed on the Poisson-inverse Gaussian as an alternative to the commonly-used Poisson-gamma and Poisson-lognormal models. We present a first full Bayesian formulation and demonstrate that the Poisson-inverse Gaussian is particularly suited for OD analysis due to desirable marginal and hierarchical properties. In addition, the integrated nested Laplace approximation (INLA) is considered as an alternative to Markov chain Monte Carlo and the two methodologies are compared under specific modeling assumptions. The case study is based on 2001 Belgian census data and focuses on a large, sparsely-distributed OD matrix containing trip information for 308 Flemish municipalities.


title: Saturating stable matchings

abstract: I relate bipartite graph matchings to stable matchings. I prove a necessary and sufficient condition for the existence of a saturating stable matching, where every agent on one side is matched, for all possible preferences. I extend my analysis to perfect stable matchings, where every agent on both sides is matched.


title: Multi-Hypothesis Interactions in Game-Theoretic Motion Planning

abstract: We present a novel method for handling uncertainty about the intentions of non-ego players in dynamic games, with application to motion planning for autonomous vehicles. Equilibria in these games explicitly account for interaction among other agents in the environment, such as drivers and pedestrians. Our method models the uncertainty about the intention of other agents by constructing multiple hypotheses about the objectives and constraints of other agents in the scene. For each candidate hypothesis, we associate a Bernoulli random variable representing the probability of that hypothesis, which may or may not be independent of the probability of other hypotheses. We leverage constraint asymmetries and feedback information patterns to incorporate the probabilities of hypotheses in a natural way. Specifically, increasing the probability associated with a given hypothesis from $0$ to $1$ shifts the responsibility of collision avoidance from the hypothesized agent to the ego agent. This method allows the generation of interactive trajectories for the ego agent, where the level of assertiveness or caution that the ego exhibits is directly related to the easy-to-model uncertainty it maintains about the scene.


title: Comparing Piezoresistive Substrates for Tactile Sensing in Dexterous Hands

abstract: While tactile skins have been shown to be useful for detecting collisions between a robotic arm and its environment, they have not been extensively used for improving robotic grasping and in-hand manipulation. We propose a novel sensor design for use in covering existing multi-fingered robot hands. We analyze the performance of four different piezoresistive materials using both fabric and anti-static foam substrates in benchtop experiments. We find that although the piezoresistive foam was designed as packing material and not for use as a sensing substrate, it performs comparably with fabrics specifically designed for this purpose. While these results demonstrate the potential of piezoresistive foams for tactile sensing applications, they do not fully characterize the efficacy of these sensors for use in robot manipulation. As such, we use a low density foam substrate to develop a scalable tactile skin that can be attached to the palm of a robotic hand. We demonstrate several robotic manipulation tasks using this sensor to show its ability to reliably detect and localize contact, as well as analyze contact patterns during grasping and transport tasks. Our project website provides details on all materials, software, and data used in the sensor development and analysis: https://sites.google.com/gcloud.utah.edu/piezoresistive-tactile-sensing/.


title: Colorado in Context: Congressional Redistricting and Competing Fairness Criteria in Colorado

abstract: In this paper, we apply techniques of ensemble analysis to understand the political baseline for Congressional representation in Colorado. We generate a large random sample of reasonable redistricting plans and determine the partisan balance of each district using returns from state-wide elections in 2018, and analyze the 2011/2012 enacted districts in this context. Colorado recently adopted a new framework for redistricting, creating an independent commission to draw district boundaries, prohibiting partisan bias and incumbency considerations, requiring that political boundaries (such as counties) be preserved as much as possible, and also requiring that mapmakers maximize the number of competitive districts. We investigate the relationships between partisan outcomes, number of counties which are split, and number of competitive districts in a plan. This paper also features two novel improvements in methodology--a more rigorous statistical framework for understanding necessary sample size, and a weighted-graph method for generating random plans which split approximately as few counties as acceptable human-drawn maps.


title: Global Position Control on Underactuated Bipedal Robots: Step-to-step Dynamics Approximation for Step Planning

abstract: Global position control for underactuated bipedal walking is a challenging problem due to the lack of actuation on the feet of the robots. In this paper, we apply the Hybrid-Linear Inverted Pendulum (H-LIP) based stepping on 3D underactuated bipedal robots for global position control. The step-to-step (S2S) dynamics of the H-LIP walking approximates the actual S2S dynamics of the walking of the robot, where the step size is considered as the input. Thus the feedback controller based on the H-LIP approximately controls the robot to behave like the H-LIP, the differences between which stay in an error invariant set. Model Predictive Control (MPC) is applied to the H-LIP for global position control in 3D. The H-LIP stepping then generates desired step sizes for the robot to track. Moreover, turning behavior is integrated with the step planning. The proposed framework is verified on the 3D underactuated bipedal robot Cassie in simulation together with a proof-of-concept experiment.


title: Extremely Low Molecular Gas Content in the Vicinity of a Red Nugget Galaxy at $z=1.91$

abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 5 observations of a galaxy at $z=1.91$, GDS24569, in search of molecular gas in its vicinity via the [C I] $^3$P$_1$-$^3$P$_0$ line. GDS24569 is a massive ($\log M_*/M_\odot=11$) passively evolving galaxy, and characterized by compact morphology with an effective radius of $\sim0.5$ kpc. We apply two blind detection algorithms to the spectral data cubes, and find no promising detection in or around GDS24569 out to projected distance of $\sim320$ kpc, while a narrow tentative line ($4.1 \sigma$) is identified at $+1200$ km/s by one of the algorithms. From the non-detection of [C I], we place a $3\sigma$ upper limit on molecular hydrogen mass, $\sim 7.1 \times 10^9 M_\odot$, which converts to an extremely low gas-to-stellar mass fraction, $< 5 \%$. We conduct a spectral energy distribution modeling by including optical-to-far-infrared data, and find a considerably high ($\sim0.1\%$) dust-to-stellar mass ratio, $\sim10$-$100\times$ higher than those of local early-type galaxies. In combination with a previous result of an insufficient number of surrounding satellite galaxies, it is suggested that GDS24569 is unlikely to experience significant size evolution via satellite mergers. We discuss possible physical mechanisms that quenched GDS24569.


title: Optimization under rare chance constraints

abstract: Chance constraints provide a principled framework to mitigate the risk of high-impact extreme events by modifying the controllable properties of a system. The low probability and rare occurrence of such events, however, impose severe sampling and computational requirements on classical solution methods that render them impractical. This work proposes a novel sampling-free method for solving rare chance constrained optimization problems affected by uncertainties that follow general Gaussian mixture distributions. By integrating modern developments in large deviation theory with tools from convex analysis and bilevel optimization, we propose tractable formulations that can be solved by off-the-shelf solvers. Our formulations enjoy several advantages compared to classical methods: their size and complexity is independent of event rarity, they do not require linearity or convexity assumptions on system constraints, and under easily verifiable conditions, serve as safe conservative approximations or asymptotically exact reformulations of the true problem. Computational experiments on linear, nonlinear and PDE-constrained problems from applications in portfolio management, structural engineering and fluid dynamics illustrate the broad applicability of our method and its advantages over classical sampling-based approaches in terms of both accuracy and efficiency.


title: Efficient magnetic superstructure optimization with $\Theta\Phi$

abstract: Simulating the incommensurate spin density waves (ISDW) states is not a simple task within the standard \emph{ab-initio} methods. Moreover, in the context of new material discovery, there is a need for fast and reliable tool capable to scan and optimize the total energy as a function of the pitch vector, thus allowing to automatize the search for new materials. In this paper we show how the ISDW can be efficiently obtained within the recently released $\Theta\Phi$ program. We illustrate this on an example of the single orbital Hubbard model and of $\gamma$-Fe, where the ISDW emerge within the mean-field approximation and by using the twisted boundary conditions. We show the excellent agreement of the $\Theta\Phi$ with the previously published ones and discuss possible extensions. Finally, we generalize the previously given framework for spin quantization axis rotation to the most general case of spin-dependent hopping matrix elements.


title: On the Geometric Orbit Property for Lorentz Manifolds

abstract: The geodesic orbit property has been studied intensively for Riemannian manifolds. Geodesic orbit spaces are homogeneous and allow simplifications of many structural questions using the Lie algebra of the isometry group. Weakly symmetric Riemannian manifolds are geodesic orbit spaces. Here we define "naturally reductive" for pseudo-Riemannian manifolds and note that they are geodesic orbit spaces. A few years ago two of the authors proved that weakly symmetric pseudo-Riemannian manifolds are geodesic orbit spaces. In particular these results apply to pseudo-Riemannian Lorentz manifolds. There our main results are Theorems 4.2 and 5.1. In the Riemannian case the nilpotent isometry group for a geodesic orbit nilmanifold is abelian or $2$-step nilpotent. Examples show that this fails dramatically in the pseudo-Riemannian case. Here we concentrate on the geodesic orbit property for Lorentz nilmanifolds $G/H$ with $G = N \rtimes H$ and $N$ nilpotent. When the metric is nondegenerate on $[\mathfrak{n},\mathfrak{n}]$, Theorem 4.2 shows that $N$ either is at most $2$-step nilpotent as in the Riemannian situation, or is $4$-step nilpotent, but cannot be $3$-step nilpotent. Examples show that these bounds are the best possible. Surprisingly, Theorem 5.1 shows that $N$ is at most $2$-step nilpotent when the metric is degenerate on $[\mathfrak{n},\mathfrak{n}]$. Both theorems give additional structural information and specialize to naturally reductive and to weakly symmetric Lorentz nilmanifolds. Key Words: Geodesic Orbit Space; Lorentz nilmanifold; Weakly Symmetric Space; Naturally Reductive Space; Pseudo-Riemannian Manifold.


title: Tamarkin-Tsygan Calculus and Chiral Poisson Cohomology

abstract: We construct and study some vertex theoretic invariants associated to Poisson varieties, specialising in the conformal weight $0$ case to the familiar package of Poisson homology and cohomology. In order to do this conceptually we sketch a version of the \emph{calculus}, in the sense of \cite{TT}, adapted to the context of vertex algebras. We obtain the standard theorems of Poisson (co)homology in this \emph{chiral} context. This is part of a larger project related to promoting non-commutative geometric structures to chiral versions of such.


title: Text Augmentation for Language Models in High Error Recognition Scenario

abstract: We examine the effect of data augmentation for training of language models for speech recognition. We compare augmentation based on global error statistics with one based on per-word unigram statistics of ASR errors and observe that it is better to only pay attention the global substitution, deletion and insertion rates. This simple scheme also performs consistently better than label smoothing and its sampled variants. Additionally, we investigate into the behavior of perplexity estimated on augmented data, but conclude that it gives no better prediction of the final error rate. Our best augmentation scheme increases the absolute WER improvement from second-pass rescoring from 1.1 % to 1.9 % absolute on the CHiMe-6 challenge.


title: Exploring the Value of Personalized Word Embeddings

abstract: In this paper, we introduce personalized word embeddings, and examine their value for language modeling. We compare the performance of our proposed prediction model when using personalized versus generic word representations, and study how these representations can be leveraged for improved performance. We provide insight into what types of words can be more accurately predicted when building personalized models. Our results show that a subset of words belonging to specific psycholinguistic categories tend to vary more in their representations across users and that combining generic and personalized word embeddings yields the best performance, with a 4.7% relative reduction in perplexity. Additionally, we show that a language model using personalized word embeddings can be effectively used for authorship attribution.


title: Forecasting Emergency Department Capacity Constraints for COVID Isolation Beds

abstract: Predicting patient volumes in a hospital setting is a well-studied application of time series forecasting. Existing tools usually make forecasts at the daily or weekly level to assist in planning for staffing requirements. Prompted by new COVID-related capacity constraints placed on our pediatric hospital's emergency department, we developed an hourly forecasting tool to make predictions over a 24 hour window. These forecasts would give our hospital sufficient time to be able to martial resources towards expanding capacity and augmenting staff (e.g. transforming wards or bringing in physicians on call). Using Gaussian Process Regressions (GPRs), we obtain strong performance for both point predictions (average R-squared: 82%) as well as classification accuracy when predicting the ordinal tiers of our hospital's capacity (average precision/recall: 82%/74%). Compared to traditional regression approaches, GPRs not only obtain consistently higher performance, but are also robust to the dataset shifts that have occurred throughout 2020. Hospital stakeholders are encouraged by the strength of our results, and we are currently working on moving our tool to a real-time setting with the goal of augmenting the capabilities of our healthcare workers.


title: Empirically Characterizing the Buffer Behaviour of Real Devices

abstract: All the routers include a buffer in order to enqueue packets waiting to be transmitted. The behaviour of the routers' buffer is of primary importance when studying network traffic, since it may modify some characteristics, as delay or jitter, and may also drop packets. As a consequence, the characterization of this buffer is interesting, especially when real-time flows are being transmitted: if the buffer characteristics are known, then different techniques can be used so as to adapt the traffic: multiplexing a number of small packets into a big one, fragmentation, etc. This work presents a preliminary study of how to determine the technical and functional characteristics of the buffer of a certain device (as e.g. behaviour, size, limits, input and output rate), or even in a remote Internet network node. Two different methodologies are considered, and tested on two real scenarios which have been implemented; real measurements permit the estimation of the buffer size, and the input and output rates, when there is physical or remote access to the "System Under Test". In case of having physical access, the maximum number of packets in the queue can be determined by counting. In contrast, if the node is remote, its buffer size has to be estimated. We have obtained accurate results in wired and wireless networks.


title: Forecasting and Analyzing the Military Expenditure of India Using Box-Jenkins ARIMA Model

abstract: The advancement in the field of statistical methodologies to economic data has paved its path towards the dire need for designing efficient military management policies. India is ranked as the third largest country in terms of military spender for the year 2019. Therefore, this study aims at utilizing the Box-Jenkins ARIMA model for time series forecasting of the military expenditure of India in forthcoming times. The model was generated on the SIPRI dataset of Indian military expenditure of 60 years from the year 1960 to 2019. The trend was analysed for the generation of the model that best fitted the forecasting. The study highlights the minimum AIC value and involves ADF testing (Augmented Dickey-Fuller) to transform expenditure data into stationary form for model generation. It also focused on plotting the residual error distribution for efficient forecasting. This research proposed an ARIMA (0,1,6) model for optimal forecasting of military expenditure of India with an accuracy of 95.7%. The model, thus, acts as a Moving Average (MA) model and predicts the steady-state exponential growth of 36.94% in military expenditure of India by 2024.


title: A Framework for Mediation Analysis with Multiple Exposures, Multivariate Mediators, and Non-Linear Response Models

abstract: Mediation analysis seeks to identify and quantify the paths by which an exposure affects an outcome. Intermediate variables which are effected by the exposure and which effect the outcome are known as mediators. There exists extensive work on mediation analysis in the context of models with a single mediator and continuous and binary outcomes. However these methods are often not suitable for multi-omic data that include highly interconnected variables measuring biological mechanisms and various types of outcome variables such as censored survival responses. In this article, we develop a general framework for causal mediation analysis with multiple exposures, multivariate mediators, and continuous, binary, and survival responses. We estimate mediation effects on several scales including the mean difference, odds ratio, and restricted mean scale as appropriate for various outcome models. Our estimation method avoids imposing constraints on model parameters such as the rare disease assumption while accommodating continuous exposures. We evaluate the framework and compare it to other methods in extensive simulation studies by assessing bias, type I error and power at a range of sample sizes, disease prevalences, and number of false mediators. Using Kidney Renal Clear Cell Carcinoma data from The Cancer Genome Atlas, we identify proteins which mediate the effect of metabolic gene expression on survival. Software for implementing this unified framework is made available in an R package (https://github.com/longjp/mediateR).


title: Rank-Based Testing for Semiparametric VAR Models: a measure transportation approach

abstract: We develop a class of tests for semiparametric vector autoregressive (VAR) models with unspecified innovation densities, based on the recent measure-transportation-based concepts of multivariate {\it center-outward ranks} and {\it signs}. We show that these concepts, combined with Le Cam's asymptotic theory of statistical experiments, yield novel testing procedures, which (a)~are valid under a broad class of innovation densities (possibly non-elliptical, skewed, and/or with infinite moments), (b)~are optimal (locally asymptotically maximin or most stringent) at selected ones, and (c) are robust against additive outliers. In order to do so, we establish a H\' ajek asymptotic representation result, of independent interest, for a general class of center-outward rank-based serial statistics. As an illustration, we consider the problems of testing the absence of serial correlation in multiple-output and possibly non-linear regression (an extension of the classical Durbin-Watson problem) and the sequential identification of the order $p$ of a vector autoregressive (VAR($p$)) model. A Monte Carlo comparative study of our tests and their routinely-applied Gaussian competitors demonstrates the benefits (in terms of size, power, and robustness) of our methodology; these benefits are particularly significant in the presence of asymmetric and leptokurtic innovation densities. A real data application concludes the paper.


title: H-chromatic symmetric functions

abstract: We introduce $H$-chromatic symmetric functions, $X_{G}^{H}$, which use the $H$-coloring of a graph $G$ to define a generalization of Stanley's chromatic symmetric functions. We say two graphs $G_1$ and $G_2$ are $H$-chromatically equivalent if $X_{G_1}^{H} = X_{G_2}^{H}$, and use this idea to study uniqueness results for $H$-chromatic symmetric functions, with a particular emphasis on the case $H$ is a complete bipartite graph. We also show that several of the classical bases of the space of symmetric functions, i.e. the monomial symmetric functions, power sum symmetric functions, and elementary symmetric functions, can be realized as $H$-chromatic symmetric functions. We end with some conjectures and open problems.


title: Non-local Optimization: Imposing Structure on Optimization Problems by Relaxation

abstract: In stochastic optimization, particularly in evolutionary computation and reinforcement learning, the optimization of a function $f: \Omega \to \mathbb{R}$ is often addressed through optimizing a so-called relaxation $\theta \in \Theta \mapsto \mathbb{E}_\theta(f)$ of $f$, where $\Theta$ resembles the parameters of a family of probability measures on $\Omega$. We investigate the structure of such relaxations by means of measure theory and Fourier analysis, enabling us to shed light on the success of many associated stochastic optimization methods. The main structural traits we derive and that allow fast and reliable optimization of relaxations are the consistency of optimal values of $f$, Lipschitzness of gradients, and convexity. We emphasize settings where $f$ itself is not differentiable or convex, e.g., in the presence of (stochastic) disturbance.


title: Surface magnetic structure investigation of a nanolaminated Mn$_2$GaC thin film using a magnetic field microscope based on Nitrogen-Vacancy centers

abstract: This work presents a magnetic field imaging method based on color centers in diamond crystal applied to a thin film of a nanolaminated Mn$_2$GaC MAX phase. Magnetic properties of the surface related structures have been described around the first order transition at 214 K by performing measurements in the temperature range between 200 K and 235 K with the surface features fading out by increasing temperature above the transition temperature. The results presented here demonstrate how Nitrogen-Vacancy center based magnetic microscopy can supplement the traditionally used set of experimental techniques, giving additional information of microscopic scale magnetic field features, and allowing to investigate the temperature dependent magnetic behavior. The additional information acquired in this way is relevant to applications where surface magnetic properties are of essence.


title: Quality of internal representation shapes learning performance in feedback neural networks

abstract: A fundamental feature of complex biological systems is the ability to form feedback interactions with their environment. A prominent model for studying such interactions is reservoir computing, where learning acts on low-dimensional bottlenecks. Despite the simplicity of this learning scheme, the factors contributing to or hindering the success of training in reservoir networks are in general not well understood. In this work, we study non-linear feedback networks trained to generate a sinusoidal signal, and analyze how learning performance is shaped by the interplay between internal network dynamics and target properties. By performing exact mathematical analysis of linearized networks, we predict that learning performance is maximized when the target is characterized by an optimal, intermediate frequency which monotonically decreases with the strength of the internal reservoir connectivity. At the optimal frequency, the reservoir representation of the target signal is high-dimensional, de-synchronized, and thus maximally robust to noise. We show that our predictions successfully capture the qualitative behaviour of performance in non-linear networks. Moreover, we find that the relationship between internal representations and performance can be further exploited in trained non-linear networks to explain behaviours which do not have a linear counterpart. Our results indicate that a major determinant of learning success is the quality of the internal representation of the target, which in turn is shaped by an interplay between parameters controlling the internal network and those defining the task.


title: Fractionally Integrated Moving Average Stable Processes With Long-Range Dependence

abstract: Long memory processes driven by L\'evy noise with finite second-order moments have been well studied in the literature. They form a very rich class of processes presenting an autocovariance function which decays like a power function. Here, we study a class of L\'evy process whose second-order moments are infinite, the so-called $\alpha$-stable processes. Based on Samorodnitsky and Taqqu (2000), we construct an isometry that allows us to define stochastic integrals concerning the linear fractional stable motion using Riemann-Liouville fractional integrals. With this construction, follows naturally an integration by parts formula. We then present a family of stationary $S\alpha S$ processes with the property of long-range dependence, using a generalized measure to investigate its dependence structure. In the end, the law of large number's result for a time's sample of the process is shown as an application of the isometry and integration by parts formula.


title: Three-Body Baryonic $ \bar{B}^{0}\rightarrow\Lambda^+_c \bar{p} \pi^0$ Decay

abstract: We study the three-body baryonic decay $ \bar{B}^{0}\rightarrow\Lambda^+_c \bar{p} \pi^0$ based on the factorization approach. The most important Feynman diagrams for this decay mode are factorizable and non-factorizable contributions, so that the former plays the key role in the decay amplitude and the latter affects the amplitude coefficient $a_2$. We calculate the decay branching ratio for some values of $a_2$ in the range of $0.16-0.40$ and we find that the value of $a_2=0.17$ leads us to the result of $(1.50\pm0.21)\times10^{-4}$ which evidently agrees with the experimental result $(1.55\pm0.17\pm0.08)\times 10^{-4}$ reported by BABAR


title: Ecole: A Gym-like Library for Machine Learning in Combinatorial Optimization Solvers

abstract: We present Ecole, a new library to simplify machine learning research for combinatorial optimization. Ecole exposes several key decision tasks arising in general-purpose combinatorial optimization solvers as control problems over Markov decision processes. Its interface mimics the popular OpenAI Gym library and is both extensible and intuitive to use. We aim at making this library a standardized platform that will lower the bar of entry and accelerate innovation in the field. Documentation and code can be found at https://www.ecole.ai.


title: Quantifying and Learning Linear Symmetry-Based Disentanglement

abstract: The definition of Linear Symmetry-Based Disentanglement (LSBD) formalizes the notion of linearly disentangled representations, but there is currently no metric to quantify LSBD. Such a metric is crucial to evaluate LSBD methods and to compare to previous understandings of disentanglement. We propose $\mathcal{D}_\mathrm{LSBD}$, a mathematically sound metric to quantify LSBD, and provide a practical implementation for $\mathrm{SO}(2)$ groups. Furthermore, from this metric we derive LSBD-VAE, a semi-supervised method to learn LSBD representations. We demonstrate the utility of our metric by showing that (1) common VAE-based disentanglement methods don't learn LSBD representations, (2) LSBD-VAE as well as other recent methods can learn LSBD representations, needing only limited supervision on transformations, and (3) various desirable properties expressed by existing disentanglement metrics are also achieved by LSBD representations.


title: Disentangling the multi-phase circumgalactic medium shared between a dwarf and a massive star-forming galaxy at z~0.4

abstract: The multi-phase circumgalactic medium (CGM) arises within the complex environment around a galaxy, or collection of galaxies, and possibly originates from a wide range of physical mechanisms. In this paper, we attempt to disentangle the origins of these multi-phase structures and present a detailed analysis of the quasar field Q0122-003 field using Keck/KCWI galaxy observations and HST/COS spectra probing the CGM. Our re-analysis of this field shows that there are two galaxies associated with the absorption. We have discovered a dwarf galaxy, G_27kpc ($M_{\star}=10^{8.7}$ M$_{\odot}$), at z=0.39863 that is 27 kpc from the quasar sightline. G_27kpc is only +21 km/s from a more massive ($M_{\star}=10^{10.5}$ M$_{\odot}$) star-forming galaxy, G_163kpc, at an impact parameter of 163 kpc. While G_163kpc is actively forming stars (SFR=6.9 M$_{\odot}$ yr$^{-1}$), G_27kpc has a low star-formation rate (SFR=$0.08\pm0.03$ M$_{\odot}$ yr$^{-1}$) and star formation surface density ($\Sigma_{SFR}=0.006$ M$_{\odot}$ kpc$^{-2}$ yr$^{-1}$), implying no active outflows. By comparing galaxy SFRs, kinematics, masses and distances from the quasar sightline to the absorption kinematics, column densities and metallicities, we have inferred the following: (1) Part of the low-ionization phase has a metallicity and kinematics consistent with being accreted onto G_27kpc. (2) The remainder of the low ionization phase has metallicities and kinematics consistent with being intragroup gas being transferred from G_27kpc to G_163kpc. (3) The high ionization phase is consistent with being produced solely by outflows originating from the massive halo of G_163kpc. Our results demonstrate the complex nature of the multi-phase CGM, especially around galaxy groups, and that detailed case-by-case studies are critical for disentangling its origins.


title: A high-order / low-order (HOLO) algorithm for preserving conservation in time-dependent low-rank transport calculations

abstract: Dynamical low-rank (DLR) approximation methods have previously been developed for time-dependent radiation transport problems. One crucial drawback of DLR is that it does not conserve important quantities of the calculation, which limits the applicability of the method. Here we address this conservation issue by solving a low-order equation with closure terms computed via a high-order solution calculated with DLR. We observe that the high-order solution well approximates the closure term, and the low-order solution can be used to correct the conservation bias in the DLR evolution. We also apply the linear discontinuous Galerkin method for the spatial discretization to obtain the asymptotic limit. We then demonstrate with the numerical results that this so-called high-order / low-order (HOLO) algorithm is conservative without sacrificing computational efficiency and accuracy.


title: Imaging Spectroscopy of CME-Associated Solar Radio Bursts

abstract: We present first results of a solar radio event observed with the Owens Valley Radio Observatory Long Wavelength Array (OVRO-LWA) at metric wavelengths. We examine a complex event consisting of multiple radio sources/bursts associated with a fast coronal mass ejection (CME) and an M2.1 GOES soft X-ray flare from 2015 September 20. Images of 9--s cadence are used to analyze the event over a 120-minute period, and solar emission is observed out to a distance of $\approx3.5\,R_\odot$, with an instantaneous bandwidth covering 22~MHz within the frequency range of 40--70~MHz. We present our results from the investigation of the radio event, focusing particularly on one burst source that exhibits outward motion, which we classify as a moving type IV burst. We image the event at multiple frequencies and use the source centroids to obtain the velocity for the outward motion. Spatial and temporal comparison with observations of the CME in white light from the LASCO(C2) coronagraph, indicates an association of the outward motion with the core of the CME. By performing graduated-cylindrical-shell (GCS) reconstruction of the CME, we constrain the density in the volume. The electron plasma frequency obtained from the density estimates do not allow us to completely dismiss plasma emission as the underlying mechanism. However, based on source height and smoothness of the emission in frequency and time, we argue that gyrosynchrotron is the more plausible mechanism. We use gyrosynchrotron spectral fitting techniques to estimate the evolving physical conditions during the outward motion of this burst source.


title: Gillet descent for connective K-theory

abstract: Using Gillet's technique of projective envelopes, we prove a homological descent theorem for the connective K-homology of schemes.


title: Domain Wall Leaky Integrate-and-Fire Neurons with Shape-Based Configurable Activation Functions

abstract: Complementary metal oxide semiconductor (CMOS) devices display volatile characteristics, and are not well suited for analog applications such as neuromorphic computing. Spintronic devices, on the other hand, exhibit both non-volatile and analog features, which are well-suited to neuromorphic computing. Consequently, these novel devices are at the forefront of beyond-CMOS artificial intelligence applications. However, a large quantity of these artificial neuromorphic devices still require the use of CMOS, which decreases the efficiency of the system. To resolve this, we have previously proposed a number of artificial neurons and synapses that do not require CMOS for operation. Although these devices are a significant improvement over previous renditions, their ability to enable neural network learning and recognition is limited by their intrinsic activation functions. This work proposes modifications to these spintronic neurons that enable configuration of the activation functions through control of the shape of a magnetic domain wall track. Linear and sigmoidal activation functions are demonstrated in this work, which can be extended through a similar approach to enable a wide variety of activation functions.

categories: cs.NE cond-mat.mes-hall cs.ET physics.app-ph


title: Face Dimensions of General-Purpose Cutting Planes for Mixed-Integer Linear Programs

abstract: Cutting planes are a key ingredient to successfully solve mixed-integer linear programs. For specific problems, their strength is often theoretically assessed by showing that they are facet-defining for the corresponding mixed-integer hull. In this paper we experimentally investigate the dimensions of faces induced by general-purpose cutting planes generated by a state-of-the-art solver. Therefore, we relate the dimension of each cutting plane to its impact in a branch-and-bound algorithm.


title: Temporal Splitting algorithms for non-stationary multiscale problems

abstract: In this paper, we study temporal splitting algorithms for multiscale problems. The exact fine-grid spatial problems typically require some reduction in degrees of freedom. Multiscale algorithms are designed to represent the fine-scale details on a coarse grid and, thus, reduce the problems' size. When solving time-dependent problems, one can take advantage of the multiscale decomposition of the solution and perform temporal splitting by solving smaller-dimensional problems, which is studied in the paper. In the proposed approach, we consider the temporal splitting based on various low dimensional spatial approximations. Because a multiscale spatial splitting gives a "good" decomposition of the solution space, one can achieve an efficient implicit-explicit temporal discretization. We present a recently developed theoretical result in our earlier work and adopt it in this paper for multiscale problems. Numerical results are presented to demonstrate the efficiency of the proposed splitting algorithm.


title: $\rm \Lambda_{c}^{+}$ production and baryon-to-meson ratios in pp and p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV at the LHC

abstract: The prompt production of the charm baryon $\rm \Lambda_{c}^{+}$ and the $\rm \Lambda_{c}^{+}/\mathrm {D^0}$ production ratios were measured at midrapidity with the ALICE detector in pp and p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$TeV. These new measurements show a clear decrease of the $\rm \Lambda_{c}^{+}/\mathrm {D^0}$ ratio with increasing transverse momentum ($p_{\rm T}$) in both collision systems in the range $2<p_{\rm T}<12$ GeV/$c$, exhibiting similarities with the light-flavour baryon-to-meson ratios ${\rm p}/\pi$ and $\Lambda/\mathrm {K^0_S}$. At low $p_{\rm T}$, predictions that include additional colour-reconnection mechanisms beyond the leading-colour approximation; assume the existence of additional higher-mass charm-baryon states; or include hadronisation via coalescence can describe the data, while predictions driven by charm-quark fragmentation processes measured in $\mathrm {e^+e^-}$ and $\mathrm {e^-p}$ collisions significantly underestimate the data. The results presented in this letter provide significant evidence that the established assumption of universality (colliding-system independence) of parton-to-hadron fragmentation is not sufficient to describe charm-baryon production in hadronic collisions at LHC energies.


title: $\mathrm{\Lambda_{c}^{+}}$ production in pp and in p-Pb collisions at $\sqrt{s_{\rm {NN}}} = 5.02$ TeV

abstract: The production cross section of prompt $\mathrm{\Lambda_{c}^{+}}$ charm baryons was measured with the ALICE detector at the LHC at midrapidity in proton-proton (pp) and proton-lead (p-Pb) collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}} = 5.02$ TeV. The $\mathrm{\Lambda_{c}^{+}}$ and $\rm {\overline{\Lambda}{}_c^-}$ baryons were reconstructed in the hadronic decay channels $\rm \Lambda_{c}^{+} \rightarrow p K^{-}\pi^{+}$ and $\rm \Lambda_{c}^{+}\to p K^{0}_{S}$ and respective charge conjugates. The measured differential cross sections as a function of transverse momentum ($p_{\rm T}$) and the $p_{\rm T}$-integrated $\mathrm{\Lambda_{c}^{+}}$ production cross section in pp and in p-Pb collisions are presented. The $\mathrm{\Lambda_{c}^{+}}$ nuclear modification factor ($R_\mathrm{pPb}$), calculated from the cross sections in pp and in p-Pb collisions, is presented and compared with the $R_\mathrm{pPb}$ of D mesons. The $\mathrm {\Lambda_{c}^{+}}/\mathrm {D^0}$ ratio is also presented and compared with the light-flavour baryon-to-meson ratios p$/\pi$ and $\Lambda /\mathrm {K^0_S}$, and measurements from other LHC experiments. The results are compared to predictions from model calculations and Monte Carlo event generators.


title: Statistical learning for change point and anomaly detection in graphs

abstract: Complex systems which can be represented in the form of static and dynamic graphs arise in different fields, e.g. communication, engineering and industry. One of the interesting problems in analysing dynamic network structures is to monitor changes in their development. Statistical learning, which encompasses both methods based on artificial intelligence and traditional statistics, can be used to progress in this research area. However, the majority of approaches apply only one or the other framework. In this paper, we discuss the possibility of bringing together both disciplines in order to create enhanced network monitoring procedures focussing on the example of combining statistical process control and deep learning algorithms. Together with the presentation of change point and anomaly detection in network data, we propose to monitor the response times of ambulance services, applying jointly the control chart for quantile function values and a graph convolutional network.


title: Single-quadrature quantum magnetometry in cavity electromagnonics

abstract: A scheme of an ultra-sensitive magnetometer in the cavity quantum electromagnonics where the intracavity microwave mode coupled to a magnonic mode via magnetic dipole interaction is proposed. It is shown that by driving both magnonic and microwave modes with external classical fields and controlling the system parameters, one can reduce the added noise of magnetic field measurement below the standard quantum limit (SQL). Surprisingly, we show that beyond the rotating wave approximation (RWA), not only the added noise can be suppressed, but also the output cavity response to the input signal can be substantially amplified in order to achieve a precise magnetic-field measurement. The estimated theoretical sensitivity of the proposed magnetic amplifier-sensor is approximately in the order of $10^{-18}T/\sqrt{Hz}$ which is competitive compared to the current state-of-the-art magnetometers like superconducting quantum interference devices (SQUIDs) and atomic magnetometers. The advantage of the proposed sensor in comparison with the other magnetometers is its high sensitivity at room temperature and sensing in a wide range of frequency up to MHz as well as its capability to signal-response amplification.


title: Controlling ferroelectric hysteresis offsets in PbTiO$_{3}$ based superlattices

abstract: Ferroelectric materials are characterized by degenerate ground states with multiple polarization directions. In a ferroelectric capacitor this should manifest as equally favourable up and down polarization states. However, this ideal behavior is rarely observed in ferroelectric thin films and superlattice devices, which generally exhibit a built-in bias which favors one polarization state over the other. Often this polarization asymmetry can be attributed to the electrodes. In this study we examine bias in PbTiO$_3$-based ferroelectric superlattices that is not due to the electrodes, but rather to the nature of the defects that form at the interfaces during growth. Using a combination of experiments and first-principles simulations, we are able to explain the sign of the observed built-in bias and its evolution with composition. Our insights allow us to design devices with zero built-in bias by controlling the composition and periodicity of the superlattices.


title: A NUT Charge Weak Gravity Conjecture from Dimensional Reduction

abstract: We analyze the constraints on four-derivative corrections to 5d Einstein-Maxwell theory from the black hole Weak Gravity Conjecture (WGC). We calculate the leading corrections to the extremal mass of asymptotically flat 5d charged solutions as well as 4d Kaluza-Klein compactifications. The WGC bounds from the latter, interpreted as 4d dyonic black holes, are found to be strictly stronger. As magnetic graviphoton charge lifts to a NUT-like charge in 5d, we argue that the logic of the WGC should apply to these topological charges as well and leads to new constraints on purely gravitational theories.


title: Attractive and repulsive dipolar interaction in bilayers of indirect excitons

abstract: We explore attractive dipolar interaction in indirect excitons (IXs). For one layer of IXs in a single pair of coupled quantum wells (CQW), the out-of-plane IX electric dipoles lead to repulsive dipolar interaction between IXs. The attractive dipolar interaction between IXs is realized in a 2-CQW heterostructure with two IX layers in two separated CQW pairs. We found both in experimental measurements and theoretical simulations that increasing density of IXs in one layer causes a monotonic energy reduction for IXs in the other layer. We also found an in-plane shift of a cloud of IXs in one layer towards a cloud of IXs in the other layer. This behaviour is qualitatively consistent with attractive dipolar interaction. The measured IX energy reduction and IX cloud shift are higher than the values given by the correlated liquid theory.


title: Analysis of an initial value problem for an extracellular and intracellular model of hepatitis C virus infection

abstract: In this paper, a mathematical analysis of the global dynamics of a viral infection model in vivo is carried out. We study the dynamics of a hepatitis C virus (HCV) model, under therapy, that considers both extracellular and intracellular levels of infection. At present most mathematical modeling of viral kinetics after treatment only addresses the process of infection of a cell by the virus and the release of virions by the cell, while the processes taking place inside the cell are not included. We prove that the solutions of the new model with positive initial values are positive, exist globally in time and are bounded. The model has two virus-free steady states. They are distinguished by the fact that viral RNA is absent inside the cells in the first state and present inside the cells in the second. There are basic reproduction numbers associated to each of these steady states. If the basic reproduction number of the first steady state is less than one then that state is asymptotically stable. If the basic reproduction number of the first steady state is greater than one and that of the second less than one then the second steady state is asymptotically stable. If both basic reproduction numbers are greater than one then we obtain various conclusions which depend on different restrictions on the parameters of the model. Under increasingly strong assumptions we prove that there is at least one positive steady state (infected equilibrium), that there is a unique positive steady state and that the positive steady state is stable. We also give a condition under which every positive solution converges to a positive steady state. This is proved by methods of Li and Muldowney. Finally, we illustrate the theoretical results by numerical simulations.


title: The Bayesian Committee Approach for Computational Physics Problems

abstract: In this work, we propose a method for efficient learning of a multi-dimensional function. This method combines the Bayesian neural networks and the query-by-committee method. A committee made of deep Bayesian neural networks not only can provide uncertainty of the prediction but also can provide the discrepancy between committee members. Both the uncertainty and the discrepancy are large in the regions where the target function varies rapidly, and therefore, both quantities can be used to guide sampling data to such regions. In this way, we can learn a function accurately with the number of queried data points much less than uniform sampling. Here we test our method with two examples. One example is to find a rare phase in a phase diagram, which is separated from other phases by a second-order phase transition. In this example, the target function is the susceptibility function, and since the divergence of the susceptibility function locates the phase diagram, the task of searching such a phase perfectly matches the advantage of our method. Another example is to learn the distribution function for Monte Carlo integration of a high-dimensional function. In both examples, we show that our method performs significantly efficiently than uniform sampling. Our method can find broad applications in computational scientific problems.

categories: physics.comp-ph cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.quant-gas


title: Evoking Places from Spaces. The application of multimodal narrative techniques in the creation of "U Modified"

abstract: Multimodal diegetic narrative tools, as applied in multimedia arts practices, possess the ability to cross the spaces that exist between the physical world and the imaginary. Within this paper we present the findings of a multidiscipline practice based research project that explored the potential of an audiovisual art performance to purposefully interact with an audience's perception of narrative place. To achieve this goal, research was undertaken to investigate the function of multimodal diegetic practices as applied in the context of a sonic art narrative. This project direction was undertaken to facilitate the transformation of previous experiences of place through the creative amalgamation and presentation of collected audio and visual footage from real world spaces. Through the presentation of multimedia relating to familiar geographical spatial features, the audience were affected to evoke memories of place and to construct and manipulate their own narrative.


title: The Key Features of the CoCalc Cloud Service Use in the Process of Mathematics and Science Teachers Training

abstract: The article describes the features of the evolution of the cloud-based systems in the context of mathematics and science teachers training. The notion of the systems of computer mathematics (SCM) is considered as the special kind of learning software. The lists of the special mathematical operations, functions and methods used within the structure of modern SCM are examined. Web-SCM properties are considered as general and specific. The cloud-based SCM are taken as the special stage of Web-SCM development. The criteria for an appropriate cloud service selection to support the process of mathematics disciplines learning are revealed. The key features of the CoCalc cloud service use are considered in concern to the proposed selection criteria: the available computing capacities, the diversity of tools for training and its control, the possibility of increasing the computing capability, the openness of the program code, the availability of tools for joint editing. The results of the experimental testing of CoCalc use in the process of mathematics teachers training are described. Besides, the discussion covering the issues of integration and compatibility with the other tools and services used for mathematics learning (on the example of Octave) is presented. The main trends of the CoCalc cloud service introduction into the training process and the prospects for further research are revealed.


title: Continuous Perception for Classifying Shapes and Weights of Garmentsfor Robotic Vision Applications

abstract: We present an approach to continuous perception for robotic laundry tasks. Our assumption is that the visual prediction of a garment's shapes and weights is possible via a neural network that learns the dynamic changes of garments from video sequences. Continuous perception is leveraged during training by inputting consecutive frames, of which the network learns how a garment deforms. To evaluate our hypothesis, we captured a dataset of 40K RGB and 40K depth video sequences while a garment is being manipulated. We also conducted ablation studies to understand whether the neural network learns the physical and dynamic properties of garments. Our findings suggest that a modified AlexNet-LSTM architecture has the best classification performance for the garment's shape and weights. To further provide evidence that continuous perception facilitates the prediction of the garment's shapes and weights, we evaluated our network on unseen video sequences and computed the 'Moving Average' over a sequence of predictions. We found that our network has a classification accuracy of 48% and 60% for shapes and weights of garments, respectively.


title: An approximation to determine the source of the WOW! Signal

abstract: In this paper it is analysed which of the thousands of stars in the WOW! Signal region could have the highest chance of being the real source of the signal, providing that it came from a star system similar to ours. A total of 66 G and K-type stars are sampled, but only one of them is identified as a potential Sun-like star considering the available information in the Gaia Archive. This candidate source, which is named 2MASS 19281982-2640123, therefore becomes an ideal target to conduct observations in the search for potentially habitable exoplanets. Another two candidate stars have a luminosity error interval that includes the luminosity of the Sun, and 14 candidates more are also identified as potential Sun-like stars, but the estimations on their luminosity were unknown.


title: Infinite degeneracy of Landau levels from the Euclidean symmetry with central extension revisited

abstract: The planar Landau system which describes the quantum mechanical motion of a charged particle in a plane with a uniform magnetic field perpendicular to the plane, is explored within pedagogical settings aimed at the beginning graduate level. The system is known to possess the Euclidean symmetry in two dimensions with central extension $\bar{E}(2)$. In this paper, we revisit the well-known energy eigenvalues of the system, known as the Landau levels, by exploiting the related $\bar{e}(2)$ symmetry algebra. Specifically, we utilize the Casimir operator and the commutation relations of the generators of the $\bar{E}(2)$ group. More importantly, an algebraic formalism on this topic based on Schwinger's oscillator model of angular momentum is also presented. The dimensions of irreducible representations of the $\bar{E}(2)$ group and their implications on the degeneracy of Landau levels is discussed.


title: 40 Gbps Readout interface STARE for the AGATA Project

abstract: The Advanced GAmma Tracking Array (AGATA) multi detector spectrometer will provide precise information for the study of the properties of the exotic nuclear matter (very unbalanced proton (Z) and neutron (N) numbers) along proton- and neutron- drip lines and of super-heavy nuclei. This is done using the latest technology of particle accelerators. The AGATA spectrometer consists of 180 high purity Germanium detectors. Each detector is segmented into 38 segments. The very harsh project requirements are to measure gamma ray energies with very high resolution (< 1x 10 -3) at a high detector counting rate (50 Kevents / sec / crystal). This results in a very high data transfer rate per crystal (5 to 8 Gbps). The 38 segments are sampled @ 100 MHz with 14 bits of resolution. The samples are continuously transferred to the CAP module which reduces the data rate from 64 Gbps to 5 Gbps. The CAP module also adds continuous monitoring data which results in total outgoing data rate of 10 Gbps. The STARE module is designed to fit between the CAP module and the computer farm. It will package the data from the CAP module and transmit it to the server farm using a 10 Gbps UDP connection with a delivery insurance mechanism implemented to ensure that all data is transferred.


title: Non-equilibrium heating path for utrafast laser-induced nucleation of skyrmion lattices

abstract: We explore the helicity-independent light-induced nucleation of skyrmion lattices in ferromagnetic cobalt-based trilayers with perpendicular magnetic anisotropy. Using Atomistic Spin Dynamics simulations, we show that a high temperature excitation followed by magnon drops and their non-equilibrium relaxation, accessed by an ultrafast laser excitation with specific duration and intensity, can lead to the generation of a skyrmion lattice stable at room temperature. The nucleation window, the topological density and the skyrmion polarity can be additionally manipulated by external magnetic fields. Our results provide insight into the non-equilibrium nature of skyrmionic excitations at non-zero temperatures and pave additional routes for their use in information technologies.


title: Guiding user annotations for units-of-measure verification

abstract: This extended abstract reports on previous work of the CamFort project in which we developed an external units-of-measure type system for Fortran code, targeted at scientists. Our approach can guide the programmer in adding specifications (type annotations) to existing code, with the aim of easing adoption on legacy code. Pertinent to the topics of the HATRA workshop, we discuss the human-aspects of the tool here. CamFort is open-source and freely available online.


title: Space-based Measurements of Neutron Lifetime: Approaches to Resolving the Neutron Lifetime Anomaly

abstract: Free neutrons have a measured lifetime of 880 s, but disagreement between existing laboratory measurements of ~10 s have persisted over many years. This uncertainty has implications for multiple physics disciplines, including standard-model particle physics and Big-Bang nucleosynthesis. Space-based neutron lifetime measurements have been shown to be feasible using existing data taken at Venus and the Moon, although the uncertainties for these measurements of tens of seconds prevent addressing the current lifetime discrepancy. We investigate the implementation of a dedicated space-based experiment that could provide a competitive and independent lifetime measurement. We considered a variety of scenarios, including measurements made from orbit about the Earth, Moon, and Venus, as well as on the surface of the Moon. For a standard-sized neutron detector, a measurement with three-second statistical precision can be obtained from Venus orbit in less than a day; a one-second statistical precision can be obtained from Venus orbit in less than a week. Similarly precise measurements in Earth orbit and on the lunar surface can be acquired in less than 40 days (three-second precision) and ~300 days (one-second precision). Systematic uncertainties that affect a space-based neutron lifetime measurement are investigated, and the feasibility of developing such an experiment is discussed.


title: Laser-driven neutron source from high temperature D-D fusion reactions

abstract: We report a laser-driven neutron source with high yield ($>10^8$/J) and high peak flux ($>10^{25}$/cm$^2$/s) derived from high-temperature deuteron-deuteron fusion reactions. The neutron yield and the fusion temperature ($\sim 200$ keV) in our experiment are respectively two orders of magnitude and one order of magnitude higher than any previous laser-induced D-D fusion reaction. The high-temperature plasma is generated from thin ($\sim 2\,\mu$m), solid-density deuterium targets, produced by a cryogenic jet, irradiated by a 140 fs, 130 J petawatt laser with an F/3 off-axis parabola and a plasma mirror achieving fast volumetric heating of the target. The fusion temperature and neutron fluxes achieved here suggest future laser experiments can take advantage of neutrons to diagnose the plasma conditions and come closer to laboratory study of astrophysically-relevant nuclear physics.


title: Subleading non linear gravitational memory effect

abstract: We apply the new method based on null geodesics for detecting gravitational memory and find the bulk memory in Newman-Unti gauge around the boundary of the conformally compactified space time. We show how we use the newly found conserved charges in the subleading orders of large-$r$ expansion of the BMS charges to define the gravitational memory at each order in the non-linearised gravitational theory. We also find the gravitational shift in the $r$ direction. It is shown that the longitudinal displacement at order $1/r$ is the relative radius change between two detectors derived by Strominger and Zhiboedov.


title: Zero-Temperature Coarsening in the Two-Dimensional Long-Range Ising Model

abstract: We investigate the nonequilibrium dynamics following a quench to zero temperature of the non-conserved Ising model with power-law decaying long-range interactions $\propto 1/r^{d+\sigma}$ in $d=2$ spatial dimensions. The zero-temperature coarsening is always of special interest among nonequilibrium processes, because often peculiar behavior is observed. We provide estimates of the nonequilibrium exponents, viz., the growth exponent $\alpha$, the persistence exponent $\theta$, and the fractal dimension $d_f$. It is found that the growth exponent $\alpha\approx 3/4$ is independent of $\sigma$ and different from $\alpha=1/2$ as expected for nearest-neighbor models. In the large $\sigma$ regime of the tunable interactions only the fractal dimension $d_f$ of the nearest-neighbor Ising model is recovered, while the other exponents differ significantly. For the persistence exponent $\theta$ this is a direct consequence of the different growth exponents $\alpha$ as can be understood from the relation $d-d_f=\theta/\alpha$; they just differ by the ratio of the growth exponents $\approx 3/2$. This relation has been proposed for annihilation processes and later numerically tested for the $d=2$ nearest-neighbor Ising model. We confirm this relation for all $\sigma$ studied, reinforcing its general validity.


title: Deep Learning-Based CSI Feedback for Beamforming in Single- and Multi-cell Massive MIMO Systems

abstract: The potentials of massive multiple-input multiple-output (MIMO) are all based on the available instantaneous channel state information (CSI) at the base station (BS). Therefore, the user in frequency-division duplexing (FDD) systems has to keep on feeding back the CSI to the BS, thereby occupying large uplink transmission resources. Recently, deep learning (DL) has achieved great success in the CSI feedback. However, the existing works just focus on improving the feedback accuracy and ignore the effects on the following modules, e.g., beamforming (BF). In this paper, we propose a DL-based CSI feedback framework for BF design, called CsiFBnet. The key idea of the CsiFBnet is to maximize the BF performance gain rather than the feedback accuracy. We apply it to two representative scenarios: single- and multi-cell systems. The CsiFBnet-s in the single-cell system is based on the autoencoder architecture, where the encoder at the user compresses the CSI and the decoder at the BS generates the BF vector. The CsiFBnet-m in the multicell system has to feed back two kinds of CSI: the desired and the interfering CSI. The entire neural networks are trained by an unsupervised learning strategy. Simulation results show the great performance improvement and complexity reduction of the CsiFBnet compared with the conventional DL-based CSI feedback methods.


title: Exploring Gender Disparities in Time to Diagnosis

abstract: Sex and gender-based healthcare disparities contribute to differences in health outcomes. We focus on time to diagnosis (TTD) by conducting two large-scale, complementary analyses among men and women across 29 phenotypes and 195K patients. We first find that women are consistently more likely to experience a longer TTD than men, even when presenting with the same conditions. We further explore how TTD disparities affect diagnostic performance between genders, both across and persistent to time, by evaluating gender-agnostic disease classifiers across increasing diagnostic information. In both fairness analyses, the diagnostic process favors men over women, contradicting the previous observation that women may demonstrate relevant symptoms earlier than men. These analyses suggest that TTD is an important yet complex aspect when studying gender disparities, and warrants further investigation.


title: Low-lying single-particle structure of 17C and the N = 14 sub-shell closure

abstract: The first investigation of the single-particle structure of the bound states of 17C, via the d(16C, p) transfer reaction, has been undertaken. The measured angular distributions confirm the spin-parity assignments of 1/2+ and 5/2+ for the excited states located at 217 and 335 keV, respectively. The spectroscopic factors deduced for these states exhibit a marked single-particle character, in agreement with shell model and particle-core model calculations, and combined with their near degeneracy in energy provide clear evidence for the absence of the N = 14 sub-shell closure. The very small spectroscopic factor found for the 3/2+ ground state is consistent with theoretical predictions and indicates that the {\nu}1d3/2 strength is carried by unbound states. With a dominant l = 0 valence neutron configuration and a very low separation energy, the 1/2+ excited state is a one-neutron halo candidate.


title: Improving Multimodal Accuracy Through Modality Pre-training and Attention

abstract: Training a multimodal network is challenging and it requires complex architectures to achieve reasonable performance. We show that one reason for this phenomena is the difference between the convergence rate of various modalities. We address this by pre-training modality-specific sub-networks in multimodal architectures independently before end-to-end training of the entire network. Furthermore, we show that the addition of an attention mechanism between sub-networks after pre-training helps identify the most important modality during ambiguous scenarios boosting the performance. We demonstrate that by performing these two tricks a simple network can achieve similar performance to a complicated architecture that is significantly more expensive to train on multiple tasks including sentiment analysis, emotion recognition, and speaker trait recognition.


title: Sketch and Scale: Geo-distributed tSNE and UMAP

abstract: Running machine learning analytics over geographically distributed datasets is a rapidly arising problem in the world of data management policies ensuring privacy and data security. Visualizing high dimensional data using tools such as t-distributed Stochastic Neighbor Embedding (tSNE) and Uniform Manifold Approximation and Projection (UMAP) became common practice for data scientists. Both tools scale poorly in time and memory. While recent optimizations showed successful handling of 10,000 data points, scaling beyond million points is still challenging. We introduce a novel framework: Sketch and Scale (SnS). It leverages a Count Sketch data structure to compress the data on the edge nodes, aggregates the reduced size sketches on the master node, and runs vanilla tSNE or UMAP on the summary, representing the densest areas, extracted from the aggregated sketch. We show this technique to be fully parallel, scale linearly in time, logarithmically in memory, and communication, making it possible to analyze datasets with many millions, potentially billions of data points, spread across several data centers around the globe. We demonstrate the power of our method on two mid-size datasets: cancer data with 52 million 35-band pixels from multiple images of tumor biopsies; and astrophysics data of 100 million stars with multi-color photometry from the Sloan Digital Sky Survey (SDSS).


title: FS-HGR: Few-shot Learning for Hand Gesture Recognition via ElectroMyography

abstract: This work is motivated by the recent advances in Deep Neural Networks (DNNs) and their widespread applications in human-machine interfaces. DNNs have been recently used for detecting the intended hand gesture through processing of surface electromyogram (sEMG) signals. The ultimate goal of these approaches is to realize high-performance controllers for prosthetic. However, although DNNs have shown superior accuracy than conventional methods when large amounts of data are available for training, their performance substantially decreases when data are limited. Collecting large datasets for training may be feasible in research laboratories, but it is not a practical approach for real-life applications. Therefore, there is an unmet need for the design of a modern gesture detection technique that relies on minimal training data while providing high accuracy. Here we propose an innovative and novel "Few-Shot Learning" framework based on the formulation of meta-learning, referred to as the FS-HGR, to address this need. Few-shot learning is a variant of domain adaptation with the goal of inferring the required output based on just one or a few training examples. More specifically, the proposed FS-HGR quickly generalizes after seeing very few examples from each class. The proposed approach led to 85.94% classification accuracy on new repetitions with few-shot observation (5-way 5-shot), 81.29% accuracy on new subjects with few-shot observation (5-way 5-shot), and 73.36% accuracy on new gestures with few-shot observation (5-way 5-shot).


title: Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

abstract: The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.


title: Assessment of weak-coupling approximations on a driven two-level system under dissipation

abstract: The standard weak-coupling approximations associated to open quantum systems have been extensively used in the description of a two-level quantum system, qubit, subjected to relatively weak dissipation compared with the qubit frequency. However, recent progress in the experimental implementations of controlled quantum systems with increased levels of on-demand engineered dissipation has motivated precision studies in parameter regimes that question the validity of the approximations, especially in the presence of time-dependent drive fields. In this paper, we address the precision of weak-coupling approximations by studying a driven qubit through the numerically exact and non-perturbative method known as the stochastic Liouville-von Neumann equation with dissipation. By considering weak drive fields and a cold Ohmic environment with a high cutoff frequency, we use the Markovian Lindblad master equation as a point of comparison for the SLED method and study the influence of the bath-induced energy shift on the qubit dynamics. We also propose a metric that may be used in experiments to map the regime of validity of the Lindblad equation in predicting the steady state of the driven qubit. In addition, we study signatures of the well-known Mollow triplet and observe its meltdown owing to dissipation in an experimentally feasible parameter regime of circuit electrodynamics. Besides shedding light on the practical limitations of the Lindblad equation, we expect our results to inspire future experimental research on engineered open quantum systems, the accurate modeling of which may benefit from non-perturbative methods.


title: Multipolar engineering of subwavelength dielectric particles for scattering enhancement

abstract: Electromagnetic scattering on subwavelength structures keeps attracting attention owing to abroad range of possible applications, where this phenomenon is in use. Fundamental limits of scattering cross-section, being well understood in spherical geometries, are overlooked in cases of low-symmetry resonators. Here, we revise the notion of superscattering and link this property with symmetry groups of the scattering potential. We demonstrate pathways to spectrally overlap several eigenmodes of a resonator in a way they interfere constructively and enhance the scattering cross-section. As a particular example, we demonstrate spectral overlapping of several electric and magnetic modes in a subwavelength entirely homogeneous ceramic resonator. The optimized structures show the excess of a dipolar scattering cross-section limit for a sphere up to a factor of four. The revealed rules, which link symmetry groups with fundamental scattering limits, allow performing and assessing designs of subwavelength supperscatterers, which can find a use in label-free imaging, compact antennas, long-range radio frequency identification, and many other fields.


title: An Optimal Rounding for Half-Integral Weighted Minimum Strongly Connected Spanning Subgraph

abstract: In the weighted minimum strongly connected spanning subgraph (WMSCSS) problem we must purchase a minimum-cost strongly connected spanning subgraph of a digraph. We show that half-integral linear program (LP) solutions for WMSCSS can be efficiently rounded to integral solutions at a multiplicative $1.5$ cost. This rounding matches a known $1.5$ integrality gap lower bound for a half-integral instance. More generally, we show that LP solutions whose non-zero entries are at least a value $f > 0$ can be rounded at a multiplicative cost of $2 - f$.


title: On scalar products and form factors by Separation of Variables: the antiperiodic XXZ model

abstract: We consider the XXZ spin-1/2 Heisenberg chain with antiperiodic boundary conditions. The inhomogeneous version of this model can be solved by Separation of Variables (SoV), and the eigenstates can be constructed in terms of Q-functions, solution of a Baxter TQ-equation, which have double periodicity compared to the periodic case. We compute in this framework the scalar products of a particular class of separate states which notably includes the eigenstates of the transfer matrix. We also compute the form factors of local spin operators, i.e. their matrix elements between two eigenstates of the transfer matrix. We show that these quantities admit determinant representations with rows and columns labelled by the roots of the Q-functions of the corresponding separate states, as in the periodic case, although the form of the determinant are here slightly different. We also propose alternative types of determinant representations written directly in terms of the transfer matrix eigenvalues.


title: Efficient Knowledge Distillation for RNN-Transducer Models

abstract: Knowledge Distillation is an effective method of transferring knowledge from a large model to a smaller model. Distillation can be viewed as a type of model compression, and has played an important role for on-device ASR applications. In this paper, we develop a distillation method for RNN-Transducer (RNN-T) models, a popular end-to-end neural network architecture for streaming speech recognition. Our proposed distillation loss is simple and efficient, and uses only the "y" and "blank" posterior probabilities from the RNN-T output probability lattice. We study the effectiveness of the proposed approach in improving the accuracy of sparse RNN-T models obtained by gradually pruning a larger uncompressed model, which also serves as the teacher during distillation. With distillation of 60% and 90% sparse multi-domain RNN-T models, we obtain WER reductions of 4.3% and 12.1% respectively, on a noisy FarField eval set. We also present results of experiments on LibriSpeech, where the introduction of the distillation loss yields a 4.8% relative WER reduction on the test-other dataset for a small Conformer model.


title: The Balloon-borne Investigation of Temperature and Speed of Electrons in the corona (BITSE): Mission Description and Preliminary Results

abstract: We report on the Balloonborne Investigation of Temperature and Speed of Electrons in the corona (BITSE) mission launched recently to observe the solar corona from about 3 Rs to 15 Rs at four wavelengths (393.5, 405.0, 398.7, and 423.4 nm). The BITSE instrument is an externally occulted single stage coronagraph developed at NASA's Goddard Space Flight Center in collaboration with the Korea Astronomy and Space Science Institute (KASI). BITSE used a polarization camera that provided polarization and total brightness images of size 1024 x 1024 pixels. The Wallops Arc Second Pointing (WASP) system developed at NASA's Wallops Flight Facility (WFF) was used for Sun-pointing. The coronagraph and WASP were mounted on a gondola provided by WFF and launched from the Fort Sumner, New Mexico station of Columbia Scientific Balloon Facility (CSBF) on September 18, 2019. BITSE obtained 17,060 coronal images at a float altitude of about 128,000 feet (39 km) over a period of about 4 hrs. BITSE flight software was based on NASA's core Flight System, which was designed to help develop flight quality software. We used EVTM (Ethernet Via Telemetry) to download science data during operations; all images were stored onboard using flash storage. At the end of the mission, all data were recovered and analyzed. Preliminary analysis shows that BITSE imaged the solar minimum corona with the equatorial streamers on the east and west limbs. The narrow streamers observed by BITSE are in good agreement with the geometric properties obtained by SOHO coronagraphs in the overlapping physical domain. In spite of the small signal-to-noise ratio (about 14) we were able to obtain the temperature and flow speed of the western steamer region in the range 4 to 7 Rs as: For the equatorial streamer on the west limb, we obtained a temperature of 1.0 +/- 0.3 MK and a flow speed of about 260 km/s with a large uncertainty interval.


title: Tree Embeddings for Hop-Constrained Network Design

abstract: Network design problems aim to compute low-cost structures such as routes, trees and subgraphs. Often, it is natural and desirable to require that these structures have small hop length or hop diameter. Unfortunately, optimization problems with hop constraints are much harder and less well understood than their hop-unconstrained counterparts. A significant algorithmic barrier in this setting is the fact that hop-constrained distances in graphs are very far from being a metric. We show that, nonetheless, hop-constrained distances can be approximated by distributions over "partial tree metrics." We build this result into a powerful and versatile algorithmic tool which, similarly to classic probabilistic tree embeddings, reduces hop-constrained problems in general graphs to hop-unconstrained problems on trees. We then use this tool to give the first poly-logarithmic bicriteria approximations for the hop-constrained variants of many classic network design problems. These include Steiner forest, group Steiner tree, group Steiner forest, buy-at-bulk network design as well as online and oblivious versions of many of these problems.


title: Storage capacity and learning capability of quantum neural networks

abstract: We study the storage capacity of quantum neural networks (QNNs) described as completely positive trace preserving (CPTP) maps, which act on an $N$-dimensional Hilbert space. We demonstrate that QNNs can store up to $N$ linearly independent pure states and provide the structure of the corresponding maps. While the storage capacity of a classical Hopfield network scales linearly with the number of neurons, we show that QNNs can store an exponential number of linearly independent states. We estimate, employing the Gardner program, the relative volume of CPTP maps with $M$ stationary states. The volume decreases exponentially with $M$ and shrinks to zero for $M\geq N+1$. We generalize our results to QNNs storing mixed states as well as input-output relations for feed-forward QNNs. Our approach opens the path to relate storage properties of QNNs to the quantum properties of the input-output states. This paper is dedicated to the memory of Peter Wittek.


title: Variation of GIT and Variation of Lagrangian Skeletons II: Quasi-Symmetric Case

abstract: Consider $(\mathbb{C}^*)^k$ acting on $\mathbb{C}^N$ satisfying certain 'quasi-symmetric' condition which produces a class of toric Calabi-Yau GIT quotient stacks. Using subcategories of $Coh([\mathbb{C}^N / (\mathbb{C}^*)^k])$ generated by line bundles whose weights are inside certain zonotope called the 'magic window', Halpern-Leistner and Sam give a combinatorial construction of equivalences between derived categories of coherent sheaves for various GIT quotients. We apply the coherent-constructible correspondence for toric varieties to the magic windows and obtain a non-characteristic deformation of Lagrangian skeletons in $\mathbb{R}^{N-k}$ parameterized by $\mathbb{R}^k$, exhibiting derived equivalences between A-models of the various phases. Moreover, by translating the magic window zonotope in $\mathbb{R}^k$, we obtain a universal skeleton over $\mathbb{R}^k \times \mathbb{R}^k \setminus \mathcal{D}$ for some fattening of hyperplane arrangements $\mathcal{D}$, and we show that the the universal skeleton induces a local system of categories over $\mathbb{R}^k \times \mathbb{R}^k \setminus \mathcal{D}$. We also connect our results to the perverse schober structure identified by \v{S}penko and Van den Bergh.


title: Making the Universe at 20 MeV

abstract: We present a testable mechanism of low-scale baryogenesis and dark matter production in which neither baryon nor lepton number are violated. Charged $D$ mesons are produced out-of-equilibrium at tens of MeV temperatures. The $D$ mesons quickly undergo CP-violating decays to charged pions, which then decay into dark-sector leptons without violating lepton number. To transfer this lepton asymmetry to the baryon asymmetry, the dark leptons scatter on additional dark-sector states charged under lepton and baryon number. Amusingly, this transfer proceeds without electroweak sphalerons, which are no longer active at such low scales. We present two example models which can achieve this transfer while remaining consistent with current limits. The required amount of CP violation in charged $D$ meson decays, while currently allowed, will be probed by colliders. Additionally, the relevant decays of charged pions to dark-sector leptons have been constrained by the PIENU and PSI experiments and will be further explored in upcoming experiments.


title: A Data-Driven Reinforcement Learning Solution Framework for Optimal and Adaptive Personalization of a Hip Exoskeleton

abstract: Robotic exoskeletons are exciting technologies for augmenting human mobility. However, designing such a device for seamless integration with the human user and to assist human movement still is a major challenge. This paper aims at developing a novel data-driven solution framework based on reinforcement learning (RL), without first modeling the human-robot dynamics, to provide optimal and adaptive personalized torque assistance for reducing human efforts during walking. Our automatic personalization solution framework includes the assistive torque profile with two control timing parameters (peak and offset timings), the least square policy iteration (LSPI) for learning the parameter tuning policy, and a cost function based on transferred work ratio. The proposed controller was successfully validated on a healthy human subject to assist unilateral hip extension in walking. The results showed that the optimal and adaptive RL controller as a new approach was feasible for tuning assistive torque profile of the hip exoskeleton that coordinated with human actions and reduced activation level of hip extensor muscle in human.


title: Modified Macdonald polynomials and the multispecies zero range process: I

abstract: In this paper we prove a new combinatorial formula for the modified Macdonald polynomials $\widetilde{H}_\lambda(X;q,t)$, motivated by connections to the theory of interacting particle systems from statistical mechanics. The formula involves a new statistic called queue inversions on fillings of tableaux. This statistic is closely related to the multiline queues which were recently used to give a formula for the Macdonald polynomials $P_\lambda(X;q,t)$. In the case $q=1$ and $X=(1,1,\dots,1)$, that formula had also been shown to compute stationary probabilities for a particle system known as the multispecies ASEP on a ring, and it is natural to ask whether a similar connection exists between the modified Macdonald polynomials and a suitable statistical mechanics model. In a sequel to this work, we demonstrate such a connection, showing that the stationary probabilities of the multispecies totally asymmetric zero-range process (mTAZRP) on a ring can be computed using tableaux formulas with the queue inversion statistic. This connection extends to arbitrary $X=(x_1,\dots, x_n)$; the $x_i$ play the role of site-dependent jump rates for the mTAZRP.


title: I Know What You Meant: Learning Human Objectives by (Under)estimating Their Choice Set

abstract: Assistive robots have the potential to help people perform everyday tasks. However, these robots first need to learn what it is their user wants them to do. Teaching assistive robots is hard for inexperienced users, elderly users, and users living with physical disabilities, since often these individuals are unable to show the robot their desired behavior. We know that inclusive learners should give human teachers credit for what they cannot demonstrate. But today's robots do the opposite: they assume every user is capable of providing any demonstration. As a result, these robots learn to mimic the demonstrated behavior, even when that behavior is not what the human really meant! Here we propose a different approach to reward learning: robots that reason about the user's demonstrations in the context of similar or simpler alternatives. Unlike prior works -- which err towards overestimating the human's capabilities -- here we err towards underestimating what the human can input (i.e., their choice set). Our theoretical analysis proves that underestimating the human's choice set is risk-averse, with better worst-case performance than overestimating. We formalize three properties to generate similar and simpler alternatives. Across simulations and a user study, our resulting algorithm better extrapolates the human's objective. See the user study here: https://youtu.be/RgbH2YULVRo


title: The effect of the WGC condition on the maximal energy extracted from black holes

abstract: In this paper, we study the Penrose process and the maximum energy extracted from the collision of two particles near the Kerr-Newman black hole with WGC condition. We consider the collision process when two particles collide in the ergosphere region of a black hole, scattering two new particles; one of them falls into the black hole and the other escapes to infinity. Our calculations also show an increasing in the energy. We also checked how much energy could be received from a black hole for particles with different spins. Results of this paper will help us to identify and study the black holes of astrophysics and the particles with different spins. Also, it states how black holes work in WGC conditions.


title: An Argument for Strong Positivity of the Decoherence Functional

abstract: We give an argument for strong positivity of the decoherence functional as the correct, physical positivity condition in formulations of quantum theory based fundamentally on the path integral. We extend to infinite systems work by Boes and Navascues that shows that the set of strongly positive quantum systems is maximal amongst sets of systems that are closed under tensor product composition. We show further that the set of strongly positive quantum systems is the unique set that is maximal amongst sets that are closed under tensor product composition.


title: Multi-wavelength high energy gas-filled fiber Raman laser spanning from 1.53 um to 2.4 um

abstract: In this work, we present a high pulse energy multi-wavelength Raman laser spanning from 1.53 um up to 2.4 um by employing the cascaded rotational stimulated Raman scattering (SRS) effect in a 5-m hydrogen (H2) -filled nested anti-resonant fiber (NARF), pumped by a linearly polarized Er/Yb fiber laser with a peak power of ~13 kW and pulse duration of ~7 ns in the C-band. The developed Raman laser has distinct lines at 1683 nm, 1868 nm, 2100 nm, and 2400 nm, with pulse energies as high as 18.25 uJ, 14.4 uJ, 14.1 uJ, and 8.2 uJ, respectively. We demonstrate how the energy in the Raman lines can be controlled by tuning the H2 pressure from 1 bar to 20 bar


title: Bayes Optimal Informer Sets for Early-Stage Drug Discovery

abstract: An important experimental design problem in early-stage drug discovery is how to prioritize available compounds for testing when very little is known about the target protein. Informer based ranking (IBR) methods address the prioritization problem when the compounds have provided bioactivity data on other potentially relevant targets. An IBR method selects an informer set of compounds, and then prioritizes the remaining compounds on the basis of new bioactivity experiments performed with the informer set on the target. We formalize the problem as a two-stage decision problem and introduce the Bayes Optimal Informer SEt (BOISE) method for its solution. BOISE leverages a flexible model of the initial bioactivity data, a relevant loss function, and effective computational schemes to resolve the two-step design problem. We evaluate BOISE and compare it to other IBR strategies in two retrospective studies, one on protein-kinase inhibition and the other on anti-cancer drug sensitivity. In both empirical settings BOISE exhibits better predictive performance than available methods. It also behaves well with missing data, where methods that use matrix completion show worse predictive performance. We provide an R implementation of BOISE at https://github.com/wiscstatman/esdd/BOISE


title: Deep learning and hand-crafted features for virus image classification

abstract: In this work, we present an ensemble of descriptors for the classification of transmission electron microscopy images of viruses. We propose to combine handcrafted and deep learning approaches for virus image classification. The set of handcrafted is mainly based on Local Binary Pattern variants, for each descriptor a different Support Vector Machine is trained, then the set of classifiers is combined by sum rule. The deep learning approach is a densenet201 pretrained on ImageNet and then tuned in the virus dataset, the net is used as features extractor for feeding another Support Vector Machine, in particular the last average pooling layer is used as feature extractor. Finally, classifiers trained on handcrafted features and classifier trained on deep learning features are combined by sum rule. The proposed fusion strongly boosts the performance obtained by each stand-alone approach, obtaining state of the art performance.


title: SEIR-Campus: Modeling Infectious Diseases on University Campuses

abstract: We introduce a Python package for modeling and studying the spread of infectious diseases using an agent-based SEIR style epidemiological model with a focus on university campuses. This document explains the epidemiological model used in the package and gives examples highlighting the ways that the package can be used.


title: Hurricane Forecasting: A Novel Multimodal Machine Learning Framework

abstract: This paper describes a novel machine learning (ML) framework for tropical cyclone intensity and track forecasting, combining multiple ML techniques and utilizing diverse data sources. Our multimodal framework, called Hurricast, efficiently combines spatial-temporal data with statistical data by extracting features with deep-learning encoder-decoder architectures and predicting with gradient-boosted trees. We evaluate our models in the North Atlantic and Eastern Pacific basins on 2016-2019 for 24-hour lead time track and intensity forecasts and show they achieve comparable mean absolute error and skill to current operational forecast models while computing in seconds. Furthermore, the inclusion of Hurricast into an operational forecast consensus model could improve over the National Hurricane Center's official forecast, thus highlighting the complementary properties with existing approaches. In summary, our work demonstrates that utilizing machine learning techniques to combine different data sources can lead to new opportunities in tropical cyclone forecasting.