title: Relating Apartness and Bisimulation

abstract: A bisimulation for a coalgebra of a functor on the category of sets can be described via a coalgebra in the category of relations, of a lifted functor. A final coalgebra then gives rise to the coinduction principle, which states that two bisimilar elements are equal. For polynomial functors, this leads to well-known descriptions. In the present paper we look at the dual notion of "apartness". Intuitively, two elements are apart if there is a positive way to distinguish them. Phrased differently: two elements are apart if and only if they are not bisimilar. Since apartness is an inductive notion, described by a least fixed point, we can give a proof system, to derive that two elements are apart. This proof system has derivation rules and two elements are apart if and only if there is a finite derivation (using the rules) of this fact. We study apartness versus bisimulation in two separate ways. First, for weak forms of bisimulation on labelled transition systems, where silent (tau) steps are included, we define an apartness notion that corresponds to weak bisimulation and another apartness that corresponds to branching bisimulation. The rules for apartness can be used to show that two states of a labelled transition system are not branching bismilar. To support the apartness view on labelled transition systems, we cast a number of well-known properties of branching bisimulation in terms of branching apartness and prove them. Next, we also study the more general categorical situation and show that indeed, apartness is the dual of bisimilarity in a precise categorical sense: apartness is an initial algebra and gives rise to an induction principle. In this analogy, we include the powerset functor, which gives a semantics to non-deterministic choice in process-theory.


title: Multi Type Mean Field Reinforcement Learning

abstract: Mean field theory provides an effective way of scaling multiagent reinforcement learning algorithms to environments with many agents that can be abstracted by a virtual mean agent. In this paper, we extend mean field multiagent algorithms to multiple types. The types enable the relaxation of a core assumption in mean field reinforcement learning, which is that all agents in the environment are playing almost similar strategies and have the same goal. We conduct experiments on three different testbeds for the field of many agent reinforcement learning, based on the standard MAgents framework. We consider two different kinds of mean field environments: a) Games where agents belong to predefined types that are known a priori and b) Games where the type of each agent is unknown and therefore must be learned based on observations. We introduce new algorithms for each type of game and demonstrate their superior performance over state of the art algorithms that assume that all agents belong to the same type and other baseline algorithms in the MAgent framework.


title: On the restricted Jordan plane in odd characteristic

abstract: In positive characteristic the Jordan plane covers a finite-dimensional Nichols algebra that was described by Cibils, Lauve and Witherspoon and we call the restricted Jordan plane. In this paper the characteristic is odd. The defining relations of the Drinfeld double of the restricted Jordan plane are presented and its simple modules are determined. A Hopf algebra that deserves the name of double of the Jordan plane is introduced and various quantum Frobenius maps are described. The finite-dimensional pre-Nichols algebras intermediate between the Jordan plane and its restricted version are classified. The defining relations of the graded dual of the Jordan plane are given.


title: Quasi-Equivalence of Width and Depth of Neural Networks

abstract: While classic studies proved that wide networks allow universal approximation, recent research and successes of deep learning demonstrate the power of deep networks. Based on a symmetric consideration, we investigate if the design of artificial neural networks should have a directional preference, and what the mechanism of interaction is between the width and depth of a network. Inspired by the De Morgan law, we address this fundamental question by establishing a quasi-equivalence between the width and depth of ReLU networks in two aspects. First, we formulate two transforms for mapping an arbitrary ReLU network to a wide network and a deep network respectively for either regression or classification so that the essentially same capability of the original network can be implemented. Then, we replace the mainstream artificial neuron type with a quadratic counterpart, and utilize the factorization and continued fraction representations of the same polynomial function to construct a wide network and a deep network, respectively. Based on our findings, a deep network has a wide equivalent, and vice versa, subject to an arbitrarily small error.


title: Breaking the $O(\sqrt n)$-Bit Barrier: Byzantine Agreement with Polylog Bits Per Party

