title: Learning by example: fast reliability-aware seismic imaging with normalizing flows

abstract: Uncertainty quantification provides quantitative measures on the reliability of candidate solutions of ill-posed inverse problems. Due to their sequential nature, Monte Carlo sampling methods require large numbers of sampling steps for accurate Bayesian inference and are often computationally infeasible for large-scale inverse problems, such as seismic imaging. Our main contribution is a data-driven variational inference approach where we train a normalizing flow (NF), a type of invertible neural net, capable of cheaply sampling the posterior distribution given previously unseen seismic data from neighboring surveys. To arrive at this result, we train the NF on pairs of low- and high-fidelity migrated images. In our numerical example, we obtain high-fidelity images from the Parihaka dataset and low-fidelity images are derived from these images through the process of demigration, followed by adding noise and migration. During inference, given shot records from a new neighboring seismic survey, we first compute the reverse-time migration image. Next, by feeding this low-fidelity migrated image to the NF we gain access to samples from the posterior distribution virtually for free. We use these samples to compute a high-fidelity image including a first assessment of the image's reliability. To our knowledge, this is the first attempt to train a conditional network on what we know from neighboring images to improve the current image and assess its reliability.


title: Learning-Based UAV Trajectory Optimization with Collision Avoidance and Connectivity Constraints

abstract: Unmanned aerial vehicles (UAVs) are expected to be an integral part of wireless networks, and determining collision-free trajectories for multiple UAVs while satisfying requirements of connectivity with ground base stations (GBSs) is a challenging task. In this paper, we first reformulate the multi-UAV trajectory optimization problem with collision avoidance and wireless connectivity constraints as a sequential decision making problem in the discrete time domain. We, then, propose a decentralized deep reinforcement learning approach to solve the problem. More specifically, a value network is developed to encode the expected time to destination given the agent's joint state (including the agent's information, the nearby agents' observable information, and the locations of the nearby GBSs). A signal-to-interference-plus-noise ratio (SINR)-prediction neural network is also designed, using accumulated SINR measurements obtained when interacting with the cellular network, to map the GBSs' locations into the SINR levels in order to predict the UAV's SINR. Numerical results show that with the value network and SINR-prediction network, real-time navigation for multi-UAVs can be efficiently performed in various environments with high success rate.


title: A fault slip model to study earthquakes due to pore pressure perturbations

abstract: The burgeoning need to sequester anthropogenic CO$_2$ for climate mitigation and the need for energy sustenance leading upto enhanced geothermal energy production has made it incredibly critical to study potential earthquakes due to fluid activity in the subsurface. These earthquakes result from reactivation of faults in the subsurface due to pore pressure perturbations. In this work, we provide a framework to model fault slip due to pore pressure change leading upto quantifying the earthquake magnitude.


title: De Sitter scalar-spinor interaction in Minkowski limit

abstract: The scalar-spinor interaction Lagrangian is presented by the Yukawa potential. In dS ambient space formalism, the interaction Lagrangian of scalar-spinor fields was obtained from a new transformation which is very similar to the guage theory. The interaction of massless minimally coupled scalar and spinor fields was investigated. The Minkowski limit of the massless minimally coupled scalar field and massive spinor field interaction in the ambient space formalism of de Sitter space time is calculated. The interaction Lagrangian and massless minimally coupled scalar field in the null curvature limit become zero and the local transformation in the null curvature limit become a constant phase transformation and the interaction in this limit become zero. The covariant derivative reduces to ordinary derivative too. Then we conclude that this interaction is due to the curvature of space time and then the massless minimally coupled scalar field may be a part of a gravitational field.


title: Profitability Analysis in Stock Investment Using an LSTM-Based Deep Learning Model

abstract: Designing robust systems for precise prediction of future prices of stocks has always been considered a very challenging research problem. Even more challenging is to build a system for constructing an optimum portfolio of stocks based on the forecasted future stock prices. We present a deep learning-based regression model built on a long-and-short-term memory network (LSTM) network that automatically scraps the web and extracts historical stock prices based on a stock's ticker name for a specified pair of start and end dates, and forecasts the future stock prices. We deploy the model on 75 significant stocks chosen from 15 critical sectors of the Indian stock market. For each of the stocks, the model is evaluated for its forecast accuracy. Moreover, the predicted values of the stock prices are used as the basis for investment decisions, and the returns on the investments are computed. Extensive results are presented on the performance of the model. The analysis of the results demonstrates the efficacy and effectiveness of the system and enables us to compare the profitability of the sectors from the point of view of the investors in the stock market.


title: A Two-Level Fourth-Order Approach For Time-Fractional Convection-Diffusion-Reaction Equation With Variable Coefficients

abstract: This paper develops a two-level fourth-order scheme for solving time-fractional convection-diffusion-reaction equation with variable coefficients subjected to suitable initial and boundary conditions. The basis properties of the new approach are investigated and both stability and error estimates of the proposed numerical scheme are deeply analyzed in the $L^{\infty}(0,T;L^{2})$-norm. The theory indicates that the method is unconditionally stable with convergence of order $O(k^{2-\frac{\lambda}{2}}+h^{4})$, where $k$ and $h$ are time step and mesh size, respectively, and $\lambda\in(0,1)$. This result suggests that the two-level fourth-order technique is more efficient than a large class of numerical techniques widely studied in the literature for the considered problem. Some numerical evidences are provided to verify the unconditional stability and convergence rate of the proposed algorithm.


title: About connection of one matrix of composite numbers with Legendre's conjecture

abstract: Scientific paper is devoted to establish connection of T-matrix - matrix of composite numbers 6h+1 v 6h-1 in special view - with Legendre's conjecture.


title: On Determinism of Game Engines used for Simulation-based Autonomous Vehicle Verification

abstract: Game engines are increasingly used as simulation platforms by the autonomous vehicle (AV) community to develop vehicle control systems and test environments. A key requirement for simulation-based development and verification is determinism, since a deterministic process will always produce the same output given the same initial conditions and event history. Thus, in a deterministic simulation environment, tests are rendered repeatable and yield simulation results that are trustworthy and straightforward to debug. However, game engines are seldom deterministic. This paper reviews and identifies the potential causes of non-deterministic behaviours in game engines. A case study using CARLA, an open-source autonomous driving simulation environment powered by Unreal Engine, is presented to highlight its inherent shortcomings in providing sufficient precision in experimental results. Different configurations and utilisations of the software and hardware are explored to determine an operational domain where the simulation precision is sufficiently low i.e.\ variance between repeated executions becomes negligible for development and testing work. Finally, a method of a general nature is proposed, that can be used to find the domains of permissible variance in game engine simulations for any given system configuration.


title: Another Proof of $e^{x/y}$ being irrational

abstract: Continued fractions are used to give an alternate proof of $e^{x/y}$ is irrational.


title: CAN Coach: Vehicular Control through Human Cyber-Physical Systems

abstract: This work addresses whether a human-in-the-loop cyber-physical system (HCPS) can be effective in improving the longitudinal control of an individual vehicle in a traffic flow. We introduce the CAN Coach, which is a system that gives feedback to the human-in-the-loop using radar data (relative speed and position information to objects ahead) that is available on the controller area network (CAN). Using a cohort of six human subjects driving an instrumented vehicle, we compare the ability of the human-in-the-loop driver to achieve a constant time-gap control policy using only human-based visual perception to the car ahead, and by augmenting human perception with audible feedback from CAN sensor data. The addition of CAN-based feedback reduces the mean time-gap error by an average of 73%, and also improves the consistency of the human by reducing the standard deviation of the time-gap error by 53%. We remove human perception from the loop using a ghost mode in which the human-in-the-loop is coached to track a virtual vehicle on the road, rather than a physical one. The loss of visual perception of the vehicle ahead degrades the performance for most drivers, but by varying amounts. We show that human subjects can match the velocity of the lead vehicle ahead with and without CAN-based feedback, but velocity matching does not offer regulation of vehicle spacing. The viability of dynamic time-gap control is also demonstrated. We conclude that (1) it is possible to coach drivers to improve performance on driving tasks using CAN data, and (2) it is a true HCPS, since removing human perception from the control loop reduces performance at the given control objective.


title: Wireless Environment as a Service Enabled by Reconfigurable Intelligent Surfaces: The RISE-6G Perspective

abstract: The design of 6th Generation (6G) wireless networks points towards flexible connect-and-compute technologies capable to support innovative services and use cases. Targeting the 2030 horizon, 6G networks are poised to pave the way for sustainable human-centered smart societies and vertical industries, such that wireless networks will be transformed into a distributed smart connectivity infrastructure, where new terminal types are embedded in the daily environment. In this context, the RISE-6G project aims at investigating innovative solutions that capitalize on the latest advances in the emerging technology of Reconfigurable Intelligent Surfaces (RISs), which offers dynamic and goal-oriented radio wave propagation control, enabling the concept of the wireless environment as a service. The project will focus on: i) the realistic modeling of RIS-assisted signal propagation, ii) the investigation of the fundamental limits of RIS-empowered wireless communications and sensing, and iii) the design of efficient algorithms for orchestrating networking RISs, in order to implement intelligent, sustainable, and dynamically programmable wireless environments enabling diverse services that go well beyond the 5G capabilities. RISE-6G will offer two unprecedented proof-of-concepts for realizing controlled wireless environments in near-future use cases.


title: A comment on instantons and their fermion zero modes in adjoint QCD_2

abstract: The adjoint 2-dimensional $QCD$ with the gauge group $SU(N)/Z_N$ admits topologically nontrivial gauge field configurations associated with nontrivial $\pi_1[SU(N)/Z_N] = Z_N$. The topological sectors are labelled by an integer $k=0,\ldots, N-1$. However, in contrast to $QED_2$ and $QCD_4$, this topology is not associated with an integral invariant like the magnetic flux or Pontryagin index. These instantons may admit fermion zero modes, but there is always an equal number of left-handed and right-handed modes, so that the Atiyah-Singer theorem, which determines in other cases the number of the modes, does not apply. The mod. 2 argument suggests that, for a generic gauge field configuration, there is either a single doublet of such zero modes or no modes whatsoever. However, the known solution of the Dirac problem for a wide class of gauge field configurations indicates the presence of $k(N-k)$ zero mode doublets in the topological sector $k$. In this note, we demonstrate in an explicit way that these modes are not robust under a generic enough deformation of the gauge background and confirm thereby the mod. 2 conjecture. The implications for the physics of this theory (screening vs. confinement issue) are briefly discussed.


