title: Practical Massively Parallel Monte-Carlo Tree Search Applied to Molecular Design

abstract: It is common practice to use large computational resources to train neural networks, as is known from many examples, such as reinforcement learning applications. However, while massively parallel computing is often used for training models, it is rarely used for searching solutions for combinatorial optimization problems. In this paper, we propose a novel massively parallel Monte-Carlo Tree Search (MP-MCTS) algorithm that works efficiently for 1,000 worker scale, and apply it to molecular design. This is the first work that applies distributed MCTS to a real-world and non-game problem. Existing work on large-scale parallel MCTS show efficient scalability in terms of the number of rollouts up to 100 workers, but suffer from the degradation in the quality of the solutions. MP-MCTS maintains the search quality at larger scale, and by running MP-MCTS on 256 CPU cores for only 10 minutes, we obtained candidate molecules having similar score to non-parallel MCTS running for 42 hours. Moreover, our results based on parallel MCTS (combined with a simple RNN model) significantly outperforms existing state-of-the-art work. Our method is generic and is expected to speed up other applications of MCTS.


title: Investment Disputes and Abnormal Volatility of Stocks

abstract: Dramatic growth of investment disputes between foreign investors and host states rises serious questions about the impact of those disputes on investors. This paper is the first to explain increased uncertainty of investors about the outcome of arbitration, which may or may not lead to compensation for damages claimed by the investor. We find robust evidence that investment disputes lead to abnormal share fluctuations of companies involved in disputes with host countries. Importantly, while a positive outcome for an investor decreases uncertainty back to original levels, we document strong increase in the volatility of companies with negative outcome for the investor. We find that several variables including size of the award, political instability, location of arbitration, country of origin of investor or public policy considerations in host country explain large portion of the investor's uncertainty.


title: Monte Carlo Simulation Techniques

abstract: Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection method, and Markov chain Monte Carlo to sample a probability distribution function, and methods for variance reduction to evaluate numerical integrals using the Monte Carlo simulation. We will also briefly introduce the quasi-Monte Carlo sampling at the end of this lecture.


title: Invariant half-spaces for rank-one perturbations

abstract: If T is a bounded linear operator acting on an infinite-dimensional Banach space, then there exists and operator F of rank at most one and arbitrarily small norm such that T-F has an invariant subspace of infinite dimension and codimension. This improves results of Tcaciuc and other authors.


title: Interpretability in PRA

abstract: In this paper from 2009 we study IL(PRA), the interpretability logic of PRA. As PRA is neither an essentially reflexive theory nor finitely axiomatizable, the two known arithmetical completeness results do not apply to PRA: IL(PRA) is not ILM or ILP. IL(PRA) does of course contain all the principles known to be part of IL(All), the interpretability logic of the principles common to all reasonable arithmetical theories. In this paper, we take two arithmetical properties of PRA and see what their consequences in the modal logic IL(PRA) are. These properties are reflected in the so-called Beklemishev Principle B, and Zambella's Principle Z, neither of which is a part of IL(All). Both principles and their interrelation are submitted to a modal study. In particular, we prove a frame condition for B. Moreover, we prove that Z follows from a restricted form of B. Finally, we give an overview of the known relationships of IL(PRA) to important other interpetability principles.


title: Structure and Design of HoloGen

abstract: Increasing popularity of augmented and mixed reality systems has seen a similar increase of interest in 2D and 3D computer generated holography (CGH). Unlike stereoscopic approaches, CGH can fully represent a light field including depth of focus, accommodation and vergence. Along with existing telecommunications, imaging, projection, lithography, beam shaping and optical tweezing applications, CGH is an exciting technique applicable to a wide array of photonic problems including full 3D representation. Traditionally, the primary roadblock to acceptance has been the significant numerical processing required to generate holograms requiring both significant expertise and significant computational power. This article discusses the structure and design of HoloGen. HoloGen is an MIT licensed application that may be used to generate holograms using a wide array of algorithms without expert guidance. HoloGen uses a Cuda C and C++ backend with a C# and Windows Presentation Framework graphical user interface. The article begins by introducing HoloGen before providing an in-depth discussion of its design and structure. Particular focus is given to the communication, data transfer and algorithmic aspects.


title: Base sizes for primitive groups with soluble stabilisers

abstract: Let $G$ be a finite primitive permutation group on a set $\Omega$ with point stabiliser $H$. Recall that a subset of $\Omega$ is a base for $G$ if its pointwise stabiliser is trivial. We define the base size of $G$, denoted $b(G,H)$, to be the minimal size of a base for $G$. Determining the base size of a group is a fundamental problem in permutation group theory, with a long history stretching back to the 19th century. Here one of our main motivations is a theorem of Seress from 1996, which states that $b(G,H) \leqslant 4$ if $G$ is soluble. In this paper we extend Seress' result by proving that $b(G,H) \leqslant 5$ for all finite primitive groups $G$ with a soluble point stabiliser $H$. This bound is best possible. We also determine the exact base size for all almost simple groups and we study random bases in this setting. For example, we prove that the probability that $4$ random elements in $\Omega$ form a base tends to $1$ as $|G|$ tends to infinity.


title: Contrastive learning of global and local features for medical image segmentation with limited annotations

abstract: A key requirement for the success of supervised deep learning is a large labeled dataset - a condition that is difficult to meet in medical image analysis. Self-supervised learning (SSL) can help in this regard by providing a strategy to pre-train a neural network with unlabeled data, followed by fine-tuning for a downstream task with limited annotations. Contrastive learning, a particular variant of SSL, is a powerful technique for learning image-level representations. In this work, we propose strategies for extending the contrastive learning framework for segmentation of volumetric medical images in the semi-supervised setting with limited annotations, by leveraging domain-specific and problem-specific cues. Specifically, we propose (1) novel contrasting strategies that leverage structural similarity across volumetric medical images (domain-specific cue) and (2) a local version of the contrastive loss to learn distinctive representations of local regions that are useful for per-pixel segmentation (problem-specific cue). We carry out an extensive evaluation on three Magnetic Resonance Imaging (MRI) datasets. In the limited annotation setting, the proposed method yields substantial improvements compared to other self-supervision and semi-supervised learning techniques. When combined with a simple data augmentation technique, the proposed method reaches within 8% of benchmark performance using only two labeled MRI volumes for training, corresponding to only 4% (for ACDC) of the training data used to train the benchmark. The code is made public at https://github.com/krishnabits001/domain_specific_cl.


title: Covariant reduction of classical Hamiltonian Field Theories: From D'Alembert to Klein-Gordon and Schr\"odinger

abstract: A novel reduction procedure for covariant classical field theories, reflecting the generalized symplectic reduction theory of Hamiltonian systems, is presented. The departure point of this reduction procedure consists in the choice of a submanifold of the manifold of solutions of the equations describing a field theory. Then, the covariance of the geometrical objects involved, will allow to define equations of motion on a reduced space. The computation of the canonical geometrical structure is performed neatly by using the geometrical framework provided by the multisymplectic description of covariant field theories. The procedure is illustrated by reducing the D'Alembert theory on a five-dimensional Minkowski space-time to a massive Klein-Gordon theory in four dimensions and, more interestingly, to the Schr\"odinger equation in 3+1 dimensions.


title: VLA/Realfast Detection of a Burst from FRB180916.J0158+65 and Tests for Periodic Activity

abstract: We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs.


title: $T\bar{T}$ deformation of chiral bosons and Chern-Simons AdS$_3$ gravity

abstract: We study the $T\bar{T}$ deformation of the chiral bosons and show the equivalence between the chiral bosons of opposite chiralities and the scalar fields at the Hamiltonian level under the deformation. We also derive the deformed Lagrangian of more generic theories which contain an arbitrary number of chiral bosons to all orders. By using these results, we derive the $T\bar{T}$ deformed boundary action of the AdS$_3$ gravity theory in the Chern-Simons formulation. We compute the deformed one-loop torus partition function, which satisfies the $T\bar{T}$ flow equation up to the one-loop order. Finally, we calculate the deformed stress tensor of a solution describing a BTZ black hole in the boundary theory, which coincides with the boundary stress tensor derived from the BTZ black hole with a finite cutoff.


title: Analytic Behaviour of Wave Function

abstract: We obtained the order $\rho$ and the type $\sigma$ of wave function with all kind of power potentials, and found that the order and the type are compatible with the condition $|\sigma|\rho =\sqrt{v_m}$ except to $m=-2$. At the same time, we ansatz that the wave function satisfaction relation $\psi(r)=f(r)\exp[g(r)]$, where $f(r),g(r)$ are polynomials, and the power of $g(r)$ is no more than the order $\rho$. Using this ansatz we can solve the Schr\"{o}dinger equation, but unfortunately sometimes we need the parameters of potential to satisfy certain relations.


title: Self-Attention Enhanced Patient Journey Understanding in Healthcare System

abstract: Understanding patients' journeys in healthcare system is a fundamental prepositive task for a broad range of AI-based healthcare applications. This task aims to learn an informative representation that can comprehensively encode hidden dependencies among medical events and its inner entities, and then the use of encoding outputs can greatly benefit the downstream application-driven tasks. A patient journey is a sequence of electronic health records (EHRs) over time that is organized at multiple levels: patient, visits and medical codes. The key challenge of patient journey understanding is to design an effective encoding mechanism which can properly tackle the aforementioned multi-level structured patient journey data with temporal sequential visits and a set of medical codes. This paper proposes a novel self-attention mechanism that can simultaneously capture the contextual and temporal relationships hidden in patient journeys. A multi-level self-attention network (MusaNet) is specifically designed to learn the representations of patient journeys that is used to be a long sequence of activities. The MusaNet is trained in end-to-end manner using the training data derived from EHRs. We evaluated the efficacy of our method on two medical application tasks with real-world benchmark datasets. The results have demonstrated the proposed MusaNet produces higher-quality representations than state-of-the-art baseline methods. The source code is available in https://github.com/xueping/MusaNet.


