title: Robust Revenue Maximization Under Minimal Statistical Information

abstract: We study the problem of multi-dimensional revenue maximization when selling $m$ items to a buyer that has additive valuations for them, drawn from a (possibly correlated) prior distribution. Unlike traditional Bayesian auction design, we assume that the seller has a very restricted knowledge of this prior: they only know the mean $\mu_j$ and an upper bound $\sigma_j$ on the standard deviation of each item's marginal distribution. Our goal is to design mechanisms that achieve good revenue against an ideal optimal auction that has full knowledge of the distribution in advance. Informally, our main contribution is a tight quantification of the interplay between the dispersity of the priors and the aforementioned robust approximation ratio. Furthermore, this can be achieved by very simple selling mechanisms. More precisely, we show that selling the items via separate price lotteries achieves an $O(\log r)$ approximation ratio where $r=\max_j(\sigma_j/\mu_j)$ is the maximum coefficient of variation across the items. To prove the result, we leverage a price lottery for the single-item case. If forced to restrict ourselves to deterministic mechanisms, this guarantee degrades to $O(r^2)$. Assuming independence of the item valuations, these ratios can be further improved by pricing the full bundle. For the case of identical means and variances, in particular, we get a guarantee of $O(\log(r/m))$ which converges to optimality as the number of items grows large. We demonstrate the optimality of the above mechanisms by providing matching lower bounds. Our tight analysis for the single-item deterministic case resolves an open gap from the work of Azar and Micali [ITCS'13]. As a by-product, we also show how one can directly use our upper bounds to improve and extend previous results related to the parametric auctions of Azar et al. [SODA'13].


title: Quantum key distribution based on the quantum eraser

abstract: Quantum information and quantum foundations are becoming popular topics for advanced undergraduate courses. Many of the fundamental concepts and applications in these two fields, such as delayed choice experiments and quantum encryption, are comprehensible to undergraduates with basic knowledge of quantum mechanics. In this paper, we show that the quantum eraser, usually used to study the duality between wave and particle properties, can also serve as a generic platform for quantum key distribution. We present a pedagogical example of an algorithm to securely share random keys using the quantum eraser platform and propose its implementation with quantum circuits.


title: Void region segmentation in ball grid array using u-net approach and synthetic data

abstract: The quality inspection of solder balls by detecting and measuring the void is important to improve the board yield issues in electronic circuits. In general, the inspection is carried out manually, based on 2D or 3D X-ray images. For high quality inspection, it is difficult to detect and measure voids accurately with high repeatability through the manual inspection and the process is time consuming. In need of high quality and fast inspection, various approaches were proposed, but, due to the various challenges like vias, reflections from the plating or vias, inconsistent lighting, noise and void-like artifacts makes these approaches difficult to work in all these challenging conditions. In recent times, deep learning approaches are providing the outstanding accuracy in various computer vision tasks. Considering the need of high quality and fast inspection, in this paper, we applied U-Net to segment the void regions in soldering balls. As it is difficult to get the annotated dataset covering all the variations of void, we proposed an approach to generated the synthetic dataset. The proposed approach is able to segment the voids and can be easily scaled to various electronic products.


title: Characterizing Inter-Layer Functional Mappings of Deep Learning Models

abstract: Deep learning architectures have demonstrated state-of-the-art performance for object classification and have become ubiquitous in commercial products. These methods are often applied without understanding (a) the difficulty of a classification task given the input data, and (b) how a specific deep learning architecture transforms that data. To answer (a) and (b), we illustrate the utility of a multivariate nonparametric estimator of class separation, the Henze-Penrose (HP) statistic, in the original as well as layer-induced representations. Given an $N$-class problem, our contribution defines the $C(N,2)$ combinations of HP statistics as a sample from a distribution of class-pair separations. This allows us to characterize the distributional change to class separation induced at each layer of the model. Fisher permutation tests are used to detect statistically significant changes within a model. By comparing the HP statistic distributions between layers, one can statistically characterize: layer adaptation during training, the contribution of each layer to the classification task, and the presence or absence of consistency between training and validation data. This is demonstrated for a simple deep neural network using CIFAR10 with random-labels, CIFAR10, and MNIST datasets.


title: Analyzing Phonetic and Graphemic Representations in End-to-End Automatic Speech Recognition

abstract: End-to-end neural network systems for automatic speech recognition (ASR) are trained from acoustic features to text transcriptions. In contrast to modular ASR systems, which contain separately-trained components for acoustic modeling, pronunciation lexicon, and language modeling, the end-to-end paradigm is both conceptually simpler and has the potential benefit of training the entire system on the end task. However, such neural network models are more opaque: it is not clear how to interpret the role of different parts of the network and what information it learns during training. In this paper, we analyze the learned internal representations in an end-to-end ASR model. We evaluate the representation quality in terms of several classification tasks, comparing phonemes and graphemes, as well as different articulatory features. We study two languages (English and Arabic) and three datasets, finding remarkable consistency in how different properties are represented in different layers of the deep neural network.


title: Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE

abstract: The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.


title: $\pi$-type Fermions and $\pi$-type KP hierarchy

abstract: In this paper, we firstly construct $\pi$-type Fermions. According to these, we define $\pi$-type Boson-Fermion correspondence which is a generalization of the classical Boson-Fermion correspondence. We can obtain $\pi$-type symmetric functions $S_\lambda^\pi$ from the $\pi$-type Boson-Fermion correspondence, analogously to the way we get the Schur functions $S_\lambda$ from the classical Boson-Fermion correspondence (which is the same thing as the Jacobi-Trudi formula). Then as a generalization of KP hierarchy, we construct the $\pi$-type KP hierarchy and obtain its tau functions.


title: H\"older regularity for the time fractional Schr\"odinger equation

abstract: In this paper, we investigate that the H\"older regularity of solutions to the time fractional Schr\"odinger equation of order $1<\alpha<2$, which interpolates between the Schr\"odinger and wave equations. This is inspired by Hirata and Miao's work which studied the fractional diffusion-wave equation. First, we give the asymptotic behavior for the oscillatory distributional kernels and their Bessel potentials by using Fourier analytic techniques. Then, the space regularity is derived by employing some results on singular Fourier multipliers. Using the asymptotic behavior for the above kernels, we prove the time regularity. Finally, we use mismatch estimates to prove the pointwise convergence to the initial data in H\"older spaces. In addition, we also prove H\"older regularity result for the Schr\"odinger equation.


title: Fundamental limits of quantum-secure covert communication over bosonic channels

abstract: We investigate the fundamental limit of quantum-secure covert communication over the lossy thermal noise bosonic channel, the quantum-mechanical model underlying many practical channels. We assume that the adversary has unlimited quantum information processing capabilities as well as access to all transmitted photons that do not reach the legitimate receiver. Given existence of noise that is uncontrolled by the adversary, the square root law (SRL) governs covert communication: up to c*sqrt{n} covert bits can be transmitted reliably in n channel uses. Attempting to surpass this limit results in detection with unity probability as n approaches infinity. Here we present the expression for c, characterizing the SRL for the bosonic channel. We also prove that discrete-valued coherent state quadrature phase shift keying (QPSK) constellation achieves the optimal c, which is the same as that achieved by a circularly-symmetric complex-valued Gaussian prior on coherent state amplitude. Finally, while binary phase shift keying (BPSK) achieves the Holevo capacity for non-covert bosonic channels in the low received signal-to-noise ratio regime, we show that it is strictly sub-optimal for covert communication.


title: A constrained pressure-temperature residual (CPTR) method for non-isothermal multiphase flow in porous media

abstract: For both isothermal and thermal petroleum reservoir simulation, the Constrained Pressure Residual (CPR) method is the industry-standard preconditioner. This method is a two-stage process involving the solution of a restricted pressure system. While initially designed for the isothermal case, CPR is also the standard for thermal cases. However, its treatment of the energy conservation equation does not incorporate heat diffusion, which is often dominant in thermal cases. In this paper, we present an extension of CPR: the Constrained Pressure-Temperature Residual (CPTR) method, where a restricted pressure-temperature system is solved in the first stage. In previous work, we introduced a block preconditioner with an efficient Schur complement approximation for a pressure-temperature system. Here, we extend this method for multiphase flow as the first stage of CPTR. The algorithmic performance of different two-stage preconditioners is evaluated for reservoir simulation test cases.


title: Tax- and expense-modified risk-minimization for insurance payment processes

abstract: We study the problem of determining risk-minimizing investment strategies for insurance payment processes in the presence of taxes and expenses. We consider the situation where taxes and expenses are paid continuously and symmetrically and introduce the concept of tax- and expense-modified risk-minimization. Risk-minimizing strategies in the presence of taxes and expenses are derived and linked to Galtchouk-Kunita-Watanabe decompositions associated with modified versions of the original payment processes. Furthermore, we show equivalence to an alternative approach involving an artificial market consisting of after-tax and after-expense assets, and we establish a type of consistency with classic risk-minimization. Finally, a case study involving classic multi-state life insurance payments in combination with a bond market exemplifies the results.


title: Synthetic dimensions and topological chiral currents in mesoscopic rings

abstract: The recently-introduced concept of "synthetic dimensions" allows for the realization of higher-dimensional topological phenomena in lower-dimensional systems. In this work we study the complementary aspect that synthetic dimensions provide a natural route to topological states in mesoscopic hybrid devices. We demonstrate this for the current induced into a closed one-dimensional Aharonov-Bohm ring by the interaction with a dynamic mesoscopic magnet. The quantization of the magnetic moment provides a synthetic dimension that complements the charge motion around the ring. We present a direct mapping that places the combined ring-magnet system into the class of quantum Hall models, and demonstrate that topological features, combined with the magnet's anisotropy, can lead to clear signatures in the persistent current of the single-particle ground state.


title: Unified Optimal Analysis of the (Stochastic) Gradient Method

abstract: In this note we give a simple proof for the convergence of stochastic gradient (SGD) methods on $\mu$-convex functions under a (milder than standard) $L$-smoothness assumption. We show that for carefully chosen stepsizes SGD converges after $T$ iterations as $O\left( LR^2 \exp \bigl[-\frac{\mu}{4L}T\bigr] + \frac{\sigma^2}{\mu T} \right)$ where $\sigma^2$ measures the variance in the stochastic noise. For deterministic gradient descent (GD) and SGD in the interpolation setting we have $\sigma^2 =0$ and we recover the exponential convergence rate. The bound matches with the best known iteration complexity of GD and SGD, up to constants.


title: Contextual One-Class Classification in Data Streams

abstract: In machine learning, the one-class classification problem occurs when training instances are only available from one class. It has been observed that making use of this class's structure, or its different contexts, may improve one-class classifier performance. Although this observation has been demonstrated for static data, a rigorous application of the idea within the data stream environment is lacking. To address this gap, we propose the use of context to guide one-class classifier learning in data streams, paying particular attention to the challenges presented by the dynamic learning environment. We present three frameworks that learn contexts and conduct experiments with synthetic and benchmark data streams. We conclude that the paradigm of contexts in data streams can be used to improve the performance of streaming one-class classifiers.


title: Weak field limit of higher dimensional massive Brans-Dicke gravity: Observational constraints

abstract: We consider higher-dimensional massive Brans-Dicke theory with Ricci-flat internal space. The background model is perturbed by a massive gravitating source which is pressureless in the external (our space) but has an arbitrary equation-of-state parameter $\Omega$ in the internal space. We obtain the exact solution of the system of linearized equations for the perturbations of the metric coefficients and scalar field. For a massless scalar field, relying on the fine-tuning between the Brans-Dicke parameter $\omega$ and $\Omega$, we demonstrate that (i) the model does not contradict gravitational tests relevant to the parameterized post-Newtonian parameter $\gamma$, and (ii) the scalar field is not ghost in the case of nonzero $|\Omega|\sim O(1)$ along with the natural value $|\omega|\sim O(1)$. In the general case of a massive scalar field, the metric coefficients acquire the Yukawa correction terms, where the Yukawa mass scale $m$ is defined by the mass of the scalar field. For the natural value $\omega\sim O(1)$, the inverse-square-law experiments impose the following restriction on the lower bound of the mass: $m\gtrsim 10^{-11}\,$GeV. The experimental constraints on $\gamma$ requires that $\Omega$ must be extremely close to $-1/2$.


title: A Simple Derivation of AMP and its State Evolution via First-Order Cancellation

abstract: We consider the linear regression problem, where the goal is to recover the vector $\boldsymbol{x}\in\mathbb{R}^n$ from measurements $\boldsymbol{y}=\boldsymbol{A}\boldsymbol{x}+\boldsymbol{w}\in\mathbb{R}^m$ under known matrix $\boldsymbol{A}$ and unknown noise $\boldsymbol{w}$. For large i.i.d. sub-Gaussian $\boldsymbol{A}$, the approximate message passing (AMP) algorithm is precisely analyzable through a state-evolution (SE) formalism, which furthermore shows that AMP is Bayes optimal in certain regimes. The rigorous SE proof, however, is long and complicated. And, although the AMP algorithm can be derived as an approximation of loop belief propagation (LBP), this viewpoint provides little insight into why large i.i.d. $\boldsymbol{A}$ matrices are important for AMP, and why AMP has a state evolution. In this work, we provide a heuristic derivation of AMP and its state evolution, based on the idea of "first-order cancellation," that provides insights missing from the LBP derivation while being much shorter than the rigorous SE proof.


title: Primordial black holes from the preheating instability

abstract: After the end of inflation, the inflaton field oscillates around a local minimum of its potential and decays into ordinary matter. These oscillations trigger a resonant instability for cosmological perturbations with wavelengths that exit the Hubble radius close to the end of inflation. In this paper, we study the formation of Primordial Black Holes (PBHs) at these enhanced scales. We find that the production mechanism can be so efficient that PBHs subsequently dominate the content of the universe and reheating proceeds from their evaporation. Observational constraints on the PBH abundance also restrict the duration of the resonant instability phase, leading to tight limits on the reheating temperature that we derive. We conclude that the production of PBHs during reheating is a generic and inevitable property of the simplest inflationary models, and does not require any fine tuning of the inflationary potential.


title: On the origins of extreme wealth inequality in the Talent vs Luck Model

abstract: While wealth distribution in the world is highly skewed and heavy-tailed, human talent - as the majority of individual features - is normally distributed. In a recent computational study by Pluchino et al [Talent vs luck: The role of randomness in success and failure, Adv. Complex Syst. 21 (03-04) (2018) 1850014], it has been shown that the combined effects of both random external factors (lucky and unlucky events) and multiplicative dynamics in capital accumulation are able to clarify this apparent contradiction. We introduce here a simplified version (STvL) of the original Talent versus Luck (TvL) model, where only lucky events are present, and verify that its dynamical rules lead to the same very large wealth inequality as the original model. We also derive some analytical approximations aimed to capture the mechanism responsible for the creation of such wealth inequality from a Gaussian-distributed talent. Under these approximations, our analysis is able to reproduce quite well the results of the numerical simulations of the simplified model in special cases. On the other hand, it also shows that the complexity of the model lies in the fact that lucky events are transformed into an increase of capital with heterogeneous rates, which yields a non-trivial generalization of the role of multiplicative processes in generating wealth inequality, whose fully generic case is still not amenable to analytical computations.


title: Lightening the Dark Matter from its Viscosity and Explanation of EDGES Anomaly

abstract: We study the visible photon production from the viscous dissipation of the dark matter (DM) fluid. The visible photon production depends on the magnitude of the dark matter viscosity and becomes important at the late times. We argue that for sufficiently large dark matter viscosity, the number of the resonantly converted visible photons becomes large which populates the Rayleigh-Jeans (RJ) tails of the Cosmic Microwave Background (CMB) radiation. Consequently, these excess visible photons possibly can explain the reported EDGES anomaly in the 21 cm signal. Further, we explore the parameter space for which the 21 cm signal can provide the region to probe the dark radiation and the DM viscosity.


title: Channel Impulse Response-based Source Localization in a Diffusion-based Molecular Communication System

abstract: This work localizes a molecular source in a diffusion based molecular communication (DbMC) system via a set of passive sensors and a fusion center. Molecular source localization finds its applications in future healthcare systems, including proactive diagnostics. In this paper, we propose two distinct methods which both utilize (the peak of) the channel impulse response measurements to uniquely localize the source, under assumption that the molecular source of interest lies within the open convex-hull of the sensor/anchor nodes. The first method is a one-shot, triangulation-based approach which estimates the unknown location of the molecular source using least-squares method. The corresponding Cramer-Rao bound (CRB) is also derived. The second method is an iterative approach, which utilizes gradient descent law to minimize a non-convex cost function. Simulation results reveal that the triangulation-based method performs very close to the CRB, for any given signal- to-noise ratio. Additionally, the gradient descent-based method converges to the true optima/source location uniformly (in less than hundred iterations).


title: Bayesian deep learning with hierarchical prior: Predictions from limited and noisy data

abstract: Datasets in engineering applications are often limited and contaminated, mainly due to unavoidable measurement noise and signal distortion. Thus, using conventional data-driven approaches to build a reliable discriminative model, and further applying this identified surrogate to uncertainty analysis remains to be very challenging. A deep learning approach is presented to provide predictions based on limited and noisy data. To address noise perturbation, the Bayesian learning method that naturally facilitates an automatic updating mechanism is considered to quantify and propagate model uncertainties into predictive quantities. Specifically, hierarchical Bayesian modeling (HBM) is first adopted to describe model uncertainties, which allows the prior assumption to be less subjective, while also makes the proposed surrogate more robust. Next, the Bayesian inference is seamlessly integrated into the DL framework, which in turn supports probabilistic programming by yielding a probability distribution of the quantities of interest rather than their point estimates. Variational inference (VI) is implemented for the posterior distribution analysis where the intractable marginalization of the likelihood function over parameter space is framed in an optimization format, and stochastic gradient descent method is applied to solve this optimization problem. Finally, Monte Carlo simulation is used to obtain an unbiased estimator in the predictive phase of Bayesian inference, where the proposed Bayesian deep learning (BDL) scheme is able to offer confidence bounds for the output estimation by analyzing propagated uncertainties. The effectiveness of Bayesian shrinkage is demonstrated in improving predictive performance using contaminated data, and various examples are provided to illustrate concepts, methodologies, and algorithms of this proposed BDL modeling technique.


title: CHOP: Bypassing Runtime Bounds Checking Through Convex Hull OPtimization

abstract: Unsafe memory accesses in programs written using popular programming languages like C/C++ have been among the leading causes for software vulnerability. Prior memory safety checkers such as SoftBound enforce memory spatial safety by checking if every access to array elements are within the corresponding array bounds. However, it often results in high execution time overhead due to the cost of executing the instructions associated with bounds checking. To mitigate this problem, redundant bounds check elimination techniques are needed. In this paper, we propose CHOP, a Convex Hull OPtimization based framework, for bypassing redundant memory bounds checking via profile-guided inferences. In contrast to existing check elimination techniques that are limited by static code analysis, our solution leverages a model-based inference to identify redundant bounds checking based on runtime data from past program executions. For a given function, it rapidly derives and updates a knowledge base containing sufficient conditions for identifying redundant array bounds checking. We evaluate CHOP on real-world applications and benchmark (such as SPEC) and the experimental results show that on average 80.12% of dynamic bounds check instructions can be avoided, resulting in improved performance up to 95.80% over SoftBound.


title: Topological Information Data Analysis

abstract: This paper presents methods that quantify the structure of statistical interactions within a given data set, and was first used in \cite{Tapia2018}. It establishes new results on the k-multivariate mutual-informations (I_k) inspired by the topological formulation of Information introduced in. In particular we show that the vanishing of all I_k for 2\leq k \leq n of n random variables is equivalent to their statistical independence. Pursuing the work of Hu Kuo Ting and Te Sun Han, we show that information functions provide co-ordinates for binary variables, and that they are analytically independent on the probability simplex for any set of finite variables. The maximal positive I_k identifies the variables that co-vary the most in the population, whereas the minimal negative I_k identifies synergistic clusters and the variables that differentiate-segregate the most the population. Finite data size effects and estimation biases severely constrain the effective computation of the information topology on data, and we provide simple statistical tests for the undersampling bias and the k-dependences following. We give an example of application of these methods to genetic expression and unsupervised cell-type classification. The methods unravel biologically relevant subtypes, with a sample size of 41 genes and with few errors. It establishes generic basic methods to quantify the epigenetic information storage and a unified epigenetic unsupervised learning formalism. We propose that higher-order statistical interactions and non identically distributed variables are constitutive characteristics of biological systems that should be estimated in order to unravel their significant statistical structure and diversity. The topological information data analysis presented here allows to precisely estimate this higher-order structure characteristic of biological systems.


title: The Combinatorics of Barrier Synchronization

abstract: In this paper we study the notion of synchronization from the point of view of combinatorics. As a first step, we address the quantitative problem of counting the number of executions of simple processes interacting with synchronization barriers. We elaborate a systematic decomposition of processes that produces a symbolic integral formula to solve the problem. Based on this procedure, we develop a generic algorithm to generate process executions uniformly at random. For some interesting sub-classes of processes we propose very efficient counting and random sampling algorithms. All these algorithms have one important characteristic in common: they work on the control graph of processes and thus do not require the explicit construction of the state-space.


title: Search for Lorentz Violation Using High-Energy Atmospheric Neutrinos In IceCube

abstract: High-energy atmospheric neutrinos observed by the IceCube Neutrino Observatory are extremely sensitive probes of Lorentz violation (LV). Here we report the result of analyzing two years of IceCube data in the search for LV. This analysis places some of the strongest constraints on LV when considering high-dimensional operators.


title: ICLab: A Global, Longitudinal Internet Censorship Measurement Platform

abstract: Researchers have studied Internet censorship for nearly as long as attempts to censor contents have taken place. Most studies have however been limited to a short period of time and/or a few countries; the few exceptions have traded off detail for breadth of coverage. Collecting enough data for a comprehensive, global, longitudinal perspective remains challenging. In this work, we present ICLab, an Internet measurement platform specialized for censorship research. It achieves a new balance between breadth of coverage and detail of measurements, by using commercial VPNs as vantage points distributed around the world. ICLab has been operated continuously since late 2016. It can currently detect DNS manipulation and TCP packet injection, and overt "block pages" however they are delivered. ICLab records and archives raw observations in detail, making retrospective analysis with new techniques possible. At every stage of processing, ICLab seeks to minimize false positives and manual validation. Within 53,906,532 measurements of individual web pages, collected by ICLab in 2017 and 2018, we observe blocking of 3,602 unique URLs in 60 countries. Using this data, we compare how different blocking techniques are deployed in different regions and/or against different types of content. Our longitudinal monitoring pinpoints changes in censorship in India and Turkey concurrent with political shifts, and our clustering techniques discover 48 previously unknown block pages. ICLab's broad and detailed measurements also expose other forms of network interference, such as surveillance and malware injection.


title: Security for Distributed Deep Neural Networks Towards Data Confidentiality & Intellectual Property Protection

abstract: Current developments in Enterprise Systems observe a paradigm shift, moving the needle from the backend to the edge sectors of those; by distributing data, decentralizing applications and integrating novel components seamlessly to the central systems. Distributively deployed AI capabilities will thrust this transition. Several non-functional requirements arise along with these developments, security being at the center of the discussions. Bearing those requirements in mind, hereby we propose an approach to holistically protect distributed Deep Neural Network (DNN) based/enhanced software assets, i.e. confidentiality of their input & output data streams as well as safeguarding their Intellectual Property. Making use of Fully Homomorphic Encryption (FHE), our approach enables the protection of Distributed Neural Networks, while processing encrypted data. On that respect we evaluate the feasibility of this solution on a Convolutional Neuronal Network (CNN) for image classification deployed on distributed infrastructures.


title: Dynamical low-rank integrator for the linear Boltzmann equation: error analysis in the diffusion limit

abstract: Dynamical low-rank algorithms are a class of numerical methods that compute low-rank approximations of dynamical systems. This is accomplished by projecting the dynamics onto a low-dimensional manifold and writing the solution directly in terms of the low-rank factors. The approach has been successfully applied to many types of differential equations. Recently, efficient dynamical low-rank algorithms have been applied to treat kinetic equations, including the Vlasov--Poisson and the Boltzmann equation, where it was demonstrated that the methods are able to capture the low-rank structure of the solution and significantly reduce numerical cost, while often maintaining high accuracy. However, no numerical analysis is currently available. In this paper, we investigate the error analysis for a dynamical low-rank algorithm applied to the multi-scale linear Boltzmann equation (a classical model in kinetic theory) to showcase the validity of the application of dynamical low-rank algorithms to kinetic theory. The equation, in its parabolic regime, is known to be rank one theoretically, and we will prove that the scheme can dynamically and automatically capture this low-rank structure. This work thus serves as the first mathematical error analysis for a dynamical low-rank approximation applied to a kinetic problem.


title: Helicity amplitudes, polarization of EM waves and Stokes parameters

abstract: It is shown that the helicity amplitudes can be used to describe and analyze the properties of the electromagnetic field in classical and in quantum theory. On the one hand they embody the relativistic content of electromagnetic theory. On the other hand they give a concise description of such experimentally important notions as polarization, the Stokes parameters and the Poincar\'e sphere.


title: Analysis of incubation time preceding the Ga-assisted nucleation and growth of GaAs nanowires on Si(111)

abstract: The incubation time preceding nucleation and growth of surface nanostructures is interesting from a fundamental viewpoint but also of practical relevance as it determines statistical properties of nanostructure ensembles such as size homogeneity. Using in situ reflection high-energy electron diffraction, we accurately deduce the incubation times for Ga-assisted GaAs nanowires grown on unpatterned Si(111) substrates by molecular beam epitaxy under different conditions. We develop a nucleation model that explains and fits very well the data. We find that, for a given temperature and Ga flux, the incubation time always increases with decreasing As flux and becomes infinite at a certain minimum flux, which is larger for higher temperature. For given As and Ga fluxes, the incubation time always increases with temperature and rapidly tends to infinity above 640 {\deg}C under typical conditions. The strong temperature dependence of the incubation time is reflected in a similar variation of the nanowire number density with temperature. Our analysis provides understanding and guidance for choosing appropriate growth conditions that avoid unnecessary material consumption, long nucleation delays, and highly inhomogeneous ensembles of nanowires. On a more general ground, the existence of a minimum flux and maximum temperature for growing surface nanostructures should be a general phenomenon pertaining for a wide range of material-substrate combinations.


title: Singular limits of the quasi-linear Kolmogorov-type equation with a source term

abstract: Existence, uniqueness and stability of kinetic and entropy solutions to the boundary value problem for the Kolmogorov-type genuinely nonlinear ultra-parabolic equation with a smooth source term is established. After this, we consider the case when the source term contains a small positive parameter and collapses to the Dirac delta-function, as this parameter tends to zero. In this case, the limiting passage from the original equation with the smooth source to the impulsive ultra-parabolic equation is fulfilled and rigorously justified. The proofs rely on the method of kinetic equation and on the compensated compactness techniques for genuinely nonlinear equations.


title: A divide-and-conquer algorithm for binary matrix completion

abstract: We propose an algorithm for low rank matrix completion for matrices with binary entries which obtains explicit binary factors. Our algorithm, which we call TBMC (\emph{Tiling for Binary Matrix Completion}), gives interpretable output in the form of binary factors which represent a decomposition of the matrix into tiles. Our approach is inspired by a popular algorithm from the data mining community called PROXIMUS: it adopts the same recursive partitioning approach while extending to missing data. The algorithm relies upon rank-one approximations of incomplete binary matrices, and we propose a linear programming (LP) approach for solving this subproblem. We also prove a $2$-approximation result for the LP approach which holds for any level of subsampling and for any subsampling pattern. Our numerical experiments show that TBMC outperforms existing methods on recommender systems arising in the context of real datasets.


title: The Secretary Recommendation Problem

abstract: In this paper we revisit the basic variant of the classical secretary problem. We propose a new approach in which we separate between an agent that evaluates the secretary performance and one that has to make the hiring decision. The evaluating agent (the sender) signals the quality of the candidate to the hiring agent (the receiver) who must make a decision. Whenever the two agents' interests are not fully aligned, this induces an information transmission (signaling) challenge for the sender. We study the sender's optimization problem subject to persuasiveness constraints of the receiver for several variants of the problem. Our results quantify the loss in performance for the sender due to online arrival. We provide optimal and near-optimal persuasive mechanisms that recover at least a constant fraction of a natural utility benchmark for the sender. The separation of evaluation and decision making can have a substantial impact on the approximation results. While in some scenarios, techniques and results closely mirror the conditions in the standard secretary problem, we also explore conditions that lead to very different characteristics.


title: Gated Multiple Feedback Network for Image Super-Resolution

abstract: The rapid development of deep learning (DL) has driven single image super-resolution (SR) into a new era. However, in most existing DL based image SR networks, the information flows are solely feedforward, and the high-level features cannot be fully explored. In this paper, we propose the gated multiple feedback network (GMFN) for accurate image SR, in which the representation of low-level features are efficiently enriched by rerouting multiple high-level features. We cascade multiple residual dense blocks (RDBs) and recurrently unfolds them across time. The multiple feedback connections between two adjacent time steps in the proposed GMFN exploits multiple high-level features captured under large receptive fields to refine the low-level features lacking enough contextual information. The elaborately designed gated feedback module (GFM) efficiently selects and further enhances useful information from multiple rerouted high-level features, and then refine the low-level features with the enhanced high-level information. Extensive experiments demonstrate the superiority of our proposed GMFN against state-of-the-art SR methods in terms of both quantitative metrics and visual quality. Code is available at https://github.com/liqilei/GMFN.


title: Arbitrarily High-order Unconditionally Energy Stable Schemes for Gradient Flow Models Using the Scalar Auxiliary Variable Approach

abstract: In this paper, we propose a novel family of high-order numerical schemes for the gradient flow models based on the scalar auxiliary variable (SAV) approach, which is named the high-order scalar auxiliary variable (HSAV) method. The newly proposed schemes could be shown to reach arbitrarily high order in time while preserving the energy dissipation law without any restriction on the time step size (i.e., unconditionally energy stable). The HSAV strategy is rather general that it does not depend on the specific expression of the effective free energy, such that it applies to a class of thermodynamically consistent gradient flow models arriving at semi-discrete high-order energy-stable schemes. We then employ the Fourier pseudospectral method for spatial discretization. The fully discrete schemes are also shown to be unconditionally energy stable. Furthermore, we present several numerical experiments on several widely-used gradient flow models, to demonstrate the accuracy, efficiency and unconditionally energy stability of the HSAV schemes. The numerical results verify that the HSAV schemes can reach the expected order of accuracy, and it allows a much larger time step size to reach the same accuracy than the standard SAV schemes.


title: Quantum-critical phase out of frustrated magnetism in a strongly correlated metal

abstract: Strange-metal phenomena often develop at the border of antiferromagnetic order in strongly correlated metals. It has been well established that they can originate from the fluctuations anchored by the point of continuous quantum phase transition out of the antiferromagnetic order, i.e., a quantum critical point. What has been unclear is how these phenomena can be associated with a potential new phase of matter at zero temperature. Here we show that magnetic frustration of the 4f-local moments in the distorted Kagome intermetallic compound CePdAl gives rise to such a paramagnetic quantum-critical phase. Moreover, we demonstrate that this phase turns into a Fermi liquid through a Mott-like crossover; in a two-dimensional parameter space of pressure and magnetic field, this crossover is linked to a line of zero-temperature 4f-electron localization-delocalization phase transitions at low and moderate pressures. Our discovery motivates a new design principle for strongly correlated metallic states with unconventional excitations that may underlie the development of such effects as high temperature superconductivity.


title: Singling Out Dynamic and Nondynamic Correlation

abstract: The correlation part of the pair density is separated into two components, one of them being predominant at short electronic ranges and the other at long ranges. The analysis of the intracular part of these components permits to classify molecular systems according to the prevailing correlation: dynamic or nondynamic. The study of the long-range asymptotics reveals the key component of the pair density that is responsible for the description of London dispersion forces and a universal decay with the interelectronic distance. The natural range-separation, the identification of the dispersion forces, and the kind of predominant correlation type that arise from this analysis are expected to be important assets in the development of new electronic structure methods in wave function, density, and reduced density-matrix functional theories.


abstract: We propose a novel concept of a Systemic Optimal Risk Transfer Equilibrium (SORTE), which is inspired by the B\"uhlmann's classical notion of an Equilibrium Risk Exchange. We provide sufficient general assumptions that guarantee existence, uniqueness, and Pareto optimality of such a SORTE. In both the B\"uhlmann and the SORTE definition, each agent is behaving rationally by maximizing his/her expected utility given a budget constraint. The two approaches differ by the budget constraints. In B\"uhlmann's definition the vector that assigns the budget constraint is given a priori. On the contrary, in the SORTE approach, the vector that assigns the budget constraint is endogenously determined by solving a systemic utility maximization. SORTE gives priority to the systemic aspects of the problem, in order to optimize the overall systemic performance, rather than to individual rationality.


title: Melody Generation using an Interactive Evolutionary Algorithm

abstract: Music generation with the aid of computers has been recently grabbed the attention of many scientists in the area of artificial intelligence. Deep learning techniques have evolved sequence production methods for this purpose. Yet, a challenging problem is how to evaluate generated music by a machine. In this paper, a methodology has been developed based upon an interactive evolutionary optimization method, with which the scoring of the generated melodies is primarily performed by human expertise, during the training. This music quality scoring is modeled using a Bi-LSTM recurrent neural network. Moreover, the innovative generated melody through a Genetic algorithm will then be evaluated using this Bi-LSTM network. The results of this mechanism clearly show that the proposed method is able to create pleasurable melodies with desired styles and pieces. This method is also quite fast, compared to the state-of-the-art data-oriented evolutionary systems.


title: Zeros of certain combinations of Eisenstein series of weight 2k, 3k, and k + l

abstract: We locate the zeros of the modular forms $E_k^2(\tau) + E_{2k}(\tau), E_k^3(\tau) + E_{3k} (\tau),$ and $E_k(\tau)E_l(\tau) +E_{k+l}(\tau),$ where $E_k(\tau)$ is the Eisenstein series for the full modular group $\text{SL}_2(\mathbb{Z})$. By utilizing work of F.K.C. Rankin and Swinnerton-Dyer, we prove that for sufficiently large $k,l$, all zeros in the standard fundamental domain are located on the lower boundary $\mathcal{A} = \{ e^{i\theta} : \pi/2 \leq \theta \leq 2\pi/3\}$.


title: Affinity-dependent bound on the spectrum of stochastic matrices

abstract: Affinity has proven to be a useful tool for quantifying the non-equilibrium character of time continuous Markov processes since it serves as a measure for the breaking of time reversal symmetry. It has recently been conjectured that the number of coherent oscillations, which is given by the ratio of imaginary and real part of the first non-trivial eigenvalue of the corresponding master matrix, is constrained by the maximum cycle affinity present in the network. In this paper, we conjecture a bound on the whole spectrum of these master matrices that constrains all eigenvalues in a fashion similar to the well known Perron-Frobenius theorem that is valid for any stochastic matrix. As in other studies that are based on affinity-dependent bounds, the limiting process that saturates the bound is given by the asymmetric random walk. For unicyclic networks, we prove that it is not possible to violate the bound by small perturbation of the asymmetric random walk and provide numerical evidence for its validity in randomly generated networks. The results are extended to multicyclic networks, backed up by numerical evidence provided by networks with randomly constructed topology and transition rates.


title: Heating Rates in Periodically Driven Strongly Interacting Quantum Many-Body Systems

abstract: We study heating rates in strongly interacting quantum lattice systems in the thermodynamic limit. Using a numerical linked cluster expansion, we calculate the energy as a function of the driving time and find a robust exponential regime. The heating rates are shown to be in excellent agreement with Fermi's golden rule. We discuss the relationship between heating rates and, within the eigenstate thermalization hypothesis, the smooth function that characterizes the off-diagonal matrix elements of the drive operator in the eigenbasis of the static Hamiltonian. We show that such a function, in nonintegrable and (remarkably) integrable Hamiltonians, can be probed experimentally by studying heating rates as functions of the drive frequency.


title: solc-verify: A Modular Verifier for Solidity Smart Contracts

abstract: We present solc-verify, a source-level verification tool for Ethereum smart contracts. Solc-verify takes smart contracts written in Solidity and discharges verification conditions using modular program analysis and SMT solvers. Built on top of the Solidity compiler, solc-verify reasons at the level of the contract source code, as opposed to the more common approaches that operate at the level of Ethereum bytecode. This enables solc-verify to effectively reason about high-level contract properties while modeling low-level language semantics precisely. The contract properties, such as contract invariants, loop invariants, and function pre- and post-conditions, can be provided as annotations in the code by the developer. This enables automated, yet user-friendly formal verification for smart contracts. We demonstrate solc-verify by examining real-world examples where our tool can effectively find bugs and prove correctness of non-trivial properties with minimal user effort.


title: Genuine multipartite correlations in Dicke Superradiance

abstract: A thorough understanding of the structure of correlations in multipartite systems is essential for the success of most quantum information processing applications. The problem, however, quickly becomes non-trivial as the size of the multipartition increases. With this motivation in mind, in this paper, we put forth a detailed study of genuine multipartite correlations (GMCs) in the Dicke model of superradiance. We compute all genuine $k$-partite correlations for Dicke states with arbitrary excitations and use these results to characterize the evolution of multipartite correlations during the superradiant dynamics. Non-trivial effects in the way correlations in Dicke states are distributed between the multiple parts are found, showing strong finite-size effects. We also employ the concept of weaving to classify how multipartite correlations scale with the number of particles in the system.


title: Optimal estimation of parameters encoded in coherent states quadratures

abstract: In the context of multiparameter quantum estimation theory, we investigate the construction of linear schemes in order to infer two classical parameters that are encoded in the quadratures of two quantum coherent states. The optimality of the scheme built on two phase-conjugate coherent states is proven with the saturation of the quantum Cram\'er--Rao bound under some global energy constraint. In a more general setting, we consider and analyze a variety of $n$-mode schemes that can be used to encode $n$ classical parameters into $n$ quantum coherent states and then estimate all parameters optimally and simultaneously.


title: Translating neural signals to text using a Brain-Machine Interface

abstract: Brain-Computer Interfaces (BCI) help patients with faltering communication abilities due to neurodegenerative diseases produce text or speech output by direct neural processing. However, practical implementation of such a system has proven difficult due to limitations in speed, accuracy, and generalizability of the existing interfaces. To this end, we aim to create a BCI system that decodes text directly from neural signals. We implement a framework that initially isolates frequency bands in the input signal encapsulating differential information regarding production of various phonemic classes. These bands then form a feature set that feeds into an LSTM which discerns at each time point probability distributions across all phonemes uttered by a subject. Finally, these probabilities are fed into a particle filtering algorithm which incorporates prior knowledge of the English language to output text corresponding to the decoded word. Performance of this model on data obtained from six patients shows encouragingly high levels of accuracy at speeds and bit rates significantly higher than existing BCI communication systems. Further, in producing an output, our network abstains from constraining the reconstructed word to be from a given bag-of-words, unlike previous studies. The success of our proposed approach, offers promise for the employment of a BCI interface by patients in unfettered, naturalistic environments.


title: Identifying the Influential Inputs for Network Output Variance Using Sparse Polynomial Chaos Expansion

abstract: Sensitivity analysis (SA) is an important aspect of process automation. It often aims to identify the process inputs that influence the process output's variance significantly. Existing SA approaches typically consider the input-output relationship as a black-box and conduct extensive random sampling from the actual process or its high-fidelity simulation model to identify the influential inputs. In this paper, an alternate, novel approach is proposed using a sparse polynomial chaos expansion-based model for a class of input-output relationships represented as directed acyclic networks. The model exploits the relationship structure by recursively relating a network node to its direct predecessors to trace the output variance back to the inputs. It, thereby, estimates the Sobol indices, which measure the influence of each input on the output variance, accurately and efficiently. Theoretical analysis establishes the validity of the model as the prediction of the network output converges in probability to the true output under certain regularity conditions. Empirical evaluation on two manufacturing processes shows that the model estimates the Sobol indices accurately with far fewer observations than a state-of-the-art Monte Carlo sampling method.


title: Magnetic instabilities in doped Fe$_2YZ$ full-Heusler thermoelectric compounds

abstract: Thermoelectricity is a promising avenue for harvesting energy but large-scale applications are still hampered by the lack of highly-efficient low-cost materials. Recently, Fe$_2YZ$ Heusler compounds were predicted theoretically to be interesting candidates with large thermoelectric power factor. Here, we show that under doping conditions compatible with thermoelectric applications, these materials are prone to an unexpected magnetic instability detrimental to their thermoelectric performance. We rationalize the physics at the origin of this instability, provide guidelines for avoiding it and discuss its impact on the thermoelectric power factor. Doing so, we also point out the shortcomings of the rigid band approximation commonly used in high-throughput theoretical searches of new thermoelectrics.


title: de Sitter Cosmology on an expanding bubble

abstract: Constructing an explicit compactification yielding a metastable de Sitter (dS) vacuum in a UV consistent string theory is an incredibly difficult open problem. Motivated by this issue, as well as the conjecture that all non-supersymmetric AdS vacua must decay, we discuss the alternative possibility of realizing an effective four-dimensional dS cosmology on a codimension-one bubble wall separating two AdS$_5$ vacua. The construction further elaborates on the scenario of arXiv:1807.01570, where the aforementioned cosmology arises due to a non-perturbative decay and is embedded in a five-dimensional bulk in a time-dependent way. In this paper we discuss the relation between this scenario and the weak gravity conjecture and further develop the details of the four-dimensional cosmology. We provide a bulk interpretation for the dS temperature as the Unruh temperature experienced by an accelerated observer riding the bubble. A source of four-dimensional matter arises from a string cloud in the bulk, and we examine the consequences for the particle mass spectrum. Furthermore, we show how effective four-dimensional Einstein gravity on the bubble is obtained from the five-dimensional Gauss equation. We conclude by outlining some implications that this paradigm will have for holography, inflation, the standard model, and black holes.


title: A Scheme for Dynamic Risk-Sensitive Sequential Decision Making

abstract: We present a scheme for sequential decision making with a risk-sensitive objective and constraints in a dynamic environment. A neural network is trained as an approximator of the mapping from parameter space to space of risk and policy with risk-sensitive constraints. For a given risk-sensitive problem, in which the objective and constraints are, or can be estimated by, functions of the mean and variance of return, we generate a synthetic dataset as training data. Parameters defining a targeted process might be dynamic, i.e., they might vary over time, so we sample them within specified intervals to deal with these dynamics. We show that: i). Most risk measures can be estimated using return variance; ii). By virtue of the state-augmentation transformation, practical problems modeled by Markov decision processes with stochastic rewards can be solved in a risk-sensitive scenario; and iii). The proposed scheme is validated by a numerical experiment.


