title: Any Target Function Exists in a Neighborhood of Any Sufficiently Wide Random Network: A Geometrical Perspective

abstract: It is known that any target function is realized in a sufficiently small neighborhood of any randomly connected deep network, provided the width (the number of neurons in a layer) is sufficiently large. There are sophisticated theories and discussions concerning this striking fact, but rigorous theories are very complicated. We give an elementary geometrical proof by using a simple model for the purpose of elucidating its structure. We show that high-dimensional geometry plays a magical role: When we project a high-dimensional sphere of radius 1 to a low-dimensional subspace, the uniform distribution over the sphere reduces to a Gaussian distribution of negligibly small covariances.


title: Covert Communication in Continuous-Time Systems

abstract: Recent works have considered the ability of transmitter Alice to communicate reliably to receiver Bob without being detected by warden Willie. These works generally assume a standard discrete-time model. But the assumption of a discrete-time model in standard communication scenarios is often predicated on its equivalence to a continuous-time model, which has not been established for the covert communications problem. Here, we consider the continuous-time channel directly and study if efficient covert communication can still be achieved. We assume that an uninformed jammer is present to assist Alice, and we consider additive white Gaussian noise (AWGN) channels between all parties. For a channel with approximate bandwidth W, we establish constructions such that O(WT) information bits can be transmitted covertly and reliably from Alice to Bob in T seconds for two separate scenarios: 1) when the path-loss between Alice and Willie is known; and 2) when the path-loss between Alice and Willie is unknown.


title: Fides: Managing Data on Untrusted Infrastructure

abstract: Significant amounts of data are currently being stored and managed on third-party servers. It is impractical for many small scale enterprises to own their private datacenters, hence renting third-party servers is a viable solution for such businesses. But the increasing number of malicious attacks, both internal and external, as well as buggy software on third-party servers is causing clients to lose their trust in these external infrastructures. While small enterprises cannot avoid using external infrastructures, they need the right set of protocols to manage their data on untrusted infrastructures. In this paper, we propose TFCommit, a novel atomic commitment protocol that executes transactions on data stored across multiple untrusted servers. To our knowledge, TFCommit is the first atomic commitment protocol to execute transactions in an untrusted environment without using expensive Byzantine replication. Using TFCommit, we propose an auditable data management system, Fides, residing completely on untrustworthy infrastructure. As an auditable system, Fides guarantees the detection of potentially malicious failures occurring on untrusted servers using tamper-resistant logs with the support of cryptographic techniques. The experimental evaluation demonstrates the scalability and the relatively low overhead of our approach that allows executing transactions on untrusted infrastructure.


title: Spectral conditions for graph rigidity in the Euclidean plane

abstract: Rigidity is the property of a structure that does not flex. It is well studied in discrete geometry and mechanics, and has applications in material science, engineering and biological sciences. A bar-and-joint framework is a pair $(G,p)$ of graph $G$ together with a map $p$ of the vertices of $G$ into the Euclidean plane. We view the edges of $(G, p)$ as bars and the vertices as universal joints. The vertices can move continuously as long as the distances between pairs of adjacent vertices are preserved. The framework is rigid if any such motion preserves the distances between all pairs of vertices. In 1970, Laman obtained a combinatorial characterization of rigid graphs in the Euclidean plane. In 1982, Lov\'asz and Yemini discovered a new characterization and proved that every $6$-connected graph is rigid. Combined with a characterization of global rigidity, their proof actually implies that every 6-connected graph is globally rigid. Consequently, if Fiedler's algebraic connectivity is greater than 5, then $G$ is globally rigid. In this paper, we improve this bound and show that for a graph $G$ with minimum degree $\delta\geq 6$, if its algebraic connectivity is greater than $2+\frac{1}{\delta -1}$, then $G$ is rigid and if its algebraic connectivity is greater than $2+\frac{2}{\delta -1}$, then $G$ is globally rigid. Our results imply that every connected regular Ramanujan graph with degree at least $8$ is globally rigid. We also prove a more general result giving a sufficient spectral condition for the existence of $k$ edge-disjoint spanning rigid subgraphs. The same condition implies that a graph contains $k$ edge-disjoint spanning $2$-connected subgraphs. This result extends previous spectral conditions for packing edge-disjoint spanning trees.


title: 75,000,000,000 Streaming Inserts/Second Using Hierarchical Hypersparse GraphBLAS Matrices

abstract: The SuiteSparse GraphBLAS C-library implements high performance hypersparse matrices with bindings to a variety of languages (Python, Julia, and Matlab/Octave). GraphBLAS provides a lightweight in-memory database implementation of hypersparse matrices that are ideal for analyzing many types of network data, while providing rigorous mathematical guarantees, such as linearity. Streaming updates of hypersparse matrices put enormous pressure on the memory hierarchy. This work benchmarks an implementation of hierarchical hypersparse matrices that reduces memory pressure and dramatically increases the update rate into a hypersparse matrices. The parameters of hierarchical hypersparse matrices rely on controlling the number of entries in each level in the hierarchy before an update is cascaded. The parameters are easily tunable to achieve optimal performance for a variety of applications. Hierarchical hypersparse matrices achieve over 1,000,000 updates per second in a single instance. Scaling to 31,000 instances of hierarchical hypersparse matrices arrays on 1,100 server nodes on the MIT SuperCloud achieved a sustained update rate of 75,000,000,000 updates per second. This capability allows the MIT SuperCloud to analyze extremely large streaming network data sets.


title: Lp-improving properties of certain singular measures on the Heisenberg group

abstract: We study Lp-improving properties as well as the type set of certain singular measures on the Heisenberg group.


title: A Review on Generative Adversarial Networks: Algorithms, Theory, and Applications

abstract: Generative adversarial networks (GANs) are a hot research topic recently. GANs have been widely studied since 2014, and a large number of algorithms have been proposed. However, there is few comprehensive study explaining the connections among different GANs variants, and how they have evolved. In this paper, we attempt to provide a review on various GANs methods from the perspectives of algorithms, theory, and applications. Firstly, the motivations, mathematical representations, and structure of most GANs algorithms are introduced in details. Furthermore, GANs have been combined with other machine learning algorithms for specific applications, such as semi-supervised learning, transfer learning, and reinforcement learning. This paper compares the commonalities and differences of these GANs methods. Secondly, theoretical issues related to GANs are investigated. Thirdly, typical applications of GANs in image processing and computer vision, natural language processing, music, speech and audio, medical field, and data science are illustrated. Finally, the future open research problems for GANs are pointed out.


title: Memory capacity of neural networks with threshold and ReLU activations

abstract: Overwhelming theoretical and empirical evidence shows that mildly overparametrized neural networks -- those with more connections than the size of the training data -- are often able to memorize the training data with $100\%$ accuracy. This was rigorously proved for networks with sigmoid activation functions and, very recently, for ReLU activations. Addressing a 1988 open question of Baum, we prove that this phenomenon holds for general multilayered perceptrons, i.e. neural networks with threshold activation functions, or with any mix of threshold and ReLU activations. Our construction is probabilistic and exploits sparsity.


title: Time-domain digital-to-analog converter for spiking neural network hardware

abstract: We propose a new digital-to-analog converter (DAC) for realizing a synapse circuit of mixed-signal spiking neural networks. We named this circuit "time-domain DAC (TDAC)". This produces weights for converting a digital input code into voltage using one current waveform. Therefore, a TDAC is more compact than a conventional DAC that comprises many current sources and resistors. Moreover, a TDAC with leak resistance reproduces biological plausible synaptic responses that are expressed as alpha functions or dual exponential equations. We will show numerical analysis results of a TDAC and circuit simulation results of a circuit designed by the TSMC 40 nm CMOS process.


title: Reinforcement Learning with Probabilistically Complete Exploration

abstract: Balancing exploration and exploitation remains a key challenge in reinforcement learning (RL). State-of-the-art RL algorithms suffer from high sample complexity, particularly in the sparse reward case, where they can do no better than to explore in all directions until the first positive rewards are found. To mitigate this, we propose Rapidly Randomly-exploring Reinforcement Learning (R3L). We formulate exploration as a search problem and leverage widely-used planning algorithms such as Rapidly-exploring Random Tree (RRT) to find initial solutions. These solutions are used as demonstrations to initialize a policy, then refined by a generic RL algorithm, leading to faster and more stable convergence. We provide theoretical guarantees of R3L exploration finding successful solutions, as well as bounds for its sampling complexity. We experimentally demonstrate the method outperforms classic and intrinsic exploration techniques, requiring only a fraction of exploration samples and achieving better asymptotic performance.


title: Gate modulation of the long-range magnetic order in a vanadium-doped WSe2 semiconductor

abstract: We demonstrate the gate-tunability of the long-range magnetic order in a p-type V-doped WSe2 monolayer using ab initio calculations. We found that at a low V-doping concentration limit, the long-range ferromagnetic order is enhanced by increasing the hole density. In contrast, the short-range antiferromagnetic order is manifested at a high electron density by full compensation of the p-type V doping concentration. The hole-mediated long-range magnetic exchange is ~70 meV, thus strongly suggesting the ferromagnetism in V-doped WSe2 at room temperature. Our findings on strong coupling between charge and spin order in V-doped WSe2 provide plenty of room for multifunctional gate-tunable spintronics.


title: Quantitative Aspects of Programming Languages and Systems over the past $2^4$ years and beyond

abstract: Quantitative aspects of computation are related to the use of both physical and mathematical quantities, including time, performance metrics, probability, and measures for reliability and security. They are essential in characterizing the behaviour of many critical systems and in estimating their properties. Hence, they need to be integrated both at the level of system modeling and within the verification methodologies and tools. Along the last two decades a variety of theoretical achievements and automated techniques have contributed to make quantitative modeling and verification mainstream in the research community. In the same period, they represented the central theme of the series of workshops entitled Quantitative Aspects of Programming Languages and Systems (QAPL) and born in 2001. The aim of this survey is to revisit such achievements and results from the standpoint of QAPL and its community.