abstract: Byzantine agreement (BA), the task of $n$ parties to agree on one of their input bits in the face of malicious agents, is a powerful primitive that lies at the core of a vast range of distributed protocols. Interestingly, in protocols with the best overall communication, the demands of the parties are highly unbalanced: the amortized cost is $\tilde O(1)$ bits per party, but some parties must send $\Omega(n)$ bits. In best known balanced protocols, the overall communication is sub-optimal, with each party communicating $\tilde O(\sqrt{n})$. In this work, we ask whether asymmetry is inherent for optimizing total communication. Our contributions in this line are as follows: 1) We define a cryptographic primitive, succinctly reconstructed distributed signatures (SRDS), that suffices for constructing $\tilde O(1)$ balanced BA. We provide two constructions of SRDS from different cryptographic and Public-Key Infrastructure (PKI) assumptions. 2) The SRDS-based BA follows a paradigm of boosting from "almost-everywhere" agreement to full agreement, and does so in a single round. We prove that PKI setup and cryptographic assumptions are necessary for such protocols in which every party sends $o(n)$ messages. 3) We further explore connections between a natural approach toward attaining SRDS and average-case succinct non-interactive argument systems (SNARGs) for a particular type of NP-Complete problems (generalizing Subset-Sum and Subset-Product). Our results provide new approaches forward, as well as limitations and barriers, towards minimizing per-party communication of BA. In particular, we construct the first two BA protocols with $\tilde O(1)$ balanced communication, offering a tradeoff between setup and cryptographic assumptions, and answering an open question presented by King and Saia (DISC'09).


title: Hybrid Solver for the Radiative Transport Equation Using Finite Volume and Discontinuous Galerkin

abstract: We propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.


title: Provably Efficient Online Hyperparameter Optimization with Population-Based Bandits

abstract: Many of the recent triumphs in machine learning are dependent on well-tuned hyperparameters. This is particularly prominent in reinforcement learning (RL) where a small change in the configuration can lead to failure. Despite the importance of tuning hyperparameters, it remains expensive and is often done in a naive and laborious way. A recent solution to this problem is Population Based Training (PBT) which updates both weights and hyperparameters in a single training run of a population of agents. PBT has been shown to be particularly effective in RL, leading to widespread use in the field. However, PBT lacks theoretical guarantees since it relies on random heuristics to explore the hyperparameter space. This inefficiency means it typically requires vast computational resources, which is prohibitive for many small and medium sized labs. In this work, we introduce the first provably efficient PBT-style algorithm, Population-Based Bandits (PB2). PB2 uses a probabilistic model to guide the search in an efficient way, making it possible to discover high performing hyperparameter configurations with far fewer agents than typically required by PBT. We show in a series of RL experiments that PB2 is able to achieve high performance with a modest computational budget.


title: Data-Driven False Data Injection Attacks Against Power Grids: A Random Matrix Approach

abstract: We address the problem of constructing false data injection (FDI) attacks that can bypass the bad data detector (BDD) of a power grid. The attacker is assumed to have access to only power flow measurement data traces (collected over a limited period of time) and no other prior knowledge about the grid. Existing related algorithms are formulated under the assumption that the attacker has access to measurements collected over a long (asymptotically infinite) time period, which may not be realistic. We show that these approaches do not perform well when the attacker has a limited number of data samples only. We design an enhanced algorithm to construct FDI attack vectors in the face of limited measurements that can nevertheless bypass the BDD with high probability. The algorithm design is guided by results from random matrix theory. Furthermore, we characterize an important trade-off between the attack's BDD-bypass probability and its sparsity, which affects the spatial extent of the attack that must be achieved. Extensive simulations using data traces collected from the MATPOWER simulator and benchmark IEEE bus systems validate our findings.

categories: cs.CR cs.IT cs.SY eess.SP eess.SY math.IT


title: Robust Multi-channel Speech Recognition using Frequency Aligned Network

abstract: Conventional speech enhancement technique such as beamforming has known benefits for far-field speech recognition. Our own work in frequency-domain multi-channel acoustic modeling has shown additional improvements by training a spatial filtering layer jointly within an acoustic model. In this paper, we further develop this idea and use frequency aligned network for robust multi-channel automatic speech recognition (ASR). Unlike an affine layer in the frequency domain, the proposed frequency aligned component prevents one frequency bin influencing other frequency bins. We show that this modification not only reduces the number of parameters in the model but also significantly and improves the ASR performance. We investigate effects of frequency aligned network through ASR experiments on the real-world far-field data where users are interacting with an ASR system in uncontrolled acoustic environments. We show that our multi-channel acoustic model with a frequency aligned network shows up to 18% relative reduction in word error rate.


title: Enhancement of shock-capturing methods via machine learning

abstract: In recent years, machine learning has been used to create data-driven solutions to problems for which an algorithmic solution is intractable, as well as fine-tuning existing algorithms. This research applies machine learning to the development of an improved finite-volume method for simulating PDEs with discontinuous solutions. Shock capturing methods make use of nonlinear switching functions that are not guaranteed to be optimal. Because data can be used to learn nonlinear relationships, we train a neural network to improve the results of a fifth-order WENO method. We post-process the outputs of the neural network to guarantee that the method is consistent. The training data consists of the exact mapping between cell averages and interpolated values for a set of integrable functions that represent waveforms we would expect to see while simulating a PDE. We demonstrate our method on linear advection of a discontinuous function, the inviscid Burgers' equation, and the 1-D Euler equations. For the latter, we examine the Shu-Osher model problem for turbulence-shockwave interactions. We find that our method outperforms WENO in simulations where the numerical solution becomes overly diffused due to numerical viscosity.

categories: physics.comp-ph cs.LG cs.NA cs.NE math.NA physics.flu-dyn


title: Link Capacity Distributions and Optimal Capacities for Competent Network Performance

abstract: This work addresses the problem of evaluating optimal link capacities of a packet-flow network for the objective of congestion minimization. We present a simple model of packet flow in networks and present a numerical approach to evaluate packet flow probability mass function at any arbitrary edge of the network for a given routing algorithm and traffic rate. We further discuss techniques of assigning optimal capacity at each edge for attaining desired minimized congestion and discuss related trade-offs. Our framework is built around the assumption of Poisson traffic, however the numerical approach fits for any general distribution of packet influx. Lastly, we define metrics of global performance of link capacities allocation and discuss the effect of network structure on capacity allocation and performance.


title: Max Min vertex cover and the size of Betti tables

abstract: Let $G$ be a finite simple graph on $n$ vertices, that contains no isolated vertices, and let $I(G) \subseteq S = K[x_1, \dots, x_n]$ be its edge ideal. In this paper, we study the pair of integers that measure the projective dimension and the regularity of $S/I(G)$. We show that if the projective dimension of $S/I(G)$ attains its minimum value $2\sqrt{n}-2$ then, with only one exception, the its regularity must be 1. We also provide a full description for the spectrum of the projective dimension of $S/I(G)$ when the regularity attains its minimum value 1.


title: High-Accuracy Multi-Node Ranging For Coherent Distributed Antenna Arrays

abstract: The design and experimental implementation of a waveform for high-accuracy inter-node ranging in a coherent distributed antenna array is presented. Based on a spectrally-sparse high-accuracy ranging waveform, the presented multi-frequency waveform enables high-accuracy ranging between multiple nodes in an array simultaneously, without interference. The waveform is based on a unique time-frequency duplexing approach combining a stepped-frequency waveform with different step cycles per node pair. The waveform also inherently includes beneficial disambiguation properties. The ambiguity function of the waveform is derived, and theoretical bounds on the ranging accuracy are obtained. Measurements were conducted in software-defined radio-based nodes in a three-element distributed array, demonstrating high-accuracy unambiguous ranging between two slave nodes and one master node.


title: Interpolating Predictors in High-Dimensional Factor Regression

abstract: This work studies finite-sample properties of the risk of the minimum-norm interpolating predictor in high-dimensional regression models. If the effective rank of the covariance matrix $\Sigma$ of the $p$ regression features is much larger than the sample size $n$, we show that the min-norm interpolating predictor is not desirable, as its risk approaches the risk of trivially predicting the response by 0. However, our detailed finite-sample analysis reveals, surprisingly, that this behavior is not present when the regression response and the features are {\it jointly} low-dimensional, following a widely used factor regression model. Within this popular model class, and when the effective rank of $\Sigma$ is smaller than $n$, while still allowing for $p \gg n$, both the bias and the variance terms of the excess risk can be controlled, and the risk of the minimum-norm interpolating predictor approaches optimal benchmarks. Moreover, through a detailed analysis of the bias term, we exhibit model classes under which our upper bound on the excess risk approaches zero, while the corresponding upper bound in the recent work arXiv:1906.11300 diverges. Furthermore, we show that the minimum-norm interpolating predictor analyzed under the factor regression model, despite being model-agnostic and devoid of tuning parameters, can have similar risk to predictors based on principal components regression and ridge regression, and can improve over LASSO based predictors, in the high-dimensional regime.


title: How to Answer Why -- Evaluating the Explanations of AI Through Mental Model Analysis

abstract: To achieve optimal human-system integration in the context of user-AI interaction it is important that users develop a valid representation of how AI works. In most of the everyday interaction with technical systems users construct mental models (i.e., an abstraction of the anticipated mechanisms a system uses to perform a given task). If no explicit explanations are provided by a system (e.g. by a self-explaining AI) or other sources (e.g. an instructor), the mental model is typically formed based on experiences, i.e. the observations of the user during the interaction. The congruence of this mental model and the actual systems functioning is vital, as it is used for assumptions, predictions and consequently for decisions regarding system use. A key question for human-centered AI research is therefore how to validly survey users' mental models. The objective of the present research is to identify suitable elicitation methods for mental model analysis. We evaluated whether mental models are suitable as an empirical research method. Additionally, methods of cognitive tutoring are integrated. We propose an exemplary method to evaluate explainable AI approaches in a human-centered way.


title: Synthetic Magnetic Resonance Images with Generative Adversarial Networks

abstract: Data augmentation is essential for medical research to increase the size of training datasets and achieve better results. In this work, we experiment three GAN architectures with different loss functions to generate new brain MRIs. The results show the importance of hyperparameter tuning and the use of mini-batch similarity layer in the Discriminator and gradient penalty in the loss function to achieve convergence with high quality and realism. Moreover, huge computation time is needed to generate indistinguishable images from the original dataset.


title: On generalized residue network for deep learning of unknown dynamical systems

abstract: We present a general numerical approach for learning unknown dynamical systems using deep neural networks (DNNs). Our method is built upon recent studies that identified the residue network (ResNet) as an effective neural network structure. In this paper, we present a generalized ResNet framework and broadly define residue as the discrepancy between observation data and prediction made by another model, which can be an existing coarse model or reduced-order model. In this case, the generalized ResNet serves as a model correction to the existing model and recovers the unresolved dynamics. When an existing coarse model is not available, we present numerical strategies for fast creation of coarse models, to be used in conjunction with the generalized ResNet. These coarse models are constructed using the same data set and thus do not require additional resources. The generalized ResNet is capable of learning the underlying unknown equations and producing predictions with accuracy higher than the standard ResNet structure. This is demonstrated via several numerical examples, including long-term prediction of a chaotic system.


title: A Repeated Game Freeway Lane Changing Model

abstract: Lane changes are complex safety and throughput critical driver actions. Most lane changing models deal with lane-changing maneuvers solely from the merging driver's standpoint and thus ignore driver interaction. To overcome this shortcoming, we develop a game-theoretical decision-making model and validate the model using empirical merging maneuver data at a freeway on-ramp. Specifically, this paper advances our repeated game model in a previous paper by using updated payoff functions. Validation results using the NGSIM empirical data show that the developed game-theoretical model provides better prediction accuracy compared to previous work, with correct predictions approximately 86 percent of the time. In addition, a sensitivity analysis demonstrates the rationality and sensitivity of the model to variations in various factors. To provide evidence of the benefits of the repeated game approach, which takes into account previous decision-making results, a case study is conducted using an agent-based simulation model. The proposed repeated game model produces superior performance to a one-shot game model, when simulating actual freeway merging behaviors. Finally, this lane change model, which captures the collective decision-making between human drivers, can be used to develop automated vehicle driving strategies.


title: Machine learning on DNA-encoded libraries: A new paradigm for hit-finding

abstract: DNA-encoded small molecule libraries (DELs) have enabled discovery of novel inhibitors for many distinct protein targets of therapeutic value through screening of libraries with up to billions of unique small molecules. We demonstrate a new approach applying machine learning to DEL selection data by identifying active molecules from a large commercial collection and a virtual library of easily synthesizable compounds. We train models using only DEL selection data and apply automated or automatable filters with chemist review restricted to the removal of molecules with potential for instability or reactivity. We validate this approach with a large prospective study (nearly 2000 compounds tested) across three diverse protein targets: sEH (a hydrolase), ER{\alpha} (a nuclear receptor), and c-KIT (a kinase). The approach is effective, with an overall hit rate of {\sim}30% at 30 {\textmu}M and discovery of potent compounds (IC50 <10 nM) for every target. The model makes useful predictions even for molecules dissimilar to the original DEL and the compounds identified are diverse, predominantly drug-like, and different from known ligands. Collectively, the quality and quantity of DEL selection data; the power of modern machine learning methods; and access to large, inexpensive, commercially-available libraries creates a powerful new approach for hit finding.


title: Homogenization of two-phase flow in porous media from Pore to Darcy Scale: A phase-field approach

abstract: We extend the two-scale expansion approach of periodic homogenization to include time scales and thus can tackle the full instationary Navier-Stokes-Cahn-Hilliard model at the pore scale as microscale. Time scale separation allows us to keep microscale dynamics, responsible e.g. for hysteresis, and arrive at a numerically tractable micro-macro model including coupled generalized Darcy's laws.


title: Singularities of transient processes in dynamics and beyond

abstract: This note is a brief review of the analysis of long transients in dynamical systems. The problem of long transients arose in many disciplines, from physical and chemical kinetic to biology and even social sciences. Detailed analysis of singularities of various `relaxation times' associated long transients with bifurcations of $\omega$-limit sets, homoclinic structures (intersections of $\alpha$- and $\omega$-limit sets) and other peculiarities of dynamics. This review was stimulated by the analysis of anomalously long transients in ecology published recently by A. Morozov and S. Petrovskii with co-authors.


title: Understanding the dynamics of message passing algorithms: a free probability heuristics

abstract: We use freeness assumptions of random matrix theory to analyze the dynamical behavior of inference algorithms for probabilistic models with dense coupling matrices in the limit of large systems. For a toy Ising model, we are able to recover previous results such as the property of vanishing effective memories and the analytical convergence rate of the algorithm.

categories: cond-mat.stat-mech cond-mat.dis-nn cs.LG stat.ML


title: Fast inference of Boosted Decision Trees in FPGAs for particle physics

abstract: We describe the implementation of Boosted Decision Trees in the hls4ml library, which allows the translation of a trained model into FPGA firmware through an automated conversion process. Thanks to its fully on-chip implementation, hls4ml performs inference of Boosted Decision Tree models with extremely low latency. With a typical latency less than 100 ns, this solution is suitable for FPGA-based real-time processing, such as in the Level-1 Trigger system of a collider experiment. These developments open up prospects for physicists to deploy BDTs in FPGAs for identifying the origin of jets, better reconstructing the energies of muons, and enabling better selection of rare signal processes.

categories: physics.comp-ph astro-ph.IM cs.LG hep-ex


title: Kinematics of a simple reciprocal model swimmer at intermediate Reynolds numbers

abstract: We computationally study the kinematics of a simple model reciprocal swimmer (asymmetric dumbbell) as a function of the Reynolds number (Re) and investigate how the onset and gradual increase of inertia impacts the swimming behavior: a reversal in the swim direction, flow directions, and the swim stroke. We divide the swim stroke into the expansion and compression of the two spheres and relate them to power and recovery strokes. We find that the switch in swim direction also corresponds to a switch in power and recovery strokes. We obtain expressions for the mean swimming velocity by collapsing the net displacement during expansion and compression under power law relationships with respect to Re, the swimmer's amplitude, and the distance between the two spheres. Analyzing the fluid flows, we see the averaged flow field during expansion always resembles a pusher and compression always a puller, but when averaged over the whole cycle, the flow that dominates is the one that occurs during the power stroke. We also relate the power and recovery strokes to the swimming efficiency during times of expansion and compression, and we find that the power stroke is, surprisingly, not always more efficient than the recovery stroke. Our results may have important implications for biology and ultimately the design of artificial swimmers.


title: Constructive Hybrid Games

abstract: Hybrid games are models which combine discrete, continuous, and adversarial dynamics. Game logic enables proving (classical) existence of winning strategies. We introduce constructive differential game logic (CdGL) for hybrid games, where proofs that a player can win the game correspond to computable winning strategies. This is the logical foundation for synthesis of correct control and monitoring code for safety-critical cyber-physical systems. Our contributions include novel static and dynamic semantics as well as soundness and consistency.


title: Quantum walks and Dirac cellular automata on a programmable trapped-ion quantum computer

abstract: The quantum walk formalism is a widely used and highly successful framework for modeling quantum systems, such as simulations of the Dirac equation, different dynamics in both the low and high energy regime, and for developing a wide range of quantum algorithms. Here we present the circuit-based implementation of a discrete-time quantum walk in position space on a five-qubit trapped-ion quantum processor. We encode the space of walker positions in particular multi-qubit states and program the system to operate with different quantum walk parameters, experimentally realizing a Dirac cellular automaton with tunable mass parameter. The quantum walk circuits and position state mapping scale favorably to a larger model and physical systems, allowing the implementation of any algorithm based on discrete-time quantum walks algorithm and the dynamics associated with the discretized version of the Dirac equation.


title: Sim2Real2Sim: Bridging the Gap Between Simulation and Real-World in Flexible Object Manipulation

abstract: This paper addresses a new strategy called Simulation-to-Real-to-Simulation (Sim2Real2Sim) to bridge the gap between simulation and real-world, and automate a flexible object manipulation task. This strategy consists of three steps: (1) using the rough environment with the estimated models to develop the methods to complete the manipulation task in the simulation; (2) applying the methods from simulation to real-world and comparing their performance; (3) updating the models and methods in simulation based on the differences between the real world and the simulation. The Plug Task from the 2015 DARPA Robotics Challenge Finals is chosen to evaluate our Sim2Real2Sim strategy. A new identification approach for building the model of the linear flexible objects is derived from real-world to simulation. The automation of the DRC plug task in both simulation and real-world proves the success of the Sim2Real2Sim strategy. Numerical experiments are implemented to validate the simulated model.


title: Weyl Superconductivity in UTe2

abstract: The search for a material platform for topological quantum computation has recently focused on unconventional superconductors. Such material systems, where the superconducting order parameter breaks a symmetry of the crystal point group, are capable of hosting novel phenomena, including emergent Majorana quasiparticles. Unique among unconventional superconductors is the recently discovered UTe2, where spin-triplet superconductivity emerges from a paramagnetic normal state. Although UTe2 could be considered a relative of a family of known ferromagnetic superconductors, the unique crystal structure of this material and experimentally suggested zero Curie temperature pose a great challenge to determining the symmetries, magnetism, and topology underlying the superconducting state. These emergent properties will determine the utility of UTe2 for future spintronics and quantum information applications. Here, we report observations of a non-zero polar Kerr effect and of two transitions in the specific heat upon entering the superconducting state, which together show that the superconductivity in UTe2 is characterized by an order parameter with two components that breaks time reversal symmetry. These data allow us to place firm constraints on the symmetries of the order parameter, which strongly suggest that UTe2 is a Weyl superconductor that hosts chiral Fermi arc surface states.


title: Obstruction to a Higman embedding theorem for residually finite groups with solvable word problem

abstract: We prove that, for a finitely generated residually finite group, having solvable word problem is not a sufficient condition to be a subgroup of a finitely presented residually finite group. The obstruction is given by a residually finite group with solvable word problem for which there is no effective method that allows, given some non-identity element, to find a morphism onto a finite group in which this element has a non-trivial image. We also prove that the depth function of this group grows faster than any recursive function.


title: Computability of finite quotients of finitely generated groups

abstract: We study systematically groups whose marked finite quotients form a recursive set. We give several definitions, and prove basic properties of this class of groups, and in particular emphasize the link between the growth of the depth function and solvability of the word problem. We give examples of infinitely presented groups whose finite quotients can be effectively enumerated. Finally, our main result is that a residually finite group can be even not recursively presented and still have computable finite quotients, and that, on the other hand, it can have solvable word problem while still not having computable finite quotients.


title: Exciton g-factors of van der Waals heterostructures from first principles calculations

abstract: External fields are a powerful tool to probe optical excitations in a material. The linear energy shift of an excitation in a magnetic field is quantified by its effective g-factor. Here we show how exciton g-factors and their sign can be determined by converged first principles calculations. We apply the method to monolayer excitons in semiconducting transition metal dichalcogenides and to interlayer excitons in MoSe$_2$/WSe$_2$ heterobilayers and obtain good agreement with recent experimental data. The precision of our method allows to assign measured g-factors of optical peaks to specific transitions in the band structure and also to specific regions of the samples. This revealed the nature of various, previously measured interlayer exciton peaks. We further show that, due to specific optical selection rules, g-factors in van der Waals heterostructures are strongly spin- and stacking-dependent. The calculation of orbital angular momenta requires the summation over hundreds of bands, indicating that for the considered two-dimensional materials the basis set size is a critical numerical issue. The presented approach can potentially be applied to a wide variety of semiconductors.


title: Dimerization and N\'eel order in different quantum spin chains through a shared loop representation

abstract: The ground states of the spin-$ S $ antiferromagnetic chain $H_\textrm{AF}$ with a projection-based interaction and the spin-$ 1/2$ XXZ-chain $ H_\textrm{XXZ} $ at anisotropy parameter $\Delta=\cosh(\lambda) $ share a common loop representation in terms of a two-dimensional functional integral which is similar to the classical planar $Q$-state Potts model at $ \sqrt Q= 2S+1 =2\cosh(\lambda)$. The multifaceted relation is used here to directly relate the distinct forms of translation symmetry breaking which are manifested in the ground states of these two models: dimerization for $H_\textrm{AF}$ at all $S> 1/2$, and N\'eel order for $ H_\textrm{XXZ} $ at $\lambda >0$. The results presented include: i) a translation to the above quantum spin systems of the results which were recently proven by Duminil-Copin-Li-Manolescu for a broad class of two-dimensional random-cluster models, and ii) a short proof of the symmetry breaking in a manner similar to the recent structural proof by Ray-Spinka of the discontinuity of the phase transition for $Q>4$. Altogether, the quantum manifestation of the change between $Q=4$ and $Q>4$ is a transition from a gapless ground state to a pair of gapped and extensively distinct ground states.

categories: math-ph cond-mat.stat-mech math.MP math.PR


title: The Voltage Regulation of a Buck Converter Using a Neural Network Predictive Controller

abstract: In this paper, a neural network predictive controller (NNPC) is proposed to control a buck converter. Conventional controllers such as proportional integral (PI) or proportional integral derivative (PID) are designed based on the linearized small signal model near the operating point. Therefore, the performance of the controller in the start up, load change, or reference change is not optimal since the system model changes by changing the operating point. The neural network predictive controller optimally controls the buck converter by following the concept of the traditional model predictive controller. The advantage of the NNPC is that the neural network system identification decreases the inaccuracy of the system model with inaccurate parameters. A NNPC with a well trained neural network can perform as an optimal controller for the buck converter. To compare the effectiveness of the traditional buck converter and the NNPC, the simulation results are provided.


title: Contradictory Structure Learning for Semi-supervised Domain Adaptation

abstract: Current adversarial adaptation methods attempt to align the cross-domain features, whereas two challenges remain unsolved: 1) the conditional distribution mismatch and 2) the bias of the decision boundary towards the source domain. To solve these challenges, we propose a novel framework for semi-supervised domain adaptation by unifying the learning of opposite structures (UODA). UODA consists of a generator and two classifiers (i.e., the source-scattering classifier and the target-clustering classifier), which are trained for contradictory purposes. The target-clustering classifier attempts to cluster the target features to improve intra-class density and enlarge inter-class divergence. Meanwhile, the source-scattering classifier is designed to scatter the source features to enhance the decision boundary's smoothness. Through the alternation of source-feature expansion and target-feature clustering procedures, the target features are well-enclosed within the dilated boundary of the corresponding source features. This strategy can make the cross-domain features to be precisely aligned against the source bias simultaneously. Moreover, to overcome the model collapse through training, we progressively update the measurement of feature's distance and their representation via an adversarial training paradigm. Extensive experiments on the benchmarks of DomainNet and Office-home datasets demonstrate the superiority of our approach over the state-of-the-art methods.


title: Convective boundary mixing in low- and intermediate-mass stars I. Core properties from pressure-mode asteroseismology

abstract: Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot' measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraints from pressure modes to infer the CBM properties of Kepler and CoRoT main-sequence and subgiant oscillators, as well as in two radial velocity targets (Procyon A and $\alpha$ Cen A). Collectively these targets allow us to identify how measurement precision, stellar spectral type, and overshoot implementation impact the asteroseismic solution. With these new measures we find that the `effective overshoot' for most stars is in line with physical expectations and calibrations from binaries and clusters. However, two F-stars in the CoRoT field (HD 49933 and HD 181906) still necessitate high overshoot in the models. Due to short mode lifetimes, mode identification can be difficult in these stars. We demonstrate that an incongruence between the radial and non-radial modes drives the asteroseismic solution to extreme structures with highly-efficient CBM as an inevitable outcome. Understanding the cause of seemingly anomalous physics for such stars is vital for inferring accurate stellar parameters from TESS data with comparable timeseries length.


title: Closing the Dequantization Gap: PixelCNN as a Single-Layer Flow

abstract: Flow models have recently made great progress at modeling ordinal discrete data such as images and audio. Due to the continuous nature of flow models, dequantization is typically applied when using them for such discrete data, resulting in lower bound estimates of the likelihood. In this paper, we introduce subset flows, a class of flows that can tractably transform finite volumes and thus allow exact computation of likelihoods for discrete data. Based on subset flows, we identify ordinal discrete autoregressive models, including WaveNets, PixelCNNs and Transformers, as single-layer flows. We use the flow formulation to compare models trained and evaluated with either the exact likelihood or its dequantization lower bound. Finally, we study multilayer flows composed of PixelCNNs and non-autoregressive coupling layers and demonstrate state-of-the-art results on CIFAR-10 for flow models trained with dequantization.


title: Temperature scaling of two-ion anisotropy in pure and mixed anisotropy systems

abstract: Magnetic anisotropy plays an essential role in information technology applications of magnetic materials, providing a means to retain the long-term stability of a magnetic state in the presence of thermal fluctuations. Anisotropy consists of a single-ion contribution stemming from the crystal structure and two-ion terms attributed to the exchange interactions between magnetic atoms. A lack of robust theory crucially limits the understanding of the temperature dependence of the anisotropy in pure two-ion and mixed single-ion and two-ion systems. Here, we use Green's function theory and atomistic Monte Carlo simulations to determine the temperature scaling of the effective anisotropy in ferromagnets in these pure and mixed cases, from saturated to vanishing magnetization. At low temperature, we find that the pure two-ion anisotropy scales with the reduced magnetization as $k(m) \sim m^{2.28}$, while the mixed scenario describes the diversity of the temperature dependence of the anisotropy observed in real materials. The temperature dependence of the mixed anisotropy determined here is expected to considerably contribute to the understanding and the control of the thermal properties of magnetic materials.


title: The Impact of Planetary Rotation Rate on the Reflectance and Thermal Emission Spectrum of Terrestrial Exoplanets Around Sun-like Stars

abstract: Robust atmospheric and radiative transfer modeling will be required to properly interpret reflected light and thermal emission spectra of terrestrial exoplanets. This will help break observational degeneracies between the numerous atmospheric, planetary, and stellar factors that drive planetary climate. Here we simulate the climates of Earth-like worlds around the Sun with increasingly slow rotation periods, from Earth-like to fully Sun-synchronous, using the ROCKE-3D general circulation model. We then provide these results as input to the Spectral Planet Model (SPM), which employs the SMART radiative transfer model to simulate the spectra of a planet as it would be observed from a future space-based telescope. We find that the primary observable effects of slowing planetary rotation rate are the altered cloud distributions, altitudes, and opacities which subsequently drive many changes to the spectra by altering the absorption band depths of biologically-relevant gas species (e.g., H2O, O2, and O3). We also identify a potentially diagnostic feature of synchronously rotating worlds in mid-infrared H2O absorption/emission lines.


title: Nullities for a class of skew-symmetric Toeplitz band matrices

abstract: For all $n > k \ge 1$, we give formulas for the nullity $N(n,k)$ of the $n \times n$ skew-symmetric Toeplitz band matrix whose first $k$ superdiagonals have all entries $1$ and whose remaining superdiagonals have all entries $0$. This is accomplished by counting the number of cycles in certain directed graphs. As an application, for each fixed integer $z\ge 0$ and large fixed $k$, we give an asymptotic formula for the percentage of $n > k$ satisfying $N(n,k)=z$. For the purpose of rapid computation, an algorithm is devised that quickly computes $N(n,k)$ even for extremely large values of $n$ and $k$.


title: Conformal infinitesimal variations of submanifolds

abstract: This paper belongs to the realm of conformal geometry and deals with Euclidean submanifolds that admit smooth variations that are infinitesimally conformal. Conformal variations of Euclidean submanifolds is a classical subject in differential geometry. In fact, already in 1917 Cartan classified parametrically the Euclidean hypersurfaces that admit nontrivial conformal variations. Our first main result is a Fundamental theorem for conformal infinitesimal variations. The second is a rigidity theorem for Euclidean submanifolds that lie in low codimension.


title: Progressive Cleaning and Mining of Uncertain Smart Water Meter Data

abstract: Several municipalities have recently installed wireless 'smart' water meters that allow functionalities such as demand response, leak alerts, identification of characteristic demand patterns, and detailed consumption analysis. To achieve these benefits, the meter data needs to be error-free, which is not necessarily available in practice, due to 'dirtiness' or 'uncertainty' of data, which is mostly unavoidable. The focus of this paper is to investigate practical solutions to mine uncertain data for reliable results and to evaluate the impact of dirty data on filters. This evaluation would eventually lead to valuable information, which can be used for educated decision making on water planning strategies. We perform a systematic study of the errors existing in a large-scale smart water meter deployments, which is helpful to better understand the nature of errors. Identifying customers contributing to a load peak is used as the main filter. The filter outputs are then combined with the domain expert knowledge to evaluate their accuracy and validity and also to look for potential errors. After discovering each error, we analyze its trails in the data and track back its source, which would eventually lead to the removal of the error or dealing with it accordingly. This procedure is applied progressively to ensure that all detectable errors are discovered and characterized in the data model. We evaluate the performance of the proposed approach using the smart water meter consumption data obtained from the City of Abbotsford, British Columbia, Canada. We present the results of both unprocessed and cleaned data and analyze, in detail, the sensitivity of the selected filter to the errors.


title: Orbit-induced spin precession as an origin of periodicity in periodically-repeating fast radio bursts

abstract: FRB 180916.J0158+65 has been found to repeatedly emit fast radio bursts with the period in roughly 16 days. We propose that such periodicity comes from orbit-induced, spin precession of the emitter, which is possibly a neutron star. Depending on the mass of the companion, the binary period ranges from several hundred to thousands of seconds. Such tight binaries have a relatively short lifetime, which does not likely come from the gravitational decay of a wide binary. We comment on the relation to GW190425 and the possibility in LISA and LIGO detections.


title: Cartesian differential categories as skew enriched categories

abstract: We exhibit the cartesian differential categories of Blute, Cockett and Seely as a particular kind of enriched category. The base for the enrichment is the category of commutative monoids -- or in a straightforward generalisation, the category of modules over a commutative rig $k$. However, the tensor product on this category is not the usual one, but rather a warping of it by a certain monoidal comonad $Q$. Thus the enrichment base is not a monoidal category in the usual sense, but rather a skew monoidal category in the sense of Szlach\'anyi. Our first main result is that cartesian differential categories are the same as categories with finite products enriched over this skew monoidal base. The comonad $Q$ involved is, in fact, an example of a differential modality. Differential modalities are a kind of comonad on a symmetric monoidal $k$-linear category with the characteristic feature that their co-Kleisli categories are cartesian differential categories. Using our first main result, we are able to prove our second one: that every small cartesian differential category admits a full, structure-preserving embedding into the cartesian differential category induced by a differential modality (in fact, a monoidal differential modality on a monoidal closed category -- thus, a model of intuitionistic differential linear logic). This resolves an important open question in this area.


title: Axisymmetric deformations of neutron stars and gravitational-wave astronomy

abstract: Einstein's theory of general relativity predicts that the only stationary configuration of an isolated black hole is the Kerr spacetime, which has a unique multipolar structure and a spherical shape when non-spinning. This is in striking contrast to the case of other self-gravitating objects, which instead can in principle have arbitrary deformations even in the static case. Here we develop a general perturbative framework to construct stationary stars with small axisymmetric deformations, and study explicitly compact stars with an intrinsic quadrupole moment. The latter can be sustained, for instance, by crust stresses or strong magnetic fields. While our framework is general, we focus on quadrupolar deformations of neutron stars induced by an anisotropic crust, which continuously connect to spherical neutron stars in the isotropic limit. Deformed neutron stars might provide a more accurate description for stellar remnants formed in supernovae and in binary mergers, and can be used to improve constraints on the neutron-star equation of state through gravitational-wave detections and through the observation of low-mass X-ray binaries. We argue that, if the (dimensionless) intrinsic quadrupole moment is of a few percent or higher, the effect of the deformation is stronger than that of tidal interactions in coalescing neutron-star binaries, and might also significantly affect the electromagnetic signal from accreting neutron stars.


title: Radial processes for sub-Riemannian Brownian motions and applications

abstract: We study the radial part of sub-Riemannian Brownian motion in the context of totally geodesic foliations. It\^o's formula is proved for the radial processes associated to Riemannian distances approximating the Riemannian one. We deduce very general stochastic completeness criteria for the sub-Riemannian Brownian motion. In the context of Sasakian foliations and H-type groups, one can push the analysis further, and taking advantage of the recently proved sub-Laplacian comparison theorems one can compare the radial processes for the sub-Riemannian distance to one-dimensional model diffusions. As a geometric application, we prove Cheng's type estimates for the Dirichlet eigenvalues of the sub-Riemannian metric balls, a result which seems to be new even in the Heisenberg group.


title: JPLink: On Linking Jobs to Vocational Interest Types

abstract: Linking job seekers with relevant jobs requires matching based on not only skills, but also personality types. Although the Holland Code also known as RIASEC has frequently been used to group people by their suitability for six different categories of occupations, the RIASEC category labels of individual jobs are often not found in job posts. This is attributed to significant manual efforts required for assigning job posts with RIASEC labels. To cope with assigning massive number of jobs with RIASEC labels, we propose JPLink, a machine learning approach using the text content in job titles and job descriptions. JPLink exploits domain knowledge available in an occupation-specific knowledge base known as O*NET to improve feature representation of job posts. To incorporate relative ranking of RIASEC labels of each job, JPLink proposes a listwise loss function inspired by learning to rank. Both our quantitative and qualitative evaluations show that JPLink outperforms conventional baselines. We conduct an error analysis on JPLink's predictions to show that it can uncover label errors in existing job posts.


title: Wolf-Rayet stars in M81 using GTC/OSIRIS: 7 new detections, analysis and classification of the full sample

abstract: We report the detection of 7 new Wolf-Rayet (WR) star locations in M81 using the Multi-Object Spectrograph of the OSIRIS instrument at Gran Telescopio Canarias. These detections are the result of a follow-up of an earlier study using the same instrumental set-up that resulted in the detection of 14 WR locations. We analyse the entire sample of 21 spectra to classify them to one of the known WR sub-types using template spectra of WR stars in the Large Magellanic Cloud (LMC), with similar metallicity to M81. Taking into consideration the dispersion in the strengths of the bumps for a given WR sub-type, we found that 19 of the 21 locations correspond to individual stars, including all the 7 new detections, of sub-types: WNL, WNE, WCE and the transitional WN/C. None of the detections correspond to WCL or WO types. The positions of these stars in the red bump vs blue bump luminosity diagram agrees well with an evolutionary path according to the Conti scenario. Based on this, we propose this diagram as a straightforward tool for spectral classification of extragalactic WR sources. The detection of individual WR stars in M81, which is at a distance of 3.6 Mpc, opens up a new environment for testing the massive star evolutionary models.


title: Impact of ImageNet Model Selection on Domain Adaptation

abstract: Deep neural networks are widely used in image classification problems. However, little work addresses how features from different deep neural networks affect the domain adaptation problem. Existing methods often extract deep features from one ImageNet model, without exploring other neural networks. In this paper, we investigate how different ImageNet models affect transfer accuracy on domain adaptation problems. We extract features from sixteen distinct pre-trained ImageNet models and examine the performance of twelve benchmarking methods when using the features. Extensive experimental results show that a higher accuracy ImageNet model produces better features, and leads to higher accuracy on domain adaptation problems (with a correlation coefficient of up to 0.95). We also examine the architecture of each neural network to find the best layer for feature extraction. Together, performance from our features exceeds that of the state-of-the-art in three benchmark datasets.


title: Global stability of a Caputo fractional SIRS model with general incidence rate

abstract: We introduce a fractional order SIRS model with non-linear incidence rate. Existence of a unique positive solution to the model is proved. Stability analysis of the disease free equilibrium and positive fixed points are investigated. Finally, a numerical example is presented.


title: Spectrum Dependent Learning Curves in Kernel Regression and Wide Neural Networks

abstract: We derive analytical expressions for the generalization performance of kernel regression as a function of the number of training samples using theoretical methods from Gaussian processes and statistical physics. Our expressions apply to wide neural networks due to an equivalence between training them and kernel regression with the Neural Tangent Kernel (NTK). By computing the decomposition of the total generalization error due to different spectral components of the kernel, we identify a new spectral principle: as the size of the training set grows, kernel machines and neural networks fit successively higher spectral modes of the target function. When data are sampled from a uniform distribution on a high-dimensional hypersphere, dot product kernels, including NTK, exhibit learning stages where different frequency modes of the target function are learned. We verify our theory with simulations on synthetic data and MNIST dataset.


title: Transformer Transducer: A Streamable Speech Recognition Model with Transformer Encoders and RNN-T Loss

abstract: In this paper we present an end-to-end speech recognition model with Transformer encoders that can be used in a streaming speech recognition system. Transformer computation blocks based on self-attention are used to encode both audio and label sequences independently. The activations from both audio and label encoders are combined with a feed-forward layer to compute a probability distribution over the label space for every combination of acoustic frame position and label history. This is similar to the Recurrent Neural Network Transducer (RNN-T) model, which uses RNNs for information encoding instead of Transformer encoders. The model is trained with the RNN-T loss well-suited to streaming decoding. We present results on the LibriSpeech dataset showing that limiting the left context for self-attention in the Transformer layers makes decoding computationally tractable for streaming, with only a slight degradation in accuracy. We also show that the full attention version of our model beats the-state-of-the art accuracy on the LibriSpeech benchmarks. Our results also show that we can bridge the gap between full attention and limited attention versions of our model by attending to a limited number of future frames.


title: Breaking Band: A Breakdown of High-performance Communication

abstract: The critical path of internode communication on large-scale systems is composed of multiple components. When a supercomputing application initiates the transfer of a message using a high-level communication routine such as an MPI_Send, the payload of the message traverses multiple software stacks, the I/O subsystem on both the host and target nodes, and network components such as the switch. In this paper, we analyze where, why, and how much time is spent on the critical path of communication by modeling the overall injection overhead and end-to-end latency of a system. We focus our analysis on the performance of small messages since fine-grained communication is becoming increasingly important with the growing trend of an increasing number of cores per node. The analytical models present an accurate and detailed breakdown of time spent in internode communication. We validate the models on Arm ThunderX2-based servers connected with Mellanox InfiniBand. This is the first work of this kind on Arm. Alongside our breakdown, we describe the methodology to measure the time spent in each component so that readers with access to precise CPU timers and a PCIe analyzer can measure breakdowns on systems of their interest. Such a breakdown is crucial for software developers, system architects, and researchers to guide their optimization efforts. As researchers ourselves, we use the breakdown to simulate the impacts and discuss the likelihoods of a set of optimizations that target the bottlenecks in today's high-performance communication.


title: Empirical Bayes for Large-scale Randomized Experiments: a Spectral Approach

abstract: Large-scale randomized experiments, sometimes called A/B tests, are increasingly prevalent in many industries. Though such experiments are often analyzed via frequentist $t$-tests, arguably such analyses are deficient: $p$-values are hard to interpret and not easily incorporated into decision-making. As an alternative, we propose an empirical Bayes approach, which assumes that the treatment effects are realized from a "true prior". This requires inferring the prior from previous experiments. Following Robbins, we estimate a family of marginal densities of empirical effects, indexed by the noise scale. We show that this family is characterized by the heat equation. We develop a spectral maximum likelihood estimate based on a Fourier series representation, which can be efficiently computed via convex optimization. In order to select hyperparameters and compare models, we describe two model selection criteria. We demonstrate our method on simulated and real data, and compare posterior inference to that under a Gaussian mixture model of the prior.


title: Borel structures on the space of left-orderings

abstract: In this paper we study the Borel structure of the space of left-orderings $\mathrm{LO}(G)$ of a group $G$ modulo the natural conjugacy action, and by using tools from descriptive set theory we find many examples of countable left-orderable groups such that the quotient space $\mathrm{LO}(G)/G$ is not standard. This answers a question of Deroin, Navas, and Rivas. We also prove that the countable Borel equivalence relation induced from the conjugacy action of $\mathbb{F}_{2}$ on $\mathrm{LO}(\mathbb{F}_{2})$ is universal, and leverage this result to provide many other examples of countable left-orderable groups $G$ such that the natural $G$-action on $\mathrm{LO}(G)$ induces a universal countable Borel equivalence relation.


title: Disjoint weighing matrices

abstract: The notion of disjoint weighing matrices is introduced as a generalization of orthogonal designs. A recursive construction along with a computer search lead to some infinite classes of disjoint weighing matrices, which in turn are shown to form commutative association schemes with 3 or 4 classes.


title: Stability and Scalability of Blockchain Systems

abstract: The blockchain paradigm provides a mechanism for content dissemination and distributed consensus on Peer-to-Peer (P2P) networks. While this paradigm has been widely adopted in industry, it has not been carefully analyzed in terms of its network scaling with respect to the number of peers. Applications for blockchain systems, such as cryptocurrencies and IoT, require this form of network scaling. In this paper, we propose a new stochastic network model for a blockchain system. We identify a structural property called \emph{one-endedness}, which we show to be desirable in any blockchain system as it is directly related to distributed consensus among the peers. We show that the stochastic stability of the network is sufficient for the one-endedness of a blockchain. We further establish that our model belongs to a class of network models, called monotone separable models. This allows us to establish upper and lower bounds on the stability region. The bounds on stability depend on the connectivity of the P2P network through its conductance and allow us to analyze the scalability of blockchain systems on large P2P networks. We verify our theoretical insights using both synthetic data and real data from the Bitcoin network.


title: High Temporal Resolution Rainfall Runoff Modelling Using Long-Short-Term-Memory (LSTM) Networks

abstract: Accurate and efficient models for rainfall runoff (RR) simulations are crucial for flood risk management. Most rainfall models in use today are process-driven; i.e. they solve either simplified empirical formulas or some variation of the St. Venant (shallow water) equations. With the development of machine-learning techniques, we may now be able to emulate rainfall models using, for example, neural networks. In this study, a data-driven RR model using a sequence-to-sequence Long-short-Term-Memory (LSTM) network was constructed. The model was tested for a watershed in Houston, TX, known for severe flood events. The LSTM network's capability in learning long-term dependencies between the input and output of the network allowed modeling RR with high resolution in time (15 minutes). Using 10-years precipitation from 153 rainfall gages and river channel discharge data (more than 5.3 million data points), and by designing several numerical tests the developed model performance in predicting river discharge was tested. The model results were also compared with the output of a process-driven model Gridded Surface Subsurface Hydrologic Analysis (GSSHA). Moreover, physical consistency of the LSTM model was explored. The model results showed that the LSTM model was able to efficiently predict discharge and achieve good model performance. When compared to GSSHA, the data-driven model was more efficient and robust in terms of prediction and calibration. Interestingly, the performance of the LSTM model improved (test Nash-Sutcliffe model efficiency from 0.666 to 0.942) when a selected subset of rainfall gages based on the model performance, were used as input instead of all rainfall gages.


title: Reentrant Kondo effect in a quantum impurity coupled to a metal-semiconductor hybrid contact

abstract: Using NRG, we show that a system containing a quantum impurity (QI), strongly coupled to a semiconductor (gap $2 \Delta$) and weakly coupled to a metal, displays a 'reentrant' Kondo stage at low temperatures. The NRG analysis of the corresponding Single Impurity Anderson Model (SIAM) shows that the reentrant stage is characterized by a second sequence of SIAM fixed points: free orbital (FO) > local moment (LM) > strong coupling (SC). In the first stage, the SC fixed point (Kondo temperature $T_{K1}$) is unstable, while the second stage exhibits a much lower Kondo temperature $T_{K2}$, associated to a stable SC fixed point. The results indicate that the reentrant Kondo screening is associated to an effective SIAM, with an effective repulsion $U_{eff}$. This low temperature effective SIAM, which we dub as 'reentrant' SIAM, behaves as a 'replica' of the high temperature (bare) SIAM. The intuitive picture that emerges is that the first Kondo state develops through impurity screening by semiconducting electrons, while the second Kondo state involves screening by metallic electrons, once the semiconducting electrons are out of reach to thermal excitations ($T < \Delta$) and only the metallic spectral weight inside the gap is available for impurity screening. In addition, we analyze a hybrid system formed by a QI `sandwiched' between an armchair graphene nanoribbon (AGNR) and a scanning tunneling microscope (STM) tip, with respective couplings set to reproduce the generic model described above. The energy gap in the AGNR can be externally tuned by an electric-field-induced Rashba spin-orbit interaction. We analyzed this system for realistic parameter values, using NRG, and concluded that the reentrant SIAM, with its associated second stage Kondo, is worthy of experimental investigation.


title: AGB winds in interacting binary stars

abstract: We perform numerical simulations to investigate the stellar wind from interacting binary stars. Our aim is to find analytical formulae describing the outflow structure. In each binary system the more massive star is in the asymptotic giant branch and its wind is driven by a combination of pulsations in the stellar surface layers and radiation pressure on dust, while the less massive star is in the main sequence. Time averages of density and outflow velocity of the stellar wind are calculated and plotted as profiles against distance from the centre of mass and colatitude angle. We find that mass is lost mainly through the outer Lagrangian point L2. The resultant outflow develops into a spiral at low distances from the binary. The outflowing spiral is quickly smoothed out by shocks and becomes an excretion disk at larger distances. This leads to the formation of an outflow structure with an equatorial density excess, which is greater in binaries with smaller orbital separation. The pole-to-equator density ratio reaches a maximum value of $\sim10^5$ at Roche-Lobe Overflow state. We also find that the gas stream leaving L2 does not form a circumbinary ring for stellar mass ratios above 0.78, when radiation pressure on dust is taken into account. Analytical formulae are obtained by curve fitting the 2-dimensional, azimuthally averaged density and outflow velocity profiles. The formulae can be used in future studies to setup the initial outflow structure in hydrodynamic simulations of common-envelope evolution and formation of planetary nebulae.


title: Adaptive semiparametric Bayesian differential equations via sequential Monte Carlo

abstract: Nonlinear differential equations (DEs) are used in a wide range of scientific problems to model complex dynamic systems. The differential equations often contain unknown parameters that are of scientific interest, which have to be estimated from noisy measurements of the dynamic system. Generally, there is no closed-form solution for nonlinear DEs, and the likelihood surface for the parameter of interest is multi-modal and very sensitive to different parameter values. We propose a Bayesian framework for nonlinear DE systems. A flexible nonparametric function is used to represent the dynamic process such that expensive numerical solvers can be avoided. A sequential Monte Carlo algorithm in the annealing framework is proposed to conduct Bayesian inference for parameters in DEs. In our numerical experiments, we use examples of ordinary differential equations and delay differential equations to demonstrate the effectiveness of the proposed algorithm. We developed an R package that is available at \url{https://github.com/shijiaw/smcDE}.


title: Multimodal Controller for Generative Models

abstract: Class-conditional generative models are crucial tools for data generation from user-specified class labels. Existing approaches for class-conditional generative models require nontrivial modifications of backbone generative architectures to model conditional information fed into the model. This paper introduces a plug-and-play module named `multimodal controller' to generate multimodal data without introducing additional learning parameters. In the absence of the controllers, our model reduces to non-conditional generative models. We test the efficacy of multimodal controllers on CIFAR10, COIL100, and Omniglot benchmark datasets. We demonstrate that multimodal controlled generative models (including VAE, PixelCNN, Glow, and GAN) can generate class-conditional images of significantly better quality when compared with conditional generative models. Moreover, we show that multimodal controlled models can also create novel modalities of images.


title: Occurrence Rates of Planets Orbiting M Stars: Applying ABC to Kepler DR25, Gaia DR2, and 2MASS Data

abstract: We present robust planet occurrence rates for Kepler planet candidates around M stars for planet radii $R_p = 0.5-4~\textrm{R}_\oplus$ and orbital periods $P = 0.5-256$ days using the approximate Bayesian computation (ABC) technique. This work incorporates the final Kepler DR25 planet candidate catalog and data products and augments them with updated stellar properties using Gaia DR2 and 2MASS PSC. We apply a set of selection criteria to select a sample of 1,746 Kepler M dwarf targets that host 89 associated planet candidates. These early M dwarfs and late K dwarfs were selected from cross-referenced targets using several photometric quality flags from Gaia DR2 and color-magnitude cuts using 2MASS magnitudes. We estimate a habitable zone occurrence rate of $f_{\textrm{M,HZ}} = 0.33^{+0.10}_{-0.12}$ for planets with $0.75-1.5$ R$_\oplus$ size. We caution that occurrence rate estimates for Kepler M stars are sensitive to the choice of prior due to the small sample of target stars and planet candidates. For example, we find an occurrence rate of $4.2^{+0.6}_{-0.6}$ or $8.4^{+1.2}_{-1.1}$ planets per M dwarf (integrating over $R_p = 0.5-4~\textrm{R}_\oplus$ and $P = 0.5-256$ days) for our two choices of prior. These occurrence rates are greater than those for FGK dwarf target when compared at the same range of orbital periods, but similar to occurrence rates when computed as a function of equivalent stellar insolation. Combining our result with recent studies of exoplanet architectures indicates that most, and potentially all, early-M dwarfs harbor planetary systems.


title: The Power of Linear Controllers in LQR Control

abstract: The Linear Quadratic Regulator (LQR) framework considers the problem of regulating a linear dynamical system perturbed by environmental noise. We compute the policy regret between three distinct control policies: i) the optimal online policy, whose linear structure is given by the Ricatti equations; ii) the optimal offline linear policy, which is the best linear state feedback policy given the noise sequence; and iii) the optimal offline policy, which selects the globally optimal control actions given the noise sequence. We fully characterize the optimal offline policy and show that it has a recursive form in terms of the optimal online policy and future disturbances. We also show that cost of the optimal offline linear policy converges to the cost of the optimal online policy as the time horizon grows large, and consequently the optimal offline linear policy incurs linear regret relative to the optimal offline policy, even in the optimistic setting where the noise is drawn i.i.d from a known distribution. Although we focus on the setting where the noise is stochastic, our results also imply new lower bounds on the policy regret achievable when the noise is chosen by an adaptive adversary.