title: Privacy-Preserving Technology to Help Millions of People: Federated Prediction Model for Stroke Prevention

abstract: Prevention of stroke with its associated risk factors has been one of the public health priorities worldwide. Emerging artificial intelligence technology is being increasingly adopted to predict stroke. Because of privacy concerns, patient data are stored in distributed electronic health record (EHR) databases, voluminous clinical datasets, which prevent patient data from being aggregated and restrains AI technology to boost the accuracy of stroke prediction with centralized training data. In this work, our scientists and engineers propose a privacy-preserving scheme to predict the risk of stroke and deploy our federated prediction model on cloud servers. Our system of federated prediction model asynchronously supports any number of client connections and arbitrary local gradient iterations in each communication round. It adopts federated averaging during the model training process, without patient data being taken out of the hospitals during the whole process of model training and forecasting. With the privacy-preserving mechanism, our federated prediction model trains over all the healthcare data from hospitals in a certain city without actual data sharing among them. Therefore, it is not only secure but also more accurate than any single prediction model that trains over the data only from one single hospital. Especially for small hospitals with few confirmed stroke cases, our federated model boosts model performance by 10%~20% in several machine learning metrics. To help stroke experts comprehend the advantage of our prediction system more intuitively, we developed a mobile app that collects the key information of patients' statistics and demonstrates performance comparisons between the federated prediction model and the single prediction model during the federated training process.


title: Improving Post Training Neural Quantization: Layer-wise Calibration and Integer Programming

abstract: Lately, post-training quantization methods have gained considerable attention, as they are simple to use, and require only a small unlabeled calibration set. This small dataset cannot be used to fine-tune the model without significant over-fitting. Instead, these methods only use the calibration set to set the activations' dynamic ranges. However, such methods always resulted in significant accuracy degradation, when used below 8-bits (except on small datasets). Here we aim to break the 8-bit barrier. To this end, we minimize the quantization errors of each layer separately by optimizing its parameters over the calibration set. We empirically demonstrate that this approach is: (1) much less susceptible to over-fitting than the standard fine-tuning approaches, and can be used even on a very small calibration set; and (2) more powerful than previous methods, which only set the activations' dynamic ranges. Furthermore, we demonstrate how to optimally allocate the bit-widths for each layer, while constraining accuracy degradation or model compression by proposing a novel integer programming formulation. Finally, we suggest model global statistics tuning, to correct biases introduced during quantization. Together, these methods yield state-of-the-art results for both vision and text models. For instance, on ResNet50, we obtain less than 1\% accuracy degradation --- with 4-bit weights and activations in all layers, but the smallest two. We open-sourced our code.


title: Symmetric contact systems of segments, pseudotriangulations and inductive constructions for corresponding surface graphs

abstract: We characterise the quotient surface graphs arising from symmetric contact systems of line segments in the plane and also from symmetric pointed pseudotriangulations in the case where the group of symmetries is generated by a translation or a rotation of finite order. These results generalise well known results of Thomassen, in the case of line segments, and of Streinu and Haas et al., in the case of pseudotriangulations. Our main tool is a new inductive characterisation of the appropriate classes of surface graphs. We also discuss some consequences of our results in the area of geometric rigidity theory.


title: Multi-view Low-rank Preserving Embedding: A Novel Method for Multi-view Representation

abstract: In recent years, we have witnessed a surge of interest in multi-view representation learning, which is concerned with the problem of learning representations of multi-view data. When facing multiple views that are highly related but sightly different from each other, most of existing multi-view methods might fail to fully integrate multi-view information. Besides, correlations between features from multiple views always vary seriously, which makes multi-view representation challenging. Therefore, how to learn appropriate embedding from multi-view information is still an open problem but challenging. To handle this issue, this paper proposes a novel multi-view learning method, named Multi-view Low-rank Preserving Embedding (MvLPE). It integrates different views into one centroid view by minimizing the disagreement term, based on distance or similarity matrix among instances, between the centroid view and each view meanwhile maintaining low-rank reconstruction relations among samples for each view, which could make more full use of compatible and complementary information from multi-view features. Unlike existing methods with additive parameters, the proposed method could automatically allocate a suitable weight for each view in multi-view information fusion. However, MvLPE couldn't be directly solved, which makes the proposed MvLPE difficult to obtain an analytic solution. To this end, we approximate this solution based on stationary hypothesis and normalization post-processing to efficiently obtain the optimal solution. Furthermore, an iterative alternating strategy is provided to solve this multi-view representation problem. The experiments on six benchmark datasets demonstrate that the proposed method outperforms its counterparts while achieving very competitive performance.


title: LSTM-TrajGAN: A Deep Learning Approach to Trajectory Privacy Protection

abstract: The prevalence of location-based services contributes to the explosive growth of individual-level trajectory data and raises public concerns about privacy issues. In this research, we propose a novel LSTM-TrajGAN approach, which is an end-to-end deep learning model to generate privacy-preserving synthetic trajectory data for data sharing and publication. We design a loss metric function TrajLoss to measure the trajectory similarity losses for model training and optimization. The model is evaluated on the trajectory-user-linking task on a real-world semantic trajectory dataset. Compared with other common geomasking methods, our model can better prevent users from being re-identified, and it also preserves essential spatial, temporal, and thematic characteristics of the real trajectory data. The model better balances the effectiveness of trajectory privacy protection and the utility for spatial and temporal analyses, which offers new insights into the GeoAI-powered privacy protection.


title: PIDA: Smooth and Stable Flight Using Stochastic Dual Simplex Algorithm and Genetic Filter

abstract: This paper presents a new Proportional-Integral-Derivative-Accelerated (PIDA) control with a derivative filter to improve quadcopter flight stability in a noisy environment. The mathematical model is derived from having an accurate model with a high level of fidelity by addressing the problems of non-linearity, uncertainties, and coupling. These uncertainties and measurement noises cause instability in flight and automatic hovering. The proposed controller associated with a heuristic Genetic Filter (GF) addresses these challenges. The tuning of the proposed PIDA controller associated with the objective of controlling is performed by Stochastic Dual Simplex Algorithm (SDSA). GF is applied to the PIDA control to estimate the observed states and parameters of quadcopters in both attitude and altitude. The simulation results show that the proposed control associated with GF has a strong ability to track the desired point in the presence of disturbances.


title: Optimization of Fuzzy Controller of a Wind Power Plant Based on the Swarm Intelligence

abstract: The article considers the problem of the optimal control of a wind power plant based on fuzzy control and automation of generating the fuzzy rule base. Fuzzy rules by experts do not always provide a maximum power output of the wind plant and fuzzy rule bases require an adjustment in the case of changing the parameters of the wind power plant or the environment. This research proposes the method for optimizing the fuzzy rules base compiled by various experts. The method is based on balancing weights of fuzzy rules into the base by the Particle Swarm Optimization algorithm. The experiment has shown that the proposed method allows forming the fuzzy rule base as an exemplary optimal base from a non-optimized set of fuzzy rules. The optimal fuzzy rule base has been taken under consideration for the concrete control loop of wind power plant and the concrete fuzzy model of the wind.


title: Reinforcement Learning as Iterative and Amortised Inference

abstract: There are several ways to categorise reinforcement learning (RL) algorithms, such as either model-based or model-free, policy-based or planning-based, on-policy or off-policy, and online or offline. Broad classification schemes such as these help provide a unified perspective on disparate techniques and can contextualise and guide the development of new algorithms. In this paper, we utilise the control as inference framework to outline a novel classification scheme based on amortised and iterative inference. We demonstrate that a wide range of algorithms can be classified in this manner providing a fresh perspective and highlighting a range of existing similarities. Moreover, we show that taking this perspective allows us to identify parts of the algorithmic design space which have been relatively unexplored, suggesting new routes to innovative RL algorithms.


title: StarCraft II Build Order Optimization using Deep Reinforcement Learning and Monte-Carlo Tree Search

abstract: The real-time strategy game of StarCraft II has been posed as a challenge for reinforcement learning by Google's DeepMind. This study examines the use of an agent based on the Monte-Carlo Tree Search algorithm for optimizing the build order in StarCraft II, and discusses how its performance can be improved even further by combining it with a deep reinforcement learning neural network. The experimental results accomplished using Monte-Carlo Tree Search achieves a score similar to a novice human player by only using very limited time and computational resources, which paves the way to achieving scores comparable to those of a human expert by combining it with the use of deep reinforcement learning.


title: The effect of adaptive mobility policy to the spread of COVID-19 in urban environment: intervention analysis of Seoul, South Korea

abstract: Although severe mobility restrictions are recognized as the key enabler to contain COVID-19, there has been few scientific studies to validate such approach, especially in urban context. This study analyzes mobility pattern changes in Seoul, South Korea that adopted adaptive approach toward mobility. Intervention analyses reveal that major mobility reduction did occur two weeks before the city's case peak. Such voluntary adjustments exhibit strong preference shift toward private mode from public transit. Large reductions occurred in non-essential and high-contact activities of shopping and dining, while work and Starbucks trips were less affected. The collective evaluation reveal that major changes in epidemiology, mobility and policy occurred simultaneously, with no lagging nor leading contributors. Our study demonstrates that collective understanding the mutual aspects among mobility, epidemiology and policy is essential. Incremental and flexible mobility restriction is not only possible but necessary, especially for a pandemic of extensive spatial and temporal scales.


title: A decoding algorithm for 2D convolutional codes over the erasure channel