title: Probabilistic Output Analyses for Deterministic Programs --- Reusing Existing Non-probabilistic Analyses

abstract: We consider reusing established non-probabilistic output analyses (either forward or backwards) that yield over-approximations of a program's pre-image or image relation, e.g., interval analyses. We assume a probability measure over the program input and present two techniques (one for forward and one for backward analyses) that both derive upper and lower probability bounds for the output events. We demonstrate the most involved technique, namely the forward technique, for two examples and compare their results to a cutting-edge probabilistic output analysis.


title: Nested-Wasserstein Self-Imitation Learning for Sequence Generation

abstract: Reinforcement learning (RL) has been widely studied for improving sequence-generation models. However, the conventional rewards used for RL training typically cannot capture sufficient semantic information and therefore render model bias. Further, the sparse and delayed rewards make RL exploration inefficient. To alleviate these issues, we propose the concept of nested-Wasserstein distance for distributional semantic matching. To further exploit it, a novel nested-Wasserstein self-imitation learning framework is developed, encouraging the model to exploit historical high-rewarded sequences for enhanced exploration and better semantic matching. Our solution can be understood as approximately executing proximal policy optimization with Wasserstein trust-regions. Experiments on a variety of unconditional and conditional sequence-generation tasks demonstrate the proposed approach consistently leads to improved performance.


title: Averaging Principles for Mixed Fast-Slow Systems Driven by Fractional Brownian Motion

abstract: We focus on fast-slow systems involving both fractional Brownian motion (fBm) and standard Brownian motion (Bm). The integral with respect to Bm is the standard Ito integral, and the integral with respect to fBm is the generalised Riemann-Stieltjes integral using the tools of fractional calculus. An averaging principle in which the fast-varying diffusion process of the fast-slow systems acts as a noise to be averaged out in the limit is established. It is shown that the slow process has a limit in the mean square sense, which is characterized by the solution of stochastic differential equations driven by fBm whose coefficients are averaged with respect to the stationary measure of the fast-varying diffusion. The implication is that one can ignore the complex original systems and concentrate on the averaged systems instead. This averaging principle paves the way for reduction of computational complexity.


title: Chaotic D1-D5 Black Hole Dynamics through Networks

abstract: This work studies dynamics controlling the transition between different microstates of two charge D1-D5 black holes by network methods, in which microstates of the system are defined as network nodes, while transitions between them are defined as edges. It is found that the eigenspectrum of this network's Laplacian matrix, which is identified with Hamiltonians of the microstate system, has completely the same Nearest-Neighbor Spacing Distribution as that of general Gaussian Orthogonal Ensemble of Random Matrices. According to the BGS, i.e. Bohigas, Giannoni and Schmit conjecture, this forms evidence for chaotic features of the D1-D5 microstate dynamics. This evidence is further strengthened by observations that inverse of the first/minimal nonzero eigenvalue of the Laplacian matrix is proportional to logarithms of the microstate number of the system. By Sekino and Susskind, this means that dynamics of the D1-D5 black hole microstates are not only chaotic, but also the fastest scrambler in nature.


title: The Herglotz wave function, the Vekua transform and the enclosure method

abstract: This paper gives applications of the enclosure method introduced by the author to typical inverse obstacle and crack scattering problems in two dimensions. Explicit extraction formulae of the convex hull of unknown polygonal sound-hard obstacles and piecewise linear cracks from the far field pattern of the scattered field at a fixed wave number and at most two incident directions are given. The main new points of this paper are: a combination of the enclosure method and the Herglotz wave function; explicit construction of the density in the Herglotz wave function by using the idea of the Vekua transform. By virtue of the construction, one can avoid any restriction on the wave number in the extraction formulae. An attempt for the case when the far field pattern is given on limited angles is also given.


title: Antibiotics, duration of infectiousness, and transmission of disease

abstract: Humans, domestic animals, orchard crops, and ornamental plants are commonly treated with antibiotics in response to bacterial infection. By curing infectious individuals, antibiotic therapy might limit the spread of contagious disease among hosts. But antibiotic suppression of within-host pathogen density might also reduce the probability that the host is otherwise removed from infectious status before recovery. When rates of both recovery and removal (isolation or mortality) depend directly on within-host density, antibiotic therapy can relax the removal rate and so increase between-host disease transmission. In this paper a deterministic within-host dynamics drives the infectious host's probability of infection transmission, as well as the host's time-dependent probability of surviving to recovery. The model varies (1) inoculum size, (2) the time elapsing between infection and initiation of therapy, (2) antibiotic efficacy, and (3) the size/susceptibility of groups encountered by an infectious host. Results identify conditions where antibiotic treatment simultaneously increases host survival and increases the expected number of new infections. That is, antibiotics might convert a rare, serious bacterial disease into a common, but treatable infection.


title: Prospect of undoped inorganic crystals at 77 Kelvin for low-mass dark matter search at Spallation Neutron Source

abstract: Investigated in this work were sensitivities of a prototype detector for the detection of low-mass dark matter particles produced at the Spallation Neutron Source at the Oak Ridge National Laboratory in two years of data taking. The presumed prototype consisted of 10 kg undoped CsI or NaI scintillation crystals directly coupled with SiPM arrays operated at 77 K. Compared to the COHERENT CsI(Na) detector, a much higher light yield was assumed for the prototype. An experiment with a cylindrical 1~kg undoped CsI crystal coupled directly to two photomultiplier tubes at about 77~K was conducted as the first step to verify the idea. A light yield of $26.0 \pm 0.4$ photoelectrons per keV electron-equivalent was achieved. This eliminated the concern of self light absorption in large crystals raised in some of the early studies.


title: Homogenization of stationary Navier--Stokes--Fourier system in domains with tiny holes

abstract: We study the homogenization of stationary compressible Navier--Stokes--Fourier system in a bounded three dimensional domain perforated with a large number of very tiny holes. Under suitable assumptions imposed on the smallness and distribution of the holes, we show that the homogenized limit system remains the same in the domain without holes.


title: Sparse Nonnegative Matrix Factorization for Multiple Local Community Detection

abstract: Local community detection consists of finding a group of nodes closely related to the seeds, a small set of nodes of interest. Such group of nodes are densely connected or have a high probability of being connected internally than their connections to other clusters in the network. Existing local community detection methods focus on finding either one local community that all seeds are most likely to be in or finding a single community for each of the seeds. However, a seed member usually belongs to multiple local overlapping communities. In this work, we present a novel method of detecting multiple local communities to which a single seed member belongs. The proposed method consists of three key steps: (1) local sampling with Personalized PageRank (PPR); (2) using the sparseness generated by a sparse nonnegative matrix factorization (SNMF) to estimate the number of communities in the sampled subgraph; (3) using SNMF soft community membership vectors to assign nodes to communities. The proposed method shows favorable accuracy performance and a good conductance when compared to state-of-the-art community detection methods by experiments using a combination of artificial and real-world networks.


title: Streaming Transformations of Infinite Ordered-Data Words

abstract: In this paper, we define streaming register transducer (SRT), a one-way, letter-to-letter, transductional machine model for transformations of infinite data words whose data domain forms a linear group. Comparing with existing data word transducers, SRT are able to perform two extra operations on the registers: a linear-order-based comparison and an additive update. We consider the transformations that can be defined by SRT and several subclasses of SRT. We investigate the expressiveness of these languages and several decision problems. Our main results include: 1) SRT are closed under union and intersection, and add-free SRT are also closed under composition; 2) SRT-definable transformations can be defined in monadic second-order (MSO) logic, but are not comparable with first-order (FO) definable transformations; 3) the functionality problem is decidable for add-free SRT, the reactivity problem and inclusion problem are decidable for deterministic add-free SRT, but none of these problems is decidable in general for SRT.


title: On Exactly $3$-Deficient-Perfect Numbers

abstract: Let $n$ and $k$ be positive integers and $\sigma(n)$ the sum of all positive divisors of $n$. We call $n$ an exactly $k$-deficient-perfect number with deficient divisors $d_1, d_2, \ldots, d_k$ if $d_1, d_2, \ldots, d_k$ are distinct proper divisors of $n$ and $\sigma (n)=2n-(d_1+d_2+\ldots + d_k)$. In this article, we show that the only odd exactly $3$-deficient-perfect number with at most two distinct prime factors is $1521=3^2 \cdot 13^2$.


title: CNN-based InSAR Denoising and Coherence Metric

abstract: Interferometric Synthetic Aperture Radar (InSAR) imagery for estimating ground movement, based on microwaves reflected off ground targets is gaining increasing importance in remote sensing. However, noise corrupts microwave reflections received at satellite and contaminates the signal's wrapped phase. We introduce Convolutional Neural Networks (CNNs) to this problem domain and show the effectiveness of autoencoder CNN architectures to learn InSAR image denoising filters in the absence of clean ground truth images, and for artefact reduction in estimated coherence through intelligent preprocessing of training data. We compare our results with four established methods to illustrate superiority of proposed method.


title: Asymptotic normality of a linear threshold estimator in fixed dimension with near-optimal rate

abstract: Linear thresholding models postulate that the conditional distribution of a response variable in terms of covariates differs on the two sides of a (typically unknown) hyperplane in the covariate space. A key goal in such models is to learn about this separating hyperplane. Exact likelihood or least squares methods to estimate the thresholding parameter involve an indicator function which make them difficult to optimize and are, therefore, often tackled by using a surrogate loss that uses a smooth approximation to the indicator. In this paper, we demonstrate that the resulting estimator is asymptotically normal with a near optimal rate of convergence: $n^{-1}$ up to a log factor, in both classification and regression thresholding models. This is substantially faster than the currently established convergence rates of smoothed estimators for similar models in the statistics and econometrics literatures. We also present a real-data application of our approach to an environmental data set where $CO_2$ emission is explained in terms of a separating hyperplane defined through per-capita GDP and urban agglomeration.


