title: A mathematical model of the COVID-19 pandemic dynamics with dependent variable infection rate: Application to the Republic of Korea

abstract: This work constructs, analyzes, and simulates a new compartmental SEIR-type model for the dynamics and potential control of the current COVID-19 pandemic. The novelty in this work is two-fold. First, the population is divided according to its compliance with disease control directives (lockdown, shelter-in-place, masks/face coverings, physical distancing, etc.) into those who fully comply and those who follow the directives partially, or are necessarily mobile (such as medical staff). This split, indirectly, reflects on the quality and consistency of these measures. This allows the assessment of the overall effectiveness of the control measures and the impact of their relaxing or tightening on the disease spread. Second, the adequate contact rate, which directly affects the infection rate, is one of the model unknowns, as it keeps track of the changes in the population behavior and the effectiveness of various disease treatment modalities via a differential inclusion. Existence, uniqueness and positivity results are proved using a nonstandard convex analysis-based approach. As a case study, the pandemic outbreak in the Republic of Korea (South Korea) is simulated. The model parameters were found by minimizing the deviation of the model prediction from the reported data over the first 100 days of the pandemic in South Korea.The simulations show that the model captures accurately the pandemic dynamics in the subsequent 75 days, which provides confidence in the model predictions and its future use. In particular, the model predicts that about 40% of the infections were not documented, which implies that asymptomatic infections contribute silently but substantially to the spread of the disease indicating that more widespread asymptomatic testing is necessary.


title: Voltage-controlled magnetic anisotropy in antiferromagnetic MgO-capped MnPt films

abstract: The magnetic anisotropy in MgO-capped MnPt films and its voltage control are studied using first-principles calculations. Sharp variation of the magnetic anisotropy with film thickness, especially in the Pt-terminated film, suggests that it may be widely tuned by adjusting the film thickness. In thick films the linear voltage control coefficient is as large as 1.5 and $-0.6$ pJ/Vm for Pt-terminated and Mn-terminated interfaces, respectively. The combination of a widely tunable magnetic anisotropy energy and a large voltage-control coefficient suggest that MgO-capped MnPt films can serve as a versatile platform for magnetic memory and antiferromagnonic applications.


title: Hyper Resolute Ultra thin Low Cost All-Dielectric Broadband Achromatic Metalenses

abstract: Metalenses offer the ground-breaking opportunity to realize highly performing low-weight, flat and ultrathin, optical elements which substantially reduce size and complexity of imaging systems. Today, a major challenge in metalenses design is still the realization of achromatic optical elements ideally focussing a broad wavelength spectrum at a single focal length. Here we present a fast and effective way to design and fabricate extremely thin all-dielectric metalenses, optimally solving achromaticity issues by means of machine learning codes. The enabling technology for fabrication is a recently developed hyper resolute two-photon direct laser writing lithography equipment. The fabricated metalenses, based on a completely flat and ultrathin design, show intriguing optical features. Overall, achromatic behavior, focal length of 1.14 mm, depth of focus of hundreds of microns and thickness of only few nanometers allow considering the design of novel and efficient imaging systems in a completely new perspective.


title: Spontaneous polarization and locomotion of an active particle with surface-mobile enzymes

abstract: We examine a mechanism of locomotion of active particles whose surface is uniformly coated with mobile enzymes. The enzymes catalyze a reaction that drives phoretic flows but their homogeneous distribution forbids locomotion by symmetry. We find that the ability of the enzymes to migrate over the surface combined with self-phoresis can lead to a spontaneous symmetry breaking instability whereby the homogeneous distribution of enzymes polarizes and the particle propels. The instability is driven by the advection of enzymes by the phoretic flows and occurs above a critical P\'eclet number. The transition to polarized motile states occurs via a supercritical or subcritical pitchfork bifurcations, the latter of which enables hysteresis and coexistence of uniform and polarized states.


title: A Survey on Security and Privacy Issues in Edge Computing-Assisted Internet of Things

abstract: Internet of Things (IoT) is an innovative paradigm envisioned to provide massive applications that are now part of our daily lives. Millions of smart devices are deployed within complex networks to provide vibrant functionalities including communications, monitoring, and controlling of critical infrastructures. However, this massive growth of IoT devices and the corresponding huge data traffic generated at the edge of the network created additional burdens on the state-of-the-art centralized cloud computing paradigm due to the bandwidth and resources scarcity. Hence, edge computing (EC) is emerging as an innovative strategy that brings data processing and storage near to the end users, leading to what is called EC-assisted IoT. Although this paradigm provides unique features and enhanced quality of service (QoS), it also introduces huge risks in data security and privacy aspects. This paper conducts a comprehensive survey on security and privacy issues in the context of EC-assisted IoT. In particular, we first present an overview of EC-assisted IoT including definitions, applications, architecture, advantages, and challenges. Second, we define security and privacy in the context of EC-assisted IoT. Then, we extensively discuss the major classifications of attacks in EC-assisted IoT and provide possible solutions and countermeasures along with the related research efforts. After that, we further classify some security and privacy issues as discussed in the literature based on security services and based on security objectives and functions. Finally, several open challenges and future research directions for secure EC-assisted IoT paradigm are also extensively provided.


title: A conditional proof of the ISP for quasinilpotent operators

abstract: The invariant subspace problem (ISP) is a well known unsolved problem in funtional analysis. While many partial results are known, the general case for complex, infinite dimensional separable Hilbert spaces is still open. It has been shown that the problem can be reduced to the case of operators which are norm limits of nilpotents. One of the most important subcases is the one of quasinilpotent operators, for which the problem has been extensively studied for many years. In this paper, we will introduce a new conjecture (supported by a heuristic argument), and we will prove conditionally that every quasinilpotent operator has a nontrivial invariant subspace. We will conclude by posing an open problem which would have deep implications regarding the ISP.


title: Evaluating Snowflake as an Indistinguishable Censorship Circumvention Tool

abstract: Tor is the most well-known tool for circumventing censorship. Unfortunately, Tor traffic has been shown to be detectable using deep-packet inspection. WebRTC is a popular web frame-work that enables browser-to-browser connections. Snowflake is a novel pluggable transport that leverages WebRTC to connect Tor clients to the Tor network. In theory, Snowflake was created to be indistinguishable from other WebRTC services. In this paper, we evaluate the indistinguishability of Snowflake. We collect over 6,500 DTLS handshakes from Snowflake, Facebook Messenger, Google Hangouts, and Discord WebRTC connections and show that Snowflake is identifiable among these applications with 100% accuracy. We show that several features, including the extensions offered and the number of packets in the handshake, distinguish Snowflake among these services. Finally, we suggest recommendations for improving identification resistance in Snowflake. We have made the dataset publicly available.


title: Optimal quantization for discrete distributions

abstract: In this paper, we first determine the optimal sets of $n$-means and the $n$th quantization errors for all $1\leq n\leq 6$ for two nonuniform discrete distributions with support the set $\{1, 2, 3, 4, 5, 6\}$. Then, for a probability distribution $P$ with support $\{\frac 1n : n\in \mathbb N\}$ associated with a mass function $f$, given by $f(x)=\frac 1 {2^k}$ if $x=\frac 1 k$ for $k\in \mathbb N$, and zero otherwise, we determine the optimal sets of $n$-means and the $n$th quantization errors for all positive integers up to $n=300$. Further, for a probability distribution $P$ with support the set $\mathbb N$ of natural number associated with a mass function $f$, given by $f(x)=\frac 1 {2^k}$ if $x=k$ for $k\in \mathbb N$, and zero otherwise, we determine the optimal sets of $n$-means and the $n$th quantization errors for all positive integers $n$. At last we discuss for a discrete distribution, if the optimal sets are given, how to obtain the probability distributions.


title: Opening practice: supporting Reproducibility and Critical spatial data science

abstract: This paper reflects on a number of trends towards a more open and reproducible approach to geographic and spatial data science over recent years. In particular it considers trends towards Big Data, and the impacts this is having on spatial data analysis and modelling. It identifies a turn in academia towards coding as a core analytic tool, and away from proprietary software tools offering 'black boxes' where the internal workings of the analysis are not revealed. It is argued that this closed form software is problematic, and considers a number of ways in which issues identified in spatial data analysis (such as the MAUP) could be overlooked when working with closed tools, leading to problems of interpretation and possibly inappropriate actions and policies based on these. In addition, this paper and considers the role that reproducible and open spatial science may play in such an approach, taking into account the issues raised. It highlights the dangers of failing to account for the geographical properties of data, now that all data are spatial (they are collected somewhere), the problems of a desire for n=all observations in data science and it identifies the need for a critical approach. This is one in which openness, transparency, sharing and reproducibility provide a mantra for defensible and robust spatial data science.


title: Interplanetary Medium as Fuel for Plasma Thrusters

abstract: We examine the use of the interplanetary substance as a possible fuel for spacecraft plasma thrusters for interplanetary flights. The solar radiation is considered as an energy source for the ionization and the acceleration of particles of captured space environment that is used as drop weight. The method of capturing the space environment, which is dependent on its density, the velocity of the spacecraft, and the power of solar radiation, defines the technical feasibility of this scheme. Both pulsed and continuous operations of the plasma accelerator are possible under some conditions as shown in the estimates.


title: A Particular Upper Expectation as Global Belief Model for Discrete-Time Finite-State Uncertain Processes

abstract: To model discrete-time finite-state uncertain processes, we argue for the use of a global belief model in the form of an upper expectation that is the most conservative one under a set of basic axioms. Our motivation for these axioms, which describe how local and global belief models should be related, is based on two possible interpretations for an upper expectation: a behavioural one similar to Walley's, and an interpretation in terms of upper envelopes of linear expectations. We show that the most conservative upper expectation satisfying our axioms, that is, our model of choice, coincides with a particular version of the game-theoretic upper expectation introduced by Shafer and Vovk. This has two important implications: it guarantees that there is a unique most conservative global belief model satisfying our axioms; and it shows that Shafer and Vovk's model can be given an axiomatic characterisation and thereby provides an alternative motivation for adopting this model, even outside their game-theoretic framework. Finally, we relate our model to the upper expectation resulting from a traditional measure-theoretic approach. We show that this measure-theoretic upper expectation also satisfies the proposed axioms, which implies that it is dominated by our model or, equivalently, the game-theoretic model. Moreover, if all local models are precise, all three models coincide.


