title: Second virial coefficients for helium-4 and helium-3 from accurate relativistic interaction potential

abstract: The second virial coefficient and the second acoustic virial coefficient for helium-4 and helium-3 are computed for a wide range of temperatures (0.5 - 1000K) using a highly accurate nonrelativistic interaction potential [M. Przybytek et al., Phys. Rev. Lett. 119, 123401 (2017)] and recalculated relativistic and quantum-electrodynamic components. The effects of the long-range retardation and of the nonadiabatic coupling of the nuclear and electronic motion are also taken into account. The results of our calculations represent at least fivefold improvement in accuracy compared to the previous ab initio work. The computed virial coefficients agree well with the most accurate recent measurements but have significantly smaller uncertainty.


title: FPT Algorithms for Finding Near-Cliques in $c$-Closed Graphs

abstract: Finding large cliques or cliques missing a few edges is a fundamental algorithmic task in the study of real-world graphs, with applications in community detection, pattern recognition, and clustering. A number of effective backtracking-based heuristics for these problems have emerged from recent empirical work in social network analysis. Given the NP-hardness of variants of clique counting, these results raise a challenge for beyond worst-case analysis of these problems. Inspired by the triadic closure of real-world graphs, Fox et al. (SICOMP 2020) introduced the notion of $c$-closed graphs and proved that maximal clique enumeration is fixed-parameter tractable with respect to $c$. In practice, due to noise in data, one wishes to actually discover "near-cliques", which can be characterized as cliques with a sparse subgraph removed. In this work, we prove that many different kinds of maximal near-cliques can be enumerated in polynomial time (and FPT in $c$) for $c$-closed graphs. We study various established notions of such substructures, including $k$-plexes, complements of bounded-degeneracy and bounded-treewidth graphs. Interestingly, our algorithms follow relatively simple backtracking procedures, analogous to what is done in practice. Our results underscore the significance of the $c$-closed graph class for theoretical understanding of social network analysis.


title: Improving the sensitivity of planar Fabry-P\'erot cavities via adaptive optics and mode filtering

abstract: Fabry-P\'erot (FP) cavities are fundamental and ubiquitous optical elements frequently used in various sensing applications. Here, we introduce a general theoretical framework to study arbitrary light-cavity mode interactions for planar FPs and show how optical aberrations, intrinsic to the interrogating beam or due to imperfect cavities, reduce optical sensitivity by exciting higher-order spatial modes in the cavity. We find that particular Zernike aberrations play a dominant role in sensitivity degradation, and that the general loss of sensitivity can be significantly recovered by appropriate wavefront correction or mode filtering. We then demonstrate our theoretical findings also experimentally and show that in practice the sensitivity of realistic planar FP sensors can be improved up to three-fold by a synergistic combination of adaptive optics and passive mode filtering.


title: Multi-stage Power Scheduling Framework for Data Center with Chilled Water Storage in Energy and Regulation Markets

abstract: Leveraging electrochemical and thermal energy storage systems has been proposed as a strategy to reduce peak power in data centers. Thermal energy storage systems, such as chilled water tanks, have gained increasing attention in data centers for load shifting due to their relatively small capital and operational costs compared to electrochemical energy storage. However, there are few studies investigating the possibility of utilizing thermal energy storage system with resources to provide ancillary services (e.g., frequency regulation) to the grid. This paper proposes a synergistic control strategy for the data center with a chilled water storage providing frequency regulation service by adjusting the chiller capacity, storage charging rate, and IT server CPU frequency. Then, a three-stage multi-market scheduling framework based on a model predictive control scheme is developed to minimize operational costs of data centers participating in both energy and regulation markets. The framework solves a power baseline scheduling problem, a regulation reserve problem, and a real-time power signal tracking problem sequentially. Simulation results show that utilizing the thermal energy storage can increase the regulation capacity bid, reduce energy costs and demand charges, and also harvest frequency regulation revenues. The proposed multi-market scheduling framework in a span of two days can reduce the operational costs up to 8.8% ($1,606.4) compared to the baseline with 0.2% (\$38.7) energy cost reduction, 6.5% (\$1,179.4) from demand reduction, and 2.1% (\$338.3) from regulation revenues.


title: Risk-Averse Equilibrium for Autonomous Vehicles in Stochastic Congestion Games

abstract: The fast-growing market of autonomous vehicles, unmanned aerial vehicles, and fleets in general necessitates the design of smart and automatic navigation systems considering the stochastic latency along different paths in the traffic network. The longstanding shortest path problem in a deterministic network, whose counterpart in a congestion game setting is Wardrop equilibrium, has been studied extensively, but it is well known that finding the notion of an optimal path is challenging in a traffic network with stochastic arc delays. In this work, we propose three classes of risk-averse equilibria for an atomic stochastic congestion game in its general form where the arc delay distributions are load dependent and not necessarily independent of each other. The three classes are risk-averse equilibrium (RAE), mean-variance equilibrium (MVE), and conditional value at risk level $\alpha$ equilibrium (CVaR$_\alpha$E) whose notions of risk-averse best responses are based on maximizing the probability of taking the shortest path, minimizing a linear combination of mean and variance of path delay, and minimizing the expected delay at a specified risky quantile of the delay distributions, respectively. We prove that for any finite stochastic atomic congestion game, the risk-averse, mean-variance, and CVaR$_\alpha$ equilibria exist. We show that for risk-averse travelers, the Braess paradox may not occur to the extent presented originally since players do not necessarily travel along the shortest path in expectation, but they take the uncertainty of travel time into consideration as well. We show through some examples that the price of anarchy can be improved when players are risk-averse and travel according to one of the three classes of risk-averse equilibria rather than the Wardrop equilibrium.


title: Updated analysis of jet quenching at RHIC and LHC within the light cone path integral approach

abstract: We present results of a detailed analysis of experimental data on the nuclear modification factor $R_{AA}$ and the flow coefficient $v_2$ for light hadrons from RHIC for $0.2$ TeV Au+Au collisions and from LHC for $2.76$ and $5.02$ TeV Pb+Pb, and $5.44$ TeV Xe+Xe collisions. We perform calculations within the light-cone path integral approach to induced gluon emission. We use running $\alpha_s$ which is frozen at low momenta at some value $\alpha_{s}^{fr}$. We find that the RHIC data support somewhat larger value of $\alpha_{s}^{fr}$. For the $\chi^2$ optimized values of $\alpha_{s}^{fr}$, the theoretical predictions are in reasonable agreement with data on $R_{AA}$ and $v_2$. Calculations made for different formation times and life-times of the QGP show that jet quenching at the RHIC and LHC energies is only weakly sensitive to the initial and final stages of the QCD matter evolution.


title: Shortest Secure Path in a Voronoi Diagram

abstract: We investigate the problem of computing the shortest secure path in a Voronoi diagram. Here, a path is secure if it is a sequence of touching Voronoi cells, where each Voronoi cell in the path has a uniform cost of being secured. Importantly, we allow inserting new sites, which in some cases leads to significantly shorter paths. We present an $O(n \log n)$ time algorithm for solving this problem in the plane, which uses a dynamic additive weighted Voronoi diagram to compute this path. The algorithm is an interesting combination of the continuous and discrete Dijkstra algorithms. We also implemented the algorithm using CGAL.


title: An Overview of Natural Language State Representation for Reinforcement Learning

abstract: A suitable state representation is a fundamental part of the learning process in Reinforcement Learning. In various tasks, the state can either be described by natural language or be natural language itself. This survey outlines the strategies used in the literature to build natural language state representations. We appeal for more linguistically interpretable and grounded representations, careful justification of design decisions and evaluation of the effectiveness of different approaches.


title: A Procedure for Developing Uncertainty-Consistent Vs Profiles from Inversion of Surface Wave Dispersion Data

abstract: Non-invasive surface wave methods have become a popular alternative to traditional invasive forms of site-characterization for inferring a site's subsurface shear wave velocity (Vs) structure. The advantage of surface wave methods over traditional forms of site characterization is that measurements made solely at the ground surface can be used routinely and economically to infer the subsurface structure of a site to depths of engineering interest (20-50 m), and much greater depths (>1 km) in some special cases. However, the quantification and propagation of uncertainties from surface wave measurements into the Vs profiles used in subsequent engineering analyses remains challenging. While this has been the focus of much work in recent years, and while considerable progress has been made, no approach for doing so has been widely accepted, leading analysts to either address the propagation of uncertainties in their own specialized manner or, worse, to ignore these uncertainties entirely. In response, this paper presents an easy-to-implement, effective, and verifiable method for developing uncertainty-consistent Vs profiles from inversion of surface wave dispersion data. We begin by examining four approaches presented in the literature for developing suites of Vs profiles meant to account for uncertainty present in the measured dispersion data. These methods are shown to be deficient in three specific ways. First, all approaches are shown to be highly sensitive to their many user-defined inversion input parameters, making it difficult/impossible for them to be performed repeatedly by different analysts. Second, the suites of inverted Vs profiles, when viewed in terms of their implied theoretical dispersion data, are shown to significantly underestimate the uncertainty present in the experimental dispersion data, though some may appear satisfactory when viewed purely qualitatively...


title: Implicit-explicit multirate infinitesimal GARK methods

abstract: This work focuses on the development of a new class of high-order accurate methods for multirate time integration of systems of ordinary differential equations. Unlike other recent work in this area, the proposed methods support mixed implicit-explicit (IMEX) treatment of the slow time scale. In addition to allowing this slow time scale flexibility, the proposed methods utilize a so-called `infinitesimal' formulation for the fast time scale through definition of a sequence of modified `fast' initial-value problems, that may be solved using any viable algorithm. We name the proposed class as implicit-explicit multirate infinitesimal generalized-structure additive Runge--Kutta (IMEX-MRI-GARK) methods. In addition to defining these methods, we prove that they may be viewed as specific instances of GARK methods, and derive a set of order conditions on the IMEX-MRI-GARK coefficients to guarantee both third and fourth order accuracy for the overall multirate method. Additionally, we provide three specific IMEX-MRI-GARK methods, two of order three and one of order four. We conclude with numerical simulations on two multirate test problems, demonstrating the methods' predicted convergence rates and comparing their efficiency against both legacy IMEX multirate schemes and recent third and fourth order implicit MRI-GARK methods.