title: A Sufficient Condition to Guarantee Non-Simultaneous Charging and Discharging of Household Battery Energy Storage

abstract: In this letter, we model the day-ahead price-based demand response of a residential household with battery energy storage and other controllable loads, as a convex optimization problem. Further using duality theory and Karush-Kuhn-Tucker optimality conditions, we derive a sufficient criterion which guarantees non-simultaneous charging and discharging of the battery energy storage, without explicitly modelling it as a constraint


title: Multilingual Transfer Learning for Code-Switched Language and Speech Neural Modeling

abstract: In this thesis, we address the data scarcity and limitations of linguistic theory by proposing language-agnostic multi-task training methods. First, we introduce a meta-learning-based approach, meta-transfer learning, in which information is judiciously extracted from high-resource monolingual speech data to the code-switching domain. The meta-transfer learning quickly adapts the model to the code-switching task from a number of monolingual tasks by learning to learn in a multi-task learning fashion. Second, we propose a novel multilingual meta-embeddings approach to effectively represent code-switching data by acquiring useful knowledge learned in other languages, learning the commonalities of closely related languages and leveraging lexical composition. The method is far more efficient compared to contextualized pre-trained multilingual models. Third, we introduce multi-task learning to integrate syntactic information as a transfer learning strategy to a language model and learn where to code-switch. To further alleviate the aforementioned issues, we propose a data augmentation method using Pointer-Gen, a neural network using a copy mechanism to teach the model the code-switch points from monolingual parallel sentences. We disentangle the need for linguistic theory, and the model captures code-switching points by attending to input words and aligning the parallel words, without requiring any word alignments or constituency parsers. More importantly, the model can be effectively used for languages that are syntactically different, and it outperforms the linguistic theory-based models.


title: Density matrices in integrable face models

abstract: Using the properties of the local Boltzmann weights of integrable interaction-round-a-face (IRF or face) models we express local operators in terms of generalized transfer matrices. This allows for the derivation of discrete functional equations for the reduced density matrices in inhomogeneous generalizations of these models. We apply these equations to study the density matrices for IRF models of various solid-on-solid type and quantum chains of non-Abelian $\bm{su(2)_3}$ or Fibonacci anyons. Similar as in the six vertex model we find that reduced density matrices for a sequence of consecutive sites can be 'factorized', i.e.\ expressed in terms of nearest-neighbour correlators with coefficients which are independent of the model parameters. Explicit expressions are provided for correlation functions on up to three neighbouring sites.

categories: cond-mat.stat-mech hep-th math-ph math.MP


title: Pairs of equiperimeter and equiareal triangles whose sides are perfect squares

abstract: In this paper we consider the problem of finding pairs of triangles whose sides are perfect squares of integers, and which have a common perimeter and common area. We find two such pairs of triangles, and prove that there exist infinitely many pairs of triangles with the specified properties.


title: A Tale of Two Lexica Testing Computational Hypotheses with Deep Convolutional Neural Networks

abstract: Gow's (2012) dual lexicon model suggests that the primary purpose of words is to mediate the mappings between acoustic-phonetic input and other forms of linguistic representation. Motivated by evidence from functional imaging, aphasia, and behavioral results, the model argues for the existence of two parallel wordform stores: the dorsal and ventral processing streams. In this paper, we tested the hypothesis that the complex, but systematic mapping between sound and articulation in the dorsal stream poses different computational pressures on feature sets than the more arbitrary mapping between sound and meaning. To test this hypothesis, we created two deep convolutional neural networks (CNNs). While the dorsal network was trained to identify individual spoken words, the ventral network was trained to map them onto semantic classes. We then extracted patterns of network activation from the penultimate level of each network and tested how well features generated by the network supported generalization to linguistic categorization associated with the dorsal versus ventral processing streams. Our preliminary results showed both models successfully learned their tasks. Secondary generalization testing showed the ventral CNN outperformed the dorsal CNN on a semantic task: concreteness classification, while the dorsal CNN outperformed the ventral CNN on articulation tasks: classification by onset phoneme class and syllable length. These results are consistent with the hypothesis that the divergent processing demands of the ventral and dorsal processing streams impose computational pressures for the development of multiple lexica.


title: Podracer architectures for scalable Reinforcement Learning

abstract: Supporting state-of-the-art AI research requires balancing rapid prototyping, ease of use, and quick iteration, with the ability to deploy experiments at a scale traditionally associated with production systems.Deep learning frameworks such as TensorFlow, PyTorch and JAX allow users to transparently make use of accelerators, such as TPUs and GPUs, to offload the more computationally intensive parts of training and inference in modern deep learning systems. Popular training pipelines that use these frameworks for deep learning typically focus on (un-)supervised learning. How to best train reinforcement learning (RL) agents at scale is still an active research area. In this report we argue that TPUs are particularly well suited for training RL agents in a scalable, efficient and reproducible way. Specifically we describe two architectures designed to make the best use of the resources available on a TPU Pod (a special configuration in a Google data center that features multiple TPU devices connected to each other by extremely low latency communication channels).


title: $\mathbb{A}^1$-connected components of classifying spaces and purity for torsors

abstract: In this paper, we study the Nisnevich sheafification $\mathcal{H}^1_{\acute{e}t}(G)$ of the presheaf associating to a smooth scheme the set of isomorphism classes of $G$-torsors, for a reductive group $G$. We show that if $G$-torsors on affine lines are extended, then $\mathcal{H}^1_{\acute{e}t}(G)$ is homotopy invariant and show that the sheaf is unramified if and only if Nisnevich-local purity holds for $G$-torsors. We also identify the sheaf $\mathcal{H}^1_{\acute{e}t}(G)$ with the sheaf of $\mathbb{A}^1$-connected components of the classifying space ${\rm B}_{\acute{e}t}G$. This establishes the homotopy invariance of the sheaves of components as conjectured by Morel. It moreover provides a computation of the sheaf of $\mathbb{A}^1$-connected components in terms of unramified $G$-torsors over function fields whenever Nisnevich-local purity holds for $G$-torsors.


title: Optimal Data Placement for Data-Sharing Scientific Workflows in Heterogeneous Edge-Cloud Computing Environments

abstract: The heterogeneous edge-cloud computing paradigm can provide a more optimal direction to deploy scientific workflows than traditional distributed computing or cloud computing environments. Due to the different sizes of scientific datasets and some of these datasets must keep private, it is still a difficult problem to finding an data placement strategy that can minimize data transmission as well as placement cost. To address this issue, this paper combines advantages of both edge and cloud computing to construct a data placement model, which can balance data transfer time and data placement cost using intelligent computation. The most difficult research challenge the model solved is to consider many constrain in this hybrid computing environments, which including shared datasets within individual and among multiple workflows across various geographical regions. According to the constructed model, the study propose a new data placement strategy named DE-DPSO-DPS, which using a discrete particle swarm optimization algorithm with differential evolution (DE-DPSO-DPA) to distribute these scientific datasets. The strategy also not only consider the characteristics such as the number and storage capacity of edge micro-datacenters, the bandwidth between different datacenters and the proportion of private datasets, but also analysis the performance of algorithm during the workflows execution. Comprehensive experiments are designed in simulated heterogeneous edge-cloud computing environments demonstrate that the data placement strategy can effectively reduce the data transmission time and placement cost as compared to traditional strategies for data-sharing scientific workflows.


title: Mixed-horizon optimal feedback control as a model of human movement

abstract: Computational optimal feedback control (OFC) models in the sensorimotor control literature span a vast range of different implementations. Among the popular algorithms, finite-horizon, receding-horizon or infinite-horizon linear-quadratic regulators (LQR) have been broadly used to model human reaching movements. While these different implementations have their unique merits, all three have limitations in simulating the temporal evolution of visuomotor feedback responses. Here we propose a novel approach - a mixed-horizon OFC - by combining the strengths of the traditional finite-horizon and the infinite-horizon controllers to address their individual limitations. Specifically, we use the infinite-horizon OFC to generate durations of the movements, which are then fed into the finite-horizon controller to generate control gains. We then demonstrate the stability of our model by performing extensive sensitivity analysis of both re-optimisation and different cost functions. Finally, we use our model to provide a fresh look to previously published studies by reinforcing the previous results, providing alternative explanations to previous studies, or generating new predictive results for prior experiments.


title: Stein variational gradient descent with local approximations

abstract: Bayesian computation plays an important role in modern machine learning and statistics to reason about uncertainty. A key computational challenge in Bayesian inference is to develop efficient techniques to approximate, or draw samples from posterior distributions. Stein variational gradient decent (SVGD) has been shown to be a powerful approximate inference algorithm for this issue. However, the vanilla SVGD requires calculating the gradient of the target density and cannot be applied when the gradient is unavailable or too expensive to evaluate. In this paper we explore one way to address this challenge by the construction of a local surrogate for the target distribution in which the gradient can be obtained in a much more computationally feasible manner. More specifically, we approximate the forward model using a deep neural network (DNN) which is trained on a carefully chosen training set, which also determines the quality of the surrogate. To this end, we propose a general adaptation procedure to refine the local approximation online without destroying the convergence of the resulting SVGD. This significantly reduces the computational cost of SVGD and leads to a suite of algorithms that are straightforward to implement. The new algorithm is illustrated on a set of challenging Bayesian inverse problems, and numerical experiments demonstrate a clear improvement in performance and applicability of standard SVGD.


