title: MOTS: Minimax Optimal Thompson Sampling

abstract: Thompson sampling is one of the most widely used algorithms for many online decision problems, due to its simplicity in implementation and superior empirical performance over other state-of-the-art methods. Despite its popularity and empirical success, it has remained an open problem whether Thompson sampling can match the minimax lower bound $\Omega(\sqrt{KT})$ for $K$-armed bandit problems, where $T$ is the total time horizon. In this paper, we solve this long open problem by proposing a variant of Thompson sampling called MOTS that adaptively clips the sampling instance of the chosen arm at each time step. We prove that this simple variant of Thompson sampling achieves the minimax optimal regret bound $O(\sqrt{KT})$ for finite time horizon $T$, as well as the asymptotic optimal regret bound for Gaussian rewards when $T$ approaches infinity. To our knowledge, MOTS is the first Thompson sampling type algorithm that achieves the minimax optimality for multi-armed bandit problems.


title: Prediction of Time-to-terminal Event (TTTE) in a Class of Joint Dynamic Models

abstract: In different areas of research, multiple recurrent competing risks (RCR) are often observed on the same observational unit. For instance, different types of cancer relapses are observed on the same patient and several types of component failures are observed in the same reliability system. When a terminal event (TE) such as death is also observed on the same unit, since the RCRs are generally informative about death, we develop joint dynamic models that simultaneously model the RCRs and the TE. A key interest of such joint dynamic modeling is to predict time-to-terminal event (TTTE) for new units that have not experienced the TE by the end of monitoring period. In this paper, we propose a simulation approach to predict TTTE which arises from a class of joint dynamic models of RCRs and TE. The proposed approach can be applied to problems in precision medicine and potentially many other settings. The simulation method makes personalized predictions of TTTE and provides an empirical predictive distribution of TTTE. Predictions of the RCR occurrences beyond a possibly random monitoring time and leading up to the TE occurrence are also produced. The approach is dynamic in that each simulated occurrence of RCR increases the amount of knowledge we obtain on an observational unit which informs the simulation of TTTE. We demonstrate the approach on a synthetic dataset and evaluate predictive accuracy of the prediction method through 5-fold cross-validation using empirical Brier score.


title: Adaptive Hyper-box Matching for Interpretable Individualized Treatment Effect Estimation

abstract: We propose a matching method for observational data that matches units with others in unit-specific, hyper-box-shaped regions of the covariate space. These regions are large enough that many matches are created for each unit and small enough that the treatment effect is roughly constant throughout. The regions are found as either the solution to a mixed integer program, or using a (fast) approximation algorithm. The result is an interpretable and tailored estimate of a causal effect for each unit.


title: The nature of a primary jet within a circumbinary disc outflow in a young stellar system

abstract: Most stars form in binaries, and both stars may grow by accreting material from a circumbinary disc onto their personal discs. We suspect that in many cases a wide molecular wind will envelope a collimated atomic jet emanating from close to an orbiting young star. This so-called Circumbinary Scenario is explored here in order to find common identifiable properties. The dynamical set up is studied with three dimensional simulations with chemistry and cooling included. We extract the properties on scales of order 100\,AU and compare to the Co-Orbital Scenario in which the wind and jet sources are in orbit. We find that the rapid orbital motion generates a wide ionised sheath around the jet core with a large opening angle at the base. This is independent of the presence of the surrounding molecular outflow. However, the atomic jet is recollimated beyond ~ 55 AU when the molecular outflow restricts the motion of the ambient medium which, in turn, confines the jet. These physical properties are relayed to the optical H-alpha imaging, providing a means of distinguishing between models. The high excitation sheath and recollimation region can be explored on these scales through the next generation of instruments. However, in general, the amount and location of the ionised material, whether in the knots or the sheath, will depend on several parameters including the orbital period, axis alignment and pulse amplitude.


title: Horizons Protect Church-Turing

abstract: The quantum-Extended Church-Turing thesis is a principle of physics as well as computer science. It asserts that the laws of physics will prevent the construction of a machine that can efficiently determine the results of any calculation which cannot be done efficiently by a quantum Turing machine (or a universal quantum circuit). In this note I will argue that an observer falling into a black hole can learn the result of such a calculation in a very short time, thereby seemingly violating the thesis. A viable reformulation requires that the thesis only applies to observers who have access to the holographic boundary of space. The properties of the horizon play a crucial a role in protecting the thesis. The arguments are closely related to, and were partially motivated by a recent paper by Bouland, Fefferman, and Vazirani, and by a question raised by Aaronson.


title: Decay and recurrence of non-Gaussian correlations in a quantum many-body system

abstract: Gaussian models provide an excellent effective description of many quantum many-body systems ranging from condensed matter systems all the way to neutron stars. Gaussian states are common at equilibrium when the interactions are weak. Recently it was proposed that they can also emerge dynamically from a non-Gaussian initial state evolving under non-interacting dynamics. In this work, we present the experimental observation of such a dynamical emergence of Gaussian correlations in a quantum many-body system. This non-equilibrium evolution is triggered by abruptly switching off the effective interaction between the observed collective degrees of freedom, while leaving the interactions between the microscopic constituents unchanged. Starting from highly non-Gaussian correlations, consistent with the sine-Gordon model, we observe a Gaussian state to emerge over time as revealed by the decay of the fourth and sixth order connected correlations in the quantum field. A description of this dynamics requires a novel mechanism for the emergence of Gaussian correlations, which is relevant for a wide class of quantum many-body systems. In our closed system with non-interacting effective degrees of freedom, we do not expect full thermalization. This memory of the initial state is confirmed by observing recurrences of non-Gaussian correlations.


title: Numerical Solution of Dynamic Portfolio Optimization with Transaction Costs

abstract: We apply numerical dynamic programming techniques to solve discrete-time multi-asset dynamic portfolio optimization problems with proportional transaction costs and shorting/borrowing constraints. Examples include problems with multiple assets, and many trading periods in a finite horizon problem. We also solve dynamic stochastic problems, with a portfolio including one risk-free asset, an option, and its underlying risky asset, under the existence of transaction costs and constraints. These examples show that it is now tractable to solve such problems.

categories: q-fin.PM econ.GN q-fin.CP q-fin.EC q-fin.MF


title: Curvature-Induced Skyrmion Mass

abstract: We investigate the propagation of magnetic skyrmions on elastically deformable geometries by employing imaginary time quantum field theory methods. We demonstrate that the Euclidean action of the problem carries information of the elements of the surface space metric, and develop a description of the skyrmion dynamics in terms of a set of collective coordinates. We reveal that novel curvature-driven effects emerge in geometries with non-constant curvature, which explicitly break the translational invariance of flat space. In particular, for a skyrmion stabilized by a curvilinear defect, an inertia term and a pinning potential are generated by the varying curvature, while both of these terms vanish in the flat-space limit.


title: RMP-SNN: Residual Membrane Potential Neuron for Enabling Deeper High-Accuracy and Low-Latency Spiking Neural Network

abstract: Spiking Neural Networks (SNNs) have recently attracted significant research interest as the third generation of artificial neural networks that can enable low-power event-driven data analytics. The best performing SNNs for image recognition tasks are obtained by converting a trained Analog Neural Network (ANN), consisting of Rectified Linear Units (ReLU), to SNN composed of integrate-and-fire neurons with "proper" firing thresholds. The converted SNNs typically incur loss in accuracy compared to that provided by the original ANN and require sizable number of inference time-steps to achieve the best accuracy. We find that performance degradation in the converted SNN stems from using "hard reset" spiking neuron that is driven to fixed reset potential once its membrane potential exceeds the firing threshold, leading to information loss during SNN inference. We propose ANN-SNN conversion using "soft reset" spiking neuron model, referred to as Residual Membrane Potential (RMP) spiking neuron, which retains the "residual" membrane potential above threshold at the firing instants. We demonstrate near loss-less ANN-SNN conversion using RMP neurons for VGG-16, ResNet-20, and ResNet-34 SNNs on challenging datasets including CIFAR-10 (93.63% top-1), CIFAR-100 (70.93% top-1), and ImageNet (73.09% top-1 accuracy). Our results also show that RMP-SNN surpasses the best inference accuracy provided by the converted SNN with "hard reset" spiking neurons using 2-8 times fewer inference time-steps across network architectures and datasets.


title: The case for algebraic biology: from research to education

abstract: Though it goes without saying that linear algebra is fundamental to mathematical biology, polynomial algebra is less visible. In this article, we will give a brief tour of four diverse biological problems where multivariate polynomials play a central role -- a subfield that is sometimes called "algebraic biology." Namely, these topics include biochemical reaction networks, Boolean models of gene regulatory networks, algebraic statistics and genomics, and place fields in neuroscience. After that, we will summarize the history of discrete and algebraic structures in mathematical biology, from their early appearances in the late 1960s to the current day. Finally, we will discuss the role of algebraic biology in the modern classroom and curriculum, including resources in the literature and relevant software. Our goal is to make this article widely accessible, reaching the mathematical biologist who knows no algebra, the algebraist who knows no biology, and especially the interested student who is curious about the synergy between these two seemingly unrelated fields.


title: Accuracy of Difference Schemes in Electromagnetic Applications: a Trefftz Analysis

abstract: The paper examines local approximation errors of finite difference schemes in electromagnetic analysis. Despite a long history of the subject, several accuracy-related issues have been overlooked and/or remain controversial. For example, conflicting claims have been made in the literature about the order of Yee-like schemes in the vicinity of slanted or curved material interfaces. Two novel practical methods for comparison of local accuracy of difference schemes are proposed: one makes use of Trefftz test matrices, and the other one relies on a new measure of approximation accuracy in scheme-exact Trefftz subspaces. One particular conclusion is that a loss of accuracy for Yee-like schemes at slanted material boundaries is unavoidable.


