title: Fermion transfer in the $\phi^4$ model with a half-BPS preserving impurity

abstract: We study a fermion field coupled to a scalar via a Yukawa term. The scalar field is the $\phi^4$ model with an impurity that preserves half of the BPS property. We analyze the spectrum of the defects of the model and collisions between them both close to the BPS regime and not. As the fermion binds to these defects, it may be transferred from one to the other, which we quantify via overlaps, known as Bogoliubov coefficients. BPS collisions are less likely to transfer the fermion between defects and can be adiabatic for non-relativistic velocities, especially for small coupling constants. Moreover, closer to the BPS limit only a small fraction of the fermion number is radiated away. In contrast, non-BPS collisions lead to more radiation in the fermion field and excitation of the fermion to higher bound states, and the result is more sensitive to the parameters.


title: The dependence of global super-rotation on planetary rotation rate

abstract: An atmosphere may be described as globally super-rotating if its total zonal angular momentum exceeds that associated with solid-body co-rotation with the underlying planet. In this paper, we discuss the dependence of global super-rotation in terrestrial atmospheres on planetary rotation rate. This dependence is revealed through analysis of global super-rotation in idealised General Circulation Model experiments with time-independent axisymmetric forcing, compared with estimates for global super-rotation in Solar System atmospheres. Axisymmetric and three-dimensional experiments are conducted. We find that the degree of global super-rotation in the three-dimensional experiments is closely related to that of the axisymmetric experiments, with some differences in detail. A scaling theory for global super-rotation in an axisymmetric atmosphere is derived from the Held-Hou model. At high rotation rate, our numerical experiments inhabit a regime where global super-rotation scales geostrophically, and we suggest that the Earth and Mars occupy this regime. At low rotation rate, our experiments occupy a regime determined by angular momentum conservation, where global super-rotation is independent of rotation rate. Global super-rotation in our experiments saturates at a value significantly lower than that achieved in the atmospheres of Venus and Titan, which instead occupy a regime where global super-rotation scales cyclostrophically. This regime can only be accessed when eddy induced up-gradient angular momentum transport is sufficiently large, which is not the case in our idealised numerical experiments. We suggest that the 'default' regime for a slowly rotating planet is the angular momentum conserving regime, characterised by mild global (and local) superrotation.


title: High-resolution Transmission Spectroscopy of MASCARA-2 b with EXPRES

abstract: We report detections of atomic species in the atmosphere of MASCARA-2 b, using the first transit observations obtained with the newly commissioned EXPRES spectrograph. EXPRES is a highly stabilised optical echelle spectrograph, designed to detect stellar reflex motions with amplitudes down to 30 cm/s, and was recently deployed at the Lowell Discovery Telescope. By analysing the transmission spectrum of the ultra-hot Jupiter MASCARA-2 b using the cross-correlation method, we confirm previous detections of Fe I, Fe II and Na I, which likely originate in the upper regions of the inflated atmosphere. In addition, we report significant detections of Mg I and Cr II. The absorption strengths change slightly with time, possibly indicating different temperatures and chemistry in the day-side and night-side terminators. Using the effective stellar line-shape variation induced by the transiting planet, we constrain the projected spin-orbit misalignment of the system to $1.6\pm3.1$ degrees, consistent with an aligned orbit. We demonstrate that EXPRES joins a suite of instruments capable of phase-resolved spectroscopy of exoplanet atmospheres.


title: Large-scale modelling and forecasting of ambulance calls in northern Sweden using spatio-temporal log-Gaussian Cox processes

abstract: Although ambulance call data typically come in the form of spatio-temporal point patterns, point process-based modelling approaches presented in the literature are scarce. In this paper, we study a unique set of Swedish spatio-temporal ambulance call data, which consist of the spatial (GPS) locations of the calls (within the four northernmost regions of Sweden) and the associated days of occurrence of the calls (January 1, 2014, to December 31, 2018). Motivated by the nature of the data, we here employ log-Gaussian Cox processes (LGCPs) for the spatio-temporal modelling and forecasting of the calls. To this end, we propose a K-means clustering based bandwidth selection method for the kernel estimation of the spatial component of the separable spatio-temporal intensity function. The temporal component of the intensity function is modelled using Poisson regression, using different calendar covariates, and the spatio-temporal random field component of the random intensity of the LGCP is fitted using the Metropolis-adjusted Langevin algorithm. Spatial hot-spots have been found in the south-eastern part of the study region, where most people in the region live and our fitted model/forecasts manage to capture this behavior quite well. Also, there is a significant association between the expected number of calls and the day-of-the-week and the season-of-the-year. A non-parametric second-order analysis indicates that LGCPs seem to be reasonable models for the data. Finally, we find that the fitted forecasts generate simulated future spatial event patterns that quite well resemble the actual future data.


title: Identifying Weakly Connected Subsystems in Building Energy Model for Effective Load Estimation in Presence of Parametric Uncertainty

abstract: It is necessary to estimate the expected energy usage of a building to determine how to reduce energy usage. The expected energy usage of a building can be reliably simulated using a Building Energy Model (BEM). Many of the numerous input parameters in a BEM are uncertain. To ensure that the building simulation is sufficiently accurate, and to better understand the impact of imprecisions in the input parameters and calculation methods, it is desirable to quantify uncertainty in the BEM throughout the modeling process. Uncertainty quantification (UQ) typically requires a large number of simulations to produce meaningful data, which, due to the vast number of input parameters and the dynamic nature of building simulation, is computationally expensive. Uncertainty Quantification (UQ) in BEM domain is thus intractable due to the size of the problem and parameters involved and hence it needs an advanced methodology for analysis. The current paper outlines a novel Weakly-Connected-Systems (WCSs) identification-based UQ framework developed to propagate the quantifiable uncertainty in the BEM. The overall approach is demonstrated on the physics-based thermal model of an actual building in Central New York.


title: Constraints on $H^\pm$ parameter space in 2HDM at $\sqrt{s}=$ 8 TeV and $\sqrt{s}=$ 13 TeV

abstract: This paper reflects the heavy Higgs scenario where the mass of charged Higgs is equal to or greater than 200 GeV. The CMS observed and expected values of upper limits on the product $\sigma_H^\pm BR(H^\pm \rightarrow tb^\mp)$, assuming $H^\pm \rightarrow tb^\mp=1$, both at 8 TeV (at integrated luminosity of 19.7 $fb^{-1}$ ) and 13 TeV (at integrated luminosity of 35.9 $fb^{-1}$ ) c.m energies are used. By comparing these expected and observed upper limits with computational values , we find out the expected and observed exclusion regions of charged Higgs parameter space ($ m_H^\pm - tan\beta $ space ) in 2HDM both at $\sqrt{s}=$8 and $\sqrt{s}=$ 13 TeV. We compare the expected and observed exclusion regions and observe that exclusion regions made by observed upper limits are always greater than the exclusion made by expected upper limits both at 8 and 13 TeV c.m energies. Only in the mass range from 200 GeV to 220 GeV the expected exclusion region is greater than the observed one only at $\sqrt{s}=$13 TeV. We also equate the exclusion regions at these two different center of mass energies and find that the expected exclusion region and observed exclusion region at $\sqrt{s}=$13 TeV are always greater than the expected exclusion region and observed exclusion region at $\sqrt{S}=$8 TeV respectively.


title: Angular momentum balance and vortex production in wall-bounded flows

abstract: To produce a vortex, a torque must be applied to the fluid. In viscous fluids, the torques that produce turbulent vortices result from the loss of symmetry of the stress tensor, once the viscous friction exceeds the shear stress resistance of the fluid. In wall-bounded flows, in particular, the turbulent vortices form in a thin layer of fluid adjacent to the wall, practically coinciding with the so-called viscous sublayer, where the viscous friction reaches the largest values. The present paper determines a vortex structure for this sublayer, consistent with the well-known linearity of the diagram of the mean streamwise velocity of this region. The analysis enables us to calculate the diameter, angular velocity, and interaxis of the vortices in the viscous sublayer in steady-state conditions. The lifting force that makes the vortices migrate from the wall towards the mainstream flow is determined, and the crucial role played by gyroscopic precession in the reorientation of the vortex axis is discussed.


title: Advanced Equivalence Checking for Quantum Circuits

abstract: Quantum computing will change the way we tackle certain problems. It promises to dramatically speed-up many chemical, financial, and machine-learning applications. However, to capitalize on those promises, complex design flows composed of steps such as compilation, decomposition, or mapping need to be employed before being able to execute a conceptual quantum algorithm on an actual device. This results in descriptions at various levels of abstraction which may significantly differ from each other. The complexity of the underlying design problems necessitates to not only provide efficient solutions for the single steps, but also to verify that the originally intended functionality is preserved throughout all levels of abstraction. This motivates methods for equivalence checking of quantum circuits. However, most existing methods are inspired by the classical realm and have merely been extended to support quantum circuits (i.e., circuits which do not only rely on 0's and 1's, but also employ superposition and entanglement). In this work, we propose an advanced methodology which takes the different paradigms of quantum circuits not only as a burden, but as an opportunity. In fact, the proposed methodology explicitly utilizes characteristics unique to quantum computing in order to overcome the shortcomings of existing approaches. We show that, by exploiting the reversibility of quantum circuits, complexity can be kept feasible in many cases. Moreover, we show that, in contrast to the classical realm, simulation is very powerful in verifying quantum circuits. Experimental evaluations confirm that the resulting methodology allows one to conduct equivalence checking dramatically faster than ever before--in many cases just a single simulation run is sufficient. An implementation of the proposed methodology is publicly available at https://iic.jku.at/eda/research/quantum_verification/.


