title: Subgradients of Minimal Time Functions without Calmness

abstract: In recent years there has been great interest in variational analysis of a class of nonsmooth functions called the minimal time function. In this paper we continue this line of research by providing new results on generalized differentiation of this class of functions, relaxing assumptions imposed on the functions and sets involved for the results. In particular, we focus on the singular subdifferential and the limiting subdifferential of this class of functions.


title: Boosting the accuracy of multi-spectral image pan-sharpening by learning a deep residual network

abstract: In the field of fusing multi-spectral and panchromatic images (Pan-sharpening), the impressive effectiveness of deep neural networks has been recently employed to overcome the drawbacks of traditional linear models and boost the fusing accuracy. However, to the best of our knowledge, existing research works are mainly based on simple and flat networks with relatively shallow architecture, which severely limited their performances. In this paper, the concept of residual learning has been introduced to form a very deep convolutional neural network to make a full use of the high non-linearity of deep learning models. By both quantitative and visual assessments on a large number of high quality multi-spectral images from various sources, it has been supported that our proposed model is superior to all mainstream algorithms included in the comparison, and achieved the highest spatial-spectral unified accuracy.


title: Whittaker models for depth zero representations of covering groups

abstract: We study the dimension of the space of Whittaker functionals for depth zero representations of covering groups. In particular, we determine such dimensions for arbitrary Brylinski-Deligne coverings of the general linear group. The results in the paper are motivated by and compatible with the work of Howard and the second author, and earlier work by Blondel.


title: Note on Evolution and Forecasting of Requirements: Communications Example

abstract: Combinatorial evolution and forecasting of system requirements is examined. The morphological model is used for a hierarchical requirements system (i.e., system parts, design alternatives for the system parts, ordinal estimates for the alternatives). A set of system changes involves changes of the system structure, component alternatives and their estimates. The composition process of the forecast is based on combinatorial synthesis (knapsack problem, multiple choice problem, hierarchical morphological design). An illustrative numerical example for four-phase evolution and forecasting of requirements to communications is described.


title: L-groups and the Langlands program for covering groups: a historical introduction

abstract: In this joint introduction to an Asterisque volume, we give a short discussion of the historical developments in the study of nonlinear covering groups, touching on their structure theory, representation theory and the theory of automorphic forms. This serves as a historical motivation and sets the scene for the papers in the volume. Our discussion is necessarily subjective and will undoubtedly leave out the contributions of many authors, to whom we apologize in earnest.


title: Exoplanet Biosignatures: Understanding Oxygen as a Biosignature in the Context of Its Environment

abstract: Here we review how environmental context can be used to interpret whether O2 is a biosignature in extrasolar planetary observations. This paper builds on the overview of current biosignature research discussed in Schwieterman et al. (2017), and provides an in-depth, interdisciplinary example of biosignature identification and observation that serves as a basis for the development of the general framework for biosignature assessment described in Catling et al., (2017). O2 is a potentially strong biosignature that was originally thought to be an unambiguous indicator for life at high-abundance. We describe the coevolution of life with the early Earth's environment, and how the interplay of sources and sinks in the planetary environment may have resulted in suppression of O2 release into the atmosphere for several billion years, a false negative for biologically generated O2. False positives may also be possible, with recent research showing potential mechanisms in exoplanet environments that may generate relatively high abundances of atmospheric O2 without a biosphere being present. These studies suggest that planetary characteristics that may enhance false negatives should be considered when selecting targets for biosignature searches. Similarly our ability to interpret O2 observed in an exoplanetary atmosphere is also crucially dependent on environmental context to rule out false positive mechanisms. We describe future photometric, spectroscopic and time-dependent observations of O2 and the planetary environment that could increase our confidence that any observed O2 is a biosignature, and help discriminate it from potential false positives. By observing and understanding O2 in its planetary context we can increase our confidence in the remote detection of life, and provide a model for biosignature development for other proposed biosignatures.


title: Detection Estimation and Grid matching of Multiple Targets with Single Snapshot Measurements

abstract: In this work, we explore the problems of detecting the number of narrow-band, far-field targets and estimating their corresponding directions from single snapshot measurements. The principles of sparse signal recovery (SSR) are used for the single snapshot detection and estimation of multiple targets. In the SSR framework, the DoA estimation problem is grid based and can be posed as the lasso optimization problem. However, the SSR framework for DoA estimation gives rise to the grid mismatch problem, when the unknown targets (sources) are not matched with the estimation grid chosen for the construction of the array steering matrix at the receiver. The block sparse recovery framework is known to mitigate the grid mismatch problem by jointly estimating the targets and their corresponding offsets from the estimation grid using the group lasso estimator. The corresponding detection problem reduces to estimating the optimal regularization parameter ($\tau$) of the lasso (in case of perfect grid-matching) or group-lasso estimation problem for achieving the required probability of correct detection ($P_c$). We propose asymptotic and finite sample test statistics for detecting the number of sources with the required $P_c$ at moderate to high signal to noise ratios. Once the number of sources are detected, or equivalently the optimal $\hat{\tau}$ is estimated, the corresponding estimation and grid matching of the DoAs can be performed by solving the lasso or group-lasso problem at $\hat{\tau}$


title: On the diffusion approximation of nonconvex stochastic gradient descent

abstract: We study the Stochastic Gradient Descent (SGD) method in nonconvex optimization problems from the point of view of approximating diffusion processes. We prove rigorously that the diffusion process can approximate the SGD algorithm weakly using the weak form of master equation for probability evolution. In the small step size regime and the presence of omnidirectional noise, our weak approximating diffusion process suggests the following dynamics for the SGD iteration starting from a local minimizer (resp.~saddle point): it escapes in a number of iterations exponentially (resp.~almost linearly) dependent on the inverse stepsize. The results are obtained using the theory for random perturbations of dynamical systems (theory of large deviations for local minimizers and theory of exiting for unstable stationary points). In addition, we discuss the effects of batch size for the deep neural networks, and we find that small batch size is helpful for SGD algorithms to escape unstable stationary points and sharp minimizers. Our theory indicates that one should increase the batch size at later stage for the SGD to be trapped in flat minimizers for better generalization.


title: Learning to Rank Using Localized Geometric Mean Metrics

abstract: Many learning-to-rank (LtR) algorithms focus on query-independent model, in which query and document do not lie in the same feature space, and the rankers rely on the feature ensemble about query-document pair instead of the similarity between query instance and documents. However, existing algorithms do not consider local structures in query-document feature space, and are fragile to irrelevant noise features. In this paper, we propose a novel Riemannian metric learning algorithm to capture the local structures and develop a robust LtR algorithm. First, we design a concept called \textit{ideal candidate document} to introduce metric learning algorithm to query-independent model. Previous metric learning algorithms aiming to find an optimal metric space are only suitable for query-dependent model, in which query instance and documents belong to the same feature space and the similarity is directly computed from the metric space. Then we extend the new and extremely fast global Geometric Mean Metric Learning (GMML) algorithm to develop a localized GMML, namely L-GMML. Based on the combination of local learned metrics, we employ the popular Normalized Discounted Cumulative Gain~(NDCG) scorer and Weighted Approximate Rank Pairwise (WARP) loss to optimize the \textit{ideal candidate document} for each query candidate set. Finally, we can quickly evaluate all candidates via the similarity between the \textit{ideal candidate document} and other candidates. By leveraging the ability of metric learning algorithms to describe the complex structural information, our approach gives us a principled and efficient way to perform LtR tasks. The experiments on real-world datasets demonstrate that our proposed L-GMML algorithm outperforms the state-of-the-art metric learning to rank methods and the stylish query-independent LtR algorithms regarding accuracy and computational efficiency.


title: Difference equations and pseudo-differential operators on $\mathbb{Z}^n$

abstract: In this paper we develop the calculus of pseudo-differential operators on the lattice $\mathbb{Z}^n$, which we can call pseudo-difference operators. An interesting feature of this calculus is that the phase space is compact so the symbol classes are defined in terms of the behaviour with respect to the lattice variable. We establish formulae for composition, adjoint, transpose, and for parametrix for the elliptic operators. We also give conditions for the $\ell^2$, weighted $\ell^2$, and $\ell^p$ boundedness of operators and for their compactness on $\ell^p$. We describe a link to the toroidal quantization on the torus $\mathbb{T}^n$, and apply it to give conditions for the membership in Schatten classes on $\ell^2(\mathbb{Z}^n)$. Furthermore, we discuss a version of Fourier integral operators on the lattice and give conditions for their $\ell^2$-boundedness. The results are applied to give estimates for solutions to difference equations on the lattice $\mathbb{Z}^n$. Moreover, we establish Garding and sharp Garding inequalities, with an application to the unique solvability of parabolic equations on the lattice $\mathbb{Z}^n$.


title: Learning to Prune Deep Neural Networks via Layer-wise Optimal Brain Surgeon

abstract: How to develop slim and accurate deep neural networks has become crucial for real- world applications, especially for those employed in embedded systems. Though previous work along this research line has shown some promising results, most existing methods either fail to significantly compress a well-trained deep network or require a heavy retraining process for the pruned deep network to re-boost its prediction performance. In this paper, we propose a new layer-wise pruning method for deep neural networks. In our proposed method, parameters of each individual layer are pruned independently based on second order derivatives of a layer-wise error function with respect to the corresponding parameters. We prove that the final prediction performance drop after pruning is bounded by a linear combination of the reconstructed errors caused at each layer. Therefore, there is a guarantee that one only needs to perform a light retraining process on the pruned network to resume its original prediction performance. We conduct extensive experiments on benchmark datasets to demonstrate the effectiveness of our pruning method compared with several state-of-the-art baseline methods.


title: Hypergroups derived from random walks on some infinite graphs

abstract: Wildberger gave a method to construct a finite hermitian discrete hypergroup from a random walk on a certain kind of graphs. In this article, we reveal that his method is applicable to a random walk on a certain kind of infinite graphs. Moreover, we make some observations of finite or infinite graphs on which a random walk produces a hermitian discrete hypergroup.


