title: Plasmonics as a fabrication tool

abstract: Plasmonic nanolithography is a nanofabrication technique which exploits surface plasmons to generate subdiffraction patterns. In this review, we highlight the key concepts in plasmonic nanolithography, review the pioneering papers and provide a classification that can be useful in determining the proper method and tool for subdiffraction patterning. Considering the ongoing developments in the field of plasmonics, lithography using surface plasmons is an extremely promising approach for development of high-throughput nanofabrication techniques.


title: High Probability Lower Bounds for the Total Variation Distance

abstract: The statistics and machine learning communities have recently seen a growing interest in classification-based approaches to two-sample testing. The outcome of a classification-based two-sample test remains a rejection decision, which is not always informative since the null hypothesis is seldom strictly true. Therefore, when a test rejects, it would be beneficial to provide an additional quantity serving as a refined measure of distributional difference. In this work, we introduce a framework for the construction of high-probability lower bounds on the total variation distance. These bounds are based on a one-dimensional projection, such as a classification or regression method, and can be interpreted as the minimal fraction of samples pointing towards a distributional difference. We further derive asymptotic power and detection rates of two proposed estimators and discuss potential uses through an application to a reanalysis climate dataset.


title: Augmented resolution of linear hyperbolic systems under nonconservative form

abstract: Hyperbolic systems under nonconservative form arise in numerous applications modeling physical processes, for example from the relaxation of more general equations (e.g. with dissipative terms). This paper reviews an existing class of augmented Roe schemes and discusses their application to linear nonconservative hyperbolic systems with source terms. We extend existing augmented methods by redefining them within a common framework which uses a geometric reinterpretation of source terms. This results in intrinsically well-balanced numerical discretizations. We discuss two equivalent formulations: (1) a nonconservative approach and (2) a conservative reformulation of the problem. The equilibrium properties of the schemes are examined and the conditions for the preservation of the well-balanced property are provided. Transient and steady state test cases for linear acoustics and hyperbolic heat equations are presented. A complete set of benchmark problems with analytical solution, including transient and steady situations with discontinuities in the medium properties, are presented and used to assess the equilibrium properties of the schemes. It is shown that the proposed schemes satisfy the expected equilibrium and convergence properties.


title: Coincidences between Calabi-Yau manifolds of Berglund-Hubsch type and Batyrev polytopes

abstract: In this article, we consider the phenomenon of complete coincidence of the key properties of pairs of Calabi-Yau manifolds realized as hypersurfaces in two different weighted projective spaces. More precisely, the first manifold in such a pair is realized as a hypersurface in a weighted projective space, and the second as a hypersurface in the orbifold of another weighted projective space. The two manifolds in each pair have the same Hodge numbers and special K\"ahler geometry on the complex structure moduli space and are associated with the same $N=2$ gauge linear sigma model. We give the explanation of this interesting coincidence using the Batyrev's correspondence between Calabi-Yau manifolds and the reflexive polyhedra.


title: Scalable robustness of interconnected systems subject to structural changes

abstract: This paper studies the robustness of large-scale interconnected systems with respect to external disturbances, focussing on their scalability properties. Specifically, a notion of scalability is introduced that asks for these robustness properties to remain unchanged under a structural change of the system, such as the addition/removal of a subsystem or a change in the interconnection structure. Both necessary and sufficient conditions, in terms of the interconnection structure and edge weights, are given under which elementary structural changes are scalable. The results are illustrated through a simple example.


title: Role of Substrate Stiffness in Tissue Spreading: Wetting Transition and Tissue Durotaxis

abstract: Living tissues undergo wetting transitions: On a surface, they can either form a droplet-like cell aggregate or spread as a monolayer of migrating cells. Tissue wetting depends not only on the chemical but also on the mechanical properties of the substrate. Here, we study the role of substrate stiffness in tissue spreading, which we describe by means of an active polar fluid model. Taking into account that cells exert larger active traction forces on stiffer substrates, we predict a tissue wetting transition at a critical substrate stiffness that decreases with tissue size. On substrates with a stiffness gradient, we find that the tissue spreads faster on the stiffer side. Further, we show that the tissue can wet the substrate on the stiffer side while dewetting from the softer side. We also show that, by means of viscous forces transmitted across the tissue, the stiffer-side interface can transiently drag the softer-side interface towards increasing stiffness, against its spreading tendency. These two effects result in directed tissue migration up the stiffness gradient. This phenomenon -- tissue durotaxis -- can thus emerge both from dewetting at the soft side and from hydrodynamic interactions between the tissue interfaces. Overall, our work unveils mechanisms whereby substrate stiffness impacts the collective migration and the active wetting properties of living tissues, which are relevant in development, regeneration, and cancer.


title: Data Comets: Designing a Visualization Tool for Analyzing Autonomous Aerial Vehicle Logs with Grounded Evaluation

abstract: Autonomous unmanned aerial vehicles are complex systems of hardware, software, and human input. Understanding this complexity is key to their development and operation. Information visualizations already exist for exploring flight logs but comprehensive analyses currently require several disparate and custom tools. This design study helps address the pain points faced by autonomous unmanned aerial vehicle developers and operators. We contribute: a spiral development process model for grounded evaluation visualization development focused on progressively broadening target user involvement and refining user goals; a demonstration of the model as part of developing a deployed and adopted visualization system; a data and task abstraction for developers and operators performing post-flight analysis of autonomous unmanned aerial vehicle logs; the design and implementation of DATA COMETS, an open-source and web-based interactive visualization tool for post-flight log analysis incorporating temporal, geospatial, and multivariate data; and the results of a summative evaluation of the visualization system and our abstractions based on in-the-wild usage. A free copy of this paper and source code are available at osf.io/h4p7g


title: Mega-COV: A Billion-Scale Dataset of 100+ Languages for COVID-19

abstract: We describe Mega-COV, a billion-scale dataset from Twitter for studying COVID-19. The dataset is diverse (covers 268 countries), longitudinal (goes as back as 2007), multilingual (comes in 100+ languages), and has a significant number of location-tagged tweets (~169M tweets). We release tweet IDs from the dataset. We also develop and release two powerful models, one for identifying whether or not a tweet is related to the pandemic (best F1=97%) and another for detecting misinformation about COVID-19 (best F1=92%). A human annotation study reveals the utility of our models on a subset of Mega-COV. Our data and models can be useful for studying a wide host of phenomena related to the pandemic. Mega-COV and our models are publicly available.


title: Two Dimensions for Organizing Immersive Analytics: Toward a Taxonomy for Facet and Position

abstract: As immersive analytics continues to grow as a discipline, so too should its underlying methodological support. Taxonomies play an important role for information visualization and human computer interaction. They provide an organization of the techniques used in a particular domain that better enable researchers to describe their work, discover existing methods, and identify gaps in the literature. Existing taxonomies in related fields do not capture or describe the unique paradigms employed in immersive analytics. We conceptualize a taxonomy that organizes immersive analytics according to two dimensions: spatial and visual presentation. Each intersection of this taxonomy represents a unique design paradigm which, when thoroughly explored, can aid in the design and research of new immersive analytic applications.


title: Well-posedness for a regularised inertial Dean-Kawasaki model for slender particles in several space dimensions

