title: On the Uniqueness of Binary Quantizers for Maximizing Mutual Information

abstract: We consider a channel with a binary input X being corrupted by a continuous-valued noise that results in a continuous-valued output Y. An optimal binary quantizer is used to quantize the continuous-valued output Y to the final binary output Z to maximize the mutual information I(X; Z). We show that when the ratio of the channel conditional density r(y) = P(Y=y|X=0)/ P(Y =y|X=1) is a strictly increasing/decreasing function of y, then a quantizer having a single threshold can maximize mutual information. Furthermore, we show that an optimal quantizer (possibly with multiple thresholds) is the one with the thresholding vector whose elements are all the solutions of r(y) = r* for some constant r* > 0. Interestingly, the optimal constant r* is unique. This uniqueness property allows for fast algorithmic implementation such as a bisection algorithm to find the optimal quantizer. Our results also confirm some previous results using alternative elementary proofs. We show some numerical examples of applying our results to channels with additive Gaussian noises.


title: Towards Secure IoT Communication with Smart Contracts in a Blockchain Infrastructure

abstract: The Internet of Things (IoT) is undergoing rapid growth in the IT industry, but, it continues to be associated with several security and privacy concerns as a result of its massive scale, decentralised topology, and resource-constrained devices. Blockchain (BC), a distributed ledger technology used in cryptocurrency has attracted significant attention in the realm of IoT security and privacy. However, adopting BC to IoT is not straightforward in most cases, due to overheads and delays caused by BC operations. In this paper, we apply a BC technology known as Hyperledgder Fabric, to an IoT network. This technology introduces an execute-order technique for transactions that separates the transaction execution from consensus, resulting in increased efficiency. We demonstrate that our proposed IoT-BC architecture is sufficiently secure with regard to fundamental security goals i.e., confidentiality, integrity, and availability. Finally, the simulation results are highlighted that shows the performance overheads associated with our approach are as minimal as those associated with the Hyperledger Fabric framework and negligible in terms of security and privacy.


title: Network-Based Analysis of Public Transportation Systems in North American Cities

abstract: A comprehensive data analysis system is implemented for the extraction of information and comparison of North American public transport systems. The system is based on network representations of the transport systems and makes use of a span of metrics and algorithms from the established properties in graph theory to complicated domain specific measurements. Due to nature of big data systems and the requirement of scalability, many heuristic optimizations and approximations have been considered in the system. Integration with other sources of data specially population density maps is also executed in the system. Formal evaluations are done on subcomponents of the system to make sure the approximations have reasonable precision. Results on comparison of four cities, San Francisco, Boston, Toronto and Los Angeles, approves that the big data approach to comparison of public transit systems can successfully reveal the underlying similarities and differences.


title: Variability in the Massive Open Cluster NGC 1817 from K2: A Rich Population of Asteroseismic Red Clump, Eclipsing Binary, and Main Sequence Pulsating Stars

abstract: We present a survey of variable stars detected in K2 Campaign 13 within the massive intermediate age ($\sim1$ Gyr) open cluster NGC 1817. We identify a complete sample of 44 red clump stars in the cluster, and have measured asteroseismic quantities ($\nu_{\rm max}$ and/or $\Delta \nu$) for 29 of them. Five stars showed suppressed dipole modes, and the occurrence rates indicate that mode suppression is unaffected by evolution through core helium burning. A subset of the giants in NGC 1817 (and in the similarly aged cluster NGC 6811) have $\nu_{\rm max}$ and $\Delta \nu$ values at or near the maximum observed for core helium burning stars, indicating they have core masses near the minimum for fully non-degenerate helium ignition. Further asteroseismic study of these stars can constrain the minimum helium core mass in red clump stars and the physics that determines this limit. Two giant stars show photometric variations on timescales similar to previously measured spectroscopic orbits. Thirteen systems in the field show eclipses, but only five are probable cluster members. We identify 32 $\delta$ Sct pulsators, 27 $\gamma$ Dor candidates, and 7 hybrids that are probable cluster members, with most new detections. We used the ensemble properties of the $\delta$ Sct stars to identify stars with possible radial pulsation modes. Among the oddities we have uncovered are: an eccentric orbit for a short-period binary containing a $\delta$ Sct pulsating star; a rare subgiant within the Hertzsprung gap showing $\delta$ Sct pulsations; and two hot $\gamma$ Dor pulsating star candidates.


title: A Conceptual Paper on SERVQUAL-Framework for Assessing Quality of Internet of Things (IoT) Services

abstract: Service quality possesses the vital prominence in usability of innovative products and services. As technological innovation has made the life synchronized and effective, Internet of Things (IoT) is matter of discussion everywhere. From users' perspective, IoT services are always embraced by various system characteristics of security and performance. A service quality model can better present the preference of such technology customers. the study intends to project theoretical model of service quality for internet of things (IoT). Based on the existing models of service quality and the literature in internet of things, a framework is proposed to conceptualize and measure service quality for internet of things.This study established the IoT-Servqual model with four dimensions (i.e., Privacy, Functionality, Efficiency, and Tangibility) of multiple service quality models. These dimensions are essential and inclined towards the users' leaning of IoT Services. This paper contributes to research on internet of things services by development of a comprehensive framework for customers' quality apprehension. This model will previse the expression of information secrecy of users related with internet of things (IoT). This research will advance understanding of service quality in modern day technology and assist firms to devise the fruitful service structure.


title: Towards a secure behavior modeling for IoT networks using Blockchain

abstract: Internet of Things (IoT) occupies a vital aspect of our everyday lives. IoT networks composed of smart-devices which communicate and transfer the information without the physical intervention of humans. Due to such proliferation and autonomous nature of IoT systems make these devices threatened and prone to a severe kind of threats. In this paper, we introduces a behavior capturing, and verification procedures in blockchain supported smart-IoT systems that can be able to show the trust-level confidence to outside networks. We defined a custom \emph{Behavior Monitor} and implement on a selected node that can extract the activity of each device and analyzes the behavior using deep machine learning strategy. Besides, we deploy Trusted Execution Technology (TEE) which can be used to provide a secure execution environment (enclave) for sensitive application code and data on the blockchain. Finally, in the evaluation phase we analyze various IoT devices data that is infected by Mirai attack. The evaluation results show the strength of our proposed method in terms of accuracy and time required for detection.


title: Single-bit Quantization Capacity of Binary-input Continuous-output Channels

abstract: We consider a channel with discrete binary input X that is corrupted by a given continuous noise to produce a continuous-valued output Y. A quantizer is then used to quantize the continuous-valued output Y to the final binary output Z. The goal is to design an optimal quantizer Q* and also find the optimal input distribution p*(X) that maximizes the mutual information I(X; Z) between the binary input and the binary quantized output. A linear time complexity searching procedure is proposed. Based on the properties of the optimal quantizer and the optimal input distribution, we reduced the searching range that results in a faster implementation algorithm. Both theoretical and numerical results are provided to illustrate our method.


title: Quantum signatures of transitions from stable fixed points to limit cycles in optomechanical systems

abstract: Optomechanical systems, due to its inherent nonlinear optomechanical coupling, owns rich nonlinear dynamics of different types of motion. The interesting question is that whether there exist some common quantum features to infer the nonlinear dynamical transitions from one type to another. In this paper, we have studied the quantum signatures of transitions from stable fixed points to limit cycles in an optomechanical phonon laser system. Our calculations show that the entanglement of stable fixed points in the long run does not change with time, however, it will oscillate periodically with time at the mechanical vibration frequency for the limit cycles. Most strikingly, the entanglement quite close to the boundary line keeps as a constant, and it is very robust to the thermal phonon noise, as strong indications of this particular classical transitions.


title: Characterization of Sr2RuO4 Josephson junctions made of epitaxial films

abstract: We have studied fundamental properties of weak-link Sr2RuO4/Sr2RuO4 Josephson junctions fabricated by making a narrow constriction on superconducting Sr2RuO4 films through laser micro-patterning. The junctions show a typical overdamped behavior with much higher critical current density, compared with those previously reported for bulk Sr2RuO4/s-wave superconductor junctions. Observed magnetic field and temperature dependences of the Josephson critical current suggest that the chiral p-wave is unlikely for the superconducting symmetry, encouraging further theoretical calculations of the Sr2RuO4/Sr2RuO4 type junctions.


title: Permutation Enhances Classical Communication Assisted by Entangled States

abstract: We give a capacity formula for the classical communication over a noisy quantum channel, when local operations and global permutations allowed in the encoding and bipartite states preshared between the sender and the receiver. The two endpoints of this formula are the Holevo capacity (without entanglement assistance) and the entanglement-assisted capacity (with unlimited entanglement assistance). What's more, we show that the capacity satisfies the strong converse property and thus the formula serves as a sharp dividing line between achievable and unachievable rates of communication. We prove that the difference between the assisted capacity and the Holevo capacity is upper bounded by the discord of formation of the preshared state. As examples, we derive analytically the classical capacity of various quantum channels of interests. Our result witnesses the power of random permutation in classical communication, whenever entanglement assistance is available.


title: Dirac neutrinos from Peccei-Quinn symmetry: two examples

abstract: We aim to explain the nature of neutrinos using Peccei-Quinn symmetry. We discuss two simple scenarios, one based on a type-II Dirac seesaw and the other in a one-loop neutrino mass generation, which solve the strong CP problem and naturally lead to Dirac neutrinos. In the first setup latest neutrino mass limit gives rise to axion which is in the reach of conventional searches. Moreover, we have both axion as well as WIMP dark mater for our second set up.


title: On Bounds and Closed Form Expressions for Capacities of Discrete Memoryless Channels with Invertible Positive Matrices

abstract: While capacities of discrete memoryless channels are well studied, it is still not possible to obtain a closed-form expression for the capacity of an arbitrary discrete memoryless channel. This paper describes an elementary technique based on Karush Kuhn Tucker (KKT) conditions to obtain (1) a good upper bound of a discrete memoryless channel having an invertible positive channel matrix and (2) a closed-form expression for the capacity if the channel matrix satisfies certain conditions related to its singular value and its Gershgorin disk.


