title: New robust ScaLAPACK routine for computing the QR factorization with column pivoting

abstract: In this note we describe two modifications of the ScaLAPACK subroutines PxGEQPF for computing the QR factorization with the Businger-Golub column pivoting. First, we resolve a subtle numerical instability in the same way as we have done it for the LAPACK subroutines xGEQPF, xGEQP3 in 2006. [LAPACK Working Note 176 (2006); ACM Trans. Math. Softw. 2008]. The problem originates in the first release of LINPACK in the 1970's: due to severe cancellations in the down-dating of partial column norms, the pivoting procedure may be in the dark completely about the true norms of the pivot column candidates. This may cause miss-pivoting, and as a result loss of the important rank revealing structure of the computed triangular factor, with severe consequences on other solvers that rely on the rank revealing pivoting. The instability is so subtle that, e.g., inserting a WRITE statement or changing the process topology can drastically change the result. Secondly, we also correct a programming error in the complex subroutines PCGEQPF, PZGEQPF, which also causes wrong pivoting because of erroneous use of PSCNRM2, PDZNRM2 for the explicit norm computation.


title: A Research Platform for Multi-Robot Dialogue with Humans

abstract: This paper presents a research platform that supports spoken dialogue interaction with multiple robots. The demonstration showcases our crafted MultiBot testing scenario in which users can verbally issue search, navigate, and follow instructions to two robotic teammates: a simulated ground robot and an aerial robot. This flexible language and robotic platform takes advantage of existing tools for speech recognition and dialogue management that are compatible with new domains, and implements an inter-agent communication protocol (tactical behavior specification), where verbal instructions are encoded for tasks assigned to the appropriate robot.


title: What You See May Not Be What You Get: UCB Bandit Algorithms Robust to {\epsilon}-Contamination

abstract: Motivated by applications of bandit algorithms in education, we consider a stochastic multi-armed bandit problem with $\varepsilon$-contaminated rewards. We allow an adversary to give arbitrary unbounded contaminated rewards with full knowledge of the past and future. We impose the constraint that for each time $t$ the proportion of contaminated rewards for any action is less than or equal to $\varepsilon$. We derive concentration inequalities for two robust mean estimators for sub-Gaussian distributions in the $\varepsilon$-contamination context. We define the $\varepsilon$-contaminated stochastic bandit problem and use our robust mean estimators to give two variants of a robust Upper Confidence Bound (UCB) algorithm, crUCB. Using regret derived from only the underlying stochastic rewards, both variants of crUCB achieve $\mathcal{O} (\sqrt{KT\log T})$ regret for small enough contamination proportions. Our simulations assume small horizons, reflecting the newly explored setting of bandits in education. We show that in certain adversarial regimes crUCB not only outperforms algorithms designed for stochastic (UCB1) and adversarial (EXP3) bandits but also those that have "best of both worlds" guarantees (EXP3++ and TsallisInf) even when our constraint on the proportion of contaminated rewards is broken.


title: Isolation and Localization of Unknown Faults Using Neural Network-Based Residuals

abstract: Localization of unknown faults in industrial systems is a difficult task for data-driven diagnosis methods. The classification performance of many machine learning methods relies on the quality of training data. Unknown faults, for example faults not represented in training data, can be detected using, for example, anomaly classifiers. However, mapping these unknown faults to an actual location in the real system is a non-trivial problem. In model-based diagnosis, physical-based models are used to create residuals that isolate faults by mapping model equations to faulty system components. Developing sufficiently accurate physical-based models can be a time-consuming process. Hybrid modeling methods combining physical-based methods and machine learning is one solution to design data-driven residuals for fault isolation. In this work, a set of neural network-based residuals are designed by incorporating physical insights about the system behavior in the residual model structure. The residuals are trained using only fault-free data and a simulation case study shows that they can be used to perform fault isolation and localization of unknown faults in the system.


title: Lehto--Virtanen-type and big Picard-type theorems for Berkovich analytic spaces

abstract: In non-archimedean setting, we establish a Lehto--Virtanen-type theorem for a morphism from the punctured Berkovich closed unit disk $\overline{\mathsf{D}}\setminus\{0\}$ in the Berkovich affine line to the Berkovich projective line $\mathsf{P}^1$ having an isolated essential singularity at the origin, and then establish a big Picard-type theorem for such an open subset $\Omega$ in the Berkovich projective space $\mathsf{P}^N$ of any dimension $N$ that the family of all morphisms from $\overline{\mathsf{D}}\setminus\{0\}$ to $\Omega$ is normal in a non-archimedean Montel's sense. As an application of the latter theorem, we see a big Brody-type hyperbolicity of the Berkovich harmonic Fatou set of an endomorphism of $\mathsf{P}^N$ of degree $>1$.


title: Entropy-driven impurity-induced nematic-isotropic transition of liquid crystals

abstract: Phase behavior of liquid crystals is of long-standing interest due to numerous applications, with one of the key issues being how the presence of impurities affects the liquid crystalline order. Here we study the orientational order of 4-cyano-4$^{'}$-pentylbiphenyl (5CB) and 4-cyano-4$^{'}$-hexylbiphenyl (6CB) nematic liquid crystals in the presence of varying concentrations of water and n-hexane molecules serving as impurities, by carrying out both fully atomistic simulations and experiments. Our results reveal that mixing of the impurities (in case of hexane) with the host liquid crystals causes a nematic-to-isotropic phase transition with hexane concentration as the control parameter while demixing (in case of water) results in only weak impurity-induced perturbations to the nematic liquid-crystalline order. We develop a coarse-grained model illustrating the general nature and entropic origin of the mixing-induced phase transition.


title: Asymptotically Weyl-Invariant Gravity

abstract: We propose a novel theory of gravity that by construction is renormalizable, evades Ostragadsky's no-go theorem, is locally scale-invariant in the high-energy limit, and equivalent to general relativity in the low-energy limit. The theory is defined by a pure f(R) = R^n action in the Palatini formalism, where the dimensionless exponent n runs from a value of two in the high-energy limit to one in the low-energy limit. We show that the proposed model contains no obvious cosmological curvature singularities. The viability of the proposed model is qualitatively assessed using several key criteria.


title: The linear stability of Reissner-Nordstr\"om spacetime: the full subextremal range

abstract: We prove the linear stability of subextremal Reissner-Nordstr\"om spacetimes as solutions to the Einstein-Maxwell equation. We make use of a novel representation of gauge-invariant quantities which satisfy a symmetric system of coupled wave equations. This system is composed of two of the three equations separately derived in previous works, where the estimates required arbitrary smallness of the charge. Here, the estimates are obtained by defining a combined energy-momentum tensor for the system in terms of the symmetric structure of the right hand sides of the equations. We obtain boundedness of the energy, Morawetz estimates and decay for the full subextremal range |Q|<M, completely in physical space. Such decay estimates, together with the estimates for the gauge-dependent quantities of the perturbations previously obtained, settle the problem of linear stability to gravitational and electromagnetic perturbations of Reissner-Nordstr\"om solution in the full subextremal range |Q|< M.


title: Forecasts on the speed of gravitational waves at high $z$

abstract: The observation of GW170817 binary neutron star (BNS) merger event has imposed strong bounds on the speed of gravitational waves (GWs) locally, inferring that the speed of GWs propagation is equal to the speed of light. Current GW detectors in operation will not be able to observe BNS merger to long cosmological distance, where possible cosmological corrections on the cosmic expansion history are expected to play an important role, specially for investigating possible deviations from general relativity. Future GW detectors designer projects will be able to detect many coalescences of BNS at high $z$, such as the third generation of the ground GW detector called Einstein Telescope (ET) and the space-based detector deci-hertz interferometer gravitational wave observatory (DECIGO). In this paper, we relax the condition $c_T/c = 1$ to investigate modified GW propagation where the speed of GWs propagation is not necessarily equal to the speed of light. Also, we consider the possibility for the running of the Planck mass corrections on modified GW propagation. We parametrize both corrections in terms of an effective GW luminosity distance and we perform a forecast analysis using standard siren events from BNS mergers, within the sensitivity predicted for the ET and DECIGO. We find at high $z$ very strong forecast bounds on the running of the Planck mass, namely $\mathcal{O}(10^{-1})$ and $\mathcal{O}(10^{-2})$ from ET and DECIGO, respectively. Possible anomalies on GW propagation are bound to $|c_T/c - 1| \leq 10^{-2} \,\,\, (10^{-2})$ from ET (DECIGO), respectively. We finally discuss the consequences of our results on modified gravity phenomenology.


title: Hydrogen as a Renewable Energy Carrier for Commercial Aircraft

abstract: This report starts with the state of the art of electrolysis and liquefaction technology, and then presents a techno-economic case study of converting LAX airport to hydrogen fuel. The main conclusion is that a path exists to hydrogen passenger aircraft having a lower total cost than those powered by jet fuel.


title: Velocity-dependent self-interacting dark matter from thermal freeze-out and tests in direct detections

abstract: A small fraction of millicharged dark matter (DM) is considered in the literature to give an interpretation about the enhanced 21-cm absorption at the cosmic dawn. Here we focus on the case that the main component of DM is self-interacting dark matter (SIDM), motivated by the small scale problems. For self interactions of SIDM being compatible from dwarf to cluster scales, velocity-dependent self interactions mediated by a light scalar $\phi$ is considered. To fermionic SIDM $\Psi$, the main annihilation mode $\Psi \bar{\Psi} \to \phi \phi$ is a $p -$wave process. The thermal transition of SIDM $\rightleftarrows \phi \rightleftarrows$ standard model (SM) particles in the early universe sets a lower bound on couplings of $\phi$ to SM particles, which has been excluded by DM direct detections, and here we consider SIDM in the thermal equilibrium via millicharged DM. For $m_\phi >$ twice millicharged DM mass, $\phi$ could decay quickly and avoid excess energy injection to the big bang nucleosynthesis. Thus, the $\phi -$SM particle couplings could be very tiny and evade DM direct detections. The picture of weakly interacting massive particle (WIMP)-nucleus scattering with contact interactions fails for SIDM-nucleus scattering with a light mediator, and a method is explored in this paper, with which a WIMP search result can be converted into the hunt for SIDM in direct detections.