title: ExoReL$^\Re$: A Bayesian Inverse Retrieval Framework For Exoplanetary Reflected Light Spectra

abstract: The high-contrast imaging technique is meant to provide insight into those planets orbiting several astronomical units from their host star. Space missions such as WFIRST, HabEx, and LUVOIR will measure reflected light spectra of cold gaseous and rocky planets. To interpret these observations we introduce ExoReL$^\Re$ (Exoplanetary Reflected Light Retrieval), a novel Bayesian retrieval framework to retrieve cloud properties and atmospheric structures from exoplanetary reflected light spectra. As a unique feature, it assumes a vertically non-uniform volume mixing ratio profile of water and ammonia, and use it to construct cloud densities. In this way, clouds and molecular mixture ratios are consistent. We apply ExoReL$^\Re$ on three test cases: two exoplanets ($\upsilon$ And e and 47 Uma b) and Jupiter. We show that we are able to retrieve the concentration of methane in the atmosphere, and estimate the position of clouds when the S/N of the spectrum is higher than 15, in line with previous works. Moreover, we described the ability of our model of giving a chemical identity to clouds, and we discussed whether or not we can observe this difference in the planetary reflection spectrum. Finally, we demonstrate how it could be possible to retrieve molecular concentrations (water and ammonia in this work) below the clouds by linking the non-uniform volume mixing ratio profile to the cloud presence. This will help to constrain the concentration of water and ammonia unseen in direct measurements.


title: Energy conditions in general relativity and quantum field theory

abstract: This review summarizes the current status of the energy conditions in general relativity and quantum field theory. We provide a historical review and a summary of technical results and applications, complemented with a few new derivations and discussions. We pay special attention to the role of the equations of motion and to the relation between classical and quantum theories. Pointwise energy conditions were first introduced as physically reasonable restrictions on matter in the context of general relativity. They aim to express e.g. the positivity of mass or the attractiveness of gravity. Perhaps more importantly, they have been used as assumptions in mathematical relativity to prove singularity theorems and the non-existence of wormholes and similar exotic phenomena. However, the delicate balance between conceptual simplicity, general validity and strong results has faced serious challenges, because all pointwise energy conditions are systematically violated by quantum fields and also by some rather simple classical fields. In response to these challenges, weaker statements were introduced, such as quantum energy inequalities and averaged energy conditions. These have a larger range of validity and may still suffice to prove at least some of the earlier results. One of these conditions, the achronal averaged null energy condition, has recently received increased attention. It is expected to be a universal property of the dynamics of all gravitating physical matter, even in the context of semiclassical or quantum gravity.


title: RODNet: Radar Object Detection Using Cross-Modal Supervision

abstract: Radar is usually more robust than the camera in severe driving scenarios, e.g., weak/strong lighting and bad weather. However, unlike RGB images captured by a camera, the semantic information from the radar signals is noticeably difficult to extract. In this paper, we propose a deep radar object detection network (RODNet), to effectively detect objects purely from the carefully processed radar frequency data in the format of range-azimuth frequency heatmaps (RAMaps). Three different 3D autoencoder based architectures are introduced to predict object confidence distribution from each snippet of the input RAMaps. The final detection results are then calculated using our post-processing method, called location-based non-maximum suppression (L-NMS). Instead of using burdensome human-labeled ground truth, we train the RODNet using the annotations generated automatically by a novel 3D localization method using a camera-radar fusion (CRF) strategy. To train and evaluate our method, we build a new dataset -- CRUW, containing synchronized videos and RAMaps in various driving scenarios. After intensive experiments, our RODNet shows favorable object detection performance without the presence of the camera.


title: Unveiling Defect-Mediated Carrier Dynamics in Monolayer Semiconductors by Spatiotemporal Microwave Imaging

abstract: The optoelectronic properties of atomically thin transition-metal dichalcogenides are strongly correlated with the presence of defects in the materials, which are not necessarily detrimental for certain applications. For instance, defects can lead to an enhanced photoconduction, a complicated process involving charge generation and recombination in the time domain and carrier transport in the spatial domain. Here, we report the simultaneous spatial and temporal photoconductivity imaging in two types of WS2 monolayers by laser-illuminated microwave impedance microscopy. The diffusion length and carrier lifetime were directly extracted from the spatial profile and temporal relaxation of microwave signals respectively. Time-resolved experiments indicate that the critical process for photo-excited carriers is the escape of holes from trap states, which prolongs the apparent lifetime of mobile electrons in the conduction band. As a result, counterintuitively, the photoconductivity is stronger in CVD samples than exfoliated monolayers with a lower defect density. Our work reveals the intrinsic time and length scales of electrical response to photo-excitation in van der Waals materials, which is essential for their applications in novel optoelectronic devices.


title: Building a larger class of graphs for efficient reconfiguration of vertex colouring

abstract: A $k$-colouring of a graph $G$ is an assignment of at most $k$ colours to the vertices of $G$ so that adjacent vertices are assigned different colours. The reconfiguration graph of the $k$-colourings, $\mathcal{R}_k(G)$, is the graph whose vertices are the $k$-colourings of $G$ and two colourings are joined by an edge in $\mathcal{R}_k(G)$ if they differ in colour on exactly one vertex. For a $k$-colourable graph $G$, we investigate the connectivity and diameter of $\mathcal{R}_{k+1}(G)$. It is known that not all weakly chordal graphs have the property that $\mathcal{R}_{k+1}(G)$ is connected. On the other hand, $\mathcal{R}_{k+1}(G)$ is connected and of diameter $O(n^2)$ for several subclasses of weakly chordal graphs such as chordal, chordal bipartite, and $P_4$-free graphs. We introduce a new class of graphs called OAT graphs that extends the latter classes and in fact extends outside the class of weakly chordal graphs. OAT graphs are built from four simple operations, disjoint union, join, and the addition of a clique or comparable vertex. We prove that if $G$ is a $k$-colourable OAT graph then $\mathcal{R}_{k+1}(G)$ is connected with diameter $O(n^2)$. Furthermore, we give polynomial time algorithms to recognize OAT graphs and to find a path between any two colourings in $\mathcal{R}_{k+1}(G)$.


title: Hopf Galois structures on extensions of degree twice an odd prime square and their associated left braces

abstract: We determine the Hopf Galois structures on a Galois field extension of degree twice an odd prime square and classify the corresponding left braces. Besides we determine the separable field extensions of degree twice an odd prime square allowing a cyclic Hopf Galois structure and the number of these structures.


title: Robust Market Making via Adversarial Reinforcement Learning

abstract: We show that adversarial reinforcement learning (ARL) can be used to produce market marking agents that are robust to adversarial and adaptively-chosen market conditions. To apply ARL, we turn the well-studied single-agent model of Avellaneda and Stoikov [2008] into a discrete-time zero-sum game between a market maker and adversary. The adversary acts as a proxy for other market participants that would like to profit at the market maker's expense. We empirically compare two conventional single-agent RL agents with ARL, and show that our ARL approach leads to: 1) the emergence of risk-averse behaviour without constraints or domain-specific penalties; 2) significant improvements in performance across a set of standard metrics, evaluated with or without an adversary in the test environment, and; 3) improved robustness to model uncertainty. We empirically demonstrate that our ARL method consistently converges, and we prove for several special cases that the profiles that we converge to correspond to Nash equilibria in a simplified single-stage game.


title: HyperEmbed: Tradeoffs Between Resources and Performance in NLP Tasks with Hyperdimensional Computing enabled Embedding of n-gram Statistics

abstract: Recent advances in Deep Learning have led to a significant performance increase on several NLP tasks, however, the models become more and more computationally demanding. Therefore, this paper tackles the domain of computationally efficient algorithms for NLP tasks. In particular, it investigates distributed representations of n-gram statistics of texts. The representations are formed using hyperdimensional computing enabled embedding. These representations then serve as features, which are used as input to standard classifiers. We investigate the applicability of the embedding on one large and three small standard datasets for classification tasks using nine classifiers. The embedding achieved on par F1 scores while decreasing the time and memory requirements by several times compared to the conventional n-gram statistics, e.g., for one of the classifiers on a small dataset, the memory reduction was 6.18 times; while train and test speed-ups were 4.62 and 3.84 times, respectively. For many classifiers on the large dataset, memory reduction was ca. 100 times and train and test speed-ups were over 100 times. Importantly, the usage of distributed representations formed via hyperdimensional computing allows dissecting strict dependency between the dimensionality of the representation and n-gram size, thus, opening a room for tradeoffs.


title: Implicitly Defined Layers in Neural Networks

abstract: In conventional formulations of multilayer feedforward neural networks, the individual layers are customarily defined by explicit functions. In this paper we demonstrate that defining individual layers in a neural network \emph{implicitly} provide much richer representations over the standard explicit one, consequently enabling a vastly broader class of end-to-end trainable architectures. We present a general framework of implicitly defined layers, where much of the theoretical analysis of such layers can be addressed through the implicit function theorem. We also show how implicitly defined layers can be seamlessly incorporated into existing machine learning libraries. In particular with respect to current automatic differentiation techniques for use in backpropagation based training. Finally, we demonstrate the versatility and relevance of our proposed approach on a number of diverse example problems with promising results.


title: A New Approach to Solar Flare Prediction

