title: Monitoring fetal electroencephalogram intrapartum: a systematic literature review

abstract: Background: Studies about the feasibility of monitoring fetal electroencephalogram (fEEG) during labor began in the early 1940s. By the 1970s, clear diagnostic and prognostic benefits from intrapartum fEEG monitoring were reported, but until today, this monitoring technology has remained a curiosity. Objectives: Our goal was to review the studies reporting the use of fEEG including the insights from interpreting fEEG patterns in response to uterine contractions during labor. We also used the most relevant information gathered from clinical studies to provide recommendations for enrollment in the unique environment of a labor and delivery unit. Data sources: PubMed. Eligibility criteria: The search strategy was: ("fetus"[MeSH Terms] OR "fetus"[All Fields] OR "fetal"[All Fields]) AND ("electroencephalography"[MeSH Terms] OR "electroencephalography"[All Fields] OR "eeg"[All Fields]) AND (Clinical Trial[ptyp] AND "humans"[MeSH Terms]). Because the landscape of fEEG research has been international, we included studies in English, French, German, and Russian. Results: From 256 screened studies, 40 studies were ultimately included in the qualitative analysis. We summarize and report features of fEEG which clearly show its potential to act as a direct biomarker of fetal brain health during delivery, ancillary to fetal heart rate monitoring. However, clinical prospective studies are needed to further establish the utility of fEEG monitoring intrapartum. We identified clinical study designs likely to succeed in bringing this intrapartum monitoring modality to the bedside. Limitations: Despite 80 years of studies in clinical cohorts and animal models, the field of research on intrapartum fEEG is still nascent and shows great promise to augment the currently practiced electronic fetal monitoring.


title: Language Representation Models for Fine-Grained Sentiment Classification

abstract: Sentiment classification is a quickly advancing field of study with applications in almost any field. While various models and datasets have shown high accuracy inthe task of binary classification, the task of fine-grained sentiment classification is still an area with room for significant improvement. Analyzing the SST-5 dataset,previous work by Munikar et al. (2019) showed that the embedding tool BERT allowed a simple model to achieve state-of-the-art accuracy. Since that paper, several BERT alternatives have been published, with three primary ones being AlBERT (Lan et al., 2019), DistilBERT (Sanh et al. 2019), and RoBERTa (Liu etal. 2019). While these models report some improvement over BERT on the popular benchmarks GLUE, SQuAD, and RACE, they have not been applied to the fine-grained classification task. In this paper, we examine whether the improvements hold true when applied to a novel task, by replicating the BERT model from Munikar et al., and swapping the embedding layer to the alternative models. Over the experiments, we found that AlBERT suffers significantly more accuracy loss than reported on other tasks, DistilBERT has accuracy loss similar to their reported loss on other tasks while being the fastest model to train, and RoBERTa reaches anew state-of-the-art accuracy for prediction on the SST-5 root level (60.2%).


title: Further $q$-series identities and conjectures relating false theta functions and characters

abstract: In this short note, a companion of [20], we discuss several families of $q$-series identities in connection to false and mock theta functions, characters of modules of vertex algebras, and "sum of tails".


title: Dynamic Coupling and Market Instability

abstract: We examine dynamic coupling and feedback effects between High Frequency Traders (HFTs) and how they can destabilize markets. We develop a general framework for modelling dynamic interaction based on recurrence relations, and use this to show how unexpected latency and feedback can trigger oscillatory instability between HFT market makers with inventory constraints. Our analysis suggests that the modelled instability is an unintentional emergent behaviour of the market that does not depend on the complexity of HFT strategies - even apparently stable strategies are vulnerable. Feedback instability can lead to substantial movements in market prices such as price spikes and crashes.


title: Assessing variable activity for Bayesian regression trees

abstract: Bayesian Additive Regression Trees (BART) are non-parametric models that can capture complex exogenous variable effects. In any regression problem, it is often of interest to learn which variables are most active. Variable activity in BART is usually measured by counting the number of times a tree splits for each variable. Such one-way counts have the advantage of fast computations. Despite their convenience, one-way counts have several issues. They are statistically unjustified, cannot distinguish between main effects and interaction effects, and become inflated when measuring interaction effects. An alternative method well-established in the literature is Sobol' indices, a variance-based global sensitivity analysis technique. However, these indices often require Monte Carlo integration, which can be computationally expensive. This paper provides analytic expressions for Sobol' indices for BART posterior samples. These expressions are easy to interpret and are computationally feasible. Furthermore, we will show a fascinating connection between first-order (main-effects) Sobol' indices and one-way counts. We also introduce a novel ranking method, and use this to demonstrate that the proposed indices preserve the Sobol'-based rank order of variable importance. Finally, we compare these methods using analytic test functions and the En-ROADS climate impacts simulator.


title: On two-weight codes

abstract: We consider $q$-ary (linear and nonlinear) block codes with exactly two distances: $d$ and $d+\delta$. Several combinatorial constructions of optimal such codes are given. In the linear (but not necessary projective) case, we prove that under certain conditions the existence of such linear $2$-weight code with $\delta > 1$ implies the following equality of great common divisors: $(d,q) = (\delta,q)$. Upper bounds for the maximum cardinality of such codes are derived by linear programming and from few-distance spherical codes. Tables of lower and upper bounds for small $q = 2,3,4$ and $q\,n < 50$ are presented.


title: Exclusive four-pion photoproduction in ultra-peripheral heavy-ion collisions at RHIC and LHC energies

abstract: We study the photoproduction of exclusive $2\pi^+2\pi^-$ mesons in ultra-peripheral heavy-ion collisions at RHIC and LHC energies. Predictions in photon-nucleus interactions are calculated for various resonances at central and forward rapidities. The recent H1 preliminary data are utilized to improve the description of the poorly known $\gamma p \to 4\pi^\pm p$ process. We present the comparisons of our results to the available STAR data at RHIC, and made predictions for LHC energies.


title: Revisiting Parameter Sharing in Multi-Agent Deep Reinforcement Learning

abstract: Parameter sharing, where each agent independently learns a policy with fully shared parameters between all policies, is a popular baseline method for multi-agent deep reinforcement learning. Unfortunately, since all agents share the same policy network, they cannot learn different policies or tasks. This issue has been circumvented experimentally by adding an agent-specific indicator signal to observations, which we term "agent indication". Agent indication is limited, however, in that without modification it does not allow parameter sharing to be applied to environments where the action spaces and/or observation spaces are heterogeneous. This work formalizes the notion of agent indication and proves that it enables convergence to optimal policies for the first time. Next, we formally introduce methods to extend parameter sharing to learning in heterogeneous observation and action spaces, and prove that these methods allow for convergence to optimal policies. Finally, we experimentally confirm that the methods we introduce function empirically, and conduct a wide array of experiments studying the empirical efficacy of many different agent indication schemes for image based observation spaces.


title: Demonstration of site-selective angular-resolved absorption spectroscopy of the $^{4}I_{15/2} \rightarrow $$^{4}I_{13/2}$ erbium transition in the monoclinic crystal Y$_2$SiO$_5$

abstract: We study the angular dependence in polarized light of the optical absorption for the Er$^{3+}$ transition $^{4}I_{15/2} \rightarrow $$^{4}I_{13/2}$ in Y$_2$SiO$_5$, revealing thus the associated anisotropy and the orientation of the related absorption principal directions in the dielectric plan perpendicular to the monoclinic axis b. The measurements are performed at low temperature. This allows us to isolate the lowest crystal field levels in the ground and excited states. We spectrally resolve and independently characterize the two yttrium substitution sites in the Y$_2$SiO$_5$ matrix. The absorption tensor components cannot be unambiguously determined yet while only considering the investigated dielectric plane. Still, measurements remarkably demonstrate that this transition of interest well resolved at low temperature is not only a magnetic-dipole allowed transition but indeed a hybrid electric-magnetic transition.


title: Blow-up rates for higher-order semilinear parabolic equations with nonlinear memory term

abstract: In this paper, we establish blow-up rates for higher-order semilinear parabolic equations with nonlocal in time nonlinearity with no positive assumption on the solution. We also give Liouville-type theorem for higher-order semilinear parabolic equation with infinite memory nonlinear term which plays the main tools to prove our blow-up rate result. Finally, we study the well-posedness of mild solutions.


title: Distributed algorithms for covering, packing and maximum weighted matching

abstract: This paper gives poly-logarithmic-round, distributed D-approximation algorithms for covering problems with submodular cost and monotone covering constraints (Submodular-cost Covering). The approximation ratio D is the maximum number of variables in any constraint. Special cases include Covering Mixed Integer Linear Programs (CMIP), and Weighted Vertex Cover (with D=2). Via duality, the paper also gives poly-logarithmic-round, distributed D-approximation algorithms for Fractional Packing linear programs (where D is the maximum number of constraints in which any variable occurs), and for Max Weighted c-Matching in hypergraphs (where D is the maximum size of any of the hyperedges; for graphs D=2). The paper also gives parallel (RNC) 2-approximation algorithms for CMIP with two variables per constraint and Weighted Vertex Cover. The algorithms are randomized. All of the approximation ratios exactly match those of comparable centralized algorithms.


title: Simultaneous Diagonalization of Incomplete Matrices and Applications

abstract: We consider the problem of recovering the entries of diagonal matrices $\{U_a\}_a$ for $a = 1,\ldots,t$ from multiple "incomplete" samples $\{W_a\}_a$ of the form $W_a=PU_aQ$, where $P$ and $Q$ are unknown matrices of low rank. We devise practical algorithms for this problem depending on the ranks of $P$ and $Q$. This problem finds its motivation in cryptanalysis: we show how to significantly improve previous algorithms for solving the approximate common divisor problem and breaking CLT13 cryptographic multilinear maps.


title: Explaining Neural Networks by Decoding Layer Activations