title: Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations

abstract: We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.


title: Refining Constructive Hybrid Games

abstract: We extend the constructive differential game logic (CdGL) of hybrid games with a refinement connective that relates two hybrid games. We use this connective to prove a folk theorem relating hybrid games to hybrid systems.


title: Thin-shell wormholes in Einstein and Einstein-Gauss-Bonnet theories of gravity

abstract: We review recent works on the possibility for eternal existence of thin-shell wormholes on Einstein and Einstein-Gauss-Bonnet gravity. We introduce thin-shell wormholes that are categorized into a class of traversable wormhole solutions. After that, we discuss stable thin-shell wormholes with negative-tension branes in Reissner-Nordstr\"om-(anti) de Sitter spacetimes in $d$ dimensional Einstein gravity. Imposing $Z_2$ symmetry, we construct and classify traversable static thin-shell wormholes in spherical, planar and hyperbolic symmetries. It is found that the spherical wormholes are stable against spherically symmetric perturbations. It is also found that some classes of wormholes in planar and hyperbolic symmetries with a negative cosmological constant are stable against perturbations preserving symmetries. In most cases, stable wormholes are found with the appropriate combination of an electric charge and a negative cosmological constant. However, as special cases, there are stable wormholes even with a vanishing cosmological constant in spherical symmetry and with a vanishing electric charge in hyperbolic symmetry. Subsequently, the existence and dynamical stability of traversable thin-shell wormholes with electrically neutral negative-tension branes is discussed in Einstein-Gauss-Bonnet theory of gravitation. We consider radial perturbations against the shell for the solutions, which have the $Z_2$ symmetry. The effect of the Gauss-Bonnet term on the stability depends on the spacetime symmetry.


