title: Anomalous atmospheric circulation favored the spread of COVID-19 in Europe

abstract: The current pandemic caused by the coronavirus SARS-CoV-2 is having negative health, social and economic consequences worldwide. In Europe, the pandemic started to develop strongly at the end of February and beginning of March 2020. It has subsequently spread over the continent, with special virulence in northern Italy and inland Spain. In this study we show that an unusual persistent anticyclonic situation prevailing in southwestern Europe during February 2020 (i.e. anomalously strong positive phase of the North Atlantic and Arctic Oscillations) could have resulted in favorable conditions, in terms of air temperature and humidity, in Italy and Spain for a quicker spread of the virus compared with the rest of the European countries. It seems plausible that the strong atmospheric stability and associated dry conditions that dominated in these regions may have favored the virus's propagation, by short-range droplet transmission as well as likely by long-range aerosol (airborne) transmission.


title: How Much Should I Pay? An Empirical Analysis on Monetary Prize in TopCoder

abstract: It is reported that task monetary prize is one of the most important motivating factors to attract crowd workers. While using expert-based methods to price Crowdsourcing tasks is a common practice, the challenge of validating the associated prices across different tasks is a constant issue. To address this issue, three different classifications of multiple linear regression, logistic regression, and K-nearest neighbor were compared to find the most accurate predicted price, using a dataset from the TopCoder website. The result of comparing chosen algorithms showed that the logistics regression model will provide the highest accuracy of 90% to predict the associated price to tasks and KNN ranked the second with an accuracy of 64% for K = 7. Also, applying PCA wouldn't lead to any better prediction accuracy as data components are not correlated.


title: Paretian social welfare relations and Baire property

abstract: We study the topological and set-theoretical nature of Paretian social welfare relations in a setting with infinite time horizon. Specifically, we answer questions posed in \citet{mathias2020} about the interplay between total welfare relations satisfying Pareto and anonymity principles with subsets of real numbers not satisfying the Baire property.


title: Assessing Discourse Relations in Language Generation from GPT-2

abstract: Recent advances in NLP have been attributed to the emergence of large-scale pre-trained language models. GPT-2, in particular, is suited for generation tasks given its left-to-right language modeling objective, yet the linguistic quality of its generated text has largely remain unexplored. Our work takes a step in understanding GPT-2's outputs in terms of discourse coherence. We perform a comprehensive study on the validity of explicit discourse relations in GPT-2's outputs under both organic generation and fine-tuned scenarios. Results show GPT-2 does not always generate text containing valid discourse relations; nevertheless, its text is more aligned with human expectation in the fine-tuned scenario. We propose a decoupled strategy to mitigate these problems and highlight the importance of explicitly modeling discourse information.


title: Simple curl-curl-conforming finite elements in two dimensions

abstract: We construct smooth finite element de Rham complexes in two space dimensions. This leads to three families of curl-curl conforming finite elements, two of which contain two existing families. The simplest triangular and rectangular finite elements have only 6 and 8 degrees of freedom, respectively. Numerical experiments for each family demonstrate the convergence and efficiency of the elements for solving the quad-curl problem.


title: Noisy Adaptive Group Testing using Bayesian Sequential Experimental Design

abstract: When the infection prevalence of a disease is low, Dorfman showed 80 years ago that testing groups of people can prove more efficient than testing people individually. Our goal in this paper is to propose new group testing algorithms that can operate in a noisy setting (tests can be mistaken) to decide adaptively (looking at past results) which groups to test next, with the goal to converge to a good detection, as quickly, and with as few tests as possible. We cast this problem as a Bayesian sequential experimental design problem. Using the posterior distribution of infection status vectors for $n$ patients, given observed tests carried out so far, we seek to form groups that have a maximal utility. We consider utilities such as mutual information, but also quantities that have a more direct relevance to testing, such as the AUC of the ROC curve of the test. Practically, the posterior distributions on $\{0,1\}^n$ are approximated by sequential Monte Carlo (SMC) samplers and the utility maximized by a greedy optimizer. Our procedures show in simulations significant improvements over both adaptive and non-adaptive baselines, and are far more efficient than individual tests when disease prevalence is low. Additionally, we show empirically that loopy belief propagation (LBP), widely regarded as the SoTA decoder to decide whether an individual is infected or not given previous tests, can be unreliable and exhibit oscillatory behavior. Our SMC decoder is more reliable, and can improve the performance of other group testing algorithms.


title: Pluto's Ultraviolet Spectrum, Surface Reflectance, and Airglow Emissions

abstract: During the New Horizons spacecraft's encounter with Pluto, the Alice ultraviolet spectrograph conducted a series of observations that detected emissions from both the interplanetary medium (IPM) and Pluto. In the direction of Pluto, the IPM was found to be 133.4$\pm$0.6R at Lyman $\alpha$, 0.24$\pm$0.02R at Lyman $\beta$, and <0.10R at He I 584{\AA}. We analyzed 3,900s of data obtained shortly before closest approach to Pluto and detect airglow emissions from H I, N I, N II, N$_2$, and CO above the disk of Pluto. We find Pluto's brightness at Lyman $\alpha$ to be $29.3\pm1.9$R, in good agreement with pre-encounter estimates. The detection of the N II multiplet at 1085{\AA} marks the first direct detection of ions in Pluto's atmosphere. We do not detect any emissions from noble gasses and place a 3$\sigma$ upper limit of 0.14 R on the brightness of the Ar I 1048{\AA} line. We compare pre-encounter model predictions and predictions from our own airglow model, based on atmospheric profiles derived from the solar occultation observed by New Horizons, to the observed brightness of Pluto's airglow. Although completely opaque at Lyman $\alpha$, Pluto's atmosphere is optically thin at wavelengths longer than 1425{\AA}. Consequently, a significant amount of solar FUV light reaches the surface, where it can participate in space weathering processes. From the brightness of sunlight reflected from Pluto, we find the surface has a reflectance factor (I/F) of 17% between 1400-1850{\AA}. We also report the first detection of an C$_3$ hydrocarbon molecule, methylacetylene, in absorption, at a column density of ~5$\times10^{15}$ cm$^{-2}$, corresponding to a column-integrated mixing ratio of $1.6\times10^{-6}$.


title: Self-consistent Nonlinear Force-free Field Reconstruction from Weighted Boundary Conditions

abstract: Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.


title: Hierarchical Low-Rank Approximation of Regularized Wasserstein Distance

abstract: Sinkhorn divergence is a measure of dissimilarity between two probability measures. It is obtained through adding an entropic regularization term to Kantorovich's optimal transport problem and can hence be viewed as an entropically regularized Wasserstein distance. Given two discrete probability vectors in the $n$-simplex and supported on two bounded spaces in ${\mathbb R}^d$, we present a fast method for computing Sinkhorn divergence when the cost matrix can be decomposed into a $d$-term sum of asymptotically smooth Kronecker product factors. The method combines Sinkhorn's matrix scaling iteration with a low-rank hierarchical representation of the scaling matrices to achieve a near-linear complexity ${\mathcal O}(n \log^3 n)$. This provides a fast and easy-to-implement algorithm for computing Sinkhorn divergence, enabling its applicability to large-scale optimization problems, where the computation of classical Wasserstein metric is not feasible. We present a numerical example related to signal processing to demonstrate the applicability of quadratic Sinkhorn divergence in comparison with quadratic Wasserstein distance and to verify the accuracy and efficiency of the proposed method.


title: Guided search for desired functional responses via Bayesian optimization of generative model: Hysteresis loop shape engineering in ferroelectrics

abstract: Advances in predictive modeling across multiple disciplines have yielded generative models capable of high veracity in predicting macroscopic functional responses of materials. Correspondingly, of interest is the inverse problem of finding the model parameter that will yield desired macroscopic responses, such as stress-strain curves, ferroelectric hysteresis loops, etc. Here we suggest and implement a Gaussian Process based methods that allow to effectively sample the degenerate parameter space of a complex non-local model to output regions of parameter space which yield desired functionalities. We discuss the specific adaptation of the acquisition function and sampling function to make the process efficient and balance the efficient exploration of parameter space for multiple possible minima and exploitation to densely sample the regions of interest where target behaviors are optimized. This approach is illustrated via the hysteresis loop engineering in ferroelectric materials, but can be adapted to other functionalities and generative models. The code is open-sourced and available at [github.com/ramav87/Ferrosim].


title: Exact solutions of the angular Teukolsky equation in particular cases

abstract: In this work, we propose a new scheme to solve the angular Teukolsky equation for the particular case: $m=0, s=0$. We first transform this equation to a confluent Heun differential equation and then construct the Wronskian determinant to calculate the eigenvalues and normalized eigenfunctions. We find that the eigenvalues for larger $l$ are approximately given by $_{0}{A_{l0}} \approx [l(l + 1) - \tau_{R}^2/2] - i\;\tau_{I}^2/2$ with an arbitrary $\tau^2=\tau_R^2 + i\,\tau_{I}^2$. The angular probability distribution (APD) for the ground state moves towards the north and south poles for $\tau_R^2>0$, but aggregates to the equator for $\tau_R^2\leq0$. However, we also notice that the APD for large angular momentum $l$ always moves towards the north and south poles , regardless the choice of $\tau^2$.


title: The Erd\H{o}s-R\'{e}nyi law of large numbers for ballistic random walk in random environment

abstract: We consider a one dimensional ballistic nearest-neighbor random walk in a random environment. We prove an Erd\H{o}s-R\'enyi strong law for the increments.


