title: Fast tree-based algorithms for DBSCAN for low-dimensional data on GPUs

abstract: DBSCAN is a well-known density-based clustering algorithm to discover arbitrary shape clusters. While conceptually simple in serial, the algorithm is challenging to efficiently parallelize on manycore GPU architectures. Common pitfalls, such as asynchronous range query calls, result in high thread execution divergence in many implementations. In this paper, we propose a new framework for GPU-accelerated DBSCAN, and describe two tree-based algorithms within that framework. Both algorithms fuse the search for neighbors with updating cluster information, but differ in their treatment of dense regions of the data. We show that the time taken to compute clusters is at most twice that of determination of the neighbors. We compare the proposed algorithms with existing CPU and GPU implementations, and demonstrate their competitiveness and performance using a fast traversal structure (bounding volume hierarchy) for low dimensional data. We also show that the memory usage can be reduced by processing object neighbors dynamically without storing them.


title: Equidistribution of rational subspaces and their shapes

abstract: To any $k$-dimensional subspace of $\mathbb Q^n$ one can naturally associate a point in the Grassmannian ${\rm Gr}_{n,k}(\mathbb R)$ and two shapes of lattices of rank $k$ and $n-k$ respectively. These lattices originate by intersecting the $k$-dimensional subspace with the lattice $\mathbb Z^n$. Using unipotent dynamics we prove simultaneous equidistribution of all of these objects under a congruence conditions when $(k,n) \neq (2,4)$.


title: Anomalous entities detection using a cascade of deep learning models

abstract: Human actions that do not conform to usual behavior are considered as anomalous and such actors are called anomalous entities. Detection of anomalous entities using visual data is a challenging problem in computer vision. This paper presents a new approach to detect anomalous entities in complex situations of examination halls. The proposed method uses a cascade of deep convolutional neural network models. In the first stage, we apply a pretrained model of human pose estimation on frames of videos to extract key feature points of body. Patches extracted from each key point are utilized in the second stage to build a densely connected deep convolutional neural network model for detecting anomalous entities. For experiments we collect a video database of students undertaking examination in a hall. Our results show that the proposed method can detect anomalous entities and warrant unusual behavior with high accuracy.


title: High pressure effects on La(O,F)BiS$_{\mathrm{2}}$ single crystal using diamond anvil cell with dual-probe diamond electrodes

abstract: The high-pressure phase of La(O,F)BiS$_{\mathrm{2}}$ exhibits the highest transition temperature among all the BiS$_{\mathrm{2}}$-based superconductors. Various studies, such as investigation of isotope effects, have been conducted to explain its superconducting mechanism. However, there are very few reports on the electrical transport properties and vibration modes of single crystalline La(O,F)BiS$_{\mathrm{2}}$ under high pressure. In this study, we developed a diamond anvil cell with dual-probe diamond electrodes to measure the electrical transport properties of La(O,F)BiS$_{\mathrm{2}}$ single crystal and Pb as a manometer at low temperature. Using the developed system, a linear decrease in the transition temperature and phonon hardening was observed under the application of pressure on La(O,F)BiS$_{\mathrm{2}}$ single crystal.


title: N-type doping of LPCVD-grown \b{eta}-Ga2O3 thin films using solid-source germanium

abstract: We report on the growth and characterization of Ge-doped \b{eta}-Ga2O3 thin films using a solid germanium source. \b{eta}-Ga2O3 thin films were grown using a low-pressure chemical vapor deposition (LPCVD) reactor with either an oxygen or gallium delivery tube. Films were grown on 6 degree offcut sapphire and (010) \b{eta}-Ga2O3 substrates with growth rates between 0.5 - 22 {\mu}m/hr. By controlling the germanium vapor pressure, a wide range of Hall carrier concentrations between 10^17 - 10^19 cm-3 were achieved. Low-temperature Hall data revealed a difference in donor incorporation depending on the reactor configuration. At low growth rates, germanium occupied a single donor energy level between 8 - 10 meV. At higher growth rates, germanium doping predominantly results in a deeper donor energy level at 85 meV. This work shows the effect of reactor design and growth regime on the kinetics of impurity incorporation. Studying donor incorporation in \b{eta}-Ga2O3 is important for the design of high-power electronic devices.


title: Improving Document-Level Sentiment Classification Using Importance of Sentences

abstract: Previous researchers have considered sentiment analysis as a document classification task, in which input documents are classified into predefined sentiment classes. Although there are sentences in a document that support important evidences for sentiment analysis and sentences that do not, they have treated the document as a bag of sentences. In other words, they have not considered the importance of each sentence in the document. To effectively determine polarity of a document, each sentence in the document should be dealt with different degrees of importance. To address this problem, we propose a document-level sentence classification model based on deep neural networks, in which the importance degrees of sentences in documents are automatically determined through gate mechanisms. To verify our new sentiment analysis model, we conducted experiments using the sentiment datasets in the four different domains such as movie reviews, hotel reviews, restaurant reviews, and music reviews. In the experiments, the proposed model outperformed previous state-of-the-art models that do not consider importance differences of sentences in a document. The experimental results show that the importance of sentences should be considered in a document-level sentiment classification task.


title: Stochastic Entry Guidance

abstract: In this paper, closed-loop entry guidance in a randomly perturbed atmosphere, using bank angle control, is posed as a stochastic optimal control problem. The entry trajectory, as well as the closed-loop controls, are both modeled as random processes with statistics determined by the entry dynamics, the entry guidance, and the probabilistic structure of altitude-dependent atmospheric density variations. The entry guidance, which is parameterized as a sequence of linear feedback gains, is designed to steer the probability distribution of the entry trajectories while satisfying bounds on the allowable control inputs and on the maximum allowable state errors. Numerical simulations of a Mars entry scenario demonstrate improved range targeting performance with approximately 50% lower 1st and 99th percentile final range errors when using the developed stochastic guidance scheme as compared to the existing Apollo final phase algorithm.


title: Coherence Protection of Electron Spin in Earth-field Range by All-optical Dynamic Decoupling

abstract: In recent years, unshielded atomic systems have been attracting researchers' attention, in which decoherence is one of the major problems, especially for high precision measurements. The nonlinear Zeeman effect and magnetic field gradient are the main decoherence sources of atomic electron spin in Earth-field range. Here, we propose a method to cancel out the two dominant broadening effects simultaneously by an all-optical dynamic decoupling approach based on Raman scattering in the 87Rb Zeeman sublevels. By adjusting the parameters of the Raman lasers, we realize spin control along an arbitrary direction. We analyze the state evolution of atomic spin under the Raman light control sequence in detail. The results show that both the nonlinear Zeeman effect and magnetic field gradient can be significantly suppressed.


title: Sequential Learning on Liver Tumor Boundary Semantics and Prognostic Biomarker Mining

abstract: The boundary of tumors (hepatocellular carcinoma, or HCC) contains rich semantics: capsular invasion, visibility, smoothness, folding and protuberance, etc. Capsular invasion on tumor boundary has proven to be clinically correlated with the prognostic indicator, microvascular invasion (MVI). Investigating tumor boundary semantics has tremendous clinical values. In this paper, we propose the first and novel computational framework that disentangles the task into two components: spatial vertex localization and sequential semantic classification. (1) A HCC tumor segmentor is built for tumor mask boundary extraction, followed by polar transform representing the boundary with radius and angle. Vertex generator is used to produce fixed-length boundary vertices where vertex features are sampled on the corresponding spatial locations. (2) The sampled deep vertex features with positional embedding are mapped into a sequential space and decoded by a multilayer perceptron (MLP) for semantic classification. Extensive experiments on tumor capsule semantics demonstrate the effectiveness of our framework. Mining the correlation between the boundary semantics and MVI status proves the feasibility to integrate this boundary semantics as a valid HCC prognostic biomarker.


title: An improvement to the vertex-splitting conjecture

abstract: For a simple graph $G$, denote by $n$, $\Delta(G)$, and $\chi'(G)$ its order, maximum degree, and chromatic index, respectively. A connected class 2 graph $G$ is edge-chromatic critical if $\chi'(G-e)<\Delta(G)+1$ for every edge $e$ of $G$. Define $G$ to be overfull if $|E(G)|>\Delta(G) \lfloor n/2 \rfloor$. Clearly, overfull graphs are class 2 and any graph obtained from a regular graph of even order by splitting a vertex is overfull. Let $G$ be an $n$-vertex connected regular class 1 graph with $\Delta(G) >n/3$. Hilton and Zhao in 1997 conjectured that if $G^*$ is obtained from $G$ by splitting one vertex of $G$ into two vertices, then $G^*$ is edge-chromatic critical, and they verified the conjecture for graphs $G$ with $\Delta(G)\ge \frac{n}{2}(\sqrt{7}-1)\approx 0.82n$. The graph $G^*$ is easily verified to be overfull, and so the hardness of the conjecture lies in showing that the deletion of every of its edge decreases the chromatic index. Except in 2002, Song showed that the conjecture is true for a special class of graphs $G$ with $\Delta(G)\ge \frac{n}{2}$, no other progress on this conjecture had been made. In this paper, we confirm the conjecture for graphs $G$ with $\Delta(G) \ge 0.75n$.


title: Fast Quantized Average Consensus over Static and Dynamic Directed Graphs

abstract: In this paper we study the distributed average consensus problem in multi-agent systems with directed communication links that are subject to quantized information flow. Specifically, we present and analyze a distributed averaging algorithm which operates exclusively with quantized values (i.e., the information stored, processed and exchanged between neighboring agents is subject to deterministic uniform quantization) and relies on event-driven updates (e.g., to reduce energy consumption, communication bandwidth, network congestion, and/or processor usage). The main idea of the proposed algorithm is that each node (i) models its initial state as two quantized fractions which have numerators equal to the node's initial state and denominators equal to one, and (ii) transmits one fraction randomly while it keeps the other stored. Then, every time it receives one or more fractions, it averages their numerators with the numerator of the fraction it stored, and then transmits them to randomly selected out-neighbors. We characterize the properties of the proposed distributed algorithm and show that its execution, on any static and strongly connected digraph, allows each agent to reach in finite time a fixed state that is equal (within one quantisation level) to the average of the initial states. We extend the operation of the algorithm to achieve finite-time convergence in the presence of a dynamic directed communication topology subject to some connectivity conditions. Finally, we provide examples to illustrate the operation, performance, and potential advantages of the proposed algorithm. We compare against state-of-the-art quantized average consensus algorithms and show that our algorithm's convergence speed significantly outperforms most existing protocols.