abstract: We present a `CLAssifier-DECoder' architecture (\emph{ClaDec}) which facilitates the comprehension of the output of an arbitrary layer in a neural network (NN). It uses a decoder to transform the non-interpretable representation of the given layer to a representation that is more similar to the domain a human is familiar with. In an image recognition problem, one can recognize what information is represented by a layer by contrasting reconstructed images of \emph{ClaDec} with those of a conventional auto-encoder(AE) serving as reference. We also extend \emph{ClaDec} to allow the trade-off between human interpretability and fidelity. We evaluate our approach for image classification using Convolutional NNs. We show that reconstructed visualizations using encodings from a classifier capture more relevant information for classification than conventional AEs. Relevant code is available at \url{https://github.com/JohnTailor/ClaDec}


title: Response and Uncertainty of the Parabolic Variance PVAR to Non-Integer Exponents of the Power Law

abstract: Oscillator fluctuations are described as the phase or frequency noise spectrum, or in terms of a wavelet variance as a function of the measurement time. The spectrum is generally approximated by the `power law,' i.e., a Laurent polynomial with integer exponents of the frequency. This article extends the domain of application of PVAR, a wavelet variance which uses the linear regression on phase data to estimate the frequency, and called `parabolic' because such regression is equivalent to a parabolic-shaped weight function applied to frequency fluctuations. In turn, PVAR is relevant in that it improves on the widely-used Modified Allan variance (MVAR) enabling the detection of the same noise processes at the same confidence level in a shorter measurement time. More specifically, we provide (i) the analytical expression of the response of the PVAR to the frequency-noise spectrum in the general case of non-integer exponents of the frequency, and (ii) a useful approximate expression of the statistical uncertainty.


title: GraFS: Graph Analytics Fusion and Synthesis

abstract: Graph analytics elicits insights from large graphs to inform critical decisions for business, safety and security. Several large-scale graph processing frameworks feature efficient runtime systems; however, they often provide programming models that are low-level and subtly different from each other. Therefore, end users can find implementation and specially optimization of graph analytics time-consuming and error-prone. This paper regards the abstract interface of the graph processing frameworks as the instruction set for graph analytics, and presents Grafs, a high-level declarative specification language for graph analytics and a synthesizer that automatically generates efficient code for five high-performance graph processing frameworks. It features novel semantics-preserving fusion transformations that optimize the specifications and reduce them to three primitives: reduction over paths, mapping over vertices and reduction over vertices. Reductions over paths are commonly calculated based on push or pull models that iteratively apply kernel functions at the vertices. This paper presents conditions, parametric in terms of the kernel functions, for the correctness and termination of the iterative models, and uses these conditions as specifications to automatically synthesize the kernel functions. Experimental results show that the generated code matches or outperforms hand-optimized code, and that fusion accelerates execution.


title: AGN anisotropic radiative feedback set by black hole spin

abstract: We consider the impact of anisotropic radiation on the active galactic nucleus (AGN) radiative dusty feedback. The radiation pattern originating from the accretion disc is determined by the central black hole (BH) spin. Here we analyse how such BH spin-induced angular dependence affects the dynamics and energetics of the radiation pressure-driven outflows, as well as AGN obscuration and BH accretion. In addition, we explore the effect of a spatially varying dust-to-gas ratio on the outflow propagation. We obtain two distinct trends for high-spin and low-spin objects, providing a direct connection between anisotropic feedback and BH spin. In the case of maximum spin, powerful quasi-spherical outflows can propagate on large scales, at all inclination angles with fairly uniform energetics. In contrast, in the case of zero spin, only weaker bipolar outflows can be driven in the polar directions. As a result, high BH spins can efficiently clear out the obscuring gas from most directions, whereas low BH spins can only remove dusty gas from the polar regions, hence also determining the overall AGN obscuration geometry. Due to such anisotropic feedback, high BH spins can prevent accretion of gas from most directions (except in the equatorial plane), while low BH spins allow inflows to proceed from a wider range of directions. This may have important implications for the BH growth in the early Universe. Anisotropic radiative dusty feedback, ruled by the BH spin, may thus play a major role in shaping AGN evolution over cosmic time.


title: On motifs in colored graphs

abstract: One of the most important concepts in biological network analysis is that of network motifs, which are patterns of interconnections that occur in a given network at a frequency higher than expected in a random network. In this work we are interested in searching and inferring network motifs in a class of biological networks that can be represented by vertex-colored graphs. We show the computational complexity for many problems related to colorful topological motifs and present efficient algorithms for special cases. We also present a probabilistic strategy to detect highly frequent motifs in vertex-colored graphs. Experiments on real data sets show that our algorithms are very competitive both in efficiency and in quality of the solutions.


title: Towards AI Forensics: Did the Artificial Intelligence System Do It?

abstract: Artificial intelligence (AI) makes decisions impacting our daily lives in an increasingly autonomous manner. Their actions might cause accidents, harm, or, more generally, violate regulations. Determining whether an AI caused a specific event and, if so, what triggered the AI's action, are key forensic questions. We provide a conceptualization of the problems and strategies for forensic investigation. We focus on AI that is potentially ``malicious by design'' and grey box analysis. Our evaluation using convolutional neural networks illustrates challenges and ideas for identifying malicious AI.


title: On the convergence of Kac-Moody Eisenstein series

abstract: Let $G$ be a representation-theoretic Kac--Moody group associated to a nonsingular symmetrizable generalized Cartan matrix. We first consider Kac-Moody analogs of Borel Eisenstein series (induced from quasicharacters on the Borel), and prove they converge almost everywhere inside the Tits cone for arbitrary spectral parameters in the Godement range. We then use this result to show the full absolute convergence everywhere inside the Tits cone (again for spectral parameters in the Godement range) for a class of Kac-Moody groups satisfying a certain combinatorial property, in particular for rank-2 hyperbolic groups.


title: Magnon antibunching in a nanomagnet

abstract: We investigate the correlations of magnons inside a nanomagnet and identify a regime of parameters where the magnons become antibunched, i.e., where there is a large probability for occupation of the single-magnon state. This antibunched state is very different from magnons at thermal equilibrium and microwave-driven coherent magnons. We further obtain the steady state analytically and describe the magnon dynamics numerically, and ascertain the stability of such antibunched magnons over a large window of magnetic anisotropy, damping and temperature. This means that the antibunched magnon state is feasible in a wide class of low-damping magnetic nanoparticles. To detect this quantum effect, we propose to transfer the quantum information of magnons to photons by magnon-photon coupling and then measure the correlations of photons to retrieve the magnon correlations. Our findings may provide a promising platform to study quantum-classical transitions and for designing a single magnon source.


title: Looking back to lower-level information in few-shot learning

abstract: Humans are capable of learning new concepts from small numbers of examples. In contrast, supervised deep learning models usually lack the ability to extract reliable predictive rules from limited data scenarios when attempting to classify new examples. This challenging scenario is commonly known as few-shot learning. Few-shot learning has garnered increased attention in recent years due to its significance for many real-world problems. Recently, new methods relying on meta-learning paradigms combined with graph-based structures, which model the relationship between examples, have shown promising results on a variety of few-shot classification tasks. However, existing work on few-shot learning is only focused on the feature embeddings produced by the last layer of the neural network. In this work, we propose the utilization of lower-level, supporting information, namely the feature embeddings of the hidden neural network layers, to improve classifier accuracy. Based on a graph-based meta-learning framework, we develop a method called Looking-Back, where such lower-level information is used to construct additional graphs for label propagation in limited data settings. Our experiments on two popular few-shot learning datasets, miniImageNet and tieredImageNet, show that our method can utilize the lower-level information in the network to improve state-of-the-art classification performance.


title: PNKH-B: A Projected Newton-Krylov Method for Large-Scale Bound-Constrained Optimization

abstract: We present PNKH-B, a projected Newton-Krylov method for iteratively solving large-scale optimization problems with bound constraints. PNKH-B is geared toward situations in which function and gradient evaluations are expensive, and the (approximate) Hessian is only available through matrix-vector products. This is commonly the case in large-scale parameter estimation, machine learning, and image processing. In each iteration, PNKH-B uses a low-rank approximation of the (approximate) Hessian to determine the search direction and construct the metric used in a projected line search. The key feature of the metric is its consistency with the low-rank approximation of the Hessian on the Krylov subspace. This renders PNKH-B similar to a projected variable metric method. We present an interior point method to solve the quadratic projection problem efficiently. Since the interior point method effectively exploits the low-rank structure, its computational cost only scales linearly with respect to the number of variables, and it only adds negligible computational time. We also experiment with variants of PNKH-B that incorporate estimates of the active set into the Hessian approximation. We prove the global convergence to a stationary point under standard assumptions. Using three numerical experiments motivated by parameter estimation, machine learning, and image reconstruction, we show that the consistent use of the Hessian metric in PNKH-B leads to fast convergence, particularly in the first few iterations. We provide our MATLAB implementation at https://github.com/EmoryMLIP/PNKH-B.


title: Hypergeometric L-functions in average polynomial time

abstract: We describe an algorithm for computing, for all primes $p \leq X$, the mod-$p$ reduction of the trace of Frobenius at $p$ of a fixed hypergeometric motive in time quasilinear in $X$. This combines the Beukers--Cohen--Mellit trace formula with average polynomial time techniques of Harvey et al.


title: Hyper-Kamiokande

abstract: A next generation water Cherenkov detector Hyper-Kamiokande to be built in Japan is described. The main goals of this project include a sensitive measurement of CP violation in neutrino oscillations, a search for proton decay and study of solar, atmospherics and astrophysical neutrinos. Key features of the Hyper-Kamiokande detector are described. The main emphasis is put on large photosensors. The recent progress in the development of near neutrino detectors is also presented.


title: Finite Quantum Instruments

abstract: This article considers quantum systems described by a finite-dimensional complex Hilbert space $H$. We first define the concept of a finite observable on $H$. We then discuss ways of combining observables in terms of convex combinations, post-processing and sequential products. We also define complementary and coexistent observables. We then introduce finite instruments and their related compatible observables. The previous combinations and relations for observables are extended to instruments and their properties are compared. We present four types of instruments; namely, identity, trivial, L\"uders and Kraus instruments. These types are used to illustrate different ways that instruments can act. We next consider joint probabilities for observables and instruments. The article concludes with a discussion of measurement models and the instruments they measure.


