title: Performance Dynamics and Success in Online Games

abstract: Online data provide a way to monitor how users behave in social systems like social networks and online games, and understand which features turn an ordinary individual into a successful one. Here, we propose to study individual performance and success in Multiplayer Online Battle Arena (MOBA) games. Our purpose is to identify those behaviors and playing styles that are characteristic of players with high skill level and that distinguish them from other players. To this aim, we study Defense of the ancient 2 (Dota 2), a popular MOBA game. Our findings highlight three main aspects to be successful in the game: (i) players need to have a warm-up period to enhance their performance in the game; (ii) having a long in-game experience does not necessarily translate in achieving better skills; but rather, (iii) players that reach high skill levels differentiate from others because of their aggressive playing strategy, which implies to kill opponents more often than cooperating with teammates, and trying to give an early end to the match.


title: Depth of Ultra High Energy Cosmic Ray Induced Air Shower Maxima Measured by the Telescope Array Black Rock and Long Ridge FADC Fluorescence Detectors and Surface Array in Hybrid Mode

abstract: The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.


title: Optical Kerr Effect in Graphene: Theoretical Analysis of the Optical Heterodyne Detection Technique

abstract: Graphene is an atomically thin two-dimensional material demonstrating strong optical nonlinearities including harmonics generation, four wave mixing, Kerr and other nonlinear effects. In this paper we theoretically analyze the optical heterodyne detection (OHD) technique of measuring the optical Kerr effect (OKE) in two-dimensional crystals and show how to relate the quantities measured in such experiments with components of the third-order conductivity tensor $\sigma^{(3)}_{\alpha\beta\gamma\delta}(\omega_1,\omega_2,\omega_3)$ of the two-dimensional crystal. Using results of a recently developed quantum theory of the third-order nonlinear electrodynamic response of graphene we analyze the frequency, charge carrier density, temperature and other dependencies of the OHD-OKE response of this material. We compare our results with a recent OHD-OKE experiment in graphene and find good agreement between the theory and experiment.


title: An Algorithm to Decompose Permutation Representations of Finite Groups: Polynomial Algebra Approach

abstract: We describe an algorithm for splitting permutation representations of finite group over fields of characteristic zero into irreducible components. The algorithm is based on the fact that the components of the invariant inner product in invariant subspaces are operators of projection into these subspaces. An important element of the algorithm is the calculation of Gr\"obner bases of polynomial ideals. A preliminary implementation of the algorithm splits representations up to dimensions of tens of thousands. Some examples of computations are given in appendix.


title: Kaleidoscopic groups: permutation groups constructed from dendrite homeomorphisms

abstract: Given a transitive permutation group, a fundamental object for studying its higher transitivity properties is the permutation action of its isotropy subgroup. We reverse this relationship and introduce a universal construction of infinite permutation groups that takes as input a given system of imprimitivity for its isotropy subgroup. This produces vast families kaleidoscopic groups. We investigate their algebraic properties, such as simplicity and oligomorphy; their homological properties, such as acyclicity or contrariwise large Schur multipliers; their topological properties, such as unique polishability. Our construction is carried out within the framework of homeomorphism groups of topological dendrites.


title: Evaluating approaches for supervised semantic labeling

abstract: Relational data sources are still one of the most popular ways to store enterprise or Web data, however, the issue with relational schema is the lack of a well-defined semantic description. A common ontology provides a way to represent the meaning of a relational schema and can facilitate the integration of heterogeneous data sources within a domain. Semantic labeling is achieved by mapping attributes from the data sources to the classes and properties in the ontology. We formulate this problem as a multi-class classification problem where previously labeled data sources are used to learn rules for labeling new data sources. The majority of existing approaches for semantic labeling have focused on data integration challenges such as naming conflicts and semantic heterogeneity. In addition, machine learning approaches typically have issues around class imbalance, lack of labeled instances and relative importance of attributes. To address these issues, we develop a new machine learning model with engineered features as well as two deep learning models which do not require extensive feature engineering. We evaluate our new approaches with the state-of-the-art.


title: Preserving of the unconditional basis property under non-self-adjoint perturbations of self-adjoint operators

abstract: Let $T$ be a self-adjoint operator in a Hilbert space $H$ with domain $\mathcal D(T)$. Assume that the spectrum of $T$ is confined in the union of disjoint intervals $\Delta_k =[\alpha_{2k-1},\alpha_{2k}]$, $k\in \mathbb{Z}$, and $$ \alpha_{2k+1}-\alpha_{2k} \geq b|\alpha_{2k+1}+\alpha_{2k}|^p\quad \text{ for some }\,b>0,\,p\in[0,1). $$ Suppose that a linear operator $B$ in $H$ is $p$-subordinated to $T$, i.e. $\mathcal D(B) \supset\mathcal D(T)$ and $\|Bx\| \leq b'\,\|Tx\|^p\|x\|^{1-p} +M\|x\| \text{\, for all } x\in \mathcal D(T)$, with some $b'>0$ and $M\geq 0$. Then the spectrum of the perturbed operator $A=T+B$ lies in the union of a rectangle in $\mathbb{C}$ and double parabola $P_{p,h} = \bigl\{\lambda \in \mathbb{C}\,\bigl|\,|\mathop{\rm Im} \lambda|\leq h|\mathop{\rm Re} \lambda|^p\bigr\}$, provided that $h>b'$. The vertical strips $\Omega_k =\{\lambda\in\mathbb{C}|\,|r_k-{\rm Re}\,\lambda|\leq \delta r_k^p\}$, $r_k =(\alpha_{2k}+\alpha_{2k+1})/2$, belong to the resolvent set of $T$, provided that $\delta <b -b'$ and $ |k|\geq N$ for $N$ large enough. For $|k|\ge N+1$, denote by $\Pi_k$ the curvilinear trapezoid formed by the lines ${\rm Re}\,\lambda = r_{k-1}$, ${\rm Re}\,\lambda = r_{k}$, and the boundary of the parabola $P_{p,h}$. Assume that $Q_0$ is the Riesz projection corresponding to the (bounded) part of the spectrum of $T$ that lies outside $\bigcup_{|k|\geq N+1}\Pi_k$. And let $Q_k$, $|k|\ge N+1$, be the Riesz projection for the part of the spectrum of $T$ confined within $\Pi_k$. Main result of the work consists in proving that the system of the invariant subspaces $Q_k(H)$, $|k|\geq N+1$, together with the invariant subspace $Q_0(H)$ forms an unconditional basis of subspaces in the space $H$.


title: $\Delta(1232)$-Resonance in the Hydrogen Spectrum

abstract: The electromagnetic excitation of the $\Delta(1232)$-resonance plays an appreciable role in the Lamb shift and hyperfine structure of muonic and electronic hydrogen. Its effect appears at the subleading order $\mathcal{O}(\alpha^5)$, together with other proton-polarizability contributions from forward two-photon exchange. We use the large-$N_c$ relations for the nucleon-to-delta transition form factors to compute the effect of the $\Delta(1232)$ in the hydrogen spectrum. We pay particular attention to a subtile difference between predictions based on a direct calculation of the two-photon exchange (or Compton scattering amplitudes) and predictions based on the $\Delta(1232)$-production photoabsorption cross sections. The mismatch is explained by studying the dispersion relations for tree-level Compton scattering off the proton in more details.


title: Noise-induced Stop-and-Go Dynamics

abstract: Stop-and-go waves are commonly observed in traffic and pedestrian flows. In traffic theory they are described by phase transitions of metastable models. The self-organization phenomenon occurs due to inertia mechanisms but requires fine tuning of the parameters. Here, a novel explanation for stop-and-go waves based on stochastic effects is presented for pedestrian dynamics. We show that the introduction of specific coloured noises in a stable microscopic model allows to describe realistic pedestrian stop-and-go behaviour without requirement of metastability and phase transition. We compare simulation results of the stochastic model to real pedestrian trajectories and discuss plausible values for the model's parameters.


title: Time domain boundary elements for dynamic contact problems

abstract: This article considers a unilateral contact problem for the wave equation. The problem is reduced to a variational inequality for the Dirichlet-to-Neumann operator for the wave equation on the boundary, which is solved in a saddle point formulation using boundary elements in the time domain. As a model problem, also a variational inequality for the single layer operator is considered. A priori estimates are obtained for Galerkin approximations both to the variational inequality and the mixed formulation in the case of a flat contact area, where the existence of solutions to the continuous problem is known. Numerical experiments demonstrate the performance of the proposed mixed method. They indicate the stability and convergence beyond flat geometries.


title: K\"ahler fibrations in quantum information theory

abstract: We discuss the fibre bundle of co-adjoint orbits of compact Lie groups, and show how it admits a compatible K\"ahler structure. The case of the unitary group allows us to reformulate the geometric framework of quantum information theory. In particular, we show that the Fisher information tensor gives rise to a structure that is sufficiently close to a K\"ahler structure to generalise some classical result on co-adjoint orbits.


title: FEAST Eigensolver for Nonlinear Eigenvalue Problems

abstract: The linear FEAST algorithm is a method for solving linear eigenvalue problems. It uses complex contour integration to calculate the eigenvectors whose eigenvalues that are located inside some user-defined region in the complex plane. This makes it possible to parallelize the process of solving eigenvalue problems by simply dividing the complex plane into a collection of disjoint regions and calculating the eigenpairs in each region independently of the eigenpairs in the other regions. In this paper we present a generalization of the linear FEAST algorithm that can be used to solve nonlinear eigenvalue problems. Like its linear progenitor, the nonlinear FEAST algorithm can be used to solve nonlinear eigenvalue problems for the eigenpairs whose eigenvalues lie in a user-defined region in the complex plane, thereby allowing for the calculation of large numbers of eigenpairs in parallel. We describe the nonlinear FEAST algorithm, and use several physically-motivated examples to demonstrate its properties.