title: The threshold bias of the clique-factor game

abstract: Let $r \ge 4$ be an integer and consider the following game on the complete graph $K_n$ for $n \in r \mathbb{Z}$: Two players, Maker and Breaker, alternately claim previously unclaimed edges of $K_n$ such that in each turn Maker claims one and Breaker claims $b \in \mathbb{N}$ edges. Maker wins if her graph contains a $K_r$-factor, that is a collection of $n/r$ vertex-disjoint copies of $K_r$, and Breaker wins otherwise. In other words, we consider a $b$-biased $K_r$-factor Maker-Breaker game. We show that the threshold bias for this game is of order $n^{2/(r+2)}$. This makes a step towards determining the threshold bias for making bounded-degree spanning graphs and extends a result of Allen et al.\ who resolved the case $r \in \{3,4\}$ up to a logarithmic factor.


title: Estimating Optimal Treatment Rules with an Instrumental Variable: A Partial Identification Learning Approach

abstract: Individualized treatment rules (ITRs) are considered a promising recipe to deliver better policy interventions. One key ingredient in optimal ITR estimation problems is to estimate the average treatment effect conditional on a subject's covariate information, which is often challenging in observational studies due to the universal concern of unmeasured confounding. Instrumental variables (IVs) are widely-used tools to infer the treatment effect when there is unmeasured confounding between the treatment and outcome. In this work, we propose a general framework of approaching the optimal ITR estimation problem when a valid IV is allowed to only partially identify the treatment effect. We introduce a novel notion of optimality called "IV-optimality". A treatment rule is said to be IV-optimal if it minimizes the maximum risk with respect to the putative IV and the set of IV identification assumptions. We derive a bound on the risk of an IV-optimal rule that illuminates when an IV-optimal rule has favorable generalization performance. We propose a classification-based statistical learning method that estimates such an IV-optimal rule, design computationally-efficient algorithms, and prove theoretical guarantees. We contrast our proposed method to the popular outcome weighted learning (OWL) approach via extensive simulations, and apply our method to study which mothers would benefit from traveling to deliver their premature babies at hospitals with high level neonatal intensive care units.


