title: On cutting blocking sets and their codes

abstract: Let PG$(r, q)$ be the $r$-dimensional projective space over the finite field ${\rm GF}(q)$. A set $\cal X$ of points of PG$(r, q)$ is a cutting blocking set if for each hyperplane $\Pi$ of PG$(r, q)$ the set $\Pi \cap \cal X$ spans $\Pi$. Cutting blocking sets give rise to saturating sets and minimal linear codes and those having size as small as possible are of particular interest. We observe that from a cutting blocking set obtained by Fancsali and Sziklai, by using a set of pairwise disjoint lines, there arises a minimal linear code whose length grows linearly with respect to its dimension. We also provide two distinct constructions: a cutting blocking set of PG$(3, q^3)$ of size $3(q+1)(q^2+1)$ as a union of three pairwise disjoint $q$-order subgeometries and a cutting blocking set of PG$(5, q)$ of size $7(q+1)$ from seven lines of a Desarguesian line spread of PG$(5, q)$. In both cases the cutting blocking sets obtained are smaller than the known ones. As a byproduct we further improve on the upper bound of the smallest size of certain saturating sets and on the minimum length of a minimal $q$-ary linear code having dimension $4$ and $6$.


title: Structural study of TATB under detonation-induced shock conditions

abstract: We explore the response of the insensitive high explosive (IHE) 1,3,5-Triamino-2,4,6- trinitrobenzene (TATB) under detonation-induced shock conditions using in-situ synchrotron X-ray diffraction in the 100 ns time scale using either a conventional or a colliding detonation drive. In all of the detonation experiments on various sizes and morphologies of TATB, we observe an extended stability of the TATB triclinic crystal structure. As the detonation front passes through the TATB, X-ray diffraction indicates a portion of the TATB exhibits a compression up to 30+ GPa, followed subsequently by a pressure release and continued decomposition over a few hundred nanoseconds. Likewise, for colliding detonation-driven shock compression of single crystals of TATB, a significant portion of the triclinic crystal structure appears to be stable up to 60+ GPa. Conversely, in similar detonations of a LLM-105 PBX, X-ray diffraction is simply indicative of decomposition without the apparent compression and slow decomposition seen in TATB. The results indicate the surprising resilience of TATB under these high-pressure, temperature and shock conditions, providing a baseline for understanding the insensitivity of TATB that is considered the industry standard for an insensitive high explosive. The results also provide intriguing information for the extended reaction zone in TATB, and the hot-spot mechanisms for initiating and propagating detonation in this uniquely insensitive explosive.


title: Measurement of Proton Quenching in a Plastic Scintillator Detector

abstract: The non-linear energy response of the plastic scintillator EJ-260 is measured with the MicroCHANDLER detector, using neutron beams of energy 5 to 27 MeV at the Triangle Universities Nuclear Laboratory. The first and second order Birks' constants are extracted from the data, and found to be $k_B = (8.70 \pm 0.93)\times 10^{-3}\ {\rm g/cm^2/MeV}$ and $k_C = (1.42 \pm 1.00) \times 10^{-5}\ {\rm (g/cm^2/MeV)^2}$. This result covers a unique energy range that is of direct relevance for fast neutron backgrounds in reactor inverse beta decay detectors. These measurements will improve the energy non-linearity modeling of plastic scintillator detectors. In particular, the updated energy response model will lead to an improvement of fast neutron modeling for detectors based on the CHANDLER reactor neutrino detector technology.


title: Model Predictive Control for Micro Aerial Vehicles: A Survey

abstract: This paper presents a review of the design and application of model predictive control strategies for Micro Aerial Vehicles and specifically multirotor configurations such as quadrotors. The diverse set of works in the domain is organized based on the control law being optimized over linear or nonlinear dynamics, the integration of state and input constraints, possible fault-tolerant design, if reinforcement learning methods have been utilized and if the controller refers to free-flight or other tasks such as physical interaction or load transportation. A selected set of comparison results are also presented and serve to provide insight for the selection between linear and nonlinear schemes, the tuning of the prediction horizon, the importance of disturbance observer-based offset-free tracking and the intrinsic robustness of such methods to parameter uncertainty. Furthermore, an overview of recent research trends on the combined application of modern deep reinforcement learning techniques and model predictive control for multirotor vehicles is presented. Finally, this review concludes with explicit discussion regarding selected open-source software packages that deliver off-the-shelf model predictive control functionality applicable to a wide variety of Micro Aerial Vehicle configurations.


title: Investigation of Charged Higgs Boson in the Bottom and Top Quark Decay Channel at the FCC-hh

abstract: After the recent discovery of a neutral Higgs boson with a mass about 125 GeV, we assess the extend of discovery potential of future circular hadron collider (FCC-hh) for a charged Higgs boson in the bottom and top quark decay channel. The charged Higgs boson can be produced through the pp->h^{-}t+X process with a subsequent decay h^{-}->b bar{t} channel. This decay channel is particularly important for studying the charged Higgs boson heavier than the top quark. We consider an extension of the standard model Higgs sector, namely two Higgs doublet model (2HDM), and perform a dedicated signal significance analysis to test this channel for the FCC-hh running at the center of mass energy of 100 TeV and the integrated luminosity of 1 ab^{-1} (initial) and 30 ab^{-1} (ultimate). We find that an important part of the parameter spaces of two Higgs doublet model are examinable at the FCC-hh.


title: Successive common envelope events from multiple planets

abstract: Many stars harbour multi-planet systems. As these stars expand late in their evolutions, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this is insufficient to eject the envelope, it may expand the star further, causing additional CE events, with the last one unbinding what remains of the envelope. This multi-planet CE scenario may have broad implications for stellar and planetary evolution across a range of systems. We develop a simplified version and show that it may be able to explain the recently observed planet WD 1856 b.


title: The $\operatorname{Ext}$-algebra of standard modules over dual extension algebras

abstract: We exhibit an isomorphism of associative algebras between the $\operatorname{Ext}$-algebra $\operatorname{Ext}_\Lambda^\ast(\Delta,\Delta)$ of standard modules over the dual extension algebra $\Lambda$ of two directed algebras $B$ and $A$ and the dual extension algebra of the $\operatorname{Ext}$-algebra $\operatorname{Ext}_B^\ast(\mathbb{L},\mathbb{L})$ with $A$. There are natural $A_\infty$-structures on these $\operatorname{Ext}$-algebras, and, under certain technical assumptions on $B$, we describe that on $\operatorname{Ext}_\Lambda^\ast(\Delta,\Delta)$ completely in terms of that on $\operatorname{Ext}_B^\ast(\mathbb{L},\mathbb{L})$. As an example, we compute these $A_\infty$-structures explicitly in the case where $B=A=K\mathbb{A}_n /(\operatorname{rad}K\mathbb{A}_n)^\ell$.


title: Multiresolution Knowledge Distillation for Anomaly Detection

abstract: Unsupervised representation learning has proved to be a critical component of anomaly detection/localization in images. The challenges to learn such a representation are two-fold. Firstly, the sample size is not often large enough to learn a rich generalizable representation through conventional techniques. Secondly, while only normal samples are available at training, the learned features should be discriminative of normal and anomalous samples. Here, we propose to use the "distillation" of features at various layers of an expert network, pre-trained on ImageNet, into a simpler cloner network to tackle both issues. We detect and localize anomalies using the discrepancy between the expert and cloner networks' intermediate activation values given the input data. We show that considering multiple intermediate hints in distillation leads to better exploiting the expert's knowledge and more distinctive discrepancy compared to solely utilizing the last layer activation values. Notably, previous methods either fail in precise anomaly localization or need expensive region-based training. In contrast, with no need for any special or intensive training procedure, we incorporate interpretability algorithms in our novel framework for the localization of anomalous regions. Despite the striking contrast between some test datasets and ImageNet, we achieve competitive or significantly superior results compared to the SOTA methods on MNIST, F-MNIST, CIFAR-10, MVTecAD, Retinal-OCT, and two Medical datasets on both anomaly detection and localization.


title: A Mathematical Dashboard for the Analysis of Italian COVID-19 Epidemic Data