title: In-Network Volumetric DDoS Victim Identification Using Programmable Commodity Switches

abstract: Volumetric distributed Denial-of-Service (DDoS) attacks have become one of the most significant threats to modern telecommunication networks. However, most existing defense systems require that detection software operates from a centralized monitoring collector, leading to increased traffic load and delayed response. The recent advent of Data Plane Programmability (DPP) enables an alternative solution: threshold-based volumetric DDoS detection can be performed directly in programmable switches to skim only potentially hazardous traffic, to be analyzed in depth at the controller. In this paper, we first introduce the BACON data structure based on sketches, to estimate per-destination flow cardinality, and theoretically analyze it. Then we employ it in a simple in-network DDoS victim identification strategy, INDDoS, to detect the destination IPs for which the number of incoming connections exceeds a pre-defined threshold. We describe its hardware implementation on a Tofino-based programmable switch using the domain-specific P4 language, proving that some limitations imposed by real hardware to safeguard processing speed can be overcome to implement relatively complex packet manipulations. Finally, we present some experimental performance measurements, showing that our programmable switch is able to keep processing packets at line-rate while performing volumetric DDoS detection, and also achieves a high F1 score on DDoS victim identification.


title: Blue supergiants as tests for stellar physics

abstract: Massive star evolution is still poorly understood, and observational tests are required to discriminate between different implementations of physical phenomenon in stellar evolution codes. By confronting stellar evolution models with observed properties of blue supergiants, such as pulsations, chemical composition and position in the Hertzsprung-Russell diagram, we aim at determining which of the criterion used for convection (Schwarzschild or Ledoux) is best able to explain the observations. We compute state-of-the-art stellar evolution models with either the Schwarzschild or the Ledoux criterion for convection. Models are for $14$ to $35\,M_\odot$ at solar or Large Magellanic Cloud metallicity. For each model, we compute the pulsation properties to know when radial modes are excited. We then compare our results with the position of blue supergiants in the Hertzsprung-Russell diagram, with their surface chemical composition, and with their variability. Our results at Large Magellanic Cloud metallicity shows only a slight preference for the Ledoux criterion over the Schwarzschild one in reproducing at the same time the observed properties of blue supergiants, even if the Schwarzschild criterion cannot be excluded at this metallicity. We check that changing the overshoot parameter at solar metallicity does not improve the situation. We also check that our models are able to reproduce the position of Galactic blue supergiant in the flux-weighted-gravity -- luminosity relation. We confirm that overall, models computed with the Ledoux criterion are slightly better in matching observations. Our results also support the idea that most Galactic $\alpha$ Cyg variables are blue supergiants of the group 2, i.e. stars that have been through a previous red supergiant phase where they have lost large amount of mass.


title: Very Lightweight Photo Retouching Network with Conditional Sequential Modulation

abstract: Photo retouching aims at improving the aesthetic visual quality of images that suffer from photographic defects, especially for poor contrast, over/under exposure, and inharmonious saturation. In practice, photo retouching can be accomplished by a series of image processing operations. As most commonly-used retouching operations are pixel-independent, i.e., the manipulation on one pixel is uncorrelated with its neighboring pixels, we can take advantage of this property and design a specialized algorithm for efficient global photo retouching. We analyze these global operations and find that they can be mathematically formulated by a Multi-Layer Perceptron (MLP). Based on this observation, we propose an extremely lightweight framework -- Conditional Sequential Retouching Network (CSRNet). Benefiting from the utilization of $1\times1$ convolution, CSRNet only contains less than 37K trainable parameters, which are orders of magnitude smaller than existing learning-based methods. Experiments show that our method achieves state-of-the-art performance on the benchmark MIT-Adobe FiveK dataset quantitively and qualitatively. In addition to achieve global photo retouching, the proposed framework can be easily extended to learn local enhancement effects. The extended model, namely CSRNet-L, also achieves competitive results in various local enhancement tasks. Codes are available at https://github.com/lyh-18/CSRNet.


title: Fair Allocation of Conflicting Items

abstract: We study fair allocation of indivisible items, where the items are furnished with a set of conflicts, and agents are not permitted to receive conflicting items. This kind of constraint captures, for example, participating in events that overlap in time, or taking on roles in the presence of conflicting interests. We demonstrate, both theoretically and experimentally, that fairness characterizations such as EF1, MMS and MNW still are applicable and useful under item conflicts. Among other existence, non-existence and computability results, we show that a $1/\Delta$-approximate MMS allocation for $n$ agents may be found in polynomial time when $n>\Delta>2$, for any conflict graph with maximum degree $\Delta$, and that, if $n > \Delta$, a 1/3-approximate MMS allocation always exists.


title: Decoupling lattice and magnetic instabilities in frustrated CuMnO$_2$

abstract: The $A$MnO$_{2}$ delafossites ($A$=Na, Cu), are model frustrated antiferromagnets, with triangular layers of Mn$^{3+}$~spins. At low temperatures ($T_{N}$=65 K), a $C2/m \rightarrow P\overline{1}$ transition is found in CuMnO$_2$, which breaks frustration and establishes magnetic order. In contrast to this clean transition, $A$=Na only shows short-range distortions at $T_N$. Here we report a systematic crystallographic, spectroscopic, and theoretical investigation of CuMnO$_2$. We show that, even in stoichiometric samples, non-zero anisotropic Cu displacements co-exist with magnetic order. Using X-ray/neutron diffraction and Raman scattering, we show that high pressures acts to decouple these degrees of freedom. This manifests as an isostuctural phase transition at $\sim$10 GPa, with a reversible collapse of the $c$-axis. This is shown to be the high pressure analog of the $c$-axis negative thermal expansion seen at ambient pressure. DFT simulations confirm that dynamical instabilities of the Cu$^{+}$ cations and edge-shared MnO$_{6}$ layers are intertwined at ambient pressure. However, high pressure selectively activates the former, before an eventual predicted re-emergence of magnetism at the highest pressures. Our results show that the lattice dynamics and local structure of CuMnO$_2$ are quantitatively different to non-magnetic Cu delafossites, and raise questions about the role of intrinsic inhomogeniety in frustrated antiferromagnets.


title: Numerics and computation in gyrokinetic simulations of electromagnetic turbulence with global particle-in-cell codes

abstract: Electromagnetic turbulence is addressed in tokamak and stellarator plasmas with the global gyrokinetic particle-in-cell codes ORB5 [E. Lanti et al, Comp. Phys. Comm, vol. 251, 107072 (2020)] and EUTERPE [V. Kornilov et al, Phys. Plasmas, vol. 11, 3196 (2004)]. The large-aspect-ratio tokamak, down-scaled ITER, and Wendelstein 7-X geometries are considered. The main goal is to increase the plasma beta, the machine size, the ion-to-electron mass ratio, as well as to include realistic-geometry features in such simulations. The associated numerical requirements and the computational cost for the cases on computer systems with massive GPU deployments are investigated. These are necessary steps to enable electromagnetic turbulence simulations in future reactor plasmas.


title: Energy Efficiency Optimization of Reconfigurable Intelligent Surfaces with Electromagnetic Field Exposure Constraints

abstract: This work considers the problem of energy efficiency maximization in a RIS-based communication link, subject to not only the conventional maximum power constraints, but also additional constraints on the maximum exposure to electromagnetic radiations of the end-users. The RIS phase shifts, the transmit beamforming, the linear receive filter, and the transmit power are jointly optimized, and two provably convergent and low-complexity algorithms are developed. One algorithm can be applied to the general system setups, but does not guarantee global optimality. The second algorithm is provably optimal in a notable special case. The numerical results show that RIS-based communications can ensure high energy efficiency while fulfilling users' exposure constraints to radio frequency emissions.


title: Noisy quantum metrology with the assistance of indefinite causal order

abstract: A generic qubit unitary operator affected by depolarizing noise is duplicated and inserted in a quantum switch process realizing a superposition of causal orders. The characterization of the resulting switched quantum channel is worked out for its action on the joint state of the probe-control qubit pair. The switched channel is then specifically investigated for the important metrological task of phase estimation on the noisy unitary operator, with the performance assessed by the Fisher information, classical or quantum. A comparison is made with conventional techniques of estimation where the noisy unitary is directly probed in a one-stage or two-stage cascade with definite order, or several uses of them with two or more qubits. In the switched channel with indefinite order, specific properties are reported, meaningful for estimation and not present with conventional techniques. It is shown that the control qubit, although it never directly interacts with the unitary, can nevertheless be measured alone for effective estimation, while discarding the probe qubit that interacts with the unitary. Also, measurement of the control qubit maintains the possibility of efficient estimation in difficult conditions where conventional estimation becomes less efficient, as with ill-configured input probes, or in blind situations when the axis of the unitary is unknown. Effective estimation by measuring the control qubit remains possible even when the input probe tends to align with the axis of the unitary, or with a fully depolarized input probe, while in these conditions conventional estimation gets inoperative. Measurement of the probe qubit of the switched channel is also shown to add useful capabilities for phase estimation. The results contribute to the analysis of switched quantum channels with indefinite order for information processing, and uncover new possibilities for qubit metrology.


title: An acceleration strategy for randomize-then-optimize sampling via deep neural networks

abstract: Randomize-then-optimize (RTO) is widely used for sampling from posterior distributions in Bayesian inverse problems. However, RTO may be computationally intensive for complexity problems due to repetitive evaluations of the expensive forward model and its gradient. In this work, we present a novel strategy to substantially reduce the computation burden of RTO by using a goal-oriented deep neural networks (DNN) surrogate approach. In particular, the training points for the DNN-surrogate are drawn from a local approximated posterior distribution, and it is shown that the resulting algorithm can provide a flexible and efficient sampling algorithm, which converges to the direct RTO approach. We present a Bayesian inverse problem governed by a benchmark elliptic PDE to demonstrate the computational accuracy and efficiency of our new algorithm (i.e., DNN-RTO). It is shown that with our algorithm, one can significantly outperform the traditional RTO.