title: CNN-based InSAR Coherence Classification

abstract: Interferometric Synthetic Aperture Radar (InSAR) imagery based on microwaves reflected off ground targets is becoming increasingly important in remote sensing for ground movement estimation. However, the reflections are contaminated by noise, which distorts the signal's wrapped phase. Demarcation of image regions based on degree of contamination ("coherence") is an important component of the InSAR processing pipeline. We introduce Convolutional Neural Networks (CNNs) to this problem domain and show their effectiveness in improving coherence-based demarcation and reducing misclassifications in completely incoherent regions through intelligent preprocessing of training data. Quantitative and qualitative comparisons prove superiority of proposed method over three established methods.


title: Test of the unitary coupled-cluster variational quantum eigensolver for a simple strongly correlated condensed-matter system

abstract: The variational quantum eigensolver has been proposed as a low-depth quantum circuit that can be employed to examine strongly correlated systems on today's noisy intermediate-scale quantum computers. We examine details associated with the factorized form of the unitary coupled-cluster variant of this algorithm. We apply it to a simple strongly correlated condensed-matter system with nontrivial behavior---the four-site Hubbard model at half filling. This work show some of the subtle issues one needs to take into account when applying this algorithm in practice, especially to condensed-matter systems.


title: Globally optimal dense and sparse spanning trees, and their applications

abstract: Finding spanning trees under various constraints is a classic problem with applications in many fields. Recently, a novel notion of "dense" ("sparse") tree, and in particular spanning tree (DST and SST respectively), is introduced as the structure that have a large (small) number of subtrees, or small (large) sum of distances between vertices. We show that finding DST and SST reduces to solving the discrete optimization problems. New and efficient approaches to find such spanning trees is achieved by imposing certain conditions on the vertex degrees which are then used to define an objective function that is minimized over all spanning trees of the graph under consideration. Solving this minimization problem exactly may be prohibitively time consuming for large graphs. Hence, we propose to use genetic algorithm (GA) which is one of well known metaheuristics methods to solve DST and SST approximately. As far as we are aware this is the first time GA has been used in this context. We also demonstrate on a number of applications that GA approach is well suited for these types of problems both in computational efficiency and accuracy of the approximate solution. Furthermore, we improve the efficiency of the proposed method by using Kruskal's algorithm in combination with GA. The application of our methods to several practical large graphs and networks is presented. Computational results show that they perform faster than previously proposed heuristic methods and produce more accurate solutions. Furthermore, the new feature of the proposed approach is that it can be applied recursively to sub-trees or spanning trees with additional constraints in order to further investigate the graphical properties of the graph and/or network. The application of this methodology on the gene network of a cancer cell led to isolating key genes in a network that were not obvious from previous studies.


title: Non-Orthogonal Multiple Access with Wireless Caching for 5G-Enabled Vehicular Networks

abstract: The proliferation of connected vehicles along with the high demand for rich multimedia services constitute key challenges for the emerging 5G-enabled vehicular networks. These challenges include, but are not limited to, high spectral efficiency and low latency requirements. Recently, the integration of cache-enabled networks with non-orthogonal multiple access (NOMA) has been shown to reduce the content delivery time and traffic congestion in wireless networks. Ac-cordingly, in this article, we envisage cache-aided NOMA as a technology facilitator for 5G-enabled vehicular networks. In particular, we present a cache-aided NOMA architecture, which can address some of the aforementioned challenges in these networks. We demonstrate that the spectral efficiency gain of the proposed architecture, which depends largely on the cached contents, significantly outperforms that of conventional vehicular networks. Finally, we provide deep insights into the challenges, opportunities, and future research trends that will enable the practical realization of cache-aided NOMA in 5G-enabled vehicular networks.


title: CSR wake expressions with two bending magnets and simulation results for the multi-turn ERL, CBETA

abstract: This paper consists of two main parts regarding Coherent Synchrotron Radiation (CSR). The first part extends the CSR theory of two particle interaction from a system of one bending magnet to two bending magnets, where the wake can leak from the first to the second. The new theory agrees well with the established simulation code Bmad. The second part of the paper presents the CSR simulation results on CBETA, the Cornell BNL Energy-Recovery-Linac (ERL) Test Accelerator, in which the magnets are so close to each other and the new extended theory becomes important. CBETA is the first multi-turn ERL with Superconducting Radio Frequency (SRF) accelerating cavities and a Fixed Field Alternating gradient (FFA) beamline. Simulations show that CSR causes phase space dilution that becomes more significant as the bunch charge and the number of recirculation passes increase. Potential ways to mitigate the CSR effects, including adding vacuum chamber shielding and increasing bunch length, are being investigated.


title: CNN-Based Real-Time Parameter Tuning for Optimizing Denoising Filter Performance

abstract: We propose a novel direction to improve the denoising quality of filtering-based denoising algorithms in real time by predicting the best filter parameter value using a Convolutional Neural Network (CNN). We take the use case of BM3D, the state-of-the-art filtering-based denoising algorithm, to demonstrate and validate our approach. We propose and train a simple, shallow CNN to predict in real time, the optimum filter parameter value, given the input noisy image. Each training example consists of a noisy input image (training data) and the filter parameter value that produces the best output (training label). Both qualitative and quantitative results using the widely used PSNR and SSIM metrics on the popular BSD68 dataset show that the CNN-guided BM3D outperforms the original, unguided BM3D across different noise levels. Thus, our proposed method is a CNN-based improvement on the original BM3D which uses a fixed, default parameter value for all images.


title: On the Joint Typicality of Permutations of Sequences of Random Variables

abstract: Permutations of correlated sequences of random variables appear naturally in a variety of applications such as graph matching and asynchronous communications. In this paper, the asymptotic statistical behavior of such permuted sequences is studied. It is assumed that a collection of random vectors is produced based on an arbitrary joint distribution, and the vectors undergo a permutation operation. The joint typicality of the resulting permuted vectors with respect to the original distribution is investigated. As an initial step, permutations of pairs of correlated random vectors are considered. It is shown that the probability of joint typicality of the permuted vectors depends only on the number and length of the disjoint cycles of the permutation. Consequently, it suffices to study typicality for a class of permutations called 'standard permutations', for which, upper-bounds on the probability of joint typicality are derived. The notion of standard permutations is extended to a class of permutation vectors called 'Bell permutation vectors'. By investigating Bell permutation vectors, upper-bounds on the probability of joint typicality of permutations of arbitrary collections of random sequences are derived.


title: A Novel Image Dehazing and Assessment Method

abstract: Images captured in hazy weather conditions often suffer from color contrast and color fidelity. This degradation is represented by transmission map which represents the amount of attenuation and airlight which represents the color of additive noise. In this paper, we have proposed a method to estimate the transmission map using haze levels instead of airlight color since there are some ambiguities in estimation of airlight. Qualitative and quantitative results of proposed method show competitiveness of the method given. In addition we have proposed two metrics which are based on statistics of natural outdoor images for assessment of haze removal algorithms.


title: Parity violation in high-energy proton-proton scattering

abstract: A new approach is proposed to describe the parity violation in the process of scattering of a polarized proton by a proton at high energies. It is shown that in the tree approximation the P-odd amplitude of proton-proton scattering is substantially smaller than the P-odd amplitude in proton-neutron scattering. The main source of P-odd asymmetry in proton-proton scattering is radiative corrections due to the charge-exchange strong interaction. We predict that the asymmetry in the inelastic cross section is substantially less than the asymmetry in the elastic cross section.


title: G2MF-WA: Geometric Multi-Model Fitting with Weakly Annotated Data

abstract: In this paper we attempt to address the problem of geometric multi-model fitting with resorting to a few weakly annotated (WA) data points, which has been sparsely studied so far. In weak annotating, most of the manual annotations are supposed to be correct yet inevitably mixed with incorrect ones. The WA data can be naturally obtained in an interactive way for specific tasks, for example, in the case of homography estimation, one can easily annotate points on the same plane/object with a single label by observing the image. Motivated by this, we propose a novel method to make full use of the WA data to boost the multi-model fitting performance. Specifically, a graph for model proposal sampling is first constructed using the WA data, given the prior that the WA data annotated with the same weak label has a high probability of being assigned to the same model. By incorporating this prior knowledge into the calculation of edge probabilities, vertices (i.e., data points) lie on/near the latent model are likely to connect together and further form a subset/cluster for effective proposals generation. With the proposals generated, the $\alpha$-expansion is adopted for labeling, and our method in return updates the proposals. This works in an iterative way. Extensive experiments validate our method and show that the proposed method produces noticeably better results than state-of-the-art techniques in most cases.


title: Reconciling the Gaussian and Whittle Likelihood with an application to estimation in the frequency domain

abstract: In time series analysis there is an apparent dichotomy between time and frequency domain methods. The aim of this paper is to draw connections between frequency and time domain methods. Our focus will be on reconciling the Gaussian likelihood and the Whittle likelihood. We derive an exact, interpretable, bound between the Gaussian and Whittle likelihood of a second order stationary time series. The derivation is based on obtaining the transformation which is biorthogonal to the discrete Fourier transform of the time series. Such a transformation yields a new decomposition for the inverse of a Toeplitz matrix and enables the representation of the Gaussian likelihood within the frequency domain. We show that the difference between the Gaussian and Whittle likelihood is due to the omission of the best linear predictions outside the domain of observation in the periodogram associated with the Whittle likelihood. Based on this result, we obtain an approximation for the difference between the Gaussian and Whittle likelihoods in terms of the best fitting, finite order autoregressive parameters. These approximations are used to define two new frequency domain quasi-likelihoods criteria. We show that these new criteria can yield a better approximation of the spectral divergence criterion, as compared to both the Gaussian and Whittle likelihoods. In simulations, we show that the proposed estimators have satisfactory finite sample properties.