title: Sr$_2$MoO$_4$ and Sr$_2$RuO$_4$: Disentangling the Roles of Hund's and van Hove Physics

abstract: Sr$_2$MoO$_4$ is isostructural to the unconventional superconductor Sr$_2$RuO$_4$ but with two electrons instead of two holes in the Mo/Ru-t$_{2g}$ orbitals. Both materials are Hund's metals, but while Sr$_2$RuO$_4$ has a van Hove singularity in close proximity to the Fermi surface, the van Hove singularity of Sr$_2$MoO$_4$ is far from the Fermi surface. By using density functional plus dynamical mean-field theory we determine the relative influence of van Hove and Hund's metal physics on the correlation properties. We show that theoretically predicted signatures of Hund's metal physics occur on the occupied side of the electronic spectrum of Sr$_2$MoO$_4$, identifying Sr$_2$MoO$_4$ as an ideal candidate system for a direct experimental confirmation of the theoretical concept of Hund's metals via photoemission spectroscopy.


title: Persistent Radio Emission from Synchrotron Heating by a Repeating Fast Radio Burst Source in a Nebula

abstract: The first repeating fast radio burst (FRB), FRB 121102, was found to be associated with a spatially coincident, persistent nonthermal radio source, but the origin of the persistent emission remains unknown. In this paper, we propose that the persistent emission is produced via synchrotron-heating process by multiple bursts of FRB 121102 in a self-absorbed synchrotron nebula. As a population of bursts of the repeating FRB absorbed by the synchrotron nebula, the energy distribution of electrons in the nebula will change significantly. As a result, the spectrum of the nebula will show a hump steadily. For the persistent emission of FRB 121102, the total energy of bursts injecting into the nebula is required to be about $3.3\times10^{49}\,\unit{erg}$, the burst injection age is over $6.7\times 10^4\,\unit{yr}$, the nebula size is $\sim0.02\,\unit{pc}$, and the electron number is about $3.2\times10^{55}$. We predict that as more bursts inject, the brightness of the nebula would be brighter than the current observation, and meanwhile, the peak frequency would become higher. Due to the synchrotron absorption of the nebula, some low-frequency bursts would be absorbed, which may explain why most bursts were detected above $\sim1~\unit{GHz}$.


title: Relating CAT(0) cubical complexes and flag simplicial complexes

abstract: Given a finite CAT(0) cubical complex, we define a flag simplicial complex associated to it, called the crossing complex. We show that the crossing complex holds much of the combinatorial information of the original cubical complex: for example, hyperplanes in the cubical complex correspond to vertex links in the crossing complex, and the crossing complex is balanced if and only if the cubical complex is cubically balanced. The most significant result is that the sets of $f$-vectors of CAT(0) cubical complexes and flag simplicial complexes are equal, up to an invertible linear transformation.


title: Hamiltonicity in Cherry-quasirandom 3-graphs

abstract: We show that for any fixed $\alpha>0$, cherry-quasirandom 3-graphs of positive density and sufficiently large order $n$ with minimum vertex degree $\alpha \binom n2$ have a tight Hamilton cycle. This solves a conjecture of Aigner-Horev and Levy.


title: Transferable Perturbations of Deep Feature Distributions

abstract: Almost all current adversarial attacks of CNN classifiers rely on information derived from the output layer of the network. This work presents a new adversarial attack based on the modeling and exploitation of class-wise and layer-wise deep feature distributions. We achieve state-of-the-art targeted blackbox transfer-based attack results for undefended ImageNet models. Further, we place a priority on explainability and interpretability of the attacking process. Our methodology affords an analysis of how adversarial attacks change the intermediate feature distributions of CNNs, as well as a measure of layer-wise and class-wise feature distributional separability/entanglement. We also conceptualize a transition from task/data-specific to model-specific features within a CNN architecture that directly impacts the transferability of adversarial examples.


title: Inclusive diffractive $\rm \eta_{c}$ production in pp, pA and AA modes at the LHC

abstract: In this paper, the inclusive Pomeron-Pomeron, Reggeon-Reggeon, Pomeron-Reggeon as well as gluon-Pomeron(-Reggeon) and photon-Pomeron(-Reggeon) interactions for the $\rm \eta_{c}$ at the LHC energies have been examined in proton-proton, proton-nucleus and nucleus-nucleus collision modes. The cross section has been computed based on the NRQCD factorization and Regge theory formalism. The cross exchange of Pomeron-Reggeon contribution is important in $\rm pp$ and $\rm pA$ modes. The Pomeron-Pomeron contribution is significant in $\rm AA$ mode. The Pomeron contribution is considerable for $\rm AA$ and $\rm pA$ modes in single diffractive process where $\rm A$ undergoes the diffractive process. Reggeon contribution is sizable in $\rm pp$ and $\rm pA$ modes where $\rm p$ only undergoes the diffractive process. The Pomeron and Reggeon contributions in photon-Pomeron and-Reggeon process remain smaller than that of gluon-Pomeron and-Reggeon processes. Our results show that the experimental study of Reggeon, Pomeron and their cross exchange can be carried out in certain kinematic windows with the specific choice of the mode at the LHC. The investigation can be useful to better constrain the Reggeon and Pomeron parton content. The inclusive process serves as the background to related exclusive processes which should be predicted.


title: Convex hulls of polynomial Julia sets

abstract: We prove P. Alexandersson's conjecture that for every complex polynomial $p$ of degree $d \geq 2$ the convex hull $H_p$ of the Julia set $J_p$ of $p$ satisfies $p^{-1}(H_p) \subset H_p$. We further prove that the equality $p^{-1}(H_p) = H_p$ is achieved only if $p$ is affinely conjugated to the Chebyshev polynomial $T_d$ of degree $d$, to $-T_d$ or a monomial $c z^d$ with $|c|=1$.


title: Foliated corona decompositions

abstract: We prove that the $L_4$ norm of the vertical perimeter of any measurable subset of the $3$-dimensional Heisenberg group $\mathbb{H}$ is at most a universal constant multiple of the (Heisenberg) perimeter of the subset. We show that this isoperimetric-type inequality is optimal in the sense that there are sets for which it fails to hold with the $L_4$ norm replaced by the $L_q$ norm for any $q<4$. This is in contrast to the $5$-dimensional setting, where the above result holds with the $L_4$ norm replaced by the $L_2$ norm. The proof of the aforementioned isoperimetric inequality introduces a new structural methodology for understanding the geometry of surfaces in $\mathbb{H}$. In previous work (2017) we showed how to obtain a hierarchical decomposition of Ahlfors-regular surfaces into pieces that are approximately intrinsic Lipschitz graphs. Here we prove that any such graph admits a foliated corona decomposition, which is a family of nested partitions into pieces that are close to ruled surfaces. Apart from the intrinsic geometric and analytic significance of these results, which settle questions posed by Cheeger-Kleiner-Naor (2009) and Lafforgue-Naor (2012), they have several noteworthy implications, including the fact that the $L_1$ distortion of a word-ball of radius $n\ge 2$ in the discrete $3$-dimensional Heisenberg group is bounded above and below by universal constant multiples of $\sqrt[4]{\log n}$; this is in contrast to higher dimensional Heisenberg groups, where our previous work showed that the distortion of a word-ball of radius $n\ge 2$ is of order $\sqrt{\log n}$.


title: About apparent superluminal drives in generic gravity theories

abstract: As is well known, there exists warp drives in GR, such as the Alcubierre bubbles, which achieve an apparent faster than light travel \cite{alcubierre}. A result due to Gao and Wald \cite{gaowald} suggests that such a travel is unlikely for GR with matter satisfying both the Null Energy and the Null Generic Conditions. There exists a generalization of this statement due to Galloway, that ensures that the Gao-Wald result is true regardless the underlying gravity model, unless there exists at least one inextendible null geodesic with achronal image in the space time (a null line). The proof of this proposition is based on techniques of causal theories, and has never been released. In the present work an independent proof of this result is presented by use of the Raychaudhuri equation, and avoiding several technical complications described along the text. Some consequences of these affirmations are discussed at last section, in particular their potential use in problems of causality.


title: Sequential Interpretability: Methods, Applications, and Future Direction for Understanding Deep Learning Models in the Context of Sequential Data

abstract: Deep learning continues to revolutionize an ever-growing number of critical application areas including healthcare, transportation, finance, and basic sciences. Despite their increased predictive power, model transparency and human explainability remain a significant challenge due to the "black box" nature of modern deep learning models. In many cases the desired balance between interpretability and performance is predominately task specific. Human-centric domains such as healthcare necessitate a renewed focus on understanding how and why these frameworks are arriving at critical and potentially life-or-death decisions. Given the quantity of research and empirical successes of deep learning for computer vision, most of the existing interpretability research has focused on image processing techniques. Comparatively, less attention has been paid to interpreting deep learning frameworks using sequential data. Given recent deep learning advancements in highly sequential domains such as natural language processing and physiological signal processing, the need for deep sequential explanations is at an all-time high. In this paper, we review current techniques for interpreting deep learning techniques involving sequential data, identify similarities to non-sequential methods, and discuss current limitations and future avenues of sequential interpretability research.


title: Fully Embedding Fast Convolutional Networks on Pixel Processor Arrays