title: Deep Representation Learning For Multimodal Brain Networks

abstract: Applying network science approaches to investigate the functions and anatomy of the human brain is prevalent in modern medical imaging analysis. Due to the complex network topology, for an individual brain, mining a discriminative network representation from the multimodal brain networks is non-trivial. The recent success of deep learning techniques on graph-structured data suggests a new way to model the non-linear cross-modality relationship. However, current deep brain network methods either ignore the intrinsic graph topology or require a network basis shared within a group. To address these challenges, we propose a novel end-to-end deep graph representation learning (Deep Multimodal Brain Networks - DMBN) to fuse multimodal brain networks. Specifically, we decipher the cross-modality relationship through a graph encoding and decoding process. The higher-order network mappings from brain structural networks to functional networks are learned in the node domain. The learned network representation is a set of node features that are informative to induce brain saliency maps in a supervised manner. We test our framework in both synthetic and real image data. The experimental results show the superiority of the proposed method over some other state-of-the-art deep brain network models.


title: Normal Reflection Subgroups of Complex Reflection Groups

abstract: We study normal reflection subgroups of complex reflection groups. Our approach leads to a refinement of a theorem of Orlik and Solomon to the effect that the generating function for fixed-space dimension over a reflection group is a product of linear factors involving generalized exponents. Our refinement gives a uniform proof and generalization of a recent theorem of the second author.


title: The magnetoelectric effect due to a semispherical capacitor surrounded by a spherical topologically insulating shell

abstract: We consider the magnetoelectric effect produced by a capacitor formed by two semispherical perfectly conducting plates subjected to a potential difference and surrounded by a spherical shell of a topologically insulating material. The modified Maxwell equations are solved in terms of coupled electric and magnetic scalar potentials using spherical coordinates and in the approximation where the effective magnetoelectric coupling is of the order of the fine structure constant. The emphasis is placed in the calculation of the magnetic field for several relevant configurations designed to enhance the possibility of measuring this field. The magnitudes we obtain fall within the sensitivities of magnetometers based upon nitrogen-vacancy centers in diamond as well as of devices using scanning SQUID magnetometry.


title: Implementing Stepped Pooled Testing for Rapid COVID-19 Detection

abstract: COVID-19, a viral respiratory pandemic, has rapidly spread throughout the globe. Large scale and rapid testing of the population is required to contain the disease, but such testing is prohibitive in terms of resources, cost and time. Recently RT-PCR based pooled testing has emerged as a promising way to boost testing efficiency. We introduce a stepped pooled testing strategy, a probability driven approach which significantly reduces the number of tests required to identify infected individuals in a large population. Our comprehensive methodology incorporates the effect of false negative and positive rates to accurately determine not only the efficiency of pooling but also it's accuracy. Under various plausible scenarios, we show that this approach significantly reduces the cost of testing and also reduces the effective false positive rate of tests when compared to a strategy of testing every individual of a population. We also outline an optimization strategy to obtain the pool size that maximizes the efficiency of pooling given the diagnostic protocol parameters and local infection conditions.


title: A Maximum Theorem for Incomplete Preferences

abstract: We extend Berge's Maximum Theorem to allow for incomplete preferences. We first provide a simple version of the Maximum Theorem for convex feasible sets and a fixed preference. Then, we show that if, in addition to the traditional continuity assumptions, a new continuity property for the domains of comparability holds, the limits of maximal elements along a sequence of decision problems are maximal elements in the limit problem. While this new continuity property for the domains of comparability is not generally necessary for optimality to be preserved by limits, we provide conditions under which it is necessary and sufficient.


title: Geometric implications of fast volume growth and capacity estimates

abstract: We obtain connectivity of annuli for a volume doubling metric measure Dirichlet space which satisfies a Poincar\'e inequality, a capacity estimate and a fast volume growth condition. This type of connectivity was introduced by Grigor'yan and Saloff-Coste in order to obtain stability results for Harnack inequalities and to study diffusions on manifolds with ends. As an application of our result, we obtain stability of the elliptic Harnack inequality under perturbations of the Dirichlet form with radial type weights.


title: The radius of comparison of the crossed product by a tracially strictly approximately inner action

abstract: Let $G$ be a finite group, let $A$ be an infinite-dimensional stably finite simple unital C*-algebra, and let $\alpha \colon G \to \operatorname{Aut} (A)$ be a tracially strictly approximately inner action of $G$ on $A$. Then the radius of comparison satisfies $\operatorname{rc} (A) \leq \operatorname{rc} \big( C^*(G, A, \alpha) \big)$ and if $C^*(G, A, \alpha)$ is simple, then $\operatorname{rc} (A) \leq \operatorname{rc} \big( C^*(G, A, \alpha) \big) \leq \operatorname{rc} (A^{\alpha})$. Further, the inclusion of $A$ in $C^*(G, A, \alpha)$ induces an isomorphism from the purely positive part of the Cuntz semigroup $\operatorname{Cu} (A)$ to its image in $\operatorname{Cu} \left(C^*(G, A, \alpha)\right)$. If $\alpha$ is strictly approximately inner, then in fact $\operatorname{Cu} (A) \to \operatorname{Cu} \left(C^*(G, A, \alpha) \right)$ is an ordered semigroup isomorphism onto its range. Also, for every finite group $G$ and for every $\eta \in \left(0, \frac{1}{\operatorname{card} (G)}\right)$, we construct a simple separable unital AH algebra $A$ with stable rank one and a strictly approximately inner action $\alpha \colon G \to \operatorname{Aut} (A)$ such that: (1) $\alpha$ is pointwise outer and doesn't have the weak tracial Rokhlin property. (2) $\operatorname{rc} (A) =\operatorname{rc} \left(C^*(G, A, \alpha)\right)= \eta$.


title: A bivariate extension of the Crouzeix-Palencia result with an application to Fr\'echet derivatives of matrix functions

abstract: A result by Crouzeix and Palencia states that the spectral norm of a matrix function $f(A)$ is bounded by $K = 1+\sqrt{2}$ times the maximum of $f$ on $W(A)$, the numerical range of $A$. The purpose of this work is to point out that this result extends to a certain notion of bivariate matrix functions; the spectral norm of $f\{A,B\}$ is bounded by $K^2$ times the maximum of $f$ on $W(A)\times W(B)$. As a special case, it follows that the spectral norm of the Fr\'echet derivative of $f(A)$ is bounded by $K^2$ times the maximum of $f^\prime$ on $W(A)$. An application to the convergence analysis of certain Krylov subspace methods and the extension to functions in more than two variables are discussed.


title: ASAP-NMS: Accelerating Non-Maximum Suppression Using Spatially Aware Priors

abstract: The widely adopted sequential variant of Non Maximum Suppression (or Greedy-NMS) is a crucial module for object-detection pipelines. Unfortunately, for the region proposal stage of two/multi-stage detectors, NMS is turning out to be a latency bottleneck due to its sequential nature. In this article, we carefully profile Greedy-NMS iterations to find that a major chunk of computation is wasted in comparing proposals that are already far-away and have a small chance of suppressing each other. We address this issue by comparing only those proposals that are generated from nearby anchors. The translation-invariant property of the anchor lattice affords generation of a lookup table, which provides an efficient access to nearby proposals, during NMS. This leads to an Accelerated NMS algorithm which leverages Spatially Aware Priors, or ASAP-NMS, and improves the latency of the NMS step from 13.6ms to 1.2 ms on a CPU without sacrificing the accuracy of a state-of-the-art two-stage detector on COCO and VOC datasets. Importantly, ASAP-NMS is agnostic to image resolution and can be used as a simple drop-in module during inference. Using ASAP-NMS at run-time only, we obtain an mAP of 44.2\%@25Hz on the COCO dataset with a V100 GPU.


title: Election Control by Manipulating Issue Significance

abstract: Integrity of elections is vital to democratic systems, but it is frequently threatened by malicious actors. The study of algorithmic complexity of the problem of manipulating election outcomes by changing its structural features is known as election control. One means of election control that has been proposed is to select a subset of issues that determine voter preferences over candidates. We study a variation of this model in which voters have judgments about relative importance of issues, and a malicious actor can manipulate these judgments. We show that computing effective manipulations in this model is NP-hard even with two candidates or binary issues. However, we demonstrate that the problem is tractable with a constant number of voters or issues. Additionally, while it remains intractable when voters can vote stochastically, we exhibit an important special case in which stochastic voting enables tractable manipulation.


title: On the existence of pairs of primitive and normal elements over finite fields

abstract: Let $\mathbb{F}_{q^n}$ be a finite field with $q^n$ elements, and let $m_1$ and $m_2$ be positive integers. Given polynomials $f_1(x), f_2(x) \in \mathbb{F}_q[x]$ with $\textrm{deg}(f_i(x)) \leq m_i$, for $i = 1, 2$, and such that the rational function $f_1(x)/f_2(x)$ belongs to a certain set which we define, we present a sufficient condition for the existence of a primitive element $\alpha \in \mathbb{F}_{q^n}$, normal over $\mathbb{F}_q$, such that $f_1(\alpha)/f_2(\alpha)$ is also primitive.


title: Coarse-grained spectral projection (CGSP): a deep learning-assisted approach to quantum unitary dynamics

abstract: We propose the coarse-grained spectral projection method (CGSP), a deep learning-assisted approach for tackling quantum unitary dynamic problems with an emphasis on quench dynamics. We show CGSP can extract spectral components of many-body quantum states systematically with sophisticated neural network quantum ansatz. CGSP exploits fully the linear unitary nature of the quantum dynamics, and is potentially superior to other quantum Monte Carlo methods for ergodic dynamics. Preliminary numerical results on 1D XXZ models with periodic boundary condition are carried out to demonstrate the practicality of CGSP.


title: Adaptive Control Plane Load Balancing in vSDN Enabled 5G Network