abstract: A stochastic PDE, describing mesoscopic fluctuations in systems of weakly interacting inertial particles of finite volume, is proposed and analysed in any finite dimension $d\in\mathbb{N}$. It is a regularised and inertial version of the Dean-Kawasaki model. A high-probability well-posedness theory for this model is developed. This theory improves significantly on the spatial scaling restrictions imposed in an earlier work of the same authors, which applied only to significantly larger particles in one dimension. The well-posedness theory now applies in $d$-dimensions when the particle-width $\epsilon$ is proportional to $N^{-1/\theta}$ for $\theta>2d$ and $N$ is the number of particles. This scaling is optimal in a certain Sobolev norm. Key tools of the analysis are fractional Sobolev spaces, sharp bounds on Bessel functions, separability of the regularisation in the $d$-spatial dimensions, and use of the Fa\`a di Bruno's formula.


title: Quadratic Hedging for Sequential Claims with Random Weights in Discrete Time

abstract: We study a quadratic hedging problem for a sequence of contingent claims with random weights in discrete time. We obtain the optimal hedging strategy explicitly in a recursive representation, without imposing the non-degeneracy (ND) condition on the model and square integrability on hedging strategies. We relate the general results to hedging under random horizon and fair pricing in the quadratic sense. We illustrate the significance of our results in an example in which the ND condition fails.


title: Detection of cellular micromotion by advanced signal processing

abstract: Cellular micromotion - a tiny movement of cell membranes on the nanometer-micrometer scale - has been proposed as a pathway for inter-cellular signal transduction and as a label-free proxy signal to neural activity. Here we harness several recent approaches of signal processing to detect such micromotion in video recordings of unlabeled cells. Our survey includes spectral filtering of the video signal, matched filtering, as well as 1D and 3D convolutional neural networks acting on pixel-wise time-domain data and a whole recording respectively.


title: Nonlinear dynamics of the fishbone-induced alpha transport on ITER

abstract: The fishbone-induced transport of alpha particles is computed for the ITER 15 MA baseline scenario, using the nonlinear hybrid Kinetic-MHD code XTOR-K. Two limit cases have been studied, in order to analyse the characteristic regimes of the fishbone instability : the weak kinetic drive limit and the strong kinetic drive limit. In both those regimes, characteristic features of the n = m = 1 fishbone instability are recovered, such as a strong up/down-chirping of the mode frequency, associated with a resonant transport of trapped and passing alpha particles. The effects of the n = m = 0 sheared poloidal and toroidal plasma rotation are taken into account in the simulations. The shear is not negligible, which implies that the fishbone mode frequency has a radial dependency, impacting the wave-particle resonance condition. Phase space hole and clump structures are observed in both nonlinear regimes, centered around the precessional and passing resonances. These structures remains attached to the resonances as the different mode frequencies chirp up and down. In the nonlinear phase, the transport of individual resonant trapped particles is identified to be linked to mode-particle synchronization. On this basis, a partial mechanism for the nonlinear coupling between particle transport and mode dominant down-chirping is proposed. The overall transport of alpha particles inside out the q = 1 surface is of order 2-5% of the initial population between the simulations. The loss of alpha power is found to be directly equal to the loss of alpha particles.


title: Particle density in diffusion-limited annihilating systems

abstract: Place an $A$-particle at each site of a graph independently with probability $p$ and otherwise place a $B$-particle. $A$- and $B$-particles perform independent continuous time random walks at rates $\lambda_A$ and $\lambda_B$, respectively, and annihilate upon colliding with a particle of opposite type. Bramson and Lebowitz studied the setting $\lambda_A = \lambda_B$ in the early 1990s. Despite recent progress, many basic questions remain unanswered for when $\lambda_A \neq \lambda_B$. For the critical case $p=1/2$ on low-dimensional integer lattices, we give a lower bound on the expected number of particles at the origin that matches physicists' predictions. For the process with $\lambda_B=0$ on the integers and the bidirected regular tree, we give sharp upper and lower bounds for the expected total occupation time of the root at and approaching criticality.


title: Triposes as a Generalization of Localic Geometric Morphisms

abstract: We discuss how triposes may be understood as generalizations of localic geometric morphisms.


title: Parametrically narrow tetraquark states containing one heavy quark and one heavy antiquark

abstract: This paper focuses on tetraquarks containing one heavy quark and one heavy antiquark in the formal limit where the heavy quark masses go to infinity. It extends the theoretical framework developed for tetraquark containing two heavy quarks based on semi-classical analysis and the Born-Oppenheimer approximation. It is shown that for sufficiently small light quark masses, parametrically narrow $q\bar{q}'Q\bar{Q}$ tetraquarks must exist in the formal limit of arbitrarily large heavy quark masses. The analysis is model-independent; it requires no assumptions about the structure of the state or the mechanism of its decay beyond what can be deduced directly from QCD. The formal analysis here is only directly applicable to tetraquarks with large angular momentum and hence does not directly apply to experimentally observed tetraquark candidates. We also discuss the condition needed for a version of the analysis to apply for low angular momenta states. Whether this condition holds can be explored via lattice calculations. Current lattice calculations suggest that the condition is not met. If these preliminary calculations are confirmed, there are strong implications for the nature of the observed tetraquark resonances containing $\bar{c}c$ or $\bar{b}b$ pairs. It raises the possibility that the observed tetraquarks exist and are narrow not because the charm and bottom quarks are extremely heavy, but because they are sufficiently light.


title: Becoming the Super Turker: Increasing Wages via a Strategy from High Earning Workers

abstract: Crowd markets have traditionally limited workers by not providing transparency information concerning which tasks pay fairly or which requesters are unreliable. Researchers believe that a key reason why crowd workers earn low wages is due to this lack of transparency. As a result, tools have been developed to provide more transparency within crowd markets to help workers. However, while most workers use these tools, they still earn less than minimum wage. We argue that the missing element is guidance on how to use transparency information. In this paper, we explore how novice workers can improve their earnings by following the transparency criteria of Super Turkers, i.e., crowd workers who earn higher salaries on Amazon Mechanical Turk (MTurk). We believe that Super Turkers have developed effective processes for using transparency information. Therefore, by having novices follow a Super Turker criteria (one that is simple and popular among Super Turkers), we can help novices increase their wages. For this purpose, we: (i) conducted a survey and data analysis to computationally identify a simple yet common criteria that Super Turkers use for handling transparency tools; (ii) deployed a two-week field experiment with novices who followed this Super Turker criteria to find better work on MTurk. Novices in our study viewed over 25,000 tasks by 1,394 requesters. We found that novices who utilized this Super Turkers' criteria earned better wages than other novices. Our results highlight that tool development to support crowd workers should be paired with educational opportunities that teach workers how to effectively use the tools and their related metrics (e.g., transparency values). We finish with design recommendations for empowering crowd workers to earn higher salaries.


title: Reputation Agent: Prompting Fair Reviews in Gig Markets

abstract: Our study presents a new tool, Reputation Agent, to promote fairer reviews from requesters (employers or customers) on gig markets. Unfair reviews, created when requesters consider factors outside of a worker's control, are known to plague gig workers and can result in lost job opportunities and even termination from the marketplace. Our tool leverages machine learning to implement an intelligent interface that: (1) uses deep learning to automatically detect when an individual has included unfair factors into her review (factors outside the worker's control per the policies of the market); and (2) prompts the individual to reconsider her review if she has incorporated unfair factors. To study the effectiveness of Reputation Agent, we conducted a controlled experiment over different gig markets. Our experiment illustrates that across markets, Reputation Agent, in contrast with traditional approaches, motivates requesters to review gig workers' performance more fairly. We discuss how tools that bring more transparency to employers about the policies of a gig market can help build empathy thus resulting in reasoned discussions around potential injustices towards workers generated by these interfaces. Our vision is that with tools that promote truth and transparency we can bring fairer treatment to gig workers.


title: Increased-confidence adversarial examples for deep learning counter-forensics

abstract: Transferability of adversarial examples is a key issue to apply this kind of attacks against multimedia forensics (MMF) techniques based on Deep Learning (DL) in a real-life setting. Adversarial example transferability, in fact, would open the way to the deployment of successful counter forensics attacks also in cases where the attacker does not have a full knowledge of the to-be-attacked system. Some preliminary works have shown that adversarial examples against CNN-based image forensics detectors are in general non-transferrable, at least when the basic versions of the attacks implemented in the most popular libraries are adopted. In this paper, we introduce a general strategy to increase the strength of the attacks and evaluate their transferability when such a strength varies. We experimentally show that, in this way, attack transferability can be largely increased, at the expense of a larger distortion. Our research confirms the security threats posed by the existence of adversarial examples even in multimedia forensics scenarios, thus calling for new defense strategies to improve the security of DL-based MMF techniques.


title: Magnetization reversal driven by spin-transfer-torque in perpendicular shape anisotropy magnetic tunnel junctions

abstract: The concept of perpendicular shape anisotropy spin-transfer torque magnetic random-access memory (PSA-STT-MRAM) consists in increasing the storage layer thickness to values comparable to the cell diameter, to induce a perpendicular shape anisotropy in the magnetic storage layer. Making use of that contribution, the downsize scalability of the STT-MRAM may be extended towards sub-20 nm technological nodes, thanks to a reinforcement of the thermal stability factor $\Delta$. Although the larger storage layer thickness improves $\Delta$, it is expected to negatively impact the writing current and switching time. Hence, optimization of the cell dimensions (diameter, thickness) is of utmost importance for attaining a sufficiently high $\Delta$ while keeping a moderate writing current. Micromagnetic simulations were carried out for different pillar thicknesses of fixed lateral size 20 nm. The switching time and the reversal mechanism were analysed as a function of the applied voltage and aspect-ratio (AR) of the storage layer. For AR $<$ 1, the magnetization reversal resembles a macrospin-like mechanism, while for AR $>$ 1 a non-coherent reversal is observed, characterized by the nucleation of a transverse domain wall at the ferromagnet/insulator interface which then propagates along the vertical axis of the pillar. It was further observed that the inverse of the switching time is linearly dependent on the applied voltage. This study was extended to sub-20 nm width with a value of $\Delta$ around 80. It was observed that the voltage necessary to reverse the magnetic layer increases as the lateral size is reduced, accompanied with a transition from macrospin-reversal to a buckling-like reversal at high aspect-ratios.


title: Neutron-$\gamma$ discrimination with organic scintillators: Intrinsic pulse shape and light yield contributions

abstract: A comparative study of the neutron-$\gamma$ Pulse Shape Discrimination (PSD) with seven organic scintillators is performed using an identical setup and digital electronics. The scintillators include plastics (EJ-299-33 and a plastic prototype), single crystals (stilbene and the recent doped $p$-terphenyl) and liquids (BC501A, NE213 and the deuterated liquid BC537). First, the overall PSD performance of the different scintillators is compared and threshold neutron energies for a given discrimination quality are determined. Then, using statistical arguments, two intrinsic contributions to the PSD capability of the scintillating materials are disentangled: the light yield and the specific pulse shapes induced by neutrons and $\gamma$-rays. This separation provides additional insight into the behaviour of organic scintillators and allows a detailed comparison of the discrimination performance of the various materials. On the basis of this analysis, limitations of current organic scintillators and of recently proposed alternative scintillators are discussed.


title: Enabling Deletion in Append-Only Blockchains (Short Summary / Work in Progress)

abstract: Conventional blockchain implementations with append-only semantics do not support deleting or overwriting data in confirmed blocks. However, many industry-relevant use cases require the ability to delete data, especially when personally identifiable information is stored or when data growth has to be constrained. Existing attempts to reconcile these contradictions compromise on core qualities of the blockchain paradigm, as they include backdoor-like approaches such as central authorities with elevated rights or usage of specialized chameleon hash algorithms in chaining of the blocks. In this technical report, we outline a novel architecture for the blockchain ledger and consensus, which uses a tree of context chains with simultaneous validity. A context chain captures the transactions of a closed group of entities and persons, thus structuring blocks in a precisely defined way. The resulting context isolation enables consensus-steered deletion of an entire context without side effects to other contextes. This architecture opens the possibility of truncation, data rollover and separation of concerns, and can help to fulfill the GDPR regulations.


title: The Orbital Pericenter Precession in the 2PN Approximation

abstract: Recent article "Revisiting the 2PN Pericenter Precession in View of Possible Future Measurements" published by Iorio (Universe, 2020) argues that calculations of the secular 2PN precession of the orbital pericenter of a binary system accomplished by Damour and Schaefer (Nuovo Cim. B, 1988) and by Kopeikin and Potapov (Astron. Rep., 1994) with different mathematical techniques are inconsistent, differ from each other and do not agree with the result obtained by Iorio (Universe, 2020). The purpose of this communication is to demonstrate that the article by Iorio is erroneous, suffers from misconceptions and cannot be trusted.


title: The nature of the companion in the Wolf-Rayet system EZ Canis Majoris

abstract: EZ Canis Majoris is a classical Wolf-Rayet star whose binary nature has been debated for decades. It was recently modeled as an eccentric binary with a periodic brightening at periastron of the emission originating in a shock heated zone near the companion. The focus of this paper is to further test the binary model and to constrain the nature of the unseen close companion by searching for emission arising in the shock-heated region. We analyze over 400 high resolution the International Ultraviolet Explorer spectra obtained between 1983 and 1995 and XMM-Newton observations obtained in 2010. The light curve and radial velocity (RV) variations were fit with the eccentric binary model and the orbital elements were constrained. We find RV variations in the primary emission lines with a semi-amplitude K$_1\sim$30 km/s in 1992 and 1995, and a second set of emissions with an anti-phase RV curve with K$_2\sim$150 km/s. The simultaneous model fit to the RVs and the light curve yields the orbital elements for each epoch. Adopting a Wolf-Rayet mass M$_1\sim$20 M$_\odot$ leads to M$_2\sim$3-5 M$_\odot$, which implies that the companion could be a late B-type star. The eccentric (e=0.1) binary model also explains the hard X-ray light curve obtained by XMM-Newton and the fit to these data indicates that the duration of maximum is shorter than the typical exposure times. The anti-phase RV variations of two emission components and the simultaneous fit to the RVs and the light curve are concrete evidence in favor of the binary nature of EZ Canis Majoris. The assumption that the emission from the shock-heated region closely traces the orbit of the companion is less certain, although it is feasible because the companion is significantly heated by the WR radiation field and impacted by the WR wind.


title: Existence of Nonnegative Solutions of Nonlinear Fractional Parabolic Inequalities

abstract: We study the existence of nontrivial nonlocal nonnegative solutions $u(x,t)$ of the nonlinear initial value problems \[ (\partial_t -\Delta)^\alpha u\geq u^\lambda \quad \text{in } \mathbb{R}^n \times\mathbb{R},\,n\geq 1 \] \[ u=0 \quad\text{in } \mathbb{R}^n \times(-\infty,0) \] and \[ C_1 u^\lambda \leq(\partial_t -\Delta)^\alpha u\leq C_2 u^\lambda \quad\text{in } \mathbb{R}^n \times\mathbb{R},\,n\geq1 \] \[ u=0 \quad\text{in } \mathbb{R}^n \times(-\infty,0), \] where $\lambda,\alpha,C_1$, and $C_2$ are positive constants with $C_1 <C_2$. We use the definition of the fractional heat operator $(\partial_t -\Delta)^\alpha$ given in [Taliaferro, 2020] and compare our results in the classical case $\alpha=1$ to known results.


title: Analytic continuation of Bethe energies and application to the thermodynamic limit of the $SL(2,\mathbb{C})$ non-compact spin chains

abstract: We consider the problem of analytically continuing energies computed with the Bethe ansatz, as posed by the study of non-compact integrable spin chains. By introducing an imaginary extensive twist in the Bethe equations, we show that one can expand the analytic continuation of energies in the scaling limit around another 'pseudo-vacuum' sitting at a negative number of Bethe roots, in the same way as around the usual pseudo-vacuum. We show that this method can be used to compute the energy levels of some states of the $SL(2,\mathbb{C})$ integrable spin chain in the infinite-volume limit, and as a proof of principle recover the ground-state value previously obtained in [1] (for the case of spins $s=0, \bar{s}=-1$) by extrapolating results in small sizes. These results represent, as far as we know, the first (partial) description of the spectrum of $SL(2,\mathbb{C})$ non-compact spin chains in the thermodynamic limit.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: Topological energy barrier for skyrmion lattice formation in MnSi

abstract: We report the direct measurement of the topological skyrmion energy barrier through a hysteresis of the skyrmion lattice in the chiral magnet MnSi. Measurements were made using small-angle neutron scattering with a custom-built resistive coil to allow for high-precision minor hysteresis loops. The experimental data was analyzed using an adapted Preisach model to quantify the energy barrier for skyrmion formation and corroborated by the minimum-energy path analysis based on atomistic spin simulations. We reveal that the skyrmion lattice in MnSi forms from the conical phase progressively in small domains, each of which consisting of hundreds of skyrmions, and with an activation barrier of several eV.


title: Kaplan-Meier type survival curves for COVID-19: a health data based decision-making tool

abstract: Countries are recording health information on the global spread of COVID-19 using different methods, sometimes changing the rules after a few days. They are all publishing the number of new individuals infected, cured and dead, along with some supplementary data. These figures are often recorded in a non-uniform manner and do not match the standard definitions of these variables. However, in this paper we show that the Kaplan-Meier curves calculated with them could provide useful information about the dynamics of the disease in different countries. Our aim is to present a robust and simple model to show certain characteristics of the evolution of the dynamic process, showing that the differences of evolution among the countries is reflected in the corresponding Kaplan-Meier-type curves. We compare the curves obtained for the most affected countries so far, proposing possible interpretations of the properties that distinguish them.


title: Strong Asymptotic Composition Theorems for Mutual Information Measures

abstract: We characterize the growth of the Sibson and Arimoto mutual informations and $\alpha$-maximal leakage, of any order that is at least unity, between a random variable and a growing set of noisy, conditionally independent and identically-distributed observations of the random variable. Each of these measures increases exponentially fast to a limit that is order- and measure-dependent, with an exponent that is order- and measure-independent.


title: An adaptive Euler-Maruyama scheme for McKean-Vlasov SDEs with super-linear growth and application to the mean-field FitzHugh-Nagumo model

abstract: In this paper, we introduce adaptive Euler-Maruyama schemes for McKean-Vlasov stochastic differential equations (SDEs) assuming only a standard monotonicity condition on the drift and diffusion coefficients but no global Lipschitz continuity in the state variable for either, while global Lipschitz continuity is required for the measure component only. We prove moment stability of the discretised processes and a strong convergence rate of $1/2$. Several numerical examples, centred around a mean-field model for FitzHugh-Nagumo neurons, illustrate that the standard uniform scheme fails and that the adaptive approach shows in most cases superior performance to tamed approximation schemes. In addition, we introduce and analyse an adaptive Milstein scheme for a certain sub-class of McKean-Vlasov SDEs with linear measure-dependence of the drift.


title: Cross-Modality Relevance for Reasoning on Language and Vision

abstract: This work deals with the challenge of learning and reasoning over language and vision data for the related downstream tasks such as visual question answering (VQA) and natural language for visual reasoning (NLVR). We design a novel cross-modality relevance module that is used in an end-to-end framework to learn the relevance representation between components of various input modalities under the supervision of a target task, which is more generalizable to unobserved data compared to merely reshaping the original representation space. In addition to modeling the relevance between the textual entities and visual entities, we model the higher-order relevance between entity relations in the text and object relations in the image. Our proposed approach shows competitive performance on two different language and vision tasks using public benchmarks and improves the state-of-the-art published results. The learned alignments of input spaces and their relevance representations by NLVR task boost the training efficiency of VQA task.


title: On Operadic Actions on Spaces of Knots and 2-Links

abstract: In the present work, we realize the space of string 2-links $\mathcal{L}$ as a free algebra over a colored operad denoted $\mathcal{SCL}$ (for "Swiss-Cheese for links"). This result extends works of Burke and Koytcheff about the quotient of $\mathcal{L}$ by its center and is compatible with Budney's freeness theorem for long knots. From an algebraic point of view, our main result refines Blaire, Burke and Koytcheff's theorem on the monoid of isotopy classes of string links. Topologically, it expresses the homotopy type of the isotopy class of a string 2-link in terms of the homotopy types of the classes of its prime factors.


title: Computer Vision Toolkit for Non-invasive Monitoring of Factory Floor Artifacts

abstract: Digitization has led to smart, connected technologies be an integral part of businesses, governments and communities. For manufacturing digitization, there has been active research and development with a focus on Cloud Manufacturing (CM) and the Industrial Internet of Things (IIoT). This work presents a computer vision toolkit (CV Toolkit) for non-invasive digitization of the factory floor in line with Industry 4.0 requirements for factory data collection. Currently, technical challenges persist towards digitization of legacy systems due to the limitation for changes in their design and sensors. This novel toolkit is developed to facilitate easy integration of legacy production machinery and factory floor artifacts with the digital and smart manufacturing environment with no requirement of any physical changes in the machines. The system developed is modular, and allows real-time monitoring of production machinery. Modularity aspect allows the incorporation of new software applications in the current framework of CV Toolkit. To allow connectivity of this toolkit with manufacturing floors in a simple, deployable and cost-effective manner, the toolkit is integrated with a known manufacturing data standard, MTConnect, to "translate" the digital inputs into data streams that can be read by commercial status tracking and reporting software solutions. The proposed toolkit is demonstrated using a mock-panel environment developed in house at the University of Cincinnati to highlight its usability.


title: Generalized Multi-view Shared Subspace Learning using View Bootstrapping

abstract: A key objective in multi-view learning is to model the information common to multiple parallel views of a class of objects/events to improve downstream learning tasks. In this context, two open research questions remain: How can we model hundreds of views per event? Can we learn robust multi-view embeddings without any knowledge of how these views are acquired? We present a neural method based on multi-view correlation to capture the information shared across a large number of views by subsampling them in a view-agnostic manner during training. To provide an upper bound on the number of views to subsample for a given embedding dimension, we analyze the error of the bootstrapped multi-view correlation objective using matrix concentration theory. Our experiments on spoken word recognition, 3D object classification and pose-invariant face recognition demonstrate the robustness of view bootstrapping to model a large number of views. Results underscore the applicability of our method for a view-agnostic learning setting.

categories: cs.LG cs.CV cs.SD eess.AS eess.SP stat.ML


title: Evaluating the Effect of Timeline Shape on Visualization Task Performance

abstract: Timelines are commonly represented on a horizontal line, which is not necessarily the most effective way to visualize temporal event sequences. However, few experiments have evaluated how timeline shape influences task performance. We present the design and results of a controlled experiment run on Amazon Mechanical Turk (n=192) in which we evaluate how timeline shape affects task completion time, correctness, and user preference. We tested 12 combinations of 4 shapes -- horizontal line, vertical line, circle, and spiral -- and 3 data types -- recurrent, non-recurrent, and mixed event sequences. We found good evidence that timeline shape meaningfully affects user task completion time but not correctness and that users have a strong shape preference. Building on our results, we present design guidelines for creating effective timeline visualizations based on user task and data types. A free copy of this paper, the evaluation stimuli and data, and code are available at https://osf.io/qr5yu/


title: Occlusion-Adaptive Deep Network for Robust Facial Expression Recognition

abstract: Recognizing the expressions of partially occluded faces is a challenging computer vision problem. Previous expression recognition methods, either overlooked this issue or resolved it using extreme assumptions. Motivated by the fact that the human visual system is adept at ignoring the occlusion and focus on non-occluded facial areas, we propose a landmark-guided attention branch to find and discard corrupted features from occluded regions so that they are not used for recognition. An attention map is first generated to indicate if a specific facial part is occluded and guide our model to attend to non-occluded regions. To further improve robustness, we propose a facial region branch to partition the feature maps into non-overlapping facial blocks and task each block to predict the expression independently. This results in more diverse and discriminative features, enabling the expression recognition system to recover even though the face is partially occluded. Depending on the synergistic effects of the two branches, our occlusion-adaptive deep network significantly outperforms state-of-the-art methods on two challenging in-the-wild benchmark datasets and three real-world occluded expression datasets.


title: Guaranteeing Reproducibility in Deep Learning Competitions

abstract: To encourage the development of methods with reproducible and robust training behavior, we propose a challenge paradigm where competitors are evaluated directly on the performance of their learning procedures rather than pre-trained agents. Since competition organizers re-train proposed methods in a controlled setting they can guarantee reproducibility, and -- by retraining submissions using a held-out test set -- help ensure generalization past the environments on which they were trained.


title: Reliable Frequency Regulation through Vehicle-to-Grid: Encoding Legislation with Robust Constraints

abstract: Problem definition: Vehicle-to-grid increases the low utilization rate of privately owned electric vehicles by making their batteries available to electricity grids. We formulate a robust optimization problem that maximizes a vehicle owner's expected profit from selling primary frequency regulation to the grid and guarantees that market commitments are met at all times for all frequency deviation trajectories in a functional uncertainty set that encodes applicable legislation. Faithfully modeling the energy conversion losses during battery charging and discharging renders this optimization problem non-convex. Methodology/results: By exploiting a total unimodularity property of the uncertainty set and an exact linear decision rule reformulation, we prove that this non-convex robust optimization problem with functional uncertainties is equivalent to a tractable linear program. Through extensive numerical experiments using real-world data, we quantify the economic value of vehicle-to-grid and elucidate the financial incentives of vehicle owners, aggregators, equipment manufacturers, and regulators. Managerial implications: We find that the prevailing penalties for non-delivery of promised regulation power are too low to incentivize vehicle owners to honor the delivery guarantees given to grid operators.


title: Serdab: An IoT Framework for Partitioning Neural Networks Computation across Multiple Enclaves

abstract: Recent advances in Deep Neural Networks (DNN) and Edge Computing have made it possible to automatically analyze streams of videos from home/security cameras over hierarchical clusters that include edge devices, close to the video source, as well as remote cloud compute resources. However, preserving the privacy and confidentiality of users' sensitive data as it passes through different devices remains a concern to most users. Private user data is subject to attacks by malicious attackers or misuse by internal administrators who may use the data in activities that are not explicitly approved by the user. To address this challenge, we present Serdab, a distributed orchestration framework for deploying deep neural network computation across multiple secure enclaves (e.g., Intel SGX). Secure enclaves provide a guarantee on the privacy of the data/code deployed inside it. However, their limited hardware resources make them inefficient when solely running an entire deep neural network. To bridge this gap, Serdab presents a DNN partitioning strategy to distribute the layers of the neural network across multiple enclave devices or across an enclave device and other hardware accelerators. Our partitioning strategy achieves up to 4.7x speedup compared to executing the entire neural network in one enclave.


title: A new approach to the Thomas-Fermi boundary-value problem

abstract: Given the Thomas-Fermi equation sqrt(x)phi''=phi*(3/2), this paper changes first the dependent variable by defining y(x)=sqrt(x phi(x)). The boundary conditions require that y(x) must vanish at the origin as sqrt(x), whereas it has a fall-off behaviour at infinity proportional to the power (1/2)(1-chi) of the independent variable x, chi being a positive number. Such boundary conditions lead to a 1-parameter family of approximate solutions in the form sqrt(x) times a ratio of finite linear combinations of integer and half-odd powers of x. If chi is set equal to 3, in order to agree exactly with the asymptotic solution of Sommerfeld, explicit forms of the approximate solution are obtained for all values of x. They agree exactly with the Majorana solution at small x, and remain very close to the numerical solution for all values of x. Remarkably, without making any use of series, our approximate solutions achieve a smooth transition from small-x to large-x behaviour. Eventually, the generalized Thomas-Fermi equation that includes relativistic, non-extensive and thermal effects is studied, finding approximate solutions at small and large x for small or finite values of the physical parameters in this equation.


title: Power Quality Monitor for Residential Voltage

abstract: A simple and inexpensive system capable of displaying relevant power quality data of a residential voltage signal is presented, and it is based on Arduino, Python and Kivy. A circuit for adjusting the voltage signal for sampling is suggested. The software of the system allows the transmission of data with any device capable of serial communication at a 2000000 Baud rate, although an Arduino NANO/UNO board is recommended A board based on an Arduino NANO R3 was used. The software plots the last 6 cycles contained in the data by default, but lets analyze any amount of cycles. It provides the RMS voltage, the peak voltage and the total harmonic distortion of the scrutinized cycles, up to the 25th harmonic. The software is also capable of calculating the fast Fourier transform of the studied cycles, and returning the corresponding plot. A sample of data captured at Maracaibo, Venezuela is shown.


title: Red-Blue Point Separation for Points on a Circle

abstract: Given a set R of red points and a set B of blue points in the plane, the Red-Blue point separation problem asks if there are at most k lines that separate R from B, that is, each cell induced by the lines of the solution is either empty or monochromatic (containing points of only one color). A common variant of the problem is when the lines are required to be axis-parallel. The problem is known to be NP-complete for both scenarios, and W[1]-hard parameterized by k in the former setting and FPT in the latter. We demonstrate a polynomial-time algorithm for the special case when the points lie on a circle. Further, we also demonstrate the W-hardness of a related problem in the axis-parallel setting, where the question is if there are p horizontal and q vertical lines that separate R from B. The hardness here is shown in the parameter p.


title: Compositional Few-Shot Recognition with Primitive Discovery and Enhancing

abstract: Few-shot learning (FSL) aims at recognizing novel classes given only few training samples, which still remains a great challenge for deep learning. However, humans can easily recognize novel classes with only few samples. A key component of such ability is the compositional recognition that human can perform, which has been well studied in cognitive science but is not well explored in FSL. Inspired by such capability of humans, to imitate humans' ability of learning visual primitives and composing primitives to recognize novel classes, we propose an approach to FSL to learn a feature representation composed of important primitives, which is jointly trained with two parts, i.e. primitive discovery and primitive enhancing. In primitive discovery, we focus on learning primitives related to object parts by self-supervision from the order of image splits, avoiding extra laborious annotations and alleviating the effect of semantic gaps. In primitive enhancing, inspired by current studies on the interpretability of deep networks, we provide our composition view for the FSL baseline model. To modify this model for effective composition, inspired by both mathematical deduction and biological studies (the Hebbian Learning rule and the Winner-Take-All mechanism), we propose a soft composition mechanism by enlarging the activation of important primitives while reducing that of others, so as to enhance the influence of important primitives and better utilize these primitives to compose novel classes. Extensive experiments on public benchmarks are conducted on both the few-shot image classification and video recognition tasks. Our method achieves the state-of-the-art performance on all these datasets and shows better interpretability.


title: Multiple-line study of molecular gas in spiral galaxy NGC 2903

abstract: Multiple molecular lines with radiative transfer modelling are a powerful tool to probe the physics of star-forming gas in galaxies. We investigate the gas properties in the centre of the spiral galaxy NGC 2903 using low-J CO lines, i.e. 12CO(1-0), 12CO(2-1), 12CO(3-2), 13CO(1-0), and HCN(1-0). We apply a nonlocal thermodynamic equilibrium radiative transfer code to derive beam-averaged molecular gas properties. We use two methods (i.e. chi-square minimizations and likelihood) to define the best model representing the observed line ratios best. The line ratio diagnostics suggest that CO gas in the centre of NGC 2903 is thinner and the dense gas fraction is similar compared to that of spirals, starburst, and early-type galaxies (ETGs), while the gas in the centre of the galaxy is warmer than that of ETGs and colder than that of starbursts. Based on the best-fitting model results, we find that the beam-averaged gas kinetic temperature is T_K = 20 K, H_2 volume number density is log(n(H_2)) = 4.2 cm-3 , and CO column number density is log(N(CO)) = 19.0 cm-2 in the centre of NGC 2903. Both methods, i.e. the line ratio diagnostics and modelling, indicate an ISM in the centre of NGC 2903 having intermediate temperature and star formation activity (also supported by star formation rates), thinner CO gas with similar dense gas fraction, and higher H_2 volume number density compared to that of spirals, ETGs, and starbursts.


title: Efficient experimental quantum fingerprinting with wavelength division multiplexing

abstract: Quantum communication complexity studies the efficiency of information communication (that is, the minimum amount of communication required to achieve a certain task) using quantum states. One representative example is quantum fingerprinting, in which the simultaneous message passing model is considered and the minimum amount of communication could be exponentially smaller than the classical fingerprinting protocol. Experimental demonstrations based on a practical quantum fingerprinting protocol where coherent states are used to construct the fingerprints, have successfully shown the superiority of quantum fingerprinting. However, as a consequence of using coherent states, the communication time in this coherent quantum fingerprinting (CQF) protocol is increased. Moreover, the minimum amount of information communicated in these experimental demonstrations is largely limited by the dark counts of the single photon detectors. Here, we propose to enhance the performance of the existing CQF protocol through applying wavelength division multiplexing (WDM) and simultaneous detection of multiple wavelength channels. We show that the new WDM-CQF protocol can potentially reduce the communication time significantly. More importantly, with the same experimental parameters, the amount of communication is much reduced in the new scheme compared with the original CQF protocol. The more wavelength channels are used, the less communication is required by WDM-CQF protocol. We also perform a proof-of-concept experimental demonstration of the new WDM-CQF protocol with 6 wavelength channels. The experimental results further validate that the new scheme with WDM not only beats the classical protocol, but also reduces the amount of communication required by the original CQF protocol by more than half.


title: Class-Incremental Learning for Semantic Segmentation Re-Using Neither Old Data Nor Old Labels

abstract: While neural networks trained for semantic segmentation are essential for perception in autonomous driving, most current algorithms assume a fixed number of classes, presenting a major limitation when developing new autonomous driving systems with the need of additional classes. In this paper we present a technique implementing class-incremental learning for semantic segmentation without using the labeled data the model was initially trained on. Previous approaches still either rely on labels for both old and new classes, or fail to properly distinguish between them. We show how to overcome these problems with a novel class-incremental learning technique, which nonetheless requires labels only for the new classes. Specifically, (i) we introduce a new loss function that neither relies on old data nor on old labels, (ii) we show how new classes can be integrated in a modular fashion into pretrained semantic segmentation models, and finally (iii) we re-implement previous approaches in a unified setting to compare them to ours. We evaluate our method on the Cityscapes dataset, where we exceed the mIoU performance of all baselines by 3.5% absolute reaching a result, which is only 2.2% absolute below the upper performance limit of single-stage training, relying on all data and labels simultaneously.


title: Laser spectroscopy of indium Rydberg atom bunches by electric field ionization

abstract: This work reports on the application of a novel electric field-ionization setup for high-resolution laser spectroscopy measurements on bunched fast atomic beams in a collinear geometry. In combination with multi-step resonant excitation to Rydberg states using pulsed lasers, the field ionization technique demonstrates increased sensitivity for isotope separation and measurement of atomic parameters over non-resonant laser ionization methods. The setup was tested at the Collinear Resonance Ionization Spectroscopy experiment at ISOLDE-CERN to perform high-resolution measurements of transitions in the indium atom from the 5s$^2$5d~$^2$D$_{5/2}$ and 5s$^2$5d~$^2$D$_{3/2}$ states to 5s$^2$($n$)p~$^2$P and 5s$^2$($n$)f~$^2$F Rydberg states, up to a principal quantum number of $n$ = 72. The extracted Rydberg level energies were used to re-evaluate the ionization potential of the indium atom to be 46670.1055(21) cm$^{-1}$. The nuclear magnetic dipole and nuclear electric quadrupole hyperfine structure constants and level isotope shifts of the 5s$^2$5d~$^2$D$_{5/2}$ and 5s$^2$5d~$^2$D$_{3/2}$ states were determined for $^{113,115}$In. The results are compared to calculations using relativistic coupled-cluster theory. A good agreement is found with the ionization potential and isotope shifts, while disagreement of hyperfine structure constants indicates an increased importance of electron correlations in these excited atomic states. With the aim of further increasing the detection sensitivity for measurements on exotic isotopes, a systematic study of the field-ionization arrangement implemented in the work was performed and an improved design was simulated and is presented. The improved design offers increased background suppression independent of the distance from field ionization to ion detection.


title: Diffusion in unimodular gravity: Analytical solutions, late-time acceleration, and cosmological constraints

abstract: Unimodular gravity is an appealing approach to address the cosmological constant problem. In this scenario, the vacuum energy density of quantum fields does not gravitate and the cosmological constant appears merely as an integration constant. Recently, it has been shown that energy diffusion that may arise in quantum gravity and in theories with spontaneous collapse is compatible with this framework by virtue of its restricted diffeomorphism invariance. New studies suggest that this phenomenon could lead to higher-order equations in the context of homogeneous and isotropic Universe, affecting the well-posedness of their Cauchy initial-value problem. In this work, we show that this issue can be circumvented by assuming an equation of state that relates the energy density to the function that characterizes the diffusion. As an application, we solve the field equations analytically for an isotropic and homogeneous Universes in a barotropic model and in the mass-proportional continuous spontaneous localization (CSL) scenario, assuming that only dark matter develops energy diffusion. Different solutions possessing phase transition from decelerated to accelerated expansion are found. We use cosmological data of type Ia Supernovae and observational Hubble data to constrain the free parameters of both models. It is found that very small but nontrivial energy nonconservation is compatible with the barotropic model. However, for the CSL model, we find that the best-fit values are not compatible with previous laboratory experiments. We comment on this fact and propose future directions to explore energy diffusion in cosmology.


title: Photonic production of the pair of $B_c$ mesons

abstract: We study the pair production of $B_c$ mesons in the photon-photon interaction in the framework of perturbative quantum chromodynamics and the relativistic quark model. The production amplitudes of a pair of pseudoscalar and vector $B_c$ mesons are constructed in the nonrelativistic approximation and taking into account relativistic effects. Relativistic corrections related to the relative motion of heavy quarks in the production amplitude, as well as in the wave function of the bound state of heavy quarks, are taken into account. Analytical expressions are constructed for the relativistic differential and total cross sections for the pair $B_c$ meson production. Based on them, numerical values of the production cross sections are obtained for various energies and scattering angles.


title: The importance of the diffuse ionized gas for interpreting galaxy spectra

abstract: Diffuse ionized gas (DIG) in galaxies can be found in early-type galaxies, in bulges of late-type galaxies, in the interarm regions of galaxy disks, and outside the plane of such disks. The emission-line spectrum of the DIG can be confused with that of a weakly active galactic nucleus. It can also bias the inference of chemical abundances and star formation rates in star forming galaxies. We discuss how one can detect and feasibly correct for the DIG contribution in galaxy spectra.


title: Estimating the integral length scale on turbulent flows from the zero crossings of the longitudinal velocity fluctuation

abstract: The integral length scale ($\mathcal{L}$) is considered to be characteristic of the largest motions of a turbulent flow, and as such, it is an input parameter in modern and classical approaches of turbulence theory and numerical simulations. Its experimental estimation, however, could be difficult in certain conditions, for instance, when the experimental calibration required to measure $\mathcal{L}$ is hard to achieve (hot-wire anemometry on large scale wind-tunnels, and field measurements), or in 'standard' facilities using active grids due to the behaviour of their velocity autocorrelation function $\rho(r)$, which does not in general cross zero. In this work, we provide two alternative methods to estimate $\mathcal{L}$ using the variance of the distance between successive zero crossings of the streamwise velocity fluctuations, thereby reducing the uncertainty of estimating $\mathcal{L}$ under similar experimental conditions. These methods are applicable to variety of situations such as active grids flows, field measurements, and large scale wind tunnels.


title: COVID-19 Contact Tracing and Privacy: Studying Opinion and Preferences

abstract: There is growing interest in technology-enabled contact tracing, the process of identifying potentially infected COVID-19 patients by notifying all recent contacts of an infected person. Governments, technology companies, and research groups alike recognize the potential for smartphones, IoT devices, and wearable technology to automatically track "close contacts" and identify prior contacts in the event of an individual's positive test. However, there is currently significant public discussion about the tensions between effective technology-based contact tracing and the privacy of individuals. To inform this discussion, we present the results of a sequence of online surveys focused on contact tracing and privacy, each with 100 participants. Our first surveys were on April 1 and 3, and we report primarily on those first two surveys, though we present initial findings from later survey dates as well. Our results present the diversity of public opinion and can inform the public discussion on whether and how to leverage technology to reduce the spread of COVID-19. We are continuing to conduct longitudinal measurements, and will update this report over time; citations to this version of the report should reference Report Version 1.0, May 8, 2020. NOTE: As of December 4, 2020, this report has been superseded by Report Version 2.0, found at arXiv:2012.01553. Please read and cite Report Version 2.0 instead.


title: Visual Analytics and Human Involvement in Machine Learning

abstract: The rapidly developing AI systems and applications still require human involvement in practically all parts of the analytics process. Human decisions are largely based on visualizations, providing data scientists details of data properties and the results of analytical procedures. Different visualizations are used in the different steps of the Machine Learning (ML) process. The decision which visualization to use depends on factors, such as the data domain, the data model and the step in the ML process. In this chapter, we describe the seven steps in the ML process and review different visualization techniques that are relevant for the different steps for different types of data, models and purposes.


title: That is a Known Lie: Detecting Previously Fact-Checked Claims

abstract: The recent proliferation of "fake news" has triggered a number of responses, most notably the emergence of several manual fact-checking initiatives. As a result and over time, a large number of fact-checked claims have been accumulated, which increases the likelihood that a new claim in social media or a new statement by a politician might have already been fact-checked by some trusted fact-checking organization, as viral claims often come back after a while in social media, and politicians like to repeat their favorite statements, true or false, over and over again. As manual fact-checking is very time-consuming (and fully automatic fact-checking has credibility issues), it is important to try to save this effort and to avoid wasting time on claims that have already been fact-checked. Interestingly, despite the importance of the task, it has been largely ignored by the research community so far. Here, we aim to bridge this gap. In particular, we formulate the task and we discuss how it relates to, but also differs from, previous work. We further create a specialized dataset, which we release to the research community. Finally, we present learning-to-rank experiments that demonstrate sizable improvements over state-of-the-art retrieval and textual similarity approaches.


title: A computational model implementing subjectivity with the 'Room Theory'. The case of detecting Emotion from Text

abstract: This work introduces a new method to consider subjectivity and general context dependency in text analysis and uses as example the detection of emotions conveyed in text. The proposed method takes into account subjectivity using a computational version of the Framework Theory by Marvin Minsky (1974) leveraging on the Word2Vec approach to text vectorization by Mikolov et al. (2013), used to generate distributed representation of words based on the context where they appear. Our approach is based on three components: 1. a framework/'room' representing the point of view; 2. a benchmark representing the criteria for the analysis - in this case the emotion classification, from a study of human emotions by Robert Plutchik (1980); and 3. the document to be analyzed. By using similarity measure between words, we are able to extract the relative relevance of the elements in the benchmark - intensities of emotions in our case study - for the document to be analyzed. Our method provides a measure that take into account the point of view of the entity reading the document. This method could be applied to all the cases where evaluating subjectivity is relevant to understand the relative value or meaning of a text. Subjectivity can be not limited to human reactions, but it could be used to provide a text with an interpretation related to a given domain ("room"). To evaluate our method, we used a test case in the political domain.


title: Artificial life properties of directed interaction combinators vs. chemlambda

abstract: We provide a framework for experimentation at https://mbuliga.github.io/quinegraphs/ic-vs-chem.html#icvschem with two artificial chemistries: directed interaction combinators (dirIC, defined in section 2) and chemlambda. We are interested if these chemistries allow for artificial life behaviour: replication, metabolism and death. The main conclusion of these experiments is that graph rewrites systems which allow conflicting rewrites are better than those which don't, as concerns their artificial life properties. This is in contradiction with the search for good graph rewrite systems for decentralized computing, where non-conflicting graph rewrite systems are historically preferred. This continues the artificial chemistry experiments with chemlambda, lambda calculus or interaction combinators, available from the entry page at https://chemlambda.github.io/index.html and described in arXiv:2003.14332.


title: TOMA: Topological Map Abstraction for Reinforcement Learning

abstract: Animals are able to discover the topological map (graph) of surrounding environment, which will be used for navigation. Inspired by this biological phenomenon, researchers have recently proposed to generate graph representation for Markov decision process (MDP) and use such graphs for planning in reinforcement learning (RL). However, existing graph generation methods suffer from many drawbacks. One drawback is that existing methods do not learn an abstraction for graphs, which results in high memory and computation cost. This drawback also makes generated graph non-robust, which degrades the planning performance. Another drawback is that existing methods cannot be used for facilitating exploration which is important in RL. In this paper, we propose a new method, called topological map abstraction (TOMA), for graph generation. TOMA can generate an abstract graph representation for MDP, which costs much less memory and computation cost than existing methods. Furthermore, TOMA can be used for facilitating exploration. In particular, we propose planning to explore, in which TOMA is used to accelerate exploration by guiding the agent towards unexplored states. A novel experience replay module called vertex memory is also proposed to improve exploration performance. Experimental results show that TOMA can outperform existing methods to achieve the state-of-the-art performance.


title: Detection thresholds in very sparse matrix completion

abstract: Let $A$ be a rectangular matrix of size $m\times n$ and $A_1$ be the random matrix where each entry of $A$ is multiplied by an independent $\{0,1\}$-Bernoulli random variable with parameter $1/2$. This paper is about when, how and why the non-Hermitian eigen-spectra of the randomly induced asymmetric matrices $A_1 (A - A_1)^*$ and $(A-A_1)^*A_1$ captures more of the relevant information about the principal component structure of $A$ than via its SVD or the eigen-spectra of $A A^*$ and $A^* A$, respectively. Hint: the asymmetry inducing randomness breaks the echo-chamber effect that cripples the SVD. We illustrate the application of this striking phenomenon on the low-rank matrix completion problem for the setting where each entry is observed with probability $d/n$, including the very sparse regime where $d$ is of order $1$, where matrix completion via the SVD of $A$ fails or produces unreliable recovery. We determine an asymptotically exact, matrix-dependent, non-universal detection threshold above which reliable, statistically optimal matrix recovery using a new, universal data-driven matrix-completion algorithm is possible. Averaging the left and right eigenvectors provably improves the recovered matrix but not the detection threshold. We define another variant of this asymmetric procedure that bypasses the randomization step and has a detection threshold that is smaller by a constant factor but with a computational cost that is larger by a polynomial factor of the number of observed entries. Both detection thresholds shatter the seeming barrier due to the well-known information theoretical limit $d \asymp \log n$ for matrix completion found in the literature.


title: CO and HCN isotopologue ratios in the outflows of AGB stars

abstract: Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abundance ratios vary radially due to chemical reactions in the outflows of AGB stars and the effect of excitation and optical depth on the resulting line intensity ratios. We find that the circumstellar 12CO/13CO can deviate from its atmospheric value by up to 25-94% and 6-60% for C- and O-type CSEs, respectively. We show that variations of the intensity of the ISRF and the gas kinetic temperature can significantly influence the CO isotopologue ratio in the outer CSEs. On the contrary, the H12CN/H13CN ratio is stable for all tested mass-loss rates. The RT modeling shows that the integrated line intensity ratio of CO of different rotational transitions varies significantly for stars with intermediate mass-loss rates due to combined chemical and excitation effects. In contrast, the excitation conditions for the both HCN isotopologues are the same. We demonstrate the importance of using the isotopologue abundance profiles from chemical models as inputs to RT models in the interpretation of isotopologue observations. Previous studies of CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN as a tracer of C isotope ratios is affected by fewer complicating problems, provided one accounts corrections for high optical depths.


title: A case-study of multi-temperature coronal jets for emerging flux MHD models

abstract: Context: Hot coronal jets are a basic observed feature of the solar atmosphere whose physical origin is still being actively debated. Aims: We study six recurrent jets occurring in the active region NOAA 12644 on April 04, 2017. They are observed in all the hot filters of AIA as well as cool surges in IRIS slit-jaw high spatial and temporal resolution images. Methods: The AIA filters allow us to study the temperature and the emission measure of the jets using the filter ratio method. We study the pre-jet phases by analyzing the intensity oscillations at the base of the jets with the wavelet technique. Results: A fine co-alignment of the AIA and IRIS data shows that the jets are initiated at the top of a canopy-like, double-chambered structure with cool emission on one side and hot emission in the other. The hot jets are collimated in the hot temperature filters, have high velocities (around 250 km/s) and accompanied by the cool surges and ejected kernels both moving at about 45 km/s. In the pre-phase of the jets, at their base we find quasi-periodic intensity oscillations in phase with small ejections; they have a period between 2 and 6 minutes and are reminiscent of acoustic or MHD waves. Conclusions: This series of jets and surges provides a good case-study to test the 2D and 3D magnetohydrodynamic (MHD) models that result from magnetic flux emergence. The double-chambered structure found in the observations corresponds to the cold and hot loop regions found in the models beneath the current sheet that contains the reconnection site. The cool surge with kernels is comparable with the cool ejection and plasmoids that naturally appear in the models.


title: Automatic Estimation of Intelligibility Measure for Consonants in Speech

abstract: In this article, we provide a model to estimate a real-valued measure of the intelligibility of individual speech segments. We trained regression models based on Convolutional Neural Networks (CNN) for stop consonants \textipa{/p,t,k,b,d,g/} associated with vowel \textipa{/A/}, to estimate the corresponding Signal to Noise Ratio (SNR) at which the Consonant-Vowel (CV) sound becomes intelligible for Normal Hearing (NH) ears. The intelligibility measure for each sound is called SNR$_{90}$, and is defined to be the SNR level at which human participants are able to recognize the consonant at least 90\% correctly, on average, as determined in prior experiments with NH subjects. Performance of the CNN is compared to a baseline prediction based on automatic speech recognition (ASR), specifically, a constant offset subtracted from the SNR at which the ASR becomes capable of correctly labeling the consonant. Compared to baseline, our models were able to accurately estimate the SNR$_{90}$~intelligibility measure with less than 2 [dB$^2$] Mean Squared Error (MSE) on average, while the baseline ASR-defined measure computes SNR$_{90}$~with a variance of 5.2 to 26.6 [dB$^2$], depending on the consonant.


title: Defect filtering for thermal expansion induced dislocations in III-V lasers on silicon

abstract: Epitaxially integrated III-V semiconductor lasers for silicon photonics have the potential to dramatically transform information networks, but currently, dislocations limit performance and reliability even in defect tolerant InAs quantum dot (QD) based lasers. Despite being below critical thickness, QD layers in these devices contain previously unexplained misfit dislocations, which facilitate non-radiative recombination. We demonstrate here that these misfit dislocations form during post-growth cooldown due to the combined effects of (1) thermal-expansion mismatch between the III-V layers and silicon and (2) precipitate and alloy hardening in the active region. By incorporating an additional sub-critical thickness, indium-alloyed misfit dislocation trapping layer, we leverage these mechanical hardening effects to our advantage, successfully displacing 95% of misfit dislocations from the QD layer in model structures. Unlike conventional dislocation mitigation strategies, the trapping layer reduces neither the number of threading dislocations nor the number of misfit dislocations. It simply shifts the position of misfit dislocations away from the QD layer, reducing the defects' impact on luminescence. In full lasers, adding a misfit dislocation trapping layer both above and below the QD active region displaces misfit dislocations and substantially improves performance: we measure a twofold reduction in lasing threshold currents and a greater than threefold increase in output power. Our results suggest that devices employing both traditional threading dislocation reduction techniques and optimized misfit dislocation trapping layers may finally lead to fully integrated, commercially viable silicon-based photonic integrated circuits.


title: Shot noise, weak convergence and diffusion approximations

abstract: Shot noise processes have been extensively studied due to their mathematical properties and their relevance in several applications. Here, we consider nonnegative shot noise processes and prove their weak convergence to L\'evy-driven Ornstein-Uhlenbeck (OU), whose features depend on the underlying jump distributions. Among others, we obtain the OU-Gamma and OU-Inverse Gaussian processes, having gamma and inverse gaussian processes as background L\'evy processes, respectively. Then, we derive the necessary conditions guaranteeing the diffusion limit to a Gaussian OU process, show that they are not met unless allowing for negative jumps happening with probability going to zero, and quantify the error occurred when replacing the shot noise with the OU process and the non-Gaussian OU processes. The results offer a new class of models to be used instead of the commonly applied Gaussian OU processes to approximate synaptic input currents, membrane voltages or conductances modelled by shot noise in single neuron modelling.


title: Deep Learning for Wireless Communications

abstract: Existing communication systems exhibit inherent limitations in translating theory to practice when handling the complexity of optimization for emerging wireless applications with high degrees of freedom. Deep learning has a strong potential to overcome this challenge via data-driven solutions and improve the performance of wireless systems in utilizing limited spectrum resources. In this chapter, we first describe how deep learning is used to design an end-to-end communication system using autoencoders. This flexible design effectively captures channel impairments and optimizes transmitter and receiver operations jointly in single-antenna, multiple-antenna, and multiuser communications. Next, we present the benefits of deep learning in spectrum situation awareness ranging from channel modeling and estimation to signal detection and classification tasks. Deep learning improves the performance when the model-based methods fail. Finally, we discuss how deep learning applies to wireless communication security. In this context, adversarial machine learning provides novel means to launch and defend against wireless attacks. These applications demonstrate the power of deep learning in providing novel means to design, optimize, adapt, and secure wireless communications.


title: Soil searching by an artificial root

abstract: We model an artificial root which grows in the soil for underground prospecting. Its evolution is described by a controlled system of two integro-partial differential equations: one for the growth of the body and the other for the elongation of the tip. At any given time, the angular velocity of the root is obtained by solving a minimization problem with state constraints. We prove the existence of solutions to the evolution problem, up to the first time where a "breakdown configuration" is reached. Some numerical simulations are performed, to test the effectiveness of our feedback control algorithm.


title: Automated Extraction of Socio-political Events from News (AESPEN): Workshop and Shared Task Report

abstract: We describe our effort on automated extraction of socio-political events from news in the scope of a workshop and a shared task we organized at Language Resources and Evaluation Conference (LREC 2020). We believe the event extraction studies in computational linguistics and social and political sciences should further support each other in order to enable large scale socio-political event information collection across sources, countries, and languages. The event consists of regular research papers and a shared task, which is about event sentence coreference identification (ESCI), tracks. All submissions were reviewed by five members of the program committee. The workshop attracted research papers related to evaluation of machine learning methodologies, language resources, material conflict forecasting, and a shared task participation report in the scope of socio-political event information collection. It has shown us the volume and variety of both the data sources and event information collection approaches related to socio-political events and the need to fill the gap between automated text processing techniques and requirements of social and political sciences.


title: Glide Symmetry Protected Higher-Order Topological Insulators from Semimetals with butterfly-like Nodal Lines

abstract: Most topological insulators discovered today in spinful systems can be transformed from topological semimetals (TSMs) with vanishing bulk gap via introducing the spin-orbit coupling (SOC), which manifests the intrinsic links between the gapped TI phases and the gapless TSMs. Recently, we have proposed a new family of TSMs in time-reversal invariant {\it spinless} systems, which host butterfly-like nodal-lines (NLs) consisting of a pair of identical concentric intersecting coplanar ellipses (CICE). In this Communication, we unveil the intrinsic link between this exotic class of nodal-line semimetals (NLSMs) and a $\mathbb{Z}_{4}$ = 2 topological crystalline insulator (TCI), by including substantial SOC. We demonstrate that in three space groups ({\it i.e.} $Pbam$ (No.55), $P4/mbm$ (No.127) and $P4_2/mbc$ (No.135)), the TCI supports a fourfold Dirac fermion on the (001) surface protected by two glide symmetries, which originates from the intertwined drumhead surface states of the CICE NLs. The higher order topology is further demonstrated by the emergence of one-dimensional helical hinge states, indicating a new higher order topological insulator protected by a glide symmetry.


title: A time splitting method for the three-dimensional linear Pauli equation

abstract: We analyze a numerical method to solve the time-dependent linear Pauli equation in three space dimensions. The Pauli equation is a semi-relativistic generalization of the Schr\"odinger equation for 2-spinors which accounts both for magnetic fields and for spin, with the latter missing in preceding numerical work on the linear magnetic Schr\"odinger equation. We use a four operator splitting in time, prove stability and convergence of the method and derive error estimates as well as meshing strategies for the case of given time-independent electromagnetic potentials, thus providing a generalization of previous results for the magnetic Schr\"odinger equation.


title: Impact of the lattice on magnetic properties and possible spin nematicity in the S=1 triangular antiferromagnet NiGa$_2$S$_4$

abstract: NiGa$_2$S$_4$ is a triangular lattice S=1 system with strong two-dimensionality of the lattice, actively discussed as a candidate to host spin-nematic order brought about by strong quadrupole coupling. Using Raman scattering spectroscopy we identify a phonon of E$_g$ symmetry which can modulate magnetic exchange $J_1$ and produce quadrupole coupling. Additionally, our Raman scattering results demonstrate a loss of local inversion symmetry on cooling, which we associate with sulfur vacancies. This will lead to disordered Dzyaloshinskii-Moriya interactions, which can prevent long range magnetic order. Using magnetic Raman scattering response we identify 160~K as a temperature of an upturn of magnetic correlations. The temperature below 160~K, but above 50~K where antiferromagnetic magnetic start to increase, is a candidate for spin-nematic regime.


title: Propagation of Wave Packets in Dispersive Media

abstract: The propagation and refraction of waves in dispersive media are considered. A primary objective is to determine whether waves refract as monochromatic waves by Snell's law with phase velocity or as wave packets. The refraction of wave packets requires two refraction laws. The wave packets refract according to Snell's law with the geometric group velocity while the wavelets within the wave packets refract according to Snell's law with phase velocity. The tests were performed using directional power spectra from gravity water waves originating in the Caribbean Sea, Gulf of Mexico, Indian Ocean, and the Southern Ocean. In all cases the waves were measured with arrays of pressure sensors. It was found that the refraction of monochromatic waves led to backtracks from a measurement site that spread out for the different wave periods and did not go back to a common source. When the refraction was large the results were especially bad. When the data were used to test the refraction of wave packets very good results were obtained. For every storm considered the backtracks of the different wave periods went back together to the wave generation area. These results offer strong support for refracting waves in dispersive media as wave packets, not as monochromatic waves. This result for wave groups was also determined in a controlled wave tank experiment of gravity water waves where a critical angle was tested. The properties of propagating wave packets are also considered and illustrated. Reflection points and critical angles are found and explained. There can be major differences between the propagation of wave packets and monochromatic waves.


title: Android Malware Clustering using Community Detection on Android Packages Similarity Network

abstract: The daily amount of Android malicious applications (apps) targeting the app repositories is increasing, and their number is overwhelming the process of fingerprinting. To address this issue, we propose an enhanced Cypider framework, a set of techniques and tools aiming to perform a systematic detection of mobile malware by building a scalable and obfuscation resilient similarity network infrastructure of malicious apps. Our approach is based on our proposed concept, namely malicious community, in which we consider malicious instances that share common features are the most likely part of the same malware family. Using this concept, we presumably assume that multiple similar Android apps with different authors are most likely to be malicious. Specifically, Cypider leverages this assumption for the detection of variants of known malware families and zero-day malicious apps. Cypider applies community detection algorithms on the similarity network, which extracts sub-graphs considered as suspicious and possibly malicious communities. Furthermore, we propose a novel fingerprinting technique, namely community fingerprint, based on a one-class machine learning model for each malicious community. Besides, we proposed an enhanced Cypider framework, which requires less memory, x650, and less time to build the similarity network, x700, compared to the original version, without affecting the fingerprinting performance of the framework. We introduce a systematic approach to locate the best threshold on different feature content vectors, which simplifies the overall detection process.


title: Discrete Bessel functions and transform

abstract: We present a straightforward discretization of the Bessel functions $J_n(x)$ to discrete counterparts $B^{(N)}_n(x_m)$, of $N$ integer orders $n$ on $N$ integer points $x_m \equiv m$, that we call discrete Bessel functions. These are built from a Bessel integral generating function, restricting the Fourier transform over the circle to $N$ points. We show that the discrete Bessel functions satisfy several linear and quadratic relations, particularly Graf's product-displacement formulas, that are exact analogues of well-known relations between the continuous functions. It is noteworthy that these discrete Bessel functions approximate very closely the values of the continuous functions in ranges $n + |m| < N$. For fixed $N$, this provides an $N$-point transform between functions of order and of position,$f_n$ and $\widetilde{f}_m$, which is efficient for the Fourier analysis of finite decaying signals.


title: Three-Photon Absorption Spectra and Bandgap Scaling in Direct-Gap Semiconductors

abstract: This paper presents three-photon absorption (3PA) measurement results for nine direct-gap semiconductors, including full 3PA spectra for ZnSe, ZnS, and GaAs. These results, along with our theory of 3PA using an 8-band Kane model (4 bands with double spin degeneracy), help to explain the significant disagreements between experiments and theory in the literature to date. 3PA in the 8-band model exhibits quantum interference between the various possible pathways that is not observed in previous 2-band theories. We present measurements of degenerate 3PA coefficients in InSb, GaAs, CdTe, CdSe, ZnTe, CdS, ZnSe, ZnO, and ZnS. We examine bandgap, Eg, scaling using 2-band tunneling and perturbation theories that show agreement with the predicted Eg^-7 dependence; however, For those semiconductors for which we measured full 3PA spectra, we observe significant discrepancies with both 2-band theories. On the other hand, our 8-band model shows excellent agreement with the spectral data. We then use our 8-band theory to predict the 3PA spectra for 15 different semiconductors in their zincblende form. These results allow prediction and interpretation of the 3PA coefficients for various narrow to wide bandgap semiconductors.


title: LiSBOA: LiDAR Statistical Barnes Objective Analysis for optimal design of LiDAR scans and retrieval of wind statistics. Part I: Theoretical framework

abstract: A LiDAR Statistical Barnes Objective Analysis (LiSBOA) for optimal design of LiDAR scans and retrieval of the velocity statistical moments is proposed. The LiSBOA represents an adaptation of the classical Barnes scheme for the statistical analysis of unstructured experimental data in N-dimensional spaces and it is a suitable technique for the evaluation over a structured Cartesian grid of the statistics of scalar fields sampled through scanning LiDARs. The LiSBOA is validated and characterized via a Monte Carlo approach applied to a synthetic velocity field. This revisited theoretical framework for the Barnes objective analysis enables the formulation of guidelines for optimal design of LiDAR experiments and efficient application of the LiSBOA for the post-processing of LiDAR measurements. The optimal design of LiDAR scans is formulated as a two cost-function optimization problem including the minimization of the percentage of the measurement volume not sampled with adequate spatial resolution and the minimization of the error on the mean of the velocity field. The optimal design of the LiDAR scans also guides the selection of the smoothing parameter and the total number of iterations to use for the Barnes scheme.


title: The variable and non-variable X-ray absorbers in Compton-thin type-II Active Galactic Nuclei

abstract: We have conducted an extensive X-ray spectral variability study of a sample of 20 Compton-thin type II galaxies using broad band spectra from XMM-Newton, Chandra, and Suzaku. The aim is to study the variability of the neutral intrinsic X-ray obscuration along the line of sight and investigate the properties and location of the dominant component of the X-ray-obscuring gas. The observations are sensitive to absorption columns of $N_{\rm H} \sim 10^{20.5-24} {\rm cm^{-2}}$ of fully- and partially-covering neutral and/or lowly-ionized gas on timescales spanning days to well over a decade. We detected variability in the column density of the full-covering absorber in 7/20 sources, on timescales of months-years, indicating a component of compact-scale X-ray-obscuring gas lying along the line of sight of each of these objects. Our results imply that torus models incorporating clouds or overdense regions should account for line of sight column densities as low as $\sim$ a few $\times 10^{21}$ cm$^{-2}$. However, 13/20 sources yielded no detection of significant variability in the full-covering obscurer, with upper limits to ${\Delta}N_{\rm H}$ spanning $10^{21-23}$ cm$^{-2}$. The dominant absorbing media in these systems could be distant, such as kpc-scale dusty structures associated with the host galaxy, or a homogeneous medium along the line of sight. Thus, we find that overall, strong variability in full-covering obscurers is not highly prevalent in Compton-thin type IIs, at least for our sample, in contrast to previous results in the literature. Finally, 11/20 sources required a partial-covering, obscuring component in all or some of their observations, consistent with clumpy near-Compton-thick compact-scale gas.


title: CP Violation in $B^\pm\to \rho^0\pi^\pm$ and $B^\pm\to \sigma\pi^\pm$ Decays

abstract: The decay amplitude of $B^+\to \pi^+\pi^-\pi^+$ in the Dalitz plot has been analyzed by the LHCb using three different approaches for the $S$-wave component. It was found that the mode with $\sigma$ (or $f_0(500)$) exhibited a CP asymmetry of 15\% in the isobar model, whereas the $f_2(1270)$ mode had a 40\% asymmetry. On the contrary, CP asymmetry for the dominant quasi-two-body decay $B^-\to\rho^0\pi^-$ was found to be consistent with zero in all three approaches, while all the existing theoretical predictions lead to a negative CP asymmetry ranging from $-7\%$ to $-45\%$. We show that the nearly vanishing CP violation in $B^-\to\rho^0\pi^-$ is understandable in the framework of QCD factorization (QCDF). It arises from the $1/m_b$ power corrections to the penguin amplitudes due to penguin annihilations and to the color-suppressed tree amplitude due to hard spectator interactions. Penguin annihilation and hard spectator interactions contribute destructively to $A_{CP}(B^-\to\rho^0\pi^-)$ to render it consistent with zero. The branching fraction and CP asymmetry in $B^-\to\sigma/f_0(500)\pi^-$ are investigated in QCDF with results in agreement with experiment.


title: A fast sparse spectral method for nonlinear integro-differential Volterra equations with general kernels

abstract: We present a sparse spectral method for nonlinear integro-differential Volterra equations based on the Volterra operator's banded sparsity structure when acting on specific Jacobi polynomial bases. The method is not restricted to convolution-type kernels of the form $K(x,y)=K(x-y)$ but instead works for general kernels at competitive speeds and with exponential convergence. We provide various numerical experiments on problems with or without known analytic solutions and comparisons with other methods.


title: Universal coherence protection in a solid-state spin qubit

abstract: Decoherence largely limits the physical realization of qubits and its mitigation is critical to quantum science. Here, we construct a robust qubit embedded in a decoherence-protected subspace, obtained by hybridizing an applied microwave drive with the ground-state electron spin of a silicon carbide divacancy defect. The qubit is protected from magnetic, electric, and temperature fluctuations, which account for nearly all relevant decoherence channels in the solid state. This culminates in an increase of the qubit's inhomogeneous dephasing time by over four orders of magnitude (to > 22 milliseconds), while its Hahn-echo coherence time approaches 64 milliseconds. Requiring few key platform-independent components, this result suggests that substantial coherence improvements can be achieved in a wide selection of quantum architectures.


title: Stochastic Learning for Sparse Discrete Markov Random Fields with Controlled Gradient Approximation Error

abstract: We study the $L_1$-regularized maximum likelihood estimator/estimation (MLE) problem for discrete Markov random fields (MRFs), where efficient and scalable learning requires both sparse regularization and approximate inference. To address these challenges, we consider a stochastic learning framework called stochastic proximal gradient (SPG; Honorio 2012a, Atchade et al. 2014,Miasojedow and Rejchel 2016). SPG is an inexact proximal gradient algorithm [Schmidtet al., 2011], whose inexactness stems from the stochastic oracle (Gibbs sampling) for gradient approximation - exact gradient evaluation is infeasible in general due to the NP-hard inference problem for discrete MRFs [Koller and Friedman, 2009]. Theoretically, we provide novel verifiable bounds to inspect and control the quality of gradient approximation. Empirically, we propose the tighten asymptotically (TAY) learning strategy based on the verifiable bounds to boost the performance of SPG.


title: Parameterization method for state-dependent delay perturbation of an ordinary differential equation

abstract: We consider state-dependent delay equations (SDDE) obtained by adding delays to a planar ordinary differential equation with a limit cycle. These situations appear in models of several physical processes, where small delay effects are added. Even if the delays are small, they are very singular perturbations since the natural phase space of an SDDE is an infinite dimensional space. We show that the SDDE admits solutions which resemble the solutions of the ODE. That is, there exist a periodic solution and a two parameter family of solutions whose evolution converges to the periodic solution. Even if the phase space of the SDDE is naturally a space of functions, we show that there are initial values which lead to solutions similar to that of the ODE. The method of proof bypasses the theory of existence, uniqueness, dependence on parameters of SDDE. We consider the class of functions of time that have a well defined behavior (e.g. periodic, or asymptotic to periodic) and derive a functional equation which imposes that they are solutions of the SDDE. These functional equations are studied using methods of functional analysis. We provide a result in "a posteriori" format: Given an approximate solution of the functional equation, which has some good condition numbers, we prove that there is true solution close to the approximate one. Thus, we can use the result to validate the results of numerical computations. The method of proof leads also to practical algorithms. In a companion paper, we present the implementation details and representative results. One feature of the method presented here is that it allows to obtain smooth dependence on parameters for the periodic solutions and their slow stable manifolds without studying the smoothness of the flow (which seems to be problematic for SDDEs, for now the optimal result on smoothness of the flow is $C^1$).


title: Constructions of highest weight modules of double Ringel-Hall algebras via functions

abstract: In [19], Zheng studied the bounded derived categories of constructible $\bar{\mathbb{Q}}_l$-sheaves on some algebraic stacks consisting of the representations of a enlarged quiver and categorified the integrable highest weight modules of the corresponding quantum group by using these categories. In this paper, we shall generalize Zheng's work to highest weight modules of a subalgebra of the double Ringel-Hall algebra associated to a quiver in a functional version.