abstract: We present a novel method of CNN inference for pixel processor array (PPA) vision sensors, designed to take advantage of their massive parallelism and analog compute capabilities. PPA sensors consist of an array of processing elements (PEs), with each PE capable of light capture, data storage and computation, allowing various computer vision processing to be executed directly upon the sensor device. The key idea behind our approach is storing network weights "in-pixel" within the PEs of the PPA sensor itself to allow various computations, such as multiple different image convolutions, to be carried out in parallel. Our approach can perform convolutional layers, max pooling, ReLu, and a final fully connected layer entirely upon the PPA sensor, while leaving no untapped computational resources. This is in contrast to previous works that only use a sensor-level processing to sequentially compute image convolutions, and must transfer data to an external digital processor to complete the computation. We demonstrate our approach on the SCAMP-5 vision system, performing inference of a MNIST digit classification network at over 3000 frames per second and over 93% classification accuracy. This is the first work demonstrating CNN inference conducted entirely upon the processor array of a PPA vision sensor device, requiring no external processing.


title: Subgoal Planning Algorithm for Autonomous Vehicle Guidance

abstract: Trained humans exhibit highly agile spatial skills, enabling them to operate vehicles with complex dynamics in demanding tasks and conditions. Prior work shows that humans achieve this performance by using strategies such as satisficing, learning hierarchical task structure, and using a library of motion primitive elements. A key aspect of efficient and versatile solutions is the decomposition of a task into a sequence of smaller tasks, represented by subgoals. The present work uses properties of constrained optimal control theory to define conditions that specify candidate subgoal states and enable this decomposition. The proposed subgoal algorithm uses graph search to determine a subgoal sequence that links a series of unconstrained motion guidance elements into a constrained solution trajectory. In simulation experiments, the subgoal guidance algorithm generates paths with higher performance and less computation time than an RRT* benchmark. Examples illustrate the robustness and versatility of this approach in multiple environment types.


title: Neural Machine Translation with Monte-Carlo Tree Search

abstract: Recent algorithms in machine translation have included a value network to assist the policy network when deciding which word to output at each step of the translation. The addition of a value network helps the algorithm perform better on evaluation metrics like the BLEU score. After training the policy and value networks in a supervised setting, the policy and value networks can be jointly improved through common actor-critic methods. The main idea of our project is to instead leverage Monte-Carlo Tree Search (MCTS) to search for good output words with guidance from a combined policy and value network architecture in a similar fashion as AlphaZero. This network serves both as a local and a global look-ahead reference that uses the result of the search to improve itself. Experiments using the IWLST14 German to English translation dataset show that our method outperforms the actor-critic methods used in recent machine translation papers.


title: The fourth moment of central values of quadratic Hecke $L$-functions in the Gaussian field

abstract: We obtain an asymptotic formula for the fourth moment of central values of a family of quadratic Hecke $L$-functions in the Gaussian field under the generalized Riemann hypothesis (GRH). We also establish lower bounds unconditionally and upper bounds under GRH for higher moments of the same family.


title: Towards Accurate and Robust Domain Adaptation under Noisy Environments

abstract: In non-stationary environments, learning machines usually confront the domain adaptation scenario where the data distribution does change over time. Previous domain adaptation works have achieved great success in theory and practice. However, they always lose robustness in noisy environments where the labels and features of examples from the source domain become corrupted. In this paper, we report our attempt towards achieving accurate noise-robust domain adaptation. We first give a theoretical analysis that reveals how harmful noises influence unsupervised domain adaptation. To eliminate the effect of label noise, we propose an offline curriculum learning for minimizing a newly-defined empirical source risk. To reduce the impact of feature noise, we propose a proxy distribution based margin discrepancy. We seamlessly transform our methods into an adversarial network that performs efficient joint optimization for them, successfully mitigating the negative influence from both data corruption and distribution shift. A series of empirical studies show that our algorithm remarkably outperforms state of the art, over 10% accuracy improvements in some domain adaptation tasks under noisy environments.


title: Analysis of the Stochastic Alternating Least Squares Method for the Decomposition of Random Tensors

abstract: Stochastic Alternating Least Squares (SALS) is a method that approximates the canonical decomposition of averages of sampled random tensors. Its simplicity and efficient memory usage make SALS an ideal tool for decomposing tensors in an online setting. We show, under mild regularization and readily verifiable assumptions on the boundedness of the data, that the SALS algorithm is globally convergent. Numerical experiments validate our theoretical findings and demonstrate the algorithm's performance and complexity.


title: Spatial-Temporal Mitosis Detection in Phase-Contrast Microscopy via Likelihood Map Estimation by 3DCNN

abstract: Automated mitotic detection in time-lapse phasecontrast microscopy provides us much information for cell behavior analysis, and thus several mitosis detection methods have been proposed. However, these methods still have two problems; 1) they cannot detect multiple mitosis events when there are closely placed. 2) they do not consider the annotation gaps, which may occur since the appearances of mitosis cells are very similar before and after the annotated frame. In this paper, we propose a novel mitosis detection method that can detect multiple mitosis events in a candidate sequence and mitigate the human annotation gap via estimating a spatiotemporal likelihood map by 3DCNN. In this training, the loss gradually decreases with the gap size between ground truth and estimation. This mitigates the annotation gaps. Our method outperformed the compared methods in terms of F1- score using a challenging dataset that contains the data under four different conditions.


title: In-Place Parallel-Partition Algorithms using Exclusive-Read-and-Write Memory

abstract: We present an in-place algorithm for the partition problem that has linear work and polylogarithmic span. The algorithm uses only exclusive read/write shared variables, and can be implemented using parallel-for-loops without any additional concurrency considerations (i.e., the algorithm is EREW). A key feature of the algorithm is that it exhibits provably optimal cache behavior, up to small-order factors. We also present a second in-place EREW algorithm for the partition problem that has linear work and span $O(\log n \cdot \log \log n)$, which is within an $O(\log\log n)$ factor of the optimal span. By using this low-span algorithm as a subroutine within the cache-friendly algorithm, we are able to obtain a single EREW algorithm that combines their theoretical guarantees: the algorithm achieves span $O(\log n \cdot \log \log n)$ and optimal cache behavior. As an immediate consequence, we also get an in-place EREW quicksort algorithm with work $O(n \log n)$, span $O(\log^2 n \cdot \log \log n)$. Whereas the standard EREW algorithm for parallel partitioning is memory-bandwidth bound on large numbers of cores, our cache-friendly algorithm is able to achieve near-ideal scaling in practice by avoiding the memory-bandwidth bottleneck. The algorithm's performance is comparable to that of the Blocked Strided Algorithm of Francis, Pannan, Frias, and Petit, which is the previous state-of-the art for parallel EREW sorting algorithms, but which lacks theoretical guarantees on its span and cache behavior.


title: Interfacial thermodynamics of spherical nanodroplets: Molecular understanding of surface tension via hydrogen bond network

abstract: Surface tension plays a ubiquitous role in phase transitions including condensation or evaporation of atmospheric liquid droplets. Especially, understanding of interfacial thermodynamics of the critical nucleus of 1 nm scale is important for molecular characterization of the activation energy barrier of nucleation. Here, we investigate surface tension of spherical nanodroplets with both molecular dynamics and density functional theory, and find that surface tension decreases appreciably below 1 nm radius, whose analytic expression is consistently derived from the classic Tolman's equation. In particular, the free energy analysis of nanodroplets shows that the change of surface tension originates dominantly from the configurational energy of interfacial molecules, which is evidenced by the increasingly disrupted hydrogen bond network as the droplet size decreases. Our result can be applied to the interface-related phenomena associated with molecular fluctuations such as biomolecule adsorption at sub-nm scale where the macroscopic thermodynamic quantities are ill-defined.


title: Tsallis meets Boltzmann: q-index for a finite ideal gas and its thermodynamic limit

abstract: Nonadditive Tsallis $q$-statistics has successfully been applied for a plethora of systems in natural sciences and other branches of knowledge. Nevertheless, its foundations have been severely criticised by some authors based on the standard additive Boltzmann-Gibbs approach thereby remaining a quite controversial subject. In order to clarify some polemical concepts, the distribution function for an ideal gas with a finite number of point particles and its $q$-index are analytically determined. The two-particle correlation function is also derived. The degree of correlation diminishes continuously with the growth of the number of particles. The ideal finite gas system is usually correlated, becomes less correlated when the number of particles grows, and is finally, fully uncorrelated when the molecular chaos regime is reached. It is also advocated that both approaches can be confronted through a careful kinetic spectroscopic experiment. The analytical results derived here suggest that Tsallis q-statistics may play a physical role more fundamental than usually discussed in the literature.


title: Difficulty Translation in Histopathology Images

abstract: The unique nature of histopathology images opens the door to domain-specific formulations of image translation models. We propose a difficulty translation model that modifies colorectal histopathology images to be more challenging to classify. Our model comprises a scorer, which provides an output confidence to measure the difficulty of images, and an image translator, which learns to translate images from easy-to-classify to hard-to-classify using a training set defined by the scorer. We present three findings. First, generated images were indeed harder to classify for both human pathologists and machine learning classifiers than their corresponding source images. Second, image classifiers trained with generated images as augmented data performed better on both easy and hard images from an independent test set. Finally, human annotator agreement and our model's measure of difficulty correlated strongly, implying that for future work requiring human annotator agreement, the confidence score of a machine learning classifier could be used as a proxy.


title: Correlation-induced insulating topological phases at charge neutrality in twisted bilayer graphene