title: Reparametrization of COM-Poisson Regression Models with Applications in the Analysis of Experimental Data

abstract: In the analysis of count data often the equidispersion assumption is not suitable, hence the Poisson regression model is inappropriate. As a generalization of the Poisson distribution, the COM-Poisson distribution can deal with under-, equi- and overdispersed count data. It is a member of the exponential family of distributions and has well known special cases. In spite of the nice properties of the COM-Poisson distribution, its location parameter does not correspond to the expectation, which complicates the interpretation of regression models. In this paper, we propose a straightforward reparametrization of the COM-Poisson distribution based on an approximation to the expectation of this distribution. The main advantage of our new parametrization is the straightforward interpretation of the regression coefficients in terms of the expectation, as usual in the context of generalized linear models. Furthermore, the estimation and inference for the new COM-Poisson regression model can be done based on the likelihood paradigm. We carried out simulation studies to verify the finite sample properties of the maximum likelihood estimators. The results from our simulation study show that the maximum likelihood estimators are unbiased and consistent for both regression and dispersion parameters. We observed that the empirical correlation between the regression and dispersion parameter estimators is close to zero, which suggests that these parameters are orthogonal. We illustrate the application of the proposed model through the analysis of three data sets with over-, under- and equidispersed count data. The study of distribution properties through a consideration of dispersion, zero-inflated and heavy tail indexes, together with the results of data analysis show the flexibility over standard approaches.


title: Communication-Efficient Search for an Approximate Closest Lattice Point

abstract: We consider the problem of finding the closest lattice point to a vector in n-dimensional Euclidean space when each component of the vector is available at a distinct node in a network. Our objectives are (i) minimize the communication cost and (ii) obtain the error probability. The approximate closest lattice point considered here is the one obtained using the nearest-plane (Babai) algorithm. Assuming a triangular special basis for the lattice, we develop communication-efficient protocols for computing the approximate lattice point and determine the communication cost for lattices of dimension n>1. Based on available parameterizations of reduced bases, we determine the error probability of the nearest plane algorithm for two dimensional lattices analytically, and present a computational error estimation algorithm in three dimensions. For dimensions 2 and 3, our results show that the error probability increases with the packing density of the lattice.


title: Discrete Autoencoders for Sequence Models

abstract: Recurrent models for sequences have been recently successful at many tasks, especially for language modeling and machine translation. Nevertheless, it remains challenging to extract good representations from these models. For instance, even though language has a clear hierarchical structure going from characters through words to sentences, it is not apparent in current language models. We propose to improve the representation in sequence models by augmenting current approaches with an autoencoder that is forced to compress the sequence through an intermediate discrete latent space. In order to propagate gradients though this discrete representation we introduce an improved semantic hashing technique. We show that this technique performs well on a newly proposed quantitative efficiency measure. We also analyze latent codes produced by the model showing how they correspond to words and phrases. Finally, we present an application of the autoencoder-augmented model to generating diverse translations.


title: Earthmover Resilience and Testing in Ordered Structures

abstract: One of the main challenges in property testing is to characterize those properties that are testable with a constant number of queries. For unordered structures such as graphs and hypergraphs this task has been mostly settled. However, for ordered structures such as strings, images, and ordered graphs, the characterization problem seems very difficult in general. In this paper, we identify a wide class of properties of ordered structures - the earthmover resilient (ER) properties - and show that the "good behavior" of such properties allows us to obtain general testability results that are similar to (and more general than) those of unordered graphs. A property P is ER if, roughly speaking, slight changes in the order of the elements in an object satisfying P cannot make this object far from P. The class of ER properties includes, e.g., all unordered graph properties, many natural visual properties of images, such as convexity, and all hereditary properties of ordered graphs and images. A special case of our results implies, building on a recent result of Alon and the authors, that the distance of a given image or ordered graph from any hereditary property can be estimated (with good probability) up to a constant additive error, using a constant number of queries.


title: Matrix Completion for Low-Observability Voltage Estimation

abstract: With the rising penetration of distributed energy resources, distribution system control and enabling techniques such as state estimation have become essential to distribution system operation. However, traditional state estimation techniques have difficulty coping with the low-observability conditions often present on the distribution system due to the paucity of sensors and heterogeneity of measurements. To address these limitations, we propose a distribution system state estimation algorithm that employs matrix completion (a tool for estimating missing values in low-rank matrices) augmented with noise-resilient power flow constraints. This method operates under low-observability conditions where standard least-squares-based methods cannot operate, and flexibly incorporates any network quantities measured in the field. We empirically evaluate our method on the IEEE 33- and 123-bus test systems, and find that it provides near-perfect state estimation performance (within 1\% mean absolute percent error) across many low-observability data availability regimes.


title: Reducing triangular systems of ODEs with rational coefficients, with applications to coupled Regge-Wheeler equations

abstract: We concisely summarize a method of finding all rational solutions to an inhomogeneous rational ODE system of arbitrary order (but solvable for its highest order terms) by converting it into a finite dimensional linear algebra problem. This method is then used to solve the problem of conclusively deciding when certain rational ODE systems in upper triangular form can or cannot be reduced to diagonal form by differential operators with rational coefficients. As specific examples, we consider systems of coupled Regge-Wheeler equations, which have naturally appeared in previous work on vector and tensor perturbations on the Schwarzschild black hole spacetime. Our systematic approach reproduces and complements identities that have been previously found by trial and error methods.


title: A NIKA view of two star-forming infrared dark clouds: Dust emissivity variations and mass concentration

abstract: The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.


title: Automatically Leveraging MapReduce Frameworks for Data-Intensive Applications

abstract: MapReduce is a popular programming paradigm for developing large-scale, data-intensive computation. Many frameworks that implement this paradigm have recently been developed. To leverage these frameworks, however, developers must become familiar with their APIs and rewrite existing code. Casper is a new tool that automatically translates sequential Java programs into the MapReduce paradigm. Casper identifies potential code fragments to rewrite and translates them in two steps: (1) Casper uses program synthesis to search for a program summary (i.e., a functional specification) of each code fragment. The summary is expressed using a high-level intermediate language resembling the MapReduce paradigm and verified to be semantically equivalent to the original using a theorem prover. (2) Casper generates executable code from the summary, using either the Hadoop, Spark, or Flink API. We evaluated Casper by automatically converting real-world, sequential Java benchmarks to MapReduce. The resulting benchmarks perform up to 48.2x faster compared to the original.


title: ExoCross: a general program for generating spectra from molecular line lists

abstract: ExoCross is a Fortran code for generating spectra (emission, absorption) and thermodynamic properties (partition function, specific heat etc.) from molecular line lists. Input is taken in several formats, including ExoMol and HITRAN formats. ExoCross is efficiently parallelized showing also a high degree of vectorization. It can work with several line profiles such as Doppler, Lorentzian and Voigt and support several broadening schemes. Voigt profiles are handled by several methods allowing fast and accurate simulations. Two of these methods are new. ExoCross is also capable of working with the recently proposed method of super-lines. It supports calculations of lifetimes, cooling functions, specific heats and other properties. ExoCross can be used to convert between different formats, such as HITRAN, ExoMol and Phoenix. It is capable of simulating non-LTE spectra using a simple two-temperature approach. Different electronic, vibronic or vibrational bands can be simulated separately using an efficient filtering scheme based on the quantum numbers.


title: Predicting Rapid Fire Growth (Flashover) Using Conditional Generative Adversarial Networks

abstract: A flashover occurs when a fire spreads very rapidly through crevices due to intense heat. Flashovers present one of the most frightening and challenging fire phenomena to those who regularly encounter them: firefighters. Firefighters' safety and lives often depend on their ability to predict flashovers before they occur. Typical pre-flashover fire characteristics include dark smoke, high heat, and rollover ("angel fingers") and can be quantified by color, size, and shape. Using a color video stream from a firefighter's body camera, we applied generative adversarial neural networks for image enhancement. The neural networks were trained to enhance very dark fire and smoke patterns in videos and monitor dynamic changes in smoke and fire areas. Preliminary tests with limited flashover training videos showed that we predicted a flashover as early as 55 seconds before it occurred.


title: Parameter Box: High Performance Parameter Servers for Efficient Distributed Deep Neural Network Training

abstract: Most work in the deep learning systems community has focused on faster inference, but arriving at a trained model requires lengthy experiments. Accelerating training lets developers iterate faster and come up with better models. DNN training is often seen as a compute-bound problem, best done in a single large compute node with many GPUs. As DNNs get bigger, training requires going distributed. Distributed deep neural network (DDNN) training constitutes an important workload on the cloud. Larger DNN models and faster compute engines shift the training performance bottleneck from computation to communication. Our experiments show existing DNN training frameworks do not scale in a typical cloud environment due to insufficient bandwidth and inefficient parameter server software stacks.We propose PBox, a balanced, scalable central PS hardware that balances compute and communication resources, and PHub, a high performance parameter server (PS) software design that provides an optimized network stack and a streamlined gradient processing pipeline to benefit common PS setups to utilize PBox. We show that in a typical cloud environment, PBox can achieve up to 3.8x speedup over state-of-the-art designs when training ImageNet. We discuss future directions of integrating PBox with programmable switches for in-network aggregation during training, leveraging the datacenter network topology to reduce bandwidth usage and localize data movement.


title: A proof of Jones' conjecture

abstract: In this paper, we prove that Wright's equation $y'(t) = - \alpha y(t-1) \{1 + y(t)\}$ has a unique slowly oscillating periodic solution for parameter values $\alpha \in (\tfrac{\pi}{2}, 1.9]$, up to time translation. This result proves Jones' Conjecture formulated in 1962, that there is a unique slowly oscillating periodic orbit for all $ \alpha > \tfrac{\pi}{2}$. Furthermore, there are no isolas of periodic solutions to Wright's equation; all periodic orbits arise from Hopf bifurcations.