abstract: In this work, we have formulated a controllerhypervisor (C-H) pair deployment approach to provide various network functions (such as differentiation between control and data signals, different translation functions etc.) with low latency. The system model deals with a real network topology and four well defined network latency matrices with an MILP model has been well constructed to optimize latency objectives. Reverse Path-Flow Mechanism (RPFM) has been proposed in order to provide feasible solutions by keeping the network load and link capacity constrained under tolerance limit. We further minimize the H-plane load by distributing the network resources based on the arrival time of SERVICE IN requests from users. Through extensive experiment, the efficacy of the proposed system has been demonstrated in terms of significant latency reduction and an e-ULL (evolved-Ultra Low Latency) experience with efficient handling of all real-time user demands. The proposed approach can be used for similar localization problems like service chain mapping in 5G-NR, baseband unit deployment in 5G C-RAN, firewall deployment in distributed CPS.


title: Generative Adversarial Stacked Autoencoders for Facial Pose Normalization and Emotion Recognition

abstract: In this work, we propose a novel Generative Adversarial Stacked Autoencoder that learns to map facial expressions, with up to plus or minus 60 degrees, to an illumination invariant facial representation of 0 degrees. We accomplish this by using a novel convolutional layer that exploits both local and global spatial information, and a convolutional layer with a reduced number of parameters that exploits facial symmetry. Furthermore, we introduce a generative adversarial gradual greedy layer-wise learning algorithm designed to train Adversarial Autoencoders in an efficient and incremental manner. We demonstrate the efficiency of our method and report state-of-the-art performance on several facial emotion recognition corpora, including one collected in the wild.


title: E$^2$Net: An Edge Enhanced Network for Accurate Liver and Tumor Segmentation on CT Scans

abstract: Developing an effective liver and liver tumor segmentation model from CT scans is very important for the success of liver cancer diagnosis, surgical planning and cancer treatment. In this work, we propose a two-stage framework for 2D liver and tumor segmentation. The first stage is a coarse liver segmentation network, while the second stage is an edge enhanced network (E$^2$Net) for more accurate liver and tumor segmentation. E$^2$Net explicitly models complementary objects (liver and tumor) and their edge information within the network to preserve the organ and lesion boundaries. We introduce an edge prediction module in E$^2$Net and design an edge distance map between liver and tumor boundaries, which is used as an extra supervision signal to train the edge enhanced network. We also propose a deep cross feature fusion module to refine multi-scale features from both objects and their edges. E$^2$Net is more easily and efficiently trained with a small labeled dataset, and it can be trained/tested on the original 2D CT slices (resolve resampling error issue in 3D models). The proposed framework has shown superior performance on both liver and liver tumor segmentation compared to several state-of-the-art 2D, 3D and 2D/3D hybrid frameworks.


title: Entanglement-Coherence and Discord-Coherence analytical relations for X states

abstract: In this work we derive analytical relations between Entanglement and Coherence as well as between Discord and Coherence, for Bell-diagonal states and for X states, evolving under the action of several noise channels: Bit Flip, Phase Damping and Depolarizing. We demonstrate that for these families, Coherence is the fundamental correlation, that is: Coherence is necessary for the presence of Entanglement and Discord.


title: $2$-blocks in strongly biconnected directed graphs

abstract: A directed graph $G=(V,E)$ is called strongly biconnected if $G$ is strongly connected and the underlying graph of $G$ is biconnected. A strongly biconnected component of a strongly connected graph $G=(V,E)$ is a maximal vertex subset $L\subseteq V$ such that the induced subgraph on $L$ is strongly biconnected. Let $G=(V,E)$ be a strongly biconnected directed graph. A $2$-edge-biconnected block in $G$ is a maximal vertex subset $U\subseteq V$ such that for any two distict vertices $v,w \in U$ and for each edge $b\in E$, the vertices $v,w$ are in the same strongly biconnected components of $G\setminus\left\lbrace b\right\rbrace $. A $2$-strong-biconnected block in $G$ is a maximal vertex subset $U\subseteq V$ of size at least $2$ such that for every pair of distinct vertices $v,w\in U$ and for every vertex $z\in V\setminus\left\lbrace v,w \right\rbrace $, the vertices $v$ and $w$ are in the same strongly biconnected component of $G\setminus \left\lbrace v,w \right\rbrace $. In this paper we study $2$-edge-biconnected blocks and $2$-strong biconnected blocks.


title: A Bijection Between Two Different Classes of Partitions Enumerated by $p_\nu(n)$

abstract: In this paper, we give a purely bijective proof that two different partition classes that are both combinatorial interpretations of the partition function $p_\nu(n)$, a partition function related to the third order mock theta function $\nu(q)$, are equinumerous. In doing so, we give a partial solution to a combinatorial problem proposed in a paper by Andrews.


title: Probing the current-phase relation of graphene Josephson junctions using microwave measurements

abstract: We perform extensive analysis of graphene Josephson junctions embedded in microwave circuits. By comparing a diffusive junction at 15 mK with a ballistic one at 15 mK and 1 K, we are able to reconstruct the current-phase relation.


title: Quantifying the transition from spiral waves to spiral wave chimeras in a lattice of self-sustained oscillators

abstract: The present work is devoted to the detailed quantification of the transition from spiral waves to spiral wave chimeras in a network of self-sustained oscillators with two-dimensional geometry. The basic elements of the networks are the van der Pol oscillator and the FitzHugh-Nagumo neuron. Both models are in the regime of relaxation oscillations. We analyze the regime by using the indices of local sensitivity which enables us to evaluate the sensitivity of each individual oscillator at finite time. Spi-ral waves are observed in both lattices when the interaction between elements have the local character. The dynamics of all the elements is regular. There are no high-sensitive regions. We have discovered that when the coupling becomes nonlocal, the features of the systems significantly changes. The oscillation regime of the spiral wave center element switches to chaotic one. Besides this, a region with high sensitivity occurs around this oscillator. Moreover, we show that the latter expands in space with elongation of the coupling range. As a result, an incoherence cluster of the spiral wave chimera is formed exactly within this high-sensitive area. Formation of this cluster is accompanied by the sharp increase in values of the maximal Lyapunov exponent to the positive region. Furthermore, we explore that the system can even switch to hyperchaotic regime, when several Lyapunov exponents becomes positive.


title: Anisotropic Neutron Stars Modelling: Constraints in Krori-Barua Spacetime

abstract: Dense nuclear matter is expected to be anisotropic due to effects such as solidification, superfluidity, strong magnetic fields, hyperons, pion-condesation. Therefore an anisotropic neutron star core seems more realistic than an ideally isotropic one. We model anisotropic neutron stars working in the Krori-Barua (KB) ansatz without preassuming an equation of state. We show that the physics of general KB solutions is encapsulated in the compactness. Imposing physical and stability requirements yields a maximum allowed compactness $2GM/Rc^2 < 0.71$ for a KB-spacetime. We further input observational data from numerous pulsars and calculate the boundary density. We focus especially on data from the LIGO/Virgo collaboration as well as recent independent measurements of mass and radius of miilisecond pulsars with white dwarf companions by the Neutron Star Interior Composition Explorer (NICER). For these data the KB-spacetime gives the same boundary density which surprisingly equals the nuclear saturation density within the data precision. Since this value designates the boundary of a neutron core, the KB-spacetime applies naturally to neutron stars. For this boundary condition we calculate a maximum mass of 4.1 solar masses.


title: Counterfactual Learning to Rank using Heterogeneous Treatment Effect Estimation

abstract: Learning-to-Rank (LTR) models trained from implicit feedback (e.g. clicks) suffer from inherent biases. A well-known one is the position bias -- documents in top positions are more likely to receive clicks due in part to their position advantages. To unbiasedly learn to rank, existing counterfactual frameworks first estimate the propensity (probability) of missing clicks with intervention data from a small portion of search traffic, and then use inverse propensity score (IPS) to debias LTR algorithms on the whole data set. These approaches often assume the propensity only depends on the position of the document, which may cause high estimation variance in applications where the search context (e.g. query, user) varies frequently. While context-dependent propensity models reduce variance, accurate estimations may require randomization or intervention on a large amount of traffic, which may not be realistic in real-world systems, especially for long tail queries. In this work, we employ heterogeneous treatment effect estimation techniques to estimate position bias when intervention click data is limited. We then use such estimations to debias the observed click distribution and re-draw a new de-biased data set, which can be used for any LTR algorithms. We conduct simulations with varying experiment conditions and show the effectiveness of the proposed method in regimes with long tail queries and sparse clicks.


title: Kazhdan-Lusztig conjecture via Zastava spaces

abstract: We deduce the Kazhdan-Lusztig conjecture on the multiplicities of simple modules over a simple complex Lie algebra in Verma modules in category O from the equivariant geometric Satake correspondence and the analysis of torus fixed points in zastava spaces. We make similar speculations for the affine Lie algebras and W-algebras.


title: Improving the Estimation of the COVID-19 Effective Reproduction Number using Nowcasting

abstract: As the interactions between people increases, the impending menace of COVID-19 outbreaks materialize, and there is an inclination to apply lockdowns. In this context, it is essential to have easy-to-use indicators for people to use as a reference. The basic reproduction number of confirmed positives, $R_t$, fulfill such a role. This document proposes a data-driven approach to nowcast $R_t$ based on previous observations' statistical behavior. As more information arrives, the method naturally becomes more precise about the final count of confirmed positives. Our method's strength is that it is based on the self-reported onset of symptoms, in contrast to other methods that use the daily report's count to infer this quantity. We show that our approach may be the foundation for determining useful epidemy tracking indicators.


title: Explosive synchronization in temporal networks: A comparative study

abstract: We present a comparative study on Explosive Synchronization (ES) in temporal networks consisting of phase oscillators. The temporal nature of the networks is modeled with two configurations: (1) oscillators are allowed to move in a closed two dimensional box such that they couple with their neighbors, (2) oscillators are static and they randomly switch their coupling partners. Configuration (1) is further studied under two possible scenarios: in the first case oscillators couple to fixed numbers of neighbors while in other they couple to all oscillators lying in their circle of vision. Under these circumstances, we monitor the degrees of temporal networks, velocities, and radius of circle of vision of the oscillators, and the probability of forming connections in order to study and compare the critical values of the coupling required to induce ES in the population of phase oscillators.