title: Minimal coloring number for $\mathbb{Z}$-colorable links II

abstract: The minimal coloring number of a $\mathbb{Z}$-colorable link is the minimal number of colors for non-trivial $\mathbb{Z}$-colorings on diagrams of the link. In this paper, we show that the minimal coloring number of any non-splittable $\mathbb{Z}$-colorable links is four. As an example, we consider the link obtained by replacing each component of the given link with several parallel strands, which we call a parallel of a link. We show that an even parallel of a link is $\mathbb{Z}$-colorable except for the case of 2 parallels with non-zero linking number. We then give a simple way to obtain a diagram which attains the minimal coloring number for such even parallels of links.


title: Understanding the nature of the intriguing source X Persei: A deep look with a Suzaku observation

abstract: We present detailed broadband timing and spectral analysis of the persistent, low luminosity and slowly spinning pulsar X Persei using a deep Suzaku observation of the source. The spectrum is unusually hard with a cyclotron resonance scattering feature (CRSF), the presence of which has been debated. By comparing the spectral models relevant for accretion powered pulsars, we have obtained the best constraint on the broadband spectral model of X Persei obtained so far. The CRSF is not confirmed in the average spectrum. We have also identified, the presence of different intensity levels in the source with distinct changes in the pulse profile and energy spectrum indicating changes in the accretion geometry. We further find evidence of a CRSF in the highest intensity levels at $\sim$ 40 keV, indicating a magnetic field strength of 3.4$\times10^{12}$ G.


title: One-dimensional phosphorus chain and two-dimensional blue phosphorene grown on Au(111) by molecular-beam epitaxy

abstract: Single layer (SL) phosphorus (phosphorene) has drawn considerable research attention recently as a two-dimensional (2D) material for application promises. It is a semiconductor showing superior transport and optical properties. Few-layer or SL black phosphorus has been successfully isolated by exfoliation from bulk crystals and extensively studied thereof for its electronic and optical properties. Blue phosphorus (blueP), an allotrope of black phosphorus where atoms are arranged in a more flat atomic configuration, has been recently suggested by theory to exist in the SL form on some substrates. In this work, we report the formation of a blueP-like epilayer on Au(111) by molecular-beam epitaxy. In particular, we uncover by scanning tunneling microscopy (STM) one-dimensional (1D) atomic chains at low coverage, which develop into more compact islands or patches of $(\sqrt{3}\times\sqrt{3})R30^\circ$ structure with increasing coverage before blueP-like islands nucleate and grow. We also note an interesting growth characteristic where the $(\sqrt{3}\times\sqrt{3})R30^\circ$ surface at intermediate coverage tends to phase-separate into locally low-coverage 1D chain and high-coverage blueP-like structures, respectively. This experiment thus not only lends a support of the recently proposed half-layer by half-layer (HLBHL) growth mechanism but also reveals the kinetic details of blueP growth processes.


title: Many-Body-Localization : Strong Disorder perturbative approach for the Local Integrals of Motion

abstract: For random quantum spin models, the strong disorder perturbative expansion of the Local Integrals of Motion (LIOMs) around the real-spin operators is revisited. The emphasis is on the links with other properties of the Many-Body-Localized phase, in particular the memory in the dynamics of the local magnetizations and the statistics of matrix elements of local operators in the eigenstate basis. Finally, this approach is applied to analyze the Many-Body-Localization transition in a toy model studied previously from the point of view of the entanglement entropy.


title: Analyzing of singlet fermionic dark matter via the updated direct detection data

abstract: We revisit the parameter space of singlet fermionic cold dark matter model in order to determine the role of the mixing angle between the standard model Higgs and new singlet one. Furthermore, we restudy the direct detection constraints with the updated and new experimental data. As an important conclusion, this model is completely excluded by recent XENON100, PandaX II and LUX data.


title: Cooling of trapped ions by resonant charge exchange

abstract: The two most widely used ion cooling methods are laser cooling and sympathetic cooling by elastic collisions (ECs). Here we demonstrate another method of cooling ions that is based on resonant charge exchange (RCE) between the trapped ion and the ultracold parent atom. Specifically, trapped Cs$^{+}$ ions are cooled by collisions with co-trapped, ultracold Cs atoms and, separately, by collisions with co-trapped, ultracold Rb atoms. We observe that the cooling of Cs$^{+}$ ions by Cs atoms is more efficient than cooling of Cs$^{+}$ ions by Rb atoms. This signals the presence of a cooling mechanism apart from the elastic ion-atom collision channel for the Cs-Cs$^{+}$ case, which is cooling by RCE. The efficiency of cooling by RCE is experimentally determined and the per-collision cooling is found to be two orders of magnitude higher than cooling by EC. The result provides the experimental basis for future studies on charge transport by electron hopping in atom-ion hybrid systems.

categories: physics.atom-ph cond-mat.quant-gas physics.chem-ph quant-ph


title: Ruelle-Pollicott Resonances of Stochastic Systems in Reduced State Space. Part II: Stochastic Hopf Bifurcation

abstract: The spectrum of the generator (Kolmogorov operator) of a diffusion process, referred to as the Ruelle-Pollicott (RP) spectrum, provides a detailed characterization of correlation functions and power spectra of stochastic systems via decomposition formulas in terms of RP resonances. Stochastic analysis techniques relying on the theory of Markov semigroups for the study of the RP spectrum and a rigorous reduction method is presented in Part I. This framework is here applied to study a stochastic Hopf bifurcation in view of characterizing the statistical properties of nonlinear oscillators perturbed by noise, depending on their stability. In light of the H\"ormander theorem, it is first shown that the geometry of the unperturbed limit cycle, in particular its isochrons, is essential to understand the effect of noise and the phenomenon of phase diffusion. In addition, it is shown that the spectrum has a spectral gap, even at the bifurcation point, and that correlations decay exponentially fast. Explicit small-noise expansions of the RP eigenvalues and eigenfunctions are then obtained, away from the bifurcation point, based on the knowledge of the linearized deterministic dynamics and the characteristics of the noise. These formulas allow one to understand how the interaction of the noise with the deterministic dynamics affect the decay of correlations. Numerical results complement the study of the RP spectrum at the bifurcation, revealing useful scaling laws. The analysis of the Markov semigroup for stochastic bifurcations is thus promising in providing a complementary approach to the more geometric random dynamical system approach. This approach is not limited to low-dimensional systems and the reduction method presented in part I is applied to a stochastic model relevant to climate dynamics in part III.


title: Self-Fulfilling Signal of an Endogenous State in Network Congestion Games

abstract: We consider the problem of coordination via signaling in network congestion games to improve social welfare deteriorated by incomplete information about traffic flow. Traditional studies on signaling, which focus on exogenous factors of congestion and ignore congestion externalities, fail to discuss the oscillations of traffic flow. To address this gap, we formulate a problem of designing a coordination signal on endogenous information about traffic flow and introduce a it self-fulfilling characteristic of a signal that guarantees an outcome flow consistent with the signal itself without causing the unwanted oscillation. An instance of the self-fulfilling signal is shown in the case of a Gaussian signal distribution. In addition, we show simple numerical examples. The results reveal how a self-fulfilling signal suppresses the oscillation and simultaneously improves social welfare through improved network efficiency.


title: Mira: A Framework for Static Performance Analysis

abstract: The performance model of an application can pro- vide understanding about its runtime behavior on particular hardware. Such information can be analyzed by developers for performance tuning. However, model building and analyzing is frequently ignored during software development until perfor- mance problems arise because they require significant expertise and can involve many time-consuming application runs. In this paper, we propose a fast, accurate, flexible and user-friendly tool, Mira, for generating performance models by applying static program analysis, targeting scientific applications running on supercomputers. We parse both the source code and binary to estimate performance attributes with better accuracy than considering just source or just binary code. Because our analysis is static, the target program does not need to be executed on the target architecture, which enables users to perform analysis on available machines instead of conducting expensive exper- iments on potentially expensive resources. Moreover, statically generated models enable performance prediction on non-existent or unavailable architectures. In addition to flexibility, because model generation time is significantly reduced compared to dynamic analysis approaches, our method is suitable for rapid application performance analysis and improvement. We present several scientific application validation results to demonstrate the current capabilities of our approach on small benchmarks and a mini application.


title: Global Guarantees for Enforcing Deep Generative Priors by Empirical Risk

abstract: We examine the theoretical properties of enforcing priors provided by generative deep neural networks via empirical risk minimization. In particular we consider two models, one in which the task is to invert a generative neural network given access to its last layer and another in which the task is to invert a generative neural network given only compressive linear observations of its last layer. We establish that in both cases, in suitable regimes of network layer sizes and a randomness assumption on the network weights, that the non-convex objective function given by empirical risk minimization does not have any spurious stationary points. That is, we establish that with high probability, at any point away from small neighborhoods around two scalar multiples of the desired solution, there is a descent direction. Hence, there are no local minima, saddle points, or other stationary points outside these neighborhoods. These results constitute the first theoretical guarantees which establish the favorable global geometry of these non-convex optimization problems, and they bridge the gap between the empirical success of enforcing deep generative priors and a rigorous understanding of non-linear inverse problems.


title: Semiparametric Efficient Empirical Higher Order Influence Function Estimators

abstract: Robins et al. (2008, 2017) applied the theory of higher order influence functions (HOIFs) to derive an estimator of the mean $\psi$ of an outcome Y in a missing data model with Y missing at random conditional on a vector X of continuous covariates; their estimator, in contrast to previous estimators, is semiparametric efficient under the minimal conditions of Robins et al. (2009b), together with an additional (non-minimal) smoothness condition on the density g of X, because the Robins et al. (2008, 2017) estimator depends on a nonparametric estimate of g. In this paper, we introduce a new HOIF estimator that has the same asymptotic properties as the original one, but does not impose any smoothness requirement on g. This is important for two reasons. First, one rarely has the knowledge about the properties of g. Second, even when g is smooth, if the dimension of X is even moderate, accurate nonparametric estimation of its density is not feasible at the sample sizes often encountered in applications. In fact, to the best of our knowledge, this new HOIF estimator remains the only semiparametric efficient estimator of $\psi$ under minimal conditions, despite the rapidly growing literature on causal effect estimation. We also show that our estimator can be generalized to the entire class of functionals considered by Robins et al. (2008) which include the average effect of a treatment on a response Y when a vector X suffices to control confounding and the expected conditional variance of a response Y given a vector X. Simulation experiments are also conducted, which demonstrate that our new estimator outperforms those of Robins et al. (2008, 2017) in finite samples, when g is not very smooth.


