title: The Impact of Damping in Second-Order Dynamical Systems with Applications to Power Grid Stability

abstract: We consider a broad class of second-order dynamical systems and study the impact of damping as a system parameter on the stability, hyperbolicity, and bifurcation in such systems. We prove a monotonic effect of damping on the hyperbolicity of the equilibrium points of the corresponding first-order system. This provides a rigorous formulation and theoretical justification for the intuitive notion that damping increases stability. To establish this result, we prove a matrix perturbation result for complex symmetric matrices with positive semidefinite perturbations to their imaginary parts, which may be of independent interest. Furthermore, we establish necessary and sufficient conditions for the breakdown of hyperbolicity of the first-order system under damping variations in terms of observability of a pair of matrices relating damping, inertia, and Jacobian matrices, and propose sufficient conditions for Hopf bifurcation resulting from such hyperbolicity breakdown. The developed theory has significant applications in the stability of electric power systems, which are one of the most complex and important engineering systems. In particular, we characterize the impact of damping on the hyperbolicity of the swing equation model which is the fundamental dynamical model of power systems, and demonstrate Hopf bifurcations resulting from damping variations.


title: Intrapersonal Parameter Optimization for Offline Handwritten Signature Augmentation

abstract: Usually, in a real-world scenario, few signature samples are available to train an automatic signature verification system (ASVS). However, such systems do indeed need a lot of signatures to achieve an acceptable performance. Neuromotor signature duplication methods and feature space augmentation methods may be used to meet the need for an increase in the number of samples. Such techniques manually or empirically define a set of parameters to introduce a degree of writer variability. Therefore, in the present study, a method to automatically model the most common writer variability traits is proposed. The method is used to generate offline signatures in the image and the feature space and train an ASVS. We also introduce an alternative approach to evaluate the quality of samples considering their feature vectors. We evaluated the performance of an ASVS with the generated samples using three well-known offline signature datasets: GPDS, MCYT-75, and CEDAR. In GPDS-300, when the SVM classifier was trained using one genuine signature per writer and the duplicates generated in the image space, the Equal Error Rate (EER) decreased from 5.71% to 1.08%. Under the same conditions, the EER decreased to 1.04% using the feature space augmentation technique. We also verified that the model that generates duplicates in the image space reproduces the most common writer variability traits in the three different datasets.


title: Experimentally Testing a Generalized Coarsening Model for Individual Bubbles in Quasi-Two-Dimensional Wet Foams

abstract: We present high-precision data for the time evolution of bubble area $A(t)$ and circularity shape parameter $C(t)$ for quasi-2d foams consisting of bubbles squashed between parallel plates. In order to fully compare with predictions by Roth et al. [Phys. Rev. E 87 2013] and Schimming et al. [Phys. Rev. E 96 2017], foam wetness is systematically varied by controlling the height of the sample above a liquid reservoir which in turn controls the radius $r$ of the inflation of the Plateau borders. For very dry foams, where the borders are very small, classic von Neumann behavior is observed where a bubble's growth rate depends only on its number $n$ of sides. For wet foams, the inflated borders impede gas exchange and cause deviations from von Neumann's law that are found to be in accord with the generalized coarsening equation. In particular, the overall growth rate varies linearly with the film height, which decrease as surface Plateau borders inflate. And, more interestingly, the deviation from $dA/dt\propto (n-6)$ von Neumann behavior grows in proportion to $nCr/\sqrt{A}$. This is highlighted definitively by data for six-sided bubbles, which are forbidden to grow or shrink except for the existence of this term. And it is tested quantitatively by variation of all four relevant quantities: $n$, $C$, $r$, and $A$


title: Mass structure of hadrons and light-front sum rules in t' Hooft model