title: Towards Practical Encrypted Network Traffic Pattern Matching for Secure Middleboxes

abstract: Network Function Virtualisation (NFV) advances the adoption of composable software middleboxes. Accordingly, cloud data centres become major NFV vendors for enterprise traffic processing. Due to the privacy concern of traffic redirection to the cloud, secure middlebox systems (e.g., BlindBox) draw much attention; they can process encrypted packets against encrypted rules directly. However, most of the existing systems supporting pattern matching based network functions require the enterprise gateway to tokenise packet payloads via sliding windows. Such tokenisation induces a considerable communication overhead, which can be over 100$\times$ to the packet size. To overcome this bottleneck, in this paper, we propose the first bandwidth-efficient encrypted pattern matching protocol for secure middleboxes. We resort to a primitive called symmetric hidden vector encryption (SHVE), and propose a variant of it, aka SHVE+, to achieve constant and moderate communication cost. To speed up, we devise encrypted filters to reduce the number of accesses to SHVE+ during matching highly. We formalise the security of our proposed protocol and conduct comprehensive evaluations over real-world rulesets and traffic dumps. The results show that our design can inspect a packet over 20k rules within 100 $\mu$s. Compared to prior work, it brings a saving of $94\%$ in bandwidth consumption.


title: Exploring two-neutron halo formation in the ground-state of $^{29}$F within a three-body model

abstract: Background$\colon$ The $^{29}$F system is located at the lower-N boundary of the "island of inversion" and is an exotic, weakly bound system. Little is known about this system beyond its two-neutron separation energy ($S_{2n}$) with large uncertainties. A similar situation is found for the low-lying spectrum of its unbound binary subsystem $^{28}$F. Purpose$\colon$ To investigate the configuration mixing, matter radius and neutron-neutron correlations in the ground state of $^{29}$F within a three-body model, exploring the possibility of $^{29}$F to be a two-neutron halo nucleus. Method$\colon$ The $^{29}$F ground-state wave function is built within the hyperspherical formalism by using an analytical transformed harmonic oscillator basis. The Gogny-Pires-Tourreil (GPT) nn interaction with central, spin-orbit and tensor terms is employed in the present calculations, together with different core$+n$ potentials constrained by the available experimental information on $^{28}$F. Results$\colon$ The $^{29}$F ground-state configuration mixing and its matter radius are computed for different choices of the $^{28}$F structure and $S_{2n}$ value. The admixture of d-waves with pf components are found to play an important role, favoring the dominance of dineutron configurations in the wave function. Our computed radii show a mild sensitivity to the $^{27}$F$+n$ potential and $S_{2n}$ values. The relative increase of the matter radius with respect to the $^{27}$F core lies in the range 0.1-0.4 fm depending upon these choices. Conclusions$\colon$ Our three-body results for $^{29}$F indicate the presence of a moderate halo structure in its ground state, which is enhanced by larger intruder components. This finding calls for an experimental confirmation.


title: Li abundances for solar twins in the open cluster M67

abstract: We determine lithium (Li) abundances for solar twins in the M67 open cluster to add valuable information about the correlation between Li depletion and stellar age and, then, better understand stellar structure and evolution. We use high resolution and good signal-to-noise ratio spectra to characterize Li depletion in three solar twins from M67, using spectral synthesis in the region of the asymmetric 6707.75 \AA \, Li I feature. The mean Li abundance value of A(Li)$=1.6\pm0.2$ dex for our sample of M67 solar twins (our three stars plus a fourth solar twin from a previous analysis in the literature) presents Li abundance expected for its age. Also, the scatter estimated from the standard deviation of the Li abundances in this work is similar to the typical scatter found in a sample of field solar twins presented in the literature.


title: Uncertainty Quantification for Materials Properties in Density Functional Theory with k-Point Density

abstract: Many computational databases emerged over the last five years that report material properties calculated with density functional theory. The properties in these databases are commonly calculated to a precision that is set by choice of the basis set and the k-point density for the Brillouin zone integration. We determine how the precision of properties obtained from the Birch equation of state for 29 transition metals and aluminum in the three common structures -- fcc, bcc, and hcp -- correlate with the k-point density and the precision of the energy. We show that the precision of the equilibrium volume, bulk modulus, and the pressure derivative of the bulk modulus correlate comparably well with the k-point density and the precision of the energy, following an approximate power law. We recommend the k-point density as the convergence parameter because it is computationally efficient, easy to use as a direct input parameter, and correlates with property precision at least as well as the energy precision. We predict that common k-point density choices in high throughput DFT databases result in precision for the volume of 0.1%, the bulk modulus of 1%, and the pressure derivative of 10%.


title: Impact of network assortativity on epidemic and vaccination behaviour

abstract: The resurgence of measles is largely attributed to the decline in vaccine adoption and the increase in mobility. Although the vaccine for measles is readily available and highly successful, its current adoption is not adequate to prevent epidemics. Vaccine adoption is directly affected by individual vaccination decisions, and has a complex interplay with the spatial spread of disease shaped by an underlying mobility (travelling) network. In this paper, we model the travelling connectivity as a scale-free network, and investigate dependencies between the network's assortativity and the resultant epidemic and vaccination dynamics. In doing so we extend an SIR-network model with game-theoretic components, capturing the imitation dynamics under a voluntary vaccination scheme. Our results show a correlation between the epidemic dynamics and the network's assortativity, highlighting that networks with high assortativity tend to suppress epidemics under certain conditions. In highly assortative networks, the suppression is sustained producing an early convergence to equilibrium. In highly disassortative networks, however, the suppression effect diminishes over time due to scattering of non-vaccinating nodes, and frequent switching between the predominantly vaccinating and non-vaccinating phases of the dynamics.


title: Starbursting [O III] emitters and quiescent [C II] emitters in the reionization era

abstract: Recent observations have successfully detected [O III] $88.3\,{\rm \mu m}$ and [C II] $157.6\,{\rm \mu m}$ lines from galaxies in the early Universe with the Atacama Large Millimeter Array (ALMA). Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present relations between the metal line emission and galaxy evolution at $z=6-15$. We find that galaxies during their starburst phases have high [O III] luminosity of $\sim 10^{42}~\rm erg~s^{-1}$. Once supernova feedback quenches star formation, [O III] luminosities rapidly decrease and continue to be zero for $\sim 100\,{\rm Myr}$. The slope of the relation between $\log{(\rm SFR/M_{\odot}~yr^{-1})}$ and $\log{(L_{\rm [O_{III}]}/{\rm L_{\odot}})}$ at $z=6-9$ is 1.03, and 1.43 for $\log{(L_{\rm [C_{II}]}/{\rm L_{\odot}})}$. As gas metallicity increases from sub-solar to solar metallicity by metal enrichment from star formation and feedback, the line luminosity ratio $L_{\rm [O_{III}]} / L_{\rm [C_{II}]}$ decreases from $\sim 10$ to $\sim 1$ because the O/C abundance ratio decreases due to carbon-rich winds from AGB stars and the mass ratio of {\sc H\,ii} to {\sc H\,i} regions decreases due to rapid recombination. Therefore, we suggest that the combination of [O III] and [C II] lines is a good probe to investigate the relative distribution of ionized and neutral gas in high-$z$ galaxies. In addition, we show that deep [C II] observations with a sensitivity of $\sim 10^{-2}~{\rm mJy~arcsec^{-2}}$ can probe the extended neutral gas disks of high-$z$ galaxies.


title: Switching dynamics of single and coupled VO2-based oscillators as elements of neural networks

abstract: In the present paper, we report on the switching dynamics of both single and coupled VO2-based oscillators, with resistive and capacitive coupling, and explore the capability of their application in oscillatory neural networks. Based on these results, we further select an adequate SPICE model to describe the modes of operation of coupled oscillator circuits. Physical mechanisms influencing the time of forward and reverse electrical switching, that determine the applicability limits of the proposed model, are identified. For the resistive coupling, it is shown that synchronization takes place at a certain value of the coupling resistance, though it is unstable and a synchronization failure occurs periodically. For the capacitive coupling, two synchronization modes, with weak and strong coupling, are found. The transition between these modes is accompanied by chaotic oscillations. A decrease in the width of the spectrum harmonics in the weak-coupling mode, and its increase in the strong-coupling one, is detected. The dependences of frequencies and phase differences of the coupled oscillatory circuits on the coupling capacitance are found. Examples of operation of coupled VO2 oscillators as a central pattern generator are demonstrated.


title: Integrality of $v$-adic multiple zeta values

abstract: In this article, we prove the integrality of $v$-adic multiple zeta values (MZVs). For any index $\mathfrak{s}\in\mathbb{N}^r$ and finite place $v\in A:=\mathbb{F}_q[\theta]$, Chang and Mishiba introduced the notion of the $v$-adic MZVs $\zeta_A(\mathfrak{s})_v$, which is a function field analogue of Furusho's $p$-adic MZVs. By estimating the $v$-adic valuation of $\zeta_A(\mathfrak{s})_v$, we show that $\zeta_A(\mathfrak{s})_v$ is a $v$-adic integer for almost all $v$. This result can be viewed as a function field analogue of the integrality of $p$-adic MZVs, which was proved by Akagi-Hirose-Yasuda and Chatzistamatiou.


title: Rigidity theorem by the minimal point of the Bergman kernel

abstract: We use the Suita conjecture (now a theorem) to prove that for any domain $\Omega \subset \mathbb{C}$ its Bergman kernel $K(\cdot, \cdot)$ satisfies $K(z_0, z_0) = \hbox{Volume}(\Omega)^{-1}$ for some $z_0 \in \Omega$ if and only if $\Omega$ is either a disk minus a (possibly empty) closed polar set or $\mathbb{C}$ minus a (possibly empty) closed polar set. When $\Omega$ is bounded with $C^{\infty}$-boundary, we provide a simple proof of this using the zero set of the Szeg\"o kernel. Finally, we show that this theorem fails to hold in $\mathbb{C}^n$ for $n > 1$ by constructing a bounded complete Reinhardt domain (with algebraic boundary) which is strongly convex and not biholomorphic to the unit ball $\mathbb{B}^n \subset \mathbb{C}^n$.