abstract: An analysis of the COVID-19 epidemic is proposed on the basis of the epiMOX dashboard (publicly accessible at https://www.epimox.polimi.it) that deals with data of the epidemic trends and outbreaks in Italy from late February 2020. Our analysis provides an immediate appreciation of the past epidemic development, together with its current trends by fostering a deeper interpretation of available data through several critical epidemic indicators. In addition, we complement the epiMOX dashboard with a predictive tool based on an epidemiological compartmental model, named SUIHTER, for the forecast on the near future epidemic evolution.


title: Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers

abstract: We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear equations of state lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.


title: Photophysics of single nitrogen-vacancy centers in nanodiamonds coupled to photonic crystal cavities

abstract: The nitrogen vacancy center in diamond in its negative charge state is a promising candidate for quantum optic experiments that require single photon emitters. Important benefits of the NV center are its high brightness and photo-stability, even at room temperature. Engineering the emission properties of NV centers with optical resonators is a widely followed approach to meet the requirements for quantum technological applications, but the effect on non-radiative decay paths is yet to be understood. Here we report on modifying the internal quantum efficiency of a single NV center in a nanodiamond coupled to a 1D photonic crystal cavity. We assess the Purcell enhancement via three independent measurement techniques and perform autocorrelation measurements at elevated excitation powers in confocal microscopy. Employing a three-level model allows us to extract the setup efficiency, individual transition rates and thus the internal quantum efficiency of our system. Combining our results, we find that the enhancement of the radiative decay rate via the Purcell effect results in an internal quantum efficiency of 90 % for cavity-coupled NV centers. Our findings will facilitate the realization of nano-scale single photon sources with near-unity internal quantum efficiencies operating at high repetition rates.


title: An application of cyberpsychology in business email compromise

abstract: This paper introduces Business Email Compromise (BEC) and why it is becoming a major issue to businesses worldwide. It also presents a case study of a BEC incident against Unatrac Holding Ltd and analyses the techniques used by the cybercriminals to defraud the company. A critical analysis of the psychological and sociotechnical impacts of BEC to both the company and employees are conducted, and potential risk mitigations strategies and recommendations are provided to prevent future attacks.


title: {\it Gaia} DR2 giants in the Galactic dust -- I. Reddening across the whole dust layer and some properties of the giant clump

abstract: We consider a complete sample of 101810 giants with Gaia DR2 parallaxes $\varpi$ within the red clump domain of the Hertzsprung-Russell diagram in the space cylinder with a radius of 700 pc around the Sun and a height of $|Z|=1800$ pc. We use the Gaia DR2 $G_\mathrm{BP}$, $G_\mathrm{RP}$ and WISE $W3$ photometry. We describe the spatial variations of the modes of the observables $G_\mathrm{BP}-G_\mathrm{RP}$, $G_\mathrm{RP}-W3$, $G_\mathrm{BP}+5+5\,\log_{10}\varpi$, $G_\mathrm{RP}+5+5\,\log_{10}\varpi$, and $W3+5+5\,\log_{10}\varpi$ by extinction and reddening in combination with linear vertical gradients of the intrinsic colours and absolute magnitudes of the red giant clump. The derived clump median absolute magnitude in $W3$ agrees with its recent literature estimates. The clump median intrinsic colours and absolute magnitudes in $G_\mathrm{BP}$ and $G_\mathrm{RP}$ are derived for the first time at a precision level of 0.01 mag. We confirm the reliability of the derived clump absolute magnitudes, intrinsic colours, and their vertical gradients by comparing them with the theoretical predictions from PARSEC, MIST and BaSTI isochrones. This leads us to the median age and [Fe/H] of the clump within $|Z|<1.7$ kpc from the Galactic mid-plane as $(2.3\pm0.5)+(3.2\pm1.6)\,|Z|$ Gyr and $(-0.08\pm0.08)-(0.16\pm0.07)\,|Z|$ dex, respectively, where $Z$ is expressed in kpc. These results agree with recent empirical and theoretical estimates. Moreover, all the models give similar age-metallicity relations by use of our results in the optical range. The derived extinctions and reddenings across the whole dust half-layer below or above the Sun converge to the reddening $E(B-V)=0.06$ mag by use of the most reliable extinction laws.


title: Analogue time machine in a photonic system

abstract: Analogue physics has successfully tackled the problems of chromodynamics, event horizons, Big Bang and Universe expansion, and many others. Here, we suggest a photonic model system for a "time machine" based on the paraxial beam approximation. We demonstrate how the closed time-like curves and the well-known grandfather paradox can be studied experimentally in this system. We show how the Novikov's self-consistency principle is realized in quantum mechanics thanks to the Heisenberg's uncertainty principle.


title: Employing distributional semantics to organize task-focused vocabulary learning

abstract: How can a learner systematically prepare for reading a book they are interested in? In this paper,we explore how computational linguistic methods such as distributional semantics, morphological clustering, and exercise generation can be combined with graph-based learner models to answer this question both conceptually and in practice. Based on the highly structured learner model and concepts from network analysis, the learner is guided to efficiently explore the targeted lexical space. They practice using multi-gap learning activities generated from the book focused on words that are central to the targeted lexical space. As such the approach offers a unique combination of computational linguistic methods with concepts from network analysis and the tutoring system domain to support learners in achieving their individual, reading task-based learning goals.


title: $\mathbb{Z}_2$-graded polynomial identities for the Jordan algebra of $2\times 2$ upper triangular matrices

abstract: Let $K$ be a field (finite or infinite) of char$(K)\neq 2$ and let $UT_n=UT_n(K)$ be the $n\times n$ upper triangular matrix algebra over $K$. If $\cdot $ is the usual product on $UT_n$ then with the new product $a\circ b=(1/2)(a\cdot b +b\cdot a)$ we have that $UT_n$ is a Jordan algebra, denoted by $UJ_n=UJ_n(K)$. In this paper, we describe the set of all $\mathbb{Z}_2$-graded polynomial identities of $UJ_2$ with any nontrivial $\mathbb{Z}_2$-grading. Moreover, we describe a linear basis for the corresponding relatively free $\mathbb{Z}_2$-graded algebra.


title: Online Orthogonal Matching Pursuit

abstract: Greedy algorithms for feature selection are widely used for recovering sparse high-dimensional vectors in linear models. In classical procedures, the main emphasis was put on the sample complexity, with little or no consideration of the computation resources required. We present a novel online algorithm: Online Orthogonal Matching Pursuit (OOMP) for online support recovery in the random design setting of sparse linear regression. Our procedure selects features sequentially, alternating between allocation of samples only as needed to candidate features, and optimization over the selected set of variables to estimate the regression coefficients. Theoretical guarantees about the output of this algorithm are proven and its computational complexity is analysed.


title: Peer-to-Peer Energy Systems for Connected Communities: A Review of Recent Advances and Emerging Challenges

abstract: After a century of relative stability of the electricity industry, extensive deployment of distributed energy resources and recent advances in computation and communication technologies have changed the nature of how we consume, trade, and apply energy. The power system is facing a transition from its traditional hierarchical structure to a more deregulated model by introducing new energy distribution models such as peer-to-peer sharing for connected communities. The proven effectiveness of P2P sharing in benefiting both prosumers and the grid has been demonstrated in many studies and pilot projects. However, there is still no extensive implementation of such sharing models in today's electricity markets. This paper aims to shed some light on this gap through a comprehensive overview of recent advances in the P2P energy system and an insightful discussion of the challenges that need to be addressed in order to establish P2P sharing as a viable energy management option in today's electricity market. To this end, in this article, we provide some background on different aspects of P2P sharing. Then, we discuss advances in P2P sharing through a systematic domain-based classification. We also review different pilot projects on P2P sharing across the globe. Finally, we identify and discuss a number of challenges that need to be addressed for scaling up P2P sharing in the electricity market followed by concluding remarks at the end of the paper.


title: The balancing number and list balancing number of some graph classes

abstract: Given a graph $G$, a 2-coloring of the edges of $K_n$ is said to contain a balanced copy of $G$ if we can find a copy of $G$ such that half of its edges is in each color class. If there exists an integer $k$ such that, for $n$ sufficiently large, every 2-coloring of $K_n$ with more than $k$ edges in each color contains a balanced copy of $G$, then we say that $G$ is balanceable. The smallest integer $k$ such that this holds is called the balancing number of $G$. In this paper, we define a more general variant of the balancing number, the list balancing number, by considering 2-list edge colorings of $K_n$, where every edge $e$ has an associated list $L(e)$ which is a nonempty subset of the color set $\{r,b\}$. In this case, edges $e$ with $L(e) = \{r,b\}$ act as jokers in the sense that their color can be chosen $r$ or $b$ as needed. In contrast to the balancing number, every graph has a list balancing number. Moreover, if the balancing number exists, then it coincides with the list balancing number. We give the exact value of the list balancing number for all cycles except for $4k$-cycles for which we give tight bounds. In addition, we give general bounds for the list balancing number of non-balanceable graphs based on the extremal number of its subgraphs, and study the list balancing number of $K_5$, which turns out to be surprisingly large.


title: On the implementation of large-scale integral operators with modern HPC solutions -- Application to 3D Marchenko imaging by least-squares inversion

abstract: Numerical integral operators of convolution type form the basis of most wave-equation-based methods for processing and imaging of seismic data. As several of these methods require the solution of an inverse problem, multiple forward and adjoint passes of the modelling operator must be performed to converge to a satisfactory solution. This work highlights the challenges that arise when implementing such operators on 3D seismic datasets and it provides insights into their usage for solving large systems of integral equations. A Python framework is presented that leverages libraries for distributed storage and computing, and provides an high-level symbolic representation of linear operators. To validate its effectiveness, the forward and adjoint implementations of a multi-dimensional convolution operator are evaluated with respect to increasing size of the kernel and number of computational resources. Our computational framework is further shown to be suitable for both classic on-premise High-Performance Computing and cloud computing architectures. An example of target-oriented imaging of a 3D synthetic dataset which comprises of two subsequent steps of seismic redatuming is finally presented. In both cases, the redatumed fields are estimated by means of least-squares inversion using the full dataset as well as spatially decimated versions of the dataset as a way to investigate the robustness of both inverse problems to spatial aliasing in the input dataset. We observe that less strict sampling requirements apply in three dimensions for these algorithms compared to their two dimensions counterparts. Whilst aliasing introduces noise in the redatumed fields, they are however deprived of the well-known spurious artefacts arising from incorrect handling of the overburden propagation in cheaper, adjoint-based redatuming techniques.


title: Using Ensemble Analysis to study the effects of Network Topology on Performance in Urban Road Networks

abstract: Studies of the effect of network structure on performance have, thus far, been restricted to examining ensembles of synthetic networks generated by canonical models from the Network Science literature, which do not plausibly represent real road networks. Furthermore disparate ranges of parameter settings for demand and supply structure are often used, making it difficult to compare one study with another. To address these deficiencies, this paper proposes an experimental framework for the design of numerical investigations to study how network performance varies with respect to the structural properties of urban road networks. The key step within this experimental framework is the generation of a spectrum of ensembles of synthetic networks, which vary with respect to a single aspect of network structure and which adhere to plausibility constraints that ensure relevance to real urban road networks. This paper then demonstrates the application of this experimental framework to an investigation of how two performance indicators; the average link Volume-to-Capacity ratio and the Price of Anarchy, vary with respect to the density of travel demand, and the size, density and connectivity of network infrastructure. In so doing, several key challenges of such investigations are identified, with respect to the computational burden of the experiments and the need for empirical data on some important aspects of network structure such as the distribution of capacity.


title: Controlling symmetries and clustered dynamics of complex networks

abstract: Symmetries are an essential feature of complex networks as they regulate how the graph collective dynamics organizes into clustered states. We here show how to control network symmetries, and how to enforce patterned states of synchronization with nodes clustered in a desired way. Our approach consists of perturbing the original network connectivity, either by adding new edges or by adding/removing links together with modifying their weights. By solving suitable optimization problems, we furthermore guarantee that changes made on the existing topology are minimal. The conditions for the stability of the enforced pattern are derived for the general case, and the performance of the method is illustrated with paradigmatic examples. Our results are relevant to all the practical situations in which coordination of the networked systems into diverse groups may be desirable, such as for teams of robots, unmanned autonomous vehicles, power grids and central pattern generators.


title: Data Driven Robust Estimation Methods for Fixed Effects Panel Data Models

abstract: The panel data regression models have gained increasing attention in different areas of research including but not limited to econometrics, environmental sciences, epidemiology, behavioral and social sciences. However, the presence of outlying observations in panel data may often lead to biased and inefficient estimates of the model parameters resulting in unreliable inferences when the least squares (LS) method is applied. We propose extensions of the M-estimation approach with a data-driven selection of tuning parameters to achieve desirable level of robustness against outliers without loss of estimation efficiency. The consistency and asymptotic normality of the proposed estimators have also been proved under some mild regularity conditions. The finite sample properties of the existing and proposed robust estimators have been examined through an extensive simulation study and an application to macroeconomic data. Our findings reveal that the proposed methods often exhibits improved estimation and prediction performances in the presence of outliers and are consistent with the traditional LS method when there is no contamination.


title: Uncorrelated Semi-paired Subspace Learning

abstract: Multi-view datasets are increasingly collected in many real-world applications, and we have seen better learning performance by existing multi-view learning methods than by conventional single-view learning methods applied to each view individually. But, most of these multi-view learning methods are built on the assumption that at each instance no view is missing and all data points from all views must be perfectly paired. Hence they cannot handle unpaired data but ignore them completely from their learning process. However, unpaired data can be more abundant in reality than paired ones and simply ignoring all unpaired data incur tremendous waste in resources. In this paper, we focus on learning uncorrelated features by semi-paired subspace learning, motivated by many existing works that show great successes of learning uncorrelated features. Specifically, we propose a generalized uncorrelated multi-view subspace learning framework, which can naturally integrate many proven learning criteria on the semi-paired data. To showcase the flexibility of the framework, we instantiate five new semi-paired models for both unsupervised and semi-supervised learning. We also design a successive alternating approximation (SAA) method to solve the resulting optimization problem and the method can be combined with the powerful Krylov subspace projection technique if needed. Extensive experimental results on multi-view feature extraction and multi-modality classification show that our proposed models perform competitively to or better than the baselines.


title: Deep exclusive electroproduction of $\pi^0$ at high $Q^2$ in the quark valence regime

abstract: We report measurements of the exclusive neutral pion electroproduction cross section off protons at large values of $x_B$ (0.36, 0.48 and 0.60) and $Q^2$ (3.1 to 8.4 GeV$^2$) obtained from Jefferson Lab Hall A experiment E12-06-014. The corresponding structure functions $d\sigma_L/dt+\epsilon d\sigma_T/dt$, $d\sigma_{TT}/dt$, $d\sigma_{LT}/dt$ and $d\sigma_{LT'}/dt$ are extracted as a function of the proton momentum transfer $t-t_{min}$. The results suggest the amplitude for transversely polarized virtual photons continues to dominate the cross-section throughout this kinematic range. The data are well described by calculations based on transversity Generalized Parton Distributions coupled to a helicity flip Distribution Amplitude of the pion, thus providing a unique way to probe the structure of the nucleon.


title: On Connectivity-Aware Distributed Mobility Models for Area Coverage in Drone Networks

abstract: Drone networks are becoming increasingly popular in recent years and they are being used in many applications such as area coverage, delivery systems, military operations, etc. Area coverage is a broad family of applications where a group of connected drones collaboratively visit the whole or parts of an area to fulfill a specific objective and is widely being researched. Accordingly, different mobility models have been designed to define the rules of movements of the participating drones. However, most of them do not consider the network connectivity which is crucial, plus many models lack the priorities and optimization strategies that are important for drone networks. Therefore within this study, three known connectivity-aware mobility models have been analyzed comparatively. Two non-connectivity-aware mobility models have further been implemented to catch the placebo effect if any. Per the detailed experiments on the mobility models, coverage rates, connectivity levels, and message traffic have been evaluated. The study shows that the Distributed Pheromone Repel (DPR) model provides a decent coverage performance, while the Connectivity-based model and the Connected Coverage model provide better connectivity and communication quality.


title: Effects of thermal inversion induced air pollution on COVID-19

abstract: Air pollution is a threat to human health, in particular since it aggravates respiratory diseases. Early COVID-19 outbreaks in Wuhan, China and Lombardy, Italy coincided with high levels of air pollution drawing attention to a potential role of particulate matter and other pollutants in infections and more severe outcomes of the new lung disease. Both air pollution and COVID-19 outcomes are driven by human mobility and economic activity leading to spurious correlations in regression estimates. We use district-level panel data from Belgium, Brazil, Germany, Italy, the UK, and the US to estimate the impact of daily variation in air pollution levels on COVID-19 infections and deaths. Using random variation in air pollution generated by thermal inversions, we rule out that changes in mobility and economic activity are driving the results. We find that a 1%-increase in air pollution levels over the three preceding weeks leads to a 1.5% increase in weekly cases. A 1%-increase in air pollution over four weeks leads to 5.1% more COVID-19 deaths. These results indicate that short-term measures to reduce air pollution can help mitigate the health damages of the virus.


title: Deep Learning in EEG: Advance of the Last Ten-Year Critical Period

abstract: Deep learning has achieved excellent performance in a wide range of domains, especially in speech recognition and computer vision. Relatively less work has been done for EEG, but there is still significant progress attained in the last decade. Due to the lack of a comprehensive and topic widely covered survey for deep learning in EEG, we attempt to summarize recent progress to provide an overview, as well as perspectives for future developments. We first briefly mention the artifacts removal for EEG signal and then introduce deep learning models that have been utilized in EEG processing and classification. Subsequently, the applications of deep learning in EEG are reviewed by categorizing them into groups such as brain-computer interface, disease detection, and emotion recognition. They are followed by the discussion, in which the pros and cons of deep learning are presented and future directions and challenges for deep learning in EEG are proposed. We hope that this paper could serve as a summary of past work for deep learning in EEG and the beginning of further developments and achievements of EEG studies based on deep learning.


title: Making mean-estimation more efficient using an MCMC trace variance approach: DynaMITE

abstract: We introduce a novel statistical measure for MCMC-mean estimation, the inter-trace variance ${\rm trv}^{(\tau_{rel})}({\cal M},f)$, which depends on a Markov chain ${\cal M}$ and a function $f:S\to [a,b]$. The inter-trace variance can be efficiently estimated from observed data and leads to a more efficient MCMC-mean estimator. Prior MCMC mean-estimators receive, as input, upper-bounds on $\tau_{mix}$ or $\tau_{rel}$, and often also the stationary variance, and their performance is highly dependent to the sharpness of these bounds. In contrast, we introduce DynaMITE, which dynamically adjusts the sample size, it is less sensitive to the looseness of input upper-bounds on $\tau_{rel}$, and requires no bound on $v_{\pi}$. Receiving only an upper-bound ${\cal T}_{rel}$ on $\tau_{rel}$, DynaMITE estimates the mean of $f$ in $\tilde{\cal{O}}\bigl(\smash{\frac{{\cal T}_{rel} R}{\varepsilon}}+\frac{\tau_{rel}\cdot {\rm trv}^{(\tau{{rel}})}}{\varepsilon^{2}}\bigr)$ steps, without a priori bounds on the stationary variance $v_{\pi}$ or the inter-trace variance ${\rm trv}^{(\tau rel)}$. Thus we depend minimally on the tightness of ${\cal T}_{mix}$, as the complexity is dominated by $\tau_{rel}\rm{trv}^{(\tau{rel})}$ as $\varepsilon \to 0$. Note that bounding $\tau_{\rm rel}$ is known to be prohibitively difficult, however, DynaMITE is able to reduce its principal dependence on ${\cal T}_{rel}$ to $\tau_{rel}$, simply by exploiting properties of the inter-trace variance. To compare our method to known variance-aware bounds, we show ${\rm trv}^{(\tau{rel})}({\cal M},f) \leq v_{\pi}$. Furthermore, we show when $f$'s image is distributed (semi)symmetrically on ${\cal M}$'s traces, we have ${\rm trv}^{({\tau{rel}})}({\cal M},f)=o(v_{\pi}(f))$, thus DynaMITE outperforms prior methods in these cases.


title: Ride-hailing Impacts on Transit Ridership: Chicago Case Study

abstract: Existing literature on the relationship between ride-hailing (RH) and transit services is limited to empirical studies that lack real-time spatial contexts. To fill this gap, we took a novel real-time geospatial analysis approach. With source data on ride-hailing trips in Chicago, Illinois, we computed real-time transit-equivalent trips for all 7,949,902 ride-hailing trips in June 2019; the sheer size of our sample is incomparable to the samples studied in existing literature. An existing Multinomial Nested Logit Model was used to determine the probability of a ride-hailer selecting a transit alternative to serve the specific O-D pair, P(Transit|CTA). We find that 31% of ride-hailing trips are replaceable, whereas 61% of trips are not replaceable. The remaining 8% lie within a buffer zone. We measured the robustness of this probability using a parametric sensitivity analysis and performed a two-tailed t-test. Our results indicate that of the four sensitivity parameters, the probability was most sensitive to the total travel time of a transit trip. The main contribution of our research is our thorough approach and fine-tuned series of real-time spatiotemporal analyses that investigate the replaceability of ride-hailing trips for public transit. The results and discussion intend to provide perspective derived from real trips and we anticipate that this paper will demonstrate the research benefits associated with the recording and release of ride-hailing data.


title: Charged Higgs Observability via Charged Higgs Pair Production at Future Lepton Collider

abstract: The observability of charged Higgs $ H^{\pm} $ has been investigated at future lepton collider by assuming type-I 2HDM, at a centre of mass energy $ \sqrt{s}=1.5$ TeV. The signal process chain is $ e^{+}e^{-} \rightarrow Z^{*}/ \gamma^{*}\rightarrow H^{+} H^{-}\rightarrow H W^{+} H W^{-}\rightarrow b \overline{b} jjb \overline{b} jj$. The process proceed through virtual gamma and Z-boson exchange in S-channel. Several benchmark points are selected and events are analyzed to reconstruct the mass of charged Higgs bosons $ H^{\pm} $. The value of $\tan\beta $ is kept relatively high to enhance the branching ratio of $ H\rightarrow b\overline{b} $ to benefit the signal processes. The main SM background processes produced is $ t\overline{t} $. Signal selection and significance efficiencies are calculated at integrated luminosities of $100 fb^{-1} $ and $ 500 fb^{-1} $. The reconstructed and corrected mass of charged Higgs bosons $ H^{\pm} $ is determined. Through analyzing the results, it is demonstrated that charged Higgs bosons can be discovered through pair production process via its bosonic decays. This study is supposed to provide the experimentalists a good way to examine the Higgs bosons beyond SM, as well as to check the validity of 2HDM models in considered parameter space.


title: Tripling the 99 % mode conversion bandwidth of gratings using a single phase shift

abstract: Achieving large bandwidth at strong mode conversion strengths is a central problem in few-mode optical waveguides. The most uncomplicated modification of a standard mode-converting grating is to incorporate a single extra gap. Such simple design has not been applied yet for broadening conversion spectrum between core modes of any few-mode waveguide e.g. few-mode fiber, possibly due to lacking knowledge of critical parameter control: (a) phase shift value at resonance wavelength, and (b) deterministic method for placement of the phase shift(s) accommodating the change in waveguide properties during grating writing. From semi-analytical considerations, we show in full generality that our design of a single-phase-shifted grating increases the bandwidth at 99 % conversion strength by more than 2.9 times, together with elucidating crucial parameter tolerances. The phase-shift placement is determinable on-the-go from the conversion spectra, thereby obviating any estimation of modified waveguide properties during grating-writing. This offers the longed on-demand in-fiber solution for dispersion control in fiber laser cavities and communication systems.


title: Ab initio study of oxygen segregation in silicon grain boundaries: the role of strain and vacancies

abstract: Multi-crystalline silicon is widely used for producing low-cost and high-efficiency solar cells. During crystal growth and device fabrication, silicon solar cells contain grain boundaries (GBs) which are preferential segregation sites for atomic impurities such as oxygen atoms. GBs can induce charge carriers recombination significantly reducing carrier lifetimes and therefore they can be detrimental for Si device performance. We studied the correlation between structural, energetic and electronic properties of {\Sigma}3{111} Si GB in the presence of vacancies, strain and multiple O segregation. The study of the structural and energetic properties of GBs in the presence of strain and vacancies gives an accurate description of the complex mechanisms that control the segregation of oxygen atoms. We analysed tensile and compressive strain and we obtained that local tensile strain around O impurities is very effective for segregation. We also studied the role of multiple O impurities in the presence of Si vacancies finding that the segregation is favorite for those structures which have restored tetrahedral covalent bonds. The presence of vacancies attract atomic impurities in order to restore the electronic stability: the interstitial impurity becomes substitutional. This analysis was the starting point to correlate the change of the electronic properties in {\Sigma}3{111}Si GBs with O impurities in the presence of strain and vacancies. For each structure we analysed the density of states and its projection on atoms and states, the band gaps, the segregation energy and their correlation in order to characterise the nature of new energy levels. Actually, knowing the origin of defined electronic states would allow the optimization of materials in order to reduce non-radiative electron-hole recombination avoiding charge and energy losses and therefore improving solar cell efficiency.


title: Differentiable Computational Geometry for 2D and 3D machine learning

abstract: With the growth of machine learning algorithms with geometry primitives, a high-efficiency library with differentiable geometric operators are desired. We present an optimized Differentiable Geometry Algorithm Library (DGAL) loaded with implementations of differentiable operators for geometric primitives like lines and polygons. The library is a header-only templated C++ library with GPU support. We discuss the internal design of the library and benchmark its performance on some tasks with other implementations.


title: On Fields of dimension one that are Galois extensions of a global or local field

abstract: Let $K$ be a global or local field, $E/K$ a Galois extension, and Br$(E)$ the Brauer group of $E$. This paper shows that if $K$ is a local field, $v$ is its natural discrete valuation, $v'$ is the valuation of $E$ extending $v$, and $q$ is the characteristic of the residue field $\widehat E$ of $(E, v')$, then Br$(E) = \{0\}$ if and only if the following conditions hold: $\widehat E$ contains as a subfield the maximal $p$-extension of $\widehat K$, for each prime $p \neq q$; $\widehat E$ is an algebraically closed field in case the value group $v'(E)$ is $q$-indivisible. When $K$ is a global field, it characterizes the fields $E$ with Br$(E) = \{0\}$, which lie in the class of tame abelian extensions of $K$. We also give a criterion that, in the latter case, for any integer $n \ge 2$, there exists an $n$-variate $E$-form of degree $n$, which violates the Hasse principle.


title: Predictive process mining by network of classifiers and clusterers: the PEDF model

abstract: In this research, a model is proposed to learn from event log and predict future events of a system. The proposed PEDF model learns based on events' sequences, durations, and extra features. The PEDF model is built by a network made of standard clusterers and classifiers, and it has high flexibility to update the model iteratively. The model requires to extract two sets of data from log files i.e., transition differences, and cumulative features. The model has one layer of memory which means that each transition is dependent on both the current event and the previous event. To evaluate the performance of the proposed model, it is compared to the Recurrent Neural Network and Sequential Prediction models, and it outperforms them. Since there is missing performance measure for event log prediction models, three measures are proposed.


title: Combined Effects of Homogenization and Singular Perturbations : A Bloch Wave Approach

abstract: We study Bloch wave homogenization of periodically heterogeneous media with fourth order singular perturbations. We recover different homogenization regimes depending on the relative strength of the singular perturbation and length scale of the periodic heterogeneity. The homogenized tensor is obtained in terms of the first Bloch eigenvalue. The higher Bloch modes do not contribute to the homogenization limit.


title: MAC for Machine Type Communications in Industrial IoT -- Part I: Protocol Design and Analysis

abstract: In this two-part paper, we propose a novel medium access control (MAC) protocol for machine-type communications in the industrial internet of things. The considered use case features a limited geographical area and a massive number of devices with sporadic data traffic and different priority types. We target at supporting the devices while satisfying their quality of service (QoS) requirements with a single access point and a single channel, which necessitates a customized design that can significantly improve the MAC performance. In Part I of this paper, we present the MAC protocol that comprises a new slot structure, corresponding channel access procedure, and mechanisms for supporting high device density and providing differentiated QoS. A key idea behind this protocol is sensing-based distributed coordination for significantly improving channel utilization. To characterize the proposed protocol, we analyze its delay performance based on the packet arrival rates of devices. The analytical results provide insights and lay the groundwork for the fine-grained scheduling with QoS guarantee as presented in Part II.


title: MAC for Machine Type Communications in Industrial IoT -- Part II: Scheduling and Numerical Results

abstract: In the second part of this paper, we develop a centralized packet transmission scheduling scheme to pair with the protocol designed in Part I and complete our medium access control (MAC) design for machine-type communications in the industrial internet of things. For the networking scenario, fine-grained scheduling that attends to each device becomes necessary, given stringent quality of service (QoS) requirements and diversified service types, but prohibitively complex for a large number of devices. To address this challenge, we propose a scheduling solution in two steps. First, we develop algorithms for device assignment based on the analytical results from Part I, when parameters of the proposed protocol are given. Then, we train a deep neural network for assisting in the determination of the protocol parameters. The two-step approach ensures the accuracy and granularity necessary for satisfying the QoS requirements and avoids excessive complexity from handling a large number of devices. Integrating the distributed coordination in the protocol design from Part I and the centralized scheduling from this part, the proposed MAC protocol achieves high performance, demonstrated through extensive simulations. For example, the results show that the proposed MAC can support 1000 devices under an aggregated traffic load of 3000 packets per second with a single channel and achieve <0.5ms average delay and <1% average collision probability among 50 high priority devices.


title: Weighted lens depth: Some applications to supervised classification

abstract: Starting with Tukey's pioneering work in the 1970's, the notion of depth in statistics has been widely extended especially in the last decade. These extensions include high dimensional data, functional data, and manifold-valued data. In particular, in the learning paradigm, the depth-depth method has become a useful technique. In this paper we extend the notion of lens depth to the case of data in metric spaces, and prove its main properties, with particular emphasis on the case of Riemannian manifolds, where we extend the concept of lens depth in such a way that it takes into account non-convex structures on the data distribution. Next we illustrate our results with some simulation results and also in some interesting real datasets, including pattern recognition in phylogenetic trees using the depth--depth approach.


title: Vanishing relaxation time dynamics of the Jordan Moore-Gibson-Thompson equation arising in nonlinear acoustics

abstract: The (third-order in time) JMGT equation \cite{Jordan2,HCP} is a nonlinear (quasi-linear) Partial Differential Equation (PDE) model introduced to describe a nonlinear propagation of sound in an acoustic medium. The important feature is that the model avoids the infinite speed of propagation paradox associated with a classical second order in time equation referred to as Westervelt equation. Replacing Fourier's law by Maxwell-Cattaneo's law gives rise to the third order in time derivative scaled by a small parameter $\tau >0$, the latter represents the thermal relaxation time parameter and is intrinsic to the medium where the dynamics occur. In this paper we provide an asymptotic analysis of the third order model when $\tau \rightarrow 0 $. It is shown that the corresponding solutions converge {\it in a strong topology of the phase space } to a limit which is the solution of Westervelt equation. In addition, rate of convergence is provided for solutions displaying higher order regularity. This addresses an open question raised in \cite{kaltev2}, where a related JMGT equation has been studied and {\it weak star } convergence of the solutions when $\tau \rightarrow 0$ has been established. Thus, our main contribution is showing {\it strong convergence on infinite time horizon,} along with related rates of convergence valid on a finite time horizon. The key to unlocking the difficulty owns to a tight control and propagation of the "smallness" of the initial data in carrying the estimates at three different topological levels. The rate of convergence allows one then to estimate the relaxation time needed for the signal to reach the target. The interest in studying this type of problems is motivated by a large array of applications arising in engineering and medical sciences.


title: Spectral shift via "lateral" perturbation

abstract: We consider a compact perturbation $H_0 = S + K_0^* K_0$ of a self-adjoint operator $S$ with an eigenvalue $\lambda^\circ$ below its essential spectrum and the corresponding eigenfunction $f$. The perturbation is assumed to be "along" the eigenfunction $f$, namely $K_0f=0$. The eigenvalue $\lambda^\circ$ belongs to the spectra of both $H_0$ and $S$. Let $S$ have $\sigma$ more eigenvalues below $\lambda^\circ$ than $H_0$; $\sigma$ is known as the spectral shift at $\lambda^\circ$. We now allow the perturbation to vary in a suitable operator space and study the continuation of the eigenvalue $\lambda^\circ$ in the spectrum of $H(K)=S + K^* K$. We show that the eigenvalue as a function of $K$ has a critical point at $K=K_0$ and the Morse index of this critical point is the spectral shift $\sigma$. A version of this theorem also holds for some non-positive perturbations.


title: Optical frequency generation using fiber Bragg grating filters for applications in portable quantum sensing

abstract: A method for the agile generation of the optical frequencies required for laser cooling and atom interferometry of rubidium is demonstrated. It relies on fiber Bragg grating technology to filter the output of an electro-optic modulator and was demonstrated in a robust, alignment-free, single-seed, frequency-doubled, telecom fiber laser system. The system was capable of 50 ns frequency switching over a ~40 GHz range, ~0.5 W output power and amplitude modulation with a ~15 ns rise/fall time and an extinction ratio of 120 $\pm$ 2 dB. The technology is ideal for enabling high-bandwidth, mobile industrial and space applications of quantum technologies.


title: Searching and Sorting with O(n^2) processors in O(1) time

abstract: The proliferation of number of processing elements (PEs) in parallel computer systems, along with the use of more extensive parallelization of algorithms causes the interprocessor communications dominate VLSI chip space. This paper proposes a new architecture to overcome this issue by using simple crosspoint switches to pair PEs instead of a complex interconnection network. Based on the cyclic permutation wiring idea described in \cite{oruc2016self}, this pairing leads to a linear crosspoint array of $n(n-1)/2$ processing elements and as many crosspoints. We demonstrate the versatility of this new parallel architecture by designing fast searching and sorting algorithms for it. In particular, we show that finding a minimum, maximum, and searching a list of $n$ elements can all be performed in $O(1)$ time with elementary logic gates with $O(n)$ fan-in, and in $O(\lg n)$ time with $O(1)$ fan-in. We further show that sorting a list of $n$ elements can also be carried out in $O(1)$ time using elementary logic gates with $O(n)$ fan-in and threshold logic gates. The sorting time increases to $O(\lg n\lg\lg n)$ if only elementary logic gates with $O(1)$ fan-in are used. The algorithm can find the maximum among $n$ elements in $O(1)$ time, and sort $n$ elements in $O(\lg n (\lg\lg n))$ time. In addition, we show how other fundamental queries can be handled within the same order of time complexities.


title: On the energy budget of the transition of a neutron star into the third family branch

abstract: Transition of a compact star into the third family for an equation of state (EoS) featuring mass twins is considered. The energy released at a baryon number conserving transition for static compact stars configurations is computed for two sets of models for comparison. The EoS of choice is the density dependent functional DD2 EoS with excluded model correction for hadronic matter which suffers a phase transition into deconfined quark matter described by a constant speed of sound approach. The two sets of EoS models feature different compact star mass onsets that maximize the energy and radius difference at the transition while simultaneously fulfilling state-of-the-art constraints from multi-messenger astronomy and empirical nuclear data. It is found that the maximal energy budget at the transitions falls in the range of $10^{49}$ - $10^{52}$ ergs.


title: Level sets of depth measures in abstract spaces

abstract: The lens depth of a point has been recently extended to general metric spaces, which is not the case for most depths. It is defined as the probability of being included in the intersection of two random balls centred at two random points X and Y, with the same radius d(X, Y). We study the consistency in Hausdorff and measure distance, of the level sets of the empirical lens depth, based on an iid sample on a general metric space. We also prove that the boundary of the empirical level sets are consistent estimators of their population counterparts, and analyze two real-life examples


title: Random matrices and controllability of dynamical systems

abstract: We introduce the concept of $\epsilon$-uncontrollability for random linear systems, i.e. linear system in which the usual matrices have been replaced by random matrices. We also estimate the $\epsilon$-uncontrollability in the case where the matrices come from the Gaussian orthogonal ensemble. Our proof utilizes tools from systems theory, probability theory and convex geometry.


title: Asymmetric spin wave dispersion due to a saturation magnetization gradient

abstract: We demonstrate using micromagnetic simulations and a theoretical model that a gradient in the saturation magnetization ($M_s$) of a perpendicularly magnetized ferromagnetic film induces a non-reciprocal spin wave propagation and, consequently an asymmetric dispersion relation. The $M_s$ gradient adds a linear potential to the spin wave equation of motion consistent with the presence of a force. We consider a transformation from an inertial reference frame in which the $M_s$ is constant to an accelerated reference frame where the resulting inertial force corresponds to the force from the $M_s$ gradient. As in the Doppler effect, the frequency shift leads to an asymmetric dispersion relation. Additionally, we show that under certain circumstances, unidirectional propagation of spin waves can be achieved which is essential for the design of magnonic circuits. Our results become more relevant in light of recent experimental works in which a suitable thermal landscape is used to dynamically modulate the saturation magnetization.


title: Convergence of energy forms on Sierpinski gaskets with added rotated triangle

abstract: We study the convergence of resistance metrics and resistance forms on a converging sequence of spaces. As an application, we study the existence and uniqueness of self-similar Dirichlet forms on Sierpinski gaskets with added rotated triangles. The fractals depend on a parameter in a continuous way. When the parameter is irrational, the fractal is not post critically finite (p.c.f.), and there are infinitely many ways that two cells intersect. In this case, we will define the Dirichlet form as a limit in some $\Gamma$-convergence sense of the Dirichlet forms on p.c.f. fractals that approximate it.


title: On the Convergence of Continuous Constrained Optimization for Structure Learning

abstract: Recently, structure learning of directed acyclic graphs (DAGs) has been formulated as a continuous optimization problem by leveraging an algebraic characterization of acyclicity. The constrained problem is solved using the augmented Lagrangian method (ALM) which is often preferred to the quadratic penalty method (QPM) by virtue of its standard convergence result that does not require the penalty coefficient to go to infinity, hence avoiding ill-conditioning. However, the convergence properties of these methods for structure learning, including whether they are guaranteed to return a DAG solution, remain unclear, which might limit their practical applications. In this work, we examine the convergence of ALM and QPM for structure learning in the linear, nonlinear, and confounded cases. We show that the standard convergence result of ALM does not hold in these settings, and demonstrate empirically that its behavior is akin to that of the QPM which is prone to ill-conditioning. We further establish the convergence guarantee of QPM to a DAG solution, under mild conditions. Lastly, we connect our theoretical results with existing approaches to help resolve the convergence issue, and verify our findings in light of an empirical comparison of them.


title: LaHAR: Latent Human Activity Recognition using LDA

abstract: Processing sequential multi-sensor data becomes important in many tasks due to the dramatic increase in the availability of sensors that can acquire sequential data over time. Human Activity Recognition (HAR) is one of the fields which are actively benefiting from this availability. Unlike most of the approaches addressing HAR by considering predefined activity classes, this paper proposes a novel approach to discover the latent HAR patterns in sequential data. To this end, we employed Latent Dirichlet Allocation (LDA), which is initially a topic modelling approach used in text analysis. To make the data suitable for LDA, we extract the so-called "sensory words" from the sequential data. We carried out experiments on a challenging HAR dataset, demonstrating that LDA is capable of uncovering underlying structures in sequential data, which provide a human-understandable representation of the data. The extrinsic evaluations reveal that LDA is capable of accurately clustering HAR data sequences compared to the labelled activities.


title: On the Overlooked Pitfalls of Weight Decay and How to Mitigate Them: A Gradient-Norm Perspective

abstract: Weight decay is a simple yet powerful regularization technique that has been very widely used in training of deep neural networks (DNNs). While weight decay has attracted much attention, previous studies fail to discover some overlooked pitfalls on large gradient norms resulted by weight decay. In this paper, we discover that, weight decay can unfortunately lead to large gradient norms at the final phase (or the terminated solution) of training, which often indicates bad convergence and poor generalization. To mitigate the gradient-norm-centered pitfalls, we present the first practical scheduler for weight decay, called the Scheduled Weight Decay (SWD) method that can dynamically adjust the weight decay strength according to the gradient norm and significantly penalize large gradient norms during training. Our experiments also support that SWD indeed mitigates large gradient norms and often significantly outperforms the conventional constant weight decay strategy for Adaptive Moment Estimation (Adam).


title: Investigation of superspreading COVID-19 outbreaks events in meat and poultry processing plants in Germany: A cross-sectional study

abstract: Since May 2020, several COVID-19 outbreaks have occurred in the German meat industry despite various protective measures, and temperature and ventilation conditions were considered as possible high-risk factors. This cross-sectional study examined meat and poultry plants to examine possible risk factors. Companies completed a self-administered questionnaire on the work environment and protective measures taken to prevent SARS-CoV-2 infection. Multivariable logistic regression analysis adjusted for the possibility to distance at least 1.5 meters, break rules, and employment status was performed to identify risk factors associated with COVID-19 cases. Twenty-two meat and poultry plants with 19,072 employees participated. The prevalence of COVID-19 in the seven plants with more than 10 cases was 12.1% and was highest in the deboning and meat cutting area with 16.1%. A subsample analysis where information on maximal ventilation rate per employee was available revealed an effect for ventilation rate (adjusted odds ratio (AOR) 0.996, 95% CI 0.993-0.999). When including temperature as an interaction term in the working area, the effect of the ventilation rate did not change. Increasing room temperatures resulted in a lower chance of obtaining a positive COVID-19 test result (AOR 0.90 95% CI 0.82-0.99), and a 0.1% greater chance of a positive COVID-19 test for the interaction term (AOR 1.001, 95% CI 1.000-1.003). Our results further indicate that climate conditions and low outdoor air flow are factors that can promote the spread of SARS-CoV-2 aerosols. A possible requirement for pandemic mitigation strategies in industrial workplace settings is to increase the ventilation rate.


title: Pseudospectra and Simultaneous Power Control

abstract: We prove that, for $M>0$ and $n,m\ge2$, we can find two $10\times10$ matrices $A$ and $B$ with identical pseudospectra such that we have simultaneously $\|A^n\|/\|B^n\|>M$ and $\|A^m\|/\|B^m\|>M$. We also prove that, under certain conditions, this result holds for two more general functions of $A$ and $B$.


title: Imbalance Robust Softmax for Deep Embeeding Learning

abstract: Deep embedding learning is expected to learn a metric space in which features have smaller maximal intra-class distance than minimal inter-class distance. In recent years, one research focus is to solve the open-set problem by discriminative deep embedding learning in the field of face recognition (FR) and person re-identification (re-ID). Apart from open-set problem, we find that imbalanced training data is another main factor causing the performance degradation of FR and re-ID, and data imbalance widely exists in the real applications. However, very little research explores why and how data imbalance influences the performance of FR and re-ID with softmax or its variants. In this work, we deeply investigate data imbalance in the perspective of neural network optimisation and feature distribution about softmax. We find one main reason of performance degradation caused by data imbalance is that the weights (from the penultimate fully-connected layer) are far from their class centers in feature space. Based on this investigation, we propose a unified framework, Imbalance-Robust Softmax (IR-Softmax), which can simultaneously solve the open-set problem and reduce the influence of data imbalance. IR-Softmax can generalise to any softmax and its variants (which are discriminative for open-set problem) by directly setting the weights as their class centers, naturally solving the data imbalance problem. In this work, we explicitly re-formulate two discriminative softmax (A-Softmax and AM-Softmax) under the framework of IR-Softmax. We conduct extensive experiments on FR databases (LFW, MegaFace) and re-ID database (Market-1501, Duke), and IR-Softmax outperforms many state-of-the-art methods.


title: Better Aggregation in Test-Time Augmentation

abstract: Test-time augmentation -- the aggregation of predictions across transformed versions of a test input -- is a common practice in image classification. Traditionally, predictions are combined using a simple average. In this paper, we present 1) experimental analyses that shed light on cases in which the simple average is suboptimal and 2) a method to address these shortcomings. A key finding is that even when test-time augmentation produces a net improvement in accuracy, it can change many correct predictions into incorrect predictions. We delve into when and why test-time augmentation changes a prediction from being correct to incorrect and vice versa. Building on these insights, we present a learning-based method for aggregating test-time augmentations. Experiments across a diverse set of models, datasets, and augmentations show that our method delivers consistent improvements over existing approaches.


