title: Reduced-Space Interior Point Methods in Power Grid Problems

abstract: Due to critical environmental issues, the power systems have to accommodate a significant level of penetration of renewable generation which requires smart approaches to the power grid control. Associated optimal control problems are large-scale nonlinear optimization problems with up to hundreds of millions of variables and constraints. The interior point methods become computationally intractable, mainly due to the solution of large linear systems. This document addresses the computational bottlenecks of the interior point method during the solution of the security constrained optimal power flow problems by applying reduced space quasi-Newton IPM, which could utilize high-performance computers due to the inherent parallelism in the adjoint method. Reduced space IPM approach and the adjoint method is a novel approach when it comes to solving the (security constrained) optimal power flow problems. These were previously used in the PDE-constrained optimization. The presented methodology is suitable for high-performance architectures due to inherent parallelism in the adjoint method during the gradient evaluation, since the individual contingency scenarios are modeled by independent set of the constraints. Preliminary evaluation of the performance and convergence is performed to study the reduced space approach.


title: Multi-task Learning for Speaker Verification and Voice Trigger Detection

abstract: Automatic speech transcription and speaker recognition are usually treated as separate tasks even though they are interdependent. In this study, we investigate training a single network to perform both tasks jointly. We train the network in a supervised multi-task learning setup, where the speech transcription branch of the network is trained to minimise a phonetic connectionist temporal classification (CTC) loss while the speaker recognition branch of the network is trained to label the input sequence with the correct label for the speaker. We present a large-scale empirical study where the model is trained using several thousand hours of labelled training data for each task. We evaluate the speech transcription branch of the network on a voice trigger detection task while the speaker recognition branch is evaluated on a speaker verification task. Results demonstrate that the network is able to encode both phonetic \emph{and} speaker information in its learnt representations while yielding accuracies at least as good as the baseline models for each task, with the same number of parameters as the independent models.


title: MCSAE: Masked Cross Self-Attentive Encoding for Speaker Embedding

abstract: In general, a self-attention mechanism has been applied for speaker embedding encoding. Previous studies focused on training the self-attention in a high-level layer, such as the last pooling layer. However, the effect of low-level features was reduced in the speaker embedding encoding. Therefore, we propose masked cross self-attentive encoding (MCSAE) using ResNet. It focuses on the features of both high-level and lowlevel layers. Based on multi-layer aggregation, the output features of each residual layer are used for the MCSAE. In the MCSAE, cross self-attention module is trained the interdependence of each input features. A random masking regularization module also applied to preventing overfitting problem. As such, the MCSAE enhances the weight of frames representing the speaker information. Then, the output features are concatenated and encoded to the speaker embedding. Therefore, a more informative speaker embedding is encoded by using the MCSAE. The experimental results showed an equal error rate of 2.63% and a minimum detection cost function of 0.1453 using the VoxCeleb1 evaluation dataset. These were improved performances compared with the previous self-attentive encoding and state-of-the-art encoding methods.


title: Convergence Guarantees for Gaussian Process Means With Misspecified Likelihoods and Smoothness

abstract: Gaussian processes are ubiquitous in machine learning, statistics, and applied mathematics. They provide a flexible modelling framework for approximating functions, whilst simultaneously quantifying uncertainty. However, this is only true when the model is well-specified, which is often not the case in practice. In this paper, we study the properties of Gaussian process means when the smoothness of the model and the likelihood function are misspecified. In this setting, an important theoretical question of practial relevance is how accurate the Gaussian process approximations will be given the difficulty of the problem, our model and the extent of the misspecification. The answer to this problem is particularly useful since it can inform our choice of model and experimental design. In particular, we describe how the experimental design and choice of kernel and kernel hyperparameters can be adapted to alleviate model misspecification.

categories: math.ST cs.LG cs.NA math.NA stat.ML stat.TH


title: Causal discrete field theory for quantum gravity

abstract: The proposed theory of causally structured discrete fields studies integer values on directed edges of a self-similar graph with a propagation rule, which we define as a set of valid combinations of integer values and edge directions around any vertex of the graph. There is an infinite countable number of variants of the theory for a given self-similar graph depending on the choice of propagation rules, some of these models can generate infinite uncountable sets of patterns. This theory takes minimum assumptions of causality, discreteness, locality, and determinism, as well as fundamental symmetries of isotropy, CPT invariance, and charge conservation. It combines the elements of cellular automata, causal sets, loop quantum gravity, and causal dynamical triangulations to become an excellent candidate to describe quantum gravity at the Planck scale. In addition to the self-consistent generation of spacetime and metrics to describe gravity and an expanding closed Universe, the theory allows for the many-worlds interpretation of quantum mechanics. We also demonstrate how to get to unitary evolution in Hilbert space for a stationary Universe with deterministic propagation.


title: Reproducibility Challenge NeurIPS 2019 Report on "Competitive Gradient Descent"