title: A novel experimental setup for an oblique impact onto an inclined granular layer

abstract: We develop an original apparatus of the granular impact experiment by which the incident angle of the solid projectile and inclination angle of the target granular layer can be systematically varied. Whereas most of the natural cratering events occur on inclined surfaces with various incident angles, there have not been any experiments on oblique impacts on an inclined target surface. To perform systematic impact experiments, a novel experimental apparatus has to be developed. Therefore, we build an apparatus for impact experiments where both the incident angle and the inclination angle can be independently varied. The projectile-injection unit accelerates a plastic ball (6~mm in diameter) up to $v_i\simeq 100$~m~s$^{-1}$ impact velocity. The barrel of the injection unit is made with a three-dimensional printer. The impact dynamics is captured by high-speed cameras to directly measure the impact velocity and incident angle. The rebound dynamics of the projectile (restitution coefficient and rebound angle) is also measured. The final crater shapes are measured using a line-laser profiler mounted on the electric stages. By scanning the surface using this system, a three-dimensional crater shape (height map) can be constructed. From the measured result, we can define and measure the characteristic quantities of the crater. The analyzed result on the restitution dynamics is presented as an example of systematic experiments using the developed system.


title: High Performance I/O For Large Scale Deep Learning

abstract: Training deep learning (DL) models on petascale datasets is essential for achieving competitive and state-of-the-art performance in applications such as speech, video analytics, and object recognition. However, existing distributed filesystems were not developed for the access patterns and usability requirements of DL jobs. In this paper, we describe AIStore, a highly scalable, easy-to-deploy storage system, and WebDataset, a standards-based storage format and library that permits efficient access to very large datasets. We compare system performance experimentally using image classification workloads and storing training data on a variety of backends, including local SSDs, single-node NFS, and two identical bare-metal clusters: HDFS and AIStore.


title: Tunable non-reciprocal quantum transport through a dissipative Aharonov-Bohm ring in ultracold atoms

abstract: We report the experimental observation of tunable, non-reciprocal quantum transport of a Bose-Einstein condensate in a momentum lattice. By implementing a dissipative Aharonov-Bohm (AB) ring in momentum space and sending atoms through it, we demonstrate a directional atom flow by measuring the momentum distribution of the condensate at different times. While the dissipative AB ring is characterized by the synthetic magnetic flux through the ring and the laser-induced loss on it, both the propagation direction and transport rate of the atom flow sensitively depend on these highly tunable parameters. We demonstrate that the non-reciprocity originates from the interplay of the synthetic magnetic flux and the laser-induced loss, which simultaneously breaks the inversion and the time-reversal symmetries. Our results open up the avenue for investigating non-reciprocal dynamics in cold atoms, and highlight the dissipative AB ring as a flexible building element for applications in quantum simulation and quantum information.


title: A simple microstructural explanation of the concavity of price impact

abstract: This article provides a simple explanation of the asymptotic concavity of the price impact of a meta-order via the microstructural properties of the market. This explanation is made more precise by a model in which the local relationship between the order flow and the fundamental price (i.e. the local price impact) is linear, with a constant slope, which makes the model dynamically consistent. Nevertheless, the expected impact on midprice from a large sequence of co-directional trades is nonlinear and asymptotically concave. The main practical conclusion of the proposed explanation is that, throughout a meta-order, the volumes at the best bid and ask prices change (on average) in favor of the executor. This conclusion, in turn, relies on two more concrete predictions, one of which can be tested, at least for large-tick stocks, using publicly available market data.


title: Vamsa: Automated Provenance Tracking in Data Science Scripts

abstract: There has recently been a lot of ongoing research in the areas of fairness, bias and explainability of machine learning (ML) models due to the self-evident or regulatory requirements of various ML applications. We make the following observation: All of these approaches require a robust understanding of the relationship between ML models and the data used to train them. In this work, we introduce the ML provenance tracking problem: the fundamental idea is to automatically track which columns in a dataset have been used to derive the features/labels of an ML model. We discuss the challenges in capturing such information in the context of Python, the most common language used by data scientists. We then present Vamsa, a modular system that extracts provenance from Python scripts without requiring any changes to the users' code. Using 26K real data science scripts, we verify the effectiveness of Vamsa in terms of coverage, and performance. We also evaluate Vamsa's accuracy on a smaller subset of manually labeled data. Our analysis shows that Vamsa's precision and recall range from 90.4% to 99.1% and its latency is in the order of milliseconds for average size scripts. Drawing from our experience in deploying ML models in production, we also present an example in which Vamsa helps automatically identify models that are affected by data corruption issues.


title: A Behavioural Theory of Recursive Algorithms

abstract: "What is an algorithm?" is a fundamental question of computer science. Gurevich's behavioural theory of sequential algorithms (aka the sequential ASM thesis) gives a partial answer by defining (non-deterministic) sequential algorithms axiomatically, without referring to a particular machine model or programming language, and showing that they are captured by (non-deterministic) sequential Abstract State Machines (nd-seq ASMs). Moschovakis pointed out that recursive algorithms such as mergesort are not covered by this theory. In this article we propose an axiomatic definition of the notion of sequential recursive algorithm which extends Gurevich's axioms for sequential algorithms by a Recursion Postulate and allows us to prove that sequential recursive algorithms are captured by recursive Abstract State Machines, an extension of nd-seq ASMs by a CALL rule. Applying this recursive ASM thesis yields a characterization of sequential recursive algorithms as finitely composed concurrent algorithms all of whose concurrent runs are partial-order runs.


title: Text Complexity Classification Based on Linguistic Information: Application to Intelligent Tutoring of ESL

abstract: The goal of this work is to build a classifier that can identify text complexity within the context of teaching reading to English as a Second Language (ESL) learners. To present language learners with texts that are suitable to their level of English, a set of features that can describe the phonological, morphological, lexical, syntactic, discursive, and psychological complexity of a given text were identified. Using a corpus of 6171 texts, which had already been classified into three different levels of difficulty by ESL experts, different experiments were conducted with five machine learning algorithms. The results showed that the adopted linguistic features provide a good overall classification performance (F-Score = 0.97). A scalability evaluation was conducted to test if such a classifier could be used within real applications, where it can be, for example, plugged into a search engine or a web-scraping module. In this evaluation, the texts in the test set are not only different from those from the training set but also of different types (ESL texts vs. children reading texts). Although the overall performance of the classifier decreased significantly (F-Score = 0.65), the confusion matrix shows that most of the classification errors are between the classes two and three (the middle-level classes) and that the system has a robust performance in categorizing texts of class one and four. This behavior can be explained by the difference in classification criteria between the two corpora. Hence, the observed results confirm the usability of such a classifier within a real-world application.


title: Full-sky searches for anisotropies in UHECR arrival directions with the Pierre Auger Observatory and the Telescope Array

abstract: The arrival directions of ultra-high-energy cosmic rays appear to be approximately isotropically distributed over the whole sky, but the last-generation UHECR detector arrays, the Pierre Auger Observatory (Auger) and the Telescope Array (TA), have detected thousands of events, allowing us to study small deviations from isotropy. So far, Auger has detected a large-scale anisotropy consistent with a $\sim 6.5\%$ dipole moment in the distribution of cosmic rays with $E > 8$~EeV, and both collaborations have reported indications for smaller-scale anisotropies at higher energies. On the other hand, neither array has full-sky coverage, the Auger field of view being limited to declinations $\delta < +45^\circ$ and the TA one to $\delta > -16^\circ$. Searches for anisotropies with full-sky coverage thus require combining data from both arrays. A working group with members from both collaborations has been established for this task. Since even a minor systematic error in the energy determinations at either of the arrays could result in a sizeable spurious anisotropy along the north--south axis, we devised a method to cross-calibrate the energy scales of the two experiments with respect to each other by using events in a declination band within the intersection of their fields of view. In this contribution, we report on both updates on blind searches for anisotropies and the first full-sky studies of flux models based on possible local extragalactic sources.


title: Effective Scaling of Blockchain Beyond Consensus Innovations and Moore's Law

abstract: As an emerging technology, blockchain has achieved great success in numerous application scenarios, from intelligent healthcare to smart cities. However, a long-standing bottleneck hindering its further development is the massive resource consumption attributed to the distributed storage and computation methods. This makes blockchain suffer from insufficient performance and poor scalability. Here, we analyze the recent blockchain techniques and demonstrate that the potential of widely-adopted consensus-based scaling is seriously limited, especially in the current era when Moore's law-based hardware scaling is about to end. We achieve this by developing an open-source benchmarking tool, called Prism, for investigating the key factors causing low resource efficiency and then discuss various topology and hardware innovations which could help to scale up blockchain. To the best of our knowledge, this is the first in-depth study that explores the next-generation scaling strategies by conducting large-scale and comprehensive benchmarking.


title: Reinforcement Learning via Fenchel-Rockafellar Duality

abstract: We review basic concepts of convex duality, focusing on the very general and supremely useful Fenchel-Rockafellar duality. We summarize how this duality may be applied to a variety of reinforcement learning (RL) settings, including policy evaluation or optimization, online or offline learning, and discounted or undiscounted rewards. The derivations yield a number of intriguing results, including the ability to perform policy evaluation and on-policy policy gradient with behavior-agnostic offline data and methods to learn a policy via max-likelihood optimization. Although many of these results have appeared previously in various forms, we provide a unified treatment and perspective on these results, which we hope will enable researchers to better use and apply the tools of convex duality to make further progress in RL.


title: Invariance of polymer partition functions under the geometric RSK correspondence

abstract: We prove that the values of discrete directed polymer partition functions involving multiple non-intersecting paths remain invariant under replacing the background weights by their images under the geometric RSK correspondence. This result is inspired by a recent and remarkable identity proved by Dauvergne, Orthmann and Vir\'ag which is recovered as the zero-temperature, semi-discrete limit of our main result.