title: BAO signatures in the 2-point angular correlations and the Hubble tension

abstract: An observational tension on estimates of the Hubble parameter, $H_0$, using early and late Universe information, is being of intense discussion in the literature. Additionally, it is of great importance to measure $H_0$ independently of CMB data and local distance ladder method. In this sense, we analyze 15 measurements of the transversal BAO scale, $\theta_{\rm BAO}$, obtained in a weakly model-dependent approach, in combination with other data sets obtained in a model-independent way, namely, Big Bang Nucleosynthesis (BBN) information, 6 gravitationally lensed quasars with measured time delays by the H0LiCOW team, and measures of cosmic chronometers (CC). We find $H_0 = 74.88_{-2.1}^{+1.9}$ km s${}^{-1}$ Mpc${}^{-1}$ and $H_0 = 72.06_{-1.3}^{+1.2}$ km s${}^{-1}$ Mpc${}^{-1}$ from $\theta_{BAO}$+BBN+H0LiCOW and $\theta_{BAO}$+BBN+CC, respectively, in fully accordance with local measurements. Moreover, we estimate the sound horizon at drag epoch, $r_{\rm d}$, independent of CMB data, and find $r_{\rm d}=144.1_{-5.5}^{+5.3}$ Mpc (from $\theta_{BAO}$+BBN+H0LiCOW) and $r_{\rm d} =150.4_{-3.3}^{+2.7}$ Mpc (from $\theta_{BAO}$+BBN+CC). In a second round of analysis, we test how the presence of a possible spatial curvature, $\Omega_k$, can influence the main results. We compare our constraints on $H_0$ and $r_{\rm d}$ with other reported values. Our results show that it is possible to use a robust compilation of transversal BAO data, $\theta_{BAO}$, jointly with other model-independent measurements, in such a way that the tension on the Hubble parameter can be alleviated.


title: Distributed Tera-Scale Similarity Search with MPI: Provably Efficient Similarity Search over billions without a Single Distance Computation

abstract: We present SLASH (Sketched LocAlity Sensitive Hashing), an MPI (Message Passing Interface) based distributed system for approximate similarity search over terabyte scale datasets. SLASH provides a multi-node implementation of the popular LSH (locality sensitive hashing) algorithm, which is generally implemented on a single machine. We show how we can append the LSH algorithm with heavy hitters sketches to provably solve the (high) similarity search problem without a single distance computation. Overall, we mathematically show that, under realistic data assumptions, we can identify the near-neighbor of a given query $q$ in sub-linear ($ \ll O(n)$) number of simple sketch aggregation operations only. To make such a system practical, we offer a novel design and sketching solution to reduce the inter-machine communication overheads exponentially. In a direct comparison on comparable hardware, SLASH is more than 10000x faster than the popular LSH package in PySpark. PySpark is a widely-adopted distributed implementation of the LSH algorithm for large datasets and is deployed in commercial platforms. In the end, we show how our system scale to Tera-scale Criteo dataset with more than 4 billion samples. SLASH can index this 2.3 terabyte data over 20 nodes in under an hour, with query times in a fraction of milliseconds. To the best of our knowledge, there is no open-source system that can index and perform a similarity search on Criteo with a commodity cluster.


title: Cabibbo anomaly versus electroweak precision tests: An exploration of extensions of the Standard Model

abstract: There is a newly emerging tension between determinations of $V_{us}$ from different sources (known as the Cabibbo anomaly), which is clearly demonstrated by the new $R(V_{us})$ observable. We explore this observable from the perspective of the Standard Model Effective Field Theory and show there is a discrepancy between $R(V_{us})$ and existing electroweak precision observables (EWPO) in a simple single operator dominated scenario. We explore all possible single particle extensions of the Standard Model that can generate the Cabibbo anomaly effect and show how they cannot reconcile the current data. We further examine the future of EWPO at the ILC or FCC-ee experiments and discuss the effect on the tension of a change in specific EW observables.


title: Tidal dissipation in evolving low-mass and solar-type stars with predictions for planetary orbital decay

abstract: We study tidal dissipation in stars with masses in the range $0.1-1.6 M_\odot$ throughout their evolution, including turbulent effective viscosity acting on equilibrium tides and inertial waves in convection zones, and internal gravity waves in radiation zones. We consider a range of stellar evolutionary models and incorporate the frequency-dependent effective viscosity acting on equilibrium tides based on the latest simulations. We compare the tidal flow and dissipation obtained with the conventional equilibrium tide, which is strictly invalid in convection zones, finding that the latter typically over-predicts the dissipation by a factor of 2-3. Dissipation of inertial waves is computed using a frequency-averaged formalism accounting for realistic stellar structure for the first time, and is the dominant mechanism for binary circularization and synchronization on the main sequence. Dissipation of gravity waves in the radiation zone assumes these waves to be fully damped (e.g.~by wave breaking), and is the dominant mechanism for planetary orbital decay. We calculate the critical planetary mass required for wave breaking as a function of stellar mass and age, and show that this mechanism predicts destruction of many hot Jupiters but probably not Earth-mass planets on the main sequence. We apply our results to compute tidal quality factors following stellar evolution, and tidal evolutionary timescales, for the orbital decay of hot Jupiters, and the spin synchronization and circularization of binary stars. We also provide predictions for shifts in transit arrival times due to tidally-driven orbital decay of hot Jupiters that may be detected with NGTS, TESS or PLATO.


title: COVID, BLM, and the polarization of US politicians on Twitter

abstract: We mapped the tweets of 520 US Congress members, focusing on analyzing their engagement with two broad topics: first, the COVID-19 pandemic, and second, the recent wave of anti-racist protest. We find that, in discussing COVID-19, Democrats frame the issue in terms of public health, while Republicans are more likely to focus on small businesses and the economy. When looking at the discourse around anti-Black violence, we find that Democrats are far more likely to name police brutality as a specific concern. In contrast, Republicans not only discuss the issue far less, but also keep their terms more general, as well as criticizing perceived protest violence.


title: Modeling and Simulation of the spread of coronavirus disease (COVID-19) in Lebanon

abstract: In this paper, we develop a probabilistic mathematical model for the spread of coronavirus disease (COVID-19). It takes into account the known special characteristics of this disease such as the existence of infectious undetected cases and the different social and infectiousness conditions of infected people. In particular, it considers the social structure and governmental measures in a country, the fraction of detected cases over the real total infected cases, and the influx of undetected infected people from outside the borders. Although the model is simple and allows a reasonable identification of its parameters, using the data provided by local authorities on this pandemic, it is also complex enough to capture the most important effects. We study the particular case of Lebanon and use its reported data to estimate the model parameters, which can be of interest for estimating the spread of COVID-19 in other countries. We show a good agreement between the reported data and the estimations given by our model. We also simulate several scenarios that help policy makers in deciding how to loosen different measures without risking a severe wave of COVID-19. We are also able to identify the main factors that lead to specific scenarios which helps in a better understanding of the spread of the virus.


title: Algebraic and symplectic viewpoint on compactifications of two-dimensional cluster varieties of finite type

abstract: In this article we explore compactifications of cluster varieties of finite type in complex dimension two. Cluster varieties can be viewed as the spec of a ring generated by theta functions and a compactification of such varieties can be given by a grading on that ring, which can be described by positive polytopes [17]. In the examples we exploit, the cluster variety can be interpreted as the complement of certain divisors in del Pezzo surfaces. In the symplectic viewpoint, they can be described via almost toric fibrations over $\R^2$ (after completion). Once identifying them as almost toric manifolds, one can symplectically view them inside other del Pezzo surfaces. So we can identify other symplectic compactifications of the same cluster variety, which we expect should also correspond to different algebraic compactifications. Both viewpoints are presented here and several compactifications have their corresponding polytopes compared. The finiteness of the cluster mutations are explored to provide cycles in the graph describing monotone Lagrangian tori in del Pezzo surfaces connected via almost toric mutation [34].


title: Sklyanin-like algebras for ($q$-)linear grids and ($q$-)para-Krawtchouk polynomials

abstract: S-Heun operators on linear and $q$-linear grids are introduced. These operators are special cases of Heun operators and are related to Sklyanin-like algebras. The Continuous Hahn and Big $q$-Jacobi polynomials are functions on which these S-Heun operators have natural actions. We show that the S-Heun operators encompass both the bispectral operators and Kalnins and Miller's structure operators. These four structure operators realize special limit cases of the trigonometric degeneration of the original Sklyanin algebra. Finite-dimensional representations of these algebras are obtained from a truncation condition. The corresponding representation bases are finite families of polynomials: the para-Krawtchouk and $q$-para-Krawtchouk ones. A natural algebraic interpretation of these polynomials that had been missing is thus obtained. We also recover the Heun operators attached to the corresponding bispectral problems as quadratic combinations of the S-Heun operators


title: On transitivity and connectedness of Cayley graphs of gyrogroups

abstract: In this work, we explore edge direction, transitivity, and connectedness of Cayley graphs of gyrogroups. More specifically, we find conditions for a Cayley graph of a gyrogroup to be undirected, transitive, and connected. We also show a relationship between the cosets of a certain type of subgyrogroups and the connected components of Cayley graphs. Some examples regarding these findings are provided.


title: A Decoupled Scheme to Solve the Mass and Momentum Conservation Equations of the Improved Darcy-Brinkman-Forchheimer Framework in Matrix Acidization