title: On the Weights of General MDS Codes

abstract: The weight spectra of MDS codes of length $ n $ and dimension $ k $ over the arbitrary alphabets are studied. For all $ q $-ary MDS codes of dimension $ k $ containing the zero codeword, it is shown that all $ k $ weights from $ n $ to $ n-k+1 $ are realized. The remaining case $ n=q+k-1 $ is also determined. Additionally, we prove that all binary MDS codes are equivalent to linear MDS codes. The proofs are combinatorial, and self contained.


title: Theory of Quantum Work in Metallic Grains

abstract: We generalize Anderson's orthogonality determinant formula to describe the statistics of work performed on generic disordered, non-interacting fermionic nanograins during quantum quenches. The energy absorbed increases linearly with time, while its variance exhibits a superdiffusive behavior due to Pauli's exclusion principle. The probability of adiabatic evolution decays as a stretched exponential. In slowly driven systems, work statistics exhibits universal features, and can be understood in terms of fermion diffusion in energy space, generated by Landau-Zener transitions. This diffusion is very well captured by a Markovian symmetrical exclusion process, with the diffusion constant identified as the energy absorption rate. The energy absorption rate shows an anomalous frequency dependence at small energies, reflecting the symmetry class of the underlying Hamiltonian. Our predictions can be experimentally verified by calorimetric measurements performed on nanoscale circuits.


title: Trajectory optimization for a class of robots belonging to Constrained Collaborative Mobile Agents (CCMA) family

abstract: We present a novel class of robots belonging to Constrained Collaborative Mobile Agents (CCMA) family which consists of ground mobile bases with non-holonomic constraints. Moreover, these mobile robots are constrained by closed-loop kinematic chains consisting of revolute joints which can be either passive or actuated. We also describe a novel trajectory optimization method which is general with respect to number of mobile robots, topology of the closed-loop kinematic chains and placement of the actuators at the revolute joints. We also extend the standalone trajectory optimization method to optimize concurrently the design parameters and the control policy. We describe various CCMA system examples, in simulation, differing in design, topology, number of mobile robots and actuation space. The simulation results for standalone trajectory optimization with fixed design parameters is presented for CCMA system examples. We also show how this method can be used for tasks other than end-effector positioning such as internal collision avoidance and external obstacle avoidance. The concurrent design and control policy optimization is demonstrated, in simulations, to increase the CCMA system workspace and manipulation capabilities. Finally, the trajectory optimization method is validated in experiments through two 4-DOF prototypes consisting of 3 tracked mobile bases.


title: Predictions of Cosmic Microwave Background Foregrounds Dust Polarization Using Velocity Gradients

abstract: The observations of fluctuations in the cosmic microwave background provide information about primordial inhomogeneities in the universe. However, the B-mode polarization of the inflationary gravitational wave is contaminated by the Galactic foreground polarized radiation arising from dust aligned by interstellar magnetic fields. To trace magnetic fields we use the Velocity Gradient Technique, which employs a modern understanding of the nature of magneto-hydrodynamics turbulent motions. In this paper, we combine the VGT with the Principal Component Analysis (PCA) to improve the accuracy of magnetic field tracing. We apply the VGT-PCA to the high-resolution neutral hydrogen data from the GALFA-HI survey to predict the polarization of dust. We report that the predicted directions of dust polarization provide good correspondence with those reported by Planck 353GHz, the alignment measure between the two AM=0.79. We show that our results statistically agree with the Planck polarization in terms of magnetic field tracing. We find that the variation of dust emission efficiency across the sky is small. Using our maps of predicted polarization we calculate the ratio of the E and B modes and show that BB/EE=0.53, which is similar to the result from Plank polarization.


title: Coset Posets of Infinite Groups

abstract: We consider the coset poset associated with the families of proper subgroups, proper subgroups of finite index, and proper normal subgroups of finite index. We investigate under which conditions those coset posets have contractible geometric realizations.


title: Disentangling Interpretable Generative Parameters of Random and Real-World Graphs

abstract: While a wide range of interpretable generative procedures for graphs exist, matching observed graph topologies with such procedures and choices for its parameters remains an open problem. Devising generative models that closely reproduce real-world graphs requires domain knowledge and time-consuming simulation. While existing deep learning approaches rely on less manual modelling, they offer little interpretability. This work approaches graph generation (decoding) as the inverse of graph compression (encoding). We show that in a disentanglement-focused deep autoencoding framework, specifically Beta-Variational Autoencoders (Beta-VAE), choices of generative procedures and their parameters arise naturally in the latent space. Our model is capable of learning disentangled, interpretable latent variables that represent the generative parameters of procedurally generated random graphs and real-world graphs. The degree of disentanglement is quantitatively measured using the Mutual Information Gap (MIG). When training our Beta-VAE model on ER random graphs, its latent variables have a near one-to-one mapping to the ER random graph parameters n and p. We deploy the model to analyse the correlation between graph topology and node attributes measuring their mutual dependence without handpicking topological properties.


title: Deep Learning for Predicting Dynamic Uncertain Opinions in Network Data

abstract: Subjective Logic (SL) is one of well-known belief models that can explicitly deal with uncertain opinions and infer unknown opinions based on a rich set of operators of fusing multiple opinions. Due to high simplicity and applicability, SL has been substantially applied in a variety of decision making in the area of cybersecurity, opinion models, trust models, and/or social network analysis. However, SL and its variants have exposed limitations in predicting uncertain opinions in real-world dynamic network data mainly in three-fold: (1) a lack of scalability to deal with a large-scale network; (2) limited capability to handle heterogeneous topological and temporal dependencies among node-level opinions; and (3) a high sensitivity with conflicting evidence that may generate counterintuitive opinions derived from the evidence. In this work, we proposed a novel deep learning (DL)-based dynamic opinion inference model while node-level opinions are still formalized based on SL meaning that an opinion has a dimension of uncertainty in addition to belief and disbelief in a binomial opinion (i.e., agree or disagree). The proposed DL-based dynamic opinion inference model overcomes the above three limitations by integrating the following techniques: (1) state-of-the-art DL techniques, such as the Graph Convolutional Network (GCN) and the Gated Recurrent Units (GRU) for modeling the topological and temporal heterogeneous dependency information of a given dynamic network; (2) modeling conflicting opinions based on robust statistics; and (3) a highly scalable inference algorithm to predict dynamic, uncertain opinions in a linear computation time. We validated the outperformance of our proposed DL-based algorithm (i.e., GCN-GRU-opinion model) via extensive comparative performance analysis based on four real-world datasets.


title: On categories of o-minimal structures

abstract: Our aim in this paper is to look at some transfer results in model theory (mainly in the context of o-minimal structures) from the category theory viewpoint.


title: Kinematical Gravitational Charge Algebra

abstract: When formulated in terms of connection and coframes, and in the time gauge, the phase space of general relativity consists of a pair of conjugate fields: the flux 2-form and the Ashtekar connection. On this phase-space, one has to impose the Gauss constraints, the vector, and scalar Hamiltonian constraints. These are respectively generating local SU(2) gauge transformations, spatial diffeomorphisms, and time diffeomorphisms. We write the Gauss and space diffeomorphism constraints as conservation laws for a set of boundary charges, representing spin and momenta, respectively. We prove that these kinematical charges generate a local Poincar\'e ISU(2) symmetry algebra. This gives strong support to the recent proposal of Poincar\'e charge networks as a new realm for discretized general relativity [Classical Quantum Gravity 36, 195014 (2019)].


title: Molecular and DNA Artificial Neural Networks via Fractional Coding

abstract: This paper considers implementation of artificial neural networks (ANNs) using molecular computing and DNA based on fractional coding. Prior work had addressed molecular two-layer ANNs with binary inputs and arbitrary weights. In prior work using fractional coding, a simple molecular perceptron that computes sigmoid of scaled weighted sum of the inputs was presented where the inputs and the weights lie between [-1, 1]. Even for computing the perceptron, the prior approach suffers from two major limitations. First, it cannot compute the sigmoid of the weighted sum, but only the sigmoid of the scaled weighted sum. Second, many machine learning applications require the coefficients to be arbitrarily positive and negative numbers that are not bounded between [-1, 1]; such numbers cannot be handled by the prior perceptron using fractional coding. This paper makes four contributions. First molecular perceptrons that can handle arbitrary weights and can compute sigmoid of the weighted sums are presented. Thus, these molecular perceptrons are ideal for regression applications and multi-layer ANNs. A new molecular divider is introduced and is used to compute sigmoid(ax) where a > 1. Second, based on fractional coding, a molecular artificial neural network (ANN) with one hidden layer is presented. Third, a trained ANN classifier with one hidden layer from seizure prediction application from electroencephalogram is mapped to molecular reactions and DNA and their performances are presented. Fourth, molecular activation functions for rectified linear unit (ReLU) and softmax are also presented.


title: Response of Vertical Velocities in Extratropical Precipitation Extremes to Climate Change

abstract: Precipitation extremes intensify in most regions in climate-model projections. Changes in vertical velocities contribute to the changes in intensity of precipitation extremes but remain poorly understood. Here, we find that mid-tropospheric vertical velocities in extratropical precipitation extremes strengthen overall in simulations of 21st-century climate change. For each extreme event, we solve the quasi-geostrophic omega equation to decompose this strengthening into different physical contributions. We first consider a dry decomposition in which latent heating is treated as an external forcing of upward motion. Much of the positive contribution to upward motion from increased latent heating is offset by negative contributions from increases in dry static stability and changes in the horizontal length scale of vertical velocities. However, taking changes in latent heating as given is a limitation when the aim is to understand changes in precipitation, since latent heating and precipitation are closely linked. Therefore, we also perform a moist decomposition of the changes in vertical velocities in which latent heating is represented through a moist static stability. In the moist decomposition, changes in moist static stability play a key role and contributions from other factors such as changes in the depth of the upward motion increase in importance. While both dry and moist decompositions are self-consistent, the moist dynamical perspective has greater potential to give insights into the causes of the dynamical contributions to changes in precipitation extremes in different regions.