title: PCOR: Private Contextual Outlier Release via Differentially Private Search

abstract: Outlier detection plays a significant role in various real world applications such as intrusion, malfunction, and fraud detection. Traditionally, outlier detection techniques are applied to find outliers in the context of the whole dataset. However, this practice neglects contextual outliers, that are not outliers in the whole dataset but in some specific neighborhoods. Contextual outliers are particularly important in data exploration and targeted anomaly explanation and diagnosis. In these scenarios, the data owner computes the following information: i) The attributes that contribute to the abnormality of an outlier (metric), ii) Contextual description of the outlier's neighborhoods (context), and iii) The utility score of the context, e.g. its strength in showing the outlier's significance, or in relation to a particular explanation for the outlier. However, revealing the outlier's context leaks information about the other individuals in the population as well, violating their privacy. We address the issue of population privacy violations in this paper, and propose a solution for the two main challenges. In this setting, the data owner is required to release a valid context for the queried record, i.e. a context in which the record is an outlier. Hence, the first major challenge is that the privacy technique must preserve the validity of the context for each record. We propose techniques to protect the privacy of individuals through a relaxed notion of differential privacy to satisfy this requirement. The second major challenge is applying the proposed techniques efficiently, as they impose intensive computation to the base algorithm. To overcome this challenge, we propose a graph structure to map the contexts to, and introduce differentially private graph search algorithms as efficient solutions for the computation problem caused by differential privacy techniques.


title: Learning to Play Soccer From Scratch: Sample-Efficient Emergent Coordination through Curriculum-Learning and Competition

abstract: This work proposes a scheme that allows learning complex multi-agent behaviors in a sample efficient manner, applied to 2v2 soccer. The problem is formulated as a Markov game, and solved using deep reinforcement learning. We propose a basic multi-agent extension of TD3 for learning the policy of each player, in a decentralized manner. To ease learning, the task of 2v2 soccer is divided in three stages: 1v0, 1v1 and 2v2. The process of learning in multi-agent stages (1v1 and 2v2) uses agents trained on a previous stage as fixed opponents. In addition, we propose using experience sharing, a method that shares experience from a fixed opponent, trained in a previous stage, for training the agent currently learning, and a form of frame-skipping, to raise performance significantly. Our results show that high quality soccer play can be obtained with our approach in just under 40M interactions. A summarized video of the resulting game play can be found in https://youtu.be/f25l1j1U9RM.


title: Non-Gaussian mechanical motion via single and multi-phonon subtraction from a thermal state

abstract: Quantum optical measurement techniques offer a rich avenue for quantum control of mechanical oscillators via cavity optomechanics. In particular, a powerful yet little explored combination utilizes optical measurements to perform heralded non-Gaussian mechanical state preparation followed by tomography to determine the mechanical phase-space distribution. Here, we experimentally perform heralded single- and multi-phonon subtraction via photon counting to a laser-cooled mechanical thermal state with a Brillouin optomechanical system at room temperature, and use optical heterodyne detection to measure the $s$-parameterized Wigner distribution of the non-Gaussian mechanical states generated. The techniques developed here advance the state-of-the-art for optics-based tomography of mechanical states and will be useful for a broad range of applied and fundamental studies that utilize mechanical quantum-state engineering and tomography.


title: Unbiased approximation of posteriors via coupled particle Markov chain Monte Carlo

abstract: Markov chain Monte Carlo (MCMC) is a powerful methodology for the approximation of posterior distributions. However, the iterative nature of MCMC does not naturally facilitate its use with modern highly parallel computation on HPC and cloud environments. Another concern is the identification of the bias and Monte Carlo error of produced averages. The above have prompted the recent development of fully ('embarrassingly') parallel unbiased Monte Carlo methodology based on coupling of MCMC algorithms. A caveat is that formulation of effective coupling is typically not trivial and requires model-specific technical effort. We propose coupling of MCMC chains deriving from sequential Monte Carlo (SMC) by considering adaptive SMC methods in combination with recent advances in unbiased estimation for state-space models. Coupling is then achieved at the SMC level and is, in principle, not problem-specific. The resulting methodology enjoys desirable theoretical properties. A central motivation is to extend unbiased MCMC to more challenging targets compared to the ones typically considered in the relevant literature. We illustrate the effectiveness of the algorithm via application to two complex statistical models: (i) horseshoe regression; (ii) Gaussian graphical models.


title: $q$-Moment Estimates for the Singular $p$-Laplace Equation and Applications

abstract: We provide $q$-moment estimates on annuli for weak solutions of the singular $p$-Laplace equation where $p$ and $q$ are conjugates. We derive $q$-uniform integrability for some critical parameter range. As a application, we derive a mass conservation as well as a weak convergence result for a larger critical parameter range. Concerning the latter point, we further provide a rate of convergence of order $t^{q-1}$ of the solution in the $q$-Wasserstein distance.


title: All-Polarization Maintaining Single-Longitudinal-Mode Fiber Laser with Ultra-High OSNR, Sub-kHz Linewidth and Extremely High Stability

abstract: An all-polarization maintaining (PM) single-longitudinal-mode (SLM) erbium-doped fiber laser (EDFL) with ultra-high optical signal-to-noise ratio (OSNR), ultra-narrow linewidth and extremely high stability is proposed and experimentally demonstrated. A double-ring passive subring resonator (DR-PSR) composed of two single-coupler fiber rings and a length of unpumped EDF-based saturable absorber filter is designed and employed in the EDFL to serve as the efficient SLM selecting element to guarantee SLM lasing with excellent output performance. The all-PM structure enables the proposed EDFL to present strong ability to resist the environment disturbance. At the pump power of 100 mW, we obtain an SLM EDFL with an ultra-high OSNR of 83 dB and an ultra-narrow linewidth of 459 Hz. For the SLM operation, the all-PM EDFL processes outstanding stability performance of both the wavelength lasing and the output power. The maximum fluctuations of the center wavelength and output power are 0.012 nm and 0.01 dB.


title: Diagnosis of information scrambling from Hamiltonian evolution under decoherence

abstract: We apply a quantum teleportation protocol based on the Hayden-Preskill thought experiment to quantify how scrambling a given quantum evolution is. It has an advantage over the direct measurement of out-of-time ordered correlators when used to diagnose the information scrambling in the presence of decoherence effects stemming from a noisy quantum device. We demonstrate the protocol by applying it to two physical systems: Ising spin chain and SU(2) lattice Yang-Mills theory. To this end, we numerically simulate the time evolution of the two theories in the Hamiltonian formalism. The lattice Yang-Mills theory is implemented with a suitable truncation of Hilbert space on the basis of the Kogut-Susskind formalism. On a two-leg ladder geometry and with the lowest nontrivial spin representations, it can be mapped to a spin chain, which we call it Yang-Mills-Ising model and is also directly applicable to future digital quantum simulations. We find that the Yang-Mills-Ising model shows the signal of information scrambling at late times.

categories: quant-ph cond-mat.quant-gas cond-mat.stat-mech hep-lat hep-th


title: An Introduction to Deep Generative Modeling

abstract: Deep generative models (DGM) are neural networks with many hidden layers trained to approximate complicated, high-dimensional probability distributions using a large number of samples. When trained successfully, we can use the DGMs to estimate the likelihood of each observation and to create new samples from the underlying distribution. Developing DGMs has become one of the most hotly researched fields in artificial intelligence in recent years. The literature on DGMs has become vast and is growing rapidly. Some advances have even reached the public sphere, for example, the recent successes in generating realistic-looking images, voices, or movies; so-called deep fakes. Despite these successes, several mathematical and practical issues limit the broader use of DGMs: given a specific dataset, it remains challenging to design and train a DGM and even more challenging to find out why a particular model is or is not effective. To help advance the theoretical understanding of DGMs, we introduce DGMs and provide a concise mathematical framework for modeling the three most popular approaches: normalizing flows (NF), variational autoencoders (VAE), and generative adversarial networks (GAN). We illustrate the advantages and disadvantages of these basic approaches using numerical experiments. Our goal is to enable and motivate the reader to contribute to this proliferating research area. Our presentation also emphasizes relations between generative modeling and optimal transport.


title: Classification of resonances and pairing effects on $NA$-scattering within the HFB framework

abstract: We analyze the properties of the scattering solutions obtained as the pole of the S- and K-matrix with the help of the Jost function framework and the Strum-Liouville theory within the Hartree-Fock-Bogoliubov(HFB) framework, and clarify the scattering solutions which can be defined as the physical state. We found that there are three types of the resonances; "{\it shape resonance}", "{\it particle-type}" and "{\it hole-type quasiparticle resonances}", and another two types of solutions are given as the independent S-matrix and K-matrix poles. The shape resonance is formed by the Hartree-Fock(HF) mean field potential, is not affected by the pairing correlation so much. The particle-type and hole-type quasiparticle resonances originate from the particle and hole states by the configuration mixing effect by pairing. All of resonance are represented by the S-matrix pole which has the corresponding K-matrix pole. Two other types of solutions are given by the independent S-matrix and K-matrix poles. These poles are formed by the HF mean field potential. The effect of pairing for the independent S-matrix pole is small, but the one for the independent K-matrix pole has the remarkable effect. The independent K-matrix pole destroys the quasiparticle resonance as it approaches to the resonance by the pairing effect. The wave function of all resonances have the characteristic structure of the metastable property. However, the metastable structure of the wave function of the quasiparticle resonance can be broken by the independent standing wave solution or the Fano effect.


