title: Seismic Negative Belt of Acoustic Metamaterials

abstract: An earthquake-proof seismic negative belt of an artificial seismic shadow zone is introduced. The belt is composed of acoustic materials which has one of the constituent parameter between density and modulus is negative effectively. It converts the velocity of the seismic wave imaginary, and then creates a stop-band for the seismic frequency range. The belt is an attenuator of a seismic wave that reduces the amplitude of the wave exponentially. Passing the belt underground, the seismic energy turns into sound and heat in air and the magnitude of the seismic wave is weakened to be defended by conventional method. Models of the negative density and negative modulus for engineering are suggested.


title: The Prototypical Young L/T-Transition Dwarf HD 203030B Likely Has Planetary Mass

abstract: Upon its discovery in 2006, the young L7.5 companion to the solar analog HD 203030 was found to be unusual in being $\approx$200 K cooler than older late-L dwarfs. HD 203030B offered the first clear indication that the effective temperature at the L-to-T spectral type transition depends on surface gravity: now a well-known characteristic of low-gravity ultra-cool dwarfs. An initial age analysis of the G8V primary star indicated that the system was 130--400 Myr old, and so the companion between 12--31 $M_{\rm Jup}$. Using moderate resolution near-infrared spectra of HD 203030B, we now find features of very low gravity comparable to those of 10--150 Myr-old L7--L8 dwarfs. We also obtained more accurate near infrared and {\sl Spitzer}/IRAC photometry, and find a $(J-K)_{\rm MKO}$ color of $2.56\pm0.13$ mag---comparable to those observed in other young planetary-mass objects---and a luminosity of log($L_{\rm bol}/L_{\odot}$)$\,=\,-4.75\pm0.04$ dex. We further reassess the evidence for the young age of the host star, HD 203030, with a more comprehensive analysis of the photometry and updated stellar activity measurements and age calibrations. Summarizing the age diagnostics for both components of the binary, we adopt an age of 100 Myr for HD 203030B and an age range of 30--150 Myr. Using cloudy evolutionary models, the new companion age range and luminosity result in a mass of 11 $M_{\rm Jup}$ with a range of 8--15 $M_{\rm Jup}$, and an effective temperature of $1040\pm50$ K.


title: Estimates for Sums of Eigenvalues of the Free Plate via the Fourier Transform

abstract: Using the Fourier transform, we obtain upper bounds for sums of eigenvalues of the free plate.


title: Synchronization in Networks of Diffusively Coupled Nonlinear Systems: Robustness Against Time-Delays

abstract: In this manuscript, we study the problem of robust synchronization in networks of diffusively time-delayed coupled nonlinear systems. In particular, we prove that, under some mild conditions on the input-output dynamics of the systems and the network topology, there always exists a unimodal region in the parameter space (coupling strength versus time-delay), such that if they belong to this region, the systems synchronize. Moreover, we show how this unimodal region scales with the network topology, which, in turn, provides useful insights on how to design the network topology to maximize robustness against time-delays. The results are illustrated by extensive simulation experiments of time-delayed coupled Hindmarsh-Rose neural chaotic oscillators.


title: Adversarial Advantage Actor-Critic Model for Task-Completion Dialogue Policy Learning

abstract: This paper presents a new method --- adversarial advantage actor-critic (Adversarial A2C), which significantly improves the efficiency of dialogue policy learning in task-completion dialogue systems. Inspired by generative adversarial networks (GAN), we train a discriminator to differentiate responses/actions generated by dialogue agents from responses/actions by experts. Then, we incorporate the discriminator as another critic into the advantage actor-critic (A2C) framework, to encourage the dialogue agent to explore state-action within the regions where the agent takes actions similar to those of the experts. Experimental results in a movie-ticket booking domain show that the proposed Adversarial A2C can accelerate policy exploration efficiently.


title: Approximating Continuous Functions by ReLU Nets of Minimal Width

abstract: This article concerns the expressive power of depth in deep feed-forward neural nets with ReLU activations. Specifically, we answer the following question: for a fixed $d_{in}\geq 1,$ what is the minimal width $w$ so that neural nets with ReLU activations, input dimension $d_{in}$, hidden layer widths at most $w,$ and arbitrary depth can approximate any continuous, real-valued function of $d_{in}$ variables arbitrarily well? It turns out that this minimal width is exactly equal to $d_{in}+1.$ That is, if all the hidden layer widths are bounded by $d_{in}$, then even in the infinite depth limit, ReLU nets can only express a very limited class of functions, and, on the other hand, any continuous function on the $d_{in}$-dimensional unit cube can be approximated to arbitrary precision by ReLU nets in which all hidden layers have width exactly $d_{in}+1.$ Our construction in fact shows that any continuous function $f:[0,1]^{d_{in}}\to\mathbb R^{d_{out}}$ can be approximated by a net of width $d_{in}+d_{out}$. We obtain quantitative depth estimates for such an approximation in terms of the modulus of continuity of $f$.

categories: stat.ML cs.CC cs.LG math.CO math.ST stat.TH


title: The Importance of Forgetting: Limiting Memory Improves Recovery of Topological Characteristics from Neural Data

abstract: We develop of a line of work initiated by Curto and Itskov towards understanding the amount of information contained in the spike trains of hippocampal place cells via topology considerations. Previously, it was established that simply knowing which groups of place cells fire together in an animal's hippocampus is sufficient to extract the global topology of the animal's physical environment. We model a system where collections of place cells group and ungroup according to short-term plasticity rules. In particular, we obtain the surprising result that in experiments with spurious firing, the accuracy of the extracted topological information decreases with the persistence (beyond a certain regime) of the cell groups. This suggests that synaptic transience, or forgetting, is a mechanism by which the brain counteracts the effects of spurious place cell activity.


title: Improving mechanical sensor performance through larger damping

abstract: Mechanical resonances are used in a wide variety of devices; from smart phone accelerometers to computer clocks and from wireless communication filters to atomic force microscope sensors. Frequency stability, a critical performance metric, is generally assumed to be tantamount to resonance quality factor (the inverse of the linewidth and of the damping). Here we show that frequency stability of resonant nanomechanical sensors can generally be made independent of quality factor. At high bandwidths, we show that quality factor reduction is completely mitigated by increases in signal to noise ratio. At low bandwidths, strikingly, increased damping leads to better stability and sensor resolution, with improvement proportional to damping. We confirm the findings by demonstrating temperature resolution of 50 \mu K at 200 Hz bandwidth. These results open the door for high performance ultrasensitive resonant sensors in gaseous or liquid environments, single cell nanocalorimetry, nanoscale gas chromatography, and atmospheric pressure nanoscale mass spectrometry.


title: Notes on Cops and Robber game on graphs

abstract: These are some personal notes about the pursuit game of Cops and Robbers that I made starting in 2007. More old and new problems (and some solutions) will be added in future versions of these notes.


title: Uniform analytic approximation of Wigner rotation matrices

abstract: We derive the leading asymptotic approximation, for low angle {\theta}, of the Wigner rotation matrix elements $d^j_{m_1m_2}(\theta)$, uniform in $j,m_1$ and $m_2$. The result is in terms of a Bessel function of integer order. We numerically investigate the error for a variety of cases and find that the approximation can be useful over a significant range of angles. This approximation has application in the partial wave analysis of wavepacket scattering.


title: Macroeconomics and FinTech: Uncovering Latent Macroeconomic Effects on Peer-to-Peer Lending

abstract: Peer-to-peer (P2P) lending is a fast growing financial technology (FinTech) trend that is displacing traditional retail banking. Studies on P2P lending have focused on predicting individual interest rates or default probabilities. However, the relationship between aggregated P2P interest rates and the general economy will be of interest to investors and borrowers as the P2P credit market matures. We show that the variation in P2P interest rates across grade types are determined by three macroeconomic latent factors formed by Canonical Correlation Analysis (CCA) - macro default, investor uncertainty, and the fundamental value of the market. However, the variation in P2P interest rates across term types cannot be explained by the general economy.


title: Some regularity and convergence results for parabolic Hamilton-Jacobi-Bellman equations in bounded domains

abstract: We study the approximation of parabolic Hamilton-Jacobi-Bellman (HJB) equations in bounded domains with strong Dirichlet boundary conditions. We work under the assumption of the existence of a sufficiently regular barrier function for the problem to obtain well-posedness and regularity of a related switching system and the convergence of its components to the HJB equation. In particular, we show existence of a viscosity solution to the switching system by a novel construction of sub- and supersolutions and application of Perron's method. Error bounds for monotone schemes for the HJB equation are then derived from estimates near the boundary, where the standard regularisation procedure for viscosity solutions is not applicable, and are found to be of the same order as known results for the whole space. We deduce error bounds for some common finite difference and truncated semi-Lagrangian schemes.


title: Dissipative extension theory for linear relations

abstract: This work is devoted to dissipative extension theory for dissipative linear relations. We give a self-consistent theory of extensions by generalizing the theory on symmetric extensions of symmetric operators. Several results on the properties of dissipative relations are proven. Finally, we deal with the spectral properties of dissipative extensions of dissipative relations and provide results concerning particular realizations of this general setting.


title: Consistency of Generalized Dynamic Principal Components in Dynamic Factor Models

abstract: We study the theoretical properties of the generalized dynamic principal components introduced in Pe\~na and Yohai (2016). In particular, we prove that when the data follows a dynamic factor model, the reconstruction provided by the procedure converges in mean square to the common part of the model as the number of series and periods diverge to infinity. The results of a simulation study support our findings.


title: Asymptotic behavior as $p\rightarrow\infty$ of least energy solutions of a $(p,q(p))$-Laplacian problem