title: Performance Evaluation of Differential Privacy Mechanisms in Blockchain based Smart Metering

abstract: The concept of differential privacy emerged as a strong notion to protect database privacy in an untrusted environment. Later on, researchers proposed several variants of differential privacy in order to preserve privacy in certain other scenarios, such as real-time cyber physical systems. Since then, differential privacy has rigorously been applied to certain other domains which has the need of privacy preservation. One such domain is decentralized blockchain based smart metering, in which smart meters acting as blockchain nodes sent their real-time data to grid utility databases for real-time reporting. This data is further used to carry out statistical tasks, such as load forecasting, demand response calculation, etc. However, in case if any intruder gets access to this data it can leak privacy of smart meter users. In this context, differential privacy can be used to protect privacy of this data. In this chapter, we carry out comparison of four variants of differential privacy (Laplace, Gaussian, Uniform, and Geometric) in blockchain based smart metering scenario. We test these variants on smart metering data and carry out their performance evaluation by varying different parameters. Experimental outcomes shows at low privacy budget ($\varepsilon$) and at low reading sensitivity value ($\delta$), these privacy preserving mechanisms provide high privacy by adding large amount of noise. However, among these four privacy preserving parameters Geometric parameters is more suitable for protecting high peak values and Laplace mechanism is more suitable for protecting low peak values at ($\varepsilon$ = 0.01).


title: On $L$-Functions of Modular Elliptic Curves and Certain $K3$ Surfaces

abstract: Inspired by Lehmer's conjecture on the nonvanishing of the Ramanujan $\tau$-function, one may ask whether an odd integer $\alpha$ can be equal to $\tau(n)$ or any coefficient of a newform $f(z)$. Balakrishnan, Craig, Ono, and Tsai used the theory of Lucas sequences and Diophantine analysis to characterize non-admissible values of newforms of even weight $k\geq 4$. We use these methods for weight $2$ and $3$ newforms and apply our results to $L$-functions of modular elliptic curves and certain $K3$ surfaces with Picard number $\ge 19$. In particular, for the complete list of weight $3$ newforms $f_\lambda(z)=\sum a_\lambda(n)q^n$ that are $\eta$-products, and for $N_\lambda$ the conductor of some elliptic curve $E_\lambda$, we show that if $|a_\lambda(n)|<100$ is odd with $n>1$ and $(n,2N_\lambda)=1$, then \begin{align*} a_\lambda(n) \in \,& \{-5,9,\pm 11,25, \pm41, \pm 43, -45,\pm47,49, \pm53,55, \pm59, \pm61, \pm 67\}\\ & \,\,\, \cup \, \{-69,\pm 71, \pm 73,75, \pm79,\pm81, \pm 83, \pm89,\pm 93 \pm 97, 99\}. \end{align*} Assuming the Generalized Riemann Hypothesis, we can rule out a few more possibilities leaving \begin{align*} a_\lambda(n) \in \{-5,9,\pm 11,25,-45,49,55,-69,75,\pm 81,\pm 93, 99\}. \end{align*}


title: Fault-tolerant measurement-free quantum error correction with multi-qubit gates

abstract: Measurement-free quantum error correction (MFQEC) offers an alternative to standard measurement-based QEC in platforms with an unconditional qubit reset gate. We revisit the question of fault tolerance (FT) for a measurement-free variant of the Steane code that leverages multi-qubit gates and redundant syndrome extraction, finding previously overlooked phase-flip errors that undermine FT. We then construct a revised MFQEC circuit that is resistant to all single-qubit errors, but which nonetheless cannot tolerate certain correlated errors. In order to investigate FT systematically, we introduce an efficient method to classically simulate MFQEC circuits with (i) Clifford gates for syndrome extraction, (ii) syndrome-controlled Pauli operations for decoding, and (iii) a Pauli noise model. We thereby find a pseudothreshold of $\sim0.7\%$ for our revised MFQEC Steane code under a restricted noise model previously considered in the literature. We then relax noise model assumptions to identify general requirements for FT with multi-qubit gates, finding that existing multi-qubit neutral atom gates are incompatible with fault-tolerant syndrome extraction in a straightforward implementation of both measurement-based and measurement-free variants of the Steane code. Decomposing multi-qubit gates into two-qubit gates similarly spoils FT. Finally, we discuss the theoretical ingredients that are necessary to recover FT for MFQEC codes, including single-shot FT and a recent proposal by Heu{\ss}en \textit{et al.}~[arXiv:2307.13296] to achieve FT by ``copying'' errors onto an ancilla register. By combining multi-qubit gates, redundant syndrome extraction, and copy-assisted FT, we construct a measurement-free and fault-tolerant variant of the Steane code with a pseudothreshold of $\sim0.1\%$.


title: Learning to Generate Customized Dynamic 3D Facial Expressions

abstract: Recent advances in deep learning have significantly pushed the state-of-the-art in photorealistic video animation given a single image. In this paper, we extrapolate those advances to the 3D domain, by studying 3D image-to-video translation with a particular focus on 4D facial expressions. Although 3D facial generative models have been widely explored during the past years, 4D animation remains relatively unexplored. To this end, in this study we employ a deep mesh encoder-decoder like architecture to synthesize realistic high resolution facial expressions by using a single neutral frame along with an expression identification. In addition, processing 3D meshes remains a non-trivial task compared to data that live on grid-like structures, such as images. Given the recent progress in mesh processing with graph convolutions, we make use of a recently introduced learnable operator which acts directly on the mesh structure by taking advantage of local vertex orderings. In order to generalize to 4D facial expressions across subjects, we trained our model using a high resolution dataset with 4D scans of six facial expressions from 180 subjects. Experimental results demonstrate that our approach preserves the subject's identity information even for unseen subjects and generates high quality expressions. To the best of our knowledge, this is the first study tackling the problem of 4D facial expression synthesis.


title: Limitations of Hartree-Fock with quantum resources

abstract: The Hartree-Fock problem provides the conceptual and mathematical underpinning of a large portion of quantum chemistry. As efforts in quantum technology aim to enhance computational chemistry algorithms, the fundamental Hartree-Fock problem is a natural target. While quantum computers and quantum simulation offer many prospects for the future of modern chemistry, the Hartree-Fock problem is not a likely candidate. We highlight this fact from a number of perspectives including computational complexity, practical examples, and the full characterization of the energy landscapes for simple systems.


title: Separable Indecomposable Continuum with Exactly One Composant

abstract: Indecomposable continua with one composant are $\textit{large}$ in the sense of being non-metrisable. We adapt the method of Smith $[18]$ to construct an example which is $\textit{small}$ in the sense of being separable.


title: Convergence and positivity of finite element methods for a haptotaxis model of tumoral invasion

abstract: In this paper, we consider a mathematical model for the invasion of host tissue by tumour cells in a $d$-dimensional bounded domain, $d\leq 3$. This model consists of a system of differential equations describing the evolution of cancer cell density, the extracellular matrix protein density and the matrix degrading enzyme concentration. We develop two fully discrete schemes for approximating the solutions based on the Finite Element (FE) method. For the first numerical scheme, we use a splitting technique to deal with the haptotaxis term, leading to introduce an equivalent system with a new variable given by the gradient of extracellular matrix. This scheme is well-posed and preserves the non-negativity of extracellular matrix and the degrading enzyme. We analyze error estimates and convergence towards regular solutions. The second numerical scheme is based on an equivalent formulation in which the cancer cell density equation is expressed in a divergence form through a suitable change of variables. This second numerical scheme preserves the non-negativity of all the discrete variables. Finally, we present some numerical simulations in agreement with the theoretical analysis.


title: Geno: A Developer Tool for Authoring Multimodal Interaction on Existing Web Applications

abstract: Supporting voice commands in applications presents significant benefits to users. However, adding such support to existing GUI-based web apps is effort-consuming with a high learning barrier, as shown in our formative study, due to the lack of unified support for creating multimodal interfaces. We present Geno---a developer tool for adding the voice input modality to existing web apps without requiring significant NLP expertise. Geno provides a high-level workflow for developers to specify functionalities to be supported by voice (intents), create language models for detecting intents and the relevant information (parameters) from user utterances, and fulfill the intents by either programmatically invoking the corresponding functions or replaying GUI actions on the web app. Geno further supports multimodal references to GUI context in voice commands (e.g. "move this [event] to next week" while pointing at an event with the cursor). In a study, developers with little NLP expertise were able to add multimodal voice command support for two existing web apps using Geno.


title: Transport in Conductors and Rectifiers: Mean-Field Redfield Equations and Non-Equilibrium Green's Functions

abstract: We derive a closed equation of motion for the one particle density matrix of a quantum system coupled to multiple baths using the Redfield master equation combined with a mean-field approximation. The steady-state solution may be found analytically with perturbation theory. Application of the method to a one-dimensional non-interacting quantum wire yields an expression for the current that reproduces the celebrated Landauer's formula. Nonlinear rectification is found for the case of a mesoscopic three-dimensional semiconductor p-n junction. The results are in good agreement with numerical simulations obtained using non-equilibrium Green's functions, supporting the validity of the Redfield equations for the description of transport.


title: Kernel Assisted Learning for Personalized Dose Finding

abstract: An individualized dose rule recommends a dose level within a continuous safe dose range based on patient level information such as physical conditions, genetic factors and medication histories. Traditionally, personalized dose finding process requires repeating clinical visits of the patient and frequent adjustments of the dosage. Thus the patient is constantly exposed to the risk of underdosing and overdosing during the process. Statistical methods for finding an optimal individualized dose rule can lower the costs and risks for patients. In this article, we propose a kernel assisted learning method for estimating the optimal individualized dose rule. The proposed methodology can also be applied to all other continuous decision-making problems. Advantages of the proposed method include robustness to model misspecification and capability of providing statistical inference for the estimated parameters. In the simulation studies, we show that this method is capable of identifying the optimal individualized dose rule and produces favorable expected outcomes in the population. Finally, we illustrate our approach using data from a warfarin dosing study for thrombosis patients.