title: A hybrid algorithm for disparity calculation from sparse disparity estimates based on stereo vision

abstract: In this paper, we have proposed a novel method for stereo disparity estimation by combining the existing methods of block based and region based stereo matching. Our method can generate dense disparity maps from disparity measurements of only 18% pixels of either the left or the right image of a stereo image pair. It works by segmenting the lightness values of image pixels using a fast implementation of K-Means clustering. It then refines those segment boundaries by morphological filtering and connected components analysis, thus removing a lot of redundant boundary pixels. This is followed by determining the boundaries' disparities by the SAD cost function. Lastly, we reconstruct the entire disparity map of the scene from the boundaries' disparities through disparity propagation along the scan lines and disparity prediction of regions of uncertainty by considering disparities of the neighboring regions. Experimental results on the Middlebury stereo vision dataset demonstrate that the proposed method outperforms traditional disparity determination methods like SAD and NCC by up to 30% and achieves an improvement of 2.6% when compared to a recent approach based on absolute difference (AD) cost function for disparity calculations [1].


title: FD-GAN: Generative Adversarial Networks with Fusion-discriminator for Single Image Dehazing

abstract: Recently, convolutional neural networks (CNNs) have achieved great improvements in single image dehazing and attained much attention in research. Most existing learning-based dehazing methods are not fully end-to-end, which still follow the traditional dehazing procedure: first estimate the medium transmission and the atmospheric light, then recover the haze-free image based on the atmospheric scattering model. However, in practice, due to lack of priors and constraints, it is hard to precisely estimate these intermediate parameters. Inaccurate estimation further degrades the performance of dehazing, resulting in artifacts, color distortion and insufficient haze removal. To address this, we propose a fully end-to-end Generative Adversarial Networks with Fusion-discriminator (FD-GAN) for image dehazing. With the proposed Fusion-discriminator which takes frequency information as additional priors, our model can generator more natural and realistic dehazed images with less color distortion and fewer artifacts. Moreover, we synthesize a large-scale training dataset including various indoor and outdoor hazy images to boost the performance and we reveal that for learning-based dehazing methods, the performance is strictly influenced by the training data. Experiments have shown that our method reaches state-of-the-art performance on both public synthetic datasets and real-world images with more visually pleasing dehazed results.


title: Effect of feedback of massive stars in the fragmentation, distribution, and kinematics of the gas in two star forming regions in the Carina Nebula

abstract: We present ALMA high spatial resolution observations towards two star forming regions located in one of the most extreme zones of star formation in the Galaxy, the Carina Nebula. One region is located at the center of the nebula and is severally affected by the stellar feedback from high-mass stars, while the other region is located further south and is less disturbed by the massive star clusters. We found that the region at the center of the nebula is forming less but more massive cores than the region located in the south, suggesting that the level of stellar feedback effectively influence the fragmentation process in clumps. Lines such as HCN, HCO$^{+}$ and SiO show abundant and complex gas distributions in both regions, confirming the presence of ionization and shock fronts. Jeans analysis suggests that the observed core masses in the region less affected by the massive stars are consistent with thermal fragmentation, but turbulent Jeans fragmentation might explain the high masses of the cores identified in the region in the center of Carina. Consistently, two different analyses in the HCO$^{+}$ line provided evidence for a higher level of turbulence in the gas more affected by the stellar feedback. The gas column density probability functions, N-PDFs, show log-normal shapes with clear transitions to power law regimes. We observed a wider N-PDF in the region at the center of the nebula, which provides further evidence for a higher level of turbulence in the material with a higher level of massive stellar feedback.


title: Finding the Sparsest Vectors in a Subspace: Theory, Algorithms, and Applications

abstract: The problem of finding the sparsest vector (direction) in a low dimensional subspace can be considered as a homogeneous variant of the sparse recovery problem, which finds applications in robust subspace recovery, dictionary learning, sparse blind deconvolution, and many other problems in signal processing and machine learning. However, in contrast to the classical sparse recovery problem, the most natural formulation for finding the sparsest vector in a subspace is usually nonconvex. In this paper, we overview recent advances on global nonconvex optimization theory for solving this problem, ranging from geometric analysis of its optimization landscapes, to efficient optimization algorithms for solving the associated nonconvex optimization problem, to applications in machine intelligence, representation learning, and imaging sciences. Finally, we conclude this review by pointing out several interesting open problems for future research.

categories: cs.LG cs.IT eess.IV math.IT math.OC stat.ML


title: Multi-switches and virtual knot invariants

abstract: In the paper we introduce a general approach how for a given virtual biquandle multi-switch $(S,V)$ on an algebraic system $X$ (from some category) and a given virtual link $L$ construct an algebraic system $X_{S,V}(L)$ (from the same category) which is an invariant of $L$. As a corollary we introduce a new quandle invariant for virtual links which generalizes previously known quandle invariants for virtual links.


title: Counter-propagating edge states in Floquet topological insulating phases

abstract: Nonequilibrium Floquet topological phases due to periodic driving are known to exhibit rich and interesting features with no static analogs. Various known topological invariants usually proposed to characterize static topological systems often fail to fully characterize Floquet topological phases. This fact has motivated extensive studies of Floquet topological phases to better understand nonequilibrium topological phases and to explore their possible applications. Here we present a theoretically simple Floquet topological insulating system that may possess an arbitrary number of counter-propagating chiral edge states. Further investigation into our system reveals another related feature by tuning the same set of system parameters, namely, the emergence of almost flat (dispersionless) edge modes. In particular, we employ two-terminal conductance and dynamical winding numbers to characterize counter-propagating chiral edge states. We further demonstrate the robustness of such edge states against symmetry preserving disorder. Finally, we identify an emergent chiral symmetry at certain sub-regimes of the Brillouin zone that can explain the presence of almost flat edge modes. Our results have exposed more interesting possibilities in Floquet topological matter.


title: Kink Dynamics in a Nonlinear Beam Model

abstract: In this paper, we study the single kink and the kink-antikink collisions of a nonlinear beam equation bearing a fourth-derivative term. We numerically explore some of the key characteristics of the single kink both in its standing wave and in its traveling wave form. A point of emphasis is the study of kink-antikink collisions, exploring the critical velocity for single-bounce (and separation) and infinite-bounce (where the kink and antikink trap each other) windows. The relevant phenomenology turns out to be dramatically different than that of the corresponding nonlinear Klein-Gordon (i.e., $\phi^4$) model. Our computations show that for small initial velocities, the kink and antikink reflect nearly elastically without colliding. For an intermediate interval of velocities, the two waves trap each other, while for large speeds a single inelastic collision between them takes place. Lastly, we briefly touch upon the use of collective coordinates (CC) method and their predictions of the relevant phenomenology. When one degree of freedom is used in the CC approach, the results match well the numerical ones for small values of initial velocity. However, for bigger values of initial velocity, it is inferred that more degrees of freedom need to be self-consistently included in order to capture the collision phenomenology.


title: Investigation of Patient-sharing Networks Using a Bayesian Network Model Selection Approach for Congruence Class Models

abstract: A Bayesian approach to conduct network model selection is presented for a general class of network models referred to as the congruence class models (CCMs). CCMs form a broad class that includes as special cases several common network models, such as the Erd\H{o}s-R\'{e}nyi-Gilbert model, stochastic block model and many exponential random graph models. Due to the range of models able to be specified as a CCM, investigators are better able to select a model consistent with generative mechanisms associated with the observed network compared to current approaches. In addition, the approach allows for incorporation of prior information. We utilize the proposed Bayesian network model selection approach for CCMs to investigate several mechanisms that may be responsible for the structure of patient-sharing networks, which are associated with the cost and quality of medical care. We found evidence in support of heterogeneity in sociality but not selective mixing by provider type nor degree.


title: Incentive-Compatible Diffusion Auctions

abstract: Diffusion auction is a new model in auction design. It can incentivize the buyers who have already joined in the auction to further diffuse the sale information to others via social relations, whereby both the seller's revenue and the social welfare can be improved. Diffusion auctions are essentially non-typical multidimensional mechanism design problems and agents' social relations are complicatedly involved with their bids. In such auctions, incentive-compatibility (IC) means it is best for every agent to honestly report her valuation and fully diffuse the sale information to all her neighbors. Existing work identified some specific mechanisms for diffusion auctions, while a general theory characterizing all incentive-compatible diffusion auctions is still missing. In this work, we identify a sufficient and necessary condition for all dominant-strategy incentive-compatible (DSIC) diffusion auctions. We formulate the monotonic allocation policies in such multidimensional problems and show that any monotonic allocation policy can be implemented in a DSIC diffusion auction mechanism. Moreover, given any monotonic allocation policy, we obtain the optimal payment policy to maximize the seller's revenue.


title: The quotient Unimodular Vector group is nilpotent

abstract: Jose-Rao introduced and studied the Special Unimodular Vector group $SUm_r(R)$ and $EUm_r(R)$, its Elementary Unimodular Vector subgroup. They proved that for $r \geq 2$, $EUm_r(R)$ is a normal subgroup of $SUm_r(R)$. The Jose-Rao theorem says that the quotient Unimodular Vector group, $SUm_r(R)/EUm_r(R)$, for $r \geq 2$, is a subgroup of the orthogonal quotient group $SO_{2(r+1)}(R)/EO_{2(r + 1)}(R)$. The latter group is known to be nilpotent by the work of Hazrat-Vavilov, following methods of A. Bak; and so is the former. In this article we give a direct proof, following ideas of A. Bak, to show that the quotient Unimodular Vector group is nilpotent of class $\leq d = \dim(R)$. We also use the Quillen-Suslin theory, inspired by A. Bak's method, to prove that if $R = A[X]$, with $A$ a local ring, then the quotient Unimodular Vector group is abelian.