abstract: We study the mass/energy structure of the bound state of hadrons in two-dimensional quantum chromodynamics in the large number of color limit (t' Hooft model). We analyze separately the contributions from the traceless and trace part of the energy-momentum tensor, and show that the masses are related to the matrix elements of the scalar charge and Coulomb energy. We derive the light-front sum rules for the scalar charge and Coulomb energy, expressed in terms of the light-front wave functions, and find that they are regular at $x=0$ without the delta function contribution. We also consider the result for the massless Goldstone boson, as well as the structure of the gravitational form factors of the bound meson states.


title: Probing for Multilingual Numerical Understanding in Transformer-Based Language Models

abstract: Natural language numbers are an example of compositional structures, where larger numbers are composed of operations on smaller numbers. Given that compositional reasoning is a key to natural language understanding, we propose novel multilingual probing tasks tested on DistilBERT, XLM, and BERT to investigate for evidence of compositional reasoning over numerical data in various natural language number systems. By using both grammaticality judgment and value comparison classification tasks in English, Japanese, Danish, and French, we find evidence that the information encoded in these pretrained models' embeddings is sufficient for grammaticality judgments but generally not for value comparisons. We analyze possible reasons for this and discuss how our tasks could be extended in further studies.


title: Chasing Your Long Tails: Differentially Private Prediction in Health Care Settings

abstract: Machine learning models in health care are often deployed in settings where it is important to protect patient privacy. In such settings, methods for differentially private (DP) learning provide a general-purpose approach to learn models with privacy guarantees. Modern methods for DP learning ensure privacy through mechanisms that censor information judged as too unique. The resulting privacy-preserving models, therefore, neglect information from the tails of a data distribution, resulting in a loss of accuracy that can disproportionately affect small groups. In this paper, we study the effects of DP learning in health care. We use state-of-the-art methods for DP learning to train privacy-preserving models in clinical prediction tasks, including x-ray classification of images and mortality prediction in time series data. We use these models to perform a comprehensive empirical investigation of the tradeoffs between privacy, utility, robustness to dataset shift, and fairness. Our results highlight lesser-known limitations of methods for DP learning in health care, models that exhibit steep tradeoffs between privacy and utility, and models whose predictions are disproportionately influenced by large demographic groups in the training data. We discuss the costs and benefits of differentially private learning in health care.


title: LiDAM: Semi-Supervised Learning with Localized Domain Adaptation and Iterative Matching

abstract: Although data is abundant, data labeling is expensive. Semi-supervised learning methods combine a few labeled samples with a large corpus of unlabeled data to effectively train models. This paper introduces our proposed method LiDAM, a semi-supervised learning approach rooted in both domain adaptation and self-paced learning. LiDAM first performs localized domain shifts to extract better domain-invariant features for the model that results in more accurate clusters and pseudo-labels. These pseudo-labels are then aligned with real class labels in a self-paced fashion using a novel iterative matching technique that is based on majority consistency over high-confidence predictions. Simultaneously, a final classifier is trained to predict ground-truth labels until convergence. LiDAM achieves state-of-the-art performance on the CIFAR-100 dataset, outperforming FixMatch (73.50% vs. 71.82%) when using 2500 labels.


title: Symplectic orbits of unimodular rows

abstract: For a smooth affine algebra $R$ of dimension $d \geq 3$ over a field $k$ and an invertible alternating matrix $\chi$ of rank $2n$, the group $Sp(\chi)$ of invertible matrices of rank $2n$ over $R$ which are symplectic with respect to $\chi$ acts on the right on the set $Um_{2n}(R)$ of unimodular rows of length $2n$ over $R$. In this paper, we prove that $Sp(\chi)$ acts transitively on $Um_{2n}(R)$ if $k$ is algebraically closed, $d! \in k^{\times}$ and $2n \geq d$.


title: HAWC and Fermi-LAT Detection of Extended Emission from the Unidentified Source 2HWC J2006+341

abstract: The discovery of the TeV point source 2HWC J2006+341 was reported in the second HAWC gamma-ray catalog. We present a follow-up study of this source here. The TeV emission is best described by an extended source with a soft spectrum. At GeV energies, an extended source is significantly detected in Fermi-LAT data. The matching locations, sizes and spectra suggest that both gamma-ray detections correspond to the same source. Different scenarios for the origin of the emission are considered and we rule out an association to the pulsar PSR J2004+3429 due to extreme energetics required, if located at a distance of 10.8 kpc.


title: A Multi-Modal Method for Satire Detection using Textual and Visual Cues

abstract: Satire is a form of humorous critique, but it is sometimes misinterpreted by readers as legitimate news, which can lead to harmful consequences. We observe that the images used in satirical news articles often contain absurd or ridiculous content and that image manipulation is used to create fictional scenarios. While previous work have studied text-based methods, in this work we propose a multi-modal approach based on state-of-the-art visiolinguistic model ViLBERT. To this end, we create a new dataset consisting of images and headlines of regular and satirical news for the task of satire detection. We fine-tune ViLBERT on the dataset and train a convolutional neural network that uses an image forensics technique. Evaluation on the dataset shows that our proposed multi-modal approach outperforms image-only, text-only, and simple fusion baselines.


title: Quantum Cycle in Relativistic Non-Commutative Space with Generalized Uncertainty Principle correction

abstract: Quantum heat cycles and quantum refrigerators are analyzed using various quantum systems as their working mediums. For example, to evaluate the efficiency and the work done of the Carnot cycle in the quantum regime, one can consider the harmonic oscillator as it's working medium. For all these well-defined working substances (which are analyzed in commutative space structure), the efficiency of the engine is not up to the mark of the Carnot efficiency. So, one inevitable question arise, can one observe a catalytic effect on the efficiency of the engines and refrigerators when the space structure is changed? In this paper, two different working substance in non-commutative spacetime with relativistic and generalized uncertainty principle corrections has been considered for the analysis of the efficiency of the heat engine cycles. The efficiency of the quantum heat engine gets a boost for higher values of the non-commutative parameter with a harmonic oscillator as the working substance. In the case of the second working medium (one-dimensional infinite potential well), the efficiency shows a constant result in the non-commutative space structure.


title: Alternative notion to intercategories: part I. A tricategory of double categories

abstract: This is a first of a series of two papers. Our motive is to tackle the question raised in B\"ohm's "The Gray Monoidal Product of Double Categories" from Applied Categorical Structures: which would be an alternative notion to intercategories of Grandis and Par\'e, so that monoids in B\"ohm's monoidal category $Dbl$ of strict double categories and double pseudo functors be an example of it? Before addressing this question we observe that although bicategories embed into pseudo double categories, this embedding is not monoidal, with the usual notion of a monoidal pseudo double category. We then prove that monoidal bicategories embed into the mentioned monoids of B\"ohm. In order to fit B\"ohm's monoid into an intercategory-type object, we start by upgrading the category $Dbl$ to a 2-category and end up rather with a tricategory $\DblPs$. We propose an alternative definition of intercategories as internal categories in this tricategory, enabling $Dbl$ to be an example of this gadget. The formal definition of a category internal to the (type of a) tricategory (of) $\DblPs$, as well as another important example of these in the literature, we leave for a subsequent paper.


title: Adaptive Testing for Specification Coverage

abstract: Ensuring correctness of cyber-physical systems (CPS) is an extremely challenging task that is in practice often addressed with simulation based testing. Formal specification languages, such as Signal Temporal Logic (STL), are used to mathematically express CPS requirements and thus render the simulation activity more systematic and principled. We propose a novel method for adaptive generation of tests with specification coverage for STL. To achieve this goal, we devise cooperative reachability games that we combine with numerical optimization to create tests that explore the system in a way that exercise various parts of the specification. To the best of our knowledge our approach is the first adaptive testing approach that can be applied directly to MATLAB\texttrademark\; Simulink/Stateflow models. We implemented our approach in a prototype tool and evaluated it on several illustrating examples and a case study from the avionics domain, demonstrating the effectiveness of adaptive testing to (1) incrementally build a test case that reaches a test objective, (2) generate a test suite that increases the specification coverage, and (3) infer what part of the specification is actually implemented.


title: Noise performance & thermalization of single electron transistors using quantum fluids

abstract: We report on low-temperature noise measurements of a single electron transistor (SET) immersed in superfluid $^4$He. The device acts as a charge sensitive electrometer able to detect the fluctuations of charged defects in close proximity to the SET. In particular, we measure telegraph switching of the electric current through the device originating from a strongly coupled individual two-level fluctuator. By embedding the device in a superfluid helium immersion cell we are able to systematically control the thermalizing environment surrounding the SET and investigate the effect of the superfluid on the SET noise performance. We find that the presence of superfluid $^4$He can strongly suppress the switching rate of the defect by cooling the surrounding phonon bath.


title: On Front-end Gain Invariant Modeling for Wake Word Spotting

abstract: Wake word (WW) spotting is challenging in far-field due to the complexities and variations in acoustic conditions and the environmental interference in signal transmission. A suite of carefully designed and optimized audio front-end (AFE) algorithms help mitigate these challenges and provide better quality audio signals to the downstream modules such as WW spotter. Since the WW model is trained with the AFE-processed audio data, its performance is sensitive to AFE variations, such as gain changes. In addition, when deploying to new devices, the WW performance is not guaranteed because the AFE is unknown to the WW model. To address these issues, we propose a novel approach to use a new feature called $\Delta$LFBE to decouple the AFE gain variations from the WW model. We modified the neural network architectures to accommodate the delta computation, with the feature extraction module unchanged. We evaluate our WW models using data collected from real household settings and showed the models with the $\Delta$LFBE is robust to AFE gain changes. Specifically, when AFE gain changes up to $\pm$12dB, the baseline CNN model lost up to relative 19.0% in false alarm rate or 34.3% in false reject rate, while the model with $\Delta$LFBE demonstrates no performance loss.


title: xVIO: A Range-Visual-Inertial Odometry Framework

abstract: xVIO is a range-visual-inertial odometry algorithm implemented at JPL. It has been demonstrated with closed-loop controls on-board unmanned rotorcraft equipped with off-the-shelf embedded computers and sensors. It can operate at daytime with visible-spectrum cameras, or at night time using thermal infrared cameras. This report is a complete technical description of xVIO. It includes an overview of the system architecture, the implementation of the navigation filter, along with the derivations of the Jacobian matrices which are not already published in the literature.


title: Entanglement dynamics in the three-dimensional Anderson model

abstract: We numerically study the entanglement dynamics of free fermions on a cubic lattice with potential disorder following a quantum quench. We focus, in particular, on the metal-insulator transition at a critical disorder strength and compare the results to the putative many-body localization (MBL) transition in interacting one-dimensional systems. We find that at the transition point the entanglement entropy grows logarithmically with time $t$ while the number entropy grows $\sim\ln\ln t$. This is exactly the same scaling recently found in the MBL phase of the Heisenberg chain with random magnetic fields suggesting that the MBL phase might be more akin to an extended critical regime with both localized and delocalized states rather than a fully localized phase. We also show that the experimentally easily accessible number entropy can be used to bound the full entanglement entropy of the Anderson model and that the critical properties at the metal-insulator transition obtained from entanglement measures are consistent with those obtained by other probes.


title: Three-particle fields as a method of description of baryons in the scattering processes

abstract: In this paper we propose a model within which the keeping of quarks inside a proton and the interaction of these quarks with quarks of another hadron can be described simultaneously. The model of a three-particle bispinor field is used for this purpose. Such a field, in comparison with the usual one-particle field, is considered not on Minkowski space, but on a subset of the tensor product of three Minkowski spaces, which is characterized by the equality of time coordinates of three events (subset of simultaneity). Such a subset cannot be distinguished by the Lorentz invariant method. But due to the multiparticle field transformation law in the transition from one inertial system to another one, the dynamic equations for a three-particle field are consistent with the principle of relativity. The interaction between quarks is introduced by requiring a local SUc(3) invariance. But in the paper a more general way to achieve this invariance than the replacement in Lagrangian of ordinary derivatives by covariance is proposed. The applied procedure leads to a two-index tensor of gauge field, both for Lorentz and internal indices. The connection of this field with the bound states of gluons - glueballs is discussed. It is shown that this glueball field can be represented as the sum of two terms: the classical field, which plays the role of the potential energy of quark interaction in the proton, and the field corresponding to quantum fluctuations around the classical field and describes the quark interaction in proton with other hadrons. A dynamic equation is obtained for the field that determines the potential energy of quark interaction. It is shown that this equation has solutions that correspond to the confinement of quarks.


title: A "head/tail" plasmon model with a Hubble law velocity profile

abstract: We present a model of a hypersonic, collimated, "single pulse" outflow, produced by an event with an ejection velocity that first grows, reaches a peak, and then decreases again to zero velocity in a finite time (simultaneously, the ejection density can have an arbitrary time-variability). We obtain a flow with a leading "head" and a trailing "tail" that for times greater than the width of the pulse develops a linear, "Hubble law" velocity vs. position. We present an analytic model for a simple pulse with a parabolic ejection velocity vs. time and time-independent mass-loss rate, and compare it to an axisymmetric gasdynamic simulation with parameters appropriate for fast knots in planetary nebulae. This "head/tail plasmon" flow might be applicable to other high-velocity clumps with "Hubble law" tails.


title: InsClustering: Instantly Clustering LiDAR Range Measures for Autonomous Vehicle

abstract: LiDARs are usually more accurate than cameras in distance measuring. Hence, there is strong interest to apply LiDARs in autonomous driving. Different existing approaches process the rich 3D point clouds for object detection, tracking and recognition. These methods generally require two initial steps: (1) filter points on the ground plane and (2) cluster non-ground points into objects. This paper proposes a field-tested fast 3D point cloud segmentation method for these two steps. Our specially designed algorithms allow instantly process raw LiDAR data packets, which significantly reduce the processing delay. In our tests on Velodyne UltraPuck, a 32 layers spinning LiDAR, the processing delay of clustering all the $360^\circ$ LiDAR measures is less than 1ms. Meanwhile, a coarse-to-fine scheme is applied to ensure the clustering quality. Our field experiments in public roads have shown that the proposed method significantly improves the speed of 3D point cloud clustering whilst maintains good accuracy.


title: Are all negatives created equal in contrastive instance discrimination?

abstract: Self-supervised learning has recently begun to rival supervised learning on computer vision tasks. Many of the recent approaches have been based on contrastive instance discrimination (CID), in which the network is trained to recognize two augmented versions of the same instance (a query and positive) while discriminating against a pool of other instances (negatives). The learned representation is then used on downstream tasks such as image classification. Using methodology from MoCo v2 (Chen et al., 2020), we divided negatives by their difficulty for a given query and studied which difficulty ranges were most important for learning useful representations. We found a minority of negatives -- the hardest 5% -- were both necessary and sufficient for the downstream task to reach nearly full accuracy. Conversely, the easiest 95% of negatives were unnecessary and insufficient. Moreover, the very hardest 0.1% of negatives were unnecessary and sometimes detrimental. Finally, we studied the properties of negatives that affect their hardness, and found that hard negatives were more semantically similar to the query, and that some negatives were more consistently easy or hard than we would expect by chance. Together, our results indicate that negatives vary in importance and that CID may benefit from more intelligent negative treatment.


title: Medial Injury/Dysfunction Induced Granulation Tissue Repair is the Pathogenesis of Atherosclerosis

abstract: Atherosclerosis, a chronic lesion of vascular wall, remains a leading cause of death and loss of life years. Classical hypotheses for atherosclerosis are long-standing mainly to explain atherogenesis. Unfortunately, these hypotheses may not explain the variation in the susceptibility to atherosclerosis. These issues are controversial over the past 150 years. Atherosclerosis from human coronary arteries was examined and triangle of media was found to be a true portraiture of cells injury in the media, and triangle of intima was a true portraiture of myofibroblast proliferation, extracellular matrix (ECM) secretion, collagen fiber formation and intimal thickening to repair media dysfunction. Myofibroblasts, ECM and lumen (intima)/vasa vasorum (VV) (adventitia) constitute granulation tissue repair. With granulation tissue hyperplasia, lots of collagen fibers (normal or denatured), foam cells and new capillaries formed. Thus, the following theory was postulated: Risk factors induce smooth muscle cells (SMCs) injury/loss, and fibrosis or structure destruction could be developed in the media, which lead to media dysfunction. Media dysfunction prompts disturbed mechanical properties of blood vessels, resulting in bigger pressure buildup in the intima and adventitia. Granulation tissues in the intima/adventitia develop to repair the injured media. Atherosclerosis, stiffening or aneurysm develops depending on media dysfunction severity and granulated tissue repair mode/degree. Nearly all characteristics of clinical atherosclerosis could be ideally interpreted with the theory. We believe that media dysfunction is a key initiator in the pathogenesis of atherosclerosis. It is expected that media dysfunction theory of atherosclerosis, should offer better understanding of the etiology for atherosclerosis.


title: Motif Learning in Knowledge Graphs Using Trajectories Of Differential Equations

abstract: Knowledge Graph Embeddings (KGEs) have shown promising performance on link prediction tasks by mapping the entities and relations from a knowledge graph into a geometric space (usually a vector space). Ultimately, the plausibility of the predicted links is measured by using a scoring function over the learned embeddings (vectors). Therefore, the capability in preserving graph characteristics including structural aspects and semantics highly depends on the design of the KGE, as well as the inherited abilities from the underlying geometry. Many KGEs use the flat geometry which renders them incapable of preserving complex structures and consequently causes wrong inferences by the models. To address this problem, we propose a neuro differential KGE that embeds nodes of a KG on the trajectories of Ordinary Differential Equations (ODEs). To this end, we represent each relation (edge) in a KG as a vector field on a smooth Riemannian manifold. We specifically parameterize ODEs by a neural network to represent various complex shape manifolds and more importantly complex shape vector fields on the manifold. Therefore, the underlying embedding space is capable of getting various geometric forms to encode complexity in subgraph structures with different motifs. Experiments on synthetic and benchmark dataset as well as social network KGs justify the ODE trajectories as a means to structure preservation and consequently avoiding wrong inferences over state-of-the-art KGE models.


title: Four-lepton Z-boson decay constraints on the SMEFT

abstract: We discuss how four-lepton decays of the Z-boson probe currently unconstrained flat directions in the parameter space of the Standard Model Effective Field Theory (SMEFT). We derive the constraints from these decays on four-lepton operators in the SMEFT and show how the LHC data for this process complements probes from neutrino-trident production. Future differential measurements with high-luminosity data can strongly constrain four-lepton operators and remove all flat directions in the four-muon sector of the SMEFT. We comment briefly on the possibility of using rare Z-decays to tau-leptons to probe untested directions in the SMEFT parameter space.


title: Applying Graph-based Deep Learning To Realistic Network Scenarios

abstract: Recent advances in Machine Learning (ML) have shown a great potential to build data-driven solutions for a plethora of network-related problems. In this context, building fast and accurate network models is essential to achieve functional optimization tools for networking. However, state-of-the-art ML-based techniques for network modelling are not able to provide accurate estimates of important performance metrics such as delay or jitter in realistic network scenarios with sophisticated queue scheduling configurations. This paper presents a new Graph-based deep learning model able to estimate accurately the per-path mean delay in networks. The proposed model can generalize successfully over topologies, routing configurations, queue scheduling policies and traffic matrices unseen during the training phase.


title: Symplectic embeddings of four-dimensional polydisks into half integer ellipsoids

abstract: We obtain new sharp obstructions to symplectic embeddings of four-dimensional polydisks $P(a,1)$ into four-dimensional ellipsoids $E(bc,c)$ when $1\le a< 2$ and $b$ is a half-integer. When $1 \leq a < 2-O(b^{-1})$ we demonstrate that $P(a,1)$ symplectically embeds into $E(bc,c)$ if and only if $a+b\le bc$. Our results show that inclusion is optimal and extend the result by Hutchings \cite{H} when $b$ is an integer. Our proof is based on a combinatorial criterion developed by Hutchings \cite{H} to obstruct symplectic embeddings.


title: The Kendall Interaction Filter for Variable Interaction Screening in Ultra High Dimensional Classification Problems

abstract: Accounting for important interaction effects can improve prediction of many statistical learning models. Identification of relevant interactions, however, is a challenging issue owing to their ultrahigh-dimensional nature. Interaction screening strategies can alleviate such issues. However, due to heavier tail distribution and complex dependence structure of interaction effects, innovative robust and/or model-free methods for screening interactions are required to better scale analysis of complex and high-throughput data. In this work, we develop a new model-free interaction screening method, termed Kendall Interaction Filter (KIF), for the classification in high-dimensional settings. The KIF method suggests a weighted-sum measure, which compares the overall to the within-cluster Kendall's $\tau$ of pairs of predictors, to select interactive couples of features. The proposed KIF measure captures relevant interactions for the clusters response-variable, handles continuous, categorical or a mixture of continuous-categorical features, and is invariant under monotonic transformations. We show that the KIF measure enjoys the sure screening property in the high-dimensional setting under mild conditions, without imposing sub-exponential moment assumptions on the features' distributions. We illustrate the favorable behavior of the proposed methodology compared to the methods in the same category using simulation studies, and we conduct real data analyses to demonstrate its utility.


title: 3D flow motions in the viscous sublayer

abstract: We employ novel digital Fresnel reflection holography to capture the 3D flows within the viscous sublayer of a smooth-wall turbulent channel flow at Re{\tau}=400. The measurements reveal unsteady and diverse flow patterns in the sublayer including nearly uniform high and low speed flows and strong small-scale (on the order of viscous wall units) spanwise meandering motions. The probability density functions (PDFs) of wall shear stresses show a clear discrepancy in high stress range with those from direct numerical simulation (DNS), which is attributed to the unresolved streamwise and spanwise motions by DNS. Moreover, the PDF of Lagrangian particle accelerations yields a stretched exponential shape like that in homogenous isotropic turbulence, indicating strong intermittency in the sublayer. We find a significant fraction of high accelerations is associated with the small-scale meandering motions. Our study helps explain the effect of sublayer-scale roughness on reducing drag and flow separation reported in the literature.


title: Disease transmission through expiratory aerosols on an urban bus

abstract: Airborne respiratory diseases such as SARS-CoV-2 (COVID-19) pose significant challenges for public transportation. Several recent outbreaks of SARS-CoV-2 indicate the high risk of transmission among passengers on public buses if special precautions are not taken. This study presents a combined experimental and numerical analysis to identify transmission mechanisms on an urban bus and assess strategies to reduce risk. The effects of the ventilation and air-conditioning systems, opening windows and doors, and wearing masks are analyzed. Specific attention is made to the transport of sub-micron and micron-size particles relevant to typical respiratory droplets. High-resolution instrumentation was used to measure size distribution and aerosol response time on a University of Michigan campus bus under these different conditions. Computational fluid dynamics was employed to measure the airflow within the bus and evaluate risk. A risk metric was adopted based on the number of particles exposed to susceptible passengers. The flow that carries these aerosols is predominantly controlled by the ventilation system, which acts to uniformly distribute the aerosol concentration throughout the bus while simultaneously diluting it with fresh air. The opening of doors and windows was found to reduce the concentration by approximately one half, albeit its benefit does not uniformly impact all passengers on the bus due to recirculation of airflow caused by entrainment through windows. Finally, it was found that well fitted surgical masks, when worn by both infected and susceptible passengers, can nearly eliminate the transmission of the disease.


title: Free energy fluctuations of the $2$-spin spherical SK model at critical temperature

abstract: We investigate the fluctuations of the free energy of the $2$-spin spherical Sherrington-Kirkpatrick model at critical temperature $\beta_c = 1$. When $\beta = 1$ we find asymptotic Gaussian fluctuations with variance $\frac{1}{6N^2} \log(N)$, confirming in the spherical case a physics prediction for the SK model with Ising spins. We furthermore prove the existence of a critical window on the scale $\beta = 1 +\alpha \sqrt{ \log(N) } N^{-1/3}$. For any $\alpha \in \mathbb{R}$ we show that the fluctuations are at most order $\sqrt{ \log(N) } / N$, in the sense of tightness. If $ \alpha \to \infty$ at any rate as $N \to \infty$ then, properly normalized, the fluctuations converge to the Tracy-Widom$_1$ distribution. If $ \alpha \to 0$ at any rate as $N \to \infty$ or $ \alpha <0$ is fixed, the fluctuations are asymptotically Gaussian as in the $\alpha=0$ case. In determining the fluctuations, we apply a recent result of Lambert and Paquette on the behavior of the Gaussian-$\beta$-ensemble at the spectral edge.


title: A correlation between the dark content of elliptical galaxies and their ellipticity

abstract: Observations indicate that the baryonic matter of galaxies is surrounded by vast dark matter halos, which nature remains unknown. This document details the analysis of the results published in MNRAS 438, 2, 1535 (2014) reporting an empirical correlation between the ellipticity of elliptical galaxies and their dark matter content. Large and homogeneous samples of elliptical galaxies for which their dark matter content is inferred were selected using different methods. Possible methodological biases in the dark mass extraction are alleviated by the multiple methods employed. Effects from galaxy peculiarities are minimized by a homogeneity requirement and further suppressed statistically. After forming homogeneous samples (rejection of galaxies with signs of interaction or dependence on their environment, of peculiar elliptical galaxies and of S0-type galaxies) a clear correlation emerges. Such a correlation is either spurious --in which case it signals an ubiquitous systematic bias in elliptical galaxy observations or their analysis-- or genuine --in which case it implies in particular that at equal luminosity, flattened medium-size elliptical galaxies are on average five times heavier than rounder ones, and that the non-baryonic matter content of medium-size round galaxies is small. It would also provides a new testing ground for models of dark matter and galaxy formation.


title: Handwriting Quality Analysis using Online-Offline Models

abstract: This work is part of an innovative e-learning project allowing the development of an advanced digital educational tool that provides feedback during the process of learning handwriting for young school children (three to eight years old). In this paper, we describe a new method for children handwriting quality analysis. It automatically detects mistakes, gives real-time on-line feedback for children's writing, and helps teachers comprehend and evaluate children's writing skills. The proposed method adjudges five main criteria shape, direction, stroke order, position respect to the reference lines, and kinematics of the trace. It analyzes the handwriting quality and automatically gives feedback based on the combination of three extracted models: Beta-Elliptic Model (BEM) using similarity detection (SD) and dissimilarity distance (DD) measure, Fourier Descriptor Model (FDM), and perceptive Convolutional Neural Network (CNN) with Support Vector Machine (SVM) comparison engine. The originality of our work lies partly in the system architecture which apprehends complementary dynamic, geometric, and visual representation of the examined handwritten scripts and in the efficient selected features adapted to various handwriting styles and multiple script languages such as Arabic, Latin, digits, and symbol drawing. The application offers two interactive interfaces respectively dedicated to learners, educators, experts or teachers and allows them to adapt it easily to the specificity of their disciples. The evaluation of our framework is enhanced by a database collected in Tunisia primary school with 400 children. Experimental results show the efficiency and robustness of our suggested framework that helps teachers and children by offering positive feedback throughout the handwriting learning process using tactile digital devices.


title: Easy, Reproducible and Quality-Controlled Data Collection with Crowdaq

abstract: High-quality and large-scale data are key to success for AI systems. However, large-scale data annotation efforts are often confronted with a set of common challenges: (1) designing a user-friendly annotation interface; (2) training enough annotators efficiently; and (3) reproducibility. To address these problems, we introduce Crowdaq, an open-source platform that standardizes the data collection pipeline with customizable user-interface components, automated annotator qualification, and saved pipelines in a re-usable format. We show that Crowdaq simplifies data annotation significantly on a diverse set of data collection use cases and we hope it will be a convenient tool for the community.


title: Non-Bayesian Social Learning on Random Digraphs with Aperiodically Varying Network Connectivity

abstract: We study non-Bayesian social learning on random directed graphs and show that under mild connectivity assumptions, all the agents almost surely learn the true state of the world asymptotically in time if the sequence of the associated weighted adjacency matrices belongs to Class $\pstar$ (a broad class of stochastic chains that subsumes uniformly strongly connected chains). We show that uniform strong connectivity, while being unnecessary for asymptotic learning, ensures that all the agents' beliefs converge to a consensus almost surely, even when the true state is not identifiable. We then provide a few corollaries of our main results, some of which apply to variants of the original update rule such as inertial non-Bayesian learning and learning via diffusion and adaptation. Others include extensions of known results on social learning. We also show that, if the network of influences is balanced in a certain sense, then asymptotic learning occurs almost surely even in the absence of uniform strong connectivity.


title: Digitalization of exceptional simple Lie algebras into matrices over complex numbers

abstract: We give the images of the adjoint representations of exceptional simple Lie algebras by matrices over complex numbers. Next, we digitalize these matrices by the use of Maxima, which is a computer algebra system. These digitalized matrices are provided by using Maxima's function. We prove that these digitalized matrices are closed by Lie bracket operations and make up simple Lie algebras. Moreover, to prove the type of the exceptional simple Lie algebra, we calculate the root system using Maxima for Lie bracket operations as matrix calculations. We show some examples of classical Lie subalgebras of these digitalized matrices.


title: Accelerated computational micromechanics

abstract: We present an approach to solving problems in micromechanics that is amenable to massively parallel calculations through the use of graphical processing units and other accelerators. The problems lead to nonlinear differential equations that are typically second order in space and first order in time. This combination of nonlinearity and nonlocality makes such problems difficult to solve in parallel. However, this combination is a result of collapsing nonlocal, but linear and universal physical laws (kinematic compatibility, balance laws), and nonlinear but local constitutive relations. We propose an operator-splitting scheme inspired by this structure. The governing equations are formulated as (incremental) variational problems, the differential constraints like compatibility are introduced using an augmented Lagrangian, and the resulting incremental variational principle is solved by the alternating direction method of multipliers. The resulting algorithm has a natural connection to physical principles, and also enables massively parallel implementation on structured grids. We present this method and use it to study two examples. The first concerns the long wavelength instability of finite elasticity, and allows us to verify the approach against previous numerical simulations. We also use this example to study convergence and parallel performance. The second example concerns microstructure evolution in liquid crystal elastomers and provides new insights into some counter-intuitive properties of these materials. We use this example to validate the model and the approach against experimental observations.


title: Product risk assessment: a Bayesian network approach

abstract: Product risk assessment is the overall process of determining whether a product, which could be anything from a type of washing machine to a type of teddy bear, is judged safe for consumers to use. There are several methods used for product risk assessment, including RAPEX, which is the primary method used by regulators in the UK and EU. However, despite its widespread use, we identify several limitations of RAPEX including a limited approach to handling uncertainty and the inability to incorporate causal explanations for using and interpreting test data. In contrast, Bayesian Networks (BNs) are a rigorous, normative method for modelling uncertainty and causality which are already used for risk assessment in domains such as medicine and finance, as well as critical systems generally. This article proposes a BN model that provides an improved systematic method for product risk assessment that resolves the identified limitations with RAPEX. We use our proposed method to demonstrate risk assessments for a teddy bear and a new uncertified kettle for which there is no testing data and the number of product instances is unknown. We show that, while we can replicate the results of the RAPEX method, the BN approach is more powerful and flexible.


title: Exit versus escape in a stochastic dynamical system of neuronal networks explains heterogenous bursting intervals

abstract: Neuronal networks can generate burst events. It remains unclear how to analyse interburst periods and their statistics. We study here the phase-space of a mean-field model, based on synaptic short-term changes, that exhibit burst and interburst dynamics and we identify that interburst corresponds to the escape from a basin of attraction. Using stochastic simulations, we report here that the distribution of the these durations do not match with the time to reach the boundary. We further analyse this phenomenon by studying a generic class of two-dimensional dynamical systems perturbed by small noise that exhibits two peculiar behaviors: 1- the maximum associated to the probability density function is not located at the point attractor, which came as a surprise. The distance between the maximum and the attractor increases with the noise amplitude $\sigma$, as we show using WKB approximation and numerical simulations. 2- For such systems, exiting from the basin of attraction is not sufficient to characterize the entire escape time, due to trajectories that can return several times inside the basin of attraction after crossing the boundary, before eventually escaping far away. To conclude, long-interburst durations are inherent properties of the dynamics and sould be expected in empirical time series.


title: Should the Ransomware be Paid?

abstract: Ransomware has emerged as one of the most concerned cyber risks in recent years, which has caused millions of dollars monetary loss over the world. It typically demands a certain amount of ransom payment within a limited timeframe to decrypt the encrypted victim's files. This paper explores whether the ransomware should be paid in a novel game-theoretic model from the perspective of Bayesian game. In particular, the new model analyzes the ransom payment strategies within the framework of incomplete information for both hacker and victim. Our results show that there exist pure and randomized Bayesian Nash equilibria under some mild conditions for the hacker and victim. The sufficient conditions that when the ransom should be paid are presented when an organization is compromised by the ransomware attack. We further study how the costs and probabilities of cracking or recovering affect the expected payoffs of the hacker and the victim in the equilibria. In particular, it is found that the backup option for computer files is not always beneficial, which actually depends on the related cost. Moreover, it is discovered that fake ransomware may be more than expected because of the potential high payoffs. Numerical examples are also presented for illustration.


title: Operator Inference and Physics-Informed Learning of Low-Dimensional Models for Incompressible Flows

abstract: Reduced-order modeling has a long tradition in computational fluid dynamics. The ever-increasing significance of data for the synthesis of low-order models is well reflected in the recent successes of data-driven approaches such as Dynamic Mode Decomposition and Operator Inference. With this work, we suggest a new approach to learning structured low-order models for incompressible flow from data that can be used for engineering studies such as control, optimization, and simulation. To that end, we utilize the intrinsic structure of the Navier-Stokes equations for incompressible flows and show that learning dynamics of the velocity and pressure can be decoupled, thus leading to an efficient operator inference approach for learning the underlying dynamics of incompressible flows. Furthermore, we show the operator inference performance in learning low-order models using two benchmark problems and compare with an intrusive method, namely proper orthogonal decomposition, and other data-driven approaches.


title: Accuracy of meteor positioning from space-and ground-based observations

abstract: Aims. The knowledge of the orbits and origins derived from meteors is important for the study of meteoroids and of the early solar system. With an increase in nano-satellite projects dedicated to Earth observations or directly to meteor observations (e.g., the Meteorix CubeSat), we investigate the stereoscopic measurement of meteor positions using a pair of cameras, one deployed in space and one on the ground, and aim to understand the accuracy and the main driving factors. This study will reveal the requirements for system setups and the geometry favorable for meteor triangulation. Methods. This Letter presents the principle of the stereoscopic measurement from space and the ground, and an error analysis. Specifically, the impacts of the resolutions of the cameras, the attitude and orbit determination accuracy of the satellite, and the geometry formed by the moving target and observers are investigated. Results. To reach a desirable positioning accuracy of 1 km it is necessary to equip the satellite with high-accuracy sensors (e.g., star tracker and GPS receiver) to perform fine attitude and orbit determination. The best accuracy can occur when the target is at an elevation of 30 deg. with respect to the ground station.


title: Critical role of parallel momentum in quantum well state couplings in multi-stacked nanofilms: an angle resolved photoemission study

abstract: We use angle resolved photoemission spectroscopy (ARPES) to investigate the coupling of electron quantum well states (QWS) in epitaxial thin Pb and Ag films. More specifically, we investigate the Ag/Si, Pb/Si, and Pb/Ag/Si systems. We found that the parallel momentum plays a very profound role determining how two adjacent quantum wells are coupled electronically across the interface. We revealed that in the Pb/Ag bimetallic system, there exist two distinctly different regimes in the energy versus momentum (E vs. k) space. In one regime the electronic states in Ag and Pb are strongly coupled resulting in a new set of quantum well states for the bi-metallic system. In the other regime the electronic states in individual metallic layers are retained in their respective regions, as if they are totally decoupled. This result is corroborated by calculations using density functional theory. We further unravel the underlying mechanism associated with the electron refraction and total internal reflection across the interface.


title: Schauder estimates for degenerate L\'evy Ornstein-Ulhenbeck operators

abstract: We establish global Schauder estimates for integro-partial differential equations (IPDE) driven by a possibly degenerate L\'evy Ornstein-Uhlenbeck operator, both in the elliptic and parabolic setting, using some suitable anisotropic H\"older spaces. The class of operators we consider is composed by a linear drift plus a L\'evy operator that is comparable, in a suitable sense, with a possibly truncated stable operator. It includes for example, the relativistic, the tempered, the layered or the Lamperti stable operators. Our method does not assume neither the symmetry of the L\'evy operator nor the invariance for dilations of the linear part of the operator. Thanks to our estimates, we prove in addition the well-posedness of the considered IPDE in suitable functional spaces. In the final section, we extend some of these results to more general operators involving non-linear, space-time dependent drifts.


title: Weakly-Supervised Aspect-Based Sentiment Analysis via Joint Aspect-Sentiment Topic Embedding

abstract: Aspect-based sentiment analysis of review texts is of great value for understanding user feedback in a fine-grained manner. It has in general two sub-tasks: (i) extracting aspects from each review, and (ii) classifying aspect-based reviews by sentiment polarity. In this paper, we propose a weakly-supervised approach for aspect-based sentiment analysis, which uses only a few keywords describing each aspect/sentiment without using any labeled examples. Existing methods are either designed only for one of the sub-tasks, neglecting the benefit of coupling both, or are based on topic models that may contain overlapping concepts. We propose to first learn <sentiment, aspect> joint topic embeddings in the word embedding space by imposing regularizations to encourage topic distinctiveness, and then use neural models to generalize the word-level discriminative information by pre-training the classifiers with embedding-based predictions and self-training them on unlabeled data. Our comprehensive performance analysis shows that our method generates quality joint topics and outperforms the baselines significantly (7.4% and 5.1% F1-score gain on average for aspect and sentiment classification respectively) on benchmark datasets. Our code and data are available at https://github.com/teapot123/JASen.


title: Fault-tolerant and Transactional Stateful Serverless Workflows (extended version)

abstract: This paper introduces Beldi, a library and runtime system for writing and composing fault-tolerant and transactional stateful serverless functions. Beldi runs on existing providers and lets developers write complex stateful applications that require fault tolerance and transactional semantics without the need to deal with tasks such as load balancing or maintaining virtual machines. Beldi's contributions include extending the log-based fault-tolerant approach in Olive (OSDI 2016) with new data structures, transaction protocols, function invocations, and garbage collection. They also include adapting the resulting framework to work over a federated environment where each serverless function has sovereignty over its own data. We implement three applications on Beldi, including a movie review service, a travel reservation system, and a social media site. Our evaluation on 1,000 AWS Lambdas shows that Beldi's approach is effective and affordable.


title: Functors of actions

abstract: In this document, we introduce a novel formalism for any field theory and apply it to the effective field theories of large-scale structure. The new formalism is based on functors of actions composing those theories. This new formalism predicts the actionic fields. We discuss our findings in a cosmological gravitology framework. We present these results with a cosmological inference approach and give guidelines on how we can choose the best candidate between those models with some latest understanding of model selection using Bayesian inference.


title: BD+14$^\circ$3061: A Luminous Yellow Post-AGB Star in the Galactic Halo

abstract: I report the discovery that the 9th-magnitude Galactic-halo star BD+14$^\circ$3061 is a member of the rare class of luminous metal-poor "yellow post-AGB" stars. Its Gaia DR2 parallax implies an absolute magnitude of $M_V=-3.44\pm0.27$, and it is a very high-velocity star moving in a retrograde Galactic orbit. BD+14$^\circ$3061 is a field analog of the half-dozen yellow PAGB stars known in Galactic globular clusters, which have closely similar absolute magnitudes. These objects are the visually brightest members of old stellar populations; their apparently narrow luminosity function makes them potentially useful as Population II standard candles. The spectral-energy distribution of BD+14$^\circ$3061 out to $22\,\mu$m shows no evidence for circumstellar dust. The star is a low-amplitude semi-regular pulsating variable, with typical periods of 30-32 days. A radial-velocity study suggests that it is a spectroscopic binary with a period of 429.6 days, making it similar to known binary yellow PAGB stars such as HD 46703 and BD+39$^\circ$4926.


title: Local Differential Privacy for Bayesian Optimization

abstract: Motivated by the increasing concern about privacy in nowadays data-intensive online learning systems, we consider a black-box optimization in the nonparametric Gaussian process setting with local differential privacy (LDP) guarantee. Specifically, the rewards from each user are further corrupted to protect privacy and the learner only has access to the corrupted rewards to minimize the regret. We first derive the regret lower bounds for any LDP mechanism and any learning algorithm. Then, we present three almost optimal algorithms based on the GP-UCB framework and Laplace DP mechanism. In this process, we also propose a new Bayesian optimization (BO) method (called MoMA-GP-UCB) based on median-of-means techniques and kernel approximations, which complements previous BO algorithms for heavy-tailed payoffs with a reduced complexity. Further, empirical comparisons of different algorithms on both synthetic and real-world datasets highlight the superior performance of MoMA-GP-UCB in both private and non-private scenarios.


title: Language Networks: a Practical Approach

abstract: This manuscript provides a short and practical introduction to the topic of language networks. This text aims at assisting researchers with no practical experience in text and/or network analysis. We provide a practical tutorial on how to model and characterize texts using network-based features. In this tutorial, we also include examples of pre-processing and network representations. A brief description of the main tasks allying network science and text analysis is also provided. A further development of this text shall include a practical description of network classification via machine learning methods.


title: Hydrodynamical simulations of protoplanetary disks including irradiation of stellar photons. I. Resolution study for Vertical Shear Instability (VSI)

abstract: In recent years hydrodynamical (HD) models have become important to describe the gas kinematics in protoplanetary disks, especially in combination with models of photoevaporation and/or magnetic-driven winds. We focus on diagnosing the the vertical extent of the VSI at 203 cells per scale height and allude at what resolution per scale height we obtain convergence. Finally, we determine the regions where EUV, FUV and X-Rays are dominant in the disk. We perform global HD simulations using the PLUTO code. We adopt a global isothermal accretion disk setup, 2.5D (2 dimensions, 3 components) which covers a radial domain from 0.5 to 5.0 and an approximately full meridional extension. We determine the 50 cells per scale height to be the lower limit to resolve the VSI. For higher resolutions, greater than 50 cells per scale height, we observe the convergence for the saturation level of the kinetic energy. We are also able to identify the growth of the `body' modes, with higher growth rate for higher resolution. Full energy saturation and a turbulent steady state is reached after 70 local orbits. We determine the location of the EUV-heated region defined by the radial column density to be 10$^{19}$ cm$^{-2}$ located at $H_\mathrm{R}\sim9.7$, and the FUV/X-Rays-heated boundary layer defined by 10$^{22}$ cm$^{-2}$ located at $H_\mathrm{R}\sim6.2$, making it necessary to introduce the need of a hot atmosphere. For the first time, we report the presence of small scale vortices in the r-Z plane, between the characteristic layers of large scale vertical velocity motions. Such vortices could lead to dust concentration, promoting grain growth. Our results highlight the importance to combine photoevaporation processes in the future high-resolution studies of the turbulence and accretion processes in disks.


title: SafetyPin: Encrypted Backups with Human-Memorable Secrets

abstract: We present the design and implementation of SafetyPin, a system for encrypted mobile-device backups. Like existing cloud-based mobile-backup systems, including those of Apple and Google, SafetyPin requires users to remember only a short PIN and defends against brute-force PIN-guessing attacks using hardware security protections. Unlike today's systems, SafetyPin splits trust over a cluster of hardware security modules (HSMs) in order to provide security guarantees that scale with the number of HSMs. In this way, SafetyPin protects backed-up user data even against an attacker that can adaptively compromise many of the system's constituent HSMs. SafetyPin provides this protection without sacrificing scalability or fault tolerance. Decentralizing trust while respecting the resource limits of today's HSMs requires a synthesis of systems-design principles and cryptographic tools. We evaluate SafetyPin on a cluster of 100 low-cost HSMs and show that a SafetyPin-protected recovery takes 1.01 seconds. To process 1B recoveries a year, we estimate that a SafetyPin deployment would need 3,100 low-cost HSMs.


title: The ping pong ball water cannon

abstract: The course "Phy Ex" was created by Yves Couder in the Paris VII university to teach experimental physics through projects. In this article, we present this teaching method through a particular project that took place in the autumn semester 2019: the ping-pong ball water cannon. In this experiment, a glass containing water and a floating table tennis ball is dropped from some height to the ground. Following the impact, the ball is ejected vertically upwards at speeds that can be several times the impact speed. We report the student team's initial dimensional and order-of-magnitude analysis, and describe the successive experimental setups that showed (1) that free flight is essential for the phenomenon to occur, (2) that the order of magnitude of the ball ejection momentum is correctly predicted by a momentum balance based on integrating the pressure impulse during impact and (3) that making the ball surface more wettable, or stirring the liquid, drastically increases the momentum transfer. The proposed explanation, confirmed by direct high-speed video observations, is that the immersion depth of the ball increases during free fall due to capillary forces or vortex depression-in the absence of buoyancy-and that the enormous excess pressure on the bottom of the ball during impact drives the ball up towards its buoyancy equilibrium. The transfered momentum is sufficient to expel the ball at high velocity, very similar to the formation of liquid jets in collapsing cavities in liquids.


title: CoRel: Seed-Guided Topical Taxonomy Construction by Concept Learning and Relation Transferring

abstract: Taxonomy is not only a fundamental form of knowledge representation, but also crucial to vast knowledge-rich applications, such as question answering and web search. Most existing taxonomy construction methods extract hypernym-hyponym entity pairs to organize a "universal" taxonomy. However, these generic taxonomies cannot satisfy user's specific interest in certain areas and relations. Moreover, the nature of instance taxonomy treats each node as a single word, which has low semantic coverage. In this paper, we propose a method for seed-guided topical taxonomy construction, which takes a corpus and a seed taxonomy described by concept names as input, and constructs a more complete taxonomy based on user's interest, wherein each node is represented by a cluster of coherent terms. Our framework, CoRel, has two modules to fulfill this goal. A relation transferring module learns and transfers the user's interested relation along multiple paths to expand the seed taxonomy structure in width and depth. A concept learning module enriches the semantics of each concept node by jointly embedding the taxonomy and text. Comprehensive experiments conducted on real-world datasets show that Corel generates high-quality topical taxonomies and outperforms all the baselines significantly.


title: Random Network Distillation as a Diversity Metric for Both Image and Text Generation

abstract: Generative models are increasingly able to produce remarkably high quality images and text. The community has developed numerous evaluation metrics for comparing generative models. However, these metrics do not effectively quantify data diversity. We develop a new diversity metric that can readily be applied to data, both synthetic and natural, of any type. Our method employs random network distillation, a technique introduced in reinforcement learning. We validate and deploy this metric on both images and text. We further explore diversity in few-shot image generation, a setting which was previously difficult to evaluate.


title: Sensitivity of BLANC to human-scored qualities of text summaries

abstract: We explore the sensitivity of a document summary quality estimator, BLANC, to human assessment of qualities for the same summaries. In our human evaluations, we distinguish five summary qualities, defined by how fluent, understandable, informative, compact, and factually correct the summary is. We make the case for optimal BLANC parameters, at which the BLANC sensitivity to almost all of summary qualities is about as good as the sensitivity of a human annotator.


title: Downscale energy fluxes in scale invariant oceanic internal wave turbulence

abstract: We analyze analytically and numerically the scale invariant stationary solution to the internal wave kinetic equation. Our analysis of the resonant energy transfers shows that the leading order contributions are given (i) by triads with extreme scale separation and (ii) by triads of waves that are quasi-colinear in the horizontal plane. The contributions from other types of triads is found to be subleading. We use the modified scale invariant limit of the Garrett and Munk spectrum of internal waves to calculate the magnitude of the energy flux towards high wave numbers in both the vertical and the horizontal directions. Our results compare favorably with the finescale parametrization of ocean mixing that was proposed in [Polzin et al. (1995)].


title: Grid-Interactive Multi-Zone Building Control Using Reinforcement Learning with Global-Local Policy Search

abstract: In this paper, we develop a grid-interactive multi-zone building controller based on a deep reinforcement learning (RL) approach. The controller is designed to facilitate building operation during normal conditions and demand response events, while ensuring occupants comfort and energy efficiency. We leverage a continuous action space RL formulation, and devise a two-stage global-local RL training framework. In the first stage, a global fast policy search is performed using a gradient-free RL algorithm. In the second stage, a local fine-tuning is conducted using a policy gradient method. In contrast to the state-of-the-art model predictive control (MPC) approach, the proposed RL controller does not require complex computation during real-time operation and can adapt to non-linear building models. We illustrate the controller performance numerically using a five-zone commercial building.


title: Convergent series of integers with missing digits

abstract: A classical theorem of Kempner states that the sum of the reciprocals of positive integers with missing decimal digits converges. This result is extended to much larger families of "missing digits" sets of positive integers with convergent harmonic series.


title: The global asymptotic structure of period mappings

abstract: This work is part of a project to construct completions of period mappings. A proper topological SBB-esque completion is constructed. The fibres of are projective varieties, and the image is a union of quasi-projective varieties; one wants to endow the topological completion with a compatible algebraic structure. This raises questions about: (i) the global geometry of the fibres; and (ii) the existence of period matrix representations on neighborhoods of such fibres over which the restricted extension is still proper. The purpose of this paper is to investigate these questions.


title: Ensemble Distillation for Structured Prediction: Calibrated, Accurate, Fast-Choose Three

abstract: Modern neural networks do not always produce well-calibrated predictions, even when trained with a proper scoring function such as cross-entropy. In classification settings, simple methods such as isotonic regression or temperature scaling may be used in conjunction with a held-out dataset to calibrate model outputs. However, extending these methods to structured prediction is not always straightforward or effective; furthermore, a held-out calibration set may not always be available. In this paper, we study ensemble distillation as a general framework for producing well-calibrated structured prediction models while avoiding the prohibitive inference-time cost of ensembles. We validate this framework on two tasks: named-entity recognition and machine translation. We find that, across both tasks, ensemble distillation produces models which retain much of, and occasionally improve upon, the performance and calibration benefits of ensembles, while only requiring a single model during test-time.


title: Factorized class $S$ theories and surface defects

abstract: It is known that some theories of class $S$ are actually factorized into multiple decoupled nontrivial four-dimensional $N=2$ theories. We propose a way of constructing examples of this phenomenon using the physics of half-BPS surface defects, and check that it works in one simple example: it correctly reproduces a known realization of two copies of $N=2$ superconformal $SU(2)$ QCD, describing this factorized theory as a class $S$ theory of type $A_3$ on a five-punctured sphere with a twist line. Separately, we also present explicit checks that the Coulomb branch of a putative factorized class $S$ theory has the expected product structure, in two examples.


title: The role of negative emissions in meeting China's 2060 carbon neutrality goal

abstract: China's pledge to reach carbon neutrality before 2060 is an ambitious goal and could provide the world with much-needed leadership on how to limit warming to +1.5C warming above pre-industrial levels by the end of the century. But the pathways that would achieve net zero by 2060 are still unclear, including the role of negative emissions technologies. We use the Global Change Analysis Model to simulate how negative emissions technologies, in general, and direct air capture (DAC) in particular, could contribute to China's meeting this target. Our results show that negative emissions could play a large role, offsetting on the order of 3 GtCO2 per year from difficult-to-mitigate sectors such as freight transportation and heavy industry. This includes up to a 1.6 GtCO2 per year contribution from DAC, constituting up to 60% of total projected negative emissions in China. But DAC, like bioenergy with carbon capture and storage and afforestation, has not yet been demonstrated at anywhere approaching the scales required to meaningfully contribute to climate mitigation. Deploying NETs at these scales will have widespread impacts on financial systems and natural resources such as water, land, and energy in China.


title: "What Are You Trying to Do?" Semantic Typing of Event Processes

abstract: This paper studies a new cognitively motivated semantic typing task, multi-axis event process typing, that, given an event process, attempts to infer free-form type labels describing (i) the type of action made by the process and (ii) the type of object the process seeks to affect. This task is inspired by computational and cognitive studies of event understanding, which suggest that understanding processes of events is often directed by recognizing the goals, plans or intentions of the protagonist(s). We develop a large dataset containing over 60k event processes, featuring ultra fine-grained typing on both the action and object type axes with very large ($10^3\sim 10^4$) label vocabularies. We then propose a hybrid learning framework, P2GT, which addresses the challenging typing problem with indirect supervision from glosses1and a joint learning-to-rank framework. As our experiments indicate, P2GT supports identifying the intent of processes, as well as the fine semantic type of the affected object. It also demonstrates the capability of handling few-shot cases, and strong generalizability on out-of-domain event processes.


title: Clustering of red and blue galaxies around high-redshift 3C radio sources as seen by the Hubble Space Telescope

abstract: To properly understand the evolution of high-redshift galaxy clusters, both passive and star-forming galaxies have to be considered. Here we study the clustering environment of 21 radio galaxies and quasars at 1<z<2.5 from the third Cambridge catalog (3C). We use optical and near-infrared Hubble Space Telescope images with a 2' field-of-view, where the filters encompass the rest-frame 4000 Angstroem break. Passive red and star-forming blue galaxies were separated in the color--magnitude diagram using a redshift-dependent cut derived from galaxy evolution models. We find that about 16 of 21 radio sources inhabit a galaxy overdensity on scales of 250 kpc (30") projected radius. The sample shows a diversity of red and blue overdensities and also sometimes a deficiency of blue galaxies in the center. The following tentative evolutionary trends are seen: extended proto-clusters with only weak overdensities at z > 1.6, red overdensities at 1.2<z<1.6, and red overdensities with an increased deficit of central blue galaxies at z<1.2. Only a few 3C sources show a blue overdensity tracing active star-formation in the cluster centers; this rarity could indicate that the powerful quasar activity may quench star-formation in the vicinity of most radio sources. The derived number of central luminous red galaxies and the radial density profiles are comparable to those found in local clusters, indicating that some 3C clusters are already mass-rich and compact.


title: One-phase free-boundary problems with degeneracy

abstract: In this paper, we study local minimizers of a degenerate version of the Alt-Caffarelli functional. Specifically, we consider local minimizers of the functional $J_{Q}(u, \Omega):= \int_{\Omega} |\nabla u|^2 + Q(x)^2\chi_{\{u>0\}}dx$ where $Q(x) = \text{dist}(x, \Gamma)^{\gamma}$ for $\gamma>0$ and $\Gamma$ a $C^{1, \alpha}$ submanifold of dimension $0 \le k \le n-1$. We show that the free boundary may be decomposed into a rectifiable set, on which we prove upper Minkowski content estimates, and a degenerate cusp set about which little can be said in general with the current techniques. Work in the theory of water waves and the Stokes wave serves as our inspiration, however the main thrust of this paper is to study the geometry of the free boundary for degenerate one-phase Bernoulli free-boundary problems in the context of local minimizers.


title: Joint Constrained Learning for Event-Event Relation Extraction

abstract: Understanding natural language involves recognizing how multiple event mentions structurally and temporally interact with each other. In this process, one can induce event complexes that organize multi-granular events with temporal order and membership relations interweaving among them. Due to the lack of jointly labeled data for these relational phenomena and the restriction on the structures they articulate, we propose a joint constrained learning framework for modeling event-event relations. Specifically, the framework enforces logical constraints within and across multiple temporal and subevent relations by converting these constraints into differentiable learning objectives. We show that our joint constrained learning approach effectively compensates for the lack of jointly labeled data, and outperforms SOTA methods on benchmarks for both temporal relation extraction and event hierarchy construction, replacing a commonly used but more expensive global inference process. We also present a promising case study showing the effectiveness of our approach in inducing event complexes on an external corpus.


title: $L^p$-Bernstein inequalities on $C^2$-domains and applications to discretization

abstract: We prove a new Bernstein type inequality in $L^p$ spaces associated with the tangential derivatives on the boundary of a general compact $C^2$-domain. We give two applications: Marcinkiewicz type inequality for discretization of $L^p$ norm and positive cubature formula. Both results are optimal in the sense that the number of function samples used has the order of the dimension of the corresponding space of algebraic polynomials.


title: Rigidity results on gradient Schouten solitons

abstract: In this paper we consider $\rho$-Einstein solitons of type $M= \left(B^n, g^{*}\right) \times (F^m,g_F)$, where $\left(B^n,g^{*}\right)$ is conformal to a pseudo-Euclidean space and invariant under the action of the pseudo-orthogonal group, and $\left(F^m,g_{F}\right)$ is an Einstein manifold. We provide all the solutions for the gradient Schouten soliton case. Moreover, in the Riemannian case, we prove that if $M= \left(B^n, g^{*}\right) \times (F^m,g_F)$ is a complete gradient Schouten soliton then $\left(B^{n},g^{*}\right)$ is isometric to $\mathbb{S}^{n-1}\times \mathbb{R}$ and $F^m$ is a compact Einstein manifold.


title: Optimal minimax bounds for time and ensemble averages for the incompressible Navier-Stokes equations

abstract: Obtaining sharp estimates for quantities involved in a given model is an integral part of the modeling process. For dynamical systems whose orbits display a complicated, perhaps chaotic, behaviour, the aim is usually to estimate time or ensemble averages of given quantities. This is the case, for instance, in turbulent flows. In this work, the aim is to present a minimax optimization formula that yields optimal bounds for time and/or ensemble averages for the two- and three-dimensional Navier-Stokes equations. The results presented here are extensions to the infinite-dimensional setting of a recent result on the finite-dimensional case given by Tobasco, Goluskin, and Doering in 2017. The optimal result occurs in the form of a minimax optimization problem and does not require knowledge of the solutions, only the law of the system. The minimax optimization problem appears in the form of a maximization over a portion of the phase space of the system and a minimization over a family of auxiliary functions made of cylindrical test functionals defined on the phase space. The function to be optimized is the desired quantity plus the duality product between the law of the system and the derivative of the auxiliary function.


title: Primitive Elements of the Hopf Algebras of Tableaux

abstract: The character theory of symmetric groups, and the theory of symmetric functions, both make use of the combinatorics of Young tableaux, such as the Robinson-Schensted algorithm, Schuetzenberger's "jeu de taquin", and evacuation. In 1995 Poirier and the second author introduced some algebraic structures, different from the plactic monoid, which induce some products and coproducts of tableaux, with homomorphisms. Their starting point are the two dual Hopf algebras of permutations, introduced by the authors in 1995. In 2006 Aguiar and Sottile studied in more detail the Hopf algebra of permutations: among other things, they introduce a new basis, by Moebius inversion in the poset of weak order, that allows them to describe the primitive elements of the Hopf algebra of permutations. In the present note, by a similar method, we determine the primitive elements of the Poirier-Reutenauer algebra of tableaux, using a partial order on tableaux defined by Taskin.


title: A Fast Gateway Placement Algorithm for Flying Networks

abstract: The ability to operate anywhere, anytime, as well as their capability to hover and carry cargo on board make Unmanned Aerial Vehicles (UAVs) suitable platforms to act as Flying Gateways (FGWs) to the Internet. The problem is the optimal placement of the FGWs within the flying network, such that the Quality of Service (QoS) offered is maximized. The literature has been focused on optimizing the placement of the Flying Access Points (FAPs), which establish high-capacity small cells to serve the users on the ground, overlooking the backhaul network design, including the FGW placement. The FGW placement problem is exacerbated in highly dynamic flying networks, where the dynamic traffic demand and the movements of the users may induce frequent changes in the placement of the FAPs. The main contribution of this paper is a fast gateway placement (F-GWP) algorithm for flying networks that determines the optimal position of a FGW. With F-GWP, backhaul communications paths with high enough capacity are established between the FAPs and the FGW, in order to accommodate the traffic demand of the users on the ground. Simulation and experimental results show F-GWP is two orders of magnitude faster than its state of the art counterpart, while ensuring the same flying network performance.


title: On $\beta=6$ Tracy-Widom distribution and the second Calogero-Painlev\'e system

abstract: The Calogero-Painlev\'e systems were introduced in 2001 by K. Takasaki as a natural generalization of the classical Painlev\'e equations to the case of the several Painlev\'e "particles" coupled via the Calogero type interactions. In 2014, I. Rumanov discovered a remarkable fact that a particular case of the second Calogero-Painlev\'e II equation describes the Tracy-Widom distribution function for the general beta-ensembles with the even values of parameter beta. Most recently, in 2017 work of M. Bertola, M. Cafasso, and V. Rubtsov, it was proven that all Calogero-Painlev\'e systems are Lax integrable, and hence their solutions admit a Riemann-Hilbert representation. This important observation has opened the door to rigorous asymptotic analysis of the Calogero-Painlev\'e equations which in turn yields the possibility of rigorous evaluation of the asymptotic behavior of the Tracy-Widom distributions for the values of beta beyond the classical $\beta =1, 2, 4.$ In this work we shall start an asymptotic analysis of the Calogero-Painlev\'e system with a special focus on the Calogero-Painlev\'e system corresponding to $\beta = 6$ Tracy-Widom distribution function. Some technical details were omitted here and will be presented in the next version of the text.


title: Evaluating Tree Explanation Methods for Anomaly Reasoning: A Case Study of SHAP TreeExplainer and TreeInterpreter

abstract: Understanding predictions made by Machine Learning models is critical in many applications. In this work, we investigate the performance of two methods for explaining tree-based models- Tree Interpreter (TI) and SHapley Additive exPlanations TreeExplainer (SHAP-TE). Using a case study on detecting anomalies in job runtimes of applications that utilize cloud-computing platforms, we compare these approaches using a variety of metrics, including computation time, significance of attribution value, and explanation accuracy. We find that, although the SHAP-TE offers consistency guarantees over TI, at the cost of increased computation, consistency does not necessarily improve the explanation performance in our case study.


title: Error-guided likelihood-free MCMC

abstract: This work presents a novel posterior inference method for models with intractable evidence and likelihood functions. Error-guided likelihood-free MCMC, or EG-LF-MCMC in short, has been developed for scientific applications, where a researcher is interested in obtaining approximate posterior densities over model parameters, while avoiding the need for expensive training of component estimators on full observational data or the tedious design of expressive summary statistics, as in related approaches. Our technique is based on two phases. In the first phase, we draw samples from the prior, simulate respective observations and record their errors $\epsilon$ in relation to the true observation. We train a classifier to distinguish between corresponding and non-corresponding $(\epsilon, \boldsymbol{\theta})$-tuples. In the second stage the said classifier is conditioned on the smallest recorded $\epsilon$ value from the training set and employed for the calculation of transition probabilities in a Markov Chain Monte Carlo sampling procedure. By conditioning the MCMC on specific $\epsilon$ values, our method may also be used in an amortized fashion to infer posterior densities for observations, which are located a given distance away from the observed data. We evaluate the proposed method on benchmark problems with semantically and structurally different data and compare its performance against the state of the art approximate Bayesian computation (ABC).


title: Approximation on slabs and uniqueness for inhomogeneous percolation with a plane of defects

abstract: Let $ \mathbb{L}^{d} = ( \mathbb{Z}^{d},\mathbb{E}^{d} ) $ be the $ d $-dimensional hypercubic lattice. We consider a model of inhomogeneous Bernoulli percolation on $ \mathbb{L}^{d} $ in which every edge inside the $ s $-dimensional hyperplane $ \mathbb{Z}^{s} \times \{ 0 \}^{d-s} $, $ 2 \leq s < d $, is open with probability $ q $ and every other edge is open with probability $ p $. We prove the uniqueness of the infinite cluster in the supercritical regime whenever $ p \neq p_{c}(d) $, where $ p_{c}(d) $ denotes the threshold for homogeneous percolation, and that the critical point $ (p,q_{c}(p)) $ can be approximated on the phase space by the critical points of slabs, for any $ p < p_{c}(d) $.


title: Nonparametric, data-based kernel interpolation for particle-tracking simulations and kernel density estimation

abstract: Traditional interpolation techniques for particle tracking include binning and convolutional formulas that use pre-determined (i.e., closed-form, parameteric) kernels. In many instances, the particles are introduced as point sources in time and space, so the cloud of particles (either in space or time) is a discrete representation of the Green's function of an underlying PDE. As such, each particle is a sample from the Green's function; therefore, each particle should be distributed according to the Green's function. In short, the kernel of a convolutional interpolation of the particle sample "cloud" should be a replica of the cloud itself. This idea gives rise to an iterative method by which the form of the kernel may be discerned in the process of interpolating the Green's function. When the Green's function is a density, this method is broadly applicable to interpolating a kernel density estimate based on random data drawn from a single distribution. We formulate and construct the algorithm and demonstrate its ability to perform kernel density estimation of skewed and/or heavy-tailed data including breakthrough curves.


title: Variational Approximation of Factor Stochastic Volatility Models

abstract: Estimation and prediction in high dimensional multivariate factor stochastic volatility models is an important and active research area because such models allow a parsimonious representation of multivariate stochastic volatility. Bayesian inference for factor stochastic volatility models is usually done by Markov chain Monte Carlo methods, often by particle Markov chain Monte Carlo, which are usually slow for high dimensional or long time series because of the large number of parameters and latent states involved. Our article makes two contributions. The first is to propose fast and accurate variational Bayes methods to approximate the posterior distribution of the states and parameters in factor stochastic volatility models. The second contribution is to extend this batch methodology to develop fast sequential variational updates for prediction as new observations arrive. The methods are applied to simulated and real datasets and shown to produce good approximate inference and prediction compared to the latest particle Markov chain Monte Carlo approaches, but are much faster.


title: The Topology of General Cosmological Models

abstract: Is the universe finite or infinite, and what shape does it have? These fundamental questions, of which relatively little is known, are typically studied within the context of the standard model of cosmology where the universe is assumed to be homogeneous and isotropic. Here we address the above questions in highly general cosmological models, with the only assumption being that the average flow of matter is irrotational. Using techniques from differential geometry, specifically extensions of the Bonnet-Myers theorem, we derive a condition which implies a finite universe and yields a bound for its diameter. Furthermore, under a weaker condition involving the interplay between curvature and diameter, together with the assumption that the universe is finite (i.e., has closed spatial slices), we provide a concise list of possible topologies. Namely, the spatial sections then would be either the ring topologies $S^1 \times S^2$, $S^1\tilde{\times}S^2$, $S^1\times\mathbb{RP}^2$, $\mathbb{RP}^3 \# \mathbb{RP}^3$, or covered by the sphere $S^3$ or torus $T^3$. In particular, under this condition the basic construction of connected sums would be ruled out (save for one), along with the plethora of topologies associated with negative curvature. These results are obtained from consequences of the geometrization of 3-manifolds, by applying a generalization of the almost splitting theorem together with a curvature formula of Ehlers and Ellis.


title: Measuring Visual Generalization in Continuous Control from Pixels

abstract: Self-supervised learning and data augmentation have significantly reduced the performance gap between state and image-based reinforcement learning agents in continuous control tasks. However, it is still unclear whether current techniques can face a variety of visual conditions required by real-world environments. We propose a challenging benchmark that tests agents' visual generalization by adding graphical variety to existing continuous control domains. Our empirical analysis shows that current methods struggle to generalize across a diverse set of visual changes, and we examine the specific factors of variation that make these tasks difficult. We find that data augmentation techniques outperform self-supervised learning approaches and that more significant image transformations provide better visual generalization \footnote{The benchmark and our augmented actor-critic implementation are open-sourced @ https://github.com/QData/dmc_remastered)


title: Collective excitations in the tetravalent lanthanide honeycomb antiferromagnet, Na2PrO3

abstract: Thermomagnetic and inelastic neutron scattering measurements on Na2PrO3 are reported. This material is an antiferromagnetic honeycomb magnet based on the tetravalent lanthanide Pr4+ and has been proposed to host dominant antiferromagnetic Kitaev interactions. These measurements reveal magnetic fluctuations in Na2PrO3 below an energy of 2 meV as well as crystal-field excitations around 230 meV. The latter energy is comparable to the scale of the spin-orbit interaction and explains both the very small effective moment of around 1.0 {\mu}B per Pr4+ and the difficulty to uncover any static magnetic scattering below the ordering transition at TN = 4.6 K. By comparing the low-energy magnetic excitations in Na2PrO3 to that of the isostructural spin-only compound, Na2TbO3, a microscopic model of exchange interactions is developed that implicates dominant and surprisingly large Heisenberg exchange interactions J approximately equal to 1.1(1) meV. Although antiferromagnetic Kitaev interactions with K less than or equal to 0.2J cannot be excluded, the inelastic neutron scattering data of Na2PrO3 is best explained with a delta = 1.24(2) easy-axis XXZ exchange anisotropy.


title: Contracts under Moral Hazard and Adverse Selection

abstract: In the classical principal-agent problem, a principal must design a contract to incentivize an agent to perform an action on behalf of the principal. We study the classical principal-agent problem in a setting where the agent can be of one of several types (affecting the outcome of actions they might take). This combines the contract theory phenomena of "moral hazard" (incomplete information about actions) with that of "adverse selection" (incomplete information about types). We examine this problem through the computational lens. We show that in this setting it is APX-hard to compute either the profit-maximizing single contract or the profit-maximizing menu of contracts (as opposed to in the absence of types, where one can efficiently compute the optimal contract). We then show that the performance of the best linear contract scales especially well in the number of types: if agent has $n$ available actions and $T$ possible types, the best linear contract achieves an $O(n\log T)$ approximation of the best possible profit. Finally, we apply our framework to prove tight worst-case approximation bounds between a variety of benchmarks of mechanisms for the principal.


title: Thermal Motors with Enhanced Performance due to Engineered Exceptional Points

abstract: A thermal current, generated by a temperature gradient between two reservoirs coupled to a carefully designed photonic or (micro-) electromechanical circuit, might induce non-conservative forces that impulse a mechanical degree of freedom to move along a closed trajectory. We show that in the limit of long - but finite - modulation periods, the extracted power and the efficiency of such autonomous motors can be maximized when an appropriately designed spatio-temporal symmetry violation is induced and when the motor operates in the vicinity of exceptional point (EP) degeneracies. These singularities appear in the spectrum of the effective non-Hermitian Hamiltonian that describes the combined circuit-reservoirs system when we judiciously tailor the coupling between them. In the photonic framework, these motors can be propelled by thermal radiation and can be utilized as autonomous self-powered microrobots, or as micro-pumps for microfluidics in biological environments. The same designs can be also implemented with electromechanical elements for harvesting ambient mechanical (or electrical) noise for powering a variety of auxiliary systems.


title: Solving Singular Control Problems in Mathematical Biology, Using PASA

abstract: In this paper, we will demonstrate how to use a nonlinear polyhedral constrained optimization solver called the Polyhedral Active Set Algorithm (PASA) for solving a general singular control problem. We present methods of discretizing a general optimal control problem that involves the use of the gradient of the Lagrangian for computing the gradient of the cost functional so that PASA can be applied. When a numerical solution contains artifacts that resemble "chattering'', a phenomenon where the control oscillates wildly along the singular region, we recommend a method of regularizing the singular control problem by adding a term to the cost functional that measures a scalar multiple of the total variation of the control, where the scalar is viewed as a tuning parameter. We then demonstrate PASA's performance on three singular control problems that give rise to different applications of mathematical biology. We also provide some exposition on the heuristics that we use in determining an appropriate size for the tuning parameter.


title: Weakly coupled alternating $S=1/2$ chains in the distorted honeycomb lattice compound Na$_2$Cu$_2$TeO$_6$

abstract: Spin-1/2 chains with alternating antiferromagnetic (AF) and ferromagnetic (FM) couplings exhibit quantum entanglement like the integer-spin Haldane chains and might be similarly utilized for quantum computations. Such alternating AF-FM chains have been proposed to be realized in the distorted honeycomb-lattice compound Na$_2$Cu$_2$TeO$_6$, but to confirm this picture a comprehensive understanding of the exchange interactions including terms outside of the idealized model is required. Here we employ neutron scattering to study the spin dynamics in Na$_2$Cu$_2$TeO$_6$ and accurately determine the coupling strengths through the random phase approximation and density functional theory (DFT) approaches. We find the AF and FM intrachain couplings are the dominant terms in the spin Hamiltonian, while the interchain couplings are AF but perturbative. This hierarchy in the coupling strengths and the alternating signs of the intrachain couplings can be understood through their different exchange paths. Our results establish Na$_2$Cu$_2$TeO$_6$ as a weakly-coupled alternating AF-FM chain compound and reveal the robustness of the gapped ground state in alternating chains under weak interchain couplings.


title: Analogical and Relational Reasoning with Spiking Neural Networks

abstract: Raven's Progressive Matrices have been widely used for measuring abstract reasoning and intelligence in humans. However for artificial learning systems, abstract reasoning remains a challenging problem. In this paper we investigate how neural networks augmented with biologically inspired spiking modules gain a significant advantage in solving this problem. To illustrate this, we first investigate the performance of our networks with supervised learning, then with unsupervised learning. Experiments on the RAVEN dataset show that the overall accuracy of our supervised networks surpass human-level performance, while our unsupervised networks significantly outperform existing unsupervised methods. Finally, our results from both supervised and unsupervised learning illustrate that, unlike their non-augmented counterparts, networks with spiking modules are able to extract and encode temporal features without any explicit instruction, do not heavily rely on training data, and generalise more readily to new problems. In summary, the results reported here indicate that artificial neural networks with spiking modules are well suited to solving abstract reasoning.


title: An Application of Dirac's Interaction Picture to Option Pricing

abstract: In this paper, the Dirac's quantum mechanical interaction picture is applied to option pricing to obtain a solution of the Black-Scholes equation in the presence of a time-dependent arbitrage bubble. In particular, for the case of a call perturbed by a square bubble, an approximate solution (valid up third order in a perturbation series) is given in terms of the three first Greeks: Delta, Gamma, and Speed. Then an exact solution is constructed in terms of all higher order $S$-derivatives of the Black-Scholes formula. It is also shown that the interacting Black-Scholes equation is invariant under a discrete transformation that interchanges the interest rate with the mean of the underlying asset and vice versa. This implies that the interacting Black-Scholes equation can be written in a 'low energy' and a 'high energy' form, in such a way that the high-interaction limit of the low energy form corresponds to the weak-interaction limit of the high energy form. One can apply a perturbative analysis to the high energy form to study the high-interaction limit of the low energy form.