title: Reachability and Shortest Paths in the Broadcast CONGEST Model

abstract: In this paper we study the time complexity of the single-source reachability problem and the single-source shortest path problem for directed unweighted graphs in the Broadcast CONGEST model. We focus on the case where the diameter $D$ of the underlying network is constant. We show that for the case where $D = 1$ there is, quite surprisingly, a very simple algorithm that solves the reachability problem in $1$(!) round. In contrast, for networks with $D = 2$, we show that any distributed algorithm (possibly randomized) for this problem requires $\Omega(\sqrt{n/ \log{n}}\,)$ rounds. Our results therefore completely resolve (up to a small polylogarithmic factor) the complexity of the single-source reachability problem for a wide range of diameters. Furthermore, we show that when $D = 1$, it is even possible to get a $3$-approximation for the all-pairs shortest path problem (for directed unweighted graphs) in just $2$ rounds. We also prove a stronger lower bound of $\Omega(\sqrt{n}\,)$ for the single-source shortest path problem for unweighted directed graphs that holds even when the diameter of the underlying network is $2$. As far as we know this is the first lower bound that achieves $\Omega(\sqrt{n}\,)$ for this problem.


title: "Bring Your Own Greedy"+Max: Near-Optimal $1/2$-Approximations for Submodular Knapsack

abstract: The problem of selecting a small-size representative summary of a large dataset is a cornerstone of machine learning, optimization and data science. Motivated by applications to recommendation systems and other scenarios with query-limited access to vast amounts of data, we propose a new rigorous algorithmic framework for a standard formulation of this problem as a submodular maximization subject to a linear (knapsack) constraint. Our framework is based on augmenting all partial Greedy solutions with the best additional item. It can be instantiated with negligible overhead in any model of computation, which allows the classic \greedy algorithm and its variants to be implemented. We give such instantiations in the offline (Greedy+Max), multi-pass streaming (Sieve+Max) and distributed (Distributed+Max) settings. Our algorithms give ($1/2-\epsilon$)-approximation with most other key parameters of interest being near-optimal. Our analysis is based on a new set of first-order linear differential inequalities and their robust approximate versions. Experiments on typical datasets (movie recommendations, influence maximization) confirm scalability and high quality of solutions obtained via our framework. Instance-specific approximations are typically in the 0.6-0.7 range and frequently beat even the $(1-1/e) \approx 0.63$ worst-case barrier for polynomial-time algorithms.


title: IoT or NoT: Identifying IoT Devices in a ShortTime Scale

abstract: In recent years the number of IoT devices in home networks has increased dramatically. Whenever a new device connects to the network, it must be quickly managed and secured using the relevant security mechanism or QoS policy. Thus a key challenge is to distinguish between IoT and NoT devices in a matter of minutes. Unfortunately, there is no clear indication of whether a device in a network is an IoT. In this paper, we propose different classifiers that identify a device as IoT or non-IoT, in a short time scale, and with high accuracy. Our classifiers were constructed using machine learning techniques on a seen (training) dataset and were tested on an unseen (test) dataset. They successfully classified devices that were not in the seen dataset with accuracy above 95%. The first classifier is a logistic regression classifier based on traffic features. The second classifier is based on features we retrieve from DHCP packets. Finally, we present a unified classifier that leverages the advantages of the other two classifiers. We focus on the home-network environment, but our classifiers are also applicable to enterprise networks.


title: A Variational Finite Element Discretization of Compressible Flow

abstract: We present a finite element variational integrator for compressible flows. The numerical scheme is derived by discretizing, in a structure preserving way, the Lie group formulation of fluid dynamics on diffeomorphism groups and the associated variational principles. Given a triangulation on the fluid domain, the discrete group of diffeomorphisms is defined as a certain subgroup of the group of linear isomorphisms of a finite element space of functions. In this setting, discrete vector fields correspond to a certain subspace of the Lie algebra of this group. This subspace is shown to be isomorphic to a Raviart-Thomas finite element space. The resulting finite element discretization corresponds to a weak form of the compressible fluid equation that doesn't seem to have been used in the finite element literature. It extends previous work done on incompressible flows and at the lowest order on compressible flows. We illustrate the conservation properties of the scheme with some numerical simulations.


title: Robustness of cluster states and surface code states against random local fields

abstract: In ideal quantum circuits, qubits are tacitly assumed to be uniformly fabricated and operated by prescribed signals. In reality, however, we must cope with different control signals to adjust individual qubits, which requires a large overhead of control circuits. Here, we theoretically investigate how random local fields affect cluster states and surface code states which constitute the key highly entangled states in quantum computation. We find similar behavior of temporal degradation of the fidelity for both cluster states and surface code states for the number of qubits up to ten. We find that the effect of local field fluctuations is greatly mitigated if the magnitude of fluctuations can be suppressed below 10 % of the energy gap $\Delta$ for both cluster states and surface code states. If the magnitude of fluctuations exceeds $\Delta/2$, the fidelity for both states deteriorates dramatically. A simple estimation based on the average fidelity up to $t\sim 2\hbar/\Delta$ shows that the maximum number of qubits that can be corrected with the 1 % error threshold is less than 31 for surface code states and 27 for cluster states. This means that the error correction should be carried out during a time shorter than $2\hbar/\Delta$.


title: On the probability distribution of the local times of diagonally operator-self-similar Gaussian fields with stationary increments

abstract: In this paper we study the local times of vector-valued Gaussian fields that are `diagonally operator-self-similar' and whose increments are stationary. Denoting the local time of such a Gaussian field around the spatial origin and over the temporal unit hypercube by $Z$, we show that there exists $\lambda\in(0,1)$ such that under some quite weak conditions, $\lim_{n\rightarrow +\infty}\frac{\sqrt[n]{\mathbb{E}(Z^n)}}{n^\lambda}$ and $\lim_{x\rightarrow +\infty}\frac{-\log \mathbb{P}(Z>x)}{x^{\frac{1}{\lambda}}}$ both exist and are strictly positive (possibly $+\infty$). Moreover, we show that if the underlying Gaussian field is `strongly locally nondeterministic', the above limits will be finite as well. These results are then applied to establish similar statements for the intersection local times of diagonally operator-self-similar Gaussian fields with stationary increments.


title: Interventional Experiment Design for Causal Structure Learning

abstract: It is known that from purely observational data, a causal DAG is identifiable only up to its Markov equivalence class, and for many ground truth DAGs, the direction of a large portion of the edges will be remained unidentified. The golden standard for learning the causal DAG beyond Markov equivalence is to perform a sequence of interventions in the system and use the data gathered from the interventional distributions. We consider a setup in which given a budget $k$, we design $k$ interventions non-adaptively. We cast the problem of finding the best intervention target set as an optimization problem which aims to maximize the number of edges whose directions are identified due to the performed interventions. First, we consider the case that the underlying causal structure is a tree. For this case, we propose an efficient exact algorithm for the worst-case gain setup, as well as an approximate algorithm for the average gain setup. We then show that the proposed approach for the average gain setup can be extended to the case of general causal structures. In this case, besides the design of interventions, calculating the objective function is also challenging. We propose an efficient exact calculator as well as two estimators for this task. We evaluate the proposed methods using synthetic as well as real data.


title: What is the Largest Sparsity Pattern that Can Be Recovered by 1-Norm Minimization?

abstract: Much of the existing literature in sparse recovery is concerned with the following question: given a sparsity pattern and a corresponding regularizer, derive conditions on the dictionary under which exact recovery is possible. In this paper, we study the opposite question: given a dictionary and the 1-norm regularizer, find the largest sparsity pattern that can be recovered. We show that such a pattern is described by a mathematical object called a "maximum abstract simplicial complex", and provide two different characterizations of this object: one based on extreme points and the other based on vectors of minimal support. In addition, we show how this new framework is useful in the study of sparse recovery problems when the dictionary takes the form of a graph incidence matrix or a partial discrete Fourier transform. In case of incidence matrices, we show that the largest sparsity pattern that can be recovered is determined by the set of simple cycles of the graph. As a byproduct, we show that standard sparse recovery can be certified in polynomial time, although this is known to be NP-hard for general matrices. In the case of the partial discrete Fourier transform, our characterization of the largest sparsity pattern that can be recovered requires the unknown signal to be real and its dimension to be a prime number.


title: Model Fusion via Optimal Transport

abstract: Combining different models is a widely used paradigm in machine learning applications. While the most common approach is to form an ensemble of models and average their individual predictions, this approach is often rendered infeasible by given resource constraints in terms of memory and computation, which grow linearly with the number of models. We present a layer-wise model fusion algorithm for neural networks that utilizes optimal transport to (soft-) align neurons across the models before averaging their associated parameters. We show that this can successfully yield "one-shot" knowledge transfer (i.e, without requiring any retraining) between neural networks trained on heterogeneous non-i.i.d. data. In both i.i.d. and non-i.i.d. settings , we illustrate that our approach significantly outperforms vanilla averaging, as well as how it can serve as an efficient replacement for the ensemble with moderate fine-tuning, for standard convolutional networks (like VGG11), residual networks (like ResNet18), and multi-layer perceptrons on CIFAR10, CIFAR100, and MNIST. Finally, our approach also provides a principled way to combine the parameters of neural networks with different widths, and we explore its application for model compression. The code is available at the following link, https://github.com/sidak/otfusion.


title: Thompson Sampling in Non-Episodic Restless Bandits

abstract: Restless bandit problems assume time-varying reward distributions of the arms, which adds flexibility to the model but makes the analysis more challenging. We study learning algorithms over the unknown reward distributions and prove a sub-linear, $O(\sqrt{T}\log T)$, regret bound for a variant of Thompson sampling. Our analysis applies in the infinite time horizon setting, resolving the open question raised by Jung and Tewari (2019) whose analysis is limited to the episodic case. We adopt their policy mapping framework, which allows our algorithm to be efficient and simultaneously keeps the regret meaningful. Our algorithm adapts the TSDE algorithm of Ouyang et al. (2017) in a non-trivial manner to account for the special structure of restless bandits. We test our algorithm on a simulated dynamic channel access problem with several policy mappings, and the empirical regrets agree with the theoretical bound regardless of the choice of the policy mapping.