title: Causal Effect Estimation and Optimal Dose Suggestions in Mobile Health

abstract: In this article, we propose novel structural nested models to estimate causal effects of continuous treatments based on mobile health data. To find the treatment regime that optimizes the expected short-term outcomes for patients, we define a weighted lag-K advantage as the value function. The optimal treatment regime is then defined to be the one that maximizes the value function. Our method imposes minimal assumptions on the data generating process. Statistical inference is provided for the estimated parameters. Simulation studies and an application to the Ohio type 1 diabetes dataset show that our method could provide meaningful insights for dose suggestions with mobile health data.


title: Two-loop mixed QCD-EW corrections to $gg \to Hg$

abstract: We compute the two-loop mixed QCD-Electroweak (QCD-EW) corrections to the production of a Higgs boson and a gluon in gluon fusion. The relevant four-point functions with internal massive propagators are expressed as multiple polylogarithms with algebraic arguments. We perform the calculation by integration over Feynman parameters and, independently, by the method of differential equations. We compute the two independent helicity amplitudes for the process and we find that they are both finite. Moreover, we observe a weight drop when all gluons have the same helicity. We also provide a simplified expression for the all-plus helicity amplitude, which is optimised for fast and reliable numerical evaluation in the physical region.


title: Machine Learning a Molecular Hamiltonian for Predicting Electron Dynamics

abstract: We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. We train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

categories: physics.comp-ph cs.LG physics.chem-ph stat.ML


title: Event Prediction in the Big Data Era: A Systematic Survey

abstract: Events are occurrences in specific locations, time, and semantics that nontrivially impact either our society or the nature, such as civil unrest, system failures, and epidemics. It is highly desirable to be able to anticipate the occurrence of such events in advance in order to reduce the potential social upheaval and damage caused. Event prediction, which has traditionally been prohibitively challenging, is now becoming a viable option in the big data era and is thus experiencing rapid growth. There is a large amount of existing work that focuses on addressing the challenges involved, including heterogeneous multi-faceted outputs, complex dependencies, and streaming data feeds. Most existing event prediction methods were initially designed to deal with specific application domains, though the techniques and evaluation procedures utilized are usually generalizable across different domains. However, it is imperative yet difficult to cross-reference the techniques across different domains, given the absence of a comprehensive literature survey for event prediction. This paper aims to provide a systematic and comprehensive survey of the technologies, applications, and evaluations of event prediction in the big data era. First, systematic categorization and summary of existing techniques are presented, which facilitate domain experts' searches for suitable techniques and help model developers consolidate their research at the frontiers. Then, comprehensive categorization and summary of major application domains are provided. Evaluation metrics and procedures are summarized and standardized to unify the understanding of model performance among stakeholders, model developers, and domain experts in various application domains. Finally, open problems and future directions for this promising and important domain are elucidated and discussed.


title: The Energy Complexity of BFS in Radio Networks

abstract: We consider a model of energy complexity in Radio Networks in which transmitting or listening on the channel costs one unit of energy and computation is free. This simplified model captures key aspects of battery-powered sensors: that battery life is most influenced by transceiver usage, and that at low transmission powers, the actual cost of transmitting and listening are very similar. The energy complexity of tasks in single-hop networks is well understood. Recent work of Chang et al. considered energy complexity in multi-hop networks and showed that $\mathsf{Broadcast}$ admits an energy-efficient protocol, by which we mean each of the $n$ nodes in the network spends $O(\text{polylog}(n))$ energy. This work left open the strange possibility that all natural problems in multi-hop networks might admit such an energy-efficient solution. In this paper we prove that the landscape of energy complexity is rich enough to support a multitude of problem complexities. Whereas $\mathsf{Broadcast}$ can be solved by an energy-efficient protocol, exact computation of $\mathsf{Diameter}$ cannot, requiring $\Omega(n)$ energy. Our main result is that $\mathsf{Breadth First Search}$ has sub-polynomial energy complexity at most $2^{O(\sqrt{\log n\log\log n})}=n^{o(1)}$; whether it admits an efficient $O(\text{polylog}(n))$-energy protocol is an open problem. Our main algorithm involves recursively solving a generalized BFS problem on a cluster graph introduced by Miller, Peng, and Xu. In this application, we make crucial use of a close relationship between distances in this cluster graph, and distances in the original network. This relationship is new and may be of independent interest.


title: Interfacial Dzyaloshinskii-Moriya interaction and spin-orbit torque in Au1-xPtx/Co bilayers with varying interfacial spin-orbit coupling

abstract: The quantitative roles of the interfacial spin-orbit coupling (SOC) in Dzyaloshinskii-Moriya interaction (DMI) and dampinglike spin-orbit torque ({\tau}DL) have remained unsettled after a decade of intensive study. Here, we report a conclusive experiment evidence that, because of the critical role of the interfacial orbital hybridization, the interfacial DMI is not necessarily a linear function of the interfacial SOC, e.g. at Au1-xPtx/Co interfaces where the interfacial SOC can be tuned significantly via strongly composition (x)-dependent spin-orbit proximity effect without varying the bulk SOC and the electronegativity of the Au1-xPtx layer. We also find that {\tau}DL in the Au1-xPtx/Co bilayers varies distinctly from the interfacial SOC as a function of x, indicating no important {\tau}DL contribution from the interfacial Rashba-Edelstein effect.


title: DBQ: A Differentiable Branch Quantizer for Lightweight Deep Neural Networks

abstract: Deep neural networks have achieved state-of-the art performance on various computer vision tasks. However, their deployment on resource-constrained devices has been hindered due to their high computational and storage complexity. While various complexity reduction techniques, such as lightweight network architecture design and parameter quantization, have been successful in reducing the cost of implementing these networks, these methods have often been considered orthogonal. In reality, existing quantization techniques fail to replicate their success on lightweight architectures such as MobileNet. To this end, we present a novel fully differentiable non-uniform quantizer that can be seamlessly mapped onto efficient ternary-based dot product engines. We conduct comprehensive experiments on CIFAR-10, ImageNet, and Visual Wake Words datasets. The proposed quantizer (DBQ) successfully tackles the daunting task of aggressively quantizing lightweight networks such as MobileNetV1, MobileNetV2, and ShuffleNetV2. DBQ achieves state-of-the art results with minimal training overhead and provides the best (pareto-optimal) accuracy-complexity trade-off.


title: Congruences for the coefficients of the Gordon and McIntosh mock theta function $\xi(q)$

abstract: Recently Gordon and McIntosh introduced the third order mock theta function $\xi(q)$ defined by $$ \xi(q)=1+2\sum_{n=1}^{\infty}\frac{q^{6n^2-6n+1}}{(q;q^6)_{n}(q^5;q^6)_{n}}. $$ Our goal in this paper is to study arithmetic properties of the coefficients of this function. We present a number of such properties, including several infinite families of Ramanujan--like congruences.


title: Reinforcement Communication Learning in Different Social Network Structures

abstract: Social network structure is one of the key determinants of human language evolution. Previous work has shown that the network of social interactions shapes decentralized learning in human groups, leading to the emergence of different kinds of communicative conventions. We examined the effects of social network organization on the properties of communication systems emerging in decentralized, multi-agent reinforcement learning communities. We found that the global connectivity of a social network drives the convergence of populations on shared and symmetric communication systems, preventing the agents from forming many local "dialects". Moreover, the agent's degree is inversely related to the consistency of its use of communicative conventions. These results show the importance of the basic properties of social network structure on reinforcement communication learning and suggest a new interpretation of findings on human convergence on word conventions.


title: Hankel Determinants of sequences related to Bernoulli and Euler Polynomials

abstract: We evaluate the Hankel determinants of various sequences related to Bernoulli and Euler numbers and special values of the corresponding polynomials. Some of these results arise as special cases of Hankel determinants of certain sums and differences of Bernoulli and Euler polynomials, while others are consequences of a method that uses the derivatives of Bernoulli and Euler polynomials. We also obtain Hankel determinants for sequences of sums and differences of powers and for generalized Bernoulli polynomials belonging to certain Dirichlet characters with small conductors. Finally, we collect and organize Hankel determinant identities for numerous sequences, both new and known, containing Bernoulli and Euler numbers and polynomials.


title: UVMBench: A Comprehensive Benchmark Suite for Researching Unified Virtual Memory in GPUs

abstract: The recent introduction of Unified Virtual Memory (UVM) in GPUs offers a new programming model that allows GPUs and CPUs to share the same virtual memory space, which shifts the complex memory management from programmers to GPU driver/ hardware and enables kernel execution even when memory is oversubscribed. Meanwhile, UVM may also incur considerable performance overhead due to tracking and data migration along with special handling of page faults and page table walk. As UVM is attracting significant attention from the research community to develop innovative solutions to these problems, in this paper, we propose a comprehensive UVM benchmark suite named UVMBench to facilitate future research on this important topic. The proposed UVMBench consists of 32 representative benchmarks from a wide range of application domains. The suite also features unified programming implementation and diverse memory access patterns across benchmarks, thus allowing thorough evaluation and comparison with current state-of-the-art. A set of experiments have been conducted on real GPUs to verify and analyze the benchmark suite behaviors under various scenarios.


title: Stoichiometry and defect superstructures in epitaxial FeSe films on SrTiO3

abstract: Cryogenic scanning tunneling microscopy is employed to investigate the stoichiometry and defects of epitaxial FeSe thin films on SrTiO3(001) substrates under various post-growth annealing conditions. Low-temperature annealing with an excess supply of Se leads to formation of Fe vacancies and superstructures, accompanied by a superconductivity (metal)-to-insulator transition in FeSe films. By contrast, high-temperature annealing could eliminate the Fe vacancies and superstructures, and thus recover the high-temperature superconducting phase of monolayer FeSe films. We also observe multilayer FeSe during low-temperature annealing, which is revealed to link with Fe vacancy formation and adatom migration. Our results document very special roles of film stoichiometry and help unravel several controversies in the properties of monolayer FeSe films.