title: Autonomous Industrial Assembly using Force, Torque, and RGB-D sensing

abstract: We present algorithms and results for a robotic manipulation system that was designed to be easily programmable and adaptable to various tasks common to industrial setting, which is inspired by the Industrial Assembly Challenge at the 2018 World Robotics Summit in Tokyo. This challenge included assembly of standard, commercially available industrial parts into 2D and 3D assemblies. We demonstrate three tasks that can be classified into "peg-in-hole" and "hole-on-peg" tasks and identify two canonical algorithms: spiral-based search and tilting insertion. Both algorithms use hand-coded thresholds in the force and torque domains to detect critical points in the assembly. After briefly summarizing the state of the art in research, we describe the strategy and approach utilized by the tested system, how it's design bears on its performance, statistics on 20 experimental trials for each task, lessons learned during the development of the system, and open research challenges that still remain.


title: Dynamic Energy Dispatch Based on Deep Reinforcement Learning in IoT-Driven Smart Isolated Microgrids

abstract: Microgrids (MGs) are small, local power grids that can operate independently from the larger utility grid. Combined with the Internet of Things (IoT), a smart MG can leverage the sensory data and machine learning techniques for intelligent energy management. This paper focuses on deep reinforcement learning (DRL)-based energy dispatch for IoT-driven smart isolated MGs with diesel generators (DGs), photovoltaic (PV) panels, and a battery. A finite-horizon Partial Observable Markov Decision Process (POMDP) model is formulated and solved by learning from historical data to capture the uncertainty in future electricity consumption and renewable power generation. In order to deal with the instability problem of DRL algorithms and unique characteristics of finite-horizon models, two novel DRL algorithms, namely, finite-horizon deep deterministic policy gradient (FH-DDPG) and finite-horizon recurrent deterministic policy gradient (FH-RDPG), are proposed to derive energy dispatch policies with and without fully observable state information. A case study using real isolated MG data is performed, where the performance of the proposed algorithms are compared with the other baseline DRL and non-DRL algorithms. Moreover, the impact of uncertainties on MG performance is decoupled into two levels and evaluated respectively.