title: The supermoduli space of genus zero SUSY curves with Ramond punctures

abstract: We give an explicit quotient construction of the supermoduli space $\mathfrak{M}_{0, n_R}$ of genus zero super Riemann surfaces with $n_R \ge 4$ Ramond punctures and prove it is a Deligne-Mumford superstack. We also make an explicit quotient construction of the moduli space of genus zero supercurves without a SUSY structure and thereby prove it is an algebraic superstack.


title: K\"unneth Formulae in Persistent Homology

abstract: The classical K\"{u}nneth formula in algebraic topology describes the homology of a product space in terms of that of its factors. In this paper, we prove K\"{u}nneth-type theorems for the persistent homology of the categorical and tensor product of filtered spaces. That is, we describe the persistent homology of these product filtrations in terms of that of the filtered components. In addition to comparing the two products, we present two applications in the setting of Vietoris-Rips complexes: one towards more efficient algorithms for product metric spaces with the maximum metric, and the other recovering persistent homology calculations on the $n$-torus.


title: How are attributes expressed in face DCNNs?

abstract: As deep networks become increasingly accurate at recognizing faces, it is vital to understand how these networks process faces. While these networks are solely trained to recognize identities, they also contain face related information such as sex, age, and pose of the face. The networks are not trained to learn these attributes. We introduce expressivity as a measure of how much a feature vector informs us about an attribute, where a feature vector can be from internal or final layers of a network. Expressivity is computed by a second neural network whose inputs are features and attributes. The output of the second neural network approximates the mutual information between feature vectors and an attribute. We investigate the expressivity for two different deep convolutional neural network (DCNN) architectures: a Resnet-101 and an Inception Resnet v2. In the final fully connected layer of the networks, we found the order of expressivity for facial attributes to be Age > Sex > Yaw. Additionally, we studied the changes in the encoding of facial attributes over training iterations. We found that as training progresses, expressivities of yaw, sex, and age decrease. Our technique can be a tool for investigating the sources of bias in a network and a step towards explaining the network's identity decisions.


title: Universal Basic Income: The Last Bullet in the Darkness

abstract: Universal Basic Income (UBI) has recently been gaining traction. Arguments exist on both sides in favor of and against it. Like any other financial tool, UBI can be useful if used with discretion. This paper seeks to clarify how UBI affects the economy, including how it can be beneficial. The key point is to regulate the rate of UBI based on the inflation rate. This should be done by an independent institution from the executive branch of the government. If implemented correctly, UBI can add a powerful tool to the Federal Reserve toolkit. UBI can be used to reintroduce inflation to the countries which suffer long-lasting deflationary environment. UBI has the potential to decrease the wealth disparity, decrease the national debt, increase productivity, and increase comparative advantage of the economy. UBI also can substitute the current welfare systems because of its transparency and efficiency. This article focuses more on the United States, but similar ideas can be implemented in other developed nations.


title: A construction of complex analytic elliptic cohomology from double free loop spaces

abstract: We construct a complex analytic version of an equivariant cohomology theory which appeared in a recent paper of Rezk, and which is roughly modeled on the Borel-equivariant cohomology of the double free loop space. The construction is defined on finite, torus-equivariant CW complexes and takes values in coherent holomorphic sheaves over the moduli stack of complex elliptic curves. Our methods involve an inverse limit construction over all finite dimensional subcomplexes of the double free loop space, following an analogous construction of Kitchloo for single free loop spaces. We show that, for any given complex elliptic curve $\mathcal{C}$, the fiber of our construction over $\mathcal{C}$ is isomorphic to Grojnowski's equivariant elliptic cohomology theory associated to $\mathcal{C}$.


title: Lattice Boltzmann simulations of droplet breakup in confined and time-dependent flows

abstract: We study droplet dynamics and breakup in generic time-dependent flows via a multicomponent Lattice Boltzmann algorithm, with emphasis on flow start up conditions. We first study droplet breakup in a confined oscillatory shear flow via two different protocols. In one set up, we start from an initially spherical droplet and turn on the flow abruptly ("shock method"); in the other protocol, we start from an initially spherical droplet as well, but we progressively increase the amplitude of the flow, by allowing the droplet to relax to the steady state for each increase in amplitude, before increasing the flow amplitude again ("relaxation method"). The two protocols are shown to produce substantially different breakup scenarios. The mismatch between these two protocols is also studied for variations in the flow topology, the degree of confinement and the inertia of the fluid. All results point to the fact that under extreme conditions of confinement the relaxation protocols can drive the droplets into metastable states, which break only for very intense flow amplitudes, but their stability is prone to external perturbations, such as an oscillatory driving force.


title: Three-phase Golay sequence and array triads

abstract: 3-phase Golay sequence and array triads are a natural generalisation of 2-phase Golay sequence pairs, yet their study has until now been largely neglected. We present exhaustive search results for 3-phase Golay sequence triads for all lengths up to 24, showing that the existence pattern is much richer than that for 2-phase Golay sequence pairs. We give an elementary proof that there is no 3-phase Golay sequence triad whose length is congruent to 4 modulo 6. We give two construction methods for 3-phase Golay array triads, and show how to project 3-phase Golay array triads to lower-dimensional Golay triads. In this way we explain much of the existence pattern found by exhaustive search.


title: X-rays from Stars and Planetary Systems

abstract: The Chandra X-ray Observatory has completed a remarkable twenty years in orbit. A large part of the science program of Chandra during this time has involved the study of stars and their planetary systems. This primer aims to give the reader a taste of the enormous range of stellar and planetary astrophysics that Chandra has enabled. Beginning with a tour of the X-ray solar system zoo, including the stunning pulsating X-ray aurorae of Jupiter, we then move on to the hot million-degree outer atmospheres of stars like our own Sun, whose X-ray emission is driven by an internal magnetic dynamo. The same emission processes are also vigorously present in the youngest stars, and we highlight some Chandra observations and results on nascent stellar and planetary systems. Chandra surveys and high-resolution spectroscopy of massive stars have provided a new window on the means by which they scavenge X-ray emission from their radiatively-driven winds, sometimes modulating this output by strong underlying stellar magnetic fields. We touch upon the evanescent X-radiation from intermediate-mass stars before arriving at the inevitable evolutionary endpoints of all but massive stars, first in energized X-ray emitting planetary nebulae, then in the slowly cooling, soft-X-ray emitting photospheres of white dwarfs. We conclude with white dwarfs in close binary systems, rejuvenated by interaction with a companion and where accretion gives play to a new range of energetic behavior even more spectacular and cataclysmic than the coruscant astrophysical road down which they have travelled.


title: Hadamard Renormalization of a 2-Dimensional Dirac Field

abstract: The Hadamard renormalization procedure is applied to a free, massive Dirac field $\psi$ on a 2 dimensional Lorentzian spacetime. This yields the state-independent divergent terms in the Hadamard bispinor $G^{(1)}(x, x') = \frac{1}{2} \left\langle \left[ \bar{\psi}(x'), \psi(x) \right] \right\rangle$ as $x$ and $x'$ are brought together along the unique geodesic connecting them. Subtracting these divergent terms within the limit assigns $G^{(1)}(x, x')$, and thus any operator expressed in terms of it, a finite value at the coincident point $x' = x$. In this limit, one obtains a quadratic operator instead of a bispinor. The procedure is thus used to assign finite values to various quadratic operators, including the stress-energy tensor. Results are presented covariantly, in a conformally-flat coordinate chart at purely spatial separations, and in the Minkowski metric. These terms can be directly subtracted from combinations of $G^{(1)}(x, x')$ - themselves obtained, for example, from a numerical simulation - to obtain finite expectation values defined in the continuum.


title: Extracting Incentives from Black-Box Decisions

abstract: An algorithmic decision-maker incentivizes people to act in certain ways to receive better decisions. These incentives can dramatically influence subjects' behaviors and lives, and it is important that both decision-makers and decision-recipients have clarity on which actions are incentivized by the chosen model. While for linear functions, the changes a subject is incentivized to make may be clear, we prove that for many non-linear functions (e.g. neural networks, random forests), classical methods for interpreting the behavior of models (e.g. input gradients) provide poor advice to individuals on which actions they should take. In this work, we propose a mathematical framework for understanding algorithmic incentives as the challenge of solving a Markov Decision Process, where the state includes the set of input features, and the reward is a function of the model's output. We can then leverage the many toolkits for solving MDPs (e.g. tree-based planning, reinforcement learning) to identify the optimal actions each individual is incentivized to take to improve their decision under a given model. We demonstrate the utility of our method by estimating the maximally-incentivized actions in two real-world settings: a recidivism risk predictor we train using ProPublica's COMPAS dataset, and an online credit scoring tool published by the Fair Isaac Corporation (FICO).


title: Dispersionless Integrable Hierarchy via Kodaira-Spencer Gravity

abstract: We explain how dispersionless integrable hierarchy in 2d topological field theory arises from the Kodaira-Spencer gravity (BCOV theory). The infinitely many commuting Hamiltonians are given by the current observables associated to the infinite abelian symmetries of the Kodaira-Spencer gravity. We describe a BV framework of effective field theories that leads to the B-model interpretation of dispersionless integrable hierarchy.


title: Constraint on the fifth force through perihelion precession of planets

abstract: The equivalence principle is important in fundamental physics. The fifth force, as a describing formalism of the equivalence principle, may indicate the property of an unknown theory. Dark matter is one of the most mysterious objects in the current natural science. It is interesting to constrain the fifth force of dark matter. We propose a new method to use perihelion precession of planets to constrain the long-range fifth force of dark matter. Due to the high accuracy of perihelion precession observation, and the large difference of matter composition between the Sun and planets, we get one of the strongest constraints on the fifth force of dark matter. In the near future, the BepiColombo mission will be capable to improve the test by another factor of ten.