abstract: All three components of the current density are required to compute the heating rate due to free magnetic energy dissipation. Here we present a first test of a new model developed to determine if the times of increases in the resistive heating rate in active region (AR) photospheres are correlated with the subsequent occurrence of M and X flares in the corona. A data driven, 3 D, non-force-free magnetohydrodynamic model restricted to the near-photospheric region is used to compute time series of the complete current density and the resistive heating rate per unit volume $(Q(t))$ in each pixel in neutral line regions (NLRs) of 14 ARs. The model is driven by time series of the magnetic field ${\bf B}$ measured by the Helioseismic & Magnetic Imager on the Solar Dynamics Observatory (SDO) satellite. Spurious Doppler periods due to SDO orbital motion are filtered out of the time series for ${\bf B}$ in every AR pixel. For each AR, the cumulative distribution function (CDF) of the values of the NLR area integral $Q_i(t)$ of $Q(t)$ is found to be a scale invariant power law distribution essentially identical to the observed CDF for the total energy released in coronal flares. This suggests that coronal flares and the photospheric $Q_i$ are correlated, and powered by the same process. The model predicts spikes in $Q_i$ with values orders of magnitude above background values. These spikes are driven by spikes in the non-force free component of the current density. The times of these spikes are plausibly correlated with times of subsequent M or X flares a few hours to a few days later. The spikes occur on granulation scales, and may be signatures of heating in horizontal current sheets. It is also found that the times of relatively large values of the rate of change of the NLR unsigned magnetic flux are also plausibly correlated with the times of subsequent M and X flares, and spikes in $Q_i$.


title: Dependency-Aware Release Planning for Software Projects using Fuzzy Graphs and Integer Programming

abstract: Software Release Planning (SRP) is to find, for the software, a subset of the requirements with the highest value while respecting the budget. The value of a requirement however may, to various degrees, depend on selecting or ignoring other requirements. However, existing SRP models ignore either Value-Related Dependencies altogether or the strengths of those dependencies even if they consider them. This paper presents an Integer Programming model for software release planning, which considers the variances of strengths of positive and negative value-related dependencies among software requirements. To capture the imprecision associated with strengths of value-related dependencies we have made use of the algebraic structure of fuzzy graphs. We have further, contributed a scalable technique for automated identification of value-related dependencies based on user preferences for software requirements. The validity of the work is verified through simulations.


title: Scaling MAP-Elites to Deep Neuroevolution

abstract: Quality-Diversity (QD) algorithms, and MAP-Elites (ME) in particular, have proven very useful for a broad range of applications including enabling real robots to recover quickly from joint damage, solving strongly deceptive maze tasks or evolving robot morphologies to discover new gaits. However, present implementations of MAP-Elites and other QD algorithms seem to be limited to low-dimensional controllers with far fewer parameters than modern deep neural network models. In this paper, we propose to leverage the efficiency of Evolution Strategies (ES) to scale MAP-Elites to high-dimensional controllers parameterized by large neural networks. We design and evaluate a new hybrid algorithm called MAP-Elites with Evolution Strategies (ME-ES) for post-damage recovery in a difficult high-dimensional control task where traditional ME fails. Additionally, we show that ME-ES performs efficient exploration, on par with state-of-the-art exploration algorithms in high-dimensional control tasks with strongly deceptive rewards.


title: Watch your Up-Convolution: CNN Based Generative Deep Neural Networks are Failing to Reproduce Spectral Distributions

abstract: Generative convolutional deep neural networks, e.g. popular GAN architectures, are relying on convolution based up-sampling methods to produce non-scalar outputs like images or video sequences. In this paper, we show that common up-sampling methods, i.e. known as up-convolution or transposed convolution, are causing the inability of such models to reproduce spectral distributions of natural training data correctly. This effect is independent of the underlying architecture and we show that it can be used to easily detect generated data like deepfakes with up to 100% accuracy on public benchmarks. To overcome this drawback of current generative models, we propose to add a novel spectral regularization term to the training optimization objective. We show that this approach not only allows to train spectral consistent GANs that are avoiding high frequency errors. Also, we show that a correct approximation of the frequency spectrum has positive effects on the training stability and output quality of generative networks.


title: Gauss and the identity function -- a tale of characterizations of the normal distribution

abstract: The normal distribution is well-known for several results that it is the only to fulfil. The aim of the present paper is to show that many of these characterizations actually follow from the fact that the derivative of the log-density of the normal distribution is the (negative) identity function. This \emph{a priori} very simple yet surprising observation allows a deeper understanding of existing characterizations and paves the way to an immediate extension to a general density $x\mapsto p(x)$ by replacing $-x$ in these results with $(\log p(x))'$.


title: Pulse reverse-engineering for strong field-matter interaction

abstract: We propose a scheme to control the evolution of a two-level quantum system in the strong coupling regime based on the idea of reverse-engineering. A coherent control field is designed to drive both closed and open two-level quantum systems along user predefined evolution trajectory without utilizing the rotating-wave approximation (RWA). As concrete examples, we show that complete population inversion, an equally weighted coherent superposition, and even oscillationlike dynamics can be achieved. As there are no limitations on the coupling strength between the control field and matter, the scheme is attractive for applications such as accelerating desired system dynamics and fast quantum information processing.


title: Structure and fluctuations of the premelted liquid film of ice at the triple point

abstract: In this paper we study the structure of the ice/vapor interface in the neighborhood of the triple point for the TIP4P/2005 model. We probe the fluctuations of the ice/film and film/vapor surfaces that separate the liquid film from the coexisting bulk phases at basal, primary prismatic and secondary prismatic planes. The results are interpreted using a coupled sine Gordon plus Interface Hamiltonian model. At large length-scales, the two bounding surfaces are correlated and behave as a single complex ice/vapor interface. For small length, on the contrary, the ice/film and film/vapor surfaces behave very much like independent ice/water and water/vapor interfaces. The study suggests that the basal facet of the TIP4P/2005 model is smooth, the prismatic facet is close to a roughening transition, and the secondary prismatic facet is rough. For the faceted basal face, our fluctuation analysis allows us to estimate the step free energy in good agreement with experiment. Our results allow for a quantitative characterization of the extent to which the adsorbed quasi-liquid layer behaves as water, and explains experimental observations which reveal similar activation energies for crystals grown in bulk vapor or bulk water.

categories: physics.chem-ph cond-mat.soft cond-mat.stat-mech physics.ao-ph


title: Models and Integral Differentials of Hyperelliptic Curves

abstract: Let $C: y^2=f(x)$ be a hyperelliptic curve of genus $g\geq 1$, defined over a complete discretely valued field $K$, with ring of integers $O_K$. Under certain conditions on $C$, mild when residue characteristic is not $2$, we explicitly construct the minimal regular model with normal crossings $\mathcal{C}/O_K$ of $C$. In the same setting we determine a basis of integral differentials of $C$, that is an $O_K$-basis for the global sections of the relative dualising sheaf $\omega_{\mathcal{C}/O_K}$.


title: QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion Quantum Monte Carlo

abstract: We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.


title: Liquid crystal nose based on chiral photonic band gap materials: principles of selective response

abstract: Novel liquid crystalline (LC) compositions are suggested and studied as elements of LC-nose. This allows for optical detection of several volatile organic compounds (VOCs). Ethanol, toluene, pyridine and acetic acid were detected by means of colorimetric and spectroscopic techniques during their diffusion inside chiral elements of LC-nose. Selectivity to different VOCs is enhanced by means of components of liquid crystal matrix with different viscosity, affinities to the solvents, and abilities to form hydrogen bonding.


title: Treecode-accelerated Green Iteration for Kohn-Sham Density Functional Theory

abstract: We present a real-space computational method called treecode-accelerated Green Iteration (TAGI) for all-electron Kohn-Sham Density Functional Theory. TAGI is based on a reformulation of the Kohn-Sham equations in which the eigenvalue problem in differential form is converted into a fixed-point problem in integral form by convolution with the modified Helmholtz Green's function. In each self-consistent field (SCF) iteration, the fixed-points are computed by Green Iteration, where the discrete convolution sums are efficiently evaluated by a GPU-accelerated barycentric Lagrange treecode. Other techniques used in TAGI include adaptive mesh refinement, Fej\'er quadrature, singularity subtraction, gradient-free eigenvalue update, and Anderson mixing to accelerate convergence of the SCF and Green Iterations. Ground state energy computations of several atoms (Li, Be, O) and small molecules (H$_2$, CO, C$_6$H$_6$) demonstrate TAGI's ability to efficiently achieve chemical accuracy.


title: Acoustic diamond resonators with ultra-small mode volumes

abstract: Quantum acoustodynamics (QAD) is a rapidly developing field of research, offering possibilities to realize and study macroscopic quantum-mechanical systems in a new range of frequencies, and implement transducers and new types of memories for hybrid quantum devices. Here we propose a novel design for a versatile diamond QAD cavity operating at GHz frequencies, exhibiting effective mode volumes of about $10^{-4}\lambda^3$. Our phononic crystal waveguide cavity implements a non-resonant analogue of the optical lightning-rod effect to localize the energy of an acoustic mode into a deeply-subwavelength volume. We demonstrate that this confinement can readily enhance the orbit-strain interaction with embedded nitrogen-vacancy (NV) centres towards the high-cooperativity regime, and enable efficient resonant cooling of the acoustic vibrations towards the ground state using a single NV. This architecture can be readily translated towards setup with multiple cavities in one- or two-dimensional phononic crystals, and the underlying non-resonant localization mechanism will pave the way to further enhance optoacoustic coupling in phoxonic crystal cavities.


title: Learning Rope Manipulation Policies Using Dense Object Descriptors Trained on Synthetic Depth Data

