title: On the general dual Orlicz-Minkowski problem

abstract: For $K\subseteq \mathbb{R}^n$ a convex body with the origin $o$ in its interior, and $\phi:\mathbb{R}^n\setminus\{o\}\rightarrow(0, \infty)$ a continuous function, define the general dual ($L_{\phi})$ Orlicz quermassintegral of $K$ by $$\mathcal{V}_\phi(K)=\int_{\mathbb{R}^n \setminus K} \phi(x)\,dx.$$ Under certain conditions on $\phi$, we prove a variational formula for the general dual ($L_{\phi})$ Orlicz quermassintegral, which motivates the definition of $\widetilde{C}_{\phi,\mathcal{V}}(K, \cdot)$, the general dual ($L_{\phi})$ Orlicz curvature measure of $K$. We pose the following general dual Orlicz-Minkowski problem: {\it Given a nonzero finite Borel measure $\mu$ defined on $S^{n-1}$ and a continuous function $\phi: \mathbb{R}^n\setminus\{o\}\rightarrow (0, \infty)$, can one find a constant $\tau>0$ and a convex body $K$ (ideally, containing $o$ in its interior), such that,} $$\mu=\tau\widetilde{C}_{\phi,\mathcal{V}}(K,\cdot)? $$ Based on the method of Lagrange multipliers and the established variational formula for the general dual ($L_{\phi})$ Orlicz quermassintegral, a solution to the general dual Orlicz-Minkowski problem is provided. In some special cases, the uniqueness of solutions is proved and the solution for $\mu$ being a discrete measure is characterized.


title: A rank-based Cram\'er-von-Mises-type test for two samples

abstract: We study a rank based univariate two-sample distribution-free test. The test statistic is the difference between the average of between-group rank distances and the average of within-group rank distances. This test statistic is closely related to the two-sample Cram\'er-von Mises criterion. They are different empirical versions of a same quantity for testing the equality of two population distributions. Although they may be different for finite samples, they share the same expected value, variance and asymptotic properties. The advantage of the new rank based test over the classical one is its ease to generalize to the multivariate case. Rather than using the empirical process approach, we provide a different easier proof, bringing in a different perspective and insight. In particular, we apply the H\'ajek projection and orthogonal decomposition technique in deriving the asymptotics of the proposed rank based statistic. A numerical study compares power performance of the rank formulation test with other commonly-used nonparametric tests and recommendations on those tests are provided. Lastly, we propose a multivariate extension of the test based on the spatial rank.


title: Explicit equations of a fake projective plane

abstract: Fake projective planes are smooth complex surfaces of general type with Betti numbers equal to those of the usual projective plane. They come in complex conjugate pairs and have been classified as quotients of the two-dimensional ball by explicitly written arithmetic subgroups. In this paper we find equations of a projective model of a conjugate pair of fake projective planes by studying the geometry of the quotient of such surface by an order seven automorphism.


title: New Physics Searches Using Precision Spectroscopy

abstract: The exceptional precision attainable using modern spectroscopic techniques provides a promising avenue to search for signatures of physics beyond the Standard Model in tiny shifts of the energy levels of atoms and molecules. We briefly review three categories of new-physics searches based in precision measurements: tests of QED using measurements of the anomalous magnetic moment of the electron and the value of the fine-structure constant, searches for time variation of the fundamental constants, and searches for a permanent electric dipole moment of an electron or atomic nucleus.


title: A Pieri-type formula and a factorization formula for sums of $K$-$k$-Schur functions

abstract: We give a Pieri-type formula for the sum of $K$-$k$-Schur functions $\sum_{\mu\le\lambda} g^{(k)}_{\mu}$ over a principal order ideal of the poset of $k$-bounded partitions under the strong Bruhat order, which sum we denote by $\widetilde{g}^{(k)}_{\lambda}$. As an application of this, we also give a $k$-rectangle factorization formula $\widetilde{g}^{(k)}_{R_t\cup\lambda}=\widetilde{g}^{(k)}_{R_t} \widetilde{g}^{(k)}_{\lambda}$ where $R_t=(t^{k+1-t})$, analogous to that of $k$-Schur functions $s^{(k)}_{R_t\cup\lambda}=s^{(k)}_{R_t}s^{(k)}_{\lambda}$.


title: Uniqueness of power of a meromorphic function with its differential polynomial

abstract: In this paper, taking the question of Zhang and L\"{u} into the background, we present one theorem which will improve and extend some recent results related to the Br\"{u}ck Conjecture.


title: Large thermoelectric response in a diluted ferroelectric system: Ba0.7Eu0.3Ti1-xNbxO3

abstract: We investigated the electrical conductivity, thermal conductivity and thermopower as a function of Nb content (x) in Ba0.7Eu0.3Ti1-xNbxO3 (x = 0.001- 0.10) in the temperature range T = 400-2 K. The substitution of Nb destabilizes the ferroelectric insulating ground state of Ba0.7Eu0.3TiO3 and transforms into a paramagnetic metal for x = 0.1. Thermopower is negative in the entire composition range (S = -613 microVolt/K at 400 K for x = 0.001) and its magnitude decreases with increasing Nb content which suggests doping of electrons into empty Ti-3d(t2g) conduction band. In this series, the dimensionless figure of merit (ZT) increases with temperature for all the compositions and the x = 0.03 composition exhibits the maximum ZT (= 0.12 at 400 K). The enhanced value of ZT is primarily due to the low thermal conductivity of samples in this series (~ 0.7 to 1 W/(m.K) at 400 K) compared to other potential high temperature n-type thermoelectric oxides such as carrier doped SrTiO3 and CaMnO3. The low thermal conductivity in our compounds most likely arises from heavy Eu2+ ion and lattice disorder introduced by Nb5+ which scatter phonons effectively.


title: Sequence-to-Sequence Prediction of Vehicle Trajectory via LSTM Encoder-Decoder Architecture

abstract: In this paper, we propose a deep learning based vehicle trajectory prediction technique which can generate the future trajectory sequence of surrounding vehicles in real time. We employ the encoder-decoder architecture which analyzes the pattern underlying in the past trajectory using the long short-term memory (LSTM) based encoder and generates the future trajectory sequence using the LSTM based decoder. This structure produces the $K$ most likely trajectory candidates over occupancy grid map by employing the beam search technique which keeps the $K$ locally best candidates from the decoder output. The experiments conducted on highway traffic scenarios show that the prediction accuracy of the proposed method is significantly higher than the conventional trajectory prediction techniques.


title: Level-zero van der Kallen modules and specialization of nonsymmetric Macdonald polynomials at $t = \infty$

abstract: Let $\lambda \in P^{+}$ be a level-zero dominant integral weight, and $w$ an arbitrary coset representative of minimal length for the cosets in $W/W_{\lambda}$, where $W_{\lambda}$ is the stabilizer of $\lambda$ in a finite Weyl group $W$. In this paper, we give a module $\mathbb{K}_{w}(\lambda)$ over the negative part of a quantum affine algebra whose graded character is identical to the specialization at $t = \infty$ of the nonsymmetric Macdonald polynomial $E_{w \lambda}(q,\,t)$ multiplied by a certain explicit finite product of rational functions of $q$ of the form $(1 - q^{-r})^{-1}$ for a positive integer $r$. This module $\mathbb{K}_{w}(\lambda)$ (called a level-zero van der Kallen module) is defined to be the quotient module of the level-zero Demazure module $V_{w}^{-}(\lambda)$ by the sum of the submodules $V_{z}^{-}(\lambda)$ for all those coset representatives $z$ of minimal length for the cosets in $W/W_{\lambda}$ such that $z > w$ in the Bruhat order $<$ on $W$.


title: Graphical Models for Non-Negative Data Using Generalized Score Matching

abstract: A common challenge in estimating parameters of probability density functions is the intractability of the normalizing constant. While in such cases maximum likelihood estimation may be implemented using numerical integration, the approach becomes computationally intensive. In contrast, the score matching method of Hyv\"arinen (2005) avoids direct calculation of the normalizing constant and yields closed-form estimates for exponential families of continuous distributions over $\mathbb{R}^m$. Hyv\"arinen (2007) extended the approach to distributions supported on the non-negative orthant $\mathbb{R}_+^m$. In this paper, we give a generalized form of score matching for non-negative data that improves estimation efficiency. We also generalize the regularized score matching method of Lin et al. (2016) for non-negative Gaussian graphical models, with improved theoretical guarantees.


title: The Mu2e Calorimeter Final Technical Design Report

abstract: Since the first version of the Mu2e TDR released at the beginning of 2015, the Mu2e Calorimeter system has undergone a long list of changes to arrive to its final design. These changes were primarily caused by two reasons: (i) the technology choice between the TDR proposed solution of BaF2 crystals readout with solar blind Avalanche Photodiodes (APDs) and the backup option of CsI crystals readout with Silicon Photomultipliers (SiPM) has been completed and (ii) the channels numbering, the mechanical system and the readout electronics were substantially modified while proceeding with engineering towards the final project. This document updates the description of the calorimeter system adding the most recent engineering drawings and tecnical progresses.


title: Stability Theorems for Group Actions on Uniform Spaces

abstract: We extend the notions of topological stability, shadowing and persistence from homeomorphisms to finitely generated group actions on uniform spaces and prove that an expansive action with either shadowing or persistence is topologically stable. Using the concept of null set of a Borel measure $\mu$, we introduce the notions of $\mu$-expansivity, $\mu$-topological stability, $\mu$-shadowing and $\mu$-persistence for finitely generated group actions on uniform spaces and show that a $\mu$-expansive action with either $\mu$-shadowing or $\mu$-persistence is $\mu$-topologically stable.


title: On an infinite limit of BGG categories O

abstract: We study a version of the BGG category O for Dynkin Borel subalgebras of root-reductive Lie algebras g, such as gl(\infty). We prove results about extension fullness and compute the higher extensions of simple modules by Verma modules. In addition, we show that our category O is Ringel self-dual and initiate the study of Koszul duality. An important tool in obtaining these results is an equivalence we establish between appropriate Serre subquotients of category O for g and category O for finite dimensional reductive subalgebras of g.