title: A Gated and Bifurcated Stacked U-Net Module for Document Image Dewarping

abstract: Capturing images of documents is one of the easiest and most used methods of recording them. These images however, being captured with the help of handheld devices, often lead to undesirable distortions that are hard to remove. We propose a supervised Gated and Bifurcated Stacked U-Net module to predict a dewarping grid and create a distortion free image from the input. While the network is trained on synthetically warped document images, results are calculated on the basis of real world images. The novelty in our methods exists not only in a bifurcation of the U-Net to help eliminate the intermingling of the grid coordinates, but also in the use of a gated network which adds boundary and other minute line level details to the model. The end-to-end pipeline proposed by us achieves state-of-the-art performance on the DocUNet dataset after being trained on just 8 percent of the data used in previous methods.


title: Quantum phase interference in a fullerene-based molecular qutrit

abstract: High spin magnetic molecules are promising candidates for quantum information processing because they intrinsically have multiple sublevels for information storage and computational operations. However, due to their susceptibility to the environment and limitation from the selection rule, the arbitrary control of the quantum state of a multilevel system on a molecular and electron spin basis has not been realized. Here we exploit the photoexcited triplet of C70 as a molecular electron spin qutrit. After the system was initialized by photoexcitation, we prepared it into representative three-level superposition states characteristic of the qutrit, measured their density matrices, and showed the interference of the quantum phases in the superposition. The interference pattern is further interpreted as a map of evolution through time under different conditions.


title: Estimation for High-Dimensional Multi-Layer Generalized Linear Model -- Part I: The Exact MMSE Estimator

abstract: This two-part work considers the minimum means square error (MMSE) estimation problem for a high dimensional multi-layer generalized linear model (ML-GLM), which resembles a feed-forward fully connected deep learning network in that each of its layer mixes up the random input with a known weighting matrix and activates the results via non-linear functions, except that the activation here is stochastic and following some random distribution. Part I of the work focuses on the exact MMSE estimator, whose implementation is long known infeasible. For this exact estimator, an asymptotic analysis on the performance is carried out using a new replica method that is refined from certain aspects. A decoupling principle is then established, suggesting that, in terms of joint input-and-estimate distribution, the original estimation problem of multiple-input multiple-output is indeed identical to a simple single-input single-output one subjected to additive white Gaussian noise (AWGN) only. The variance of the AWGN is further shown to be determined by some coupled equations, whose dependency on the weighting and activation is given explicitly and analytically. Comparing to existing results, this paper is the first to offer a decoupling principle for the ML-GLM estimation problem. To further address the implementation issue of an exact solution, Part II proposes an approximate estimator, ML-GAMP, whose per-iteration complexity is as low as GAMP, while its asymptotic MSE (if converged) is as optimal as the exact MMSE estimator.


title: Estimation for High-Dimensional Multi-Layer Generalized Linear Model -- Part II: The ML-GAMP Estimator

abstract: This is Part II of a two-part work on the estimation for a multi-layer generalized linear model (ML-GLM) in large system limits. In Part I, we had analyzed the asymptotic performance of an exact MMSE estimator, and obtained a set of coupled equations that could characterize its MSE performance. To work around the implementation difficulty of the exact estimator, this paper continues to propose an approximate solution, ML-GAMP, which could be derived by blending a moment-matching projection into the Gaussian approximated loopy belief propagation. The ML-GAMP estimator is then shown to enjoy a great simplicity in its implementation, where its per-iteration complexity is as low as GAMP. Further analysis on its asymptotic performance also reveals that, in large system limits, its dynamical MSE behavior is fully characterized by a set of simple one-dimensional iterating equations, termed state evolution (SE). Interestingly, this SE of ML-GAMP share exactly the same fixed points with an exact MMSE estimator whose fixed points were obtained in Part I via a replica analysis. Given the Bayes-optimality of the exact implementation, this proposed estimator (if converged) is optimal in the MSE sense.


title: Over then Under Tangles

abstract: Over-then-Under (OU) tangles are oriented tangles whose strands travel through all of their over crossings before any under crossings. In this paper we discuss the idea of gliding: an algorithm by which any tangle diagram could be brought to OU form. Unfortunately, the algorithm is flawed. However, by analyzing cases in which it does succeed we obtain a braid classification result, which we also extend to virtual braids, and provide a Mathematica implementation. We discuss other instances of successful "gliding ideas" which appear in the literature - sometimes in disguise - such as the Drinfel'd double construction, Enriquez's work on quantization of Lie bialgebras, and Audoux and Meilhan's classification of welded homotopy links,


title: Wireless Performance Evaluation of Building Layouts: Closed-Form Computation of Figures of Merit

abstract: This paper presents a part of our ground-breaking work on evaluation of buildings in terms of wireless friendliness in the building-design stage. The main goal is to devise construction practices that provide for a good performance of wireless networks deployed in buildings. In this paper, the interference gain (IG) and power gain (PG) are defined as two figures of merit (FoM) of the wireless performance of buildings. The FoMs bridge the gap between building design and wireless communications industries. An approach to derive exact closed-form equations for these FoMs is proposed for the first time. The derived analytic expressions facilitate straightforward and more computationally efficient numerical evaluation of the proposed FoMs as compared to Monte Carlo simulations for well-known indoor propagation models. It is shown that the derived closed-form expression can be readily employed to evaluate the impact of building properties, such as the sizes and the aspect ratios (ARs) of rooms, on the wireless performance. The proposed approach sheds light to architects on evaluation and design of wireless-friendly building layouts.


title: Sizable spin-transfer torque in Bi/Ni80Fe20 bilayer film

abstract: The search for efficient spin conversion in Bi has attracted great attention in spin-orbitronics. In the present work, we employ spin-torque ferromagnetic resonance to investigate spin conversion in Bi/Ni80Fe20(Py) bilayer films with continuously varying Bi thickness. In contrast with previous studies, sizable spin-transfer torque (i.e., a sizable spin-conversion effect) is observed in Bi/Py bilayer film. Considering the absence of spin conversion in Bi/yttrium-iron-garnet bilayers and the enhancement of spin conversion in Bi-doped Cu, the present results indicate the importance of material combinations to generate substantial spin-conversion effects in Bi.


title: Analytical Derivation of Second-Order Deflection in Equatorial Plane of a Radially Moving Kerr-Newman Black Hole

abstract: In this work, we base on the second-order post-Minkowskian equations of motion, and apply an iterative technique to analytically derive the gravitational deflection of the relativistic particles in the equatorial plane of Kerr-Newman black hole with a radial (or longitudinal) and constant velocity. We find that the kinematically correctional effects on the second-order contributions to the deflection can be expressed into a very compact form, which are valid for both the massive particle and photon. Our result reduces to the previous formulation for the first-order deflection caused by a moving Schwarzschild black hole when the second-order contributions are dropped.


title: A New Method for Detecting Axion With Cylindrical Superconductor

abstract: We propose a method for searching dark matter axion in axion-photon conversion. We consider a superconductor of cylindrical shape under strong magnetic field. The dark matter axion generates oscillating electric field which induces oscillating superconducting current in the surface of the superconductor. The current gives rise to dipole radiation with the frequency $m_a/2\pi$ given by axion mass $m_a$. We show that the radiation flux generated by the current is of the order of $10^{-18}$W under magnetic field $\sim 3$T in the case of radius $\sim 1$cm and length $\sim 10$cm of the cylindrical superconductor. %The superconducting current flows in the surface to the depth $\lambda\simeq 5\times 10^{-6}$cm. The large amount of the radiation flux arises because Cooper pairs with large number density ( $\sim 10^{22}/\rm cm^{3}$ ) is present in the superconductor. With the high detection sensitivity, we can simultaneously search wide bandwidth of the radio frequency with existing radio telescope.


title: MINI-Net: Multiple Instance Ranking Network for Video Highlight Detection

abstract: We address the weakly supervised video highlight detection problem for learning to detect segments that are more attractive in training videos given their video event label but without expensive supervision of manually annotating highlight segments. While manually averting localizing highlight segments, weakly supervised modeling is challenging, as a video in our daily life could contain highlight segments with multiple event types, e.g., skiing and surfing. In this work, we propose casting weakly supervised video highlight detection modeling for a given specific event as a multiple instance ranking network (MINI-Net) learning. We consider each video as a bag of segments, and therefore, the proposed MINI-Net learns to enforce a higher highlight score for a positive bag that contains highlight segments of a specific event than those for negative bags that are irrelevant. In particular, we form a max-max ranking loss to acquire a reliable relative comparison between the most likely positive segment instance and the hardest negative segment instance. With this max-max ranking loss, our MINI-Net effectively leverages all segment information to acquire a more distinct video feature representation for localizing the highlight segments of a specific event in a video. The extensive experimental results on three challenging public benchmarks clearly validate the efficacy of our multiple instance ranking approach for solving the problem.


title: Migration and Refugee Crisis: a Critical Analysis of Online Public Perception

abstract: The migration rate and the level of resentments towards migrants are an important issue in modern civilisation. The infamous EU refugee crisis caught many countries unprepared, leading to sporadic and rudimentary containment measures that, in turn, led to significant public discourse. Decades of offline data collected via traditional survey methods have been utilised earlier to understand public opinion to foster peaceful coexistence. Capturing and understanding online public opinion via social media is crucial towards a joint strategic regulation spanning safety, rights of migrants and cordial integration for economic prosperity. We present a analysis of opinions on migrants and refugees expressed by the users of a very popular social platform, Twitter. We analyse sentiment and the associated context of expressions in a vast collection of tweets related to the EU refugee crisis. Our study reveals a marginally higher proportion of negative sentiments vis-a-vis migrants and a large proportion of the negative sentiments is more reflected among the ordinary users. Users with many followers and non-governmental organisations (NGO) tend to tweet favourably about the topic, offsetting the distribution of negative sentiment. We opine that they can be encouraged to be more proactive in neutralising negative attitudes that may arise concerning similar incidences.