abstract: Robotic manipulation of deformable 1D objects such as ropes, cables, and hoses is challenging due to the lack of high-fidelity analytic models and large configuration spaces. Furthermore, learning end-to-end manipulation policies directly from images and physical interaction requires significant time on a robot and can fail to generalize across tasks. We address these challenges using interpretable deep visual representations for rope, extending recent work on dense object descriptors for robot manipulation. This facilitates the design of interpretable and transferable geometric policies built on top of the learned representations, decoupling visual reasoning and control. We present an approach that learns point-pair correspondences between initial and goal rope configurations, which implicitly encodes geometric structure, entirely in simulation from synthetic depth images. We demonstrate that the learned representation -- dense depth object descriptors (DDODs) -- can be used to manipulate a real rope into a variety of different arrangements either by learning from demonstrations or using interpretable geometric policies. In 50 trials of a knot-tying task with the ABB YuMi Robot, the system achieves a 66% knot-tying success rate from previously unseen configurations. See https://tinyurl.com/rope-learning for supplementary material and videos.


title: A GPU-Accelerated Barycentric Lagrange Treecode

abstract: We present an MPI + OpenACC implementation of the kernel-independent barycentric Lagrange treecode (BLTC) for fast summation of particle interactions on GPUs. The distributed memory parallelization uses recursive coordinate bisection for domain decomposition and MPI remote memory access to build locally essential trees on each rank. The particle interactions are organized into target batch/source cluster interactions which efficiently map onto the GPU; target batching provides an outer level of parallelism, while the direct sum form of the barycentric particle-cluster approximation provides an inner level of parallelism. The GPU-accelerated BLTC performance is demonstrated on several test cases up to 1~billion particles interacting via the Coulomb potential and Yukawa potential.


title: Computing Lipschitz Constants for Hydraulic Models of Water Distribution Networks

abstract: Drinking water distribution networks (WDN) are large-scale, dynamic systems spanning large geographic areas. Water networks include various components such as junctions, reservoirs, tanks, pipes, pumps, and valves. Hydraulic models for these components depicting mass and energy balance form nonlinear algebraic differential equations (NDAE). While control theoretic studies have been thoroughly explored for other complex infrastructure such as power and transportation systems, little is understood or even investigated for feedback control and state estimation problems for the NDAE models of WDN. The objective of this letter is to showcase a complete NDAE model of WDN followed by computing Lipschitz constants of the vector-valued nonlinearity in that model. The computation of Lipschitz constants of hydraulic models is crucial as it paves the way to apply a plethora of control-theoretic studies for water system applications. In particular, the computation of Lipschitz constant is explored through closed-form, analytical expressions as well as via numerical methods. Case studies reveal how such computations fare against each other for various water networks.


title: Effect of receiver orientation on resource allocation in optical wireless systems

abstract: Optical wireless communication (OWC) systems have been the subject of a significant amount of interest as they can be used in sixth generation (6G) wireless communication to provide high data rates and support multiple users simultaneously. This paper investigates the impact of receiver orientation on resource allocation in optical wireless systems, using a wavelength division multiple access (WDMA) scheme. Three different systems that have different receiver orientations are examined in this work. Each of these systems considers 8 simultaneous users in two scenarios. WDMA is utilised to support multiple users and is based on four wavelengths offered by Red, Yellow, Green and Blue (RYGB) LDs for each AP. An angle diversity receiver (ADR) is used in each system with different orientations. The optimised resource allocations in terms of wavelengths and access point (AP) is obtained by using a mixed-integer linear programming (MILP) model. The channel bandwidth and SINR are determined in the two scenarios in all systems. The results show that a change in the orientation of the receiver can affect the level of channel bandwidth and SINR. However, SINRs in both scenarios for all users are above the threshold (15.6 dB). The SINR obtained can support t data rate of 5.7 Gbps in both scenarios in all systems.


title: Alfnoor: A Retrieval Simulation of the Ariel Target List

abstract: In this work, we present Alfnoor, a dedicated tool optimised for population studies of exoplanet atmospheres. Alfnoor combines the latest version of the retrieval algorithm TauREx 3, with the instrument noise simulator ArielRad and enables the simultaneous retrieval analysis of a large sample of exo-atmospheres. We applied this tool to the Ariel list of planetary candidates and focus on hydrogen dominated, cloudy atmospheres observed in transit with the Tier-2 mode (medium Ariel resolution). As a first experiment, we randomised the abundances - ranging from 10$^{-7}$ to 10$^{-2}$ - of the trace gases, which include H$_2$O, CH$_4$, CO, CO$_2$ and NH$_3$. This exercise allowed to estimate the detection limits for Ariel Tier-2 and Tier-3 modes when clouds are present. In a second experiment, we imposed an arbitrary trend between a chemical species and the effective temperature of the planet. A last experiment was run requiring molecular abundances being dictated by equilibrium chemistry at a certain temperature. Our results demonstrate the ability of Ariel Tier-2 and Tier-3 surveys to reveal trends between the chemistry and associated planetary parameters. Future work will focus on eclipse data, on atmospheres heavier than hydrogen and will be applied also to other observatories.


title: Deciphering the 3-D Orion Nebula-I: Expanding Shells in the Huygens Region

abstract: Based on imaging and spectroscopic data, we develop a 3-D model for the Huygens Region of the Orion Nebula. Theta1OriC , the hottest star in the Trapezium, is surrounded by a wind-blown Central Bubble that opens SW into the Extended Orion Nebula. Outside of this feature lies a layer of ionized gas at about 0.4 pc from Theta1OriC. Both of these features are moving rapidly away from Theta1OriC with an expansion age for the Central Bubble of only 15,000 yrs.


title: IsoRAN: Isolation and Scaling for 5G RANvia User-Level Data Plane Virtualization

abstract: 5G presents a unique set of challenges for cellular network architecture. The architecture needs to be versatile in order to handle a variety of use cases. While network slicing has been proposed as a way to provide such versatility, it is also important to ensure that slices do not adversely interfere with each other. In other words, isolation among network slices is needed. Additionally, the large number of use cases also implies a large number of users, making it imperative that 5G architectures scale efficiently. In this paper we propose IsoRAN, which provides isolation and scaling along with the flexibility needed for 5G architecture. In IsoRAN, users are processed by daemon threads in the Cloud Radio Access Network (CRAN) architecture. Our design allows users from different use cases to be executed, in a distributed manner, on the most efficient hardware to ensure that the Service Level Agreements (SLAs) are met while minimising power consumption. Our experiments show that IsoRAN handles users with different SLA while providing isolation to reduce interference. This increased isolation reduces the drop rate for different users from 42% to nearly 0% in some cases. Finally, we run large scale simulations on real traces to show the benefits for power consumption and cost reduction scale while increasing the number of base stations.


title: Characterisation of polyhedral products with finite generalised Postnikov decomposition

abstract: A generalised Postnikov tower for a space $X$ is a tower of principal fibrations with fibres generalised Eilenberg-MacLane spaces, whose inverse limit is weakly homotopy equivalent to $X$. In this paper we give a characterisation of a polyhedral product $Z_K(X,A)$ whose universal cover either admits a generalised Postnikov tower of finite length, or is a homotopy retract of a space admitting such a tower. We also include $p$-local and rational versions of the theorem. We end with a group theoretic application.


title: ${\cal N}=1$ conformal duals of gauged $E_n$ MN models

abstract: We suggest three new ${\cal N}=1$ conformal dual pairs. First, we argue that the ${\cal N}=2$ $E_6$ Minahan-Nemeschansky (MN) theory with a $USp(4)$ subgroup of the $E_6$ global symmetry conformally gauged with an ${\cal N}=1$ vector multiplet and certain additional chiral multiplet matter resides at some cusp of the conformal manifold of an $SU(2)^5$ quiver gauge theory. Second, we argue that the ${\cal N}=2$ $E_7$ MN theory with an $SU(2)$ subgroup of the $E_7$ global symmetry conformally gauged with an ${\cal N}=1$ vector multiplet and certain additional chiral multiplet matter resides at some cusp of the conformal manifold of a conformal ${\cal N}=1$ $USp(4)$ gauge theory. Finally, we claim that the ${\cal N}=2$ $E_8$ MN theory with a $USp(4)$ subgroup of the $E_8$ global symmetry conformally gauged with an ${\cal N}=1$ vector multiplet and certain additional chiral multiplet matter resides at some cusp of the conformal manifold of an ${\cal N}=1$ $Spin(7)$ conformal gauge theory. We argue for the dualities using a variety of non-perturbative techniques including anomaly and index computations. The dualities can be viewed as ${\cal N}=1$ analogues of ${\cal N}=2$ Argyres-Seiberg/Argyres-Wittig duals of the $E_n$ MN models. We also briefly comment on an ${\cal N}=1$ version of the Schur limit of the superconformal index.


title: Kashiwara-Vergne and dihedral bigraded Lie algebras in mould theory

abstract: We introduce the Kashiwara-Vergne bigraded Lie algebra associated with a finite abelian group and give its mould theoretic reformulation. By using the mould theory, we show that it includes Goncharov's dihedral Lie algebra, which generalizes the result of Raphael and Schneps.


title: Molecule formation in dust-poor irradiated jets I. Stationary disk winds