title: Quantifying the Value of Lateral Views in Deep Learning for Chest X-rays

abstract: Most deep learning models in chest X-ray prediction utilize the posteroanterior (PA) view due to the lack of other views available. PadChest is a large-scale chest X-ray dataset that has almost 200 labels and multiple views available. In this work, we use PadChest to explore multiple approaches to merging the PA and lateral views for predicting the radiological labels associated with the X-ray image. We find that different methods of merging the model utilize the lateral view differently. We also find that including the lateral view increases performance for 32 labels in the dataset, while being neutral for the others. The increase in overall performance is comparable to the one obtained by using only the PA view with twice the amount of patients in the training set.


title: The microscopic relationships between triangular arbitrage and cross-currency correlations in a simple agent based model of foreign exchange markets

abstract: Foreign exchange rates movements exhibit significant cross-correlations even on very short time-scales. The effect of these statistical relationships become evident during extreme market events, such as flash crashes.In this scenario, an abrupt price swing occurring on a given market is immediately followed by anomalous movements in several related foreign exchange rates. Although a deep understanding of cross-currency correlations would be clearly beneficial for conceiving more stable and safer foreign exchange markets, the microscopic origins of these interdependencies have not been extensively investigated. We introduce an agent-based model which describes the emergence of cross-currency correlations from the interactions between market makers and an arbitrager. Our model qualitatively replicates the time-scale vs. cross-correlation diagrams observed in real trading data, suggesting that triangular arbitrage plays a primary role in the entanglement of the dynamics of different foreign exchange rates. Furthermore, the model shows how the features of the cross-correlation function between two foreign exchange rates, such as its sign and value, emerge from the interplay between triangular arbitrage and trend-following strategies.