title: Face flips in origami tessellations

abstract: Given a flat-foldable origami crease pattern $G=(V,E)$ (a straight-line drawing of a planar graph on a region of the plane) with a mountain-valley (MV) assignment $\mu:E\to\{-1,1\}$ indicating which creases in $E$ bend convexly (mountain) or concavely (valley), we may \emph{flip} a face $F$ of $G$ to create a new MV assignment $\mu_F$ which equals $\mu$ except for all creases $e$ bordering $F$, where we have $\mu_F(e)=-\mu(e)$. In this paper we explore the configuration space of face flips for a variety of crease patterns $G$ that are tilings of the plane, proving examples where $\mu_F$ results in a MV assignment that is either never, sometimes, or always flat-foldable for various choices of $F$. We also consider the problem of finding, given two foldable MV assignments $\mu_1$ and $\mu_2$ of a given crease pattern $G$, a minimal sequence of face flips to turn $\mu_1$ into $\mu_2$. We find polynomial-time algorithms for this in the cases where $G$ is either a square grid or the Miura-ori, and show that this problem is NP-hard in the case where $G$ is the triangle lattice.


title: A systematically study of thermal width of heavy quarkonia in a finite temperature magnetized background from holography

abstract: By simulating the finite temperatures magnetized background in the RHIC and LHC energies, we systematically study the characteristics of thermal widths and potentials of heavy quarkonia. It is found that the magnetic field has less influence on the real potential, but has a significant influence on the imaginary potential, especially in the low deconfined temperature. Extracted from the effect of thermal worldsheet fluctuations about the classical configuration, the thermal width of $\Upsilon(1s)$ in the finite temperature magnetized background is investigated. It is found that at the low deconfined temperature the magnetic field can generate a significant thermal fluctuation of the thermal width of $\Upsilon(1s)$, but with the increase of temperature, the effect of magnetic field on the thermal width becomes less important, which means the effect of high temperature completely exceeds that of magnetic field and magnetic field become less important at high temperature. The thermal width decreases with the increasing rapidity at the finite temperature magnetized background. It is also observed that the effect of the magnetic field on the thermal width when dipole moving parallel to the magnetic field direction are larger than that moving perpendicular to the magnetic field direction, which implies that the magnetic field tends to enhance thermal fluctuation when dipole moving parallel to the direction of magnetic field. The thermal width of $\Upsilon(1S)$ hardly changes with the increasing temperature when dipole moving perpendicular to the magnetic field. But when dipole moving parallel to the magnetic field, the thermal width at low temperature is obviously larger than that at high temperature.


title: Model Predictive Tracking Control for Invariant Systems on Matrix Lie Groups via Stable Embedding into Euclidean Spaces

abstract: For controller design for systems on manifolds embedded in Euclidean space, it is convenient to utilize a theory that requires a single global coordinate system on the ambient Euclidean space rather than multiple local charts on the manifold or coordinate-free tools from differential geometry. In this article, we apply such a theory to design model predictive tracking controllers for systems whose dynamics evolve on manifolds and illustrate its efficacy with the fully actuated rigid body attitude control system.


title: Band-limited Features in the Primordial Power Spectrum Do Not Resolve the Hubble Tension

abstract: For a standard $\Lambda$CDM universe with a power-law primordial power spectrum, the discrepancy between early- and late-universe measurements of the Hubble constant continued to grow and recently reached $5.3\sigma$. During inflation, local features in the inflationary potential often lead to band-limited features in the primordial power spectrum, hence breaking the power-law assumption in the derivation of the Hubble tension. We investigate whether such inflationary "glitches" can ease the Hubble tension. The recently released \emph{Planck} temperature and polarization data and the 2019 SH0ES+H0LiCOW joint constraint on the Hubble constant are combined to drive a blind Daubechies wavelet signal search in the primordial power spectrum, up to a resolution $\Delta \ln k\sim 0.1$. We find no significant detection of any features beyond power law. With 64 more degrees of freedom injected in the primordial power spectrum, the Hubble tension persists at a $4.9\sigma$ level.


title: A Banach space mixed formulation for the unsteady Brinkman-Forchheimer equations

abstract: We propose and analyze a mixed formulation for the Brinkman-Forchheimer equations for unsteady flows. Our approach is based on the introduction of a pseudostress tensor related to the velocity gradient, leading to a mixed formulation where the pseudostress tensor and the velocity are the main unknowns of the system. We establish existence and uniqueness of a solution to the weak formulation in a Banach space setting, employing classical results on nonlinear monotone operators and a regularization technique. We then present well-posedness and error analysis for semidiscrete continuous-in-time and fully discrete finite element approximations on simplicial grids with spatial discretization based on the Raviart-Thomas spaces of degree $k$ for the pseudostress tensor and discontinuous piecewise polynomial elements of degree $k$ for the velocity and backward Euler time discretization. We provide several numerical results to confirm the theoretical rates of convergence and illustrate the performance and flexibility of the method for a range of model parameters.


title: Optic-Net: A Novel Convolutional Neural Network for Diagnosis of Retinal Diseases from Optical Tomography Images

abstract: Diagnosing different retinal diseases from Spectral Domain Optical Coherence Tomography (SD-OCT) images is a challenging task. Different automated approaches such as image processing, machine learning and deep learning algorithms have been used for early detection and diagnosis of retinal diseases. Unfortunately, these are prone to error and computational inefficiency, which requires further intervention from human experts. In this paper, we propose a novel convolution neural network architecture to successfully distinguish between different degeneration of retinal layers and their underlying causes. The proposed novel architecture outperforms other classification models while addressing the issue of gradient explosion. Our approach reaches near perfect accuracy of 99.8% and 100% for two separately available Retinal SD-OCT data-set respectively. Additionally, our architecture predicts retinal diseases in real time while outperforming human diagnosticians.


title: Dynamics of bacterial populations under the feast-famine cycles

abstract: Bacterial populations in natural conditions are expected to experience stochastic environmental fluctuations, and in addition, environments are affected by bacterial activities since they consume substrates and excrete various chemicals. We here study possible outcomes of population dynamics and evolution under the repeated cycle of substrate-rich conditions and starvation, called the "feast-famine cycle", by a simple stochastic model with the trade-off relationship between the growth rate and the growth yield or the death rate. In the model, the feast (substrate-rich) period is led by a stochastic substrate addition event, while the famine (starvation) period is evoked because bacteria use the supplied substrate. Under the repeated feast-famine cycle, the bacterial population tends to increase the growth rate, even though that tends to decrease the total population size due to the trade-off. Analysis of the model shows that the ratio between the growth rate and the death rate becomes the effective fitness of the population. Hence, the functional form of the trade-off between the growth and death rate determines if the bacterial population eventually goes extinct as an evolutionary consequence. We then show that the increase of the added substrate in the feast period can drive the extinction faster. Overall, the model sheds light on non-trivial possible outcomes under repeated feast-famine cycles.


title: Structure-preserving Interpolatory Model Reduction for Port-Hamiltonian Differential-Algebraic Systems

abstract: We examine interpolatory model reduction methods that are well-suited for treating large scale port-Hamiltonian differential-algebraic systems in a way that is able to preserve and indeed, take advantage of the underlying structural features of the system. We introduce approaches that incorporate regularization together with prudent selection of interpolation data. We focus on linear time-invariant systems and present a systematic treatment of a variety of model classes that include combinations of index-$1$ and index-$2$ systems, describing in particular how constraints may be represented in the transfer function and then preserved with interpolatory methods. We propose an algorithm to generate effective interpolation data and illustrate its effectiveness via two numerical examples.


title: On the Fibonacci $(p,r)$-cubes

abstract: In this paper, first it is shown that the "FSibonacci $(p,r)$-cube"(denoted as $I\Gamma_{n}^{(p,r)}$) studied in many papers, such as \cite{OZY}, \cite{K1}, \cite{OZ}, \cite{KR} and \cite{JZ}, is a new topological structure different from the original one (denoted as $O\Gamma_{n}^{(p,r)}$) presented by Egiazarian and Astola $\cite{EA}$. Then some topological properties of $I\Gamma_{n}^{(p,r)}$ and $O\Gamma_{n}^{(p,r)}$ are studied, including the recursive structure of them, the cubes $O\Gamma_{n}^{(p,r)}$ which are partial cubes and median graphs, some distance invariants of $I\Gamma_{n}^{(p,r)}$ and $O\Gamma_{n}^{(p,r)}$, and the maximum and minimum degree of these two types of cubes. Finally, several problems and conjectures on $I\Gamma_{n}^{(p,r)}$ and $O\Gamma_{n}^{(p,r)}$ are listed


title: Regression for Copula-linked Compound Distributions with Applications in Modeling Aggregate Insurance Claims

abstract: In actuarial research, a task of particular interest and importance is to predict the loss cost for individual risks so that informative decisions are made in various insurance operations such as underwriting, ratemaking, and capital management. The loss cost is typically viewed to follow a compound distribution where the summation of the severity variables is stopped by the frequency variable. A challenging issue in modeling such outcome is to accommodate the potential dependence between the number of claims and the size of each individual claim. In this article, we introduce a novel regression framework for compound distributions that uses a copula to accommodate the association between the frequency and the severity variables, and thus allows for arbitrary dependence between the two components. We further show that the new model is very flexible and is easily modified to account for incomplete data due to censoring or truncation. The flexibility of the proposed model is illustrated using both simulated and real data sets. In the analysis of granular claims data from property insurance, we find substantive negative relationship between the number and the size of insurance claims. In addition, we demonstrate that ignoring the frequency-severity association could lead to biased decision-making in insurance operations.


title: Computer-inspired concept for high-dimensional multipartite quantum gates