title: Anomalous Goos-H\"anchen shift in the Floquet scattering of Dirac fermions

abstract: We study the inelastic scattering of two-dimensional massless Dirac fermions by an inhomogeneous time-dependent driving field. As a physical realization we consider a monolayer graphene normally illuminated with a circularly polarized laser of frequency $\Omega$ in a given region. The interaction Hamiltonian introduced by the laser, being periodic in time, can be treated with the Floquet method which naturally leads to a multi-channel scattering problem. We analyze planar and circular geometries of the interface separating the irradiated and non-irradiated regions and find that there is an anomalous Goos-H\"anchen shift in the inelastic channel. The latter is independent of the amplitude of the driving while its sign is determined by the polarization of the laser field. We related this shift with the appearance of topological edge states between two illuminated regions of opposite chiralities.


title: Ultra-light scalar saving the 3+1 neutrino scheme from the cosmological bounds

abstract: The LSND and MiniBooNE results as well as the reactor and Gallium anomalies seem to indicate the presence of a sterile neutrino with a mass of $\sim 1$ eV mixed with active neutrinos. Such sterile neutrino can be produced in the early universe before the neutrino decoupling, leading to a contribution to the effective number of neutrinos ($N_{eff}$) as well as to a contribution to the sum of neutrino masses which are in tension with cosmological observations. We propose a scenario to relax this tension by a Yukawa coupling of the sterile neutrinos to ultra-light scalar particles which contribute to the dark matter in the background. The coupling induces an effective mass for $\nu_s$ which prevents its production in the early universe. We discuss the implications for the upcoming KATRIN experiment and future relic neutrino search experiments such as PTOLEMY. We also briefly comment on certain non-renormalizable forms of interaction between $\nu_s$ and the scalar and their consequences for the $\nu_s$ production in the early universe.