title: A QCA Layout Design Methodology. Part I

abstract: Quantum-dot Cellular Automata (QCA) as a nanoscale transistor-less device technology offers distinguishing advantages over the limitations of CMOS circuits. While more than 2 decades of design progress have been achieved with QCA, a comprehensive composition approach for the layout design in this technology is needed. In this study, the Priority-Phased Decomposition-Driven (PPDD) QCA logic design methodology is proposed. In this methodology, after partitioning combinational and sequential circuits into primary-level priority blocks including 2:1 MUX and XOR gates, and secondary-level priority blocks comprising multi-input majority gates, there are three streamlined approaches which are developed to compose the desired QCA circuit using blocks from each priority level.


title: The Intriguing Properties of Model Explanations

abstract: Linear approximations to the decision boundary of a complex model have become one of the most popular tools for interpreting predictions. In this paper, we study such linear explanations produced either post-hoc by a few recent methods or generated along with predictions with contextual explanation networks (CENs). We focus on two questions: (i) whether linear explanations are always consistent or can be misleading, and (ii) when integrated into the prediction process, whether and how explanations affect the performance of the model. Our analysis sheds more light on certain properties of explanations produced by different methods and suggests that learning models that explain and predict jointly is often advantageous.


title: A Distributed-Memory Algorithm for Computing a Heavy-Weight Perfect Matching on Bipartite Graphs

abstract: We design and implement an efficient parallel algorithm for finding a perfect matching in a weighted bipartite graph such that weights on the edges of the matching are large. This problem differs from the maximum weight matching problem, for which scalable approximation algorithms are known. It is primarily motivated by finding good pivots in scalable sparse direct solvers before factorization. Due to the lack of scalable alternatives, distributed solvers use sequential implementations of maximum weight perfect matching algorithms, such as those available in MC64. To overcome this limitation, we propose a fully parallel distributed memory algorithm that first generates a perfect matching and then iteratively improves the weight of the perfect matching by searching for weight-increasing cycles of length four in parallel. For most practical problems the weights of the perfect matchings generated by our algorithm are very close to the optimum. An efficient implementation of the algorithm scales up to 256 nodes (17,408 cores) on a Cray XC40 supercomputer and can solve instances that are too large to be handled by a single node using the sequential algorithm.


title: Personalized Survival Prediction with Contextual Explanation Networks

abstract: Accurate and transparent prediction of cancer survival times on the level of individual patients can inform and improve patient care and treatment practices. In this paper, we design a model that concurrently learns to accurately predict patient-specific survival distributions and to explain its predictions in terms of patient attributes such as clinical tests or assessments. Our model is flexible and based on a recurrent network, can handle various modalities of data including temporal measurements, and yet constructs and uses simple explanations in the form of patient- and time-specific linear regression. For analysis, we use two publicly available datasets and show that our networks outperform a number of baselines in prediction while providing a way to inspect the reasons behind each prediction.


title: Operational Markov condition for quantum processes

abstract: We derive a necessary and sufficient condition for a quantum process to be Markovian which coincides with the classical one in the relevant limit. Our condition unifies all previously known definitions for quantum Markov processes by accounting for all potentially detectable memory effects. We then derive a family of measures of non-Markovianity with clear operational interpretations, such as the size of the memory required to simulate a process, or the experimental falsifiability of a Markovian hypothesis.


title: Enabling Low-power Duplex Visible Light Communication

abstract: The new generation of LED-based illuminating infrastructures has enabled a "dual-paradigm" where LEDs are used for both illumination and communication purposes. The ubiquity of lighting makes visible light communication (VLC) well suited for communication with mobile devices and sensor nodes in indoor environment. Existing research on VLC has primarily been focused on advancing the performance of one-way communication. In this paper, we present Retro-VLC, a low-power duplex VLC system that enables a mobile device to perform bi-directional communication with the illuminating LEDs over the same light carrier. The design features a retro-reflector fabric that backscatters light, an LCD shutter that modulates information bits on the backscattered light carrier, and several low-power optimization techniques. We have prototyped the Reader system and made a few battery-free tag devices. Experimental results show that the tag can achieve a 10kbps downlink speed and 0.5kbps uplink speed over a distance of 2.4m. We also outline several potential applications of the proposed Retro-VLC system.


title: Subgraph counts for dense random graphs with specified degrees

abstract: We prove two estimates for the expectation of the exponential of a complex function of a random permutation or subset. Using this theory, we find asymptotic expressions for the expected number of copies and induced copies of a given graph in a uniformly random graph with degree sequence $(d_1,\ldots,d_n)$ as $n \rightarrow \infty$. We also determine the expected number of spanning trees in this model. The range of degrees covered includes $d_j = \lambda n + O(n^{1/2+\varepsilon})$ for some $\lambda$ bounded away from $0$ and $1$.


title: Hardy's paradox according to non-classical semantics

abstract: In the paper, using the language of spin-half particles, Hardy's paradox is examined within different semantics: a partial one, a many-valued one, and one defined as a set of weak values of projection operators. As it is shown in this paper, any of such non-classical semantics can resolve Hardy's paradox.


title: A brief survey on singularities of geodesic flows in smooth signature changing metrics on 2-surfaces

abstract: We present a survey on generic singularities of geodesic flows in smooth signature changing metrics (often called pseudo-Riemannian) in dimension 2. Generically, a pseudo-Riemannian metric on a 2-manifold $S$ changes its signature (degenerates) along a curve $S_0$, which locally separates $S$ into a Riemannian ($R$) and a Lorentzian ($L$) domain. The geodesic flow does not have singularities over $R$ and $L$, and for any point $q \in R \cup L$ and every tangential direction $p$ there exists a unique geodesic passing through the point $q$ with the direction $p$. On the contrary, geodesics cannot pass through a point $q \in S_0$ in arbitrary tangential directions, but only in some admissible directions; the number of admissible directions is 1 or 2 or 3. We study this phenomenon and the local properties of geodesics near $q \in S_0$.


title: Spatiotemporal intermittency and localized dynamic fluctuations upon approaching the glass transition

abstract: We introduce a new and robust approach for characterizing spatially and temporally heterogeneous behavior within a system based on the evolution of dynamic fuctuations once averaged over different space lengths and time scales. We apply it to investigate the dynamics in two canonical systems as the glass transition is approached: simulated Lennard-Jones liquids and experimental dense colloidal suspensions. We find that in both cases the onset of glassines is marked by spatially localized dynamic fluctuations originating in regions of correlated mobile particles. By removing the trivial system size dependence of the fluctuations we show that such regions contain tens to hundreds of particles for time scales corresponding to maximally non-Gaussian dynamics.


title: Model-assisted inference for treatment effects using regularized calibrated estimation with high-dimensional data

abstract: Consider the problem of estimating average treatment effects when a large number of covariates are used to adjust for possible confounding through outcome regression and propensity score models. The conventional approach of model building and fitting iteratively can be difficult to implement, depending on ad hoc choices of what variables are included. In addition, uncertainty from the iterative process of model selection is complicated and often ignored in subsequent inference about treatment effects. We develop new methods and theory to obtain not only doubly robust point estimators for average treatment effects, which remain consistent if either the propensity score model or the outcome regression model is correctly specified, but also model-assisted confidence intervals, which are valid when the propensity score model is correctly specified but the outcome regression model may be misspecified. With a linear outcome model, the confidence intervals are doubly robust, that is, being also valid when the outcome model is correctly specified but the propensity score model may be misspecified. Our methods involve regularized calibrated estimators with Lasso penalties, but carefully chosen loss functions, for fitting propensity score and outcome regression models. We provide high-dimensional analysis to establish the desired properties of our methods under comparable conditions to previous results, which give valid confidence intervals when both the propensity score and outcome regression are correctly specified. We present a simulation study and an empirical application which confirm the advantages of the proposed methods compared with related methods based on regularized maximum likelihood estimation.


title: On the crossroads of enumerative geometry and geometric representation theory

abstract: The subjects in the title are interwoven in many different and very deep ways. I recently wrote several expository accounts [64-66] that reflect a certain range of developments, but even in their totality they cannot be taken as a comprehensive survey. In the format of a 30-page contribution aimed at a general mathematical audience, I have decided to illustrate some of the basic ideas in one very interesting example - that of HilbpC2, nq, hoping to spark the curiosity of colleagues in those numerous fields of study where one should expect applications.


title: Transformation Autoregressive Networks

abstract: The fundamental task of general density estimation $p(x)$ has been of keen interest to machine learning. In this work, we attempt to systematically characterize methods for density estimation. Broadly speaking, most of the existing methods can be categorized into either using: \textit{a}) autoregressive models to estimate the conditional factors of the chain rule, $p(x_{i}\, |\, x_{i-1}, \ldots)$; or \textit{b}) non-linear transformations of variables of a simple base distribution. Based on the study of the characteristics of these categories, we propose multiple novel methods for each category. For example we proposed RNN based transformations to model non-Markovian dependencies. Further, through a comprehensive study over both real world and synthetic data, we show for that jointly leveraging transformations of variables and autoregressive conditional models, results in a considerable improvement in performance. We illustrate the use of our models in outlier detection and image modeling. Finally we introduce a novel data driven framework for learning a family of distributions.


title: A consequence of the growth of rotation sets for families of diffeomorphisms of the torus