abstract: Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlations is fundamental to shed light on the nature of the superconducting and correlated insulating states observed in the TBG experiments. Here, we use unbiased, sign-problem-free quantum Monte Carlo simulations to solve an effective interacting lattice model for TBG at charge neutrality. Besides the usual cluster Hubbard-like repulsion, this model also contains an assisted hopping interaction that emerges due to the non-trivial topological properties of TBG. Such a non-local interaction fundamentally alters the phase diagram at charge neutrality, gapping the Dirac cones even for infinitesimally small interaction. As the interaction strength increases, a sequence of different correlated insulating phases emerge, including a quantum valley Hall state with topological edge states, an intervalley-coherent insulator, and a valence bond solid. The charge-neutrality correlated insulating phases discovered here provide the sought-after reference states needed for a comprehensive understanding of the insulating states at integer fillings and the proximate superconducting states of TBG.


title: Towards Data-Efficient Learning: A Benchmark for COVID-19 CT Lung and Infection Segmentation

abstract: Purpose: Accurate segmentation of lung and infection in COVID-19 CT scans plays an important role in the quantitative management of patients. Most of the existing studies are based on large and private annotated datasets that are impractical to obtain from a single institution, especially when radiologists are busy fighting the coronavirus disease. Furthermore, it is hard to compare current COVID-19 CT segmentation methods as they are developed on different datasets, trained in different settings, and evaluated with different metrics. Methods: To promote the development of data-efficient deep learning methods, in this paper, we built three benchmarks for lung and infection segmentation based on 70 annotated COVID-19 cases, which contain current active research areas, e.g., few-shot learning, domain generalization, and knowledge transfer. For a fair comparison among different segmentation methods, we also provide standard training, validation and testing splits, evaluation metrics and, the corresponding code. Results: Based on the state-of-the-art network, we provide more than 40 pre-trained baseline models, which not only serve as out-of-the-box segmentation tools but also save computational time for researchers who are interested in COVID-19 lung and infection segmentation. We achieve average Dice Similarity Coefficient (DSC) scores of 97.3\%, 97.7\%, and 67.3\% and average Normalized Surface Dice (NSD) scores of 90.6\%, 91.4\%, and 70.0\% for left lung, right lung, and infection, respectively. Conclusions: To the best of our knowledge, this work presents the first data-efficient learning benchmark for medical image segmentation and the largest number of pre-trained models up to now. All these resources are publicly available, and our work lays the foundation for promoting the development of deep learning methods for efficient COVID-19 CT segmentation with limited data.


title: Interpretation of Deep Temporal Representations by Selective Visualization of Internally Activated Nodes

abstract: Recently deep neural networks demonstrate competitive performances in classification and regression tasks for many temporal or sequential data. However, it is still hard to understand the classification mechanisms of temporal deep neural networks. In this paper, we propose two new frameworks to visualize temporal representations learned from deep neural networks. Given input data and output, our algorithm interprets the decision of temporal neural network by extracting highly activated periods and visualizes a sub-sequence of input data which contributes to activate the units. Furthermore, we characterize such sub-sequences with clustering and calculate the uncertainty of the suggested type and actual data. We also suggest Layer-wise Relevance from the output of a unit, not from the final output, with backward Monte-Carlo dropout to show the relevance scores of each input point to activate units with providing a visual representation of the uncertainty about this impact.


title: Intelligent Reflecting Surface Assisted Anti-Jamming Communications: A Fast Reinforcement Learning Approach

abstract: Malicious jamming launched by smart jammers can attack legitimate transmissions, which has been regarded as one of the critical security challenges in wireless communications. With this focus, this paper considers the use of an intelligent reflecting surface (IRS) to enhance anti-jamming communication performance and mitigate jamming interference by adjusting the surface reflecting elements at the IRS. Aiming to enhance the communication performance against a smart jammer, an optimization problem for jointly optimizing power allocation at the base station (BS), and reflecting beamforming at the IRS is formulated while considering quality of service (QoS) requirements of legitimate users. As the jamming model and jamming behavior are dynamic and unknown, a fuzzy win or learn fast-policy hill-climbing (WoLFPHC) learning approach is proposed to jointly optimize the anti-jamming power allocation and reflecting beamforming strategy, where WoLFPHC is capable of quickly achieving the optimal policy without the knowledge of the jamming model, and fuzzy state aggregation can represent the uncertain environment states as aggregate states. Simulation results demonstrate that the proposed anti-jamming learning-based approach can efficiently improve both the IRS-assisted system rate and transmission protection level compared with existing solutions.


title: Data-Driven Construction of Data Center Graph of Things for Anomaly Detection

abstract: Data center (DC) contains both IT devices and facility equipment, and the operation of a DC requires a high-quality monitoring (anomaly detection) system. There are lots of sensors in computer rooms for the DC monitoring system, and they are inherently related. This work proposes a data-driven pipeline (ts2graph) to build a DC graph of things (sensor graph) from the time series measurements of sensors. The sensor graph is an undirected weighted property graph, where sensors are the nodes, sensor features are the node properties, and sensor connections are the edges. The sensor node property is defined by features that characterize the sensor events (behaviors), instead of the original time series. The sensor connection (edge weight) is defined by the probability of concurrent events between two sensors. A graph of things prototype is constructed from the sensor time series of a real data center, and it successfully reveals meaningful relationships between the sensors. To demonstrate the use of the DC sensor graph for anomaly detection, we compare the performance of graph neural network (GNN) and existing standard methods on synthetic anomaly data. GNN outperforms existing algorithms by a factor of 2 to 3 (in terms of precision and F1 score), because it takes into account the topology relationship between DC sensors. We expect that the DC sensor graph can serve as the infrastructure for the DC monitoring system since it represents the sensor relationships.


title: Forecast analysis of the epidemics trend of COVID-19 in the United States by a generalized fractional-order SEIR model

abstract: In this paper, a generalized fractional-order SEIR model is proposed, denoted by SEIQRP model, which has a basic guiding significance for the prediction of the possible outbreak of infectious diseases like COVID-19 and other insect diseases in the future. Firstly, some qualitative properties of the model are analyzed. The basic reproduction number $R_{0}$ is derived. When $R_{0}<1$, the disease-free equilibrium point is unique and locally asymptotically stable. When $R_{0}>1$, the endemic equilibrium point is also unique. Furthermore, some conditions are established to ensure the local asymptotic stability of disease-free and endemic equilibrium points. The trend of COVID-19 spread in the United States is predicted. Considering the influence of the individual behavior and government mitigation measurement, a modified SEIQRP model is proposed, defined as SEIQRPD model. According to the real data of the United States, it is found that our improved model has a better prediction ability for the epidemic trend in the next two weeks. Hence, the epidemic trend of the United States in the next two weeks is investigated, and the peak of isolated cases are predicted. The modified SEIQRP model successfully capture the development process of COVID-19, which provides an important reference for understanding the trend of the outbreak.


title: Isolated Photonic Flatband with the Effective Magnetic Flux in A Synthetic Space including the Frequency Dimension

abstract: The exploration of flatband in photonics is fundamentally important, aiming to control the localization of light for potential applications in optical communications. We study the flatband physics in a synthetic space including the frequency axis of light. A ring-resonator array is used to construct a synthetic Lieb-type lattice, where the modulation phase distribution supports a locally non-zero effective magnetic flux pointing into the synthetic space. We find that the flatband is isolated from other dispersive bands and the light can still be localized even with the perturbation from the group velocity dispersion of the waveguide. Our work points out a route towards manipulating the localization performance of a wavepacket of light along both the spatial dimension and the frequency dimension, which holds the potential implication for controlling the storage of optical information in dispersive materials.


title: Wide viewing-angle holographic display based on enhanced-NA Fresnel hologram

abstract: The viewing-angle enlargement of a holographic image is a crucial factor for realizing the holographic display. The numerical aperture (NA) of digital hologram other than a pixel specification has been known to determine the angular field extent of image. Here, we provide a valid foundation for the dependence of viewing angle on the hologram numerical aperture by investigating mathematically the internal structure of the sampled point spread function showing a self-similarity of its modulating curves. The enhanced-NA Fresnel hologram reconstructs the images at a viewing angle larger than a diffraction angle by a hologram pixel pitch where its angle value is expressed in terms of the NA of whole hologram aperture, which is systematically observed by optical hologram imaging. Finally, we found that the aliased replica noises generated in the enhanced-NA Fresnel diffraction regime are effectively suppressed within the diffraction scope by a digitized pixel. This characteristic enables us to overcome the image reduction and to remove the interference of high-order images, which leads to the wide viewing-angle holographic display. Optical experiments are shown to be consistent with the results of numerical simulation.


title: Optically detected magnetic resonance in neutral silicon vacancy centers in diamond via bound exciton states

abstract: Neutral silicon vacancy (SiV0) centers in diamond are promising candidates for quantum networks because of their excellent optical properties and long spin coherence times. However, spin-dependent fluorescence in such defects has been elusive due to poor understanding of the excited state fine structure and limited off-resonant spin polarization. Here we report the realization of optically detected magnetic resonance and coherent control of SiV0 centers at cryogenic temperatures, enabled by efficient optical spin polarization via previously unreported higher-lying excited states. We assign these states as bound exciton states using group theory and density functional theory. These bound exciton states enable new control schemes for SiV0 as well as other emerging defect systems.


title: Wireless VR/Haptic Open Platform for Multimodal Teleoperation

abstract: With emerging trends in the fifth generation and robotics, the Internet of Skills will enable us to deliver skills or expertise anywhere over the Internet. In this paper, we propose a wireless connected virtual reality and haptic communication open platform to show the proof of concept for multimodal teleoperation systems in real-time. We focus on a practical implementation with commercial products to facilitate the access and modification of the system. The performance of the system is measured in terms of system latency and user-centric metrics.