title: Ordinal Imitative Dynamics

abstract: This paper introduces an evolutionary dynamics based on imitate the better realization (IBR) rule. Under this rule, agents in a population game imitate the strategy of a randomly chosen opponent whenever the opponent`s realized payoff is higher than their own. Such behavior generates an ordinal mean dynamics which is polynomial in strategy utilization frequencies. We demonstrate that while the dynamics does not possess Nash stationarity or payoff monotonicity, under it pure strategies iteratively strictly dominated by pure strategies are eliminated and strict equilibria are locally stable. We investigate the relationship between the dynamics based on the IBR rule and the replicator dynamics. In trivial cases, the two dynamics are topologically equivalent. In Rock-Paper-Scissors games we conjecture that both dynamics exhibit the same types of behavior, but the partitions of the game set do not coincide. In other cases, the IBR dynamics exhibits behaviors that are impossible under the replicator dynamics.


title: An MPC Approach to Transient Control of Liquid-Propellant Rocket Engines

abstract: The current context of launchers reusability requires the improvement of control algorithms for their liquid-propellant rocket engines. Their transient phases are generally still performed in open loop. In this paper, it is aimed at enhancing the control performance and robustness during the fully continuous phase of the start-up transient of a generic gas-generator cycle. The main control goals concern end-state tracking in terms of combustion-chamber pressure and chambers mixture ratios, as well as the verification of a set of hard operational constraints. A controller based on a nonlinear preprocessor and on linear MPC (Model-Predictive Control) has been synthesised, making use of nonlinear state-space models of the engine. The former generates the full-state reference to be tracked while the latter achieves the aforementioned goals with sufficient accuracy and verifying constraints for the required pressure levels. Robustness considerations are included in the MPC algorithm via an epigraph formulation of the minimax robust optimisation problem, where a finite set of perturbation scenarios is considered.


title: Reconstruction under outliers for Fourier-sparse functions

abstract: We consider the problem of learning an unknown $f$ with a sparse Fourier spectrum in the presence of outlier noise. In particular, the algorithm has access to a noisy oracle for (an unknown) $f$ such that (i) the Fourier spectrum of $f$ is $k$-sparse; (ii) at any query point $x$, the oracle returns $y$ such that with probability $1-\rho$, $|y-f(x)| \le \epsilon$. However, with probability $\rho$, the error $y-f(x)$ can be arbitrarily large. We study Fourier sparse functions over both the discrete cube $\{0,1\}^n$ and the torus $[0,1)$ and for both these domains, we design efficient algorithms which can tolerate any $\rho<1/2$ fraction of outliers. We note that the analogous problem for low-degree polynomials has recently been studied in several works~[AK03, GZ16, KKP17] and similar algorithmic guarantees are known in that setting. While our main results pertain to the case where the location of the outliers, i.e., $x$ such that $|y-f(x)|>\epsilon$ is randomly distributed, we also study the case where the outliers are adversarially located. In particular, we show that over the torus, assuming that the Fourier transform satisfies a certain \emph{granularity} condition, there is a sample efficient algorithm to tolerate $\rho =\Omega(1)$ fraction of outliers and further, that this is not possible without such a granularity condition. Finally, while not the principal thrust, our techniques also allow us non-trivially improve on learning low-degree functions $f$ on the hypercube in the presence of adversarial outlier noise. Our techniques combine a diverse array of tools from compressive sensing, sparse Fourier transform, chaining arguments and complex analysis.


title: Learning to Optimize Domain Specific Normalization for Domain Generalization

abstract: We propose a simple but effective multi-source domain generalization technique based on deep neural networks by incorporating optimized normalization layers that are specific to individual domains. Our approach employs multiple normalization methods while learning separate affine parameters per domain. For each domain, the activations are normalized by a weighted average of multiple normalization statistics. The normalization statistics are kept track of separately for each normalization type if necessary. Specifically, we employ batch and instance normalizations in our implementation to identify the best combination of these two normalization methods in each domain. The optimized normalization layers are effective to enhance the generalizability of the learned model. We demonstrate the state-of-the-art accuracy of our algorithm in the standard domain generalization benchmarks, as well as viability to further tasks such as multi-source domain adaptation and domain generalization in the presence of label noise.


title: A Conformance Checking-based Approach for Drift Detection in Business Processes

abstract: Real life business processes change over time, in both planned and unexpected ways. The detection of these changes is crucial for organizations to ensure that the expected and the real behavior are as similar as possible. These changes over time are called concept drift and its detection is a big challenge in process mining since the inherent complexity of the data makes difficult distinguishing between a change and an anomalous execution. In this paper, we present C2D2 (Conformance Checking-based Drift Detection), a new approach to detect sudden control-flow changes in the process models from event traces. C2D2 combines discovery techniques with conformance checking methods to perform an offline detection. Our approach has been validated with a synthetic benchmarking dataset formed by 68 logs, showing an improvement in the accuracy while maintaining a minimum delay in the drift detection.


title: The long term impact of Daylight Saving Time regulations in daily life at several circles of latitude

abstract: We analyze large scale (N~10000) time use surveys in United States, Spain, Italy, France and Great Britain to ascertain seasonal variations in the sleep/wake cycle and the labor cycle after daylight saving time -- summer time arrangements -- regulations have stood up for at least forty years. That is, not the usual search for short-term effects of the biannual transitions of DST but a search for how industrialized societies have responded to DST regulations. Results show that the labor cycle is equally distributed through seasons, an everyday experience which is a major outcome of DST. The sleep/wake cycle displays disturbance punctuated by solar events: sunrise, sunset and noon. Sleep onset in non-employees in a week-end -- the group more prone to free choices -- delays in summer opposing to the regulation. Yet, sleep offset in this group advances in summer, despite clock time already advanced in the spring transition. The best explanation for these findings is that size and directionality of DST regulations are close to an optimum. The role of latitude is discussed after sleep/wake cycle in Spain, Italy and United States exhibit fewer statistically significant excursions that the sleep/wake cycle in Great Britain and France, which could be indicating preference for time synchronization as latitude increases. The consequences of a seasonal regulation of clocks vs the choice of permanent summer time or permanent winter time is sketched from a previous report on human activity.


title: Nonlinear Electrophoresis of Highly Charged Nonpolarizable Particles

abstract: Nonlinear field dependence of electrophoresis in high fields has been investigated theoretically, yet experimental studies have failed to reach consensus on the effect. In this work, we present a systematic study on the nonlinear electrophoresis of highly charged submicron particles in applied electric fields of up to several kV/cm. First, the particles are characterized in the low-field regime at different salt concentrations and the surface charge density is estimated. Subsequently, we use microfluidic channels and video tracking to systematically characterize the nonlinear response over a range of field strengths. Using velocity measurements on the single particle level, we prove that nonlinear effects are present at electric fields and surface charge densities that are accessible in practical conditions. Finally, we show that nonlinear behavior leads to unexpected particle trapping in channels.


title: Multiple Knapsack-Constrained Monotone DR-Submodular Maximization on Distributive Lattice --- Continuous Greedy Algorithm on Median Complex ---

abstract: We consider a problem of maximizing a monotone DR-submodular function under multiple order-consistent knapsack constraints on a distributive lattice. Since a distributive lattice is used to represent a dependency constraint, the problem can represent a dependency constrained version of a submodular maximization problem on a set. We propose a $1 - 1/e$ approximation algorithm for this problem. To achieve this result, we generalize the continuous greedy algorithm to distributive lattices: We choose a median complex as a continuous relaxation of a distributive lattice and define the multilinear extension on it. We show that the median complex admits special curves, named uniform linear motions, such that the multilinear extension of a DR-submodular function is concave along a positive uniform linear motion, which is a key property of the continuous greedy algorithm.


title: Revisiting Biorthogonal Polynomials. An $LU$ factorization discussion

abstract: The Gauss-Borel or $LU$ factorization of Gram matrices of bilinear forms is the pivotal element in the discussion of the theory of biorthogonal polynomials. The construction of biorthogonal families of polynomials and its second kind functions, of the spectral matrices modeling the multiplication by the independent variable $x$, the Christoffel-Darboux kernel and its projection properties, are discussed from this point of view. Then, the Hankel case is presented and different properties, specific of this case, as the three terms relations, Heine formulas, Gauss quadrature and the Christoffel-Darboux formula are given. The classical orthogonal polynomial of Hermite, Laguerre and Jacobi type are discussed and characterized within this scheme. Finally, it is shown who this approach is instrumental in the derivation of Christoffel formulas for general Christoffel and Geronimus perturbations of the bilinear forms.


title: Deep Learning Techniques for Improving Digital Gait Segmentation

abstract: Wearable technology for the automatic detection of gait events has recently gained growing interest, enabling advanced analyses that were previously limited to specialist centres and equipment (e.g., instrumented walkway). In this study, we present a novel method based on dilated convolutions for an accurate detection of gait events (initial and final foot contacts) from wearable inertial sensors. A rich dataset has been used to validate the method, featuring 71 people with Parkinson's disease (PD) and 67 healthy control subjects. Multiple sensors have been considered, one located on the fifth lumbar vertebrae and two on the ankles. The aims of this study were: (i) to apply deep learning (DL) techniques on wearable sensor data for gait segmentation and quantification in older adults and in people with PD; (ii) to validate the proposed technique for measuring gait against traditional gold standard laboratory reference and a widely used algorithm based on wavelet transforms (WT); (iii) to assess the performance of DL methods in assessing high-level gait characteristics, with focus on stride, stance and swing related features. The results showed a high reliability of the proposed approach, which achieves temporal errors considerably smaller than WT, in particular for the detection of final contacts, with an inter-quartile range below 70 ms in the worst case. This study showes encouraging results, and paves the road for further research, addressing the effectiveness and the generalization of data-driven learning systems for accurate event detection in challenging conditions.


title: Boundary integral formulations of eigenvalue problems for elliptic differential operators with singular interactions and their numerical approximation by boundary element methods

abstract: In this paper the discrete eigenvalues of elliptic second order differential operators in $L^2(\mathbb{R}^n)$, $n \in \mathbb{N}$, with singular $\delta$- and $\delta'$-interactions are studied. We show the self-adjointness of these operators and derive equivalent formulations for the eigenvalue problems involving boundary integral operators. These formulations are suitable for the numerical computations of the discrete eigenvalues and the corresponding eigenfunctions by boundary element methods. We provide convergence results and show numerical examples.

categories: math.SP cs.NA math-ph math.AP math.MP math.NA


title: Parabolic problems in generalized Sobolev spaces

abstract: We consider a general inhomogeneous parabolic initial-boundary value problem for a $2b$-parabolic differential equation given in a finite multidimensional cylinder. We investigate the solvability of this problem in some generalized anisotropic Sobolev spaces. They are parametrized with a pair of positive numbers $s$ and $s/(2b)$ and with a function $\varphi:[1,\infty)\to(0,\infty)$ that varies slowly at infinity. The function parameter $\varphi$ characterizes subordinate regularity of distributions with respect to the power regularity given by the number parameters. We prove that the operator corresponding to this problem is an isomorphism on appropriate pairs of these spaces. As an application, we give a theorem on the local regularity of the generalized solution to the problem. We also obtain sharp sufficient conditions under which chosen generalized derivatives of the solution are continuous on a given set.


title: No-dimensional Tverberg Theorems and Algorithms

abstract: Tverberg's theorem states that for any $k \ge 2$ and any set $P \subset \mathbb{R}^d$ of at least $(d + 1)(k - 1) + 1$ points in $d$ dimensions, we can partition $P$ into $k$ subsets whose convex hulls have a non-empty intersection. The associated search problem of finding the partition lies in the complexity class $\text{CLS} = \text{PPAD} \cap \text{PLS}$, but no hardness results are known. In the colorful Tverberg theorem, the points in $P$ have colors, and under certain conditions, $P$ can be partitioned into colorful sets, in which each color appears exactly once and whose convex hulls intersect. To date, the complexity of the associated search problem is unresolved. Recently, Adiprasito, Barany, and Mustafa gave a no-dimensional Tverberg theorem, in which the convex hulls may intersect in an approximate fashion. This relaxes the requirement on the cardinality of $P$. The argument is constructive, but does not result in a polynomial-time algorithm. We present a deterministic algorithm that finds for any $n$-point set $P \subset \mathbb{R}^d$ and any $k \in \{2, \dots, n\}$ in $O(nd \lceil{\log k}\rceil)$ time a $k$-partition of $P$ such that there is a ball of radius $O\left((k/\sqrt{n})\mathrm{diam(P)}\right)$ that intersects the convex hull of each set. Given that this problem is not known to be solvable exactly in polynomial time, our result provides a remarkably efficient and simple new notion of approximation. Our main contribution is to generalize Sarkaria's method to reduce the Tverberg problem to the Colorful Caratheodory problem (in the simplified tensor product interpretation of Barany and Onn) and to apply it algorithmically. It turns out that this not only leads to an alternative algorithmic proof of a no-dimensional Tverberg theorem, but it also generalizes to other settings such as the colorful variant of the problem.


title: Simulation and Control of a Nonsmooth Cahn-Hilliard Navier-Stokes System

abstract: We are concerned with the simulation and control of a two phase flow model governed by a coupled Cahn-Hilliard Navier-Stokes system involving a nonsmooth energy potential. We establish the existence of optimal solutions and present two distinct approaches to derive suitable stationarity conditions for the bilevel problem, namely C- and strong stationarity. Moreover, we demonstrate the numerical realization of these concepts at the hands of two adaptive solution algorithms relying on a specifically developed goal-oriented error estimator. In addition, we present a model order reduction approach using proper orthogonal decomposition (POD-MOR) in order to replace high-fidelity models by low order surrogates. In particular, we combine POD with space-adapted snapshots and address the challenges which are the consideration of snapshots with different spatial resolutions and the conservation of a solenoidal property.


title: UW-BHI at MEDIQA 2019: An Analysis of Representation Methods for Medical Natural Language Inference

abstract: Recent advances in distributed language modeling have led to large performance increases on a variety of natural language processing (NLP) tasks. However, it is not well understood how these methods may be augmented by knowledge-based approaches. This paper compares the performance and internal representation of an Enhanced Sequential Inference Model (ESIM) between three experimental conditions based on the representation method: Bidirectional Encoder Representations from Transformers (BERT), Embeddings of Semantic Predications (ESP), or Cui2Vec. The methods were evaluated on the Medical Natural Language Inference (MedNLI) subtask of the MEDIQA 2019 shared task. This task relied heavily on semantic understanding and thus served as a suitable evaluation set for the comparison of these representation methods.


title: Approximation in the mean by rational functions II

abstract: For $1\le t < \infty$, a compact subset $K\subset\mathbb C$, and a finite positive measure $\mu$ supported on $K$, $R^t(K, \mu)$ denotes the closure in $L^t(\mu)$ of rational functions with poles off $K$. Conway and Yang (2019) introduced the concept of non-removable boundary $\mathcal F$ and removable set $\mathcal R = K\setminus \mathcal F$ for $R^t(K, \mu)$. We continue the previous work and obtain structural results for $R^t(K, \mu)$. Assume that $S_\mu$, the multiplication by $z$ on $R^t(K, \mu)$, is pure ($R^t(K, \mu)$ does not have $L^t$ summand). Let $H^\infty_{\mathcal R}(A_{\mathcal R})$ be the weak$^*$ closure in $L^\infty (A_{\mathcal R})$ of the functions that are bounded analytic off compact subsets of $\mathcal F$, where $A_{\mathcal R}$ denotes the area measure restricted to $\mathcal R$. $\mathcal R$ is $\gamma$-connected ($\gamma$ denotes analytic capacity) if for any two disjoint open set $G_1$ and $G_2$ with $\mathcal R \subset G_1 \cup G_2 ~\gamma-a.a.$, then $\mathcal R \subset G_1 ~\gamma-a.a.$ or $\mathcal R \subset G_2 ~\gamma-a.a.$. We prove: (1) $R^t(K, \mu)$ contains no non-trivial characterization functions if and only if the removable set $\mathcal R$ is $\gamma$-connected. (2) There is an isometric isomorphism and a weak$^*$ homeomorphism from $R^t(K, \mu)\cap L^\infty(\mu )$ onto $H^\infty_{\mathcal R}(A_{\mathcal R })$.


title: Betti numbers of symmetric shifted ideals

abstract: We introduce a new class of monomial ideals which we call symmetric shifted ideals. Symmetric shifted ideals are fixed by the natural action of the symmetric group and, within the class of monomial ideals fixed by this action, they can be considered as an analogue of stable monomial ideals within the class of monomial ideals. We show that a symmetric shifted ideal has linear quotients and compute its (equivariant) graded Betti numbers. As an application of this result, we obtain several consequences for graded Betti numbers of symbolic powers of defining ideals of star configurations.


title: Early-Time and Late-Time Quantum Chaos

abstract: We show the relation between the Heisenberg averaging of regularized 2-point out-of-time ordered correlation function and the 2-point spectral form factor in bosonic quantum mechanics. The generalization to all even-point is also discussed. We also do the direct extension from the bosonic quantum mechanics to the non-interacting scalar field theory. Finally, we find that the coherent state and large-$N$ approaches are useful in the late-time study. We find that the computation of the coherent state can be simplified by the Heisenberg averaging. Therefore, this provides a simplified way to probe the late-time quantum chaos through a coherent state. The large-$N$ result is also comparable to the $N=3$ numerical result in the large-$N$ quantum mechanics. This can justify that large-$N$ technique in bosonic quantum mechanics can probe the late time, not the early time. Because the quantitative behavior of large-$N$ can be captured from the $N=3$ numerical result, the realization in experiments should be possible.


title: Fundamental Group Schemes of Hilbert Scheme of $n$ Points on a Smooth Projective Surface

abstract: Let $k$ be an algebraically closed field of characteristic $p > 3$. Let $X$ be an irreducible smooth projective surface over $k$. Fix an integer $n \geq 1$ and let ${\mathcal{H}{\it ilb}}_X^n$ be the Hilbert scheme parameterizing effective $0$-cycles of length $n$ on $X$. The aim of the present article is to find the $S$-fundamental group scheme and Nori's fundamental group scheme of the Hilbert scheme $\mathcal{H}{\it ilb}_X^n$.


title: An application of survival analysis to disruption prediction via Random Forests

abstract: One of the most pressing challenges facing the fusion community is adequately mitigating or, even better, avoiding disruptions of tokamak plasmas. However, before this can be done, disruptions must first be predicted with sufficient warning time to actuate a response. The established field of survival analysis provides a convenient statistical framework for time-to-event (i.e. time-to-disruption) studies. This paper demonstrates the integration of an existing disruption prediction machine learning algorithm with the Kaplan-Meier estimator of survival probability. Specifically discussed are the implied warning times from binary classification of disruption databases and the interpretation of output signals from Random Forest algorithms trained and tested on these databases. This survival analysis approach is applied to both smooth and noisy test data to highlight important features of the survival and hazard functions. In addition, this method is applied to three Alcator C-Mod plasma discharges and compared to a threshold-based scheme for triggering alarms. In one case, both techniques successfully predict the disruption; although, in another, neither warns of the impending disruption with enough time to mitigate. For the final discharge, the survival analysis approach could avoid the false alarm triggered by the threshold method. Limitations of this analysis and opportunities for future work are also presented.


title: Evolution of the Informational Complexity of Contemporary Western Music

abstract: We measure the complexity of songs in the Million Song Dataset (MSD) in terms of pitch, timbre, loudness, and rhythm to investigate their evolution from 1960 to 2010. By comparing the Billboard Hot 100 with random samples, we find that the complexity of popular songs tends to be more narrowly distributed around the mean, supporting the idea of an inverted U-shaped relationship between complexity and hedonistic value. We then examine the temporal evolution of complexity, reporting consistent changes across decades, such as a decrease in average loudness complexity since the 1960s, and an increase in timbre complexity overall but not for popular songs. We also show, in contrast to claims that popular songs sound more alike over time, that they are not more similar than they were 50 years ago in terms of pitch or rhythm, although similarity in timbre shows distinctive patterns across eras and similarity in loudness has been increasing. Finally, we show that musical genres can be differentiated by their distinctive complexity profiles.


title: Quench Dynamics in 1D Optomechanical Arrays

abstract: Non-equilibrium dynamics induced by rapid changes of external parameters is relevant for a wide range of scenarios across many domains of physics. For waves in spatially periodic systems, quenches will alter the bandstructure and generate new excitations. In the case of topological bandstructures, defect modes at boundaries can be generated or destroyed when quenching through a topological phase transition. Here, we demonstrate that optomechanical arrays are a promising platform for studying such dynamics, as their bandstructure can be tuned temporally by a control laser. We study the creation of nonequilibrium optical and mechanical excitations in 1D arrays, including a bosonic version of the Su-Schrieffer-Heeger model. These ideas can be transferred to other systems such as driven nonlinear cavity arrays.


title: An Attention Mechanism for Musical Instrument Recognition

abstract: While the automatic recognition of musical instruments has seen significant progress, the task is still considered hard for music featuring multiple instruments as opposed to single instrument recordings. Datasets for polyphonic instrument recognition can be categorized into roughly two categories. Some, such as MedleyDB, have strong per-frame instrument activity annotations but are usually small in size. Other, larger datasets such as OpenMIC only have weak labels, i.e., instrument presence or absence is annotated only for long snippets of a song. We explore an attention mechanism for handling weakly labeled data for multi-label instrument recognition. Attention has been found to perform well for other tasks with weakly labeled data. We compare the proposed attention model to multiple models which include a baseline binary relevance random forest, recurrent neural network, and fully connected neural networks. Our results show that incorporating attention leads to an overall improvement in classification accuracy metrics across all 20 instruments in the OpenMIC dataset. We find that attention enables models to focus on (or `attend to') specific time segments in the audio relevant to each instrument label leading to interpretable results.