abstract: \[ \left\{ \begin{array} [c]{lll} -\left( \Delta_{p}+\Delta_{q(p)}\right) u=\lambda_{p}\left\vert u(x_{u})\right\vert ^{p-2}u(x_{u})\delta_{x_{u}} & \mathrm{in} & \Omega\\ u=0 & \mathrm{on} & \partial\Omega, \end{array} \right. \] where $x_{u}$ is the (unique) maximum point of $\left\vert u\right\vert ,$ $\delta_{x_{u}}$ is the Dirac delta distribution supported at $x_{u},$ \[ \lim_{p\rightarrow\infty}\frac{q(p)}{p}=Q\in\left\{ \begin{array} [c]{lll} (0,1) & \mathrm{if} & N<q(p)<p\\ (1,\infty) & \mathrm{if} & N<p<q(p) \end{array} \right. \] and $\lambda_{p}>0$ is such that \[ \min\left\{ \frac{\left\Vert \nabla u\right\Vert _{\infty}}{\left\Vert u\right\Vert _{\infty}}:0\not \equiv u\in W^{1,\infty}(\Omega)\cap C_{0}(\overline{\Omega})\right\} \leq\lim_{p\rightarrow\infty}(\lambda _{p})^{\frac{1}{p}}<\infty. \]


title: Affine highest weight categories and quantum affine Schur-Weyl duality of Dynkin quiver types

abstract: For a Dynkin quiver $Q$ (of type ADE), we consider a central completion of the convolution algebra of the equivariant K-group of a certain Steinberg type graded quiver variety. We observe that it is affine quasi-hereditary and prove that its category of finite-dimensional modules is identified with a block of Hernandez-Leclerc's monoidal category $\mathcal{C}_Q$ of modules over the quantum loop algebra $U_q(L\mathfrak{g})$ via Nakajima's homomorphism. As an application, we show that Kang-Kashiwara-Kim's generalized quantum affine Schur-Weyl duality functor gives an equivalence between the category of finite-dimensional modules over the quiver Hecke algebra associated with $Q$ and Hernandez-Leclerc's category $\mathcal{C}_Q$, assuming the simpleness of some poles of normalized R-matrices for type E.


title: Quasi-local charges and the Generalized Gibbs Ensemble in the Lieb-Liniger model

abstract: We consider the construction of a generalized Gibbs ensemble composed of complete bases of conserved charges in the repulsive Lieb-Liniger model. We will show that it is possible to construct these bases with varying locality as well as demonstrating that such constructions are always possible provided one has in hand at least one complete basis set of charges. This procedure enables the construction of bases of charges that possess well defined, finite expectation values given an arbitrary initial state. We demonstrate the use of these charges in the context of two different quantum quenches: a quench where the strength of the interactions in a one-dimensional gas is switched suddenly from zero to some finite value and the release of a one dimensional cold atomic gas from a confining parabolic trap. While we focus on the Lieb-Liniger model in this paper, the principle of the construction of these charges applies to all integrable models, both in continuum and lattice form.


title: Collider Detection of Dark Matter Electromagnetic Anapole Moments

abstract: Dark matter that interacts with the Standard Model by exchanging photons through higher multipole interactions occurs in a wide range of both strongly as well as weakly coupled hidden sector models. We study the collider detection prospects of these candidates, with a focus on Majorana dark matter that couples through the anapole moment. The study is conducted at the effective field theory level with the mono-$Z$ signature incorporating varying levels of systematic uncertainties at the high-luminosity LHC. The projected collider reach on the anapole moment is then compared to the reach coming from direct detection experiments like LZ. Finally, the analysis is applied to a weakly coupled completion with leptophilic dark matter.


title: A Generalized Model for Light Transport in Scintillators

abstract: Transported light in the medium usually shows as an exponential decay tendency. In the DAMPE strip scintillators, however, the phenomenon of light attenuation as the hit position approaches the end of the scintillator can not be described by the simple exponential decay model. The spread angle of PMT relative to hit position is distance-dependent, so the larger the angle, the larger the proportion $\eta$ of emitted light to becomes the effective input light. We consider the contribution of the spread angle, and propose a generalized model: $f(x)=I_0\eta_{a}(x)e^{-x/\lambda}+I_0\eta_{b}(x)e^{-(2L-x)/\lambda }$. The model well describes the light attenuation in the scintillator, reducing the maximum deviation of the sample from the fit function from 29\% to below 2\%. Moreover, our model contains most of the traditional models, so the experimental data that traditional models can fit and our models fit well.


title: Correlation between Multivariate Datasets, from Inter-Graph Distance computed using Graphical Models Learnt With Uncertainties

abstract: We present a method for simultaneous Bayesian learning of the correlation matrix and graphical model of a multivariate dataset, along with uncertainties in each, to subsequently compute distance between the learnt graphical models of a pair of datasets, using a new metric that approximates an uncertainty-normalised Hellinger distance between the posterior probabilities of the graphical models given the respective dataset; correlation between the pair of datasets is then computed as a corresponding affinity measure. We achieve a closed-form likelihood of the between-columns correlation matrix by marginalising over the between-row matrices. This between-columns correlation is updated first, given the data, and the graph is then updated, given the partial correlation matrix that is computed given the updated correlation, allowing for learning of the 95$\%$ Highest Probability Density credible regions of the correlation matrix and graphical model of the data. Difference made to the learnt graphical model, by acknowledgement of measurement noise, is demonstrated on a small simulated dataset, while the large human disease-symptom network--with $>8,000$ nodes--is learnt using real data. Data on vino-chemical attributes of Portuguese red and white wine samples are employed to learn with-uncertainty graphical model of each dataset, and subsequently, the distance between these learnt graphical models.


title: Aleman-Richter-Sundberg's Theorem On $P^t(\mu )$-Spaces

abstract: Let $\nu$ be a finite complex measure with support in $\bar {\mathbb D}$ and let $\mathcal C\nu$ denote the Cauchy transform of $\nu .$ Suppose that $\nu$ annihilates polynomials in complex variable $z$ and $\nu |_{\partial \mathbb D} = hm,$ where $m$ is the normalized Lebesgue measure on $\partial {\mathbb D}$. We show that, for $\epsilon_0 > 0,$ $m$-almost all $e^{i\theta}\in \partial {\mathbb D},$ and $a > 0,$ when $r$ tends to 1, there exists $E_r \subset B(re^{i\theta}, \frac{1-r}{4})$ with analytic capacity $\gamma (E_r) < \epsilon_0 \frac{1-r}{4}$ such that $|\mathcal C\nu (\lambda) - e^{-i\theta}h(e^{i\theta}) | \le a$ area-almost all $\lambda \in B (re^{i\theta}, \frac{1-r}{4} ) \setminus E_r .$ Using this result, we provide an alternative proof of Aleman-Richter-Sundberg's Theorem on nontangential limits in $P^t(\mu )$-Spaces and the index of invariant subspaces.


title: A Positive-Definite Energy Functional for Axially Symmetric Maxwell's Equations on Kerr-de Sitter Black Hole Spacetimes

abstract: We prove that there exists a phase space of canonical variables for the initial value problem for axially symmetric Maxwell fields propagating in Kerr-de Sitter black hole spacetimes such that their motion is restricted to the level sets of a positive-definite Hamiltonian, despite the ergo-region.


title: Optimal Control of Connected and Automated Vehicles at Roundabouts: An Investigation in a Mixed-Traffic Environment

abstract: Connectivity and automation in vehicles provide the most intriguing opportunity for enabling users to better monitor transportation network conditions and make better operating decisions to improve safety and reduce pollution, energy consumption, and travel delays. This paper investigates the implications of optimally coordinating vehicles that are wirelessly connected to each other in roundabouts to achieve a smooth traffic flow. We apply an optimization framework and an analytical solution that allows optimal coordination of vehicles for merging in such traffic scenario. The effectiveness of the proposed approach is validated through simulation and it is shown that fully coordinated vehicles reduce total travel time by 51% and fuel consumption by 35%. Furthermore, we study the influence of vehicle coordination in a mixed-traffic environment and compare the network performance under different market penetration rates of connected and automated vehicles (CAVs). For this particular study with near-capacity demand, due to extremely unstable traffic, the results show that even with high penetration of CAVs (e.g., 80%), travel time and fuel savings are much less than a network of CAVs.


title: Early pathogen replacement in a model of Influenza and Respiratory Syncytial Virus with partial vaccination. A computational study

abstract: In this paper, we carry out a computational study using the spectral decomposition of the fluctuations of a two-pathogen epidemic model around its deterministic attractor, i.e., steady state or limit cycle, to examine the role of partial vaccination and between-host pathogen interaction on early pathogen replacement during seasonal epidemics of influenza and respiratory syncytial virus.


title: Sequential Adaptive Detection for In-Situ Transmission Electron Microscopy (TEM)

abstract: We develop new efficient online algorithms for detecting transient sparse signals in TEM video sequences, by adopting the recently developed framework for sequential detection jointly with online convex optimization [1]. We cast the problem as detecting an unknown sparse mean shift of Gaussian observations, and develop adaptive CUSUM and adaptive SSRS procedures, which are based on likelihood ratio statistics with post-change mean vector being online maximum likelihood estimators with $\ell_1$. We demonstrate the meritorious performance of our algorithms for TEM imaging using real data.


title: Effective Tensor Sketching via Sparsification

abstract: In this paper, we investigate effective sketching schemes via sparsification for high dimensional multilinear arrays or tensors. More specifically, we propose a novel tensor sparsification algorithm that retains a subset of the entries of a tensor in a judicious way, and prove that it can attain a given level of approximation accuracy in terms of tensor spectral norm with a much smaller sample complexity when compared with existing approaches. In particular, we show that for a $k$th order $d\times\cdots\times d$ cubic tensor of {\it stable rank} $r_s$, the sample size requirement for achieving a relative error $\varepsilon$ is, up to a logarithmic factor, of the order $r_s^{1/2} d^{k/2} /\varepsilon$ when $\varepsilon$ is relatively large, and $r_s d /\varepsilon^2$ and essentially optimal when $\varepsilon$ is sufficiently small. It is especially noteworthy that the sample size requirement for achieving a high accuracy is of an order independent of $k$. To further demonstrate the utility of our techniques, we also study how higher order singular value decomposition (HOSVD) of large tensors can be efficiently approximated via sparsification.


title: Properties of Caratheodory measure hyperbolic universal covers of compact Kahler manifolds

abstract: This article explores some properties of universal covers of compact Kahler manifolds, under the assumption of Caratheodory measure hyperbolicity. In particular, by comparing invariant volume forms, an inequality is established between the volume of canonical bundle of a compact Kahler manifolds and the Caratheodory measure of its universal cover (similar result as in [Kikuta 10]). Using similar method, an inequality is established between the restricted volume of canonical bundle of a compact Kahler manifolds and the restricted Caratheodory measure of its covering, solving a conjecture in [Kikuta 13].


title: Strong Edge Magnetism and Tunable Energy Gaps in Zigzag Graphene Quantum Dots with High Stability

abstract: Graphene is a nonmagnetic semimetal and cannot be directly used as electronic or spintronic devices. We demonstrate that graphene quantum dots (GQDs) can exhibit strong edge magnetism and tunable energy gaps due to the presence of localized edge states. By using large-scale first principle density functional theory (DFT) calculations and detailed analysis based on model Hamiltonians, we can show that the zigzag edge states in GQDs become much stronger and more localized as the system size increases. The enhanced edge states induce strong electron-electron interactions along the edges of GQDs, ultimately resulting in a magnetic phase transition from nonmagnetic to intra-edge ferromagnetic and inter-edge antiferromagnetic, when the diameter is larger than 4.5 nm. Our analysis shows that the inter-edge superexchange interaction of antiferromagnetic states between two nearest-neighbor zigzag edges in GQDs is much stronger than that exists between two parallel zigzag edges in GQDs and graphene nanoribbons. Furthermore, such strong and localized edge states also induce GQDs semiconducting with tunable energy gaps, mainly controlled by adjusting the system size. Our results show that the quantum confinement effect, inter-edge superexchange (antiferromagnetic), and intra-edge direct exchange (ferromagnetic) interactions are crucial for the electronic and magnetic properties of zigzag GQDs at the nanoscale.


title: A generalized parsing framework for Abstract Grammars

abstract: This technical report presents a general framework for parsing a variety of grammar formalisms. We develop a grammar formalism, called an Abstract Grammar, which is general enough to represent grammars at many levels of the hierarchy, including Context Free Grammars, Minimalist Grammars, and Generalized Context-free Grammars. We then develop a single parsing framework which is capable of parsing grammars which are at least up to GCFGs on the hierarchy. Our parsing framework exposes a grammar interface, so that it can parse any particular grammar formalism that can be reduced to an Abstract Grammar.


title: Stochastic Linear Quadratic Optimal Control with General Control Domain

abstract: This paper considers the stochastic linear quadratic optimal control problem in which the control domain is nonconvex. By the functional analysis and convex perturbation methods, we establish a novel maximum principle. The application of the proposed maximum principle is illustrated through a work-out example.


title: Algorithmic learning of probability distributions from random data in the limit

abstract: We study the problem of identifying a probability distribution for some given randomly sampled data in the limit, in the context of algorithmic learning theory as proposed recently by Vinanyi and Chater. We show that there exists a computable partial learner for the computable probability measures, while by Bienvenu, Monin and Shen it is known that there is no computable learner for the computable probability measures. Our main result is the characterization of the oracles that compute explanatory learners for the computable (continuous) probability measures as the high oracles. This provides an analogue of a well-known result of Adleman and Blum in the context of learning computable probability distributions. We also discuss related learning notions such as behaviorally correct learning and orther variations of explanatory learning, in the context of learning probability distributions from data.


title: Characterizing the structural diversity of complex networks across domains

abstract: The structure of complex networks has been of interest in many scientific and engineering disciplines over the decades. A number of studies in the field have been focused on finding the common properties among different kinds of networks such as heavy-tail degree distribution, small-worldness and modular structure and they have tried to establish a theory of structural universality in complex networks. However, there is no comprehensive study of network structure across a diverse set of domains in order to explain the structural diversity we observe in the real-world networks. In this paper, we study 986 real-world networks of diverse domains ranging from ecological food webs to online social networks along with 575 networks generated from four popular network models. Our study utilizes a number of machine learning techniques such as random forest and confusion matrix in order to show the relationships among network domains in terms of network structure. Our results indicate that there are some partitions of network categories in which networks are hard to distinguish based purely on network structure. We have found that these partitions of network categories tend to have similar underlying functions, constraints and/or generative mechanisms of networks even though networks in the same partition have different origins, e.g., biological processes, results of engineering by human being, etc. This suggests that the origin of a network, whether it's biological, technological or social, may not necessarily be a decisive factor of the formation of similar network structure. Our findings shed light on the possible direction along which we could uncover the hidden principles for the structural diversity of complex networks.


title: A probabilistic cellular automata model for the dynamics of a population driven by logistic growth and weak Allee effect

abstract: We propose and investigate a one-parameter probabilistic mixture of one-dimensional elementary cellular automata under the guise of a model for the dynamics of a single-species unstructured population with nonoverlapping generations in which individuals have smaller probability of reproducing and surviving in a crowded neighbourhood but also suffer from isolation and dispersal. Remarkably, the first-order mean field approximation to the dynamics of the model yields a cubic map containing terms representing both logistic and weak Allee effects. The model has a single absorbing state devoid of individuals, but depending on the reproduction and survival probabilities can achieve a stable population. We determine the critical probability separating these two phases and find that the phase transition between them is in the directed percolation universality class of critical behaviour.


title: The Exact Solution to Rank-1 L1-norm TUCKER2 Decomposition

abstract: We study rank-1 {L1-norm-based TUCKER2} (L1-TUCKER2) decomposition of 3-way tensors, treated as a collection of $N$ $D \times M$ matrices that are to be jointly decomposed. Our contributions are as follows. i) We prove that the problem is equivalent to combinatorial optimization over $N$ antipodal-binary variables. ii) We derive the first two algorithms in the literature for its exact solution. The first algorithm has cost exponential in $N$; the second one has cost polynomial in $N$ (under a mild assumption). Our algorithms are accompanied by formal complexity analysis. iii) We conduct numerical studies to compare the performance of exact L1-TUCKER2 (proposed) with standard HOSVD, HOOI, GLRAM, PCA, L1-PCA, and TPCA-L1. Our studies show that L1-TUCKER2 outperforms (in tensor approximation) all the above counterparts when the processed data are outlier corrupted.


