title: Stochastic solutions of Stefan problems

abstract: This work deals with the one-dimensional Stefan problem with a general time-dependent boundary condition at the fixed boundary. Stochastic solutions are obtained using discrete random walks, and the results are compared with analytic formulae when they exist, otherwise with numerical solutions from a finite difference method. The innovative part is to model the moving boundary with a random walk method. The results show statistical convergence for many random walkers when $\Delta x \rightarrow 0$. Stochastic methods are very competitive in large domains in higher dimensions and has the advantages of generality and ease of implementation. The stochastic method suffers from that longer execution times are required for increased accuracy. Since the code is easily adapted for parallel computing, it is possible to speed up the calculations. Regarding applications for Stefan problems, they have historically been used to model the dynamics of melting ice, and we give such an example here where the fixed boundary condition follows data from observed day temperatures at \"{O}rebro airport. Nowadays, there are a large range of examples of applications, such as climate models, the diffusion of lithium-ions in lithium-ion batteries and modelling steam chambers for petroleum extraction.


title: An Improved and Parallel Version of a Scalable Algorithm for Analyzing Time Series Data

abstract: Today, very large amounts of data are produced and stored in all branches of society including science. Mining these data meaningfully has become a considerable challenge and is of the broadest possible interest. The size, both in numbers of observations and dimensionality thereof, requires data mining algorithms to possess time complexities with both variables that are linear or nearly linear. One such algorithm, see Comput. Phys. Commun. 184, 2446-2453 (2013), arranges observations into a sequence called the progress index. The progress index steps through distinct regions of high sampling density sequentially. By means of suitable annotations, it allows a compact representation of the behavior of complex systems, which is encoded in the original data set. The only essential parameter is a notion of distance between observations. Here, we present the shared memory parallelization of the key step in constructing the progress index, which is the calculation of an approximation of the minimum spanning tree of the complete graph of observations. We demonstrate that excellent parallel efficiencies are obtained for up to 72 logical (CPU) cores. In addition, we introduce three conceptual advances to the algorithm that improve its controllability and the interpretability of the progress index itself.


title: Persona2vec: A Flexible Multi-role Representations Learning Framework for Graphs

abstract: Graph embedding techniques, which learn low-dimensional representations of a graph, are achieving state-of-the-art performance in many graph mining tasks. Most existing embedding algorithms assign a single vector to each node, implicitly assuming that a single representation is enough to capture all characteristics of the node. However, across many domains, it is common to observe pervasively overlapping community structure, where most nodes belong to multiple communities, playing different roles depending on the contexts. Here, we propose persona2vec, a graph embedding framework that efficiently learns multiple representations of nodes based on their structural contexts. Using link prediction-based evaluation, we show that our framework is significantly faster than the existing state-of-the-art model while achieving better performance.


title: CRISP: A Probabilistic Model for Individual-Level COVID-19 Infection Risk Estimation Based on Contact Data

abstract: We present CRISP (COVID-19 Risk Score Prediction), a probabilistic graphical model for COVID-19 infection spread through a population based on the SEIR model where we assume access to (1) mutual contacts between pairs of individuals across time across various channels (e.g., Bluetooth contact traces), as well as (2) test outcomes at given times for infection, exposure and immunity tests. Our micro-level model keeps track of the infection state for each individual at every point in time, ranging from susceptible, exposed, infectious to recovered. We develop both a Monte Carlo EM as well as a message passing algorithm to infer contact-channel specific infection transmission probabilities. Our Monte Carlo algorithm uses Gibbs sampling to draw samples of the latent infection status of each individual over the entire time period of analysis, given the latent infection status of all contacts and test outcome data. Experimental results with simulated data demonstrate our CRISP model can be parametrized by the reproduction factor $R_0$ and exhibits population-level infectiousness and recovery time series similar to those of the classical SEIR model. However, due to the individual contact data, this model allows fine grained control and inference for a wide range of COVID-19 mitigation and suppression policy measures. Moreover, the block-Gibbs sampling algorithm is able to support efficient testing in a test-trace-isolate approach to contain COVID-19 infection spread. To the best of our knowledge, this is the first model with efficient inference for COVID-19 infection spread based on individual-level contact data; most epidemic models are macro-level models that reason over entire populations. The implementation of CRISP is available in Python and C++ at https://github.com/zalandoresearch/CRISP.


title: How Smart is the Grid?