title: Better Sample -- Random Subset Sum in $2^{0.255n}$ and its Impact on Decoding Random Linear Codes

abstract: We propose a new heuristic algorithm for solving random subset sum instances $a_1, \ldots, a_n, t \in \mathbb{Z}_{2^n}$, which play a crucial role in cryptographic constructions. Our algorithm is search tree-based and solves the instances in a divide-and-conquer method using the representation method. From a high level perspective, our algorithm is similar to the algorithm of Howgrave-Graham-Joux (HGJ) and Becker-Coron-Joux (BCJ), but instead of enumerating the initial lists we sample candidate solutions. So whereas HGJ and BCJ are based on combinatorics, our analysis is stochastic. Our sampling technique introduces variance that increases the amount of representations and gives our algorithm more optimization flexibility. This results in the remarkable and natural property that we improve with increasing search tree depth. Whereas BCJ achieves the currently best known (heuristic) run time $2^{0.291n}$ for random subset sum, we improve (heuristically) down to $2^{0.255n}$ using a search tree of depth at least $13$. We also apply our subset algorithm to the decoding of random binary linear codes, where we improve the best known run time of the Becker-Joux-May-Meurer algorithm from $2^{0.048n}$ in the half distance decoding setting down to $2^{0.042n}$.


title: Particle Acceleration and Heating in Regions of Magnetic Flux Emergence