title: Generalized Forward-Backward Splitting with Penalization for Monotone Inclusion Problems

abstract: We introduce a generalized forward-backward splitting method with penalty term for solving monotone inclusion problems involving the sum of a finite number of maximally monotone operators and the normal cone to the nonempty set of zeros of another maximal monotone operator. We show weak ergodic convergence of the generated sequence of iterates to a solution of the considered monotone inclusion problem, provided the condition corresponded to the Fitzpatrick function of the operator describing the set of the normal cone is fulfilled. Under strong monotonicity of an operator, we show strong convergence of the iterates. Furthermore, we utilize the proposed method for minimizing a large-scale hierarchical minimization problem concerning the sum of differentiable and nondifferentiable convex functions subject to the set of minima of another differentiable convex function. We illustrate the functionality of the method through numerical experiments addressing constrained elastic net and generalized Heron location problems.


title: Bistable cooling in optomechanical system

abstract: A scheme is presented to optimize the optomechanical cooling of mechanical resonator in instability regime. Based on the stability analysis, we uncovered a distinct bistable effect of photons and phonons, which can be used to realize a strong nonlinear effect even in the single-photon weak coupling regime. Considering the experimental realization, we investigate the sideband cooling in bistable regime with and without quantum nonlinearity. It is shown that the fluctuation of the steady state phonons can be excellently suppressed at a rather low level due to the anti-rotating-wave effect, and it does not require high quality factor of the cavity. Our scheme offers a new perspective for optimizing the sideband cooling of mechanical resonators in the weak coupling regime.


title: Tumor Classification and Segmentation of MR Brain Images

abstract: The diagnosis and segmentation of tumors using any medical diagnostic tool can be challenging due to the varying nature of this pathology. Magnetic Reso- nance Imaging (MRI) is an established diagnostic tool for various diseases and disorders and plays a major role in clinical neuro-diagnosis. Supplementing this technique with automated classification and segmentation tools is gaining importance, to reduce errors and time needed to make a conclusive diagnosis. In this paper a simple three-step algorithm is proposed; (1) identification of patients that present with tumors, (2) automatic selection of abnormal slices of the patients, and (3) segmentation and detection of the tumor. Features were extracted by using discrete wavelet transform on the normalized images and classified by support vector machine (for step (1)) and random forest (for step (2)). The 400 subjects were divided in a 3:1 ratio between training and test with no overlap. This study is novel in terms of use of data, as it employed the entire T2 weighted slices as a single image for classification and a unique combination of contralateral approach with patch thresholding for segmentation, which does not require a training set or a template as is used by most segmentation studies. Using the proposed method, the tumors were segmented accurately with a classification accuracy of 95% with 100% specificity and 90% sensitivity.


title: An Innovations Approach to Viterbi Decoding of Convolutional Codes

abstract: We introduce the notion of innovations for Viterbi decoding of convolutional codes. First we define a kind of innovation corresponding to the received data, i.e., the input to a Viterbi decoder. Then the structure of a Scarce-State-Transition (SST) Viterbi decoder is derived in a natural manner. It is shown that the newly defined innovation is just the input to the main decoder in an SST Viterbi decoder and generates the same syndrome as the original received data does. A similar result holds for Quick-Look-In (QLI) codes as well. In this case, however, the precise innovation is not defined. We see that this innovation-like quantity is related to the linear smoothed estimate of the information. The essence of innovations approach to a linear filtering problem is first to whiten the observed data, and then to treat the resulting simpler white-noise observations problem. In our case, this corresponds to the reduction of decoding complexity in the main decoder in an SST Viterbi decoder. We show the distributions related to the main decoder (i.e., the input distribution and the state distribution in the code trellis for the main decoder) are much biased under moderately noisy conditions. We see that these biased distributions actually lead to the complexity reduction in the main decoder. Furthermore, it is shown that the proposed innovations approach can be extended to maximum-likelihood (ML) decoding of block codes as well.


title: Deep Forward and Inverse Perceptual Models for Tracking and Prediction

abstract: We consider the problems of learning forward models that map state to high-dimensional images and inverse models that map high-dimensional images to state in robotics. Specifically, we present a perceptual model for generating video frames from state with deep networks, and provide a framework for its use in tracking and prediction tasks. We show that our proposed model greatly outperforms standard deconvolutional methods and GANs for image generation, producing clear, photo-realistic images. We also develop a convolutional neural network model for state estimation and compare the result to an Extended Kalman Filter to estimate robot trajectories. We validate all models on a real robotic system.


title: A Gagliardo-Nirenberg-type inequality and its applications to decay estimates for solutions of a degenerate parabolic equation

abstract: We establish a Gagliardo-Nirenberg-type inequality in $\mathbb{R}^n$ for functions which decay fast as $|x|\to\infty$. We use this inequality to derive upper bounds for the decay rates of solutions of a degenerate parabolic equation. Moreover, we show that these upper bounds, hence also the Gagliardo-Nirenberg-type inequality, are sharp in an appropriate sense.


title: Lost but not Forgotten: Intracluster Light in Galaxy Groups and Clusters