title: The existence of primitive normal elements of quadratic forms over finite fields

abstract: For $q=3^r$ ($r>0$), denote by $\mathbb{F}_q$ the finite field of order $q$ and for a positive integer $m\geq2$, let $\mathbb{F}_{q^m}$ be its extension field of degree $m$. We establish a sufficient condition for existence of a primitive normal element $\alpha$ such that $f(\alpha)$ is a primitive element, where $f(x)= ax^2+bx+c$, with $a,b,c\in \mathbb{F}_{q^m}$ satisfying $b^2\neq ac$ in $\Fm$ except for at most 9 exceptional pairs $(q,m)$.


title: An Atmospheric Cerenkov Telescope Simulation System

abstract: A detailed numerical procedure has been developed to simulate the mechanical configurations and optical properties of Imaging Atmospheric Cerenkov Telescope systems. To test these procedures a few existing ACT arrays are simulated. First results from these simulations are presented.


title: Dietary Restriction of Amino Acids for Cancer Therapy

abstract: Biosyntheses of proteins, nucleotides and fatty acids, are essential for the malignant proliferation and survival of cancer cells. Cumulating research findings show that amino acid restrictions are potential strategies for cancer interventions. Meanwhile, dietary strategies are popular among cancer patients. However, there is still lacking solid rationale to clarify what is the best strategy, why and how it is. Here, integrated analyses and comprehensive summaries for the abundances, signalling and functions of amino acids in proteomes, metabolism, immunity and food compositions, suggest that, intermittent fasting or intermittent dietary lysine restriction with normal maize as an intermittent staple food for days or weeks, might have the value and potential for cancer prevention or therapy. Moreover, dietary supplements were also discussed for cancer cachexia including dietary immunomodulatory.


title: MOEA/D with Random Partial Update Strategy

abstract: Recent studies on resource allocation suggest that some subproblems are more important than others in the context of the MOEA/D, and that focusing on the most relevant ones can consistently improve the performance of that algorithm. These studies share the common characteristic of updating only a fraction of the population at any given iteration of the algorithm. In this work we investigate a new, simpler partial update strategy, in which a random subset of solutions is selected at every iteration. The performance of the MOEA/D using this new resource allocation approach is compared experimentally against that of the standard MOEA/D-DE and the MOEA/D with relative improvement-based resource allocation. The results indicate that using the MOEA/D with this new partial update strategy results in improved HV and IGD values, and a much higher proportion of non-dominated solutions, particularly as the number of updated solutions at every iteration is reduced.


title: Effects of magnon-magnon interactions in a noncollinear magnet induced by combination of a symmetric and an antisymmetric exchange interaction

abstract: We study magnon-magnon interactions and their effects in a spiral magnet induced by combination of an antiferromagnetic Heisenberg interaction and a Dzyaloshinsky-Moriya interaction. We show that the main effect of magnon-magnon interactions on low-energy magnons is to renormalize the coefficient of energy dispersion. This could explain why some experiments are consistent with the noninteracting theory. We also show that although the magnon-magnon interactions induce the pair amplitude for low-energy magnons, its effect on the excitation energy is negligible. This suggests that for magnons the finite pair amplitude does not necessarily accompany the pair condensation.


title: Chemistry and physics of a low-metallicity hot core in the Large Magellanic Cloud

abstract: We present the results of 0.1-pc-scale observations in 250 GHz and 350GHz towards a newly-discovered hot molecular core in a nearby low-metallicity galaxy, the Large Magellanic Cloud (LMC), with the Atacama Large Millimeter/submillimeter Array. A variety of C/N/O/Si/S-bearing molecules are detected towards the high-mass young stellar object, ST16. A rotating protostellar envelope is for the first time detected outside our Galaxy by SO2 and 34SO lines. An outflow cavity is traced by CCH and CN. The isotope abundance of sulfur in the source is estimated to be 32S/34S = 17 and 32S/33S = 53 based on SO, SO2, and CS isotopologues, suggesting that both 34S and 33S are overabundant in the LMC. Rotation diagram analyses show that the source is associated with hot gas (>100K) traced by high-excitation lines of CH3OH and SO2, as well as warm gas (~50K) traced by CH3OH, SO2, 34SO, OCS, CH3CN lines. A comparison of molecular abundances between LMC and Galactic hot cores suggests that organic molecules (e.g., CH3OH, a classical hot core tracer) show a large abundance variation in low metallicity, where the present source is classified into an organic-poor hot core. Our astrochemical simulations suggest that different grain temperature during the initial ice-forming stage would contribute to the chemical differentiation. In contrast, SO2 shows similar abundances within all the known LMC hot cores and the typical abundance roughly scales with the LMC's metallicity. Nitrogen-bearing molecules are generally less abundant in LMC hot cores, except for NO. The present results suggest that chemical compositions of hot cores do not always simply scale with the metallicity.


title: Adaptive Dithering Using Curved Markov-Gaussian Noise in the Quantized Domain for Mapping SDR to HDR Image

abstract: High Dynamic Range (HDR) imaging is gaining increased attention due to its realistic content, for not only regular displays but also smartphones. Before sufficient HDR content is distributed, HDR visualization still relies mostly on converting Standard Dynamic Range (SDR) content. SDR images are often quantized, or bit depth reduced, before SDR-to-HDR conversion, e.g. for video transmission. Quantization can easily lead to banding artefacts. In some computing and/or memory I/O limited environment, the traditional solution using spatial neighborhood information is not feasible. Our method includes noise generation (offline) and noise injection (online), and operates on pixels of the quantized image. We vary the magnitude and structure of the noise pattern adaptively based on the luma of the quantized pixel and the slope of the inverse-tone mapping function. Subjective user evaluations confirm the superior performance of our technique.


title: Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