abstract: An open question in quantum optics is how to manipulate and control complex quantum states in an experimentally feasible way. Here we present concepts for transformations of high-dimensional multi-photonic quantum systems. The proposals rely on two new ideas: (I) a novel high-dimensional quantum non-demolition measurement, (II) the encoding and decoding of the entire quantum transformation in an ancillary state for sharing the necessary quantum information between the involved parties. Many solutions can readily be performed in laboratories around the world, and identify important pathways for experimental research in the near future. The concept has been found using the computer algorithm Melvin for designing computer-inspired quantum experiments. This demonstrates that computer algorithms can inspire new ideas in science, which is a widely unexplored potential.


title: An Image Segmentation Model Based on a Variational Formulation

abstract: Starting from a variational formulation, we present a model for image segmentation that employs both region statistics and edge information. This combination allows for improved flexibility, making the proposed model suitable to process a wider class of images than purely region-based and edge-based models. We perform several simulations with real images that attest to the versatility of the model. We also show another set of experiments on images with certain pathologies that suggest opportunities for improvement.


title: The space of tree-based phylogenetic networks

abstract: Phylogenetic networks are generalizations of phylogenetic trees that allow the representation of reticulation events such as horizontal gene transfer or hybridization, and can also represent uncertainty in inference. A subclass of these, tree-based phylogenetic networks, have been introduced to capture the extent to which reticulate evolution nevertheless broadly follows tree-like patterns. Several important operations that change a general phylogenetic network have been developed in recent years, and are important for allowing algorithms to move around spaces of networks; a vital ingredient in finding an optimal network given some biological data. A key such operation is the Nearest Neighbor Interchange, or NNI. While it is already known that the space of unrooted phylogenetic networks is connected under NNI, it has been unclear whether this also holds for the subspace of tree-based networks. In this paper we show that the space of unrooted tree-based phylogenetic networks is indeed connected under the NNI operation. We do so by explicitly showing how to get from one such network to another one without losing tree-basedness along the way. Moreover, we introduce some new concepts, for instance ``shoat networks'', and derive some interesting aspects concerning tree-basedness. Last, we use our results to derive an upper bound on the size of the space of tree-based networks.


title: eCNN: A Block-Based and Highly-Parallel CNN Accelerator for Edge Inference

abstract: Convolutional neural networks (CNNs) have recently demonstrated superior quality for computational imaging applications. Therefore, they have great potential to revolutionize the image pipelines on cameras and displays. However, it is difficult for conventional CNN accelerators to support ultra-high-resolution videos at the edge due to their considerable DRAM bandwidth and power consumption. Therefore, finding a further memory- and computation-efficient microarchitecture is crucial to speed up this coming revolution. In this paper, we approach this goal by considering the inference flow, network model, instruction set, and processor design jointly to optimize hardware performance and image quality. We apply a block-based inference flow which can eliminate all the DRAM bandwidth for feature maps and accordingly propose a hardware-oriented network model, ERNet, to optimize image quality based on hardware constraints. Then we devise a coarse-grained instruction set architecture, FBISA, to support power-hungry convolution by massive parallelism. Finally,we implement an embedded processor---eCNN---which accommodates to ERNet and FBISA with a flexible processing architecture. Layout results show that it can support high-quality ERNets for super-resolution and denoising at up to 4K Ultra-HD 30 fps while using only DDR-400 and consuming 6.94W on average. By comparison, the state-of-the-art Diffy uses dual-channel DDR3-2133 and consumes 54.3W to support lower-quality VDSR at Full HD 30 fps. Lastly, we will also present application examples of high-performance style transfer and object recognition to demonstrate the flexibility of eCNN.


title: Continuum limit for discrete NLS with memory effect

abstract: We consider a discrete nonlinear Schr\"odinger equation with memory effect on the lattice $h\mathbb{Z}$ with mesh-size $h>0$. As $h\rightarrow 0$, we prove that solutions to this discrete equation converge strongly in $L^2$ to the solution to a continuous NLS-type equation with a memory effect, and compute the precise rate of convergence. In the process, we extend some recent ideas proposed by Hong and Yang in order to exploit a smoothing effect. This approach could therefore be adapted to tackle continuum limits of more general dispersive equations that require working in similar spaces.


title: Neutrino oscillations in supernovae: angular moments and fast instabilities