title: The Hal\'asz-Sz\'ekely Barycenter

abstract: We introduce a notion of barycenter of a probability measure related to the symmetric mean of a collection of nonnegative real numbers. Our definition is inspired by the work of Hal\'asz and Sz\'ekely, who in 1976 proved a law of large numbers for symmetric means. We study analytic properties of this Hal\'asz-Sz\'ekely barycenter. We establish fundamental inequalities that relate the symmetric mean of a list of nonnegative real numbers with the barycenter of the measure uniformly supported on these points. As consequence, we go on to establish an ergodic theorem stating that the symmetric means of a sequence of dynamical observations converges to the Hal\'asz-Sz\'ekely barycenter of the corresponding distribution.


title: A Simplified Multifractal Model for Self-Similar Traffic Flows in High-Speed Computer Networks Revisited

abstract: In the context of the simulations carried out using a simplified multifractal model that is proposed to give an explanation to the locality phenomenon that appears in the estimation of the Hurst exponent in the second-order stationary series that represent the self-similar traffic flows in high-speed computer networks, its formulation is perfected to reduce the variability in the singularity limits and it is demonstrated through by its wavelet variant that this modification leads to a higher resolution in the interval of interest under study.


title: Molecular machines for quantum error correction

abstract: Inspired by biological molecular machines we explore the idea of an active quantum robot whose purpose is delaying decoherence. A conceptual model capable of partially protecting arbitrary logical qubit states against single physical qubit errors is presented. Implementation of an instance of that model - the entanglement qubot - is proposed using laser-dressed Rydberg atoms. Dynamics of the system is studied using stochastic wavefunction methods.


title: Near-zero Downtime Recovery from Transient-error-induced Crashes

abstract: Due to the system scaling, transient errors caused by external noises, e.g., heat fluxes and particle strikes, have become a growing concern for the current and upcoming extreme-scale high-performance-computing (HPC) systems. However, since such errors are still quite rare as compared to no-fault cases, desirable solutions call for low/no-overhead systems that do not compromise the performance under no-fault conditions and also allow very fast fault recovery to minimize downtime. In this paper, we present IterPro, a light-weight compiler-assisted resilience technique to quickly and accurately recover processes from transient-error-induced crashes. IterPro repairs the corrupted process states on-the-fly upon occurrences of errors, enabling applications to continue their executions instead of being terminated. IterPro also exploits side effects introduced by induction variable based code optimization techniques to improve its recovery capability. To this end, two new code transformation passes are introduced to expose the side effects for resilience purposes. We evaluated IterPro with 4 scientific workloads as well as the NPB benchmarks suite. During their normal execution, IterPro incurs almost zero runtime overhead and a small, fixed 27MB memory overhead. Meanwhile, IterPro can recover on an average 83.55% of crash-causing errors within dozens of milliseconds with negligible downtime. With such an effective recovery mechanism, IterPro could tremendously mitigate the overheads and resource requirements of the resilience subsystem in future extreme-scale systems.


title: All longest cycles intersect in partial 3-trees

abstract: We show that all longest cycles intersect in 2-connected partial 3-trees.


title: Iterative Shrinking for Referring Expression Grounding Using Deep Reinforcement Learning

abstract: In this paper, we are tackling the proposal-free referring expression grounding task, aiming at localizing the target object according to a query sentence, without relying on off-the-shelf object proposals. Existing proposal-free methods employ a query-image matching branch to select the highest-score point in the image feature map as the target box center, with its width and height predicted by another branch. Such methods, however, fail to utilize the contextual relation between the target and reference objects, and lack interpretability on its reasoning procedure. To solve these problems, we propose an iterative shrinking mechanism to localize the target, where the shrinking direction is decided by a reinforcement learning agent, with all contents within the current image patch comprehensively considered. Beside, the sequential shrinking process enables to demonstrate the reasoning about how to iteratively find the target. Experiments show that the proposed method boosts the accuracy by 4.32% against the previous state-of-the-art (SOTA) method on the RefCOCOg dataset, where query sentences are long and complex, with many targets referred by other reference objects.


title: High-dimensional Quantum Dynamics Study on Excitation-Specific Surface Scattering including Lattice Effects of a Five-Atoms Surface Cell

abstract: In this work high-dimensional (21D) quantum dynamics calculations on mode-specific surface scattering of a carbon monoxide molecule on a copper (100) surface with lattice effects of a five-atom surface cell are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. We employ a surface model in which five surface atoms near the impact site are treated as fully flexible quantum particles while all other more distant atoms are kept at fixed locations. To efficiently perform the 21D ML-MCTDH wavepacket propagation, the potential energy surface is transferred to canonical polyadic decomposition form with the aid of a Monte Carlo based method. Excitation-specific sticking probabilities of CO on Cu(100) are computed and lattice effects caused by the flexible surface atoms are demonstrated by comparison with sticking probabilities computed for a rigid surface. The dependence of the sticking probability of the initial state of the system is studied, and it is found that the sticking probability is reduced, when the surface atom on the impact site is initially vibrationally excited.


title: Urban Epidemic Hazard Index for Chinese Cities: Why Did Small Cities Become Epidemic Hotspots?

abstract: Multiple small- to middle-scale cities, mostly located in northern China, became epidemic hotspots during the second wave of the spread of COVID-19 in early 2021. Despite qualitative discussions of potential social-economic causes, it remains unclear how this pattern could be accounted for from a quantitative approach. Through the development of an urban epidemic hazard index (EpiRank), we came up with a mathematical explanation for this phenomenon. The index is constructed from epidemic simulations on a multi-layer transportation network model on top of local SEIR transmission dynamics, which characterizes intra- and inter-city compartment population flow with a detailed mathematical description. Essentially, we argue that these highlighted cities possess greater epidemic hazards due to the combined effect of large regional population and small inter-city transportation. The proposed index, dynamic and applicable to different epidemic settings, could be a useful indicator for the risk assessment and response planning of urban epidemic hazards in China; the model framework is modularized and can be adapted for other nations without much difficulty.

categories: physics.soc-ph cs.SI econ.GN math.DS q-fin.EC


title: Rotation curves with the multistate Scalar Field Dark Matter model

abstract: We use the concept of co-added rotation curves of Salucci et al. to investigate the properties of axi-symmetric multistate Scalar Field Dark Matter halos in low surface brightness galaxies and dwarf disc galaxies. We fit their rotation curves in two-state configurations and we find that all of these can be well fitted with a particle mass $\mu \sim (10^{-23} - 10^{-24})\rm{eV}/c^2$. Comparing our results with the standard cosmological model, the well-known $\Lambda$-cold dark matter, by using the Bayesian information criterion and the Akaike information criterion, we found that our two-state model seemed to be preferred.


title: Exponential Modalities and Complementarity (extended abstract)

abstract: The exponential modalities of linear logic have been used by various authors to model infinite-dimensional quantum systems. This paper explains how these modalities can also give rise to the complementarity principle of quantum mechanics. The paper uses a formulation of quantum systems based on dagger-linear logic, whose categorical semantics lies in mixed unitary categories, and a formulation of measurement therein. The main result exhibits a complementary system as the result of measurements on free exponential modalities. Recalling that, in linear logic, exponential modalities have two distinct but dual components, ! and ?, this shows how these components under measurement become "compacted" into the usual notion of complementary Frobenius algebras from categorical quantum mechanics.


title: Swell induced vibrations of a thickening ice shelf over a shoaling seabed

abstract: A solution method is developed for a model of ice shelf vibrations in response to ocean waves, in which the ice shelf thickness and seabed beneath the ice shelf vary over distance, and the ice shelf/sub--ice--shelf cavity are connected to the open ocean. The method combines a decomposition of the ice shelf motion into free modes of vibration, a finite element method for the cavity water motion, and a non-local operator to connect to the open ocean. An investigation is conducted into the effects of ice shelf thickening, seabed shoaling and the grounding-line conditions on ice shelf vibrations, induced by regular incident waves in the swell regime. Further, results are given for ice shelf vibrations in response to irregular incident waves, and ocean-to-ice-shelf transfer functions are derived. The findings add to evidence that ice shelves experience appreciable vibrations in response to swell, and that ice shelf thickening and seabed shoaling can have a considerable influence on predictions of how ice shelves respond to ocean waves.


title: Generative Transition Mechanism to Image-to-Image Translation via Encoded Transformation

abstract: In this paper, we revisit the Image-to-Image (I2I) translation problem with transition consistency, namely the consistency defined on the conditional data mapping between each data pairs. Explicitly parameterizing each data mappings with a transition variable $t$, i.e., $x \overset{t(x,y)}{\mapsto}y$, we discover that existing I2I translation models mainly focus on maintaining consistency on results, e.g., image reconstruction or attribute prediction, named result consistency in our paper. This restricts their generalization ability to generate satisfactory results with unseen transitions in the test phase. Consequently, we propose to enforce both result consistency and transition consistency for I2I translation, to benefit the problem with a closer consistency between the input and output. To benefit the generalization ability of the translation model, we propose transition encoding to facilitate explicit regularization of these two {kinds} of consistencies on unseen transitions. We further generalize such explicitly regularized consistencies to distribution-level, thus facilitating a generalized overall consistency for I2I translation problems. With the above design, our proposed model, named Transition Encoding GAN (TEGAN), can poss superb generalization ability to generate realistic and semantically consistent translation results with unseen transitions in the test phase. It also provides a unified understanding of the existing GAN-based I2I transition models with our explicitly modeling of the data mapping, i.e., transition. Experiments on four different I2I translation tasks demonstrate the efficacy and generality of TEGAN.


title: Efficient Topology Design Algorithms for Power Grid Stability