title: On the attractor zero of a sequence of polynomials

abstract: The main perpose of this paper is to sudy the roots of a familly of polynomials that arise from a linear recurrences associated to Pascal's triangle and their zero attractor, using an analytical methods based on conformal mappings.


title: Nuclear global spaces of ultradifferentiable functions in the matrix weighted setting

abstract: We prove that the Hermite functions are an absolute Schauder basis for many global weighted spaces of ultradifferentiable functions in the matrix weighted setting and we determine also the corresponding coefficient spaces, thus extending previous work by Langenbruch. As a consequence we give very general conditions for these spaces to be nuclear. In particular, we obtain the corresponding results for spaces defined by weight functions.


title: Fitting the Search Space of Weight-sharing NAS with Graph Convolutional Networks

abstract: Neural architecture search has attracted wide attentions in both academia and industry. To accelerate it, researchers proposed weight-sharing methods which first train a super-network to reuse computation among different operators, from which exponentially many sub-networks can be sampled and efficiently evaluated. These methods enjoy great advantages in terms of computational costs, but the sampled sub-networks are not guaranteed to be estimated precisely unless an individual training process is taken. This paper owes such inaccuracy to the inevitable mismatch between assembled network layers, so that there is a random error term added to each estimation. We alleviate this issue by training a graph convolutional network to fit the performance of sampled sub-networks so that the impact of random errors becomes minimal. With this strategy, we achieve a higher rank correlation coefficient in the selected set of candidates, which consequently leads to better performance of the final architecture. In addition, our approach also enjoys the flexibility of being used under different hardware constraints, since the graph convolutional network has provided an efficient lookup table of the performance of architectures in the entire search space.


title: PySINDy: A Python package for the Sparse Identification of Nonlinear Dynamics from Data

abstract: PySINDy is a Python package for the discovery of governing dynamical systems models from data. In particular, PySINDy provides tools for applying the sparse identification of nonlinear dynamics (SINDy) (Brunton et al. 2016) approach to model discovery. In this work we provide a brief description of the mathematical underpinnings of SINDy, an overview and demonstration of the features implemented in PySINDy (with code examples), practical advice for users, and a list of potential extensions to PySINDy. Software is available at https://github.com/dynamicslab/pysindy.


title: Automated System Performance Testing at MongoDB

abstract: Distributed Systems Infrastructure (DSI) is MongoDB's framework for running fully automated system performance tests in our Continuous Integration (CI) environment. To run in CI it needs to automate everything end-to-end: provisioning and deploying multi-node clusters, executing tests, tuning the system for repeatable results, and collecting and analyzing the results. Today DSI is MongoDB's most used and most useful performance testing tool. It runs almost 200 different benchmarks in daily CI, and we also use it for manual performance investigations. As we can alert the responsible engineer in a timely fashion, all but one of the major regressions were fixed before the 4.2.0 release. We are also able to catch net new improvements, of which DSI caught 17. We open sourced DSI in March 2020.


title: Organ at Risk Segmentation for Head and Neck Cancer using Stratified Learning and Neural Architecture Search

abstract: OAR segmentation is a critical step in radiotherapy of head and neck (H&N) cancer, where inconsistencies across radiation oncologists and prohibitive labor costs motivate automated approaches. However, leading methods using standard fully convolutional network workflows that are challenged when the number of OARs becomes large, e.g. > 40. For such scenarios, insights can be gained from the stratification approaches seen in manual clinical OAR delineation. This is the goal of our work, where we introduce stratified organ at risk segmentation (SOARS), an approach that stratifies OARs into anchor, mid-level, and small & hard (S&H) categories. SOARS stratifies across two dimensions. The first dimension is that distinct processing pipelines are used for each OAR category. In particular, inspired by clinical practices, anchor OARs are used to guide the mid-level and S&H categories. The second dimension is that distinct network architectures are used to manage the significant contrast, size, and anatomy variations between different OARs. We use differentiable neural architecture search (NAS), allowing the network to choose among 2D, 3D or Pseudo-3D convolutions. Extensive 4-fold cross-validation on 142 H&N cancer patients with 42 manually labeled OARs, the most comprehensive OAR dataset to date, demonstrates that both pipeline- and NAS-stratification significantly improves quantitative performance over the state-of-the-art (from 69.52% to 73.68% in absolute Dice scores). Thus, SOARS provides a powerful and principled means to manage the highly complex segmentation space of OARs.


title: Consistent evaluation of continuum scale properties of Graphene

abstract: We handshake statistical mechanics with continuum mechanics to develop a methodology for consistent evaluation of the continuum scale properties of graphene. The scope is kept limited to elastic modulus, $E$, which has been reported to vary between 0.912 TPa to 7 TPa, Poisson's ratio, $\nu$, which has been reported to vary from being negative to a value as large as 0.46, and effective thickness, $q$, whose value varies between 0.75 \AA to 3.41 \AA. Such a large scatter arises due to inconsistent evaluation of these properties and making assumptions that may not be valid at atomistic scales. Our methodology combines three separate methods -- uniaxial tension, equibiaxial tension, and flexural out-of-plane free vibrations of simply supported sheets, which, when used in tandem in MD, can provide consistent values of $E, \nu$ and $q$. The only assumption made in the present study is the validity of the continuum scale thin plate vibration equation to represent the free vibrations of a long graphene sheet. Our results suggest that -- (i) graphene is auxetic with its Poisson's ratio increasing with increasing temperature, (ii) with increasing temperature, $E$ decreases, and (iii) the effective thickness increases with temperature.


title: Renormalizing individual performance metrics for cultural heritage management of sports records

abstract: Individual performance metrics are commonly used to compare players from different eras. However, such cross-era comparison is often biased due to significant changes in success factors underlying player achievement rates (e.g. performance enhancing drugs and modern training regimens). Such historical comparison is more than fodder for casual discussion among sports fans, as it is also an issue of critical importance to the multi-billion dollar professional sport industry and the institutions (e.g. Hall of Fame) charged with preserving sports history and the legacy of outstanding players and achievements. To address this cultural heritage management issue, we report an objective statistical method for renormalizing career achievement metrics, one that is particularly tailored for common seasonal performance metrics, which are often aggregated into summary career metrics -- despite the fact that many player careers span different eras. Remarkably, we find that the method applied to comprehensive Major League Baseball and National Basketball Association player data preserves the overall functional form of the distribution of career achievement, both at the season and career level. As such, subsequent re-ranking of the top-50 all-time records in MLB and the NBA using renormalized metrics indicates reordering at the local rank level, as opposed to bulk reordering by era. This local order refinement signals time-independent mechanisms underlying annual and career achievement in professional sports, meaning that appropriately renormalized achievement metrics can be used to compare players from eras with different season lengths, team strategies, rules -- and possibly even different sports.


title: A Proof of Gr\"unbaum's Lower Bound Conjecture for general polytopes

abstract: In 1967, Gr\"unbaum conjectured that any $d$-dimensional polytope with $d+s\leq 2d$ vertices has at least \[\phi_k(d+s,d) = {d+1 \choose k+1 }+{d \choose k+1 }-{d+1-s \choose k+1 } \] $k$-faces. We prove this conjecture and also characterize the cases in which equality holds.


title: An averaging principle for fractional stochastic differential equations with L\'evy noise

abstract: This paper is devoted to the study of an averaging principle for fractional stochastic differential equations in Rnwith L\'evy motion, using an integral transform method. We obtain a time-averaged equation under suitable assumptions. Furthermore, we show that the solutions of averaged equation approach the solutions of the original equation. Our results in this paper provide better understanding for effective approximation of fractional dynamical systems with non-Gaussian L\'evy noise.


title: Outlier detection at the parcel-level in wheat and rapeseed crops using multispectral and SAR time series

abstract: This paper studies the detection of anomalous crop development at the parcel-level based on an unsupervised outlier detection technique. The experimental validation is conducted on rapeseed and wheat parcels located in Beauce (France). The proposed methodology consists of four sequential steps: 1) preprocessing of synthetic aperture radar (SAR) and multispectral images acquired using Sentinel-1 and Sentinel-2 satellites, 2) extraction of SAR and multispectral pixel-level features, 3) computation of parcel-level features using zonal statistics and 4) outlier detection. The different types of anomalies that can affect the studied crops are analyzed and described. The different factors that can influence the outlier detection results are investigated with a particular attention devoted to the synergy between Sentinel-1 and Sentinel-2 data. Overall, the best performance is obtained when using jointly a selection of Sentinel-1 and Sentinel-2 features with the isolation forest algorithm. The selected features are VV and VH backscattering coefficients for Sentinel-1 and 5 Vegetation Indexes for Sentinel-2 (among us, the Normalized Difference Vegetation Index and two variants of the Normalized Difference Water). When using these features with an outlier ratio of 10%, the percentage of detected true positives (i.e., crop anomalies) is equal to 94.1% for rapeseed parcels and 95.5% for wheat parcels.


title: Fully-Dynamic Graph Sparsifiers Against an Adaptive Adversary