abstract: The Multi-region Relaxed MHD (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The Stepped Pressure Equilibrium Code (SPEC) [S.R. Hudson et al., Phys. Plasmas 19, 112502 (2012)] was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad-Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Finally, we provide an example equilibrium with a 3D helical field structure as the favoured lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.


title: Versatile Direct Writing of Dopants in a Solid State Host Through Recoil Implantation

abstract: Modifying material properties at the nanoscale is crucially important for devices in nanoelectronics, nanophotonics and quantum information. Optically active defects in wide band gap materials, for instance, are vital constituents for the realisation of quantum technologies. Yet, the introduction of atomic defects through direct ion implantation remains a fundamental challenge. Herein, we establish a universal method for material doping by exploiting one of the most fundamental principles of physics - momentum transfer. As a proof of concept, we direct-write arrays of emitters into diamond via momentum transfer from a Xe+ focused ion beam (FIB) to thin films of the group IV dopants pre-deposited onto a diamond surface. We conclusively show that the technique, which we term knock-on doping, can yield ultra-shallow dopant profiles localized to the top 5 nm of the target surface, and use it to achieve sub-50 nm lateral resolution. The knock-on doping method is cost-effective, yet very versatile, powerful and universally suitable for applications such as electronic and magnetic doping of atomically thin materials and engineering of near-surface states of semiconductor devices.


title: Meso-scale Instability Triggered by Dust Feedback in Dusty Rings: Origin and Observational Implications

abstract: High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations.


title: Photonic nanobeam cavities with nano-pockets for efficient integration of fluorescent nanoparticles

abstract: Integrating fluorescent nanoparticles with high-Q, small mode volume cavities is indispensable for nanophotonics and quantum technologies. To date, nanoparticles have largely been coupled to evanescent fields of cavity modes, which limits the strength of the interaction. Here, we developed both a cavity design and a fabrication method that enable efficient coupling between a fluorescent nanoparticle and a cavity optical mode. The design consists of a fishbone-shaped, one-dimensional photonic crystal cavity with a nano-pocket located at the electric field maximum of the fundamental optical mode. Furthermore, the presence of a nanoparticle inside the pocket reduces the mode volume substantially and induces subwavelength light confinement. Our approach opens exciting pathways to achieve strong light confinement around fluorescent nanoparticles for applications in energy, sensing, lasing and quantum technologies.


title: Deep learning generates custom-made logistic regression models for explaining how breast cancer subtypes are classified

abstract: Differentiating the intrinsic subtypes of breast cancer is crucial for deciding the best treatment strategy. Deep learning can predict the subtypes from genetic information more accurately than conventional statistical methods, but to date, deep learning has not been directly utilized to examine which genes are associated with which subtypes. To clarify the mechanisms embedded in the intrinsic subtypes, we developed an explainable deep learning model called a point-wise linear (PWL) model that generates a custom-made logistic regression for each patient. Logistic regression, which is familiar to both physicians and medical informatics researchers, allows us to analyze the importance of the feature variables, and the PWL model harnesses these practical abilities of logistic regression. In this study, we show that analyzing breast cancer subtypes is clinically beneficial for patients and one of the best ways to validate the capability of the PWL model. First, we trained the PWL model with RNA-seq data to predict PAM50 intrinsic subtypes and applied it to the 41/50 genes of PAM50 through the subtype prediction task. Second, we developed a deep enrichment analysis method to reveal the relationships between the PAM50 subtypes and the copy numbers of breast cancer. Our findings showed that the PWL model utilized genes relevant to the cell cycle-related pathways. These preliminary successes in breast cancer subtype analysis demonstrate the potential of our analysis strategy to clarify the mechanisms underlying breast cancer and improve overall clinical outcomes.


title: The Parallelism Motifs of Genomic Data Analysis

abstract: Genomic data sets are growing dramatically as the cost of sequencing continues to decline and small sequencing devices become available. Enormous community databases store and share this data with the research community, but some of these genomic data analysis problems require large scale computational platforms to meet both the memory and computational requirements. These applications differ from scientific simulations that dominate the workload on high end parallel systems today and place different requirements on programming support, software libraries, and parallel architectural design. For example, they involve irregular communication patterns such as asynchronous updates to shared data structures. We consider several problems in high performance genomics analysis, including alignment, profiling, clustering, and assembly for both single genomes and metagenomes. We identify some of the common computational patterns or motifs that help inform parallelization strategies and compare our motifs to some of the established lists, arguing that at least two key patterns, sorting and hashing, are missing.


title: Connection independent formulation of general relativity

abstract: A connection-independent formulation of general relativity is presented, in which the dynamics does not depend on the choice of connection. The gravity action in this formulation includes one additional scalar term in addition to the Einstein-Hilbert action. No conditions on the connection are imposed. Nevertheless, this formulation yields the Einstein equations, without adding the Gibbons-Hawking-York term even when a manifold has a boundary. Furthermore, this formulation yields a unified description of general relativity, teleparallel gravity, and symmetric teleparallel gravity.


title: Magnetisms of spinor alkali and alkaline-earth atoms in optical lattices

abstract: We theoretically investigate zero-temperature magnetic ordering of mixtures of spin-1 (alkali atoms) and spin-0 (alkaline-earth atoms) bosons in a three-dimensional optical lattice. With the single-mode approximation for the spin-1 bosons, we obtain an effective Bose-Hubbard model for describing the heteronuclear mixtures in optical lattices. By controlling the interspecies interactions between alkali and alkaline-earth atoms, we map out complete phase diagrams of the system with both positive and negative spin-dependent interactions for spin-1 atoms, based on bosonic dynamical mean-field theory. We find that the spin-1 components feature ferromagnetic and nematic insulating phases, in addition to the superfluid, depending on spin-dependent interactions. Take the spin-1 alkali bosons as spin $\uparrow$, and spin-0 alkaline-earth bosons as spin $\downarrow$, we observe that the system favors ferromagnetic insulator at filling $n=1$, and unorder insulator at $n=2$. Interestingly, we observe a two-step Mott-insulating-superfluid phase transition, as a result of mass imbalance between alkali and alkaline-earth atoms.


title: Torsion classes in the equivariant Chow groups of algebraic tori

abstract: We give an example of an algebraic torus $T$ such that the group $\operatorname{CH}^2(BT)_{\operatorname{tors}}$ is non-trivial. This answers a question of Blinstein and Merkurjev.


title: Nonequilibrium superfield and lattice Weyl transform transport approach to quantum Hall effect

abstract: We derive the topological Chern number of the integer quantum Hall effect in electrical conductivity, using Buot's superfield and lattice Weyl transform nonequilibrium quantum transport formalism. The method is naturally straightforward, appropriate for treating nonequilibrium systems acted on by external electromagnetic fields. We have identified the topological invariant in the effective (~p; ~q; E; t)-phase space via the nonequilibrium quantum transport equation, generally not to first-order in electric field but to first-order in the gradient expansion. We have also derive the Kubo current-current correlation for the Hall current as a by-product of our new transport approach. The Berry curvature related to orbital magnetic moment is also calculated.


title: On the directions determined by a Cartesian product in an affine Galois plane

abstract: We prove that the number of directions contained in a set of the form $A \times B \subset AG(2,p)$, where $p$ is prime, is at least $|A||B| - \min\{|A|,|B|\} + 2$. Here $A$ and $B$ are subsets of $GF(p)$ each with at least two elements and $|A||B| <p$. This bound is tight for an infinite class of examples. Our main tool is the use of the R\'edei polynomial with Sz\H{o}nyi's extension. As an application of our main result, we obtain an upper bound on the clique number of a Paley graph, matching the current best bound obtained recently by Hanson and Petridis.


title: Nov\'{a}k's conjecture on cyclic Steiner triple systems and its generalization

abstract: Nov\'{a}k conjectured in 1974 that for any cyclic Steiner triple systems of order $v$ with $v\equiv 1\pmod{6}$, it is always possible to choose one block from each block orbit so that the chosen blocks are pairwise disjoint. We consider the generalization of this conjecture to cyclic $(v,k,\lambda)$-designs with $1 \leq \lambda \leq k-1$. Superimposing multiple copies of a cyclic symmetric design shows that the generalization cannot hold for all $v$, but we conjecture that it holds whenever $v$ is sufficiently large compared to $k$. We confirm that the generalization of the conjecture holds when $v$ is prime and $\lambda=1$ and also when $\lambda \leq (k-1)/2$ and $v$ is sufficiently large compared to $k$. As a corollary, we show that for any $k \geq 3$, with the possible exception of finitely many composite orders $v$, every cyclic $(v,k,1)$-design without short orbits is generated by a $(v,k,1)$-disjoint difference family.


title: The impact of anisotropy on ITER scenarios

abstract: We report on the impact of anisotropy to tokamak plasma configuration and stability. Our focus is on analysis of the impact of anisotropy on ITER pre-fusion power operation 5~MA, $B=1.8$~T ICRH scenarios. To model ITER scenarios remapping tools are developed to distinguish the impact of pressure anisotropy from the change in magnetic geometry caused by an anisotropy-modified current profile. The remappings iterate the anisotropy-modified current profile to produce the same $q$ profile with matched thermal energy. The analysis is a step toward equilibria that are kinetically self-consistent for a prescribed scenario. We find characteristic detachment of flux surfaces from pressure surfaces, and an outboard (inboard) shift of peak density for $T_{\parallel}>T_\perp$ ( $T_{\parallel}<T_\perp$). Differences in the poloidal current profile are evident, albeit not as pronounced as for the spherical tokamak. We find that the incompressional continuum is largely unchanged in the presence of anisotropy, and the mode structure of gap modes is largely unchanged. The compressional branch however exhibits significant differences in the continuum. We report on the implication of these modifications.


title: The redshift and star formation mode of AzTEC2: a pair of massive galaxies at $z=4.63$

abstract: We combine observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) to assess the redshift and to study the star formation conditions in AzTEC2: one of the brightest sub-millimeter galaxies (SMGs) in the COSMOS field ($S_{\rm 1.1mm}=10.5\pm1.4$mJy). Our high-resolution observations confirm that AzTEC2 splits into two components (namely AzTEC2-A and AzTEC2-B) for which we detect [C$\,$II] and $^{12}$CO(5$\to$4) line emission, implying a redshift of $4.626\pm0.001$ ($4.633\pm0.001$) for AzTEC2-A (AzTEC2-B) and ruling out previous associations with a galaxy at $z\sim1$. We use the $^{12}$CO(5$\to$4) line emission and adopt typical SMG-like gas excitation conditions to estimate the molecular gas mass, which is $M_{\rm gas}(\alpha_{\rm CO}/2.5)=2.1\pm0.4 \times10^{11}{\rm M}_\odot$ for AzTEC2-A, and a factor four lower for AzTEC2-B. With the infrared-derived star formation rate of AzTEC2-A ($1920\pm100 \,M_\odot{\rm \, yr}^{-1}$) and AzTEC2-B ($710\pm 35\,M_\odot{\rm \,yr}^{-1}$), they both will consume their current gas reservoir within $(30-200)\,$Myr. We find evidence of a rotation-dominated [C$\,$II] disk in AzTEC2-A, with a de-projected rotational velocity of $v_{\rm rot}(i=39^\circ)=660\pm130{\rm \,km\,s}^{-1}$, velocity dispersion $\lesssim100{\rm \,km\,s}^{-1}$, and dynamical mass of $M_{\rm dyn}(i=39^\circ)=2.6^{+1.2}_{-0.9}\times10^{11}\,M_\odot$. We propose that an elevated gas accretion rate from the cosmic web might be the main driver of the intense levels of star formation in AzTEC2-A, which might be further enhanced by gravitational torques induced by its minor companion (AzTEC2-B). These results strengthen the picture whereby the population of single-dish selected SMGs is rather heterogeneous, including a population of pairs of massive, highly-active galaxies in a pre-coalescence phase.


title: Convergence rate analysis of a sequential convex programming method with line search for a class of constrained difference-of-convex optimization problems

abstract: In this paper, we study the sequential convex programming method with monotone line search (SCP$_{ls}$) in [46] for a class of difference-of-convex (DC) optimization problems with multiple smooth inequality constraints. The SCP$_{ls}$ is a representative variant of moving-ball-approximation-type algorithms [6,10,13,54] for constrained optimization problems. We analyze the convergence rate of the sequence generated by SCP$_{ls}$ in both nonconvex and convex settings by imposing suitable Kurdyka-Lojasiewicz (KL) assumptions. Specifically, in the nonconvex settings, we assume that a special potential function related to the objective and the constraints is a KL function, while in the convex settings we impose KL assumptions directly on the extended objective function (i.e., sum of the objective and the indicator function of the constraint set). A relationship between these two different KL assumptions is established in the convex settings under additional differentiability assumptions. We also discuss how to deduce the KL exponent of the extended objective function from its Lagrangian in the convex settings, under additional assumptions on the constraint functions. Thanks to this result, the extended objectives of some constrained optimization models such as minimizing $\ell_1$ subject to logistic/Poisson loss are found to be KL functions with exponent $\frac12$ under mild assumptions. To illustrate how our results can be applied, we consider SCP$_{ls}$ for minimizing $\ell_{1-2}$ [60] subject to residual error measured by $\ell_2$ norm/Lorentzian norm [21]. We first discuss how the various conditions required in our analysis can be verified, and then perform numerical experiments to illustrate the convergence behaviors of SCP$_{ls}$.


title: Radius of Starlikeness of Certain Analytic Functions

abstract: This paper studies analytic functions $f$ defined on the open unit disk of the complex plane for which $f/g$ and $(1+z)g/z$ are both functions with positive real part for some analytic function $g$. We determine radius constants of these functions to belong to classes of strong starlike functions, starlike functions of order $\alpha$, parabolic starlike functions, as well as to the classes of starlike functions associated with lemniscate of Bernoulli, cardioid, lune, reverse lemniscate, sine function, exponential function and a particular rational function. The results obtained are sharp.