abstract: The dynamic response of power grids to small disturbances influences their overall stability. This paper examines the effect of network topology on the linearized time-invariant dynamics of electric power systems. The proposed framework utilizes ${\cal H}_2$-norm based stability metrics to study the optimal placement of lines on existing networks as well as the topology design of new networks. The design task is first posed as an NP-hard mixed-integer nonlinear program (MINLP) that is exactly reformulated as a mixed-integer linear program (MILP) using McCormick linearization. To improve computation time, graph-theoretic properties are exploited to derive valid inequalities (cuts) and tighten bounds on the continuous optimization variables. Moreover, a cutting plane generation procedure is put forth that is able to interject the MILP solver and augment additional constraints to the problem on-the-fly. The efficacy of our approach in designing optimal grid topologies is demonstrated through numerical tests on the IEEE 39-bus network.


title: An efficient algorithm for deciding vanishing of Schubert polynomial coefficients

abstract: Schubert polynomials form a basis of all polynomials and appear in the study of cohomology rings of flag manifolds. The vanishing problem for Schubert polynomials asks if a coefficient of a Schubert polynomial is zero. We give a tableau criterion to solve this problem, from which we deduce the first polynomial time algorithm. These results are obtained from new characterizations of the Schubitope, a generalization of the permutahedron defined for any subset of the n x n grid. In contrast, we show that computing these coefficients explicitly is #P-complete.


title: A Learning-Based Computational Impact Time Guidance

abstract: This paper investigates the problem of impact-time-control and proposes a learning-based computational guidance algorithm to solve this problem. The proposed guidance algorithm is developed based on a general prediction-correction concept: the exact time-to-go under proportional navigation guidance with realistic aerodynamic characteristics is estimated by a deep neural network and a biased command to nullify the impact time error is developed by utilizing the emerging reinforcement learning techniques. The deep neural network is augmented into the reinforcement learning block to resolve the issue of sparse reward that has been observed in typical reinforcement learning formulation. Extensive numerical simulations are conducted to support the proposed algorithm.


title: An Identity for Expectations and Characteristic Function of Matrix Variate Skew-normalDistribution with Applications to Associated Stochastic Orderings

abstract: We establish an identity for E f (Y) -E f (X), when X and Y both have matrix variateskew-normal distributions and the function f fulfills some weak conditions. Thecharacteristic function of matrix variate skew normal distribution is then derived. Finally,we make use of it to derive some necessary and sucient conditions for the comparisonof matrix variate skew-normal distributions under six di erent orders, such as usualstochastic order, convex order, increasing convex order, upper orthant order,directionally convex order and supermodular order.


title: Continuous body 3-D reconstruction of limbless animals

abstract: Limbless animals such as snakes, limbless lizards, worms, eels, and lampreys move their slender, long bodies in three dimensions to traverse diverse environments. Accurately quantifying their continuous body's 3-D shape and motion is important for understanding body-environment interactions in complex terrain, but this is difficult to achieve (especially for local orientation and rotation). Here, we describe an interpolation method to quantify continuous body 3-D position and orientation. We simplify the body as an elastic rod and apply a backbone optimization method to interpolate continuous body shape between end constraints imposed by tracked markers. Despite over-simplifying the biomechanics, our method achieves a higher interpolation accuracy (~50% error) in both 3-D position and orientation compared with the widely-used cubic B-spline interpolation method. Beyond snakes traversing large obstacles as demonstrated, our method applies to other long, slender, limbless animals and continuum robots. We provide codes and demo files for easy application of our method.


title: Dynamic Pricing and Learning under the Bass Model

abstract: We consider a novel formulation of the dynamic pricing and demand learning problem, where the evolution of demand in response to posted prices is governed by a stochastic variant of the popular Bass model with parameters $\alpha, \beta$ that are linked to the so-called "innovation" and "imitation" effects. Unlike the more commonly used i.i.d. and contextual demand models, in this model the posted price not only affects the demand and the revenue in the current round but also the future evolution of demand, and hence the fraction of potential market size $m$ that can be ultimately captured. In this paper, we consider the more challenging incomplete information problem where dynamic pricing is applied in conjunction with learning the unknown parameters, with the objective of optimizing the cumulative revenues over a given selling horizon of length $T$. Equivalently, the goal is to minimize the regret which measures the revenue loss of the algorithm relative to the optimal expected revenue achievable under the stochastic Bass model with market size $m$ and time horizon $T$. Our main contribution is the development of an algorithm that satisfies a high probability regret guarantee of order $\tilde O(m^{2/3})$; where the market size $m$ is known a priori. Moreover, we show that no algorithm can incur smaller order of loss by deriving a matching lower bound. Unlike most regret analysis results, in the present problem the market size $m$ is the fundamental driver of the complexity; our lower bound in fact, indicates that for any fixed $\alpha, \beta$, most non-trivial instances of the problem have constant $T$ and large $m$. We believe that this insight sets the problem of dynamic pricing under the Bass model apart from the typical i.i.d. setting and multi-armed bandit based models for dynamic pricing, which typically focus only on the asymptotics with respect to time horizon $T$.


title: Virtual Reality Sickness Mitigation Methods: A Comparative Study in a Racing Game

abstract: Using virtual reality (VR) head-mounted displays (HMDs) can induce VR sickness. VR sickness can cause strong discomfort, decrease users' presence and enjoyment, especially in games, shorten the duration of the VR experience, and can even pose health risks. Previous research has explored different VR sickness mitigation methods by adding visual effects or elements. Field of View (FOV) reduction, Depth of Field (DOF) blurring, and adding a rest frame into the virtual environment are examples of such methods. Although useful in some cases, they might result in information loss. This research is the first to compare VR sickness, presence, workload to complete a search task, and information loss of these three VR sickness mitigation methods in a racing game with two levels of control. To do this, we conducted a mixed factorial user study (N = 32) with degree of control as the between-subjects factor and the VR sickness mitigation techniques as the within-subjects factor. Participants were required to find targets with three difficulty levels while steering or not steering a car in a virtual environment. Our results show that there are no significant differences in VR sickness, presence and workload among these techniques under two levels of control in our VR racing game. We also found that changing FOV dynamically or using DOF blur effects would result in information loss while adding a target reticule as a rest frame would not.


title: Multivariate tail covariance for generalized skew-elliptical distributions

abstract: In this paper, the multivariate tail covariance (MTCov) for generalized skew-elliptical distributions is considered. Some special cases for this distribution, such as generalized skew-normal, generalized skew student-t, generalized skew-logistic and generalized skew-Laplace distributions, are also considered. In order to test the theoretical feasibility of our results, the MTCov for skewed and non skewed normal distributions are computed and compared. Finally, we give a special formula of the MTCov for generalized skew-elliptical distributions.


title: Rainbow Independent Sets in Cycles

abstract: For a given class ${\cal C}$ of graphs and given integers $m \le n$, let $f_{\cal C}(n,m)$ be the minimal number $k$ such that every $k$ independent $n$-sets in any graph belonging to ${\cal C}$ have a (possibly partial) rainbow independent $m$-set. In this paper, we consider the case ${\cal C}=\{C_{2s+1}\}$ and show that $f_{C_{2s+1}}(s, s) = s$. Our result is a special case of the conjecture (Conjecture 2.9) proposed by Aharoni et al in \cite{Aharoni}.


title: Persistence and Smooth Dependence on Parameters of Periodic Orbits in Functional Differential Equations Close to an ODE or an Evolutionary PDE

abstract: We consider functional differential equations(FDEs) which are perturbations of smooth ordinary differential equations(ODEs). The FDE can involve multiple state-dependent delays or distributed delays (forward or backward). We show that, under some mild assumptions, if the ODE has a nondegenerate periodic orbit, then the FDE has a smooth periodic orbit. Moreover, we get smooth dependence of the periodic orbit and its frequency on parameters with high regularity. The result also applies to FDEs which are perturbations of some evolutionary partial differential equations(PDEs). The proof consists in solving functional equations satisfied by the parameterization of the periodic orbit and the frequency using a fixed point approach. We do not need to consider the smoothness of the evolution or even the phase space of the FDEs.


title: A New Metric on Symmetric Group and Applications to Block Permutation Codes

abstract: Permutation codes have received a great attention due to various applications. For different applications, one needs permutation codes under different metrics. The generalized Cayley metric was introduced by Chee and Vu [4] and this metric includes several other metrics as special cases. However, the generalized Cayley metric is not easily computable in general. Therefore the block permutation metric was introduced by Yang et al. [22] as the generalized Cayley metric and the block permutation metric have the same magnitude. However, the block permutation metric lacks the symmetry property which restricts more advanced algebraic tools to be involved. In this paper, by introducing a novel metric closely related to the block permutation metric, we build a bridge between some advanced algebraic methods and codes in the block permutation metric. More specifically, based on some techniques from algebraic function fields originated in [19], we give an algebraic-geometric construction of codes in the novel metric with reasonably good parameters. By observing a trivial relation between the novel metric and block permutation metric, we then produce non-systematic codes in block permutation metric that improve all known results given in [21, 22]. More importantly, based on our non-systematic codes, we provide an explicit and systematic construction of codes in block permutation metric which improves the systematic result shown in [22]. In the end, we demonstrate that our codes in the novel metric itself have reasonably good parameters by showing that our construction beats the corresponding Gilbert-Varshamov bound.


title: Target-Normal Single Spin Asymmetries Measured with Positrons

abstract: Two-photon exchange and the larger class of hadronic box diagrams are difficult to calculate without a large degree of model-dependence. At the same time, these processes are significant radiative corrections in parity-violating electron scattering, in neutron decay, and may even be responsible for the proton's form factor ratio discrepancy. New kinds of experimental data are needed to help constrain models and guide future box-diagram calculations. The target-normal single spin asymmetry, $A_n$, formed with an unpolarized beam scattering from a target that is polarized normal to the scattering plane, is sensitive to the imaginary part of the two-photon exchange amplitude, and can provide a valuable constraint. A measurement with both electrons and positrons can reduce sources of experimental error, and distinguish between the effects of two-photon exchange and those of time-reversal symmetry violation. This article describes a proposed experiment in Hall A, using the new Super Big-Bite Spectrometer that can cover a momentum transfer range in the critical zone of uncertainty between where hadronic calculations and those based on partonic degrees of freedom are expected to be accurate.


title: Controlling fingering instabilities in Hele-Shaw flows in the presence of wetting film effects