abstract: This is a report for reproducibility challenge of NeurlIPS 2019 on the paper Competitive Gradient Descent (Schafer et al., 2019). The paper introduces a novel algorithm for the numerical computation of Nash equilibria of competitive two-player games. It avoids oscillatory and divergent behaviours seen in alternating gradient descent. The purpose of this report is to critically examine the reproducibility of the work by (Schafer et al., 2019), within the framework of the NeurIPS 2019 Reproducibility Challenge. The experiments replicated in this report confirms the results of the original study. Moreover, this project offers a Python (Pytorch based) implementation of the proposed CGD algorithm which can be found at the following public git repository: (https://github.com/GopiKishan14/Reproducibility_Challenge_NeurIPS_2019)


title: Decremental SSSP in Weighted Digraphs: Faster and Against an Adaptive Adversary

abstract: Given a dynamic digraph $G = (V,E)$ undergoing edge deletions and given $s\in V$ and $\epsilon>0$, we consider the problem of maintaining $(1+\epsilon)$-approximate shortest path distances from $s$ to all vertices in $G$ over the sequence of deletions. Even and Shiloach (J.~ACM'$81$) give a deterministic data structure for the exact version of the problem with total update time $O(mn)$. Henzinger et al. (STOC'$14$, ICALP'$15$) give a Monte Carlo data structure for the approximate version with improved total update time $ O(mn^{0.9 + o(1)}\log W)$ where $W$ is the ratio between the largest and smallest edge weight. A drawback of their data structure is that they only work against an oblivious adversary, meaning that the sequence of deletions needs to be fixed in advance. This limits its application as a black box inside algorithms. We present the following $(1+\epsilon)$-approximate data structures: (1) the first data structure is Las Vegas and works against an adaptive adversary; it has total expected update time $\tilde O(m^{2/3}n^{4/3})$ for unweighted graphs and $\tilde O(m^{3/4}n^{5/4}\log W)$ for weighted graphs, (2) the second data structure is Las Vegas and assumes an oblivious adversary; it has total expected update time $\tilde O(\sqrt m n^{3/2})$ for unweighted graphs and $\tilde O(m^{2/3}n^{4/3}\log W)$ for weighted graphs, (3) the third data structure is Monte Carlo and is correct w.h.p.~against an oblivious adversary; it has total expected update time $\tilde O((mn)^{7/8}\log W) = \tilde O(mn^{3/4}\log W)$. Each of our data structures can be queried at any stage of $G$ in constant worst-case time; if the adversary is oblivious, a query can be extended to also report such a path in time proportional to its length. Our update times are faster than those of Henzinger et al.~for all graph densities.


title: Lattice-based Improvements for Voice Triggering Using Graph Neural Networks

abstract: Voice-triggered smart assistants often rely on detection of a trigger-phrase before they start listening for the user request. Mitigation of false triggers is an important aspect of building a privacy-centric non-intrusive smart assistant. In this paper, we address the task of false trigger mitigation (FTM) using a novel approach based on analyzing automatic speech recognition (ASR) lattices using graph neural networks (GNN). The proposed approach uses the fact that decoding lattice of a falsely triggered audio exhibits uncertainties in terms of many alternative paths and unexpected words on the lattice arcs as compared to the lattice of a correctly triggered audio. A pure trigger-phrase detector model doesn't fully utilize the intent of the user speech whereas by using the complete decoding lattice of user audio, we can effectively mitigate speech not intended for the smart assistant. We deploy two variants of GNNs in this paper based on 1) graph convolution layers and 2) self-attention mechanism respectively. Our experiments demonstrate that GNNs are highly accurate in FTM task by mitigating ~87% of false triggers at 99% true positive rate (TPR). Furthermore, the proposed models are fast to train and efficient in parameter requirements.


title: A note on causation versus correlation

abstract: Recently, it has been shown that the causality and information flow between two time series can be inferred in a rigorous and quantitative sense, and, besides, the resulting causality can be normalized. A corollary that follows is, in the linear limit, causation implies correlation, while correlation does not imply causation. Now suppose there is an event $A$ taking a harmonic form (sine/cosine), and it generates through some process another event $B$ so that $B$ always lags $A$ by a phase of $\pi/2$. Here the causality is obviously seen, while by computation the correlation is, however, zero. This seemingly contradiction is rooted in the fact that a harmonic system always leaves a single point on the Poincar\'e section; it does not add information. That is to say, though the absolute information flow from $A$ to $B$ is zero, i.e., $T_{A\to B}=0$, the total information increase of $B$ is also zero, so the normalized $T_{A\to B}$, denoted as $\tau_{A\to B}$, takes the form of $\frac 0 0$. By slightly perturbating the system with some noise, solving a stochastic differential equation, and letting the perturbation go to zero, it can be shown that $\tau_{A\to B}$ approaches 100\%, just as one would have expected.


title: The Clouds in Asynchronous Cellular Automata

abstract: This article introduces the notion of clouds in asynchronous cellular automata (ACAs). We show that the cloud behaviour of ACAs has similarity with natural clouds across the sky, election model of parliamentary democratic system, and electron cloud around nucleus. These systems, therefore, can be modelled by the ACAs.


title: Slow Fluors for Effective Separation of Cherenkov Light in Liquid Scintillators

abstract: The timing and spectral characteristics of four highly efficient, slow fluors are presented for liquid scintillator solutions using linear alkylbenzene (LAB) as the primary solvent. The mixtures exhibit high light yields, but with rise times of several ns or more and decay times on the order of tens of ns. Consequently, such liquid scintillator mixtures can be used for effective separation of Cherenkov and scintillation components based on timing in large scale liquid scintillation detectors. Such a separation, showing high light yield and directional information, is demonstrated here on a bench-top scale for electrons with energies extending below 1 MeV. This could have significant consequences for the future development of such detectors for measurements of solar neutrinos and neutrinoless double beta decay as well as providing good directional information for elastic scattering events from supernovae neutrinos and reactor anti-neutrinos, amongst others.


title: Singlets in the tensor product of an arbitrary number of Adjoint representations of SU(3)

abstract: We propose a set of 4 recurrence relations whose linear combination gives the number of group invariants, equivalently the dimension of the invariant subspace, in the tensor product of an arbitrary number of adjoint representations of the SU(3) Lie Group.


title: Stochastic Analysis of an Adaptive Cubic Regularisation Method under Inexact Gradient Evaluations and Dynamic Hessian Accuracy

abstract: We here adapt an extended version of the adaptive cubic regularisation method with dynamic inexact Hessian information for nonconvex optimisation in [3] to the stochastic optimisation setting. While exact function evaluations are still considered, this novel variant inherits the innovative use of adaptive accuracy requirements for Hessian approximations introduced in [3] and additionally employs inexact computations of the gradient. Without restrictions on the variance of the errors, we assume that these approximations are available within a sufficiently large, but fixed, probability and we extend, in the spirit of [18], the deterministic analysis of the framework to its stochastic counterpart, showing that the expected number of iterations to reach a first-order stationary point matches the well known worst-case optimal complexity. This is, in fact, still given by O(epsilon^(-3/2)), with respect to the first-order epsilon tolerance. Finally, numerical tests on nonconvex finite-sum minimisation confirm that using inexact first and second-order derivatives can be beneficial in terms of the computational savings.


title: A New Arc-Routing Algorithm Applied to Winter Road Maintenance

abstract: This paper studies large scale instances of a fairly general arc-routing problem as well as incorporate practical constraints in particular coming from the scheduling problem of the winter road maintenance (e.g. different priorities for and methods of road maintenance). We develop a new algorithm based on a bin-packing heuristic which is well-scalable and able to solve road networks on thousands of crossroads and road segments in few minutes. Since it is impossible to find an optimal solution for such a large instances to compare it with a result of our algorithm, we also develop techniques to compute lower bounds which are based on Integer Linear Programming and Lazy Constraints.


title: Variational Autoencoders for Opponent Modeling in Multi-Agent Systems

abstract: Multi-agent systems exhibit complex behaviors that emanate from the interactions of multiple agents in a shared environment. In this work, we are interested in controlling one agent in a multi-agent system and successfully learn to interact with the other agents that have fixed policies. Modeling the behavior of other agents (opponents) is essential in understanding the interactions of the agents in the system. By taking advantage of recent advances in unsupervised learning, we propose modeling opponents using variational autoencoders. Additionally, many existing methods in the literature assume that the opponent models have access to opponent's observations and actions during both training and execution. To eliminate this assumption, we propose a modification that attempts to identify the underlying opponent model using only local information of our agent, such as its observations, actions, and rewards. The experiments indicate that our opponent modeling methods achieve equal or greater episodic returns in reinforcement learning tasks against another modeling method.


title: Medical image reconstruction with image-adaptive priors learned by use of generative adversarial networks

abstract: Medical image reconstruction is typically an ill-posed inverse problem. In order to address such ill-posed problems, the prior distribution of the sought after object property is usually incorporated by means of some sparsity-promoting regularization. Recently, prior distributions for images estimated using generative adversarial networks (GANs) have shown great promise in regularizing some of these image reconstruction problems. In this work, we apply an image-adaptive GAN-based reconstruction method (IAGAN) to reconstruct high fidelity images from incomplete medical imaging data. It is observed that the IAGAN method can potentially recover fine structures in the object that are relevant for medical diagnosis but may be oversmoothed in reconstructions with traditional sparsity-promoting regularization.


title: The long time behavior and the rate of convergence of symplectic convex algorithms obtained via splitting discretizations of inertial damping systems

abstract: In this paper we propose new numerical algorithms in the setting of unconstrained optimization problems and we study the rate of convergence in the iterates of the objective function. Furthermore, our algorithms are based upon splitting and symplectic methods and they preserve the energy properties of the inherent continuous dynamical system that contains a Hessian perturbation. At the same time, we show that Nesterov gradient method is equivalent to a Lie-Trotter splitting applied to a Hessian driven damping system. Finally, some numerical experiments are presented in order to validate the theoretical results.


title: Audio-Visual Decision Fusion for WFST-based and seq2seq Models

abstract: Under noisy conditions, speech recognition systems suffer from high Word Error Rates (WER). In such cases, information from the visual modality comprising the speaker lip movements can help improve the performance. In this work, we propose novel methods to fuse information from audio and visual modalities at inference time. This enables us to train the acoustic and visual models independently. First, we train separate RNN-HMM based acoustic and visual models. A common WFST generated by taking a special union of the HMM components is used for decoding using a modified Viterbi algorithm. Second, we train separate seq2seq acoustic and visual models. The decoding step is performed simultaneously for both modalities using shallow fusion while maintaining a common hypothesis beam. We also present results for a novel seq2seq fusion without the weighing parameter. We present results at varying SNR and show that our methods give significant improvements over acoustic-only WER.


title: On quantum ensembles of quantum classifiers

abstract: Quantum machine learning seeks to exploit the underlying nature of a quantum computer to enhance machine learning techniques. A particular framework uses the quantum property of superposition to store sets of parameters, thereby creating an ensemble of quantum classifiers that may be computed in parallel. The idea stems from classical ensemble methods where one attempts to build a stronger model by averaging the results from many different models. In this work, we demonstrate that a specific implementation of the quantum ensemble of quantum classifiers, called the accuracy-weighted quantum ensemble, can be fully dequantised. On the other hand, the general quantum ensemble framework is shown to contain the well-known Deutsch-Jozsa algorithm that notably provides a quantum speedup and creates the potential for a useful quantum ensemble to harness this computational advantage.


title: Smart Induction for Isabelle/HOL (System Description)

abstract: Proof assistants offer tactics to facilitate inductive proofs. However, it still requires human ingenuity to decide what arguments to pass to those induction tactics. To automate this process, we present smart_induct for Isabelle/HOL. Given an inductive problem in any problem domain, smart_induct lists promising arguments for the induct tactic without relying on a search. Our evaluation demonstrated smart_induct produces valuable recommendations across problem domains.


title: Modeling single-phase permeability in uniform grain packs

abstract: Accurate estimation of single-phase permeability (k) has broad application in numerous areas, particularly modeling flow and transport in porous materials. Various techniques have been proposed in the literature to estimate k from other medium's properties, such as porosity, grain and/or pore size distribution, and pore connectivity. Among them critical path analysis (CPA) from statistical physics, first developed to model fluid flow in media with broad conductance distributions, has been successfully applied to heterogeneous soils and rocks. However, its application to uniform sphere and/or glass bead packs that represent homogeneous porous media with narrow conductance distributions needs to be investigated. In this study, we invoke concepts from CPA and estimate k from average grain diameter and formation factor in uniform sand and glass bead packs. By comparing theory with eight datasets including 105 packs from the literature, we demonstrate that CPA estimates permeability in homogeneous media accurately. We also compare our CPA-based model estimations with those from the Kozeny-Carman, Revil and Cathles, and RGPZ models. Results indicate that the CPA approach estimates k more precisely than other three models studied here.


title: On the Cooling Flow Problem in the Gaseous Halo of the Milky Way

abstract: Theoretical and observational arguments suggest that there is a large amount of hot ($\sim 10^6$ K), diffuse gas residing in the Milky Way's halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $\beta$ model. Typical mass inflow rates in our simulations range from $\sim 5 M_{\odot}$ yr$^{-1}$ to $\sim 60 M_{\odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies.


title: Measured Lightcurves and Rotational Periods of 3122 Florence, 3830 Trelleborg, and (131077) 2000 YH105

abstract: We determined the rotational periods of 3122 Florence, 3830 Trelleborg, and (131077) 2000 YH105 with the Harvard Clay Telescope and KeplerCam at the Fred L. Whipple Observatory. We found the rotational periods to be 2.3580 $\pm$ 0.0015 h, 17.059 $\pm$ 0.017 h, and 1.813 $\pm$ 0.00003 h, respectively. Our measurement of 3122 Florence's period agrees with Warner (2016), who reported 2.3580 $\pm$ 0.0002 h.


title: A stochastic output-feedback MPC scheme for distributed systems

abstract: In this paper, we present a novel stochastic output-feedback MPC scheme for distributed systems with additive process and measurement noise. The chance constraints are treated with the concept of probabilistic reachable sets, which, under an unimodality assumption on the disturbance distributions are guaranteed to be satisfied in closed-loop. By conditioning the initial state of the optimization problem on feasibility, the fundamental property of recursive feasibility is ensured. Closed-loop chance constraint satisfaction, recursive feasibility and convergence to an asymptotic average cost bound are proven. The paper closes with a numerical example of three interconnected subsystems, highlighting the chance constraint satisfaction and average cost compared to a centralized setting.


title: A kind of quaternary sequences of period $2p^mq^n$ and their linear complexity

abstract: Sequences with high linear complexity have wide applications in cryptography. In this paper, a new class of quaternary sequences over $\mathbb{F}_4$ with period $2p^mq^n$ is constructed using generalized cyclotomic classes. Results show that the linear complexity of these sequences attains the maximum.


title: The relevance of learning quantum physics from the perspective of the secondary school student: A case study

abstract: Studying quantum physics in upper secondary school is now a standard practice (Stadermann et al., 2019). But given the context of science education, with low recruitment numbers in higher education and poor attitudes towards science, it remains a question whether students find the learning of quantum physics relevant. In this study, we explore how students perceive the importance of quantum physics and technology and whether their perception changes after an intervention, namely the "Quantum Rules!" visit. We also aim to understand if they overall feel that learning quantum is relevant or not. In order to answer these questions, we followed a mixed-methods approach, combining both questionnaires and interviews. The quantitative analysis showed that the "Quantum Rules!" intervention has a positive effect on students' perception of the relevance of quantum physics and technology, especially regarding how important they feel quantum science is for society. Nevertheless, the qualitative information revealed that although students may find quantum physics and technology important for society, that does not necessarily mean that they find learning quantum physics relevant. We found that students believe the latter is relevant to them only if they find it interesting. We therefore rediscover the common expression "important, but not for me," and we further propose that this perception derives from students not seeing the societal relevance of learning quantum physics.


title: Evolution of quasiperiodic structures in non-ideal hydrodynamic description of phase transitions

abstract: Various phase transitions could have taken place in the early Universe, and may occur in the course of heavy-ion collisions and supernova explosions, in proto-neutron stars, cold compact stars, and in the condensed matter at terrestrial conditions. Most generally, the dynamics of the density and temperature at first- and second-order phase transitions can be described with the help of the equations of non-ideal hydrodynamics. In the given work some novel solutions are found describing the evolution of quasiperiodic structures that are formed in the course of the phase transitions. Although this consideration is very general, particular examples of quark-hadron and nuclear liquid-gas first-order phase transitions to the uniform k_0=0 state and pion condensate second-order phase transition to non-uniform k_0\neq 0 state in dense baryon matter are considered.


title: Predicting Large-Chern-Number Phases in a Shaken Optical Dice Lattice

abstract: With respect to the quantum anomalous Hall effect (QAHE), the detection of topological nontrivial large-Chern-number phases is an intriguing subject. Motivated by recent research on Floquet topological phases, this study proposes a periodic driving protocol to engineer large-Chern-number phases using QAHE. Herein, spinless ultracold fermionic atoms are studied in a two-dimensional optical dice lattice with nearest-neighbor hopping and a $\Lambda$/V-type sublattice potential subjected to a circular driving force. Results suggest that large-Chern-number phases exist with Chern numbers equal to $C=-2$, which is consistent with the edge-state energy spectra.


title: Finite time extinction for the 1D stochastic porous medium equation with transport noise

abstract: We establish finite time extinction with probability one for weak solutions of the Cauchy-Dirichlet problem for the 1D stochastic porous medium equation with Stratonovich transport noise and compactly supported smooth initial datum. Heuristically, this is expected to hold because Brownian motion has average spread rate $O(t^\frac{1}{2})$ whereas the support of solutions to the deterministic PME grows only with rate $O(t^{\frac{1}{m{+}1}})$. The rigorous proof relies on a contraction principle up to time-dependent shift for Wong-Zakai type approximations, the transformation to a deterministic PME with two copies of a Brownian path as the lateral boundary, and techniques from the theory of viscosity solutions.


title: Extremely slow non-equilibrium monopole dynamics in classical spin ice

abstract: We report on the non-equilibrium monopole dynamics in the classical spin ice Dy$_2$Ti$_2$O$_7$ detected by means of high-resolution magnetostriction measurements. Significant lattice changes occur at the transition from the kagome-ice to the saturated-ice phase, visible in the longitudinal and transverse magnetostriction. A hysteresis opening at temperatures below 0.6 K suggests a first-order transition between the kagome and saturated state. Extremely slow lattice relaxations, triggered by changes of the magnetic field, were observed. These lattice-relaxation effects result from non-equilibrium monopole formation or annihilation processes. The relaxation times extracted from our experiment are in good agreement with theoretical predictions with decay constants of the order of $10{^4}$ s at 0.3 K.


title: Reconfigurable Intelligent Surface for MISO Systems with Proportional Rate Constraints

abstract: This paper investigates the system spectral efficiency (SE) in reconfigurable intelligent surface (RIS)-aided multiuser multiple-input single-output (MISO) systems, where RIS can reconfigure the propagation environment via a large number of controllable and intelligent phase shifters. In order to explore the system SE performance behavior with user proportional fairness for such a system, an optimization problem is formulated to maximize the SE by jointly considering the power allocation at the base station (BS) and phase shift at the RIS, under nonlinear proportional rate fairness constraints. To solve the nonconvex optimization problem, an effective solution is developed, which capitalizes on an iterative algorithm with closed-form expressions, i.e., alternatively optimizing the transmit power at the BS and the reflecting phase shift at the RIS. Numerical simulations are provided to validate the theoretical analysis and assess the performance of the proposed alternative algorithm.

categories: cs.IT cs.NI cs.PF eess.SP math.IT math.OC


title: About Convergence and Order of Convergence of some Fractional Derivatives

abstract: In this paper we establish some convergence results for Riemann-Liouville, Caputo, and Caputo-Fabrizio fractional operators when the order of differentiation approaches one. We consider some errors given by $\left|\left| D^{1-\al}f -f'\right|\right|_p$ for p=1 and $p=\infty$ and we prove that for both Caputo and Caputo Fabrizio operators the order of convergence is a positive real r, 0<r<1. Finally, we compare the speed of convergence between Caputo and Caputo-Fabrizio operators obtaining that they a related by the Digamma function.


title: Nonlinear spline approximation in $\operatorname{BMO}(\mathbb R)$

abstract: We study nonlinear approximation in $\operatorname{BMO}$ from splines generated by a hierarchy of B-splines over regular multilevel nested partitions of $\mathbb R$. Companion Jackson and Bernstein estimates are established that allow to completely characterize the associated approximation spaces.


title: Counterintuitive example on relation between ZT and thermoelectric efficiency

abstract: The thermoelectric figure of merit ZT, which is defined using electrical conductivity, Seebeck coefficient, thermal conductivity, and absolute temperature T, has been widely used as a simple estimator of the conversion efficiency of a thermoelectric heat engine. When material properties are constant or slowly varying with T, a higher ZT ensures a higher maximum conversion efficiency of thermoelectric materials. However, as material properties can vary strongly with T, efficiency predictions based on ZT can be inaccurate, especially for wide-temperature applications. Moreover, although ZT values continue to increase, there has been no investigation of the relationship between ZT and the efficiency in the higher ZT regime. In this paper, we report a counterintuitive situation by comparing two materials: although one material has a higher ZT value over the whole operational temperature range, its maximum conversion efficiency is smaller than that of the other. This indicates that, for material comparisons, the evaluation of exact efficiencies as opposed to a simple comparison of the ZTs is necessary in certain cases.


title: Axion Constraints from Quiescent Soft Gamma-ray Emission from Magnetars

abstract: Axion-like-particles (ALPs) emitted from the core of a magnetar can convert to photons in its magnetosphere. The resulting photon flux is sensitive to the product of $(i)$ the ALP-nucleon coupling $G_{an}$ which controls the production cross section in the core and $(ii)$ the ALP-photon coupling $g_{a\gamma \gamma}$ which controls the conversion in the magnetosphere. We study such emissions in the soft-gamma-ray range (300 keV to 1 MeV), where the ALP spectrum peaks and astrophysical backgrounds from resonant Compton upscattering are expected to be suppressed. Using published quiescent soft-gamma-ray flux upper limits in 5 magnetars obtained with $CGRO$ COMPTEL and $INTEGRAL$ SPI/IBIS/ISGRI, we put limits on the product of the ALP-nucleon and ALP-photon couplings. We also provide a detailed study of the dependence of our results on the magnetar core temperature. We further show projections of our result for future $Fermi$-GBM observations. Our results motivate a program of studying quiescent soft-gamma-ray emission from magnetars with the $Fermi$-GBM.


title: A note on nonradial nodal solutions to the H\'enon problem in the disc

abstract: In this paper we consider some nodal solutions of the H\'enon problem in the unit disc with Dirichlet boundary conditions and we show that they are quasiradial, that is to say they are nonradial, they have two nodal regions and their nodal line does not touch the boundary of the disc.


title: Einselection from incompatible decoherence channels

abstract: Decoherence of quantum systems from entanglement with an unmonitored environment is to date the most compelling explanation of the emergence of a classical picture from a quantum world. While it is well understood for a single Lindblad operator, the role in the einselection process of a complex system-environment interaction remains to be clarified. In this paper, we analyze an open quantum dynamics inspired by CQED experiments with two non-commuting Lindblad operators modeling decoherence in the number basis and dissipative decoherence in the coherent state basis. We study and solve exactly the problem using quantum trajectories and phase-space techniques. The einselection optimization problem, which we consider to be about finding states minimizing the variation of some entanglement witness at a given energy, is studied numerically. We show that Fock states remain the most robust states to decoherence up to a critical coupling.


title: All-inclusive interacting dark sector cosmologies

abstract: In this paper we explore possible extensions of Interacting Dark Energy cosmologies, where Dark Energy and Dark Matter interact non-gravitationally with one another. In particular, we focus on the neutrino sector, analyzing the effect of both neutrino masses and the effective number of neutrino species. We consider the Planck 2018 legacy release data combined with several other cosmological probes, finding no evidence for new physics in the dark radiation sector. The current neutrino constraints from cosmology should be therefore regarded as robust, as they are not strongly dependent on the dark sector physics, once all the available observations are combined. Namely, we find a total neutrino mass $M_\nu<0.15$ eV and a number of effective relativistic degrees of freedom of $N_{\rm eff}=3.03^{+0.33}_{-0.33}$, both at 95\% CL, which are close to those obtained within the $\Lambda$CDM cosmology, $M_\nu<0.12$ eV and $N_{\rm eff}=3.00^{+0.36}_{-0.35}$ for the same data combination.


title: H-OWAN: Multi-distorted Image Restoration with Tensor 1x1 Convolution

abstract: It is a challenging task to restore images from their variants with combined distortions. In the existing works, a promising strategy is to apply parallel "operations" to handle different types of distortion. However, in the feature fusion phase, a small number of operations would dominate the restoration result due to the features' heterogeneity by different operations. To this end, we introduce the tensor 1x1 convolutional layer by imposing high-order tensor (outer) product, by which we not only harmonize the heterogeneous features but also take additional non-linearity into account. To avoid the unacceptable kernel size resulted from the tensor product, we construct the kernels with tensor network decomposition, which is able to convert the exponential growth of the dimension to linear growth. Armed with the new layer, we propose High-order OWAN for multi-distorted image restoration. In the numerical experiments, the proposed net outperforms the previous state-of-the-art and shows promising performance even in more difficult tasks.


title: Pressure fluctuations, viscosity, and Brownian motion

abstract: Brownian motion occurs in a variety of fluids, from rare gases to liquids. The Langevin equation, describing friction and agitation forces in statistical balance, is one of the most successful ways to treat the phenomenon. In rare gases, it is appropriate to model both friction and agitation in terms of independent molecular impacts with the particle. But in relatively dense fluids, such as water and air at standard temperature and pressure, the mean free path between collisions of fluid molecules is much smaller than the size of the Brownian particle, and the friction is normally treated as a mesoscopic viscous effect described by Stokes' Law which treats the fluid as continuous. The appropriateness of using Stokes' Law will be discussed in terms of recent experimental research in the ballistic or "coasting" phase of motion occurring at a very short time scale. Given the mesoscopic nature of the friction force for relatively dense fluids, we should expect the agitation force to also be mesoscopic. But it is often unrealistically modeled as uncorrelated individual impacts. It has been suggested occasionally that mesoscopic pressure fluctuations are appropriate for denser fluids. The purpose of this paper is to model friction as a result of mesoscopic pressure fluctuations. First, the simple random walk will be used to approximate the time and space scales below which ballistic motion begins and diffusive motion ends. Following that, pressure fluctuations and the associated time scale will be introduced to explain Brownian motion. As successful as the pressure fluctuation model is for many fluids, it will be shown to fail for fluids like glycerin that have viscosities a thousand times and more that of water.


title: On the site percolation threshold of circle packings and planar graphs

abstract: A circle packing is a collection of disks with disjoint interiors in the plane. It naturally defines a graph by tangency. It is shown that there exists $p>0$ such that the following holds for every circle packing: If each disk is retained with probability $p$ independently, then the probability that there is a path of retained disks connecting the origin to infinity is zero. The following conclusions are derived using results on circle packings of planar graphs: (i) Site percolation with parameter $p$ has no infinite connected component on recurrent simple plane triangulations, or on Benjamini--Schramm limits of finite simple planar graphs. (ii) Site percolation with parameter $1-p$ has an infinite connected component on transient simple plane triangulations with bounded degree. These results lend support to recent conjectures of Benjamini. Extensions to graphs formed from the packing of shapes other than disks, in the plane and in higher dimensions, are presented. Several conjectures and open questions are discussed.


title: Superfluid helium as the condensate of bound atoms pairs

abstract: Based upon the hypothesis of simultaneous participation of helium particles in atom(atom) and pair(pair) motions, the proof of Feynman and Bogoliubov formulas validity for elementary excitations spectrum in superfluid helium is carried out with the double mass of particle.


title: Evaluating the Progress of Deep Learning for Visual Relational Concepts

abstract: Convolutional Neural Networks (CNNs) have become the state of the art method for image classification in the last ten years. Despite the fact that they achieve superhuman classification accuracy on many popular datasets, they often perform much worse on more abstract image classification tasks. We will show that these difficult tasks are linked to relational concepts from cognitive psychology and that despite progress over the last few years, such relational reasoning tasks still remain difficult for current neural network architectures. We will review deep learning research that is linked to relational concept learning, even if it was not originally presented from this angle. Reviewing the current literature, we will argue that some form of attention will be an important component of future systems to solve relational tasks. In addition, we will point out the shortcomings of currently used datasets, and we will recommend steps to make future datasets more relevant for testing systems on relational reasoning.


title: The Case for Repeatable Analysis with Energy Economy Optimization Models

abstract: Energy economy optimization (EEO) models employ formal search techniques to explore the future decision space over several decades in order to deliver policy-relevant insights. EEO models are a critical tool for decision-makers who must make near-term decisions with long-term effects in the face of large future uncertainties. While the number of model-based analyses proliferates, insufficient attention is paid to transparency in model development and application. Given the complex, data-intensive nature of EEO models and the general lack of access to source code and data, many of the assumptions underlying model-based analysis are hidden from external observers. This paper discusses the simplifications and subjective judgments involved in the model building process, which cannot be fully articulated in journal papers, reports, or model documentation. In addition, we argue that for all practical purposes, EEO model-based insights cannot be validated through comparison to real world outcomes. As a result, modelers are left without credible metrics to assess a model's ability to deliver reliable insight. We assert that EEO models should be discoverable through interrogation of publicly available source code and data. In addition, third parties should be able to run a specific model instance in order to independently verify published results. Yet a review of twelve EEO models suggests that in most cases, replication of model results is currently impossible. We provide several recommendations to help develop and sustain a software framework for repeatable model analysis.


title: Prospects for large N gauge theories on the lattice

abstract: I will review recent progress on addressing large N gauge theories on the lattice. The focus will be put on the use of large N volume independence as an effective tool to compute non-perturbative dynamics at, otherwise unreachable, large number of colours. A selection of results will be presented and future prospects and challenges for the study of large N QCD and various extensions will be discussed.


title: A New Angle on Lattice Sieving for the Number Field Sieve

abstract: Lattice sieving in two or more dimensions has proven to be an indispensable practical aid in integer factorization and discrete log computations involving the number field sieve. The main contribution of this article is to show that a different method of lattice enumeration in three dimensions will provide a significant speedup. We use the successive minima and shortest vectors of the lattice instead of transition vectors to iterate through lattice points. We showcase the new method by a record computation in a 133-bit subgroup of $\mathbb{F}_{p^6}$, with $p^6$ having 423 bits. Our overall timing nearly $3$ times faster than the previous record of a 132-bit subgroup in a 422-bit field. The approach generalizes to dimensions 4 or more, overcoming a key obstruction to the implementation of the tower number field sieve.


title: Identifying trending coefficients with an ensemble Kalman filter

abstract: This paper extends the ensemble Kalman filter (EnKF) for inverse problems to identify trending model coefficients. This is done by repeatedly inflating the ensemble while maintaining the mean of the particles. As a benchmark serves a classic EnKF and a recursive least squares (RLS) on the example of identifying a force model in milling, which changes due to the progression of tool wear. For a proper comparison, the true values are simulated and augmented with white Gaussian noise. The results demonstrate the feasibility of the approach for dynamic identification while still achieving good accuracy in the static case. Further, the inflated EnKF shows a remarkably insensitivity on the starting set but a less smooth convergence compared to the classic EnKF.


title: Analytical Modeling of the Path-Loss for Reconfigurable Intelligent Surfaces -- Anomalous Mirror or Scatterer ?

abstract: Reconfigurable intelligent surfaces (RISs) are an emerging field of research in wireless communications. A fundamental component for analyzing and optimizing RIS-empowered wireless networks is the development of simple but sufficiently accurate models for the power scattered by an RIS. By leveraging the general scalar theory of diffraction and the Huygens-Fresnel principle, we introduce simple formulas for the electric field scattered by an RIS that is modeled as a sheet of electromagnetic material of negligible thickness. The proposed approach allows us to identify the conditions under which an RIS of finite size can or cannot be approximated as an anomalous mirror. Numerical results are illustrated to confirm the proposed approach.

categories: physics.class-ph cs.ET eess.SP physics.optics


title: The Voltage Regulation of Boost Converters Using Dual Heuristic Programming

abstract: In this paper, a dual heuristic programming controller is proposed to control a boost converter. Conventional controllers such as proportional integral derivative (PID) or proportional integral (PI) are designed based on the linearized small-signal model near the operating point. Therefore, the performance of the controller during start up, load change, or input voltage variation is not optimal since the system model changes by varying the operating point. The dual heuristic programming controller optimally controls the boost converter by following the approximate dynamic programming. The advantage of the DHP is that the neural network based characteristic of the proposed controller enables boost converters to easily cope with large disturbances. A DHP with a well trained critic and action networks can perform as an optimal controller for the boost converter. To compare the effectiveness of the traditional PI based and the DHP boost converter, the simulation results are provided.


title: On $\phi^3$ Theory Above Six Dimensions

abstract: We study the scalar $\phi^3$ theory above six dimensions. The beta function $\beta(g)=-\epsilon g-\frac{3}{4}g^3$ in $d=6-2\epsilon$ dimensions has a UV fixed point when $\epsilon<0$. Like the $O(N)$ vector models above four dimensions, such a fixed point observed perturbatively in fact corresponds to a pair of complex CFTs separated by a branch cut. Using both the numerical bootstrap method and Gliozzi's fusion rule truncation method, we argue that the fixed points of the $\phi^3$ theory above six dimensions exist.


title: Smart Resource Management for Data Streaming using an Online Bin-packing Strategy

abstract: Data stream processing frameworks provide reliable and efficient mechanisms for executing complex workflows over large datasets. A common challenge for the majority of currently available streaming frameworks is efficient utilization of resources. Most frameworks use static or semi-static settings for resource utilization that work well for established use cases but lead to marginal improvements for unseen scenarios. Another pressing issue is the efficient processing of large individual objects such as images and matrices typical for scientific datasets. HarmonicIO has proven to be a good solution for streams of relatively large individual objects, as demonstrated in a benchmark comparison with the Spark and Kafka streaming frameworks. We here present an extension of the HarmonicIO framework based on the online bin-packing algorithm, to allow for efficient utilization of resources. Based on a real world use case from large-scale microscopy pipelines, we compare results of the new system to Spark's auto-scaling mechanism.


title: Data integration and prediction models of photovoltaic production from Brazilian northeastern

abstract: All productive branches of society need an estimate to be able to control their expenses well. In the energy business, electric utilities use this information to control the power flow in the grid. For better energy production estimation of photovoltaic systems, it is necessary to join multiples geospatial and meteorological variables. This work proposes the creation of a satellite data integration platform, with production estimation models, base stations measurement and actual production capacity. This work presents statistical, probabilistic and artificial intelligence models that generate spatial and temporal production estimates that could improve production gains as well as facilitate the monitoring and supervision of new enterprises are presented.


title: Radial velocity constraints on the long-period transiting planet Kepler-1625 b with CARMENES

abstract: The star Kepler-1625 recently attracted considerable attention when an analysis of the stellar photometric time series from the Kepler mission was interpreted as showing evidence of a large exomoon around the transiting Jupiter-sized planet candidate Kepler-1625b. We aim to detect the radial velocity (RV) signal imposed by Kepler-1625b (and its putative moon) on the host star or, as the case may be, determine an upper limit on the mass of the transiting object. We took a total of 22 spectra of Kepler-1625 using CARMENES, 20 of which were useful. Observations were spread over a total of seven nights between October 2017 and October 2018, covering $125\%$ of one full orbit of Kepler-1625b. We used the automatic Spectral Radial Velocity Analyser (SERVAL) pipeline to deduce the stellar RVs and uncertainties. Then we fitted the RV curve model of a single planet on a Keplerian orbit to the observed RVs using a $\chi^2$ minimisation procedure. We derive upper limits on the mass of Kepler-1625b under the assumption of a single planet on a circular orbit. In this scenario, the $1\,\sigma$, $2\,\sigma$, and $3\,\sigma$ confidence upper limits for the mass of Kepler-1625b are $2.90\,M_{\rm J}$, $7.15\,M_{\rm J}$, and $11.60\,M_{\rm J}$, respectively. We present strong evidence for the planetary nature of Kepler-1625b, making it the 10th most long-period confirmed planet known today. Our data does not answer the question about a second, possibly more short-period planet that could be responsible for the observed transit timing variation of Kepler-1625b.


title: Uniform error bounds of a time-splitting spectral method for the long-time dynamics of the nonlinear Klein-Gordon equation with weak nonlinearity

abstract: We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein--Gordon equation (NKGE) with weak power-type nonlinearity and $O(1)$ initial data, while the nonlinearity strength is characterized by $\varepsilon^{p}$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$, for the long-time dynamics up to the time at $O(\varepsilon^{-\beta})$ with $0 \leq \beta \leq p$. In fact, when $0 < \varepsilon \ll 1$, the problem is equivalent to the long-time dynamics of NKGE with small initial data and $O(1)$ nonlinearity strength, while the amplitude of the initial data (and the solution) is at $O(\varepsilon)$. By reformulating the NKGE into a relativistic nonlinear Schr\"{o}dinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at $O(h^{m}+\varepsilon^{p-\beta}\tau^2)$ of the TSFP method with $h$ mesh size, $\tau$ time step and $m\ge2$ depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at $O(\varepsilon^{-\beta})$ and uniformly valid for $\varepsilon\in(0,1]$. Especially, the error bounds are uniformly at the second order rate for the large time step $\tau = O(\varepsilon^{-(p-\beta)/2})$ in the parameter regime $0\le\beta <p$. Numerical results are reported to confirm our error bounds in the long-time regime. Finally, the TSFP method and its error bounds are extended to a highly oscillatory complex NKGE which propagates waves with wavelength at $O(1)$ in space and $O(\varepsilon^{\beta})$ in time and wave velocity at $O(\varepsilon^{-\beta})$.


title: A geometric construction of representations of the Berezin-Toeplitz quantization

abstract: For a K\"ahler manifold $X$ equipped with a prequantum line bundle $L$, we give a geometric construction of a family of representations of the Berezin-Toeplitz deformation quantization algebra $(C^\infty(X)[[\hbar]],\star_{BT})$ parametrized by points $z_0 \in X$. The key idea is to use peak sections to suitably localize the Hilbert spaces $H^{0}\left(X,L^{\otimes m}\right)$ around $z_{0}$ in the large volume limit.


title: Is Your Quantum Program Bug-Free?

abstract: Quantum computers are becoming more mainstream. As more programmers are starting to look at writing quantum programs, they face an inevitable task of debugging their code. How should the programs for quantum computers be debugged? In this paper, we discuss existing debugging tactics, used in developing programs for classic computers, and show which ones can be readily adopted. We also highlight quantum-computer-specific debugging issues and list novel techniques that are needed to address these issues. The practitioners can readily apply some of these tactics to their process of writing quantum programs, while researchers can learn about opportunities for future work.


title: Efficient and robust chiral resolution by composite pulses

abstract: We introduce a method for detection of chiral molecules using sequences of three pulses driving a closed-loop three-state quantum system. The left- and right-handed enantiomers have identical optical properties (transition frequencies and transition dipole moments) with the only difference being the sign of one of the couplings. We identify twelve different sequences of resonant pulses for which chiral resolution with perfect contrast occurs. In all of them the first and third pulses are $\pi/2$-pulses and the middle pulse is a $\pi$-pulse. In addition, one of the three pulses must have a phase shift of $\pi/2$ with respect to the other two. The simplicity of the proposed chiral resolution technique allows for straightforward extensions to more efficient and more robust implementations by replacing the single $\pi/2$ and $\pi$-pulses by composite pulses. We present specific examples of chiral resolution by composite pulses which compensate errors in the pulse areas and the detuning of the driving fields.


title: A scale-dependent notion of effective dimension

abstract: We introduce a notion of "effective dimension" of a statistical model based on the number of cubes of size $1/\sqrt{n}$ needed to cover the model space when endowed with the Fisher Information Matrix as metric, $n$ being the number of observations. The number of observations fixes a natural scale or resolution. The effective dimension is then measured via the spectrum of the Fisher Information Matrix regularized using this natural scale.


title: Model category structures on multicomplexes

abstract: We present a family of model structures on the category of multicomplexes. There is a cofibrantly generated model structure in which the weak equivalences are the morphisms inducing an isomorphism at a fixed stage of an associated spectral sequence. Corresponding model structures are given for truncated versions of multicomplexes, interpolating between bicomplexes and multicomplexes. For a fixed stage of the spectral sequence, the model structures on all these categories are shown to be Quillen equivalent.


title: Cyclic abelian varieties over finite fields in ordinary isogeny classes

abstract: Given an abelian variety $A$ defined over a finite field $k$, we say that $A$ is "cyclic" if its group $A(k)$ of rational points is cyclic. In this paper we give a bijection between cyclic abelian varieties of an ordinary isogeny class $\mathcal{A}$ with Weil polynomial $f_{\mathcal{A}}$ and some classes of matrices with integer coefficients and having $f_{\mathcal{A}}$ as characteristic polynomial.


title: Algorithms for new types of fair stable matchings

abstract: We study the problem of finding "fair" stable matchings in the Stable Marriage problem with Incomplete lists (SMI). For an instance $I$ of SMI there may be many stable matchings, providing significantly different outcomes for the sets of men and women. We introduce two new notions of fairness in SMI. Firstly, a regret-equal stable matching minimises the difference in ranks of a worst-off man and a worst-off woman, among all stable matchings. Secondly, a min-regret sum stable matching minimises the sum of ranks of a worst-off man and a worst-off woman, among all stable matchings. We present two new efficient algorithms to find stable matchings of these types. Firstly, the Regret-Equal Degree Iteration Algorithm finds a regret-equal stable matching in $O(d_0 nm)$ time, where $d_0$ is the absolute difference in ranks between a worst-off man and a worst-off woman in the man-optimal stable matching, $n$ is the number of men or women, and $m$ is the total length of all preference lists. Secondly, the Min-Regret Sum Algorithm finds a min-regret sum stable matching in $O(d_s m)$ time, where $d_s$ is the difference in the ranks between a worst-off man in each of the woman-optimal and man-optimal stable matchings. Experiments to compare several types of fair optimal stable matchings were conducted and show that the Regret-Equal Degree Iteration Algorithm produces matchings that are competitive with respect to other fairness objectives. On the other hand, existing types of "fair" stable matchings did not provide as close an approximation to regret-equal stable matchings.


title: Compact recurrent neural networks for acoustic event detection on low-energy low-complexity platforms

abstract: Outdoor acoustic events detection is an exciting research field but challenged by the need for complex algorithms and deep learning techniques, typically requiring many computational, memory, and energy resources. This challenge discourages IoT implementation, where an efficient use of resources is required. However, current embedded technologies and microcontrollers have increased their capabilities without penalizing energy efficiency. This paper addresses the application of sound event detection at the edge, by optimizing deep learning techniques on resource-constrained embedded platforms for the IoT. The contribution is two-fold: firstly, a two-stage student-teacher approach is presented to make state-of-the-art neural networks for sound event detection fit on current microcontrollers; secondly, we test our approach on an ARM Cortex M4, particularly focusing on issues related to 8-bits quantization. Our embedded implementation can achieve 68% accuracy in recognition on Urbansound8k, not far from state-of-the-art performance, with an inference time of 125 ms for each second of the audio stream, and power consumption of 5.5 mW in just 34.3 kB of RAM.


title: On the behavior of extreme $d$-dimensional spatial quantiles under minimal assumptions

abstract: "Spatial" or "geometric" quantiles are the only multivariate quantiles coping with both high-dimensional data and functional data, also in the framework of multiple-output quantile regression. This work studies spatial quantiles in the finite-dimensional case, where the spatial quantile $\mu_{\alpha,u}(P)$ of the distribution $P$ taking values in $\mathbb{R}^d $ is a point in $\mathbb{R}^d$ indexed by an order $\alpha\in[0,1)$ and a direction $u$ in the unit sphere $\mathcal{S}^{d-1}$ of $\mathbb{R}^d$ --- or equivalently by a vector $\alpha u$ in the open unit ball of $\mathbb{R}^d$. Recently, Girard and Stupfler (2017) proved that (i) the extreme quantiles $\mu_{\alpha,u}(P)$ obtained as $\alpha\to 1$ exit all compact sets of $\mathbb{R}^d$ and that (ii) they do so in a direction converging to $u$. These results help understanding the nature of these quantiles: the first result is particularly striking as it holds even if $P$ has a bounded support, whereas the second one clarifies the delicate dependence of spatial quantiles on $u$. However, they were established under assumptions imposing that $P$ is non-atomic, so that it is unclear whether they hold for empirical probability measures. We improve on this by proving these results under much milder conditions, allowing for the sample case. This prevents using gradient condition arguments, which makes the proofs very challenging. We also weaken the well-known sufficient condition for uniqueness of finite-dimensional spatial quantiles.


title: Asteroseismology of luminous red giants with Kepler I: Long Period Variables with radial and non-radial modes

abstract: While long period variables (LPVs) have been extensively investigated, especially with MACHO and OGLE data for the Magellanic Clouds, there still exist open questions in their pulsations regarding the excitation mechanisms, radial order and angular degree assignment. Here, we perform asteroseismic analyses on LPVs observed by the 4-year Kepler mission. Using a cross-correlation method, we detect unambiguous pulsation ridges associated with radial fundamental modes ($n=1$) and overtones ($n\geqslant2$), where the radial order assignment is made by using theoretical frequencies and observed frequencies. Our results confirm that the amplitude variability seen in semiregulars is consistent with oscillations being solar-like. We identify that the dipole modes, $l=1$, are dominant in the radial orders of $3\leq n \leq6$, and that quadrupole modes, $l=2$, are dominant in the first overtone $n=2$. A test of seismic scaling relations using Gaia DR2 parallaxes reveals the possibility that the relations break down when $\nu_{\rm max}$ $\lesssim$ 3 $\mu$Hz (R $\gtrsim$ 40 R$_{\odot}$, or log $\rm L/L_{\odot}$ $\gtrsim$ 2.6). Our homogeneous measurements of pulsation amplitude and period for 3213 LPVs will be very valuable for probing effects of pulsation on mass loss, in particular in those stars with periods around 60 days, which has been argued as a threshold of substantial pulsation-triggered mass loss.


title: 2MASS J15491331-3539118: a new low-mass wide companion of the GQ Lup system

abstract: Substellar companions at wide separation around stars hosting planets or brown dwarfs (BDs) yet close enough for their formation in the circumstellar disc are of special interest. In this letter we report the discovery of a wide (projected separation $\sim$16.0arcsec, or 2400 AU, and position angle 114.61$^\circ$) companion of the GQ Lup A-B system, most likely gravitationally bound to it. A VLT/X-Shooter spectrum shows that this star, 2MASS J15491331-3539118, is a bonafide low-mass ($\sim$0.15 M$_\odot$) young stellar object (YSO) with stellar and accretion/ejection properties typical of Lupus YSOs of similar mass, and with kinematics consistent with that of the GQ Lup A-B system. A possible scenario for the formation of the triple system is that GQ Lup A and 2MASS J15491331-3539118 formed by fragmentation of a turbulent core in the Lup I filament, while GQ Lup B, the BD companion of GQ Lup A at 0.7arcsec, formed in situ by the fragmentation of the circumprimary disc. The recent discoveries that stars form along cloud filaments would favour the scenario of turbulent fragmentation for the formation of GQ Lup A and 2MASS J15491331-3539118.


title: An Analysis of the Shapes of Interstellar Extinction Curves. VIII. The Optical Extinction Structure

abstract: New HST/STIS optical spectra were obtained for a sample of early type stars with existing IUE UV spectra. These data were used to construct optical extinction curves whose general properties are discussed elsewhere. In this paper, we identify extinction features in the curves that are wider than diffuse interstellar bands (DIBs) but narrower than the well known broad band variability. This intermediate scale structure, or ISS, contains distinct features whose peaks can contribute a few percent to 20% of the total extinction. Most of the ISS variation can be captured by three principal components. We model the ISS with three Drude profiles and show that their strengths and widths vary from one sight line to another, but their central positions are stable, near 4370, 4870 and 6300A. The Very Broad Structure, VBS, in optical curves appears to be a minimum between the 4870 and 6300A absorption peaks. We find relations among the fit parameters and provide a physical interpretation of them in terms of a simplistic grain model. Finally, we note that the strengths of the 4370 and 4870A features are correlated to the strength of the 2175A UV bump, but that the 6300A feature is not, and that none of the ISS features are related to R(V). However, we verify that the broad band curvature of the continuous optical extinction is strongly related to R(V).


title: Provably Secure Isolation for Interruptible Enclaved Execution on Small Microprocessors: Extended Version

abstract: Computer systems often provide hardware support for isolation mechanisms like privilege levels, virtual memory, or enclaved execution. Over the past years, several successful software-based side-channel attacks have been developed that break, or at least significantly weaken the isolation that these mechanisms offer. Extending a processor with new architectural or micro-architectural features, brings a risk of introducing new such side-channel attacks. This paper studies the problem of extending a processor with new features without weakening the security of the isolation mechanisms that the processor offers. We propose to use full abstraction as a formal criterion for the security of a processor extension, and we instantiate that criterion to the concrete case of extending a microprocessor that supports enclaved execution with secure interruptibility of these enclaves. This is a very relevant instantiation as several recent papers have shown that interruptibility of enclaves leads to a variety of software-based side-channel attacks. We propose a design for interruptible enclaves, and prove that it satisfies our security criterion. We also implement the design on an open-source enclave-enabled microprocessor, and evaluate the cost of our design in terms of performance and hardware size.


title: Extremal Area of Polygons, sliding along a Circle

abstract: We determine all critical configurations for the Area function on polygons with vertices on a circle or an ellipse. For isolated critical points we compute their Morse index, resp index of the gradient vector field. We relate the computation at an isolated degenerate point to an eigenvalue question about combinations. In the even dimensional case non-isolated singularities occur as `zigzag trains'.


title: Unsupervised Anomaly Detection for X-Ray Images

abstract: Obtaining labels for medical (image) data requires scarce and expensive experts. Moreover, due to ambiguous symptoms, single images rarely suffice to correctly diagnose a medical condition. Instead, it often requires to take additional background information such as the patient's medical history or test results into account. Hence, instead of focusing on uninterpretable black-box systems delivering an uncertain final diagnosis in an end-to-end-fashion, we investigate how unsupervised methods trained on images without anomalies can be used to assist doctors in evaluating X-ray images of hands. Our method increases the efficiency of making a diagnosis and reduces the risk of missing important regions. Therefore, we adopt state-of-the-art approaches for unsupervised learning to detect anomalies and show how the outputs of these methods can be explained. To reduce the effect of noise, which often can be mistaken for an anomaly, we introduce a powerful preprocessing pipeline. We provide an extensive evaluation of different approaches and demonstrate empirically that even without labels it is possible to achieve satisfying results on a real-world dataset of X-ray images of hands. We also evaluate the importance of preprocessing and one of our main findings is that without it, most of our approaches perform not better than random. To foster reproducibility and accelerate research we make our code publicly available at https://github.com/Valentyn1997/xray


title: Helicity proxies from linear polarisation of solar active regions

abstract: The alpha effect is believed to play a key role in the generation of the solar magnetic field. A fundamental test for its significance in the solar dynamo is to look for magnetic helicity of opposite signs in the two hemispheres, and at small and large scales. However, measuring magnetic helicity is compromised by the inability to fully infer the magnetic field vector from observations of solar spectra, caused by what is known as the "pi ambiguity" of spectropolarimetric observations. We decompose linear polarisation into parity-even and parity-odd E and B polarisations, which are not affected by the "pi ambiguity". Furthermore, we study whether the correlations of spatial Fourier spectra of B and parity-even quantities such as E or temperature T are a robust proxy for magnetic helicity of solar magnetic fields. We analyse polarisation measurements of active regions observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics observatory. Theory predicts the magnetic helicity of active regions to have, statistically, opposite signs in the two hemispheres. We then compute the parity-odd E B and T B correlations, and test for systematic preference of their sign based on the hemisphere of the active regions. We find that: (i) E B and T B correlations are a reliable proxy for magnetic helicity, when computed from linear polarisation measurements away from spectral line cores, and (ii) E polarisation reverses its sign close to the line core. Our analysis reveals Faraday rotation to not have a significant influence on the computed parity-odd correlations. The EB decomposition of linear polarisation appears to be a good proxy for magnetic helicity independent of the "pi ambiguity". This allows us to routinely infer magnetic helicity directly from polarisation measurements.


title: Q = 0 order in quantum kagome Heisenberg antiferromagnet

abstract: We have studied the nearest neighbor Heisenberg model with added Dzyaloshinskii-Moriya interaction using Schwinger boson mean field theory considering in-plane component as well as out-of-plane component. Motivated by experimental result of Vesigniete that the ground state is in Q = 0 long range order state, we first looked at the classical ground state of the model and consider the mean field ansatz which mimics the classical ground state in the large S limit. We have obtained the ground state phase diagram of this model and calculated properties of different phases. We have also studied the above model numerically using exact diagonalization up to a system size N = 30. We have compared the obtained results from these two approaches. Our results are in agreement with the experimental result of the Vesigniete.


title: A new hybrid radiative transfer method for massive star formation

abstract: Frequency-dependent/hybrid approaches for stellar irradiation are of primary importance in numerical simulations of massive star formation. We seek to compare outflow and accretion mechanisms in star formation simulations. We investigate the accuracy of a hybrid radiative transfer method using the gray M1 closure relation for proto-stellar irradiation and gray flux-limited diffusion (FLD) for photons emitted everywhere else. We have coupled the FLD module of the adaptive-mesh refinement code Ramses with Ramses-RT, which is based on the M1 closure relation. Our hybrid (M1+FLD) method takes an average opacity at the stellar temperature for the M1 module, instead of the local environmental radiation field. We have tested this approach in radiative transfer tests of disks irradiated by a star for three levels of optical thickness and compared the temperature structure with RADMC-3D and MCFOST. We applied it to a radiation-hydrodynamical simulation of massive star formation. Our tests validate our hybrid approach for determining the temperature structure of an irradiated disk in the optically-thin and moderately optically-thick regimes and the most optically-thick test shows the limitation of our approach. The optically-thick setups highlight the ability of the hybrid method to partially capture the self-shielding in the disk while the FLD cannot. The radiative acceleration is 100 times greater with the hybrid method. It consistently leads to about +50% more extended and wider-angle radiative outflows in the massive star formation simulation. We obtain a $17.6 M_\odot$ at $t{\simeq}0.7 \tau_\mathrm{ff}$, while the accretion phase is ongoing. Finally, despite the use of refinement to resolve the radiative cavities, no Rayleigh-Taylor instability appears in our simulations, and we justify their absence by physical arguments based on the entropy gradient. (abridged)


title: Model-independent constraints on Type Ia supernova light-curve hyper-parameters and reconstructions of the expansion history of the Universe

abstract: We reconstruct the expansion history of the Universe using type Ia supernovae (SN Ia) in a manner independent of any cosmological model assumptions. To do so, we implement a non-parametric iterative smoothing method on the Joint Light-curve Analysis (JLA) data while exploring the SN Ia light-curve hyper-parameter space by Monte Carlo Markov Chain sampling. We test to see how the posteriors of these hyper-parameters depend on cosmology, whether using different dark energy models or reconstructions shift these posteriors. Our constraints on the SN Ia light-curve hyper-parameters from our model-independent analysis are very consistent with the constraints from using different parametrizations of the equation of state of dark energy, namely the flat $\Lambda\mathrm{CDM}$ cosmology, the Chevallier-Polarski-Linder (CPL) model, and the Phenomenologically Emergent Dark Energy (PEDE) model. This implies that the distance moduli constructed from the JLA data are mostly independent of the cosmological models. We also studied that the possibility the light-curve parameters evolve with redshift and our results show consistency with no evolution. The reconstructed expansion history of the Universe and dark energy properties also seem to be in good agreement with the expectations of the standard $\Lambda\mathrm{CDM}$ model. However, our results also indicate that the data still allow for considerable flexibility in the expansion history of the Universe.


title: Cross-Layer Scheduling and Beamforming in Smart-Grid Powered Cellular Networks With Heterogeneous Energy Coordination

abstract: User scheduling, beamforming and energy coordination are investigated in smart-grid powered cellular networks (SGPCNs), where the base stations are powered by a smart grid and natural renewable energy sources. Heterogeneous energy coordination is considered in SGPCNs, namely energy merchandizing with the smart grid and energy exchanging among the base stations. A long-term grid-energy expenditure minimization problem with proportional-rate constraints is formulated for SGPCNs. Since user scheduling is coupled with the beamforming vectors, the formulated problem is challenging to handle via standard convex optimization methods. In practice, the beamforming vectors need to be updated over each slot according to the channel variations. User scheduling needs to be updated over several slots (frame) since the frequent scheduling of user equipment can cause reliability issues. Therefore, the Lyapunov optimization method is used to decouple the problem. A practical two-scale algorithm is proposed to schedule users at each frame, and obtain the beamforming vectors and amount of exchanged natural renewable energy at each slot. We prove that the proposed two-scale algorithm can asymptotically achieve the optimal solutions via tuning a control parameter. Numerical results verify the performance of the proposed two-scale algorithm.


title: Observation of a potential future sensitivity limitation from ground motion at LIGO Hanford

abstract: A first detection of terrestrial gravity noise in gravitational-wave detectors is a formidable challenge. With the help of environmental sensors, it can in principle be achieved before the noise becomes dominant by estimating correlations between environmental sensors and the detector. The main complication is to disentangle different coupling mechanisms between the environment and the detector. In this paper, we analyze the relations between physical couplings and correlations that involve ground motion and LIGO strain data h(t) recorded during its second science run in 2016 and 2017. We find that all noise correlated with ground motion was more than an order of magnitude lower than dominant low-frequency instrument noise, and the dominant coupling over part of the spectrum between ground and h(t) was residual coupling through the seismic-isolation system. We also present the most accurate gravitational coupling model so far based on a detailed analysis of data from a seismic array. Despite our best efforts, we were not able to unambiguously identify gravitational coupling in the data, but our improved models confirm previous predictions that gravitational coupling might already dominate linear ground-to-h(t) coupling over parts of the low-frequency, gravitational-wave observation band.


title: Toeplitz kernels and the backward shift

abstract: In this paper we study the kernels of Toeplitz operators on both the scalar and the vector-valued Hardy space for $ 1 < p < \infty $. We show existence of a minimal kernel of any element of the vector-valued Hardy space and we determine a symbol for the corresponding Toeplitz operator. In the scalar case we give an explicit description of a maximal function for a given Toeplitz kernel. In the vectorial case we show not all Toeplitz kernels have a maximal function and in the case of $p=2$ we find the exact conditions for when a Toeplitz kernel has a maximal function. For both the scalar and vector-valued Hardy space we study the minimal Toeplitz kernel containing multiple elements of the Hardy space, we also find an equivalent condition for a function in the Smirnov class to be cyclic for the backwards shift.


title: On the RLWE/PLWE equivalence for cyclotomic number fields

abstract: We study the equivalence between the Ring Learning With Errors and Polynomial Learning With Errors problems for cyclotomic number fields,namely: we prove that both problems are equivalent via a polynomial noise increase as long as the number of distinct primes dividing the conductor is kept constant. We refine our bound in the case where the conductor is divisible by at most three primes and we give an asymptotic subexponential formula for the condition number of the attached Vandermonde matrix valid for arbitrary degree.


title: Discriminating between and within (semi)continuous classes of both Tweedie and geometric Tweedie models

abstract: In both Tweedie and geometric Tweedie models, the common power parameter $p\notin(0,1)$ works as an automatic distribution selection. It mainly separates two subclasses of semicontinuous ($1<p<2$) and positive continuous ($p\geq 2$) distributions. Our paper centers around exploring diagnostic tools based on the maximum likelihood ratio test and minimum Kolmogorov-Smirnov distance methods in order to discriminate very close distributions within each subclass of these two models according to values of $p$. Grounded on the unique equality of variation indices, we also discriminate the gamma and geometric gamma distributions with $p=2$ in Tweedie and geometric Tweedie families, respectively. Probabilities of correct selection for several combinations of dispersion parameters, means and sample sizes are examined by simulations. We thus perform a numerical comparison study to assess the discrimination procedures in these subclasses of two families. Finally, semicontinuous ($1<p\leq 2$) distributions in the broad sense are significantly more distinguishable than the over-varied continuous ($p>2$) ones; and two datasets for illustration purposes are investigated.


title: Privacy-Preserving Gaussian Process Regression -- A Modular Approach to the Application of Homomorphic Encryption

abstract: Much of machine learning relies on the use of large amounts of data to train models to make predictions. When this data comes from multiple sources, for example when evaluation of data against a machine learning model is offered as a service, there can be privacy issues and legal concerns over the sharing of data. Fully homomorphic encryption (FHE) allows data to be computed on whilst encrypted, which can provide a solution to the problem of data privacy. However, FHE is both slow and restrictive, so existing algorithms must be manipulated to make them work efficiently under the FHE paradigm. Some commonly used machine learning algorithms, such as Gaussian process regression, are poorly suited to FHE and cannot be manipulated to work both efficiently and accurately. In this paper, we show that a modular approach, which applies FHE to only the sensitive steps of a workflow that need protection, allows one party to make predictions on their data using a Gaussian process regression model built from another party's data, without either party gaining access to the other's data, in a way which is both accurate and efficient. This construction is, to our knowledge, the first example of an effectively encrypted Gaussian process.


title: On the structure of the boundary layer in a Beklemishev diamagnetic bubble

abstract: The article provides a kinetic description of the plasma equilibrium in the Beklemishev diamagnetic trap, where the traditional approach based on the theory of magnetic drifts is not applicable, since the ions move in a substantially non-circular orbit, the diameter of which is approximately equal to the diameter of the diamagnetic bubble. The ion distribution function was found in the collisionless approximation, neglecting the diamagnetic electron current. The radial profile of the magnetic field, the plasma density, the current density, and the components of the pressure tensor are calculated. It was found that the width of the boundary layer in the diamagnetic bubble varies from 6 to 8 Larmor radii calculated by a vacuum magnetic field. An adiabatic invariant is calculated that replaces the magnetic moment, which is not conserved in the diamagnetic bubble. The criterion of absolute confinement is formulated and the plasma equilibrium is found for the case when the adiabatic invariant is not conserved and only ions whose velocity satisfies the criterion of absolute confinement are trapped.


title: Chromospheric activity in bright contact binary stars

abstract: Studying chromospheric activity of contact binaries is an important way of revealing the magnetic activity processes of these systems. An efficient but somewhat neglected method for that is to follow the changes of the H$\alpha$ line profiles via optical spectroscopy. Our goal was to perform a comprehensive analysis based on the optical spectral signs of chromospheric activity on the largest sample of contact binaries to date. We collected echelle spectra on 12 bright contact binaries and derived new radial velocity curves from our observations. For quantifying the apparent chromospheric activity levels of the systems, we subtracted self-constructed synthetic spectra from the observed ones and measured the equivalent widths of the residual H$\alpha$-profiles at each observed epoch. Our well-sampled data set allowed us to study the short-term variations of chromospheric activity levels as well as to search for correlations between them and some basic physical parameters of the systems. Fitting the radial velocity curves, we re-determined the mass ratios and systemic velocities of all observed objects. We found that chromospheric activity levels of the studied systems show various changes during the orbital revolution: we see either flat, or one-peaked, or two-peaked distributions of equivalent width vs. the orbital phase. The first case means that the activity level is probably constant, while the latter two cases suggest the presence of one or two active longitudes at the stellar surfaces. Our correlation diagrams show that mean chromospheric activity levels may be in connection with orbital periods, B$-$V color indices, inverse Rossby numbers, and temperature differences of the components. At the same time, no clear trend is visible with respect to mass ratios, inclinations and fill-out factors of the systems. A- and W-type contact binaries in our sample show similar distributions.


title: Explicit Solutions to Fractional Stefan-like problems for Caputo and Riemann-Liouville Derivatives

abstract: Two fractional two-phase Stefan-like problems are considered by using Riemann-Liouville and Caputo derivatives of order $\alpha \in (0, 1)$ verifying that they coincide with the same classical Stefan problem at the limit case when $\alpha=1$. For both problems, explicit solutions in terms of the Wright functions are presented. Even though the similarity of the two solutions, a proof that they are different is also given. The convergence when $\alpha \nearrow 1$ of the one and the other solutions to the same classical solution is given. Numerical examples for the dimensionless version of the problem are also presented and analyzed.


title: Accelerating the evaluation of inspiral-merger-ringdown waveforms with adapted grids

abstract: This paper presents an algorithm to accelerate the evaluation of inspiral-merger-ringdown waveform models for gravitational wave data analysis. While the idea can also be applied in the time domain, here we focus on the frequency domain, which is most typically used to reduced computational cost in gravitational wave data analysis. Our work extends the idea of multibanding, which has been developed to accelerate frequency domain waveforms, to include the merger and ringdown and spherical harmonics beyond the dominant quadrupole spherical harmonic. The original method is based on a heuristic algorithm based on the inspiral to de-refine the equi-spaced frequency grid used for data analysis where a coarser grid is sufficient for accurate evaluation of a waveform model. Here we use a different criterion, based on the local interpolation error, which is more flexible and can easily be adapted to general waveforms, if their phenomenology is understood. We discuss our implementation in the LIGO Algorithm Library for the PhenomXHM frequency domain model, and report the acceleration in different parts of the parameter space of compact binary systems.


title: ScreenTrack: Using a Visual History of a Computer Screen to Retrieve Documents and Web Pages

abstract: Computers are used for various purposes, so frequent context switching is inevitable. In this setting, retrieving the documents, files, and web pages that have been used for a task can be a challenge. While modern applications provide a history of recent documents for users to resume work, this is not sufficient to retrieve all the digital resources relevant to a given primary document. The histories currently available do not take into account the complex dependencies among resources across applications. To address this problem, we tested the idea of using a visual history of a computer screen to retrieve digital resources within a few days of their use through the development of ScreenTrack. ScreenTrack is software that captures screenshots of a computer at regular intervals. It then generates a time-lapse video from the captured screenshots and lets users retrieve a recently opened document or web page from a screenshot after recognizing the resource by its appearance. A controlled user study found that participants were able to retrieve requested information more quickly with ScreenTrack than under the baseline condition with existing tools. A follow-up study showed that the participants used ScreenTrack to retrieve previously used resources and to recover the context for task resumption.


title: Focusing on Bandwidth: Achromatic Metalens Limits

abstract: Metalenses have shown great promise in their ability to function as ultracompact optical systems for focusing and imaging. Remarkably, several designs have been recently demonstrated that operate over a large range of frequencies with minimized chromatic aberrations, potentially paving the way for ultrathin achromatic optics. Here, we derive fundamental bandwidth limits that apply to broadband optical metalenses, regardless of their implementation. Specifically, we discuss how the product between achievable time delay and bandwidth is limited in any time-invariant system, and we apply well-established bounds on this product to a general focusing system. We then show that all metalenses designed thus far obey the appropriate bandwidth limit. The derived physical bounds provide a useful metric to compare and assess the performance of different devices, and offer fundamental insight into how to design better broadband metalenses.


title: On Learning Vehicle Detection in Satellite Video

abstract: Vehicle detection in aerial and satellite images is still challenging due to their tiny appearance in pixels compared to the overall size of remote sensing imagery. Classical methods of object detection very often fail in this scenario due to violation of implicit assumptions made such as rich texture, small to moderate ratios between image size and object size. Satellite video is a very new modality which introduces temporal consistency as inductive bias. Approaches for vehicle detection in satellite video use either background subtraction, frame differencing or subspace methods showing moderate performance (0.26 - 0.82 $F_1$ score). This work proposes to apply recent work on deep learning for wide-area motion imagery (WAMI) on satellite video. We show in a first approach comparable results (0.84 $F_1$) on Planet's SkySat-1 LasVegas video with room for further improvement.


title: q-difference equations associated with the Rubin's q-difference operator $\partial_{q}$

abstract: The aim of this paper is to prove the existence and uniqueness of solutions of the following $q$- Cauchy problem of second order linear $q$-difference problem associated with the Rubin's $q$- difference operator $\partial_q$ in a neighborhood of zero \begin{equation} \left\{ \begin{array}{cc} q\,a_0(x)\, \partial_q^2y(qx)\, +\ ,a_1(x)\,\partial_qy(x)\, + \,a_2(x)y(x) &\; = \;b(x),\quad \hbox{if $y$ is odd;}\\ q\,a_0(x) \partial_q^2y(qx)\, + \,q\,a_1(x)\partial_qy(qx)\, + \,a_2(x)y(x)&\; = \;b(x),\quad \hbox{if $y$ is even,} \end{array} \right. \end{equation} with the initial conditions \begin{equation} \partial_{q}^{i-1}y(0)= b_{i};\quad b_{i} \in{\mathbb{C}},\; i=1,2 \end{equation} where $a_i$, $i=0,1,2$, and $b$ are defined, continuous at zero and bounded on an interval $I$ containing zero such that $a_0(x)\neq 0$ for all $x\in I$. Then, as application of the main results, we study the second order homogenous linear $q$- difference equations as well as the $q$-Wronskian associated with the Rubin's $q$-difference operator $\partial_q$. Finally, we construct a fundamental set of solutions for the second order linear homogeneous $q$-difference equations in the cases when the coefficients are constants and $a_1(x)=0$ for all $x\in I$.


title: RPCA-Based High Resolution Through-the-Wall Human Motion Detection and Classification

abstract: Radar based assisted living has received great amount of research interest in recent years. By employing the micro-Doppler features of indoor human motions, accurate recognition and classification of different types of movements become possible. Whereas, most of the existing works are focused only on free space detection, the literature on detection and recognition of human motions in through-the-wall scenarios is still in its infancy. As can be anticipated, the wall media and indoor static non-human targets would cause clutters and significantly corrupt the motion information of human subjects behind wall. However, no relevant work is reported to effectively handle this problem. In the present work, we aim to fill the gap and propose to use a low center-frequency ultra-wideband (UWB) radar system to probe the behind wall scene. Then, a Robust Principal Component Analysis (RPCA) based subspace decomposition technique, as its first reported implementation, is employed not only to remove the stationary clutters in raw range slow-time map but also to mitigate the multipath effects in the time-frequency map. Onsite experiments of detecting human motions behind a single layer of concrete wall is carried out to investigate the performance of the technique. Lastly, a two dimensional (2D)-PCA algorithm-based motion classification is provided to further verify the effectiveness of the proposed technique. Classification result shows that an enhanced recognition capability can be achieved using the proposed technique in detection and classification of indoor human motions.


title: Plastic Instabilities in Charged Granular Systems: Competition between Elasticity and Electrostatics

abstract: Electrostatic theory preserves charges, but allows dipolar excitations. Elasticity theory preserves dipoles, but allows quadrupolar (Eshelby like) plastic events. Charged amorphous granular systems are interesting in their own right; here we focus on their plastic instabilities and examine their mechanical response to external strain and to external electric field, to expose the competition between elasticity and electrostatics. In this paper a generic model is offered, its mechanical instabilities are examined and a theoretical analysis is presented. Plastic instabilities are discussed as saddle-node bifurcations that can be fully understood in terms of eigenvalues and eigenfunctions of the relevant Hessian matrix. This system exhibits moduli that describe how electric polarization and stress are influenced by strain and electric field. Theoretical expression for these moduli are offered and compared to the measurements in numerical simulations.


title: Optimal Augmented-Channel Puncturing for Low-Complexity Soft-Output MIMO Detectors

abstract: We propose a computationally-efficient soft-output detector for multiple-input multiple-output channels based on augmented channel puncturing in order to reduce tree processing complexity. The proposed detector, dubbed augmented WL detector (AWLD), employs a punctured channel with a special structure derived by triangulizing the original channel in augmented form, followed by Gaussian elimination. We prove that these punctured channels are optimal in maximizing the lower-bound on the achievable information rate (AIR) based on a newly proposed mismatched detection model. We show that the AWLD decomposes into a minimum mean-square error (MMSE) prefilter and channel-gain compensation stages, followed by a regular unaugmented WL detector (WLD). It attains the same performance as the existing AIR partial marginalization (AIR-PM) detector, but with much simpler processing.