abstract: Two-dimensional (2D) convolutional codes are a generalization of (1D) convolutional codes, which are very appropriate for transmission over an erasure channel. In this paper, we present a decoding algorithm for 2D convolutional codes over this kind of channel by reducing the decoding process to several decoding steps with 1D convolutional codes. Moreover, we provide constructions of 2D convolutional codes that are specially taylored to our decoding algorithm.


title: The United Nations Sustainable Development Goals in Systems Engineering: Eliciting sustainability requirements

abstract: This paper discusses a PhD research project testing the hypothesis that using the United Nations Sustainable Development Goals(SDG) as explicit inputs to drive the Software Requirements Engineering process will result in requirements with improved sustainability benefits. The research has adopted the Design Science Research Method (DSRM) [21] to test a process named SDG Assessment for Requirements Elicitation (SDGARE). Three DSRM cycles are being used to test the hypothesis in safety-critical, highprecision, software-intensive systems in aerospace and healthcare. Initial results from the first two DSRM cycles support the hypothesis. However, these cycles are in a plan-driven (waterfall) development context and future research agenda would be a similar application in an Agile development context.


title: Neural Path Features and Neural Path Kernel : Understanding the role of gates in deep learning

abstract: Rectified linear unit (ReLU) activations can also be thought of as 'gates', which, either pass or stop their pre-activation input when they are 'on' (when the pre-activation input is positive) or 'off' (when the pre-activation input is negative) respectively. A deep neural network (DNN) with ReLU activations has many gates, and the on/off status of each gate changes across input examples as well as network weights. For a given input example, only a subset of gates are 'active', i.e., on, and the sub-network of weights connected to these active gates is responsible for producing the output. At randomised initialisation, the active sub-network corresponding to a given input example is random. During training, as the weights are learnt, the active sub-networks are also learnt, and potentially hold very valuable information. In this paper, we analytically characterise the role of active sub-networks in deep learning. To this end, we encode the on/off state of the gates of a given input in a novel 'neural path feature' (NPF), and the weights of the DNN are encoded in a novel 'neural path value' (NPV). Further, we show that the output of network is indeed the inner product of NPF and NPV. The main result of the paper shows that the 'neural path kernel' associated with the NPF is a fundamental quantity that characterises the information stored in the gates of a DNN. We show via experiments (on MNIST and CIFAR-10) that in standard DNNs with ReLU activations NPFs are learnt during training and such learning is key for generalisation. Furthermore, NPFs and NPVs can be learnt in two separate networks and such learning also generalises well in experiments.


title: Optimal speed in Thoroughbred horse racing

abstract: The objective of this work is to provide a mathematical analysis on how a Thoroughbred horse should regulate its speed over the course of a race to optimize performance. Because Thoroughbred horses are not capable of running the whole race at top speed, determining what pace to set and when to unleash the burst of speed is essential. Our model relies on mechanics, energetics (both aerobic and anaerobic) and motor control. It is a system of coupled ordinary differential equations on the velocity, the propulsive force and the anaerobic energy, that leads to an optimal control problem that we solve. In order to identify the parameters meaningful for Thoroughbred horses, we use velocity data on races in Chantilly (France) provided by France Galop, the French governing body of flat horse racing in France. Our numerical simulations of performance optimization then provide the optimal speed along the race, the oxygen uptake evolution in a race, as well as the energy or the propulsive force. It also predicts how the horse has to change its effort and velocity according to the topography (altitude and bending) of the track.


title: LimeOut: An Ensemble Approach To Improve Process Fairness

abstract: Artificial Intelligence and Machine Learning are becoming increasingly present in several aspects of human life, especially, those dealing with decision making. Many of these algorithmic decisions are taken without human supervision and through decision making processes that are not transparent. This raises concerns regarding the potential bias of these processes towards certain groups of society, which may entail unfair results and, possibly, violations of human rights. Dealing with such biased models is one of the major concerns to maintain the public trust. In this paper, we address the question of process or procedural fairness. More precisely, we consider the problem of making classifiers fairer by reducing their dependence on sensitive features while increasing (or, at least, maintaining) their accuracy. To achieve both, we draw inspiration from "dropout" techniques in neural based approaches, and propose a framework that relies on "feature drop-out" to tackle process fairness. We make use of "LIME Explanations" to assess a classifier's fairness and to determine the sensitive features to remove. This produces a pool of classifiers (through feature dropout) whose ensemble is shown empirically to be less dependent on sensitive features, and with improved or no impact on accuracy.


title: COVID-ABS: An Agent-Based Model of COVID-19 Epidemic to Simulate Health and Economic Effects of Social Distancing Interventions

abstract: The COVID-19 pandemic due to the SARS-CoV-2 coronavirus has directly impacted the public health and economy worldwide. To overcome this problem, countries have adopted different policies and non-pharmaceutical interventions for controlling the spread of the virus. This paper proposes the COVID-ABS, a new SEIR (Susceptible-Exposed-Infected-Recovered) agent-based model that aims to simulate the pandemic dynamics using a society of agents emulating people, business and government. Seven different scenarios of social distancing interventions were analyzed, with varying epidemiological and economic effects: (1) do nothing, (2) lockdown, (3) conditional lockdown, (4) vertical isolation, (5) partial isolation, (6) use of face masks, and (7) use of face masks together with 50% of adhesion to social isolation. In the impossibility of implementing scenarios with lockdown, which present the lowest number of deaths and highest impact on the economy, scenarios combining the use of face masks and partial isolation can be the more realistic for implementation in terms of social cooperation. The COVID-ABS model was implemented in Python programming language, with source code publicly available. The model can be easily extended to other societies by changing the input parameters, as well as allowing the creation of a multitude of other scenarios. Therefore, it is a useful tool to assist politicians and health authorities to plan their actions against the COVID-19 epidemic.


title: Endpoints for randomized controlled clinical trials for COVID-19 treatments

abstract: Introduction: Endpoint choice for randomized controlled trials of treatments for COVID-19 is complex. A new disease brings many uncertainties, but trials must start rapidly. COVID-19 is heterogeneous, ranging from mild disease that improves within days to critical disease that can last weeks and can end in death. While improvement in mortality would provide unquestionable evidence about clinical significance of a treatment, sample sizes for a study evaluating mortality are large and may be impractical. Furthermore, patient states in between "cure" and "death" represent meaningful distinctions. Clinical severity scores have been proposed as an alternative. However, the appropriate summary measure for severity scores has been the subject of debate, particularly in relating to the uncertainty about the time-course of COVID-19. Outcomes measured at fixed time-points may risk missing the time of clinical benefit. An endpoint such as time-to-improvement (or recovery), avoids the timing problem. However, some have argued that power losses will result from reducing the ordinal scale to a binary state of "recovered" vs "not recovered." Methods: We evaluate statistical power for possible trial endpoints for COVID-19 treatment trials using simulation models and data from two recent COVID-19 treatment trials. Results: Power for fixed-time point methods depends heavily on the time selected for evaluation. Time-to-improvement (or recovery) analyses do not specify a time-point. Time-to-event approaches have reasonable statistical power, even when compared to a fixed time-point method evaluated at the optimal time. Discussion: Time-to-event analyses methods have advantages in the COVID-19 setting, unless the optimal time for evaluating treatment effect is known in advance. Even when the optimal time is known, a time-to-event approach may increase power for interim analyses.


title: Analytic Continuation and Reciprocity Relation for Collinear Splitting in QCD

abstract: It is well-known that direct analytic continuation of DGLAP evolution kernel (splitting functions) from space-like to time-like kinematics breaks down at three loops. We identify the origin of this breakdown as splitting functions are not analytic function of external momenta. However, splitting functions can be constructed from square of (generalized) splitting amplitudes. We establish the rule of analytic continuation for splitting amplitudes, and use them to determine the analytic continuation of certain holomorphic and anti-holomorphic part of splitting functions and transverse-momentum dependent distributions. In this way we derive the time-like splitting functions at three loops without ambiguity. We also propose a reciprocity relation for singlet splitting functions, and provide non-trivial evidence that it holds in QCD at least through three loops.


title: Regime transitions in stratified shear flows: the link between horizontal and inclined ducts

abstract: We present the analytical solution for the two-dimensional velocity and density fields within an approximation for laminar stratified inclined duct (SID) flows where diffusion dominates over inertia in the along-channel momentum equation but it is negligible in the density transport equation. We refer to this approximation as the hydrostatic/gravitational/viscous in momentum and advective in density (HGV-A) approximation due to the leading balances in the governing equations. The analytical solution is valid for laminar flows in a two-layer configuration in the limit of long ducts. Under such conditions, the non-dimensional volume flux is given by the Froude number $Fr^* =Re_g/(A\,K)$ with $Re_g$ the gravitational Reynolds number, $A$ the aspect ratio of the duct, and $K$ a geometrical parameter that depends on the tilt of the duct and is obtained from the analytical solution. The analytical solution in the HGV-A approximation is validated against results from laboratory experiments, and allows us to gain new insight into the dynamics and properties of SID flows. Most importantly, constant values of $Fr^*$ describe, in both horizontal and inclined ducts, the transitions between increasingly turbulent flow regimes: from laminar flow, to interfacial waves, to intermittent turbulence and sustained turbulence.


title: On the existence and the uniqueness of the solution to a fluid-structure interaction problem

abstract: In this paper we consider the linearized version of a system of partial differential equations arising from a fluid-structure interaction model. We prove the existence and the uniqueness of the solution under natural regularity assumptions.


title: Quotients of monoid extensions and their interplay with Baer sums