title: Segmentation of the Myocardium on Late-Gadolinium Enhanced MRI based on 2.5 D Residual Squeeze and Excitation Deep Learning Model

abstract: Cardiac left ventricular (LV) segmentation from short-axis MRI acquired 10 minutes after the injection of a contrast agent (LGE-MRI) is a necessary step in the processing allowing the identification and diagnosis of cardiac diseases such as myocardial infarction. However, this segmentation is challenging due to high variability across subjects and the potential lack of contrast between structures. Then, the main objective of this work is to develop an accurate automatic segmentation method based on deep learning models for the myocardial borders on LGE-MRI. To this end, 2.5 D residual neural network integrated with a squeeze and excitation blocks in encoder side with specialized convolutional has been proposed. Late fusion has been used to merge the output of the best trained proposed models from a different set of hyperparameters. A total number of 320 exams (with a mean number of 6 slices per exam) were used for training and 28 exams used for testing. The performance analysis of the proposed ensemble model in the basal and middle slices was similar as compared to intra-observer study and slightly lower at apical slices. The overall Dice score was 82.01% by our proposed method as compared to Dice score of 83.22% obtained from the intra observer study. The proposed model could be used for the automatic segmentation of myocardial border that is a very important step for accurate quantification of no-reflow, myocardial infarction, myocarditis, and hypertrophic cardiomyopathy, among others.


title: A note on the wave equation on black hole spacetimes with small non-decaying first order terms

abstract: We present an elementary physical space argument to establish local integrated decay estimates for the perturbed wave equation $\Box_g \phi = \epsilon \beta^a \partial_a \phi$ on the exterior of the Schwarzschild geometry $(\mathcal{M},g)$. Here $\beta$ is a regular vectorfield on $\mathcal{M}$ decaying suitably in space but not necessarily in time. The proof is formulated to cover also perturbations of the Regge--Wheeler equation.


title: Algorithmic approaches to selecting control clones in DNA array hybridization experiments

abstract: We study the problem of selecting control clones in DNA array hybridization experiments. The problem arises in the OFRG method for analyzing microbial communities. The OFRG method performs classification of rRNA gene clones using binary fingerprints created from a series of hybridization experiments, where each experiment consists of hybridizing a collection of arrayed clones with a single oligonucleotide probe. This experiment produces analog signals, one for each clone, which then need to be classified, that is, converted into binary values 1 and 0 that represent hybridization and non-hybridization events. In addition to the sample rRNA gene clones, the array contains a number of control clones needed to calibrate the classification procedure of the hybridization signals. These control clones must be selected with care to optimize the classification process. We formulate this as a combinatorial optimization problem called Balanced Covering. We prove that the problem is NP-hard, and we show some results on hardness of approximation. We propose approximation algorithms based on randomized rounding and we show that, with high probability, our algorithms approximate well the optimum solution. The experimental results confirm that the algorithms find high quality control clones. The algorithms have been implemented and are publicly available as part of the software package called CloneTools.


title: Exact spectral function of a Tonks-Girardeau gas in a lattice

abstract: The single-particle spectral function of a strongly correlated system is an essential ingredient to describe its dynamics and transport properties. We develop a general method to calculate the exact spectral function of a strongly interacting one-dimensional Bose gas in the Tonks-Girardeau regime, valid for any type of confining potential, and apply it to bosons on a lattice to obtain the full spectral function, at all energy and momentum scales. We find that it displays three main singularity lines. The first two can be identified as the analogs of Lieb-I and Lieb-II modes of a uniform fluid; the third one, instead, is specifically due to the presence of the lattice. We show that the spectral function displays a power-law behaviour close to the Lieb-I and Lieb-II singularities, as predicted by the non-linear Luttinger liquid description, and obtain the exact exponents. In particular, the Lieb-II mode shows a divergence in the spectral function, differently from what happens in the dynamical structure factor, thus providing a route to probe it in experiments with ultracold atoms.


title: Spicule jets in the solar atmosphere modeled with resistive MHD and thermal conduction

abstract: Using numerical simulations, we study the effects of magnetic resistivity and thermal conductivity in the dynamics and properties of solar jets with characteristics of Type II spicules and cool coronal jets. The dynamic evolution of the jets is governed by the resistive MHD equations with thermal conduction along the magnetic field lines on a 2.5D slice. The magnetic field configuration consists of two symmetric neighboring loops with opposite polarity, used to support reconnection and followed by the plasma jet formation. In total 10 simulations were carried out with different values of resistivity and thermal conductivity, that produce jets with different morphological and thermal properties we quantify. We find that an increase in magnetic resistivity does not produce significant effects on the morphology, velocity and temperature of the jets. However, thermal conductivity affects both temperature and morphology of the jets. In particular, thermal conductivity causes jets to reach greater heights and increases the temperature of the jet-apex. Also, heat flux maps indicate the jet-apex and corona interchange energy more efficiently than the body of the jet. These results could potentially open a new avenue for plasma diagnostics in the Sun's atmosphere.


title: In-depth Real-World Evaluation of NB-IoT Module Energy Consumption

abstract: Narrowband Internet of Things (NB-IoT) is a recent addition to the 3GPP standards offering low power wide area networking (LP-WAN) for a massive amount of IoT devices communicating at low data rates in the licensed bands. As the number of deployed NB-IoT devices worldwide increases, the question of energy-efficient real-world NB-IoT device and network operation that would maximize the device battery lifetime becomes crucial for NB-IoT service popularity and adoption. In this paper, we present a detailed energy consumption of an NB-IoT module obtained from a custom-designed NB-IoT platform from which a large amount of high time-resolution data is collected. We perform a detailed energy consumption analysis of each NB-IoT data transfer phase and discuss both the device and the network-side configurations that may affect the module energy consumption in each of the phases. Preliminary results and conclusions are put forth along with the discussion of ongoing and future study plans.


title: The Orbit Type Stratification OF the Moduli Space of Higgs Bundles

abstract: The moduli space of Higgs bundles can be constructed as a quotient of an infinite-dimensional space and hence admits an orbit type decomposition. In this paper, we show that the orbit type decomposition is a complex Whitney stratification such that each stratum is a complex symplectic submanifold and hence admits a complex Poisson bracket. Moreover, these Poisson brackets glue to a Poisson bracket on the structure sheaf of the moduli space so that the moduli space is a stratified complex symplectic space


title: Group testing with nested pools

abstract: In order to identify the infected individuals of a population, their samples are divided in equally sized groups called pools and a single laboratory test is applied to each pool. Individuals whose samples belong to pools that test negative are declared healthy, while each pool that tests positive is divided into smaller, equally sized pools which are tested in the next stage. In the $(k+1)$-th stage all remaining samples are tested. If $p<1-3^{-1/3}$, we minimize the expected number of tests per individual as a function of the number $k+1$ of stages, and of the pool sizes in the first $k$ stages. We show that for each $p\in (0, 1-3^{-1/3})$ the optimal choice is one of four possible schemes, which are explicitly described. We conjecture that for each $p$, the optimal choice is one of the two sequences of pool sizes $(3^k\text{ or }3^{k-1}4,3^{k-1},\dots,3^2,3 )$, with a precise description of the range of $p$'s where each is optimal. The conjecture is supported by overwhelming numerical evidence for $p>2^{-51}$. We also show that the cost of the best among the schemes $(3^k,\dots,3)$ is of order $O\big(p\log(1/p)\big)$, comparable to the information theoretical lower bound $p\log_2(1/p)+(1-p)\log_2(1/(1-p))$, the entropy of a Bernoulli$(p)$ random variable.

categories: math.ST physics.data-an q-bio.QM stat.AP stat.TH


title: Stability and gravitational collapse of neutron stars with realistic equations of state

abstract: We discuss the stability and construct dynamical configurations describing the gravitational collapse of unstable neutron stars with realistic equations of state compatible with the recent LIGO-Virgo constraints. Unlike other works that consider the collapse of a stellar configuration without a priori knowledge if it is stable or unstable, we first perform a complete analysis on stellar stability for such equations of state. Negative values of the squared frequency of the fundamental mode indicate us radial instability with respect to the collapse of the unstable star to a black hole. We find numerical solutions corresponding to the temporal and radial behavior during the evolution of the collapse for certain relevant physical quantities such as mass, luminosity, energy density, pressure, heat flow, temperature and quantities that describe bulk viscous processes. Our results show that the equation of state undergoes abrupt changes close to the moment of event horizon formation as a consequence of dissipative effects. During the collapse process all energy conditions are respected, which implies that our model is physically acceptable.


title: Using Hypotheses as a Debugging Aid

abstract: As developers debug, developers formulate hypotheses about the cause of the defect and gather evidence to test these hypotheses. To better understand the role of hypotheses in debugging, we conducted two studies. In a preliminary study, we found that, even with the benefit of modern internet resources, incorrect hypotheses can cause developers to investigate irrelevant information and block progress. We then conducted a controlled experiment where 20 developers debugged and recorded their hypotheses. We found that developers have few hypotheses, two per defect. Having a correct hypothesis early strongly predicted later success. We also studied the impact of two debugging aids: fault locations and potential hypotheses. Offering fault locations did not help developers formulate more correct hypotheses or debug more successfully. In contrast, offering potential hypotheses made developers six times more likely to succeed. These results demonstrate the potential of future debugging tools that enable finding and sharing relevant hypotheses.


title: Unveiling the molecular mechanism of SARS-CoV-2 main protease inhibition from 92 crystal structures