title: Semi-compactness and multiple oscillating pulses in kink scattering

abstract: In this work we consider model of asymmetric kinks, where the behavior of the solution in one side is different from the other side. Also, the models depend of an integer $n$ and, with the increase of $n$, the constructed kink assumes a hybrid character: a compactlike profile on one side and a kinklike profile on the other side. We investigate numerically the kink-antikink and antikink-kink dynamics, with the aim to understand the effect of the transition of the usual kink to the semi-compacton structure. The kink-antikink process shows the formation of one-bounce windows for small values of $n$. The increase of $n$ favors the breaking this structure and the appearance of oscillatory modes. For antikink-kink collisions we report the appearance of two-bounce windows for small values of the parameter. We also found an intricate structure of two-oscillation windows.


title: Nonsingular Parameterization for Modeling Translational Motion Using Euler Parameters

abstract: A parameterization is described for quantifying translational motion of a point in three-dimensional Euclidean space. The parameterization is similar to well-known parameterizations such as spherical coordinates in that both position and velocity are decoupled into magnitude and orientation components. Unlike these standard parameterizations, where principal rotation sequences are employed, the method presented in this research employs Euler parameters. By using Euler parameters instead of Euler angles, singularities and trigonometric functions are removed from the equations of motion. The parameterization is demonstrated on two examples, where it is found that the new parameterization offers both mathematical and computational advantages over other commonly used parameterizations.