title: Naturality of ${\rm SL}_3$ quantum trace maps for surfaces

abstract: Fock-Goncharov's moduli spaces $\mathscr{X}_{{\rm PGL}_3,\frak{S}}$ of framed ${\rm PGL}_3$-local systems on punctured surfaces $\frak{S}$ provide prominent examples of cluster $\mathscr{X}$-varieties and higher Teichm\"uller spaces. In a previous paper of the author (arXiv:2011.14765), building on the works of others, the so-called ${\rm SL}_3$ quantum trace map is constructed for each triangulable punctured surface $\frak{S}$ and an ideal triangulation $\Delta$ of $\frak{S}$, as a homomorphism from the stated ${\rm SL}_3$-skein algebra of the surface to a quantum torus algebra that deforms the ring of Laurent polynomials in the cube-roots of the cluster coordinate variables for the cluster $\mathscr{X}$-chart for $\mathscr{X}_{{\rm PGL}_3,\frak{S}}$ associated to $\Delta$. We develop quantum mutation maps between special subalgebras of the cube-root quantum torus algebras for different triangulations and show that the ${\rm SL}_3$ quantum trace maps are natural, in the sense that they are compatible under these quantum mutation maps. As an application, the quantum ${\rm SL}_3$-${\rm PGL}_3$ duality map constructed in the previous paper is shown to be independent of the choice of an ideal triangulation.


title: Modeling the dominance of the gradient drift or Kelvin-Helmholtz instability in sheared ionospheric E x B flows

abstract: Studies have shown that in sheared $\mathbf{E}\times\mathbf{B}$ flows in an inhomogeneous ionospheric plasma, the gradient drift (GDI) or the Kelvin-Helmholtz (KHI) instability may grow. This work examines the conditions that cause one of these instabilities to dominate over the other using a novel model to study localized ionospheric instabilities. The effect of collisions with neutral particles plays an important role in the instability development. It is found that the KHI is dominant in low collisionality regimes, the GDI is dominant in high collisionality regimes, and there exists an intermediate region in which both instabilities exist in tandem. For low collisionality cases in which the velocity shear is sufficiently far from the density gradient, the GDI is found to grow as a secondary instability extending from the KHI vortices. The inclusion of a neutral wind driven electric field in the direction of the velocity shear does not impact the dominance of either instability. Using data from empirical ionospheric models, two altitude limits are found. For altitudes above the higher limit, the KHI is dominant. For altitudes below the lower limit, the GDI is dominant. In the intermediate region, both instabilities grow together. Increasing the velocity shear causes both limits to be lower in altitude. This implies that for ionospheric phenomena whose density and velocity gradients span large altitude ranges, such as subauroral polarization streams, the instabilities observed by space-based and ground-based observation instruments could be significantly different.


title: Defocused travelling-fringes in scanning triple-Laue x-ray interferometry

abstract: The measurement of the silicon lattice parameter by a separate-crystal triple-Laue x-ray interferometer is a key step for the kilogram realisation by counting atoms. Since the measurement accuracy is approaching nine significant digits, a reliable model of the interferometer operation is demanded to quantify or exclude systematic errors. This paper investigates both analytically and experimentally the effect of defocus (a difference between the splitter-to-mirror distance on the one hand and the analyser-to-mirror one on the other) on the phase of the interference fringes and the measurement of the lattice parameter.


title: Traffic models and traffic-jam transition in quantum ($N$+1)-level systems

abstract: We propose a model to implement and simulate different traffic-flow conditions in terms of quantum graphs hosting an ($N$+1)-level dot at each site, which allows us to keep track of the type and of the destination of each vehicle. By implementing proper Lindbladian local dissipators, we derive the master equations that describe the traffic flow in our system. To show the versatility and the reliability of our technique, we employ it to model different types of traffic flow (the symmetric three-way roundabout and the three-road intersection). Eventually, we successfully compare our predictions with results from classical models.


title: Fermat functional equations over Riemann surfaces

abstract: We investigate the existence of non-trivial holomorphic and meromorphic solutions of Fermat functional equations over an open Riemann surface $S$. When $S$ is hyperbolic, we prove that any $k$-term Fermat functional equation always exists non-trivial holomorphic and meromorphic solution. When $S$ is a general open Riemann surface, we prove that every non-trivial holomorphic or meromorphic solution satisfies a growth condition, provided that the power exponents of the equations are bigger than some certain positive integers.


title: Optimal design for kernel interpolation: applications to uncertainty quantification

abstract: The paper is concerned with classic kernel interpolation methods, in addition to approximation methods that are augmented by gradient measurements. To apply kernel interpolation using radial basis functions (RBFs) in a stable way, we propose a type of quasi-optimal interpolation points, searching from a large set of \textit{candidate} points, using a procedure similar to designing Fekete points or power function maximizing points that use pivot from a Cholesky decomposition. The proposed quasi-optimal points result in a smaller condition number, and thus mitigates the instability of the interpolation procedure when the number of points becomes large. Applications to parametric uncertainty quantification are presented, and it is shown that the proposed interpolation method can outperform sparse grid methods in many interesting cases. We also demonstrate the new procedure can be applied to constructing gradient-enhanced Gaussian process emulators.


title: Nonlocal cross-diffusion systems for multi-species populations and networks

abstract: Nonlocal cross-diffusion systems on the torus, arising in population dynamics and neuroscience, are analyzed. The global existence of weak solutions, the weak-strong uniqueness, and the localization limit are proved. The kernels are assumed to be positive definite and in detailed balance. The proofs are based on entropy estimates coming from Shannon-type and Rao-type entropies, while the weak-strong uniqueness result follows from the relative entropy method. The existence and uniqueness theorems hold for nondifferentiable kernels. The associated local cross-diffusion system, derived in the localization limit, is also discussed.


title: Analysis of optimal portfolio on finite and small time horizons for a stochastic volatility market model

abstract: In this paper, we consider the portfolio optimization problem in a financial market under a general utility function. Empirical results suggest that if a significant market fluctuation occurs, invested wealth tends to have a notable change from its current value. We consider an incomplete stochastic volatility market model, that is driven by both a Brownian motion and a jump process. At first, we obtain a closed-form formula for an approximation to the optimal portfolio in a small-time horizon. This is obtained by finding the associated Hamilton-Jacobi-Bellman integro-differential equation and then approximating the value function by constructing appropriate super-solution and sub-solution. It is shown that the true value function can be obtained by sandwiching the constructed super-solution and sub-solution. We also prove the accuracy of the approximation formulas. Finally, we provide a procedure for generating a close-to-optimal portfolio for a finite time horizon.


title: Online and Offline Reinforcement Learning by Planning with a Learned Model

abstract: Learning efficiently from small amounts of data has long been the focus of model-based reinforcement learning, both for the online case when interacting with the environment and the offline case when learning from a fixed dataset. However, to date no single unified algorithm could demonstrate state-of-the-art results in both settings. In this work, we describe the Reanalyse algorithm which uses model-based policy and value improvement operators to compute new improved training targets on existing data points, allowing efficient learning for data budgets varying by several orders of magnitude. We further show that Reanalyse can also be used to learn entirely from demonstrations without any environment interactions, as in the case of offline Reinforcement Learning (offline RL). Combining Reanalyse with the MuZero algorithm, we introduce MuZero Unplugged, a single unified algorithm for any data budget, including offline RL. In contrast to previous work, our algorithm does not require any special adaptations for the off-policy or offline RL settings. MuZero Unplugged sets new state-of-the-art results in the RL Unplugged offline RL benchmark as well as in the online RL benchmark of Atari in the standard 200 million frame setting.


title: LyST: a Scalar-Tensor Theory of Gravity on Lyra Manifold

abstract: We present a scalar-tensor theory of gravity on a torsion-free and metric compatible Lyra manifold. This is obtained by generalizing the concept of physical reference frame by considering a scale function defined over the manifold. The choice of a specific frame induces a local base, naturally non-holonomic, whose structure constants give rise to extra terms in the expression of the connection coefficients and in the expression for the covariant derivative. In the Lyra manifold, transformations between reference frames involving both coordinates and scale change the transformation law of tensor fields, when compared to those of the Riemann manifold. From a direct generalization of the Einstein-Hilbert minimal action coupled with a matter term, it was possible to build a Lyra invariant action, which gives rise to the associated Lyra Scalar-Tensor theory of gravity (LyST), with field equations for $g_{\mu\nu}$ and $\phi$. These equations have a well-defined Newtonian limit, from which it can be seen that both metric and scale play a role in the description gravitational interaction. We present a spherically symmetric solution for the LyST gravity field equations. It dependent on two parameters $m$ and $r_{L}$, whose physical meaning is carefully investigated. We highlight the properties of LyST spherically symmetric line element and compare it to Schwarzchild solution.


title: Count Network Autoregression

abstract: We consider network autoregressive models for count data with a non-random neighborhood structure. The main methodological contribution is the development of conditions that guarantee stability and valid statistical inference for such models. We consider both cases of fixed and increasing network dimension and we show that quasi-likelihood inference provides consistent and asymptotically normally distributed estimators. The work is complemented by simulation results and a data example.


title: Adversarial autoencoders and adversarial LSTM for improved forecasts of urban air pollution simulations

abstract: This paper presents an approach to improve the forecast of computational fluid dynamics (CFD) simulations of urban air pollution using deep learning, and most specifically adversarial training. This adversarial approach aims to reduce the divergence of the forecasts from the underlying physical model. Our two-step method integrates a Principal Components Analysis (PCA) based adversarial autoencoder (PC-AAE) with adversarial Long short-term memory (LSTM) networks. Once the reduced-order model (ROM) of the CFD solution is obtained via PCA, an adversarial autoencoder is used on the principal components time series. Subsequentially, a Long Short-Term Memory network (LSTM) is adversarially trained on the latent space produced by the PC-AAE to make forecasts. Once trained, the adversarially trained LSTM outperforms a LSTM trained in a classical way. The study area is in South London, including three-dimensional velocity vectors in a busy traffic junction.