abstract: With Hubble Space Telescope imaging, we investigate the progenitor population and formation mechanisms of the intracluster light (ICL) for 23 galaxy groups and clusters ranging from 3$\times10^{13}<$M$_{500,c}$ [M$_\odot$]$<9\times10^{14}$ at 0.29$<$z$<$0.89. The color gradients of the BCG+ICL become bluer with increasing radius out to 53-100 kpc for all but one system, suggesting that violent relaxation after major mergers with the BCG cannot be the dominant source of ICL. For clusters the BCG+ICL luminosity at r$<$100 kpc (0.08-0.13 r$_{500,c}$) is 1.2-3.5$\times 10^{12}$L$_\odot$; for the groups, BCG+ICL luminosities within 100 kpc (0.17-0.23 r$_{500,c}$) range between 0.7-1.3$\times 10^{12}$ L$_\odot$. The BCG+ICL stellar mass in the inner 100 kpc increases with total cluster mass as M$_\bigstar\propto$M$_{500,c}$$^{0.37\pm0.05}$. This steep slope implies that the BCG+ICL is a higher fraction of the total mass in groups than in clusters. The BCG+ICL luminosities and stellar masses are too large for the ICL stars to come from the dissolution of dwarf galaxies alone, implying instead that the ICL grows from the stripping of more massive galaxies. Using the colors of cluster members from the CLASH sample, we place conservative lower limits on the luminosities of galaxies from which the ICL could originate. We find that at 10 kpc the ICL has a color similar to massive, passive cluster galaxies ($>10^{11.6}$ M$_\odot$), while by 100 kpc this colour is equivalent to that of a 10$^{10}$ M$_\odot$ galaxy. Additionally, we find 75% of the total BCG+ICL luminosity is consistent in color of galaxies with L$>$0.2 L$_*$ (log(M$_\bigstar$[M$_\odot$])$>$10.4), assuming conservatively that these galaxies are completely disrupted. We conclude that tidal stripping of massive galaxies is the likely source of the intracluster light from 10-100 kpc (0.008-0.23 r$_{500,c}$) for galaxy groups and clusters.


title: Vanishing Ideals Of Affine Sets Parameterized By Odd Cycles

abstract: Let K be a finite field. Let X* be a subset of the affine space Kn, which is parameterized by odd cycles. In this paper we give an explicit Gr\"obner basis for the vanishing ideal, I(X*), of X*. We give an explicit formula for the regularity of I(X*) and finally if X* is parameterized by an odd cycle of length k, we show that the Hilbert function of the vanishing ideal of X* can be written as linear combination of Hilbert functions of degenerate torus.


title: Rate-optimal Meta Learning of Classification Error

abstract: Meta learning of optimal classifier error rates allows an experimenter to empirically estimate the intrinsic ability of any estimator to discriminate between two populations, circumventing the difficult problem of estimating the optimal Bayes classifier. To this end we propose a weighted nearest neighbor (WNN) graph estimator for a tight bound on the Bayes classification error; the Henze-Penrose (HP) divergence. Similar to recently proposed HP estimators [berisha2016], the proposed estimator is non-parametric and does not require density estimation. However, unlike previous approaches the proposed estimator is rate-optimal, i.e., its mean squared estimation error (MSEE) decays to zero at the fastest possible rate of $O(1/M+1/N)$ where $M,N$ are the sample sizes of the respective populations. We illustrate the proposed WNN meta estimator for several simulated and real data sets.


title: The spin-1/2 Heisenberg antiferromagnet on the pyrochlore lattice: An exact diagonalisation study

abstract: We present exact diagonalisation calculations for the spin-1/2 nearest neighbour antiferromagnet on the pyrochlore lattice. We study a section of the lattice in the $[111]$ direction and analyse the Hamiltonian of the breathing pyrochlore system with two coupling constants $J_1$ and $J_2$ for tetrahedra of different orientations and investigate the evolution of the system from the limit of disconnected tetrahedra ($J_2 = 0$) to a correlated state at $J_1=J_2$. We evaluate the low energy spectrum, two and four spin correlations, and spin chirality correlations for a system size of up to 36 sites. The model shows a fast decay of spin correlations and we confirm the presence of several singlet excitations below the lowest magnetic excitation. We find chirality correlations near $J_1=J_2$ to be small at the length scales available at this system size. Evaluation of dimer-dimer correlations and analysis of the nature of the entanglement of the tetrahedral unit shows that the triplet sector of the tetrahedron contributes significantly to the ground state entanglement at $J_1 = J_2$.


title: Nebula: F0 Estimation and Voicing Detection by Modeling the Statistical Properties of Feature Extractors

abstract: A F0 and voicing status estimation algorithm for high quality speech analysis/synthesis is proposed. This problem is approached from a different perspective that models the behavior of feature extractors under noise, instead of directly modeling speech signals. Under time-frequency locality assumptions, the joint distribution of extracted features and target F0 can be characterized by training a bank of Gaussian mixture models (GMM) on artificial data generated from Monte-Carlo simulations. The trained GMMs can then be used to generate a set of conditional distributions on the predicted F0, which are then combined and post-processed by Viterbi algorithm to give a final F0 trajectory. Evaluation on CSTR and CMU Arctic speech databases shows that the proposed method, trained on fully synthetic data, achieves lower gross error rates than state-of-the-art methods.


title: Data-based prediction and causality inference of nonlinear dynamics

abstract: Natural systems are typically nonlinear and complex, and it is of great interest to be able to reconstruct a system in order to understand its mechanism, which can not only recover nonlinear behaviors but also predict future dynamics. Due to the advances of modern technology, big data becomes increasingly accessible and consequently the problem of reconstructing systems from measured data or time series plays a central role in many scientific disciplines. In recent decades, nonlinear methods rooted in state space reconstruction have been developed, and they do not assume any model equations but can recover the dynamics purely from the measured time series data. In this review, the development of state space reconstruction techniques will be introduced and the recent advances in systems prediction and causality inference using state space reconstruction will be presented. Particularly, the cutting-edge method to deal with short-term time series data will be focused. Finally, the advantages as well as the remaining problems in this field are discussed.


title: Learning Motion Predictors for Smart Wheelchair using Autoregressive Sparse Gaussian Process

abstract: Constructing a smart wheelchair on a commercially available powered wheelchair (PWC) platform avoids a host of seating, mechanical design and reliability issues but requires methods of predicting and controlling the motion of a device never intended for robotics. Analog joystick inputs are subject to black-box transformations which may produce intuitive and adaptable motion control for human operators, but complicate robotic control approaches; furthermore, installation of standard axle mounted odometers on a commercial PWC is difficult. In this work, we present an integrated hardware and software system for predicting the motion of a commercial PWC platform that does not require any physical or electronic modification of the chair beyond plugging into an industry standard auxiliary input port. This system uses an RGB-D camera and an Arduino interface board to capture motion data, including visual odometry and joystick signals, via ROS communication. Future motion is predicted using an autoregressive sparse Gaussian process model. We evaluate the proposed system on real-world short-term path prediction experiments. Experimental results demonstrate the system's efficacy when compared to a baseline neural network model.


title: Suspension and simple optical characterization of two-dimensional membranes

abstract: We report on a method for suspending two-dimensional crystal materials in an electronic circuit using an only photoresists and solvents. Graphene and NbSe${}_2$ are suspended tens of nanometers above metal electrodes with clamping diameters of several microns. The optical cavity formed from the membrane/air/metal structures enables a quick method to measure the number of layers and the gap separation using comparisons between the expected colour and optical microscope images. This characterization technique can be used with just an illuminated microscope with a digital camera which makes it adaptable to environments where other means of characterization are not possible, such as inside nitrogen glove boxes used in handling oxygen-sensitive materials.


title: Existence of Kirillov-Reshetikhin crystals of type $G_2^{(1)}$ and $D_4^{(3)}$

abstract: In this paper we prove that every Kirillov-Reshetikhin module of type $G_2^{(1)}$ and $D_4^{(3)}$ has a crystal pseudobase (crystal base modulo signs), by applying the criterion for the existence of a crystal pseudobase due to Kang et al.


title: High-cadence, High-resolution Spectroscopic Observations of Herbig Stars HD 98922 and V1295 Aquila

abstract: Recent observational work has indicated that mechanisms for accretion and outflow in Herbig Ae/Be star-disk systems may differ from magnetospheric accretion (MA) as it is thought to occur in T Tauri star-disk systems. In this work, we assess the temporal evolution of spectral lines probing accretion and mass loss in Herbig Ae/Be systems and test for consistency with the MA paradigm. For two Herbig Ae/Be stars, HD 98922 (B9e) and V1295 Aql (A2e), we have gathered multi-epoch (~years) and high-cadence (~minutes) high-resolution optical spectra to probe a wide range of kinematic processes. Employing a line equivalent width evolution correlation metric introduced here, we identify species co-evolving (indicative of common line origin) via novel visualization. We interferometrically constrain often problematically degenerate parameters, inclination and inner disk radius, allowing us to focus on the structure of the wind, magnetosphere, and inner gaseous disk in radiative transfer models. Over all timescales sampled, the strongest variability occurs within the blueshifted absorption components of the Balmer series lines; the strength of variability increases with the cadence of the observations. Finally, high-resolution spectra allow us to probe substructure within the Balmer series' blueshifted absorption components: we observe static, low-velocity features and time-evolving features at higher velocities. Overall, we find the observed line morphologies and variability are inconsistent with a scaled-up T Tauri MA scenario. We suggest that as magnetic field structure and strength change dramatically with increasing stellar mass from T Tauri to Herbig Ae/Be stars, so too may accretion and outflow processes.


title: The Kontsevich--Zorich cocycle over Veech--McMullen family of symmetric translation surfaces

abstract: We describe the Kontsevich--Zorich cocycle over an affine invariant orbifold coming from a (cyclic) covering construction inspired by works of Veech and McMullen. In particular, using the terminology in a recent paper of Filip, we show that all cases of Kontsevich--Zorich monodromies of $SU(p,q)$ type are realized by appropriate covering constructions.


title: Resolution and Relevance Trade-offs in Deep Learning