title: Crystallization Mechanism Tuned Phase-Change Materials: Quantum Effect on Te-Terminated Boundary

abstract: While phase-change materials (PCMs) composed of chalcogenide have different crystallization mechanisms (CM), such as nucleation-dominated Ge2Sb2Te5 (GST) and growth-dominated GeTe (GT), revealing the essential reason of CM as well as the tuned properties is still a long-standing issue. Here, we remarkably find the distinct stability of Te-terminated (111) boundaries (TTB) in different systems, which provides a path to understand the difference in CM. It stems from the quantum effect of molecular orbital theory: the optimal local chemical composition results in the formation of TTB without dangling bonds (DB) in GST but with DB in GT, where DB destabilizes boundary due to its distorted local environment mismatching Oh symmetry of p orbitals. Moreover, the inner vacancy concentration in GST is alterable and controlled by TTB, manifested by the absence of cubic-to-hexagonal transition in carbon-doped GST of small grains and minimized inner vacancy. Finally, the charge transport property (CTP) is controlled by boundary via changing the density of charge or hole nearby as well as vacancy. These findings open the door to tune CTP by CM, which is necessary for achieving low-power and ultrafast devices.


title: LINDT: Tackling Negative Federated Learning with Local Adaptation

abstract: Federated Learning (FL) is a promising distributed learning paradigm, which allows a number of data owners (also called clients) to collaboratively learn a shared model without disclosing each client's data. However, FL may fail to proceed properly, amid a state that we call negative federated learning (NFL). This paper addresses the problem of negative federated learning. We formulate a rigorous definition of NFL and analyze its essential cause. We propose a novel framework called LINDT for tackling NFL in run-time. The framework can potentially work with any neural-network-based FL systems for NFL detection and recovery. Specifically, we introduce a metric for detecting NFL from the server. On occasion of NFL recovery, the framework makes adaptation to the federated model on each client's local data by learning a Layer-wise Intertwined Dual-model. Experiment results show that the proposed approach can significantly improve the performance of FL on local data in various scenarios of NFL.