title: CGC/saturation approach: an impact-parameter dependent model for diffraction production in DIS

abstract: In the paper we discussed the evolution equations for diffractive production in the framework of CGC/saturation approach, and found the analytical solutions for several kinematic regions. The most impressive features of these solutions are, that diffractive production does not exibit geometric scaling behaviour i.e. being a function of one variable. Based on these solutions, we suggest an impact parameter dependent saturation model, which is suitable for describing diffraction production both deep in the saturation region, and in the vicinity of the saturation scale. Using the model we attempted to fit the combined data on diffraction production from H1 and ZEUS collaborations. We found that we are able describe both $x_\pom$ and $\beta$ dependence, as well as $Q$ behavior of the measured cross sections. In spite of the sufficiently large $\chi^2/d.o.f.$ we believe that our description provides an initial impetus to find a fit of the experimental data, based on the solution of the CGC/saturation equation, rather than on describing the diffraction system in simplistic manner, assuming that only quark-antiquark pair and one extra gluons, are produced.


title: Resolving Two Conjectures on Staircase Encodings and Boundary Grids of $132$ and $123$-avoiding permutations

abstract: This paper analyzes relations between pattern avoidance of certain permutations and graphs on staircase grids and boundary grids, and proves two conjectures posed by Bean, Tannock, and Ulfarsson (2015). More specifically, this paper enumerates a certain family of staircase encodings and proves that the downcore graph, a certain graph established on the boundary grid, is pure if and only if the permutation corresponding to the boundary grid avoids the classical patterns 123 and 2143.


title: Design and status of the Mu2e crystal calorimeter

abstract: The Mu2e experiment at Fermilab searches for the charged-lepton flavour violating (CLFV) conversion of a negative muon into an electron in the field of an aluminum nucleus, with a distinctive signature of a mono-energetic electron of energy slightly below the muon rest mass (104.967 MeV). The Mu2e goal is to improve by four orders of magnitude the search sensitivity with respect to the previous experiments. Any observation of a CLFV signal will be a clear indication of new physics. The Mu2e detector is composed of a tracker, an electro- magnetic calorimeter and an external veto for cosmic rays surrounding the solenoid. The calorimeter plays an important role in providing particle identification capabilities, a fast online trigger filter, a seed for track reconstruction while working in vacuum, in the presence of 1 T axial magnetic field and in an harsh radiation environment. The calorimeter requirements are to provide a large acceptance for 100 MeV electrons and reach at these energies: (a) a time resolution better than 0.5 ns; (b) an energy resolution < 10% and (c) a position resolution of 1 cm. The calorimeter design consists of two disks, each one made of 674 undoped CsI crystals read by two large area arrays of UV-extended SiPMs. We report here the construction and test of the Module-0 prototype. The Module-0 has been exposed to an electron beam in the energy range around 100 MeV at the Beam Test Facility in Frascati. Preliminary results of timing and energy resolution at normal incidence are shown. A discussion of the technical aspects of the calorimeter engineering is also reported in this paper.


title: Optimal stopping, randomized stopping and singular control with partial information flow

abstract: The purpose of this paper is two-fold: We extend the well-known relation between optimal stopping and randomized stopping of a given stochastic process to a situation where the available information flow is a filtration with no a priori assumed relation to the filtration of the process. We call these problems optimal stopping and randomized stopping with general information. Following an idea of Krylov [K] we introduce a special singular stochastic control problem with general information and show that this is also equivalent to the partial information optimal stopping and randomized stopping problems. Then we show that the solution of this singular control problem can be expressed in terms of partial information variational inequalities.


title: Quantitative Predictions in Quantum Decision Theory

abstract: Quantum Decision Theory, advanced earlier by the authors, and illustrated for lotteries with gains, is generalized to the games containing lotteries with gains as well as losses. The mathematical structure of the approach is based on the theory of quantum measurements, which makes this approach relevant both for the description of decision making of humans and the creation of artificial quantum intelligence. General rules are formulated allowing for the explicit calculation of quantum probabilities representing the fraction of decision makers preferring the considered prospects. This provides a method to quantitatively predict decision-maker choices, including the cases of games with high uncertainty for which the classical expected utility theory fails. The approach is applied to experimental results obtained on a set of lottery gambles with gains and losses. Our predictions, involving no fitting parameters, are in very good agreement with experimental data. The use of quantum decision making in game theory is described. A principal scheme of creating quantum artificial intelligence is suggested.


title: Isotropic LQC and LQC--inspired Models with a massless scalar field as Generalised Brans--Dicke theories

abstract: We explore whether generalised Brans -- Dicke theories, which have a scalar field $\Phi$ and a function $\omega(\Phi)$, can be the effective actions leading to the effective equations of motion of the LQC and the LQC--inspired models, which have a massless scalar field $\sigma$ and a function $f(m) \;$. We find that this is possible for isotropic cosmology. We relate the pairs $(\sigma, f)$ and $(\Phi, \omega)$ and, using examples, illustrate these relations. We find that near the bounce of the LQC evolutions for which $f(m) = sin \; m$, the corresponding field $\Phi \to 0$ and the function $\omega(\Phi) \propto \Phi^2 \;$. We also find that the class of generalised Brans -- Dicke theories, which we had found earlier to lead to non singular isotropic evolutions, may be written as an LQC--inspired model. The relations found here in the isotropic cases do not apply to the anisotropic cases, which perhaps require more general effective actions.


title: Spatial modelling with R-INLA: A review

abstract: Coming up with Bayesian models for spatial data is easy, but performing inference with them can be challenging. Writing fast inference code for a complex spatial model with realistically-sized datasets from scratch is time-consuming, and if changes are made to the model, there is little guarantee that the code performs well. The key advantages of R-INLA are the ease with which complex models can be created and modified, without the need to write complex code, and the speed at which inference can be done even for spatial problems with hundreds of thousands of observations. R-INLA handles latent Gaussian models, where fixed effects, structured and unstructured Gaussian random effects are combined linearly in a linear predictor, and the elements of the linear predictor are observed through one or more likelihoods. The structured random effects can be both standard areal model such as the Besag and the BYM models, and geostatistical models from a subset of the Mat\'ern Gaussian random fields. In this review, we discuss the large success of spatial modelling with R-INLA and the types of spatial models that can be fitted, we give an overview of recent developments for areal models, and we give an overview of the stochastic partial differential equation (SPDE) approach and some of the ways it can be extended beyond the assumptions of isotropy and separability. In particular, we describe how slight changes to the SPDE approach leads to straight-forward approaches for non-stationary spatial models and non-separable space-time models.


title: Mixing angle and decay constants of $J^P=1^+$ heavy-light mesons

abstract: The mass spectra, mixing angle and decay constants of the $J^P=1^+$ heavy-light mesons are systematically studied within the framework of the Bethe-Salpeter equation (BSE). The full $1^+$ Salpeter wave function is given for the first time. The mixing between the $1^{+-}$ and $1^{++}$ in the $1^+$ heavy-light systems are automatically determined by the dynamics in the equation without any man-made mixing. The results indicate that in a rigorous study there exists the phenomenon of mixing angle inversion or mass inversion within $1^{+}$ heavy-light doublet, which is sensitive to the $s$-quark mass for the charmed mesons and $u$- or $d$-quark masses for the bottomed mesons. This inversion phenomenon can answer the question of why we have confused mixing angles in the literature and partly explain the lower mass of $D_{s1}(2460)$ compared to that of $D_{s1}(2536)$. The decay constants are also presented and can behave as a good quantity to distinguish the $1^+$ doublet in heavy-light mesons. This study indicates that the light-quark mass may play an important role in deciding the mass order, mixing angle, and decay constant relation between the $\ket{j_l=\frac{3}{2}}$ and $\ket{j_l=\frac{1}{2}}$ heavy-light mesons.


title: Dynamic spin injection into a quantum well coupled to a spin-split bound state

abstract: We present a theoretical analysis of dynamic spin injection due to spin-dependent tunneling between a quantum well (QW) and a bound state split in spin projection due to an exchange interaction or external magnetic field. We focus on the impact of Coulomb correlations at the bound state on spin polarization and sheet density kinetics of the charge carriers in the QW. The theoretical approach is based on kinetic equations for the electron occupation numbers taking into account high order correlation functions for the bound state electrons. It is shown that the on-site Coulomb repulsion leads to an enhanced dynamic spin polarization of the electrons in the QW and a delay in the carriers tunneling into the bound state. The interplay of these two effects leads to non-trivial dependence of the spin polarization degree, which can be probed experimentally using time-resolved photoluminescence experiments. It is demonstrated that the influence of the Coulomb interactions can be controlled by adjusting the relaxation rates. These findings open a new way of studying the Hubbard-like electron interactions experimentally.


title: On the well-posedness of a multiscale mathematical model for Lithium-ion batteries

abstract: We consider the mathematical treatment of a system of nonlinear partial differential equations based on a model, proposed in 1972 by J. Newman, in which the coupling between the Lithium concentration, the phase potentials and temperature in the electrodes and the electrolyte of a Lithium battery cell is considered. After introducing some functional spaces well-adapted to our framework we obtain some rigorous results showing the well-posedness of the system, first for some short time and then, by considering some hypothesis on the nonlinearities, globally in time. As far as we know, this is the first result in the literature proving existence in time of the full Newman model, which follows previous results by the third author in 2016 regarding a simplified case.


title: K2SUPERSTAMP: The release of calibrated mosaics for the {\em Kepler/K2} Mission