abstract: The interaction between emerging and pre-existing magnetic fields in the solar atmosphere can trigger several dynamic phenomena, such as eruptions and jets. A key element during this interaction is the formation of large scale current sheets and, eventually, their fragmentation that leads to the creation of a strongly turbulent environment. In this paper, we study the kinetic aspects of the interaction (reconnection) between emerging and ambient magnetic fields. We show that the statistical properties of the spontaneously fragmented and fractal electric fields are responsible for the efficient heating and acceleration of charged particles, which form a power law tail at high energies on sub-second time scales. A fraction of the energized particles escapes from the acceleration volume, with a super-hot component with temperature close to $150\,$MK, and with a power law high energy tail with index between -2 and -3. We estimate the transport coefficients in energy space from the dynamics of the charged particles inside the fragmented and fractal electric fields, and the solution of a fractional transport equation, as appropriate for a strongly turbulent plasma, agrees with the test particle simulations. We also show that the acceleration mechanism is not related to Fermi acceleration, and the Fokker Planck equation is inconsistent and not adequate as a transport model. Finally, we address the problem of correlations between spatial transport and transport in energy space. Our results confirm the observations reported for high energy particles (hard X-rays, type III bursts and solar energetic particles) during the emission of solar jets.


title: Quantum jumps and attractors of the Maxwell-Schr\"odinger equations