abstract: In this paper, the interfacial motion between two immiscible viscous fluids in the confined geometry of a Hele-Shaw cell is studied. We consider the influence of a thin wetting film trailing behind the displaced fluid, which dynamically affects the pressure drop at the fluid-fluid interface by introducing a nonlinear dependence on the interfacial velocity. In this framework, two cases of interest are analyzed: The injection-driven flow (expanding evolution), and the lifting plate flow (shrinking evolution). In particular, we investigate the possibility of controlling the development of fingering instabilities in these two different Hele-Shaw setups when wetting effects are taken into account. By employing linear stability theory, we find the proper time-dependent injection rate $Q(t)$ and the time-dependent lifting speed ${\dot b}(t)$ required to control the number of emerging fingers during the expanding and shrinking evolution, respectively. Our results indicate that the consideration of wetting leads to an increase in the magnitude of $Q(t)$ [and ${\dot b}(t)$] in comparison to the non-wetting strategy. Moreover, a spectrally accurate boundary integral approach is utilized to examine the validity and effectiveness of the controlling protocols at the fully nonlinear regime of the dynamics and confirms that the proposed injection and lifting schemes are feasible strategies to prescribe the morphologies of the resulting patterns in the presence of the wetting film.


title: Queer dual equivalence graphs

abstract: We introduce a new paradigm for proving the Schur $P$-positivity. Generalizing dual equivalence, we give an axiomatic definition for a family of involutions on a set of objects to be a queer dual equivalence, and we prove whenever such a family exists, the fundamental quasisymmetric generating function is Schur $P$-positive. In contrast with shifted dual equivalence, the queer dual equivalence involutions restrict to a dual equivalence when the queer involution is omitted. We highlight the difference between these two generalization with a new application to the product of Schur $P$-functions.


title: Random nanowire networks: Identification of a current-carrying subset of wires using a modified wall follower algorithm

abstract: We mimic random nanowire networks by the homogeneous, isotropic, and random deposition of conductive zero-width sticks onto an insulating substrate. The number density (the number of objects per unit area of the surface) of these sticks is supposed to exceed the percolation threshold, i.e., the system under consideration is a conductor. To identify any current-carrying part (the backbone) of the percolation cluster, we have proposed and implemented a modification of the well-known wall follower algorithm -- one type of maze solving algorithm. The advantage of the modified algorithm is its identification of the whole backbone without visiting all the edges. The complexity of the algorithm depends significantly on the structure of the graph and varies from $O\left(\sqrt{N_\text{V}}\right)$ to $\Theta(N_\text{V})$. The algorithm has been applied to backbone identification in networks with different number densities of conducting sticks. We have found that (i) for number densities of sticks above the percolation threshold, the strength of the percolation cluster quickly approaches unity as the number density of the sticks increases; (ii) simultaneously, the percolation cluster becomes identical to its backbone plus simplest dead ends, i.e., edges that are incident to vertices of degree 1. This behavior is consistent with the presented analytical evaluations.


title: Rapid-response radio observations of short GRB 181123B with the Australia Telescope Compact Array

abstract: We introduce the Australia Telescope Compact Array (ATCA) rapid-response mode by presenting the first successful trigger on the short-duration gamma-ray burst (GRB) 181123B. Early-time radio observations of short GRBs may provide vital insights into the radio afterglow properties of Advanced LIGO- and Virgo-detected gravitational wave events, which will in turn inform follow-up strategies to search for counterparts within their large positional uncertainties. The ATCA was on target within 12.6 hr post-burst, when the source had risen above the horizon. While no radio afterglow was detected during the 8.3 hr observation, we obtained force-fitted flux densities of $7 \pm 12$ and $15 \pm 11~\mu$Jy at 5.5 and 9 GHz, respectively. Afterglow modelling of GRB 181123B showed that the addition of the ATCA force-fitted radio flux densities to the Swift X-ray Telescope detections provided more stringent constraints on the fraction of thermal energy in the electrons (log$\epsilon_e = -0.75^{+0.39}_{-0.40}$ rather than log$\epsilon_e = -1.13^{+0.82}_{-1.2}$ derived without the inclusion of the ATCA values), which is consistent with the range of typical $\epsilon_e$ derived from GRB afterglow modelling. This allowed us to predict that the forward shock may have peaked in the radio band $\sim10$ days post-burst, producing detectable radio emission $\gtrsim3-4$ days post-burst. Overall, we demonstrate the potential for extremely rapid radio follow-up of transients and the importance of triggered radio observations for constraining GRB blast wave properties, regardless of whether there is a detection, via the inclusion of force-fitted radio flux densities in afterglow modelling efforts.


title: An atomic perspective on the serpentine-chlorite solid-state transformation

abstract: Serpentine minerals are important components of metamorphic rocks and promising geo-materials for nanotechnology. Lizardite, the most abundant serpentine mineral, can be transformed into chlorite during metamorphism. This intriguing phase transformation should affect the deformation behavior during aseismic creep and slow slip at the base of the subduction zone, but has not been understood structurally and chemically at the atomic scale. Here we visualized cations and oxygen atoms using the state-of-the-art low-dose scanning transmission electron microscopy and found that restructuring mainly involves the synergistic migration of tetrahedral cations and oxygen anions, coupled with the migration of octahedral trivalent cations into the brucite-like interlayer. Further, we show that different serpentine polytypes result in distinct regular interstratifications of serpentine and chlorite. Our results clarify the long-standing puzzle of how solid-state layer silicate transformations occur and lead to long-period ordered structures.


title: Combing a double helix

abstract: Combing hair involves brushing away the topological tangles in a collective curl. Using a combination of experiment and computation, we study this problem that naturally links topology, geometry and mechanics. Observations show that the dominant interactions in hair are those of a two-body nature, corresponding to a braided homochiral double helix. Using this minimal model, we study the detangling of an elastic double helix via a single stiff tine that moves along it, leaving two untangled filaments in its wake. Our results quantify how the forces of detangling are correlated with the magnitude and spatial extent of the link density, a topological quantity, that propagates ahead of the tine. This in turn provides a measure of the maximum characteristic length of a single combing stroke, and thus the trade-offs between comfort, efficiency and speed of combing in the many-body problem on a head of hair.


title: Testing f(R) gravity with the simulated data of gravitational waves from the Einstein Telescope

abstract: In this paper we analyze the implications of gravitational waves (GWs) as standard sirens on the modified gravity models by using the third-generation gravitational wave detector, i.e., the Einstein Telescope. Two viable models in $f(R)$ theories within the Palatini formalism are considered in our analysis ($f_{1}(\mathcal{R})=\mathcal{R}-\frac{\beta}{\mathcal{R}^{n}}$ and $f_{2}(\mathcal{R})=\mathcal{R}+\alpha\ln{\mathcal{R}}-\beta$), with the combination of simulated GW data and the latest electromagnetic (EM) observational data (including the recently released Pantheon type Ia supernovae sample, the cosmic chronometer data, and baryon acoustic oscillation distance measurements). Our analysis reveals that the standard sirens GWs, which provide an independent and complementary alternative to current experiments, could effectively eliminate the degeneracies among parameters in the two modified gravity models. In addition, we thoroughly investigate the nature of geometrical dark energy in the modified gravity theories with the assistance of $Om(z)$ and statefinder diagnostic analysis. The present analysis makes it clear-cut that the simplest cosmological constant model is still the most preferred by the current data. However, the combination of future naturally improved GW data most recent EM observations will reveal the consistency or acknowledge the tension between the $\Lambda$CDM model and modified gravity theories.


title: Enhancing Medical Image Registration via Appearance Adjustment Networks

abstract: Deformable image registration is fundamental for many medical image analyses. A key obstacle for accurate image registration lies in image appearance variations such as the variations in texture, intensities, and noise. These variations are readily apparent in medical images, especially in brain images where registration is frequently used. Recently, deep learning-based registration methods (DLRs), using deep neural networks, have shown computational efficiency that is several orders of magnitude faster than traditional optimization-based registration methods (ORs). DLRs rely on a globally optimized network that is trained with a set of training samples to achieve faster registration. DLRs tend, however, to disregard the target-pair-specific optimization inherent in ORs and thus have degraded adaptability to variations in testing samples. This limitation is severe for registering medical images with large appearance variations, especially since few existing DLRs explicitly take into account appearance variations. In this study, we propose an Appearance Adjustment Network (AAN) to enhance the adaptability of DLRs to appearance variations. Our AAN, when integrated into a DLR, provides appearance transformations to reduce the appearance variations during registration. In addition, we propose an anatomy-constrained loss function through which our AAN generates anatomy-preserving transformations. Our AAN has been purposely designed to be readily inserted into a wide range of DLRs and can be trained cooperatively in an unsupervised and end-to-end manner. We evaluated our AAN with three state-of-the-art DLRs on three well-established public datasets of 3D brain magnetic resonance imaging (MRI). The results show that our AAN consistently improved existing DLRs and outperformed state-of-the-art ORs on registration accuracy, while adding a fractional computational load to existing DLRs.


title: Universal Undersampled MRI Reconstruction

abstract: Deep neural networks have been extensively studied for undersampled MRI reconstruction. While achieving state-of-the-art performance, they are trained and deployed specifically for one anatomy with limited generalization ability to another anatomy. Rather than building multiple models, a universal model that reconstructs images across different anatomies is highly desirable for efficient deployment and better generalization. Simply mixing images from multiple anatomies for training a single network does not lead to an ideal universal model due to the statistical shift among datasets of various anatomies, the need to retrain from scratch on all datasets with the addition of a new dataset, and the difficulty in dealing with imbalanced sampling when the new dataset is further of a smaller size. In this paper, for the first time, we propose a framework to learn a universal deep neural network for undersampled MRI reconstruction. Specifically, anatomy-specific instance normalization is proposed to compensate for statistical shift and allow easy generalization to new datasets. Moreover, the universal model is trained by distilling knowledge from available independent models to further exploit representations across anatomies. Experimental results show the proposed universal model can reconstruct both brain and knee images with high image quality. Also, it is easy to adapt the trained model to new datasets of smaller size, i.e., abdomen, cardiac and prostate, with little effort and superior performance.