abstract: Currently, there is no effective antiviral drugs nor vaccine for coronavirus disease 2019 (COVID-19) caused by acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Due to its high conservativeness and low similarity with human genes, SARS-CoV-2 main protease (M$^{\text{pro}}$) is one of the most favorable drug targets. However, the current understanding of the molecular mechanism of M$^{\text{pro}}$ inhibition is limited by the lack of reliable binding affinity ranking and prediction of existing structures of M$^{\text{pro}}$-inhibitor complexes. This work integrates mathematics and deep learning (MathDL) to provide a reliable ranking of the binding affinities of 92 SARS-CoV-2 M$^{\text{pro}}$ inhibitor structures. We reveal that Gly143 residue in M$^{\text{pro}}$ is the most attractive site to form hydrogen bonds, followed by Cys145, Glu166, and His163. We also identify 45 targeted covalent bonding inhibitors. Validation on the PDBbind v2016 core set benchmark shows the MathDL has achieved the top performance with Pearson's correlation coefficient ($R_p$) being 0.858. Most importantly, MathDL is validated on a carefully curated SARS-CoV-2 inhibitor dataset with the averaged $R_p$ as high as 0.751, which endows the reliability of the present binding affinity prediction. The present binding affinity ranking, interaction analysis, and fragment decomposition offer a foundation for future drug discovery efforts.


title: Local Algebraic Effect Theories

abstract: Algebraic effects are computational effects that can be described with a set of basic operations and equations between them. As many interesting effect handlers do not respect these equations, most approaches assume a trivial theory, sacrificing both reasoning power and safety. We present an alternative approach where the type system tracks equations that are observed in subparts of the program, yielding a sound and flexible logic, and paving a way for practical optimizations and reasoning tools.


title: BeCAPTCHA: Behavioral Bot Detection using Touchscreen and Mobile Sensors benchmarked on HuMIdb

abstract: In this paper we study the suitability of a new generation of CAPTCHA methods based on smartphone interactions. The heterogeneous flow of data generated during the interaction with the smartphones can be used to model human behavior when interacting with the technology and improve bot detection algorithms. For this, we propose BeCAPTCHA, a CAPTCHA method based on the analysis of the touchscreen information obtained during a single drag and drop task in combination with the accelerometer data. The goal of BeCAPTCHA is to determine whether the drag and drop task was realized by a human or a bot. We evaluate the method by generating fake samples synthesized with Generative Adversarial Neural Networks and handcrafted methods. Our results suggest the potential of mobile sensors to characterize the human behavior and develop a new generation of CAPTCHAs. The experiments are evaluated with HuMIdb (Human Mobile Interaction database), a novel multimodal mobile database that comprises 14 mobile sensors acquired from 600 users. HuMIdb is freely available to the research community.


title: Modeling cell crawling strategies with a bistable model: From amoeboid to fan-shaped cell motion

abstract: Eukaryotic cell motility involves a complex network of interactions between biochemical components and mechanical processes. The cell employs this network to polarize and induce shape changes that give rise to membrane protrusions and retractions, ultimately leading to locomotion of the entire cell body. The combination of a nonlinear reaction-diffusion model of cell polarization, noisy bistable kinetics, and a dynamic phase field for the cell shape permits us to capture the key features of this complex system to investigate several motility scenarios, including amoeboid and fan-shaped forms as well as intermediate states with distinct displacement mechanisms. We compare the numerical simulations of our model to live cell imaging experiments of motile {\it Dictyostelium discoideum} cells under different developmental conditions. The dominant parameters of the mathematical model that determine the different motility regimes are identified and discussed.


title: Gelfand type problems involving the 1-Laplacian operator