abstract: Deep learning has been successfully applied to various tasks, but its underlying mechanism remains unclear. Neural networks associate similar inputs in the visible layer to the same state of hidden variables in deep layers. The fraction of inputs that are associated to the same state is a natural measure of similarity and is simply related to the cost in bits required to represent these inputs. The degeneracy of states with the same information cost provides instead a natural measure of noise and is simply related the entropy of the frequency of states, that we call relevance. Representations with minimal noise, at a given level of similarity (resolution), are those that maximise the relevance. A signature of such efficient representations is that frequency distributions follow power laws. We show, in extensive numerical experiments, that deep neural networks extract a hierarchy of efficient representations from data, because they i) achieve low levels of noise (i.e. high relevance) and ii) exhibit power law distributions. We also find that the layer that is most efficient to reliably generate patterns of training data is the one for which relevance and resolution are traded at the same price, which implies that frequency distribution follows Zipf's law.


title: Statistical Studies of Solar White-Light Flares and Comparisons with Superflares on Solar-type Stars

abstract: Recently, many superflares on solar-type stars have been discovered as white-light flares (WLFs). The statistical study found a correlation between their energies ($E$) and durations ($\tau$): $\tau \propto E^{0.39}$ (Maehara et al. 2017 $EP\& S$, 67, 59), similar to those of solar hard/soft X-ray flares: $\tau \propto E^{0.2-0.33}$. This indicates a universal mechanism of energy release on solar and stellar flares, i.e., magnetic reconnection. We here carried out a statistical research on 50 solar WLFs observed with \textit{SDO}/HMI and examined the correlation between the energies and durations. As a result, the $E$--$\tau$ relation on solar WLFs ($\tau \propto E^{0.38}$) is quite similar to that on stellar superflares ($\tau \propto E^{0.39}$). However, the durations of stellar superflares are one order of magnitude shorter than those expected from solar WLFs. We present the following two interpretations for the discrepancy. (1) In solar flares, the cooling timescale of WLFs may be longer than the reconnection one, and the decay time of solar WLFs can be elongated by the cooling effect. (2) The distribution can be understood by applying a scaling law ($\tau \propto E^{1/3}B^{-5/3}$) derived from the magnetic reconnection theory. In this case, the observed superflares are expected to have 2-4 times stronger magnetic field strength than solar flares.


title: Geometry of spin coherent states

abstract: Spin states of maximal projection along some direction in space are called (spin) coherent, and are, in many aspects, the "most classical" available. For any spin $s$, the spin coherent states form a 2-sphere in the projective Hilbert space $\mathbb{P}$ of the system. We address several questions regarding that sphere, in particular its possible intersections with complex lines. We also find that, like Dali's iconic clocks, it extends in all possible directions in $\mathbb{P}$. We give a simple expression for the Majorana constellation of the linear combination of two coherent states, and use Mason's theorem to give a lower bound on the number of distinct stars of a linear combination of two arbitrary spin-$s$ states. Finally, we plot the image of the spin coherent sphere, assuming light in $\mathbb{P}$ propagates along Fubini-Study geodesics. We argue that, apart from their intrinsic geometric interest, such questions translate into statements experimentalists might find useful.


title: The incompatibility of crossing number and bridge number for knot diagrams

abstract: We define and compare several natural ways to compute the bridge number of a knot diagram. We study bridge numbers of crossing number minimizing diagrams, as well as the behavior of diagrammatic bridge numbers under the connected sum operation. For each notion of diagrammatic bridge number considered, we find crossing number minimizing knot diagrams which fail to minimize bridge number. Furthermore, we construct a family of minimal crossing diagrams for which the difference between diagrammatic bridge number and the actual bridge number of the knot grows to infinity.


title: Gaussian Approximation of the Distribution of Strongly Repelling Particles on the Unit Circle

abstract: In this paper, we consider a strongly-repelling model of $n$ ordered particles $\{e^{i \theta_j}\}_{j=0}^{n-1}$ with the density $p({\theta_0},\cdots, \theta_{n-1})=\frac{1}{Z_n} \exp \left\{-\frac{\beta}{2}\sum_{j \neq k} \sin^{-2} \left( \frac{\theta_j-\theta_k}{2}\right)\right\}$, $\beta>0$. Let $\theta_j=\frac{2 \pi j}{n}+\frac{x_j}{n^2}+const$ such that $\sum_{j=0}^{n-1}x_j=0$. Define $\zeta_n \left( \frac{2 \pi j}{n}\right) =\frac{x_j}{\sqrt{n}}$ and extend $\zeta_n$ piecewise linearly to $[0, 2 \pi]$. We prove the functional convergence of $\zeta_n(t)$ to $\zeta(t)=\sqrt{\frac{2}{\beta}} \mathfrak{Re} \left( \sum_{k=1}^{\infty} \frac{1}{k} e^{ikt} Z_k \right)$, where $Z_k$ are i.i.d. complex standard Gaussian random variables.


title: Modelling the cosmic spectral energy distribution and extragalactic background light over all time

abstract: We present a phenomological model of the Cosmic Spectral Energy Distribution (CSED) and the integrated galactic light (IGL) over all cosmic time. This model, based on an earlier model by Driver et al. (2013), attributes the cosmic star formation history to two processes -- firstly, chaotic clump accretion and major mergers, resulting in the early-time formation of bulges and secondly, cold gas accretion, resulting in late-time disc formation. Under the assumption of a Universal Chabrier initial mass function, we combine the Bruzual & Charlot (2003) stellar libraries, the Charlot & Fall (2000) dust attenuation prescription and template spectra for emission by dust and active galactic nuclei to predict the CSED -- pre- and post-dust attenuation -- and the IGL throughout cosmic time. The phenomological model, as constructed, adopts a number of basic axioms and empirical results and has minimal free parameters. We compare the model output, as well as predictions from the semi-analytic model GALFORM to recent estimates of the CSED out to $z=1$. By construction, our empirical model reproduces the full energy output of the Universe from the ultraviolet to the far-infrared extremely well. We use the model to derive predictions of the stellar and dust mass densities, again finding good agreement. We find that GALFORM predicts the CSED for $z < 0.3$ in good agreement with the observations. This agreement becomes increasingly poor towards $z = 1$, when the model CSED is $\sim$50 per cent fainter. The latter is consistent with the model underpredicting the cosmic star formation history. As a consequence, GALFORM predicts a $\sim$30 per cent fainter IGL.


title: Homogenization of the variational principle for discrete random maps

abstract: We consider homogenization of random surfaces and study the variational principle for graph homomorphisms from subsets of $\mathbb{Z}^m$ into $\mathbb{Z}$, where the underlying uniform measure is perturbed by a random field. Motivated by the theories of random walks in random potentials, we assume that random field is stationary, ergodic, and bounded in $L^1$ . We show that the variational principle holds in probability and that the entropy functional homogenizes, i.e.\ is independent of the values taken by the random field. The main ingredients in the argument are the existence of the quenched surface tension, the equivalence of the quenched and the annealed surface tension, and robustness of the surface tension under change in boundary data. These ingredients are deduced by a combination of a superadditive ergodic theorem and combinatorial results, especially the Kirszbraun theorem.


title: Refined Heinz Mean Operator Inequality

abstract: It is shown that if $A,B\in \mathbb{B}\left( \mathcal{H} \right)$ be positive operators, then \begin{equation*} \begin{aligned} A\#B&\le \frac{1}{1-2\mu }{A^{\frac{1}{2}}}{{F}_{\mu }}\left( {A^{-\frac{1}{2}}}B{A^{-\frac{1}{2}}} \right){A^{\frac{1}{2}}}\\ & \le \frac{1}{2}\left[ A\#B+{{H}_{\mu }}\left( A,B \right) \right]\\ & \le \frac{1}{2}\left[ \frac{1}{1-2\mu }{A^{\frac{1}{2}}} {F_{\mu }}\left( {A^{-\frac{1}{2}}}B{A^{-\frac{1}{2}}} \right){A^{\frac{1}{2}}}+{H_{\mu }}\left( A,B \right) \right]\\ & \le \cdots \le \frac{1}{{{2}^{n}}}A\#B+\frac{{{2}^n}-1}{{2^n}}{H_\mu }\left( A,B \right)\\ & \le \frac{1}{{{2}^{n}}\left( 1-2\mu \right)}{A^{\frac{1}{2}}}{F_{\mu }}\left( {A^{-\frac{1}{2}}}B{A^{-\frac{1}{2}}} \right){A^{\frac{1}{2}}}+\frac{{{2}^{n}}-1}{{{2}^{n}}}{H_\mu }\left( A,B \right)\\ & \le \frac{1}{{2^{n+1}}}A\#B+\frac{{{2}^{n+1}}-1}{{2^{n+1}}}{H_\mu }\left( A,B \right)\\ & \le \cdots \le {H_\mu }\left( A,B \right). \end{aligned} \end{equation*} for each $\mu \in \left[ 0,1 \right]\backslash \left\{ \frac{1}{2} \right\}$. As an application, we present several inequalities for unitarily invariant norms. Our results are refinements of some existing inequalities.


title: Improving Social Media Text Summarization by Learning Sentence Weight Distribution

abstract: Recently, encoder-decoder models are widely used in social media text summarization. However, these models sometimes select noise words in irrelevant sentences as part of a summary by error, thus declining the performance. In order to inhibit irrelevant sentences and focus on key information, we propose an effective approach by learning sentence weight distribution. In our model, we build a multi-layer perceptron to predict sentence weights. During training, we use the ROUGE score as an alternative to the estimated sentence weight, and try to minimize the gap between estimated weights and predicted weights. In this way, we encourage our model to focus on the key sentences, which have high relevance with the summary. Experimental results show that our approach outperforms baselines on a large-scale social media corpus.


title: Spectral geometries on a compact metric space

abstract: The notion of a spectral geometry on a compact metric space X is introduced. This notion serves as a discrete approximation of X motivated by the notion of a spectral triple from non-commutative geometry. A set of axioms charaterising spectral geometries is given. Bounded deformations of spectral geometries are studied and the relationship between the dimension of a spectral geometry and more traditional dimensions of metric spaces is investigated.


title: Shallow Discourse Parsing with Maximum Entropy Model

abstract: In recent years, more research has been devoted to studying the subtask of the complete shallow discourse parsing, such as indentifying discourse connective and arguments of connective. There is a need to design a full discourse parser to pull these subtasks together. So we develop a discourse parser turning the free text into discourse relations. The parser includes connective identifier, arguments identifier, sense classifier and non-explicit identifier, which connects with each other in pipeline. Each component applies the maximum entropy model with abundant lexical and syntax features extracted from the Penn Discourse Tree-bank. The head-based representation of the PDTB is adopted in the arguments identifier, which turns the problem of indentifying the arguments of discourse connective into finding the head and end of the arguments. In the non-explicit identifier, the contextual type features like words which have high frequency and can reflect the discourse relation are introduced to improve the performance of non-explicit identifier. Compared with other methods, experimental results achieve the considerable performance.