title: RT3D: Achieving Real-Time Execution of 3D Convolutional Neural Networks on Mobile Devices

abstract: Mobile devices are becoming an important carrier for deep learning tasks, as they are being equipped with powerful, high-end mobile CPUs and GPUs. However, it is still a challenging task to execute 3D Convolutional Neural Networks (CNNs) targeting for real-time performance, besides high inference accuracy. The reason is more complex model structure and higher model dimensionality overwhelm the available computation/storage resources on mobile devices. A natural way may be turning to deep learning weight pruning techniques. However, the direct generalization of existing 2D CNN weight pruning methods to 3D CNNs is not ideal for fully exploiting mobile parallelism while achieving high inference accuracy. This paper proposes RT3D, a model compression and mobile acceleration framework for 3D CNNs, seamlessly integrating neural network weight pruning and compiler code generation techniques. We propose and investigate two structured sparsity schemes i.e., the vanilla structured sparsity and kernel group structured (KGS) sparsity that are mobile acceleration friendly. The vanilla sparsity removes whole kernel groups, while KGS sparsity is a more fine-grained structured sparsity that enjoys higher flexibility while exploiting full on-device parallelism. We propose a reweighted regularization pruning algorithm to achieve the proposed sparsity schemes. The inference time speedup due to sparsity is approaching the pruning rate of the whole model FLOPs (floating point operations). RT3D demonstrates up to 29.1$\times$ speedup in end-to-end inference time comparing with current mobile frameworks supporting 3D CNNs, with moderate 1%-1.5% accuracy loss. The end-to-end inference time for 16 video frames could be within 150 ms, when executing representative C3D and R(2+1)D models on a cellphone. For the first time, real-time execution of 3D CNNs is achieved on off-the-shelf mobiles.


title: Object-Aware Centroid Voting for Monocular 3D Object Detection

abstract: Monocular 3D object detection aims to detect objects in a 3D physical world from a single camera. However, recent approaches either rely on expensive LiDAR devices, or resort to dense pixel-wise depth estimation that causes prohibitive computational cost. In this paper, we propose an end-to-end trainable monocular 3D object detector without learning the dense depth. Specifically, the grid coordinates of a 2D box are first projected back to 3D space with the pinhole model as 3D centroids proposals. Then, a novel object-aware voting approach is introduced, which considers both the region-wise appearance attention and the geometric projection distribution, to vote the 3D centroid proposals for 3D object localization. With the late fusion and the predicted 3D orientation and dimension, the 3D bounding boxes of objects can be detected from a single RGB image. The method is straightforward yet significantly superior to other monocular-based methods. Extensive experimental results on the challenging KITTI benchmark validate the effectiveness of the proposed method.


title: Permutation-based tests for discontinuities in event studies

abstract: We propose using a permutation test to detect discontinuities in an underlying economic model at a known cutoff point. Relative to the existing literature, we show that this test is well suited for event studies based on time-series data. The test statistic measures the distance between the empirical distribution functions of observed data in two local subsamples on the two sides of the cutoff. Critical values are computed via a standard permutation algorithm. Under a high-level condition that the observed data can be coupled by a collection of conditionally independent variables, we establish the asymptotic validity of the permutation test, allowing the sizes of the local subsamples to be either be fixed or grow to infinity. In the latter case, we also establish that the permutation test is consistent. We demonstrate that our high-level condition can be verified in a broad range of problems in the infill asymptotic time-series setting, which justifies using the permutation test to detect jumps in economic variables such as volatility, trading activity, and liquidity. These potential applications are illustrated in an empirical case study for selected FOMC announcements during the ongoing COVID-19 pandemic.


title: Noise and decoherence induced by gravitons

abstract: We study quantum noise and decoherence induced by gravitons. We derive a Langevin equation of geodesic deviation in the presence of gravitons. The amplitude of noise correlations tells us that large squeezing is necessary to detect the noise. We also consider the decoherence of spatial superpositions of two massive particles caused by gravitons in the vacuum state and find that gravitons could give a relevant contribution to the decoherence. The decoherence induced by gravitons would offer new vistas to test quantum gravity in tabletop experiments.


title: A refinement on a theorem of Z. Janko

abstract: We say that a subgroup $H$ is isolated in a group $G$ if for each $x\in G$ we have either $x\in H$ or $\langle x\rangle \cap H={1}$. Z. Janko, in his paper [J. Algebra, 465(2016), 41--61], determined certain classes of finite nonabelian $p$-groups which possess some isolated subgroups. In this note, a theorem of his paper is refined.


title: Global well-posedness for the fractional Boussinesq-Coriolis system with stratification in a framework of Fourier-Besov type

abstract: We establish the global well-posedness of the 3D fractional Boussinesq-Coriolis system with stratification in a framework of Fourier type, namely spaces of Fourier-Besov type with underlying space being Morrey spaces (FBM-spaces, for short). Under suitable conditions and rescaled density fluctuation, the result is uniform with respect to the Coriolis and stratification parameters. We cover the critical case of the dissipation, namely half-Laplacian, in which the nonlocal dissipation has the same differential order as the nonlinearity and balances critically the scaling of the quadratic nonlinearities. As a byproduct, considering trivial initial temperature and null stratification, we also obtain well-posed results in FBM-spaces for the fractional Navier-Stokes-Coriolis system as well as for the Navier-Stokes equations with critical dissipation. Moreover, since small conditions are taken in the weak norm of FBM-spaces, we can consider some initial data with arbitrarily large $H^{s}$-norms, $s\geq0$.


title: Seeing the Un-Scene: Learning Amodal Semantic Maps for Room Navigation

abstract: We introduce a learning-based approach for room navigation using semantic maps. Our proposed architecture learns to predict top-down belief maps of regions that lie beyond the agent's field of view while modeling architectural and stylistic regularities in houses. First, we train a model to generate amodal semantic top-down maps indicating beliefs of location, size, and shape of rooms by learning the underlying architectural patterns in houses. Next, we use these maps to predict a point that lies in the target room and train a policy to navigate to the point. We empirically demonstrate that by predicting semantic maps, the model learns common correlations found in houses and generalizes to novel environments. We also demonstrate that reducing the task of room navigation to point navigation improves the performance further.


title: Emergent topology under slow non-adiabatic quantum dynamics

abstract: Characterization of equilibrium topological quantum phases by non-equilibrium quench dynamics provides a novel and efficient scheme in detecting topological invariants defined in equilibrium. Nevertheless, most of the previous studies have focused on the ideal sudden quench regime. Here we provide a generic non-adiabatic protocol of slowly quenching the system Hamiltonian, and investigate the non-adiabatic dynamical characterization scheme of topological phase. The {\it slow} quench protocol is realized by introducing a Coulomb-like Landau-Zener problem, and it can describe, in a unified way, the crossover from sudden quench regime (deep non-adiabatic limit) to adiabatic regime. By analytically obtaining the final state vector after non-adiabatic evolution, we can calculate the time-averaged spin polarization and the corresponding topological spin texture. We find that the topological invariants of the post-quench Hamiltonian are characterized directly by the values of spin texture on the band inversion surfaces. Compared to the sudden quench regime, where one has to take an additional step to calculate the {\it gradients} of spin polarization, this non-adiabatic characterization provides a {\it minimal} scheme in detecting the topological invariants. Our findings are not restricted to 1D and 2D topological phases under Coulomb-like quench protocol, but are also valid for higher dimensional system or different quench protocol.


title: pH-Sensitive Ultra-thin Oxide Liquid Metal System

abstract: A complex experiment was performed to unravel the simplest particle,Higgs Boson, and we performed a simple experiment to address a complex system. We tested the pH response of liquid metal EGaInSn in the form of a pendant drop and measured the sensitivity of 92.96 mV.pH-1 in the pH range from 4 to 10. We derived a unified Nernst equation to explain high pH sensitivity and coincidentally rediscovered Marcel Pourbaix illuminating work on the pH-Potential Diagram. The surface potential in the sensor probe is originated from a spontaneous electrochemical reaction purely driven by thermodynamics, rendering to the lowest system energy possible. Our findings have a great scientific significance, which could redefine the conventional concept of the ion sensing mechanism in a solid-state electrochemical sensor.


title: Bayesian EWMA and CUSUM Control Charts Under Different Loss Functions

abstract: The Exponentially Weighted Moving Average (EWMA) and Cumulative Sum (CUSUM) control charts have been used in profile monitoring to track drift shifts that occur in a monitored process. We construct Bayesian EWMA and Bayesian CUSUM charts informed by posterior and posterior predictive distributions using different loss functions, prior distributions, and likelihood distributions. A simulation study is performed, and the performance of the charts are evaluated via average run length (ARL), standard deviation of the run length (SDRL), average time to signal (ATS), and standard deviation of time to signal (SDTS). A sensitivity analysis is conducted using choices for the smoothing parameter, out-of-control shift size, and hyper-parameters of the distribution. Based on obtained results, we provide recommendations for use of the Bayesian EWMA and Bayesian CUSUM control charts.


title: Estimating heterogeneous effects of continuous exposures using Bayesian tree ensembles: revisiting the impact of abortion rates on crime

abstract: In estimating the causal effect of a continuous exposure or treatment, it is important to control for all confounding factors. However, most existing methods require parametric specification for how control variables influence the outcome or generalized propensity score, and inference on treatment effects is usually sensitive to this choice. Additionally, it is often the goal to estimate how the treatment effect varies across observed units. To address this gap, we propose a semiparametric model using Bayesian tree ensembles for estimating the causal effect of a continuous treatment of exposure which (i) does not require a priori parametric specification of the influence of control variables, and (ii) allows for identification of effect modification by pre-specified moderators. The main parametric assumption we make is that the effect of the exposure on the outcome is linear, with the steepness of this relationship determined by a nonparametric function of the moderators, and we provide heuristics to diagnose the validity of this assumption. We apply our methods to revisit a 2001 study of how abortion rates affect incidence of crime.


title: Pure metric geometry: introductory lectures