title: Poisson-bracket formulation of the dynamics of fluids of deformable particles

abstract: Using the Poisson bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes. We explicitly consider the example of a dense biological tissue as described by the Vertex model energy, where cell shape has been proposed as a structural order parameter for a liquid-solid transition. The hydrodynamic model of biological tissue proposed here captures the coupling of cell shape to flow, and provides a starting point for modeling the rheology of dense tissue.


title: Fast Decentralized Linear Functions Over Edge Fluctuating Graphs

abstract: Implementing linear transformations is a key task in the decentralized signal processing framework, which performs learning tasks on data sets distributed over multi-node networks. That kind of network can be represented by a graph. Recently, some decentralized methods have been proposed to compute linear transformations by leveraging the notion of graph shift operator, which captures the local structure of the graph. However, existing approaches have some drawbacks such as considering some special instances of linear transformations, or reducing the family of transformations by assuming that a shift matrix is given such that a subset of its eigenvectors spans the subspace of interest. In contrast, this paper develops a framework for computing a wide class of linear transformations in a decentralized fashion by relying on the notion of graph shift operator. The main goal of the proposed method is to compute the desired linear transformation in a small number of iterations. To this end, a set of successive graph shift operators is employed, then, a new optimization problem is proposed whose goal is to compute the desired transformation as fast as possible. In addition, usually, the topology of the networks, especially the wireless sensor networks, change randomly because of node failures or random links. In this paper, the effect of edge fluctuations on the performance of the proposed method is studied. To deal with the negative effect of edge fluctuations, an online kernel-based method is proposed which enables nodes to estimate the missed values with their at hand information. The proposed method can also be employed to sparsify the network graph or reduce the number of local exchanges between nodes, which saves sensors power in the wireless sensor networks.