abstract: In this paper we consider $C^\infty $-generic families of area-preserving diffeomorphisms of the torus homotopic to the identity and their rotation sets. Let $f_t:\rm{T^2\rightarrow T^2}$ be such a family, $\widetilde{f}_t:\rm I\negthinspace R^2 \rightarrow \rm I\negthinspace R^2$ be a fixed family of lifts and $\rho (\widetilde{f}_t)$ be their rotation sets, which we assume to have interior for $t$ in a certain open interval $I.$ We also assume that some rational point $(\frac pq,\frac rq)\in \partial \rho (\widetilde{f}_{\overline{t}})$ for a certain parameter $\overline{t}\in I$ and we want to understand consequences of the following hypothesis: For all $t>\overline{t},$ $t\in I,$ $(\frac pq,\frac rq)\in int(\partial \rho (\widetilde{f}_t)).$ Under these very natural assumptions, we prove that there exists a $f_{\overline{t}}^q$-fixed hyperbolic saddle $P_{\overline{t}}$ such that its rotation vector is $(\frac pq,\frac rq)$ and, there exists a sequence $t_i>\overline{t},$ $t_i\rightarrow \overline{t},$ such that if $P_t$ is the continuation of $P_{\overline{t}}$ with the parameter, then $W^u(\widetilde{P}_{t_i})$ (the unstable manifold) has quadratic tangencies with $W^s(\widetilde{P}_{t_i})+(c,d)$ (the stable manifold translated by $(c,d)),$ where $\widetilde{P}_{t_i}$ is any lift of $P_{t_i}$ to the plane, in other words, $\widetilde{P}_{t_i}$ is a fixed point for $(\widetilde{f}_{t_i})^q-(p,r),$ and $(c,d)\neq (0,0)$ are certain integer vectors such that $W^u(\widetilde{P}_{\overline{t}})$ do not intersect $W^s(\widetilde{P}_{\overline{t}})+(c,d).$ And these tangencies become transverse as $t$ increases.


title: Learning to Emulate an Expert Projective Cone Scheduler

abstract: Projective cone scheduling defines a large class of rate-stabilizing policies for queueing models relevant to several applications. While there exists considerable theory on the properties of projective cone schedulers, there is little practical guidance on choosing the parameters that define them. In this paper, we propose an algorithm for designing an automated projective cone scheduling system based on observations of an expert projective cone scheduler. We show that the estimated scheduling policy is able to emulate the expert in the sense that the average loss realized by the learned policy will converge to zero. Specifically, for a system with $n$ queues observed over a time horizon $T$, the average loss for the algorithm is $O(\ln(T)\sqrt{\ln(n)/T})$. This upper bound holds regardless of the statistical characteristics of the system. The algorithm uses the multiplicative weights update method and can be applied online so that additional observations of the expert scheduler can be used to improve an existing estimate of the policy. This provides a data-driven method for designing a scheduling policy based on observations of a human expert. We demonstrate the efficacy of the algorithm with a simple numerical example and discuss several extensions.


title: $N^*$ Production from $e^+e^-$ Annihilations

abstract: Up to now, the $N^*$ production from $e^+e^-$ annihilations has been studied only around charmonium region. Charmonium decays to $N^*$s are analogous to (time-like) EM form factors in that the charm quark annihilation provides a nearly pointlike (ggg) current. Complementary to other sources, such as $\pi N$, $eN$ and $\gamma N$ reactions, this new source for $N^*$ spectroscopy has a few advantages, such as an isospin filter and a low spin filter. The experimental results on $N^*$ from $e^+e^-$ annihilations and their phenomenological implications are reviewed. Possible new sources on $N^*$ production from $e^+e^-$ annihilations are discussed.


title: Object Detection in Videos by High Quality Object Linking

abstract: Compared with object detection in static images, object detection in videos is more challenging due to degraded image qualities. An effective way to address this problem is to exploit temporal contexts by linking the same object across video to form tubelets and aggregating classification scores in the tubelets. In this paper, we focus on obtaining high quality object linking results for better classification. Unlike previous methods that link objects by checking boxes between neighboring frames, we propose to link in the same frame. To achieve this goal, we extend prior methods in following aspects: (1) a cuboid proposal network that extracts spatio-temporal candidate cuboids which bound the movement of objects; (2) a short tubelet detection network that detects short tubelets in short video segments; (3) a short tubelet linking algorithm that links temporally-overlapping short tubelets to form long tubelets. Experiments on the ImageNet VID dataset show that our method outperforms both the static image detector and the previous state of the art. In particular, our method improves results by 8.8% over the static image detector for fast moving objects.


title: Evidence for rotational to vibrational evolution along the yrast line in the odd-A rare-earth nuclei

abstract: The phase transition of nuclei to increasing angular momentum (or spin) and excitation energy is one of the most fundamental topics of nuclear structure research. The odd-N nuclei with A equal 160 are widely considered belonging to the well-deformed region, and their excitation spectra are energetically favored to exhibit the rotational characteristics. In the present work, however, there is evidence indicating that the nuclei can evolve from rotation to vibration along the yrast lines while increasing spin. The simple method, named as E-Gamma Over Spin (E-GOS) curves, would be used to discern the evolution from rotational to vibrational structure in nuclei as a function of spin. In addition, in order to get the insight into the rotational-like properties of nuclei, theoretical calculations have been performed for the yrast bands of the odd-A rare-earth nuclei using the total Routhian surfaces (TRS) model. The TRS plots indicate that the 165Yb and 157Dy isotopes have stable prolate shapes at low spin states. At higher rotational frequency (larger than 0.50 MeV), a distinct decrease in the quadrupole deformation is predicted by the calculations, and the isotopes becomes rigid in the gamma deformation.


title: On the global stability of departure time user equilibrium: A Lyapunov approach

abstract: In (Jin, 2018), a new day-to-day dynamical system was proposed for drivers' departure time choice at a single bottleneck. Based on three behavioral principles, the nonlocal departure and arrival times choice problems were converted to the local scheduling payoff choice problem, whose day-to-day dynamics are described by the Lighthill-Whitham-Richards (LWR) model on an imaginary road of increasing scheduling payoff. Thus the departure time user equilibrium (DTUE), the arrival time user equilibrium (ATUE), and the scheduling payoff user equilibrium (SPUE) are uniquely determined by the stationary state of the LWR model, which was shown to be locally, asymptotically stable with analysis of the discrete approximation of the LWR model and through a numerical example. In this study attempt to analytically prove the global stability of the SPUE, ATUE, and DTUE. We first generalize the conceptual models for arrival time and scheduling payoff choices developed in (Jin, 2018) for a single bottleneck with a generalized scheduling cost function, which includes the cost of the free-flow travel time. Then we present the LWR model for the day-to-day dynamics for the scheduling payoff choice as well as the SPUE. We further formulate a new optimization problem for the SPUE and demonstrate its equivalent to the optimization problem for the ATUE in (Iryo and Yoshii, 2007). Finally we show that the objective functions in the two optimization formulations are equal and can be used as the potential function for the LWR model and prove that the stationary state of the LWR model, and therefore, the SPUE, DTUE, and ATUE, are globally, asymptotically stable, by using Lyapunov's second method. Such a globally stable behavioral model can provide more efficient departure time and route choice guidance for human drivers and connected and autonomous vehicles in more complicated networks.


title: Manhattan Curves for Hyperbolic Surfaces with Cusps

abstract: In this paper, we study an interesting curve, so-called the Manhattan curve, associated with a pair of boundary-preserving Fuchsian representations of a (non-compact) surface, especially representations corresponding to Riemann surfaces with cusps. Using Thermodynamic Formalism (for countable Markov shifts), we prove the analyticity of the Manhattan curve. Moreover, we derive several dynamical and geometric rigidity results, which generalize results of Marc Burger and Richard Sharp for convex-cocompact Fuchsian representations.


title: Robustness of classification ability of spiking neural networks

abstract: It is well-known that the robustness of artificial neural networks (ANNs) is important for their wide ranges of applications. In this paper, we focus on the robustness of the classification ability of a spiking neural network which receives perturbed inputs. Actually, the perturbation is allowed to be arbitrary styles. However, Gaussian perturbation and other regular ones have been rarely investigated. For classification problems, the closer to the desired point, the more perturbed points there are in the input space. In addition, the perturbation may be periodic. Based on these facts, we only consider sinusoidal and Gaussian perturbations in this paper. With the SpikeProp algorithm, we perform extensive experiments on the classical XOR problem and other three benchmark datasets. The numerical results show that there is not significant reduction in the classification ability of the network if the input signals are subject to sinusoidal and Gaussian perturbations.


title: Regularity and continuity of the multilinear strong maximal operators

abstract: Let $m\ge 1$, in this paper, our object of investigation is the regularity and and continuity properties of the following multilinear strong maximal operator $${\mathscr{M}}_{\mathcal{R}}(\vec{f})(x)=\sup_{\substack{R \ni x R\in\mathcal{R}}}\prod\limits_{i=1}^m\frac{1}{|R|}\int_{R}|f_i(y)|dy,$$ where $x\in\mathbb{R}^d$ and $\mathcal{R}$ denotes the family of all rectangles in $\mathbb{R}^d$ with sides parallel to the axes. When $m=1$, denote $\mathscr{M}_{\mathcal{R}}$ by $\mathcal {M}_{\mathcal{R}}$.Then, $\mathcal {M}_{\mathcal{R}}$ coincides with the classical strong maximal function initially studied by Jessen, Marcinkiewicz and Zygmund. We showed that ${\mathscr{M}}_{\mathcal{R}}$ is bounded and continuous from the Sobolev spaces $W^{1,p_1}(\mathbb{R}^d)\times \cdots\times W^{1,p_m}(\mathbb{R}^d)$ to $W^{1,p} (\mathbb{R}^d)$, from the Besov spaces $B_{s}^{p_1,q} (\mathbb{R}^d)\times\cdots\times B_s^{p_m,q}(\mathbb{R}^d)$ to $B_s^{p,q}(\mathbb{R}^d)$, from the Triebel-Lizorkin spaces $F_{s}^{p_1,q}(\mathbb{R}^d)\times\cdots\times F_s^{p_m,q}(\mathbb{R}^d)$ to $F_s^{p,q}(\mathbb{R}^d)$. As a consequence, we further showed that ${\mathscr{M}}_{\mathcal{R}}$ is bounded and continuous from the fractional Sobolev spaces $W^{s,p_1}(\mathbb{R}^d)\times \cdots\times W^{s,p_m}(\mathbb{R}^d)$ to $W^{s,p}(\mathbb{R}^d)$ for $0<s\leq 1$ and $1<p<\infty$. As an application, we obtain a weak type inequality for the Sobolev capacity, which can be used to prove the $p$-quasicontinuity of $\mathscr{M}_{\mathcal{R}}$. The discrete type of the strong maximal operators has also been considered. We showed that this discrete type of the maximal operators enjoys somewhat unexpected regularity properties.