abstract: We describe the release of a new High Level Science Product (HLSP) available at the MAST archive. The HLSP, called K2Superstamp, consists of a series of FITS images for four open star clusters observed by the K2 Mission using so-called "superstamp" pixel masks: M35, the $\sim$150 Myr old open cluster observed during K2 Campaign 0, M67, the solar-age, solar-metallicity benchmark cluster observed during Campaign 5, Ruprecht 147, the $\sim$3 Gyr-old open cluster observed during Campaign 7, and the Lagoon Nebula (M8/NGC 6530), the high-mass star-forming region observed during Campaign 9. While the data for these regions have long been served on MAST, until now they were only available as a disconnected set of smaller Target Pixel Files (TPFs) because the spacecraft stored these observations in small chunks. As a result, these regions have hitherto been ignored by many lightcurve and planet search pipelines. With this new release, we have stitched these TPFs together into spatially contiguous FITS images (one per cadence) to make their scientific analysis easier. In addition, each image has been fit with an accurate WCS solution so that you may locate any object of interest via its right ascension and declination. We describe here the process of stitching and astrometric calibration.


title: Stochastic Chebyshev Gradient Descent for Spectral Optimization

abstract: A large class of machine learning techniques requires the solution of optimization problems involving spectral functions of parametric matrices, e.g. log-determinant and nuclear norm. Unfortunately, computing the gradient of a spectral function is generally of cubic complexity, as such gradient descent methods are rather expensive for optimizing objectives involving the spectral function. Thus, one naturally turns to stochastic gradient methods in hope that they will provide a way to reduce or altogether avoid the computation of full gradients. However, here a new challenge appears: there is no straightforward way to compute unbiased stochastic gradients for spectral functions. In this paper, we develop unbiased stochastic gradients for spectral-sums, an important subclass of spectral functions. Our unbiased stochastic gradients are based on combining randomized trace estimators with stochastic truncation of the Chebyshev expansions. A careful design of the truncation distribution allows us to offer distributions that are variance-optimal, which is crucial for fast and stable convergence of stochastic gradient methods. We further leverage our proposed stochastic gradients to devise stochastic methods for objective functions involving spectral-sums, and rigorously analyze their convergence rate. The utility of our methods is demonstrated in numerical experiments.


title: Lamellar ordering, droplet formation and phase inversion in exotic active emulsions

abstract: We study numerically the behaviour of a mixture of a passive isotropic fluid and an active polar gel, in the presence of a surfactant favouring emulsification. Focussing on parameters for which the underlying free energy favours the lamellar phase in the passive limit, we show that the interplay between nonequilibrium and thermodynamic forces creates a range of multifarious exotic emulsions. When the active component is contractile (e.g., an actomyosin solution), moderate activity enhances the efficiency of lamellar ordering, whereas strong activity favours the creation of passive droplets within an active matrix. For extensile activity (occurring, e.g., in microtubule-motor suspensions), instead, we observe an emulsion of spontaneously rotating droplets of different size. By tuning the overall composition, we can create high internal phase emulsions, which undergo sudden phase inversion when activity is switched off. Therefore, we find that activity provides a single control parameter to design composite materials with a strikingly rich range of morphologies.


title: Convergence of Online Mirror Descent

abstract: In this paper we consider online mirror descent (OMD) algorithms, a class of scalable online learning algorithms exploiting data geometric structures through mirror maps. Necessary and sufficient conditions are presented in terms of the step size sequence $\{\eta_t\}_{t}$ for the convergence of an OMD algorithm with respect to the expected Bregman distance induced by the mirror map. The condition is $\lim_{t\to\infty}\eta_t=0, \sum_{t=1}^{\infty}\eta_t=\infty$ in the case of positive variances. It is reduced to $\sum_{t=1}^{\infty}\eta_t=\infty$ in the case of zero variances for which the linear convergence may be achieved by taking a constant step size sequence. A sufficient condition on the almost sure convergence is also given. We establish tight error bounds under mild conditions on the mirror map, the loss function, and the regularizer. Our results are achieved by some novel analysis on the one-step progress of the OMD algorithm using smoothness and strong convexity of the mirror map and the loss function.


title: Comments on $^"$Cosmic evolution in Brans-Dicke chameleon cosmology$^"$

abstract: The authors of Ref. \cite{1-2}, investigated cosmic evolution in an external interacting model of scalar tensor gravity namely Brans Dicke chameleon scenario. The procedure of this work contains novelties but, it shall be observed from this comment, their report faces three fundamental drawbacks. One of them concerns the energy conservation equation, and the other two flaws are about mathematical mistakes. In the scalar tensor gravity models, a well-known method, in order to obtain conservation equation one must combine Friedmann equations with modified Klein-Gordon equation. But in Ref. \cite{1-2}, by virtue of the mentioned approach one would not be able to achieve a correct result for conservation equation. In addition, one can readily realize that their mathematical mistakes lead to different results compared to the present report.


title: Geostatistical methods for disease mapping and visualization using data from spatio-temporally referenced prevalence surveys

abstract: In this paper we set out general principles and develop geostatistical methods for the analysis of data from spatio-temporally referenced prevalence surveys. Our objective is to provide a tutorial guide that can be used in order to identify parsimonious geostatistical models for prevalence mapping. A general variogram-based Monte Carlo procedure is proposed to check the validity of the modelling assumptions. We describe and contrast likelihood-based and Bayesian methods of inference, showing how to account for parameter uncertainty under each of the two paradigms. We also describe extensions of the standard model for disease prevalence that can be used when stationarity of the spatio-temporal covariance function is not supported by the data. We discuss how to define predictive targets and argue that exceedance probabilities provide one of the most effective ways to convey uncertainty in prevalence estimates. We describe statistical software for the visualization of spatio-temporal predictive summaries of prevalence through interactive animations. Finally, we illustrate an application to historical malaria prevalence data from 1334 surveys conducted in Senegal between 1905 and 2014.


title: Anomaly Detection using One-Class Neural Networks

abstract: We propose a one-class neural network (OC-NN) model to detect anomalies in complex data sets. OC-NN combines the ability of deep networks to extract a progressively rich representation of data with the one-class objective of creating a tight envelope around normal data. The OC-NN approach breaks new ground for the following crucial reason: data representation in the hidden layer is driven by the OC-NN objective and is thus customized for anomaly detection. This is a departure from other approaches which use a hybrid approach of learning deep features using an autoencoder and then feeding the features into a separate anomaly detection method like one-class SVM (OC-SVM). The hybrid OC-SVM approach is sub-optimal because it is unable to influence representational learning in the hidden layers. A comprehensive set of experiments demonstrate that on complex data sets (like CIFAR and GTSRB), OC-NN performs on par with state-of-the-art methods and outperformed conventional shallow methods in some scenarios.


title: On Finding Dense Common Subgraphs

abstract: We study the recently introduced problem of finding dense common subgraphs: Given a sequence of graphs that share the same vertex set, the goal is to find a subset of vertices $S$ that maximizes some aggregate measure of the density of the subgraphs induced by $S$ in each of the given graphs. Different choices for the aggregation function give rise to variants of the problem that were studied recently. We settle many of the questions left open by previous works, showing NP-hardness, hardness of approximation, non-trivial approximation algorithms, and an integrality gap for a natural relaxation.


title: Experimental demonstration of three-dimensional broadband underwater acoustic carpet cloak

abstract: We present the design, architecture and detailed performance of a three-dimensional (3D) underwater acoustic carpet cloak (UACC). The proposed system of the 3D UACC is an octahedral pyramid which is composed of periodical steel strips. This underwater acoustic device, placed over the target to hide, is able to manipulate the scattered wavefront to mimic a reflecting plane. The effectiveness of the prototype is experimentally demonstrated in an anechoic tank. The measured acoustic pressure distributions show that the 3D UACC can work in all directions in a wide frequency range. This experimental verification of 3D device paves the way for guidelines on future practical applications.


title: Generalized Bessel multipliers in Hilbert spaces

abstract: The notation of generalized Bessel multipliers is obtained by a bounded operator on $\ell^2$ which is inserted between the analysis and synthesis operators. We show that various properties of generalized multipliers are closely related to their parameters, in particular it will be shown that the membership of generalized Bessel multiplier in the certain operator classes requires that its symbol belongs in the same classes, in special sense. Also, we give some examples to illustrate our results. As we shall see, generalized multipliers associated with Riesz bases are well-behaved, more precisely in this case multipliers can be easily composed and inverted. Special attention is devoted to the study of invertible generalized multipliers. Sufficient and/or necessary conditions for invertibility are determined. Finally, the behavior of these operators under perturbations is discussed.


title: The variation of invariant graphs in forced systems

abstract: In skew-product systems with contractive factors, all orbits asymptotically approach the graph of the so-called sync function; hence, the corresponding regularity properties primarily matter. In the literature, sync function Lipschitz continuity and differentiability have been proved to hold depending on the derivative of the base reciprocal, if not on its Lyapunov exponent. However, forcing topological features can also impact the sync function regularity. Here, we estimate the total variation of sync functions generated by one-dimensional Markov maps. A sharp condition for bounded variation is obtained depending on parameters, that involves the Markov map topological entropy. The results are illustrated with examples.


title: Orbitally limited pair-density wave phase of multilayer superconductors

abstract: We investigate the magnetic field dependence of an ideal superconducting vortex lattice in the parity-mixed pair-density wave phase of multilayer superconductors within a circular cell Ginzburg-Landau approach. In multilayer systems, due to local inversion symmetry breaking, a Rashba spin-orbit coupling is induced at the outer layers. This combined with a perpendicular paramagnetic (Pauli) limiting magnetic field stabilizes a staggered layer dependent pair-density wave phase in the superconducting singlet channel. The high-field pair-density wave phase is separated from the low-field BCS phase by a first-order phase transition. The motivating guiding question in this paper is: what is the minimal necessary Maki parameter $\alpha_M$ for the appearance of the pair-density wave phase of a superconducting trilayer system? To address this problem we generalize the circular cell method for the regular flux-line lattice of a type-II superconductor to include paramagnetic depairing effects. Then, we apply the model to the trilayer system, where each of the layers are characterized by Ginzburg-Landau parameter $\kappa_0$, and a Maki parameter $\alpha_M$. We find that when the spin-orbit Rashba interaction compares to the superconducting condensation energy, the orbitally limited pair-density wave phase stabilizes for Maki parameters $\alpha_M> 10$.