abstract: Designing dynamic graph algorithms against an adaptive adversary is a major goal in the field of dynamic graph algorithms. While a few such algorithms are known for spanning trees, matchings, and single-source shortest paths, very little was known for an important primitive like graph sparsifiers. The challenge is how to approximately preserve so much information about the graph (e.g., all-pairs distances and all cuts) without revealing the algorithms' underlying randomness to the adaptive adversary. In this paper we present the first non-trivial efficient adaptive algorithms for maintaining spanners and cut sparisifers. These algorithms in turn imply improvements over existing algorithms for other problems. Our first algorithm maintains a polylog$(n)$-spanner of size $\tilde O(n)$ in polylog$(n)$ amortized update time. The second algorithm maintains an $O(k)$-approximate cut sparsifier of size $\tilde O(n)$ in $\tilde O(n^{1/k})$ amortized update time, for any $k\ge1$, which is polylog$(n)$ time when $k=\log(n)$. The third algorithm maintains a polylog$(n)$-approximate spectral sparsifier in polylog$(n)$ amortized update time. The amortized update time of both algorithms can be made worst-case by paying some sub-polynomial factors. Prior to our result, there were near-optimal algorithms against oblivious adversaries (e.g. Baswana et al. [TALG'12] and Abraham et al. [FOCS'16]), but the only non-trivial adaptive dynamic algorithm requires $O(n)$ amortized update time to maintain $3$- and $5$-spanner of size $O(n^{1+1/2})$ and $O(n^{1+1/3})$, respectively [Ausiello et al. ESA'05]. Our results are based on two novel techniques. The first technique, is a generic black-box reduction that allows us to assume that the graph undergoes only edge deletions and, more importantly, remains an expander with almost-uniform degree. The second technique we call proactive resampling. [...]


title: A Weighted Population Update Rule for PACO Applied to the Single Machine Total Weighted Tardiness Problem

abstract: In this paper a new population update rule for population based ant colony optimization (PACO) is proposed. PACO is a well known alternative to the standard ant colony optimization algorithm. The new update rule allows to weight different parts of the solutions. PACO with the new update rule is evaluated for the example of the single machine total weighted tardiness problem (SMTWTP). This is an $\mathcal{NP}$-hard optimization problem where the aim is to schedule jobs on a single machine such that their total weighted tardiness is minimized. PACO with the new population update rule is evaluated with several benchmark instances from the OR-Library. Moreover, the impact of the weights of the jobs on the solutions in the population and on the convergence of the algorithm are analyzed experimentally. The results show that PACO with the new update rule has on average better solution quality than PACO with the standard update rule.


title: Projection-Cost-Preserving Sketches: Proof Strategies and Constructions

abstract: In this note we illustrate how common matrix approximation methods, such as random projection and random sampling, yield projection-cost-preserving sketches, as introduced in [FSS13, CEM+15]. A projection-cost-preserving sketch is a matrix approximation which, for a given parameter $k$, approximately preserves the distance of the target matrix to all $k$-dimensional subspaces. Such sketches have applications to scalable algorithms for linear algebra, data science, and machine learning. Our goal is to simplify the presentation of proof techniques introduced in [CEM+15] and [CMM17] so that they can serve as a guide for future work. We also refer the reader to [CYD19], which gives a similar simplified exposition of the proof covered in Section 2.


title: Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

abstract: We study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu $L$-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped $t$-$t^\prime$-$J$ model for a set of $t^\prime$ values. We find evidence for the breakdown of spin-charge separation as $|t^\prime|$ increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.


title: Analyzing the discrepancy principle for kernelized spectral filter learning algorithms

abstract: We investigate the construction of early stopping rules in the nonparametric regression problem where iterative learning algorithms are used and the optimal iteration number is unknown. More precisely, we study the discrepancy principle, as well as modifications based on smoothed residuals, for kernelized spectral filter learning algorithms including gradient descent. Our main theoretical bounds are oracle inequalities established for the empirical estimation error (fixed design), and for the prediction error (random design). From these finite-sample bounds it follows that the classical discrepancy principle is statistically adaptive for slow rates occurring in the hard learning scenario, while the smoothed discrepancy principles are adaptive over ranges of faster rates (resp. higher smoothness parameters). Our approach relies on deviation inequalities for the stopping rules in the fixed design setting, combined with change-of-norm arguments to deal with the random design setting.


title: Symmetry of relative equilibria with one dominant and four infinitesimal point vortices

abstract: We investigate the symmetry of point vortices with one dominant vortex and four vortices with infinitesimal circulations in the (1+4)-vortex problem, a subcase of the five-vortex problem. The four infinitesimal vortices inscribe quadrilaterals in the unit circle with the dominant vortex at the origin. We consider symmetric configurations which have one degree of spacial freedom, namely the (1+N)-gon, kites, rectangles, and trapezoids with three equal sides. We show there is only one possible rectangular configuration (up to rotation and ordering of the vortices) and one possible trapezoid with three equal sides (up to rotation and ordering), while there are parametrically defined families of kites. Additionally we consider the (1+4)-gon and show that the infinitesimal vortices must have equal circulations on opposite corners of the square. The proofs are heavily dependent on techniques from algebraic geometry and require the use of a computer to calculate Grobner bases.


title: Absorption by deformed black holes

abstract: Alternative theories of gravity and the parameterized deviation approach allow black hole solutions to have additional parameters beyond mass, charge and angular momentum. Matter fields could be, in principle, affected by the additional parameters of these solutions. We compute the absorption cross section of massless spin-0 waves by static Konoplya-Zhidenko black holes, characterized by a deformation parameter introduced in the mass term, and compare it with the well-known absorption of a Schwarzschild black hole with the same mass. We compare our numerical results with the sinc approximation in the high-frequency limit, finding excellent agreement.


title: Scaling the training of particle classification on simulated MicroBooNE events to multiple GPUs

abstract: Measurements in Liquid Argon Time Projection Chamber (LArTPC) neutrino detectors, such as the MicroBooNE detector at Fermilab, feature large, high fidelity event images. Deep learning techniques have been extremely successful in classification tasks of photographs, but their application to LArTPC event images is challenging, due to the large size of the events. Events in these detectors are typically two orders of magnitude larger than images found in classical challenges, like recognition of handwritten digits contained in the MNIST database or object recognition in the ImageNet database. Ideally, training would occur on many instances of the entire event data, instead of many instances of cropped regions of interest from the event data. However, such efforts lead to extremely long training cycles, which slow down the exploration of new network architectures and hyperparameter scans to improve the classification performance. We present studies of scaling a LArTPC classification problem on multiple architectures, spanning multiple nodes. The studies are carried out on simulated events in the MicroBooNE detector. We emphasize that it is beyond the scope of this study to optimize networks or extract the physics from any results here. Institutional computing at Pacific Northwest National Laboratory and the SummitDev machine at Oak Ridge National Laboratory's Leadership Computing Facility have been used. To our knowledge, this is the first use of state-of-the-art Convolutional Neural Networks for particle physics and their attendant compute techniques onto the DOE Leadership Class Facilities. We expect benefits to accrue particularly to the Deep Underground Neutrino Experiment (DUNE) LArTPC program, the flagship US High Energy Physics (HEP) program for the coming decades.

categories: physics.comp-ph cs.LG physics.data-an stat.ML


title: Parallelization Techniques for Verifying Neural Networks

abstract: Inspired by recent successes with parallel optimization techniques for solving Boolean satisfiability, we investigate a set of strategies and heuristics that aim to leverage parallel computing to improve the scalability of neural network verification. We introduce an algorithm based on partitioning the verification problem in an iterative manner and explore two partitioning strategies, that work by partitioning the input space or by case splitting on the phases of the neuron activations, respectively. We also introduce a highly parallelizable pre-processing algorithm that uses the neuron activation phases to simplify the neural network verification problems. An extensive experimental evaluation shows the benefit of these techniques on both existing benchmarks and new benchmarks from the aviation domain. A preliminary experiment with ultra-scaling our algorithm using a large distributed cloud-based platform also shows promising results.


title: Structuring Communities for Sharing Human Digital Memories in a Social P2P Network

abstract: A community is sub-network inside P2P networks that partition the network into groups of similar peers to improve performance by reducing network traffic and high search query success rate. Large communities are common in online social networks than traditional file-sharing P2P networks because many people capture huge amounts of data through their lives. This increases the number of hosts bearing similar data in the network and hence increases the size of communities. This article presents a Memory Thread-based Communities for our entity-based social P2P network that partition the network into groups of peers sharing data belonging to an entity - person, place, object or interest, having its own digital memory or be a part another memory. These connected peers having further similarities by organizing the network using linear orderings. A Memory-Thread is the collection of digital memories having a common reference key and organized according to some form of correlation. The simulation results show an increase in network performance for the proposed scheme along with a decrease in network overhead and higher query success rate compared to other similar schemes. The network maintains its performance even while the network traffic and size increase.


title: Spring Festival points the way to cleaner air in China

abstract: Human migration during the Chinese Spring Festival (SF) is the largest collective human activity of its kind in the modern era-involving about one-tenth of the world population and over six percent of the earth's land surface area. The festival results in a drop of air pollutant emissions that causes dramatic changes of atmospheric composition over China's most polluted regions. Based on satellite and in-situ measurements for the years 2005-2019 over 50 cities in eastern China, we find that the atmospheric NO2 pollution dropped by ~40% during the SF week, and fine particulate matter (PM2.5) decreased by ~30% in the following week, reflecting the effectiveness of precursor emission controls on the mitigation of secondary PM2.5 formation. However, although human activity and emissions are at the lowest level, air pollution over eastern China during the SF still far exceeds that over other worldwide pollution hotspots. Our analyses suggest that measures based solely on end-of-pipe controls and industry upgrades may not suffice to meet air quality goals. Further cleaning of the air in China depends fundamentally on sustainable advances in both heavy industry upgrades and clean energy transition.


title: Adversarial Attack on Deep Learning-Based Splice Localization