title: Effect of stand-off distance and spatial resolution on the pressure impact of near-wall vapor bubble collapses

abstract: We consider the collapse behavior of cavitation bubbles near walls under high ambient pressure conditions. Generic configurations with different stand-off distances are investigated by numerical simulation using a fully compressible two-phase flow solver including phase change. The results show that the stand-off distance has significant effects on collapse dynamics, micro-jet formation, rebound, and maximum wall pressure. A relation between cavitation induced material damage and corresponding collapse mechanisms is obtained from pressure-impact data at the wall. We analyze the resolution dependence of collapse and rebound and the observed maximum pressure distributions. The comparison of the results on six different grid resolutions shows that main collapse features are already captured on the coarsest resolution, while the peak pressures are strongly resolution dependent.


title: A review and evaluation of secondary school accountability in England: Statistical strengths, weaknesses, and challenges for 'Progress 8'

abstract: School performance measures are published annually in England to hold schools to account and to support parental school choice. This article reviews and evaluates the Progress 8 secondary school accountability system for state-funded schools. We assess the statistical strengths and weaknesses of Progress 8 relating to: choice of pupil outcome attainment measures; potential adjustments for pupil input attainment and background characteristics; decisions around which schools and pupils are excluded from the measure; presentation of Progress 8 to users, choice of statistical model, and calculation of statistical uncertainty; and issues related to the volatility of school performance over time, including scope for reporting multi-year averages. We then discuss challenges for Progress 8 raised by the COVID-19 pandemic. Six simple recommendations follow to improve Progress 8 and school accountability in England.


title: Chaos and turbulence in clouds

abstract: Spatial buoyancy-helical distributed chaos (turbulence) in the cumulus, stratiform, stratocumulus, cirrus and cirrus mammatus clouds have been studied using results of direct numerical simulations and measurements in cloudy atmosphere. It is shown that in the considered cases the moments of helicity distribution dominate the kinetic energy spectra both in the middle of the clouds and in the cloud-top regions.


title: A single-qubit position verification protocol that is secure against multi-qubit attacks

abstract: The position of a device or agent is an important security credential in today's society, both online and in the real world. Unless in direct proximity, however, the secure verification of a position is impossible without further assumptions. This is true classically, but also in any future quantum-equipped communications infrastructure. We show in this work that minimal quantum resources, in the form of a single qubit, combined with classical communication are sufficient to thwart quantum adversaries that pretend to be at a specific position and have the ability to coordinate their action with entanglement. More precisely, we show that the adversaries using an increasing amount of entanglement can be combatted solely by increasing the number of classical bits used in the protocol. The presented protocols are noise-robust and within reach of current quantum technology.


title: Stabilisation of the complex double integrator by means of a saturated linear feedback

abstract: Consider the saturated complex double integrator, i.e., the linear control system $\dot x=Ax+B\sigma(u)$, where $x\in\R^4$, $u\in\R$, $B\in\R^4$, the $4\times 4$ matrix $A$ is not diagolizable and admits a non zero purely imaginary eigenvalue of multiplicity two, the pair $(A,B)$ is controllable and $\sigma:\R\to\R$ is a saturation function. We prove that there exists a linear feedback $u=K^Tx$ such that the resulting closed loop system given by $\dot x=Ax+B\sigma(K^Tx)$ is globally asymptotically stable with respect to the origin.


title: Learning and Planning in Complex Action Spaces

abstract: Many important real-world problems have action spaces that are high-dimensional, continuous or both, making full enumeration of all possible actions infeasible. Instead, only small subsets of actions can be sampled for the purpose of policy evaluation and improvement. In this paper, we propose a general framework to reason in a principled way about policy evaluation and improvement over such sampled action subsets. This sample-based policy iteration framework can in principle be applied to any reinforcement learning algorithm based upon policy iteration. Concretely, we propose Sampled MuZero, an extension of the MuZero algorithm that is able to learn in domains with arbitrarily complex action spaces by planning over sampled actions. We demonstrate this approach on the classical board game of Go and on two continuous control benchmark domains: DeepMind Control Suite and Real-World RL Suite.


title: Lifetime Optimization of Dense Wireless Sensor Networks Using Continuous Ring-sector Model

abstract: Wireless sensor networks (WSNs) are becoming increasingly utilized in applications that require remote collection of data on environmental conditions. In particular dense WSNs are emerging as an important sensing platforms for the Internet of Things (IoT). WSNs are able to generate huge volumes of raw data, which require network structuring and efficient collaboration between nodes to ensure efficient transmission. In order to reduce the amount of data carried in the network, data aggregation is used in WSNs to define a policy of data fusion and compression. In this paper, we investigate a model for data aggregation in a dense {WSN} with a single sink. The model divides a circular coverage region centered at the sink into patches which are intersections of sectors of concentric rings, and data in each patch is aggregated at a single node before transmission. Nodes only communicate with other nodes in the same sector. Based on these assumptions, we formulate a linear programming problem to maximize system lifetime by minimizing the maximum proportionate energy consumption over all nodes. Under a wide variety of conditions, the optimal solution employs two transmissions mechanisms: direct transmission, in which nodes send information directly to the sink; and stepwise transmission, in which nodes transmit information to adjacent nodes. An exact formula is given for the proportionate energy consumption rate of the network. Asymptotic forms of this exact solution are also derived, and are verified to agree with the linear programming solution. We investigate three strategies for improving system lifetime: nonuniform energy and information density; iterated compression; and modifications of rings. We conclude that iterated compression has the biggest effect in increasing system lifetime.


title: Astronomy Education in Covid Times

abstract: We shall describe the various activities done by us in Covid Times including outreach and educational workshops in Physics and Astronomy. We shall discuss the caveats in virtual teaching of Astronomy and the lessons learnt in the process.


title: Port Parameter Extraction Based Self Consistent Coupled EM-Circuit FEM Solvers

abstract: Self consistent solution to electromagnetic (EM)-circuit systems is of significant interest for a number of applications. This has resulted in exhaustive research on means to couple them. In time domain, this typically involves a tight integration (or coupling) with field and non-linear circuit solvers. This is in stark contrast to coupled analysis of linear/weakly non-linear circuits and EM systems in frequency domain. Here, one typically extracts equivalent port parameters that are then fed into the circuit solver. Such an approach has several advantages; (a) the number of ports is typically smaller than the number of degrees of freedom, resulting in cost savings; (b) is circuit agnostic; (c) can be integrated with a variety of device models. Port extraction is tantamount to obtaining impulse response of the linear EM system. In time domain, the deconvolution required to effect this is unstable. Recently, a novel approach was developed for time domain integral equations to overcome this bottleneck. We extend this approach to time domain finite element method, and demonstrate its utility via a number of examples; significantly, we demonstrate that self consistent solutions obtained using either a fully coupled or port extraction is identical to the desired precision for non-linear circuit systems. This is shown within a nodal network. We also demonstrate integration of port extracted data directly with drift diffusion equation to model device physics.


title: Detecting False Data Injection Attacks in Smart Grids with Modeling Errors: A Deep Transfer Learning Based Approach

abstract: Most traditional false data injection attack (FDIA) detection approaches rely on a key assumption, i.e., the power system can be accurately modeled. However, the transmission line parameters are dynamic and cannot be accurately known during operation and thus the involved modeling errors should not be neglected. In this paper, an illustrative case has revealed that modeling errors in transmission lines significantly weaken the detection effectiveness of conventional FDIA approaches. To tackle this issue, we propose an FDIA detection mechanism from the perspective of transfer learning. Specifically, the simulated power system is treated as a source domain, which provides abundant simulated normal and attack data. The real world's running system whose transmission line parameters are unknown is taken as a target domain where sufficient real normal data are collected for tracking the latest system states online. The designed transfer strategy that aims at making full use of data in hand is divided into two optimization stages. In the first stage, a deep neural network (DNN) is built by simultaneously optimizing several well-designed objective terms with both simulated data and real data, and then it is fine-tuned via real data in the second stage. Several case studies on the IEEE 14-bus and 118-bus systems verify the effectiveness of the proposed mechanism.


title: SFE-Net: EEG-based Emotion Recognition with Symmetrical Spatial Feature Extraction

abstract: Emotion recognition based on EEG (electroencephalography) has been widely used in human-computer interaction, distance education and health care. However, the conventional methods ignore the adjacent and symmetrical characteristics of EEG signals, which also contain salient information related to emotion. In this paper, a spatial folding ensemble network (SFE-Net) is presented for EEG feature extraction and emotion recognition. Firstly, for the undetected area between EEG electrodes, an improved Bicubic-EEG interpolation algorithm is developed for EEG channels information completion, which allows us to extract a wider range of adjacent space features. Then, motivated by the spatial symmetric mechanism of human brain, we fold the input EEG channels data with five different symmetrical strategies, which enable the proposed network to extract the information of space features of EEG signals more effectively. Finally, a 3DCNN-based spatial, temporal extraction, and a multi-voting strategy of ensemble learning are integrated to model a new neural network. With this network, the spatial features of different symmetric folding signals can be extracted simultaneously, which greatly improves the robustness and accuracy of emotion recognition. The experimental results on DEAP and SEED datasets show that the proposed algorithm has comparable performance in terms of recognition accuracy.


title: Signal Processing and Machine Learning Techniques for Terahertz Sensing: An Overview

abstract: Following the recent progress in Terahertz (THz) signal generation and radiation methods, joint THz communications and sensing applications are shaping the future of wireless systems. Towards this end, THz spectroscopy is expected to be carried over user equipment devices to identify material and gaseous components of interest. THz-specific signal processing techniques should complement this re-surged interest in THz sensing for efficient utilization of the THz band. In this paper, we present an overview of these techniques, with an emphasis on signal pre-processing (standard normal variate normalization, min-max normalization, and Savitzky-Golay filtering), feature extraction (principal component analysis, partial least squares, t-distributed stochastic neighbor embedding, and nonnegative matrix factorization), and classification techniques (support vector machines, k-nearest neighbor, discriminant analysis, and naive Bayes). We also address the effectiveness of deep learning techniques by exploring their promising sensing capabilities at the THz band. Lastly, we investigate the performance and complexity trade-offs of the studied methods in the context of joint communications and sensing; we motivate the corresponding use-cases, and we present few future research directions in the field.