abstract: Cosetal extensions of monoids generalise extensions of groups, special Schreier extensions of monoids and Leech's normal extensions of groups by monoids. They share a number of properties with group extensions, including a notion of Baer sum when the kernel is abelian. However, unlike group extensions (with fixed kernel and cokernel) there may be nontrivial morphisms between them. We explore the structure of the category of cosetal extensions and relate it to an analogue of second cohomology groups. Finally, the order structure and additive structures are combined to give an indexed family of inverse semigroups of extensions. These in turn can be combined into an inverse category.


title: Subgraph Neural Networks

abstract: Deep learning methods for graphs achieve remarkable performance on many node-level and graph-level prediction tasks. However, despite the proliferation of the methods and their success, prevailing Graph Neural Networks (GNNs) neglect subgraphs, rendering subgraph prediction tasks challenging to tackle in many impactful applications. Further, subgraph prediction tasks present several unique challenges: subgraphs can have non-trivial internal topology, but also carry a notion of position and external connectivity information relative to the underlying graph in which they exist. Here, we introduce SubGNN, a subgraph neural network to learn disentangled subgraph representations. We propose a novel subgraph routing mechanism that propagates neural messages between the subgraph's components and randomly sampled anchor patches from the underlying graph, yielding highly accurate subgraph representations. SubGNN specifies three channels, each designed to capture a distinct aspect of subgraph topology, and we provide empirical evidence that the channels encode their intended properties. We design a series of new synthetic and real-world subgraph datasets. Empirical results for subgraph classification on eight datasets show that SubGNN achieves considerable performance gains, outperforming strong baseline methods, including node-level and graph-level GNNs, by 19.8% over the strongest baseline. SubGNN performs exceptionally well on challenging biomedical datasets where subgraphs have complex topology and even comprise multiple disconnected components.


title: Provability and interpretability logics with restricted realizations

abstract: The provability logic of a theory T is the set of modal formulas, which under any arithmetical realization are provable in T . We slightly modify this notion by requiring the arithmetical realizations to come from a specified set $\Gamma$. We make an analogous modification for interpretability logics. This is a paper from 2012. We first studied provability logics with restricted realizations, and show that for various natural candidates of theory T and restriction set $\Gamma$, where each sentence in $\Gamma$ has a well understood (meta)-mathematical content in T, the result is the logic of linear frames. However, for the theory Primitive Recursive Arithmetic (PRA), we define a fragment that gives rise to a more interesting provability logic, by capitalizing on the well-studied relationship between PRA and I$\Sigma_1$. We then study interpretability logics, obtaining some upper bounds for IL(PRA), whose characterization remains a major open question in interpretability logic. Again this upper bound is closely relatively to linear frames. The technique is also applied to yield the non-trivial result that IL(PRA) $\subset$ ILM.


title: Infinite attention: NNGP and NTK for deep attention networks

abstract: There is a growing amount of literature on the relationship between wide neural networks (NNs) and Gaussian processes (GPs), identifying an equivalence between the two for a variety of NN architectures. This equivalence enables, for instance, accurate approximation of the behaviour of wide Bayesian NNs without MCMC or variational approximations, or characterisation of the distribution of randomly initialised wide NNs optimised by gradient descent without ever running an optimiser. We provide a rigorous extension of these results to NNs involving attention layers, showing that unlike single-head attention, which induces non-Gaussian behaviour, multi-head attention architectures behave as GPs as the number of heads tends to infinity. We further discuss the effects of positional encodings and layer normalisation, and propose modifications of the attention mechanism which lead to improved results for both finite and infinitely wide NNs. We evaluate attention kernels empirically, leading to a moderate improvement upon the previous state-of-the-art on CIFAR-10 for GPs without trainable kernels and advanced data preprocessing. Finally, we introduce new features to the Neural Tangents library (Novak et al., 2020) allowing applications of NNGP/NTK models, with and without attention, to variable-length sequences, with an example on the IMDb reviews dataset.


title: Exact posterior distributions of wide Bayesian neural networks

abstract: Recent work has shown that the prior over functions induced by a deep Bayesian neural network (BNN) behaves as a Gaussian process (GP) as the width of all layers becomes large. However, many BNN applications are concerned with the BNN function space posterior. While some empirical evidence of the posterior convergence was provided in the original works of Neal (1996) and Matthews et al. (2018), it is limited to small datasets or architectures due to the notorious difficulty of obtaining and verifying exactness of BNN posterior approximations. We provide the missing theoretical proof that the exact BNN posterior converges (weakly) to the one induced by the GP limit of the prior. For empirical validation, we show how to generate exact samples from a finite BNN on a small dataset via rejection sampling.


title: The Randers metrics of weakly isotropic scalar curvature

abstract: In this paper, we study the Randers metrics of weakly isotropic scalar curvature. We prove that a Randers metric of weakly isotropic scalar curvature must be of isotropic $S$-curvature. Further, we prove that a conformally flat Randers metric of weakly isotropic scalar curvature is either Minkowskian or Riemannian.


title: AstroDance: Engaging Deaf and Hard-of-Hearing Students in Astrophysics via Multimedia Performances

abstract: The dynamics of gravitating astrophysical systems such as black holes and neutron stars are fascinatingly complex, offer some of nature's most spectacular phenomena, and capture the public's imagination in ways that few subjects can. Here, we describe {\it AstroDance}, a multi-media project to engage deaf and hard-of-hearing (DHH) students in astronomy and gravitational physics. {\it AstroDance} incorporates multiple means of representation of scientific concepts and was performed primarily for secondary and post-secondary audiences at $\sim$20 venues in the northeastern US prior to the historic first detection of gravitational waves. As part of the {\it AstroDance} project, we surveyed $\sim$1000 audience members roughly split evenly between hearing and DHH audience members. While both groups reported statistically equivalent high-rates of enjoyment of the performance, the DHH group reported an increase in how much they learned about science at a statistically significant rate compared to the hearing audience. Our findings suggest that multi-sensory approaches benefit both hearing and deaf audiences and enable accessible participation for broader groups.

categories: physics.ed-ph astro-ph.EP astro-ph.GA astro-ph.SR


title: Methods of Stellar Space-Density Analyses: A Retrospective Review

abstract: Stellar space density analyses, once a very active area of astronomical research, involved transforming counts of stars to given limiting magnitudes in selected areas of the sky into graphs of the number of stars per cubic parsec as a function of distance. Several methods of computing the transformation have been published, varying from manual spreadsheets to computer programs based on very sophisticated mathematical techniques for deconvolving integral equations. This paper revisits these techniques and compares the performance of seven of them published between 1913 and 2003 as applied to both simulated and real data.


title: On the Largest Common Subtree of Random Leaf-Labeled Binary Trees

abstract: The size of the largest common subtree (maximum agreement subtree) of two independent uniform random binary trees on $n$ leaves is known to be between orders $n^{1/8}$ and $n^{1/2}$. By a construction based on recursive splitting and analyzable by standard "stochastic fragmentation" methods, we improve the lower bound to order $n^\beta$ for $\beta = \frac{\sqrt{3} - 1}{2} = 0.366$. Improving the upper bound remains a challenging problem.


title: Cauchy-Szeg\"o operator, quaternionic Siegel upper half space, commutator, weighted Morrey space

abstract: In the setting of quaternionic Heisenberg group $\mathscr H^{n-1}$, we characterize the boundedness and compactness of commutator $[b,\mathcal C]$ for the Cauchy--Szeg\"o operator $\mathcal C$ on the weighted Morrey space $L_w^{p,\,\kappa}(\mathscr H^{n-1})$ with $p\in(1, \infty)$, $\kappa\in(0, 1)$ and $w\in A_p(\mathscr H^{n-1}).$ More precisely, we prove that $[b,\mathcal C]$ is bounded on $L_w^{p,\,\kappa}(\mathscr H^{n-1})$ if and only if $b\in {\rm BMO}(\mathscr H^{n-1})$. And $[b,\mathcal C]$ is compact on $L_w^{p,\,\kappa}(\mathscr H^{n-1})$ if and only if $b\in {\rm VMO}(\mathscr H^{n-1})$.


title: MOSQUITO-NET: A deep learning based CADx system for malaria diagnosis along with model interpretation using GradCam and class activation maps

abstract: Malaria is considered one of the deadliest diseases in today world which causes thousands of deaths per year. The parasites responsible for malaria are scientifically known as Plasmodium which infects the red blood cells in human beings. The parasites are transmitted by a female class of mosquitos known as Anopheles. The diagnosis of malaria requires identification and manual counting of parasitized cells by medical practitioners in microscopic blood smears. Due to the unavailability of resources, its diagnostic accuracy is largely affected by large scale screening. State of the art Computer-aided diagnostic techniques based on deep learning algorithms such as CNNs, with end to end feature extraction and classification, have widely contributed to various image recognition tasks. In this paper, we evaluate the performance of custom made convnet Mosquito-Net, to classify the infected and uninfected cells for malaria diagnosis which could be deployed on the edge and mobile devices owing to its fewer parameters and less computation power. Therefore, it can be wildly preferred for diagnosis in remote and countryside areas where there is a lack of medical facilities.


title: Full classification of dynamics for one-dimensional continuous time Markov chains with polynomial transition rates

abstract: This paper provides full classification of dynamics for continuous time Markov chains (CTMCs) on the non-negative integers with polynomial transition rate functions. Such stochastic processes are abundant in applications, in particular in biology. More precisely, for CTMCs of bounded jumps, we provide necessary and sufficient conditions in terms of calculable parameters for explosivity, recurrence vs transience, certain absorption, positive recurrence vs null recurrence, and implosivity. Simple sufficient conditions for exponential ergodicity of stationary distributions and quasi-stationary distributions as well as existence and non-existence of moments of hitting times are also obtained. Similar simple sufficient conditions for the aforementioned dynamics together with their opposite dynamics are established for CTMCs with unbounded jumps. The results generalize respective criteria for birth-death processes by Karlin and McGregor in the 1960s. Finally, we apply our results to stochastic reaction networks, an extended class of branching processes, a general bursty single-cell stochastic gene expression model, and population processes, none of which are birth-death processes. The approach is based on a mixture of Lyapunov-Foster type results, semimartingale approach, as well as estimates of stationary measures.