title: On the c-concavity with respect to the quadratic cost on a manifold

abstract: Pushing a little forward an approach proposed by Villani, we are going to prove that in the Riemannian setting the condition $\nabla^2 f< g$ implies that $f$ is $c$-concave with respect to the quadratic cost as soon as it has a sufficiently small $C^1$-norm. From this, we deduce a sufficient condition for the optimality of transport maps.


title: Efficient Sparse-Winograd Convolutional Neural Networks

abstract: Convolutional Neural Networks (CNNs) are computationally intensive, which limits their application on mobile devices. Their energy is dominated by the number of multiplies needed to perform the convolutions. Winograd's minimal filtering algorithm (Lavin, 2015) and network pruning (Han et al., 2015) can reduce the operation count, but these two methods cannot be directly combined $-$ applying the Winograd transform fills in the sparsity in both the weights and the activations. We propose two modifications to Winograd-based CNNs to enable these methods to exploit sparsity. First, we move the ReLU operation into the Winograd domain to increase the sparsity of the transformed activations. Second, we prune the weights in the Winograd domain to exploit static weight sparsity. For models on CIFAR-10, CIFAR-100 and ImageNet datasets, our method reduces the number of multiplications by $10.4\times$, $6.8\times$ and $10.8\times$ respectively with loss of accuracy less than $0.1\%$, outperforming previous baselines by $2.0\times$-$3.0\times$. We also show that moving ReLU to the Winograd domain allows more aggressive pruning.


title: Node Centralities and Classification Performance for Characterizing Node Embedding Algorithms

abstract: Embedding graph nodes into a vector space can allow the use of machine learning to e.g. predict node classes, but the study of node embedding algorithms is immature compared to the natural language processing field because of a diverse nature of graphs. We examine the performance of node embedding algorithms with respect to graph centrality measures that characterize diverse graphs, through systematic experiments with four node embedding algorithms, four or five graph centralities, and six datasets. Experimental results give insights into the properties of node embedding algorithms, which can be a basis for further research on this topic.


title: Linear-Time Algorithm for Long LCF with $k$ Mismatches

abstract: In the Longest Common Factor with $k$ Mismatches (LCF$_k$) problem, we are given two strings $X$ and $Y$ of total length $n$, and we are asked to find a pair of maximal-length factors, one of $X$ and the other of $Y$, such that their Hamming distance is at most $k$. Thankachan et al. show that this problem can be solved in $\mathcal{O}(n \log^k n)$ time and $\mathcal{O}(n)$ space for constant $k$. We consider the LCF$_k$($\ell$) problem in which we assume that the sought factors have length at least $\ell$, and the LCF$_k$($\ell$) problem for $\ell=\Omega(\log^{2k+2} n)$, which we call the Long LCF$_k$ problem. We use difference covers to reduce the Long LCF$_k$ problem to a task involving $m=\mathcal{O}(n/\log^{k+1}n)$ synchronized factors. The latter can be solved in $\mathcal{O}(m \log^{k+1}m)$ time, which results in a linear-time algorithm for Long LCF$_k$. In general, our solution to LCF$_k$($\ell$) for arbitrary $\ell$ takes $\mathcal{O}(n + n \log^{k+1} n/\sqrt{\ell})$ time.


title: Hamiltonian Zoo for systems with one degree of freedom

abstract: We present alternative explicit forms of the standard Hamiltonian for systems with one degree of freedom. This new class of infinite Hamiltonians is called Newton-equivalent Hamiltonian zoo, producing the same equation of motion. These Hamiltonians are directly solved from the Hamilton's equations and come with extra-parameters, which are interpreted as scaling factors for the time evolution on phase space. Moreover, each Hamiltonian in the zoo can be used as a generating function for a Hamiltonian hierarchy.


title: Inductive Framework for Multi-Aspect Streaming Tensor Completion with Side Information

abstract: Low rank tensor completion is a well studied problem and has applications in various fields. However, in many real world applications the data is dynamic, i.e., new data arrives at different time intervals. As a result, the tensors used to represent the data grow in size. Besides the tensors, in many real world scenarios, side information is also available in the form of matrices which also grow in size with time. The problem of predicting missing values in the dynamically growing tensor is called dynamic tensor completion. Most of the previous work in dynamic tensor completion make an assumption that the tensor grows only in one mode. To the best of our Knowledge, there is no previous work which incorporates side information with dynamic tensor completion. We bridge this gap in this paper by proposing a dynamic tensor completion framework called Side Information infused Incremental Tensor Analysis (SIITA), which incorporates side information and works for general incremental tensors. We also show how non-negative constraints can be incorporated with SIITA, which is essential for mining interpretable latent clusters. We carry out extensive experiments on multiple real world datasets to demonstrate the effectiveness of SIITA in various different settings.


title: Strong convergence rates of semi-discrete splitting approximations for stochastic Allen--Cahn equation

abstract: This article analyzes an explicit temporal splitting numerical scheme for the stochastic Allen-Cahn equation driven by additive noise, in a bounded spatial domain with smooth boundary in dimension $d\le 3$. The splitting strategy is combined with an exponential Euler scheme of an auxiliary problem. When $d=1$ and the driving noise is a space-time white noise, we first show some a priori estimates of this splitting scheme. Using the monotonicity of the drift nonlinearity, we then prove that under very mild assumptions on the initial data, this scheme achieves the optimal strong convergence rate $\OO(\delta t^{\frac 14})$. When $d\le 3$ and the driving noise possesses some regularity in space, we study exponential integrability properties of the exact and numerical solutions. Finally, in dimension $d=1$, these properties are used to prove that the splitting scheme has a strong convergence rate $\OO(\delta t)$.


title: Estimation of the linear fractional stable motion

abstract: In this paper we investigate the parametric inference for the linear fractional stable motion in high and low frequency setting. The symmetric linear fractional stable motion is a three-parameter family, which constitutes a natural non-Gaussian analogue of the scaled fractional Brownian motion. It is fully characterised by the scaling parameter $\sigma>0$, the self-similarity parameter $H \in (0,1)$ and the stability index $\alpha \in (0,2)$ of the driving stable motion. The parametric estimation of the model is inspired by the limit theory for stationary increments L\'evy moving average processes that has been recently studied in \cite{BLP}. More specifically, we combine (negative) power variation statistics and empirical characteristic functions to obtain consistent estimates of $(\sigma, \alpha, H)$. We present the law of large numbers and some fully feasible weak limit theorems.


title: Quantum metamaterials: entanglement of spin and orbital angular momentum of a single photon

abstract: Metamaterials have been a major research area for more than two decades now, involving artificial structures with predesigned electromagnetic properties constructed from deep subwavelength building blocks. They have been used to demonstrate a wealth of fascinating phenomena ranging from negative refractive index and epsilon-near-zero to cloaking, emulations of general relativity effects, and super-resolution imaging, to name a few. In the past few years, metamaterials have been suggested as a new platform for quantum optics, and several pioneering experiments have already been carried out with single photons. Here, we employ a dielectric metasurface to generate entanglement between spin and orbital angular momentum of single photons. We demonstrate experimentally the generation of the four Bell states by utilizing the geometric phase arising from the photonic spin-orbit interaction. These are the first experiments with entangled states with metasurfaces, and as such they are paving the way to the new area of quantum metamaterials.


title: Efficient Gradual Typing

abstract: Gradual typing combines static and dynamic typing in the same program. One would hope that the performance in a gradually typed language would range between that of a dynamically typed language and a statically typed language. Existing implementations of gradually typed languages have not achieved this goal due to overheads associated with runtime casts. Takikawa et al. (2016) report up to 100$\times$ slowdowns for partially typed programs. In this paper we present a compiler, named Grift, for evaluating implementation techniques for gradual typing. We take a straightforward but surprisingly unexplored implementation approach for gradual typing, that is, ahead-of-time compilation to native assembly code with carefully chosen runtime representations and space-efficient coercions. Our experiments show that this approach achieves performance on par with OCaml on statically typed programs and performance between that of Gambit and Racket on untyped programs. On partially typed code, the geometric mean ranges from 0.42$\times$ to 2.36$\times$ that of (untyped) Racket across the benchmarks. We implement casts using the coercions of Siek, Thiemann, and Wadler (2015). This technique eliminates all catastrophic slowdowns without introducing significant overhead. Across the benchmarks, coercions range from 15% slower (fft) to almost 2$\times$ faster (matmult) than regular casts. We also implement the monotonic references of Siek et al. (2015). Monotonic references eliminate all overhead in statically typed code, and for partially typed code, they are faster than proxied references, sometimes up to 1.48$\times$.


title: On the entropy norm on the group of diffeomorphisms of closed oriented surface

abstract: We prove that the entropy norm on the group of diffeomorphisms of a closed orientable surface of positive genus is unbounded.


title: On the structure of affine flat group schemes over discrete valuation rings, II

abstract: In the 1st part of this work [DHdS18], we studied affine group schemes over a discrete valuation ring (DVR) by means of Neron blowups. We also showed how to apply these findings to throw light on the group schemes coming from Tannakian categories of D-modules. In the present work, we follow up this theme. We show that a certain class of affine group schemes of "infinite type", Neron blowups of formal subgroups, are quite typical. We also explain how these group schemes appear naturally in Tannakian categories of D-modules. To conclude, we isolate a Tannakian property of affine group schemes, named prudence, which allows one to verify if the underlying ring of functions is a free module over the base ring. This is then successfully applied to obtain a general result on the structure of differential Galois groups over complete DVRs.


title: Neutral gas properties of Lyman continuum emitting galaxies: column densities and covering fractions from UV absorption lines