title: Gate Tunable Supercurrent in Josephson Junctions Based on Bi2Te3 Topological Insulator Thin Films

abstract: We report transport measurements on Josephson junctions consisting of Bi2Te3 topological insulator (TI) thin films contacted by superconducting Nb electrodes. For a device with junction length L = 134 nm, the critical supercurrent Ic can be modulated by an electrical gate which tunes the carrier type and density of the TI film. Ic can reach a minimum when the TI is near the charge neutrality regime with the Fermi energy lying close to the Dirac point of the surface state. In the p-type regime the Josephson current can be well described by a short ballistic junction model. In the n-type regime the junction is ballistic at 0.7 K < T < 3.8 K while for T < 0.7 K the diffusive bulk modes emerge and contribute a larger Ic than the ballistic model. We attribute the lack of diffusive bulk modes in the p-type regime to the formation of p-n junctions. Our work provides new clues for search of Majorana zero mode in TI-based superconducting devices.


title: Tidal Effects on the Radial Velocities of V723 Mon: Additional Evidence for a Dark $3\,M_\odot$ Companion

abstract: Jayasinghe et al. (2021) identified a dark $\approx 3\,M_\odot$ companion on a nearly edge-on $\approx 60\,\mathrm{day}$ orbit around the red giant star V723 Monoceros as a black hole candidate in the mass gap. This scenario was shown to explain most of the data presented by Jayasinghe et al. (2021), except for periodic radial velocity (RV) residuals from the circular Keplerian model. Here we show that the RV residuals are explained by orbital phase-dependent distortion of the absorption line profile associated with changing visible fractions of the approaching and receding sides of the red giant star, whose surface is tidally deformed by and rotating synchronously with the dark companion. Our RV model constrains the companion mass $M_\bullet = 2.95\pm0.17\,M_\odot$ and orbital inclination $i=82.9^{+7.0}_{-3.3}\,\mathrm{deg}$ (medians and 68.3% highest density intervals of the marginal posteriors) adopting the radius of the red giant $24.0\pm0.9\,R_\odot$ as constrained from its SED and distance. The analysis provides independent support for the companion mass from ellipsoidal variations and the limits on the companion's luminosity from the absence of eclipses both derived by Jayasinghe et al. (2021). We also show that a common scheme to evaluate the tidal RV signal as the flux-weighted mean of the surface velocity field can significantly underestimate its amplitude for RVs measured with a cross-correlation technique, and present a modified prescription that directly models the distorted line profile and its effects on the measured RVs. The formulation will be useful for estimating the component masses and inclinations in other similar binaries.


title: Sequential Importance Sampling With Corrections For Partially Observed States

abstract: We consider an evolving system for which a sequence of observations is being made, with each observation revealing additional information about current and past states of the system. We suppose each observation is made without error, but does not fully determine the state of the system at the time it is made. Our motivating example is drawn from invasive species biology, where it is common to know the precise location of invasive organisms that have been detected by a surveillance program, but at any time during the program there are invaders that have not been detected. We propose a sequential importance sampling strategy to infer the state of the invasion under a Bayesian model of such a system. The strategy involves simulating multiple alternative states consistent with current knowledge of the system, as revealed by the observations. However, a difficult problem that arises is that observations made at a later time are invariably incompatible with previously simulated states. To solve this problem, we propose a two-step iterative process in which states of the system are alternately simulated in accordance with past observations, then corrected in light of new observations. We identify criteria under which such corrections can be made while maintaining appropriate importance weights.


title: Searching for axion-like particles at future electron-positron colliders

abstract: We investigate the prospects for discovering axion-like particles (ALPs) via a light-by-light (LBL) scattering at two colliders, the future circular collider (FCC-ee) and circular electron-positron collider (CEPC). The promising sensitivities to the effective ALP-photon coupling $g_{a\gamma\gamma}$ are obtained. Our numerical results show that the FCC-ee and CEPC might be more sensitive to the ALPs with mass 2 GeV $\sim$ 10 GeV than the LHC and CLIC.


title: Outdoor sub-THz Position Location and Tracking using Field Measurements at 142 GHz

abstract: Future sub-THz cellular deployments may be utilized to complement the coverage of the global positioning system (GPS) and provide centimeter-level accuracy. In this work, we use measurement data at 142 GHz to test a map-based position location algorithm in an outdoor urban microcell (UMi) environment. We utilize an extended Kalman filter (EKF) to track the position of the user equipment (UE) along a rectangular track, with the transmitter-receiver separation distances varying from 24.3 m to 52.8 m. The position and velocity of the UE are tracked by the EKF, with measurements of the angle of arrival and time of flight information obtained along an outdoor track, to provide a mean accuracy of 24.8 cm at 142 GHz, over 34 UE locations, using a single base station in line-of-sight and non-line-of-sight.


title: Prediction of 5-year Progression-Free Survival in Advanced Nasopharyngeal Carcinoma with Pretreatment PET/CT using Multi-Modality Deep Learning-based Radiomics

abstract: Objective: Deep Learning-based Radiomics (DLR) has achieved great success in medical image analysis and has been considered a replacement for conventional radiomics that relies on handcrafted features. In this study, we aimed to explore the capability of DLR for the prediction of 5-year Progression-Free Survival (PFS) in Nasopharyngeal Carcinoma (NPC) using pretreatment PET/CT. Methods: A total of 257 patients (170/87 in internal/external cohorts) with advanced NPC (TNM stage III or IVa) were enrolled. We developed an end-to-end multi-modality DLR model, in which a 3D convolutional neural network was optimized to extract deep features from pretreatment PET/CT images and predict the probability of 5-year PFS. TNM stage, as a high-level clinical feature, could be integrated into our DLR model to further improve the prognostic performance. To compare conventional radiomics and DLR, 1456 handcrafted features were extracted, and optimal conventional radiomics methods were selected from 54 cross-combinations of 6 feature selection methods and 9 classification methods. In addition, risk group stratification was performed with clinical signature, conventional radiomics signature, and DLR signature. Results: Our multi-modality DLR model using both PET and CT achieved higher prognostic performance than the optimal conventional radiomics method. Furthermore, the multi-modality DLR model outperformed single-modality DLR models using only PET or only CT. For risk group stratification, the conventional radiomics signature and DLR signature enabled significant differences between the high- and low-risk patient groups in both internal and external cohorts, while the clinical signature failed in the external cohort. Conclusion: Our study identified potential prognostic tools for survival prediction in advanced NPC, suggesting that DLR could provide complementary values to the current TNM staging.


title: Learning to Find Usages of Library Functions in Optimized Binaries

abstract: Much software, whether beneficent or malevolent, is distributed only as binaries, sans source code. Absent source code, understanding binaries' behavior can be quite challenging, especially when compiled under higher levels of compiler optimization. These optimizations can transform comprehensible, "natural" source constructions into something entirely unrecognizable. Reverse engineering binaries, especially those suspected of being malevolent or guilty of intellectual property theft, are important and time-consuming tasks. There is a great deal of interest in tools to "decompile" binaries back into more natural source code to aid reverse engineering. Decompilation involves several desirable steps, including recreating source-language constructions, variable names, and perhaps even comments. One central step in creating binaries is optimizing function calls, using steps such as inlining. Recovering these (possibly inlined) function calls from optimized binaries is an essential task that most state-of-the-art decompiler tools try to do but do not perform very well. In this paper, we evaluate a supervised learning approach to the problem of recovering optimized function calls. We leverage open-source software and develop an automated labeling scheme to generate a reasonably large dataset of binaries labeled with actual function usages. We augment this large but limited labeled dataset with a pre-training step, which learns the decompiled code statistics from a much larger unlabeled dataset. Thus augmented, our learned labeling model can be combined with an existing decompilation tool, Ghidra, to achieve substantially improved performance in function call recovery, especially at higher levels of optimization.


title: Data augmentation by morphological mixup for solving Raven's Progressive Matrices

abstract: Raven's Progressive Matrices (RPMs) are frequently used in testing human's visual reasoning ability. Recent advances of RPM-like datasets and solution models partially address the challenges of visually understanding the RPM questions and logically reasoning the missing answers. In view of the poor generalization performance due to insufficient samples in RPM datasets, we propose an effective scheme, namely Candidate Answer Morphological Mixup (CAM-Mix). CAM-Mix serves as a data augmentation strategy by gray-scale image morphological mixup, which regularizes various solution methods and overcomes the model overfitting problem. By creating new negative candidate answers semantically similar to the correct answers, a more accurate decision boundary could be defined. By applying the proposed data augmentation method, a significant and consistent performance improvement is achieved on various RPM-like datasets compared with the state-of-the-art models.


title: The pointwise stabilities of piecewise linear finite element method on non-obtuse tetrahedral meshes of nonconvex polyhedra

abstract: Let $\Omega$ be a Lipschitz polyhedral (can be nonconvex) domain in $\mathbb{R}^{3}$, and $V_{h}$ denotes the finite element space of continuous piecewise linear polynomials. On non-obtuse quasi-uniform tetrahedral meshes, we prove that the finite element projection $R_{h}u$ of $u \in H^{1}(\Omega) \cap C(\overline{\Omega})$ (with $R_{h} u$ interpolating $u$ at the boundary nodes) satisfies \begin{align*} \Vert R_{h} u\Vert_{L^{\infty}(\Omega)} \leq C \vert \log h \vert \Vert u\Vert_{L^{\infty}(\Omega)}. \end{align*} If we further assume $u \in W^{1,\infty}(\Omega)$, then \begin{align*} \Vert R_{h} u\Vert_{W^{1, \infty}(\Omega)} \leq C \vert \log h \vert \Vert u\Vert_{W^{1, \infty}(\Omega)}. \end{align*}


title: A review of flywheel energy storage systems: state of the art and opportunities