title: Explicit Mean-Square Error Bounds for Monte-Carlo and Linear Stochastic Approximation

abstract: This paper concerns error bounds for recursive equations subject to Markovian disturbances. Motivating examples abound within the fields of Markov chain Monte Carlo (MCMC) and Reinforcement Learning (RL), and many of these algorithms can be interpreted as special cases of stochastic approximation (SA). It is argued that it is not possible in general to obtain a Hoeffding bound on the error sequence, even when the underlying Markov chain is reversible and geometrically ergodic, such as the M/M/1 queue. This is motivation for the focus on mean square error bounds for parameter estimates. It is shown that mean square error achieves the optimal rate of $O(1/n)$, subject to conditions on the step-size sequence. Moreover, the exact constants in the rate are obtained, which is of great value in algorithm design.

categories: math.PR cs.LG cs.SY eess.SY math.OC math.ST stat.ML stat.TH


title: Learning Hyperspectral Feature Extraction and Classification with ResNeXt Network

abstract: The Hyperspectral image (HSI) classification is a standard remote sensing task, in which each image pixel is given a label indicating the physical land-cover on the earth's surface. The achievements of image semantic segmentation and deep learning approaches on ordinary images have accelerated the research on hyperspectral image classification. Moreover, the utilization of both the spectral and spatial cues in hyperspectral images has shown improved classification accuracy in hyperspectral image classification. The use of only 3D Convolutional Neural Networks (3D-CNN) to extract both spatial and spectral cues from Hyperspectral images results in an explosion of parameters hence high computational cost. We propose network architecture called the MixedSN that utilizes the 3D convolutions to modeling spectral-spatial information in the early layers of the architecture and the 2D convolutions at the top layers which majorly deal with semantic abstraction. We constrain our architecture to ResNeXt block because of their performance and simplicity. Our model drastically reduced the number of parameters and achieved comparable classification performance with state-of-the-art methods on Indian Pine (IP) scene dataset, Pavia University scene (PU) dataset, Salinas (SA) Scene dataset, and Botswana (BW) dataset.