abstract: The processes allowing the escape of ionizing photons from galaxies into the intergalactic medium are poorly known. To understand how Lyman continuum (LyC) photons escape galaxies, we constrain the HI covering fractions and column densities using ultraviolet HI and metal absorption lines of 18 star-forming galaxies which have Lyman series observations. Nine of these galaxies are confirmed LyC emitters. We fit the stellar continuum, dust attenuation, metal, and HI properties to consistently determine the UV attenuation, as well as the column densities and covering factors of neutral hydrogen and metals. We use synthetic interstellar absorption lines to explore the systematics of our measurements. Then we apply our method to the observed UV spectra of low-redshift and z-2 galaxies. The observed HI lines are found to be saturated in all galaxies. An indirect approach using OI column densities and the observed O/H abundances yields HI column densities of 18.6 to 20 cm-2. These columns are too high to allow the escape of ionizing photons. We find that the known LyC leakers have HI covering fractions less than unity. Ionizing photons escape through optically thin holes/channels in a clumpy interstellar medium. Our simulations confirm that the HI covering fractions are accurately recovered. The SiII and HI covering fractions scale linearly, in agreement with observations from stacked Lyman break galaxy spectra at z-3. Thus, with an empirical correction, the SiII absorption lines can also be used to determine the HI coverage. Finally, we show that a consistent fitting of dust attenuation, continuum and absorption lines is required to properly infer the covering fraction of neutral gas and subsequently to infer the escape fraction of ionizing radiation. These measurements can estimate the LyC escape fraction, as we demonstrate in a companion paper.


title: Long-Term Oscillations of Sunspots and a Special Class of Artifacts in SOHO(MDI) and SDO(HMI) Data

abstract: A specific type of artifacts, that originate due to displacement of the image of a moving object along the digital (pixel) matrix of receiver are analyzed in detail. The criteria of appearance and the influence of these artifacts on the study of long-term oscillations of sunspots are deduced. The obtained criteria suggest us methods for reduction or even elimination of these artifacts. It is shown that the use of integral parameters can be very effective against the artifact distortions. The simultaneous observations of sunspot magnetic field and ultraviolet intensity of the umbra have given the same periods for the long-term oscillations. In this way the real physical nature of the oscillatory process, which is independent of the artifacts have been confirmed again. A number of examples considered here confirm the dependence between the periods of main mode of the sunspot magnetic field long-term oscillations and its strength. The dependence was derived earlier from both the observations and the theoretical model of the shallow sunspot. The anti-phase behavior of time variations of sunspot umbra area and magnetic field of the sunspot demonstrates that the umbra of sunspot moves in long-term oscillations as a whole: all its points oscillate with the same phase.


title: A new discontinuous Galerkin spectral element method for elastic waves with physically motivated numerical fluxes

abstract: The discontinuous Galerkin (DG) method is an established method for computing approximate solutions of partial differential equations in many applications. Unlike continuous finite elements, in DG methods, numerical fluxes are used to enforce inter-element conditions, and internal and external physical boundary conditions. However, for certain problems such as elastic wave propagation in complex media, where several wave types and wave speeds are simultaneously present, a standard numerical flux may not be compatible with the physical boundary conditions. If surface or interface waves are present, this incompatibility may lead to numerical instabilities. We present a stable and arbitrary order accurate DG method for elastic waves with a physically motivated numerical flux. Our numerical flux is compatible with all well-posed, internal and external, boundary conditions, including linear and nonlinear frictional constitutive equations for modelling spontaneously propagating shear ruptures in elastic solids and dynamic earthquake rupture processes. We present numerical experiments in one and two space dimensions verifying high order accuracy and asymptotic numerical stability, and demonstrating potentials for modelling complex nonlinear frictional problems in elastic solids.


title: Smooth maps compatible with simplicial structures and inverse images

abstract: As a higher dimensional version of the theory of Morse functions, there have been various studies of smooth manifolds using generic smooth maps. As fundamental results, in these studies, they have found that inverse images of such maps often restrict the types of the source manifolds. For example, if a generic map such that the inverse image of a regular value is not null-cobordant, then the homology group of its {\it Reeb} space, which is defined as the space of all the connected components of inverse images of the map and a fundamental tool in the theory of generic smooth maps, is known to be non-trivial. In this paper, we show similar results in new appropriate situations. These works are regarded as extensions of works by Hiratuka and Saeki in 2013--4.


title: Space-efficient Feature Maps for String Alignment Kernels

abstract: String kernels are attractive data analysis tools for analyzing string data. Among them, alignment kernels are known for their high prediction accuracies in string classifications when tested in combination with SVM in various applications. However, alignment kernels have a crucial drawback in that they scale poorly due to their quadratic computation complexity in the number of input strings, which limits large-scale applications in practice. We address this need by presenting the first approximation for string alignment kernels, which we call space-efficient feature maps for edit distance with moves (SFMEDM), by leveraging a metric embedding named edit sensitive parsing (ESP) and feature maps (FMs) of random Fourier features (RFFs) for large-scale string analyses. The original FMs for RFFs consume a huge amount of memory proportional to the dimension d of input vectors and the dimension D of output vectors, which prohibits its large-scale applications. We present novel space-efficient feature maps (SFMs) of RFFs for a space reduction from O(dD) of the original FMs to O(d) of SFMs with a theoretical guarantee with respect to concentration bounds. We experimentally test SFMEDM on its ability to learn SVM for large-scale string classifications with various massive string data, and we demonstrate the superior performance of SFMEDM with respect to prediction accuracy, scalability and computation efficiency.


title: Efficient Gaussian Process Classification Using Polya-Gamma Data Augmentation

abstract: We propose a scalable stochastic variational approach to GP classification building on Polya-Gamma data augmentation and inducing points. Unlike former approaches, we obtain closed-form updates based on natural gradients that lead to efficient optimization. We evaluate the algorithm on real-world datasets containing up to 11 million data points and demonstrate that it is up to two orders of magnitude faster than the state-of-the-art while being competitive in terms of prediction performance.


title: Spurious Valleys in Two-layer Neural Network Optimization Landscapes

abstract: Neural networks provide a rich class of high-dimensional, non-convex optimization problems. Despite their non-convexity, gradient-descent methods often successfully optimize these models. This has motivated a recent spur in research attempting to characterize properties of their loss surface that may explain such success. In this paper, we address this phenomenon by studying a key topological property of the loss: the presence or absence of spurious valleys, defined as connected components of sub-level sets that do not include a global minimum. Focusing on a class of two-layer neural networks defined by smooth (but generally non-linear) activation functions, we identify a notion of intrinsic dimension and show that it provides necessary and sufficient conditions for the absence of spurious valleys. More concretely, finite intrinsic dimension guarantees that for sufficiently overparametrised models no spurious valleys exist, independently of the data distribution. Conversely, infinite intrinsic dimension implies that spurious valleys do exist for certain data distributions, independently of model overparametrisation. Besides these positive and negative results, we show that, although spurious valleys may exist in general, they are confined to low risk levels and avoided with high probability on overparametrised models.


title: Connection between diphoton and triboson channels in new physics searches

abstract: Diphoton channel provides a clean signature in searches for new physics. In this paper, we discuss a connection between the diphoton channel ($\gamma\gamma$) and triboson channels ($Z\gamma\gamma$, $ZZ\gamma$, $WW\gamma$) imposed by the $SU(2)_{L}\times U(1)_{Y}$ symmetry of the Standard Model (SM) in certain classes of models. To illustrate this idea we choose a simple model that has all these channels. In this model, the same physics can give rise to $\gamma+$MET instead of $\gamma\gamma$ and 2 bosons plus missing energy instead of 3-boson channels. We analyze existing constraints and previous searches and show that channels $WW\gamma$ and especially $Z\gamma+$MET have a potential to discover new physics at the LHC.


title: How local in time is the no-arbitrage property under capital gains taxes ?

abstract: In frictionless financial markets, no-arbitrage is a local property in time. This means that a discrete time model is arbitrage-free if and only if there does not exist a one-period-arbitrage. With capital gains taxes, this equivalence fails. For a model with a linear tax and one non-shortable risky stock, we introduce the concept of robust local no-arbitrage (RLNA) as the weakest local condition which guarantees dynamic no-arbitrage. Under a sharp dichotomy condition, we prove (RLNA). Since no-one-period-arbitrage is necessary for no-arbitrage, the latter is sandwiched between two local conditions, which allows us to estimate its non-locality. Furthermore, we construct a stock price process such that two long positions in the same stock hedge each other. This puzzling phenomenon that cannot occur in arbitrage-free frictionless markets (or markets with proportional transaction costs) is used to show that no-arbitrage alone does not imply the existence of an equivalent separating measure if the probability space is infinite. Finally, we show that the model with a linear tax on capital gains can be written as a model with proportional transaction costs by introducing several fictitious securities.


title: Quantifying tolerance of a nonlocal multi-qudit state to any local noise

abstract: We present a general approach for quantifying tolerance of a nonlocal N-partite state to any local noise under different classes of quantum correlation scenarios with arbitrary numbers of settings and outcomes at each site. This allows us to derive new precise bounds in d and N on noise tolerances for: (i) an arbitrary nonlocal N-qudit state; (ii) the N-qudit Greenberger-Horne-Zeilinger (GHZ) state; (iii) the N-qubit W state and the N-qubit Dicke states, and to analyse asymptotics of these precise bounds for large N and d.


title: On energy stable discontinuous Galerkin spectral element approximations of the perfectly matched layer for the wave equation

abstract: We develop a provably energy stable discontinuous Galerkin spectral element method (DGSEM) approximation of the perfectly matched layer (PML) for the three and two space dimensional (3D and 2D) linear acoustic wave equations, in first order form, subject to well-posed linear boundary conditions. First, using the well-known complex coordinate stretching, we derive an efficient un-split modal PML for the 3D acoustic wave equation. Second, we prove asymptotic stability of the continuous PML by deriving energy estimates in the Laplace space, for the 3D PML in a heterogeneous acoustic medium, assuming piece-wise constant PML damping. Third, we develop a DGSEM for the wave equation using physically motivated numerical flux, with penalty weights, which are compatible with all well-posed, internal and external, boundary conditions. When the PML damping vanishes, by construction, our choice of penalty parameters yield an upwind scheme and a discrete energy estimate analogous to the continuous energy estimate. Fourth, to ensure numerical stability when PML damping is present, it is necessary to systematically extend the numerical numerical fluxes, and the inter-element and boundary procedures, to the PML auxiliary differential equations. This is critical for deriving discrete energy estimates analogous to the continuous energy estimates. Finally, we propose a procedure to compute PML damping coefficients such that the PML error converges to zero, at the optimal convergence rate of the underlying numerical method. Numerical experiments are presented in 2D and 3D corroborating the theoretical results.