abstract: Recent ALMA observations suggest that the highest velocity part of molecular protostellar jets are launched from the dust-sublimation regions of the accretion disks (<0.3 au). However, formation and survival of molecules in inner protostellar disk winds, in the presence of a harsh FUV radiation field and the absence of dust, remain unexplored. We aim at determining if simple molecules can be synthesized and spared in fast and collimated dust-free disk winds or if a fraction of dust is necessary to explain the observed molecular abundances. This work is based on the Paris-Durham shock code designed to model irradiated environments. Fundamental properties of the dust-free chemistry are investigated from single point models. A laminar 1D disk wind model is then built using a parametric flow geometry. This model includes time-dependent chemistry and the attenuation of the radiation field by gas-phase photoprocesses. We show that a small fraction of H2 (< 1e-2), primarily formed through the H- route, can efficiently initiate molecule synthesis such as CO and SiO above TK ~ 800 K. The attenuation of the radiation field by atomic species (eg. C, Si, S) proceeds through continuum self-shielding. This process ensures efficient formation of CO, OH, SiO, H2O through neutral-neutral reactions, and the survival of these molecules. Class 0 dust-free winds with high mass-loss rates ($\dot{M}_w >$ 2e-6 Msun/yr) are predicted to be rich in molecules if warm (TK > 800 K). The molecular content of disk winds is very sensitive to the presence of dust and a mass-fraction of surviving dust as small as 1e-5 significantly increases the H2O and SiO abundances. Chemistry of high-velocity jets is a powerful tool to probe their content in dust and uncover their launching point. Models of internal shocks are required to fully exploit the current (sub-)millimeter observations and prepare future JWST observations.


title: New Examples of Minimal Non-Strongly-Perfect Graphs

abstract: A graph is strongly perfect if every induced subgraph H has a stable set that meets every nonempty maximal clique of H. The characterization of strongly perfect graphs by a set of forbidden induced subgraphs is not known. Here we provide several new minimal non-strongly-perfect graphs.


title: Generalized Gumbel-Softmax Gradient Estimator for Generic Discrete Random Variables

abstract: Estimating the gradients of stochastic nodes in stochastic computational graphs is one of the crucial research questions in the deep generative modeling community, which enables the gradient descent optimization on neural network parameters. Stochastic gradient estimators of discrete random variables are widely explored, for example, Gumbel-Softmax reparameterization trick for Bernoulli and categorical distributions. Meanwhile, other discrete distribution cases such as the Poisson, geometric, binomial, multinomial, negative binomial, etc. have not been explored. This paper proposes a generalized version of the Gumbel-Softmax estimator, which is able to reparameterize generic discrete distributions, not restricted to the Bernoulli and the categorical. The proposed estimator utilizes the truncation of discrete random variables, the Gumbel-Softmax trick, and a special form of linear transformation. Our experiments consist of (1) synthetic examples and applications on VAE, which show the efficacy of our methods; and (2) topic models, which demonstrate the value of the proposed estimation in practice.


title: On Emergent Communication in Competitive Multi-Agent Teams

abstract: Several recent works have found the emergence of grounded compositional language in the communication protocols developed by mostly cooperative multi-agent systems when learned end-to-end to maximize performance on a downstream task. However, human populations learn to solve complex tasks involving communicative behaviors not only in fully cooperative settings but also in scenarios where competition acts as an additional external pressure for improvement. In this work, we investigate whether competition for performance from an external, similar agent team could act as a social influence that encourages multi-agent populations to develop better communication protocols for improved performance, compositionality, and convergence speed. We start from Task & Talk, a previously proposed referential game between two cooperative agents as our testbed and extend it into Task, Talk & Compete, a game involving two competitive teams each consisting of two aforementioned cooperative agents. Using this new setting, we provide an empirical study demonstrating the impact of competitive influence on multi-agent teams. Our results show that an external competitive influence leads to improved accuracy and generalization, as well as faster emergence of communicative languages that are more informative and compositional.


title: Distributed velocity-constrained consensus of discrete-time multi-agent systems with nonconvex constraints, switching topologies, and delays

abstract: In this paper, a distributed velocity-constrained consensus problem is studied for discrete-time multi-agent systems, where each agent's velocity is constrained to lie in a nonconvex set. A distributed constrained control algorithm is proposed to enable all agents to converge to a common point using only local information. {The gains of the algorithm for all agents need not to be the same or predesigned and can be adjusted by each agent itself based on its own and neighbors' information.} It is shown that the algorithm is robust to arbitrarily bounded communication delays and arbitrarily switching communication graphs provided that the union of the graphs has directed spanning trees among each certain time interval. The analysis approach is based on multiple novel model transformations, proper control parameter selections, boundedness analysis of state-dependent stochastic matrices, exploitation of the convexity of stochastic matrices, and the joint connectivity of the communication graphs. Numerical examples are included to illustrate the theoretical results.


title: Risk Projection for Time-to-event Outcome Leveraging Summary Statistics With Source Individual-level Data

abstract: Predicting risks of chronic diseases has become increasingly important in clinical practice. When a prediction model is developed in a given source cohort, there is often a great interest to apply the model to other cohorts. However, due to potential discrepancy in baseline disease incidences between different cohorts and shifts in patient composition, the risk predicted by the original model often under- or over-estimates the risk in the new cohort. The remedy of such a poorly calibrated prediction is needed for proper medical decision-making. In this article, we assume the relative risks of predictors are the same between the two cohorts, and propose a novel weighted estimating equation approach to re-calibrating the projected risk for the targeted population through updating the baseline risk. The recalibration leverages the knowledge about the overall survival probabilities for the disease of interest and competing events, and the summary information of risk factors from the targeted population. The proposed re-calibrated risk estimators gain efficiency if the risk factor distributions are the same for both the source and target cohorts, and are robust with little bias if they differ. We establish the consistency and asymptotic normality of the proposed estimators. Extensive simulation studies demonstrate that the proposed estimators perform very well in terms of robustness and efficiency in finite samples. A real data application to colorectal cancer risk prediction also illustrates that the proposed method can be used in practice for model recalibration.


title: A Distributed Pipeline for Scalable, Deconflicted Formation Flying

abstract: Reliance on external localization infrastructure and centralized coordination are main limiting factors for formation flying of vehicles in large numbers and in unprepared environments. While solutions using onboard localization address the dependency on external infrastructure, the associated coordination strategies typically lack collision avoidance and scalability. To address these shortcomings, we present a unified pipeline with onboard localization and a distributed, collision-free motion planning strategy that scales to a large number of vehicles. Since distributed collision avoidance strategies are known to result in gridlock, we also present a decentralized task assignment solution to deconflict vehicles. We experimentally validate our pipeline in simulation and hardware. The results show that our approach for solving the optimization problem associated with motion planning gives solutions within seconds in cases where general purpose solvers fail due to high complexity. In addition, our lightweight assignment strategy leads to successful and quicker formation convergence in 96-100% of all trials, whereas indefinite gridlocks occur without it for 33-50% of trials. By enabling large-scale, deconflicted coordination, this pipeline should help pave the way for anytime, anywhere deployment of aerial swarms.


title: q-VAE for Disentangled Representation Learning and Latent Dynamical Systems

abstract: A variational autoencoder (VAE) derived from Tsallis statistics called q-VAE is proposed. In the proposed method, a standard VAE is employed to statistically extract latent space hidden in sampled data, and this latent space helps make robots controllable in feasible computational time and cost. To improve the usefulness of the latent space, this paper focuses on disentangled representation learning, e.g., $\beta$-VAE, which is the baseline for it. Starting from a Tsallis statistics perspective, a new lower bound for the proposed q-VAE is derived to maximize the likelihood of the sampled data, which can be considered an adaptive $\beta$-VAE with deformed Kullback-Leibler divergence. To verify the benefits of the proposed q-VAE, a benchmark task to extract the latent space from the MNIST dataset was performed. The results demonstrate that the proposed q-VAE improved disentangled representation while maintaining the reconstruction accuracy of the data. In addition, it relaxes the independency condition between data, which is demonstrated by learning the latent dynamics of nonlinear dynamical systems. By combining disentangled representation, the proposed q-VAE achieves stable and accurate long-term state prediction from the initial state and the action sequence. The dataset for hexapod walking is available on IEEE Dataport, doi: https://dx.doi.org/10.21227/99af-jw71.


title: Hypergraph Motifs: Concepts, Algorithms, and Discoveries

abstract: Hypergraphs naturally represent group interactions, which are omnipresent in many domains: collaborations of researchers, co-purchases of items, joint interactions of proteins, to name a few. In this work, we propose tools for answering the following questions in a systematic manner: (Q1) what are structural design principles of real-world hypergraphs? (Q2) how can we compare local structures of hypergraphs of different sizes? (Q3) how can we identify domains which hypergraphs are from? We first define hypergraph motifs (h-motifs), which describe the connectivity patterns of three connected hyperedges. Then, we define the significance of each h-motif in a hypergraph as its occurrences relative to those in properly randomized hypergraphs. Lastly, we define the characteristic profile (CP) as the vector of the normalized significance of every h-motif. Regarding Q1, we find that h-motifs' occurrences in 11 real-world hypergraphs from 5 domains are clearly distinguished from those of randomized hypergraphs. In addition, we demonstrate that CPs capture local structural patterns unique to each domain, and thus comparing CPs of hypergraphs addresses Q2 and Q3. Our algorithmic contribution is to propose MoCHy, a family of parallel algorithms for counting h-motifs' occurrences in a hypergraph. We theoretically analyze their speed and accuracy, and we show empirically that the advanced approximate version MoCHy-A+ is up to 25X more accurate and 32X faster than the basic approximate and exact versions, respectively.


title: Chimera States in Hybrid Coupled Neuron Populations