title: Weighted Community Detection and Data Clustering Using Message Passing

abstract: Grouping objects into clusters based on similarities or weights between them is one of the most important problems in science and engineering. In this work, by extending message passing algorithms and spectral algorithms proposed for unweighted community detection problem, we develop a non-parametric method based on statistical physics, by mapping the problem to Potts model at the critical temperature of spin glass transition and applying belief propagation to solve the marginals corresponding to the Boltzmann distribution. Our algorithm is robust to over-fitting and gives a principled way to determine whether there are significant clusters in the data and how many clusters there are. We apply our method to different clustering tasks and use extensive numerical experiments to illustrate the advantage of our method over existing algorithms. In the community detection problem in weighted and directed networks, we show that our algorithm significantly outperforms existing algorithms. In the clustering problem when the data was generated by mixture models in the sparse regime we show that our method works to the theoretical limit of detectability and gives accuracy very close to that of the optimal Bayesian inference. In the semi-supervised clustering problem, our method only needs several labels to work perfectly in classic datasets. Finally, we further develop Thouless-Anderson-Palmer equations which reduce heavily the computation complexity in dense-networks but gives almost the same performance as belief propagation.


title: Linking Phase Transitions and Quantum Entanglement at Arbitrary Temperature

abstract: In this work, we establish a general theory of phase transitions and quantum entanglement in the equilibrium state at arbitrary temperatures. First, we derived a set of universal functional relations between the matrix elements of two-body reduced density matrix of the canonical density matrix and the Helmholtz free energy of the equilibrium state, which implies that the Helmholtz free energy and its derivatives are directly related to entanglement measures because any entanglement measures are defined as a function of the reduced density matrix. Then we show that the first order phase transitions are signaled by the matrix elements of reduced density matrix while the second order phase transitions are witnessed by the first derivatives of the reduced density matrix elements. Near second order phase transition point, we show that the first derivative of the reduced density matrix elements present universal scaling behaviors. Finally we establish a theorem which connects the phase transitions and entanglement at arbitrary temperatures. Our general results are demonstrated in an experimentally relevant many-body spin model.


title: Evaluation of Parallel Tempering to Accelerate Bayesian Parameter Estimation in Systems Biology

abstract: Models of biological systems often have many unknown parameters that must be determined in order for model behavior to match experimental observations. Commonly-used methods for parameter estimation that return point estimates of the best-fit parameters are insufficient when models are high dimensional and under-constrained. As a result, Bayesian methods, which treat model parameters as random variables and attempt to estimate their probability distributions given data, have become popular in systems biology. Bayesian parameter estimation often relies on Markov Chain Monte Carlo (MCMC) methods to sample model parameter distributions, but the slow convergence of MCMC sampling can be a major bottleneck. One approach to improving performance is parallel tempering (PT), a physics-based method that uses swapping between multiple Markov chains run in parallel at different temperatures to accelerate sampling. The temperature of a Markov chain determines the probability of accepting an unfavorable move, so swapping with higher temperatures chains enables the sampling chain to escape from local minima. In this work we compared the MCMC performance of PT and the commonly-used Metropolis-Hastings (MH) algorithm on six biological models of varying complexity. We found that for simpler models PT accelerated convergence and sampling, and that for more complex models, PT often converged in cases MH became trapped in non-optimal local minima. We also developed a freely-available MATLAB package for Bayesian parameter estimation called PTempEst (http://github.com/RuleWorld/ptempest), which is closely integrated with the popular BioNetGen software for rule-based modeling of biological systems.


title: Characterizations for inner functions in certain function spaces

abstract: For $\frac12<p<\infty$, $0<q<\infty$ and a certain two-sided doubling weight $\omega$, we characterize those inner functions $\Theta$ for which $$\|\Theta'\|_{A^{p,q}_\omega}^q=\int_0^1 \left(\int_0^{2\pi} |\Theta'(re^{i\theta})|^p d\theta\right)^{q/p} \omega(r)\,dr<\infty.$$ Then we show a modified version of this result for $p\ge q$. Moreover, two additional characterizations for inner functions whose derivative belongs to the Bergman space $A_\omega^{p,p}$ are given.


title: Strain engineering of the silicon-vacancy center in diamond

abstract: We control the electronic structure of the silicon-vacancy (SiV) color-center in diamond by changing its static strain environment with a nano-electro-mechanical system. This allows deterministic and local tuning of SiV optical and spin transition frequencies over a wide range, an essential step towards multi-qubit networks. In the process, we infer the strain Hamiltonian of the SiV revealing large strain susceptibilities of order 1 PHz/strain for the electronic orbital states. We identify regimes where the spin-orbit interaction results in a large strain suseptibility of order 100 THz/strain for spin transitions, and propose an experiment where the SiV spin is strongly coupled to a nanomechanical resonator.


title: A Flexible Procedure for Mixture Proportion Estimation in Positive-Unlabeled Learning

abstract: Positive--unlabeled (PU) learning considers two samples, a positive set P with observations from only one class and an unlabeled set U with observations from two classes. The goal is to classify observations in U. Class mixture proportion estimation (MPE) in U is a key step in PU learning. Blanchard et al. [2010] showed that MPE in PU learning is a generalization of the problem of estimating the proportion of true null hypotheses in multiple testing problems. Motivated by this idea, we propose reducing the problem to one dimension via construction of a probabilistic classifier trained on the P and U data sets followed by application of a one--dimensional mixture proportion method from the multiple testing literature to the observation class probabilities. The flexibility of this framework lies in the freedom to choose the classifier and the one--dimensional MPE method. We prove consistency of two mixture proportion estimators using bounds from empirical process theory, develop tuning parameter free implementations, and demonstrate that they have competitive performance on simulated waveform data and a protein signaling problem.


title: A study of deformation localization in nonlinear elastic lattices

abstract: The paper investigates localized deformation patterns resulting from the onset of instabilities in lattice structures. The study is motivated by previous observations on discrete hexagonal lattices, where the onset of non-uniform, quasi-static deformation patterns was associated with the loss of convexity of the interaction potential, and where a variety of localized deformations were found depending on loading configuration, lattice parameters and boundary conditions. These observations are here conducted on other lattice structures, with the goal of identifying models of reduced complexity that are able to provide insight into the key parameters that govern the onset of instability-induced localization. To this end, we first consider a two-dimensional square lattice consisting of point masses connected by in-plane axial springs and vertical ground springs. Results illustrate that depending on the choice of spring constants and their relative values, the lattice exhibits in-plane or out-of plane instabilities leading to folding and unfolding. This model is further simplified by considering the one-dimensional case of a spring-mass chain sitting on an elastic foundation. A bifurcation analysis of this lattice identifies the stable and unstable branches and illustrates its hysteretic and loading path-dependent behaviors. Finally, the lattice is further reduced to a minimal four mass model which undergoes a folding/unfolding process qualitatively similar to the same process in the central part of a longer chain, helping our understanding of localization in more complex systems. In contrast to the widespread assumption that localization is induced by defects or imperfections in a structure, this work illustrates that such phenomena can arise in perfect lattices as a consequence of the mode-shapes at the bifurcation points.


title: On conormal and oblique derivative problem for elliptic equations with Dini mean oscillation coefficients

abstract: We show that weak solutions to conormal derivative problem for elliptic equations in divergence form are continuously differentiable up to the boundary provided that the mean oscillations of the leading coefficients satisfy the Dini condition, the lower order coefficients satisfy certain suitable conditions, and the boundary is locally represented by a $C^1$ function whose derivatives are Dini continuous. We also prove that strong solutions to oblique derivative problem for elliptic equations in nondivergence form are twice continuously differentiable up to the boundary if the mean oscillations of coefficients satisfy the Dini condition and the boundary is locally represented by a $C^1$ function whose derivatives are double Dini continuous. This in particular extends a result of M. V. Safonov (Comm. Partial Differential Equations 20:1349--1367, 1995)


title: Splitting Localization and Prediction Numbers

abstract: In this paper the work done by Newelski and Roslanowski is revisited to solve a question done by Blass about one of the possible evasion and prediction numbers. This led to define a variation of the $k$-localization property (the $(k+1)^{\omega}$-localization property) and the use of a forcing notion with accelerating trees.


title: Multiple Accounts Detection on Facebook Using Semi-Supervised Learning on Graphs

abstract: In social networks, a single user may create multiple accounts to spread his / her opinions and to influence others, by actively comment on different news pages. It would be beneficial to both social networks and their communities, to demote such abnormal activities, and the first step is to detect those accounts. However, the detection is challenging, because these accounts may have very realistic names and reasonable activity patterns. In this paper, we investigate three different approaches, and propose using graph embedding together with semi-supervised learning, to predict whether a pair of accounts are created by the same user. We carry out extensive experimental analyses to understand how changes in the input data and algorithmic parameters / optimization affect the prediction performance. We also discover that local information have higher importance than the global ones for such prediction, and point out the threshold leading to the best results. We test the proposed approach with 6700 Facebook pages from the Middle East, and achieve the averaged accuracy at 0.996 and AUC (area under curve) at 0.952 for users with the same name; with the U.S. 2016 election dataset, we obtain the best AUC at 0.877 for users with different names.


title: Yet another UFO in the X-ray spectrum of a high-z lensed QSO