title: Malliavin calculus for the stochastic Cahn-Hilliard / Allen Cahn equation with unbounded noise diffusion

abstract: The stochastic partial differential equation analyzed in this work, is motivated by a simplified mesoscopic physical model for phase separation. It describes pattern formation due to adsorption and desorption mechanisms involved in surface processes, in the presence of a stochastic driving force. This equation is a combination of Cahn-Hilliard and Allen-Cahn type operators with a multiplicative, white, space-time noise of unbounded diffusion. We apply Malliavin calculus, in order to investigate the existence of a density for the stochastic solution $u$. In dimension one, according to the regularity result in \cite{AKM}, $u$ admits continuous paths a.s. Using this property, and inspired by a method proposed in \cite{CW1}, we construct a modified approximating sequence for $u$, which properly treats the new second order Allen-Cahn operator. Under a localization argument, we prove that the Malliavin derivative of $u$ exists locally, and that the law of $u$ is absolutely continuous, establishing thus that a density exists.


title: Current-density implementation for calculating flexoelectric coefficients

abstract: The flexoelectric effect refers to polarization induced in an insulator when a strain gradient is applied. We have developed a first-principles methodology based on density-functional perturbation theory to calculate the elements of the bulk, clamped-ion flexoelectric tensor. In order to determine the transverse and shear components directly from a unit cell calculation, we calculate the current density induced by the adiabatic atomic displacements of a long-wavelength acoustic phonon. Previous implementations based on the charge-density response required supercells to capture these components. Our density-functional-theory implementation requires the development of an expression for the current density that is valid for the case of nonlocal pseudopotentials, and long-wavelength phonon perturbations. We benchmark our methodology on simple systems of isolated noble gas atoms, and apply it to calculate the clamped-ion flexoelectric constants for a variety of technologically important cubic oxides. We also discuss some technical issues that are associated with the definition of current density in a nonlocal pseudopotential context, and their relevance to the calculation of macroscopic response properties of crystals.


title: Orbits in elementary, power-law galaxy bars: 1. Occurence and role of single loops

abstract: Orbits in galaxy bars are generally complex, but simple closed loop orbits play an important role in our conceptual understanding of bars. Such orbits are found in some well-studied potentials, provide a simple model of the bar in themselves, and may generate complex orbit families. The precessing, power ellipse (p-ellipse) orbit approximation provides accurate analytic orbit orbits in symmetric galaxy potentials. It remains useful for finding and fitting simple loop orbits in the frame of a rotating bar with bar-like and symmetric power-law potentials. Second order perturbation theory yields two or fewer simple loop solutions in these potentials. Numerical integrations in the parameter space neighborhood of perturbation solutions reveal zero or one ac- tual loops in a range of such potentials with rising rotation curves. These loops are embedded in a small parameter region of similar, but librating orbits, which have a subharmonic frequency superimposed on the basic loop. These loops and their librat- ing companions support annular bars. Solid bars can be produced in more complex potentials, as shown by an example with power-law indices varying with radius. The power-law potentials can be viewed as the elementary constituents of more complex potentials. Numerical integrations also reveal interesting classes of orbits with multiple loops. In two-dimensional, self-gravitating bars, with power-law potentials, single loop orbits are very rare. This result suggests that gas bars or oval distortions are unlikely to be long-lived, and that complex orbits or three-dimensional structure must support self-gravitating stellar bars.


title: Center-of-Mass-Based Grasp Pose Adaptation Using 3D Range and Force/Torque Sensing

abstract: Lifting objects, whose mass may produce high wrist torques that exceed the hardware strength limits, could lead to unstable grasps or serious robot damage. This work introduces a new Center-of-Mass (CoM)-based grasp pose adaptation method, for picking up objects using a combination of exteroceptive 3D perception and proprioceptive force/torque sensor feedback. The method works in two iterative stages to provide reliable and wrist torque efficient grasps. Initially, a geometric object CoM is estimated from the input range data. In the first stage, a set of hand-size handle grasps are localized on the object and the closest to its CoM is selected for grasping. In the second stage, the object is lifted using a single arm, while the force and torque readings from the sensor on the wrist are monitored. Based on these readings, a displacement to the new CoM estimation is calculated. The object is released and the process is repeated until the wrist torque effort is minimized. The advantage of our method is the blending of both exteroceptive (3D range) and proprioceptive (force/torque) sensing for finding the grasp location that minimizes the wrist effort, potentially improving the reliability of the grasping and the subsequent manipulation task. We experimentally validate the proposed method by executing a number of tests on a set of objects that include handles, using the humanoid robot WALK-MAN.


title: Primordial black holes and uncertainties in the choice of the window function

abstract: Primordial black holes (PBHs) can be produced by the perturbations that exit the horizon during inflationary phase. While inflation models predict the power spectrum of the perturbations in Fourier space, the PBH abundance depends on the probability distribution function (PDF) of density perturbations in real space. In order to estimate the PBH abundance in a given inflation model, we must relate the power spectrum in Fourier space to the PDF in real space by coarse-graining the perturbations with a window function. However, there are uncertainties on what window function should be used, which could change the relation between the PBH abundance and the power spectrum. This is particularly important in considering PBHs with mass $30 M_\odot$ that account for the LIGO events because the required power spectrum is severely constrained by the observations. In this paper, we investigate how large influence the uncertainties on the choice of a window function have over the power spectrum required for LIGO PBHs. As a result, it is found that the uncertainties significantly affect the prediction for the stochastic gravitational waves (GWs) induced by the second order effect of the perturbations. In particular, the pulsar timing array constraints on the produced GWs could disappear for the real-space top-hat window function.


title: Training Big Random Forests with Little Resources

abstract: Without access to large compute clusters, building random forests on large datasets is still a challenging problem. This is, in particular, the case if fully-grown trees are desired. We propose a simple yet effective framework that allows to efficiently construct ensembles of huge trees for hundreds of millions or even billions of training instances using a cheap desktop computer with commodity hardware. The basic idea is to consider a multi-level construction scheme, which builds top trees for small random subsets of the available data and which subsequently distributes all training instances to the top trees' leaves for further processing. While being conceptually simple, the overall efficiency crucially depends on the particular implementation of the different phases. The practical merits of our approach are demonstrated using dense datasets with hundreds of millions of training instances.


title: A quasilinear approach to fully nonlinear parabolic (S)PDEs on $\mathbb{R}^d$

abstract: We study the Cauchy problem for fully nonlinear (stochastic) parabolic partial differential equations. We provide both in deterministic and stochastic case the existence of a maximal defined solution for the problem and we provide suitable blow-up criterion. The key idea is the use of the paradifferential operator calculus in order to reduce the fully nonlinear problem into an abstract quasilinear (stochastic) parabolic equation. This allows us to use some recent results on abstract quasilinear (stochastic) evolution equations in Banach spaces. To do so, we analyse the properties of the paradifferential operator, in light of known results on the boundedness of the $\mathcal{H}^{\infty}$-calculus for pseudodifferential operator. Finally, we extend the theory just developed to cover high order fully nonlinear parabolic (S)PDEs.


title: Can quantum mechanics be considered consistent? a discussion of Frauchinger and Renner's argument

abstract: We discuss the argument proposed in Ref.~\cite{Frauchiger-Renner}, and show that it does not particularly illustrate any inconsistency in quantum mechanics, but rather the well known difficulty often described as the \textquotedblleft shifty split\textquotedblright: the exact point at which the von Neumann reduction postulate should be applied is ill defined. This is the origin of the famous Schr\"{o}dinger's cat or Wigner's friend paradoxes. We investigate the argument of Ref.~\cite{Frauchiger-Renner} and show that it combines statements obtained by different agents assuming very different positions of the shifty split, and therefore applying the reduction postulate in different ways. This results in the introduction of several different state vectors, while such descriptions are considered as incompatible in standard quantum mechanics. To our knowledge, no interpretation of quantum mechanics includes this possibility; the argument thus refers to a new form of quantum mechanics that should be specified more precisely.


title: Semiclassical approach to finite temperature quantum annealing with trapped ions

abstract: Recently it has been demonstrated that an ensemble of trapped ions may serve as a quantum annealer for the number-partitioning problem [Nature Comm. DOI: 10.1038/ncomms11524]. This hard computational problem may be addressed employing a tunable spin glass architecture. Following the proposal of the trapped ions annealer, we study here its robustness against thermal effects, that is, we investigate the role played by thermal phonons. For the efficient description of the system, we use a semiclassical approach, and benchmark it against the exact quantum evolution. The aim is to understand better and characterize how the quantum device approaches a solution of, an otherwise, difficult to solve NP-hard problem.


title: HybridSVD: When Collaborative Information is Not Enough

abstract: We propose a new hybrid algorithm that allows incorporating both user and item side information within the standard collaborative filtering technique. One of its key features is that it naturally extends a simple PureSVD approach and inherits its unique advantages, such as highly efficient Lanczos-based optimization procedure, simplified hyper-parameter tuning and a quick folding-in computation for generating recommendations instantly even in highly dynamic online environments. The algorithm utilizes a generalized formulation of the singular value decomposition, which adds flexibility to the solution and allows imposing the desired structure on its latent space. Conveniently, the resulting model also admits an efficient and straightforward solution for the cold start scenario. We evaluate our approach on a diverse set of datasets and show its superiority over similar classes of hybrid models.


title: Visual-Only Recognition of Normal, Whispered and Silent Speech