abstract: In this paper, the theory of Gelfand problems is adapted to the 1--Laplacian setting. Concretely, we deal with the following problem \begin{equation*} \left\{\begin{array}{cc} -\Delta_1u=\lambda f(u) &\hbox{in }\Omega\,;\\[2mm] u=0 &\hbox{on }\partial\Omega\,; \end{array} \right. \end{equation*} where $\Omega\subset\mathbb{R}^N$ ($N\ge1$) is a domain, $\lambda \geq 0$ and $f\>:\>[0,+\infty[\to]0,+\infty[$ is any continuous increasing and unbounded function with $f(0)>0$. It is proved the existence of a threshold $\lambda^*=\frac{h(\Omega)}{f(0)}$ (being $h(\Omega)$ the Cheeger constant of $\Omega$) such that there exists no solution when $\lambda>\lambda^*$ and the trivial function is always a solution when $\lambda\le\lambda^*$. The radial case is analyzed in more detail showing the existence of multiple solutions (even singular) as well as the behaviour of solutions to problems involving the $p$--Laplacian as $p$ tends to 1, which allows us to identify proper solutions through an extra condition.


title: Orthogonal Projections on Hyperplanes Intertwined With Unitaries

abstract: Fix a point in a finite-dimensional complex vector space and consider the sequence of iterates of this point under the composition of a unitary map with the orthogonal projection on the hyperplane orthogonal to the starting point. We prove that, generically, the series of the squared norms of these iterates sums to the dimension of the underlying space. This leads us to construct a (device-dependent) dimension witness for quantum systems which involves the probabilities of obtaining certain strings of outcomes in a sequential yes-no measurement. The exact formula for this series in non-generic cases is provided as well as its analogue in the real case.


title: Rydberg spectrum of a single trapped Ca$^+$ ion: A Floquet analysis

abstract: We compute the Rydberg spectrum of a single Ca$^+$ ion in a Paul trap by incorporating various internal and external coupling terms of the ion to the trap in the Hamiltonian. The coupling terms include spin-orbit coupling in Ca$^+$, charge (electron and ionic core) coupling to the radio frequency and static fields, ion-electron coupling in the Paul trap, and ion center-of-mass coupling. The electronic Rydberg states are precisely described by a one-electron model potential for e$^-$+Ca$^{2+}$, and accurate eigenenergies, quantum defect parameters, and static and tensor polarizabilities for a number of excited Rydberg states are obtained. The time-periodic rf Hamiltonian is expanded in the Floquet basis, and the trapping-field-broadened Rydberg lines are compared with recent observations of Ca$^+(23P)$ and Ca$^+(52F)$ Rydberg lines.


title: The Maximum Surface Area Polyhedron with Five Vertices Inscribed in the Sphere $\mathbb{S}^2$

abstract: This article focuses on the problem of analytically determining the optimal placement of five points on the unit sphere $\mathbb{S}^2$ so that the surface area of the convex hull of the points is maximized. It is shown that the optimal polyhedron has a trigonal bipyramidal structure with two vertices placed at the north and south poles and the other three vertices forming an equilateral triangle inscribed in the equator. This result confirms a conjecture of Akkiraju, who conducted a numerical search for the maximizer. As an application to crystallography, the surface area discrepancy is considered as a measure of distortion between an observed coordination polyhedron and an ideal one. The main result yields a formula for the surface area discrepancy of any coordination polyhedron with five vertices.


title: Microrheology of semiflexible filament solutions based on relaxation simulations

abstract: We present an efficient computational methodology to obtain the viscoelastic response of dilute solutions of semiflexible filaments. By considering an approach based on the fluctuation-dissipation theorem, we were able to evaluate the dynamical properties of probe particles immersed in solutions of semiflexible filaments from relaxation simulations with a relatively low computational cost and higher precision in comparison to those based on stochastic dynamics. We used a microrheological approach to obtain the complex shear modulus and the complex viscosity of the solution through its compliance which was obtained directly from the dynamical properties of a probe particle attached to an effective medium described by a mesoscopic model, i.e., an effective filament model (EFM). The relaxation simulations were applied to assess the effects of the bending energy on the viscoelasticity of semiflexible filament solutions and our methodology was validated by comparing the numerical results to experimental data on DNA and collagen solutions.

categories: cond-mat.soft cond-mat.mtrl-sci physics.bio-ph physics.flu-dyn


title: Thermodynamically consistent gradient elasticity with an internal variable

abstract: The role of thermodynamics in continuum mechanics and the derivation of the proper constitutive relations is a discussed subject of Rational Mechanics. The classical literature did not use the accumulated knowledge of thermostatics and was very critical with the heuristic methods of irreversible thermodynamics. In this paper, a small strain gradient elasticity theory is constructed with memory effects and dissipation. The method is nonequilibrium thermodynamics with internal variables; therefore, the constitutive relations are compatible with thermodynamics by construction. Thermostatic Gibbs relation is introduced for elastic bodies with a single tensorial internal variable. The thermodynamic potentials are first-order weakly nonlocal, and the entropy production is calculated. Then the constitutive functions and the evolution equation of the internal variable is constructed. The second law analysis has shown a contribution of gradient terms to the stress, also without dissipation.


title: A Comparative Study of Long and Short GRBs. II. A Multi-wavelength Method to distinguish Type II (massive star) and Type I (compact star) GRBs

abstract: Gamma Ray Burst (GRBs) are empirically classified as long-duration GRBs (LGRBs, $>$ 2s) and short-duration GRBs (SGRBs, $<$ 2s). Physically they can be grouped into two distinct progenitor categories: those originating from collapse of massive stars (also known as Type II) and those related to mergers of compact stars (also known as Type I). Even though most LGRBs are Type II and most SGRBs are Type I, the duration criterion is not always reliable to determine the physical category of a certain GRB. Based on our previous comprehensive study of the multi-wavelength properties of long and short GRBs, here we utilize the Naive Bayes method to physically classify GRBs as Type I and Type II GRBs based on multi-wavelength criteria. It results in 0.5\% training error rate and 1\% test error rate. Moreover, there is a gap [-1.2, -0.16] in the distribution of the posterior Odds, $\log O({\rm II:I})$, the Type II to Type I probability ratio. Therefore, we propose to use ${\cal O} = \log O({\rm II:I})+0.7$ as the parameter to classify GRBs into Type I ($<0$) or Type II ($>0$). The only confirmed Type I GRB, GRB 170817A, has log $O({\rm II:I})=-10$. According to this criterion, the supernova-less long GRBs 060614 and 060505 belong to Type I, and two controversial short GRBs 090426 and 060121 belong to Type II.


title: Generalized Grad-Shafranov equation for non-axisymmetric MHD equilibria

abstract: The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad-Shafranov equation, which is a single elliptic PDE in two space dimensions. By revealing a hidden symmetry, we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad-Shafranov equation. In contrast to solutions of the classical Grad-Shafranov equation, solutions of he generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad-Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.


title: Deep Learning for Portfolio Optimization

abstract: We adopt deep learning models to directly optimise the portfolio Sharpe ratio. The framework we present circumvents the requirements for forecasting expected returns and allows us to directly optimise portfolio weights by updating model parameters. Instead of selecting individual assets, we trade Exchange-Traded Funds (ETFs) of market indices to form a portfolio. Indices of different asset classes show robust correlations and trading them substantially reduces the spectrum of available assets to choose from. We compare our method with a wide range of algorithms with results showing that our model obtains the best performance over the testing period, from 2011 to the end of April 2020, including the financial instabilities of the first quarter of 2020. A sensitivity analysis is included to understand the relevance of input features and we further study the performance of our approach under different cost rates and different risk levels via volatility scaling.


title: Deletion and contraction in configuration spaces of graphs

abstract: The aim of this article is to provide space level maps between configuration spaces of graphs that are predicted by algebraic manipulations of cellular chains. More explicitly, we consider edge contraction and half-edge deletion, and identify the homotopy cofibers in terms of configuration spaces of simpler graphs. The construction's main benefit lies in making the operations functorial - in particular, graph minors give rise to compatible maps at the level of fundamental groups as well as generalized (co)homology theories. As applications we provide a long exact sequence for half-edge deletion in any generalized cohomology theory, compatible with cohomology operations such as the Steenrod and Adams operations, allowing for inductive calculations in this general context. We also show that the generalized homology of unordered configuration spaces is finitely generated as a representation of the opposite graph minor category.


title: Nanosensitive optical coherence tomography to assess wound healing within the cornea

abstract: Optical Coherence Tomography (OCT) is a non-invasive depth resolved optical imaging modality, that enables high resolution, cross-sectional imaging in biological tissues and materials at clinically relevant depths. Though OCT offers high resolution imaging, the best ultra-high-resolution OCT systems are limited to imaging structural changes with a resolution of one micron on a single B-scan within very limited depth. Nanosensitive OCT (nsOCT) is a recently developed technique that is capable of providing enhanced sensitivity of OCT to structural changes. Improving the sensitivity of OCT to detect structural changes at the nanoscale level, to a depth typical for conventional OCT, could potentially improve the diagnostic capability of OCT in medical applications. In this paper, we demonstrate the capability of nsOCT to detect structural changes deep in the rat cornea following superficial corneal injury.


title: Self-generating lower bounds and continuation for the Boltzmann equation

abstract: For the spatially inhomogeneous, non-cutoff Boltzmann equation posed in the whole space $\mathbb R^3_x$, we establish pointwise lower bounds that appear instantaneously even if the initial data contains vacuum regions. Our lower bounds depend only on the initial data and upper bounds for the mass and energy densities of the solution. As an application, we improve the weakest known continuation criterion for large-data solutions, by removing the assumptions of mass bounded below and entropy bounded above.


title: The Manufacturing Data and Machine Learning Platform: Enabling Real-time Monitoring and Control of Scientific Experiments via IoT

abstract: IoT devices and sensor networks present new opportunities for measuring, monitoring, and guiding scientific experiments. Sensors, cameras, and instruments can be combined to provide previously unachievable insights into the state of ongoing experiments. However, IoT devices can vary greatly in the type, volume, and velocity of data they generate, making it challenging to fully realize this potential. Indeed, synergizing diverse IoT data streams in near-real time can require the use of machine learning (ML). In addition, new tools and technologies are required to facilitate the collection, aggregation, and manipulation of sensor data in order to simplify the application of ML models and in turn, fully realize the utility of IoT devices in laboratories. Here we will demonstrate how the use of the Argonne-developed Manufacturing Data and Machine Learning (MDML) platform can analyze and use IoT devices in a manufacturing experiment. MDML is designed to standardize the research and operational environment for advanced data analytics and AI-enabled automated process optimization by providing the infrastructure to integrate AI in cyber-physical systems for in situ analysis. We will show that MDML is capable of processing diverse IoT data streams, using multiple computing resources, and integrating ML models to guide an experiment.


title: On the orbits associated with the Collatz conjecture

abstract: This article is based upon previous work by Sousa Ramos and his collaborators. They first prove that the existence of only one orbit associated with the Collatz conjecture is equivalent to the determinant of each matrix of a certain sequence of matrices to have the same value. These matrices are called Collatz matrices. The second step in their work would be to calculate this determinant for each of the Collatz matrices. Having calculated this determinant for the first few terms of the sequence of matrices, their plan was to prove the determinant of the current term equals the determinant of the previous one. Unfortunately, they could not prove it for the cases where the dimensions of the matrices are 26+54l or 44+54l, where l is a positive integer. In the current article we improve on these results.


title: Reconstruction and Optimization of Coherent Synthesis by Fourier Optics Based Genetic Algorithm

abstract: We present a numerical method for the reconstruction and optimization of complex field synthesis using coherent pulse combination systems. A genetic algorithm utilizing a Fourier optics based propagation method is developed for accurate convergence and modeling of near and far field distributions, achieving better than $\pi/10$ phase accuracy in reconstructed input parameters.


title: Measurement of the cosmic muon annual and diurnal flux variation with the COSINE-100 detector

abstract: We report measurements of annual and diurnal modulations of the cosmic-ray muon rate in the Yangyang underground laboratory (Y2L) using 952 days of COSINE-100 data acquired between September 2016 and July 2019. A correlation of the muon rate with the atmospheric temperature is observed and its amplitude on the muon rate is determined. The effective atmospheric temperature and muon rate variations are positively correlated with a measured effective temperature coefficient of $\alpha_{T}$ = 0.80 $\pm$ 0.11. This result is consistent with a model of meson production in the atmosphere. We also searched for a diurnal modulation in the underground muon rate by comparing one-hour intervals. No significant diurnal modulation of the muon rate was observed.


title: A six-octave optical frequency comb from a scalable few-cycle erbium fiber laser

abstract: A compact and robust coherent laser light source that provides spectral coverage from the ultraviolet to infrared is desirable for numerous applications, including heterodyne super resolution imaging[1], broadband infrared microscopy[2], protein structure determination[3], and standoff atmospheric trace-gas detection[4]. Addressing these demanding measurement problems, laser frequency combs[5] combine user-defined spectral resolution with sub-femtosecond timing and waveform control to enable new modalities of high-resolution, high-speed, and broadband spectroscopy[6-9]. In this Letter we introduce a scalable source of near-single-cycle, 0.56 MW pulses generated from robust and low-noise erbium fiber (Er:fiber) technology, and we use it to generate a frequency comb that spans six octaves from the ultraviolet (350 nm) to mid-infrared (22500 nm). The high peak power allows us to exploit the second-order nonlinearities in infrared-transparent, nonlinear crystals (LiNbO$_3$, GaSe, and CSP) to provide a robust source of phase-stable infrared ultra-short pulses with simultaneous spectral brightness exceeding that of an infrared synchrotron[10]. Additional cascaded second-order nonlinearities in LiNbO$_3$ lead to comb generation with four octaves of simultaneous coverage (0.350 to 5.6 $\mu$m). With a comb-tooth linewidth of 10 kHz at 193 THz, we realize a notable spectral resolving power exceeding 10$^{10}$ across 0.86 PHz of bandwidth. We anticipate that this compact and accessible technology will open new opportunities for multi-band precision spectroscopy, coherent microscopy, ultra-high sensitivity nanoscopy, astronomical spectroscopy, and precision carrier-envelope phase (CEP) stable strong field phenomena.


title: HAADF-STEM block-scanning strategy for local measurement of strain at the nanoscale

abstract: Lattice strain measurement of nanoscale semiconductor devices is crucial for the semiconductor industry as strain substantially improves the electrical performance of transistors. High resolution scanning transmission electron microscopy (HR-STEM) imaging is an excellent tool that provides spatial resolution at the atomic scale and strain information by applying Geometric Phase Analysis or image fitting procedures. However, HR-STEM images regularly suffer from scanning distortions and sample drift during image acquisition. In this paper, we propose a new scanning strategy that drastically reduces artefacts due to drift and scanning distortion, along with extending the field of view. The method allows flexible tuning of the spatial resolution and decouples the choice of field of view from the need for local atomic resolution. It consists of the acquisition of a series of independent small subimages containing an atomic resolution image of the local lattice. All subimages are then analysed individually for strain by fitting a nonlinear model to the lattice images. The obtained experimental strain maps are quantitatively benchmarked against the Bessel diffraction technique. We demonstrate that the proposed scanning strategy approaches the performance of the diffraction technique while having the advantage that it does not require specialized diffraction cameras.


title: Band structure and inter-tube optical transitions in double-walled carbon nanotubes

abstract: Usually, in optical spectra of double-walled carbon nanotubes (DWCNTs) weak van der Waals coupling between the layers leads only to a small shift of transition energies with respect to their values in pristine single-walled nanotubes. However, recent results have shown that the Rayleigh spectrum of the DWCNT (12,11)@(17,16) contains additional peaks. Using the tight binding approximation, we demonstrate that in specific DWCNTs the interlayer coupling can slightly modify the band structure of pristine nanotubes in such a way that the unconventional inter-tube electronic transitions become possible and additional peaks in the DWCNT optical spectrum appear. Using the known experimental data on 118 optical transitions in DWCNTs, in addition to the recently published case we reveal six more DWCNTs with inter-tube transitions and obtain geometrical selection rules permitting them. In few dozens of DWCNTs our approach yields the energies of electronic transitions close to the experimentally observed ones and may be useful for structural identification of this-type nanotubes.


title: Feynman-Kac formula for iterated derivatives of the parabolic Anderson model

abstract: The purpose of this paper is to establish a Feynman-Kac formula for the moments of the iterated Malliavin derivatives of the solution to the parabolic Anderson model in terms of pinned Brownian motions. As an application, we obtain estimates for the moments of the iterated derivatives of the solution.


title: Extraction of many-body Chern number from a single wave function

abstract: The quantized Hall conductivity of integer and fractional quantum Hall (IQH and FQH) states is directly related to a topological invariant, the many-body Chern number. The conventional calculation of this invariant in interacting systems requires a family of many-body wave functions parameterized by twist angles in order to calculate the Berry curvature. In this paper, we demonstrate how to extract the Chern number given a single many-body wave function, without knowledge of the Hamiltonian. For FQH states, our method requires one additional integer invariant as input: the number of $2\pi$ flux quanta, $s$, that must be inserted to obtain a topologically trivial excitation. As we discuss, $s$ can be obtained in principle from the degenerate set of ground state wave functions on the torus, without knowledge of the Hamiltonian. We perform extensive numerical simulations involving IQH and FQH states to validate these methods.

categories: cond-mat.str-el hep-th math-ph math.MP quant-ph


title: Mixing and Combustion in a Laminar Shear Layer with Imposed Counterflow

abstract: Three-dimensional laminar flow structures with mixing, chemical reaction, normal strain, and shear strain qualitatively representative of turbulent combustion at the small scales are analyzed. A mixing layer is subjected to counterflow in the transverse y- and z-directions. Both non-reactive and reactive flows are examined. Reduction of the three-dimensional boundary-layer equations to a one-dimensional similar form is obtained allowing for heat and mass diffusion with variations in density and properties. In steady configurations, a set of ODEs governs the three velocity components as well as the scalar-field variables. A flamelet model for individual diffusion flames with combined shear and normal strain is developed. Another model with solution in similar form is obtained for a configuration with a dominant diffusion flame and a weaker fuel-rich premixed flame. Results for the velocity and scalar fields are found for ranges of Damkohler number Da, normal strain rate due to the counterflow, streamwise-velocity ratio across the mixing layer, Prandtl number, and Mach number. For the flamelet model, a conserved scalar is cast as the independent variable to give an alternative description of the results. The imposed normal strain decreases mixing-layer thickness and increases scalar gradients and transport rates. There is indication of diffusion control for partially premixed flames in the multi-branched flame situation. The enhancement of the mixing and combustion rates by imposed normal strain on a shear layer can be very substantial. Also, the imposition of shear strain and thereby vorticity on the counterflow can be substantial indicating the need for flamelet models with both shear strain and normal strain.


title: Hashing algorithms, optimized mappings and massive parallelization of multiconfigurational methods for bosons

abstract: Numerical routines for Fock states indexing and to handle creation and annihilation operators in the spanned multiconfigurational space are developed. From the combinatorial problem of fitting particles in a truncated basis of individual particle states, which defines the spanned multiconfigurational space, a hashing function is provided based on a metric to sort all possible configurations, which refers to sets of occupation numbers required in the definition of Fock states. Despite the hashing function unambiguously relates the configuration to the coefficient index of the many-particle state expansion in the Fock basis, averages of creation and annihilation operators can be a highly demanding computation, especially when they are embedded in a time-dependent problem. Therefore, improvements in the conversion between configurations after the action of creation and annihilation operators are thoroughly inspected, highlighting the advantages and additional memory consumption. We also exploit massive parallel processors from graphics processor units with CUDA to improve a routine to act with the many-body Hamiltonian matrix on the spanned multiconfigurational space, which demonstrated quantitatively the scalability of the problem. The improvements shown here seem promising especially for calculations involving a large number of particles, in which case, the optimized CUDA code provided a drastic performance gain of roughly fifty times faster than a single core processor. The codes were consistently tested with an application to the Lieb-Liniger gas, evaluating the ground state and comparing with the analytical solution.


title: Performance analysis and optimization of power systems with spatially correlated noise

abstract: Based on stochastic differential equations (SDEs), we analyse the overall performance of heterogeneous power systems network, subject to spatially distributed and correlated noise with random initial conditions. We determine bounds on the H_2 norm of the heterogeneous system based on a closed-form of the norm of the homogeneous power system. Then, we formulate possible scenarios for performance optimization and link these to applications for network design and control problems in power systems. Our results are corroborated by numerical simulations from Kundur's four-machine two-area network after adaption to our setup.


title: Phone Features Improve Speech Translation

abstract: End-to-end models for speech translation (ST) more tightly couple speech recognition (ASR) and machine translation (MT) than a traditional cascade of separate ASR and MT models, with simpler model architectures and the potential for reduced error propagation. Their performance is often assumed to be superior, though in many conditions this is not yet the case. We compare cascaded and end-to-end models across high, medium, and low-resource conditions, and show that cascades remain stronger baselines. Further, we introduce two methods to incorporate phone features into ST models. We show that these features improve both architectures, closing the gap between end-to-end models and cascades, and outperforming previous academic work -- by up to 9 BLEU on our low-resource setting.


title: Dark Matter Halo Sparsity of Modified Gravity Scenarios

abstract: Modified Gravity (MG) scenarios have been advocated to account for the dark energy phenomenon in the universe. These models predict departures from General Relativity on large cosmic scales that can be tested through a variety of probes such as observations of galaxy clusters among others. Here, we investigate the imprint of MG models on the internal mass distribution of cluster-like halos as probed by the dark matter halo sparsity. To this purpose we perform a comparative analysis of the properties of the halo sparsity using N-body simulation halo catalogs of a standard flat $\Lambda$CDM model and MG scenarios from the DUSTGRAIN-pathfinder simulation suite. We find that the onset of the screening mechanism leaves a distinct signature in the redshift evolution of the ensemble average halos sparsity. Measurements of the sparsity of galaxy clusters from currently available mass estimates are unable to test MG models due to the large uncertainties on the cluster masses. We show that this should be possible in the future provided large cluster samples with cluster masses determined to better than $30\%$ accuracy level.


title: Microscopic models for Kitaev's sixteenfold way of anyon theories

abstract: In two dimensions, the topological order described by $\mathbb{Z}_2$ gauge theory coupled to free or weakly interacting fermions with a nonzero spectral Chern number $\nu$ is classified by $\nu \; \mathrm{mod}\; 16$ as predicted by Kitaev [Ann. Phys. 321, 2 (2006)]. Here we provide a systematic and complete construction of microscopic models realizing this so-called sixteenfold way of anyon theories. These models are defined by $\Gamma$ matrices satisfying the Clifford algebra, enjoy a global $\mathrm{SO}(\nu)$ symmetry, and live on either square or honeycomb lattices depending on the parity of $\nu$. We show that all these models are exactly solvable by using a Majorana representation and characterize the topological order by calculating the topological spin of an anyonic quasiparticle and the ground-state degeneracy. The possible relevance of the $\nu=2$ and $\nu=3$ models to materials with Kugel-Khomskii-type spin-orbital interactions is discussed.


title: Detection of the diffuse HI emission in the Circumgalactic Medium of NGC 891 and NGC 4565

abstract: We present detections of 21-cm emission from neutral hydrogen (HI) in the circumgalactic medium (CGM) of the local edge-on galaxies NGC 891 and NGC 4565 using the Robert C. Byrd Green Bank Telescope (GBT). With our 5$\sigma$ sensitivity of $8.2 \times 10^{16}$ cm$^{-2}$ calculated over a 20 km s$^{-1}$ channel, we achieve $>5\sigma$ detections out to $90-120$ kpc along the minor axes. The velocity width of the CGM emission is as large as that of the disk $\approx 500$ km s$^{-1}$, indicating the existence of a diffuse component permeating the halo. We compare our GBT measurements with interferometric data from the Westerbork Synthesis Radio Telescope (WSRT). The WSRT maps the HI emission from the disk at high S/N but has limited surface brightness sensitivity at the angular scales probed with the GBT. After convolving the WSRT data to the spatial resolution of the GBT (FWHM = 9.1$'$), we find that the emission detected by the WSRT accounts for $48^{+15}_{-25}$% ($58^{+4}_{-18}$%) of the total flux recovered by the GBT from the CGM of NGC 891(NGC 4565). The existence of significant GBT-only flux suggests the presence of a large amount of diffuse, low column density HI emission in the CGM. For reasonable assumptions, the extended diffuse HI could account for $5.2\pm0.9$% and $2.0\pm0.8$% of the total HI emission of NGC 891 and NGC 4565.


title: ProTuner: Tuning Programs with Monte Carlo Tree Search

abstract: We explore applying the Monte Carlo Tree Search (MCTS) algorithm in a notoriously difficult task: tuning programs for high-performance deep learning and image processing. We build our framework on top of Halide and show that MCTS can outperform the state-of-the-art beam-search algorithm. Unlike beam search, which is guided by greedy intermediate performance comparisons between partial and less meaningful schedules, MCTS compares complete schedules and looks ahead before making any intermediate scheduling decision. We further explore modifications to the standard MCTS algorithm as well as combining real execution time measurements with the cost model. Our results show that MCTS can outperform beam search on a suite of 16 real benchmarks.


title: Sorting probability of Catalan posets

abstract: We show that the sorting probability of the Catalan poset} $P_n$ satisfies $\delta(P_n)= O \bigl(n^{-5/4}\bigr)$.


title: Lepton dark matter portal in the inert Zee model

abstract: The inert Zee model is an extension of the Zee model for neutrino masses. This new model explains the dark matter relic abundance, generates a one-loop neutrino masses and forbids tree-level Higgs-mediated flavor changing neutral currents. Although the dark matter phenomenology of the model is similar to that of the inert doublet model, the presence of new vector-like fermions opens the lepton portal as a new dark matter annihilation channel. We study the impact of such a new portal in the low mass regime and show the parameter space allowed by direct and indirect searches of dark matter. Remarkably, we show that the region for $m_{H^0} \lesssim$ 70 GeV is recovered for $\lambda_L \lesssim 10^{-3}$. We also show that future experiments like LZ and DARWIN could test a large region of the parameter space of the model.


title: Great SCO2T! Rapid tool for carbon sequestration science, engineering, and economics

abstract: CO2 capture and storage (CCS) technology is likely to be widely deployed in coming decades in response to major climate and economics drivers: CCS is part of every clean energy pathway that limits global warming to 2C or less and receives significant CO2 tax credits in the United States. These drivers are likely to stimulate capture, transport, and storage of hundreds of millions or billions of tonnes of CO2 annually. A key part of the CCS puzzle will be identifying and characterizing suitable storage sites for vast amounts of CO2. We introduce a new software tool called SCO2T (Sequestration of CO2 Tool, pronounced "Scott") to rapidly characterizing saline storage reservoirs. The tool is designed to rapidly screen hundreds of thousands of reservoirs, perform sensitivity and uncertainty analyses, and link sequestration engineering (injection rates, reservoir capacities, plume dimensions) to sequestration economics (costs constructed from around 70 separate economic inputs). We describe the novel science developments supporting SCO2T including a new approach to estimating CO2 injection rates and CO2 plume dimensions as well as key advances linking sequestration engineering with economics. Next, we perform a sensitivity and uncertainty analysis of geology combinations (including formation depth, thickness, permeability, porosity, and temperature) to understand the impact on carbon sequestration. Through the sensitivity analysis we show that increasing depth and permeability both can lead to increased CO2 injection rates, increased storage potential, and reduced costs, while increasing porosity reduces costs without impacting the injection rate (CO2 is injected at a constant pressure in all cases) by increasing the reservoir capacity.


title: Stochasticity and heterogeneity in the transmission dynamics of SARS-CoV-2

abstract: SARS-CoV-2 causing COVID-19 disease has moved rapidly around the globe, infecting millions and killing hundreds of thousands. The basic reproduction number, which has been widely used and misused to characterize the transmissibility of the virus, hides the fact that transmission is stochastic, is dominated by a small number of individuals, and is driven by super-spreading events (SSEs). The distinct transmission features, such as high stochasticity under low prevalence, and the central role played by SSEs on transmission dynamics, should not be overlooked. Many explosive SSEs have occurred in indoor settings stoking the pandemic and shaping its spread, such as long-term care facilities, prisons, meat-packing plants, fish factories, cruise ships, family gatherings, parties and night clubs. These SSEs demonstrate the urgent need to understand routes of transmission, while posing an opportunity that outbreak can be effectively contained with targeted interventions to eliminate SSEs. Here, we describe the potential types of SSEs, how they influence transmission, and give recommendations for control of SARS-CoV-2.


title: Multiple resolution residual network for automatic thoracic organs-at-risk segmentation from CT

abstract: We implemented and evaluated a multiple resolution residual network (MRRN) for multiple normal organs-at-risk (OAR) segmentation from computed tomography (CT) images for thoracic radiotherapy treatment (RT) planning. Our approach simultaneously combines feature streams computed at multiple image resolutions and feature levels through residual connections. The feature streams at each level are updated as the images are passed through various feature levels. We trained our approach using 206 thoracic CT scans of lung cancer patients with 35 scans held out for validation to segment the left and right lungs, heart, esophagus, and spinal cord. This approach was tested on 60 CT scans from the open-source AAPM Thoracic Auto-Segmentation Challenge dataset. Performance was measured using the Dice Similarity Coefficient (DSC). Our approach outperformed the best-performing method in the grand challenge for hard-to-segment structures like the esophagus and achieved comparable results for all other structures. Median DSC using our method was 0.97 (interquartile range [IQR]: 0.97-0.98) for the left and right lungs, 0.93 (IQR: 0.93-0.95) for the heart, 0.78 (IQR: 0.76-0.80) for the esophagus, and 0.88 (IQR: 0.86-0.89) for the spinal cord.


title: Challenges in Combating COVID-19 Infodemic -- Data, Tools, and Ethics

abstract: While the COVID-19 pandemic continues its global devastation, numerous accompanying challenges emerge. One important challenge we face is to efficiently and effectively use recently gathered data and find computational tools to combat the COVID-19 infodemic, a typical information overloading problem. Novel coronavirus presents many questions without ready answers; its uncertainty and our eagerness in search of solutions offer a fertile environment for infodemic. It is thus necessary to combat the infodemic and make a concerted effort to confront COVID-19 and mitigate its negative impact in all walks of life when saving lives and maintaining normal orders during trying times. In this position paper of combating the COVID-19 infodemic, we illustrate its need by providing real-world examples of rampant conspiracy theories, misinformation, and various types of scams that take advantage of human kindness, fear, and ignorance. We present three key challenges in this fight against the COVID-19 infodemic where researchers and practitioners instinctively want to contribute and help. We demonstrate that these three challenges can and will be effectively addressed by collective wisdom, crowdsourcing, and collaborative research.


title: First-Passage Time Distributions in Two-State Protein Folding Kinetics: Exploring the Native-Like States vs Overcoming the Free Energy Barrier

abstract: Using a beta-hairpin protein as a representative example of two-state folders, we studied how the exploration of native-like states affects the folding kinetics. It has been found that the first-passage time (FPT) distributions are essentially single-exponential not only for the times to overcome the free energy barrier that separates unfolded and native-like states but also for the times to find the native state among the native-like ones. If the protein explores native-like states for a time much longer than the time to overcome the free energy barrier, which was found to be characteristic of high temperatures, the resulting FPT distribution to reach the native state remains close to exponential but the mean FPT (MFPT) is determined not by the height of the free energy barrier but by the time to explore native-like states. The mean time to overcome the free energy barrier is found to be in reasonable agreement with the Kramers rate formula and generally far shorter than the MFPT to reach the native state. The time to find the native state among native-like ones increases with temperature, which explains the known U-shape dependence of the MFPTs on temperature.


title: Third order optical nonlinearity of three dimensional massless Dirac fermions

abstract: We present analytic expressions for the electronic contributions to the linear conductivity $\sigma^{(1)}_{3d}(\omega)$ and the third order optical conductivity $\sigma^{(3)}_{3d}(\omega_1,\omega_2,\omega_3)$ of three dimensional massless Dirac fermions, the quasi-particles relevant for the low energy excitation of topological Dirac semimetals and Weyl semimetals. Although there is no gap for massless Dirac fermions, a finite chemical potential $\mu$ can lead to an effective gap parameter, which plays an important role in the qualitative features of interband optical transitions. For gapless linear dispersion in three dimension, the imaginary part of the linear conductivity diverges as a logarithmic function of the cutoff energy, while the real part is linear with photon frequency $\omega$ as $\hbar\omega>2|\mu|$. The third order conductivity exhibits features very similar to those of two dimensional Dirac fermions, i.e., graphene, but with the amplitude for a single Dirac cone generally two orders of magnitude smaller in three dimension than in two dimension. There are many resonances associated with the chemical potential induced gap parameters, and divergences associated with the intraband transitions. The details of the third order conductivity are discussed for third harmonic generation, the Kerr effect and two-photon carrier injection, parametric frequency conversion, and two-color coherent current injection. Although the expressions we derive are limited to the clean limit at zero temperature, the generalization to include phenomenological relaxation processes at finite temperature is straightforward and is presented.


title: Learning Secured Modulation With Deep Adversarial Neural Networks

abstract: Growing interest in utilizing the wireless spectrum by heterogeneous devices compels us to rethink the physical layer security to protect the transmitted waveform from an eavesdropper. We propose an end-to-end symmetric key neural encryption and decryption algorithm with a modulation technique, which remains undeciphered by an eavesdropper, equipped with the same neural network and trained on the same dataset as the intended users. We solve encryption and modulation as a joint problem for which we map the bits to complex analog signals, without adhering to any particular encryption algorithm or modulation technique. We train to cooperatively learn encryption and decryption algorithms between our trusted pair of neural networks, while eavesdropper's model is trained adversarially on the same data to minimize the error. We introduce a discrete activation layer with a defined gradient to combat noise in a lossy channel. Our results show that a trusted pair of users can exchange data bits in both clean and noisy channels, where a trained adversary cannot decipher the data.


title: An ENAS Based Approach for Constructing Deep Learning Models for Breast Cancer Recognition from Ultrasound Images

abstract: Deep Convolutional Neural Networks (CNN) provides an "end-to-end" solution for image pattern recognition with impressive performance in many areas of application including medical imaging. Most CNN models of high performance use hand-crafted network architectures that require expertise in CNNs to utilise their potentials. In this paper, we applied the Efficient Neural Architecture Search (ENAS) method to find optimal CNN architectures for classifying breast lesions from ultrasound (US) images. Our empirical study with a dataset of 524 US images shows that the optimal models generated by using ENAS achieve an average accuracy of 89.3%, surpassing other hand-crafted alternatives. Furthermore, the models are simpler in complexity and more efficient. Our study demonstrates that the ENAS approach to CNN model design is a promising direction for classifying ultrasound images of breast lesions.


title: Instabilities in neutron-star postmerger remnants

abstract: Using nonlinear, fully relativistic, simulations we investigate the dynamics and gravitational-wave signature associated with instabilities in neutron star postmerger remnants. For simplified models of the remnant we establish the presence of an instability in stars with moderate $T/|W|$, the ratio between the kinetic and the gravitational potential energies. Detailed analysis of the density oscillation pattern reveals a \emph{local} instability in the inner region of the more realistic differential rotation profile. We apply Rayleigh's inflection theorem and Fj{\o}rtoft's theorem to analyze the stability criteria concluding that this \emph{inner local} instability originates from a shear instability close to the peak of the angular velocity profile, and that it later evolves into a fast-rotating m=2 oscillation pattern. We discuss the importance of the presence of a corotation point in the fluid, its connection with the shear instability, and comparisons to the Rossby wave and Papaloizou Pringle instabilities considered in the wider literature.


title: Measuring the mass of the supermassive black hole of the lenticular galaxy NGC 4546

abstract: Most galaxies with a well-structured bulge host a supermassive black hole (SMBH) in their centre. Stellar kinematics models applied to adaptive optics (AO) assisted integral field unit observations are well-suited to measure the SMBH mass ($M_{BH}$) and also the total mass-to-light ratio [$(M/L)_{TOT}$] and possible anisotropies in the stellar velocity distribution in the central region of galaxies. In this work, we used new AO assisted Near-Infrared Integral Field Spectrometer (NIFS) observations and also photometric data from the Hubble Space Telescope Legacy Archive of the galaxy NGC 4546 in order to determine its SMBH mass. To do this, we applied the Jeans Anisotropic Modelling (JAM) method to fit the average second velocity moment in the line of sight $(\overline{v^2_{los}})$ of the stellar structure. In addition, we also obtained $(M/L)_{TOT}$ and the classical anisotropy parameter $\beta_z$=1--($\sigma_z$/$\sigma_R$)$^2$ for this object within a field of view of 200$\times$200 pc$^2$. Maps of the stellar radial velocity and of the velocity dispersion were built for this galaxy using the penalized pixel fitting ({\sc ppxf}) technique. We applied the Multi Gaussian Expansion procedure to fit the stellar brightness distribution. Using JAM, the best-fitting model for $\overline{v^2_{los}}$ of the stellar structure was obtained with $(M/L)_{TOT}$ = 4.34$\pm$0.07 (Johnson's R band), $M_{BH}$ = (2.56$\pm$0.16)$\times$10$^8$ M$_\odot$ and $\beta_z$ = --0.015$\pm$0.03 (3$\sigma$ confidence level). With these results, we found that NGC 4546 follows the $M_{BH}$ $\times$ $\sigma$ relation. We also measured the central velocity dispersion within a radius of 1 arcsec of this object as $\sigma_c$ = 241$\pm$2 km s$^{-1}$.


title: A quantitative framework for exploring exit strategies from the COVID-19 lockdown

abstract: Following the highly restrictive measures adopted by many countries for combating the current pandemic, the number of individuals infected by SARS-CoV-2 and the associated number of deaths is steadily decreasing. This fact, together with the impossibility of maintaining the lockdown indefinitely, raises the crucial question of whether it is possible to design an exit strategy based on quantitative analysis. Guided by rigorous mathematical results, we show that this is indeed possible: we present a robust numerical algorithm which can compute the cumulative number of deaths that will occur as a result of increasing the number of contacts by a given multiple, using as input only the most reliable of all data available during the lockdown, namely the cumulative number of deaths.


title: HAWC J2227+610 and its association with G106.3+2.7, a new potential Galactic PeVatron

abstract: We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time.


title: The dust-to-gas ratio, size distribution, and dust fall-back fraction of comet 67P/Churyumov-Gerasimenko: Inferences from linking the optical and dynamical properties of the inner comae

abstract: In this work, we present results that simultaneously constrain the dust size distribution, dust-to-gas ratio, fraction of dust re-deposition, and total mass production rates for comet 67P/Churyumov-Gerasimenko. We use a 3D Direct Simulation Monte Carlo (DSMC) gas dynamics code to simulate the inner gas coma of the comet for the duration of the Rosetta mission. The gas model is constrained by ROSINA/COPS data. Further, we simulate for different epochs the inner dust coma using a 3D dust dynamics code including gas drag and the nucleus' gravity. Using advanced dust scattering properties these results are used to produce synthetic images that can be compared to the OSIRIS data set. These simulations allow us to constrain the properties of the dust coma and the total gas and dust production rates. We determined a total volatile mass loss of $(6.1 \pm 1.5) \cdot 10^9$~kg during the 2015 apparition. Further, we found that power-laws with $q=3.7^{+0.57}_{-0.078}$ are consistent with the data. This results in a total of $5.1^{+6.0}_{-4.9}\cdot10^9$~kg of dust being ejected from the nucleus surface, of which $4.4^{+4.9}_{-4.2}\cdot10^9$~kg escape to space and $6.8^{+11}_{-6.8}\cdot10^8$~kg (or an equivalent of $14^{+22}_{-14}$~cm over the smooth regions) is re-deposited on the surface. This leads to a dust-to-gas ratio of $0.73^{+1.3}_{-0.70}$ for the escaping material and $0.84^{+1.6}_{-0.81}$ for the ejected material. We have further found that the smallest dust size must be strictly smaller than $\sim30\mu$m and nominally even smaller than $\sim12\mu$m.


title: Modular Politics: Toward a Governance Layer for Online Communities

abstract: Governance in online communities is an increasingly high-stakes challenge, and yet many basic features of offline governance legacies--juries, political parties, term limits, and formal debates, to name a few--are not in the feature-sets of the software most community platforms use. Drawing on the paradigm of Institutional Analysis and Development, this paper proposes a strategy for addressing this lapse by specifying basic features of a generalizable paradigm for online governance called Modular Politics. Whereas classical governance typologies tend to present a choice among wholesale ideologies, such as democracy or oligarchy, Modular Politics would enable platform operators and their users to build bottom-up governance processes from computational components that are modular and composable, highly versatile in their expressiveness, portable from one context to another, and interoperable across platforms. This kind of approach could implement pre-digital governance systems as well as accelerate innovation in uniquely digital techniques. As diverse communities share and connect their components and data, governance could occur through a ubiquitous network layer. To that end, this paper proposes the development of an open standard for networked governance.


title: On the Difficulty of Membership Inference Attacks

abstract: Recent studies propose membership inference (MI) attacks on deep models, where the goal is to infer if a sample has been used in the training process. Despite their apparent success, these studies only report accuracy, precision, and recall of the positive class (member class). Hence, the performance of these attacks have not been clearly reported on negative class (non-member class). In this paper, we show that the way the MI attack performance has been reported is often misleading because they suffer from high false positive rate or false alarm rate (FAR) that has not been reported. FAR shows how often the attack model mislabel non-training samples (non-member) as training (member) ones. The high FAR makes MI attacks fundamentally impractical, which is particularly more significant for tasks such as membership inference where the majority of samples in reality belong to the negative (non-training) class. Moreover, we show that the current MI attack models can only identify the membership of misclassified samples with mediocre accuracy at best, which only constitute a very small portion of training samples. We analyze several new features that have not been comprehensively explored for membership inference before, including distance to the decision boundary and gradient norms, and conclude that deep models' responses are mostly similar among train and non-train samples. We conduct several experiments on image classification tasks, including MNIST, CIFAR-10, CIFAR-100, and ImageNet, using various model architecture, including LeNet, AlexNet, ResNet, etc. We show that the current state-of-the-art MI attacks cannot achieve high accuracy and low FAR at the same time, even when the attacker is given several advantages. The source code is available at https://github.com/shrezaei/MI-Attack.


title: Scale-free spanning trees: complexity, bounds and algorithms

abstract: We introduce and study the general problem of finding a most "scale-free-like" spanning tree of a connected graph. It is motivated by a particular problem in epidemiology, and may be useful in studies of various dynamical processes in networks. We employ two possible objective functions for this problem and introduce the corresponding algorithmic problems termed $m$-SF and $s$-SF Spanning Tree problems. We prove that those problems are APX- and NP-hard, respectively, even in the classes of cubic, bipartite and split graphs. We study the relations between scale-free spanning tree problems and the max-leaf spanning tree problem, which is the classical algorithmic problem closest to ours. For split graphs, we explicitly describe the structure of optimal spanning trees and graphs with extremal solutions. Finally, we propose two Integer Linear Programming formulations and two fast heuristics for the $s$-SF Spanning Tree problem, and experimentally assess their performance using simulated and real data.


title: Graph-based Proprioceptive Localization Using a Discrete Heading-Length Feature Sequence Matching Approach

abstract: Proprioceptive localization refers to a new class of robot egocentric localization methods that do not rely on the perception and recognition of external landmarks. These methods are naturally immune to bad weather, poor lighting conditions, or other extreme environmental conditions that may hinder exteroceptive sensors such as a camera or a laser ranger finder. These methods depend on proprioceptive sensors such as inertial measurement units (IMUs) and/or wheel encoders. Assisted by magnetoreception, the sensors can provide a rudimentary estimation of vehicle trajectory which is used to query a prior known map to obtain location. Named as graph-based proprioceptive localization (GBPL), we provide a low cost fallback solution for localization under challenging environmental conditions. As a robot/vehicle travels, we extract a sequence of heading-length values for straight segments from the trajectory and match the sequence with a pre-processed heading-length graph (HLG) abstracted from the prior known map to localize the robot under a graph-matching approach. Using the information from HLG, our location alignment and verification module compensates for trajectory drift, wheel slip, or tire inflation level. We have implemented our algorithm and tested it in both simulated and physical experiments. The algorithm runs successfully in finding robot location continuously and achieves localization accurate at the level that the prior map allows (less than 10m).


title: Detecting Scatteredly-Distributed, Small, andCritically Important Objects in 3D OncologyImaging via Decision Stratification

abstract: Finding and identifying scatteredly-distributed, small, and critically important objects in 3D oncology images is very challenging. We focus on the detection and segmentation of oncology-significant (or suspicious cancer metastasized) lymph nodes (OSLNs), which has not been studied before as a computational task. Determining and delineating the spread of OSLNs is essential in defining the corresponding resection/irradiating regions for the downstream workflows of surgical resection and radiotherapy of various cancers. For patients who are treated with radiotherapy, this task is performed by experienced radiation oncologists that involves high-level reasoning on whether LNs are metastasized, which is subject to high inter-observer variations. In this work, we propose a divide-and-conquer decision stratification approach that divides OSLNs into tumor-proximal and tumor-distal categories. This is motivated by the observation that each category has its own different underlying distributions in appearance, size and other characteristics. Two separate detection-by-segmentation networks are trained per category and fused. To further reduce false positives (FP), we present a novel global-local network (GLNet) that combines high-level lesion characteristics with features learned from localized 3D image patches. Our method is evaluated on a dataset of 141 esophageal cancer patients with PET and CT modalities (the largest to-date). Our results significantly improve the recall from $45\%$ to $67\%$ at $3$ FPs per patient as compared to previous state-of-the-art methods. The highest achieved OSLN recall of $0.828$ is clinically relevant and valuable.