title: Closed loop application of electroadhesion for increased precision in texture rendering

abstract: Tactile displays based on friction modulation offer wide-bandwidth forces rendered directly on the fingertip. However, due to a number of touch conditions (e.g., normal force, skin hydration) that result in variations in the friction force and the strength of modulation effect, the precision of the force rendering remains limited. In this paper we demonstrate a closed-loop electroadhesion method for precise playback of friction force profiles on a human finger and we apply this method to the tactile rendering of several textiles encountered in everyday life.


title: HybridPose: 6D Object Pose Estimation under Hybrid Representations

abstract: We introduce HybridPose, a novel 6D object pose estimation approach. HybridPose utilizes a hybrid intermediate representation to express different geometric information in the input image, including keypoints, edge vectors, and symmetry correspondences. Compared to a unitary representation, our hybrid representation allows pose regression to exploit more and diverse features when one type of predicted representation is inaccurate (e.g., because of occlusion). Different intermediate representations used by HybridPose can all be predicted by the same simple neural network, and outliers in predicted intermediate representations are filtered by a robust regression module. Compared to state-of-the-art pose estimation approaches, HybridPose is comparable in running time and accuracy. For example, on Occlusion Linemod dataset, our method achieves a prediction speed of 30 fps with a mean ADD(-S) accuracy of 47.5%, representing a state-of-the-art performance. The implementation of HybridPose is available at https://github.com/chensong1995/HybridPose.


title: MW-GAN: Multi-Warping GAN for Caricature Generation with Multi-Style Geometric Exaggeration

abstract: Given an input face photo, the goal of caricature generation is to produce stylized, exaggerated caricatures that share the same identity as the photo. It requires simultaneous style transfer and shape exaggeration with rich diversity, and meanwhile preserving the identity of the input. To address this challenging problem, we propose a novel framework called Multi-Warping GAN (MW-GAN), including a style network and a geometric network that are designed to conduct style transfer and geometric exaggeration respectively. We bridge the gap between the style and landmarks of an image with corresponding latent code spaces by a dual way design, so as to generate caricatures with arbitrary styles and geometric exaggeration, which can be specified either through random sampling of latent code or from a given caricature sample. Besides, we apply identity preserving loss to both image space and landmark space, leading to a great improvement in quality of generated caricatures. Experiments show that caricatures generated by MW-GAN have better quality than existing methods.


title: Attention over Parameters for Dialogue Systems

abstract: Dialogue systems require a great deal of different but complementary expertise to assist, inform, and entertain humans. For example, different domains (e.g., restaurant reservation, train ticket booking) of goal-oriented dialogue systems can be viewed as different skills, and so does ordinary chatting abilities of chit-chat dialogue systems. In this paper, we propose to learn a dialogue system that independently parameterizes different dialogue skills, and learns to select and combine each of them through Attention over Parameters (AoP). The experimental results show that this approach achieves competitive performance on a combined dataset of MultiWOZ, In-Car Assistant, and Persona-Chat. Finally, we demonstrate that each dialogue skill is effectively learned and can be combined with other skills to produce selective responses.


title: Spatial Applications of Topological Data Analysis: Cities, Snowflakes, Random Structures, and Spiders Spinning Under the Influence

abstract: Spatial networks are ubiquitous in social, geographical, physical, and biological applications. To understand the large-scale structure of networks, it is important to develop methods that allow one to directly probe the effects of space on structure and dynamics. Historically, algebraic topology has provided one framework for rigorously and quantitatively describing the global structure of a space, and recent advances in topological data analysis (TDA) have given scholars a new lens for analyzing network data. In this paper, we study a variety of spatial networks -- including both synthetic and natural ones -- using novel topological methods that we recently developed for analyzing spatial networks. We demonstrate that our methods are able to capture meaningful quantities, with specifics that depend on context, in spatial networks and thereby provide useful insights into the structure of those networks, including a novel approach for characterizing them based on their topological structures. We illustrate these ideas with examples of synthetic networks and dynamics on them, street networks in cities, snowflakes, and webs spun by spiders under the influence of various psychotropic substances.

categories: cs.SI cond-mat.soft cs.CG math.AT physics.soc-ph


title: Behavioural Theory of Reflective Algorithms I: Reflective Sequential Algorithms

abstract: We develop a behavioural theory of reflective sequential algorithms (RSAs), i.e. sequential algorithms that can modify their own behaviour. The theory comprises a set of language-independent postulates defining the class of RSAs, an abstract machine model, and the proof that all RSAs are captured by this machine model. As in Gurevich's behavioural theory for sequential algorithms RSAs are sequential-time, bounded parallel algorithms, where the bound depends on the algorithm only and not on the input. Different from the class of sequential algorithms every state of an RSA includes a representation of the algorithm in that state, thus enabling linguistic reflection. Bounded exploration is preserved using terms as values. The model of reflective sequential abstract state machines (rsASMs) extends sequential ASMs using extended states that include an updatable representation of the main ASM rule to be executed by the machine in that state. Updates to the representation of ASM signatures and rules are realised by means of a sophisticated tree algebra.


title: Making spaces wild (simply-connected case)

abstract: We attach copies of the circle to points of a countable dense subset $D$ of a separable metric space $X$ and construct an earring space $E(X,D)$. We show that the fundamental group of $E(X,D)$ is isomorphic to a subgroup of the Hawaiian earring group, if the space $X$ is simply-connected and locally simply-connected. In addition if the space $X$ is locally path-connected, the space $X$ can be recovered from the fundamental group of $E(X,D)$.


title: Poset subdivisions and the mixed cd-index

abstract: The cd-index is an invariant of Eulerian posets expressed as a polynomial in noncommuting variables c and d. It determines another invariant, the h-polynomial. In this paper, we study the relative setting, that of subdivisions of posets. We introduce the mixed cd-index, an invariant of strong formal subdivisions of posets, which determines the mixed h-polynomial introduced by the second author with Stapledon. The mixed cd-index is a polynomial in noncommuting variables c',d',c,d, and e and is defined in terms of the local cd-index of Karu. Here, use is made of the decomposition theorem for the cd-index. We extend the proof of the decomposition theorem, originally due to Ehrenborg-Karu, to the class of strong formal subdivisions. We also compute the mixed cd-index in a number of examples.


title: On the error in the two-term Weyl formula for the Dirichlet Laplacian

abstract: We study the optimality of the remainder term in the two-term Weyl law for the Dirichlet Laplacian within the class of Lipschitz regular subsets of $\mathbb{R}^d$. In particular, for the short-time asymptotics of the trace of the heat kernel we prove that the error term cannot be made quantitatively better than little-$o$ of the second term.


title: Null controllability and inverse source problem for stochastic Grushin equation with boundary degeneracy and singularity

abstract: In this paper, we consider a null controllability and an inverse source problem for stochastic Grushin equation with boundary degeneracy and singularity. We construct two special weight functions to establish two Carleman estimates for the whole stochastic Grushin operator with singular potential by a weighted identity method. One is for the backward stochastic Grushin equation with singular weight function. We then apply it to prove the null controllability for stochastic Grushin equation for any $T$ and any degeneracy $\gamma>0$, when our control domain touches the degeneracy line $\{x=0\}$. In order to study the inverse source problem of determining two kinds of sources simultaneously, we prove the other Carleman estimate, which is for the forward stochastic Grushin equation with regular weight function. Based on this Carleman estimate, we obtain the uniqueness of the inverse source problem.


title: Phase Transitions for the Information Bottleneck in Representation Learning

abstract: In the Information Bottleneck (IB), when tuning the relative strength between compression and prediction terms, how do the two terms behave, and what's their relationship with the dataset and the learned representation? In this paper, we set out to answer these questions by studying multiple phase transitions in the IB objective: $\text{IB}_\beta[p(z|x)] = I(X; Z) - \beta I(Y; Z)$ defined on the encoding distribution p(z|x) for input $X$, target $Y$ and representation $Z$, where sudden jumps of $dI(Y; Z)/d \beta$ and prediction accuracy are observed with increasing $\beta$. We introduce a definition for IB phase transitions as a qualitative change of the IB loss landscape, and show that the transitions correspond to the onset of learning new classes. Using second-order calculus of variations, we derive a formula that provides a practical condition for IB phase transitions, and draw its connection with the Fisher information matrix for parameterized models. We provide two perspectives to understand the formula, revealing that each IB phase transition is finding a component of maximum (nonlinear) correlation between $X$ and $Y$ orthogonal to the learned representation, in close analogy with canonical-correlation analysis (CCA) in linear settings. Based on the theory, we present an algorithm for discovering phase transition points. Finally, we verify that our theory and algorithm accurately predict phase transitions in categorical datasets, predict the onset of learning new classes and class difficulty in MNIST, and predict prominent phase transitions in CIFAR10.


title: In Leon's company, it seemed that anything might be possible

abstract: Memorial sessions celebrated Leon Lederman, Helen Edwards, and Burton Richter at the April 2019 Meeting of the American Physical Society in Denver. In the session entitled Honoring Leon Lederman, Sally Dawson gave an overview of Leon's scientific career, Marge Bardeen reviewed his work in science education, and I spoke of his years as Fermilab Director. This essay is drawn from my lecture.


title: Convexification for an Inverse Parabolic Problem

abstract: A convexification-based numerical method for a Coefficient Inverse Problem for a parabolic PDE is presented. The key element of this method is the presence of the so-called Carleman Weight Function in the numerical scheme. Convergence analysis ensures the global convergence of this method, as opposed to the local convergence of the conventional least squares minimization techniques. Numerical results demonstrate a good performance.


title: Completely determined Borel sets and measurability