title: Locally harmonic Maass forms and periods of meromorphic modular forms

abstract: We investigate a new family of locally harmonic Maass forms which correspond to periods of modular forms. They transform like negative weight modular forms and are harmonic apart from jump singularities along infinite geodesics. Our main result is an explicit splitting of the new locally harmonic Maass forms into a harmonic part and a locally polynomial part that captures the jump singularities. As an application, we obtain finite rational formulas for suitable linear combinations of periods of meromorphic modular forms associated to positive definite binary quadratic forms.


title: ChestX-Det10: Chest X-ray Dataset on Detection of Thoracic Abnormalities

abstract: Instance level detection of thoracic diseases or abnormalities are crucial for automatic diagnosis in chest X-ray images. Most existing works on chest X-rays focus on disease classification and weakly supervised localization. In order to push forward the research on disease classification and localization on chest X-rays. We provide a new benchmark called ChestX-Det10, including box-level annotations of 10 categories of disease/abnormality of $\sim$ 3,500 images. The annotations are located at https://github.com/Deepwise-AILab/ChestX-Det10-Dataset.


title: On exact solutions, conservation laws and invariant analysis of the generalized Rosenau-Hyman equation

abstract: In this paper, the nonlinear Rosenau-Hyman equation with time dependent variable coefficients is considered for investigating its invariant properties, exact solutions and conservation laws. Using Lie classical method, we derive symmetries admitted by considered equation. Symmetry reductions are performed for each components of optimal set. Also nonclassical approach is employed on considered equation to find some additional supplementary symmetries and corresponding symmetry reductions are performed. Later three kinds of exact solutions of considered equation are presented graphically for different parameters. In addition, local conservation laws are constructed for considered equation by multiplier approach.


title: XRayGAN: Consistency-preserving Generation of X-ray Images from Radiology Reports

abstract: To effectively train medical students to become qualified radiologists, a large number of X-ray images collected from patients with diverse medical conditions are needed. However, due to data privacy concerns, such images are typically difficult to obtain. To address this problem, we develop methods to generate view-consistent, high-fidelity, and high-resolution X-ray images from radiology reports to facilitate radiology training of medical students. This task is presented with several challenges. First, from a single report, images with different views (e.g., frontal, lateral) need to be generated. How to ensure consistency of these images (i.e., make sure they are about the same patient)? Second, X-ray images are required to have high resolution. Otherwise, many details of diseases would be lost. How to generate high-resolutions images? Third, radiology reports are long and have complicated structure. How to effectively understand their semantics to generate high-fidelity images that accurately reflect the contents of the reports? To address these three challenges, we propose an XRayGAN composed of three modules: (1) a view consistency network that maximizes the consistency between generated frontal-view and lateral-view images; (2) a multi-scale conditional GAN that progressively generates a cascade of images with increasing resolution; (3) a hierarchical attentional encoder that learns the latent semantics of a radiology report by capturing its hierarchical linguistic structure and various levels of clinical importance of words and sentences. Experiments on two radiology datasets demonstrate the effectiveness of our methods. To our best knowledge, this work represents the first one generating consistent and high-resolution X-ray images from radiology reports. The code is available at https://github.com/UCSD-AI4H/XRayGAN.


title: Artificial Musical Intelligence: A Survey

abstract: Computers have been used to analyze and create music since they were first introduced in the 1950s and 1960s. Beginning in the late 1990s, the rise of the Internet and large scale platforms for music recommendation and retrieval have made music an increasingly prevalent domain of machine learning and artificial intelligence research. While still nascent, several different approaches have been employed to tackle what may broadly be referred to as "musical intelligence." This article provides a definition of musical intelligence, introduces a taxonomy of its constituent components, and surveys the wide range of AI methods that can be, and have been, brought to bear in its pursuit, with a particular emphasis on machine learning methods.


title: Dimits shift, avalanche-like bursts, and Solitary propagating structures in the two-field Flux-Balanced Hasegawa-Wakatani model for plasma edge turbulence

abstract: We show that the recently introduced two-field flux-balanced Hasegawa-Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa-Wakatani models, which confirms the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multiscale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.


title: Sparse HP Filter: Finding Kinks in the COVID-19 Contact Rate

abstract: In this paper, we estimate the time-varying COVID-19 contact rate of a Susceptible-Infected-Recovered (SIR) model. Our measurement of the contact rate is constructed using data on actively infected, recovered and deceased cases. We propose a new trend filtering method that is a variant of the Hodrick-Prescott (HP) filter, constrained by the number of possible kinks. We term it the $\textit{sparse HP filter}$ and apply it to daily data from five countries: Canada, China, South Korea, the UK and the US. Our new method yields the kinks that are well aligned with actual events in each country. We find that the sparse HP filter provides a fewer kinks than the $\ell_1$ trend filter, while both methods fitting data equally well. Theoretically, we establish risk consistency of both the sparse HP and $\ell_1$ trend filters. Ultimately, we propose to use time-varying $\textit{contact growth rates}$ to document and monitor outbreaks of COVID-19.


title: Electromagnetic Waves in Hot and Dense Media

abstract: It is known that the finite temperature and density (FTD) corrections increase the electron mass, charge and modify the properties of the emitted radiation. All the signals, travelling through the astronomical bodies, carry over the information of their origin and bring minor details about the structure and composition of the source. It has been noticed that temperatures of the early universe add physically measureable mass to electron and large chemical potential lead to an increase in mass as well. However, the QED coupling is slowly increased with temperature but is decreased with increasing chemical potential in superdense systems. Due to the strong relationship between the properties of electromagnetic signals and composition of the material of the propagating objects, we propose to incorporate the modified signal properties to analyze the astronomical data and its interpretation. Existence of different phases of stellar matter at the quantum scale indicates frequent phase transitions for rapidly rotating superdense materials. This information is related to the fact that the allowed range of temperature and chemical potential for a physical system is determined by the signal properties and properties changes with the material as well..


title: The navigation problems and the curvature properties on conic Kropina manifolds

abstract: In this paper, we study navigation problems on conic Kropina manifolds. Let $F(x, y)$ be a conic Kropina metric on an $n$-dimensional manifold $M$ and $V$ be a conformal vector field on $(M, F)$ with $F(x, - V_{x})\leq 1$. Let $\widetilde{F}= \widetilde{F} (x,y)$ be the solution of the navigation problem with navigation data $(F, V)$. We prove that $\widetilde{F}$ must be either a Randers metric or a Kropina metric. Then we establish the relationships between some curvature properties of $F$ and the corresponding properties of the new metric $\widetilde{F}$, which involve S-curvature, flag curvature and Ricci curvature.


title: On the radiation-loud and radiation-quiet anti-glitches in magnetar 1E 2259+586

abstract: The radiation-loud and radiation-quiet anti-glitches in magnetar 1E 2259+586 are discussed together. It is shown that the 2019 radiation-quiet anti-glitch in 1E 2259+586 is not constraining, because its amplitude is 2.5 smaller than the previous radiation-loud anti-glitch. We note that magnetars with high X-ray luminosity and low rotational energy loss rate are more likely to show anti-glitches. Anti-glitchers and radio-emitters of magnetars occupy opposite corner on the luminosity vs. rotational energy loss rate diagram. By introducing the X-ray luminosity parameter and the activity parameter, we try to characterize this trend quantitatively. A general picture of magnetar anti-glitch is given.


title: Differentially-private Federated Neural Architecture Search

abstract: Neural architecture search, which aims to automatically search for architectures (e.g., convolution, max pooling) of neural networks that maximize validation performance, has achieved remarkable progress recently. In many application scenarios, several parties would like to collaboratively search for a shared neural architecture by leveraging data from all parties. However, due to privacy concerns, no party wants its data to be seen by other parties. To address this problem, we propose federated neural architecture search (FNAS), where different parties collectively search for a differentiable architecture by exchanging gradients of architecture variables without exposing their data to other parties. To further preserve privacy, we study differentially-private FNAS (DP-FNAS), which adds random noise to the gradients of architecture variables. We provide theoretical guarantees of DP-FNAS in achieving differential privacy. Experiments show that DP-FNAS can search highly-performant neural architectures while protecting the privacy of individual parties. The code is available at https://github.com/UCSD-AI4H/DP-FNAS


title: Gradient Amplification: An efficient way to train deep neural networks

abstract: Improving performance of deep learning models and reducing their training times are ongoing challenges in deep neural networks. There are several approaches proposed to address these challenges one of which is to increase the depth of the neural networks. Such deeper networks not only increase training times, but also suffer from vanishing gradients problem while training. In this work, we propose gradient amplification approach for training deep learning models to prevent vanishing gradients and also develop a training strategy to enable or disable gradient amplification method across several epochs with different learning rates. We perform experiments on VGG-19 and resnet (Resnet-18 and Resnet-34) models, and study the impact of amplification parameters on these models in detail. Our proposed approach improves performance of these deep learning models even at higher learning rates, thereby allowing these models to achieve higher performance with reduced training time.


title: An electromagnetic extension of the Schwarzschild interior solution and the corresponding Buchdahl limit