title: Exploration of Spanish Olive Oil Quality with a Miniaturized Low-Cost Fluorescence Sensor and Machine Learning Techniques

abstract: Extra virgin olive oil (EVOO) is the highest quality of olive oil and is characterized by highly beneficial nutritional properties. The large increase in both consumption and fraud, for example through adulteration, creates new challenges and an increasing demand for developing new quality assessment methodologies that are easier and cheaper to perform. As of today, the determination of olive oil quality is performed by producers through chemical analysis and organoleptic evaluation. The chemical analysis requires the advanced equipment and chemical knowledge of certified laboratories, and has therefore a limited accessibility. In this work a minimalist, portable and low-cost sensor is presented, which can perform olive oil quality assessment using fluorescence spectroscopy. The potential of the proposed technology is explored by analyzing several olive oils of different quality levels, EVOO, virgin olive oil (VOO), and lampante olive oil (LOO). The spectral data were analyzed using a large number of machine learning methods, including artificial neural networks. The analysis performed in this work demonstrates the possibility of performing classification of olive oil in the three mentioned classes with an accuracy of 100$\%$. These results confirm that this minimalist low-cost sensor has the potential of substituting expensive and complex chemical analysis.


title: 10,000 km Straight-line Transmission using a Real-time Software-defined GPU-Based Receiver

abstract: Real-time operation of a software-defined, GPU-based optical receiver is demonstrated over a 100-span straight-line optical link. Performance of minimum-phase Kramers-Kronig 4-, 8-, 16-, 32-, and 64-QAM signals are evaluated at various distances.


title: A Non-sequential Approach to Deep User Interest Model for CTR Prediction

abstract: Click-Through Rate (CTR) prediction plays an important role in many industrial applications, and recently a lot of attention is paid to the deep interest models which use attention mechanism to capture user interests from historical behaviors. However, most current models are based on sequential models which truncate the behavior sequences by a fixed length, thus have difficulties in handling very long behavior sequences. Another big problem is that sequences with the same length can be quite different in terms of time, carrying completely different meanings. In this paper, we propose a non-sequential approach to tackle the above problems. Specifically, we first represent the behavior data in a sparse key-vector format, where the vector contains rich behavior info such as time, count and category. Next, we enhance the Deep Interest Network to take such rich information into account by a novel attention network. The sparse representation makes it practical to handle large scale long behavior sequences. Finally, we introduce a multidimensional partition framework to mine behavior interactions. The framework can partition data into custom designed time buckets to capture the interactions among information aggregated in different time buckets. Similarly, it can also partition the data into different categories and capture the interactions among them. Experiments are conducted on two public datasets: one is an advertising dataset and the other is a production recommender dataset. Our models outperform other state-of-the-art models on both datasets.


title: SetConv: A New Approach for Learning from Imbalanced Data

abstract: For many real-world classification problems, e.g., sentiment classification, most existing machine learning methods are biased towards the majority class when the Imbalance Ratio (IR) is high. To address this problem, we propose a set convolution (SetConv) operation and an episodic training strategy to extract a single representative for each class, so that classifiers can later be trained on a balanced class distribution. We prove that our proposed algorithm is permutation-invariant despite the order of inputs, and experiments on multiple large-scale benchmark text datasets show the superiority of our proposed framework when compared to other SOTA methods.


title: Energy-Efficient 3D Deployment of Aerial Access Points in a UAV Communication System

abstract: In this letter, we propose an energy-efficient 3-dimensional placement of multiple aerial access points (AAPs), in the desired area, acting as flying base stations for uplink communication from a set of ground user equipment (UE). The globally optimal energy-efficient vertical position of AAPs is derived analytically by considering the inter-cell interference and AAP energy consumption. The horizontal position of AAPs which maximize the packing density of the AAP coverage area are determined using a novel regular polygon-based AAP placement algorithm. We also determine the maximum number of non-interfering AAPs that can be placed in the desired area. The effect of the AAP energy consumption on the optimal placement and the analytic findings are verified via numerical simulations.


title: Differential chaos shift keying-based wireless power transfer

abstract: In this work, we investigate differential chaos shift keying (DCSK), a communication-based waveform, in the context of wireless power transfer (WPT). Particularly, we present a DCSK-based WPT architecture, that employs an analog correlator at the receiver in order to boost the energy harvesting (EH) performance. By taking into account the nonlinearities of the EH process, we derive closed-form analytical expressions for the peak-to-average-power-ratio of the received signal as well as the harvested power. Nontrivial design insights are provided, where it is shown how the parameters of the transmitted waveform affects the EH performance. Furthermore, it is demonstrated that the employment of a correlator at the receiver achieves significant EH gains in DCSK-based WPT systems.


title: A Novel Solution of Using Mixed Reality in Bowel and Oral and Maxillofacial Surgical Telepresence: 3D Mean Value Cloning algorithm

abstract: Background and aim: Most of the Mixed Reality models used in the surgical telepresence are suffering from discrepancies in the boundary area and spatial-temporal inconsistency due to the illumination variation in the video frames. The aim behind this work is to propose a new solution that helps produce the composite video by merging the augmented video of the surgery site and the virtual hand of the remote expertise surgeon. The purpose of the proposed solution is to decrease the processing time and enhance the accuracy of merged video by decreasing the overlay and visualization error and removing occlusion and artefacts. Methodology: The proposed system enhanced the mean value cloning algorithm that helps to maintain the spatial-temporal consistency of the final composite video. The enhanced algorithm includes the 3D mean value coordinates and improvised mean value interpolant in the image cloning process, which helps to reduce the sawtooth, smudging and discolouration artefacts around the blending region. Results: As compared to the state of the art solution, the accuracy in terms of overlay error of the proposed solution is improved from 1.01mm to 0.80mm whereas the accuracy in terms of visualization error is improved from 98.8% to 99.4%. The processing time is reduced to 0.173 seconds from 0.211 seconds. Conclusion: Our solution helps make the object of interest consistent with the light intensity of the target image by adding the space distance that helps maintain the spatial consistency in the final merged video.


title: Probing Negative Sampling Strategies to Learn GraphRepresentations via Unsupervised Contrastive Learning

abstract: Graph representation learning has long been an important yet challenging task for various real-world applications. However, their downstream tasks are mainly performed in the settings of supervised or semi-supervised learning. Inspired by recent advances in unsupervised contrastive learning, this paper is thus motivated to investigate how the node-wise contrastive learning could be performed. Particularly, we respectively resolve the class collision issue and the imbalanced negative data distribution issue. Extensive experiments are performed on three real-world datasets and the proposed approach achieves the SOTA model performance.


title: Rigidity of proper holomorphic maps among generalized balls with Levi-degenerate boundaries

abstract: In this paper we studied a broader type of generalized balls which are domains on the complex projective with possibly Levi-degenerate boundaries. We proved rigidity theorems for proper holomorphic mappings among them by exploring the structure of the moduli spaces of projective linear subspaces, which generalized some earlier results for the ordinary generalized balls with Levi-nondegenerate boundaries.


title: Integrable modulation, curl forces and parametric Kapitza equation with trapping and escaping

abstract: In this present communication the integrable modulation problem has been applied to study parametric extension of the Kapitza rotating shaft problem, which is a protypical example of curl force as formulated by Berry and Shukla in [J. Phys. A 45 305201 (2012)] associated with simple saddle potential. The integrable modulation problems yield parametric time dependent integrable systems. The Hamiltonian and first integrals of the linear and nonlinear parametric Kapitza equation (PKE) associated with simple and monkey saddle potentials have been given. The construction has been illustrated by choosing $\omega(t) = a + b \cos t$ and that maps to Mathieu type equations, which yield Mathieu extension of PKE. We study the dynamics of these equations. The most interesting finding is the mixed mode of particle trapping and escaping via the heteroclinic orbits depicted with the parametric Mathieu-Kapitza equation which are absent in the case of non parametric cases.


title: $(g-2)_{\mu}$ anomaly within the left-right symmetric model

abstract: The muon anomalous magnetic moment is discussed within the model based on the $SU(2)_L\times SU(2)_R\times U(1)_{B-L}$-gauge group. The contributions of the neutrinos having the dipole magnetic moment (DMM) are calculated. The cases of the Majorana and Dirac neutrino are considered. It is shown that to account for the $(g-2)_{\mu}$ anomaly DMMs connected with heavy neutrinos should be on the order of $\mbox{few}\times10^{-8}\mu_B$, that is at variance with the theoretically predicted values.


title: Does our universe conform with the existence of a universal maximum energy-density $\rho^{uni}_{max}$ ?

abstract: Recent astronomical observations of high redshift quasars, dark matter-dominated galaxies, mergers of neutron stars, glitch phenomena in pulsars, cosmic microwave background and experimental data from hadronic colliders do not rule out, but they even support the hypothesis that the energy-density in our universe most likely is upper-limited by $\rho^{uni}_{max},$ which is predicted to lie between $2$ to $3$ the nuclear density $\rho_0.$ Quantum fluids in the cores of massive NSs with $\rho \approx \rho^{uni}_{max}$ reach the maximum compressibility state, where they become insensitive to further compression by the embedding spacetime and undergo a phase transition into the purely incompressible gluon-quark superfluid state. A direct correspondence between the positive energy stored in the embedding spacetime and the degree of compressibility and superfluidity of the trapped matter is proposed. In this paper relevant observation signatures that support the maximum density hypothesis are reviewed, a possible origin of $\rho^{uni}_{max}$ is proposed and finally the consequences of this scenario on the spacetime's topology of the universe as well as on the mechanisms underlying the growth rate and power of the high redshift QSOs are discussed.