title: Demo: A Reinforcement Learning-based Flexible Duplex System for B5G with Sub-6 GHz

abstract: In this paper, we propose a reinforcement learning-based flexible duplex system for B5G with Sub-6 GHz. This system combines full-duplex radios and dynamic spectrum access to maximize the spectral efficiency. We verify this method's feasibility by implementing an FPGA-based real-time testbed. In addition, we compare the proposed algorithm with the result derived from the numerical analysis through system-level evaluations.


title: On the Formation of SS433

abstract: SS433 is an extraordinary X-ray binary which is ejecting bipolar jets with $26\%$ of the speed of light. Associated with the supernova-like shell W50, SS433 is also probably one of the youngest X-ray binaries with an age of $\lesssim 10^5$ yr. However, the masses of the two components in SS433 and even the nature of the compact object are still under debate. In this work, assuming that the compact object is a black hole, we employ a binary population synthesis method to study the formation of SS433. We use previous estimates of the age of W50 and the duration of the jet activity to constrain the evolutionary history. Our calculations suggest that SS433 likely harboured a Hertzsprung gap star at the beginning of the current Roche-lobe overflow phase. The masses of the black hole and the optical/donor star in SS433 suggested by the simulations are around $8\ M_{\odot}$ and $24\ M_{\odot}$, respectively. Future measurement of the donor mass and radius can help infer not only the orgin of the binary but also the nature of the nebula W50.


title: ALMA ACA and Nobeyama observations of two Orion cores in deuterated molecular lines

abstract: We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.


title: A theory of type-II minimally modified gravity

abstract: We propose a modified gravity theory that propagates only two local gravitational degrees of freedom and that does not have an Einstein frame. According to the classification in JCAP 01 (2019) 017 [arXiv:1810.01047 [gr-qc]], this is a type-II minimally modified gravity theory. The theory is characterized by the gravitational constant $G_{\rm N}$ and a function $V(\phi)$ of a non-dynamical auxiliary field $\phi$ that plays the role of dark energy. Once one fixes a homogeneous and isotropic cosmological background, the form of $V(\phi)$ is determined and the theory no longer possesses a free parameter or a free function, besides $G_{\rm N}$. For $V'(\phi) = 0$ the theory reduces to general relativity (GR) with $G_N$ being the Newton's constant and $V=const.$ being the cosmological constant. For $V'(\phi) \ne 0$, it is shown that gravity behaves differently from GR but that GR with $G_{\rm N}$ being the Newton's constant is recovered for weak gravity at distance and time scales sufficiently shorter than the scale associated with $V(\phi)$. Therefore this theory provides the simplest framework of cosmology in which deviations from GR can be tested by observational data.


title: Hamiltonian Monte Carlo using an adjoint-differentiated Laplace approximation: Bayesian inference for latent Gaussian models and beyond

abstract: Gaussian latent variable models are a key class of Bayesian hierarchical models with applications in many fields. Performing Bayesian inference on such models can be challenging as Markov chain Monte Carlo algorithms struggle with the geometry of the resulting posterior distribution and can be prohibitively slow. An alternative is to use a Laplace approximation to marginalize out the latent Gaussian variables and then integrate out the remaining hyperparameters using dynamic Hamiltonian Monte Carlo, a gradient-based Markov chain Monte Carlo sampler. To implement this scheme efficiently, we derive a novel adjoint method that propagates the minimal information needed to construct the gradient of the approximate marginal likelihood. This strategy yields a scalable differentiation method that is orders of magnitude faster than state of the art differentiation techniques when the hyperparameters are high dimensional. We prototype the method in the probabilistic programming framework Stan and test the utility of the embedded Laplace approximation on several models, including one where the dimension of the hyperparameter is $\sim$6,000. Depending on the cases, the benefits can include an alleviation of the geometric pathologies that frustrate Hamiltonian Monte Carlo and a dramatic speed-up.


title: Dynamic Predictions of Postoperative Complications from Explainable, Uncertainty-Aware, and Multi-Task Deep Neural Networks

abstract: Accurate prediction of postoperative complications can inform shared decisions regarding prognosis, preoperative risk-reduction, and postoperative resource use. We hypothesized that multi-task deep learning models would outperform random forest models in predicting postoperative complications, and that integrating high-resolution intraoperative physiological time series would result in more granular and personalized health representations that would improve prognostication compared to preoperative predictions. In a longitudinal cohort study of 56,242 patients undergoing 67,481 inpatient surgical procedures at a university medical center, we compared deep learning models with random forests for predicting nine common postoperative complications using preoperative, intraoperative, and perioperative patient data. Our study indicated several significant results across experimental settings that suggest the utility of deep learning for capturing more precise representations of patient health for augmented surgical decision support. Multi-task learning improved efficiency by reducing computational resources without compromising predictive performance. Integrated gradients interpretability mechanisms identified potentially modifiable risk factors for each complication. Monte Carlo dropout methods provided a quantitative measure of prediction uncertainty that has the potential to enhance clinical trust. Multi-task learning, interpretability mechanisms, and uncertainty metrics demonstrated potential to facilitate effective clinical implementation.


title: Finding compositional inverses of permutations from the AGW criterion

abstract: Permutation polynomials and their compositional inverses have wide applications in cryptography, coding theory, and combinatorial designs. Motivated by several previous results on finding compositional inverses of permutation polynomials of different forms, we propose a general method for finding these inverses of permutation polynomials constructed by the AGW criterion. As a result, we have reduced the problem of finding the compositional inverse of such a permutation polynomial over a finite field to that of finding the inverse of a bijection over a smaller set. We demonstrate our method by interpreting several recent known results, as well as by providing new explicit results on more classes of permutation polynomials in different types. In addition, we give new criteria for these permutation polynomials being involutions. Explicit constructions are also provided for all involutory criteria.


title: Differentiating through Log-Log Convex Programs

abstract: We show how to efficiently compute the derivative (when it exists) of the solution map of log-log convex programs (LLCPs). These are nonconvex, nonsmooth optimization problems with positive variables that become convex when the variables, objective functions, and constraint functions are replaced with their logs. We focus specifically on LLCPs generated by disciplined geometric programming, a grammar consisting of a set of atomic functions with known log-log curvature and a composition rule for combining them. We represent a parametrized LLCP as the composition of a smooth transformation of parameters, a convex optimization problem, and an exponential transformation of the convex optimization problem's solution. The derivative of this composition can be computed efficiently, using recently developed methods for differentiating through convex optimization problems. We implement our method in CVXPY, a Python-embedded modeling language and rewriting system for convex optimization. In just a few lines of code, a user can specify a parametrized LLCP, solve it, and evaluate the derivative or its adjoint at a vector. This makes it possible to conduct sensitivity analyses of solutions, given perturbations to the parameters, and to compute the gradient of a function of the solution with respect to the parameters. We use the adjoint of the derivative to implement differentiable log-log convex optimization layers in PyTorch and TensorFlow. Finally, we present applications to designing queuing systems and fitting structured prediction models.


title: Forecasting in Non-stationary Environments with Fuzzy Time Series

abstract: In this paper we introduce a Non-Stationary Fuzzy Time Series (NSFTS) method with time varying parameters adapted from the distribution of the data. In this approach, we employ Non-Stationary Fuzzy Sets, in which perturbation functions are used to adapt the membership function parameters in the knowledge base in response to statistical changes in the time series. The proposed method is capable of dynamically adapting its fuzzy sets to reflect the changes in the stochastic process based on the residual errors, without the need to retraining the model. This method can handle non-stationary and heteroskedastic data as well as scenarios with concept-drift. The proposed approach allows the model to be trained only once and remain useful long after while keeping reasonable accuracy. The flexibility of the method by means of computational experiments was tested with eight synthetic non-stationary time series data with several kinds of concept drifts, four real market indices (Dow Jones, NASDAQ, SP500 and TAIEX), three real FOREX pairs (EUR-USD, EUR-GBP, GBP-USD), and two real cryptocoins exchange rates (Bitcoin-USD and Ethereum-USD). As competitor models the Time Variant fuzzy time series and the Incremental Ensemble were used, these are two of the major approaches for handling non-stationary data sets. Non-parametric tests are employed to check the significance of the results. The proposed method shows resilience to concept drift, by adapting parameters of the model, while preserving the symbolic structure of the knowledge base.


title: Requirements and Technologies Towards UAM: Communication, Navigation, and Surveillance

abstract: Urban air mobility (UAM) is a concept for creating an airborne transportation system that operates in urban settings with an on-board pilot and/or remote pilot in command (RPIC), or with a fully autonomous architecture. Although the passenger traffic will be mostly in and near urban environments, UAM is also being considered for air cargo, perhaps between cities. Such capability is pushing the current communication, navigation and surveillance (CNS) / air traffic management (ATM) systems that were not designed to support these types of aviation scenarios. The UAM aircraft will be operating in a congested environment, where CNS and ATM systems need to provide integrity, robustness, security, and very high availability for safety of UAM operations while evolving. As UAM is under research by academia and government agencies, the industry is driving technology towards aircraft prototypes. Critical UAM requirements are derived from command and control (C2) (particularly for RPIC scenario), data connectivity for passengers and flight systems, unmanned aircraft systems (UAS) to UAS communication to avoid collision, and data exchange for positioning and surveillance. In this paper, we study connectivity challenges and present requirements towards a robust UAM architecture through its concept of operations. In addition, we review the existing/potential CNS technologies towards UAM, i.e., 3rd generation partnership project (3GPP) fifth generation (5G) new radio (NR), navigation detect & avoid (DAA), and satellite systems and present conclusions on a future road-map for UAM CNS architecture.