abstract: We discuss domestic affairs of metric spaces, keeping away from any extra structure. Topics include universal spaces, injective spaces, Hausdorff and Gromov--Hausdorff convergences, and ultralimits.


title: Quantum Immortality and Non-Classical Logic

abstract: The Everett Box is a device in which an observer and a lethal quantum apparatus are isolated from the rest of the universe. On a regular basis, successive trials occur, in each of which an automatic measurement of a quantum superposition inside the apparatus either causes instant death or does nothing to the observer. From the observer's perspective, the chances of surviving $m$ trials monotonically decreases with increasing $m$. As a result, if the observer is still alive for sufficiently large $m$ she must reject any interpretation of quantum mechanics which is not the many-worlds interpretation (MWI), since surviving $m$ trials becomes vanishingly unlikely in a single world, whereas a version of the observer will necessarily survive in the branching MWI universe. Here we ask whether this conclusion still holds if rather than a classical understanding of limits built on classical logic we instead require our physics to satisfy a computability requirement by investigating the Everett Box in a model of a computational universe running on a variety of constructive logic, Recursive Constructive Mathematics. We show that although the standard Everett argument rejecting non-MWI interpretations is no longer valid, we can show that Everett's conclusion still holds within a computable universe. Thus we argue that Everett's argument is strengthened and any counter-argument must be strengthened, since it holds not only in classical logic (with embedded notions of continuity and infinity) but also in a computable logic.


title: Adiabatic mass loss in binary stars. III. From the base of the red giant branch to the tip of asymptotic giant branch

abstract: The distinguishing feature of the evolution of close binary stars is the role played by the mass exchange between the component stars. Whether the mass transfer is dynamically stable is one of the essential questions in binary evolution. In the limit of extremely rapid mass transfer, the response of a donor star in an interacting binary becomes asymptotically one of adiabatic expansion. We use the adiabatic mass loss model to systematically survey the thresholds for dynamical timescale mass transfer over the entire span of possible donor star evolutionary states. We also simulate mass loss process with isentropic envelopes, the specific entropy of which is fixed to be that at the base of the convective envelope, to artificially mimic the effect of such mass loss in superadiabatic surface convection regions, where the adiabatic approximation fails. We illustrate the general adiabatic response of 3.2 Msun donor stars at different evolutionary stages. We extend our study to a grid of donor stars with different masses (from 0.1 to 100 Msun with Z = 0.02) and at different evolutionary stages. We proceed to present our criteria for dynamically unstable mass transfer in both tabular and graphical forms. For red giant branch and asymptotic giant branch donors in systems with such mass ratios, they may have convective envelopes deep enough to evolve into common envelopes on a thermal timescale, if the donor star overfills its outer Lagrangian radius. Our results show that the red giant branch and asymptotic giant branch stars tend to be more stable than previously believed, and this may be helpful to explain the abundance of observed post-AGB binary stars with an orbital period of around 1000 days.


title: A polynomial time 12-approximation algorithm for restricted Santa Claus problem

abstract: In this paper, we consider the restricted case of the problem and improve the current best approximation ratio by presenting a polynomial time 12-approximation algorithm using linear programming and semi-definite programming. Our algorithm starts by solving the configuration LP and uses the optimum value to get a 12-gap instance. This is then followed by the well-known clustering technique of Bansal and Sviridenko\cite{bansal}. We then apply the analysis of Asadpour \textit{et al.} \cite{AFS,AFS2} to show that the clustered instance has an integer solution which is at least $\frac{1}{6}$ times the best possible value, which was computed by solving the configuration LP. To find this solution, we formulate a problem called the Extended Assignment Problem, and formulate it as an LP. We then, show that the associated polytope is integral and gives us an fractional solution of value at least $\frac{1}{6}$ times the optimum. From this solution we find a solution to a new quadratic program that we introduce to select one machine from each cluster, and then we show that the resulting instance has an Assignment LP fractional solution of value at least $\frac{1}{6}$ times the optimum. We then use the well known rounding technique due to Bezakova and Dani \cite{bezakova} on the 12-gap instance to get our 12-approximate solution.


title: Joint Transmit Power and Placement Optimization for URLLC-enabled UAV Relay Systems

abstract: This letter considers an unmanned aerial vehicle (UAV)-enabled relay communication system for delivering latency-critical messages with ultra-high reliability, where the relay is operating under amplifier-and-forward (AF) mode. We aim to jointly optimize the UAV location and power to minimize decoding error probability while guaranteeing the latency constraints. Both the free-space channel model and three-dimensional (3-D) channel model are considered. For the first model, we propose a low-complexity iterative algorithm to solve the problem, while globally optimal solution is derived for the case when the signal-to-noise ratio (SNR) is extremely high. For the second model, we also propose a low-complexity iterative algorithm to solve the problem. Simulation results confirm the performance advantages of our proposed algorithms.


title: Testing goodness-of-fit and conditional independence with approximate co-sufficient sampling

abstract: Goodness-of-fit (GoF) testing is ubiquitous in statistics, with direct ties to model selection, confidence interval construction, conditional independence testing, and multiple testing, just to name a few applications. While testing the GoF of a simple (point) null hypothesis provides an analyst great flexibility in the choice of test statistic while still ensuring validity, most GoF tests for composite null hypotheses are far more constrained, as the test statistic must have a tractable distribution over the entire null model space. A notable exception is co-sufficient sampling (CSS): resampling the data conditional on a sufficient statistic for the null model guarantees valid GoF testing using any test statistic the analyst chooses. But CSS testing requires the null model to have a compact (in an information-theoretic sense) sufficient statistic, which only holds for a very limited class of models; even for a null model as simple as logistic regression, CSS testing is powerless. In this paper, we leverage the concept of approximate sufficiency to generalize CSS testing to essentially any parametric model with an asymptotically-efficient estimator; we call our extension "approximate CSS" (aCSS) testing. We quantify the finite-sample Type I error inflation of aCSS testing and show that it is vanishing under standard maximum likelihood asymptotics, for any choice of test statistic. We apply our proposed procedure both theoretically and in simulation to a number of models of interest to demonstrate its finite-sample Type I error and power.


title: Multi-label Contrastive Predictive Coding

abstract: Variational mutual information (MI) estimators are widely used in unsupervised representation learning methods such as contrastive predictive coding (CPC). A lower bound on MI can be obtained from a multi-class classification problem, where a critic attempts to distinguish a positive sample drawn from the underlying joint distribution from $(m-1)$ negative samples drawn from a suitable proposal distribution. Using this approach, MI estimates are bounded above by $\log m$, and could thus severely underestimate unless $m$ is very large. To overcome this limitation, we introduce a novel estimator based on a multi-label classification problem, where the critic needs to jointly identify multiple positive samples at the same time. We show that using the same amount of negative samples, multi-label CPC is able to exceed the $\log m$ bound, while still being a valid lower bound of mutual information. We demonstrate that the proposed approach is able to lead to better mutual information estimation, gain empirical improvements in unsupervised representation learning, and beat a current state-of-the-art knowledge distillation method over 10 out of 13 tasks.


title: Physical properties and evolution of (Sub-)millimeter selected galaxies in the galaxy formation simulation Shark

abstract: We thoroughly explore the properties of (sub)-millimeter (mm) selected galaxies (SMGs) in the Shark semi-analytic model of galaxy formation. Compared to observations, the predicted number counts at wavelengths (lambda) 0.6-2mm and redshift distributions at 0.1-2mm, agree well. At the bright end (>1mJy), Shark galaxies are a mix of mergers and disk instabilities. These galaxies display a stacked FUV-to-FIR spectrum that agrees well with observations. We predict that current optical/NIR surveys are deep enough to detect bright (>1mJy) lambda=0.85-2mm-selected galaxies at z<5, but too shallow to detect counterparts at higher redshift. A James Webb Space Telescope 10,000s survey should detect all counterparts for galaxies with $S_{\rm 0.85mm}>0.01$mJy. We predict SMG's disks contribute significantly (negligibly) to the rest-frame UV (IR). We investigate the 0<z<6 evolution of the intrinsic properties of >1mJy lambda=0.85-2mm-selected galaxies finding their: (i) stellar masses are $>10^{10.2}M_{\odot}$, with the 2mm ones tracing the most massive galaxies ($>10^{11}M_{\odot}$); (ii) specific star formation rates (SFR) are mildly (~3-10x) above the main sequence (MS); (iii) host halo masses are $\gtrsim 10^{12.3}M_{\odot}$, with 2mm galaxies tracing the most massive halos (proto-clusters); (iv) SMGs have lower dust masses ($\approx 10^{8}M_{\odot}$), higher dust temperatures ($\approx 40-45$K) and higher rest-frame V-band attenuation (>1.5) than MS galaxies; (v) sizes decrease with redshift, from 4kpc at z=1 to <1kpc at z=4; (vi) the Carbon Monoxide line spectra of $S_{\rm 0.85mm}>1$mJy sources peak at 4->3. Finally, we study the contribution of SMGs to the molecular gas and cosmic SFR density at 0<z<10, finding that >1mJy sources make a negligible contribution at z>3 and z>5, respectively, suggesting current observations have unveiled the majority of the star formation at 0<z<10.


title: Self-Loop Uncertainty: A Novel Pseudo-Label for Semi-Supervised Medical Image Segmentation

abstract: Witnessing the success of deep learning neural networks in natural image processing, an increasing number of studies have been proposed to develop deep-learning-based frameworks for medical image segmentation. However, since the pixel-wise annotation of medical images is laborious and expensive, the amount of annotated data is usually deficient to well-train a neural network. In this paper, we propose a semi-supervised approach to train neural networks with limited labeled data and a large quantity of unlabeled images for medical image segmentation. A novel pseudo-label (namely self-loop uncertainty), generated by recurrently optimizing the neural network with a self-supervised task, is adopted as the ground-truth for the unlabeled images to augment the training set and boost the segmentation accuracy. The proposed self-loop uncertainty can be seen as an approximation of the uncertainty estimation yielded by ensembling multiple models with a significant reduction of inference time. Experimental results on two publicly available datasets demonstrate the effectiveness of our semi-supervied approach.