abstract: We consider the reverse math strength of the statement $\mathsf{C\text-DM}$:"Every completely determined Borel set is measurable." Over $\mathsf{WWKL}_0$, we obtain the following results analogous to the previously studied category case. $\mathsf{C\text-DM}$ lies strictly between $\mathsf{ATR}_0$ and $\mathsf{L}_{\omega_1,\omega}\text-\mathsf{CA}$. Whenever $M\subseteq 2^\omega$ is the second-order part of an $\omega$-model of $\mathsf{C\text-DM}$, then for every $Z \in M$, there is a $R \in M$ such that $R$ is $\Delta^1_1$-random relative to $Z$. On the other hand, without $\mathsf{WWKL}_0$, all sets have measure zero (as measured according to $\mathsf{C\text-DM}$), and it follows vacuously that $\neg \mathsf{WWKL}_0$ implies $\mathsf{C\text-DM}$ over $\mathsf{RCA}_0$.


title: A unique continuation property for a class of parabolic differential inequalities in a bounded domain

abstract: This article is concerned with the unique continuation property of a forward differential inequality abstracted from parabolic equations proposed on a convex domain $\Omega$ prescribed with some regularity and growth conditions. Our result shows that the value of the solutions can be determined uniquely by its value on an arbitrary open subset $\omega$ in $\Omega$ at any given positive time $T$. We also derive the quantitative nature of this unique continuation, that is, the estimate of a $L^2({\Omega})$ norm of the initial data on $\Omega$, which is majorized by that of solution on the bounded open subset $\omega$ at terminal moment $t = T$.


title: Boxicity of Circulant Graph $G_k^d$

abstract: The boxicity of a graph $G$, denoted by $box(G)$, is the least positive integer $\ell$ such that $G$ can be isomorphic to the intersection graph of a family of boxes in Euclidean $\ell$-space, where box in an Euclidean $\ell$-space is the Cartesian product of $\ell$ closed intervals on the real line. Let $k$ and $d$ be two positive integers with $k\geq 2d$. The circulant graph $G_k^d$ is the graph with vertices set $V(G_k^d)=\{a_0, a_1,\ldots, a_{k-1}\}$ and edge set $E(G_k^d)=\{a_i a_j | \ d\leq |i-j|\leq k-d\}$. Denote $\chi(G)$ the chromatic number of a graph $G$. In \cite{Aki} Akira Kamibeppu proved that $box(G_k^d)\leq \chi(G_k^d)$ for some class of circulant graph $G_k^d$ and raised the question that the same result holds for all circulant graph. In this short note, we prove that $box(G_k^d)\leq \chi(G_k^d)$, for all $k$ and $d$ with $k\geq 2d$. This include all circulant graph $G_k^d$. Our proof is very simple and short. This answer the above question.


title: Stochastic Geometry Analysis of Sojourn Time in Multi-Tier Cellular Networks

abstract: Impact of mobility will be increasingly important in future generation wireless services and the related challenges will need to be addressed. Sojourn time, the time duration that a mobile user stays within a cell, is a mobility-aware parameter that can significantly impact the performance of mobile users and it can also be exploited to improve resource allocation and mobility management methods in the network. In this paper, we derive the distribution and mean of the sojourn time in multi-tier cellular networks, where spatial distribution of base stations (BSs) in each tier follows an independent homogeneous Poisson point process (PPP). To obtain the sojourn time distribution in multi-tier cellular networks with maximum biased averaged received power association, we derive the linear contact distribution function and chord length distribution of each tier. We also study the relation between mean sojourn time and other mobility-related performance metrics. We show that the mean sojourn time is inversely proportional to the handoff rate, and the complementary cumulative distribution function (CCDF) of sojourn time is bounded from above by the complement of the handoff probability. Moreover, we study the impact of user velocity and network parameters on the sojourn time.


title: Discovering Nonlinear Relations with Minimum Predictive Information Regularization

abstract: Identifying the underlying directional relations from observational time series with nonlinear interactions and complex relational structures is key to a wide range of applications, yet remains a hard problem. In this work, we introduce a novel minimum predictive information regularization method to infer directional relations from time series, allowing deep learning models to discover nonlinear relations. Our method substantially outperforms other methods for learning nonlinear relations in synthetic datasets, and discovers the directional relations in a video game environment and a heart-rate vs. breath-rate dataset.


title: From Open Set to Closed Set: Supervised Spatial Divide-and-Conquer for Object Counting

abstract: Visual counting, a task that aims to estimate the number of objects from an image/video, is an open-set problem by nature, i.e., the number of population can vary in [0, inf) in theory. However, collected data and labeled instances are limited in reality, which means that only a small closed set is observed. Existing methods typically model this task in a regression manner, while they are prone to suffer from an unseen scene with counts out of the scope of the closed set. In fact, counting has an interesting and exclusive property---spatially decomposable. A dense region can always be divided until sub-region counts are within the previously observed closed set. We therefore introduce the idea of spatial divide-and-conquer (S-DC) that transforms open-set counting into a closed-set problem. This idea is implemented by a novel Supervised Spatial Divide-and-Conquer Network (SS-DCNet). Thus, SS-DCNet can only learn from a closed set but generalize well to open-set scenarios via S-DC. SS-DCNet is also efficient. To avoid repeatedly computing sub-region convolutional features, S-DC is executed on the feature map instead of on the input image. We provide theoretical analyses as well as a controlled experiment on toy data, demonstrating why closed-set modeling makes sense. Extensive experiments show that SS-DCNet achieves the state-of-the-art performance. Code and models are available at: https://tinyurl.com/SS-DCNet.


title: IRS-based Wireless Jamming Attacks: When Jammers can Attack without Power

abstract: This paper proposes to use Intelligent Reflecting Surface (IRS) as a green jammer to attack a legitimate communication without using any internal energy to generate jamming signals. In particular, the IRS is used to intelligently reflect the signals from the legitimate transmitter to the legitimate receiver (LR) to guarantee that the received signals from direct and reflecting links can be added destructively, which thus diminishes the Signal-to-Interference-plus-Noise Ratio (SINR) at the LR. To minimize the received signal power at the LR, we consider the joint optimization of magnitudes of reflection coefficients and discrete phase shifts at the IRS. Based on the block coordinate descent, semidefinite relaxation, and Gaussian randomization techniques, the solution can be obtained efficiently. Through simulation results, we show that by using the IRS-based jammer, we can reduce the signal power received at the LR by up to 99\%. Interestingly, the performance of the proposed IRS-based jammer is even better than that of the conventional active jamming attacks in some scenarios.


title: Indoor Localization System of ROS mobile robot based on Visible Light Communication

abstract: In this paper, an indoor robot localization system based on Robot Operating System (ROS) and visible light communication (VLC) is presented. On the basis of our previous work, we innovatively designed a VLC localization and navigation package based on Robot Operating System (ROS), which contains the LED-ID detection and recognition method, the video target tracking algorithm and the double-lamp positioning algorithm. This package exploited the principle of double-lamp positioning and the loose coupling characteristics of the ROS system, which is implemented by loosely coupled ROS nodes. Consequently, this paper combines ROS and VLC, aiming at promoting the application of VLC positioning in mature robotic systems. Moreover, it pushed forward the development of localization application based on VLC technology and lays a foundation for transplanting to other ROS robot platforms. Experimental results show that the proposed system can provide indoor localization within 1 cm and possesses a good real-time performance which takes only 0.4 seconds for one-time positioning. And if a high-performance laptop is used, the single positioning time can be reduced to 0.08 seconds. Therefore, this study confirms the practical application and the superior performance of VLC technology in ROS robot, showing the great potential of VLC localization. T he video demo of the proposed robot positioning system based on VLC can be seen in *


title: Quantum Advantage for Shared Randomness Generation