title: Phase diagram reconstruction of the Bose-Hubbard model with a Restricted Boltzmann Machine wavefunction

abstract: Recently, the use of neural quantum states for describing the ground state of many- and few-body problems has been gaining popularity because of their high expressivity and ability to handle intractably large Hilbert spaces. In particular, methods based on variational Monte Carlo have proven to be successful in describing the physics of bosonic systems such as the Bose-Hubbard (BH) model. However, this technique has not been systematically tested on the parameter space of the BH model, particularly at the boundary between the Mott insulator and superfluid phases. In this work, we evaluate the capabilities of variational Monte Carlo with a trial wavefunction given by a Restricted Boltzmann Machine to reproduce the quantum ground state of the BH model on several points of its parameter space. To benchmark the technique, we compare its results to the ground state found through exact diagonalization for small one-dimensional chains. In general, we find that the learned ground state correctly estimates many observables, reproducing to a high degree the phase diagram for the first Mott lobe and part of the second one. However, we find that the technique is challenged whenever the system transitions between excitation manifolds, as the ground state is not learned correctly at these boundaries. We improve the quality of the results produced by the technique by proposing a method to discard noisy probabilities learned in the ground state.


title: High Energy Modular Bootstrap, Global Symmetries and Defects

abstract: We derive Cardy-like formulas for the growth of operators in different sectors of unitary $2$ dimensional CFT in the presence of topological defect lines by putting an upper and lower bound on the number of states with scaling dimension in the interval $[\Delta-\delta,\Delta+\delta]$ for large $\Delta$ at fixed $\delta$. Consequently we prove that given any unitary modular invariant $2$D CFT symmetric under finite global symmetry $G$ (acting faithfully), all the irreducible representations of $G$ appear in the spectra of the untwisted sector; the growth of states is Cardy like and proportional to the ''square'' of the dimension of the irrep. In the Schwarzian limit, the result matches onto that of JT gravity with a bulk gauge theory. If the symmetry is non-anomalous, the result applies to any sector twisted by a group element. For $c>1$, the statements are true for Virasoro primaries. Furthermore, the results are applicable to large c CFTs. We also extend our results for the continuous $U(1)$ group.


title: Rainbow matchings for 3-uniform hypergraphs

abstract: K\"{u}hn, Osthus, and Treglown and, independently, Khan proved that if $H$ is a $3$-uniform hypergraph with $n$ vertices such that $n\in 3\mathbb{Z}$ and large, and $\delta_1(H)>{n-1\choose 2}-{2n/3\choose 2}$, then $H$ contains a perfect matching. In this paper, we show that for $n\in 3\mathbb{Z}$ sufficiently large, if $F_1, \ldots, F_{n/3}$ are 3-uniform hypergrapghs with a common vertex set and $\delta_1(F_i)>{n-1\choose 2}-{2n/3\choose 2}$ for $i\in [n/3]$, then $\{F_1,\dots, F_{n/3}\}$ admits a rainbow matching, i.e., a matching consisting of one edge from each $F_i$. This is done by converting the rainbow matching problem to a perfect matching problem in a special class of uniform hypergraphs.


title: Response to Comment on "Low-frequency lattice phonons in halide perovskites explain high defect tolerance toward electron-hole recombination"

abstract: Recently we proposed that defect tolerance in the hybrid perovskites is due to their characteristic low-frequency lattice phonon modes that decrease the non-adiabatic coupling and weaken the overlap between the free carrier and defect states [Sci. Adv. 6 7, eaaw7453 (2020)]. Kim and Walsh disagree with the interpretation and argue that there are flaws in our employed methodology. Herein we address their concerns and show that their conclusions are not valid due to misunderstandings of nonadiabatic transition.


title: Interface-Induced Sign Reversal of the Anomalous Hall Effect in Magnetic Topological Insulator Heterostructures

abstract: The Berry phase picture provides important insights into the electronic properties of condensed matter systems. The intrinsic anomalous Hall (AH) effect can be understood as a consequence of non-zero Berry curvature in momentum space. The realization of the quantum anomalous Hall effect provided conclusive evidence for the intrinsic mechanism of the AH effect in magnetic topological insulators (TIs). Here we fabricated magnetic TI/TI heterostructures and found both the magnitude and sign of the AH effect in the magnetic TI layer can be altered by tuning the TI thickness and/or the electric gate voltage. The sign change of the AH effect with increasing TI thickness is attributed to the charge transfer across the TI and magnetic TI layers, consistent with first-principles calculations. By fabricating the magnetic TI/TI/magnetic TI sandwich heterostructures with different dopants, we created an artificial topological Hall (TH) effect-like feature in Hall traces. This artificial TH effect is induced by the superposition of two AH effects with opposite signs instead of the formation of chiral spin textures in the samples. Our study provides a new route to engineer the Berry curvature in magnetic topological materials that may lead to potential technological applications.


title: A better bound on the size of rainbow matchings

abstract: Aharoni and Howard conjectured that, for positive integers $n,k,t$ with $n\ge k$ and $n\ge t$, if $F_1,\ldots, F_t\subseteq {[n]\choose k}$ such that $|F_i|>{n\choose k}-{n-t+1\choose k}$ for $i\in [t]$ then there exist $e_i\in F_i$ for $i\in [t]$ such that $e_1,\ldots,e_t$ are pairwise disjoint. Huang, Loh, and Sudakov proved this conjecture for $t<n/(3k^2)$. In this paper, we show that this conjecture holds for $t\le n/(2k)$ and $n$ sufficiently large.


title: Temperature dependence of the anomalous Nernst effect in Ni$_{2}$MnGa shape memory alloy

abstract: We report a detailed investigation of the Ni$_{2}$MnGa shape memory alloy through magnetic, electronic, and thermal measurements. Our measurements of the anomalous Nernst effect (ANE) reveal that this technique is very sensitive to the onset of the pre-martensitic transition in sharp contrast to other transport measurements. With the ANE being sensitive to changes at the Fermi surface, we infer on the link between the structural modulations and the modulation of the Fermi surface via its nesting features, with the magnetic field induced strain being the mediating mechanism.


title: Automatic Textual Evidence Mining in COVID-19 Literature

abstract: We created this EVIDENCEMINER system for automatic textual evidence mining in COVID-19 literature. EVIDENCEMINER is a web-based system that lets users query a natural language statement and automatically retrieves textual evidence from a background corpora for life sciences. It is constructed in a completely automated way without any human effort for training data annotation. EVIDENCEMINER is supported by novel data-driven methods for distantly supervised named entity recognition and open information extraction. The named entities and meta-patterns are pre-computed and indexed offline to support fast online evidence retrieval. The annotation results are also highlighted in the original document for better visualization. EVIDENCEMINER also includes analytic functionalities such as the most frequent entity and relation summarization.


title: Counterexamples to conjectures by Gross, Mansour and Tucker on partial-dual genus polynomials of ribbon graphs

abstract: Gross, Mansour and Tucker introduced the partial-dual orientable genus polynomial and the partial-dual Euler genus polynomial. They computed these two partial-dual genus polynomials of four families of ribbon graphs, posed some research problems and made some conjectures. In this paper, we introduce the notion of signed sequences of bouquets and obtain the partial-dual Euler genus polynomials for all ribbon graphs with the number of edges less than 4 and the partial-dual orientable genus polynomials for all orientable ribbon graphs with the number of edges less than 5 in terms of signed sequences. We check all the conjectures and find a counterexample to the Conjecture 3.1 in their paper: There is no orientable ribbon graph having a non-constant partial-dual genus polynomial with only one non-zero coefficient. Motivated by this counterexample, we further find an infinite family of counterexamples to the conjecture. Moreover, we find a counterexample to the Conjecture 5.3 in their paper: The partial-dual Euler-genus polynomial for any non-orientable ribbon graph is interpolating.


title: SLSpy: Python-Based System-Level Controller Synthesis Framework

abstract: Synthesizing controllers for large, complex, and distributed systems is a challenging task. Numerous proposed methods exist in the literature, but it is difficult for practitioners to apply them -- most proposed synthesis methods lack ready-to-use software implementations, and existing proprietary components are too rigid to extend to general systems. To address this gap, we develop SLSpy, a framework for controller synthesis, comparison, and testing. SLSpy implements a highly extensible software framework which provides two essential workflows: synthesis and simulation. The workflows are built from five conceptual components that can be customized to implement a wide variety of synthesis algorithms and disturbance tests. SLSpy comes pre-equipped with a workflow for System Level Synthesis (SLS), which enables users to easily and freely specify desired design objectives and constraints. We demonstrate the effectiveness of SLSpy using two examples that have been described in the literature but do not have ready-to-use implementations. We open-source SLSpy to facilitate future controller synthesis research and practical usage.


title: Spontaneous scalarization of Gauss-Bonnet black holes surrounded by massive scalar fields

abstract: For massless scalar fields, a relation $\Delta_{n}=\frac{\sqrt{3}}{2}\pi$ for $n\rightarrow \infty$ was observed in the scalar-Gauss-Bonnet theory. In the present paper, we extend the discussion by including a nonzero scalar field mass. For massive scalar fields, we show that the relation $\Delta_{n}=\frac{\sqrt{3}}{2}\pi$ for $n\rightarrow \infty$ still holds. We demonstrate this relation with both analytical and numerical methods. The analytical analysis implies that this relation may be a very universal behavior.


title: System Level Synthesis via Dynamic Programming