abstract: Matrix acidization simulation is a challenging task in the study of flows in porous media, due to the changing porosity in the procedure. The improved DBF framework is one model to do this simulation, and its numerical scheme discretises the mass and momentum conservation equations together to form a pressure-velocity linear system. However, this linear system can only be solved by direct solvers to solve for pressure and velocity simultaneously, since zeros appear in the diagonal of the coefficient matrix. Considering the large-scale attribute of matrix acidization simulation, the solving time of direct solvers is not intolerant. Thus, a decoupled scheme is proposed in this work to decouple the coupled pressure-velocity linear system into two independent linear systems: one is to solve for pressure, and the other one is to solve for velocity. Both of the new linear systems can be solved by parallel and iterative solvers, which guarantees the large-scale simulation can be finished in a reasonable time period. A numerical experiment is carried out to demonstrate the correctness of the decoupled scheme and its higher computing efficiency.


title: Coloring the normalized Laplacian for oriented hypergraphs

abstract: The independence number, coloring number and related parameters are investigated in the setting of oriented hypergraphs using the spectrum of the normalized Laplace operator. For the independence number, both an inertia--like bound and a ratio--like bound are shown. A Sandwich Theorem involving the clique number, the vector chromatic number and the coloring number is proved, as well as a lower bound for the vector chromatic number in terms of the smallest and the largest eigenvalue of the normalized Laplacian. In addition, spectral partition numbers are studied in relation to the coloring number.


title: Multi-Level Temporal Pyramid Network for Action Detection

abstract: Currently, one-stage frameworks have been widely applied for temporal action detection, but they still suffer from the challenge that the action instances span a wide range of time. The reason is that these one-stage detectors, e.g., Single Shot Multi-Box Detector (SSD), extract temporal features only applying a single-level layer for each head, which is not discriminative enough to perform classification and regression. In this paper, we propose a Multi-Level Temporal Pyramid Network (MLTPN) to improve the discrimination of the features. Specially, we first fuse the features from multiple layers with different temporal resolutions, to encode multi-layer temporal information. We then apply a multi-level feature pyramid architecture on the features to enhance their discriminative abilities. Finally, we design a simple yet effective feature fusion module to fuse the multi-level multi-scale features. By this means, the proposed MLTPN can learn rich and discriminative features for different action instances with different durations. We evaluate MLTPN on two challenging datasets: THUMOS'14 and Activitynet v1.3, and the experimental results show that MLTPN obtains competitive performance on Activitynet v1.3 and outperforms the state-of-the-art approaches on THUMOS'14 significantly.


title: Bayesian causal inference for count potential outcomes

abstract: The literature for count modeling provides useful tools to conduct causal inference when outcomes take non-negative integer values. Applied to the potential outcomes framework, we link the Bayesian causal inference literature to statistical models for count data. We discuss the general architectural considerations for constructing the predictive posterior of the missing potential outcomes. Special considerations for estimating average treatment effects are discussed, some generalizing certain relationships and some not yet encountered in the causal inference literature.


title: Dressing a new integrable boundary of the nonlinear Schr\"odinger equation

abstract: We further develop the method of dressing the boundary for the focusing nonlinear Schr\"odinger equation (NLS) on the half-line to include the new boundary condition presented by Zambon. Additionally, the foundation to compare the solutions to the ones produced by the mirror-image technique is laid by explicitly computing the change of scattering data under the Darboux transformation. In particular, the developed method is applied to insert pure soliton solutions.


title: SafePILCO: a software tool for safe and data-efficient policy synthesis

abstract: SafePILCO is a software tool for safe and data-efficient policy search with reinforcement learning. It extends the known PILCO algorithm, originally written in MATLAB, to support safe learning. We provide a Python implementation and leverage existing libraries that allow the codebase to remain short and modular, which is appropriate for wider use by the verification, reinforcement learning, and control communities.


title: SemEval-2020 Task 10: Emphasis Selection for Written Text in Visual Media

abstract: In this paper, we present the main findings and compare the results of SemEval-2020 Task 10, Emphasis Selection for Written Text in Visual Media. The goal of this shared task is to design automatic methods for emphasis selection, i.e. choosing candidates for emphasis in textual content to enable automated design assistance in authoring. The main focus is on short text instances for social media, with a variety of examples, from social media posts to inspirational quotes. Participants were asked to model emphasis using plain text with no additional context from the user or other design considerations. SemEval-2020 Emphasis Selection shared task attracted 197 participants in the early phase and a total of 31 teams made submissions to this task. The highest-ranked submission achieved 0.823 Matchm score. The analysis of systems submitted to the task indicates that BERT and RoBERTa were the most common choice of pre-trained models used, and part of speech tag (POS) was the most useful feature. Full results can be found on the task's website.


title: Harnessing Interpretable and Unsupervised Machine Learning to Address Big Data from Modern X-ray Diffraction

abstract: The information content of crystalline materials becomes astronomical when collective electronic behavior and their fluctuations are taken into account. In the past decade, improvements in source brightness and detector technology at modern x-ray facilities have allowed a dramatically increased fraction of this information to be captured. Now, the primary challenge is to understand and discover scientific principles from big data sets when a comprehensive analysis is beyond human reach. We report the development of a novel unsupervised machine learning approach, XRD Temperature Clustering (X-TEC), that can automatically extract charge density wave (CDW) order parameters and detect intra-unit cell (IUC) ordering and its fluctuations from a series of high-volume X-ray diffraction (XRD) measurements taken at multiple temperatures. We apply X-TEC to XRD data on a quasi-skutterudite family of materials, (Ca$_x$Sr$_{1-x}$)$_3$Rh$_4$Sn$_{13}$, where a quantum critical point arising from charge order is observed as a function of Ca concentration. We further apply X-TEC to XRD data on the pyrochlore metal, Cd$_2$Re$_2$O$_7$, to investigate its two much debated structural phase transitions and uncover the Goldstone mode accompanying them. We demonstrate how unprecedented atomic scale knowledge can be gained when human researchers connect the X-TEC results to physical principles. Specifically, we extract from the X-TEC-revealed selection rule that the Cd and Re displacements are approximately equal in amplitude, but out of phase. This discovery reveals a previously unknown involvement of $5d^2$ Re, supporting the idea of an electronic origin to the structural order. Our approach can radically transform XRD experiments by allowing in-operando data analysis and enabling researchers to refine experiments by discovering interesting regions of phase space on-the-fly.


title: Total Variation Diminishing (TVD) method for Elastohydrodynamic Lubrication (EHL) problem on Parallel Computers

abstract: In this article, we offer a novel numerical approach for the solution of elastohydrodynamic lubrication line and point contact problems using a class of total variation diminishing (TVD) schemes on parallel computers. A direct parallel approach is presented by introducing a novel solver named as projected alternate quadrant interlocking factorization (PAQIF) by solving discrete variational inequality. For one-dimensional EHL case, we use weighted change in Newton-Raphson approximation to compute the Jacobian matrix in the form of a banded matrix by dividing two subregions on the whole computation domain. Such subregion matrices are assembled by measuring the ratio of diffusive coefficient and discrete grid length on the domain of the interest. The banded matrix is then processed to parallel computers for solving discrete linearized complementarity system using PAQIF algorithm. The idea is easily extended in two-dimensional EHL case by taking appropriate splitting in x and y alternating directions respectively. Numerical experiments are performed and analyzed to validate the performance of computed solution on serial and parallel computers.


title: Learning a natural-language to LTL executable semantic parser for grounded robotics

abstract: Children acquire their native language with apparent ease by observing how language is used in context and attempting to use it themselves. They do so without laborious annotations, negative examples, or even direct corrections. We take a step toward robots that can do the same by training a grounded semantic parser, which discovers latent linguistic representations that can be used for the execution of natural-language commands. In particular, we focus on the difficult domain of commands with a temporal aspect, whose semantics we capture with Linear Temporal Logic, LTL. Our parser is trained with pairs of sentences and executions as well as an executor. At training time, the parser hypothesizes a meaning representation for the input as a formula in LTL. Three competing pressures allow the parser to discover meaning from language. First, any hypothesized meaning for a sentence must be permissive enough to reflect all the annotated execution trajectories. Second, the executor -- a pretrained end-to-end LTL planner -- must find that the observe trajectories are likely executions of the meaning. Finally, a generator, which reconstructs the original input, encourages the model to find representations that conserve knowledge about the command. Together these ensure that the meaning is neither too general nor too specific. Our model generalizes well, being able to parse and execute both machine-generated and human-generated commands, with near-equal accuracy, despite the fact that the human-generated sentences are much more varied and complex with an open lexicon. The approach presented here is not specific to LTL: it can be applied to any domain where sentence meanings can be hypothesized and an executor can verify these meanings, thus opening the door to many applications for robotic agents.


title: Mapping Local Interstellar Medium With Diffuse Interstellar Bands

abstract: With the use of the data from archives, we studied the correlations between the equivalent widths of four diffuse interstellar bands (4430$\r{A}$, 5780$\r{A}$, 5797$\r{A}$, 6284$\r{A}$) and properties of the target stars (colour excess values, distances and Galactic coordinates). Many different plots of the diffuse interstellar bands and their maps were produced and further analysed. There appears to be a structure in the plot of equivalent widths of 5780$\r{A}$ DIB (and 6284$\r{A}$ DIB) against the Galactic $x$-coordinate. The structure is well defined below $\sim150$ m$\r{A}$ and within $|x|<250$ pc, peaking around $x=170$ pc. We argue that the origin of this structure is not a statistical fluctuation. Splitting the data in the Galactic longitude into several subregions improves or lowers the well known linear relation between the equivalent widths and the colour excess, which was expected. However, some of the lines of sight display drastically different behaviour. The region within $150^\circ<l<200^\circ$ shows scatter in the correlation plots with the colour excess for all of the four bands with correlation coefficients $\textrm{R}<0.58$. We suspect that the variation of physical conditions in the nearby molecular clouds could be responsible. Finally, the area $250^\circ<l<300^\circ$ displays (from the statistical point of view) significantly lower values of equivalent widths than the other regions -- this tells us that there is either a significant underabundance of carriers (when compared with the other regions) or that this has to be a result of an observational bias.


title: Criteria for the less-equal-relation between partial Lov\'asz-vectors of digraphs