abstract: Sharing correlated random variables is a resource for a number of information theoretic tasks such as privacy amplification, simultaneous message passing, secret sharing and many more. In this article, we show that to establish such a resource called shared randomness, quantum systems provide an advantage over their classical counterpart. Precisely, we show that appropriate albeit fixed measurements on a shared two-qubit state can generate correlations which cannot be obtained from any possible state on two classical bits. In a resource theoretic set-up, this feature of quantum systems can be interpreted as an advantage in winning a two players co-operative game, which we call the `non-monopolize social subsidy' game. It turns out that the quantum states leading to the desired advantage must possess non-classicality in the form of quantum discord. On the other hand, while distributing such sources of shared randomness between two parties via noisy channels, quantum channels with zero capacity as well as with classical capacity strictly less than unity perform more efficiently than the perfect classical channel. Protocols presented here are noise-robust and hence should be realizable with state-of-the-art quantum devices.


title: Statistical Inference for High-Dimensional Matrix-Variate Factor Model

abstract: This paper considers the estimation and inference of the low-rank components in high-dimensional matrix-variate factor models, where each dimension of the matrix-variates ($p \times q$) is comparable to or greater than the number of observations ($T$). We propose an estimation method called $\alpha$-PCA that preserves the matrix structure and aggregates mean and contemporary covariance through a hyper-parameter $\alpha$. We develop an inferential theory, establishing consistency, the rate of convergence, and the limiting distributions, under general conditions that allow for correlations across time, rows, or columns of the noise. We show both theoretical and empirical methods of choosing the best $\alpha$, depending on the use-case criteria. Simulation results demonstrate the adequacy of the asymptotic results in approximating the finite sample properties. The $\alpha$-PCA compares favorably with the existing ones. Finally, we illustrate its applications with a real numeric data set and two real image data sets. In all applications, the proposed estimation procedure outperforms previous methods in the power of variance explanation using out-of-sample 10-fold cross-validation.


title: IMLI: An Incremental Framework for MaxSAT-Based Learning of Interpretable Classification Rules

abstract: The wide adoption of machine learning in the critical domains such as medical diagnosis, law, education had propelled the need for interpretable techniques due to the need for end users to understand the reasoning behind decisions due to learning systems. The computational intractability of interpretable learning led practitioners to design heuristic techniques, which fail to provide sound handles to tradeoff accuracy and interpretability. Motivated by the success of MaxSAT solvers over the past decade, recently MaxSAT-based approach, called MLIC, was proposed that seeks to reduce the problem of learning interpretable rules expressed in Conjunctive Normal Form (CNF) to a MaxSAT query. While MLIC was shown to achieve accuracy similar to that of other state of the art black-box classifiers while generating small interpretable CNF formulas, the runtime performance of MLIC is significantly lagging and renders approach unusable in practice. In this context, authors raised the question: Is it possible to achieve the best of both worlds, i.e., a sound framework for interpretable learning that can take advantage of MaxSAT solvers while scaling to real-world instances? In this paper, we take a step towards answering the above question in affirmation. We propose IMLI: an incremental approach to MaxSAT based framework that achieves scalable runtime performance via partition-based training methodology. Extensive experiments on benchmarks arising from UCI repository demonstrate that IMLI achieves up to three orders of magnitude runtime improvement without loss of accuracy and interpretability.


title: The roles of individual advanced LIGO detectors on revealing the neutron star properties

abstract: In this work we re-analyze the data of GW170817, the first binary neutron star (BNS) merger event, in two ways, including the parameterized equation of state (EoS) method and gravitational wave (GW) parameter estimation analysis. Three sets of data, including the combined LIGO/Virgo (HLV) data, the LIGO-Hanford (H1) data and the LIGO-Livingston (L1) data, have been analyzed, respectively. We confirm the bimodal probability distribution of the tidal deformability of GW170817 with the HLV data and resolve the different contributions of H1 and L1 detectors. Our simulation reveals a tendency of ``overestimating" the tidal parameter with the decreasing signal-to-noise ratio (SNR). Such a tendency is likely caused by the dependence of the result on the prior (which is not well understood and has been assumed to be widely-distributed) in the low SNR case. Though this effect is interesting, the slightly higher SNR of H1 than L1 at frequencies above $\sim 100$ Hz is most-likely not the main reason for the lower tidal parameter region favored by H1. In view of the large fluctuation of the expected tide parameter in the case of low SNR, the different probability distributions favored by the L1 and H1 data of GW170817 are reasonable. We have also explored the dependence of expected difference of SNRs of H1 and L1 detectors on the source locations in the O3 and full sensitivity runs. Such maps are helpful in evaluating the relative powers of individual detectors on measuring the tidal deformabilities for the new double neutron star merger events.


title: Fast and robust multiplane single molecule localization microscopy using deep neural network

abstract: Single molecule localization microscopy is widely used in biological research for measuring the nanostructures of samples smaller than the diffraction limit. This study uses multifocal plane microscopy and addresses the 3D single molecule localization problem, where lateral and axial locations of molecules are estimated. However, when we multifocal plane microscopy is used, the estimation accuracy of 3D localization is easily deteriorated by the small lateral drifts of camera positions. We formulate a 3D molecule localization problem along with the estimation of the lateral drifts as a compressed sensing problem, A deep neural network was applied to accurately and efficiently solve this problem. The proposed method is robust to the lateral drifts and achieves an accuracy of 20 nm laterally and 50 nm axially without an explicit drift correction.


title: Causal Mosaic: Cause-Effect Inference via Nonlinear ICA and Ensemble Method

abstract: We address the problem of distinguishing cause from effect in bivariate setting. Based on recent developments in nonlinear independent component analysis (ICA), we train nonparametrically general nonlinear causal models that allow non-additive noise. Further, we build an ensemble framework, namely Causal Mosaic, which models a causal pair by a mixture of nonlinear models. We compare this method with other recent methods on artificial and real world benchmark datasets, and our method shows state-of-the-art performance.


title: Machine-learning classifiers for logographic name matching in public health applications: approaches for incorporating phonetic, visual, and keystroke similarity in large-scale probabilistic record linkage

abstract: Approximate string-matching methods to account for complex variation in highly discriminatory text fields, such as personal names, can enhance probabilistic record linkage. However, discriminating between matching and non-matching strings is challenging for logographic scripts, where similarities in pronunciation, appearance, or keystroke sequence are not directly encoded in the string data. We leverage a large Chinese administrative dataset with known match status to develop logistic regression and Xgboost classifiers integrating measures of visual, phonetic, and keystroke similarity to enhance identification of potentially-matching name pairs. We evaluate three methods of leveraging name similarity scores in large-scale probabilistic record linkage, which can adapt to varying match prevalence and information in supporting fields: (1) setting a threshold score based on predicted quality of name-matching across all record pairs; (2) setting a threshold score based on predicted discriminatory power of the linkage model; and (3) using empirical score distributions among matches and nonmatches to perform Bayesian adjustment of matching probabilities estimated from exact-agreement linkage. In experiments on holdout data, as well as data simulated with varying name error rates and supporting fields, a logistic regression classifier incorporated via the Bayesian method demonstrated marked improvements over exact-agreement linkage with respect to discriminatory power, match probability estimation, and accuracy, reducing the total number of misclassified record pairs by 21% in test data and up to an average of 93% in simulated datasets. Our results demonstrate the value of incorporating visual, phonetic, and keystroke similarity for logographic name matching, as well as the promise of our Bayesian approach to leverage name-matching within large-scale record linkage.


title: An Efficient Gradient Projection Method for Structural Topology Optimization

abstract: This paper presents an efficient gradient projection-based method for structural topological optimization problems characterized by a nonlinear objective function which is minimized over a feasible region defined by bilateral bounds and a single linear equality constraint. The specialty of the constraints type, as well as heuristic engineering experiences are exploited to improve the scaling scheme, projection, and searching step. In detail, gradient clipping and a modified projection of searching direction under certain condition are utilized to facilitate the efficiency of the proposed method. Besides, an analytical solution is proposed to approximate this projection with negligible computation and memory costs. Furthermore, the calculation of searching steps is largely simplified. Benchmark problems, including the MBB, the force inverter mechanism, and the 3D cantilever beam are used to validate the effectiveness of the method. The proposed method is implemented in MATLAB which is open-sourced for educational usage.


title: Two_Generalizations_for_Quadratic_Residue_Codes_over_Finite_Fields

abstract: It's well known that the quadratic residue code over finite fields is an interesting class of cyclic codes for its higher minimum distance. Let $g$ be a positive integer and $p,p_{1},\ldots, p_{g}$ be distinct odd primes, the present paper generalizes the constructions for the quadratic residue code with length $p$ to be the length $n=p_{1}\cdots p_{g}$, and to be the case $m$-th residue codes with length $p$ over finite fields, where $m\geq 2$ is a positive integer. Furthermore, a criterion for that these codes are self-orthogonal or complementary dual is obtained, and then the corresponding counting formula are given. In particular, the minimum distance of all 24 quaternary quadratic residue codes $[15,8]$ are determined.


title: Reanalysis of Variance Reduced Temporal Difference Learning

abstract: Temporal difference (TD) learning is a popular algorithm for policy evaluation in reinforcement learning, but the vanilla TD can substantially suffer from the inherent optimization variance. A variance reduced TD (VRTD) algorithm was proposed by Korda and La (2015), which applies the variance reduction technique directly to the online TD learning with Markovian samples. In this work, we first point out the technical errors in the analysis of VRTD in Korda and La (2015), and then provide a mathematically solid analysis of the non-asymptotic convergence of VRTD and its variance reduction performance. We show that VRTD is guaranteed to converge to a neighborhood of the fixed-point solution of TD at a linear convergence rate. Furthermore, the variance error (for both i.i.d.\ and Markovian sampling) and the bias error (for Markovian sampling) of VRTD are significantly reduced by the batch size of variance reduction in comparison to those of vanilla TD. As a result, the overall computational complexity of VRTD to attain a given accurate solution outperforms that of TD under Markov sampling and outperforms that of TD under i.i.d.\ sampling for a sufficiently small conditional number.


title: Half-integer quantized anomalous thermal Hall effect in the Kitaev material $\alpha$-RuCl$_3$

abstract: Heat transport mediated by Majorana edge modes in a magnetic insulator leads to a half-integer thermal quantum Hall conductance, which has recently been reported for the two-dimensional honeycomb material $\alpha$-RuCl$_3$. While the conventional electronic Hall effect requires a perpendicular magnetic field, we find that this is not the case in $\alpha$-RuCl$_3$. Strikingly, the thermal Hall plateau appears even for a magnetic field with no out-of-plane components. The field-angular variation of the quantized thermal Hall conductance has the same sign structure of the topological Chern number, which is either $\pm$1, as the Majorana band structure of the pure Kitaev spin liquid. This observation of a half-integer anomalous thermal Hall effect firmly establishes that the Kitaev interaction is primarily responsible and that the non-Abelian topological order associated with fractionalization of the local magnetic moments persists even in the presence of non-Kitaev interactions in $\alpha$-RuCl$_3$.


title: Regularization via Structural Label Smoothing

abstract: Regularization is an effective way to promote the generalization performance of machine learning models. In this paper, we focus on label smoothing, a form of output distribution regularization that prevents overfitting of a neural network by softening the ground-truth labels in the training data in an attempt to penalize overconfident outputs. Existing approaches typically use cross-validation to impose this smoothing, which is uniform across all training data. In this paper, we show that such label smoothing imposes a quantifiable bias in the Bayes error rate of the training data, with regions of the feature space with high overlap and low marginal likelihood having a lower bias and regions of low overlap and high marginal likelihood having a higher bias. These theoretical results motivate a simple objective function for data-dependent smoothing to mitigate the potential negative consequences of the operation while maintaining its desirable properties as a regularizer. We call this approach Structural Label Smoothing (SLS). We implement SLS and empirically validate on synthetic, Higgs, SVHN, CIFAR-10, and CIFAR-100 datasets. The results confirm our theoretical insights and demonstrate the effectiveness of the proposed method in comparison to traditional label smoothing.


title: Transparency in a periodic chain of quantum emitters strongly coupled to a waveguide

abstract: We demonstrate the emergence of transparent behavior in a chain of periodically spaced non-identical quantums emitters coupled to a waveguide, in the special case when the inter-atomic separation is a half-integral multiple of the resonant wavelength, i.e. $kL$ is an integral multiple of $\pi$, with $k$ being the spatial frequency and $L$ the spatial periodicity. When equal but opposite frequency detunings are assigned in pairs to a system of even number of atoms, perfect transmission ensues. When the chain size is odd, a similar assignment leads to the disappearance of collective effects as the odd atom determines the spectral behavior. We also manifest the robustness of these features against dissipative effects and show, how the spectral behavior hinges significantly on the relative detunings between the atoms as compared to the decay rate. A key distinction from the phenomenon of Electromagnetically Induced Transparency (EIT) is that in the waveguide case, the presence of an intrinsic waveguide mediated phase coupling between the atoms strongly affects the transport properties. Furthermore, while reciprocity in single-photon transport does not generally hold due to the phase coupling, we observe an interesting exception for $kL = n\pi$ at which the waveguide demonstrates reciprocal behavior with regard to both the transmission and reflection coefficients.


title: Monte Carlo Tree Search for Generating Interactive Data Analysis Interfaces

abstract: Interactive tools like user interfaces help democratize data access for end-users by hiding underlying programming details and exposing the necessary widget interface to users. Since customized interfaces are costly to build, automated interface generation is desirable. SQL is the dominant way to analyze data and there already exists logs to analyze data. Previous work proposed a syntactic approach to analyze structural changes in SQL query logs and automatically generates a set of widgets to express the changes. However, they do not consider layout usability and the sequential order of queries in the log. We propose to adopt Monte Carlo Tree Search(MCTS) to search for the optimal interface that accounts for hierarchical layout as well as the usability in terms of how easy to express the query log.


title: Hamiltonian formulation of dust cloud collapse

abstract: We consider the gravitational collapse of self-gravitating spherical dust cloud in the Hamiltonian formalism. We address both homogeneous and inhomogeneous cases. Our novel derivation of the Hamiltonian of the system is based on the improved variational principle that was proposed in \cite{KMM}. The present derivation differs from usual treatments due to the presence of an extra boundary term added to the Hilbert action. As expected, the standard equations of motion are retrieved. However, differently from other treatments, the total Hamiltonian obtained with our procedure in the Schwarzschild time-gauge is identical to the total mass of the system as measured from infinity, as it would be expected. Implications for the quantization of the system are suggested.


title: Dynamics of dioecious population with different fitness of genotypes

abstract: In this paper we study dynamical systems generated by an evolution operator of a dioecious population. This evolution operator is a six-parametric, non-linear operator mapping $[0,1]^2$ to itself. We find all fixed points and under some conditions on parameters we give limit points of trajectories constructed by iterations of the evolution operator.


title: Gravitational reduction of the wave function based on Bohmian quantum potential

abstract: In objective gravitational reduction of the wave function of a quantum system, the classical limit of the system is obtained in terms of the objective properties of the system. On the other hand, in Bohmian quantum mechanics the usual criterion for getting classical limit is the vanishing of the quantum potential or the quantum force of the system, which suffers from the lack of an objective description. In this regard, we investigated the usual criterion of getting the classical limit of a free particle in Bohmian quantum mechanics. Then we argued that how it is possible to have an objective gravitational classical limit related to the Bohmian mechanical concepts like quantum potential or quantum force. Also we derived a differential equation related to the wave function reduction. An interesting connection will be made between Bohmian concepts and gravitational concepts.


title: Optimal Transmission of Multi-Quality Tiled 360 VR Video by Exploiting Multicast Opportunities

abstract: In this paper, we would like to investigate fundamental impacts of multicast opportunities on efficient transmission of a 360 VR video to multiple users in the cases with and without transcoding at each user. We establish a novel mathematical model that reflects the impacts of multicast opportunities on the average transmission energy in both cases and the transcoding energy in the case with user transcoding, and facilitates the optimal exploitation of transcoding-enabled multicast opportunities. In the case without user transcoding, we optimize the transmission resource allocation to minimize the average transmission energy by exploiting natural multicast opportunities. The problem is nonconvex. We transform it to an equivalent convex problem and obtain an optimal solution using standard convex optimization techniques. In the case with user transcoding, we optimize the transmission resource allocation and the transmission quality level selection to minimize the weighted sum of the average transmission energy and the transcoding energy by exploiting both natural and transcoding-enabled multicast opportunities. The problem is a challenging mixed discrete-continuous optimization problem. We transform it to a Difference of Convex (DC) programming problem and obtain a suboptimal solution using a DC algorithm. Finally, numerical results demonstrate the importance of effective exploitation of transcoding-enabled multicast opportunities in the case with user transcoding.


title: Nuclear magnetic moments of francium 207-213 from precision hyperfine comparisons

abstract: We report a fourfold improvement in the determination of nuclear magnetic moments for neutron-deficient isotopes of francium-207--213, reducing the uncertainties from 2% for most isotopes to 0.5%. These are found by comparing our high-precision calculations of hyperfine structure constants for the ground states with experimental values. In particular, we show the importance of a careful modeling of the Bohr-Weisskopf effect, which arises due to the finite nuclear magnetization distribution. This effect is particularly large in Fr and until now has not been modeled with sufficiently high accuracy. An improved understanding of the nuclear magnetic moments and Bohr-Weisskopf effect are crucial for benchmarking the atomic theory required in precision tests of the standard model, in particular atomic parity violation studies, that are underway in francium.


title: On optimal orientations of complete tripartite graphs

abstract: Given a connected and bridgeless graph $G$, let $\mathscr{D}(G)$ be the family of strong orientations of $G$. The orientation number of $G$ is defined to be $\bar{d}(G):=min\{d(D)|D\in \mathscr{D}(G)\}$, where $d(D)$ is the diameter of the digraph $D$. In this paper, we focus on the orientation number of complete tripartite graphs. We prove a conjecture raised by Rajasekaran and Sampathkumar. Specifically, for $q\ge p\ge 3$, if $\bar{d}(K(2,p,q))=2$, then $q\le{{p}\choose{\lfloor{p/2}\rfloor}}$. We also present some sufficient conditions on $p$ and $q$ for $\bar{d}(K(p,p,q))=2$.


title: Congruence lattices of ideals in categories and (partial) semigroups

abstract: This paper presents a unified framework for determining the congruences on a number of monoids and categories of transformations, diagrams, matrices and braids, and on all their ideals. The key theoretical advances present an iterative process of stacking certain normal subgroup lattices on top of each other to successively build congruence lattices of a chain of ideals. This is applied to several specific categories of: transformations; order/orientation preserving/reversing transformations; partitions; planar/annular partitions; Brauer, Temperley--Lieb and Jones partitions; linear and projective linear transformations; and partial braids. Special considerations are needed for certain small ideals, and technically more intricate theoretical underpinnings for the linear and partial braid categories.


title: On Cross-intersecting Sperner Families

abstract: Two sets $\mathscr{A}$ and $\mathscr{B}$ are said to be cross-intersecting if $X\cap Y\neq\emptyset$ for all $X\in\mathscr{A}$ and $Y\in\mathscr{B}$. Given two cross-intersecting Sperner families (or antichains) $\mathscr{A}$ and $\mathscr{B}$ of $\mathbb{N}_n$, we prove that $|\mathscr{A}|+|\mathscr{B}|\le 2{{n}\choose{\lceil{n/2}\rceil}}$ if $n$ is odd, and $|\mathscr{A}|+|\mathscr{B}|\le {{n}\choose{n/2}}+{{n}\choose{(n/2)+1}}$ if $n$ is even. Furthermore, all extremal and almost-extremal families for $\mathscr{A}$ and $\mathscr{B}$ are determined.


title: Federated Learning for Localization: A Privacy-Preserving Crowdsourcing Method

abstract: Received Signal Strength (RSS) fingerprint-based localization has attracted a lot of research effort and cultivated many commercial applications of location-based services due to its low cost and ease of implementation. Many studies are exploring the use of deep learning (DL) algorithms for localization. DL's ability to extract features and to classify autonomously makes it an attractive solution for fingerprint-based localization. These solutions require frequent retraining of DL models with vast amounts of measurements. Although crowdsourcing is an excellent way to gather immense amounts of data, it jeopardizes the privacy of participants, as it requires to collect labeled data at a centralized server. Recently, federated learning has emerged as a practical concept in solving the privacy preservation issue of crowdsourcing participants by performing model training at the edge devices in a decentralized manner; the participants do not expose their data anymore to a centralized server. This paper presents a novel method utilizing federated learning to improve the accuracy of RSS fingerprint-based localization while preserving the privacy of the crowdsourcing participants. Employing federated learning allows ensuring \emph{preserving the privacy of user data} while enabling an adequate localization performance with experimental data captured in real-world settings. The proposed method improved localization accuracy by 1.8 meters when used as a booster for centralized learning and achieved satisfactory localization accuracy when used standalone.


title: Automated Pavement Crack Segmentation Using U-Net-based Convolutional Neural Network

abstract: Automated pavement crack image segmentation is challenging because of inherent irregular patterns, lighting conditions, and noise in images. Conventional approaches require a substantial amount of feature engineering to differentiate crack regions from non-affected regions. In this paper, we propose a deep learning technique based on a convolutional neural network to perform segmentation tasks on pavement crack images. Our approach requires minimal feature engineering compared to other machine learning techniques. We propose a U-Net-based network architecture in which we replace the encoder with a pretrained ResNet-34 neural network. We use a "one-cycle" training schedule based on cyclical learning rates to speed up the convergence. Our method achieves an F1 score of 96% on the CFD dataset and 73% on the Crack500 dataset, outperforming other algorithms tested on these datasets. We perform ablation studies on various techniques that helped us get marginal performance boosts, i.e., the addition of spatial and channel squeeze and excitation (SCSE) modules, training with gradually increasing image sizes, and training various neural network layers with different learning rates.


title: Electrical switching and oscillations in vanadium dioxide

abstract: We have studied electrical switching with S-shaped I-V characteristics in two-terminal MOM devices based on vanadium dioxide thin films. The switching effect is associated with the metal-insulator phase transition. Relaxation oscillations are observed in circuits with VO2-based switches. Dependences of the oscillator critical frequency Fmax, threshold power and voltage, as well as the time of current rise, on the switching structure size are obtained by numerical simulation. The empirical dependence of the threshold voltage on the switching region dimensions and film thickness is found. It is shown that, for the VO2 channel sizes of 10*10 nm, Fmax can reach the value of 300 MHz at a film thickness of ~20 nm. Next, it is shown that oscillatory neural networks can be implemented on the basis of coupled VO2 oscillators. For the weak capacitive coupling, we revealed the dependence of the phase difference upon synchronization on the coupling capacitance value. When the switches are scaled down, the limiting time of synchronization is reduced to Ts ~13 {\mu}s, and the number of oscillation periods for the entering to the synchronization mode remains constant, Ns ~ 17. In the case of weak thermal coupling in the synchronization mode, we observe in-phase behavior of oscillators, and there is a certain range of parameters of the supply current, in which the synchronization effect becomes possible. With a decrease in dimensions, a decrease in the thermal coupling action radius is observed, which can vary in the range from 0.5 to 50 {\mu}m for structures with characteristic dimensions of 0.1 to 5 {\mu}m, respectively. Thermal coupling may have a promising effect for realization of a 3D integrated oscillatory neural network.


title: Quantum Algorithms for the Most Frequently String Search, Intersection of Two String Sequences and Sorting of Strings Problems

abstract: We study algorithms for solving three problems on strings. The first one is the Most Frequently String Search Problem. The problem is the following. Assume that we have a sequence of $n$ strings of length $k$. The problem is finding the string that occurs in the sequence most often. We propose a quantum algorithm that has a query complexity $\tilde{O}(n \sqrt{k})$. This algorithm shows speed-up comparing with the deterministic algorithm that requires $\Omega(nk)$ queries. The second one is searching intersection of two sequences of strings. All strings have the same length $k$. The size of the first set is $n$ and the size of the second set is $m$. We propose a quantum algorithm that has a query complexity $\tilde{O}((n+m) \sqrt{k})$. This algorithm shows speed-up comparing with the deterministic algorithm that requires $\Omega((n+m)k)$ queries. The third problem is sorting of $n$ strings of length $k$. On the one hand, it is known that quantum algorithms cannot sort objects asymptotically faster than classical ones. On the other hand, we focus on sorting strings that are not arbitrary objects. We propose a quantum algorithm that has a query complexity $O(n (\log n)^2 \sqrt{k})$. This algorithm shows speed-up comparing with the deterministic algorithm (radix sort) that requires $\Omega((n+d)k)$ queries, where $d$ is a size of the alphabet.


title: A Compressed Particle-Hole Symmetric Pfaffian State for $\nu = 5/2$ Quantum Hall Effect

abstract: A recent thermal Hall conductance experiment [Banerjee et al., Nature {\bf559}, 205 (2018)] for $\nu = 5/2$ fractional quantum Hall system appears to rule out both the Pfaffian and anti-Pfaffian and be in favor of the PH-Pfaffian topological order, while the existing numerical results without disorder have shown otherwise. In this paper we offer a possible resolution by proposing a new state, termed compressed PH-Pfaffian state by "compressing" the PH-Pfaffian state with two flux quanta removed to create two abelian Laughlin type quasiparticles of the maximum avoidance from one another (or of the maximum number of zeros). The compressed PH-Pfaffian state is not particle-hole symmetric but possesses the PH-Pfaffian topological order. In spherical geometry, the compressed PH-Pfaffian state has the same magnetic flux number $N_{\phi}= 2N-3$ as the Pfaffian state, allowing a direct numerical comparison between the two states. Results of exact diagonalization of finite disorder-free systems in the second Landau level show that, by increasing the short range component of the Coulomb interaction, the ground state undergoes a phase transition from the Pfaffian state to the compressed PH-Pfaffian state before further entering into a gapless state. The low energy gapped excited states result from the breakup of the abelian Laughlin type quasiparticle into two non-abelian quasiparticles.


title: High-Performance Statistical Computing in the Computing Environments of the 2020s

abstract: Technological advances in the past decade, hardware and software alike, have made access to high-performance computing (HPC) easier than ever. We review these advances from a statistical computing perspective. Cloud computing makes access to supercomputers affordable. Deep learning software libraries make programming statistical algorithms easy and enable users to write code once and run it anywhere -- from a laptop to a workstation with multiple graphics processing units (GPUs) or a supercomputer in a cloud. Highlighting how these developments benefit statisticians, we review recent optimization algorithms that are useful for high-dimensional models and can harness the power of HPC. Code snippets are provided to demonstrate the ease of programming. We also provide an easy-to-use distributed matrix data structure suitable for HPC. Employing this data structure, we illustrate various statistical applications including large-scale positron emission tomography and $\ell_1$-regularized Cox regression. Our examples easily scale up to an 8-GPU workstation and a 720-CPU-core cluster in a cloud. As a case in point, we analyze the onset of type-2 diabetes from the UK Biobank with 200,000 subjects and about 500,000 single nucleotide polymorphisms using the HPC $\ell_1$-regularized Cox regression. Fitting this half-million-variate model takes less than 45 minutes and reconfirms known associations. To our knowledge, this is the first demonstration of the feasibility of penalized regression of survival outcomes at this scale.


title: Scalable Hybrid HMM with Gaussian Process Emission for Sequential Time-series Data Clustering

abstract: Hidden Markov Model (HMM) combined with Gaussian Process (GP) emission can be effectively used to estimate the hidden state with a sequence of complex input-output relational observations. Especially when the spectral mixture (SM) kernel is used for GP emission, we call this model as a hybrid HMM-GPSM. This model can effectively model the sequence of time-series data. However, because of a large number of parameters for the SM kernel, this model can not effectively be trained with a large volume of data having (1) long sequence for state transition and 2) a large number of time-series dataset in each sequence. This paper proposes a scalable learning method for HMM-GPSM. To effectively train the model with a long sequence, the proposed method employs a Stochastic Variational Inference (SVI) approach. Also, to effectively process a large number of data point each time-series data, we approximate the SM kernel using Reparametrized Random Fourier Feature (R-RFF). The combination of these two techniques significantly reduces the training time. We validate the proposed learning method in terms of its hidden-sate estimation accuracy and computation time using large-scale synthetic and real data sets with missing values.


title: Context-Aware Design of Cyber-Physical Human Systems (CPHS)

abstract: Recently, it has been widely accepted by the research community that interactions between humans and cyber-physical infrastructures have played a significant role in determining the performance of the latter. The existing paradigm for designing cyber-physical systems for optimal performance focuses on developing models based on historical data. The impacts of context factors driving human system interaction are challenging and are difficult to capture and replicate in existing design models. As a result, many existing models do not or only partially address those context factors of a new design owing to the lack of capabilities to capture the context factors. This limitation in many existing models often causes performance gaps between predicted and measured results. We envision a new design environment, a cyber-physical human system (CPHS) where decision-making processes for physical infrastructures under design are intelligently connected to distributed resources over cyberinfrastructure such as experiments on design features and empirical evidence from operations of existing instances. The framework combines existing design models with context-aware design-specific data involving human-infrastructure interactions in new designs, using a machine learning approach to create augmented design models with improved predictive powers.


title: The role of the GPD E in the determination of the parton angular momenta

abstract: For the generalized parton distribution (GPD) E, extracted from hard exclusive processes and the Pauli form factor of the nucleon, the second moments for quarks and gluons at $t=0$ are calculated and their scale dependence analyzed. The parton angular momenta are subsequently determined and compared with other results.


title: FedDANE: A Federated Newton-Type Method

abstract: Federated learning aims to jointly learn statistical models over massively distributed remote devices. In this work, we propose FedDANE, an optimization method that we adapt from DANE, a method for classical distributed optimization, to handle the practical constraints of federated learning. We provide convergence guarantees for this method when learning over both convex and non-convex functions. Despite encouraging theoretical results, we find that the method has underwhelming performance empirically. In particular, through empirical simulations on both synthetic and real-world datasets, FedDANE consistently underperforms baselines of FedAvg and FedProx in realistic federated settings. We identify low device participation and statistical device heterogeneity as two underlying causes of this underwhelming performance, and conclude by suggesting several directions of future work.


title: A water-obstacle separation and refinement network for unmanned surface vehicles

abstract: Obstacle detection by semantic segmentation shows a great promise for autonomous navigation in unmanned surface vehicles (USV). However, existing methods suffer from poor estimation of the water edge in the presence of visual ambiguities, poor detection of small obstacles and high false-positive rate on water reflections and wakes. We propose a new deep encoder-decoder architecture, a water-obstacle separation and refinement network (WaSR), to address these issues. Detection and water edge accuracy are improved by a novel decoder that gradually fuses inertial information from IMU with the visual features from the encoder. In addition, a novel loss function is designed to increase the separation between water and obstacle features early on in the network. Subsequently, the capacity of the remaining layers in the decoder is better utilised, leading to a significant reduction in false positives and increased true positives. Experimental results show that WaSR outperforms the current state-of-the-art by a large margin, yielding a 14% increase in F-measure over the second-best method.


title: Interpretation of normal-deformed bands in 167Lu

abstract: With the aim to get a general understanding of rotational bands in the deformed rare-earth region or in deformed nuclei in general, the observed normal-deformed rotational structures in $^{167}$Lu are interpreted within the unpaired and paired cranked Nilsson-Strutinsky formalisms, CNS and CNSB. Particular attention is devoted to the band crossings. For this nucleus with the Fermi surface high up in the $h_{11/2}$ shell, we conclude that except for the paired AB and BC crossings in configurations with an even and odd number of $i_{13/2}$ neutrons, respectively, the observed band crossings can be understood within the unpaired formalism. Especially, it means that above the AB and BC crossings, the evolution with spin is described as a gradual alignment of the spin vectors of the particles outside closed shells. Consequently, the configurations can be characterized by the number of particles occupying open $j$-shells or groups of $j$-shells. In the present study, we revise the interpretation of some experimental bands and also the nature of the crossings while some previous configuration assignments are confirmed.