abstract: System Level Synthesis (SLS) parametrization facilitates controller synthesis for large, complex, and distributed systems by incorporating system level constraints (SLCs) into a convex SLS problem and mapping its solution to stable controller design. Solving the SLS problem at scale efficiently is challenging, and current attempts take advantage of special system or controller structures to speed up the computation in parallel. However, those methods do not generalize as they rely on the specific system/controller properties. We argue that it is possible to solve general SLS problems more efficiently by exploiting the structure of SLS constraints. In particular, we derive dynamic programming (DP) algorithms to solve SLS problems. In addition to the plain SLS without any SLCs, we extend DP to tackle infinite horizon SLS approximation and entrywise linear constraints, which form a superclass of the locality constraints. Comparing to convex program solver and naive analytical derivation, DP solves SLS 4 to 12 times faster and scales with little computation overhead. We also quantize the cost of synthesizing a controller that stabilizes the system in a finite horizon through simulations.


title: Attenuation of Several Common Building Materials in Millimeter-Wave Frequency Bands: 28, 73 and 91 GHz

abstract: Future cellular systems will make use of millimeter wave (mmWave) frequency bands. Many users in these bands are located indoors, i.e., inside buildings, homes, and offices. Typical building material attenuations in these high frequency ranges are of interest for link budget calculations. In this paper, we report on a collaborative measurement campaign to find the attenuation of several typical building materials in three potential mmWave bands (28, 73, 91 GHz). Using directional antennas, we took multiple measurements at multiple locations using narrow-band and wide-band signals, and averaged out residual small-scale fading effects. Materials include clear glass, drywall (plasterboard), plywood, acoustic ceiling tile, and cinder blocks. Specific attenuations range from approximately 0.5 dB/cm for ceiling tile at 28 GHz to approximately 19 dB/cm for clear glass at 91 GHz.


title: DWT-GBT-SVD-based Robust Speech Steganography

abstract: Steganography is a method that can improve network security and make communications safer. In this method, a secret message is hidden in content like audio signals that should not be perceptible by listening to the audio or seeing the signal waves. Also, it should be robust against different common attacks such as noise and compression. In this paper, we propose a new speech steganography method based on a combination of Discrete Wavelet Transform, Graph-based Transform, and Singular Value Decomposition (SVD). In this method, we first find voiced frames based on energy and zero-crossing counts of the frames and then embed a binary message into voiced frames. Experimental results on the NOIZEUS database show that the proposed method is imperceptible and also robust against Gaussian noise, re-sampling, re-quantization, high pass filter, and low pass filter. Also, it is robust against MP3 compression and scaling for watermarking applications.


title: The Ingredients of Real-World Robotic Reinforcement Learning

abstract: The success of reinforcement learning for real world robotics has been, in many cases limited to instrumented laboratory scenarios, often requiring arduous human effort and oversight to enable continuous learning. In this work, we discuss the elements that are needed for a robotic learning system that can continually and autonomously improve with data collected in the real world. We propose a particular instantiation of such a system, using dexterous manipulation as our case study. Subsequently, we investigate a number of challenges that come up when learning without instrumentation. In such settings, learning must be feasible without manually designed resets, using only on-board perception, and without hand-engineered reward functions. We propose simple and scalable solutions to these challenges, and then demonstrate the efficacy of our proposed system on a set of dexterous robotic manipulation tasks, providing an in-depth analysis of the challenges associated with this learning paradigm. We demonstrate that our complete system can learn without any human intervention, acquiring a variety of vision-based skills with a real-world three-fingered hand. Results and videos can be found at https://sites.google.com/view/realworld-rl/


title: Exploiting Defenses against GAN-Based Feature Inference Attacks in Federated Learning

abstract: As a decentralized model training method, federated learning is designed to integrate the isolated data islands and protect data privacy. Recent studies, however, have demonstrated that the Generative Adversarial Network (GAN) based attacks can be used in federated learning to learn the distribution of the victim's private dataset and accordingly reconstruct human-distinguishable images. In this paper, we exploit defenses against GAN-based attacks in federated learning, and propose a framework, Anti-GAN, to prevent attackers from learning the real distribution of the victim's data. The core idea of Anti-GAN is to corrupt the visual features of the victim's private training images, such that the images restored by the attacker are indistinguishable to human eyes. Specifically, in Anti-GAN, the victim first projects the personal dataset onto a GAN's generator, then mixes the fake images generated by the generator with the real images to obtain the training dataset, which will be fed into the federated model for training. We redesign the structure of the victim's GAN to encourage it to learn the classification features (instead of the visual features) of the real images. We further introduce an unsupervised task to the GAN model for obfuscating the visual features of the generated images. The experiments demonstrate that Anti-GAN can effectively prevent the attacker from learning the distribution of the private images, meanwhile causing little harm to the accuracy of the federated model.


title: Strong Backward uniqueness for sublinear parabolic equations

abstract: In this paper, we establish strong backward uniqueness for solutions to sublinear parabolic equations of the type (1.1). The proof of our main result Theorem 1.1 is achieved by means of a new Carleman estimate and a Weiss type monotonicity that are tailored for such parabolic sublinear operators.


title: Bayesian Learning of Probabilistic Dipole Inversion for Quantitative Susceptibility Mapping

abstract: A learning-based posterior distribution estimation method, Probabilistic Dipole Inversion (PDI), is proposed to solve quantitative susceptibility mapping (QSM) inverse problem in MRI with uncertainty estimation. A deep convolutional neural network (CNN) is used to represent the multivariate Gaussian distribution as the approximated posterior distribution of susceptibility given the input measured field. In PDI, such CNN is firstly trained on healthy subjects dataset with labels by maximizing the posterior Gaussian distribution loss function as used in Bayesian deep learning. When tested on new dataset without any label, PDI updates the pre-trained network in an unsupervised fashion by minimizing the KL divergence between the approximated posterior distribution represented by CNN and the true posterior distribution given the likelihood distribution from known physical model and prior distribution. Based on our experiments, PDI provides additional uncertainty estimation compared to the conventional MAP approach, meanwhile addressing the potential discrepancy issue of CNN when test data deviates from training dataset.


title: Single- and Multi-Objective Evolutionary Algorithms for the Knapsack Problem with Dynamically Changing Constraints

abstract: Evolutionary algorithms are bio-inspired algorithms that can easily adapt to changing environments. Recent results in the area of runtime analysis have pointed out that algorithms such as the (1+1)~EA and Global SEMO can efficiently reoptimize linear functions under a dynamic uniform constraint. Motivated by this study, we investigate single- and multi-objective baseline evolutionary algorithms for the classical knapsack problem where the capacity of the knapsack varies over time. We establish different benchmark scenarios where the capacity changes every $\tau$ iterations according to a uniform or normal distribution. Our experimental investigations analyze the behavior of our algorithms in terms of the magnitude of changes determined by parameters of the chosen distribution, the frequency determined by $\tau$, and the class of knapsack instance under consideration. Our results show that the multi-objective approaches using a population that caters for dynamic changes have a clear advantage on many benchmarks scenarios when the frequency of changes is not too high. Furthermore, we demonstrate that the diversity mechanisms used in popular evolutionary multi-objective algorithms such as NSGA-II and SPEA2 do not necessarily result in better performance and even lead to inferior results compared to our simple multi-objective approaches.


title: Adaptive model selection in photonic reservoir computing by reinforcement learning

abstract: Photonic reservoir computing is an emergent technology toward beyond-Neumann computing. Although photonic reservoir computing provides superior performance in environments whose characteristics are coincident with the training datasets for the reservoir, the performance is significantly degraded if these characteristics deviate from the original knowledge used in the training phase. Here, we propose a scheme of adaptive model selection in photonic reservoir computing using reinforcement learning. In this scheme, a temporal waveform is generated by different dynamic source models that change over time. The system autonomously identifies the best source model for the task of time series prediction using photonic reservoir computing and reinforcement learning. We prepare two types of output weights for the source models, and the system adaptively selected the correct model using reinforcement learning, where the prediction errors are associated with rewards. We succeed in adaptive model selection when the source signal is temporally mixed, having originally been generated by two different dynamic system models, as well as when the signal is a mixture from the same model but with different parameter values. This study paves the way for autonomous behavior in photonic artificial intelligence and could lead to new applications in load forecasting and multi-objective control, where frequent environment changes are expected.


title: How to Return to Normalcy: Fast and Comprehensive Contact Tracing of COVID-19 through Proximity Sensing Using Mobile Devices

abstract: We outline a contact-tracing strategy based on proximity sensing using mobile devices. We discuss what an ideal system should look like and what it can do. We show that, when adopted sufficiently broadly, such a contact-tracing strategy can bring COVID-19 under complete control, end the need of social distancing, and return the society to full normalcy. We also review some of the challenges faced by the current generation of proximity-sensing technologies, including Bluetooth Low Energy used by phones, and consider both interim and longer-term solutions. Our main contribution is that we reason through why such a contact-tracing strategy is likely to achieve the stated goal of returning to full normalcy. Using probabilistic models, we show that universal adoption is not necessary to achieve the stated goal, thus there is some room for exceptions; however, the adoption rate needs to be very high, e.g., above $95\%$ depending on the disease parameters. With more vigilance in disease surveillance to detect mild cases earlier, the number may be brought down to about $90\%$. The results call for deployment effort to be led by public authorities at the state or federal level so that the required adoption rate can be reached and the tracing coverage is wide enough to be relevant for disease control.


title: Supersolidity in an elongated dipolar condensate