title: Heat-blanketed convection and its implications for the continental lithosphere

abstract: Earth's continents are characterized by a strong enrichment in long-lived radioactive isotopes. Recent estimates suggest that they contribute to 33\% of the heat released at the surface of the Earth, while occupying less than 1\% of the mantle. This distinctive feature has profound implications for the underlying mantle by impacting its thermal structure and heat transfer. However, the effects of a continental crust enriched in heat-producing elements on the underlying mantle have not yet been systematically investigated. Here, we conduct a preliminary investigation by considering a simplified convective system consisting in a mixed heated fluid where all the internal heating is concentrated in a top layer of thickness $d_{HL}$ (referred to as "heat-blanketed convection"). We perform 24 numerical simulations in 3D Cartesian geometry for four specific set-ups and various values of $d_{HL}$. Our results suggest that the effects of the heated layer strongly depend on its thickness relative to the thickness of the thermal boundary layer ($\delta_{TBL}$) in the homogeneous heating case ($d_{HL} = 1.0$). More specifically, for $d_{HL} > \delta_{TBL}$, the effects induced by the heated layer are quite modest, while, for $d_{HL} < \delta_{TBL}$, the properties of the convective system are strongly altered as $d_{HL}$ decreases. In particular, the surface heat flux and convective vigour are significantly enhanced for very thin heated layers compared to the case $d_{HL} = 1.0$. The vertical distribution of heat producing elements may therefore play a key role on mantle dynamics. For Earth, the presence of continents should however not affect significantly the surface heat flux, and thus the Earth's cooling rate.


title: Learnable Boundary Guided Adversarial Training