abstract: Thanks to the unique advantages such as long life cycles, high power density, minimal environmental impact, and high power quality such as fast response and voltage stability, the flywheel/kinetic energy storage system (FESS) is gaining attention recently. There is noticeable progress made in FESS, especially in utility, large-scale deployment for the electrical grid, and renewable energy applications. This paper gives a review of the recent developments in FESS technologies. Due to the highly interdisciplinary nature of FESSs, we survey different design approaches, choices of subsystems, and the effects on performance, cost, and applications. This review focuses on the state of the art of FESS technologies, especially for those who have been commissioned or prototyped. We also highlighted the opportunities and potential directions for the future development of FESS technologies.


title: A Scavenger Hunt for Service Robots

abstract: Creating robots that can perform general-purpose service tasks in a human-populated environment has been a longstanding grand challenge for AI and Robotics research. One particularly valuable skill that is relevant to a wide variety of tasks is the ability to locate and retrieve objects upon request. This paper models this skill as a Scavenger Hunt (SH) game, which we formulate as a variation of the NP-hard stochastic traveling purchaser problem. In this problem, the goal is to find a set of objects as quickly as possible, given probability distributions of where they may be found. We investigate the performance of several solution algorithms for the SH problem, both in simulation and on a real mobile robot. We use Reinforcement Learning (RL) to train an agent to plan a minimal cost path, and show that the RL agent can outperform a range of heuristic algorithms, achieving near optimal performance. In order to stimulate research on this problem, we introduce a publicly available software stack and associated website that enable users to upload scavenger hunts which robots can download, perform, and learn from to continually improve their performance on future hunts.


title: DeepSeagrass Dataset

abstract: We introduce a dataset of seagrass images collected by a biologist snorkelling in Moreton Bay, Queensland, Australia, as described in our publication: arXiv:2009.09924. The images are labelled at the image-level by collecting images of the same morphotype in a folder hierarchy. We also release pre-trained models and training codes for detection and classification of seagrass species at the patch level at https://github.com/csiro-robotics/deepseagrass.


title: Uncertainty-aware Incremental Learning for Multi-organ Segmentation

abstract: Most existing approaches to train a unified multi-organ segmentation model from several single-organ datasets require simultaneously access multiple datasets during training. In the real scenarios, due to privacy and ethics concerns, the training data of the organs of interest may not be publicly available. To this end, we investigate a data-free incremental organ segmentation scenario and propose a novel incremental training framework to solve it. We use the pretrained model instead of its own training data for privacy protection. Specifically, given a pretrained $K$ organ segmentation model and a new single-organ dataset, we train a unified $K+1$ organ segmentation model without accessing any data belonging to the previous training stages. Our approach consists of two parts: the background label alignment strategy and the uncertainty-aware guidance strategy. The first part is used for knowledge transfer from the pretained model to the training model. The second part is used to extract the uncertainty information from the pretrained model to guide the whole knowledge transfer process. By combing these two strategies, more reliable information is extracted from the pretrained model without original training data. Experiments on multiple publicly available pretrained models and a multi-organ dataset MOBA have demonstrated the effectiveness of our framework.


title: Constraints on High-Order Gravitational Theories from Weak Cosmic Censorship

abstract: Weak cosmic censorship conjecture (WCCC) is a basic principle that guarantees the predictability of spacetime and should be valid in any classical theories. One critical scientific question is whether the WCCC can serve as a constraint to the gravitational theories. To explore this question, we perform the first-order Sorce-Wald's gedanken experiments to test the WCCC in the higher-order gravitational theories and find that there exists a destruction condition $S_\text{ext}'(r_h)<0$ for the extremal black holes. To show the power of this condition, we evaluate the constraints given by WCCC in the quadratic and cubic gravitational theories. Our investigation makes an essential step toward applying WCCC to constrain the modified gravitational theories, and opens a new avenue to judge which theory is reasonable.


title: Effect of electron-irradiation on layered quantum materials

abstract: Technological advancement towards the quantum era requires secure communication, quantum computation, and ultra-sensitive sensing capabilities. Layered quantum materials (LQMs) have remarkable optoelectronic and quantum properties that can usher us into the quantum era. Electron microscopy is the tool of choice for measuring these LQMs at an atomic and nanometer scale. On the other hand, electron-irradiation of LQMs can modify various material properties, including the creation of structural defects. We review different types of structural defects, as well as electron elastic- and inelastic-scattering induced processes. Controlled modification of optoelectronic and quantum properties of LQMs using electron-irradiation, including creation of single photon emitters is discussed. Protection of electron-irradiation induced damage of LQMs via encapsulation by other layered materials is encouraged. We finally give insights into challenges and opportunities, including creating novel structures using an electron beam.


title: Physics of Superluminous Supernovae

abstract: Understanding how massive stars die as supernovae is a crucial question in modern astrophysics. Supernovae are powerful stellar explosions and key drivers in the cosmic baryonic cycles by injecting their explosion energy and heavy elements to the interstellar medium that forms new stars. After decades of effort, astrophysicists have built up a stand model for the explosion mechanism of massive stars. However, this model is challenged by new kinds of stellar explosions discovered in the recent transit surveys. In particular, the new population called superluminous supernovae, which are a hundred times brighter than typical supernovae, is revolutionizing our understanding of supernovae. New studies suggest the superluminous supernovae are associated with the unusual demise of very massive stars and their extreme supernovae powered by the radioactive isotopes or compact objects formed after the explosion. Studying these supernovae fills a gap of knowledge between the death of massive stars and their explosions; furthermore, we may apply their intense luminosity to light up the distant universe. This paper aims to provide a timely review of superluminous supernovae physics, focusing on the latest development of their theoretical models.


title: Self-supervised Regularization for Text Classification

abstract: Text classification is a widely studied problem and has broad applications. In many real-world problems, the number of texts for training classification models is limited, which renders these models prone to overfitting. To address this problem, we propose SSL-Reg, a data-dependent regularization approach based on self-supervised learning (SSL). SSL is an unsupervised learning approach which defines auxiliary tasks on input data without using any human-provided labels and learns data representations by solving these auxiliary tasks. In SSL-Reg, a supervised classification task and an unsupervised SSL task are performed simultaneously. The SSL task is unsupervised, which is defined purely on input texts without using any human-provided labels. Training a model using an SSL task can prevent the model from being overfitted to a limited number of class labels in the classification task. Experiments on 17 text classification datasets demonstrate the effectiveness of our proposed method.


title: Stabilized Medical Image Attacks

abstract: Convolutional Neural Networks (CNNs) have advanced existing medical systems for automatic disease diagnosis. However, a threat to these systems arises that adversarial attacks make CNNs vulnerable. Inaccurate diagnosis results make a negative influence on human healthcare. There is a need to investigate potential adversarial attacks to robustify deep medical diagnosis systems. On the other side, there are several modalities of medical images (e.g., CT, fundus, and endoscopic image) of which each type is significantly different from others. It is more challenging to generate adversarial perturbations for different types of medical images. In this paper, we propose an image-based medical adversarial attack method to consistently produce adversarial perturbations on medical images. The objective function of our method consists of a loss deviation term and a loss stabilization term. The loss deviation term increases the divergence between the CNN prediction of an adversarial example and its ground truth label. Meanwhile, the loss stabilization term ensures similar CNN predictions of this example and its smoothed input. From the perspective of the whole iterations for perturbation generation, the proposed loss stabilization term exhaustively searches the perturbation space to smooth the single spot for local optimum escape. We further analyze the KL-divergence of the proposed loss function and find that the loss stabilization term makes the perturbations updated towards a fixed objective spot while deviating from the ground truth. This stabilization ensures the proposed medical attack effective for different types of medical images while producing perturbations in small variance. Experiments on several medical image analysis benchmarks including the recent COVID-19 dataset show the stability of the proposed method.


title: Towards Artefact-based Requirements Engineering for Data-Centric Systems

abstract: Many modern software-intensive systems employ artificial intelligence / machine-learning (AI/ML) components and are, thus, inherently data-centric. The behaviour of such systems depends on typically large amounts of data processed at run-time rendering such non-deterministic systems as complex. This complexity growth affects our understanding on needs and practices in Requirements Engineering (RE). There is, however, still little guidance on how to handle requirements for such systems effectively: What are, for example, typical quality requirements classes? What modelling concepts do we rely on or which levels of abstraction do we need to consider? In fact, how to integrate such concepts into approaches for a more traditional RE still needs profound investigations. In this research preview paper, we report on ongoing efforts to establish an artefact-based RE approach for the development of datacentric systems (DCSs). To this end, we sketch a DCS development process with the newly proposed requirements categories and data-centric artefacts and briefly report on an ongoing investigation of current RE challenges in industry developing data-centric systems.


title: Simultaneous Conjugacy Classes of Finite $p$-groups of rank $\leq 5$

abstract: For a finite group $G$, we consider the problem of counting simultaneous conjugacy classes of $n$-tuples and simultaneous conjugacy classes of commuting $n$-tuples in $G$. Let $\alpha_{G,n}$ denote the number of simultaneous conjugacy classes of $n$-tuples, and $\beta_{G,n}$ the number of simultaneous conjugacy classes of commuting $n$-tuples in $G$. The generating functions $A_G(t) = \sum_{n\geq 0} \alpha_{G,n}t^n,$ and $B_G(t) = \sum_{n\geq 0} \beta_{G,n}t^n$ are rational functions of $t$. This paper concern studied of normalized functions $A_G(t/|G|)$ and $B_G(t/|G|)$ for finite $p$-groups of rank at most $5$.


title: A new type of spectral mapping theorem for quantum walks with a moving shift on graphs

abstract: The conventional spectral mapping theorem for quantum walks can only be applied for walks employing a shift operator whose square is the identity. This theorem gives most of the eigenvalues of the time evolution $U$ by lifting the eigenvalues of an induced self-adjoint matrix $T$ onto the unit circle on the complex plane. We acquire a new spectral mapping theorem for the Grover walk with a shift operator whose cube is the identity on finite graphs. Moreover, graphs we can consider for a quantum walk with such a shift operator is characterized by a triangulation. We call these graphs triangulable graphs in this paper. One of the differences between our spectral mapping theorem and the conventional one is that lifting the eigenvalues of $T-1/2$ onto the unit circle gives most of the eigenvalues of $U$.