abstract: Ultra-fast outflows (UFO) appear to be common in local active galactic nuclei (AGN) and may be powerful enough ($\dot{E}_{kin}$$\geq$1\% of L$_{bol}$) to effectively quench the star formation in their host galaxies. To test feedback models based on AGN outflows, it is mandatory to investigate UFOs near the peak of AGN activity, that is, at high-z where only a few studies are available to date. UFOs produce Fe resonant absorption lines measured above $\approx$7 keV. The most critical problem in detecting such features in distant objects is the difficulty in obtaining X-ray data with sufficient signal-to-noise. We therefore selected a distant QSO that gravitational lensing made bright enough for these purposes, the z=2.64 QSO MG J0414+0534, and observed it with XMM-Newton for $\approx$78 ks.} The X-ray spectrum of MG J0414+0534 is complex and shows signatures of cold absorption (N$_{H}\approx$4$\times$10$^{22}$ cm$^{-2}$) and of the presence of an iron emission line (E$\approx$6.4 keV, EW$=$95$\pm$53 eV) consistent with it originating in the cold absorber. Our main result, however, is the robust detection (more than 5$\sigma$) of an absorption line at E$_{int}\approx$9.2 keV (E$_{obs}\approx$2.5 keV observer frame). If interpreted as due to FeXXVI, it implies gas outflowing at $v_{out}\approx$0.3c. To our knowledge, this is the first detection of an UFO in a radio-loud quasar at z$\geq$1.5. We estimated that the UFO mechanical output is $\dot{E}_{kin}$$\approx$2.5$L_{bol}$ with $\dot{p}_{out}/\dot{p}_{rad}\approx$17 indicating that it is capable of installing significant feedback between the super-massive black hole (SMBH) and the bulge of the host galaxy. We argue that this also suggests a magnetic driving origin of the UFO.


title: A time-optimal algorithm for solving (block-)tridiagonal linear systems of dimension N on a distributed computer of N nodes

abstract: We are concerned with the fastest possible direct numerical solution algorithm for a thin-banded or tridiagonal linear system of dimension $N$ on a distributed computing network of $N$ nodes that is connected in a binary communication tree. Our research is driven by the need for faster ways of numerically solving discretized systems of coupled one-dimensional black-box boundary-value problems. Our paper presents two major results: First, we provide an algorithm that achieves the optimal parallel time complexity for solving a tridiagonal linear system and thin-banded linear systems. Second, we prove that it is impossible to improve the time complexity of this method by any polynomial degree. To solve a system of dimension $m\cdot N$ and bandwidth $m \in \Omega(N^{1/6})$ on $2 \cdot N-1$ computing nodes, our method needs time complexity $\mathcal{O}(\log(N)^2 \cdot m^3)$.


title: Olivine-Carbonate Mineralogy of Jezero Crater

abstract: A well-preserved, ancient delta deposit, in combination with ample exposures of Mg-carbonate rich materials, make Jezero Crater in Nili Fossae a compelling astrobiological site and a top candidate for future landed missions to Mars. We use CRISM observations to characterize the surface mineralogy of the crater and surrounding watershed. We have identified a three- endmember sequence of olivine-bearing lithologies that we hypothesize are distinguished by their Mg content. We find that Mg-carbonates are consistently identified in association with one of the olivine-bearing lithologies, although that lithology is not fully carbonatized. Surprisingly, this lithology contains relatively Fe-rich olivine. We address a range of formation scenarios, including the possibility that these olivine and carbonate associations are indicators of serpentinization on early Mars. These deposits provide an opportunity for deepening our understanding of early Mars by revealing the thermal history of the martian interior and potentially changes in its tectonic regime with time.


title: Fu-Yau Hessian Equations

abstract: We solve the Fu-Yau equation for arbitrary dimension and arbitrary slope $\alpha'$. Actually we obtain at the same time a solution of the open case $\alpha'>0$, an improved solution of the known case $\alpha'<0$, and solutions for a family of Hessian equations which includes the Fu-Yau equation as a special case. The method is based on the introduction of a more stringent ellipticity condition than the usual $\Gamma_k$ admissible cone condition, and which can be shown to be preserved by precise estimates with scale.


title: Well-posedness of a non-local abstract Cauchy problem with a singular integral

abstract: A non-local abstract Cauchy problem with a singular integral is studied, which is a closed system of two evolution equations for a real-valued function and a function-valued function. By proposing an appropriate Banach space, the well-posedness of the evolution system is proved under some boundedness and smoothness conditions on the coefficient functions. Furthermore, an isomorphism is established to extend the result to a partial integro-differential equation with singular convolution kernel, which is a generalized form of the stationary Wigner equation. Our investigation considerably improves the understanding of the open problem concerning the well-posedness of the stationary Wigner equation with inflow boundary conditions.


title: Bi-scalar integrable CFT at any dimension

abstract: We propose a $D$-dimensional generalization of $4D$ bi-scalar conformal quantum field theory recently introduced by G\"{u}rdogan and one of the authors as a strong-twist double scaling limit of $\gamma$-deformed $\mathcal{N}=4$ SYM theory. Similarly to the $4D$ case, this D-dimensional CFT is also dominated by "fishnet" Feynman graphs and is integrable in the planar limit. The dynamics of these graphs is described by the integrable conformal $SO(D+1,1)$ spin chain. In $2D$ it is the analogue of L. Lipatov's $SL(2,\mathbb{C})$ spin chain for the Regge limit of $QCD$, but with the spins $s=1/4$ instead of $s=0$. Generalizing recent $4D$ results of Grabner, Gromov, Korchemsky and one of the authors to any $D$ we compute exactly, at any coupling, a four point correlation function, dominated by the simplest fishnet graphs of cylindric topology, and extract from it exact dimensions of R-charge 2 operators with any spin and some of their OPE structure constants.


title: Bayesian Model Testing of Ellipsoidal Variations on Stars due to Hot Jupiters

abstract: A massive planet closely orbiting its host star creates tidal forces that distort the typically spherical stellar surface. These distortions, known as ellipsoidal variations, result in changes in the photometric flux emitted by the star, which can be detected within the data from the Kepler Space Telescope. Currently, there exist several models describing such variations and their effect on the photometric flux. By using Bayesian model testing in conjunction with the Bayesian-based exoplanet characterization software package EXONEST, the most probable representation for ellipsoidal variations was determined for synthetic data and the confirmed hot Jupiter exoplanet Kepler-13Ab.The most preferred model for ellipsoidal variations observed in the Kepler-13 light curve was determined to be EVIL-MC. Among the trigonometric models, the Modified Kane & Gelino model provided the best representation of ellipsoidal variations for the Kepler-13 system and may serve as a fast alternative to the more computationally intensive EVIL-MC. The computational feasibility of directly modeling the ellipsoidal variations of a star are examined and future work is outlined. Providing a more accurate model of ellipsoidal variations is expected to result in better planetary mass estimations.


title: Antenna Selection for Large-Scale MIMO Systems with Low-Resolution ADCs

abstract: One way to reduce the power consumption in large-scale multiple-input multiple-output (MIMO) systems is to employ low-resolution analog-to-digital converters (ADCs). In this paper, we investigate antenna selection for large-scale MIMO receivers with low-resolution ADCs, thereby providing more flexibility in resolution and number of ADCs. To incorporate quantization effects, we generalize an existing objective function for a greedy capacity-maximization antenna selection approach. The derived objective function offers an opportunity to select an antenna with the best tradeoff between the additional channel gain and increase in quantization error. Using the generalized objective function, we propose an antenna selection algorithm based on a conventional antenna selection algorithm without an increase in overall complexity. Simulation results show that the proposed algorithm outperforms the conventional algorithm in achievable capacity for the same number of antennas.


title: Open3D: A Modern Library for 3D Data Processing

abstract: Open3D is an open-source library that supports rapid development of software that deals with 3D data. The Open3D frontend exposes a set of carefully selected data structures and algorithms in both C++ and Python. The backend is highly optimized and is set up for parallelization. Open3D was developed from a clean slate with a small and carefully considered set of dependencies. It can be set up on different platforms and compiled from source with minimal effort. The code is clean, consistently styled, and maintained via a clear code review mechanism. Open3D has been used in a number of published research projects and is actively deployed in the cloud. We welcome contributions from the open-source community.


title: Over-representation of Extreme Events in Decision-Making: A Rational Metacognitive Account

abstract: The Availability bias, manifested in the over-representation of extreme eventualities in decision-making, is a well-known cognitive bias, and is generally taken as evidence of human irrationality. In this work, we present the first rational, metacognitive account of the Availability bias, formally articulated at Marr's algorithmic level of analysis. Concretely, we present a normative, metacognitive model of how a cognitive system should over-represent extreme eventualities, depending on the amount of time available at its disposal for decision-making. Our model also accounts for two well-known framing effects in human decision-making under risk---the fourfold pattern of risk preferences in outcome probability (Tversky & Kahneman, 1992) and in outcome magnitude (Markovitz, 1952)---thereby providing the first metacognitively-rational basis for those effects. Empirical evidence, furthermore, confirms an important prediction of our model. Surprisingly, our model is unimaginably robust with respect to its focal parameter. We discuss the implications of our work for studies on human decision-making, and conclude by presenting a counterintuitive prediction of our model, which, if confirmed, would have intriguing implications for human decision-making under risk. To our knowledge, our model is the first metacognitive, resource-rational process model of cognitive biases in decision-making.


title: A characterization of nonnegativity relative to proper cones

abstract: Let $A$ be an $m \times n$ matrix with real entries. Given two proper cones $K_1$ and $K_2$ in $\mathbb{R}^n$ and $\mathbb{R}^m$, respectively, we say that $A$ is nonnegative if $A(K_1) \subseteq K_2$. $A$ is said to be semipositive if there exists a $x \in K_1^\circ$ such that $Ax \in K_2^\circ$. We prove that $A$ is nonnegative if and only if $A+B$ is semipositive for every semipositive matrix $B$. Applications of the above result are also brought out.