abstract: Previous adversarial training raises model robustness under the compromise of accuracy on natural data. In this paper, we reduce natural accuracy degradation. We use the model logits from one clean model to guide learning of another one robust model, taking into consideration that logits from the well trained clean model embed the most discriminative features of natural data, {\it e.g.}, generalizable classifier boundary. Our solution is to constrain logits from the robust model that takes adversarial examples as input and makes it similar to those from the clean model fed with corresponding natural data. It lets the robust model inherit the classifier boundary of the clean model. Moreover, we observe such boundary guidance can not only preserve high natural accuracy but also benefit model robustness, which gives new insights and facilitates progress for the adversarial community. Finally, extensive experiments on CIFAR-10, CIFAR-100, and Tiny ImageNet testify to the effectiveness of our method. We achieve new state-of-the-art robustness on CIFAR-100 without additional real or synthetic data with auto-attack benchmark \footnote{\url{https://github.com/fra31/auto-attack}}. Our code is available at \url{https://github.com/dvlab-research/LBGAT}.


title: (3+1)D topological orders with only a $\mathbb{Z}_2$-charged particle

abstract: There is exactly one bosonic (3+1)-dimensional topological order whose only nontrivial particle is an emergent boson: pure $\mathbb{Z}_2$ gauge theory. There are exactly two (3+1)-dimensional topological orders whose only nontrivial particle is an emergent fermion: pure "spin-$\mathbb{Z}_2$" gauge theory, in which the dynamical field is a spin structure; and an anomalous version thereof. I give three proofs of this classification, varying from hands-on to abstract. Along the way, I provide a detailed study of the braided fusion $2$-category $\mathcal{Z}_{(1)}(\Sigma \mathbf{SVec})$ of string and particle operators in pure spin-$\mathbb{Z}_2$ gauge theory.


title: Quantum Phase Space Description of a Cosmological Minimal Massive Bigravity Model

abstract: Bimetric gravity theories describes gravitational interactions in the presence of an extra spin-2 field. The Hassan-Rosen (HR) nonlinear massive minimal bigravity theory is a ghost-free bimetric theory formulated with respect a flat, dynamical reference metric. In this work the deformation quantization formalism is applied to a HR cosmological model in the minisuperspace. The quantization procedure is performed explicitly for quantum cosmology in the minisuperspace. The Friedmann-Lema\^itre-Robertson-Walker (FLRW) model with flat metrics is worked out and the computation of the Wigner functions for theHartle-Hawking, Vilenkin and Linde wavefunctions are done numerically and, in the Hartle-Hawking case, also analytically. From the stability analysis in the quantum minisuper phase space it is found an interpretation of the mass of graviton as an emergent cosmological constant and as a measure of the deviation of classical behavior of the Wigner functions. Also, from the Hartle-Hawking case, an interesting relation between the curvature and the mass of graviton in a cusp catastrophe surface is discussed


title: The Selectivity and Competition of the Mind's Eye in Visual Perception

abstract: Research has shown that neurons within the brain are selective to certain stimuli. For example, the fusiform face area (FFA) region is known by neuroscientists to selectively activate when people see faces over non-face objects. However, the mechanisms by which the primary visual system directs information to the correct higher levels of the brain are currently unknown. In our work, we mimic several high-level neural mechanisms of perception by creating a novel computational model that incorporates lateral and top down feedback in the form of hierarchical competition. Not only do we show that these elements can help explain the information flow and selectivity of high level areas within the brain, we also demonstrate that these neural mechanisms provide the foundation of a novel classification framework that rivals traditional supervised learning in computer vision. Additionally, we present both quantitative and qualitative results that demonstrate that our generative framework is consistent with neurological themes and enables simple, yet robust category level classification.


title: Iterated collapsing phenomenon on $G_2$-manifolds

abstract: We propose a new collapsing mechanism for $G_2$-metrics, with the generic region admitting a circle bundle structure over a K3 fibration over a Riemann surface. The adiabatic description involves a weighted version of the maximal submanifold equation. In a local smooth setting we prove the existence of formal power series solutions, and the problem of compactification is discussed at a heuristic level. The $Spin(7)$ analogue is also discussed.


title: Properties of Galaxies in Cosmic Filaments around the Virgo Cluster

abstract: We present the properties of galaxies in filaments around the Virgo cluster with respect to their vertical distance from the filament spine using the NASA-Sloan Atlas catalog. The filaments are mainly composed of low-mass, blue dwarf galaxies. We observe that the g - r color of galaxies becomes blue and stellar mass decreases with increasing vertical filament distance. The galaxies were divided into higher-mass (log h^2M_* > 8) and lower-mass (log h^2M_* < 8) subsamples. We also examine the g - r color, stellar mass, Halpha equivalent width (EW(Halpha)), near-ultraviolet (NUV) - r color, and HI fraction distributions of the two subsamples against the vertical distance. The lower-mass galaxies exhibit a negative g - r color gradient, whereas higher-mass galaxies have a flat g - r color distribution. We observe a negative EW(Halpha) gradient for higher-mass galaxies, whereas lower-mass galaxies show no distinct EW(Halpha) variation. In contrast, the NUV - r color distribution of higher-mass galaxies shows no strong trend, whereas the lower-mass galaxies show a negative NUV - r color gradient. We do not witness clear gradients of HI fraction in either the higher- or lower-mass subsamples. We propose that the negative color and stellar mass gradients of galaxies can be explained by mass assembly from past galaxy mergers at different vertical filament distances. In addition, galaxy interactions might be responsible for the contrasting features of EW(Halpha) and NUV - r color distributions between the higher- and lower-mass subsamples. The HI fraction distributions of the two subsamples suggest that ram-pressure stripping and gas accretion could be ignorable processes in the Virgo filaments.


title: Quantifying model uncertainty for the observed non-Gaussian data by the Hellinger distance

abstract: Mathematical models for complex systems under random fluctuations often certain uncertain parameters. However, quantifying model uncertainty for a stochastic differential equation with an $\alpha$-stable L\'evy process is still lacking. Here, we propose an approach to infer all the uncertain non-Gaussian parameters and other system parameters by minimizing the Hellinger distance over the parameter space. The Hellinger distance measures the similarity between an empirical probability density of non-Gaussian observations and a solution (as a probability density) of the associated nonlocal Fokker-Planck equation. Numerical experiments verify that our method is feasible for estimating single and multiple parameters. Meanwhile, we find an optimal estimation interval of the estimated parameters. This method is beneficial for extracting governing dynamical system models under non-Gaussian fluctuations, as in the study of abrupt climate changes in the Dansgaard-Oeschger events.


title: Quantum tricriticality of chiral-coherent phase in quantum Rabi triangle

abstract: The interplay of interactions, symmetries and gauge fields usually leads to intriguing quantum many-body phases. To explore the nature of emerging phases, we study a quantum Rabi triangle system as an elementary building block for synthesizing an artificial magnetic field. We develop an analytical approach to study the rich phase diagram and the associated quantum criticality. Of particular interest is the emergence of a chiral-coherent phase, which breaks both the $\mathbb{Z}_2$ and the chiral symmetry. In this chiral phase, photons flow unidirectionally and the chirality can be tuned by the artificial gauge field, exhibiting a signature of broken time-reversal symmetry. The finite-frequency scaling analysis further confirms the associated phase transition to be in the universality class of the Dicke model. This model can simulate a broad range of physical phenomena of light-matter coupling systems, and may have an application in future developments of various quantum information technologies.


title: Application of the mixing length theory to assess the generation of melt in internally heated systems

abstract: The effect of melting in planetary mantles plays a key role in their thermo-chemical evolution. Because of the laterally heterogeneous nature of melting, 3D numerical simulations are in principle necessary prohibiting us from exploring wide ranges of conditions. To overcome this issue, we propose a new analytical framework allowing to estimate the amount and depths of melting in a 1D analytical model for a simplified convective system. To do so, we develop an approach, partly based on an extended version of the mixing length theory, able to estimate the distribution of the hottest temperatures in natural systems. The approach involves several free parameters that are calibrated by fitting 3D numerical simulations. We demonstrate that our algorithm provides the melting profile at steady-state as well as the long-term evolution in fairly good agreement with the ones obtained in 3D numerical simulations. We then apply our framework for a wide variety of planetary sizes and heating rates. We find that an increase in planetary size increases the depth of melting for small planets but decreases for large planets. This change in the trend is caused by the pressure dependence of the solidus.


title: Stochastic optimization with decision-dependent distributions

abstract: Stochastic optimization problems often involve data distributions that change in reaction to the decision variables. This is the case for example when members of the population respond to a deployed classifier by manipulating their features so as to improve the likelihood of being positively labeled. Recent works on performative prediction have identified an intriguing solution concept for such problems: find the decision that is optimal with respect to the static distribution that the decision induces. Continuing this line of work, we show that typical stochastic algorithms -- originally designed for static problems -- can be applied directly for finding such equilibria with little loss in efficiency. The reason is simple to explain: the main consequence of the distributional shift is that it corrupts algorithms with a bias that decays linearly with the distance to the solution. Using this perspective, we obtain sharp convergence guarantees for popular algorithms, such as stochastic gradient, clipped gradient, proximal point, and dual averaging methods, along with their accelerated and proximal variants. In realistic applications, deployment of a decision rule is often much more expensive than sampling. We show how to modify the aforementioned algorithms so as to maintain their sample efficiency while performing only logarithmically many deployments.


title: Acoustic Cavitation Rheometry

abstract: Characterization of soft materials is challenging due to their high compliance and the strain-rate dependence of their mechanical properties. The inertial microcavitation-based high strain-rate rheometry (IMR) method [Estrada et al., J. Mech. Phys. Solids, 2018, 112, 291-317] combines laser-induced cavitation measurements with a model for the bubble dynamics to measure local properties of polyacrylamide hydrogel under high strain-rates from $10^3$ to $10^8$ s$^{-1}$. While promising, laser-induced cavitation involves plasma formation and optical breakdown during nucleation, a process that could alter local material properties before measurements are obtained. In the present study, we extend the IMR method to another means to generate cavitation, namely high-amplitude focused ultrasound, and apply the resulting acoustic-cavitation-based IMR to characterize the mechanical properties of agarose hydrogels. Material properties including viscosity, elastic constants, and a stress-free bubble radius are inferred from bubble radius histories in $0.3\%$ and $1\%$ agarose gels. An ensemble-based data assimilation is used to further help interpret the obtained estimates. The resulting parameter distributions are consistent with available measurements of agarose gel properties and with expected trends related to gel concentration and high strain-rate loading. Our findings demonstrate the utility of applying IMR and data assimilation methods with single-bubble acoustic cavitation data for measurement of viscoelastic properties.


title: Renormalons beyond the Borel plane

abstract: The renormalon singularities are a known source of the divergent behavior of asymptotic perturbative series from field theoretical models. These singularities live in the Borel plane and are responsible for ambiguities in the physical plane. We propose that field theories can have renormalons beyond the usual first Borel plane. We show an example with a scalar field theory where, considering a chain of cat's eyes diagrams, the model develops a Gevrey-3 asymptotic series.


title: Distance-based Data Cleaning: A Survey (Technical Report)

abstract: With the rapid development of the internet technology, dirty data are commonly observed in various real scenarios, e.g., owing to unreliable sensor reading, transmission and collection from heterogeneous sources. To deal with their negative effects on downstream applications, data cleaning approaches are designed to preprocess the dirty data before conducting applications. The idea of most data cleaning methods is to identify or correct dirty data, referring to the values of their neighbors which share the same information. Unfortunately, owing to data sparsity and heterogeneity, the number of neighbors based on equality relationship is rather limited, especially in the presence of data values with variances. To tackle this problem, distance-based data cleaning approaches propose to consider similarity neighbors based on value distance. By tolerance of small variants, the enriched similarity neighbors can be identified and used for data cleaning tasks. At the same time, distance relationship between tuples is also helpful to guide the data cleaning, which contains more information and includes the equality relationship. Therefore, distance-based technology plays an important role in the data cleaning area, and we also have reason to believe that distance-based data cleaning technology will attract more attention in data preprocessing research in the future. Hence this survey provides a classification of four main data cleaning tasks, i.e., rule profiling, error detection, data repair and data imputation, and comprehensively reviews the state of the art for each class.


title: Gonogo: An R Implementation of Test Methods to Perform, Analyze and Simulate Sensitivity Experiments

abstract: This work provides documentation for a suite of R functions contained in gonogo.R. The functions provide sensitivity testing practitioners and researchers with an ability to conduct, analyze and simulate various sensitivity experiments involving binary responses and a single stimulus level (e.g., drug dosage, drop height, velocity, etc.). Included are the modern Neyer and 3pod adaptive procedures, as well as the Bruceton and Langlie. The latter two benchmark procedures are capable of being performed according to generalized up-down transformed-response rules. Each procedure is designated phase-one of a three-phase experiment. The goal of phase-one is to achieve overlapping data. The two additional (and optional) refinement phases utilize the D-optimal criteria and the Robbins-Monro-Joseph procedure. The goals of the two refinement phases are to situate testing in the vicinity of the median and tails of the latent response distribution, respectively.


title: Bayesian Nonparametric Estimation for Point Processes with Spatial Homogeneity: A Spatial Analysis of NBA Shot Locations

abstract: Basketball shot location data provide valuable summary information regarding players to coaches, sports analysts, fans, statisticians, as well as players themselves. Represented by spatial points, such data are naturally analyzed with spatial point process models. We present a novel nonparametric Bayesian method for learning the underlying intensity surface built upon a combination of Dirichlet process and Markov random field. Our method has the advantage of effectively encouraging local spatial homogeneity when estimating a globally heterogeneous intensity surface. Posterior inferences are performed with an efficient Markov chain Monte Carlo (MCMC) algorithm. Simulation studies show that the inferences are accurate and that the method is superior compared to the competing methods. Application to the shot location data of $20$ representative NBA players in the 2017-2018 regular season offers interesting insights about the shooting patterns of these players. A comparison against the competing method shows that the proposed method can effectively incorporate spatial contiguity into the estimation of intensity surfaces.


title: A difference method for solving the nonlinear $q$-factional differential equations on time scale

abstract: The $q$-fractional differential equation usually describe the physics process imposed on the time scale set $T_q$. In this paper, we first propose a difference formula for discretizing the fractional $q$-derivative $^cD_q^\alpha x(t)$ on the time scale set $T_q$ with order $0<\alpha<1$ and scale index $0<q<1$. We establish a rigours truncation error boundness and prove that this difference formula is unconditionally stable. Then, we consider the difference method for solving the initial problem of $q$-fractional differential equation: $^cD_q^\alpha x(t)=f(t,x(t))$ on the time scale set. We prove the unique existence and stability of the difference solution and give the convergence analysis. Numerical experiments show the effectiveness and high accuracy of the proposed difference method.


title: Time-Symmetry Breaking in Hamiltonian Mechanics. III. A Memoir for Douglas James Henderson [1934-2020]

abstract: Following Berni Alder [1] and Francis Ree [2], Douglas Henderson was the third of Bill's California coworkers from the 1960s to die in 2020. Motivated by Doug's death we undertook better to understand Lyapunov instability and the breaking of time symmetry in continuum and atomistic simulations. Here we have chosen to extend our explorations of an interesting pair of nonequilibrium systems, the steady shockwave and the unsteady rarefaction wave. We eliminate the need for boundary potentials by simulating the collisions of pairs of mirror-images projectiles. The resulting shock and rarefaction structures are respectively the results of the compression and the expansion of simple fluids. Shockwaves resulting from compression have a steady structure while the rarefaction fans resulting from free expansions continually broaden. We model these processes using classical molecular dynamics and Eulerian fluid mechanics in two dimensions. Although molecular dynamics is time-reversible the reversed simulation of a steady shockwave compression soon results in an unsteady rarefaction fan, violating the microscopic time symmetry of the motion equations but in agreement with the predictions of macroscopic Navier-Stokes fluid mechanics. The explanations for these results are an interesting combination of two (irreversible) instabilities, Lyapunov and Navier-Stokes.


title: On InstaHide, Phase Retrieval, and Sparse Matrix Factorization

abstract: In this work, we examine the security of InstaHide, a scheme recently proposed by [Huang, Song, Li and Arora, ICML'20] for preserving the security of private datasets in the context of distributed learning. To generate a synthetic training example to be shared among the distributed learners, InstaHide takes a convex combination of private feature vectors and randomly flips the sign of each entry of the resulting vector with probability 1/2. A salient question is whether this scheme is secure in any provable sense, perhaps under a plausible hardness assumption and assuming the distributions generating the public and private data satisfy certain properties. We show that the answer to this appears to be quite subtle and closely related to the average-case complexity of a new multi-task, missing-data version of the classic problem of phase retrieval. Motivated by this connection, we design a provable algorithm that can recover private vectors using only the public vectors and synthetic vectors generated by InstaHide, under the assumption that the private and public vectors are isotropic Gaussian.


title: Crystalline cohomology over general bases

abstract: Building on ideas of Berthelot, we develop a crystalline cohomology formalism over divided power rings $(A, I_0, \eta)$ for any ring $A$, allowing $\mathbf{Z}$-flat $A$. For a smooth $A$-scheme $Y$ and a closed subscheme $X$ of $Y$ for which $\eta$ extends to $I_0 \mathscr{O}_X$, a (quasi-coherent) crystal $\mathscr{F}$ on $(X/A)_{\rm{cris}}$ is equivalent to a specific type of module with integrable $A$-linear connection over a certain completion $D_{Y,\eta}(X)^{\wedge}$ (called "pd-adic") of the divided power envelope $D_{Y,\eta}(X)$ of $Y$ along $X$ (with divided power structure $\delta$) Our main result, building on ideas of Bhatt and de Jong for $\mathbf{Z}/(p^e)$-schemes (where pd-adic completion has no effect), is a natural isomorphism between ${\rm{R}}\Gamma((X/A)_{\rm{cris}}, \mathscr{F})$ and the Zariski hypercohomology of the pd-adically completed de Rham complex $\mathscr{F} \widehat{\otimes} \widehat{\Omega}^*_{D_{Y,\eta}(X)^{\wedge}/A,\delta}$ arising from the module with integrable connection over $D_{Y,\eta}(X)^{\wedge}$ associated to $\mathscr{F}$. By a variant of the same methods, we obtain a representative of the complex $\mathscr{F} \widehat{\otimes} \widehat{\Omega}^*_{D_{Y,\eta}(X)^{\wedge}/A,\delta}$ in the derived category of sheaves of $A$-modules on $X$ in terms of a \v{C}ech-Alexander construction. When $\mathscr{F}=\mathscr{O}_{X/A}$, our comparison theorem implies that in the derived category of sheaves of $A$-modules on $X$, the pd-adic completion of $\Omega^*_{D_{Y,\eta}(X)/A,\delta}$ functorially depends only on $X$. Over $\mathbf{Q}$-algebras $A$, so pd-adic completion becomes ideal-adic completion, this recovers a result of Hartshorne.


title: CoMatch: Semi-supervised Learning with Contrastive Graph Regularization

abstract: Semi-supervised learning has been an effective paradigm for leveraging unlabeled data to reduce the reliance on labeled data. We propose CoMatch, a new semi-supervised learning method that unifies dominant approaches and addresses their limitations. CoMatch jointly learns two representations of the training data, their class probabilities and low-dimensional embeddings. The two representations interact with each other to jointly evolve. The embeddings impose a smoothness constraint on the class probabilities to improve the pseudo-labels, whereas the pseudo-labels regularize the structure of the embeddings through graph-based contrastive learning. CoMatch achieves state-of-the-art performance on multiple datasets. It achieves substantial accuracy improvements on the label-scarce CIFAR-10 and STL-10. On ImageNet with 1% labels, CoMatch achieves a top-1 accuracy of 66.0%, outperforming FixMatch by 12.6%. Furthermore, CoMatch achieves better representation learning performance on downstream tasks, outperforming both supervised learning and self-supervised learning. Code and pre-trained models are available at https://github.com/salesforce/CoMatch.