title: Semiparametric estimation in the normal variance-mean mixture model

abstract: In this paper we study the problem of statistical inference on the parameters of the semiparametric variance-mean mixtures. This class of mixtures has recently become rather popular in statistical and financial modelling. We design a semiparametric estimation procedure that first estimates the mean of the underlying normal distribution and then recovers nonparametrically the density of the corresponding mixing distribution. We illustrate the performance of our procedure on simulated and real data.


title: Lyapunov optimization for non-generic one-dimensional expanding Markov maps

abstract: For a non-generic, yet dense subset of $C^1$ expanding Markov maps of the interval we prove the existence of uncountably many Lyapunov optimizing measures which are ergodic, fully supported and have positive entropy. These measures are equilibrium states for some H\"older continuous potentials. We also prove the existence of another non-generic dense subset for which the optimizing measure is unique and supported on a periodic orbit. A key ingredient is a new $C^1$ perturbation lemma which allows us to interpolate between expanding Markov maps and the shift map on a finite number of symbols.


title: A note on the Morse index of 2k-ended phase transitions in R^2

abstract: We show that the Morse index of every 2k-ended solution of the Allen-Cahn equation in R^2 is >= k-1. This bound is expected to be sharp.


title: The detection of an extremely bright fast radio burst in a phased array feed survey

abstract: We report the detection of an ultra-bright fast radio burst (FRB) from a modest, 3.4-day pilot survey with the Australian Square Kilometre Array Pathfinder. The survey was conducted in a wide-field fly's-eye configuration using the phased-array-feed technology deployed on the array to instantaneously observe an effective area of $160$ deg$^2$, and achieve an exposure totaling $13200$ deg$^2$ hr. We constrain the position of FRB 170107 to a region $8'\times8'$ in size (90% containment) and its fluence to be $58\pm6$ Jy ms. The spectrum of the burst shows a sharp cutoff above $1400$ MHz, which could be either due to scintillation or an intrinsic feature of the burst. This confirms the existence of an ultra-bright ($>20$ Jy ms) population of FRBs.


title: The Spatial Range of Conformity

abstract: Properties of galaxies like their absolute magnitude and their stellar mass content are correlated. These correlations are tighter for close pairs of galaxies, which is called galactic conformity. In hierarchical structure formation scenarios, galaxies form within dark matter halos. To explain the amplitude and the spatial range of galactic conformity two--halo terms or assembly bias become important. With the scale dependent correlation coefficients the amplitude and the spatial range of conformity are determined from galaxy and halo samples. The scale dependent correlation coefficients are introduced as a new descriptive statistic to quantify the correlations between properties of galaxies or halos, depending on the distances to other galaxies or halos. These scale dependent correlation coefficients can be applied to the galaxy distribution directly. Neither a splitting of the sample into subsamples, nor an a priori clustering is needed. This new descriptive statistic is applied to galaxy catalogues derived from the Sloan Digital Sky Survey III and to halo catalogues from the MultiDark simulations. In the galaxy sample the correlations between absolute Magnitude, velocity dispersion, ellipticity, and stellar mass content are investigated. The correlations of mass, spin, and ellipticity are explored in the halo samples. Both for galaxies and halos a scale dependent conformity is confirmed. Moreover the scale dependent correlation coefficients reveal a signal of conformity out to 40Mpc and beyond. The halo and galaxy samples show a differing amplitude and range of conformity.


title: A discrete Schrodinger equation via optimal transport on graphs

abstract: In 1966, Edward Nelson presented an interesting derivation of the Schrodinger equation using Brownian motion. Recently, this derivation is linked to the theory of optimal transport, which shows that the Schrodinger equation is a Hamiltonian system on the probability density manifold equipped with the Wasserstein metric. In this paper, we consider similar matters on a finite graph. By using discrete optimal transport and its corresponding Nelson's approach, we derive a discrete Schrodinger equation on a finite graph. The proposed system is quite different from the commonly referred discretized Schrodinger equations. It is a system of nonlinear ordinary differential equations (ODEs) with many desirable properties. Several numerical examples are presented to illustrate the properties.


title: Counting the solutions of $\lambda_1 x_1^{k_1}+\cdots +\lambda_t x_t^{k_t}\equiv c\bmod{n}$

abstract: Given a polynomial $Q(x_1,\cdots, x_t)=\lambda_1 x_1^{k_1}+\cdots +\lambda_t x_t^{k_t}$, for every $c\in \mathbb{Z}$ and $n\geq 2$, we study the number of solutions $N_J(Q;c,n)$ of the congruence equation $Q(x_1,\cdots, x_t)\equiv c\bmod{n}$ in $(\mathbb{Z}/n\mathbb{Z})^t$ such that $x_i\in (\mathbb{Z}/n\mathbb{Z})^\times$ for $i\in J\subseteq I= \{1,\cdots, t\}$. We deduce formulas and an algorithm to study $N_J(Q; c,p^a)$ for $p$ any prime number and $a\geq 1$ any integer. As consequences of our main results, we completely solve: the counting problem of $Q(x_i)=\sum\limits_{i\in I}\lambda_i x_i$ for any prime $p$ and any subset $J$ of $I$; the counting problem of $Q(x_i)=\sum\limits_{i\in I}\lambda_i x^2_i$ in the case $t=2$ for any $p$ and $J$, and the case $t$ general for any $p$ and $J$ satisfying $\min\{v_p(\lambda_i)\mid i\in I\}=\min\{v_p(\lambda_i)\mid i\in J\}$; the counting problem of $Q(x_i)=\sum\limits_{i\in I}\lambda_i x^k_i$ in the case $t=2$ for any $p\nmid k$ and any $J$, and in the case $t$ general for any $p\nmid k$ and $J$ satisfying $\min\{v_p(\lambda_i)\mid i\in I\}=\min\{v_p(\lambda_i)\mid i\in J\}$.


title: Union of Intersections (UoI) for Interpretable Data Driven Discovery and Prediction

abstract: The increasing size and complexity of scientific data could dramatically enhance discovery and prediction for basic scientific applications. Realizing this potential, however, requires novel statistical analysis methods that are both interpretable and predictive. We introduce Union of Intersections (UoI), a flexible, modular, and scalable framework for enhanced model selection and estimation. Methods based on UoI perform model selection and model estimation through intersection and union operations, respectively. We show that UoI-based methods achieve low-variance and nearly unbiased estimation of a small number of interpretable features, while maintaining high-quality prediction accuracy. We perform extensive numerical investigation to evaluate a UoI algorithm ($UoI_{Lasso}$) on synthetic and real data. In doing so, we demonstrate the extraction of interpretable functional networks from human electrophysiology recordings as well as accurate prediction of phenotypes from genotype-phenotype data with reduced features. We also show (with the $UoI_{L1Logistic}$ and $UoI_{CUR}$ variants of the basic framework) improved prediction parsimony for classification and matrix factorization on several benchmark biomedical data sets. These results suggest that methods based on the UoI framework could improve interpretation and prediction in data-driven discovery across scientific fields.


title: Achiral tilted domain walls in perpendicularly magnetized nanowires

abstract: Perpendicularly magnetized nanowires exhibit distinct domain wall types depending on the geometry. Wide wires contain Bloch walls, narrow wires N\'eel walls. Here, the transition region is investigated by direct imaging of the wall structure using high-resolution spin-polarized scanning electron microscopy. An achiral intermediate wall type is discovered that is unpredicted by established theoretical models. With the help of micromagnetic simulations, the formation of this novel wall type is explained.


title: Universal Behaviors of Speed of Sound from Holography

abstract: We investigate the speed of sound in $(d+1)$-dimensional field theory by studying its dual $(d+2)$-dimensional gravity theory from gauge/gravity correspondence. Instead of the well known conformal limit $c_s^2 \rightarrow 1/d$ at high temperature, we reveal two more universal quantities in various limits: $c_s^2 \rightarrow (d-1)/16\pi$ at low temperature and $c_s^2 \rightarrow (d-1)/16\pi d$ at large chemical potential.


title: Predictions of homogeneous nucleation rates for $n$-alkanes accounting for the diffuse phase interface and capillary waves

abstract: Homogeneous droplet nucleation has been studied for almost a century but has not yet been fully understood. In this work, we used the density gradient theory (DGT) and considered the influence of capillary waves (CW) on the predicted size-dependent surface tensions and nucleation rates for selected $n$-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) and compared to the classical nucleation theory and the Peng--Robinson EoS. It was found that the critical clusters are practically free of CW because they are so small that even the smallest CW wavelengths do not fit into their finite dimensions. The CW contribute to the entropy of the system and thus decrease the surface tension. A correction for the effect of CW on the surface tension is presented. The effect of the different EoSs is relatively small because by a fortuitous coincidence their predictions are similar in the relevant range of critical cluster sizes. The difference of the DGT predictions to the classical nucleation theory computations is important but not decisive. Of the effects investigated, the most pronounced is the suppression of the CW which causes a sizable decrease of the predicted nucleation rates. The major difference between experimental nucleation rate data and theoretical predictions remains in the temperature dependence. For normal alkanes, this discrepancy is much stronger than observed, e.g., for water. The theoretical corrections developed here have a minor influence on the temperature dependency. We provide empirical equations correcting the predicted nucleation rates to values comparable with experiments.


title: Eye Tracker Accuracy: Quantitative Evaluation of the Invisible Eye Center Location