abstract: We wish to construct a model for charged star as a generalization of the uniform density Schwarzschild interior solution. We employ the Vaidya and Tikekar ansatz [{\it Astrophys. Astron.} {\bf 3} (1982) 325] for one of the metric potentials and electric field is chosen in such a way that when it is switched off the metric reduces to the Schwarzschild. This relates charge distribution to the Vaidya-Tikekar parameter, $k$, indicating deviation form sphericity of three dimensional space when embedded into four dimensional Euclidean space. The model is examined against all the physical conditions required for a relativistic charged fluid sphere as an interior to a charged star. We also obtain and discuss charged analogue of the Buchdahl compactness bound.


title: Uncertainty in Gradient Boosting via Ensembles

abstract: For many practical, high-risk applications, it is essential to quantify uncertainty in a model's predictions to avoid costly mistakes. While predictive uncertainty is widely studied for neural networks, the topic seems to be under-explored for models based on gradient boosting. However, gradient boosting often achieves state-of-the-art results on tabular data. This work examines a probabilistic ensemble-based framework for deriving uncertainty estimates in the predictions of gradient boosting classification and regression models. We conducted experiments on a range of synthetic and real datasets and investigated the applicability of ensemble approaches to gradient boosting models that are themselves ensembles of decision trees. Our analysis shows that ensembles of gradient boosting models successfully detect anomalous inputs while having limited ability to improve the predicted total uncertainty. Importantly, we also propose a concept of a virtual ensemble to get the benefits of an ensemble via only one gradient boosting model, which significantly reduces complexity.


title: A self-similar solution to time-fractional Stefan problem

abstract: We derive the fractional version of one-phase one-dimensional Stefan model. We assume that the diffusive flux is given by the time-fractional Riemann-Liouville derivative, i.e. we impose the memory effect in the examined model. Furthermore, we find a special solution to this problem.


title: Distribution-free binary classification: prediction sets, confidence intervals and calibration

abstract: We study three notions of uncertainty quantification -- calibration, confidence intervals and prediction sets -- for binary classification in the distribution-free setting, that is without making any distributional assumptions on the data. With a focus towards calibration, we establish a 'tripod' of theorems that connect these three notions for score-based classifiers. A direct implication is that distribution-free calibration is only possible, even asymptotically, using a scoring function whose level sets partition the feature space into at most countably many sets. Parametric calibration schemes such as variants of Platt scaling do not satisfy this requirement, while nonparametric schemes based on binning do. To close the loop, we derive distribution-free confidence intervals for binned probabilities for both fixed-width and uniform-mass binning. As a consequence of our 'tripod' theorems, these confidence intervals for binned probabilities lead to distribution-free calibration. We also derive extensions to settings with streaming data and covariate shift.

categories: stat.ML cs.AI cs.LG math.ST stat.ME stat.TH


title: Using Computational Essays to Scaffold Professional Physics Practice

abstract: This article describes a curricular innovation designed to help students experience authentic physics inquiry with an emphasis on computational modeling and scientific communication. The educational design centers on a new type of assignment called a computational essay, which was developed and implemented over the course of two semesters of an intermediate electricity and magnetism course at the University of Oslo, Norway. We describe the motivation, learning goals, and scaffolds used in the computational essay project, with the intention that other educators will be able to replicate and adapt our design. We also report on initial findings from this implementation, including key features of student-written computational essays, student reflections on the inquiry process, and self-reported conceptual and attitudinal development. Based on these findings, we argue that computational essays can serve a key role in introducing students to open-ended, inquiry-based work and setting the foundation for future computational research and studies.


title: Carrier-envelope phase effects in Laser Wakefield Acceleration with near-single-cycle pulses

abstract: Driving laser wakefield acceleration with extremely short, near single-cycle laser pulses is crucial to the realisation of an electron source that can operate at kHz-repetition rate while relying on modest laser energy. It is also interesting from a fundamental point of view, as the ponderomotive approximation is no longer valid for such short pulses. Through particle-in-cell simulations, we show how the plasma response becomes asymmetric in the plane of laser polarization, and dependent on the carrier-envelope phase (CEP) of the laser pulse. For the case of self-injection, this in turn strongly affects the initial conditions of injected electrons, causing collective betatron oscillations of the electron beam. As a result, the electron beam pointing, electron energy spectrum and the direction of emitted betatron radiation become CEP-dependent. For injection in a density gradient the effect on beam pointing is reduced and the electron energy spectrum is CEP-independent, as electron injection is mostly longitudinal and mainly determined by the density gradient. Our results highlight the importance of controlling the CEP in this regime for producing stable and reproducible relativistic electron beams and identify how CEP effects may be observed in experiments. In the future, CEP control may become an additional tool to control the energy spectrum or pointing of the accelerated electron beam.


title: Approximation of the inverse scattering Steklov eigenvalues and the inverse spectral problem

abstract: In this paper, we consider the numerical approximation of the Steklov eigenvalue problem that arises in inverse acoustic scattering. The underlying scattering problem is for an inhomogeneous isotropic medium. These eigenvalues have been proposed to be used as a target signature since they can be recovered from the scattering data. A Galerkin method is studied where the basis functions are the Neumann eigenfunctions of the Laplacian. Error estimates for the eigenvalues and eigenfunctions are proven by appealing to Weyl's Law. We will test this method against separation of variables in order to validate the theoretical convergence. We also consider the inverse spectral problem of estimating/recovering the refractive index from the knowledge of the Steklov eigenvalues. Since the eigenvalues are monotone with respect to a real-valued refractive index implies that they can be used for non-destructive testing. Some numerical examples are provided for the inverse spectral problem.


title: Eigenvalue asymptotics for polynomially compact pseudodifferenial operators and applications

abstract: We fnd the asymptotics of eigenvalues of polynomially compact zero order pseudodiferential operators, the motivating example being the Neumann- Poincare operator in linear elasticity.


title: Towards a Neural Graphics Pipeline for Controllable Image Generation

abstract: In this paper, we leverage advances in neural networks towards forming a neural rendering for controllable image generation, and thereby bypassing the need for detailed modeling in conventional graphics pipeline. To this end, we present Neural Graphics Pipeline (NGP), a hybrid generative model that brings together neural and traditional image formation models. NGP decomposes the image into a set of interpretable appearance feature maps, uncovering direct control handles for controllable image generation. To form an image, NGP generates coarse 3D models that are fed into neural rendering modules to produce view-specific interpretable 2D maps, which are then composited into the final output image using a traditional image formation model. Our approach offers control over image generation by providing direct handles controlling illumination and camera parameters, in addition to control over shape and appearance variations. The key challenge is to learn these controls through unsupervised training that links generated coarse 3D models with unpaired real images via neural and traditional (e.g., Blinn- Phong) rendering functions, without establishing an explicit correspondence between them. We demonstrate the effectiveness of our approach on controllable image generation of single-object scenes. We evaluate our hybrid modeling framework, compare with neural-only generation methods (namely, DCGAN, LSGAN, WGAN-GP, VON, and SRNs), report improvement in FID scores against real images, and demonstrate that NGP supports direct controls common in traditional forward rendering. Code is available at http://geometry.cs.ucl.ac.uk/projects/2021/ngp.


title: A brief overview of black hole-neutron star mergers

abstract: Of the three main types of binaries detectable through ground-based gravitational wave observations, black hole-neutron star (BHNS) mergers remain the most elusive. While candidates BHNS exist in the triggers released during the third observing run of the Advanced LIGO/Virgo collaboration, no detection has been confirmed so far. As for binary neutron star systems, BHNS binaries allow us to explore a wide range of physical processes, including the neutron star equation of state, nucleosynthesis, stellar evolution, high-energy astrophysics, and the expansion of the Universe. Here, we review some of the main features of BHNS systems: the distinction between disrupting and non-disrupting binaries, the types of outflows that BHNS mergers can produce, and the information that can be extracted from the observation of their gravitational wave and electromagnetic signals. We also emphasize that for the most likely binary parameters, BHNS mergers seem less likely to power electromagnetic signals than binary neutron star systems. Finally, we discuss some of the issues that still limit our ability to model and interpret electromagnetic signals from BHNS binaries.


title: Likelihood-Free Inference with Deep Gaussian Processes

abstract: In recent years, surrogate models have been successfully used in likelihood-free inference to decrease the number of simulator evaluations. The current state-of-the-art performance for this task has been achieved by Bayesian Optimization with Gaussian Processes (GPs). While this combination works well for unimodal target distributions, it is restricting the flexibility and applicability of Bayesian Optimization for accelerating likelihood-free inference more generally. We address this problem by proposing a Deep Gaussian Process (DGP) surrogate model that can handle more irregularly behaved target distributions. Our experiments show how DGPs can outperform GPs on objective functions with multimodal distributions and maintain a comparable performance in unimodal cases. This confirms that DGPs as surrogate models can extend the applicability of Bayesian Optimization for likelihood-free inference (BOLFI), while adding computational overhead that remains negligible for computationally intensive simulators.


title: Gas, dust, and star formation in the positive AGN feedback candidate 4C 41.17 at z=3.8

abstract: We present new, spatially resolved [CI]1-0, [CI]2-1, CO(7-6), and dust continuum observations of 4C 41.17 at $z=3.8$ obtained with the IRAM NOEMA interferometer. This is one of the best-studied radio galaxies in this epoch and is arguably the best candidate of jet-triggered star formation at high redshift currently known in the literature. 4C 41.17 shows a narrow ridge of dust continuum extending over 15 kpc near the radio jet axis. Line emission is found within the galaxy in the region with signatures of positive feedback. Using the [CI]1-0 line as a molecular gas tracer, and multifrequency observations of the far-infrared dust heated by star formation, we find a total gas mass of $7.6\times 10^{10}$ M$_{\odot}$, which is somewhat greater than that previously found from CO(4-3). The gas mass surface density of $10^3$ M$_{\odot}$ yr$^{-1}$ pc$^{-2}$ and the star formation rate surface density of 10 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$ were derived over the 12 kpc$\times$8 kpc area, where signatures of positive feedback have previously been found. These densities are comparable to those in other populations of massive, dusty star-forming galaxies in this redshift range, suggesting that the jet does not currently enhance the efficiency with which stars form from the gas. This is consistent with expectations from simulations, whereby radio jets may facilitate the onset of star formation in galaxies without boosting its efficiency over longer timescales, in particular after the jet has broken out of the interstellar medium, as is the case in 4C 41.17.