title: GAN Vocoder: Multi-Resolution Discriminator Is All You Need

abstract: Several of the latest GAN-based vocoders show remarkable achievements, outperforming autoregressive and flow-based competitors in both qualitative and quantitative measures while synthesizing orders of magnitude faster. In this work, we hypothesize that the common factor underlying their success is the multi-resolution discriminating framework, not the minute details in architecture, loss function, or training strategy. We experimentally test the hypothesis by evaluating six different generators paired with one shared multi-resolution discriminating framework. For all evaluative measures with respect to text-to-speech syntheses and for all perceptual metrics, their performances are not distinguishable from one another, which supports our hypothesis.


title: Column randomization and almost-isometric embeddings

abstract: The matrix $A:\mathbb{R}^n \to \mathbb{R}^m$ is $(\delta,k)$-regular if for any $k$-sparse vector $x$, $$ \left| \|Ax\|_2^2-\|x\|_2^2\right| \leq \delta \sqrt{k} \|x\|_2^2. $$ We show that if $A$ is $(\delta,k)$-regular for $1 \leq k \leq 1/\delta^2$, then by multiplying the columns of $A$ by independent random signs, the resulting random ensemble $A_\epsilon$ acts on an arbitrary subset $T \subset \mathbb{R}^n$ (almost) as if it were gaussian, and with the optimal probability estimate: if $\ell_*(T)$ is the gaussian mean-width of $T$ and $d_T=\sup_{t \in T} \|t\|_2$, then with probability at least $1-2\exp(-c(\ell_*(T)/d_T)^2)$, $$ \sup_{t \in T} \left| \|A_\epsilon t\|_2^2-\|t\|_2^2 \right| \leq C\left(\Lambda d_T \delta\ell_*(T)+(\delta \ell_*(T))^2 \right), $$ where $\Lambda=\max\{1,\delta^2\log(n\delta^2)\}$. This estimate is optimal for $0<\delta \leq 1/\sqrt{\log n}$.


title: Fast Statistical Leverage Score Approximation in Kernel Ridge Regression

abstract: Nystr\"om approximation is a fast randomized method that rapidly solves kernel ridge regression (KRR) problems through sub-sampling the n-by-n empirical kernel matrix appearing in the objective function. However, the performance of such a sub-sampling method heavily relies on correctly estimating the statistical leverage scores for forming the sampling distribution, which can be as costly as solving the original KRR. In this work, we propose a linear time (modulo poly-log terms) algorithm to accurately approximate the statistical leverage scores in the stationary-kernel-based KRR with theoretical guarantees. Particularly, by analyzing the first-order condition of the KRR objective, we derive an analytic formula, which depends on both the input distribution and the spectral density of stationary kernels, for capturing the non-uniformity of the statistical leverage scores. Numerical experiments demonstrate that with the same prediction accuracy our method is orders of magnitude more efficient than existing methods in selecting the representative sub-samples in the Nystr\"om approximation.


title: Transfer matrix study of the Anderson transition in non-Hermitian systems

abstract: In this paper, we present in-depth transfer matrix analyses of the Anderson transition in three non-Hermitian (NH) systems, NH Anderson, U(1) and Peierls models, that belong to NH class AI$^{\dagger}$ or NH class A. We first argue a general validity of the transfer matrix analysis, and clarify the symmetry properties of the Lyapunov exponents, scattering ($S$) matrix and two-terminal conductance in these NH models. The unitarity of the $S$ matrix is violated in NH systems, where the two-terminal conductance can take arbitrarily large values. Nonetheless, we show that the transposition symmetry of a Hamiltonian leads to the symmetric nature of the $S$ matrix as well as the reciprocal symmetries of the Lyapunov exponents and conductance in certain ways in these NH models. Finite size scaling data are fitted by polynomial functions, from which we determine the critical exponent $\nu$ at different single-particle energies and system parameters of the NH models. We conclude that the critical exponents of the NH class AI$^{\dagger}$ and the NH class A are $\nu=1.19 \pm 0.01$ and $\nu=1.00 \pm 0.04$, respectively. In the NH models, a distribution of the two-terminal conductance is not Gaussian. Instead, it contains small fractions of huge conductance values, which come from rare-event states with huge transmissions amplified by on-site NH disorders. Nonetheless, a geometric mean of the conductance enables the finite-size scaling analysis. We show that the critical exponents obtained from the conductance analysis are consistent with those from the localization length in these three NH models.


title: Electrostatic solitary waves in the Earth's bow shock: nature, properties, lifetimes and origin

abstract: We present a statistical analysis of more than two thousand bipolar electrostatic solitary waves (ESW) collected from ten quasi-perpendicular Earth's bow shock crossings by Magnetospheric Multiscale spacecraft. We developed and implemented a correction procedure for reconstruction of actual electric fields, velocities, and other properties of ESW from measurements, whose spatial scales are typically comparable with or smaller than spatial distance between voltage-sensitive probes. We determined the optimal ratio between frequency response factors of axial and spin plane antennas to be around 1.65/1.8. We found that more than 95\% of the ESW in the Earth's bow shock are of negative polarity and present an in depth analysis of properties of these ESW. They have spatial scales of about 10--100 m that is within a range of $\lambda_{D}$ to $10\lambda_{D}$, amplitudes typically below a few Volts that is below 0.1 of local electron temperature, and velocities below a few hundreds km/s in spacecraft and plasma rest frames that is on the order of local ion-acoustic speed. The spatial scales of ESW are distinctly correlated with local Debye length $\lambda_{D}$. ESW with amplitudes of 5--30 V or 0.1--0.3 Te have the occurrence rate of a few percent. The ESW have electric fields generally oblique to local magnetic field and propagate highly oblique to shock normal ${\bf N}$; more than 80\% of ESW propagate within 30$^{\circ}$ of the shock plane. In the shock plane, ESW typically propagate within a few tens of degrees of local magnetic field projection ${\bf B}_{\rm LM}$ onto the shock plane and preferentially opposite to ${\bf N}\times {\bf B}_{\rm LM}$. We argue that the ESW of negative polarity are ion phase space holes produced in a nonlinear stage of ion-ion ion-streaming instabilities. We estimated lifetimes of the ion holes to be 10--100 ms, or 1--10 km in terms of spatial distance.


title: Fast and Efficient Bit-Level Precision Tuning

abstract: In this article, we introduce a new technique for precision tuning. This problem consists of finding the least data types for numerical values such that the result of the computation satisfies some accuracy requirement. State of the art techniques for precision tuning use a try and fail approach. They change the data types of some variables of the program and evaluate the accuracy of the result. Depending on what is obtained, they change more or less data types and repeat the process. Our technique is radically different. Based on semantic equations, we generate an Integer Linear Problem (ILP) from the program source code. Basically, this is done by reasoning on the most significant bit and the number of significant bits of the values which are integer quantities. The integer solution to this problem, computed in polynomial time by a (real) linear programming solver, gives the optimal data types at the bit level. A finer set of semantic equations is also proposed which does not reduce directly to an ILP problem. So we use policy iteration to find the solution. Both techniques have been implemented and we show that our results encompass the results of state of the art tools.


title: A deep learning based known plaintext attack method for chaotic cryptosystem

abstract: In this paper, we propose a known-plaintext attack (KPA) method based on deep learning for traditional chaotic encryption scheme. We employ the convolutional neural network to learn the operation mechanism of chaotic cryptosystem, and accept the trained network as the final decryption system. To evaluate the attack performance of different networks on different chaotic cryptosystem, we adopt two neural networks to perform known-plaintext attacks on two distinct chaotic encryption schemes. The experimental results demonstrate the potential of deep learning-based method for known-plaintext attack against chaotic cryptosystem. Different from the previous known-plaintext attack methods, which were usually limited to a specific chaotic cryptosystem, a neural network can be applied to the cryptanalysis of various chaotic cryptosystems with deep learning-based approach, while several different networks can be designed for the cryptanalysis of chaotic cryptosystems. This paper provides a new idea for the cryptanalysis of chaotic image encryption algorithm.


title: On the Generalization Power of Overfitted Two-Layer Neural Tangent Kernel Models

abstract: In this paper, we study the generalization performance of min $\ell_2$-norm overfitting solutions for the neural tangent kernel (NTK) model of a two-layer neural network with ReLU activation that has no bias term. We show that, depending on the ground-truth function, the test error of overfitted NTK models exhibits characteristics that are different from the "double-descent" of other overparameterized linear models with simple Fourier or Gaussian features. Specifically, for a class of learnable functions, we provide a new upper bound of the generalization error that approaches a small limiting value, even when the number of neurons $p$ approaches infinity. This limiting value further decreases with the number of training samples $n$. For functions outside of this class, we provide a lower bound on the generalization error that does not diminish to zero even when $n$ and $p$ are both large.


title: ModelingToolkit: A Composable Graph Transformation System For Equation-Based Modeling

abstract: Getting good performance out of numerical equation solvers requires that the user has provided stable and efficient functions representing their model. However, users should not be trusted to write good code. In this manuscript we describe ModelingToolkit (MTK), a symbolic equation-based modeling system which allows for composable transformations to generate stable, efficient, and parallelized model implementations. MTK blurs the lines of traditional symbolic computing by acting directly on a user's numerical code. We show the ability to apply graph algorithms for automatically parallelizing and performing index reduction on code written for differential-algebraic equation (DAE) solvers, "fixing" the performance and stability of the model without requiring any changes to on the user's part. We demonstrate how composable model transformations can be combined with automated data-driven surrogate generation techniques, allowing machine learning methods to generate accelerated approximate models within an acausal modeling framework. These reduced models are shown to outperform the Dymola Modelica compiler on an HVAC model by 590x at 3\% error. Together, this demonstrates MTK as a system for bringing the latest research in graph transformations directly to modeling applications.