abstract: Here we study the emergence of chimera states, a recently reported phenomenon referring to the coexistence of synchronized and unsynchronized dynamical units, in a population of Morris-Lecar neurons which are coupled by both electrical and chemical synapses, constituting a hybrid synaptic architecture, as in actual brain connectivity. This scheme consists of a nonlocal network where the nearest neighbor neurons are coupled by electrical synapses, while the synapses from more distant neurons are of the chemical type. We demonstrate that peculiar dynamical behaviors, including chimera state and traveling wave, exist in such a hybrid coupled neural system, and analyze how the relative abundance of chemical and electrical synapses affects the features of chimera and different synchrony states (i.e. incoherent, traveling wave and coherent) and the regions in the space of relevant parameters for their emergence. Additionally, we show that, when the relative population of chemical synapses increases further, a new intriguing chaotic dynamical behavior appears above the region for chimera states. This is characterized by the coexistence of two distinct synchronized states with different amplitude, and an unsynchronized state, that we denote as a chaotic amplitude chimera. We also discuss about the computational implications of such state.


title: Equity-Based Incentives, Production/Service Functions And Game Theory

abstract: EBIs/ESOs substantially change the traditional production/service function because ESOs/EBIs can have different psychological effects(motivation or de-motivation), and can create intangible capital and different economic payoffs. Although Game Theory is flawed, it can be helpful in describing interactions in ESO/EBIs transactions. ESOs/EBIs involve two-stage games and there are no perfect Nash Equilibria for the two sub-games. The large number of actual and potential participants in these games significantly complicates resolution of equilibria and increases the dynamism of the games given that players are more sensitive to other peoples moves in such games. This article: a) analyzes how ESOs/EBIs affect traditional assumptions of production functions (in both the manufacturing and service sectors), b) analyzes ESOs/EBIs transactions using game theory concepts, c) illustrates some of the limitations of game theory.


title: Universal sieve-based strategies for efficient estimation using machine learning tools

abstract: Suppose that we wish to estimate a finite-dimensional summary of one or more function-valued features of an underlying data-generating mechanism under a nonparametric model. One approach to estimation is by plugging in flexible estimates of these features. Unfortunately, in general, such estimators may not be asymptotically efficient, which often makes these estimators difficult to use as a basis for inference. Though there are several existing methods to construct asymptotically efficient plug-in estimators, each such method either can only be derived using knowledge of efficiency theory or is only valid under stringent smoothness assumptions. Among existing methods, sieve estimators stand out as particularly convenient because efficiency theory is not required in their construction, their tuning parameters can be selected data adaptively, and they are universal in the sense that the same fits lead to efficient plug-in estimators for a rich class of estimands. Inspired by these desirable properties, we propose two novel universal approaches for estimating function-valued features that can be analyzed using sieve estimation theory. Compared to traditional sieve estimators, these approaches are valid under more general conditions on the smoothness of the function-valued features by utilizing flexible estimates that can be obtained, for example, using machine learning.


title: SeMemNN: A Semantic Matrix-Based Memory Neural Network for Text Classification

abstract: Text categorization is the task of assigning labels to documents written in a natural language, and it has numerous real-world applications including sentiment analysis as well as traditional topic assignment tasks. In this paper, we propose 5 different configurations for the semantic matrix-based memory neural network with end-to-end learning manner and evaluate our proposed method on two corpora of news articles (AG news, Sogou news). The best performance of our proposed method outperforms the baseline VDCNN models on the text classification task and gives a faster speed for learning semantics. Moreover, we also evaluate our model on small scale datasets. The results show that our proposed method can still achieve better results in comparison to VDCNN on the small scale dataset. This paper is to appear in the Proceedings of the 2020 IEEE 14th International Conference on Semantic Computing (ICSC 2020), San Diego, California, 2020.


title: Time-frequency Analysis of two-wavelet theory in Weinstein setting

abstract: In this paper, we introduce the notion of Weinstein two-wavelet and we define the two-wavelet localization operators in the setting of the Weinstein theory. Then we give a host of sufficient conditions for the boundedness and compactness of the two-wavelet localization operator on $L^{p}_{\alpha}(\mathbb{R}^{d+1}_+)$ for all $1\leq p\leq \infty$, in terms of properties of the symbol $\sigma$ and the functions $\varphi$ and $\psi$. In the end, we study some typical examples of the Weinstein two-wavelet localization operators.


title: Applications of deep learning in stock market prediction: recent progress

abstract: Stock market prediction has been a classical yet challenging problem, with the attention from both economists and computer scientists. With the purpose of building an effective prediction model, both linear and machine learning tools have been explored for the past couple of decades. Lately, deep learning models have been introduced as new frontiers for this topic and the rapid development is too fast to catch up. Hence, our motivation for this survey is to give a latest review of recent works on deep learning models for stock market prediction. We not only category the different data sources, various neural network structures, and common used evaluation metrics, but also the implementation and reproducibility. Our goal is to help the interested researchers to synchronize with the latest progress and also help them to easily reproduce the previous studies as baselines. Base on the summary, we also highlight some future research directions in this topic.


title: Designing a Bonus-Malus system reflecting the claim size under the dependent frequency-severity model

abstract: In auto insurance, a Bonus-Malus System (BMS) is commonly used as a posteriori risk classification mechanism to set the premium for the next contract period based on a policyholder's claim history. Even though recent literature reports evidence of a significant dependence between frequency and severity, the current BMS practice is to use a frequency-based transition rule while ignoring severity information. Although Oh et al. (2019) claim that the frequency-driven BMS transition rule can accommodate the dependence between frequency and severity, their proposal is only a partial solution, as the transition rule still completely ignores the claim severity and is unable to penalize large claims. In this study, we propose to use the BMS with a transition rule based on both frequency and size of claim, based on the bivariate random effect model, which conveniently allows dependence between frequency and severity. We analytically derive the optimal relativities under the proposed BMS framework and show that the proposed BMS outperforms the existing frequency-driven BMS. Later numerical experiments are also provided using both hypothetical and actual datasets in order to assess the effect of various dependencies on the BMS risk classification and confirm our theoretical findings.


title: Multi-Microphone Complex Spectral Mapping for Speech Dereverberation

abstract: This study proposes a multi-microphone complex spectral mapping approach for speech dereverberation on a fixed array geometry. In the proposed approach, a deep neural network (DNN) is trained to predict the real and imaginary (RI) components of direct sound from the stacked reverberant (and noisy) RI components of multiple microphones. We also investigate the integration of multi-microphone complex spectral mapping with beamforming and post-filtering. Experimental results on multi-channel speech dereverberation demonstrate the effectiveness of the proposed approach.


title: An application benchmark for fermionic quantum simulations

abstract: It is expected that the simulation of correlated fermions in chemistry and material science will be one of the first practical applications of quantum processors. Given the rapid evolution of quantum hardware, it is increasingly important to develop robust benchmarking techniques to gauge the capacity of quantum hardware specifically for the purpose of fermionic simulation. Here we propose using the one-dimensional Fermi-Hubbard model as an application benchmark for variational quantum simulations on near-term quantum devices. Since the one-dimensional Hubbard model is both strongly correlated and exactly solvable with the Bethe ansatz, it provides a reference ground state energy that a given device with limited coherence will be able to approximate up to a maximal size. The length of the largest chain that can be simulated provides an effective fermionic length. We use variational quantum eigensolver to approximate the ground state energy values of Fermi-Hubbard instances and show how the fermionic length benchmark can be used in practice to assess the performance of bounded-depth devices in a scalable fashion.


title: Hyperbolic 3-manifolds with large kissing number

abstract: In this article we construct a sequence $\{M_i\}$ of non compact finite volume hyperbolic $3$-manifolds whose kissing number grows at least as $\mathrm{vol}(M_i)^{\frac{31}{27}-\epsilon}$ for any $\epsilon>0$. This extends a previous result due to Schmutz in dimension $2$.


title: Visualizing and Understanding Large-Scale Assessments in Mathematics through Dimensionality Reduction

abstract: In this paper, we apply the Logistic PCA (LPCA) as a dimensionality reduction tool for visualizing patterns and characterizing the relevance of mathematics abilities from a given population measured by a large-scale assessment. We establish an equivalence of parameters between LPCA, Inner Product Representation (IPR) and the two paramenter logistic model (2PL) from the Item Response Theory (IRT). This equivalence provides three complemetary ways of looking at data that assists professionals in education to perform in-context interpretations. Particularly, we analyse the data collected from SPAECE, a large-scale assessment in Mathematics that has been applied yearly in the public educational system of the state of Cear\'a, Brazil. As the main result, we show that the the poor performance of examinees in the end of middle school is primarily caused by their disabilities in number sense.


title: New dynamics of the classical and nonlocal Gross-Pitaevskii equation with a parabolic potential

abstract: Solutions of the classical and nonlocal Gross-Pitaevskii (GP) equation with a parabolic potential and a gain term are derived by using a second order nonisospectral Ablowitz-Kaup-Newell-Segur system and reduction technique of double Wronskians. Solutions of the classical GP equation show typical space-time localized characteristics. An interesting dynamics, solitons carrying an oscillating wave, are found with mathematical analysis and illustrations. Solutions of some nonlocal cases are also illustrated.


title: Region adaptive graph fourier transform for 3d point clouds

abstract: We introduce the Region Adaptive Graph Fourier Transform (RA-GFT) for compression of 3D point cloud attributes. The RA-GFT is a multiresolution transform, formed by combining spatially localized block transforms. We assume the points are organized by a family of nested partitions represented by a rooted tree. At each resolution level, attributes are processed in clusters using block transforms. Each block transform produces a single approximation (DC) coefficient, and various detail (AC) coefficients. The DC coefficients are promoted up the tree to the next (lower resolution) level, where the process can be repeated until reaching the root. Since clusters may have a different numbers of points, each block transform must incorporate the relative importance of each coefficient. For this, we introduce the $\mathbf{Q}$-normalized graph Laplacian, and propose using its eigenvectors as the block transform. The RA-GFT achieves better complexity-performance trade-offs than previous approaches. In particular, it outperforms the Region Adaptive Haar Transform (RAHT) by up to 2.5 dB, with a small complexity overhead.