title: A Superconducting Phase Shifter and Traveling Wave Kinetic Inductance Parametric Amplifier for W-Band Astronomy

abstract: The W-Band ($75-110\; \mathrm{GHz}$) sky contains a plethora of information about star formation, galaxy evolution and the cosmic microwave background. We have designed and fabricated a dual-purpose superconducting circuit to facilitate the next generation of astronomical observations in this regime by providing proof-of-concept for both a millimeter-wave low-loss phase shifter, which can operate as an on-chip Fourier transform spectrometer (FTS) and a traveling wave kinetic inductance parametric amplifier (TKIP). Superconducting transmission lines have a propagation speed that depends on the inductance in the line which is a combination of geometric inductance and kinetic inductance in the superconductor. The kinetic inductance has a non-linear component with a characteristic current, $I_*$, and can be modulated by applying a DC current, changing the propagation speed and effective path length. Our test circuit is designed to measure the path length difference or phase shift, $\Delta \phi$, between two symmetric transmission lines when one line is biased with a DC current. To provide a measurement of $\Delta\phi$, a key parameter for optimizing a high gain W-Band TKIP, and modulate signal path length in FTS operation, our $3.6 \times 2.5\; \mathrm{cm}$ chip employs a pair of $503\; \mathrm{mm}$ long NbTiN inverted microstrip lines coupled to circular waveguide ports through radial probes. For a line of width $3\; \mathrm{\mu m}$ and film thickness $20\; \mathrm{nm}$, we predict $\Delta\phi\approx1767\; \mathrm{rad}$ at $90\; \mathrm{GHz}$ when biased at close to $I_*$. We have fabricated a prototype with $200\; \mathrm{nm}$ thick Nb film and the same line length and width. The predicted phase shift for our prototype is $\Delta\phi\approx30\; \mathrm{rad}$ at $90\; \mathrm{GHz}$ when biased at close to $I_*$ for Nb.


title: Local and global strong solutions to the stochastic incompressible Navier-Stokes equations in critical Besov space

abstract: Considering the stochastic Navier-Stokes system in $\mathbb{R}^d$ forced by a multiplicative white noise, we establish the local existence and uniqueness of the strong solution when the initial data take values in the critical space $\dot{B}_{p,r}^{\frac{d}{p}-1}(\mathbb{R}^d)$. The proof is based on the contraction mapping principle, stopping time and stochastic estimates. Then we prove the global existence of strong solutions in probability if the initial data are sufficiently small, which contain a class of highly oscillating "large" data.


title: Building models for extended radio sources: implications for Epoch of Reionisation science

abstract: We test the hypothesis that limitations in the sky model used to calibrate an interferometric radio telescope, where the model contains extended radio sources, will generate bias in the Epoch of Reionisation (EoR) power spectrum. The information contained in a calibration model about the spatial and spectral structure of an extended source is incomplete because a radio telescope cannot sample all Fourier components. Application of an incomplete sky model to calibration of EoR data will imprint residual error in the data, which propagates forward to the EoR power spectrum. This limited information is studied in the context of current and future planned instruments and surveys at EoR frequencies, such as the Murchison Widefield Array (MWA), Giant Metrewave Radio Telescope (GMRT) and the Square Kilometre Array (SKA1-Low). For the MWA EoR experiment, we find that both the additional short baseline $uv$-coverage of the compact EoR array, and the additional long baselines provided by TGSS and planned MWA expansions, are required to obtain sufficient information on all relevant scales. For SKA1-Low, arrays with maximum baselines of 49~km and 65~km yield comparable performance at 50~MHz and 150~MHz, while 39~km, 14~km and 4~km arrays yield degraded performance.


title: Complementarity and the pathological statistics of the quantum impossible

abstract: Standard quantum physics prevents the existence of a joint statistics for complementary observables. Nevertheless, a joint distribution for complementary observables can be derived from their imperfect simultaneous measurement followed by a suitable data inversion. Quantumness is reflected then in the pathologies of the inferred distribution. We apply this program to the paradigmatic example of complementarity: the wave-particle duality in a Young interferometer.


title: Entropy production for partially observed harmonic systems

abstract: The probability distribution of the total entropy production in the non-equilibrium steady state follows a symmetry relation called the fluctuation theorem. When a certain part of the system is masked or hidden, it is difficult to infer the exact estimate of the total entropy production. Entropy produced from the observed part of the system shows a significant deviation from the steady state fluctuation theorem. This deviation occurs due to the interaction between the observed and the masked part of the system. A naive guess would be that the deviation from the steady state fluctuation theorem may disappear in the limit of small interaction between both parts of the system. In contrast, we investigate the entropy production of a particle in a harmonically coupled Brownian particle system (say, particle A and B) in a heat reservoir at a constant temperature. The system is maintained in the non-equilibrium steady state using stochastic driving. When the coupling between particle A and B is infinitesimally weak, the deviation from the steady state fluctuation theorem for the entropy production of a partial system of a coupled system is studied. Furthermore, we consider a harmonically confined system (i.e., a harmonically coupled system of particle A and B in harmonic confinement). In the weak coupling limit, the entropy produced by the partial system (e.g., particle A) of the coupled system in a harmonic trap satisfies the steady state fluctuation theorem. Numerical simulations are performed to support the analytical results. Part of these results was announced in a recent letter, Europhys. Lett. 115, 60003 (2016).


title: Secular Instabilities of Keplerian Stellar Discs

abstract: We present idealized models of a razor-thin, axisymmetric, Keplerian stellar disc around a massive black hole, and study non-axisymmetric secular instabilities in the absence of either counter-rotation or loss cones. These discs are prograde mono-energetic waterbags, whose phase space distribution functions are constant for orbits within a range of eccentricities (e) and zero outside this range. The linear normal modes of waterbags are composed of sinusoidal disturbances of the edges of distribution function in phase space. Waterbags which include circular orbits (polarcaps) have one stable linear normal mode for each azimuthal wavenumber m. The m = 1 mode always has positive pattern speed and, for polarcaps consisting of orbits with e < 0.9428, only the m = 1 mode has positive pattern speed. Waterbags excluding circular orbits (bands) have two linear normal modes for each m, which can be stable or unstable. We derive analytical expressions for the instability condition, pattern speeds, growth rates and normal mode structure. Narrow bands are unstable to modes with a wide range in m. Numerical simulations confirm linear theory and follow the non-linear evolution of instabilities. Long-time integration suggests that instabilities of different m grow, interact non-linearly and relax collisionlessly to a coarse-grained equilibrium with a wide range of eccentricities.


title: Global Rank Estimation

abstract: In real world complex networks, the importance of a node depends on two important parameters: 1. characteristics of the node, and 2. the context of the given application. The current literature contains several centrality measures that have been defined to measure the importance of a node based on the given application requirements. These centrality measures assign a centrality value to each node that denotes its importance index. But in real life applications, we are more interested in the relative importance of the node that can be measured using its centrality rank based on the given centrality measure. To compute the centrality rank of a node, we need to compute the centrality value of all the nodes and compare them to get the rank. This process requires the entire network. So, it is not feasible for real-life applications due to the large size and dynamic nature of real world networks. In the present project, we aim to propose fast and efficient methods to estimate the global centrality rank of a node without computing the centrality value of all nodes. These methods can be further extended to estimate the rank without having the entire network. The proposed methods use the structural behavior of centrality measures, sampling techniques, or the machine learning models. In this work, we also discuss how to apply these methods for degree and closeness centrality rank estimation.


title: Generating Natural Adversarial Examples

abstract: Due to their complex nature, it is hard to characterize the ways in which machine learning models can misbehave or be exploited when deployed. Recent work on adversarial examples, i.e. inputs with minor perturbations that result in substantially different model predictions, is helpful in evaluating the robustness of these models by exposing the adversarial scenarios where they fail. However, these malicious perturbations are often unnatural, not semantically meaningful, and not applicable to complicated domains such as language. In this paper, we propose a framework to generate natural and legible adversarial examples that lie on the data manifold, by searching in semantic space of dense and continuous data representation, utilizing the recent advances in generative adversarial networks. We present generated adversaries to demonstrate the potential of the proposed approach for black-box classifiers for a wide range of applications such as image classification, textual entailment, and machine translation. We include experiments to show that the generated adversaries are natural, legible to humans, and useful in evaluating and analyzing black-box classifiers.


title: Coarse-Graining Open Markov Processes

abstract: Coarse-graining is a standard method of extracting a simple Markov process from a more complicated one by identifying states. Here we extend coarse-graining to open Markov processes. An "open" Markov process is one where probability can flow in or out of certain states called "inputs" and "outputs". One can build up an ordinary Markov process from smaller open pieces in two basic ways: composition, where we identify the outputs of one open Markov process with the inputs of another, and tensoring, where we set two open Markov processes side by side. In previous work, Fong, Pollard and the first author showed that these constructions make open Markov processes into the morphisms of a symmetric monoidal category. Here we go further by constructing a symmetric monoidal double category where the 2-morphisms are ways of coarse-graining open Markov processes. We also extend the already known "black-boxing" functor from the category of open Markov processes to our double category. Black-boxing sends any open Markov process to the linear relation between input and output data that holds in steady states, including nonequilibrium steady states where there is a nonzero flow of probability through the process. To extend black-boxing to a functor between double categories, we need to prove that black-boxing is compatible with coarse-graining.


title: A Sequential Matching Framework for Multi-turn Response Selection in Retrieval-based Chatbots

abstract: We study the problem of response selection for multi-turn conversation in retrieval-based chatbots. The task requires matching a response candidate with a conversation context, whose challenges include how to recognize important parts of the context, and how to model the relationships among utterances in the context. Existing matching methods may lose important information in contexts as we can interpret them with a unified framework in which contexts are transformed to fixed-length vectors without any interaction with responses before matching. The analysis motivates us to propose a new matching framework that can sufficiently carry the important information in contexts to matching and model the relationships among utterances at the same time. The new framework, which we call a sequential matching framework (SMF), lets each utterance in a context interacts with a response candidate at the first step and transforms the pair to a matching vector. The matching vectors are then accumulated following the order of the utterances in the context with a recurrent neural network (RNN) which models the relationships among the utterances. The context-response matching is finally calculated with the hidden states of the RNN. Under SMF, we propose a sequential convolutional network and sequential attention network and conduct experiments on two public data sets to test their performance. Experimental results show that both models can significantly outperform the state-of-the-art matching methods. We also show that the models are interpretable with visualizations that provide us insights on how they capture and leverage the important information in contexts for matching.