abstract: Regarding image forensics, researchers have proposed various approaches to detect and/or localize manipulations, such as splices. Recent best performing image-forensics algorithms greatly benefit from the application of deep learning, but such tools can be vulnerable to adversarial attacks. Due to the fact that most of the proposed adversarial example generation techniques can be used only on end-to-end classifiers, the adversarial robustness of image-forensics methods that utilize deep learning only for feature extraction has not been studied yet. Using a novel algorithm capable of directly adjusting the underlying representations of patches we demonstrate on three non end-to-end deep learning-based splice localization tools that hiding manipulations of images is feasible via adversarial attacks. While the tested image-forensics methods, EXIF-SC, SpliceRadar, and Noiseprint, rely on feature extractors that were trained on different surrogate tasks, we find that the formed adversarial perturbations can be transferable among them regarding the deterioration of their localization performance.


title: On the Approximate Nearest Neighbor Queries among Curves under the Fr\'echet Distance

abstract: Approximate near-neighbors search (\textsc{ANNS}) is a long-studied problem in computational geometry. %that has received considerable attention by researchers in the community. In this paper, we revisit the problem and propose the first data structure for curves under the (continuous) Fr\'echet distance in $\Reals^d$. Given a set $\P$ of $n$ curves of size at most $m$ each in $\Reals^d$, and a real fixed $\delta>0$, we aim to preprocess $\P$ into a data structure so that for any given query curve $Q$ of size $k$, we can efficiently report all curves in $\P$ whose Fr\'echet distances to $Q$ are at most $\delta$. In the case that $k$ is given in the preprocessing stage, for any $\eps>0$ we propose a deterministic data structure whose space is $n \cdot O\big(\max\big\{\big(\frac{\sqrt{d}}{\eps}\big)^{kd}, \big(\frac{\D\sqrt{d}}{\eps^2}\big)^{kd}\big\}\big)$ that can answer \textsc{$(1+\eps)\delta$-ANNS} queries in $O(kd)$ query time, where $\D$ is the diameter of $\P$. Considering $k$ as part of the query slightly changes the space to $n \cdot O\big(\frac{1}{\eps}\big)^{md} $ with $O(kd)$ query time within an approximation factor of $5+\eps$. We show that our generic data structure for ANNS can give an alternative treatment of the approximate subtrajectory range searching problem studied by de Berg et al. [8]. We also revisit the time-window data structure for spatial density maps in [6]. Given $\theta>0$, and $n$ time-stamped points spread over $m$ regions in a map, for any query window $W$, we propose a data structure of size $O(n/\eps^2)$ and construction time $O((n+m)/\eps^2)$ that can approximately return the regions containing at least $\theta$ points whose times are within $W$ in $O(1)$ query time.


title: Signatures of multiple charge excitations in RIXS spectra of metals

abstract: We study how multiple charge excitations appear in the resonant inelastic x-ray scattering (RIXS) spectra of metals. The single excitations in the problem are the plasmons and electron-hole pairs, and multi-excitation processes are usually neglected. However, at small momentum transfer the multi-excitation contributions may dominate the signal and one needs to understand how to interpret the data. In particular, we demonstrate how to "decode" the total multi-excitation intensity and extract the plasmon dispersion. While our calculations are based on the random phase approximation, which does not allow to obtain quantitatively precise results in the entire region of parameters, we expect them to capture semi-qualitatively all features expected for charged Fermi-liquid states, including universal and singular properties of the RIXS spectra.


title: Rational cubic fourfolds with associated singular K3 surfaces

abstract: Generalizing a recent construction of Yang and Yu, we study to what extent one can intersect Hassett's Noether-Lefschetz divisors $\mathcal{C}_d$ in the moduli space of cubic fourfolds $\mathcal{C}$. In particular, we exhibit arithmetic conditions on 20 indexes $d_1,\dots, d_{20}$ that assure that the divisors $\mathcal{C}_{d_1},\dots,\mathcal{C}_{d_{20}}$ all intersect one another. This allows us to produce examples of rational cubic fourfolds with an associated K3 surface with rank 20 N\'eron-Severi group, i.e. a singular K3 surface.


title: Response of a plate in turbulent channel flow: Analysis of fluid-solid coupling

abstract: The paper performs simulation of a rectangular plate excited by turbulent channel flow at friction Reynolds numbers of 180 and 400. The fluid-structure interaction is assumed to be one-way coupled, i.e, the fluid affects the solid and not vice versa. We solve the incompressible Navier Stokes equations using finite volume direct numerical simulation in the fluid domain. In the solid domain, we solve the dynamic linear elasticity equations using a time-domain finite element method. The obtained plate averaged displacement spectra collapse in the low frequency region in outer scaling. However, the high frequency spectral levels do not collapse in inner units. This spectral behavior is reasoned using theoretical arguments. We further study the sources of plate excitation using a novel formulation. This formulation expresses the average displacement spectra of the plate as an integrated contribution from the fluid sources within the channel. Analysis of the net displacement source reveals that at the plate natural frequencies, the contribution of the fluid sources to the plate excitation peaks in the buffer layer. The corresponding wall-normal width is found to be $\approx 0.75\delta$. We analyze the decorrelated features of the sources using spectral Proper Orthogonal Decomposition (POD) of the net displacement source. We enforce the orthogonality of the modes in an inner product with a symmetric positive definite kernel. The dominant spectral POD mode contributes to the entire plate excitation. The contribution of the remaining modes from the different wall-normal regions undergo destructive interference resulting in zero net contribution. The envelope of the dominant mode further shows that the location and width of the contribution depend on inner and outer units, respectively.


title: Bourgain-Brezis-Mironescu approach in metric spaces with Euclidean tangents

abstract: In the setting of metric measure spaces satisfying the doubling condition and the $(1,p)$-Poincar\'e inequality, we prove a metric analogue of the Bourgain-Brezis-Mironescu formula for functions in the Sobolev space $W^{1,p}(X,d,\nu)$, under the assumption that for $\nu$-a.e. point the tangent space in the Gromov-Hausdorff sense is Euclidean with fixed dimension $N$.


title: Can You Put it All Together: Evaluating Conversational Agents' Ability to Blend Skills

abstract: Being engaging, knowledgeable, and empathetic are all desirable general qualities in a conversational agent. Previous work has introduced tasks and datasets that aim to help agents to learn those qualities in isolation and gauge how well they can express them. But rather than being specialized in one single quality, a good open-domain conversational agent should be able to seamlessly blend them all into one cohesive conversational flow. In this work, we investigate several ways to combine models trained towards isolated capabilities, ranging from simple model aggregation schemes that require minimal additional training, to various forms of multi-task training that encompass several skills at all training stages. We further propose a new dataset, BlendedSkillTalk, to analyze how these capabilities would mesh together in a natural conversation, and compare the performance of different architectures and training schemes. Our experiments show that multi-tasking over several tasks that focus on particular capabilities results in better blended conversation performance compared to models trained on a single skill, and that both unified or two-stage approaches perform well if they are constructed to avoid unwanted bias in skill selection or are fine-tuned on our new task.


title: Reducing Commutativity Verification to Reachability with Differencing Abstractions

abstract: Commutativity of data structure methods is of ongoing interest, with roots in the database community. In recent years commutativity has been shown to be a key ingredient to enabling multicore concurrency in contexts such as parallelizing compilers, transactional memory, speculative execution and, more broadly, software scalability. Despite this interest, it remains an open question as to how a data structure's commutativity specification can be verified automatically from its implementation. In this paper, we describe techniques to automatically prove the correctness of method commutativity conditions from data structure implementations. We introduce a new kind of abstraction that characterizes the ways in which the effects of two methods differ depending on the order in which the methods are applied, and abstracts away effects of methods that would be the same regardless of the order. We then describe a novel algorithm that reduces the problem to reachability, so that off-the-shelf program analysis tools can perform the reasoning necessary for proving commutativity. Finally, we describe a proof-of-concept implementation and experimental results, showing that our tool can verify commutativity of data structures such as a memory cell, counter, two-place Set, array-based stack, queue, and a rudimentary hash table. We conclude with a discussion of what makes a data structure's commutativity provable with today's tools and what needs to be done to prove more in the future.


title: An efficient algorithm for graph Laplacian optimization based on effective resistances

abstract: In graph signal processing, data samples are associated to vertices on a graph, while edge weights represent similarities between those samples. We propose a convex optimization problem to learn sparse well connected graphs from data. We prove that each edge weight in our solution is upper bounded by the inverse of the distance between data features of the corresponding nodes. We also show that the effective resistance distance between nodes is upper bounded by the distance between nodal data features. Thus, our proposed method learns a sparse well connected graph that encodes geometric properties of the data. We also propose a coordinate minimization algorithm that, at each iteration, updates an edge weight using exact minimization. The algorithm has a simple and low complexity implementation based on closed form expressions.


title: A Subterranean Virtual Cave World for Gazebo based on the DARPA SubT Challenge

abstract: Subterranean environments with lots of obstacles, including narrow passages, large voids, rock falls and absence of illumination were always challenging for control, navigation, and perception of mobile robots. The limited availability and access to such environments restricts the development pace of capabilities for robotic platforms to autonomously accomplish tasks in such challenging areas. The Subterranean Challenge is a competition focusing on bringing robotic exploration a step closer to real life applications for man-made underground tunnels, urban areas and natural cave networks, envisioning advanced assistance tools for first responders and disaster relief agencies. The challenge offers a software-based virtual part to showcase technologies in autonomy perception, networking and mobility for such areas. Thus, the presented open-source virtual world aims to become a test-bed for evaluating the developed algorithms and software and to foster mobile robotics developments.