title: Probing single stop production at the FCC-hh/SPPC

abstract: Top squark (stop) is a crucial part of supersymmetric models (SUSY) to understand the naturalness problem. Other than the traditional stop pair production, the single production via electroweak interaction provides signals with distinctive features which could help confirm the existence of the top squark. In this paper, we investigate the observability of stop through the mono-top channel of the single stop production at the future proton-proton colliders, FCC-hh and SPPC, in a simplified Minimal Supersymmetric Standard Model (MSSM). With the integrated luminosity of 3000 $\text{fb}^{-1}$, we can probe the stop with mass up to 3.25 TeV by the mono-top channel at $5\sigma$ level. Considering the systematic uncertainty of 10%, the exclusion limit for stop mass can be reached at about 1.5 TeV. Exclusion limits on stop mass and higgsino mass parameter $\mu$ are also presented.


title: Robust Instability Analysis with Application to Neuronal Dynamics

abstract: This paper is concerned with robust instability analysis of linear feedback systems subject to a dynamic uncertainty. The work is motivated by, and provides a basic foundation for, a more challenging problem of analyzing persistence of oscillations in nonlinear dynamical systems. We first formalize the problem for SISO LTI systems by introducing a notion of the robust instability radius (RIR). We provide a method for calculating the RIR exactly for a certain class of systems and show that it works well for a class of second order systems. This result is applied to the FitzHugh-Nagumo model for neuronal dynamics, and the effectiveness is confirmed by numerical simulations, where we properly care for the change of the equilibrium point.


title: Exactly solvable two-terminal heat engine with asymmetric Onsager coefficients: Origin of the power-efficiency bound

abstract: An engine producing a finite power at the ideal (Carnot) efficiency is a dream engine, which is not prohibited by the thermodynamic second law. Some years ago, a two-terminal heat engine with {\em asymmetric} Onsager coefficients in the linear response regime was suggested by Benenti, Saito, and Casati [Phys. Rev. Lett. {\bf 106}, 230602 (2011)], as a prototypical system to make such a dream come true with non-divergent system parameter values. However, such a system has never been realized in spite of many trials. Here, we introduce an exactly solvable two-terminal Brownian heat engine with the asymmetric Onsager coefficients in the presence of a Lorenz (magnetic) force. Nevertheless, we show that the dream engine regime cannot be accessible even with the asymmetric Onsager coefficients, due to an instability keeping the engine from reaching its steady state. This is consistent with recent trade-off relations between the engine power and efficiency, where the (cyclic) steady-state condition is implicitly presumed. We conclude that the inaccessibility to the dream engine originates from the steady-state constraint on the engine.


title: Weeping and Gnashing of Teeth: Teaching Deep Learning in Image and Video Processing Classes

abstract: In this rather informal paper and talk I will discuss my own experiences, feelings, and evolution as an Image Processing and Digital Video educator trying to navigate the Deep Learning revolution. I will discuss my own ups and downs of trying to deal with extremely rapid technological changes, and how I have reacted to, and dealt with consequent dramatic changes in the relevance of the topics I've taught for three decades. I have arranged the discussion in terms of the stages, over time, of my progression dealing with these sea changes.


title: Semantic sensor fusion: from camera to sparse lidar information

abstract: To navigate through urban roads, an automated vehicle must be able to perceive and recognize objects in a three-dimensional environment. A high-level contextual understanding of the surroundings is necessary to plan and execute accurate driving maneuvers. This paper presents an approach to fuse different sensory information, Light Detection and Ranging (lidar) scans and camera images. The output of a convolutional neural network (CNN) is used as classifier to obtain the labels of the environment. The transference of semantic information between the labelled image and the lidar point cloud is performed in four steps: initially, we use heuristic methods to associate probabilities to all the semantic classes contained in the labelled images. Then, the lidar points are corrected to compensate for the vehicle's motion given the difference between the timestamps of each lidar scan and camera image. In a third step, we calculate the pixel coordinate for the corresponding camera image. In the last step we perform the transfer of semantic information from the heuristic probability images to the lidar frame, while removing the lidar information that is not visible to the camera. We tested our approach in the Usyd Dataset \cite{usyd_dataset}, obtaining qualitative and quantitative results that demonstrate the validity of our probabilistic sensory fusion approach.


title: Type I Attack for Generative Models

abstract: Generative models are popular tools with a wide range of applications. Nevertheless, it is as vulnerable to adversarial samples as classifiers. The existing attack methods mainly focus on generating adversarial examples by adding imperceptible perturbations to input, which leads to wrong result. However, we focus on another aspect of attack, i.e., cheating models by significant changes. The former induces Type II error and the latter causes Type I error. In this paper, we propose Type I attack to generative models such as VAE and GAN. One example given in VAE is that we can change an original image significantly to a meaningless one but their reconstruction results are similar. To implement the Type I attack, we destroy the original one by increasing the distance in input space while keeping the output similar because different inputs may correspond to similar features for the property of deep neural network. Experimental results show that our attack method is effective to generate Type I adversarial examples for generative models on large-scale image datasets.


title: Large-Scale Shrinkage Estimation under Markovian Dependence

abstract: We consider the problem of simultaneous estimation of a sequence of dependent parameters that are generated from a hidden Markov model. Based on observing a noise contaminated vector of observations from such a sequence model, we consider simultaneous estimation of all the parameters irrespective of their hidden states under square error loss. We study the roles of statistical shrinkage for improved estimation of these dependent parameters. Being completely agnostic on the distributional properties of the unknown underlying Hidden Markov model, we develop a novel non-parametric shrinkage algorithm. Our proposed method elegantly combines \textit{Tweedie}-based non-parametric shrinkage ideas with efficient estimation of the hidden states under Markovian dependence. Based on extensive numerical experiments, we establish superior performance our our proposed algorithm compared to non-shrinkage based state-of-the-art parametric as well as non-parametric algorithms used in hidden Markov models. We provide decision theoretic properties of our methodology and exhibit its enhanced efficacy over popular shrinkage methods built under independence. We demonstrate the application of our methodology on real-world datasets for analyzing of temporally dependent social and economic indicators such as search trends and unemployment rates as well as estimating spatially dependent Copy Number Variations.


title: A Deep Learning Method for Complex Human Activity Recognition Using Virtual Wearable Sensors

abstract: Sensor-based human activity recognition (HAR) is now a research hotspot in multiple application areas. With the rise of smart wearable devices equipped with inertial measurement units (IMUs), researchers begin to utilize IMU data for HAR. By employing machine learning algorithms, early IMU-based research for HAR can achieve accurate classification results on traditional classical HAR datasets, containing only simple and repetitive daily activities. However, these datasets rarely display a rich diversity of information in real-scene. In this paper, we propose a novel method based on deep learning for complex HAR in the real-scene. Specially, in the off-line training stage, the AMASS dataset, containing abundant human poses and virtual IMU data, is innovatively adopted for enhancing the variety and diversity. Moreover, a deep convolutional neural network with an unsupervised penalty is proposed to automatically extract the features of AMASS and improve the robustness. In the on-line testing stage, by leveraging advantages of the transfer learning, we obtain the final result by fine-tuning the partial neural network (optimizing the parameters in the fully-connected layers) using the real IMU data. The experimental results show that the proposed method can surprisingly converge in a few iterations and achieve an accuracy of 91.15% on a real IMU dataset, demonstrating the efficiency and effectiveness of the proposed method.


title: Localising Faster: Efficient and precise lidar-based robot localisation in large-scale environments

abstract: This paper proposes a novel approach for global localisation of mobile robots in large-scale environments. Our method leverages learning-based localisation and filtering-based localisation, to localise the robot efficiently and precisely through seeding Monte Carlo Localisation (MCL) with a deep-learned distribution. In particular, a fast localisation system rapidly estimates the 6-DOF pose through a deep-probabilistic model (Gaussian Process Regression with a deep kernel), then a precise recursive estimator refines the estimated robot pose according to the geometric alignment. More importantly, the Gaussian method (i.e. deep probabilistic localisation) and non-Gaussian method (i.e. MCL) can be integrated naturally via importance sampling. Consequently, the two systems can be integrated seamlessly and mutually benefit from each other. To verify the proposed framework, we provide a case study in large-scale localisation with a 3D lidar sensor. Our experiments on the Michigan NCLT long-term dataset show that the proposed method is able to localise the robot in 1.94 s on average (median of 0.8 s) with precision 0.75~m in a large-scale environment of approximately 0.5 km2.


title: Privacy-preserving Learning via Deep Net Pruning

abstract: This paper attempts to answer the question whether neural network pruning can be used as a tool to achieve differential privacy without losing much data utility. As a first step towards understanding the relationship between neural network pruning and differential privacy, this paper proves that pruning a given layer of the neural network is equivalent to adding a certain amount of differentially private noise to its hidden-layer activations. The paper also presents experimental results to show the practical implications of the theoretical finding and the key parameter values in a simple practical setting. These results show that neural network pruning can be a more effective alternative to adding differentially private noise for neural networks.


title: Quasi-periodic pulsation detected in Lyman-alpha emission during solar flares