title: Tensor Regression Meets Gaussian Processes

abstract: Low-rank tensor regression, a new model class that learns high-order correlation from data, has recently received considerable attention. At the same time, Gaussian processes (GP) are well-studied machine learning models for structure learning. In this paper, we demonstrate interesting connections between the two, especially for multi-way data analysis. We show that low-rank tensor regression is essentially learning a multi-linear kernel in Gaussian processes, and the low-rank assumption translates to the constrained Bayesian inference problem. We prove the oracle inequality and derive the average case learning curve for the equivalent GP model. Our finding implies that low-rank tensor regression, though empirically successful, is highly dependent on the eigenvalues of covariance functions as well as variable correlations.


title: Socialbots supporting human rights

abstract: Socialbots, or non-human/algorithmic social media users, have recently been documented as competing for information dissemination and disruption on online social networks. Here we investigate the influence of socialbots in Mexican Twitter in regards to the "Tanhuato" human rights abuse report. We analyze the applicability of the BotOrNot API to generalize from English to Spanish tweets and propose adaptations for Spanish-speaking bot detection. We then use text and sentiment analysis to compare the differences between bot and human tweets. Our analysis shows that bots actually aided in information proliferation among human users. This suggests that taxonomies classifying bots should include non-adversarial roles as well. Our study contributes to the understanding of different behaviors and intentions of automated accounts observed in empirical online social network data. Since this type of analysis is seldom performed in languages different from English, the proposed techniques we employ here are also useful for other non-English corpora.