abstract: Ancient Romans called 'urbs' the set of buildings and infrastructures, and 'civitas' the Roman citizens. Today instead, while the society is surfing the digital tsunami, 'urbs' and 'civitas' tend to become much closer, almost merging, that we might attempt to condensate these into a single concept: 'smart grid'. Internet of things, artificial intelligence, blockchain, quantum cryptography is only a few of the paradigms that are likely to contribute to determining the final portrait of the future smart grid. However, to understand the effective sustainability of complex grids, specific tools are required. To this end, in this article, a systematic review of the emerging paradigms is presented, identifying intersectoral synergies and limitations with respect to the `smart grid' concept. Further, a taxonomic framework for assessing the level of sustainability of the grid is proposed. Finally, from the scenario portrayed, a set of issues involving engineering, regulation, security, and social frameworks have been derived in a theoretical fashion. The findings are likely to suggest the urgent need for multidisciplinary cooperation to wisely address engineering and ontological challenges gravitating around the smart grid concept.


title: A Machine Learning System for Retaining Patients in HIV Care

abstract: Retaining persons living with HIV (PLWH) in medical care is paramount to preventing new transmissions of the virus and allowing PLWH to live normal and healthy lifespans. Maintaining regular appointments with an HIV provider and taking medication daily for a lifetime is exceedingly difficult. 51% of PLWH are non-adherent with their medications and eventually drop out of medical care. Current methods of re-linking individuals to care are reactive (after a patient has dropped-out) and hence not very effective. We describe our system to predict who is most at risk to drop-out-of-care for use by the University of Chicago HIV clinic and the Chicago Department of Public Health. Models were selected based on their predictive performance under resource constraints, stability over time, as well as fairness. Our system is applicable as a point-of-care system in a clinical setting as well as a batch prediction system to support regular interventions at the city level. Our model performs 3x better than the baseline for the clinical model and 2.3x better than baseline for the city-wide model. The code has been released on github and we hope this methodology, particularly our focus on fairness, will be adopted by other clinics and public health agencies in order to curb the HIV epidemic.


title: Gradient Boosting Application in Forecasting of Performance Indicators Values for Measuring the Efficiency of Promotions in FMCG Retail

abstract: In the paper, a problem of forecasting promotion efficiency is raised. The authors propose a new approach, using the gradient boosting method for this task. Six performance indicators are introduced to capture the promotion effect. For each of them, within predefined groups of products, a model was trained. A description of using these models for forecasting and optimising promotion efficiency is provided. Data preparation and hyperparameters tuning processes are also described. The experiments were performed for three groups of products from a large grocery company.


title: Transforming unstructured voice and text data into insight for paramedic emergency service using recurrent and convolutional neural networks

abstract: Paramedics often have to make lifesaving decisions within a limited time in an ambulance. They sometimes ask the doctor for additional medical instructions, during which valuable time passes for the patient. This study aims to automatically fuse voice and text data to provide tailored situational awareness information to paramedics. To train and test speech recognition models, we built a bidirectional deep recurrent neural network (long short-term memory (LSTM)). Then we used convolutional neural networks on top of custom-trained word vectors for sentence-level classification tasks. Each sentence is automatically categorized into four classes, including patient status, medical history, treatment plan, and medication reminder. Subsequently, incident reports were automatically generated to extract keywords and assist paramedics and physicians in making decisions. The proposed system found that it could provide timely medication notifications based on unstructured voice and text data, which was not possible in paramedic emergencies at present. In addition, the automatic incident report generation provided by the proposed system improves the routine but error-prone tasks of paramedics and doctors, helping them focus on patient care.


title: Adherence to Personal Health Devices: A Case Study in Diabetes Management

abstract: Personal health devices can enable continuous monitoring of health parameters. However, the benefit of these devices is often directly related to the frequency of use. Therefore, adherence to personal health devices is critical. This paper takes a data mining approach to study continuous glucose monitor use in diabetes management. We evaluate two independent datasets from a total of 44 subjects for 60 - 270 days. Our results show that: 1) missed target goals (i.e. suboptimal outcomes) is a factor that is associated with wearing behavior of personal health devices, and 2) longer duration of non-adherence, identified through missing data or data gaps, is significantly associated with poorer outcomes. More specifically, we found that up to 33% of data gaps occurred when users were in abnormal blood glucose categories. The longest data gaps occurred in the most severe (i.e. very low / very high) glucose categories. Additionally, subjects with poorly-controlled diabetes had longer average data gap duration than subjects with well-controlled diabetes. This work contributes to the literature on the design of context-aware systems that can leverage data-driven approaches to understand factors that influence non-wearing behavior. The results can also support targeted interventions to improve health outcomes.


title: AI Research Considerations for Human Existential Safety (ARCHES)

abstract: Framed in positive terms, this report examines how technical AI research might be steered in a manner that is more attentive to humanity's long-term prospects for survival as a species. In negative terms, we ask what existential risks humanity might face from AI development in the next century, and by what principles contemporary technical research might be directed to address those risks. A key property of hypothetical AI technologies is introduced, called \emph{prepotence}, which is useful for delineating a variety of potential existential risks from artificial intelligence, even as AI paradigms might shift. A set of \auxref{dirtot} contemporary research \directions are then examined for their potential benefit to existential safety. Each research direction is explained with a scenario-driven motivation, and examples of existing work from which to build. The research directions present their own risks and benefits to society that could occur at various scales of impact, and in particular are not guaranteed to benefit existential safety if major developments in them are deployed without adequate forethought and oversight. As such, each direction is accompanied by a consideration of potentially negative side effects.


title: A refined dynamic finite-strain shell theory for incompressible hyperelastic materials: equations and two-dimensional shell virtual work principle

abstract: Based on previous work for the static problem, in this paper we first derive one form of dynamic finite-strain shell equations for incompressible hyperelastic materials that involve three shell constitutive relations. In order to single out the bending effect as well as to reduce the number of shell constitutive relations, a further refinement is performed, which leads to a refined dynamic finite-strain shell theory with only two shell constitutive relations (deducible from the given three-dimensional (3D) strain energy function) and some new insights are also deduced. By using the weak formulation of the shell equations and the variation of the 3D Lagrange functional, boundary conditions and the two-dimensional (2D) shell virtual work principle are derived. As a benchmark problem, we consider the extension and inflation of an arterial segment. The good agreement between the asymptotic solution based on the shell equations and that from the 3D exact one gives verification of the former. The refined shell theory is also applied to study the plane-strain vibrations of a pressurized artery, and the effects of the axial pre-stretch, pressure and fibre angle on the vibration frequencies are investigated in detail.


title: Machine Learning Systems for Intelligent Services in the IoT: A Survey

abstract: Machine learning (ML) technologies are emerging in the Internet of Things (IoT) to provision intelligent services. This survey moves beyond existing ML algorithms and cloud-driven design to investigate the less-explored systems, scaling and socio-technical aspects for consolidating ML and IoT. It covers the latest developments (up to 2020) on scaling and distributing ML across cloud, edge, and IoT devices. With a multi-layered framework to classify and illuminate system design choices, this survey exposes fundamental concerns of developing and deploying ML systems in the rising cloud-edge-device continuum in terms of functionality, stakeholder alignment and trustworthiness.


title: Network visualizations with Pyvis and VisJS

abstract: Pyvis is a Python module that enables visualizing and interactively manipulating network graphs in the Jupyter notebook, or as a standalone web application. Pyvis is built on top of the powerful and mature VisJS JavaScript library, which allows for fast and responsive interactions while also abstracting away the low-level JavaScript and HTML. This means that elements of the rendered graph visualization, such as node/edge attributes can be specified within Python and shipped to the JavaScript layer for VisJS to render. This declarative approach makes it easy to quickly explore graph visualizations and investigate data relationships. In addition, Pyvis is highly customizable so that colors, sizes, and hover tooltips can be assigned to the rendered graph. The network graph layout is controlled by a front-end physics engine that is configurable from a Python interface, allowing for the detailed placement of the graph elements. In this paper, we outline use cases for Pyvis with specific examples to highlight key features for any analysis workflow. A brief overview of Pyvis' implementation describes how the Python front-end binding uses simple Pyvis calls.


title: General method to calculate the elastic deformation and X-ray diffraction properties of bent crystal wafers

abstract: Curved single crystals are widely employed in spectrometer designs in the hard X-ray regime. Due to their large solid angle coverage and focusing properties, toroidally bent crystals are extremely useful in applications where the output of photons is low. Spherically bent crystals, a subgroup of toroidally bent crystals, particularly have found their way in many instruments at synchtrotrons and free electron laser lightsource end-stations but also in the re-emerging field of high-resolution laboratory-scale X-ray spectroscopy. A solid theoretical understanding of the diffraction properties of such crystals is essential when aiming for optimal spectrometer performance. In this work, we present a general method to calculate the internal stress and strain fields of toroidally bent crystals and how to apply it to predict their diffraction properties. Solutions are derived and discussed for circular and rectangular spherically bent wafers due to their prevalence in contemporary instrumentation.


title: Hedging in games: Faster convergence of external and swap regrets

abstract: We consider the setting where players run the Hedge algorithm or its optimistic variant to play an $n$-action game repeatedly for $T$ rounds. 1) For two-player games, we show that the regret of optimistic Hedge decays at $\tilde{O}( 1/T ^{5/6} )$, improving the previous bound $O(1/T^{3/4})$ by Syrgkanis, Agarwal, Luo and Schapire (NIPS'15) 2) In contrast, we show that the convergence rate of vanilla Hedge is no better than $\tilde{\Omega}(1/ \sqrt{T})$, addressing an open question posted in Syrgkanis, Agarwal, Luo and Schapire (NIPS'15). For general m-player games, we show that the swap regret of each player decays at rate $\tilde{O}(m^{1/2} (n/T)^{3/4})$ when they combine optimistic Hedge with the classical external-to-internal reduction of Blum and Mansour (JMLR'07). The algorithm can also be modified to achieve the same rate against itself and a rate of $\tilde{O}(\sqrt{n/T})$ against adversaries. Via standard connections, our upper bounds also imply faster convergence to coarse correlated equilibria in two-player games and to correlated equilibria in multiplayer games.


title: Morphology and numerical characteristics of epidemic curves for SARS-Cov-II using Moyal distribution

abstract: In this paper, it is shown that the Moyal distribution is an excelent tool to study the SARS-Cov-II (Covid-19) epidemiological associated curves and its propagation. The Moyal parameters give all the information to describe the form and the impact of the illness outbreak in the different affected countries and its global impact. We checked that the Moyal distribution can accurately fit the daily report of {\it{new confirmed cases of infected people}} (NCC) per country, in that places where the contagion is reaching their final phase, describing the beginning, the most intense phase and the descend of the contagion, simultaneously . In order to achieve the purpose of this work, it is important to work with a complete and well compilated set of the data to be used to fit the curves. Data from European countries like France, Spain, Italy Belgium, Sweden, United Kingdom, Denmark and others like USA and China, have been used. Also, the correlation between the parameters of the Moyal distribution fitting and the general public health measures imposed in each country, have been discussed. A relation between those policies and the features of the Moyal distribution, in terms of their parameters and critical points, is shown; from that, it can be seen that the knowledge of the time evolution of the epidemiological curve, their critical points, superposition properties and rates of the rising and the ending, could help to find a way to estimate the efficiency of social distancing measures, imposed in each country, and anticipate the different phases of the pandemic.


title: White-dwarf asteroseismology with the Kepler space telescope

abstract: In the course of their evolution, white-dwarf stars go through at least one phase of variability in which the global pulsations they undergo allow astronomers to peer into their interiors, this way making possible to shed light on their deep inner structure and evolutionary stage by means of asteroseismology. The study of pulsating white dwarfs has witnessed substantial progress in the last decade, and this has been so largely thanks to the arrival of continuous observations of unprecedented quality from space, like those of the CoRoT, Kepler, and TESS missions. This, along with the advent of new detailed thoretical models and the development of improved asteroseismological techniques, has helped to unravel the internal chemical structure of many pulsating white dwarfs, and, at the same time, has opened new questions that challenge theoreticians. In particular, uninterrupted monitoring of white-dwarf stars for months has allowed discovering phenomena impossible to detect with ground-based observations, despite admirable previous efforts like the Whole Earth Telescope (WET). Here, we start by reviewing the essential properties of white-dwarf and pre-white dwarf stars and their pulsations, and then, we go through the different families of pulsating objects known to date. Finally, we review the most outstanding findings about pulsating white dwarfs and pre-white dwarfs made possible with the unprecedented-quality observations of the Kepler space telescope, although we envisage that future analyzes of space data from this mission that still await to be examined may reveal new secrets of these extremely interesting variable stars.


title: Consistent analytical solution of the time-dependent Schr\"odinger equation for nanoscale circuits with laser-assisted quantum tunneling

abstract: It is now common practice to solve the Schr\"odinger equation to estimate the tunneling current between two metal electrodes at specified potentials, or the transmission through a potential barrier by assuming an incident, reflected, and transmitted wave. However, we suggest that these methods may not be appropriate for nanoscale circuits. The electron man-free path may be as long as 68.2 nm in metallic elements so we consider the possibility that quantum effects may occur throughout a nanoscale circuit, including the connections. Analytical methods are presented for modeling the coherent transfer of the wavefunction through a closed circuit.


title: Stable characters from permutation patterns

abstract: For a fixed permutation $\sigma \in S_k$, let $N_{\sigma}$ denote the function which counts occurrences of $\sigma$ as a pattern in permutations from $S_n$. We study the expected value (and $d$-th moments) of $N_{\sigma}$ on conjugacy classes of $S_n$ and prove that the irreducible character support of these class functions stabilizes as $n$ grows. This says that there is a single polynomial in the variables $n, m_1, \ldots, m_{dk}$ which computes these moments on any conjugacy class (of cycle type $1^{m_1}2^{m_2}\cdots$) of any symmetric group. This result generalizes results of Hultman and of Gill, who proved the cases $(d,k)=(1,2)$ and $(1,3)$ using ad hoc methods. Our proof is, to our knowledge, the first application of partition algebras to the study of permutation patterns.


title: Duality and symmetry of complexity over complete intersections via exterior homology

abstract: We study homological properties of a locally complete intersection ring by importing facts from homological algebra over exterior algebras. One application is showing that the thick subcategories of the bounded derived category of a locally complete intersection ring are self-dual under Grothendieck duality. This was proved by Stevenson when the ring is a quotient of a regular ring modulo a regular sequence; we offer two independent proofs in the more general setting. Second, we use these techniques to supply new proofs that complete intersections possess symmetry of complexity.


title: Thinking Taxonomically about Fake Accounts: Classification, False Dichotomies, and the Need for Nuance

abstract: It is often said that war creates a fog in which it becomes difficult to discern friend from foe on the battlefield. In the ongoing war on fake accounts, conscious development of taxonomies of the phenomenon has yet to occur, resulting in much confusion on the digital battlefield about what exactly a fake account is. This paper intends to address this problem, not by proposing a taxonomy of fake accounts, but by proposing a systematic way to think taxonomically about the phenomenon. Specifically, we examine fake accounts through both a combined philosophical and computer science-based perspective. Through these lenses, we deconstruct narrow binary thinking about fake accounts, both in the form of general false dichotomies and specifically in relation to the Facebook's conceptual framework "Coordinated Inauthentic Behavior" (CIB). We then address the false dichotomies by constructing a more complex way of thinking taxonomically about fake accounts.


title: A Notion of Individual Fairness for Clustering

abstract: A common distinction in fair machine learning, in particular in fair classification, is between group fairness and individual fairness. In the context of clustering, group fairness has been studied extensively in recent years; however, individual fairness for clustering has hardly been explored. In this paper, we propose a natural notion of individual fairness for clustering. Our notion asks that every data point, on average, is closer to the points in its own cluster than to the points in any other cluster. We study several questions related to our proposed notion of individual fairness. On the negative side, we show that deciding whether a given data set allows for such an individually fair clustering in general is NP-hard. On the positive side, for the special case of a data set lying on the real line, we propose an efficient dynamic programming approach to find an individually fair clustering. For general data sets, we investigate heuristics aimed at minimizing the number of individual fairness violations and compare them to standard clustering approaches on real data sets.


title: The weight distributions of linear sets in PG(1,q^5)

abstract: In this paper, we study the weight distributions of $\mathbb{F}_q$-linear sets in $\mathrm{PG}(1,q^5)$. Our main theorem proves that a linear set $S$ of rank $5$, which is not scattered has the following weight distribution for its points with weight larger than 1: (i) one point of weight $4$ or $5$, (ii) one point of weight $3$ and $0$, $q$, $q^2$ points of weight two, (iii) $s$ points of weight $2$ where $s\in [q-2\sqrt{q}+1,q+2\sqrt{q}+1]\cup\{2q,2q+1,2q+2,3q,3q+1,q^2+1\}$. In particular, there are no $2$-clubs in $\mathrm{PG}(1,q^5)$.


title: A Novel Navigation System for an Autonomous Mobile Robot in an Uncertain Environment

abstract: In this paper, we developed a new navigation system, which detects obstacles in a sliding window with an adaptive threshold clustering algorithm, classifies the detected obstacles with a decision tree, heuristically predicts potential collision and finds optimal path with a simplified Mophin algorithm. This system has the merits of optimal free-collision path, small memory size and less computing complexity, compared with the state of the arts in robot navigation. The experiments on simulation and a robot for eight scenarios demonstrate that the robot can effectively and efficiently avoid potential collisions with any static or dynamic obstacles in its surrounding environment.


title: Fast neutrino cooling of nuclear pasta in neutron stars: molecular dynamics simulations

abstract: The direct URCA process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crust of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct URCA reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct URCA reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified URCA process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars


title: InFocus: A spatial coding technique to mitigate misfocus in near field LoS beamforming

abstract: Phased arrays, commonly used in IEEE 802.11ad and 5G radios, are capable of focusing radio frequency signals in a specific direction or a spatial region. Beamforming achieves such directional or spatial concentration of signals and enables phased array-based radios to achieve high data rates. Designing beams for millimeter wave and terahertz communication using massive phased arrays, however, is challenging due to hardware constraints and the wide bandwidth in these systems. For example, beams which are optimal at the center frequency may perform poor in wideband communication systems where the radio frequencies differ substantially from the center frequency. The poor performance in such systems is due to differences in the optimal beamformers corresponding to distinct radio frequencies within the wide bandwidth. Such a mismatch leads to a misfocus effect in near field systems and the beam squint effect in far field systems. In this paper, we investigate the misfocus effect and propose InFocus, a low complexity technique to construct beams that are well suited for massive wideband phased arrays. The beams are constructed using a carefully designed frequency modulated waveform in the spatial dimension. For the special case of beamforming along the boresight of an array, this waveform is analogous to the frequency modulated continuous wave (FMCW) chirp signal in radar. InFocus mitigates beam misfocus and beam squint when applied to near field and far field systems. Simulation results indicate that InFocus enables massive wideband phased array-based radios to achieve higher data rates than comparable beamforming solutions.


title: The Kuranishi map for vector bundles on certain products of curves

abstract: We describe deformations of vector bundles on surfaces that are a product of two smooth projective curves. We explicitly describe the Kuranishi map around unstable vector bundles and compare the homologies of the Kuranishi spaces of stable and unstable deformations.


title: The fractional derivative of the Dirac delta function and new results on the inverse Laplace transform of irrational functions

abstract: Motivated from studies on anomalous diffusion, we show that the memory function $M(t)$ of complex materials, that their creep compliance follows a power law, $J(t)\sim t^q$ with $q\in \mathbb{R}^+$, is the fractional derivative of the Dirac delta function, $\frac{\mathrm{d}^q\delta(t-0)}{\mathrm{d}t^q}$ with $q\in \mathbb{R}^+$. This leads to the finding that the inverse Laplace transform of $s^q$ for any $q\in \mathbb{R}^+$ is the fractional derivative of the Dirac delta function, $\frac{\mathrm{d}^q\delta(t-0)}{\mathrm{d}t^q}$. This result, in association with the convolution theorem, makes possible the calculation of the inverse Laplace transform of $\frac{s^q}{s^{\alpha}\mp\lambda}$ where $\alpha<q\in\mathbb{R}^+$ which is the fractional derivative of order $q$ of the Rabotnov function $\varepsilon_{\alpha-1}(\pm\lambda, t)=t^{\alpha-1}E_{\alpha, \alpha}(\pm\lambda t^{\alpha})$. The fractional derivative of order $q\in \mathbb{R}^+$ of the Rabotnov function, $\varepsilon_{\alpha-1}(\pm\lambda, t)$ produces singularities which are extracted with a finite number of fractional derivatives of the Dirac delta function depending on the strength of $q$ in association with the recurrence formula of the two-parameter Mittag-Leffler function.


title: Summarising Big Data: Common GitHub Dataset for Software Engineering Challenges

abstract: In open-source software development environments; textual, numerical and relationship-based data generated are of interest to researchers. Various data sets are available for this data, which is frequently used in areas such as software engineering and natural language processing. However, since these data sets contain all the data in the environment, the problem arises in the terabytes of data processing. For this reason, almost all of the studies using GitHub data use filtered data according to certain criteria. In this context, using a different data set in each study makes a comparison of the accuracy of the studies quite difficult. In order to solve this problem, a common dataset was created and shared with the researchers, which would allow us to work on many software engineering problems.


title: Heterogeneous Effects of Job Displacement on Earnings

abstract: This paper considers how the effect of job displacement varies across different individuals. In particular, our interest centers on features of the distribution of the individual-level effect of job displacement. Identifying features of this distribution is particularly challenging -- e.g., even if we could randomly assign workers to be displaced or not, many of the parameters that we consider would not be point identified. We exploit our access to panel data, and our approach relies on comparing outcomes of displaced workers to outcomes the same workers would have experienced if they had not been displaced and if they maintained the same rank in the distribution of earnings as they had before they were displaced. Using data from the Displaced Workers Survey, we find that displaced workers earn about $157 per week less, on average, than they would have earned if they had not been displaced. We also find that there is substantial heterogeneity. We estimate that 42% of workers have higher earnings than they would have had if they had not been displaced and that a large fraction of workers have experienced substantially more negative effects than the average effect of displacement. Finally, we also document major differences in the distribution of the effect of job displacement across education levels, sex, age, and counterfactual earnings levels. Throughout the paper, we rely heavily on quantile regression. First, we use quantile regression as a flexible (yet feasible) first step estimator of conditional distributions and quantile functions that our main results build on. We also use quantile regression to study how covariates affect the distribution of the individual-level effect of job displacement.


title: Scalability in Computing and Robotics

abstract: Efficient engineered systems require scalability. A scalable system has increasing performance with increasing system size. In an ideal case, the increase in performance (e.g., speedup) corresponds to the number of units that are added to the system. However, if multiple units work on the same task, then coordination among these units is required. This coordination can introduce overheads with an impact on system performance. The coordination costs can lead to sublinear improvement or even diminishing performance with increasing system size. However, there are also systems that implement efficient coordination and exploit collaboration of units to attain superlinear improvement. Modeling the scalability dynamics is key to understanding efficient systems. Known laws of scalability, such as Amdahl's law, Gustafson's law, and Gunther's Universal Scalability Law, are minimalistic phenomenological models that explain a rich variety of system behaviors through concise equations. While useful to gain general insights, the phenomenological nature of these models may limit the understanding of the underlying dynamics, as they are detached from first principles that could explain coordination overheads among units. Through a decentralized system approach, we propose a general model based on generic interactions between units that is able to describe, as specific cases, any general pattern of scalability included by previously reported laws. The proposed general model of scalability is built on first principles, or at least on a microscopic description of interaction between units, and therefore has the potential to contribute to a better understanding of system behavior and scalability. We show that this model can be applied to a diverse set of systems, such as parallel supercomputers, robot swarms, or wireless sensor networks, creating a unified view on interdisciplinary design for scalability.


title: Degenerate Competing Three-Particle Systems

abstract: We study systems of three interacting particles, in which drifts and variances are assigned by rank. These systems are "degenerate": the variances corresponding to one or two ranks can vanish, so the corresponding ranked motions become ballistic rather than diffusive. Depending on which ranks are allowed to "go ballistic" the systems exhibit markedly different behavior which we study in some detail. Also studied are stability properties for the resulting planar process of gaps between successive ranks.


title: An extension of a cubic 2-connected plane graph G to a hamiltonian plane graph contained in G^{2}

abstract: Let $G$ be a simple cubic 2-connected plane graph. For every $2$-factor $X$ of $G$ having $n$-components there exists a simple hamiltonian plane graph $J \subset G^{2}$ such that $|E(J)|= |E(G)| + 2n -2$ and $\Delta(J) \leqslant 5$.


title: Multi-Fidelity High-Order Gaussian Processes for Physical Simulation

abstract: The key task of physical simulation is to solve partial differential equations (PDEs) on discretized domains, which is known to be costly. In particular, high-fidelity solutions are much more expensive than low-fidelity ones. To reduce the cost, we consider novel Gaussian process (GP) models that leverage simulation examples of different fidelities to predict high-dimensional PDE solution outputs. Existing GP methods are either not scalable to high-dimensional outputs or lack effective strategies to integrate multi-fidelity examples. To address these issues, we propose Multi-Fidelity High-Order Gaussian Process (MFHoGP) that can capture complex correlations both between the outputs and between the fidelities to enhance solution estimation, and scale to large numbers of outputs. Based on a novel nonlinear coregionalization model, MFHoGP propagates bases throughout fidelities to fuse information, and places a deep matrix GP prior over the basis weights to capture the (nonlinear) relationships across the fidelities. To improve inference efficiency and quality, we use bases decomposition to largely reduce the model parameters, and layer-wise matrix Gaussian posteriors to capture the posterior dependency and to simplify the computation. Our stochastic variational learning algorithm successfully handles millions of outputs without extra sparse approximations. We show the advantages of our method in several typical applications.


title: Pixel-Wise Motion Deblurring of Thermal Videos

abstract: Uncooled microbolometers can enable robots to see in the absence of visible illumination by imaging the "heat" radiated from the scene. Despite this ability to see in the dark, these sensors suffer from significant motion blur. This has limited their application on robotic systems. As described in this paper, this motion blur arises due to the thermal inertia of each pixel. This has meant that traditional motion deblurring techniques, which rely on identifying an appropriate spatial blur kernel to perform spatial deconvolution, are unable to reliably perform motion deblurring on thermal camera images. To address this problem, this paper formulates reversing the effect of thermal inertia at a single pixel as a Least Absolute Shrinkage and Selection Operator (LASSO) problem which we can solve rapidly using a quadratic programming solver. By leveraging sparsity and a high frame rate, this pixel-wise LASSO formulation is able to recover motion deblurred frames of thermal videos without using any spatial information. To compare its quality against state-of-the-art visible camera based deblurring methods, this paper evaluated the performance of a family of pre-trained object detectors on a set of images restored by different deblurring algorithms. All evaluated object detectors performed systematically better on images restored by the proposed algorithm rather than any other tested, state-of-the-art methods.


title: The ALPINE-ALMA [CII] Survey: The nature, luminosity function and star formation history of dusty galaxies up to z~6

abstract: We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 14 per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [CII] emitters at z=5). The six HST and near-IR dark galaxies with mid-IR counterpart contribute for about 17 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3.


title: Architectural Blueprints: The 4+1 View Model of Software Architecture

abstract: This article presents a model for describing the architecture of software-intensive systems, based on the use of multiple, concurrent views. This use of multiple views allows to address separately the concerns of the various stakeholders of the architecture: end-user, developers, systems engineers, project managers, etc., and to handle separately the functional and non functional requirements. Each of the five views is described, together with a notation to capture it. The views are designed using an architecture-centered, scenario-driven, iterative development process.


title: Physics Informed Deep Kernel Learning

abstract: Deep kernel learning is a promising combination of deep neural networks and nonparametric function learning. However, as a data driven approach, the performance of deep kernel learning can still be restricted by scarce or insufficient data, especially in extrapolation tasks. To address these limitations, we propose Physics Informed Deep Kernel Learning (PI-DKL) that exploits physics knowledge represented by differential equations with latent sources. Specifically, we use the posterior function sample of the Gaussian process as the surrogate for the solution of the differential equation, and construct a generative component to integrate the equation in a principled Bayesian hybrid framework. For efficient and effective inference, we marginalize out the latent variables in the joint probability and derive a collapsed model evidence lower bound (ELBO), based on which we develop a stochastic model estimation algorithm. Our ELBO can be viewed as a nice, interpretable posterior regularization objective. On synthetic datasets and real-world applications, we show the advantage of our approach in both prediction accuracy and uncertainty quantification.


title: Retakh's Motzkin paths and some combinatorial comments

abstract: Dyck paths where peaks are only allowed on level 1 and on even-indexed levels, were introduced by Retakh and analysed by Zeilberger, with assistance from Ekhad. We add some combinatorial comments to the enumeration, which involves Motzkin numbers, in particular, about the average height of such objects.


title: Estimates of dissipation of wave energy by sea ice for a field experiment in the Southern Ocean, using model/data inversion

abstract: A model-data inversion is applied to a very large observational dataset collected in the Southern Ocean north of the Ross Sea during late autumn to early winter, producing estimates of the frequency-dependent rate of dissipation by sea ice. The modeling platform is WAVEWATCH III(R) which accounts for non-stationarity, advection, wave generation, and other relevant processes. The resulting 9477 dissipation profiles are co-located with other variables such as ice thickness to quantify correlations which might be exploited in later studies to improve predictions. Mean dissipation profiles from the inversion are fitted to simple binomials. Variability about the mean profile is not small, but the binomials show remarkable qualitative similarity to prior observation-based estimates of dissipation, and the power dependence is consistent with at least three theoretical models, one of which assumes that dissipation is dominated by turbulence generated by shear at the ice-water interface.


title: Radiative Decay of Dark Exciton Related Emission in a Sandwiched Monolayer WSe2 Revealed by Room Temperature Micro and Nano Photoluminescence

abstract: TMDCs have attracted a lot of attention in recent years due to their unique indirect to direct band gap transition from bulk to monolayer thickness. Strong confinement in the out-of-plane direction enhances the Coulomb potential between the charged particles (e-h pairs) and thus increases the exciton binding energy dramatically. The lattice inversion asymmetry in a monolayer creates two non-equivalent (but degenerate in energy) band edges protected by time reversal polarisation via pseudo-spin. However, the presence of strong spin-orbit coupling in the valence band and weak spin-splitting in the conduction band results in the lowest lying exciton in WX2 (X = S, Se) being spin forbidden and optically dark. Because of their long life times, dark excitons (XD) are highly attractive for quantum optics and optoelectronic applications. To date studying XD emission is limited to cryogenic temperature or required very complex experimental configurations to observe them at room temperature (RT). Here, we demonstrate a novel approach of radiative decay of XD related emissions in 1L-WSe2 studied by micro and nano PL at RT. 1L-WSe2 flakes were sandwiched by noble metal (Au or Ag) substrates and PDMS nano-patches providing a strong local out-of-plane dipole moment with respect to the 2D plane. This strong dipole moment not only enhances the XD in WSe2, it also produces bound excitons due to extrinsic charge defects visible at RT. The spatial distributions of these XD related emissions were studied by TEPL with a spatial resolution < 10 nm confirming the confinement of these excitons within the PDMS nano-patches. Finally, by removing the nano-patches from the top of the flakes we are able to recover the bright excitons in the 1L-WSe2. Our approach paves the way for deep understanding and to harness excitonic properties in low dimensional semiconductors, thus offering a platform towards quantum optics.


title: The Structure of U.S. College Networks on Facebook

abstract: Anecdotally, social connections made in university have life-long impact. Yet knowledge of social networks formed in college remains episodic, due in large part to the difficulty and expense involved in collecting a suitable dataset for comprehensive analysis. To advance and systematize insight into college social networks, we describe a dataset of the largest online social network platform used by college students in the United States. We combine de-identified and aggregated Facebook data with College Scorecard data, campus-level information provided by U.S. Department of Education, to produce a dataset covering the 2008-2015 entry year cohorts for 1,159 U.S. colleges and universities, spanning 7.6 million students. To perform the difficult task of comparing these networks of different sizes we develop a new methodology. We compute features over sampled ego-graphs, train binary classifiers for every pair of graphs, and operationalize distance between graphs as predictive accuracy. Social networks of different year cohorts at the same school are structurally more similar to one another than to cohorts at other schools. Networks from similar schools have similar structures, with the public/private and graduation rate dimensions being the most distinguishable. We also relate school types to specific outcomes. For example, students at private schools have larger networks that are more clustered and with higher homophily by year. Our findings may help illuminate the role that colleges play in shaping social networks which partly persist throughout people's lives.


title: A Framework for Neural Network Pruning Using Gibbs Distributions

abstract: Modern deep neural networks are often too large to use in many practical scenarios. Neural network pruning is an important technique for reducing the size of such models and accelerating inference. Gibbs pruning is a novel framework for expressing and designing neural network pruning methods. Combining approaches from statistical physics and stochastic regularization methods, it can train and prune a network simultaneously in such a way that the learned weights and pruning mask are well-adapted for each other. It can be used for structured or unstructured pruning and we propose a number of specific methods for each. We compare our proposed methods to a number of contemporary neural network pruning methods and find that Gibbs pruning outperforms them. In particular, we achieve a new state-of-the-art result for pruning ResNet-56 with the CIFAR-10 dataset.


title: The soft-gluon limit and the infrared enhancement of the quark-gluon vertex

abstract: The Schwinger-Dyson quark equation (SDE) combined with results from lattice simulation for the propagators are used to obtain information on the quark-gluon vertex, taking into account the recent full QCD lattice results for the soft-gluon limit. Its inclusion leads to a clear enhancement of the infrared quark-gluon vertex. We also find that the relative contribution of the quark-ghost kernel to the quark-gluon vertex in the infrared region does not follow the rules from the perturbative analysis of the ultraviolet region. This shows that for QCD the intuition based on perturbation theory does not apply to the full momentum range. The framework developed in the current work provides analytical expressions for all the longitudinal components of vertex taken into account.


title: Modulation Format Dependent, Closed-Form Formula for Estimating Nonlinear Interference in S+C+L Band Systems

abstract: A closed-form formula for the nonlinear interference estimation of arbitrary modulation formats in ultra-wideband transmission systems is presented. Enabled by the proposed approach, the formula is applied to the entire S+C+L band (20 THz) and validated by numerical simulations with experimentally measured fibre data.


title: Making Convolutions Resilient via Algorithm-Based Error Detection Techniques

abstract: The ability of Convolutional Neural Networks (CNNs) to accurately process real-time telemetry has boosted their use in safety-critical and high-performance computing systems. As such systems require high levels of resilience to errors, CNNs must execute correctly in the presence of hardware faults. Full duplication provides the needed assurance but incurs a prohibitive 100% overhead. Algorithmic techniques are known to offer low-cost solutions, but the practical feasibility and performance of such techniques have never been studied for CNN deployment platforms (e.g., TensorFlow or TensorRT on GPUs). In this paper, we focus on algorithmically verifying Convolutions, which are the most resource-demanding operations in CNNs. We use checksums to verify convolutions, adding a small amount of redundancy, far less than full-duplication. We first identify the challenges that arise in employing Algorithm-Based Error Detection (ABED) for Convolutions in optimized inference platforms that fuse multiple network layers and use reduced-precision operations, and demonstrate how to overcome them. We propose and evaluate variations of ABED techniques that offer implementation complexity, runtime overhead, and coverage trade-offs. Results show that ABED can detect all transient hardware errors that might otherwise corrupt output and does so while incurring low runtime overheads (6-23%), offering at least 1.6X throughput to workloads compared to full duplication.


title: Medidas de distanciamento social e mobilidade na Am\'erica do Sul durante a pandemia por COVID-19: Condi\c{c}\~oes necess\'arias e suficientes?

abstract: In a scenario where there is no vaccine for COVID-19, non-pharmaceutical interventions are necessary to contain the spread of the virus and the collapse of the health system in the affected regions. One of these measures is social distancing, which aims to reduce interactions in the community by closing public and private establishments that involve crowds of people. The lockdown presupposes a drastic reduction in community interactions, representing a more extreme measure of social distancing. Based on geolocation data provided by Google for six categories of physical spaces, this article identifies the variations in the circulation of people in South America for different types of social distancing measures adopted during the COVID-19 pandemic. In this study, population mobility trends for a group of countries between February 15, 2020 and May 16, 2020 were analyzed. To summarize these trends in a single metric, a general circulation index was created, and to identify regional mobility patterns, descriptive analyzes of spatial autocorrelation (global and local Moran index) were used. The first hypothesis of this study is that countries with a lockdown decree can achieve greater success in reducing the mobility of the population, and the second hypothesis is that Argentina, Brazil and Colombia have regional mobility patterns. The first hypothesis was partially confirmed (considering 10 countries in South America), and the results obtained in the spatial analyzes confirmed the second hypothesis. In general, the observed data shows that less rigid lockdown or social distancing measures are necessary, however, they are not sufficient to achieve a significant reduction in the circulation of people during the pandemic.


title: 5* Knowledge Graph Embeddings with Projective Transformations

abstract: Performing link prediction using knowledge graph embedding models has become a popular approach for knowledge graph completion. Such models employ a transformation function that maps nodes via edges into a vector space in order to measure the likelihood of the links. While mapping the individual nodes, the structure of subgraphs is also transformed. Most of the embedding models designed in Euclidean geometry usually support a single transformation type - often translation or rotation, which is suitable for learning on graphs with small differences in neighboring subgraphs. However, multi-relational knowledge graphs often include multiple sub-graph structures in a neighborhood (e.g. combinations of path and loop structures), which current embedding models do not capture well. To tackle this problem, we propose a novel KGE model (5*E) in projective geometry, which supports multiple simultaneous transformations - specifically inversion, reflection, translation, rotation, and homothety. The model has several favorable theoretical properties and subsumes the existing approaches. It outperforms them on the most widely used link prediction benchmarks


title: Langevin dynamic for the 2D Yang-Mills measure

abstract: We define a natural state space and Markov process associated to the stochastic Yang-Mills heat flow in two dimensions. To accomplish this we first introduce a space of distributional connections for which holonomies along sufficiently regular curves (Wilson loop observables) and the action of an associated group of gauge transformations are both well-defined and satisfy good continuity properties. The desired state space is obtained as the corresponding space of orbits under this group action and is shown to be a Polish space when equipped with a natural Hausdorff metric. To construct the Markov process we show that the stochastic Yang-Mills heat flow takes values in our space of connections and use the "DeTurck trick" of introducing a time dependent gauge transformation to show invariance, in law, of the solution under gauge transformations. Our main tool for solving for the Yang-Mills heat flow is the theory of regularity structures and along the way we also develop a "basis-free" framework for applying the theory of regularity structures in the context of vector-valued noise - this provides a conceptual framework for interpreting several previous constructions and we expect this framework to be of independent interest.


title: Object Segmentation Without Labels with Large-Scale Generative Models

abstract: The recent rise of unsupervised and self-supervised learning has dramatically reduced the dependency on labeled data, providing effective image representations for transfer to downstream vision tasks. Furthermore, recent works employed these representations in a fully unsupervised setup for image classification, reducing the need for human labels on the fine-tuning stage as well. This work demonstrates that large-scale unsupervised models can also perform a more challenging object segmentation task, requiring neither pixel-level nor image-level labeling. Namely, we show that recent unsupervised GANs allow to differentiate between foreground/background pixels, providing high-quality saliency masks. By extensive comparison on standard benchmarks, we outperform existing unsupervised alternatives for object segmentation, achieving new state-of-the-art.


title: Synchronous Robotic Framework

abstract: We present a synchronous robotic testbed called SyROF that allows fast implementation of robotic swarms. Our main goal is to lower the entry barriers to cooperative-robot systems for undergraduate and graduate students. The testbed provides a high-level programming environment that allows the implementation of Timed Input/Output Automata (TIOA). SyROF offers the following unique characteristics: 1) a transparent mechanism to synchronize robot maneuvers, 2) a membership service with a failure detector, and 3) a transparent service to provide common knowledge in every round. These characteristics are fundamental to simplifying the implementation of robotic swarms. The software is organized in five layers: The lower layer consists of a real-time publish-subscribe system that allows efficient communication between tasks. The next layer is an implementation of a Kalman filter to estimate the position, orientation, and speed of the robot. The third layer consists of a synchronizer that synchronously executes the robot maneuvers, provides common knowledge to all the active participants, and handles failures. The fifth layer consists of the programming environment.


title: Coronal-Jet-Producing Minifilament Eruptions as a Possible Source of Parker Solar Probe (PSP) Switchbacks

abstract: The Parker Solar Probe (PSP) has observed copious rapid magnetic field direction changes in the near-Sun solar wind. These features have been called "switchbacks," and their origin is a mystery. But their widespread nature suggests that they may be generated by a frequently occurring process in the Sun's atmosphere. We examine the possibility that the switchbacks originate from coronal jets. Recent work suggests that many coronal jets result when photospheric magnetic flux cancels, and forms a small-scale "minifilament" flux rope that erupts and reconnects with coronal field. We argue that the reconnected erupting minifilament flux rope can manifest as an outward propagating Alfv\'enic fluctuation that steepens into an increasingly compact disturbance as it moves through the solar wind. Using previous observed properties of coronal jets that connect to coronagraph-observed white-light jets (a.k.a. "narrow CMEs"), along with typical solar wind speed values, we expect the coronal-jet-produced disturbances to traverse near-perihelion PSP in ~<25 min, with a velocity of ~400 km/s. To consider further the plausibility of this idea, we show that a previously studied series of equatorial latitude coronal jets, originating from the periphery of an active region, generate white-light jets in the outer corona (seen in STEREO/COR2 coronagraph images; 2.5---15 R_sun), and into the inner heliosphere (seen in STEREO/Hi1 heliospheric imager images; 15---84 R_sun). Thus it is tenable that disturbances put onto open coronal magnetic field lines by coronal-jet-producing erupting minifilament flux ropes can propagate out to PSP space and appear as switchbacks.


title: Beyond Triplet Loss: Meta Prototypical N-tuple Loss for Person Re-identification

abstract: Person Re-identification (ReID) aims at matching a person of interest across images. In convolutional neural network (CNN) based approaches, loss design plays a vital role in pulling closer features of the same identity and pushing far apart features of different identities. In recent years, triplet loss achieves superior performance and is predominant in ReID. However, triplet loss considers only three instances of two classes in per-query optimization (with an anchor sample as query) and it is actually equivalent to a two-class classification. There is a lack of loss design which enables the joint optimization of multiple instances (of multiple classes) within per-query optimization for person ReID. In this paper, we introduce a multi-class classification loss, i.e., N-tuple loss, to jointly consider multiple (N) instances for per-query optimization. This in fact aligns better with the ReID test/inference process, which conducts the ranking/comparisons among multiple instances. Furthermore, for more efficient multi-class classification, we propose a new meta prototypical N-tuple loss. With the multi-class classification incorporated, our model achieves the state-of-the-art performance on the benchmark person ReID datasets.


title: Deep Stock Predictions

abstract: Forecasting stock prices can be interpreted as a time series prediction problem, for which Long Short Term Memory (LSTM) neural networks are often used due to their architecture specifically built to solve such problems. In this paper, we consider the design of a trading strategy that performs portfolio optimization using the LSTM stock price prediction for four different companies. We then customize the loss function used to train the LSTM to increase the profit earned. Moreover, we propose a data driven approach for optimal selection of window length and multi-step prediction length, and consider the addition of analyst calls as technical indicators to a multi-stack Bidirectional LSTM strengthened by the addition of Attention units. We find the LSTM model with the customized loss function to have an improved performance in the training bot over a regressive baseline such as ARIMA, while the addition of analyst call does improve the performance for certain datasets.


title: On the stabilization of relative trace formulae: descent and the fundamental lemma

abstract: Motivated by the study of periods of automorphic forms and relative trace formulae, we develop the theory of descent necessary to study orbital integrals arising in the fundamental lemma for a general class of symmetric spaces over a $p$-adic field $F$. More precisely, we prove that a connected symmetric space over $F$ enjoys a notion of topological Jordan decomposition, which may be of independent interest, and establish a relative version of a lemma of Kazhdan that played a crucial role in the proof of the Langlands-Shelstad fundamental lemma. As our main application, we use these results to prove the endoscopic fundamental lemma for the unit element of the Hecke algebra for the symmetric space associated to unitary Friedberg-Jacquet periods.


title: On travelling wave solutions of a model of a liquid film flowing down a fibre

abstract: Existence of non-negative weak solutions is shown for a full curvature thin-film model of a liquid thin film flowing down a vertical fibre. The proof is based on the application of a priori estimates derived for energy-entropy functionals. Long-time behaviour of these weak solutions is analysed and, under some additional constraints for the model parameters and initial values, convergence towards a travelling wave solution is obtained. Numerical studies of energy minimizers and travelling waves are presented to illustrate analytical results.


title: Spoon or slide? The non-linear matter power spectrum in the presence of massive neutrinos

abstract: Numerical simulations of massive neutrino cosmologies consistently find a spoon-like feature in the non-linear matter power spectrum ratios of cosmological models that differ only in the neutrino mass fraction f_N. Typically, the ratio approaches unity at low wave numbers k, decreases by ~ 10 f_N at k ~ 1 h/Mpc, and turns up again at large k. Using the halo model of large-scale structure, we show that this spoon feature originates in the transition from the two-halo power spectrum to the one-halo power spectrum. The former's sensitivity to f_N rises with k, while that of the latter decreases with k. The presence of this spoon feature is robust with respect to different choices of the halo mass function and the halo density profile, and does not require any parameter tuning within the halo model. We demonstrate that a standard halo model calculation is already able to predict the depth, width, and position of this spoon as well as its evolution with redshift z with remarkable accuracy. Predictions at z >= 1 can be further improved using non-linear perturbative inputs.


title: Implicit Class-Conditioned Domain Alignment for Unsupervised Domain Adaptation

abstract: We present an approach for unsupervised domain adaptation---with a strong focus on practical considerations of within-domain class imbalance and between-domain class distribution shift---from a class-conditioned domain alignment perspective. Current methods for class-conditioned domain alignment aim to explicitly minimize a loss function based on pseudo-label estimations of the target domain. However, these methods suffer from pseudo-label bias in the form of error accumulation. We propose a method that removes the need for explicit optimization of model parameters from pseudo-labels directly. Instead, we present a sampling-based implicit alignment approach, where the sample selection procedure is implicitly guided by the pseudo-labels. Theoretical analysis reveals the existence of a domain-discriminator shortcut in misaligned classes, which is addressed by the proposed implicit alignment approach to facilitate domain-adversarial learning. Empirical results and ablation studies confirm the effectiveness of the proposed approach, especially in the presence of within-domain class imbalance and between-domain class distribution shift.


title: The Norton algebra of a $Q$-polynomial distance-regular graph

abstract: We consider the Norton algebra associated with a $Q$-polynomial primitive idempotent of the adjacency matrix for a distance-regular graph. We obtain a formula for the Norton algebra product that we find attractive.


title: Machine Learning Automatically Detects COVID-19 using Chest CTs in a Large Multicenter Cohort

abstract: Objectives: To investigate machine-learning classifiers and interpretable models using chest CT for detection of COVID-19 and differentiation from other pneumonias, ILD and normal CTs. Methods: Our retrospective multi-institutional study obtained 2096 chest CTs from 16 institutions (including 1077 COVID-19 patients). Training/testing cohorts included 927/100 COVID-19, 388/33 ILD, 189/33 other pneumonias, and 559/34 normal (no pathologies) CTs. A metric-based approach for classification of COVID-19 used interpretable features, relying on logistic regression and random forests. A deep learning-based classifier differentiated COVID-19 via 3D features extracted directly from CT attenuation and probability distribution of airspace opacities. Results: Most discriminative features of COVID-19 are percentage of airspace opacity and peripheral and basal predominant opacities, concordant with the typical characterization of COVID-19 in the literature. Unsupervised hierarchical clustering compares feature distribution across COVID-19 and control cohorts. The metrics-based classifier achieved AUC=0.83, sensitivity=0.74, and specificity=0.79 of versus respectively 0.93, 0.90, and 0.83 for the DL-based classifier. Most of ambiguity comes from non-COVID-19 pneumonia with manifestations that overlap with COVID-19, as well as mild COVID-19 cases. Non-COVID-19 classification performance is 91% for ILD, 64% for other pneumonias and 94% for no pathologies, which demonstrates the robustness of our method against different compositions of control groups. Conclusions: Our new method accurately discriminates COVID-19 from other types of pneumonia, ILD, and no pathologies CTs, using quantitative imaging features derived from chest CT, while balancing interpretability of results and classification performance, and therefore may be useful to facilitate diagnosis of COVID-19.


title: Robustness of slow contraction to cosmic initial conditions

abstract: We present numerical relativity simulations of cosmological scenarios in which the universe is smoothed and flattened by undergoing a phase of slow contraction and test their sensitivity to a wide range of initial conditions. Our numerical scheme enables the variation of all freely specifiable physical quantities that characterize the initial spatial hypersurface, such as the initial shear and spatial curvature contributions as well as the initial field and velocity distributions of the scalar that drives the cosmological evolution. In particular, we include initial conditions that are far outside the perturbative regime of the well-known attractor scaling solution. We complement our numerical results by analytically performing a complete dynamical systems analysis and show that the two approaches yield consistent results.


title: Student difficulties with the corrections to the energy spectrum of the hydrogen atom for the intermediate field Zeeman effect

abstract: We discuss an investigation of student difficulties with the corrections to the energy spectrum of the hydrogen atom for the intermediate field Zeeman effect using the degenerate perturbation theory. The investigation was carried out in advanced quantum mechanics courses by administering free-response and multiple-choice questions and conducting individual interviews with students. We find that students share many common difficulties related to relevant physics concepts. In particular, students often struggled with mathematical sense-making in this context of quantum mechanics which requires interpretation of the implications of degeneracy in the unperturbed energy spectrum and how the Zeeman perturbation will impact the splitting of the energy levels. We discuss how the common difficulties often arise from the fact that applying linear algebra concepts correctly in this context with degeneracy in the energy spectrum is challenging for students.


title: In Proximity of ReLU DNN, PWA Function, and Explicit MPC

abstract: Rectifier (ReLU) deep neural networks (DNN) and their connection with piecewise affine (PWA) functions is analyzed. The paper is an effort to find and study the possibility of representing explicit state feedback policy of model predictive control (MPC) as a ReLU DNN, and vice versa. The complexity and architecture of DNN has been examined through some theorems and discussions. An approximate method has been developed for identification of input-space in ReLU net which results a PWA function over polyhedral regions. Also, inverse multiparametric linear or quadratic programs (mp-LP or mp-QP) has been studied which deals with reconstruction of constraints and cost function given a PWA function.


title: Determination and estimation of optimal quarantine duration for infectious diseases with application to data analysis of COVID-19

abstract: Quarantine measure is a commonly used non-pharmaceutical intervention during the outbreak of infectious diseases. A key problem for implementing quarantine measure is to determine the duration of quarantine. In this paper, a policy with optimal quarantine duration is developed. The policy suggests different quarantine durations for every individual with different characteristic. The policy is optimal in the sense that it minimizes the average quarantine duration of uninfected people with the constraint that the probability of symptom presentation for infected people attains the given value closing to 1. The optimal solution for the quarantine duration is obtained and estimated by some statistic methods with application to analyzing COVID-19 data.


title: Universal Vector Neural Machine Translation With Effective Attention

abstract: Neural Machine Translation (NMT) leverages one or more trained neural networks for the translation of phrases. Sutskever introduced a sequence to sequence based encoder-decoder model which became the standard for NMT based systems. Attention mechanisms were later introduced to address the issues with the translation of long sentences and improving overall accuracy. In this paper, we propose a singular model for Neural Machine Translation based on encoder-decoder models. Most translation models are trained as one model for one translation. We introduce a neutral/universal model representation that can be used to predict more than one language depending on the source and a provided target. Secondly, we introduce an attention model by adding an overall learning vector to the multiplicative model. With these two changes, by using the novel universal model the number of models needed for multiple language translation applications are reduced.


title: Global asymptotic behavior of solutions to a class of Kirchhoff equations

abstract: In this paper, a parabolic type Kirchhoff equation and its stationary counterpart are considered. For the evolution problem, the precise decay rates of the weak solution and of the corresponding energy functional are derived. For the stationary problem, a ground-state solution is obtained by applying Lagrange multiplier method. Moreover, the asymptotic behaviors of the general global solutions are also described. These results extend some recent ones obtained in [Threshold results for the existence of global and blow-up solutions to Kirchhoff equations with arbitrary initial energy, Computers and Mathematics with Applications, 75(2018), 3283-3297] by Han and Li.


title: Blow-up phenomena for a reaction diffusion equation with special diffusion process

abstract: This paper is concerned with the blow-up property of solutions to an initial boundary value problem for a reaction diffusion equation with special diffusion processes. It is shown, under certain conditions on the initial data, that the solutions to this problem blow up in finite time, by combining Hardy inequality, "moving" potential well methods with some differential inequalities. Moreover, the upper and lower bounds for the blow-up time are also derived when blow-up occurs.


title: Lifespan of solutions to a damped fourth-order wave equation with logarithmic nonlinearity

abstract: This paper is devoted to the lifespan of solutions to a damped fourth-order wave equation with logarithmic nonlinearity $$u_{tt}+\Delta^2u-\Delta u-\omega\Delta u_t+\alpha(t)u_t=|u|^{p-2}u\ln|u|.$$ Finite time blow-up criteria for solutions at both lower and high initial energy levels are established, and an upper bound for the blow-up time is given for each case. Moreover, by constructing a new auxiliary functional and making full use of the strong damping term, a lower bound for the blow-up time is also derived.


title: The Hitchhiker's Guide to the All-Interval 12-Tone Rows

abstract: This article revisits the generation, classification and categorization of all-intervals 12-tone series (AIS). Inspired by the seminal work of Morris and Starr in 1974 (Morris and Starr, The Structure of All-Interval Series 1974), it expands their analysis using complex network theory and provides composers and theorists with the re-ordering scheme that links all AISs together by chains of relations.


title: Staggered explicit-implicit time-discretization for elastodynamics with dissipative internal variables

abstract: An extension of the two-step staggered time discretization of linear elastodynamics in stress-velocity form to systems involving internal variables subjected to a possibly non-linear dissipative evolution is proposed. The original scheme is thus enhanced by another step for the internal variables which, in general, is implicit, although even this step might be explicit if no spatial gradients of the internal variables are involved. Using an abstract Banach-space formulation, a-priori estimates and convergence are proved under a CFL condition. The developed three-step scheme finds applications in various problems of continuum mechanics at small strain. Here, we consider in particular plasticity, viscoelasticity (creep), diffusion in poroelastic media, and damage.


title: Few-Shot Generative Conversational Query Rewriting

abstract: Conversational query rewriting aims to reformulate a concise conversational query to a fully specified, context-independent query that can be effectively handled by existing information retrieval systems. This paper presents a few-shot generative approach to conversational query rewriting. We develop two methods, based on rules and self-supervised learning, to generate weak supervision data using large amounts of ad hoc search sessions, and to fine-tune GPT-2 to rewrite conversational queries. On the TREC Conversational Assistance Track, our weakly supervised GPT-2 rewriter improves the state-of-the-art ranking accuracy by 12%, only using very limited amounts of manual query rewrites. In the zero-shot learning setting, the rewriter still gives a comparable result to previous state-of-the-art systems. Our analyses reveal that GPT-2 effectively picks up the task syntax and learns to capture context dependencies, even for hard cases that involve group references and long-turn dependencies.


title: Reinforcement Learning-Based Joint Self-Optimisation Method for the Fuzzy Logic Handover Algorithm in 5G HetNets

abstract: 5G heterogeneous networks (HetNets) can provide higher network coverage and system capacity to the user by deploying massive small base stations (BSs) within the 4G macro system. However, the large-scale deployment of small BSs significantly increases the complexity and workload of network maintenance and optimisation. The current handover (HO) triggering mechanism A3 event was designed only for mobility management in the macro system. Directly implementing A3 in 5G-HetNets may degrade the user mobility robustness. Motivated by the concept of self-organisation networks (SON), this study developed a self-optimised triggering mechanism to enable automated network maintenance and enhance user mobility robustness in 5G-HetNets. The proposed method integrates the advantages of subtractive clustering and Q-learning frameworks into the conventional fuzzy logic-based HO algorithm (FLHA). Subtractive clustering is first adopted to generate a membership function (MF) for the FLHA to enable FLHA with the self-configuration feature. Subsequently, Q-learning is utilised to learn the optimal HO policy from the environment as fuzzy rules that empower the FLHA with a self-optimisation function. The FLHA with SON functionality also overcomes the limitations of the conventional FLHA that must rely heavily on professional experience to design. The simulation results show that the proposed self-optimised FLHA can effectively generate MF and fuzzy rules for the FLHA. By comparing with conventional triggering mechanisms, the proposed approach can decrease the HO, ping-pong HO, and HO failure ratios by approximately 91%, 49%, and 97.5% while improving network throughput and latency by 8% and 35%, respectively.


title: RGB-D-E: Event Camera Calibration for Fast 6-DOF Object Tracking

abstract: Augmented reality devices require multiple sensors to perform various tasks such as localization and tracking. Currently, popular cameras are mostly frame-based (e.g. RGB and Depth) which impose a high data bandwidth and power usage. With the necessity for low power and more responsive augmented reality systems, using solely frame-based sensors imposes limits to the various algorithms that needs high frequency data from the environement. As such, event-based sensors have become increasingly popular due to their low power, bandwidth and latency, as well as their very high frequency data acquisition capabilities. In this paper, we propose, for the first time, to use an event-based camera to increase the speed of 3D object tracking in 6 degrees of freedom. This application requires handling very high object speed to convey compelling AR experiences. To this end, we propose a new system which combines a recent RGB-D sensor (Kinect Azure) with an event camera (DAVIS346). We develop a deep learning approach, which combines an existing RGB-D network along with a novel event-based network in a cascade fashion, and demonstrate that our approach significantly improves the robustness of a state-of-the-art frame-based 6-DOF object tracker using our RGB-D-E pipeline.


title: Measurement of the Proton-Air Cross Section with Telescope Array's Black Rock Mesa and Long Ridge Fluorescence Detectors, and Surface Array in Hybrid Mode

abstract: Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $\sigma^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $\sigma^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb.


title: A Random Matrix Analysis of Random Fourier Features: Beyond the Gaussian Kernel, a Precise Phase Transition, and the Corresponding Double Descent

abstract: This article characterizes the exact asymptotics of random Fourier feature (RFF) regression, in the realistic setting where the number of data samples $n$, their dimension $p$, and the dimension of feature space $N$ are all large and comparable. In this regime, the random RFF Gram matrix no longer converges to the well-known limiting Gaussian kernel matrix (as it does when $N \to \infty$ alone), but it still has a tractable behavior that is captured by our analysis. This analysis also provides accurate estimates of training and test regression errors for large $n,p,N$. Based on these estimates, a precise characterization of two qualitatively different phases of learning, including the phase transition between them, is provided; and the corresponding double descent test error curve is derived from this phase transition behavior. These results do not depend on strong assumptions on the data distribution, and they perfectly match empirical results on real-world data sets.


title: HausaMT v1.0: Towards English-Hausa Neural Machine Translation

abstract: Neural Machine Translation (NMT) for low-resource languages suffers from low performance because of the lack of large amounts of parallel data and language diversity. To contribute to ameliorating this problem, we built a baseline model for English-Hausa machine translation, which is considered a task for low-resource language. The Hausa language is the second largest Afro-Asiatic language in the world after Arabic and it is the third largest language for trading across a larger swath of West Africa countries, after English and French. In this paper, we curated different datasets containing Hausa-English parallel corpus for our translation. We trained baseline models and evaluated the performance of our models using the Recurrent and Transformer encoder-decoder architecture with two tokenization approaches: standard word-level tokenization and Byte Pair Encoding (BPE) subword tokenization.


title: Can Synthetic Data Improve Object Detection Results for Remote Sensing Images?

abstract: Deep learning approaches require enough training samples to perform well, but it is a challenge to collect enough real training data and label them manually. In this letter, we propose the use of realistic synthetic data with a wide distribution to improve the performance of remote sensing image aircraft detection. Specifically, to increase the variability of synthetic data, we randomly set the parameters during rendering, such as the size of the instance and the class of background images. In order to make the synthetic images more realistic, we then refine the synthetic images at the pixel level using CycleGAN with real unlabeled images. We also fine-tune the model with a small amount of real data, to obtain a higher accuracy. Experiments on NWPU VHR-10, UCAS-AOD and DIOR datasets demonstrate that the proposed method can be applied for augmenting insufficient real data.


title: Reproducible and Portable Workflows for Scientific Computing and HPC in the Cloud

abstract: The increasing availability of cloud computing services for science has changed the way scientific code can be developed, deployed, and run. Many modern scientific workflows are capable of running on cloud computing resources. Consequently, there is an increasing interest in the scientific computing community in methods, tools, and implementations that enable moving an application to the cloud and simplifying the process, and decreasing the time to meaningful scientific results. In this paper, we have applied the concepts of containerization for portability and multi-cloud automated deployment with industry-standard tools to three scientific workflows. We show how our implementations provide reduced complexity to portability of both the applications themselves, and their deployment across private and public clouds. Each application has been packaged in a Docker container with its dependencies and necessary environment setup for production runs. Terraform and Ansible have been used to automate the provisioning of compute resources and the deployment of each scientific application in a Multi-VM cluster. Each application has been deployed on the AWS and Aristotle Cloud Federation platforms. Variation in data management constraints, Multi-VM MPI communication, and embarrassingly parallel instance deployments were all explored and reported on. We thus present a sample of scientific workflows that can be simplified using the tools and our proposed implementation to deploy and run in a variety of cloud environments.


title: Photometry and high-resolution spectroscopy of comet 21P/Giacobini-Zinner during its 2018 apparition

abstract: We report on photometry and high resolution spectroscopy of the chemically peculiar Jupiter-family Comet (hereafter JFC) 21P/Giacobini-Zinner. Comet 21P is a well known member of the carbon-chain depleted family but displays also a depletion of amines. We monitored continuously the comet over more than seven months with the two TRAPPIST telescopes (TN and TS), covering a large heliocentric distance range from 1.60 au inbound to 2.10 au outbound with a perihelion at 1.01 au on September 10, 2018. We computed and followed the evolution of the dust (represented by Af$\rho$) and gas production rates of the daughter species OH, NH, CN, C$_3$, and C$_2$ and their relative abundances to OH and to CN over the comet orbit. We compared them to those measured in the previous apparitions. The activity of the comet and its water production rate reached a maximum of (3.72$\pm$0.07)$\times$10$^{28}$ molec/s on August 17, 2018 (r$_h$=1.07 au), 24 days before perihelion. The peak value of A(0)f$\rho$ was reached on the same date (1646$\pm$13) cm in the red filter. The abundance ratios of the various species are remarkably constant over a large range of heliocentric distances, before and after perihelion, showing a high level of homogeneity of the ices in the surface of the nucleus. The behaviour and level of the activity of the comet is also remarkably similar over the last five orbits. About the coma dust colour, 21P shows reflectively gradients similar to JFCs. We obtained a high resolution spectrum of 21P with UVES at ESO VLT one week after perihelion. Using the CN B-X (0,0) violet band, we measured $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N isotopic ratios of 100$\pm$10 and 145$\pm$10, respectively, both in very good agreement with what is usually found in comets.


title: Deep learning to estimate the physical proportion of infected region of lung for COVID-19 pneumonia with CT image set

abstract: Utilizing computed tomography (CT) images to quickly estimate the severity of cases with COVID-19 is one of the most straightforward and efficacious methods. Two tasks were studied in this present paper. One was to segment the mask of intact lung in case of pneumonia. Another was to generate the masks of regions infected by COVID-19. The masks of these two parts of images then were converted to corresponding volumes to calculate the physical proportion of infected region of lung. A total of 129 CT image set were herein collected and studied. The intrinsic Hounsfiled value of CT images was firstly utilized to generate the initial dirty version of labeled masks both for intact lung and infected regions. Then, the samples were carefully adjusted and improved by two professional radiologists to generate the final training set and test benchmark. Two deep learning models were evaluated: UNet and 2.5D UNet. For the segment of infected regions, a deep learning based classifier was followed to remove unrelated blur-edged regions that were wrongly segmented out such as air tube and blood vessel tissue etc. For the segmented masks of intact lung and infected regions, the best method could achieve 0.972 and 0.757 measure in mean Dice similarity coefficient on our test benchmark. As the overall proportion of infected region of lung, the final result showed 0.961 (Pearson's correlation coefficient) and 11.7% (mean absolute percent error). The instant proportion of infected regions of lung could be used as a visual evidence to assist clinical physician to determine the severity of the case. Furthermore, a quantified report of infected regions can help predict the prognosis for COVID-19 cases which were scanned periodically within the treatment cycle.


title: Handover Rate and Sojourn Time Analysis in Mobile Drone-Assisted Cellular Networks

abstract: To improve capacity and overcome some of the limitations of cellular wireless networks, drones with aerial base stations can be deployed to assist the terrestrial cellular wireless networks. The mobility of drones allows flexible network reconfiguration to adapt to dynamic traffic and channel conditions. However, this is achieved at the expense of more handovers since even a static user may experience a handover when the drones are mobile. In this letter, we provide an exact analysis of the handover rate and sojourn time (time between two subsequent handovers) for a network of drone base stations. We also show that among different speed distributions with the same mean, the handover rate is minimum when all drone base stations move with same speed.


title: A functional-data approach to the Argo data

abstract: The Argo data is a modern oceanography dataset that provides unprecedented global coverage of temperature and salinity measurements in the upper 2,000 meters of depth of the ocean. We study the Argo data from the perspective of functional data analysis (FDA). We develop spatio-temporal functional kriging methodology for mean and covariance estimation to predict temperature and salinity at a fixed location as a smooth function of depth. By combining tools from FDA and spatial statistics, including smoothing splines, local regression, and multivariate spatial modeling and prediction, our approach provides advantages over current methodology that consider pointwise estimation at fixed depths. Our approach naturally leverages the irregularly-sampled data in space, time, and depth to fit a space-time functional model for temperature and salinity. The developed framework provides new tools to address fundamental scientific problems involving the entire upper water column of the oceans such as the estimation of ocean heat content, stratification, and thermohaline oscillation. For example, we show that our functional approach yields more accurate ocean heat content estimates than ones based on discrete integral approximations in pressure. Further, using the derivative function estimates, we obtain a new product of a global map of the mixed layer depth, a key component in the study of heat absorption and nutrient circulation in the oceans. The derivative estimates also reveal evidence for density inversions in areas distinguished by mixing of particularly different water masses.


title: CLAIMED: A CLAssification-Incorporated Minimum Energy Design to explore a multivariate response surface with feasibility constraints

abstract: Motivated by the problem of optimization of force-field systems in physics using large-scale computer simulations, we consider exploration of a deterministic complex multivariate response surface. The objective is to find input combinations that generate output close to some desired or "target" vector. In spite of reducing the problem to exploration of the input space with respect to a one-dimensional loss function, the search is nontrivial and challenging due to infeasible input combinations, high dimensionalities of the input and output space and multiple "desirable" regions in the input space and the difficulty of emulating the objective function well with a surrogate model. We propose an approach that is based on combining machine learning techniques with smart experimental design ideas to locate multiple good regions in the input space.


title: Near-Optimal Confidence Sequences for Bounded Random Variables

abstract: Many inference problems, such as sequential decision problems like A/B testing, adaptive sampling schemes like bandit selection, are often online in nature. The fundamental problem for online inference is to provide a sequence of confidence intervals that are valid uniformly over the growing-into-infinity sample sizes. To address this question, we provide a near-optimal confidence sequence for bounded random variables by utilizing Bentkus' concentration results. We show that it improves on the existing approaches that use the Cram{\'e}r-Chernoff technique such as the Hoeffding, Bernstein, and Bennett inequalities. The resulting confidence sequence is confirmed to be favorable in both synthetic coverage problems and an application to adaptive stopping algorithms.

categories: math.ST cs.AI cs.LG stat.AP stat.ML stat.TH