title: Sparsity in Max-Plus Algebra and Systems

abstract: We study sparsity in the max-plus algebraic setting. We seek both exact and approximate solutions of the max-plus linear equation with minimum cardinality of support. In the former case, the sparsest solution problem is shown to be equivalent to the minimum set cover problem and, thus, NP-complete. In the latter one, the approximation is quantified by the $\ell_{1}$ residual error norm, which is shown to have supermodular properties under some convex constraints, called lateness constraints. Thus, greedy approximation algorithms of polynomial complexity can be employed for both problems with guaranteed bounds of approximation. We also study the sparse recovery problem and present conditions, under which, the sparsest exact solution solves it. Through multi-machine interactive processes, we describe how the present framework could be applied to two practical discrete event systems problems: resource optimization and structure-seeking system identification. We also show how sparsity is related to the pruning problem. Finally, we present a numerical example of the structure-seeking system identification problem and we study the performance of the greedy algorithm via simulations.


title: Cross-type Biomedical Named Entity Recognition with Deep Multi-Task Learning

abstract: Motivation: State-of-the-art biomedical named entity recognition (BioNER) systems often require handcrafted features specific to each entity type, such as genes, chemicals and diseases. Although recent studies explored using neural network models for BioNER to free experts from manual feature engineering, the performance remains limited by the available training data for each entity type. Results: We propose a multi-task learning framework for BioNER to collectively use the training data of different types of entities and improve the performance on each of them. In experiments on 15 benchmark BioNER datasets, our multi-task model achieves substantially better performance compared with state-of-the-art BioNER systems and baseline neural sequence labeling models. Further analysis shows that the large performance gains come from sharing character- and word-level information among relevant biomedical entities across differently labeled corpora.


title: Big polynomial rings and Stillman's conjecture

abstract: The purpose of this paper is to prove that certain limits of polynomial rings are themselves polynomial rings, and show how this observation can be used to deduce some interesting results in commutative algebra. In particular, we give two new proofs of Stillman's conjecture. The first is similar to that of Ananyan-Hochster, though more streamlined; in particular, it establishes the existence of small subalgebras. The second proof is completely different, and relies on a recent noetherianity result of Draisma.


title: Spin portal to dark matter

abstract: In this work we study the possibility that dark matter fields transform in the $(1,0)\oplus(0,1)$ representation of the Homogeneous Lorentz Group. In an effective theory approach, we study the lowest dimension interacting terms of dark matter with standard model fields, assuming that dark matter fields transform as singlets under the standard model gauge group. There are three dimension-four operators, two of them yielding a Higgs portal to dark matter. The third operator couple the photon and $Z^0$ fields to the higher multipoles of dark matter, yielding a \textit{spin portal} to dark matter. For dark matter ($D$) mass below a half of the $Z^0$ mass, the decays $Z^0\to \bar{D} D$ and $H\to \bar{D}D$ are kinematically allowed and contribute to the invisible widths of the $Z^0$ and $H$. We calculate these decays and use experimental results on these invisible widths to constrain the values of the low energy constants finding in general that effects of the spin portal can be more important that those of the Higgs portal. We calculate the dark matter relic density in our formalism, use the constraints on the low energy constants from the $Z^0$ and $H$ invisible widths and compare our results with the measured relic density, finding that dark matter with a $(1,0)\oplus(0,1)$ space-time structure must have a mass $M>43 ~ GeV$.


title: Algorithms for the Greater Good! On Mental Modeling and Acceptable Symbiosis in Human-AI Collaboration

abstract: Effective collaboration between humans and AI-based systems requires effective modeling of the human in the loop, both in terms of the mental state as well as the physical capabilities of the latter. However, these models can also open up pathways for manipulating and exploiting the human in the hopes of achieving some greater good, especially when the intent or values of the AI and the human are not aligned or when they have an asymmetrical relationship with respect to knowledge or computation power. In fact, such behavior does not necessarily require any malicious intent but can rather be borne out of cooperative scenarios. It is also beyond simple misinterpretation of intents, as in the case of value alignment problems, and thus can be effectively engineered if desired. Such techniques already exist and pose several unresolved ethical and moral questions with regards to the design of autonomy. In this paper, we illustrate some of these issues in a teaming scenario and investigate how they are perceived by participants in a thought experiment.


title: Stop Coannihilation in the CMSSM and SubGUT Models

abstract: Stop coannihilation may bring the relic density of heavy supersymmetric dark matter particles into the range allowed by cosmology. The efficiency of this process is enhanced by stop-antistop annihilations into the longitudinal (Goldstone) modes of the $W$ and $Z$ bosons, as well as by Sommerfeld enhancement of stop annihilations and the effects of bound states. Since the couplings of the stops to the Goldstone modes are proportional to the trilinear soft supersymmetry-breaking $A$-terms, these annihilations are enhanced when the $A$-terms are large. However, the Higgs mass may be reduced below the measured value if the $A$-terms are too large. Unfortunately, the interpretation of this constraint on the stop coannihilation strip is clouded by differences between the available Higgs mass calculators. For our study, we use as our default calculator FeynHiggs 2.13.0, the most recent publicly available version of this code. Exploring the CMSSM parameter space, we find that along the stop coannihilation strip the masses of the stops are severely split by the large $A$-terms. This suppresses the Higgs mass drastically for $\mu$ and $A_0 > 0$, whilst the extent of the stop coannihilation strip is limited for $A_0 < 0$ and either sign of $\mu$. However, in sub-GUT models, reduced renormalization-group running mitigates the effect of the large $A$-terms, allowing larger LSP masses to be consistent with the Higgs mass calculation. We give examples where the dark matter particle mass may reach $\gtrsim 8$ TeV.


title: ReNN: Rule-embedded Neural Networks

abstract: The artificial neural network shows powerful ability of inference, but it is still criticized for lack of interpretability and prerequisite needs of big dataset. This paper proposes the Rule-embedded Neural Network (ReNN) to overcome the shortages. ReNN first makes local-based inferences to detect local patterns, and then uses rules based on domain knowledge about the local patterns to generate rule-modulated map. After that, ReNN makes global-based inferences that synthesizes the local patterns and the rule-modulated map. To solve the optimization problem caused by rules, we use a two-stage optimization strategy to train the ReNN model. By introducing rules into ReNN, we can strengthen traditional neural networks with long-term dependencies which are difficult to learn with limited empirical dataset, thus improving inference accuracy. The complexity of neural networks can be reduced since long-term dependencies are not modeled with neural connections, and thus the amount of data needed to optimize the neural networks can be reduced. Besides, inferences from ReNN can be analyzed with both local patterns and rules, and thus have better interpretability. In this paper, ReNN has been validated with a time-series detection problem.


title: On a $C$-integrable equation for second sound propagation in heated dielectrics

abstract: An exactly solvable model in heat conduction is considered. The $C$-integrable (i.e., change-of-variables-integrable) equation for second sound (i.e., heat wave) propagation in a thin, rigid dielectric heat conductor uniformly heated on its lateral side by a surrounding medium under the Stefan--Boltzmann law is derived. A simple change-of-variables transformation is shown to exactly map the nonlinear governing partial differential equation to the classical linear telegrapher's equation. In a one-dimensional context, known integral-transform solutions of the latter are adapted to construct exact solutions relevant to heat transfer applications: (i) the initial-value problem on an infinite domain (the real line), and (ii) the initial-boundary-value problem on a semi-infinite domain (the half-line). Possible "second law violations" and restrictions on the $C$-transformation are noted for some sets of parameters.


title: Decoding and mapping task states of the human brain via deep learning

abstract: Support vector machine (SVM) based multivariate pattern analysis (MVPA) has delivered promising performance in decoding specific task states based on functional magnetic resonance imaging (fMRI) of the human brain. Conventionally, the SVM-MVPA requires careful feature selection/extraction according to expert knowledge. In this study, we propose a deep neural network (DNN) for directly decoding multiple brain task states from fMRI signals of the brain without any burden for feature handcrafts. We trained and tested the DNN classifier using task fMRI data from the Human Connectome Project's S1200 dataset (N=1034). In tests to verify its performance, the proposed classification method identified seven tasks with an average accuracy of 93.7%. We also showed the general applicability of the DNN for transfer learning to small datasets (N=43), a situation encountered in typical neuroscience research. The proposed method achieved an average accuracy of 89.0% and 94.7% on a working memory task and a motor classification task, respectively, higher than the accuracy of 69.2% and 68.6% obtained by the SVM-MVPA. A network visualization analysis showed that the DNN automatically detected features from areas of the brain related to each task. Without incurring the burden of handcrafting the features, the proposed deep decoding method can classify brain task states highly accurately, and is a powerful tool for fMRI researchers.


title: Structured Memory based Deep Model to Detect as well as Characterize Novel Inputs

abstract: While deep learning has pushed the boundaries in various machine learning tasks, the current models are still far away from replicating many functions that a normal human brain can do. Explicit memorization based deep architecture have been recently proposed with the objective to understand and predict better. In this work, we design a system that involves a primary learner and an adjacent representational memory bank which is organized using a comparative learner. This spatially forked deep architecture with a structured memory can simultaneously predict and reason about the nature of an input, which may even belong to a category never seen in the training data, by relating it with the memorized past representations at the higher layers. Characterizing images of unseen object classes in both synthetic and real world datasets is used as an example to showcase the operational success of the proposed framework.


title: Quasiperiodic granular chains and Hofstadter butterflies