title: Contract-connection:An efficient communication protocol for Distributed Ledger Technology

abstract: Distributed Ledger Technology (DLT) is promising to become the foundation of many decentralised systems. However, the unbalanced and unregulated network layout contributes to the inefficiency of DLT especially in the Internet of Things (IoT) environments, where nodes connect to only a limited number of peers. The data communication speed globally is unbalanced and does not live up to the constraints of efficient real-time distributed systems. In this paper, we introduce a new communication protocol, which enables nodes to calculate the tradeoff between connecting/disconnecting a peer in a completely decentralised manner. The network layout globally is continuously re-balancing and optimising along with nodes adjusting their peers. This communication protocol weakened the inequality of the communication network. The experiment suggests this communication protocol is stable and efficient.


title: The Difference and Unity of Irregular LQ Control and Standard LQ Control and Its Solution

abstract: Irregular linear quadratic control (LQ, was called Singular LQ) has been a long-standing problem since 1970s. This paper will show that an irregular LQ control (deterministic) is solvable (for arbitrary initial value) if and only if the LQ cost can be rewritten as a regular one by changing the terminal cost $x'(T)Hx(T)$ to $x'(T)[H+P_1(T)]x(T)$, while the optimal controller can achieve $P_1(T)x(T)=0$ at the same time. In other words, the irregular controller (if exists) needs to do two things at the same time, one thing is to minimize the cost and the other is to achieve the terminal constraint $P_1(T)x(T)=0$, which clarifies the essential difference of irregular LQ from the standard LQ control where the controller is to minimize the cost only. With this breakthrough, we further study the irregular LQ control for stochastic systems with multiplicative noise. A sufficient solving condition and the optimal controller is presented based on Riccati equations.


title: An Efficient Framework for Automated Screening of Clinically Significant Macular Edema

abstract: The present study proposes a new approach to automated screening of Clinically Significant Macular Edema (CSME) and addresses two major challenges associated with such screenings, i.e., exudate segmentation and imbalanced datasets. The proposed approach replaces the conventional exudate segmentation based feature extraction by combining a pre-trained deep neural network with meta-heuristic feature selection. A feature space over-sampling technique is being used to overcome the effects of skewed datasets and the screening is accomplished by a k-NN based classifier. The role of each data-processing step (e.g., class balancing, feature selection) and the effects of limiting the region-of-interest to fovea on the classification performance are critically analyzed. Finally, the selection and implication of operating point on Receiver Operating Characteristic curve are discussed. The results of this study convincingly demonstrate that by following these fundamental practices of machine learning, a basic k-NN based classifier could effectively accomplish the CSME screening.


title: Strategy-Proof Spectrum Allocation among Multiple Operators for Demand Varying Wireless Networks

abstract: To address the exponentially increasing data rate demands of end users, necessitates efficient spectrum allocation among co-existing operators in licensed and unlicensed spectrum bands to cater to the temporal and spatial variations of traffic in the wireless network. In this paper, we address the spectrum allocation problem among non-cooperative operators via auctions. The classical Vickrey-Clarke-Groves (VCG) approach provides the framework for a strategy-proof and social welfare maximizing auction at high computational complexity, which makes it infeasible for practical implementation. We propose sealed bid auction mechanisms for spectrum allocation which are computationally tractable and hence applicable for allocating spectrum by performing auctions in short durations as per the dynamic load variations of the network. We establish that the proposed algorithm is strategy-proof for uniform demand. Furthermore, for non-uniform demand we propose an algorithm that satisfies weak strategy-proofness. We also consider non-linear increase in the marginal valuations with demand. Simulation results are presented to exhibit the performance comparison of the proposed algorithms with VCG and other existing mechanisms.


title: Bicomplex frames

abstract: The main purpose is to introduce the so-called bicomplex (bc)-frames which is a special extension to bicomplex infinite Hilbert spaces of the classical frames. The crucial result is the characterization of bc-frames in terms of their idempotent components, giving rise to generalization of certain results to bc-frames. Although the extension is natural, many basic properties satisfied by classical frames do not remain valid for bc-frames, unless we restrict ourself to complex-valued Hilbert space on bicomplex numbers. By benefiting from insight provided by the classical frame theory, we discuss the construction of bc-frame operator and Weyl--Heisenberg bc-frames and we provide some new ones which are appropriate for bc-frames.


title: TPFA Finite Volume Approximation of Wasserstein Gradient Flows

abstract: Numerous infinite dimensional dynamical systems arising in different fields have been shown to exhibit a gradient flow structure in the Wasserstein space. We construct Two Point Flux Approximation Finite Volume schemes discretizing such problems which preserve the variational structure and have second order accuracy in space. We propose an interior point method to solve the discrete variational problem, providing an efficient and robust algorithm. We present two applications to test the scheme and show its order of convergence.


title: Distributed State Estimation over Time-Varying Graphs: Exploiting the Age-of-Information

abstract: We study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. In fact, we also establish finite-time convergence based on our approach. Finally, we develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.


title: On the global shape of convex functions on locally convex spaces

abstract: In the recent paper \cite{Aza:19} D Azagra studies the global shape of continuous convex functions defined on a Banach space $X$. More precisely, when $X$ is separable, it is shown that for every continuous convex function $f:X\rightarrow\mathbb{R}$ there exist a unique closed linear subspace $Y$ of $X$, a continuous function $h:X/Y\rightarrow\mathbb{R}$ with the property that $\lim_{t\rightarrow\infty}h(u+tv)=\infty$ for all $u,v\in X/Y$, $v\neq0$, and $x^{\ast}\in X^{\ast}$ such that $f=h\circ\pi+x^{\ast}$, where $\pi :X\rightarrow X/Y$ is the natural projection. Our aim is to characterize those proper lower semi\-continuous convex functions defined on a locally convex space which have the above representation. In particular, we show that the continuity of the function $f$ and the completeness of $X$ can be removed from the hypothesis of Azagra's theorem.


title: The topological magnon bands in the Flux state in Sashtry-Sutherland lattice

abstract: We investigate low energy magnon excitations above the non-collinear flux state and non-coplanar canted flux state in a Heisenberg anti-ferromagnet with Dzyaloshinskii-Moriya interaction~(DMI) on a Sashtry-Sutherland lattice.While previous studies have shown the presence of topological magnetic excitation in the dimer and ferromagnetic phases on the Shastry-Sutherland lattice, our results establish the non-trivial topology of magnons in the anti-ferromagnetic flux and canted flux states. Our results uncover the existence of a multitude of topological phase transitions in the magnon sector -- evidenced by the changing Chern numbers of the single magnon bands -- as the Hamiltonian parameters are varied, even when the ground state remains unchanged. The thermal Hall conductivity is calculated and its derivative is shown to exhibit a logarithmic divergence at the phase transitions, independent of the type of band touching involved. This may provide a useful means to identify the energy at which the transition occurs. Finally, we propose the way to realize the studied model in a practical material.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.other


title: Global climate modeling of Saturn's atmosphere. Part IV: stratospheric equatorial oscillation

abstract: The Composite InfraRed Spectrometer (CIRS) on board Cassini revealed an equatorial oscillation of stratospheric temperature, reminiscent of the Earth's Quasi-Biennial Oscillation (QBO), as well as anomalously high temperatures under Saturn's rings. To better understand these predominant features of Saturn's atmospheric circulation in the stratosphere, we have extended towards higher altitudes the DYNAMICO-Saturn global climate model (GCM), already used in a previous publication to study the tropospheric dynamics, jets formation and planetary-scale waves activity. Firstly, we study the higher model top impact on the tropospheric zonal jets and kinetic energy distribution. Raising the model top prevents energy and enstrophy accumulation at tropopause levels. The reference GCM simulation with 1/2$^{\circ}$ latitude/longitude resolution and a raised model top exhibits a QBO-like oscillation produced by resolved planetary-scale waves. However, the period is more irregular and the downward propagation faster than observations. Furthermore, compared to the CIRS observation retrievals, the modeled QBO-like oscillation underestimates by half both the amplitude of temperature anomalies at the equator and the vertical characteristic length of this equatorial oscillation. This QBO-like oscillation is mainly driven by westward-propagating waves; a significant lack of eastward wave-forcing explains a fluctuating eastward phase of the QBO-like oscillation. At 20$^{\circ}$N and 20$^{\circ}$S latitudes, the DYNAMICO-Saturn GCM exhibits several strong seasonal eastward jets, alternatively in the northern and southern hemisphere. These jets are correlated with the rings' shadowing. Using a GCM simulation without rings' shadowing, we show its impact on Saturn's stratospheric dynamics. Both residual-mean circulation and eddy forcing are impacted by rings' shadowing.


title: The Laplacian on some self-conformal fractals and Weyl's asymptotics for its eigenvalues: A survey of the analytic aspects

abstract: This article surveys the analytic aspects of the author's recent studies on the construction and analysis of a "geometrically canonical" Laplacian on circle packing fractals invariant with respect to certain Kleinian groups (i.e., discrete groups of M\"{o}bius transformations on the Riemann sphere $\widehat{\mathbb{C}}=\mathbb{C}\cup\{\infty\}$), including the classical Apollonian gasket and some round Sierpi\'{n}ski carpets. The main result on Weyl's asymptotics for its eigenvalues is of the same form as that by Oh and Shah [Invent. Math. 187 (2012), 1--35, Theorem 1.4] on the asymptotic distribution of the circles in a very large class of such fractals.


title: Approximation Algorithms and LP Relaxations for Scheduling Problems Related to Min-Sum Set Cover