categories: astro-ph.GA astro-ph.HE cond-mat.supr-con gr-qc hep-ph hep-th


title: Towards Single Atom Computing via High Harmonic Generation

abstract: The development of alternative platforms for computing has been a longstanding goal for physics, and represents a particularly pressing concern as conventional transistors approach the limit of miniaturization. A potential alternatice paradigm is that of reservoir computing, which leverages unknown, but highly non-linear transformations of input-data to perform computations. This has the advantage that many physical systems exhibit precisely the type of non-linear input-output relationships necessary for them to function as reservoirs. Consequently, the quantum effects which obstruct the further development of silicon electronics become an advantage for a reservoir computer. Here we demonstrate that even the most basic constituents of matter - atoms - can act as a reservoir for optical computers, thanks to the phenomenon of High Harmonic Generation (HHG). A prototype single-atom computer for classification problems is proposed, where parameters of the classification model are mapped to optical elements. We numerically demonstrate that this `all-optical' computer can successfully classify data with an accuracy that is strongly dependent on dynamical non-linearities. This may pave the way for the development of petahertz information processing platforms.


title: {\delta}-CLUE: Diverse Sets of Explanations for Uncertainty Estimates

abstract: To interpret uncertainty estimates from differentiable probabilistic models, recent work has proposed generating Counterfactual Latent Uncertainty Explanations (CLUEs). However, for a single input, such approaches could output a variety of explanations due to the lack of constraints placed on the explanation. Here we augment the original CLUE approach, to provide what we call $\delta$-CLUE. CLUE indicates $\it{one}$ way to change an input, while remaining on the data manifold, such that the model becomes more confident about its prediction. We instead return a $\it{set}$ of plausible CLUEs: multiple, diverse inputs that are within a $\delta$ ball of the original input in latent space, all yielding confident predictions.


title: Modeling the dynamics of language change: logistic regression, Piotrowski's law, and a handful of examples in Polish

abstract: The study discusses modeling diachronic processes by logistic regression. The phenomenon of nonlinear changes in language was first observed by Raimund Piotrowski (hence labelled as Piotrowski's law), even if actual linguistic evidence usually speaks against using the notion of a "law" in this context. In our study, we apply logistic regression models to 9 changes which occurred between 15th and 18th century in the Polish language. The attested course of the majority of these changes closely follow the expected values, which proves that the language change might indeed resemble a nonlinear phase change scenario. We also extend the original Piotrowski's approach by proposing polynomial logistic regression for these cases which can hardly be described by its standard version. Also, we propose to consider individual language change cases jointly, in order to inspect their possible collinearity or, more likely, their different dynamics in the function of time. Last but not least, we evaluate our results by testing the influence of the subcorpus size on the model's goodness-of-fit.


title: Finding Concept-specific Biases in Form--Meaning Associations

abstract: This work presents an information-theoretic operationalisation of cross-linguistic non-arbitrariness. It is not a new idea that there are small, cross-linguistic associations between the forms and meanings of words. For instance, it has been claimed (Blasi et al., 2016) that the word for "tongue" is more likely than chance to contain the phone [l]. By controlling for the influence of language family and geographic proximity within a very large concept-aligned, cross-lingual lexicon, we extend methods previously used to detect within language non-arbitrariness (Pimentel et al., 2019) to measure cross-linguistic associations. We find that there is a significant effect of non-arbitrariness, but it is unsurprisingly small (less than 0.5% on average according to our information-theoretic estimate). We also provide a concept-level analysis which shows that a quarter of the concepts considered in our work exhibit a significant level of cross-linguistic non-arbitrariness. In sum, the paper provides new methods to detect cross-linguistic associations at scale, and confirms their effects are minor.


title: Terrain assessment for precision agriculture using vehicle dynamic modelling

abstract: Advances in precision agriculture greatly rely on innovative control and sensing technologies that allow service units to increase their level of driving automation while ensuring at the same time high safety standards. This paper deals with automatic terrain estimation and classification that is performed simultaneously by an agricultural vehicle during normal operations. Vehicle mobility and safety, and the successful implementation of important agricultural tasks including seeding, ploughing, fertilising and controlled traffic depend or can be improved by a correct identification of the terrain that is traversed. The novelty of this research lies in that terrain estimation is performed by using not only traditional appearance-based features, that is colour and geometric properties, but also contact-based features, that is measuring physics-based dynamic effects that govern the vehicleeterrain interaction and that greatly affect its mobility. Experimental results obtained from an all-terrain vehicle operating on different surfaces are presented to validate the system in the field. It was shown that a terrain classifier trained with contact features was able to achieve a correct prediction rate of 85.1%, which is comparable or better than that obtained with approaches using traditional feature sets. To further improve the classification performance, all feature sets were merged in an augmented feature space, reaching, for these tests, 89.1% of correct predictions.


title: Regarding the domain of non-symmetric and, possibly, degenerate Ornstein--Uhlenbeck operators in separable Banach spaces

abstract: Let $X$ be a separable Banach space and let $Q:X^*\rightarrow X$ be a linear, bounded, non-negative and symmetric operator and let $A:D(A)\subseteq X\rightarrow X$ be the infinitesimal generator of a strongly continuous semigroup of contractions on $X$. We consider the abstract Wiener space $(X,\mu_\infty,H_\infty)$ where $\mu_\infty$ is a centred non-degenerate Gaussian measure on $X$ with covariance operator defined, at least formally, as \begin{align*} Q_\infty=\int_0^{+\infty} e^{sA}Qe^{sA^*}ds, \end{align*} and $H_\infty$ is the Cameron--Martin space associated to $\mu_\infty$. Let $H$ be the reproducing kernel Hilbert space associated with $Q$ with inner product $[\cdot,\cdot]_H$. We assume that the operator $Q_\infty A^*:D(A^*)\subseteq X^*\rightarrow X$ extends to a bounded linear operator $B\in \mathcal L(H)$ which satisfies $B+B^*=-{\rm Id}_H$, where ${\rm Id}_H$ denotes the identity operator on $H$. Let $D$ and $D^2$ be the first and second order Fr\'echet derivative operators, we denote by $D_H$ and $D^2_H$ the closure in $L^2(X,\mu_\infty)$ of the operators $QD$ and $QD^2$ and by $W^{1,2}_H(X,\mu_\infty)$ and and $W^{2,2}_H(X,\mu_\infty)$ their domains in $L^2(X,\mu_\infty)$, respectively,. Furthermore, we denote by $D_{A_\infty}$ the closure of the operator $Q_\infty A^*D$ and by $W^{1,2}_{A_\infty}(X,\mu_\infty)$ its domain in $L^2(X,\mu_\infty)$. We characterize the domain of the operator $L$, associated to the bilinear form \begin{align*} (u,v)\mapsto-\int_{X}[BD_Hu,D_Hv]_Hd\mu_\infty, \qquad u,v\in W^{1,2}_H(X,\mu_\infty), \end{align*} in $L^2(X,\mu_\infty)$. More precisely, we prove that $D(L)$ coincides, up to an equivalent remorming, with a subspace of $W^{2,2}_H(X,\mu_\infty)\cap W^{1,2}_{A_\infty}(X,\mu_\infty)$. We stress that we are able to treat the case when $L$ is degenerate and non-symmetric.


title: List Message Passing Decoding of Non-binary Low-Density Parity-Check Codes

abstract: A decoding algorithm for $q$-ary low-density parity-check codes over the $q$-ary symmetric channel is introduced. The exchanged messages are lists of symbols from $\Fq$. A density evolution analysis for maximum list sizes $1$ and $2$ is developed. Thresholds for selected regular low-density parity-check code ensembles are computed showing gains with respect to a similar algorithm in the literature. Finite-length simulation results confirm the asymptotic analysis.


title: Strongly zero product determined Banach algebras

abstract: $C^*$-algebras, group algebras, and the algebra $\mathcal{A}(X)$ of approximable operators on a Banach space $X$ having the bounded approximation property are known to be zero product determined. We are interested in giving a quantitative estimate of this property by finding, for each Banach algebra $A$ of the above classes, a constant $\alpha$ with the property that for every continuous bilinear functional $\varphi\colon A \times A\to\mathbb{C}$ there exists a continuous linear functional $\xi$ on $A$ such that \[ \sup_{\Vert a\Vert=\Vert b\Vert=1}\vert\varphi(a,b)-\xi(ab)\vert\le \alpha\sup_{\mathclap{\substack{\Vert a\Vert=\Vert b\Vert=1, \\ ab=0}}}\vert\varphi(a,b)\vert. \]


title: On the stable eigenvalues of perturbed anharmonic oscillators in dimension two

abstract: We study the asymptotic behavior of the spectrum of a quantum system which is a perturbation of a spherically symmetric anharmonic oscillator in dimension 2. We prove that a large part of its eigenvalues can be obtained by Bohr-Sommerfeld quantization rule applied to the normal form Hamiltonian and also admit an asymptotic expansion at infinity. The proof is based on the generalization to the present context of the normal form approach developed in [BLM20b] (see also [PS10]) for the particular case of $T^d$.


title: Non-Debye relaxations: smeared time evolution, memory effects, and the Laplace exponents

abstract: The non-Debye, \textit{i.e.,} non-exponential, behavior characterizes a large plethora of dielectric relaxation phenomena. Attempts to find their theoretical explanation are dominated either by considerations rooted in the stochastic processes methodology or by the so-called \textsl{fractional dynamics} based on equations involving fractional derivatives which mimic the non-local time evolution and as such may be interpreted as describing memory effects. Using the recent results coming from the stochastic approach we link memory functions with the Laplace (characteristic) exponents of infinitely divisible probability distributions and show how to relate the latter with experimentally measurable spectral functions characterizing relaxation in the frequency domain. This enables us to incorporate phenomenological knowledge into the evolution laws. To illustrate our approach we consider the standard Havriliak-Negami and Jurlewicz-Weron-Stanislavsky models for which we derive well-defined evolution equations. Merging stochastic and fractional dynamics approaches sheds also new light on the analysis of relaxation phenomena which description needs going beyond using the single evolution pattern. We determine sufficient conditions under which such description is consistent with general requirements of our approach.


title: Lunar Dust Particles Blown By Lander Engine Exhaust in Rarefied and Compressible Flow

abstract: Previously we presented a numerical model that predicts trajectories of lunar dust, soil, and gravel blown by the engine exhaust of a lunar lander. The model uses the gas density, velocity vector field, and temperature predicted by computational fluid dynamics (CFD) or Direct Simulation Monte Carlo (DSMC) simulations to compute the forces and accelerations acting on the regolith particles, one particle at a time (ignoring particle collisions until more advanced models are developed). Here we present significant improvements to the model, including the implementation of particle drag and lift formulas to account for the rarefaction and compressibility of the flow. It turns out that the drag force is reduced due to the rarefaction, but the lift is increased due to several effects such as particle rotation. A data matrix of particle sizes, engine thrusts (descent and ascent values for Altair), horizontal and vertical starting distances, and lander height above ground, have been tested using the latest version of the software. These results suggest that the previously reported 3 degree trajectory angle limit can be exceeded in several cases by as much as a factor of five. Particles that originate at a height of 1 cm above the surface from an outer crater rim can be propelled to angles of 5 degrees or more. Particles that start at 10 cm above the surface can be ejected with angles of up to 15 degrees. Mechanisms responsible for placing particles at starting heights above the surface may include the kinetics of horizontal collisions, as suggested by Discrete Element Method (DEM) simulations. We also present results showing the distance particles travel and their impact velocities. We then use the model to evaluate the effectiveness of berms or other methods to block the spray of soil at a lunar landing site.


title: Decomposing hypergraphs into cycle factors

abstract: A famous result by R\"odl, Ruci\'nski, and Szemer\'edi guarantees a (tight) Hamilton cycle in $k$-uniform hypergraphs $H$ on $n$ vertices with minimum $(k-1)$-degree $\delta_{k-1}(H)\geq (1/2+o(1))n$, thereby extending Dirac's result from graphs to hypergraphs. For graphs, much more is known; each graph on $n$ vertices with $\delta(G)\geq (1/2+o(1))n$ contains $(1-o(1))r$ edge-disjoint Hamilton cycles where $r$ is the largest integer such that $G$ contains a spanning $2r$-regular subgraph, which is clearly asymptotically optimal. This was proved by Ferber, Krivelevich, and Sudakov answering a question raised by K\"uhn, Lapinskas, and Osthus. We extend this result to hypergraphs; every $k$-uniform hypergraph $H$ on $n$ vertices with $\delta_{k-1}(H)\geq (1/2+o(1))n$ contains $(1-o(1))r$ edge-disjoint (tight) Hamilton cycles where $r$ is the largest integer such that $H$ contains a spanning subgraph with each vertex belonging to $kr$ edges. In particular, this yields an asymptotic solution to a question of Glock, K\"uhn, and Osthus. In fact, our main result applies to approximately vertex-regular $k$-uniform hypergraphs with a weak quasirandom property and provides approximate decompositions into cycle factors without too short cycles.


title: On the asymptotic behavior of cocycles over flows

abstract: In 1968, V.I. Oseledets formulated the question of convergence in the Birkhoff theorem and the multiplicative ergodic theorem for measurable cocycles over flows under the condition of integrability for each individual t. A.M. Stepin and the author established (2016) the convergence along subsets of time of density 1. In this note, we show that moreover the convergence is fulfilled modulo time subsets of finite measure.


title: On the Use of Linguistic Features for the Evaluation of Generative Dialogue Systems

abstract: Automatically evaluating text-based, non-task-oriented dialogue systems (i.e., `chatbots') remains an open problem. Previous approaches have suffered challenges ranging from poor correlation with human judgment to poor generalization and have often required a gold standard reference for comparison or human-annotated data. Extending existing evaluation methods, we propose that a metric based on linguistic features may be able to maintain good correlation with human judgment and be interpretable, without requiring a gold-standard reference or human-annotated data. To support this proposition, we measure and analyze various linguistic features on dialogues produced by multiple dialogue models. We find that the features' behaviour is consistent with the known properties of the models tested, and is similar across domains. We also demonstrate that this approach exhibits promising properties such as zero-shot generalization to new domains on the related task of evaluating response relevance.