abstract: Our goal is the discussion of the problem of mathematical interpretation of basic postulates (or `principles') of Quantum Mechanics: transitions to quantum stationary orbits, the wave-particle duality, and the probabilistic interpretation, in the context of semiclassical self-consistent Maxwell--Schr\"odinger equations. We discuss possible relations of these postulates to the theory of attractors of Hamiltonian nonlinear PDEs and to a new general mathematical conjecture on global attractors of G-invariant nonlinear Hamiltonian partial differential equations with a Lie symmetry group G. This conjecture is inspired by our results on global attractors of nonlinear Hamiltonian PDEs obtained since 1990 for a list of model equations with three basic symmetry groups: the trivial group, the group of translations, and the unitary group U(1). We present sketchy these results.


title: Deep Learning for Spacecraft Pose Estimation from Photorealistic Rendering

abstract: On-orbit proximity operations in space rendezvous, docking and debris removal require precise and robust 6D pose estimation under a wide range of lighting conditions and against highly textured background, i.e., the Earth. This paper investigates leveraging deep learning and photorealistic rendering for monocular pose estimation of known uncooperative spacecrafts. We first present a simulator built on Unreal Engine 4, named URSO, to generate labeled images of spacecrafts orbiting the Earth, which can be used to train and evaluate neural networks. Secondly, we propose a deep learning framework for pose estimation based on orientation soft classification, which allows modelling orientation ambiguity as a mixture of Gaussians. This framework was evaluated both on URSO datasets and the ESA pose estimation challenge. In this competition, our best model achieved 3rd place on the synthetic test set and 2nd place on the real test set. Moreover, our results show the impact of several architectural and training aspects, and we demonstrate qualitatively how models learned on URSO datasets can perform on real images from space.


title: On the 3-D Placement of Airborne Base Stations Using Tethered UAVs

abstract: One of the main challenges slowing the deployment of airborne base stations (BSs) using unmanned aerial vehicles (UAVs) is the limited on-board energy and flight time. One potential solution to such problem, is to provide the UAV with power supply through a tether that connects the UAV to the ground. In this paper, we study the optimal placement of tethered UAVs (TUAVs) to minimize the average path-loss between the TUAV and a receiver located on the ground. Given that the tether has a maximum length, and the launching point of the TUAV (the starting point of the tether) is placed on a rooftop, the TUAV is only allowed to hover within a specific hovering region. Beside the maximum tether length, this hovering region also depends on the heights of the buildings surrounding the rooftop, which requires the inclination angle of the tether not to be below a given minimum value, in order to avoid tangling and ensure safety. We first formulate the optimization problem for such setup and provide some useful insights on its solution. Next, we derive upper and lower bounds for the optimal values of the tether length and inclination angle. We also propose a suboptimal closed-form solution for the tether length and its inclination angle that is based on maximizing the line-of-sight probability. Finally, we derive the probability distribution of the minimum inclination angle of the tether length. We show that its mean value varies depending on the environment from 10 degrees in suburban environments to 31 degrees in high rise urban environments. Our numerical results show that the derived upper and lower bounds on the optimal values of the tether length and inclination angle lead to tight suboptimal values of the average path-loss that are only 0-3 dBs above the minimum value.


title: Shock-accelerated cosmic rays and streaming instability in the adaptive mesh refinement code Ramses

abstract: Cosmic rays (CRs) are thought to play a dynamically important role in several key aspects of galaxy evolution, including the structure of the interstellar medium, the formation of galactic winds, and the non-thermal pressure support of halos. We introduce a numerical model solving for the CR streaming instability and acceleration of CRs at shocks with a fluid approach in the adaptive mesh refinement code ramses. CR streaming is solved with a diffusion approach and its anisotropic nature is naturally captured. We introduce a shock finder for the ramses code that automatically detects shock discontinuities in the flow. Shocks are the loci for CR injection, and their efficiency of CR acceleration is made dependent on the upstream magnetic obliquity according to the diffuse shock acceleration mechanism. We show that the shock finder accurately captures shock locations and estimates the shock Mach number for several problems. The obliquity-dependent injection of CRs in the Sedov solution leads to situations where the supernova bubble exhibits large polar caps (homogeneous background magnetic field), or a patchy structure of the CR distribution (inhomogeneous background magnetic field). Finally, we combine both accelerated CRs with streaming in a simple turbulent interstellar medium box, and show that the presence of CRs significantly modifies the structure of the gas.


title: Role of zero point energy in promoting ice formation in a spherical drop of water

abstract: We demonstrate that the Lifshitz interaction energy (excluding the self-energies of the inner and outer spherical regions) for three concentric spherical dielectric media can be evaluated easily using the immense computation power in recent processors relative to those of a few decades ago. As a prototype, we compute the Lifshitz interaction energy for a spherical shell of water immersed in water vapor of infinite extent while enclosing a spherical ball of ice inside the shell, such that two concentric spherical interfaces are formed: one between solid ice and liquid water and the other between liquid water and gaseous vapor. We evaluate the Lifshitz interaction energy for the above configuration at the triple point of water when the solid, liquid, and gaseous states of water coexist, and, thus, extend the analysis of Elbaum and Schick in Phys. Rev. Lett. 66 (1991) 1713 to spherical configurations. We find that, when the Lifshitz energy contributes dominantly to the total energy of this system, which is often the case when electrostatic interactions are absent, a drop of water surrounded by vapor of infinite extent is not stable at the triple point. This instability, that is a manifestation of the quantum fluctuations in the medium, will promote formation of ice in water, which will then grow in size indefinitely. This is a consequence of the finding here that the Lifshitz energy is minimized for large (micrometer size) radius of the ice ball and small (nanometer size) thickness of the water shell surrounding the ice. These results might be relevant to the formation of hail in thunderclouds. These results are tentative in that the self-energies are omitted; surface tension and nucleation energy are not considered.


title: Interactive Verifiable Polynomial Evaluation

abstract: Cloud computing platforms have created the possibility for computationally limited users to delegate demanding tasks to strong but untrusted servers. Verifiable computing algorithms help build trust in such interactions by enabling the server to provide a proof of correctness of his results which the user can check very efficiently. In this paper, we present a doubly-efficient interactive algorithm for verifiable polynomial evaluation. Unlike the mainstream literature on verifiable computing, the soundness of our algorithm is information-theoretic and cannot be broken by a computationally unbounded server. By relying on basic properties of error correcting codes, our algorithm enforces a dishonest server to provide false results to problems which become progressively easier to verify. After roughly $\log d$ rounds, the user can verify the response of the server against a look-up table that has been pre-computed during an initialization phase. For a polynomial of degree $d$, we achieve a user complexity of $O(d^{\epsilon})$, a server complexity of $O(d^{1+\epsilon})$, a round complexity of $O(\log d)$ and an initialization complexity of $O(d^{1+\epsilon})$.


title: A Bayesian Approach for Analyzing Data on the Stiefel Manifold

abstract: Directional data emerges in a wide array of applications, ranging from atmospheric sciences to medical imaging. Modeling such data, however, poses unique challenges by virtue of their being constrained to non-Euclidean spaces like manifolds. Here, we present a unified Bayesian framework for inference on the Stiefel manifold using the Matrix Langevin distribution. Specifically, we propose a novel family of conjugate priors and establish a number of theoretical properties relevant to statistical inference. %Importantly, these include the propriety of these priors and concentration characterization. Conjugacy enables the translation of these properties to their corresponding posteriors, which we exploit to develop the posterior inference scheme. For the implementation of the posterior computation, including the posterior sampling, we adopt a novel computational procedure for evaluating the hypergeometric function of matrix arguments that appears as normalization constants in the relevant densities.