abstract: Purpose. We present a new method to evaluate the accuracy of an eye tracker based eye localization system. Measuring the accuracy of an eye tracker's primary intention, the estimated point of gaze, is usually done with volunteers and a set of fixation points used as ground truth. However, verifying the accuracy of the location estimate of a volunteer's eye center in 3D space is not easily possible. This is because the eye center is an intangible point hidden by the iris. Methods. We evaluate the eye location accuracy by using an eye phantom instead of eyes of volunteers. For this, we developed a testing stage with a realistic artificial eye and a corresponding kinematic model, which we trained with {\mu}CT data. This enables us to precisely evaluate the eye location estimate of an eye tracker. Results. We show that the proposed testing stage with the corresponding kinematic model is suitable for such a validation. Further, we evaluate a particular eye tracker based navigation system and show that this system is able to successfully determine the eye center with sub-millimeter accuracy. Conclusions. We show the suitability of the evaluated eye tracker for eye interventions, using the proposed testing stage and the corresponding kinematic model. The results further enable specific enhancement of the navigation system to potentially get even better results.


title: Spin Currents of Charged Dirac Particles in Rotating Coordinates

abstract: The semiclassical Boltzmann transport equation of charged, massive fermions in a rotating frame of reference, in the presence of external electromagnetic fields is solved in the relaxation time approach to establish the distribution function up to linear order in the electric field in rotating coordinates, centrifugal force and the derivatives. The spin and spin current densities are calculated by means of this distribution function at zero temperature up to the first order. It is shown that the nonequilibrium part of the distribution function yields the spin Hall effect for fermions constrained to move in a plane perpendicular to the angular velocity and magnetic field. Moreover it yields an analogue of Ohm's law for spin currents whose resistivity depends on the external magnetic field and the angular velocity of the rotating frame. Spin current densities in three-dimensional systems are also established.

categories: math-ph cond-mat.mes-hall hep-th math.MP


title: Leptoquark models for the $B$-physics anomalies

abstract: The $B$-physics experiments at LHCb, BaBar and Belle hint towards deviations from Lepton Flavor Universality in both the tree-level and loop-induced $B$ meson semileptonic decays. We propose a leptoquark model with light right-handed neutrinos which can accommodate both $R_K^\mathrm{exp}<R_K^\mathrm{SM}$ and $R_{D^{(\ast)}}^\mathrm{exp}>R_{D^{(\ast)}}^\mathrm{SM}$. We discuss several of its predictions which can be tested in modern day experiments. We also comment on the recent finding at LHCb, namely $R_{K^\ast}^{\mathrm{exp}}<R_{K^\ast}^\mathrm{SM}$.


title: Improved Clustering with Augmented k-means

abstract: Identifying a set of homogeneous clusters in a heterogeneous dataset is one of the most important classes of problems in statistical modeling. In the realm of unsupervised partitional clustering, k-means is a very important algorithm for this. In this technical report, we develop a new k-means variant called Augmented k-means, which is a hybrid of k-means and logistic regression. During each iteration, logistic regression is used to predict the current cluster labels, and the cluster belonging probabilities are used to control the subsequent re-estimation of cluster means. Observations which can't be firmly identified into clusters are excluded from the re-estimation step. This can be valuable when the data exhibit many characteristics of real datasets such as heterogeneity, non-sphericity, substantial overlap, and high scatter. Augmented k-means frequently outperforms k-means by more accurately classifying observations into known clusters and / or converging in fewer iterations. We demonstrate this on both simulated and real datasets. Our algorithm is implemented in Python and will be available with this report.


title: The structure of random automorphisms of countable structures