title: Real-time interactive 4D-STEM phase-contrast imaging from electron event representation data

abstract: The arrival of direct electron detectors (DED) with high frame-rates in the field of scanning transmission electron microscopy has enabled many experimental techniques that require collection of a full diffraction pattern at each scan position, a field which is subsumed under the name four dimensional-scanning transmission electron microscopy (4D-STEM). DED frame rates approaching 100 kHz require data transmission rates and data storage capabilities that exceed commonly available computing infrastructure. Current commercial DEDs allow the user to make compromises in pixel bit depth, detector binning or windowing to reduce the per-frame file size and allow higher frame rates. This change in detector specifications requires decisions to be made before data acquisition that may reduce or lose information that could have been advantageous during data analysis. The 4D Camera, a DED with 87 kHz frame-rate developed at Lawrence Berkeley National Laboratory, reduces the raw data to a linear-index encoded electron event representation (EER). Here we show with experimental data from the 4D Camera that linear-index encoded EER and its direct use in 4D-STEM phase contrast imaging methods enables real-time, interactive phase-contrast from large-area 4D-STEM datasets. We detail the computational complexity advantages of the EER and the necessary computational steps to achieve real-time interactive ptychography and center-of-mass differential phase contrast using commonly available hardware accelerators.


title: Negative Differential Resistance in Carbon-Based Nanostructures

abstract: Nonlinear electrical properties, such as negative differential resistance (NDR), are essential in numerous electrical circuits, including memristors. Several physical origins have been proposed to lead to the NDR phenomena in semiconductor devices in the last more than half a century. Here, we report NDR behavior formation in randomly oriented graphene-like nanostructures up to 37 K and high on-current density up to 10^5 A/cm^2. Our modeling of the current-voltage characteristics, including the self-heating effects, suggests that strong temperature dependence of the low-bias resistance is responsible for the nonlinear electrical behavior. Our findings open opportunities for the practical realization of the on-demand NDR behavior in nanostructures of 2D and 3D material-based devices via heat management in the conducting films and the underlying substrates.


title: On the Impact of Random Seeds on the Fairness of Clinical Classifiers

abstract: Recent work has shown that fine-tuning large networks is surprisingly sensitive to changes in random seed(s). We explore the implications of this phenomenon for model fairness across demographic groups in clinical prediction tasks over electronic health records (EHR) in MIMIC-III -- the standard dataset in clinical NLP research. Apparent subgroup performance varies substantially for seeds that yield similar overall performance, although there is no evidence of a trade-off between overall and subgroup performance. However, we also find that the small sample sizes inherent to looking at intersections of minority groups and somewhat rare conditions limit our ability to accurately estimate disparities. Further, we find that jointly optimizing for high overall performance and low disparities does not yield statistically significant improvements. Our results suggest that fairness work using MIMIC-III should carefully account for variations in apparent differences that may arise from stochasticity and small sample sizes.


title: Deep imagination is a close to optimal policy for planning in large decision trees under limited resources

abstract: Many decisions involve choosing an uncertain course of actions in deep and wide decision trees, as when we plan to visit an exotic country for vacation. In these cases, exhaustive search for the best sequence of actions is not tractable due to the large number of possibilities and limited time or computational resources available to make the decision. Therefore, planning agents need to balance breadth (exploring many actions at each level of the tree) and depth (exploring many levels in the tree) to allocate optimally their finite search capacity. We provide efficient analytical solutions and numerical analysis to the problem of allocating finite sampling capacity in one shot to large decision trees. We find that in general the optimal policy is to allocate few samples per level so that deep levels can be reached, thus favoring depth over breadth search. In contrast, in poor environments and at low capacity, it is best to broadly sample branches at the cost of not sampling deeply, although this policy is marginally better than deep allocations. Our results provide a theoretical foundation for the optimality of deep imagination for planning and show that it is a generally valid heuristic that could have evolved from the finite constraints of cognitive systems.


title: Comments on diffraction experiments at high-field on uranium dioxide in the ordered state

abstract: D. J. Antonio et al. report an x-ray diffraction experiment on uranium dioxide at high-magnetic field and low temperature. The authors have apparently not realized that the diffraction data actually shows unambiguously the presence of a rhombohedral distortion induced by the field at low temperature and the presence of more than one domain. Our note correctly analyses their data. The crystallographic analysis we perform is based on simple arguments with respect to plane spacings in slightly distorted cubic materials. Similar studies have been reported in many materials over the last 50 years. Our analysis explains quantitatively the observations, the presence of the two peaks (the appearance of which the authors regard as unexpected) and their different variations with respect to the applied field, which is a simple consequence of the conservation of atomic volume.


title: Constraint-coupled Optimization with Unknown Costs: A Distributed Primal Decomposition Approach

abstract: In this paper, we present a distributed algorithm for solving convex, constraint-coupled, optimization problems over peer-to-peer networks. We consider a network of processors that aim to cooperatively minimize the sum of local cost functions, subject to individual constraints and to global coupling constraints. The major assumption of this work is that the cost functions are unknown and must be learned online. We propose a fully distributed algorithm, based on a primal decomposition approach, that uses iteratively refined data-driven estimations of the cost functions over the iterations. The algorithm is scalable and maintains private information of agents. We prove that, asymptotically, the distributed algorithm provides the optimal solution of the problem even though the true cost functions are never used within the algorithm. The analysis requires an in-depth exploration of the primal decomposition approach and shows that the distributed algorithm can be thought of as an epsilon-subgradient method applied to a suitable reformulation of the original problem. Finally, numerical computations cor- roborate the theoretical findings and show the efficacy of the proposed approach.