title: Traveling wave solutions for a class of discrete diffusive SIR epidemic model

abstract: This paper is concerned with the conditions of existence and nonexistence of traveling wave solutions (TWS) for a class of discrete diffusive epidemic models. We find that the existence of TWS is determined by the so-called basic reproduction number and the critical wave speed: When the basic reproduction number R0 greater than 1, there exists a critical wave speed c* > 0, such that for each c >= c * the system admits a nontrivial TWS and for c < c* there exists no nontrivial TWS for the system. In addition, the boundary asymptotic behaviour of TWS is obtained by constructing a suitable Lyapunov functional and employing Lebesgue dominated convergence theorem. Finally, we apply our results to two discrete diffusive epidemic models to verify the existence and nonexistence of TWS.


title: Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets

abstract: The superposition of atomic potentials (SAP) approach has recently been shown to be a simple and efficient way to initialize electronic structure calculations [S. Lehtola, J. Chem. Theory Comput. 15, 1593 (2019)]. Here, we study the differences between effective potentials from fully numerical density functional and optimized effective potential calculations for fixed configurations. We find that the differences are small, overall, and choose exchange-only potentials at the local density approximation level of theory computed on top of Hartree-Fock densities as a good compromise. The differences between potentials arising from different atomic configurations are also found to be small at this level of theory. Furthermore, we discuss the efficient Gaussian-basis implementation of SAP via error function fits to fully numerical atomic radial potentials. The guess obtained from the fitted potentials can be easily implemented in any Gaussian-basis quantum chemistry code in terms of two-electron integrals. Fits covering the whole periodic table from H to Og are reported for non-relativistic as well as fully relativistic four-component calculations that have been carried out with fully numerical approaches.


title: Observation of flat bands due to band hybridization in 3d-electron heavy-fermion compound CaCu3Ru4O12

abstract: We report angle-resolved photoemission spectroscopy and first-principles numerical calculations for the band structure evolution of the 3d heavy-fermion compound CaCu3Ru4O12. Below 200 K, we observed an emergent hybridization gap between the Cu 3d electron-like band and the Ru 4d hole-like band and the resulting flat band features near the Fermi energy centered around the Brillouin zone corner. Our results confirm the non-Kondo nature of CaCu3Ru4O12, in which the Cu 3dxy electrons are less correlated and not in the Kondo limit. Comparison between theory and experiment also suggests that other mechanism such as nonlocal interactions or spin fluctuations beyond the local dynamical mean-field theory may be needed in order to give a quantitative explanation of the peculiar properties in this material.


title: Poisson Kernel Avoiding Self-Smoothing in Graph Convolutional Networks

abstract: Graph convolutional network (GCN) is now an effective tool to deal with non-Euclidean data, such as social networks in social behavior analysis, molecular structure analysis in the field of chemistry, and skeleton-based action recognition. Graph convolutional kernel is one of the most significant factors in GCN to extract nodes' feature, and some improvements of it have reached promising performance theoretically and experimentally. However, there is limited research about how exactly different data types and graph structures influence the performance of these kernels. Most existing methods used an adaptive convolutional kernel to deal with a given graph structure, which still not reveals the internal reasons. In this paper, we started from theoretical analysis of the spectral graph and studied the properties of existing graph convolutional kernels. While taking some designed datasets with specific parameters into consideration, we revealed the self-smoothing phenomenon of convolutional kernels. After that, we proposed the Poisson kernel that can avoid self-smoothing without training any adaptive kernel. Experimental results demonstrate that our Poisson kernel not only works well on the benchmark dataset where state-of-the-art methods work fine, but also is evidently superior to them in synthetic datasets.


title: A Time Delay Dynamic System with External Source for the Local Outbreak of 2019-nCoV

abstract: How to model the 2019 CoronaVirus (2019-nCov) spread in China is one of the most urgent and interesting problems in applied mathematics. In this paper, we propose a novel time delay dynamic system with external source to describe the trend of local outbreak for the 2019-nCoV. The external source is introduced in the newly proposed dynamic system, which can be considered as the suspected people travel to different areas. The numerical simulations exhibit the dynamic system with the external source is more reliable than the one without it, and the rate of isolation is extremely important for controlling the increase of cumulative confirmed people of 2019-nCoV. Based on our numerical simulation results with the public data, we suggest that the local government should have some more strict measures to maintain the rate of isolation. Otherwise the local cumulative confirmed people of 2019-nCoV might be out of control.


title: Long-Range Gesture Recognition Using Millimeter Wave Radar

abstract: Millimeter wave (mmWave) based gesture recognition technology provides a good human computer interaction (HCI) experience. Prior works focus on the close-range gesture recognition, but fall short in range extension, i.e., they are unable to recognize gestures more than one meter away from considerable noise motions. In this paper, we design a long-range gesture recognition model which utilizes a novel data processing method and a customized artificial Convolutional Neural Network (CNN). Firstly, we break down gestures into multiple reflection points and extract their spatial-temporal features which depict gesture details. Secondly, we design a CNN to learn changing patterns of extracted features respectively and output the recognition result. We thoroughly evaluate our proposed system by implementing on a commodity mmWave radar. Besides, we also provide more extensive assessments to demonstrate that the proposed system is practical in several real-world scenarios.


title: Equivalence relations and $L^p$ distances between time series with application to the Black Summer Australian bushfires

abstract: This paper introduces a new framework of algebraic equivalence relations between time series and new distance metrics between them, then applies these to investigate the Australian ``Black Summer'' bushfire season of 2019-2020. First, we introduce a general framework for defining equivalence between time series, heuristically intended to be equivalent if they differ only up to noise. Our first specific implementation is based on using change point algorithms and comparing statistical quantities such as mean or variance in stationary segments. We thus derive the existence of such equivalence relations on the space of time series, such that the quotient spaces can be equipped with a metrizable topology. Next, we illustrate specifically how to define and compute such distances among a collection of time series and perform clustering and additional analysis thereon. Then, we apply these insights to analyze air quality data across New South Wales, Australia, during the 2019-2020 bushfires. There, we investigate structural similarity with respect to this data and identify locations that were impacted anonymously by the fires relative to their location. This may have implications regarding the appropriate management of resources to avoid gaps in the defense against future fires.


title: Phonon spectrum of underdoped $\text{HgBa}_2\text{CuO}_{4+\delta}$ investigated by neutron scattering

abstract: The cuprates exhibit a prominent charge-density-wave (CDW) instability with wavevector along [100], i.e., the Cu-O bond direction. Whereas CDW order is most prominent at moderate doping and low temperature, there exists increasing evidence for dynamic charge correlations throughout a large portion of the temperature-doping phase diagram. In particular, signatures of incipient charge order have been observed as phonon softening and/or broadening near the CDW wavevector approximately half-way through the Brillouin zone. Most of this work has focused on moderately-doped cuprates, for which the CDW order is robust, or on optimally-doped samples, for which the superconducting transition temperature ($T_c$) attains its maximum. Here we present a time-of-flight neutron scattering study of phonons in simple-tetragonal $\text{HgBa}_2\text{CuO}_{4+\delta}$ ($T_c = 55$ K) at a low doping level where prior work showed the CDW order to be weak. We employ and showcase a new software-based technique that mines the large number of measured Brillouin zones for useful data in order to improve accuracy and counting statistics. Density-functional theory has not provided an accurate description of phonons in $\text{HgBa}_2\text{CuO}_{4+\delta}$, yet we find the right set of parameters to qualitatively reproduce the data. The notable exception is a dispersion minimum in the longitudinal Cu-O bond-stretching branch along [100]. This discrepancy suggests that, while CDW order is weak, there exist significant dynamic charge correlations in the optic phonon range at low doping, near the edge of the superconducting dome.


title: Distribution free testing for linear regression. Extension to general parametric regression

abstract: Recently a distribution free approach for testing parametric hypotheses based on unitary transformations has been suggested in \cite{Khm13, Khm16, Khm17} and further studied in \cite{Ngu17} and \cite{Rob19}. In this note we show that the transformation takes extremely simple form in distribution free testing of linear regression. Then we extend it to general parametric regression with vector-valued covariates.


title: ML Estimation and MAP Estimation for Device Activities in Grant-Free Random Access with Interference

abstract: Device activity detection is one main challenge in grant-free random access, which is recently proposed to support massive access for massive machine-type communications (mMTC). Existing solutions fail to consider interference generated by massive Internet of Things (IoT) devices, or important prior information on device activities and interference. In this paper, we consider device activity detection at an access point (AP) in the presence of interference generated by massive devices from other cells. We consider the joint maximum likelihood (ML) estimation and the joint maximum a posterior probability (MAP) estimation of both the device activities and interference powers, jointly utilizing tools from probability, stochastic geometry and optimization. Each estimation problem is a difference of convex (DC) programming problem, and a coordinate descent algorithm is proposed to obtain a stationary point. The proposed ML estimation extends the existing ML estimation by considering the estimation of interference powers together with the estimation of device activities. The proposed MAP estimation further enhances the proposed ML estimation by exploiting prior distributions of device activities and interference powers. Numerical results show the substantial gains of the proposed joint estimation designs, and reveal the importance of explicit consideration of interference and the value of prior information in device activity detection.