abstract: We study quasiperiodicity-induced localization of waves in strongly precompressed granular chains. We propose three different setups, inspired by the Aubry--Andr\'e (AA) model, of quasiperiodic chains; and we use these models to compare the effects of on-site and off-site quasiperiodicity in nonlinear lattices. When there is purely on-site quasiperiodicity, which we implement in two different ways, we show for a chain of spherical particles that there is a localization transition (as in the original AA model). However, we observe no localization transition in a chain of cylindrical particles in which we incorporate quasiperiodicity in the distribution of contact angles between adjacent cylinders by making the angle periodicity incommensurate with that of the chain. For each of our three models, we compute the Hofstadter spectrum and the associated Minkowski--Bouligand fractal dimension, and we demonstrate that the fractal dimension decreases as one approaches the localization transition (when it exists). Finally, in a suite of numerical computations, we demonstrate both localization and that there exist regimes of ballistic, superdiffusive, diffusive, and subdiffusive transport. Our models provide a flexible set of systems to study quasiperiodicity-induced analogs of Anderson phenomena in granular chains that one can tune controllably from weakly to strongly nonlinear regimes.

categories: nlin.PS cond-mat.mtrl-sci math-ph math.DS math.MP nlin.CD


title: Self-consistent Analytic Solutions in Twisted $\mathbb{C}P^{N-1}$ Model in the Large-$N$ Limit

abstract: We construct self-consistent analytic solutions in the ${\mathbb C}P^{N-1}$ model in the large-$N$ limit, in which more than one Higgs scalar component take values inside a single or multiple soliton on an infinite space or on a ring, or around boundaries of a finite interval.


title: Electrostatic steepening of whistler waves

abstract: We present surprising observations by the NASA Van Allen Probes spacecraft of whistler waves with substantial electric field power at harmonics of the whistler wave fundamental frequency. The wave power at harmonics is due to nonlinearly steepened whistler electrostatic field that becomes possible in the two-temperature electron plasma due to whistler wave coupling to the electron-acoustic mode. The simulation and analytical estimates show that the steepening takes a few tens of milliseconds. The hydrodynamic energy cascade to higher frequencies facilitates efficient energy transfer from cyclotron resonant electrons, driving the whistler waves, to lower energy electrons.


title: Burnside groups and $n$-moves for links

abstract: Let $n$ be a positive integer. M. K. Dabkowski and J. H. Przytycki introduced the $n$th Burnside group of links which is preserved by $n$-moves, and proved that for any odd prime $p$ there exist links which are not equivalent to trivial links up to $p$-moves by using their $p$th Burnside groups. This gives counterexamples for the Montesinos-Nakanishi $3$-move conjecture. In general, it is hard to distinguish $p$th Burnside groups of a given link and a trivial link. We give a necessary condition for which $p$th Burnside groups are isomorphic to those of trivial links. The necessary condition gives us an efficient way to distinguish $p$th Burnside groups of a given link and a trivial link. As an application, we show that there exist links, each of which is not equivalent to a trivial link up to $p$-moves for any odd prime $p$.


title: NGS Based Haplotype Assembly Using Matrix Completion

abstract: We apply matrix completion methods for haplotype assembly from NGS reads to develop the new HapSVT, HapNuc, and HapOPT algorithms. This is performed by applying a mathematical model to convert the reads to an incomplete matrix and estimating unknown components. This process is followed by quantizing and decoding the completed matrix in order to estimate haplotypes. These algorithms are compared to the state-of-the-art algorithms using simulated data as well as the real fosmid data. It is shown that the SNP missing rate and the haplotype block length of the proposed HapOPT are better than those of HapCUT2 with comparable accuracy in terms of reconstruction rate and switch error rate. A program implementing the proposed algorithms in MATLAB is freely available at https://github.com/smajidian/HapMC.


title: Magnetoelectric and Raman spectroscopic studies of single-crystalline MnCr2O4

abstract: MnCr2O4 that exhibits spin frustration and complex spiral spin order is of great interest from both fundamental as well as application-oriented perspectives. Unlike CoCr2O4 whose ground state presents the coexistence of commensurate spiral spin order (CSSO) and ferroelectric order, MnCr2O4 shows no multiferroicity. One reason is that the spiral spin order is highly sensitive to the oxygen concentration in MnCr2O4. Here, we have successfully grown high-quality single-crystalline MnCr2O4 by the chemical vapor transport method. We observe a new first-order magnetic transition from the incommensurate spiral spin order (ICSSO) at 19.4 K to the CSSO at 17.4 K. This magnetic transition is verified by magnetization, specific heat, and magnetoelectric measurements, which also confirm that the ground state exhibits the coexistence of the CSSO and magnetoelectricity below 17.4 K. Interestingly, the temperature evolution of Raman spectra between 5.4 and 300 K suggests that the structure remains the same. We also find that the phase-transition temperature of the CSSO decreases as applied magnetic field increases up to 45 kOe.


title: Accelerating recurrent neural network language model based online speech recognition system

abstract: This paper presents methods to accelerate recurrent neural network based language models (RNNLMs) for online speech recognition systems. Firstly, a lossy compression of the past hidden layer outputs (history vector) with caching is introduced in order to reduce the number of LM queries. Next, RNNLM computations are deployed in a CPU-GPU hybrid manner, which computes each layer of the model on a more advantageous platform. The added overhead by data exchanges between CPU and GPU is compensated through a frame-wise batching strategy. The performance of the proposed methods evaluated on LibriSpeech test sets indicates that the reduction in history vector precision improves the average recognition speed by 1.23 times with minimum degradation in accuracy. On the other hand, the CPU-GPU hybrid parallelization enables RNNLM based real-time recognition with a four times improvement in speed.


title: An iterative support shrinking algorithm for $\ell_{p}$-$\ell_{q}$ minimization

abstract: We present an iterative support shrinking algorithm for $\ell_{p}$-$\ell_{q}$ minimization~($0 <p < 1 \leq q < \infty $). This algorithm guarantees the nonexpensiveness of the signal support set and can be easily implemented after being proximally linearized. The subproblem can be very efficiently solved due to its convexity and reducing size along iteration. We prove that the iterates of the algorithm globally converge to a stationary point of the $\ell_{p}$-$\ell_{q}$ objective function. In addition, we show a lower bound theory for the iteration sequence, which is more practical than the lower bound results for local minimizers in the literature.


title: New characterizations of freeness for hyperplane arrangements

abstract: In this article we describe two new characterizations of freeness for hyperplane arrangements via the study of the generic initial ideal and of the sectional matrix of the Jacobian ideal of arrangements.


title: Assessment of interaction-strength interpolation formulas for gold and silver clusters

abstract: The performance of functionals based on the idea of interpolating between the weak and the strong-interaction limits the global adiabatic-connection integrand is carefully studied for the challenging case of noble-metal clusters. Different interpolation formulas are considered and various features of this approach are analyzed. It is found that these functionals, when used as a correlation correction to Hartree-Fock, are quite robust for the description of atomization energies, while performing less well for ionization potentials. Future directions that can be envisaged from this study and a previous one on main group chemistry are discussed.


title: Fast Power system security analysis with Guided Dropout

abstract: We propose a new method to efficiently compute load-flows (the steady-state of the power-grid for given productions, consumptions and grid topology), substituting conventional simulators based on differential equation solvers. We use a deep feed-forward neural network trained with load-flows precomputed by simulation. Our architecture permits to train a network on so-called "n-1" problems, in which load flows are evaluated for every possible line disconnection, then generalize to "n-2" problems without retraining (a clear advantage because of the combinatorial nature of the problem). To that end, we developed a technique bearing similarity with "dropout", which we named "guided dropout".


title: An infinite family of subcubic graphs with unbounded packing chromatic number

abstract: Recently, Balogh, Kostochka and Liu in [Packing chromatic number of cubic graphs, Discrete Math.~341 (2018) 474--483] answered in negative the question that was posed in several earlier papers whether the packing chromatic number is bounded in the class of graphs with maximum degree $3$. In this note, we present an explicit infinite family of subcubic graphs with unbounded packing chromatic number.


title: A State-of-the-Art of Semantic Change Computation

abstract: This paper reviews the state-of-the-art of semantic change computation, one emerging research field in computational linguistics, proposing a framework that summarizes the literature by identifying and expounding five essential components in the field: diachronic corpus, diachronic word sense characterization, change modelling, evaluation data and data visualization. Despite the potential of the field, the review shows that current studies are mainly focused on testifying hypotheses proposed in theoretical linguistics and that several core issues remain to be solved: the need for diachronic corpora of languages other than English, the need for comprehensive evaluation data for evaluation, the comparison and construction of approaches to diachronic word sense characterization and change modelling, and further exploration of data visualization techniques for hypothesis justification.


title: Condition for confinement in non-Abelian gauge theories

abstract: We show that a criterion for confinement, based on the BRST invariance, holds in four dimensions, by solving a non-Abelian gauge theory with a set of exact solutions. The confinement condition we consider was obtained by Kugo and Ojima some decades ago. The current understanding of gauge theories permits us to apply the techniques straightforwardly for checking the validity of this criterion. In this way, we are able to show that the non-Abelian gauge theory is confining and that confinement is rooted in the BRST invariance and asymptotic freedom.


title: Change point analysis in non-stationary processes - a mass excess approach

abstract: This paper considers the problem of testing if a sequence of means $(\mu_t)_{t =1,\ldots ,n }$ of a non-stationary time series $(X_t)_{t =1,\ldots ,n }$ is stable in the sense that the difference of the means $\mu_1$ and $\mu_t$ between the initial time $t=1$ and any other time is smaller than a given level, that is $ | \mu_1 - \mu_t | \leq c $ for all $t =1,\ldots ,n $. A test for hypotheses of this type is developed using a biascorrected monotone rearranged local linear estimator and asymptotic normality of the corresponding test statistic is established. As the asymptotic variance depends on the location and order of the critical roots of the equation $| \mu_1 - \mu_t | = c$ a new bootstrap procedure is proposed to obtain critical values and its consistency is established. As a consequence we are able to quantitatively describe relevant deviations of a non-stationary sequence from its initial value. The results are illustrated by means of a simulation study and by analyzing data examples.