title: Bloch-type spaces and extended Ces\`{a}ro operators in the unit ball of a complex Banach space

abstract: Let $\mathbb{B}$ be the unit ball of a complex Banach space $X$. In this paper, we will generalize the Bloch-type spaces and the little Bloch-type spaces to the open unit ball $\mathbb{B}$ by using the radial derivative. Next, we define an extended Ces\`{a}ro operator $T_{\varphi}$ with holomorphic symbol $\varphi$ and characterize those $\varphi$ for which $T_{\varphi}$ is bounded between the Bloch-type spaces and the little Bloch-type spaces. We also characterize those $\varphi$ for which $T_{\varphi}$ is compact between the Bloch-type spaces and the little Bloch-type spaces under some additional assumption on the symbol $\varphi$. When $\mathbb{B}$ is the open unit ball of a finite dimensional complex Banach space $X$, this additional assumption is automatically satisfied.


title: Little-Parks oscillations with half-quantum fluxoid features in Sr2RuO4 micro rings

abstract: In a micro ring of a superconductor with a spin-triplet equal-spin pairing state, a fluxoid, a combined object of magnetic flux and circulating supercurrent, can penetrate as half-integer multiples of the flux quantum. A candidate material to investigate such half-quantum fluxoids is Sr$_\mathsf{2}$RuO$_\mathsf{4}$. We fabricated Sr$_\mathsf{2}$RuO$_\mathsf{4}$ micro rings using single crystals and measured their resistance behavior under magnetic fields controlled with a three-axes vector magnet. Proper Little-Parks oscillations in the magnetovoltage as a function of an axially applied field, associated with fluxoid quantization are clearly observed, for the first time using bulk single crystalline superconductors. We then performed magnetovoltage measurements with additional in-plane magnetic fields. By carefully analyzing both the voltages $V_+$ ($V_-$) measured at positive (negative) current, we find that, above an in-plane threshold field of about 10 mT, the magnetovoltage maxima convert to minima. We interpret this behavior as the peak splitting expected for the half-quantum fluxoid states.


title: Generation and control of noncollinear magnetism by supercurrent

abstract: When superconductivity couples with noncollinear spin textures, rich physics arises; for instance, singlet Cooper pairs can be converted to triplet pairs, and topological superconductors can be realized. For their applications, the controllability of noncollinear magnetism is a crucial issue. Here, we propose that a supercurrent can induce and control noncollinear magnetic orders in a correlated metal on top of a singlet superconductor. We show that the magnetic instability in the correlated metal is enhanced by the proximity effect of supercurrents, which leads to phase transitions from a paramagnetic state to noncollinear magnetic phases with helical or vortexlike spin textures. Furthermore, these magnetic orders can be switched by the direction of the supercurrent. We also discuss the effect of the Rashba spin-orbit coupling and the experimental realization.


title: Grammar Induction for Minimalist Grammars using Variational Bayesian Inference : A Technical Report

abstract: The following technical report presents a formal approach to probabilistic minimalist grammar parameter estimation. We describe a formalization of a minimalist grammar. We then present an algorithm for the application of variational Bayesian inference to this formalization.


title: ChainerMN: Scalable Distributed Deep Learning Framework

abstract: One of the keys for deep learning to have made a breakthrough in various fields was to utilize high computing powers centering around GPUs. Enabling the use of further computing abilities by distributed processing is essential not only to make the deep learning bigger and faster but also to tackle unsolved challenges. We present the design, implementation, and evaluation of ChainerMN, the distributed deep learning framework we have developed. We demonstrate that ChainerMN can scale the learning process of the ResNet-50 model to the ImageNet dataset up to 128 GPUs with the parallel efficiency of 90%.


title: Variations of the cop and robber game on graphs

abstract: We prove new theoretical results about several variations of the cop and robber game on graphs. First, we consider a variation of the cop and robber game which is more symmetric called the cop and killer game. We prove for all $c < 1$ that almost all random graphs are stalemate for the cop and killer game, where each edge occurs with probability $p$ such that $\frac{1}{n^{c}} \le p \le 1-\frac{1}{n^{c}}$. We prove that a graph can be killer-win if and only if it has exactly $k\ge 3$ triangles or none at all. We prove that graphs with multiple cycles longer than triangles permit cop-win and killer-win graphs. For $\left(m,n\right)\neq\left(1,5\right)$ and $n\geq4$, we show that there are cop-win and killer-win graphs with $m$ $C_n$s. In addition, we identify game outcomes on specific graph products. Next, we find a generalized version of Dijkstra's algorithm that can be applied to find the minimal expected capture time and the minimal evasion probability for the cop and gambler game and other variations of graph pursuit. Finally, we consider a randomized version of the killer that is similar to the gambler. We use the generalization of Dijkstra's algorithm to find optimal strategies for pursuing the random killer. We prove that if $G$ is a connected graph with maximum degree $d$, then the cop can win with probability at least $\frac{\sqrt d}{1+\sqrt d}$ after learning the killer's distribution. In addition, we prove that this bound is tight only on the $\left(d+1\right)$-vertex star, where the killer takes the center with probability $\frac1{1+\sqrt d}$ and each of the other vertices with equal probabilities.


title: On polarization of vector light beams: origin of Berry phase

abstract: When generalized from plane waves to general vector beams, the notion of polarization described by the Stokes parameters turns out to be defined in a momentum-associated system that is fixed by the so-called Stratton vector. As the true intrinsic degree of freedom in the language of quantum mechanics, the polarization of light beams in any fixed momentum-associated system is able to characterize their vectorial feature in the laboratory reference system. The Stratton vector is therefore the degree of freedom to distinguish the vectorial feature of light beams that have the same "polarization". Such an observable effect of the Stratton vector helps to understand why plane waves of the same helicity and the same momentum can be different by a Berry phase. This might be the first time to reveal the physical origin of the Berry phase.


title: Spatio-temporal interaction model for crowd video analysis

abstract: We present an unsupervised approach to analyze crowd at various levels of granularity $-$ individual, group and collective. We also propose a motion model to represent the collective motion of the crowd. The model captures the spatio-temporal interaction pattern of the crowd from the trajectory data captured over a time period. Furthermore, we also propose an effective group detection algorithm that utilizes the eigenvectors of the interaction matrix of the model. We also show that the eigenvalues of the interaction matrix characterize various group activities such as being stationary, walking, splitting and approaching. The algorithm is also extended trivially to recognize individual activity. Finally, we discover the overall crowd behavior by classifying a crowd video in one of the eight categories. Since the crowd behavior is determined by its constituent groups, we demonstrate the usefulness of group level features during classification. Extensive experimentation on various datasets demonstrates a superlative performance of our algorithms over the state-of-the-art methods.


title: Necessary Condition for Steerability of Arbitrary Two-Qubit States with Loss

abstract: Einstein-Podolsky-Rosen steering refers to the quantum phenomenon whereby the state of a system held by one party can be "steered" into different states at the will of another, distant, party by performing different local measurements. Although steering has been demonstrated in a number of experiments involving qubits, the question of which two-qubit states are steerable remains an open theoretical problem. Here, we derive a necessary condition for any two-qubit state to be steerable when the steering party suffers from a given probability of qubit loss. Our main result finds application in one-way steering demonstrations that rely upon loss. Specifically, we apply it to a recent experiment on one-way steering with projective measurements and POVMs, reported by Wollmann et. al. [Phys. Rev. Lett., 116, 160403 (2016)].


title: Collimated Ultra-Bright Gamma-Rays from a PW-Laser-Driven Wire Wiggler

abstract: It is shown by three-dimensional QED particle-in-cell simulation that as a laser pulse of 2.5 PW and 20 fs propagates along a sub-wavelength-wide solid wire, directional synchrotron $\gamma-$rays along the wire surface can be efficiently generated. With 8\% energy conversion from the pulse, the $\gamma-$rays contains $10^{12}$ photons between 5 and 500 MeV within 10 fs duration, corresponding to peak brilliance of $10^{27}$ photons ${\rm s^{-1}~ mrad^{-2}~ mm^{-2}}$ per 0.1\% bandwidth. The brilliance and photon energy are respectively 2 and 3 orders of magnitude higher than the highest values of synchrotron radiation facilities. The radiation is attributed to the generation of nC, GeV electron beams well guided along the wire surface and their wiggling motion in strong electrostatic and magnetostatic fields induced at the high-density-wire surface. In particular, these quasistatic fields are so strong that QED effects already play a significant role for the $\gamma-$ray radiation. With the laser power $P_0$ ranging from 0.5 PW to 5 PW available currently, this scheme can robustly produce $\gamma-$rays peaked at $1^\circ$ with few-mrad divergence and the photon energy and number roughly scales with $P_0$ and $P_0^{3/2}$, respectively. Our scheme embraces both the merits of high directionality comparable to those based upon laser wakefield acceleration and high charge comparable to those based upon laser-solid interaction.


title: Kitaev-Heisenberg Hamiltonian for High-Spin $d^7$ Mott Insulators

abstract: In the search for quantum spin liquids, candidate materials for the Kitaev model and its extensions have been intensively explored during the past decade, as the models realize the exact quantum spin liquids in the ground state. Thus far, insulating magnets in the low-spin $d^5$ electron configuration under the strong spin-orbit coupling have been studied for realizing the Kitaev-type bond-dependent anisotropic interactions between the spin-orbital entangled Kramers doublets. To extend the candidates, here we investigate the systems in a high-spin $d^7$ electron configuration, whose ground state is described by the spin-orbital entangled Kramers doublet. By the second- order perturbation in terms of the $t_{2g}$-$t_{2g}$ and $t_{2g}$-$e_g$ hoppings, we show that the effective spin model possesses the anisotropic Kitaev interactions as well as the isotropic Heisenberg ones. While the Kitaev interaction is always ferromagnetic, the Heisenberg interaction can become either ferromagnetic or antiferromagnetic depending on the Coulomb interactions and the crystalline electric fields. We also derive the effective model for the low-spin $d^5$ electron configuration within the same perturbation scheme, in which the Kitaev interaction becomes both ferromagnetic and antiferromagnetic, while the Heisenberg one always ferromagnetic. Referring to the previous study for the Kitaev-Heisenberg model, we find that the quantum spin liquid phase exists in the reasonable parameter region in both $d^7$ and $d^5$ cases, while the former has a richer structure of the phase diagram. We discuss the advantages of the $d^7$ case in comparison with the $d^5$ case. Our results indicate that the high-spin $d^7$ state provides another platform for the Kitaev-type quantum spin liquid.


title: Scalar one-loop Feynman integrals with complex internal masses revisited

abstract: In this paper, we study systematically scalar one-loop two-, three-, and four-point Feynman integrals with complex internal masses. Our analytic results presented in this report are valid for both real and complex internal masses. The calculations are then implemented into a {\tt Mathematica} (version $9$) package and {\tt FORTRAN} program. Our program is cross-checked numerically with {\tt LoopTools} (version $2.14$) in real as well as complex internal masses. We find a perfect agreement between our results and {\tt LoopTools} for all cases. Additionally, this work is applied for evaluating scalar one-loop Feynman integrals developed leading Landau singularities which may appear in real scattering processes at colliders. Last but not least, the method used in this report can also extend to evaluate tensor one-loop integrals. Therefore, this may open a new approach which can solve the inverse Gram determinant problem analytically.


title: Image Patch Matching Using Convolutional Descriptors with Euclidean Distance

abstract: In this work we propose a neural network based image descriptor suitable for image patch matching, which is an important task in many computer vision applications. Our approach is influenced by recent success of deep convolutional neural networks (CNNs) in object detection and classification tasks. We develop a model which maps the raw input patch to a low dimensional feature vector so that the distance between representations is small for similar patches and large otherwise. As a distance metric we utilize L2 norm, i.e. Euclidean distance, which is fast to evaluate and used in most popular hand-crafted descriptors, such as SIFT. According to the results, our approach outperforms state-of-the-art L2-based descriptors and can be considered as a direct replacement of SIFT. In addition, we conducted experiments with batch normalization and histogram equalization as a preprocessing method of the input data. The results confirm that these techniques further improve the performance of the proposed descriptor. Finally, we show promising preliminary results by appending our CNNs with recently proposed spatial transformer networks and provide a visualisation and interpretation of their impact.


title: Orbital tori for non-axisymmetric galaxies

abstract: Our Galaxy's bar makes the Galaxy's potential distinctly non-axisymmetric. All orbits are affected by non-axisymmetry, and significant numbers are qualitatively changed by being trapped at a resonance with the bar. Orbital tori are used to compute these effects. Thick-disc orbits are no less likely to be trapped by corotation or a Lindblad resonance than thin-disc orbits. Perturbation theory is used to create non-axisymmetric orbital tori from standard axisymmetric tori, and both trapped and untrapped orbits are recovered to surprising accuracy. Code is added to the TorusModeller library that makes it as easy to manipulate non-axisymmetric tori as axisymmetric ones. The augmented TorusModeller is used to compute the velocity structure of the solar neighbourhood for bars of different pattern speeds and a simple action-based distribution function. The technique developed here can be applied to any non-axisymmetric potential that is stationary in a rotating from - hence also to classical spiral structure.


title: Testing the distance duality relation using type Ia supernovae and ultracompact radio sources

abstract: We test the possible deviation of the cosmic distance duality relation $D_A(z)(1+z)^2/D_L(z)\equiv 1$ using the standard candles/rulers in a fully model-independent manner. Type-Ia supernovae are used as the standard candles to derive the luminosity distance $D_L(z)$, and ultra-compact radio sources are used as the standard rulers to obtain the angular diameter distance $D_A(z)$. We write the deviation of distance duality relation as $D_A(z)(1+z)^2/D_L(z)=\eta(z)$. Specifically, we use two parameterizations of $\eta(z)$, i.e. $\eta_1(z)=1+\eta_0 z$ and $\eta_2(z)=1+\eta_0 z/(1+z)$. The parameter $\eta_0$ is obtained using the Markov chain Monte Carlo methods by comparing $D_L(z)$ and $D_A(z)$ at the same redshift. The best-fitting results are $\eta_0=-0.06\pm 0.05$ and $-0.18\pm 0.16$ for the first and second parameterizations, respectively. Our results depend on neither the cosmological models, nor the matter contents or the curvature of the universe.


title: Asymptotic behavior in time of solution to the nonlinear Schr"odinger equation with higher order anisotropic dispersion

abstract: We consider the asymptotic behavior in time of solutions to the nonlinear Schr"odinger equation with fourth order anisotropic dispersion (4NLS) which describes the propagation of ultrashort laser pulses in a medium with anomalous time-dispersion in the presence of fourth-order time-dispersion. We prove existence of a solution to (4NLS) which scatters to a solution of the linearized equation of (4NLS).


title: Flat plumbing basket and contact structure

abstract: A flat plumbing basket is a Seifert surface consisting of a disk and bands contained in distinct pages of the disk open book decomposition of the 3-sphere. In this paper, we examine close connections between flat plumbing baskets and the contact structure supported by the open book. As an application we give lower bounds for the flat plumbing basket numbers and determine the flat plumbing basket numbers for various knots and links, including the torus links.


title: Kinematics of the Interstellar Vagabond 1I/'Oumuamua (A/2017 U1)

abstract: The initial Galactic velocity vector for the recently discovered hyperbolic asteroid 1I/'Oumuamua (A/2017 U1) is calculated for before its encounter with our solar system. The latest orbit (JPL-13) shows that 'Oumuamua has eccentricity > 1 at 944\sigma\, significance (1.19936 +- 0.00021), i.e. clearly unbound. Assuming no non-gravitational forces, the object's inbound Galactic velocity was U, V, W = -11.457, -22.395, -7.746 (+-0.009, +-0.009, +-0.011) km/s (U towards Galactic center), with total heliocentric speed 26.32 +- 0.01 km/s. When the velocity is compared to the local stars, 'Oumuamua can be ruled out as co-moving with any of the dozen nearest systems, i.e. it does not appear to be associated with any local exo-Oort clouds (most notably that of the Alpha Centauri triple system). 'Oumuamua's velocity is within 5 km/s of the median Galactic velocity of the stars in the solar neighborhood (<25 pc), and within 2 km/s of the mean velocity of the local M dwarfs. Its velocity appears to be statistically "too" typical for a body whose velocity was drawn from the Galactic velocity distribution of the local stars (i.e. less than 1 in 500 field stars in the solar neighborhood would have a velocity so close to the median UVW velocity). In the Local Standard of Rest frame (circular Galactic motion), 'Oumuamua is remarkable for showing both negligible radial (U) and vertical (W) motion, while having a slightly sub-Keplerian circular velocity (V; by ~11 km/s). These calculations strengthen the interpretation that A/2017 U1 has a distant extrasolar origin, but not among the very nearest stars. Any formation mechanism for this interstellar asteroid should account for the coincidence of 'Oumuamua's velocity being so close to the LSR.


title: Coulomb drag transistor via graphene/MoS2 heterostructures

abstract: Two-dimensional (2D) heterointerfaces often provide extraordinary carrier transport as exemplified by superconductivity or excitonic superfluidity. Recently, double-layer graphene separated by few-layered boron nitride demonstrated the Coulomb drag phenomenon: carriers in the active layer drag the carriers in the passive layer. Here, we propose a new switching device operating via Coulomb drag interaction at a graphene/MoS2 (GM) heterointerface. The ideal van der Waals distance allows strong coupling of the interlayer electron-hole pairs, whose recombination is prevented by the Schottky barrier formed due to charge transfer at the heterointerface. This device exhibits a high carrier mobility (up to ~3,700 cm^2V^-1s^-1) even at room temperature, while maintaining a high on/off current ratio (~10^8), outperforming those of individual layers. In the electron-electron drag regime, graphene-like Shubnikov-de Haas oscillations are observed at low temperatures. Our Coulomb drag transistor could provide a shortcut for the practical application of quantum-mechanical 2D heterostructures at room temperature.


title: Continuity of Gevrey-H\"ormander pseudo-differential operators on modulation spaces

abstract: Let $s\ge 1$, $\omega ,\omega_0\in \mathscr P_{E,s}^0$, $a\in \Gamma _{s}^{(\omega_0)}$, and let $\mathscr B$ be a suitable invariant quasi-Banach function space, Then we prove that the pseudo-differential operator $\operatorname{Op} (a)$ is continuous from $M(\omega_0\omega ,\mathscr B )$ to $M(\omega ,\mathscr B )$.


title: The Values of the Riemann Zeta-Function on Discrete Sets

abstract: We study the values taken by the Riemann zeta-function $\zeta$ on discrete sets. We show that infinite vertical arithmetic progressions are uniquely determined by the values of $\zeta$ taken on this set. Moreover, we prove a joint discrete universality theorem for $\zeta$ with respect to certain permutations of the set of positive integers. Finally, we study a generalization of the classical denseness theorems for $\zeta$.