abstract: Finite digraphs $R$ and $S$ are studied with $\# {\cal H}(G,R) \leq \# {\cal H}(G,S)$ for every finite digraph $G \in \mathfrak{ D }'$, where ${\cal H}(G,H)$ is the set of order homomorphisms from $G$ to $H$ and $\mathfrak{ D }'$ is a class of finite digraphs. It is shown that for several classes $\mathfrak{ D }'$ of digraphs and $R \in \mathfrak{ D }'$, the relation $\# {\cal H}(G,R) \leq \# {\cal H}(G,S)$ for every $G \in \mathfrak{ D }'$ is implied by the relation $\# {\cal S}(G,R) \leq \# {\cal S}(G,S)$ for every $G \in \mathfrak{ D }'$, where ${\cal S}(G,H)$ is the set of homomorphisms from $G$ to $H$ mapping all proper arcs of $G$ to proper arcs of $H$. Under an application-oriented regularity condition, the two relations are even equivalent. A method is developed for the rearrangement of a digraph $R$, resulting in a digraph $S$ with $\# {\cal H}(G,R) \leq \# {\cal H}(G,S)$ for every digraph $G$. The method is applied in constructing pairs of partially ordered sets $R$ and $S$ with $\# {\cal H}(P,R) \leq \# {\cal H}(P,S)$ for every partially ordered set $P$. The main part of the results holds also for undirected graphs.


title: The Generalized OTOC from Supersymmetric Quantum Mechanics: Study of Random Fluctuations from Eigenstate Representation of Correlation Functions

abstract: The concept of out-of-time-ordered correlation (OTOC) function is treated as a very strong theoretical probe of quantum randomness, using which one can study both chaotic and non-chaotic phenomena in the context of quantum statistical mechanics. In this paper, we define a general class of OTOC, which can perfectly capture quantum randomness phenomena in a better way. Further we demonstrate an equivalent formalism of computation using a general time independent Hamiltonian having well defined eigenstate representation for integrable supersymmetric quantum systems. We found that one needs to consider two new correlators apart from the usual one to have a complete quantum description. To visualize the impact of the given formalism we consider the two well known models viz. Harmonic Oscillator and one dimensional potential well within the framework of supersymmetry. For the Harmonic Oscillator case we obtain similar periodic time dependence but dissimilar parameter dependences compared to the results obtained from both micro-canonical and canonical ensembles in quantum mechanics without supersymmetry. On the other hand, for one dimensional potential well problem we found significantly different time scale and the other parameter dependence compared to the results obtained from non-supersymmetric quantum mechanics. Finally, to establish the consistency of the prescribed formalism in the classical limit, we demonstrate the phase space averaged version of the classical version of OTOCs from a model independent Hamiltonian along with the previously mentioned these well cited models.

categories: hep-th cond-mat.dis-nn gr-qc nlin.CD quant-ph


title: Scanning electron diffraction tomography of strain

abstract: Strain engineering is used to obtain desirable materials properties in a range of modern technologies. Direct nanoscale measurement of the three-dimensional strain tensor field within these materials has however been limited by a lack of suitable experimental techniques and data analysis tools. Scanning electron diffraction has emerged as a powerful tool for obtaining two-dimensional maps of strain components perpendicular to the incident electron beam direction. Extension of this method to recover the full three-dimensional strain tensor field has been restricted though by the absence of a formal framework for tensor tomography using such data. Here, we show that it is possible to reconstruct the full non-symmetric strain tensor field as the solution to an ill-posed tensor tomography inverse problem. We then demonstrate the properties of this tomography problem both analytically and computationally, highlighting why incorporating precession to perform scanning precession electron diffraction may be important. We establish a general framework for non-symmetric tensor tomography and demonstrate computationally its applicability for achieving strain tomography with scanning precession electron diffraction data.


title: Empirical Likelihood Estimation for Linear Regression Models with AR(p) Error Terms

abstract: Linear regression models are useful statistical tools to analyze data sets in several different fields. There are several methods to estimate the parameters of a linear regression model. These methods usually perform under normally distributed and uncorrelated errors with zero mean and constant variance. However, for some data sets error terms may not satisfy these or some of these assumptions. If error terms are correlated, such as the regression models with autoregressive (AR(p)) error terms, the Conditional Maximum Likelihood (CML) under normality assumption or the Least Square (LS) methods are often used to estimate the parameters of interest. For CML estimation a distributional assumption on error terms is needed to carry on estimation, but, in practice, such distributional assumptions on error terms may not be plausible. Therefore, in such cases some alternative distribution free methods are needed to conduct the parameter estimation. In this paper, we propose to estimate the parameters of a linear regression model with AR(p) error term using the Empirical Likelihood (EL) method, which is one of the distribution free estimation methods. A small simulation study and a numerical example are provided to evaluate the performance of the proposed estimation method over the CML method. The results of simulation study show that the proposed estimators based on EL method are remarkably better than the estimators obtained from the CML method in terms of mean squared errors (MSE) and bias in almost all the simulation configurations. These findings are also confirmed by the results of the numerical and real data examples.


title: COVID-19: What If Immunity Wanes?

abstract: Using a simple economic model in which social-distancing reduces contagion, we study the implications of waning immunity for the epidemiological dynamics and social activity. If immunity wanes, we find that COVID-19 likely becomes endemic and that social-distancing is here to stay until the discovery of a vaccine or cure. But waning immunity does not necessarily change optimal actions on the onset of the pandemic. Decentralized equilibria are virtually independent of waning immunity until close to peak infections. For centralized equilibria, the relevance of waning immunity decreases in the probability of finding a vaccine or cure, the costs of infection (e.g., infection-fatality rate), and the presence of other NPIs that lower contagion (e.g., quarantining and mask use). In simulations calibrated to July 2020, our model suggests that waning immunity is virtually unimportant for centralized equilibria until at least 2021. This provides vital time for individuals and policymakers to learn about immunity against SARS-CoV-2 before it becomes critical.

categories: econ.GN physics.soc-ph q-bio.PE q-fin.EC


title: More than softer-when-brighter: the X-Ray powerlaw spectral variability in NGC 4051

abstract: The powerlaw X-ray spectra of active galactic nuclei at moderate to high accretion rates normally appear softer when they brighten, for which the underlying mechanisms are yet unclear. Utilizing XMM-Newton observations and excluding photons $<$ 2 keV to avoid contamination from the soft excess, in this work we scrutinize the powerlaw spectral variability of NCG 4051 from two new aspects. We first find that a best-fit "softer-when-brighter" relation is statistically insufficient to explain the observed spectral variabilities, and intervals deviated from the empirical relation are clearly visible in the light curve of 2 -- 4 keV/4 -- 10 keV count rate ratio. The deviations are seen not only between but also within individual XMM-Newton exposures, consistent with random variations of the corona geometry or inner structure (with timescales as short as $\sim$ 1 ks), in addition to those behind the smooth "softer-when-brighter" trend. We further find the "softer-when-brighter" trend gradually weakens with the decreasing timescale (from $\sim$ 100 ks down to 0.5 ks). These findings indicate that the powerlaw spectral slope is not solely determined by its brightness. We propose a two-tier geometry, including flares/nano-flares on top of the inner disc and an embedding extended corona (heated by the flares, in analogy to solar corona) to explain the observations together with other observational clues in literature. Rapid spectral variabilities could be due to individual flares/nano-flares, while slow ones are driven by the variations in the global activity of inner disc region (akin to the variation of solar activity, but not the accretion rate) accompanied with heating/cooling and inflation/contraction of the extended corona.


title: Physics-Based Dexterous Manipulations with Estimated Hand Poses and Residual Reinforcement Learning

abstract: Dexterous manipulation of objects in virtual environments with our bare hands, by using only a depth sensor and a state-of-the-art 3D hand pose estimator (HPE), is challenging. While virtual environments are ruled by physics, e.g. object weights and surface frictions, the absence of force feedback makes the task challenging, as even slight inaccuracies on finger tips or contact points from HPE may make the interactions fail. Prior arts simply generate contact forces in the direction of the fingers' closures, when finger joints penetrate virtual objects. Although useful for simple grasping scenarios, they cannot be applied to dexterous manipulations such as in-hand manipulation. Existing reinforcement learning (RL) and imitation learning (IL) approaches train agents that learn skills by using task-specific rewards, without considering any online user input. In this work, we propose to learn a model that maps noisy input hand poses to target virtual poses, which introduces the needed contacts to accomplish the tasks on a physics simulator. The agent is trained in a residual setting by using a model-free hybrid RL+IL approach. A 3D hand pose estimation reward is introduced leading to an improvement on HPE accuracy when the physics-guided corrected target poses are remapped to the input space. As the model corrects HPE errors by applying minor but crucial joint displacements for contacts, this helps to keep the generated motion visually close to the user input. Since HPE sequences performing successful virtual interactions do not exist, a data generation scheme to train and evaluate the system is proposed. We test our framework in two applications that use hand pose estimates for dexterous manipulations: hand-object interactions in VR and hand-object motion reconstruction in-the-wild.


title: HoliCity: A City-Scale Data Platform for Learning Holistic 3D Structures

abstract: We present HoliCity, a city-scale 3D dataset with rich structural information. Currently, this dataset has 6,300 real-world panoramas of resolution $13312 \times 6656$ that are accurately aligned with the CAD model of downtown London with an area of more than 20 km$^2$, in which the median reprojection error of the alignment of an average image is less than half a degree. This dataset aims to be an all-in-one data platform for research of learning abstracted high-level holistic 3D structures that can be derived from city CAD models, e.g., corners, lines, wireframes, planes, and cuboids, with the ultimate goal of supporting real-world applications including city-scale reconstruction, localization, mapping, and augmented reality. The accurate alignment of the 3D CAD models and panoramas also benefits low-level 3D vision tasks such as surface normal estimation, as the surface normal extracted from previous LiDAR-based datasets is often noisy. We conduct experiments to demonstrate the applications of HoliCity, such as predicting surface segmentation, normal maps, depth maps, and vanishing points, as well as test the generalizability of methods trained on HoliCity and other related datasets. HoliCity is available at https://holicity.io.