title: A micromagnetic study of the switching dynamics of the BiFeO$_3$/CoFe heterojunction

abstract: The switching dynamics of a single-domain BiFeO3/CoFe heterojunction is modeled and key parameters such as interface exchange coupling coefficient are extracted from experimental results. The lower limit of the magnetic order response time of CoFe in the BiFeO3/CoFe heterojunction is theoretically quantified to be on to the order of 100 ps. Our results indicate that the switching behavior of CoFe in the BiFeO3/CoFe heterojunction is dominated by the rotation of the Neel vector in BiFeO3 rather than the unidirectional exchange bias at the interface. We also quantify the magnitude of the interface exchange coupling coefficient J_int to be 0.32 pJ/m by comparing our simulation results with the giant magnetoresistance (GMR) curves and the magnetic hysteresis loop in the experiments. To the best of our knowledge, this is the first time that J_int is extracted quantitatively from experiments. Furthermore, we demonstrate that the switching success rate and the thermal stability of the BiFeO3/CoFe heterojunction can be improved by reducing the thickness of CoFe and increasing the length to width aspect ratio of the BiFeO3/CoFe heterojunction. Our theoretical model provides a comprehensive framework to study the magnetoelectric properties and the manipulation of the magnetic order of CoFe in the BiFeO3/CoFe heterojunction.


title: Counterexamples to the Low-Degree Conjecture

abstract: A conjecture of Hopkins (2018) posits that for certain high-dimensional hypothesis testing problems, no polynomial-time algorithm can outperform so-called "simple statistics", which are low-degree polynomials in the data. This conjecture formalizes the beliefs surrounding a line of recent work that seeks to understand statistical-versus-computational tradeoffs via the low-degree likelihood ratio. In this work, we refute the conjecture of Hopkins. However, our counterexample crucially exploits the specifics of the noise operator used in the conjecture, and we point out a simple way to modify the conjecture to rule out our counterexample. We also give an example illustrating that (even after the above modification), the symmetry assumption in the conjecture is necessary. These results do not undermine the low-degree framework for computational lower bounds, but rather aim to better understand what class of problems it is applicable to.


title: Multi wavelength cross-correlation analysis of the simulated cosmic web

abstract: We used magneto-hydrodynamical cosmological simulations to investigate the cross-correlation between different observables (i.e. X-ray emission, Sunyaev-Zeldovich signal at 21 cm, HI temperature decrement, diffuse synchrotron emission and Faraday Rotation) as a probe of the diffuse matter distribution in the cosmic web. We adopt an uniform and simplistic approach to produce synthetic observations at various wavelengths, and we compare the detection chances of different combinations of observables correlated with each other and with the underlying galaxy distribution in the volume. With presently available surveys of galaxies and existing instruments, the best chances to detect the diffuse gas in the cosmic web outside of halos is by cross-correlating the distribution of galaxies with Sunyaev-Zeldovich observations. We also find that the cross-correlation between the galaxy network and the radio emission or the Faraday Rotation can already be used to limit the amplitude of extragalactic magnetic fields, well outside of the cluster volume usually explored by existing radio observations, and to probe the origin of cosmic magnetism with the future generation of radio surveys.


title: Narrowband and passband composite pulses for variable rotations

abstract: We develop a systematic approach to derive narrowband (NB) and passband (PB) composite sequences which can produce any pre-selected transition probability with any desired accuracy. The NB composite pulses are derived by successive cancellation of derivative terms in the propagator. The PB pulses are built by combining two half-$\pi$ NB pulses. Both the NB and PB pulses feature vanishing wings on both sides of the central region around the target transition probability, and the PB pulses offer also a flat (broadband) top around the target probability. The NB sequences are particularly suitable for sensing and metrology applications, while the PB sequences can be very useful for qubit control in the case of tightly spaced qubits for they ensure both selectivity (suppressed cross talk to neighbors) and robustness (suppressed laser pointing instability). Other possible applications include selective (NB), or both selective and robust (PB), control of close transitions in frequency space.


title: Long-range miniature ceramic RFID tags

abstract: Radio frequency identification (RFID) is a mature technology, which allows performing contactless data readout via wireless communication links. While communication protocols in this field are subject to regulations, there is a room of opportunities to improve the hardware realization of antennas devices, which support the technology. In particular, readout range extension and miniaturization of passive RFID tags is an important objective with far going impact on retail, security, IoT, and many others. Here we introduce a new concept of high-permittivity ceramic tag, which relies on different physical principles. Instead of using conduction currents in metallic wires to drive electronic chips and radiate electromagnetic waves, high permittivity components rely on an efficient excitation of displacement currents. Those are efficiently converted to conduction currents, powering a memory chip. The practical aspect of this approach is improved robustness to environmental fluctuations, footprint reduction, and readout range extension. In particular, our high permittivity ceramic ({\epsilon} ~ 100) elements have demonstrated a 25% reading range improvement in respect to commercial tags. In case, when state of the art readers and RFID chips are used, the readout distances of the developed ceramic tags are approaching 23 m and could be further extended with improved matching circuits.


title: Multiplicity for a Group Sequential Trial with Biomarker Subpopulations

abstract: Biomarker subpopulations have become increasingly important for drug development in targeted therapies. The use of biomarkers has the potential to facilitate more effective outcomes by guiding patient selection appropriately, thus enhancing the benefit-risk profile and improving trial power. Studying a broad population simultaneously with a more targeted one allows the trial to determine the population for which a treatment is effective and allows a goal of making approved regulatory labeling as inclusive as is appropriate. We examine new methods accounting for the complete correlation structure in group sequential designs with hypotheses in nested subgroups. The designs provide full control of family-wise Type I error rate. This extension of previous methods accounting for either group sequential design or correlation between subgroups improves efficiency (power or sample size) over a typical Bonferroni approach for testing nested populations.


title: Dynamical $\ell$-boson stars: generic stability and evidence for non-spherical solutions

abstract: $\ell$-boson stars are static, spherical, multi-field self-gravitating solitons. They are asymptotically flat, finite energy solutions of Einstein's gravity minimally coupled to an odd number of massive, complex scalar fields. A previous study assessed the stability of $\ell$-boson stars under spherical perturbations, finding that there are both stable and unstable branches of solutions, as for single-field boson stars ($\ell=0$). In this work we probe the stability of $\ell$-boson stars against non-spherical perturbations by performing numerical evolutions of the Einstein-Klein-Gordon system, with a 3D code. For the timescales explored, the $\ell$-boson stars belonging to the spherical stable branch do not exhibit measurable growing modes. We find, however, evidence of zero modes; that is, non-spherical perturbations that neither grow nor decay. This suggests the branching off towards a larger family of equilibrium solutions: we conjecture that $\ell$-boson stars are the enhanced isometry point of a larger family of static (and possibly stationary), non-spherical multi-field self-gravitating solitons.


title: Estimating and Projecting Air Passenger Traffic during the COVID-19 Coronavirus Outbreak and its Socio-Economic Impact

abstract: The main focus of this study is to collect and prepare data on air passengers traffic worldwide with the scope of analyze the impact of travel ban on the aviation sector. Based on historical data from January 2010 till October 2019, a forecasting model is implemented in order to set a reference baseline. Making use of airplane movements extracted from online flight tracking platforms and on-line booking systems, this study presents also a first assessment of recent changes in flight activity around the world as a result of the COVID-19 pandemic. To study the effects of air travel ban on aviation and in turn its socio-economic, several scenarios are constructed based on past pandemic crisis and the observed flight volumes. It turns out that, according to this hypothetical scenarios, in the first Quarter of 2020 the impact of aviation losses could have negatively reduced World GDP by 0.02% to 0.12% according to the observed data and, in the worst case scenarios, at the end of 2020 the loss could be as high as 1.41-1.67% and job losses may reach the value of 25-30 millions. Focusing on EU27, the GDP loss may amount to 1.66-1.98% by the end of 2020 and the number of job losses from 4.2 to 5 millions in the worst case scenarios. Some countries will be more affected than others in the short run and most European airlines companies will suffer from the travel ban.


title: Algebraic Relations among Goss's Zeta Values on Elliptic Curves

abstract: In 2007 Chang and Yu determined all the algebraic relations among Goss's zeta values for the rational function field - these are also known as the Carlitz zeta values. Goss raised the problem about algebraic relations among Goss's zeta values for a general base ring A, but very little is known. In this paper we develop a general method and determine all algebraic relations among Goss's zeta values for the base ring A attached to an elliptic curve over a finite field. To our knowledge, these are the first non-trivial solutions of Goss's problem for a base ring whose class number is strictly greater than 1.


title: Stability of radical-functionalized gold surfaces by self-assembly and on-surface chemistry

abstract: We have investigated the radical functionalization of gold surfaces with a derivative of the perchlorotriphenylmethyl (PTM) radical, using two methods: by chemisorption from the radical solution and by on surface chemical derivatization from a precursor. We have investigated the obtained self-assembled monolayers by photon-energy dependent X-ray photoelectron spectroscopy. Our results show that the molecules were successfully anchored on the surfaces. The monolayers are characterized by air and beam stability unprecedented for films of organic radicals. Over very long beam exposure we observed a dynamic nature of the radical-Au complex. The results clearly indicate that (mono)layers of PTM; radical derivatives have the necessary stability to stand device applications.


title: A Unified Calibration Framework for 21 cm Cosmology