abstract: In order to understand the structure of the `typical' element of an automorphism group, one has to study how large the conjugacy classes of the group are. When typical is meant in the sense of Baire category, a complete description of the size of the conjugacy classes has been given by Kechris and Rosendal. Following Dougherty and Mycielski we investigate the measure theoretic dual of this problem, using Christensen's notion of Haar null sets. When typical means random, that is, almost every with respect to this notion of Haar null sets, the behavior of the automorphisms is entirely different from the Baire category case. In this paper, we generalize the theorems of Dougherty and Mycielski about $S_\infty$ to arbitrary automorphism groups of countable structures isolating a new model theoretic property, the Cofinal Strong Amalgamation Property. As an application we show that a large class of automorphism groups can be decomposed into the union of a meager and a Haar null set.


title: Learning Robust Object Recognition Using Composed Scenes from Generative Models

abstract: Recurrent feedback connections in the mammalian visual system have been hypothesized to play a role in synthesizing input in the theoretical framework of analysis by synthesis. The comparison of internally synthesized representation with that of the input provides a validation mechanism during perceptual inference and learning. Inspired by these ideas, we proposed that the synthesis machinery can compose new, unobserved images by imagination to train the network itself so as to increase the robustness of the system in novel scenarios. As a proof of concept, we investigated whether images composed by imagination could help an object recognition system to deal with occlusion, which is challenging for the current state-of-the-art deep convolutional neural networks. We fine-tuned a network on images containing objects in various occlusion scenarios, that are imagined or self-generated through a deep generator network. Trained on imagined occluded scenarios under the object persistence constraint, our network discovered more subtle and localized image features that were neglected by the original network for object classification, obtaining better separability of different object classes in the feature space. This leads to significant improvement of object recognition under occlusion for our network relative to the original network trained only on un-occluded images. In addition to providing practical benefits in object recognition under occlusion, this work demonstrates the use of self-generated composition of visual scenes through the synthesis loop, combined with the object persistence constraint, can provide opportunities for neural networks to discover new relevant patterns in the data, and become more flexible in dealing with novel situations.


title: Magnetic field-induced enhancement of the nitrogen-vacancy fluorescence quantum yield

abstract: The nitrogen-vacancy (NV) centre in diamond is a unique optical defect that is used in many applications today and methods to enhance its fluorescence brightness are highly sought after. We observed experimentally an enhancement of the NV quantum yield by up to 7% in bulk diamond caused by an external magnetic field relative to the field-free case. This observation is rationalised phenomenologically in terms of a magnetic field dependence of the NV excited state triplet-to-singlet transition rate. The theoretical model is in good qualitative agreement with the experimental results at low excitation intensities. Our results significantly contribute to our fundamental understanding of the photophysical properties of the NV defect in diamond and may enable novel NV centre-based magnetometry techniques.


title: Galaxy And Mass Assembly: the evolution of the cosmic spectral energy distribution from z = 1 to z = 0

abstract: We present the evolution of the Cosmic Spectral Energy Distribution (CSED) from $z = 1 - 0$. Our CSEDs originate from stacking individual spectral energy distribution fits based on panchromatic photometry from the Galaxy and Mass Assembly (GAMA) and COSMOS datasets in ten redshift intervals with completeness corrections applied. Below $z = 0.45$, we have credible SED fits from 100 nm to 1 mm. Due to the relatively low sensitivity of the far-infrared data, our far-infrared CSEDs contain a mix of predicted and measured fluxes above $z = 0.45$. Our results include appropriate errors to highlight the impact of these corrections. We show that the bolometric energy output of the Universe has declined by a factor of roughly four -- from $5.1 \pm 1.0$ at $z \sim 1$ to $1.3 \pm 0.3 \times 10^{35}~h_{70}$~W~Mpc$^{-3}$ at the current epoch. We show that this decrease is robust to cosmic variance, SED modelling and other various types of error. Our CSEDs are also consistent with an increase in the mean age of stellar populations. We also show that dust attenuation has decreased over the same period, with the photon escape fraction at 150~nm increasing from $16 \pm 3$ at $z \sim 1$ to $24 \pm 5$ per cent at the current epoch, equivalent to a decrease in $A_\mathrm{FUV}$ of 0.4~mag. Our CSEDs account for $68 \pm 12$ and $61 \pm 13$ per cent of the cosmic optical and infrared backgrounds respectively as defined from integrated galaxy counts and are consistent with previous estimates of the cosmic infrared background with redshift.


title: Finite-size scaling of out-of-time-ordered correlators at late times

abstract: Chaotic dynamics in quantum many-body systems scrambles local information so that at late times it can no longer be accessed locally. This is reflected quantitatively in the out-of-time-ordered correlator of local operators, which is expected to decay to zero with time. However, for systems of finite size, out-of-time-ordered correlators do not decay exactly to zero and in this paper we show that the residual value can provide useful insights into the chaotic dynamics. When energy is conserved, the late-time saturation value of the out-of-time-ordered correlator of generic traceless local operators scales as an inverse polynomial in the system size. This is in contrast to the inverse exponential scaling expected for chaotic dynamics without energy conservation. We provide both analytical arguments and numerical simulations to support this conclusion.

categories: quant-ph cond-mat.stat-mech cond-mat.str-el hep-th


title: Rao-Blackwellized Particle Smoothing as Message Passing

abstract: In this manuscript the fixed-lag smoothing problem for conditionally linear Gaussian state-space models is investigated from a factor graph perspective. More specifically, after formulating Bayesian smoothing for an arbitrary state-space model as forward-backward message passing over a factor graph, we focus on the above mentioned class of models and derive a novel Rao-Blackwellized particle smoother for it. Then, we show how our technique can be modified to estimate a point mass approximation of the so called joint smoothing distribution. Finally, the estimation accuracy and the computational requirements of our smoothing algorithms are analysed for a specific state-space model.


title: Toric Varieties vs. Horofunction Compactifications of Polyhedral Norms

abstract: We establish a natural and geometric 1-1 correspondence between projective toric varieties of dimension $n$ and horofunction compactifications of $\mathbb{R}^n$ with respect to rational polyhedral norms. For this purpose, we explain a topological model of toric varieties. Consequently, toric varieties in algebraic geometry, normed spaces in convex analysis, and horofunction compactifications in metric geometry are directly and explicitly related.


title: Classification Using Proximity Catch Digraphs (Technical Report)

abstract: We employ random geometric digraphs to construct semi-parametric classifiers. These data-random digraphs are from parametrized random digraph families called proximity catch digraphs (PCDs). A related geometric digraph family, class cover catch digraph (CCCD), has been used to solve the class cover problem by using its approximate minimum dominating set. CCCDs showed relatively good performance in the classification of imbalanced data sets, and although CCCDs have a convenient construction in $\mathbb{R}^d$, finding minimum dominating sets is NP-hard and its probabilistic behaviour is not mathematically tractable except for $d=1$. On the other hand, a particular family of PCDs, called \emph{proportional-edge} PCDs (PE-PCDs), has mathematical tractable minimum dominating sets in $\mathbb{R}^d$; however their construction in higher dimensions may be computationally demanding. More specifically, we show that the classifiers based on PE-PCDs are prototype-based classifiers such that the exact minimum number of prototypes (equivalent to minimum dominating sets) are found in polynomial time on the number of observations. We construct two types of classifiers based on PE-PCDs. One is a family of hybrid classifiers depend on the location of the points of the training data set, and another type is a family of classifiers solely based on class covers. We assess the classification performance of our PE-PCD based classifiers by extensive Monte Carlo simulations, and compare them with that of other commonly used classifiers. We also show that, similar to CCCD classifiers, our classifiers are relatively better in classification in the presence of class imbalance.


title: The Knaster-Tarski theorem versus monotone nonexpansive mappings

abstract: Let $X$ be a partially ordered set with the property that each family of order intervals of the form $[a,b],[a,\rightarrow )$ with the finite intersection property has a nonempty intersection. We show that every directed subset of $X$ has a supremum. Then we apply the above result to prove that if $X$ is a topological space with a partial order $\preceq $ for which the order intervals are compact, $\mathcal{F}$ a nonempty commutative family of monotone maps from $X$ into $X$ and there exists $c\in X$ such that $c\preceq Tc$ for every $T\in \mathcal{F}$, then the set of common fixed points of $\mathcal{F}$ is nonempty and has a maximal element. The result, specialized to the case of Banach spaces gives a general fixed point theorem that drops almost all assumptions from the recent results in this area. An application to the theory of integral equations of Urysohn's type is also given.


title: Pairing in spherical nuclei: quasi-particle random phase approximation calculations with the Gogny interaction

abstract: We investigate the effects of the pairing in spherical nuclei. We use the same finite-range interaction of Gogny type in the three steps of our approach, Hartree-Fock, Bardeen, Cooper and Schrieffer, and quasi-particle random phase approximation calculations. We study electric and magnetic dipole, quadrupole and octuple excitations in oxygen and calcium isotopes and also in isotopes with 20 neutrons. We investigate the pairing effects on single particle energies and occupation probabilities, on the excitation energies, $B$-values and collectivity of low-lying states including the isoscalar electric dipole and the magnetic dipole excitations, and also the giant resonances. The inclusion of the pairing increases the values of the excitation energies in all the cases we have studied. In general, the effects of the pairing are too small to remarkably improve the agreement with the available experimental data.


title: Multi-output Polynomial Networks and Factorization Machines

abstract: Factorization machines and polynomial networks are supervised polynomial models based on an efficient low-rank decomposition. We extend these models to the multi-output setting, i.e., for learning vector-valued functions, with application to multi-class or multi-task problems. We cast this as the problem of learning a 3-way tensor whose slices share a common basis and propose a convex formulation of that problem. We then develop an efficient conditional gradient algorithm and prove its global convergence, despite the fact that it involves a non-convex basis selection step. On classification tasks, we show that our algorithm achieves excellent accuracy with much sparser models than existing methods. On recommendation system tasks, we show how to combine our algorithm with a reduction from ordinal regression to multi-output classification and show that the resulting algorithm outperforms simple baselines in terms of ranking accuracy.


title: Skew Howe duality and random rectangular Young tableaux

abstract: We consider the decomposition into irreducible components of the external power $\Lambda^p(\mathbb{C}^m\otimes \mathbb{C}^n)$ regarded as a $\operatorname{GL}_m\times\operatorname{GL}_n$-module. Skew Howe duality implies that the Young diagrams from each pair $(\lambda,\mu)$ which contributes to this decomposition turn out to be conjugate to each other, i.e.~$\mu=\lambda'$. We show that the Young diagram $\lambda$ which corresponds to a randomly selected irreducible component $(\lambda,\lambda')$ has the same distribution as the Young diagram which consists of the boxes with entries $\leq p$ of a random Young tableau of rectangular shape with $m$ rows and $n$ columns. This observation allows treatment of the asymptotic version of this decomposition in the limit as $m,n,p\to\infty$ tend to infinity.


title: A copula approach for dependence modeling in multivariate nonparametric time series

abstract: This paper is concerned with modeling the dependence structure of two (or more) time-series in the presence of a (possible multivariate) covariate which may include past values of the time series. We assume that the covariate influences only the conditional mean and the conditional variance of each of the time series but the distribution of the standardized innovations is not influenced by the covariate and is stable in time. The joint distribution of the time series is then determined by the conditional means, the conditional variances and the marginal distributions of the innovations, which we estimate nonparametrically, and the copula of the innovations, which represents the dependency structure. We consider a nonparametric as well as a semiparametric estimator based on the estimated residuals. We show that under suitable assumptions these copula estimators are asymptotically equivalent to estimators that would be based on the unobserved innovations. The theoretical results are illustrated by simulations and a real data example.


title: Guide Actor-Critic for Continuous Control

abstract: Actor-critic methods solve reinforcement learning problems by updating a parameterized policy known as an actor in a direction that increases an estimate of the expected return known as a critic. However, existing actor-critic methods only use values or gradients of the critic to update the policy parameter. In this paper, we propose a novel actor-critic method called the guide actor-critic (GAC). GAC firstly learns a guide actor that locally maximizes the critic and then it updates the policy parameter based on the guide actor by supervised learning. Our main theoretical contributions are two folds. First, we show that GAC updates the guide actor by performing second-order optimization in the action space where the curvature matrix is based on the Hessians of the critic. Second, we show that the deterministic policy gradient method is a special case of GAC when the Hessians are ignored. Through experiments, we show that our method is a promising reinforcement learning method for continuous controls.


title: An Equilibration Based A Posteriori Error Estimate for the Biharmonic Equation and Two Finite Element Methods

abstract: We develop an a posteriori error estimator for the Interior Penalty Discontinuous Galerkin approximation of the biharmonic equation with continuous finite elements. The error bound is based on the two-energies principle and requires the computation of an equilibrated moment tensor. The natural space for the moment tensor consists of symmetric tensor fields with continuous normal-normal components. It is known from the Hellan-Herrmann-Johnson (HHJ) mixed formulation. We propose a construction that is totally local. The procedure can also be applied to the original HHJ formulation, which directly provides an equilibrated moment tensor.


title: Stray field signatures of N\'eel textured skyrmions in Ir/Fe/Co/Pt multilayer films

abstract: Skyrmions are nanoscale spin configurations with topological properties that hold great promise for spintronic devices. Here, we establish their N\'eel texture, helicity, and size in Ir/Fe/Co/Pt multilayer films by constructing a multipole expansion to model their stray field signatures and applying it to magnetic force microscopy (MFM) images. Furthermore, the demonstrated sensitivity to inhomogeneity in skyrmion properties, coupled with a unique capability to estimate the pinning force governing dynamics, portends broad applicability in the burgeoning field of topological spin textures.

categories: physics.app-ph cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el


title: View-Invariant Recognition of Action Style Self-Dissimilarity

abstract: Self-similarity was recently introduced as a measure of inter-class congruence for classification of actions. Herein, we investigate the dual problem of intra-class dissimilarity for classification of action styles. We introduce self-dissimilarity matrices that discriminate between same actions performed by different subjects regardless of viewing direction and camera parameters. We investigate two frameworks using these invariant style dissimilarity measures based on Principal Component Analysis (PCA) and Fisher Discriminant Analysis (FDA). Extensive experiments performed on IXMAS dataset indicate remarkably good discriminant characteristics for the proposed invariant measures for gender recognition from video data.


title: Topological computation of some Stokes phenomena on the affine line

abstract: Let $\mathcal M$ be a holonomic algebraic $\mathcal D$-module on the affine line, regular everywhere including at infinity. Malgrange gave a complete description of the Fourier-Laplace transform $\widehat{\mathcal M}$, including its Stokes multipliers at infinity, in terms of the quiver of $\mathcal M$. Let $F$ be the perverse sheaf of holomorphic solutions to $\mathcal M$. By the irregular Riemann-Hilbert correspondence, $\widehat{\mathcal M}$ is determined by the enhanced Fourier-Sato transform $F^\curlywedge$ of $F$. Our aim here is to recover Malgrange's result in a purely topological way, by computing $F^\curlywedge$ using Borel-Moore cycles. In this paper, we also consider some irregular $\mathcal M$'s, like in the case of the Airy equation, where our cycles are related to steepest descent paths.


title: On the numerical solution of the elastodynamic problem by a boundary integral equation method

abstract: A numerical method for the Dirichlet initial boundary value problem for the elastic equation in the exterior and unbounded region of a smooth closed simply connected 2-dimensional domain, is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and a boundary integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the time-depended problem to a sequence of stationary boundary value problems, which are solved by a boundary layer approach resulting to a sequence of boundary integral equations of the first kind. The numerical discretization and solution are obtained by a trigonometrical quadrature method. Numerical results are included.


title: Quantum Transport in Weyl Semimetal Thin Films in the Presence of Spin-Orbit Coupled Impurities

abstract: Topological semimetals have been at the forefront of experimental and theoretical attention in condensed matter physics. Among these, recently discovered Weyl semimetals have a dispersion described by a three-dimensional Dirac cone, which is at the root of exotic physics such as the chiral anomaly in magnetotransport. In a time reversal symmetric (TRS) Weyl semimetal film, the confinement gap gives the quasiparticles a mass, while TRS is preserved by having an even number of valleys with opposite masses. The film can be tuned through a topological phase transition by a gate electric field. In this work, we present a theoretical study of the quantum corrections to the conductivity of a topological semimetal thin film, which is governed by the complex interplay of the chiral band structure, mass term, and scalar and spin-orbit scattering. We study scalar and spin-orbit scattering mechanisms on the same footing, demonstrating that they have a strong qualitative and quantitative impact on the conductivity correction. We show that, due to the spin structure of the matrix Green's functions, terms linear in the extrinsic spin-orbit scattering are present in the Bloch and momentum relaxation times, whereas previous works had identified corrections starting from the second order. In the limit of small quasiparticle mass, the terms linear in the impurity spin-orbit coupling lead to a potentially observable density dependence in the weak antilocalization correction. Moreover, when the mass term is around 30 percent of the linear Dirac terms, the system is in the unitary symmetry class with zero quantum correction and switching the extrinsic spin-orbit scattering drives the system to the weak antilocalization. We discuss the crossover between the weak localization and weak antilocalization regimes in terms of the singlet and triplet Cooperon channels, tuning the spin-orbit scattering strength.


title: Nonconvex homogenization for one-dimensional controlled random walks in random potential

abstract: We consider a finite horizon stochastic optimal control problem for nearest-neighbor random walk $\{X_i\}$ on the set of integers. The cost function is the expectation of exponential of the path sum of a random stationary and ergodic bounded potential plus $\theta X_n$. The random walk policies are measurable with respect to the random potential, and are adapted, with their drifts uniformly bounded in magnitude by a parameter $\delta\in[0,1]$. Under natural conditions on the potential, we prove that the normalized logarithm of the optimal cost function converges. The proof is constructive in the sense that we identify asymptotically optimal policies given the value of the parameter $\delta$, as well as the law of the potential. It relies on correctors from large deviation theory as opposed to arguments based on subadditivity which do not seem to work except when $\delta = 0$. The Bellman equation associated to this control problem is a second-order Hamilton-Jacobi (HJ) stochastic partial difference equation with a separable random Hamiltonian which is nonconvex in $\theta$ unless $\delta = 0$. We prove that this equation homogenizes under linear initial data to a first-order HJ deterministic partial differential equation. When $\delta = 0$, the effective Hamiltonian is the tilted free energy of random walk in random potential and it is convex in $\theta$. In contrast, when $\delta = 1$, the effective Hamiltonian is piecewise linear and nonconvex in $\theta$. Finally, when $\delta \in (0,1)$, the effective Hamiltonian is expressed completely in terms of the tilted free energy for the $\delta=0$ case and its convexity/nonconvexity in $\theta$ is characterized by a simple inequality involving $\delta$ and the magnitude of the potential, thereby marking two qualitatively distinct control regimes.


title: Feedback inhibition shapes emergent computational properties of cortical microcircuit motifs

abstract: Cortical microcircuits are very complex networks, but they are composed of a relatively small number of stereotypical motifs. Hence one strategy for throwing light on the computational function of cortical microcircuits is to analyze emergent computational properties of these stereotypical microcircuit motifs. We are addressing here the question how spike-timing dependent plasticity (STDP) shapes the computational properties of one motif that has frequently been studied experimentally: interconnected populations of pyramidal cells and parvalbumin-positive inhibitory cells in layer 2/3. Experimental studies suggest that these inhibitory neurons exert some form of divisive inhibition on the pyramidal cells. We show that this data-based form of feedback inhibition, which is softer than that of winner-take-all models that are commonly considered in theoretical analyses, contributes to the emergence of an important computational function through STDP: The capability to disentangle superimposed firing patterns in upstream networks, and to represent their information content through a sparse assembly code.


title: AIXIjs: A Software Demo for General Reinforcement Learning

abstract: Reinforcement learning is a general and powerful framework with which to study and implement artificial intelligence. Recent advances in deep learning have enabled RL algorithms to achieve impressive performance in restricted domains such as playing Atari video games (Mnih et al., 2015) and, recently, the board game Go (Silver et al., 2016). However, we are still far from constructing a generally intelligent agent. Many of the obstacles and open questions are conceptual: What does it mean to be intelligent? How does one explore and learn optimally in general, unknown environments? What, in fact, does it mean to be optimal in the general sense? The universal Bayesian agent AIXI (Hutter, 2005) is a model of a maximally intelligent agent, and plays a central role in the sub-field of general reinforcement learning (GRL). Recently, AIXI has been shown to be flawed in important ways; it doesn't explore enough to be asymptotically optimal (Orseau, 2010), and it can perform poorly with certain priors (Leike and Hutter, 2015). Several variants of AIXI have been proposed to attempt to address these shortfalls: among them are entropy-seeking agents (Orseau, 2011), knowledge-seeking agents (Orseau et al., 2013), Bayes with bursts of exploration (Lattimore, 2013), MDL agents (Leike, 2016a), Thompson sampling (Leike et al., 2016), and optimism (Sunehag and Hutter, 2015). We present AIXIjs, a JavaScript implementation of these GRL agents. This implementation is accompanied by a framework for running experiments against various environments, similar to OpenAI Gym (Brockman et al., 2016), and a suite of interactive demos that explore different properties of the agents, similar to REINFORCEjs (Karpathy, 2015). We use AIXIjs to present numerous experiments illustrating fundamental properties of, and differences between, these agents.


title: A novel algorithmic approach to Bayesian Logic Regression

abstract: Logic regression was developed more than a decade ago as a tool to construct predictors from Boolean combinations of binary covariates. It has been mainly used to model epistatic effects in genetic association studies, which is very appealing due to the intuitive interpretation of logic expressions to describe the interaction between genetic variations. Nevertheless logic regression has (partly due to computational challenges) remained less well known than other approaches to epistatic association mapping. Here we will adapt an advanced evolutionary algorithm called GMJMCMC (Genetically modified Mode Jumping Markov Chain Monte Carlo) to perform Bayesian model selection in the space of logic regression models. After describing the algorithmic details of GMJMCMC we perform a comprehensive simulation study that illustrates its performance given logic regression terms of various complexity. Specifically GMJMCMC is shown to be able to identify three-way and even four-way interactions with relatively large power, a level of complexity which has not been achieved by previous implementations of logic regression. We apply GMJMCMC to reanalyze QTL mapping data for Recombinant Inbred Lines in \textit{Arabidopsis thaliana} and from a backcross population in \textit{Drosophila} where we identify several interesting epistatic effects. The method is implemented in an R package which is available on github.


title: The role of spatial scale in joint optimisations of generation and transmission for European highly renewable scenarios

abstract: The effects of the spatial scale on the results of the optimisation of transmission and generation capacity in Europe are quantified under a 95% CO2 reduction compared to 1990 levels, interpolating between one-node-per-country solutions and many-nodes-per-country. The trade-offs that come with higher spatial detail between better exposure of transmission bottlenecks, exploitation of sites with good renewable resources (particularly wind power) and computational limitations are discussed. It is shown that solutions with no grid expansion beyond today's capacities are only around 20% more expensive than with cost-optimal grid expansion.


title: The heat and work of quantum thermodynamic processes with quantum coherence

abstract: Energy is often partitioned into heat and work by two independent paths corresponding to the change in the eigenenergies or the probability distributions of a quantum system. The discrepancies of the heat and work for various quantum thermodynamic processes have not been well characterized in literature. Here we show how the work in quantum machines is differentially related to isochoric, isothermal, and adiabatic processes. We prove that the energy exchanges during the quantum isochoric and isothermal processes are simply depending on the change in the eigenenergies or the probability distributions. However, for a time-dependent system in a non-adiabatic quantum evolution, the transitions between the different quantum states representing the quantum coherence can affect the essential thermodynamic properties, and thus the general definitions of the heat and work should be clarified with respect to the microscopic generic time-dependent system. By integrating the coherence effects in the exactly-solvable dynamics of quantum-spin precession, the internal energy is rigorously transferred as the work in the thermodynamic adiabatic process. The present study demonstrates that quantum adiabatic process is sufficient but not necessary for thermodynamic adiabatic process.


title: Dayside and nightside magnetic field responses at 780 km altitude to dayside reconnection

abstract: During southward interplanetary magnetic field, dayside reconnection will drive the Dungey cycle in the magnetosphere, which is manifested as a two-cell convection pattern in the ionosphere. We address the response of the ionospheric convection to changes in the dayside reconnection rate by examining magnetic field perturbations at 780 km altitude. The Active Magnetosphere and Planetary Electrodynamics Response Experiment data products derived from the Iridium constellation provide global maps of the magnetic field perturbations. Cluster data just upstream of the Earth's bow shock have been used to estimate the dayside reconnection rate. By using a statistical model where the magnetic field can respond on several time scales, we confirm previous reports of an almost immediate response both near noon and near midnight combined with a 10-20 min reconfiguration time of the two-cell convection pattern. The response of the ionospheric convection has been associated with the expansion of the polar cap boundary in the Cowley-Lockwood paradigm. In the original formulation of this paradigm the expansion spreads from noon to midnight in 15-20 min. However, also an immediate global response has been shown to be consistent with the paradigm when the previous dayside reconnection history is considered. In this paper we present a new explanation for how the immediate response can be accommodated in the Cowley-Lockwood paradigm. The new explanation is based on how MHD waves propagate in the magnetospheric lobes when newly reconnected open flux tubes are added to the lobes, and the magnetopause flaring angle increases.


title: Modal series expansions for plane gravitational waves

abstract: Propagation of gravitational disturbances at the speed of light is one of the key predictions of the General Theory of Relativity. This result is now backed indirectly by the observations of the behaviour of the ephemeris of binary pulsar systems. These new results have in- creased the interest in the mathematical theory of gravitational waves in the last decades and several mathematical approaches have been developed for a better understanding of the solutions. In this paper we develop a modal series expansion technique in which solutions can be built for plane waves from a seed integrable function. The convergence of these series is proven by the Raabe- Duhamel criteria and we show that these solutions are characterized by a well-defined and finite curvature tensor and also a finite energy content.


title: The grasshopper problem

abstract: We introduce and physically motivate the following problem in geometric combinatorics, originally inspired by analysing Bell inequalities. A grasshopper lands at a random point on a planar lawn of area one. It then jumps once, a fixed distance $d$, in a random direction. What shape should the lawn be to maximise the chance that the grasshopper remains on the lawn after jumping? We show that, perhaps surprisingly, a disc shaped lawn is not optimal for any $d>0$. We investigate further by introducing a spin model whose ground state corresponds to the solution of a discrete version of the grasshopper problem. Simulated annealing and parallel tempering searches are consistent with the hypothesis that for $ d < \pi^{-1/2}$ the optimal lawn resembles a cogwheel with $n$ cogs, where the integer $n$ is close to $ \pi ( \arcsin ( \sqrt{\pi} d /2 ) )^{-1}$. We find transitions to other shapes for $d \gtrsim \pi^{-1/2}$.

categories: cond-mat.stat-mech math-ph math.CO math.MP quant-ph


title: Exchange-mediated dynamic screening in the integer quantum Hall regime

abstract: We study many-body interaction effects in the spatially-resolved filling factor ($\nu$) distribution for higher Landau levels (LLs) via self-consistent Hartree-Fock simulations in the integer quantum Hall (IQH) regime. Our results indicate a strong, interaction-induced tendency to avoid the simultaneous existence of partially filled spin-up and spin-down LLs. Rather, we find that such partially filled LLs consist of coexisting regions of full and empty LLs. At the boundaries between the regions of full and empty LLs, we observe edge stripes of nearly constant $\nu$ close to half-odd filling. This suggests that the exchange interaction induces a behavior similar to a Hund's rule for the occupation of the spin split LLs. The screening of the disorder and edge potential appears significantly reduced as compared to static Thomas-Fermi screening. Our results are consistent with a local, lateral $\nu$ dependence of the exchange-enhanced spin splitting. Hence, on quantum-coherent length scales as probed here, the electron system of the IQH effect behaves similar to a non-interacting single particle system - not because of the absence, but rather due to the dominance of many-body effects.


title: Construction of strongly regular Cayley graphs based on three-valued Gauss periods

abstract: In this paper, we give a construction of strongly regular Cayley graphs on the additive groups of finite fields based on three-valued Gauss periods. As consequences, we obtain two infinite families and one sporadic example of new strongly regular Cayley graphs. This construction can be viewed as a generalization of that of strongly regular Cayley graphs obtained in \cite{BLMX}.


title: Calibration-Free Relaxation-Based Multi-Color Magnetic Particle Imaging

abstract: Magnetic Particle Imaging (MPI) is a novel imaging modality with important applications such as angiography, stem cell tracking, and cancer imaging. Recently, there have been efforts to increase the functionality of MPI via multi-color imaging methods that can distinguish the responses of different nanoparticles, or nanoparticles in different environmental conditions. The proposed techniques typically rely on extensive calibrations that capture the differences in the harmonic responses of the nanoparticles. In this work, we propose a method to directly estimate the relaxation time constant of the nanoparticles from the MPI signal, which is then used to generate a multi-color relaxation map. The technique is based on the underlying mirror symmetry of the adiabatic MPI signal when the same region is scanned back and forth. We validate the proposed method via extensive simulations, and via experiments on our in-house Magnetic Particle Spectrometer (MPS) setup at 550 Hz and our in-house MPI scanner at 9.7 kHz. Our results show that nanoparticles can be successfully distinguished with the proposed technique, without any calibration or prior knowledge about the nanoparticles.


title: Variations for Some Painlev\'e Equations

abstract: This paper first discusses irreducibility of a Painlev\'e equation $P$. We explain how the Painlev\'e property is helpful for the computation of special classical and algebraic solutions. As in a paper of Morales-Ruiz we associate an autonomous Hamiltonian $\mathbb{H}$ to a Painlev\'e equation $P$. Complete integrability of $\mathbb{H}$ is shown to imply that all solutions to $P$ are classical (which includes algebraic), so in particular $P$ is solvable by ''quadratures''. Next, we show that the variational equation of $P$ at a given algebraic solution coincides with the normal variational equation of $\mathbb{H}$ at the corresponding solution. Finally, we test the Morales-Ramis theorem in all cases $P_{2}$ to $P_{5}$ where algebraic solutions are present, by showing how our results lead to a quick computation of the component of the identity of the differential Galois group for the first two variational equations. As expected there are no cases where this group is commutative.


title: Moduli spaces of curves in tropical varieties

abstract: We describe a framework to construct tropical moduli spaces of rational stable maps to a smooth tropical hypersurface or curve. These moduli spaces will be tropical cycles of the expected dimension, corresponding to virtual fundamental classes in algebraic geometry. As we focus on the combinatorial aspect, we take the weights on certain basic 0-dimensional local combinatorial curve types as input data, and give a compatibility condition in dimension 1 to ensure that this input data glues to a global well-defined tropical cycle. As an application, we construct such moduli spaces for the case of lines in surfaces, and in a subsequent paper for stable maps to a curve.


title: Rational CFTs on Riemann surfaces

abstract: The partition function of rational conformal field theories (CFTs) on Riemann surfaces is expected to satisfy ODEs of Gauss-Manin type. We investigate the case of hyperelliptic surfaces and derive the ODE system for the $(2,5)$ minimal model.


title: JHelioviewer - Time-dependent 3D visualisation of solar and heliospheric data

abstract: Context. Solar observatories are providing the world-wide community with a wealth of data, covering large time ranges, multiple viewpoints, and returning large amounts of data. In particular, the large volume of SDO data presents challenges: it is available only from a few repositories, and full-disk, full-cadence data for reasonable durations of scientific interest are difficult to download practically due to their size and download data rates available to most users. From a scientist's perspective this poses three problems: accessing, browsing and finding interesting data as efficiently as possible. Aims. To address these challenges, we have developed JHelioviewer, a visualisation tool for solar data based on the JPEG2000 compression standard and part of the open source ESA/NASA Helioviewer Project. Since the first release of JHelioviewer, the scientific functionality of the software has been extended significantly, and the objective of this paper is to highlight these improvements. Methods. The JPEG2000 standard offers useful new features that facilitate the dissemination and analysis of high-resolution image data and offers a solution to the challenge of efficiently browsing petabyte-scale image archives. The JHelioviewer software is open source, platform independent and extendable via a plug-in architecture. Results. With JHelioviewer, users can visualise the Sun for any time period between September 1991 and today. They can perform basic image processing in real time, track features on the Sun and interactively overlay magnetic field extrapolations. The software integrates solar event data and a time line display. As a first step towards supporting science planning of the upcoming Solar Orbiter mission, JHelioviewer offers a virtual camera model that enables users to set the vantage point to the location of a spacecraft or celestial body at any given time.


title: Renormalized quark-antiquark Hamiltonian induced by a gluon mass ansatz in heavy-flavor QCD

abstract: In response to the growing need for theoretical tools that can be used in QCD to describe and understand the dynamics of gluons in hadrons in the Minkowski space-time, the renormalization group procedure for effective particles (RGPEP) is shown in the simplest available context of heavy quarkonia to exhibit a welcome degree of universality in the first approximation it yields once one assumes that beyond perturbation theory gluons obtain effective mass. Namely, in the second-order terms, the Coulomb potential with Breit-Fermi spin couplings in the effective quark-antiquark component of a heavy quarkonium, is corrected in one-flavor QCD by a spin-independent harmonic oscillator term that does not depend on the assumed effective gluon mass or the choice of the RGPEP generator. The new generator we use here is much simpler than the ones used before and has the advantage of being suitable for studies of the effective gluon dynamics at higher orders than the second and beyond the perturbative expansion.


title: Almost Gorenstein Hibi rings

abstract: In this paper, we state criteria of a Hibi ring to be level, non-Gorenstein and almost Gorenstein and to be non-level and almost Gorenstein in terms of the structure of the partially ordered set defining the Hibi ring. We also state a criterion of a ladder determinantal ring defined by 2-minors to be non-Gorenstein and almost Gorenstein in terms of the shape of the ladder.


title: An infrared view of AGN feedback in a type-2 quasar: the case of the Teacup galaxy

abstract: We present near-infrared integral field spectroscopy data obtained with VLT/SINFONI of "the Teacup galaxy". The nuclear K-band (1.95-2.45 micron) spectrum of this radio-quiet type-2 quasar reveals a blueshifted broad component of FWHM~1600-1800 km/s in the hydrogen recombination lines (Pa$\alpha$, Br$\delta$, and Br$\gamma$) and also in the coronal line [Si VI]$\lambda$1.963 micron. Thus the data confirm the presence of the nuclear ionized outflow previously detected in the optical and reveal its coronal counterpart. Both the ionized and coronal nuclear outflows are resolved, with seeing-deconvolved full widths at half maximum of 1.1$\pm$0.1 and 0.9$\pm$0.1 kpc along PA$\sim$72-74 deg. This orientation is almost coincident with the radio axis (PA=77 deg), suggesting that the radio jet could have triggered the nuclear outflow. In the case of the H$_2$ lines we do not require a broad component to reproduce the profiles, but the narrow lines are blueshifted by ~50 km/s on average from the galaxy systemic velocity. This could be an indication of the presence of a nuclear molecular outflow, although the bulk of the H$_2$ emission in the inner ~2 arcsec (~3 kpc) of the galaxy follows a rotation pattern. We find evidence for kinematically disrupted gas (FWHM>250 km/s) at up to 5.6 kpc from the AGN, which can be naturally explained by the action of the outflow. The narrow component of [Si VI] is redshifted with respect to the systemic velocity, unlike any other emission line in the K-band spectrum. This indicates that the region where the coronal lines are produced is not co-spatial with the narrow line region.


title: Computer vision-based food calorie estimation: dataset, method, and experiment

abstract: Computer vision has been introduced to estimate calories from food images. But current food image data sets don't contain volume and mass records of foods, which leads to an incomplete calorie estimation. In this paper, we present a novel food image data set with volume and mass records of foods, and a deep learning method for food detection, to make a complete calorie estimation. Our data set includes 2978 images, and every image contains corresponding each food's annotation, volume and mass records, as well as a certain calibration reference. To estimate calorie of food in the proposed data set, a deep learning method using Faster R-CNN first is put forward to detect the food. And the experiment results show our method is effective to estimate calories and our data set contains adequate information for calorie estimation. Our data set is the first released food image data set which can be used to evaluate computer vision-based calorie estimation methods.


title: Dissipatively driven strongly interacting bosons in a gauge field

abstract: The interplay between dissipation, interactions and gauge fields opens the possibility to rich emerging physics. Here we focus on a set-up in which the system is coupled at its extremities to two different baths which impose a current. We then study the system's response to a gauge field depending on the filling. We show that while the current induced by the baths has a marked dependence on the magnetic field at low fillings which is significantly reduced close to half-filling. We explain the interplay between interactions, gauge field and dissipation by studying the system's energy spectrum at the different fillings. This interplay also results in the emergence of negative differential conductivity. For this study we have developed a number-conserving treatment which allows a numerical exact treatment of fairly large system sizes, and which can be extended to a large class of systems.


title: Resonant transport and near-field effects in photonic glasses

abstract: A fundamental quantity in multiple scattering is the transport mean free path the inverse of which describes the scattering strength of a sample. In this paper, we emphasize the importance of an appropriate description of the effective refractive index $n_{\mathrm{eff}}$ in multiple light scattering to accurately describe the light transport in dense photonic glasses. Using $n_{\mathrm{eff}}$ as calculated by the energy-density coherent potential approximation we are able to predict the transport mean free path of monodisperse photonic glasses. This model without any fit parameter is in qualitative agreement with numerical simulations and in fair quantitative agreement with spectrally resolved coherent backscattering measurements on new specially synthesized polystyrene photonic glasses. These materials exhibit resonant light scattering perturbed by strong near-field coupling, all captured within the model. Our model might be used to maximize the scattering strength of high index photonic glasses, which are a key in the search for Anderson localization of light in three dimensions.


title: On the Efficient Simulation of the Left-Tail of the Sum of Correlated Log-normal Variates

abstract: The sum of Log-normal variates is encountered in many challenging applications such as in performance analysis of wireless communication systems and in financial engineering. Several approximation methods have been developed in the literature, the accuracy of which is not ensured in the tail regions. These regions are of primordial interest wherein small probability values have to be evaluated with high precision. Variance reduction techniques are known to yield accurate, yet efficient, estimates of small probability values. Most of the existing approaches, however, have considered the problem of estimating the right-tail of the sum of Log-normal random variables (RVS). In the present work, we consider instead the estimation of the left-tail of the sum of correlated Log-normal variates with Gaussian copula under a mild assumption on the covariance matrix. We propose an estimator combining an existing mean-shifting importance sampling approach with a control variate technique. The main result is that the proposed estimator has an asymptotically vanishing relative error which represents a major finding in the context of the left-tail simulation of the sum of Log-normal RVs. Finally, we assess by various simulation results the performances of the proposed estimator compared to existing estimators.


title: Wide subcategories are semistable

abstract: For an arbitrary finite dimensional algebra $\Lambda$, we prove that any wide subcategory of $\mathsf{mod} \Lambda$ satisfying a certain finiteness condition is $\theta$-semistable for some stability condition $\theta$. More generally, we show that wide subcategories of $\mathsf{mod} \Lambda$ associated with two-term presilting complexes of $\Lambda$ are semistable. This provides a complement for Ingalls-Thomas-type bijections for finite dimensional algebras.


title: Kinodynamic Planning on Constraint Manifolds

abstract: This paper presents a motion planner for systems subject to kinematic and dynamic constraints. The former appear when kinematic loops are present in the system, such as in parallel manipulators, in robots that cooperate to achieve a given task, or in situations involving contacts with the environment. The latter are necessary to obtain realistic trajectories, taking into account the forces acting on the system. The kinematic constraints make the state space become an implicitly-defined manifold, which complicates the application of common motion planning techniques. To address this issue, the planner constructs an atlas of the state space manifold incrementally, and uses this atlas both to generate random states and to dynamically simulate the steering of the system towards such states. The resulting tools are then exploited to construct a rapidly-exploring random tree (RRT) over the state space. To the best of our knowledge, this is the first randomized kinodynamic planner for implicitly-defined state spaces. The test cases presented in this paper validate the approach in significantly-complex systems.


title: $K^0-\bar{K}^0$ mixing in the minimal flavor-violating two-Higgs-doublet models

abstract: The two-Higgs-doublet model (2HDM), as one of the simplest extensions of the Standard Model (SM), is obtained by adding another scalar doublet to the SM, and is featured by a pair of charged scalars, which could affect many low-energy processes. In the "Higgs basis" for a generic 2HDM, only one scalar doublet gets a nonzero vacuum expectation value and, under the criterion of minimal flavor violation, the other one is fixed to be either color-singlet or color-octet, which are named as the type-III and the type-C 2HDM, respectively. In this paper, we study the charged-scalar effects of these two models on the $K^0-\bar{K}^0$ mixing, an ideal process to probe New Physics (NP) beyond the SM. Firstly, we perform a complete one-loop computation of the box diagrams relevant to the $K^0-\bar{K}^0$ mixing, keeping the mass and momentum of the external strange quark up to the second order. Together with the up-to-date theoretical inputs, we then give a detailed phenomenological analysis, in the cases of both real and complex Yukawa couplings of the charged scalars to quarks. The parameter spaces allowed by the current experimental data on the mass difference $\Delta m_K$ and the CP-violating parameter $\epsilon_K$ are obtained and the differences between these two 2HDMs are investigated, which are helpful to distinguish them from each other from a phenomenological point of view.


title: Weakening of the diamagnetic shielding in FeSe$_{1-x}$S$_x$ at high pressures

abstract: The superconducting transition of FeSe$_{1-x}$S$_x$ with three distinct sulphur concentrations $x$ was studied under hydrostatic pressure up to $\sim$70 kbar via bulk AC susceptibility. The pressure dependence of the superconducting transition temperature ($T_c$) features a small dome-shaped variation at low pressures for $x=0.04$ and $x=0.12$, followed by a more substantial $T_c$ enhancement to a value of around 30 K at moderate pressures. In $x=0.21$, a similar overall pressure dependence of $T_c$ is observed, except that the small dome at low pressures is flattened. For all three concentrations, a significant weakening of the diamagnetic shielding is observed beyond the pressure around which the maximum $T_c$ of 30 K is reached near the verge of pressure-induced magnetic phase. This observation points to a strong competition between the magnetic and high-$T_c$ superconducting states at high pressure in this system.


title: Dynamics Based 3D Skeletal Hand Tracking

abstract: Tracking the full skeletal pose of the hands and fingers is a challenging problem that has a plethora of applications for user interaction. Existing techniques either require wearable hardware, add restrictions to user pose, or require significant computation resources. This research explores a new approach to tracking hands, or any articulated model, by using an augmented rigid body simulation. This allows us to phrase 3D object tracking as a linear complementarity problem with a well-defined solution. Based on a depth sensor's samples, the system generates constraints that limit motion orthogonal to the rigid body model's surface. These constraints, along with prior motion, collision/contact constraints, and joint mechanics, are resolved with a projected Gauss-Seidel solver. Due to camera noise properties and attachment errors, the numerous surface constraints are impulse capped to avoid overpowering mechanical constraints. To improve tracking accuracy, multiple simulations are spawned at each frame and fed a variety of heuristics, constraints and poses. A 3D error metric selects the best-fit simulation, helping the system handle challenging hand motions. Such an approach enables real-time, robust, and accurate 3D skeletal tracking of a user's hand on a variety of depth cameras, while only utilizing a single x86 CPU core for processing.


title: Speed and fluctuations for some driven dimer models

abstract: We consider driven dimer models on the square and honeycomb graphs, starting from a stationary Gibbs measure. Each model can be thought of as a two dimensional stochastic growth model of an interface, belonging to the anisotropic KPZ universality class. We use a combinatorial approach to determine the speed of growth and show logarithmic growth in time of the variance of the height function fluctuations.


title: From optimal transport to generative modeling: the VEGAN cookbook

abstract: We study unsupervised generative modeling in terms of the optimal transport (OT) problem between true (but unknown) data distribution $P_X$ and the latent variable model distribution $P_G$. We show that the OT problem can be equivalently written in terms of probabilistic encoders, which are constrained to match the posterior and prior distributions over the latent space. When relaxed, this constrained optimization problem leads to a penalized optimal transport (POT) objective, which can be efficiently minimized using stochastic gradient descent by sampling from $P_X$ and $P_G$. We show that POT for the 2-Wasserstein distance coincides with the objective heuristically employed in adversarial auto-encoders (AAE) (Makhzani et al., 2016), which provides the first theoretical justification for AAEs known to the authors. We also compare POT to other popular techniques like variational auto-encoders (VAE) (Kingma and Welling, 2014). Our theoretical results include (a) a better understanding of the commonly observed blurriness of images generated by VAEs, and (b) establishing duality between Wasserstein GAN (Arjovsky and Bottou, 2017) and POT for the 1-Wasserstein distance.


title: Near-Feasible Stable Matchings with Budget Constraints

abstract: We consider the matching with contracts framework of Hatfield and Milgrom when one side (a firm or hospital) can make monetary transfers (offer wages) to the other (a worker or doctor). In a standard model, monetary transfers are not restricted. However, we assume that each hospital has a fixed budget; that is, the total amount of wages allocated by each hospital to the doctors is constrained. With this constraint, stable matchings may fail to exist and checking for the existence is hard. To deal with the nonexistence, we focus on near-feasible matchings that can exceed each hospital budget by a certain amount, and We introduce a new concept of compatibility. We show that the compatibility condition is a sufficient condition for the existence of a near-feasible stable matching in the matching with contracts framework. Under a slight restriction on hospitals' preferences, we provide mechanisms that efficiently return a near-feasible stable matching with respect to the actual amount of wages allocated by each hospital. By sacrificing strategy-proofness, the best possible bound of budget excess is achieved.