abstract: Silent speech interfaces have been recently proposed as a way to enable communication when the acoustic signal is not available. This introduces the need to build visual speech recognition systems for silent and whispered speech. However, almost all the recently proposed systems have been trained on vocalised data only. This is in contrast with evidence in the literature which suggests that lip movements change depending on the speech mode. In this work, we introduce a new audiovisual database which is publicly available and contains normal, whispered and silent speech. To the best of our knowledge, this is the first study which investigates the differences between the three speech modes using the visual modality only. We show that an absolute decrease in classification rate of up to 3.7% is observed when training and testing on normal and whispered, respectively, and vice versa. An even higher decrease of up to 8.5% is reported when the models are tested on silent speech. This reveals that there are indeed visual differences between the 3 speech modes and the common assumption that vocalized training data can be used directly to train a silent speech recognition system may not be true.


title: Elementary Quotient Completions, Church's Thesis, and Partioned Assemblies

abstract: Hyland's effective topos offers an important realizability model for constructive mathematics in the form of a category whose internal logic validates Church's Thesis. It also contains a boolean full sub-quasitopos of "assemblies" where only a restricted form of Church's Thesis survives. In the present paper we compare the effective topos and the quasitopos of assemblies each as the elementary quotient completions of a Lawvere doctrine based on the partitioned assemblies. In that way we can explain why the two forms of Church's Thesis each category satisfies differ by the way each is inherited from specific properties of the doctrine which determines the elementary quotient completion.


title: An Overview of Physical Layer Security with Finite-Alphabet Signaling

abstract: Providing secure communications over the physical layer with the objective of achieving perfect secrecy without requiring a secret key has been receiving growing attention within the past decade. The vast majority of the existing studies in the area of physical layer security focus exclusively on the scenarios where the channel inputs are Gaussian distributed. However, in practice, the signals employed for transmission are drawn from discrete signal constellations such as phase shift keying and quadrature amplitude modulation. Hence, understanding the impact of the finite-alphabet input constraints and designing secure transmission schemes under this assumption is a mandatory step towards a practical implementation of physical layer security. With this motivation, this article reviews recent developments on physical layer security with finite-alphabet inputs. We explore transmit signal design algorithms for single-antenna as well as multi-antenna wiretap channels under different assumptions on the channel state information at the transmitter. Moreover, we present a review of the recent results on secure transmission with discrete signaling for various scenarios including multi-carrier transmission systems, broadcast channels with confidential messages, cognitive multiple access and relay networks. Throughout the article, we stress the important behavioral differences of discrete versus Gaussian inputs in the context of the physical layer security. We also present an overview of practical code construction over Gaussian and fading wiretap channels, and we discuss some open problems and directions for future research.


title: Towards Ultra-High Performance and Energy Efficiency of Deep Learning Systems: An Algorithm-Hardware Co-Optimization Framework

abstract: Hardware accelerations of deep learning systems have been extensively investigated in industry and academia. The aim of this paper is to achieve ultra-high energy efficiency and performance for hardware implementations of deep neural networks (DNNs). An algorithm-hardware co-optimization framework is developed, which is applicable to different DNN types, sizes, and application scenarios. The algorithm part adopts the general block-circulant matrices to achieve a fine-grained tradeoff between accuracy and compression ratio. It applies to both fully-connected and convolutional layers and contains a mathematically rigorous proof of the effectiveness of the method. The proposed algorithm reduces computational complexity per layer from O($n^2$) to O($n\log n$) and storage complexity from O($n^2$) to O($n$), both for training and inference. The hardware part consists of highly efficient Field Programmable Gate Array (FPGA)-based implementations using effective reconfiguration, batch processing, deep pipelining, resource re-using, and hierarchical control. Experimental results demonstrate that the proposed framework achieves at least 152X speedup and 71X energy efficiency gain compared with IBM TrueNorth processor under the same test accuracy. It achieves at least 31X energy efficiency gain compared with the reference FPGA-based work.


title: RadialGAN: Leveraging multiple datasets to improve target-specific predictive models using Generative Adversarial Networks

abstract: Training complex machine learning models for prediction often requires a large amount of data that is not always readily available. Leveraging these external datasets from related but different sources is therefore an important task if good predictive models are to be built for deployment in settings where data can be rare. In this paper we propose a novel approach to the problem in which we use multiple GAN architectures to learn to translate from one dataset to another, thereby allowing us to effectively enlarge the target dataset, and therefore learn better predictive models than if we simply used the target dataset. We show the utility of such an approach, demonstrating that our method improves the prediction performance on the target domain over using just the target dataset and also show that our framework outperforms several other benchmarks on a collection of real-world medical datasets.


title: Using 3D Hahn Moments as A Computational Representation of ATS Drugs Molecular Structure

abstract: The campaign against drug abuse is fought by all countries, most notably on ATS drugs. The technical limitations of the current test kits to detect new brand of ATS drugs present a challenge to law enforcement authorities and forensic laboratories. Meanwhile, new molecular imaging devices which allowed mankind to characterize the physical 3D molecular structure have been recently introduced, and it can be used to remedy the limitations of existing drug test kits. Thus, a new type of 3D molecular structure representation technique should be developed to cater the 3D molecular structure acquired physically using these molecular imaging devices. One of the applications of image processing methods to represent a 3D image is 3D moments, and this study formulates a new 3D moments technique, namely 3D Hahn moments, to represent the 3D molecular structure of ATS drugs. The performance of the proposed technique was analysed using drug chemical structures obtained from UNODC for the ATS drugs, while non-ATS drugs are obtained randomly from ChemSpider database. The evaluation shows the technique is qualified to be further explored in the future works to be fully compatible with ATS drug identification domain.


title: Sums, products and ratios along the edges of a graph

abstract: In their seminal paper Erd\H{o}s and Szemer\'edi formulated conjectures on the size of sumset and product set of integers. The strongest form of their conjecture is about sums and products along the edges of a graph. In this paper we show that this strong form of the Erd\H{o}s-Szemer\'edi conjecture does not hold. We give upper and lower bounds on the cardinalities of sumsets, product sets and ratio sets along the edges of graphs.


title: Predicting B Cell Receptor Substitution Profiles Using Public Repertoire Data