abstract: Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero-crossing in $n_{\nu_e} - n_{\bar{\nu}_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using Fornax simulation data, we show that fast collective oscillations qualitatively depend on how forward-peaked the neutrino angular distributions are.


title: Hardware/Software Codesign for Training/Testing Multiple Neural Networks on Multiple FPGAs

abstract: Most neural network designs for FPGAs are inflexible. In this paper, we propose a flexible VHDL structure that would allow any neural network to be implemented on multiple FPGAs. Moreover, the VHDL structure allows for testing as well as training multiple neural networks. The VHDL design consists of multiple processor groups. There are two types of processor groups: Mini Vector Machine Processor Group and Activation Processor Group. Each processor group consists of individual Mini Vector Machines and Activation Processor. The Mini Vector Machines apply vector operations to the data, while the Activation Processors apply activation functions to the data. A ring buffer was implemented to connect the various processor groups.


title: Graphene Quantum Dot as a Probe for DNA Nucleobase Detection: A First Principle Study

abstract: In this study, structural stability, electronic, optical and vibrational properties of DNA nucleobase adsorbed Graphene Quantum Dot (GQD) has been investigated using density functional theory. Based on state-of-art electronic structure calculations, we predict order of GQD sensitivity for DNA nucleobase as Thymine > Cytosine > Guanine > Adenine. An interaction of GQD with DNA nucleobase leads to modulation in electronic energy gap. Our calculated UV/vis and IR vibrational spectra show unique spectral band that can be used fingerprints in next generation DNA sequencing diagnostics.


title: A Coding-free Software Framework of Developing Web Data Management Systems

abstract: More and more enterprises recently intend to deploy data management systems in the cloud. Due to the professionalism of software development, it has still been difficult for non-programmers to develop this kind of systems, even a small one. However, the development of SaaS brings forth the more feasibility of coding-free software development than before. Based on the SaaS architecture, this paper presents a set of theory and method for coding-free construction of a data management system, on which our contributions involve in a practical application platform, a set of construction method and a set of interface on data exchange. By abstracting the common features of data management systems, we design a universal web platform to quickly generate and publish customized system instances. Moreover, we propose a kind of method to develop a data management system using a specific requirements table in spreadsheet. The corresponding platform maps the requirements table into a system instance through parsing the table model and implementing the objective system in the running stage. Finally, we implement the proposed framework and deploy it on web. The empirical result demonstrates the feasibility and availability of the coding-free method in developing web data management systems.


title: Fast Fourier Sparsity Testing

abstract: A function $f : \mathbb{F}_2^n \to \mathbb{R}$ is $s$-sparse if it has at most $s$ non-zero Fourier coefficients. Motivated by applications to fast sparse Fourier transforms over $\mathbb{F}_2^n$, we study efficient algorithms for the problem of approximating the $\ell_2$-distance from a given function to the closest $s$-sparse function. While previous works (e.g., Gopalan et al. SICOMP 2011) study the problem of distinguishing $s$-sparse functions from those that are far from $s$-sparse under Hamming distance, to the best of our knowledge no prior work has explicitly focused on the more general problem of distance estimation in the $\ell_2$ setting, which is particularly well-motivated for noisy Fourier spectra. Given the focus on efficiency, our main result is an algorithm that solves this problem with query complexity $\mathcal{O}(s)$ for constant accuracy and error parameters, which is only quadratically worse than applicable lower bounds.


title: Constraining the magnetic field structure in collisionless relativistic shocks with a radio afterglow polarization upper limit in GW170817

abstract: Gamma-ray burst afterglows arise from relativistic collisionless shocks in which the postshock tangled magnetic field $\vec{B}$ is produced by the two-stream and/or Weibel instabilities on plasma skin-depth scales $(c/\omega_p)$. The field is expected to be oriented predominantly within the shock plane ($B_{\perp}$; transverse to the shock normal, $\hat{n}_{\rm{sh}}$), and is often approximated to be completely within it ($B_\parallel\equiv\hat{n}_{\rm{sh}}\,\cdot\,\vec{B}=0$). Current 2D/3D particle-in-cell simulations are limited to short timescales and box sizes $\lesssim10^4(c/\omega_p)\ll R/\Gamma_{\rm{sh}}$ much smaller than the shocked region's comoving width, and cannot probe the asymptotic downstream $\vec{B}$ structure. We constrain the latter using the linear polarization upper limit, $|\Pi|<12\%$, on the radio afterglow of GW170817/GRB170817A. Afterglow polarization depends on the jet's angular structure, our viewing angle, and the $\vec{B}$ structure. In GW170817/GRB170817A the latter can be tightly constrained since the former two are constrained from observations. We model $\vec{B}$ as an isotropic field in 3D that is stretched along $\hat{n}_{\rm{sh}}$ by a factor $\xi\equiv{}B_{\parallel}/B_{\perp}$, whose initial value $\xi_f\equiv{}B_{\parallel,f}/B_{\perp,f}$ describes the field that survives downstream on plasma scales $\ll{}R/\Gamma_{\rm{sh}}$. We calculate $\Pi(\xi_f)$ by integrating over the shocked volume for core-dominated structured jets, with a local Blandford-McKee self-similar radial profile. We find that independent of the exact jet structure, $\vec{B}$ has a finite, but initially sub-dominant, parallel component: $0.57\lesssim\xi_f\lesssim0.89$, making it less anisotropic. While this motivates numerical studies of the asymptotic $\vec{B}$ structure in relativistic collisionless shocks, it may be consistent with turbulence amplified magnetic field.


title: Overall variation of the $H_2O$ masers around W Hydrae in 28 years

abstract: In this paper, we present the distribution of $H_2O$ masers associated with the semi-regular variable star W Hydrae (W Hya). We have collected the radio interferometric data of the maser distribution taken with the Very Large Array (VLA), the Kashima-Nobeyama InterFErometer (KNIFE), the Multi-Element Radio Link Network (MERLIN), the VLBI Exploration of Radio Astrometry (VERA), and the combined array of the Korean VLBI Network (KVN) and VERA (KaVA) in order to trace the maser distribution variation in two decades. Even though differences in the sensitivities and angular resolutions of the interferometric observations should be taken into account, we attempt to find possible correlation of the maser distribution with the stellar light curve. Our failure in the measurement of the annual parallax of the masers with VERA is likely caused by the properties of the maser features, which have been spatially resolved by the synthesized beam and survived for only half a year or less. No dependence of the maser spot flux density on its size is found in the KNIFE data, suggesting that maser spot size is determined by the physical boundary, as is expected for a clump affected by outward propagation of a stellar pulsation shock wave, rather than the (spherical) geometry of maser beaming in the maser gas clump.


title: On graphic elementary lifts of graphic matroids

abstract: Zaslavsky introduced the concept of lifted-graphic matroid. For binary matroids, a binary elementary lift can be defined in terms of the splitting operation. In this paper, we give a method to get a forbidden-minor characterization for the class of graphic matroids whose all lifted-graphic matroids are also graphic using the splitting operation.


title: Periodicity in the cohomology of finite general linear groups via q-divided powers

abstract: We show that $\bigoplus_{n \ge 0} {\mathrm H}^t({\bf GL}_n({\bf F}_q), {\bf F}_\ell)$ canonically admits the structure of a module over the $q$-divided power algebra (assuming $q$ is invertible in ${\bf F}_{\ell}$), and that, as such, it is free and (for $q \neq 2$) generated in degrees $\le t$. As a corollary, we show that the cohomology of a finitely generated ${\bf VI}$-module in non-describing characteristic is eventually periodic in $n$. We apply this to obtain a new result on the cohomology of unipotent Specht modules.


title: Analyzing User Activities Using Vector Space Model in Online Social Networks

abstract: The increasing popularity of internet, wireless technologies and mobile devices has led to the birth of mass connectivity and online interaction through Online Social Networks (OSNs) and similar environments. OSN reflects a social structure consist of a set of individuals and different types of ties like connections, relationships, interactions etc among them and helps its users to connect with their friends and common interest groups, share views and to pass information. Now days the users choose OSN sites as a most preferred place for sharing their updates, different views, posting photographs and would like to make it available for others for viewing, rating and making comments. The current paper aims to explore and analyze the association between the objects (like photographs, posts etc) and its viewers (friends, acquaintances etc) for a given user and to find activity relationship among them by using the TF-IDF scheme of Vector Space Model. After vectorization the vector data has been presented through a weighted graph with various properties.


title: Deep Markov Chain Monte Carlo

abstract: We propose a new computationally efficient sampling scheme for Bayesian inference involving high dimensional probability distributions. Our method maps the original parameter space into a low-dimensional latent space, explores the latent space to generate samples, and maps these samples back to the original space for inference. While our method can be used in conjunction with any dimension reduction technique to obtain the latent space, and any standard sampling algorithm to explore the low-dimensional space, here we specifically use a combination of auto-encoders (for dimensionality reduction) and Hamiltonian Monte Carlo (HMC, for sampling). To this end, we first run an HMC to generate some initial samples from the original parameter space, and then use these samples to train an auto-encoder. Next, starting with an initial state, we use the encoding part of the autoencoder to map the initial state to a point in the low-dimensional latent space. Using another HMC, this point is then treated as an initial state in the latent space to generate a new state, which is then mapped to the original space using the decoding part of the auto-encoder. The resulting point can be treated as a Metropolis-Hasting (MH) proposal, which is either accepted or rejected. While the induced dynamics in the parameter space is no longer Hamiltonian, it remains time reversible, and the Markov chain could still converge to the canonical distribution using a volume correction term. Dropping the volume correction step results in convergence to an approximate but reasonably accurate distribution. The empirical results based on several high-dimensional problems show that our method could substantially reduce the computational cost of Bayesian inference.


title: Radiomic Feature Stability Analysis based on Probabilistic Segmentations

abstract: Identifying image features that are robust with respect to segmentation variability and domain shift is a tough challenge in radiomics. So far, this problem has mainly been tackled in test-retest analyses. In this work we analyze radiomics feature stability based on probabilistic segmentations. Based on a public lung cancer dataset, we generate an arbitrary number of plausible segmentations using a Probabilistic U-Net. From these segmentations, we extract a high number of plausible feature vectors for each lung tumor and analyze feature variance with respect to the segmentations. Our results suggest that there are groups of radiomic features that are more (e.g. statistics features) and less (e.g. gray-level size zone matrix features) robust against segmentation variability. Finally, we demonstrate that segmentation variance impacts the performance of a prognostic lung cancer survival model and propose a new and potentially more robust radiomics feature selection workflow.


title: What are temporal correlations

abstract: Space and time are crucial twins in physics. In quantum mechanics, spatial correlations already reveal nonclassical features, such as entanglement, and have bred many quantum technologies. However, the nature of quantum temporal correlations still remains in vague. In this Letter, based on the entangled-history formalism, we prove rigorously that temporal correlations are equivalent to spatial correlations. The effect of temporal correlations corresponds to a quantum channel. The resulting quantifications and classifications of quantum temporal correlations are illustrated in a natural way. Our proposed procedures also show how to determine temporal correlations completely.


title: Bayesian Neural Decoding Using A Diversity-Encouraging Latent Representation Learning Method

abstract: It is well established that temporal organization is critical to memory, and that the ability to temporally organize information is fundamental to many perceptual, cognitive, and motor processes. While our understanding of how the brain processes the spatial context of memories has advanced considerably, our understanding of their temporal organization lags far behind. In this paper, we propose a new approach for elucidating the neural basis of complex behaviors and temporal organization of memories. More specifically, we focus on neural decoding - the prediction of behavioral or experimental conditions based on observed neural data. In general, this is a challenging classification problem, which is of immense interest in neuroscience. Our goal is to develop a new framework that not only improves the overall accuracy of decoding, but also provides a clear latent representation of the decoding process. To accomplish this, our approach uses a Variational Auto-encoder (VAE) model with a diversity-encouraging prior based on determinantal point processes (DPP) to improve latent representation learning by avoiding redundancy in the latent space. We apply our method to data collected from a novel rat experiment that involves presenting repeated sequences of odors at a single port and testing the rats' ability to identify each odor. We show that our method leads to substantially higher accuracy rate for neural decoding and allows to discover novel biological phenomena by providing a clear latent representation of the decoding process.


title: Automatic Differentiation of Some First-Order Methods in Parametric Optimization

abstract: We aim at computing the derivative of the solution to a parametric optimization problem with respect to the involved parameters. For a class broader than that of strongly convex functions, this can be achieved by automatic differentiation of iterative minimization algorithms. If the iterative algorithm converges pointwise, then we prove that the derivative sequence also converges pointwise to the derivative of the minimizer with respect to the parameters. Moreover, we provide convergence rates for both sequences. In particular, we prove that the accelerated convergence rate of the Heavy-ball method compared to Gradient Descent also accelerates the derivative computation. An experiment with L2-Regularized Logistic Regression validates the theoretical results.


title: Generalization Bounds for Neural Networks via Approximate Description Length

abstract: We investigate the sample complexity of networks with bounds on the magnitude of its weights. In particular, we consider the class \[ H=\left\{W_t\circ\rho\circ \ldots\circ\rho\circ W_{1} :W_1,\ldots,W_{t-1}\in M_{d, d}, W_t\in M_{1,d}\right\} \] where the spectral norm of each $W_i$ is bounded by $O(1)$, the Frobenius norm is bounded by $R$, and $\rho$ is the sigmoid function $\frac{e^x}{1+e^x}$ or the smoothened ReLU function $ \ln (1+e^x)$. We show that for any depth $t$, if the inputs are in $[-1,1]^d$, the sample complexity of $H$ is $\tilde O\left(\frac{dR^2}{\epsilon^2}\right)$. This bound is optimal up to log-factors, and substantially improves over the previous state of the art of $\tilde O\left(\frac{d^2R^2}{\epsilon^2}\right)$. We furthermore show that this bound remains valid if instead of considering the magnitude of the $W_i$'s, we consider the magnitude of $W_i - W_i^0$, where $W_i^0$ are some reference matrices, with spectral norm of $O(1)$. By taking the $W_i^0$ to be the matrices at the onset of the training process, we get sample complexity bounds that are sub-linear in the number of parameters, in many typical regimes of parameters. To establish our results we develop a new technique to analyze the sample complexity of families $H$ of predictors. We start by defining a new notion of a randomized approximate description of functions $f:X\to\mathbb{R}^d$. We then show that if there is a way to approximately describe functions in a class $H$ using $d$ bits, then $d/\epsilon^2$ examples suffices to guarantee uniform convergence. Namely, that the empirical loss of all the functions in the class is $\epsilon$-close to the true loss. Finally, we develop a set of tools for calculating the approximate description length of classes of functions that can be presented as a composition of linear function classes and non-linear functions.


title: Terahertz dual phase gradient metasurface: high-efficiency binary-channel spoof surface plasmon excitation

abstract: Spoof surface plasmon meta-couplers are compact antennas which link propagating waves and surface waves. However, most of them are designed with a fixed phase gradient and channel for the incident waves with specific polarization, which limits their further applications in multichannel scenarios. In this Letter, we propose, to the best of our knowledge, a new method that combines the Brillouin folds theory with the Generalized Snell Law. We demonstrate that when the phase gradient of the metasurface is large enough, Brillouin folds effect will occur, which will create dual phase gradient space in a single metasurface. With this method, we design two novel terahertz meta-couplers with functionalities of symmetrical and asymmetrical binary-channel/bidirectional SSP excitation. Furthermore, finite element method (FEM) simulations are performed to demonstrate their functionalities. Considering the orthogonality of the incident waves, there can be a total of four independent space channels to excite SSP on one metasurface. This work may open up new routes in multi-channel SSP meta-couplers and multi-beam surface wave antennas.


title: Quantum-Inspired Classical Algorithms for Singular Value Transformation

abstract: A recent breakthrough by Tang (STOC 2019) showed how to "dequantize" the quantum algorithm for recommendation systems by Kerenidis and Prakash (ITCS 2017). The resulting algorithm, classical but "quantum-inspired", efficiently computes a low-rank approximation of the users' preference matrix. Subsequent works have shown how to construct efficient quantum-inspired algorithms for approximating the pseudo-inverse of a low-rank matrix as well, which can be used to (approximately) solve low-rank linear systems of equations. In the present paper, we pursue this line of research and develop quantum-inspired algorithms for a large class of matrix transformations that are defined via the singular value decomposition of the matrix. In particular, we obtain classical algorithms with complexity polynomially related (in most parameters) to the complexity of the best quantum algorithms for singular value transformation recently developed by Chakraborty, Gily\'{e}n and Jeffery (ICALP 2019) and Gily\'{e}n, Su, Low and Wiebe (STOC19).


title: A Novel Self-Supervised Re-labeling Approach for Training with Noisy Labels

abstract: The major driving force behind the immense success of deep learning models is the availability of large datasets along with their clean labels. Unfortunately, this is very difficult to obtain, which has motivated research on the training of deep models in the presence of label noise and ways to avoid over-fitting on the noisy labels. In this work, we build upon the seminal work in this area, Co-teaching and propose a simple, yet efficient approach termed mCT-S2R (modified co-teaching with self-supervision and re-labeling) for this task. First, to deal with significant amount of noise in the labels, we propose to use self-supervision to generate robust features without using any labels. Next, using a parallel network architecture, an estimate of the clean labeled portion of the data is obtained. Finally, using this data, a portion of the estimated noisy labeled portion is re-labeled, before resuming the network training with the augmented data. Extensive experiments on three standard datasets show the effectiveness of the proposed framework.


title: Phase Transitions in Genome-wide Association Studies and Categorical Variable Screenings

abstract: Motivated by genome-wide association screening studies (GWAS), we study high-dimensional marginal screenings of categorical variables where test statistics have approximate chi-square distributions. We characterize four new phase transitions in high-dimensional chi-square models, and derive the signal sizes necessary and sufficient for statistical procedures to simultaneously control false discovery (in terms of family-wise error rate or false discovery rate) and missed detection (in terms of family-wise non-discovery rate or false non-discovery rate) in large dimensions. Remarkably, degrees of freedom in the chi-square distributions do not affect the boundaries in all four phase transitions. Several well-known procedures are shown to attain these boundaries. Two new phase transitions are also identified in the Gaussian location model under one-sided alternatives. We then elucidate on the nature of signal sizes in association tests by characterizing its relationship with marginal frequencies, odds ratio, and sample sizes in $2\times2$ contingency tables. This allows us to illustrate an interesting manifestation of the phase transition phenomena in genome-wide association studies (GWAS). We also show, perhaps surprisingly, that given total sample sizes, balanced designs in such association studies rarely deliver optimal power for detecting the effects of rare genetic variants.


title: The Largest-Z-ratio-First algorithm is 0.8531-approximate for scheduling unreliable jobs on m parallel machines

abstract: In this paper we analyze the worst-case performance of a greedy algorithm called Largest-Z-ratio-First for the problem of scheduling unreliable jobs on m parallel machines. Each job is characterized by a success probability and a reward earned in the case of success. In the case of failure, the jobs subsequently sequenced on that machine cannot be performed. The objective is to maximize the expected reward. We show the algorithm provides an approximation ratio of approximately 0.853196, and that the bound is tight.


title: Global optimization-based dimer method for finding saddle points

abstract: Searching saddle points on the potential energy surface is a challenging problem in the rare event. When there exist multiple saddle points, sampling different initial guesses are needed in most dimer-type methods in order to find distinct saddle points. In this paper, we present a novel global optimization-based dimer method (GOD) to efficiently search saddle points by coupling ant colony optimization (ACO) algorithm with optimization-based shrinking dimer (OSD) method. In particular, we apply OSD method as a local search algorithm for saddle points and construct a pheromone function in ACO to update the global population. By applying a two-dimensional example and a benchmark problem of seven-atom island on the (111) surface of an FCC crystal, we demonstrate that GOD shows a significant improvement in computational efficiency compared with OSD method. Our algorithm offers a new framework to open up possibilities of adopting other global optimization methods to search saddle points.


title: Contour Sparse Representation with SDD Features for Object Recognition

abstract: Slope difference distribution (SDD) is computed for the one-dimensional curve. It is not only robust to calculate the partitioning point to separate the curve logically, but also robust to calculate the clustering center of each part of the separated curve. SDD has been proposed for image segmentation and it outperforms all existing image segmentation methods. For verification purpose, we have made the Matlab codes of comparing SDD method with existing image segmentation methods freely available at Matlab Central. The contour of the object is similar to the histogram of the image. Thus, feature detection by SDD from the contour of the object is also feasible. In this letter, SDD features are defined and they form the sparse representation of the object contour. The reference model of each object is built based on the SDD features and then model matching is used for on line object recognition. The experimental results are very encouraging. For the gesture recognition, SDD achieved 100% accuracy for two public datasets: the NUS dataset and the near-infrared dataset. For the object recognition, SDD achieved 100% accuracy for the Kimia 99 dataset.


title: Beyond 5G: Leveraging Cell Free TDD Massive MIMO using Cascaded Deep learning

abstract: This paper deals with the calibration of Time Division Duplexing (TDD) reciprocity in an Orthogonal Frequency Division Multiplexing (OFDM) based Cell Free Massive MIMO system where the responses of the (Radio Frequency) RF chains render the end to end channel non-reciprocal, even though the physical wireless channel is reciprocal. We further address the non-availability of the uplink channel estimates at locations other than pilot subcarriers and propose a single-shot solution to estimate the downlink channel at all subcarriers from the uplink channel at selected pilot subcarriers. We propose a cascade of two Deep Neural Networks (DNN) to achieve the objective. The proposed method is easily scalable and removes the need for relative reciprocity calibration based on the cooperation of antennas, which usually introduces dependency in Cell Free Massive MIMO systems.


title: Residue formula for an obstruction to Coupled K\"ahler-Einstein metrics

abstract: We obtain a residue formula for an obstruction to the existence of coupled K\"ahler-Einstein metrics described by Futaki-Zhang. We apply it to an example studied separately by Futaki and Hultgren which is a toric Fano manifold with reductive automorphism, does not admit a K\"ahler-Einstein metric but still admits coupled K\"ahler-Einstein metrics.


title: A Conservative Finite Element Method for the Incompressible Euler Equations with Variable Density

abstract: We construct a finite element discretization and time-stepping scheme for the incompressible Euler equations with variable density that exactly preserves total mass, total squared density, total energy, and pointwise incompressibility. The method uses Raviart-Thomas or Brezzi-Douglas-Marini finite elements to approximate the velocity and discontinuous polynomials to approximate the density and pressure. To achieve exact preservation of the aforementioned conserved quantities, we exploit a seldom-used weak formulation of the momentum equation and a second-order time-stepping scheme that is similar, but not identical, to the midpoint rule. We also describe and prove stability of an upwinded version of the method. We present numerical examples that demonstrate the order of convergence of the method.


title: Thermodynamics and statistical physics of quasiparticles within the quark-gluon plasma model

abstract: We consider thermodynamic properties of a quark-gluon plasma related to quasiparticles having the internal structure. For this purpose, we employ a possible analogy between quantum chromodynamics and non-Abelian Proca-Dirac-Higgs theory. The influence of characteristic sizes of the quasiparticles on such thermodynamic properties of the quark-gluon plasma like the internal energy and pressure is studied. Sizes of the quasiparticles are taken into account in the spirit of the van der Waals equation but we take into consideration that the quasiparticles have different sizes, and the average value of these sizes depends on temperature. It is shown that this results in a change in the internal energy and pressure of the quark-gluon plasma. Also, we show that, when the temperature increases, the average value of characteristic sizes of the quasiparticles increases as well. This leads to the occurrence of a phase transition at the temperature at which the volume occupied by the quasiparticles is compared with the volume occupied by the plasma.


title: Variational Phase Field Formulations of Polarization and Phase Transition in Ferroelectric Thin Films

abstract: Electric field plays an important role in ferroelectric phase transition. There have been numerous phase field formulations attempting to account for electrostatic interactions subject to different boundary conditions. In this paper, we develop new variational forms of the phase field electrostatic energy and the relaxation dynamics of the polarization vector that involves a hybrid representation in both real and Fourier variables. The new formulations avoid ambiguities appeared in earlier studies and lead to much more effective ways to perform variational studies and numerical simulations. Computations of polarization switching in a single domain by applying the new formulations are provided as illustrative examples.


title: Multi-layered Kelvin-Helmholtz Instability in the Solar Corona

abstract: The Kelvin-Helmholtz (KH) instability is commonly found in many astrophysical, laboratory, and space plasmas. It could mix plasma components of different properties and convert dynamic fluid energy from large scale structure to smaller ones. In this study, we combined the ground-based New Vacuum Solar Telescope (NVST) and the Solar Dynamic Observatories (SDO) / Atmospheric Imaging Assembly (AIA) to observe the plasma dynamics associated with active region 12673 on 09 September 2017. In this multi-temperature view, we identified three adjacent layers of plasma flowing at different speeds, and detected KH instabilities at their interfaces. We could unambiguously track a typical KH vortex and measure its motion. We found that the speed of this vortex suddenly tripled at a certain stage. This acceleration was synchronized with the enhancements in emission measure and average intensity of the 193 \AA{} data. We interpret this as evidence that KH instability triggers plasma heating. The intriguing feature in this event is that the KH instability observed in the NVST channel was nearly complementary to that in the AIA 193 \AA{}. Such a multi-thermal energy exchange process is easily overlooked in previous studies, as the cold plasma component is usually not visible in the extreme ultraviolet channels that are only sensitive to high temperature plasma emissions. Our finding indicates that embedded cold layers could interact with hot plasma as invisible matters. We speculate that this process could occur at a variety of length scales and could contribute to plasma heating.


title: Blois 2019: highlights and outlook

abstract: This is an idiosyncratic account of the main results presented at the 31st Rencontres de Blois, which took place from June 2nd to June 7th, 2019 in the Castle of Blois, France.