title: One idea and two proofs of the KMT theorems

abstract: Two proofs of the Koml\'os-Major-Tusn\'ady embedding theorems, one for the uniform empirical process and one for the simple symmetric random walk, are given. More precisely, what are proved are the univariate coupling results needed in the proofs, such as Tusn\'{a}dy's lemma. These proofs are modifications of existing proof architectures, one combinatorial (the original proof with many modifications, due to Cs\"{o}rg\~o, R\'{e}v\'{e}sz, Bretagnolle, Massart, Dudley, Carter, Pollard etc.) and one analytical (due to Sourav Chatterjee). There is one common idea to both proofs: we compare binomial and hypergeometric distributions among themselves, rather than with the Gaussian distribution. In the combinatorial approach, this involves comparing Binomial(n,1/2) distribution with the Binomial(4n,1/2) distribution, which mainly involves comparison between the corresponding binomial coefficients. In the analytical approach, this reduces Chatterjee's method to coupling nearest neighbour Markov chains on integers so that they stay close.


title: Rejoinder: On nearly assumption-free tests of nominal confidence interval coverage for causal parameters estimated by machine learning

abstract: This is the rejoinder to the discussion by Kennedy, Balakrishnan and Wasserman on the paper "On nearly assumption-free tests of nominal confidence interval coverage for causal parameters estimated by machine learning" published in Statistical Science.


title: The power spectrum on small scales: Robust constraints and comparing PBH methodologies

abstract: We compare primordial black hole (PBH) constraints on the power spectrum and mass distributions using the traditional Press Schechter formalism, peaks theory, and a recently developed version of peaks theory relevant to PBHs. We show that, provided the PBH formation criteria and the power spectrum smoothing are treated consistently, the constraints only vary by $\sim$10\% between methods (a difference that will become increasingly important with better data). Our robust constraints from PBHs take into account the effects of critical collapse, the non-linear relation between $\zeta$ and $\delta$, and the shift from the PBH mass to the power spectrum peak scale. We show that these constraints are remarkably similar to the pulsar timing array (PTA) constraints impacting the black hole masses detected by the LIGO and Virgo, but that the $\mu$-distortion constraints rule out supermassive black hole (SMBH) formation and potentially even the much lighter mass range of $\sim$(1-100) $\mathrm{M}_\odot$ that LIGO/Virgo probes.


title: A Novel Tampering Attack on AES Cores with Hardware Trojans

abstract: The implementation of cryptographic primitives in integrated circuits (ICs) continues to increase over the years due to the recent advancement of semiconductor manufacturing and reduction of cost per transistors. The hardware implementation makes cryptographic operations faster and more energy-efficient. However, various hardware attacks have been proposed aiming to extract the secret key in order to undermine the security of these primitives. In this paper, we focus on the widely used advanced encryption standard (AES) block cipher and demonstrate its vulnerability against tampering attack. Our proposed attack relies on implanting a hardware Trojan in the netlist by an untrusted foundry, which can design and implement such a Trojan as it has access to the design layout and mask information. The hardware Trojan's activation modifies a particular round's input data by preventing the effect of all previous rounds' key-dependent computation. We propose to use a sequential hardware Trojan to deliver the payload at the input of an internal round for achieving this modification of data. All the internal subkeys, and finally, the secret key can be computed from the observed ciphertext once the Trojan is activated. We implement our proposed tampering attack with a sequential hardware Trojan inserted into a 128-bit AES design from OpenCores benchmark suite and report the area overhead to demonstrate the feasibility of the proposed tampering attack.


title: Nuclear and dark matter heating in massive white dwarf stars

abstract: Recently, Cheng et al. identified a number of massive white dwarfs (WD) that appear to have an additional heat source providing a luminosity near $\approx 10^{-3}L_\odot$ for multiple Gyr. In this paper we explore heating from electron capture and pycnonuclear reactions. We also explore heating from dark matter annihilation. WD stars appear to be too small to capture enough dark matter for this to be important. Finally, if dark matter condenses to very high densities inside a WD this could ignite nuclear reactions. We calculate the enhanced central density of a WD in the gravitational potential of a very dense dark matter core. While this might start a supernova, it seems unlikely to provide modest heating for a long time. We conclude that electron capture, pycnonuclear, and dark matter reactions are unlikely to provide significant heating in the massive WD that Cheng considers.


title: Foatic actions of the symmetric group and fixed-point homomesy

abstract: We study maps on the set of permutations of n generated by the R\'enyi-Foata map intertwined with other dihedral symmetries (of a permutation considered as a 0-1 matrix). Iterating these maps leads to dynamical systems that in some cases exhibit interesting orbit structures, e.g., every orbit size being a power of two, and homomesic statistics (ones which have the same average over each orbit). In particular, the number of fixed points (aka 1-cycles) of a permutation appears to be homomesic with respect to three of these maps, even in one case where the orbit structures are far from nice. For the most interesting such "Foatic" action, we give a heap analysis and recursive structure that allows us to prove the fixed-point homomesy and orbit properties, but two other cases remain conjectural.


title: NuSTAR hard X-ray spectra of radio galaxies

abstract: NuSTAR observatory, with its 3 - 78 keV broadband spectral coverage, enables the detections of the high-energy cutoff in a number of active galaxies, including several individual radio loud ones. In this work we present systematic and uniform analyses of 55 NuSTAR spectra for a large sample of 28 radio galaxies, 20 of which are FR II galaxies. We perform spectral fitting to measure the high energy cut-off $E_{cut}$, photon index $\Gamma$, reflection factor R and Fe K$\alpha$ line equivalent width. Measurements of $E_{cut}$ are given for 13 sources, and lower limits for the rest. We find those $E_{cut}$ non-detections could primarily be attributed to the obviously smaller net photon counts in their spectra. This indicates that the NuSTAR spectra of the majority of our sample are dominated by the thermal corona emission, and the $E_{cut}$ distribution of the sample is indistinguishable from that of a radio quiet one in literature. The flatter NuSTAR spectra we observed, comparing with radio quiet sources, are thus unlikely due to jet contamination. The radio galaxies also show weaker X-ray reflection (both in R and Fe K$\alpha$ line EW) comparing with radio quiet ones. Combining with the radio quiet sample we see a correlation between R and EW, but with considerably large scatter. Notably, the radio loud and quiet sources appear to follow a common $\Gamma$ - R correlation trend, supporting the outflowing corona model for both populations in which higher bulk outflowing velocity yields weaker reflection and flatter X-ray slope.


title: Active Brownian Motion in two-dimensions under Stochastic Resetting

abstract: We study the position distribution of an active Brownian particle (ABP) in the presence of stochastic resetting in two spatial dimensions. We consider three different resetting protocols : (I) where both position and orientation of the particle are reset, (II) where only the position is reset, and (III) where only the orientation is reset with a certain rate $r.$ We show that in the first two cases the ABP reaches a stationary state. Using a renewal approach, we calculate exactly the stationary marginal position distributions in the limiting cases when the resetting rate $r$ is much larger or much smaller than the rotational diffusion constant $D_R$ of the ABP. We find that, in some cases, for a large resetting rate, the position distribution diverges near the resetting point; the nature of the divergence depends on the specific protocol. For the orientation resetting, there is no stationary state, but the motion changes from a ballistic one at short-times to a diffusive one at late times. We characterize the short-time non-Gaussian marginal position distributions using a perturbative approach.


title: Low Scale String Theory Benchmarks for Hidden Photon Dark Matter Interpretations of the XENON1T Anomaly

abstract: An excess of low-energy electronic recoil events over known backgrounds was recently observed in the XENON1T detector, where $285$ events are observed compared to an expected $232 \pm 15$ events from the background-only fit to the data in the energy range 1-7 keV. This could be due to the beta decay of an unexpected tritium component, or possibly to new physics. One plausible new physics explanation for the excess is absorption of hidden photon dark matter relics with mass around $2.8$ keV and kinetic mixing of about $10^{-15}$, which can also explain cooling excesses in horizontal-branch (HB) stars. Such small gauge boson masses and couplings can naturally arise from type-IIB low scale string theory. We provide a fit of the XENON1T excess in terms of a minimal low scale type-IIB string theory parameter space and present some benchmark points which provide a good fit to the data. It is also demonstrated how the required transformation properties of the massless spectrum are obtained in intersecting D-brane models.


title: On equations over direct powers of algebraic structures

abstract: We study systems of equations over graphs, posets and matroids. We give the criteria, when a direct power of such algebraic structures is equationally Noetherian. Moreover we prove that any direct power of a finite algebraic structure is weakly equationally Noetherian.


title: Multi-Stage Optimized Machine Learning Framework for Network Intrusion Detection

abstract: Cyber-security garnered significant attention due to the increased dependency of individuals and organizations on the Internet and their concern about the security and privacy of their online activities. Several previous machine learning (ML)-based network intrusion detection systems (NIDSs) have been developed to protect against malicious online behavior. This paper proposes a novel multi-stage optimized ML-based NIDS framework that reduces computational complexity while maintaining its detection performance. This work studies the impact of oversampling techniques on the models' training sample size and determines the minimal suitable training sample size. Furthermore, it compares between two feature selection techniques, information gain and correlation-based, and explores their effect on detection performance and time complexity. Moreover, different ML hyper-parameter optimization techniques are investigated to enhance the NIDS's performance. The performance of the proposed framework is evaluated using two recent intrusion detection datasets, the CICIDS 2017 and the UNSW-NB 2015 datasets. Experimental results show that the proposed model significantly reduces the required training sample size (up to 74%) and feature set size (up to 50%). Moreover, the model performance is enhanced with hyper-parameter optimization with detection accuracies over 99% for both datasets, outperforming recent literature works by 1-2% higher accuracy and 1-2% lower false alarm rate.


title: FitsGeo -- Python package for PHITS geometry development and visualization

abstract: An easy way to define and visualize geometry for PHITS input files introduced. Suggested FitsGeo Python package helps to define surfaces as Python objects and manipulate them conveniently. VPython assists to view defined geometry interactively which boosts geometry development and helps with complicated cases. Every class that sets the surface object has methods with some extra properties. As well as geometry generation for PHITS input, additional modules developed for material and cell definition. Any user with a very basic knowledge of Python can define the geometry in a convenient way and use it in further research related to particle transport.


title: Topology in cyber research

abstract: We give an idiosyncratic overview of applications of topology to cyber research, spanning the analysis of variables/assignments and control flow in computer programs, a brief sketch of topological data analysis in one dimension, and the use of sheaves to analyze wireless networks. The text is from a chapter in the forthcoming book Mathematics in Cyber Research, to be published by Taylor and Francis.


title: Electrodynamics in flat spacetime of six dimensions

abstract: We consider the dynamics of a classical charge in flat spacetime of six dimensions. The mass shell relation of a free charge admits nonlinear oscillations. Having analyzed the problem of on eigenvalues and eigenvectors of Faraday tensor, we establish the algebraic structure of electromagnetic field in 6D. We elaborate the classification scheme based on three field's invariants. Using the basic algebraic properties of the electromagnetic field tensor we analyze the motion of a charge in constant electromagnetic field. Its world line is a combination of hyperbolic and circular orbits which lie in three mutually orthogonal sheets of two dimensions. Within the braneworld scenario, we project the theory on the de Sitter space of four dimensions. Actually, as it turns out, spins of elementary particles themselves are manifestations of extra dimensions.


title: White-box Induction From SVM Models: Explainable AI with Logic Programming

abstract: We focus on the problem of inducing logic programs that explain models learned by the support vector machine (SVM) algorithm. The top-down sequential covering inductive logic programming (ILP) algorithms (e.g., FOIL) apply hill-climbing search using heuristics from information theory. A major issue with this class of algorithms is getting stuck in a local optimum. In our new approach, however, the data-dependent hill-climbing search is replaced with a model-dependent search where a globally optimal SVM model is trained first, then the algorithm looks into support vectors as the most influential data points in the model, and induces a clause that would cover the support vector and points that are most similar to that support vector. Instead of defining a fixed hypothesis search space, our algorithm makes use of SHAP, an example-specific interpreter in explainable AI, to determine a relevant set of features. This approach yields an algorithm that captures SVM model's underlying logic and outperforms %GG: the FOIL algorithm --> other ILP algorithms other ILP algorithms in terms of the number of induced clauses and classification evaluation metrics. This paper is under consideration for publication in the journal of "Theory and practice of logic programming".


title: The operational reality of quantum nonlocality

abstract: Does the remote measurement-disturbance of the quantum state of a system $B$ by measurement on system $A$ entangled with $B$, constitute a real disturbance -- i.e., an objective alteration -- of $B$ in an operational sense? Employing information theoretic criteria motivated by operational considerations alone, we argue that the disturbance in question is real for a subset of steerable correlations. This result highlights the distinction between quantum no-signaling and relativistic signal-locality. It furthermore suggests a natural reason why a convex operational theory should be non-signaling: namely, to ensure the consistency between the properties of reduced systems and those of single systems.


title: The Automatic Learning for the Rapid Classification of Events (ALeRCE) Alert Broker

abstract: We introduce the Automatic Learning for the Rapid Classification of Events (ALeRCE) broker, an astronomical alert broker designed to provide a rapid and self--consistent classification of large etendue telescope alert streams, such as that provided by the Zwicky Transient Facility (ZTF) and, in the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). ALeRCE is a Chilean--led broker run by an interdisciplinary team of astronomers and engineers, working to become intermediaries between survey and follow--up facilities. ALeRCE uses a pipeline which includes the real--time ingestion, aggregation, cross--matching, machine learning (ML) classification, and visualization of the ZTF alert stream. We use two classifiers: a stamp--based classifier, designed for rapid classification, and a light--curve--based classifier, which uses the multi--band flux evolution to achieve a more refined classification. We describe in detail our pipeline, data products, tools and services, which are made public for the community (see \url{https://alerce.science}). Since we began operating our real--time ML classification of the ZTF alert stream in early 2019, we have grown a large community of active users around the globe. We describe our results to date, including the real--time processing of $9.7\times10^7$ alerts, the stamp classification of $1.9\times10^7$ objects, the light curve classification of $8.5\times10^5$ objects, the report of 3088 supernova candidates, and different experiments using LSST-like alert streams. Finally, we discuss the challenges ahead to go from a single-stream of alerts such as ZTF to a multi--stream ecosystem dominated by LSST.


title: Bipolaron liquids at strong Peierls electron-phonon couplings

abstract: We use the Density Matrix Renormalization Group method to study a one-dimensional chain with Peierls electron-phonon coupling, which describes the modulation of the electron hopping by lattice distortions. We demonstrate that this system is stable against phase separation in the dilute density limit. We only find phase separation numerically for large couplings for which the linear approximation for the electron-phonon coupling becomes invalid; this behavior can be stabilized in a narrow sliver of the physical parameter space if the dispersion of the phonons is carefully tuned. These results indicate that in the dilute electron density limit, Peierls bipolaron liquids are generically stable, unlike in other models of electron-phonon coupling. We show that this behavior extends to finite carrier concentrations of up to quarter filling. This stability of low-density, light-mass bipolaron liquids in the Peierls model opens a path to high-$T_c$ superconductivity based on a bipolaronic mechanism, in higher dimensions.


title: X-ray Searches for Axions from Super Star Clusters

abstract: Axions may be produced in abundance inside stellar cores and then convert into observable X-rays in the Galactic magnetic fields. We focus on the Quintuplet and Westerlund 1 super star clusters, which host large numbers of hot, young stars including Wolf-Rayet stars; these stars produce axions efficiently through the axion-photon coupling. We use Galactic magnetic field models to calculate the expected X-ray flux locally from axions emitted from these clusters. We then combine the axion model predictions with archival Nuclear Spectroscopic Telescope Array (NuSTAR) data from 10 - 80 keV to search for evidence of axions. We find no significant evidence for axions and constrain the axion-photon coupling $g_{a\gamma\gamma} \lesssim 3.6 \times 10^{-12}$ GeV$^{-1}$ for masses $m_a \lesssim 5 \times 10^{-11}$ eV at 95\% confidence.


title: Two-Dimensional Simulations of Internal Gravity Waves in The Radiation Zones of Intermediate-Mass Stars

abstract: Intermediate-mass main sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective-radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. We have studied the propagation of these IGWs using two-dimensional fully-non-linear hydrodynamical simulations with realistic stellar reference states from the one-dimensional stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations which can be an effect of nonlinearities. When the same generation spectrum is applied to 3 solar mass models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.


title: Shortcuts to Squeezed Thermal States

abstract: Squeezed state in harmonic systems can be generated through a variety of techniques, including varying the oscillator frequency or using nonlinear two-photon Raman interaction. We focus on these two techniques to drive an initial thermal state into a final squeezed thermal state with controlled squeezing parameters -- amplitude and phase -- in arbitrary time. The protocols are designed through reverse engineering for both unitary and open dynamics. Control of the dissipation is achieved using stochastic processes, readily implementable via, e.g., continuous quantum measurements. Importantly, this allows controlling the state entropy and can be used for fast thermalization. The developed protocols are thus suited to generate squeezed thermal states at controlled temperature in arbitrary time.


title: Selected Topics in the Large Quantum Number Expansion

abstract: In this review we study quantum field theories and conformal field theories with global symmetries in the limit of large charge for some of the generators of the symmetry group. At low energy the sectors of the theory with large charge are described by a hybrid form of Goldstone's theorem, involving its relativistic and non-relativistic forms. The associated effective field theory in the infrared allows the computation of anomalous dimensions, and operator product expansion coefficients in a well defined expansion in inverse powers of the global charge. This applies even when the initial theory does not have a reliable semiclassical approximation. The large quantum number expansion complements, and may provide an alternative approach to the bootstrap and numerical treatments. We will present some general features of the symmetry breaking patterns and the low-energy effective actions, and a fairly large number of examples exhibiting the salient features of this method.


title: Alert Classification for the ALeRCE Broker System: The Real-time Stamp Classifier

abstract: We present a real-time stamp classifier of astronomical events for the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker. The classifier is based on a convolutional neural network, trained on alerts ingested from the Zwicky Transient Facility (ZTF). Using only the \textit{science, reference} and \textit{difference} images of the first detection as inputs, along with the metadata of the alert as features, the classifier is able to correctly classify alerts from active galactic nuclei, supernovae (SNe), variable stars, asteroids and bogus classes, with high accuracy ($\sim$94\%) in a balanced test set. In order to find and analyze SN candidates selected by our classifier from the ZTF alert stream, we designed and deployed a visualization tool called SN Hunter, where relevant information about each possible SN is displayed for the experts to choose among candidates to report to the Transient Name Server database. From June 26th 2019 to February 28th 2021, we have reported 6846 SN candidates to date (11.8 candidates per day on average), of which 971 have been confirmed spectroscopically. Our ability to report objects using only a single detection means that 70\% of the reported SNe occurred within one day after the first detection. ALeRCE has only reported candidates not otherwise detected or selected by other groups, therefore adding new early transients to the bulk of objects available for early follow-up. Our work represents an important milestone toward rapid alert classifications with the next generation of large etendue telescopes, such as the Vera C. Rubin Observatory.


title: Synchronization and Non-Markovianity in open quantum systems

abstract: Detuned systems can spontaneously achieve a synchronous dynamics and display robust quantum correlations in different local and global dissipation regimes. Beyond the Markovian limit, information backflow from the environment becomes a crucial mechanism whose interplay with spontaneous synchronization is unknown. Considering a model of two coupled qubits, one of which interacts with a dissipative environment, we show that non-Markovianity is highly detrimental for the emergence of synchronization, for the latter can be delayed and hindered because of the presence of information backflow. The results are obtained considering both a master equation approach and a collision model based on repeated interactions, which represents a very versatile tool to tailor the desired kind of environment.


title: Alert Classification for the ALeRCE Broker System: The Light Curve Classifier

abstract: We present the first version of the ALeRCE (Automatic Learning for the Rapid Classification of Events) broker light curve classifier. ALeRCE is currently processing the Zwicky Transient Facility (ZTF) alert stream, in preparation for the Vera C. Rubin Observatory. The ALeRCE light curve classifier uses variability features computed from the ZTF alert stream, and colors obtained from AllWISE and ZTF photometry. We apply a Balanced Random Forest algorithm with a two-level scheme, where the top level classifies each source as periodic, stochastic, or transient, and the bottom level further resolves each of these hierarchical classes, amongst 15 total classes. This classifier corresponds to the first attempt to classify multiple classes of stochastic variables (including core- and host-dominated active galactic nuclei, blazars, young stellar objects, and cataclysmic variables) in addition to different classes of periodic and transient sources, using real data. We created a labeled set using various public catalogs (such as the Catalina Surveys and {\em Gaia} DR2 variable stars catalogs, and the Million Quasars catalog), and we classify all objects with $\geq6$ $g$-band or $\geq6$ $r$-band detections in ZTF (868,371 sources as of 2020/06/09), providing updated classifications for sources with new alerts every day. For the top level we obtain macro-averaged precision and recall scores of 0.96 and 0.99, respectively, and for the bottom level we obtain macro-averaged precision and recall scores of 0.57 and 0.76, respectively. Updated classifications from the light curve classifier can be found at the \href{http://alerce.online}{ALeRCE Explorer website}.

categories: astro-ph.IM astro-ph.GA astro-ph.HE astro-ph.SR


title: Complete parameter inference for GW150914 using deep learning

abstract: The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-network conditional density estimator to model posterior probability distributions over the full 15-dimensional space of binary black hole system parameters, given detector strain data from multiple detectors. We use the method of normalizing flows---specifically, a neural spline normalizing flow---which allows for rapid sampling and density estimation. Training the network is likelihood-free, requiring samples from the data generative process, but no likelihood evaluations. Through training, the network learns a global set of posteriors: it can generate thousands of independent posterior samples per second for any strain data consistent with the prior and detector noise characteristics used for training. By training with the detector noise power spectral density estimated at the time of GW150914, and conditioning on the event strain data, we use the neural network to generate accurate posterior samples consistent with analyses using conventional sampling techniques.


title: Observation of eccentric binary black hole mergers with second and third generation gravitational wave detector networks

abstract: [Abridged] We introduce an improved version of the Eccentric, Non-spinning, Inspiral-Gaussian-process Merger Approximant (ENIGMA) waveform model. We find that this ready-to-use model can: (i) produce physically consistent signals when sampling over 1M samples chosen over the $m_{\{1,\,2\}}\in[5M_{\odot},\,50M_{\odot}]$ parameter space, and the entire range of binary inclination angles; (ii) produce waveforms within 0.04 seconds from an initial gravitational wave frequency $f_{\textrm{GW}} =15\,\textrm{Hz}$ and at a sample rate of 8192 Hz; and (iii) reproduce the physics of quasi-circular mergers. We utilize ENIGMA to compute the expected signal-to-noise ratio (SNR) distributions of eccentric binary black hole mergers assuming the existence of second and third generation gravitational wave detector networks that include the twin LIGO detectors, Virgo, KAGRA, LIGO-India, a LIGO-type detector in Australia, Cosmic Explorer, and the Einstein Telescope. In the context of advanced LIGO-type detectors, we find that the SNR of eccentric mergers is always larger than quasi-circular mergers for systems with $e_0\leq0.4$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$, even if the timespan of eccentric signals is just a third of quasi-circular systems with identical total mass and mass-ratio. For Cosmic Explorer-type detector networks, we find that eccentric mergers have similar SNRs than quasi-circular systems for $e_0\leq0.3$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$. Systems with $e_0\sim0.5$ at $f_{\textrm{GW}} =10\,\textrm{Hz}$ have SNRs that range between 50%-90% of the SNR produced by quasi-circular mergers, even if these eccentric signals are just between a third to a tenth the length of quasi-circular systems. For Einstein Telescope-type detectors, we find that eccentric mergers have similar SNRs than quasi-circular systems for $e_0\leq0.4$ at $f_{\textrm{GW}} =5\,\textrm{Hz}$.


title: Spectro-polarimetric observations of the CIZA J2242.8+5301 northern radio relic: no evidence of high-frequency steepening

abstract: Observations of radio relics at very high frequency (>10 GHz) can help to understand how particles age and are (re-)accelerated in galaxy cluster outskirts and how magnetic fields are amplified in these environments. In this work, we present new single-dish 18.6 GHz Sardinia Radio Telescope and 14.25 GHz Effelsberg observations of the well known northern radio relic of CIZA J2242.8+5301. We detected the relic which shows a length of $\sim$1.8 Mpc and a flux density equal to $\rm S_{14.25\,GHz}=(9.5\pm3.9)\,mJy$ and $\rm S_{18.6\,GHz}=(7.67\pm0.90)\,mJy$ at 14.25 GHz and 18.6 GHz respectively. The resulting best-fit model of the relic spectrum from 145 MHz to 18.6 GHz is a power-law spectrum with spectral index $\alpha=1.12\pm0.03$: no evidence of steepening has been found in the new data presented in this work. For the first time, polarisation properties have been derived at 18.6 GHz, revealing an averaged polarisation fraction of $\sim40\%$ and a magnetic field aligned with the 'filaments' or 'sheets' of the relic.


title: 3D Network Model for Strong Topological Insulator Transitions

abstract: We construct a three-dimensional (3D), time-reversal symmetric generalization of the Chalker-Coddington network model for the integer quantum Hall transition. The novel feature of our network model is that in addition to a weak topological insulator phase already accommodated by the network model framework in the pre-existing literature, it hosts strong topological insulator phases as well. We unambiguously demonstrate that strong topological insulator phases emerge as intermediate phases between a trivial insulator phase and a weak topological phase. Additionally, we found a non-local transformation that relates a trivial insulator phase and a weak topological phase in our network model. Remarkably, strong topological phases are mapped to themselves under this transformation. We show that upon adding sufficiently strong disorder the strong topological insulator phases undergo phase transitions into a metallic phase. We numerically determine the critical exponent of the insulator-metal transition. Our network model explicitly shows how a semi-classical percolation picture of topological phase transitions in 2D can be generalized to 3D and opens up a new venue for studying 3D topological phase transitions.


title: Controlling a quantum point junction on the surface of an antiferromagnetic topological insulator

abstract: The abstract notion of topology has led to profound insights into real materials. Notably, the surface and edges of topological materials can host physics, such as unidirectional charge or spin transport, that is unavailable in isolated one- and two-dimensional systems. However, to fully control the mixing and interference of edge-state wave functions, one needs robust and tunable junctions. We propose to achieve this control using an antiferromagnetic topological insulator that supports two distinct types of gapless unidirectional channels on its surface, one from antiferromagnetic domain walls and the other from single-height steps. The distinct geometric nature of these edge modes allows them to intersect robustly to form quantum point junctions, and their presence at the surface makes them subject to control by magnetic and electrostatic tips like those used in scanning probe microscopes. Prospects for realizing such junctions are encouraged by recent material candidate proposals, potentially leading to exciting applications in quantum computing and sensing.


title: Relativistic Bell Test within Quantum Reference Frames

abstract: A still widely debated question in the field of relativistic quantum information is whether entanglement and the degree of violation of Bell's inequalities for massive relativistic particles are frame independent or not. At the core of this question is the effect that spin gets entangled with the momentum degree of freedom at relativistic velocities. Here, we show that Bell's inequalities for a pair of particles can be maximally violated in a special-relativistic regime, even without any post-selection of the momentum of the particles. To this end, we use the methodology of quantum reference frames, which allows us to transform the problem to the rest frame of a particle, whose state can be in a superposition of relativistic momenta from the viewpoint of the laboratory frame. We show that, when the relative motion of two particles is non-collinear, the optimal measurements for violation of Bell's inequalities in the laboratory frame involve "coherent Wigner rotations". Moreover, the degree of violation of Bell's inequalities is independent of the choice of the quantum reference frame. Our results open up the possibility of extending entanglement-based quantum communication protocols to relativistic regimes.


title: Point Spectrum of Periodic Operators on Universal Covering Trees

abstract: For any multi-graph $G$ with edge weights and vertex potential, and its universal covering tree $\mathcal{T}$, we completely characterize the point spectrum of operators $A_{\mathcal{T}}$ on $\mathcal{T}$ arising as pull-backs of local, self-adjoint operators $A_{G}$ on $G$. This builds on work of Aomoto, and includes an alternative proof of the necessary condition for point spectrum he derived in (Aomoto, 1991). Our result gives a finite time algorithm to compute the point spectrum of $A_{\mathcal{T}}$ from the graph $G$, and additionally allows us to show that this point spectrum is contained in the spectrum of $A_{G}$. Finally, we prove that typical pull-back operators have a spectral delocalization property: the set of edge weight and vertex potential parameters of $A_{G}$ giving rise to $A_{\mathcal{T}}$ with purely absolutely continuous spectrum is open and its complement has large codimension.


title: Leading order corrections to the quantum extremal surface prescription

abstract: We show that a na\"{i}ve application of the quantum extremal surface (QES) prescription can lead to paradoxical results and must be corrected at leading order. The corrections arise when there is a second QES (with strictly larger generalized entropy at leading order than the minimal QES), together with a large amount of highly incompressible bulk entropy between the two surfaces. We trace the source of the corrections to a failure of the assumptions used in the replica trick derivation of the QES prescription, and show that a more careful derivation correctly computes the corrections. Using tools from one-shot quantum Shannon theory (smooth min- and max-entropies), we generalize these results to a set of refined conditions that determine whether the QES prescription holds. We find similar refinements to the conditions needed for entanglement wedge reconstruction (EWR), and show how EWR can be reinterpreted as the task of one-shot quantum state merging (using zero-bits rather than classical bits), a task gravity is able to achieve optimally efficiently.


title: The chemical properties of the Milky Way's on-bar and off-bar regions: evidence for inhomogeneous star formation history in the bulge

abstract: Numerous studies of integrated starlight, stellar counts, and kinematics have confirmed that the Milky Way is a barred galaxy. However, far fewer studies have investigated the bar's stellar population properties, which carry valuable independent information regarding the bar's formation history. Here we conduct a detailed analysis of chemical abundance distributions ([Fe/H] and [Mg/Fe]) in the on-bar and off-bar regions to study the azimuthal variation of star formation history (SFH) in the inner Galaxy. We find that the on-bar and off-bar stars at Galactocentric radii 3 $< r_{\rm GC}<$ 5 kpc have remarkably consistent [Fe/H] and [Mg/Fe] distribution functions and [Mg/Fe]--[Fe/H] relation, suggesting a common SFH shared by the long bar and the disc. In contrast, the bar and disc at smaller radii (2 $< r_{\rm GC} <$ 3 kpc) show noticeable differences, with relatively more very metal-rich ([Fe/H]~0.4) stars but fewer solar abundance stars in the bar. Given the three-phase star formation history proposed for the inner Galaxy in Lian et al. (2020b), these differences could be explained by the off-bar disc having experienced either a faster early quenching process or recent metal-poor gas accretion. Vertical variations of the abundance distributions at small $r_{\rm GC}$ suggest a wider vertical distribution of low-$\alpha$ stars in the bar, which may serve as chemical evidence for vertical heating through the bar buckling process. The lack of such vertical variations outside the bulge may then suggest a lack of vertical heating in the long bar.


title: The Effective Action of Superrotation Modes

abstract: Starting from an analysis of four-dimensional asymptotically flat gravity in first order formulation, we show that superrotation reparametrization modes are governed by an Alekseev--Shatashvili action on the celestial sphere. This two-dimensional conformal theory describes spontaneous symmetry breaking of Virasoro superrotations together with the explicit symmetry breaking of more general Diff$(\mathcal{S}^2)$ superrotations. We arrive at this result by first reformulating the asymptotic field equations and symmetries of the radiative vacuum sector in terms of a Chern--Simons theory at null infinity, and subsequently performing a Hamiltonian reduction of this theory onto the celestial sphere.


title: scGNN: scRNA-seq Dropout Imputation via Induced Hierarchical Cell Similarity Graph

abstract: Single-cell RNA sequencing provides tremendous insights to understand biological systems. However, the noise from dropout can corrupt the downstream biological analysis. Hence, it is desirable to impute the dropouts accurately. In this work, we propose a simple and powerful dropout imputation method (scGNN) by applying a bottlenecked Graph Convolutional Neural Network on an induced hierarchical cell similarity graph. We show scGNN has competitive performance against state-of-the-art baselines across three datasets and can improve downstream analysis.


title: COVID-19 in differential diagnosis of online symptom assessments

abstract: The COVID-19 pandemic has magnified an already existing trend of people looking for healthcare solutions online. One class of solutions are symptom checkers, which have become very popular in the context of COVID-19. Traditional symptom checkers, however, are based on manually curated expert systems that are inflexible and hard to modify, especially in a quickly changing situation like the one we are facing today. That is why all COVID-19 existing solutions are manual symptom checkers that can only estimate the probability of this disease and cannot contemplate alternative hypothesis or come up with a differential diagnosis. While machine learning offers an alternative, the lack of reliable data does not make it easy to apply to COVID-19 either. In this paper we present an approach that combines the strengths of traditional AI expert systems and novel deep learning models. In doing so we can leverage prior knowledge as well as any amount of existing data to quickly derive models that best adapt to the current state of the world and latest scientific knowledge. We use the approach to train a COVID-19 aware differential diagnosis model that can be used for medical decision support both for doctors or patients. We show that our approach is able to accurately model new incoming data about COVID-19 while still preserving accuracy on conditions that had been modeled in the past. While our approach shows evident and clear advantages for an extreme situation like the one we are currently facing, we also show that its flexibility generalizes beyond this concrete, but very important, example.


title: Fisher Information Field: an Efficient and Differentiable Map for Perception-aware Planning

abstract: Considering visual localization accuracy at the planning time gives preference to robot motion that can be better localized and thus has the potential of improving vision-based navigation, especially in visually degraded environments. To integrate the knowledge about localization accuracy in motion planning algorithms, a central task is to quantify the amount of information that an image taken at a 6 degree-of-freedom pose brings for localization, which is often represented by the Fisher information. However, computing the Fisher information from a set of sparse landmarks (i.e., a point cloud), which is the most common map for visual localization, is inefficient. This approach scales linearly with the number of landmarks in the environment and does not allow the reuse of the computed Fisher information. To overcome these drawbacks, we propose the first dedicated map representation for evaluating the Fisher information of 6 degree-of-freedom visual localization for perception-aware motion planning. By formulating the Fisher information and sensor visibility carefully, we are able to separate the rotational invariant component from the Fisher information and store it in a voxel grid, namely the Fisher information field. This step only needs to be performed once for a known environment. The Fisher information for arbitrary poses can then be computed from the field in constant time, eliminating the need of costly iterating all the 3D landmarks at the planning time. Experimental results show that the proposed Fisher information field can be applied to different motion planning algorithms and is at least one order-of-magnitude faster than using the point cloud directly. Moreover,the proposed map representation is differentiable, resulting in better performance than the point cloud when used in trajectory optimization algorithms.


title: Stochastic Optimization and Learning for Two-Stage Supplier Problems

abstract: The main focus of this paper is radius-based (supplier) clustering in the two-stage stochastic setting with recourse, where the inherent stochasticity of the model comes in the form of a budget constraint. In addition to the standard (homogeneous) setting where all clients must be within a distance $R$ of the nearest facility, we provide results for the more general problem where the radius demands may be inhomogeneous (i.e., different for each client). We also explore a number of variants where additional constraints are imposed on the first-stage decisions, specifically matroid and multi-knapsack constraints, and provide results for these settings. We derive results for the most general distributional setting, where there is only black-box access to the underlying distribution. To accomplish this, we first develop algorithms for the polynomial scenarios setting; we then employ a novel scenario-discarding variant of the standard Sample Average Approximation (SAA) method, which crucially exploits properties of the restricted-case algorithms. We note that the scenario-discarding modification to the SAA method is necessary in order to optimize over the radius.


title: Optimal Combination of Linear and Spectral Estimators for Generalized Linear Models

abstract: We study the problem of recovering an unknown signal $\boldsymbol x$ given measurements obtained from a generalized linear model with a Gaussian sensing matrix. Two popular solutions are based on a linear estimator $\hat{\boldsymbol x}^{\rm L}$ and a spectral estimator $\hat{\boldsymbol x}^{\rm s}$. The former is a data-dependent linear combination of the columns of the measurement matrix, and its analysis is quite simple. The latter is the principal eigenvector of a data-dependent matrix, and a recent line of work has studied its performance. In this paper, we show how to optimally combine $\hat{\boldsymbol x}^{\rm L}$ and $\hat{\boldsymbol x}^{\rm s}$. At the heart of our analysis is the exact characterization of the joint empirical distribution of $(\boldsymbol x, \hat{\boldsymbol x}^{\rm L}, \hat{\boldsymbol x}^{\rm s})$ in the high-dimensional limit. This allows us to compute the Bayes-optimal combination of $\hat{\boldsymbol x}^{\rm L}$ and $\hat{\boldsymbol x}^{\rm s}$, given the limiting distribution of the signal $\boldsymbol x$. When the distribution of the signal is Gaussian, then the Bayes-optimal combination has the form $\theta\hat{\boldsymbol x}^{\rm L}+\hat{\boldsymbol x}^{\rm s}$ and we derive the optimal combination coefficient. In order to establish the limiting distribution of $(\boldsymbol x, \hat{\boldsymbol x}^{\rm L}, \hat{\boldsymbol x}^{\rm s})$, we design and analyze an Approximate Message Passing (AMP) algorithm whose iterates give $\hat{\boldsymbol x}^{\rm L}$ and approach $\hat{\boldsymbol x}^{\rm s}$. Numerical simulations demonstrate the improvement of the proposed combination with respect to the two methods considered separately.

categories: stat.ML cs.IT cs.LG math.IT math.ST stat.TH


title: A $4/3$-Approximation Algorithm for the Minimum $2$-Edge Connected Multisubgraph Problem in the Half-Integral Case

abstract: Given a connected undirected graph $\bar{G}$ on $n$ vertices, and non-negative edge costs $c$, the 2ECM problem is that of finding a $2$-edge~connected spanning multisubgraph of $\bar{G}$ of minimum cost. The natural linear program (LP) for 2ECM, which coincides with the subtour LP for the Traveling Salesman Problem on the metric closure of $\bar{G}$, gives a lower bound on the optimal cost. For instances where this LP is optimized by a half-integral solution $x$, Carr and Ravi (1998) showed that the integrality gap is at most $\frac43$: they show that the vector $\frac43 x$ dominates a convex combination of incidence vectors of $2$-edge connected spanning multisubgraphs of $\bar{G}$. We present a simpler proof of the result due to Carr and Ravi by applying an extension of Lov\'{a}sz's splitting-off theorem. Our proof naturally leads to a $\frac43$-approximation algorithm for half-integral instances. Given a half-integral solution $x$ to the LP for 2ECM, we give an $O(n^2)$-time algorithm to obtain a $2$-edge connected spanning multisubgraph of $\bar{G}$ whose cost is at most $\frac43 c^T x$.


title: Variations on Glauberman's ZJ Theorem

abstract: We give a new proof of Glauberman's ZJ Theorem, in a form that clarifies the choices involved and offers more choices than classical treatments. In particular, we introduce two new ZJ-type subgroups of a $p$-group $S$, that contain $ZJ_r(S)$ and $ZJ_o(S)$ respectively and are often strictly larger.