abstract: We consider single-machine scheduling problems that are natural generalizations or variations of the min-sum set cover problem and the min-sum vertex cover problem. For each of these problems, we give new approximation algorithms. Some of these algorithms rely on time-indexed LP relaxations. We show how a variant of alpha-point scheduling leads to the best-known approximation ratios, including a guarantee of 4 for an interesting special case of the so-called generalized min-sum set cover problem. We also make explicit the connection between the greedy algorithm for min-sum set cover and the concept of Sidney decomposition for precedence-constrained single-machine scheduling, and show how this leads to a 4-approximation algorithm for single-machine scheduling with so-called bipartite OR-precedence constraints.


title: Gauss-Bonnet Inflation after Planck2018

abstract: We study the primordial perturbations and reheating process in the models where the Gauss-Bonnet term is non-minimally coupled to the canonical and non-canonical (DBI and tachyon) scalar fields. We consider several potentials and Gauss-Bonnet coupling terms as power-law, dilaton-like, $\cosh$-type, E-model and T-model. To seek the observational viability of these models, we study the scalar perturbations numerically and compare the results with the Planck2018 TT, TE, EE+lowE+lensing+BK14+BAO joint data at $68\%$ CL and $95\%$ CL. We also study the tensor perturbations in confrontation with the Planck2018 TT, TE, EE+lowE+lensing+BK14+BAO+ LIGO$\&$Virgo2016 joint data at $68\%$ CL and $95\%$ CL. In this regard, we obtain some constraints on the Gauss-Bonnet coupling parameter $\beta$. Another important process in the early universe is the reheating phase after inflation which is necessary to reheat the universe for subsequent evolution. In this regard, we study the reheating process in these models and find some expressions for the e-folds number and temperature during that era. Considering that from Planck TT,TE,EE+lowEB+lensing data and BICEP2/Keck Array 2014, based on the $\Lambda$CDM$+r+\frac{dn_{s}}{d\ln k}$ model, we have $n_{s}=0.9658\pm 0.0038$ and $r<0.072$, we obtain some constraints on the e-folds number and temperature. From the values of the e-folds number and the effective equation of state and also the observationally viable value of the scalar spectral index, we explore the capability of the models in explaining the reheating phase.


title: Some inequalities for Chebyshev polynomials

abstract: Askey and Gasper (1976) proved a trigonometric inequality which improves another trigonometric inequality found by M. S. Robertson (1945). Here these inequalities are reformulated in terms of the Chebyshev polynomial of the first kind $T_n$ and then put into a one-parametric family of inequalities. The extreme value of the parameter is found for which these inequalities hold true. As a step towards the proof of this result we establish the following complement to the finite increment theorem specialized to $T_n^{\prime}$: $$ T_n^{\prime}(1)-T_n^{\prime}(x)\geq (1-x)\,T_n^{\prime\prime}(x)\,,\qquad x\in [0,1]\,. $$ By a known expansion formula, this property is extended for the class of ultraspherical polynomials $P_n^{(\lambda)}$, $\lambda\geq 1$.


title: eXTP perspectives for the $\nu$MSM sterile neutrino dark matter model

abstract: We discuss the potential of the eXTP X-ray telescope, in particular its Spectroscopic Focusing Array (SFA), Large Area Detector (LAD) and Wide Field Monitor (WFM) for the detection of a signal from keV-scale decaying dark matter. We show that the sensitivity of the eXTP is sufficient to improve existing constraints on the mixing angle of the neutrino Minimal extension of the Standard Model ($\nu$MSM) by a factor of 5-10 within the dark matter mass range 2-50 keV, assuming 1% level of systematic uncertainty. We assert that the eXTP will be able to probe previously inaccessible range of $\nu$MSM parameters and serve as a precursor for the Athena mission in decaying dark matter searches.


title: Operator product expansion coefficients in the exact renormalization group formalism

abstract: We study how to compute the operator product expansion coefficients in the exact renormalization group formalism. After discussing possible strategies, we consider some examples explicitly, within the $\epsilon$-expansions, for the Wilson-Fisher fixed points of the real scalar theory in $d=4-\epsilon$ dimensions and the Lee-Yang model in $d=6-\epsilon$ dimensions. Finally we discuss how our formalism may be extended beyond perturbation theory.


title: BlockHouse: Blockchain-based Distributed Storehouse System

abstract: We propose in this paper BlockHouse, a decentralized/P2P storage system fully based on private blockchains. Each participant can rent his unused storage in order to host data of other members. This system uses a dual Smart Contract and Proof of Retrievability system to automatically check at a fixed frequency if the file is still hosted. In addition to transparency, the blockchain allows a better integration with all payments associated to this type of system ( regular payments, sequestration to ensure good behaviors of users, ...). Except the data transferred between the client and the server, all the actions go through a smart contract in the blockchain in order to log, pay and secure the entire storage process.


title: The Thue-Morse and Rudin-Shapiro sequences at primes in principal number fields

abstract: We consider a numeration system in the ring of integers ${\mathcal O}_K$ of a number field, which we assume to be principal. We prove that the property of being a prime in ${\mathcal O}_K$ is decorrelated from two fundamental examples of automatic sequences relative to the chosen numeration system: the Thue-Morse and the Rudin-Shapiro sequences. This is an analogue, in ${\mathcal O}_K$, of results of Mauduit-Rivat which were concerned with the case $K={\mathbb Q}$.


title: Rest frame vacuum of the Proca field on the de Sitter expanding universe

abstract: The general plane wave solutions of the Proca field in conformal charts of the de Sitter expanding universe are derived for arbitrary polarizations showing how the frequencies can be separated in rest frames, defining thus the rest frame vacuum of this field.


title: Numerical simulation of ionospheric depletions resulting from rocket launches using a general circulation model

abstract: Rocket exhaust plumes have been observed to cause large-scale depletions of ionospheric plasmas ("ionospheric holes"). In the F-region, charge exchange reactions occur between O^+ ions and exhaust species such as H$_2$O, H$_2$, and CO$_2$ to form ions which then undergo rapid dissociative recombination. The Global Ionosphere-Thermosphere Model (GITM) was extended to include these chemical reactions and appropriate source terms to represent rocket exhaust plumes. The resulting model was applied to ionospheric depletions resulting from the launches of Jason-3 and FORMOSAT-5 on SpaceX Falcon 9 rockets from Vandenberg Air Force Base. Outputs from the model were compared with GNSS, ionosonde, and satellite Langmuir probe measurements. Simulation indicated that the FORMOSAT-5 launch resulted in a far larger and longer-lived ionospheric depletion than the Jason-3 launch, consistent with observations.


title: Coded Caching with Polynomial Subpacketization

abstract: Consider a centralized caching network with a single server and $K$ users. The server has a database of $N$ files with each file being divided into $F$ packets ($F$ is known as subpacketization), and each user owns a local cache that can store $\frac{M}{N}$ fraction of the $N$ files. We construct a family of centralized coded caching schemes with polynomial subpacketization. Specifically, given $M$, $N$ and an integer $n\geq 0$, we construct a family of coded caching schemes for any $(K,M,N)$ caching system with $F=O(K^{n+1})$. More generally, for any $t\in\{1,2,\cdots,K-2\}$ and any integer $n$ such that $0\leq n\leq t$, we construct a coded caching scheme with $\frac{M}{N}=\frac{t}{K}$ and $F\leq K\binom{\left(1-\frac{M}{N}\right)K+n}{n}$.


title: The viable f(G) gravity models via reconstruction from the observations

abstract: We reconstruct the viable f(G) gravity models from the observations and provide the analytic solutions that well describe our numerical results. In order to avoid unphysical challenges that occur during the numerical reconstruction, we generalize f(G) models into f(GA), which is the simple extension of f(G) models with the introduction of a constant A parameter. We employ several observational data together with the stability condition, which reads d2f/dG2 > 0 and must be satisfied in the late-time evolution of the universe, to give proper initial conditions for solving the perturbation equation. As a result, we obtain the analytic functions that match the numerical solutions. Furthermore, it might be interesting if one can find the physical origin of those analytic solutions and its cosmological implications.


title: BAASH: Enabling Blockchain-as-a-Service on High-Performance Computing Systems

abstract: The state-of-the-art approach to manage blockchains is to process blocks of transactions in a shared-nothing environment. Although blockchains have the potential to provide various services for high-performance computing (HPC) systems, HPC will not be able to embrace blockchains before the following two missing pieces become available: (i) new consensus protocols being aware of the shared-storage architecture in HPC, and (ii) new fault-tolerant mechanisms compensating for HPC's programming model---the message passing interface (MPI)---that is vulnerable for blockchain-like workloads. To this end, we design a new set of consensus protocols crafted for the HPC platforms and a new fault-tolerance subsystem compensating for the failures caused by faulty MPI processes. Built on top of the new protocols and fault-tolerance mechanism, a prototype system is implemented and evaluated with two million transactions on a 500-core HPC cluster, showing $6\times$, $12\times$, and $75\times$ higher throughput than Hyperldeger, Ethereum, and Parity, respectively.


title: Segment blockchain: A size reduced storage mechanism for blockchain

abstract: The exponential growth of the blockchain size has become a major contributing factor that hinders the decentralisation of blockchain and its potential implementations in data-heavy applications. In this paper, we propose segment blockchain, an approach that segmentises blockchain and enables nodes to only store a copy of one blockchain segment. We use \emph{PoW} as a membership threshold to limit the number of nodes taken by an Adversary---the Adversary can only gain at most $n/2$ of nodes in a network of $n$ nodes when it has $50\%$ of the calculation power in the system (the Nakamoto blockchain security threshold). A segment blockchain system fails when an Adversary stores all copies of a segment, because the Adversary can then leave the system, causing a permanent loss of the segment. We theoretically prove that segment blockchain can sustain a $(AD/n)^m$ failure probability when the Adversary has no more than $AD$ number of nodes and every segment is stored by $m$ number of nodes. The storage requirement is mostly shrunken compared to the traditional design and therefore making the blockchain more suitable for data-heavy applications.