abstract: Calibration precision is currently a limiting systematic in 21 cm cosmology experiments. While there are innumerable calibration approaches, most can be categorized as either `sky-based,' relying on an extremely accurate model of astronomical foreground emission, or `redundant,' requiring a precisely regular array with near-identical antenna response patterns. Both of these classes of calibration are inflexible to the realities of interferometric measurement. In practice, errors in the foreground model, antenna position offsets, and beam response inhomogeneities degrade calibration performance and contaminate the cosmological signal. Here we show that sky-based and redundant calibration can be unified into a highly general and physically motivated calibration framework based on a Bayesian statistical formalism. Our new framework includes sky and redundant calibration as special cases but can additionally support relaxing the rigid assumptions implicit in those approaches. Furthermore, we present novel calibration techniques such as redundant calibration for arrays with no redundant baselines, representing an alternative calibration method for imaging arrays such as the MWA Phase I. These new calibration approaches could mitigate systematics and reduce calibration error, thereby improving the precision of cosmological measurements.


title: Global Analysis of the Gravitational Wave Signal from Galactic Binaries

abstract: Galactic ultra compact binaries are expected to be the dominant source of gravitational waves in the milli-Hertz frequency band. Of the tens of millions of galactic binaries with periods shorter than an hour, it is estimated that a few tens of thousand will be resolved by the future Laser Interferometer Space Antenna (LISA). The unresolved remainder will be the main source of ``noise'' between 1-3 milli-Hertz. Typical galactic binaries are millions of years from merger, and consequently their signals will persist for the the duration of the LISA mission. Extracting tens of thousands of overlapping galactic signals and characterizing the unresolved component is a central challenge in LISA data analysis, and a key contribution to arriving at a global solution that simultaneously fits for all signals in the band. Here we present an end-to-end analysis pipeline for galactic binaries that uses trans-dimensional Bayesian inference to develop a time-evolving catalog of sources as data arrive from the LISA constellation.


title: Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

abstract: We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. We analyze the effect of the system size on the ice Ih-liquid free energy differences and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.


title: Application of Progressive Hedging to Var Expansion Planning Under Uncertainty

abstract: This paper describes the application of a Progressive Hedging (PH) algorithm to the least-cost var planning under uncertainty. The method PH is a scenario-based decomposition technique for solving stochastic programs, i.e., it decomposes a large scale stochastic problem into s deterministic subproblems and couples the decision from the s subproblems to form a solution for the original stochastic problem. The effectiveness and computational performance of the proposed methodology will be illustrated with var planning studies for the IEEE 24-bus system (5 operating scenarios), the 200-bus Bolivian system (1,152 operating scenarios) and the 1,600-bus Colombian system (180 scenarios).


title: Lidar for Autonomous Driving: The principles, challenges, and trends for automotive lidar and perception systems

abstract: Autonomous vehicles rely on their perception systems to acquire information about their immediate surroundings. It is necessary to detect the presence of other vehicles, pedestrians and other relevant entities. Safety concerns and the need for accurate estimations have led to the introduction of Light Detection and Ranging (LiDAR) systems in complement to the camera or radar-based perception systems. This article presents a review of state-of-the-art automotive LiDAR technologies and the perception algorithms used with those technologies. LiDAR systems are introduced first by analyzing the main components, from laser transmitter to its beam scanning mechanism. Advantages/disadvantages and the current status of various solutions are introduced and compared. Then, the specific perception pipeline for LiDAR data processing, from an autonomous vehicle perspective is detailed. The model-driven approaches and the emerging deep learning solutions are reviewed. Finally, we provide an overview of the limitations, challenges and trends for automotive LiDARs and perception systems.


title: Loss aversion and the welfare ranking of policy interventions

abstract: This paper develops theoretical criteria and econometric methods to rank policy interventions in terms of welfare when individuals are loss-averse. Our new criterion for "loss aversion-sensitive dominance" defines a weak partial ordering of the distributions of policy-induced gains and losses. It applies to the class of welfare functions which model individual preferences with non-decreasing and loss-averse attitudes towards changes in outcomes. We also develop new statistical methods to test loss aversion-sensitive dominance in practice, using nonparametric plug-in estimates; these allow inference to be conducted through a special resampling procedure. Since point-identification of the distribution of policy-induced gains and losses may require strong assumptions, we extend our comparison criteria, test statistics, and resampling procedures to the partially-identified case. We illustrate our methods with a simple empirical application to the welfare comparison of alternative income support programs in the US.


title: Study of Four-Dimensional DOA and Polarisation Estimation with Crossed-dipole and Tripole Arrays

abstract: Electromagnetic (EM) vector sensor arrays can track both the polarisation and direction of arrival (DOA) of the impinging signals. For linear crossed-dipole arrays, as shown by our analysis, due to inherent limitation of the structure, it can only track one DOA parameter and two polarisation parameters. For full four-dimensional (4-D, 2 DOA and 2 polarization parameters) estimation, we could extend the linear crossed-dipole array to the planar case. In this paper, instead of extending the array geometry, we replace the crossed-dipoles by tripoles and construct a linear tripole array. It is proved that such a structure can estimate the 2-D DOA and 2-D polarisation information effectively in general and a dimension-reduction based MUSIC algorithm is developed so that the 4-D estimation problem can be simplified into two separate 2-D estimation problems, significantly reducing the computational complexity of the solution. The Cramr-Rao Bound (CRB) is also derived as a reference for algorithm performance. A brief comparison between the planar crossed-dipole array and the linear tripole array is performed at last, showing that although the planar structure has a better performance, it is achieved at the cost of increased physical size.


title: Growth and Characterization of Homoepitaxial $\beta$-Ga$_2$O$_3$ Layers

abstract: $\beta$-Ga$_2$O$_3$ is a next-generation ultra wide bandgap semiconductor (E$_g$ = 4.8 eV to 4.9 eV) that can be homoepitaxially grown on commercial substrates, enabling next-generation power electronic devices among other important applications. Analyzing the quality of deposited homoepitaxial layers used in such devices is challenging, in part due to the large probing depth in traditional x-ray diffraction (XRD) and also due to the surface-sensitive nature of atomic force microscopy (AFM). Here, a combination of evanescent grazing-incidence skew asymmetric XRD and AFM are investigated as an approach to effectively characterize the quality of homoepitaxial $\beta$-Ga$_2$O$_3$ layers grown by molecular beam epitaxy at a variety of Ga/O flux ratios. Accounting for both structure and morphology, optimal films are achieved at a Ga/O ratio of $\sim$1.15, a conclusion that would not be possible to achieve by either XRD or AFM methods alone. Finally, fabricated Schottky barrier diodes with thicker homoepitaxial layers are characterized by $J-V$ and $C-V$ measurements, revealing an unintentional doping density of 4.3 $\times$ 10$^{16}$ cm$^{-3}$ - 2 $\times$ 10$^{17}$ cm$^{-3}$ in the epilayer. These results demonstrate the importance of complementary measurement methods for improving the quality of the $\beta$-Ga$_2$O$_3$ homoepitaxial layers used in power electronic and other devices.


title: Model-driven Engineering of Safety and Security Systems: A Systematic Mapping Study

abstract: This paper presents a systematic mapping study on the model-driven engineering of safety and security concerns in systems. Integrated modeling and development of both safety and security concerns is an emerging field of research. Our mapping study provides an overview of the current state-of-the-art in this field. Through a rigorous and systematic process, this study carefully selected 95 publications out of 17,927 relevant papers published between 1992 and 2018. This paper then proposes and answers several relevant research questions about frequently used methods, development stages where these concerns are typically investigated in, or application domains. Additionally, we identify the community's preference for publication venues and trends.


title: Leveraging the Fisher randomization test using confidence distributions: inference, combination and fusion learning

abstract: The flexibility and wide applicability of the Fisher randomization test (FRT) makes it an attractive tool for assessment of causal effects of interventions from modern-day randomized experiments that are increasing in size and complexity. This paper provides a theoretical inferential framework for FRT by establishing its connection with confidence distributions Such a connection leads to development of (i) an unambiguous procedure for inversion of FRTs to generate confidence intervals with guaranteed coverage, (ii) generic and specific methods to combine FRTs from multiple independent experiments with theoretical guarantees and (iii) new insights on the effect of size of the Monte Carlo sample on the results of FRT. Our developments pertain to finite sample settings but have direct extensions to large samples. Simulations and a case example demonstrate the benefit of these new developments.


title: Fork-Resilient Cross-Blockchain Transactions through Algebraic Topology

abstract: The cross-blockchain transaction (CBT) serves as a cornerstone for the next-generation, blockchain-based data management systems. However, state-of-the-art CBT models do not address the effect of the possible local fork suspension that might invalidate the entire CBT. This paper takes an algebraic-topological approach to abstract the blockchains and their transactions into simplicial complexes and shows that CBTs cannot complete in either a \textit{committed} or an \textit{aborted} status by a $t$-resilient message-passing protocol. This result implies that a more sophisticated model is in need to support CBTs and, thus, sheds light on the future blockchain designs.


title: Quantifying uncertainties in neutron-alpha scattering with chiral nucleon-nucleon and three-nucleon forces