abstract: B cells develop high affinity receptors during the course of affinity maturation, a cyclic process of mutation and selection. At the end of affinity maturation, a number of cells sharing the same ancestor (i.e. in the same "clonal family") are released from the germinal center, their amino acid frequency profile reflects the allowed and disallowed substitutions at each position. These clonal-family-specific frequency profiles, called "substitution profiles", are useful for studying the course of affinity maturation as well as for antibody engineering purposes. However, most often only a single sequence is recovered from each clonal family in a sequencing experiment, making it impossible to construct a clonal-family-specific substitution profile. Given the public release of many high-quality large B cell receptor datasets, one may ask whether it is possible to use such data in a prediction model for clonal-family-specific substitution profiles. In this paper, we present the method "Substitution Profiles Using Related Families" (SPURF), a penalized tensor regression framework that integrates information from a rich assemblage of datasets to predict the clonal-family-specific substitution profile for any single input sequence. Using this framework, we show that substitution profiles from similar clonal families can be leveraged together with simulated substitution profiles and germline gene sequence information to improve prediction. We fit this model on a large public dataset and validate the robustness of our approach on an external dataset. Furthermore, we provide a command-line tool in an open-source software package (https://github.com/krdav/SPURF) implementing these ideas and providing easy prediction using our pre-fit models.


title: Systematic elimination of Stokes divergences emanating from complex phase space caustics

abstract: Stokes phenomenon refers to the fact that the asymptotic expansion of complex functions can differ in different regions of the complex plane, and that beyond the so-called Stokes lines has an unphysical divergence. An important special case is when the Stokes lines emanate from phase space caustics of a complex trajectory manifold. In this case, symmetry determines that to second order there is a double coverage of the space, one portion of which is unphysical. Building on the seminal but laconic findings of Adachi, we show that the deviation from second order can be used to rigorously determine the Stokes lines and therefore the region of the space that should be removed. The method has applications to wavepacket reconstruction from complex valued classical trajectories. With a rigorous method in hand for removing unphysical divergences, we demonstrate excellent wavepacket reconstruction for the Morse, Quartic, Coulomb and Eckart systems.


title: A study of the HI and optical properties of Low Surface Brightness galaxies: spirals, dwarfs and irregulars

abstract: We present a study of the HI and optical properties of nearby ($z$ $\le$ 0.1) Low Surface Brightness galaxies (LSBGs). We started with a literature sample of $\sim$900 LSBGs and divided them into three morphological classes: spirals, irregulars and dwarfs. Of these, we could use $\sim$490 LSBGs to study their HI and stellar masses, colours and colour magnitude diagrams, and local environment, compare them with normal, High Surface Brightness (HSB) galaxies and determine the differences between the three morphological classes. We found that LSB and HSB galaxies span a similar range in HI and stellar masses, and have a similar $M_{\rm HI}$/$M_{\star}$--$M_{\star}$ relationship. Among the LSBGs, as expected, the spirals have the highest average HI and stellar masses, both of about 10$^{9.8}$ $M_\odot$. The LSGBs' ($g$--$r$) integrated colour is nearly constant as function of HI mass for all classes. In the colour magnitude diagram, the spirals are spread over the red and blue regions whereas the irregulars and dwarfs are confined to the blue region. The spirals also exhibit a steeper slope in the $M_{\rm HI}$/$M_{\star}$--$M_{\star}$ plane. Within their local environment we confirmed that LSBGs are more isolated than HSB galaxies, and LSB spirals more isolated than irregulars and dwarfs. Kolmogorov-Smirnov statistical tests on the HI mass, stellar mass and number of neighbours indicates that the spirals are a statistically different population from the dwarfs and irregulars. This suggests that the spirals may have different formation and HI evolution than the dwarfs and irregulars.


title: The many-faceted light curves of young disk-bearing stars in Upper Sco and $\rho$ Oph observed by $K2$ Campaign 2

abstract: The $K2$ Mission has photometrically monitored thousands of stars at high precision and cadence in a series of $\sim$80-day campaigns focused on sections of the ecliptic plane. During its second campaign, $K2$ targeted over 1000 young stellar objects (YSOs) in the $\sim$1-3 Myr $\rho$ Ophiuchus and 5-10 Myr Upper Scorpius regions. From this set, we have carefully vetted photometry from {\em WISE} and {\em Spitzer} to identify those YSOs with infrared excess indicative of primordial circumstellar disks. We present here the resulting comprehensive sample of 288 young disk-bearing stars from B through M spectral types and analysis of their associated $K2$ light curves. Using statistics of periodicity and symmetry, we categorize each light curves into eight different variability classes, notably including "dippers" (fading events), "bursters" (brightening events), stochastic, and quasi-periodic types. Nearly all (96\%) of disk-bearing YSOs are identified as variable at 30-minute cadence with the sub-1\% precision of {\em K2}. Combining our variability classifications with (circum)stellar properties, we find that the bursters, stochastic sources, and the largest amplitude quasi-periodic stars have larger infrared colors, and hence stronger circumstellar disks. They also tend to have larger H$\alpha$ equivalent widths, indicative of higher accretion rates. The dippers, on the other hand, cluster toward moderate infrared colors and low H$\alpha$. Using resolved disk observations, we further find that the latter favor high inclinations, apart from a few notable exceptions with close to face-on disks. These observations support the idea that YSO time domain properties are dependent on several factors including accretion rate and view angle.


title: Emergence of oscillatory behaviors for excitable systems with noise and mean-field interaction, a slow-fast dynamics approach

abstract: We consider the long-time dynamics of a general class of nonlinear Fokker-Planck equations, describing the large population behavior of mean-field interacting units. The main motivation of this work concerns the case where the individual dynamics is excitable, i.e. when each isolated dynamics rests in a stable state, whereas a sufficiently strong perturbation induces a large excursion in the phase space. We address the question of the emergence of oscillatory behaviors induced by noise and interaction in such systems. We tackle this problem by considering this model as a slow-fast system (the mean value of the process giving the slow dynamics) in the regime of small individual dynamics and by proving the existence of a positively stable invariant manifold, whose slow dynamics is at first order the dynamics of a single individual averaged with a Gaussian kernel. We consider applications of this result to Stuart-Landau, FitzHugh-Nagumo and Cucker-Smale oscillators.

categories: math.AP math-ph math.DS math.MP math.PR nlin.AO


title: Magnetoelectric effect and orbital magnetization in skyrmion crystals: Detection and characterization of skyrmions

abstract: Skyrmions are small magnetic quasiparticles, which are uniquely characterized by their topological charge and their helicity. In this Rapid Communication, we show via calculations how both properties can be determined without relying on real-space imaging. The orbital magnetization and topological Hall conductivity measure the arising magnetization due to the circulation of electrons in the bulk and the occurrence of topologically protected edge channels due to the emergent field of a skyrmion crystal. Both observables quantify the topological Hall effect and distinguish skyrmions from antiskyrmions by sign. Additionally, we predict a magnetoelectric effect in skyrmion crystals, which is the generation of a magnetization (polarization) by application of an electric (magnetic) field. This effect is quantified by spin toroidization and magnetoelectric polarizability. The dependence of the transverse magnetoelectric effect on the skyrmion helicity fits that of the classical toroidal moment of the spin texture and allows to differentiate skyrmion helicities: it is largest for Bloch skyrmions and zero for Neel skyrmions. We predict distinct features of the four observables that can be used to detect and characterize skyrmions in experiments.


title: Improved TDNNs using Deep Kernels and Frequency Dependent Grid-RNNs

abstract: Time delay neural networks (TDNNs) are an effective acoustic model for large vocabulary speech recognition. The strength of the model can be attributed to its ability to effectively model long temporal contexts. However, current TDNN models are relatively shallow, which limits the modelling capability. This paper proposes a method of increasing the network depth by deepening the kernel used in the TDNN temporal convolutions. The best performing kernel consists of three fully connected layers with a residual (ResNet) connection from the output of the first to the output of the third. The addition of spectro-temporal processing as the input to the TDNN in the form of a convolutional neural network (CNN) and a newly designed Grid-RNN was investigated. The Grid-RNN strongly outperforms a CNN if different sets of parameters for different frequency bands are used and can be further enhanced by using a bi-directional Grid-RNN. Experiments using the multi-genre broadcast (MGB3) English data (275h) show that deep kernel TDNNs reduces the word error rate (WER) by 6% relative and when combined with the frequency dependent Grid-RNN gives a relative WER reduction of 9%.


title: New spinor classes on the Graf-Clifford algebra

abstract: Pinor and spinor fields are sections of the subbundles whose fibers are the representation spaces of the Clifford algebra of the forms, equipped with the Graf product. In this context, pinors and spinors are here considered and the geometric generalized Fierz identities provide the necessary framework to derive and construct new spinor classes on the space of smooth sections of the exterior bundle, endowed with the Graf product, for prominent specific signatures, whose applications are discussed.


title: Magnetoanisotropic spin-triplet Andreev reflection in ferromagnet-Ising superconductor junctions

abstract: We theoretically study the electronic transport through a ferromagnet-Ising superconductor junction. A tight-binding Hamiltonian describing the Ising superconductor is presented. Then by combing the non-equilibrium Green's function method, the expressions of Andreev reflection coefficient and conductance are obtained. A strong magnetoanisotropic spin-triplet Andreev reflection is shown, and the magnetoanisotropic period is $\pi$ instead of $2\pi$ as in the conventional magnetoanisotropic system. We demonstrate a significant increase of the spin-triplet Andreev reflection for the single-band Ising superconductor. Furthermore, the dependence of the Andreev reflection on the incident energy and incident angle are also investigated. A complete Andreev reflection can occur when the incident energy is equal to the superconductor gap, regardless of the Fermi energy (spin polarization) of the ferromagnet. For the suitable oblique incidence, the spin-triplet Andreev reflection can be strongly enhanced. In addition, the conductance spectroscopies of both zero bias and finite bias are studied, and the influence of gate voltage, exchange energy, and spin-orbit coupling on the conductance spectroscopy are discussed in detail. The conductance reveals a strong magnetoanisotropy with period $\pi$ as the Andreev reflection coefficient. When the magnetization direction is parallel to the junction plane, a large conductance peak always emerges at the superconductor gap. This work offers a comprehensive and systematic study of the spin-triplet Andreev reflection, and has underlying application of $\pi$-periodic spin valve in spintronics.


title: Solving Large-Scale Minimum-Weight Triangulation Instances to Provable Optimality

abstract: We consider practical methods for the problem of finding a minimum-weight triangulation (MWT) of a planar point set, a classic problem of computational geometry with many applications. While Mulzer and Rote proved in 2006 that computing an MWT is NP-hard, Beirouti and Snoeyink showed in 1998 that computing provably optimal solutions for MWT instances of up to 80,000 uniformly distributed points is possible, making use of clever heuristics that are based on geometric insights. We show that these techniques can be refined and extended to instances of much bigger size and different type, based on an array of modifications and parallelizations in combination with more efficient geometric encodings and data structures. As a result, we are able to solve MWT instances with up to 30,000,000 uniformly distributed points in less than 4 minutes to provable optimality. Moreover, we can compute optimal solutions for a vast array of other benchmark instances that are not uniformly distributed, including normally distributed instances (up to 30,000,000 points), all point sets in the TSPLIB (up to 85,900 points), and VLSI instances with up to 744,710 points. This demonstrates that from a practical point of view, MWT instances can be handled quite well, despite their theoretical difficulty.


title: Sim-to-Real Optimization of Complex Real World Mobile Network with Imperfect Information via Deep Reinforcement Learning from Self-play

abstract: Mobile network that millions of people use every day is one of the most complex systems in the world. Optimization of mobile network to meet exploding customer demand and reduce capital/operation expenditures poses great challenges. Despite recent progress, application of deep reinforcement learning (DRL) to complex real world problem still remains unsolved, given data scarcity, partial observability, risk and complex rules/dynamics in real world, as well as the huge reality gap between simulation and real world. To bridge the reality gap, we introduce a Sim-to-Real framework to directly transfer learning from simulation to real world via graph convolutional neural network (CNN) - by abstracting partially observable mobile network into graph, then distilling domain-variant irregular graph into domain-invariant tensor in locally Euclidean space as input to CNN -, domain randomization and multi-task learning. We use a novel self-play mechanism to encourage competition among DRL agents for best record on multiple tasks via simulated annealing, just like athletes compete for world record in decathlon. We also propose a decentralized multi-agent, competitive and cooperative DRL method to coordinate the actions of multi-cells to maximize global reward and minimize negative impact to neighbor cells. Using 6 field trials on commercial mobile networks, we demonstrate for the first time that a DRL agent can successfully transfer learning from simulation to complex real world problem with imperfect information, complex rules/dynamics, huge state/action space, and multi-agent interactions, without any training in the real world.


title: Superthermal photon bunching in terms of simple probability distributions

abstract: We analyze the second-order photon autocorrelation function $g^{(2)}$ with respect to the photon probability distribution and discuss the generic features of a distribution that result in superthermal photon bunching ($g^{(2)}>2$). Superthermal photon bunching has been reported for a number of optical microcavity systems that exhibit processes like superradiance or mode competition. We show that a superthermal photon number distribution cannot be constructed from the principle of maximum entropy, if only the intensity and the second-order autocorrelation are given. However, for bimodal systems an unbiased superthermal distribution can be constructed from second-order correlations and the intensities alone. Our findings suggest modeling superthermal single-mode distributions by a mixture of a thermal and a lasing like state and thus reveal a generic mechanism in the photon probability distribution responsible for creating superthermal photon bunching. We relate our general considerations to a physical system, a (single-emitter) bimodal laser, and show that its statistics can be approximated and understood within our proposed model. Furthermore the excellent agreement of the statistics of the bimodal laser and our model reveal that the bimodal laser is an ideal source of bunched photons, in the sense that it can generate statistics that contain no other features but the superthermal bunching.