title: Quantum-Limited Squeezed Light Detection with a Camera

abstract: We present a technique for squeezed light detection based on direct imaging of the displaced-squeezed-vacuum state using a CCD camera. We show that the squeezing parameter can be accurately estimated using only the first two moments of the recorded pixel-to-pixel photon fluctuation statistics, with accuracy that rivals that of the standard squeezing detection methods such as a balanced homodyne detection. Finally, we numerically simulate the camera operation, reproducing the noisy experimental results with low signal samplings and confirming the theory with high signal samplings.


title: Extenders under ZF and constructibility of rank-to-rank embeddings

abstract: Assume ZF (without the Axiom of Choice). Let $j:V_\varepsilon\to V_\delta$ be a non-trivial $\in$-cofinal $\Sigma_1$-elementary embedding, where $\varepsilon,\delta$ are limit ordinals. We prove some restrictions on the constructibility of $j$ from $V_\delta$, mostly focusing on the case $\varepsilon=\delta$. In particular, if $\varepsilon=\delta$ and $j\in L(V_\delta)$ then $\delta$ has cofinality $\omega$. However, assuming ZFC+I$_3$, with the appropriate $\varepsilon=\delta$, one can force to get such $j\in L(V^{V[G]}_\delta)$. Assuming Dependent Choice and that $\delta$ has cofinality $\omega$ (but not assuming $V=L(V_\delta)$), and $j:V_\delta\to V_\delta$ is $\Sigma_1$-elementary, we show that there are "perfectly many" such $j$, with none being "isolated". Assuming a proper class of weak Lowenheim-Skolem cardinals, we also give a first-order characterization of critical points of embeddings $j:V\to M$ with $M$ transitive. The main results rely on a development of extenders under ZF (which is most useful given such wLS cardinals).


title: A Review of 1D Convolutional Neural Networks toward Unknown Substance Identification in Portable Raman Spectrometer

abstract: Raman spectroscopy is a powerful analytical tool with applications ranging from quality control to cutting edge biomedical research. One particular area which has seen tremendous advances in the past decade is the development of powerful handheld Raman spectrometers. They have been adopted widely by first responders and law enforcement agencies for the field analysis of unknown substances. Field detection and identification of unknown substances with Raman spectroscopy rely heavily on the spectral matching capability of the devices on hand. Conventional spectral matching algorithms (such as correlation, dot product, etc.) have been used in identifying unknown Raman spectrum by comparing the unknown to a large reference database. This is typically achieved through brute-force summation of pixel-by-pixel differences between the reference and the unknown spectrum. Conventional algorithms have noticeable drawbacks. For example, they tend to work well with identifying pure compounds but less so for mixture compounds. For instance, limited reference spectra inaccessible databases with a large number of classes relative to the number of samples have been a setback for the widespread usage of Raman spectroscopy for field analysis applications. State-of-the-art deep learning methods (specifically convolutional neural networks CNNs), as an alternative approach, presents a number of advantages over conventional spectral comparison algorism. With optimization, they are ideal to be deployed in handheld spectrometers for field detection of unknown substances. In this study, we present a comprehensive survey in the use of one-dimensional CNNs for Raman spectrum identification. Specifically, we highlight the use of this powerful deep learning technique for handheld Raman spectrometers taking into consideration the potential limit in power consumption and computation ability of handheld systems.


title: Pion-photon transition form factor in LCSR and tests of asymptotics

abstract: We study the pion-photon transition form factor (TFF) $F^{\gamma*\gamma\pi^0}(Q^2)$ using a state-of-the art implementation of light cone sum rules (LCSRs) within fixed-order QCD perturbation theory. The spectral density in the dispersion relation includes all currently known radiative corrections up to the next-to-next-to-leading-order (NNLO) and all twist contributions up to order six. Predictions for the TFF are obtained for various pion distribution amplitudes (DAs) of twist two, including two-loop evolution which accounts for heavy-quark mass thresholds. The influence of the main theoretical uncertainties is quantified in order to enable a more realistic comparison with the data. The characteristics of various pion DAs are analyzed in terms of the conformal coefficients $a_2$ and $a_4$ in comparison with the $1\sigma$ and $2\sigma$ error regions of the data and the most recent lattice constraints on $a_2$ with NLO and NNLO accuracy. Our results provide more stringent bounds on the variation of the pion DA and illuminate the corresponding asymptotic behavior of the calculated TFF.


title: Effects of dynamical noises on Majorana bound states

abstract: The nonlocal nature of unpaired Majorana bound states (MBSs) in topological superconductors can be exploited to create topologically protected qubits and perform gate operations fault-tolerantly via braidings. However, the time-dependent noises induced by coupling to an environment which is inevitable in any realistic system could spoil the topological protection. In this work, we study the effects of various dynamical noises such as Lorentzian, thermal, and quantum point contact on the MBSs in the recently proposed one-dimensional topological superconductors. We begin by investigating the Kitaev p-wave superconductors and examine the effects of long-range hopping and pairing on the transition rate of MBSs. We found that, especially, the long-range pairings significantly reduce the transition rate of bound states. Then, we consider the recently discovered topological superconducting nanowires and magnetic chains. Our findings are consequential for the recent attempts to manipulate MBSs. In particular, for the latter two experimentally realized systems we argue how low magnetic/Zeeman fields and strong spin-orbit coupling make the MBSs more robust to noises.


title: Dynamical Strengthening of Covalent and Non-Covalent Molecular Interactions by Nuclear Quantum Effects at Finite Temperature

abstract: Nuclear quantum effects (NQE) tend to generate delocalized molecular dynamics due to the inclusion of the zero point energy and its coupling with the anharmonicities in interatomic interactions. Here, we present evidence that NQE often enhance electronic interactions and, in turn, can result in dynamical molecular stabilization at finite temperature. The underlying physical mechanism promoted by NQE depends on the particular interaction under consideration. First, the effective reduction of interatomic distances between functional groups within a molecule can enhance the $n\to\pi^*$ interaction by increasing the overlap between molecular orbitals or by strengthening electrostatic interactions between neighboring charge densities. Second, NQE can localize methyl rotors by temporarily changing molecular bond orders and leading to the emergence of localized transient rotor states. Third, for noncovalent van der Waals interactions the strengthening comes from the increase of the polarizability given the expanded average interatomic distances induced by NQE. The implications of these boosted interactions include counterintuitive hydroxyl--hydroxyl bonding, hindered methyl rotor dynamics, and molecular stiffening which generates smoother free-energy surfaces. Our findings yield new insights into the versatile role of nuclear quantum fluctuations in molecules and materials.

categories: physics.chem-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.comp-ph quant-ph


title: Uniqueness of STFT phase retrieval for bandlimited functions

abstract: We consider the problem of phase retrieval from magnitudes of short-time Fourier transform (STFT) measurements. It is well-known that signals are uniquely determined (up to global phase) by their STFT magnitude when the underlying window has an ambiguity function that is nowhere vanishing. It is less clear, however, what can be said in terms of unique phase-retrievability when the ambiguity function of the underlying window vanishes on some of the time-frequency plane. In this short note, we demonstrate that by considering signals in Paley-Wiener spaces, it is possible to prove new uniqueness results for STFT phase retrieval. Among those, we establish a first uniqueness theorem for STFT phase retrieval from magnitude-only samples in a real-valued setting.


title: Sharp Estimates for Blowing Down Functions in a Denjoy-Carleman Class

abstract: If F is an infinitely differentiable function whose composition with a blowing-up belongs to a Denjoy-Carleman class C_M (determined by a log convex sequence M=(M_k)), then F, in general, belongs to a larger shifted class C_N, where N_k = M_2k; i.e., there is a loss of regularity. We show that this loss of regularity is sharp. In particular, loss of regularity of Denjoy-Carleman classes is intrinsic to arguments involving resolution of singularities.


title: Minimax Lower Bounds for Transfer Learning with Linear and One-hidden Layer Neural Networks

abstract: Transfer learning has emerged as a powerful technique for improving the performance of machine learning models on new domains where labeled training data may be scarce. In this approach a model trained for a source task, where plenty of labeled training data is available, is used as a starting point for training a model on a related target task with only few labeled training data. Despite recent empirical success of transfer learning approaches, the benefits and fundamental limits of transfer learning are poorly understood. In this paper we develop a statistical minimax framework to characterize the fundamental limits of transfer learning in the context of regression with linear and one-hidden layer neural network models. Specifically, we derive a lower-bound for the target generalization error achievable by any algorithm as a function of the number of labeled source and target data as well as appropriate notions of similarity between the source and target tasks. Our lower bound provides new insights into the benefits and limitations of transfer learning. We further corroborate our theoretical finding with various experiments.


title: Flavor dependence of the thermal dissociations of vector and axial-vector mesons