abstract: Background: Modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Purpose: Starting from NN potentials up to fifth order (N4LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (cE and cD) contribute to the overall uncertainty budget of many-body calculations of neutron-He elastic scattering. Methods: We compute structure and reaction observables for three-, four- and five-nucleon systems within the ab initio frameworks of the no-core shell model an no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte-Carlo.


title: A Simple, General, and GPU Friendly Method for Computing Dual Mesh and Iso-Surfaces of Adaptive Mesh Refinement (AMR) Data

abstract: We propose a novel approach to extracting crack-free iso-surfaces from Structured AMR data that is more general than previous techniques, is trivially simple to implement, requires no information other than the list of AMR cells, and works, in particular, for different AMR formats including octree AMR, block-structured AMR with arbitrary level differences at level boundaries, and AMR data that consist of individual cells without any existing grid structure. We describe both the technique itself and a CUDA-based GPU implementation of this technique, and evaluate it on several non-trivial AMR data sets.


title: Learning-to-Rank with BERT in TF-Ranking

abstract: This paper describes a machine learning algorithm for document (re)ranking, in which queries and documents are firstly encoded using BERT [1], and on top of that a learning-to-rank (LTR) model constructed with TF-Ranking (TFR) [2] is applied to further optimize the ranking performance. This approach is proved to be effective in a public MS MARCO benchmark [3]. Our first two submissions achieve the best performance for the passage re-ranking task [4], and the second best performance for the passage full-ranking task as of April 10, 2020 [5]. To leverage the lately development of pre-trained language models, we recently integrate RoBERTa [6] and ELECTRA [7]. Our latest submissions improve our previously state-of-the-art re-ranking performance by 4.3% [8], and achieve the third best performance for the full-ranking task [9] as of June 8, 2020. Both of them demonstrate the effectiveness of combining ranking losses with BERT representations for document ranking.


title: Self-Interacting Dark Matter and the Delay of Super-Massive Black Hole Growth

abstract: Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.


title: Automorphisms of shift spaces and the Higman--Thompson groups: the one-sided case

abstract: Let $1 \le r < n$ be integers. We give a proof that the group $\mathop{\mathrm{Aut}}({X_{n}^{\mathbb{N}}, \sigma_{n}})$ of automorphisms of the one-sided shift on $n$ letters embeds naturally as a subgroup $\mathcal{H}_{n}$ of the outer automorphism group $\mathop{\mathrm{Out}}({G_{n,r}})$ of the Higman-Thompson group $G_{n,r}$. From this, we can represent the elements of $\mathop{\mathrm{Aut}}({X_{n}^{\mathbb{N}}, \sigma_{n}})$ by finite state non-initial transducers admitting a very strong synchronizing condition. Let $H \in \mathcal{H}_{n}$ and write $|H|$ for the number of states of the minimal transducer representing $H$. We show that $H$ can be written as a product of at most $|H|$ torsion elements. This result strengthens a similar result of Boyle, Franks and Kitchens, where the decomposition involves more complex torsion elements and also does not support practical \textit{a priori} estimates of the length of the resulting product. We also explore the number of foldings of de Bruijn graphs and give a counting result for these for word length $2$ and alphabet size $n$. Finally, we offer new proofs of some known results about $\mathop{\mathrm{Aut}}({X_{n}^{\mathbb{N}}, \sigma_{n}})$.


title: Lessons on early structure formation from a mature galaxy cluster observed at cosmic noon

abstract: We demonstrate a new approach of indirectly constraining both early star and structure formation via mature galaxy clusters at cosmic noon (z~2), using the cluster XLSSC 122 as an example. With the standard Press Schechter formalism, we infer a rapid evolution of the star formation efficiency (the ratio of stellar to halo mass) from 10^-4 to 0.01 during z~20-13, based on the age distribution of stars in post-starburst galaxies of XLSSC 122, measured by HST photometry assuming no dust extinction. Here, we consider all low-mass haloes, including minihaloes, that host the first stars and galaxies (5*10^5 Msun < M_halo < 10^10 Msun). We also place new constraints on fuzzy dark matter models of m_a < 5*10^-21 eV/c^2 for the ultra-light boson mass, from the abundance of galaxies with star formation at z > 13 in XLSSC 122. Our exploratory results are consistent with existing constraints. More comprehensive results will be obtained if our approach is extended to a large sample of clusters or field post-starburst galaxies at cosmic noon, with improved modelling of halo and stellar populations.


title: The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

abstract: We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.


title: An optimal pointwise Morrey-Sobolev inequality

abstract: Let $\Omega$ be a bounded, smooth domain of $\mathbb{R}^{N},$ $N\geq1.$ For each $p>N$ we study the optimal function $s=s_{p}$ in the pointwise inequality \[ \left\vert v(x)\right\vert \leq s(x)\left\Vert \nabla v\right\Vert _{L^{p}(\Omega)},\quad\forall\,(x,v)\in\overline{\Omega}\times W_{0}% ^{1,p}(\Omega). \] We show that $s_{p}\in C_{0}^{0,1-(N/p)}(\overline{\Omega})$ and that $s_{p}$ converges pointwise to the distance function to the boundary, as $p\rightarrow\infty.$ Moreover, we prove that if $\Omega$ is convex, then $s_{p}$ is concave and has a unique maximum point.


title: Evolving Ultralight Scalars into Non-Linearity with Lagrangian Perturbation Theory

abstract: Many models of high energy physics suggest that the cosmological dark sector consists of not just one, but a spectrum of ultralight scalar particles with logarithmically distributed masses. To study the potential signatures of low concentrations of ultralight axion (also known as fuzzy) dark matter, we modify Lagrangian perturbation theory (LPT) by distinguishing between trajectories of different dark matter species. We further adapt LPT to include the effects of a quantum pressure, which is necessary to generate correct initial conditions for ultralight axion simulations. Based on LPT, our modified scheme is extremely efficient on large scales and it can be extended to an arbitrary number of particle species at very little computational cost. This allows for computation of self-consistent initial conditions in mixed dark matter models. Additionally, we find that shell-crossing is delayed for ultralight particles and that the deformation tensor extracted from LPT can be used to identify the range of redshifts and scales for which the Madelung formalism of fuzzy dark matter is a reliable approximation.


title: ETC: Encoding Long and Structured Inputs in Transformers

abstract: Transformer models have advanced the state of the art in many Natural Language Processing (NLP) tasks. In this paper, we present a new Transformer architecture, Extended Transformer Construction (ETC), that addresses two key challenges of standard Transformer architectures, namely scaling input length and encoding structured inputs. To scale attention to longer inputs, we introduce a novel global-local attention mechanism between global tokens and regular input tokens. We also show that combining global-local attention with relative position encodings and a Contrastive Predictive Coding (CPC) pre-training objective allows ETC to encode structured inputs. We achieve state-of-the-art results on four natural language datasets requiring long and/or structured inputs.


title: Probing variation of the fine-structure constant in runaway dilaton models using Strong Gravitational Lensing and Type Ia Supernovae

abstract: In order to probe a possible time variation of the fine-structure constant ($\alpha$), we propose a new method based on Strong Gravitational Lensing and Type Ia Supernovae observations. By considering a class of dilaton runaway models, where $\frac{\Delta \alpha}{\alpha}= - \gamma \ln{(1+z)}$, we obtain constraints on $\frac{\Delta \alpha}{\alpha}$ at the level $\gamma \approx 10^{-2}$ ($\gamma$ captures the physical properties of the model). Since the data set covers the redshift range $0.075 \leq z \leq 2.2649$, the constraints derived here provide independent bounds on a possible time variation of $\alpha$ at low, intermediate and high redshifts.


title: Two phases with different domain wall networks and a reentrant phase transition in bilayer graphene under strain

abstract: The analytical two-chain Frenkel-Kontorova model is used to describe domain wall networks in bilayer graphene upon biaxial stretching of one of the layers. We show that the commensurate-incommensurate phase transition leading to formation of a regular triangular domain wall network at the relative biaxial elongation of $3.0\cdot 10^{-3}$ is followed by the transition to another incommensurate phase with a striped network at the elongation of $3.7\cdot 10^{-3}$. The reentrant transition to the phase with a triangular domain wall network is predicted for the elongation $\sim10^{-2}$.


title: Critical exponent for semi-linear structurally damped wave equation of derivative type

abstract: Main purpose of this paper is to study the following semi-linear structurally damped wave equation with nonlinearity of derivative type: $$u_{tt}- \Delta u+ \mu(-\Delta)^{\sigma/2} u_t= |u_t|^p,\quad u(0,x)= u_0(x),\quad u_t(0,x)=u_1(x),$$ with $\mu>0$, $n\geq1$, $\sigma \in (0,2]$ and $p>1$. In particular, we are going to prove the non-existence of global weak solutions by using a new test function and suitable sign assumptions on the initial data in both the subcritical case and the critical case.


title: *-Autonomous Envelopes and Conservativity

abstract: We prove 2-categorical conservativity for any {0,T}-free fragment of MALL over its corresponding intuitionistic version: that is, that the universal map from a closed symmetric monoidal category to the *-autonomous category that it freely generates is fully faithful, and similarly for other doctrines. This implies that linear logics and graphical calculi for *-autonomous categories can also be interpreted canonically in closed symmetric monoidal categories. In particular, every closed symmetric monoidal category can be fully embedded in a *-autonomous category, preserving both tensor products and internal-homs. In fact, we prove this directly first with a Yoneda-style embedding (an enhanced "Hyland envelope" that can be regarded as a polycategorical form of Day convolution), and deduce 2-conservativity afterwards from Hyland-Schalk double gluing and a technique of Lafont. The same is true for other fragments of *-autonomous structure, such as linear distributivity, and the embedding can be enhanced to preserve any desired family of nonempty limits and colimits.


title: Network-Aware Optimization of Distributed Learning for Fog Computing

abstract: Fog computing promises to enable machine learning tasks to scale to large amounts of data by distributing processing across connected devices. Two key challenges to achieving this goal are heterogeneity in devices compute resources and topology constraints on which devices can communicate with each other. We address these challenges by developing the first network-aware distributed learning optimization methodology where devices optimally share local data processing and send their learnt parameters to a server for aggregation at certain time intervals. Unlike traditional federated learning frameworks, our method enables devices to offload their data processing tasks to each other, with these decisions determined through a convex data transfer optimization problem that trades off costs associated with devices processing, offloading, and discarding data points. We analytically characterize the optimal data transfer solution for different fog network topologies, showing for example that the value of offloading is approximately linear in the range of computing costs in the network. Our subsequent experiments on testbed datasets we collect confirm that our algorithms are able to improve network resource utilization substantially without sacrificing the accuracy of the learned model. In these experiments, we also study the effect of network dynamics, quantifying the impact of nodes entering or exiting the network on model learning and resource costs.


title: Lax-Sato formulation of the Novikov-Veselov Hierarchy

abstract: We construct a hierarchy of pairwise commuting flows $d/dt_{i,n}$ indexed by $i \in \{1,2 \}$ and $n \in \mathbb{Z}_{\geq 0}$ on triples $(\mathcal{L}_1, \mathcal{L}_2, \mathcal{H})$ where $\partial_1$ and $\partial_2$ are two commuting derivations, $\partial_i \mathcal{L}_i$ is a self-adjoint pseudodifferential operator in $\partial_i$ and $\mathcal{H}$ is the formal Schr\"{o}dinger operator $\mathcal{H}=\partial_1 \partial_2 +u$. $\mathcal{L}_1, \mathcal{L}_2$ and $\mathcal{H}$ are coupled by the relations $\mathcal{H} \mathcal{L}_i+\mathcal{L}_i^* \mathcal{H}=0$. We show that the flows $d/dt_{1,n}+d/dt_{2,n}$ commute with the involution $(\mathcal{L}_1, \mathcal{L}_2, \mathcal{H}) \mapsto (\mathcal{L}_2, \mathcal{L}_1, \mathcal{H})$ and that the first equation of this reduced hierarchy is the Novikov-Veselov equation.


title: Robust Optimization with Decision-Dependent Information Discovery

abstract: Robust optimization is a popular paradigm for modeling and solving two- and multi-stage decision-making problems affected by uncertainty. In many real-world applications, the time of information discovery is decision-dependent and the uncertain parameters only become observable after an often costly investment. Yet, most of the literature assumes that uncertain parameters can be observed for free and that the sequence in which they are revealed is independent of the decision-maker's actions. To fill this gap. we consider two- and multi-stage robust optimization problems in which part of the decision variables control the time of information discovery. Thus, information can be discovered (at least in part) by making strategic exploratory investments in previous stages. We propose a novel dynamic formulation of the problem and prove its correctness. We leverage our model to provide a solution method inspired from the K-adaptability approximation. We reformulate the problem as a finite mixed-integer (resp. bilinear) program if none (resp. some of the) decision variables are real-valued. This finite program is solvable with off-the-shelf solvers. We generalize our approach to the minimization of piecewise linear convex functions. We demonstrate the effectiveness of our method in terms of interpretability, optimality, and speed on synthetic instances of the Pandora box problem, the preference elicitation problem, the best box problem, and the R&D project portfolio optimization problem. Finally, we evaluate it on an instance of the active preference elicitation problem used to recommend kidney allocation policies to policy-makers at the United Network for Organ Sharing.


title: Is the anti-filarial drug diethylcarbamazine useful to treat COVID-19?

abstract: The SARS-CoV-2 virus has resulted in a devastating pandemic of COVID-19. Exploring compounds that could offer a breakthrough in treatment is the need of the hour. Re-positioning cheap, freely available and safe drugs is a priority. The paper proposes evidence for the potential use of diethylcarbamazine (DEC) in the treatment of COVID-19. DEC has inhibitory effects on arachidonic acid metabolism to prostaglandins, little known anti-viral effects on animal retroviruses and demonstrated anti-inflammatory actions in animal models of lung inflammation indicating the need to explore this hypothesis further. We believe this is the first time DEC is being proposed to treat COVID-19.


title: Orbit: Probabilistic Forecast with Exponential Smoothing

abstract: Time series forecasting is an active research topic in academia as well as industry. Although we see an increasing amount of adoptions of machine learning methods in solving some of those forecasting challenges, statistical methods remain powerful while dealing with low granularity data. This paper introduces a refined Bayesian exponential smoothing model with the help of probabilistic programming languages including Stan. Our model refinements include additional global trend, transformation for multiplicative form, noise distribution and choice of priors. A benchmark study is conducted on a rich set of time-series data sets for our models along with other well-known time series models.


title: First integrals on step-two and step-three nilpotent Lie groups

abstract: The goal of this paper is the study of algebraic relations on the Lie algebra of first integrals of the geodesic flow on nilpotent Lie groups equipped with a left-invariant metric. It is proved that the isometry algebra of the $k$-step nilpotent Lie group, $k=2,3$, gives rise to a isomorphic family of first integrals for the geodesic flow. Also invariant first integrals are analyzed and new involution conditions are shown. Finally it is proved that in low dimensions complete families of first integrals can be constructed with Killing vector fields and symmetric Killing 2-tensor fields. This holds for k-step nilpotent Lie algebras of dimension $m\leq 5$ and $k=2,3$. The situation in dimension six is also studied.


title: Convergence for global curve diffusion flows

abstract: In this note we establish exponentially fast smooth convergence for global curve diffusion flows, and discuss open problems relating embeddedness to global existence (Giga's conjecture) and the shape of Type I singularities (Chou's conjecture).


title: BReG-NeXt: Facial Affect Computing Using Adaptive Residual Networks With Bounded Gradient

abstract: This paper introduces BReG-NeXt, a residual-based network architecture using a function wtih bounded derivative instead of a simple shortcut path (a.k.a. identity mapping) in the residual units for automatic recognition of facial expressions based on the categorical and dimensional models of affect. Compared to ResNet, our proposed adaptive complex mapping results in a shallower network with less numbers of training parameters and floating point operations per second (FLOPs). Adding trainable parameters to the bypass function further improves fitting and training the network and hence recognizing subtle facial expressions such as contempt with a higher accuracy. We conducted comprehensive experiments on the categorical and dimensional models of affect on the challenging in-the-wild databases of AffectNet, FER2013, and Affect-in-Wild. Our experimental results show that our adaptive complex mapping approach outperforms the original ResNet consisting of a simple identity mapping as well as other state-of-the-art methods for Facial Expression Recognition (FER). Various metrics are reported in both affect models to provide a comprehensive evaluation of our method. In the categorical model, BReG-NeXt-50 with only 3.1M training parameters and 15 MFLOPs, achieves 68.50% and 71.53% accuracy on AffectNet and FER2013 databases, respectively. In the dimensional model, BReG-NeXt achieves 0.2577 and 0.2882 RMSE value on AffectNet and Affect-in-Wild databases, respectively.


title: Continuous selections, prime number and a covering type property

abstract: Let $(X,\tau)$ be a Hausdorff space and $n\in\omega$. We prove that if $X$ admits a continuous selection over $\mathcal{F}_{n}(X)$ (nonempty subsets of $X$ of cardinality at most $n$), then for every $n\leq m\leq 2n$ such that $m$ is not a prime number, $X$ admits a continuous selection over $[X]^m$ (subsets of $X$ of cardinality $m$). As a consequence of this, a space $X$ admits a continuous selection for every natural number if and only if the same is true for every prime number. For Hausdorff spaces $(X,\tau)$ which admit continuous selections over $[X]^2$, we characterize the existence of continuous selections over $[X]^n$ for $n\geq 2$, in terms of a covering-type property.


title: The Geometric Airy Curve Flow on R^n

abstract: Langer and Perline proved that if x is a solution of the geometric Airy curve flow on R^n then there exists a parallel normal frame along x(. ,t) for each t such that the corresponding principal curvatures satisfy the (n-1) component modified KdV (vmKdV_n). They also constructed higher order curve flows whose principal curvatures are solutions of the higher order flows in the vmKdV_n soliton hierarchy. In this paper, we write down a Poisson structure on the space of curves in R^n parametrized by the arc-length, show that the geometric Airy curve flow is Hamiltonian, write down a sequence of commuting Hamiltonians, and construct Backlund transformations and explicit soliton solutions.


title: Enhanced principle component method for fringe removal in cold atom images

abstract: Many powerful imaging techniques for cold atoms are based on determining the optical density by comparing a beam image having passed through the atom cloud to a reference image taken under similar conditions with no atoms. In practice the beam profile typically contains interference fringes whose phase is not stable between camera exposures. To reduce the error of these fringes in the computed optical density, an algorithm based on principle component analysis (PCA) is often employed. However, PCA is general purpose and not tailored to the specific case of interference fringes. Here we demonstrate an algorithm that takes advantage of the Fourier-space structure of interference fringes to further reduce the residual fringe signatures in the optical density.


title: Design and Control of Roller Grasper V2 for In-Hand Manipulation

abstract: The ability to perform in-hand manipulation still remains an unsolved problem; having this capability would allow robots to perform sophisticated tasks requiring repositioning and reorienting of grasped objects. In this work, we present a novel non-anthropomorphic robot grasper with the ability to manipulate objects by means of active surfaces at the fingertips. Active surfaces are achieved by spherical rolling fingertips with two degrees of freedom (DoF) -- a pivoting motion for surface reorientation -- and a continuous rolling motion for moving the object. A further DoF is in the base of each finger, allowing the fingers to grasp objects over a range of size and shapes. Instantaneous kinematics was derived and objects were successfully manipulated both with a custom handcrafted control scheme as well as one learned through imitation learning, in simulation and experimentally on the hardware.