abstract: We investigated the quasi-periodic pulsation (QPP) in Lyman-alpha, X-ray and extreme-ultraviolet (EUV) emissions during two solar flares, i.e., an X-class (SOL2012-01-27T) and a C-class (SOL2016-02-08T). The full-disk Lyman-alpha and X-Ray flux during these solar flares were recorded by the EUV Sensor and X-Ray Sensor on board the Geostationary Operational Environmental Satellite. The {\deg}are regions were located from the EUV images measured by the Atmospheric Imaging Assembly. The QPP could be identified as a series of regular and periodic peaks in the light curves, and its quasi-periodicity was determined from the global wavelet and Fourier power spectra. A quasi-periodicity at about 3 minutes is detected during the impulsive phase of the X-class flare, which could be explained as the acoustic wave in the chromosphere (e.g., Milligan et al. 2017). Interestingly, a quasi-periodicity at roughly 1 minute is discovered during the entire evolutionary phases of solar flares, including the precursor, impulsive, and gradual phases. This is the first report of 1-minute QPP in the Lyman-alpha emission during solar flares, in particular during the flare precursor. It may be interpreted as a self-oscillatory regime of the magnetic reconnection, such as magnetic dripping.


title: Physics-informed machine learning for composition-process-property alloy design: shape memory alloy demonstration

abstract: Machine learning (ML) is shown to predict new alloys and their performances in a high dimensional, multiple-target-property design space that considers chemistry, multi-step processing routes, and characterization methodology variations. A physics-informed featured engineering approach is shown to enable otherwise poorly performing ML models to perform well with the same data. Specifically, previously engineered elemental features based on alloy chemistries are combined with newly engineered heat treatment process features. The new features result from first transforming the heat treatment parameter data as it was previously recorded using nonlinear mathematical relationships known to describe the thermodynamics and kinetics of phase transformations in alloys. The ability of the ML model to be used for predictive design is validated using blind predictions. Composition - process - property relationships for thermal hysteresis of shape memory alloys (SMAs) with complex microstructures created via multiple melting-homogenization-solutionization-precipitation processing stage variations are captured, in addition to the mean transformation temperatures of the SMAs. The quantitative models of hysteresis exhibited by such highly processed alloys demonstrate the ability for ML models to design for physical complexities that have challenged physics-based modeling approaches for decades.


title: Kinetic simulations of nonrelativistic perpendicular shocks of young supernova remnants. III. Magnetic reconnection

abstract: Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfv\'enic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.


title: Safeguarded Learned Convex Optimization

abstract: Applications abound in which optimization problems must be repeatedly solved, each time with new (but similar) data. Analytic optimization algorithms can be hand-designed to provably solve these problems in an iterative fashion. On one hand, data-driven algorithms can "learn to optimize" (L2O) with much fewer iterations and similar cost per iteration as general-purpose optimization algorithms. On the other hand, unfortunately, many L2O algorithms lack converge guarantees. To fuse the advantages of these approaches, we present a Safe-L2O framework. Safe-L2O updates incorporate a safeguard to guarantee convergence for convex problems with proximal and/or gradient oracles. The safeguard is simple and computationally cheap to implement, and it is activated only when the data-driven L2O updates would perform poorly or appear to diverge. This yields the numerical benefits of employing machine learning to create rapid L2O algorithms while still guaranteeing convergence. Our numerical examples show convergence of Safe-L2O algorithms, even when the provided data is not from the distribution of training data.


title: Modern cities emerge as 'super-cells' where enclosed industrial systems are hotspots of goods and services

abstract: Prevailing hypotheses recognize cities as 'super-organisms' which both provides organizing principles for cities and fills the scalar gap in the hierarchical living system between ecosystems and the entire planet. However, most analogies between the traits of organisms and cities are inappropriate making the super-organism model impractical as a means to acquire new knowledge. Using a cluster analysis of 15 traits of cities and other living systems, we found that modern cities are more similar to eukaryotic cells than to multicellular organisms. Enclosed industrial systems, such as factories and greenhouses, dominate modern cities and are analogous to organelles as hotspots that provide high-flux goods and services. Therefore, we propose a 'super-cell city model' as more appropriate than the super-organism model. In addition to the theoretical significance, our model also recognizes enclosed industrial systems as functional components that improve the vitality and sustainability of cities.


title: Constraining the Neutrino Mass with the Drifting Coefficient of the Field Cluster Mass Function

abstract: A new diagnostics to break the degeneracy between the total neutrino mass ($M_{\nu}$) and the primordial power spectrum amplitude ($\sigma_{8}$) by using the drifting coefficient of the field cluster mass function is presented. Analyzing the data from the Cosmological Massive Neutrino Simulations, we first determine the numerical mass functions of the field clusters at various redshifts. Then, we compare the numerical results with the analytical model characterized by a single parameter called the drifting coefficient which measures the drifts of the collapse density threshold, $\delta_{c}$, from the Einstein-de Sitter spherical value, $\delta_{sc}$, at a given mass scale. It is found that the analytic model for the field cluster mass function is found to work excellently even in the presence of massive neutrinos and that its drifting coefficient evolves differently in the cosmologies with different values of $M_{\nu}$. At low redshifts ($z\lesssim 0.3$) the more massive neutrinos drift $\delta_{c}$ further from $\delta_{sc}$, while the opposite trend is found at higher redshifts ($z\gtrsim0.3$). Speculating that this distinct redshift-dependent effect of massive neutrinos on the drifting coefficient of the field cluster mass function might help break the $\sigma_{8}$-$M_{\nu}$ degeneracy, we also show that the sensitivity of this new diagnostics to $M_{\nu}$ is high enough to discriminate the case of $M_{\nu}=0.1\,{\rm eV}$ from that of massless neutrinos.


title: Solitary waves in the resonant nonlinear Schr\"odinger equation: stability and dynamical properties

abstract: The stability and dynamical properties of the so-called resonant nonlinear Schr\"odinger (RNLS) equation, are considered. The RNLS is a variant of the nonlinear Schr\"odinger (NLS) equation with the addition of a perturbation used to describe wave propagation in cold collisionless plasmas. We first examine the modulational stability of plane waves in the RNLS model, identifying the modifications of the associated conditions from the NLS case. We then move to the study of solitary waves with vanishing and nonzero boundary conditions. Interestingly the RNLS, much like the usual NLS, exhibits both dark and bright soliton solutions depending on the relative signs of dispersion and nonlinearity. The corresponding existence, stability and dynamics of these solutions are studied systematically in this work.


title: Asymptotic behavior of branching diffusion processes in periodic media

abstract: We study the asymptotic behavior of branching diffusion processes in periodic media. For a super-critical branching process, we distinguish two types of behavior for the normalized number of particles in a bounded domain, depending on the distance of the domain from the region where the bulk of the particles is located. At distances that grow linearly in time, we observe intermittency (i.e., the $k$-th moment dominates the $k$-th power of the first moment for some $k$), while, at distances that grow sub-linearly in time, we show that all the moments converge. A key ingredient in our analysis is a sharp estimate of the transition kernel for the branching process, valid up to linear in time distances from the location of the initial particle.


title: No-go Theorem of Purification

abstract: The Shannon's bound for compression is one of the key restrictions for the compression of quantum information. Here we show that the unitarity of the compression operation imposes new bounds on the compression that are more limiting than Shannon's compression bound. This translates to a no-go theorem for the purification of quantum states. For a specific case of a two-qubit system, our results indicate that it is not possible to distill purity beyond the maximum of the individual purities. We show that this restriction results in the cooling limit of the heat-bath algorithmic cooling techniques. We formalize the limitations imposed by the unitarity of the compression operation in two theorems and use the theorems to show that the limitations of unitarity lead to the cooling limit of heat-bath algorithmic cooling. To this end, we introduce a new optimal cooling technique and show that without the limitations of the unitary operations, the new cooling technique would have exceeded the limit of Heat-bath algorithmic cooling. This work opens up new avenues to understanding the limits of dynamic cooling.


title: Efficient statistical validation with edge cases to evaluate Highly Automated Vehicles

abstract: The widescale deployment of Autonomous Vehicles (AV) seems to be imminent despite many safety challenges that are yet to be resolved. It is well known that there are no universally agreed Verification and Validation (VV) methodologies to guarantee absolute safety, which is crucial for the acceptance of this technology. Existing standards focus on deterministic processes where the validation requires only a set of test cases that cover the requirements. Modern autonomous vehicles will undoubtedly include machine learning and probabilistic techniques that require a much more comprehensive testing regime due to the non-deterministic nature of the operating design domain. A rigourous statistical validation process is an essential component required to address this challenge. Most research in this area focuses on evaluating system performance in large scale real-world data gathering exercises (number of miles travelled), or randomised test scenarios in simulation. This paper presents a new approach to compute the statistical characteristics of a system's behaviour by biasing automatically generated test cases towards the worst case scenarios, identifying potential unsafe edge cases.We use reinforcement learning (RL) to learn the behaviours of simulated actors that cause unsafe behaviour measured by the well established RSS safety metric. We demonstrate that by using the method we can more efficiently validate a system using a smaller number of test cases by focusing the simulation towards the worst case scenario, generating edge cases that correspond to unsafe situations.


title: Ising-based Consensus Clustering on Specialized Hardware

abstract: The emergence of specialized optimization hardware such as CMOS annealers and adiabatic quantum computers carries the promise of solving hard combinatorial optimization problems more efficiently in hardware. Recent work has focused on formulating different combinatorial optimization problems as Ising models, the core mathematical abstraction used by a large number of these hardware platforms, and evaluating the performance of these models when solved on specialized hardware. An interesting area of application is data mining, where combinatorial optimization problems underlie many core tasks. In this work, we focus on consensus clustering (clustering aggregation), an important combinatorial problem that has received much attention over the last two decades. We present two Ising models for consensus clustering and evaluate them using the Fujitsu Digital Annealer, a quantum-inspired CMOS annealer. Our empirical evaluation shows that our approach outperforms existing techniques and is a promising direction for future research.