abstract: We present a theory for the emergence of a supersolid state in a cigar-shaped dipolar quantum Bose gas. Our approach is based on a reduced three-dimensional (3D) theory, where the condensate wavefunction is decomposed into an axial field and a transverse part described variationally. This provides an accurate fully 3D description that is specific to the regime of current experiments and efficient to compute. We apply this theory to understand the phase diagram for a gas in an infinite tube potential. We find that the supersolid transition has continuous and discontinuous regions as the averaged density varies. We develop two simplified analytic models to characterize the phase diagram and elucidate the roles of quantum droplets and of the roton excitation.


title: Weak compactness criteria in $L_1$ space in terms of Orlicz function

abstract: In this paper, we provide a direct proof for the equivalence of K.M. Chong's and De la Vall\'{e}e Poussin's criteria of weak compactness of a subset $K$ of $L_1(0,1)$ in terms of some Orlicz function. Furthermore, we discuss the equivalence in $L_1(0, \infty)$.


title: Power spectrum of halo intrinsic alignments in simulations

abstract: We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $\Lambda$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter of galaxy or halo density field. We find that the cross- and auto-power spectra $P_{\delta E}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($\Omega_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Thus our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be used to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.


title: A Phase Transition in Arrow's Theorem

abstract: Arrow's Theorem concerns a fundamental problem in social choice theory: given the individual preferences of members of a group, how can they be aggregated to form rational group preferences? Arrow showed that in an election between three or more candidates, there are situations where any voting rule satisfying a small list of natural "fairness" axioms must produce an apparently irrational intransitive outcome. Furthermore, quantitative versions of Arrow's Theorem in the literature show that when voters choose rankings in an i.i.d.\ fashion, the outcome is intransitive with non-negligible probability. It is natural to ask if such a quantitative version of Arrow's Theorem holds for non-i.i.d.\ models. To answer this question, we study Arrow's Theorem under a natural non-i.i.d.\ model of voters inspired by canonical models in statistical physics; indeed, a version of this model was previously introduced by Raffaelli and Marsili in the physics literature. This model has a parameter, temperature, that prescribes the correlation between different voters. We show that the behavior of Arrow's Theorem in this model undergoes a striking phase transition: in the entire high temperature regime of the model, a Quantitative Arrow's Theorem holds showing that the probability of paradox for any voting rule satisfying the axioms is non-negligible; this is tight because the probability of paradox under pairwise majority goes to zero when approaching the critical temperature, and becomes exponentially small in the number of voters beyond it. We prove this occurs in another natural model of correlated voters and conjecture this phenomena is quite general.


title: Multi-IF : An Approach to Anomaly Detection in Self-Driving Systems

abstract: Autonomous driving vehicles (ADVs) are implemented with rich software functions and equipped with many sensors, which in turn brings broad attack surface. Moreover, the execution environment of ADVs is often open and complex. Hence, ADVs are always at risk of safety and security threats. This paper proposes a fast method called Multi-IF, using multiple invocation features of system calls to detect anomalies in self-driving systems. Since self-driving functions take most of the computation resources and upgrade frequently, Multi-IF is designed to work under such resource constraints and support frequent updates. Given the collected sequences of system calls, the combination of different syntax patterns is used to analyze and construct feature vectors of those sequences. By taking the feature vectors as inputs, one-class support vector machine is adopted to determine whether the current sequence of system calls is abnormal, which is trained with the feature vectors from the normal sequences. The evaluations on both simulated and real data prove that the proposed method is effective in identifying the abnormal behavior after minutes of feature extraction and training. Further comparisons with the existing methods on the ADFA-LD data set also validate that the proposed approach achieves a higher accuracy with less time overhead.


title: Fixed points of compositions of nonexpansive mappings: finitely many linear reflectors

abstract: Nonexpansive mappings play a central role in modern optimization and monotone operator theory because their fixed points can describe solutions to optimization or critical point problems. It is known that when the mappings are sufficiently "nice", then the fixed point set of the composition coincides with the intersection of the individual fixed point sets. In this paper, we explore the situation for compositions of linear reflectors. We provide positive results, upper bounds, and limiting examples. We also discuss classical reflectors in the Euclidean plane.


title: Nonlinear nonhomogeneous singular problems

abstract: We consider a nonlinear Dirichlet problem driven by a nonhomogeneous differential operator with a growth of order $(p-1)$ near $+\infty$ and with a reaction which has the competing effects of a parametric singular term and a $(p-1)$-superlinear perturbation which does not satisfy the usual Ambrosetti-Rabinowitz condition. Using variational tools, together with suitable truncation and strong comparison techniques, we prove a "bifurcation-type" theorem that describes the set of positive solutions as the parameter $\lambda$ moves on the positive semiaxis. We also show that for every $\lambda>0$, the problem has a smallest positive solution $u^*_\lambda$ and we demonstrate the monotonicity and continuity properties of the map $\lambda\mapsto u^*_\lambda$.


title: Viral air load due to sedimenting and evaporating droplets produced by speaking

abstract: The effect of evaporation on droplet sedimentation times is crucial for estimating the risk of infection from virus-containing airborne droplets. For droplet radii in the range 100 nm < R < 60 {\mu}m, evaporation can be described in the stagnant-flow approximation and is diffusion limited. Analytical equations are presented for the droplet evaporation rate, the time-dependent droplet size and the sedimentation time, including the significant effect of evaporation cooling. Evaporation makes the time for large droplets to sediment much longer and thus significantly increases the viral air load. Using recent estimates for SARS-CoV-2 concentrations in sputum and droplet production rates while speaking, a single infected person that constantly speaks without a mouth cover produces a total air load of more than 10^4 virions. In a mid-size closed room, this leads to a viral inhalation frequency of at least 2.5 per minute. Low relative humidity, as encountered inside buildings in winter and in airliners, speeds evaporation and thus keeps initially larger droplets suspended in air.


title: A Batch Normalized Inference Network Keeps the KL Vanishing Away

abstract: Variational Autoencoder (VAE) is widely used as a generative model to approximate a model's posterior on latent variables by combining the amortized variational inference and deep neural networks. However, when paired with strong autoregressive decoders, VAE often converges to a degenerated local optimum known as "posterior collapse". Previous approaches consider the Kullback Leibler divergence (KL) individual for each datapoint. We propose to let the KL follow a distribution across the whole dataset, and analyze that it is sufficient to prevent posterior collapse by keeping the expectation of the KL's distribution positive. Then we propose Batch Normalized-VAE (BN-VAE), a simple but effective approach to set a lower bound of the expectation by regularizing the distribution of the approximate posterior's parameters. Without introducing any new model component or modifying the objective, our approach can avoid the posterior collapse effectively and efficiently. We further show that the proposed BN-VAE can be extended to conditional VAE (CVAE). Empirically, our approach surpasses strong autoregressive baselines on language modeling, text classification and dialogue generation, and rivals more complex approaches while keeping almost the same training time as VAE.


title: Anisotropy driven reversal of magnetisation in Blume-Capel ferromagnet: A Monte Carlo study

abstract: The two dimensional Spin-1 Blume-Capel ferromagnet is studied by Monte Carlo simulation with Metropolis algorithm. Starting from initial ordered spin configuration the reversal of magnetisation is investigated in presence of a magnetic field ($h$) applied in the opposite direction. The variations of the reversal time with the strength of single site anisotropy are investigated in details. The exponential dependence was observed. The systematic variations of the mean reversal time with positive and negative anisotropy was found. The mean macroscopic reversal time was observed to be linearly dependent on a suitably defined microscopic reversal time. The saturated magnetisation $M_f$ after the reversal was noticed to be dependent of the strength of anisotropy $D$. An interesting scaling relation was obtained, $|M_f| \sim |h|^{\beta}f(D|h|^{-\alpha})$ with the scaling function of the form $f(x)= \frac{1}{1+e^{(x-a)/b}}$. The temporal evolution of density of $S_i^z=0$ (surrounded by all $S_i^z=+1$) is observed to be exponentially decaying. The growth of mean density of $S_i^z=-1$ has been fitted in a function $\rho_{-1}(t) \sim \frac{1}{a+e^{(t_c-t)/c}}$. The characteristic time shows $t_c \sim e^{-rD}$ and a crossover in the rate of exponential falling is observed at $D=1.5$. The metastable volume fraction has been found to obey the Avrami's law.


title: Binary MIMO Detection via Homotopy Optimization and Its Deep Adaptation

abstract: In this paper we consider maximum-likelihood (ML) MIMO detection under one-bit quantized observations and binary symbol constellations. This problem is motivated by the recent interest in adopting coarse quantization in massive MIMO systems--as an effective way to scale down the hardware complexity and energy consumption. Classical MIMO detection techniques consider unquantized observations, and many of them are not applicable to the one-bit MIMO case. We develop a new non-convex optimization algorithm for the one-bit ML MIMO detection problem, using a strategy called homotopy optimization. The idea is to transform the ML problem into a sequence of approximate problems, from easy (convex) to hard (close to ML), and with each problem being a gradual modification of its previous. Then, our attempt is to iteratively trace the solution path of these approximate problems. This homotopy algorithm is well suited to the application of deep unfolding, a recently popular approach for turning certain model-based algorithms into data-driven, and performance enhanced, ones. While our initial focus is on one-bit MIMO detection, the proposed technique also applies naturally to the classical unquantized MIMO detection. We performed extensive simulations and show that the proposed homotopy algorithms, both non-deep and deep, have satisfactory bit-error probability performance compared to many state-of-the-art algorithms. Also, the deep homotopy algorithm has attractively low computational complexity.