abstract: The in-medium behavior of ground-state $q\bar{q}$ mesons, where $q \in \{u,d,s,c\}$, in vector and axial-vector channels is studied based on the spectral analysis for mesonic correlators at finite temperature and zero chemical potential. We first compute the correlators by solving the quark gap equations and the inhomogeneous Bethe-Salpeter equations in the rainbow-ladder approximation. Using a phenomenological ansatz, the spectral functions are extracted by fitting the correlators. By analyzing the evolution of the spectral functions with the temperature, we obtain the dissociation temperatures of mesons and discuss their relations to the critical temperature of the chiral symmetry restoration. The results show a pattern of the flavor dependence of the thermal dissociation of the mesons.


title: Particle Physics at Accelerators in the United States and Asia

abstract: Particle physics experiments in the United States and Asia have greatly contributed to the understanding of elementary particles and their interactions. With the recent discovery of the Higgs boson at CERN, interest in the development of next-generation colliders has been rekindled. A linear electron-positron collider in Japan and a circular collider in China have been proposed for precision studies of the Higgs boson. In addition to the Higgs programme, new accelerator-based long-baseline neutrino mega-facilities are being built in the United States and Japan. Here, we outline the present status of key particle physics programmes at accelerators and future plans in the United States and Asia that largely complement approaches being explored in the European Strategy for Particle Physics Update. We encourage the pursuit of this global approach, reaching beyond regional boundaries for optimized development and operations of major accelerator facilities worldwide, to ensure active and productive future of the field.


title: Review of COVID-19 Antibody Therapies

abstract: Under the global health emergency caused by coronavirus disease 2019 (COVID-19), efficient and specific therapies are urgently needed. Compared with traditional small-molecular drugs, antibody therapies are relatively easy to develop and as specific as vaccines in targeting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and thus attract much attention in the past few months. This work reviews seven existing antibodies for SARS-CoV-2 spike (S) protein with three-dimensional (3D) structures deposited in the Protein Data Bank. Five antibody structures associated with SARS-CoV are evaluated for their potential in neutralizing SARS-CoV-2. The interactions of these antibodies with the S protein receptor-binding domain (RBD) are compared with those of angiotensin-converting enzyme 2 (ACE2) and RBD complexes. Due to the orders of magnitude in the discrepancies of experimental binding affinities, we introduce topological data analysis (TDA), a variety of network models, and deep learning to analyze the binding strength and therapeutic potential of the aforementioned fourteen antibody-antigen complexes. The current COVID-19 antibody clinical trials, which are not limited to the S protein target, are also reviewed.


title: Process tomography of Robust Dynamical Decoupling in Superconducting Qubits

abstract: Dynamical decoupling is a technique that protects qubits against noise. The ability to preserve quantum coherence in the presence of noise is essential for the development of quantum devices. Here the Rigetti quantum computing platform was used to test different dynamical decoupling sequences. The performance of the sequences was characterized by Quantum Process Tomography and analyzed using the quantum channels formalism. It is shown that dynamical decoupling can reduce qubit dephasing but cannot protect against spontaneous emission. Furthermore, from process tomography results, it was also possible to conclude that the action of dynamical decoupling cannot be understood as a simple modification of the qubit coherence time. It is also shown here that the performance of dynamical decoupling on the Rigetti's qubits is limited by pulse imperfections. However, the performance can be improved using robust dynamical decoupling, i.e. sequences that are robust against experimental imperfections. The sequences tested here outperformed previous dynamical decoupling sequences tested in the same platform.


title: On Calder\'on's inverse inclusion problem with smooth shapes by a single partial boundary measurement

abstract: We are concerned with the Calder\'on inverse inclusion problem, where one intends to recover the shape of an inhomogeneous conductive inclusion embedded in a homogeneous conductivity by the associated boundary measurements. We consider the highly challenging case with a single partial boundary measurement, which constitutes a long-standing open problem in the literature. It is shown in several existing works that corner singularities can help to resolve the uniqueness and stability issues for this inverse problem. In this paper, we show that the corner singularity can be relaxed to be a certain high-curvature condition and derive a novel local unique determination result. To our best knowledge, this is the first (local) uniqueness result in determining a conductive inclusions with general smooth shapes by a single (partial) boundary measurement.


title: Lightweight Collaborative Anomaly Detection for the IoT using Blockchain

abstract: Due to their rapid growth and deployment, the Internet of things (IoT) have become a central aspect of our daily lives. Unfortunately, IoT devices tend to have many vulnerabilities which can be exploited by an attacker. Unsupervised techniques, such as anomaly detection, can be used to secure these devices in a plug-and-protect manner. However, anomaly detection models must be trained for a long time in order to capture all benign behaviors. Furthermore, the anomaly detection model is vulnerable to adversarial attacks since, during the training phase, all observations are assumed to be benign. In this paper, we propose (1) a novel approach for anomaly detection and (2) a lightweight framework that utilizes the blockchain to ensemble an anomaly detection model in a distributed environment. Blockchain framework incrementally updates a trusted anomaly detection model via self-attestation and consensus among the IoT devices. We evaluate our method on a distributed IoT simulation platform, which consists of 48 Raspberry Pis. The simulation demonstrates how the approach can enhance the security of each device and the security of the network as a whole.


title: Low-Rank Parity-Check Codes over Galois Rings

abstract: Low-rank parity-check (LRPC) are rank-metric codes over finite fields, which have been proposed by Gaborit et al. (2013) for cryptographic applications. Inspired by a recent adaption of Gabidulin codes to certain finite rings by Kamche et al. (2019), we define and study LRPC codes over Galois rings - a wide class of finite commutative rings. We give a decoding algorithm similar to Gaborit et al.'s decoder, based on simple linear-algebraic operations. We derive an upper bound on the failure probability of the decoder, which is significantly more involved than in the case of finite fields. The bound depends only on the rank of an error, i.e., is independent of its free rank. Further, we analyze the complexity of the decoder. We obtain that there is a class of LRPC codes over a Galois ring that can decode roughly the same number of errors as a Gabidulin code with the same code parameters, but faster than the currently best decoder for Gabidulin codes. However, the price that one needs to pay is a small failure probability, which we can bound from above.


title: Random walks on randomly evolving graphs

abstract: A random walk is a basic stochastic process on graphs and a key primitive in the design of distributed algorithms. One of the most important features of random walks is that, under mild conditions, they converge to a stationary distribution in time that is at most polynomial in the size of the graph. This fundamental property, however, only holds if the graph does not change over time, while on the other hand many distributed networks are inherently dynamic, and their topology is subjected to potentially drastic changes. In this work we study the mixing (i.e., converging) properties of random walks on graphs subjected to random changes over time. Specifically, we consider the edge-Markovian random graph model: for each edge slot, there is a two-state Markov chain with transition probabilities $p$ (add a non-existing edge) and $q$ (remove an existing edge). We derive several positive and negative results that depend on both the density of the graph and the speed by which the graph changes. We show that if $p$ is very small (i.e., below the connectivity threshold of Erd\H{o}s-R\'{e}nyi random graphs), random walks do not mix (fast). When $p$ is larger, instead, we observe the following behavior: if the graph changes slowly over time (i.e., $q$ is small), random walks enjoy strong mixing properties that are comparable to the ones possessed by random walks on static graphs; however, if the graph changes too fast (i.e., $q$ is large), only coarse mixing properties are preserved.


title: Restriction of Scalars Chabauty and the $S$-unit equation

abstract: Given a smooth, proper, geometrically integral curve $X$ of genus $g$ with Jacobian $J$ over a number field $K$, Chabauty's method is a $p$-adic technique to bound $\# X(K)$ when $\mathrm{rank}\ J(K) < g$. We study limitations of a variant called `Restriction of Scalars Chabauty' (RoS Chabauty). RoS Chabauty typically bounds $\# X(K)$ when $\mathrm{rank}\ J(K) \leq [K:\mathbb{Q}] (g - 1)$, but fails in the presence of a subgroup obstruction, a high-rank subgroup scheme of $\mathrm{Res}_{K/\mathbb{Q}} J$ which intersects the image of $\mathrm{Res}_{K/\mathbb{Q}} X$ in higher-than-expected dimension. We define BCP obstructions, which are certain subgroup obstructions arising from the geometry of $X$. BCP obstructions explain all known examples where RoS Chabauty fails to bound $\# X(K)$. We also extend RoS Chabauty to compute $S$-integral points on affine curves. Suppose $K$ does not contain a CM-subfield. We present a $p$-adic algorithm which conjecturally computes solutions to the $S$-unit equation $x+y = 1$ for $x,y \in \mathcal{O}_{K,S}^{\times}$ by using RoS Chabauty to compute $S$-integral points on certain genus $0$ affine curves. As evidence the algorithm succeeds, we prove that all but one of these curves have no subgroup obstructions and that the remaining curve has no BCP obstructions to RoS Chabauty. In contrast, under a generalized Leopoldt conjecture, we prove that analogous methods using classical Chabauty cannot bound solutions to the $S$-unit equation when $[K:\mathbb Q] \geq 3$ and $K$ is not totally real.


title: Security Considerations for Internet of Things: A Survey

abstract: Interconnecting "things" and devices that take the form of wearables, sensors, actuators, mobiles, computers, meters, or even vehicles is a critical requirement for the current era. These internetworked connections are serving the emerging applications home and building automation, smart cities and infrastructure, smart industries, and smart-everything. However, the security of these connected Internet of things (IoT) plays a centric role with no margin for error. After a review of the relevant, online literature on the topic and after looking at the market trends and developments, one can notice that there are still concerns with regard to security in IoT products and services. This paper is focusing on a survey on IoT security and aims to highlight the most significant problems related to safety and security in the IoT ecosystems. This survey identifies the general threat and attack vectors against IoT devices while highlighting the flaws and weak points that can lead to breaching the security. Furthermore, this paper presents solutions for remediation of the compromised security, as well as methods for risk mitigation, with prevention and improvement suggestions.