title: On smooth square-free numbers in arithmetic progressions

abstract: A. Booker and C. Pomerance (2017) have shown that any residue class modulo a prime $p\ge 11$ can be represented by a positive $p$-smooth square-free integer $s = p^{O(\log p)}$ with all prime factors up to $p$ and conjectured that in fact one can find such $s$ with $s = p^{O(1)}$. Using bounds on double Kloosterman sums due to M. Z. Garaev (2010) we prove this conjecture in a stronger form $s \le p^{3/2 + o(1)}$ and also consider more general versions of this question replacing $p$-smoothness of $s$ by the stronger condition of $p^{\alpha}$-smoothness. Using bounds on multiplicative character sums and a sieve method, we also show that we can represent all residue classes by a positive square-free integer $s\le p^{2+o(1)}$ which is $p^{1/(4e^{ /2})+o(1)}$-smooth. Additionally, we obtain stronger results for almost all primes $p$.


title: Search for excited states in $^{25}$O

abstract: Theoretical calculations suggest the presence of low-lying excited states in $^{25}$O. Previous experimental searches by means of proton knockout on $^{26}$F produced no evidence for such excitations. We search for excited states in $^{25}$O using the ${ {}^{24}\text{O} (d,p) {}^{25}\text{O} }$ reaction. The theoretical analysis of excited states in unbound $^{25,27}$O is based on the configuration interaction approach that accounts for couplings to the scattering continuum. We use invariant-mass spectroscopy to measure neutron-unbound states in $^{25}$O. For the theoretical approach, we use the complex-energy Gamow Shell Model and Density Matrix Renormalization Group method with a finite-range two-body interaction optimized to the bound states and resonances of $^{23-26}$O, assuming a core of $^{22}$O. We predict energies, decay widths, and asymptotic normalization coefficients. Our calculations in a large $spdf$ space predict several low-lying excited states in $^{25}$O of positive and negative parity, and we obtain an experimental limit on the relative cross section of a possible ${ {J}^{\pi} = {1/2}^{+} }$ state with respect to the ground-state of $^{25}$O at $\sigma_{1/2+}/\sigma_{g.s.} = 0.25_{-0.25}^{+1.0}$. We also discuss how the observation of negative parity states in $^{25}$O could guide the search for the low-lying negative parity states in $^{27}$O. Previous experiments based on the proton knockout of $^{26}$F suffered from the low cross sections for the population of excited states in $^{25}$O because of low spectroscopic factors. In this respect, neutron transfer reactions carry more promise.


title: A search for transit timing variations and orbital decay in WASP-46b

abstract: We present 12 new transit observations of the exoplanet WASP-46b obtained with the 1.54-m telescope at Estaci\'on Astrof\'isica de Bosque Alegre (EABA, Argentina) and the 0.40-m Horacio Ghielmetti and 2.15-m Jorge Sahade telescopes at Complejo Astron\'omico El Leoncito (CASLEO, Argentina). We analyse them together with 37 light curves from the literature to re-determine the physical parameters and search for additional planets via transit timing variations (TTVs). We consider the 31 transits with uncertainties in their mid-transit times ($e_\mathrm{T_{0}}$) $<$ 1 minute, to perform the first homogeneous study of TTVs for the system, finding a dispersion of $\sigma= 1.66$ minutes over a 6 year baseline. Since no periodic variations are found, our interpretation for this relatively high value of $\sigma$ is that the stellar activity could be affecting the measured mid-transit times. This value of dispersion allows us to rule out the presence of additional bodies with masses larger than 2.3, 4.6, 7, and 9.3 $M_{\mathrm{\bigoplus}}$ at the first-order mean-motion resonances 2:1, 3:2, 4:3, and 5:4 with the transiting planet, respectively. Despite the 6 year baseline and a typical light curve precision of $2 \times 10^{-3}$, we find that we cannot significantly demonstrate a slow decrease of the orbital period of WASP-46b. We place a lower limit of Q$_{\star}$ > $7 \times 10^{3}$ on the tidal quality factor and determine that an additional 6 year baseline is required to rule out Q$_{\star}$ < $10^{5}$.


title: Perfect Fluids

abstract: We develop a new theory of perfect fluids with translation and rotation symmetry, which is also applicable in the absence of any type of boost symmetry. It involves introducing a new fluid variable, the kinetic mass density, which is needed to define the most general energy-momentum tensor for perfect fluids. Our theory leads to corrections to the Euler equations for perfect fluids that might be observable in hydrodynamic fluid experiments. We also derive new expressions for the speed of sound in perfect fluids. Our theory reduces to the known perfect fluid models when boost symmetry is present. It can also be adapted to (non-relativistic) scale invariant fluids with critical exponent $z$. We show that perfect fluids cannot have Schr\"odinger symmetry unless $z=2$. For generic values of $z$ there can be fluids with Lifshitz symmetry, and as a concrete example, we work out in detail the thermodynamics and fluid description of an ideal gas of Lifshitz particles and compute the speed of sound for the classical and quantum Lifshitz gasses.


title: Computing twisted KLV polynomials

abstract: In order to compute Hermitian forms on representations of real reductive groups, in the unequal rank case, it is necessary to compute twisted Kazhdan-Lusztig-Vogan polynomials. These were defined by Lusztig and Vogan (Quasisplit Hecke algebras and Symmetric Spaces, Duke, 2014) and discussed further by Adams and Vogan (Parameters for twisted representations, 2015). These notes contain the details necessary to go from what is in those papers to an explicit algorithm. This algorithm has been implemented in the Atlas of Lie Groups and Representations software.


title: A resource theory of quantum memories and their faithful verification with minimal assumptions

abstract: We provide a complete set of game-theoretic conditions equivalent to the existence of a transformation from one quantum channel into another one, by means of classically correlated pre/post processing maps only. Such conditions naturally induce tests to certify that a quantum memory is capable of storing quantum information, as opposed to memories that can be simulated by measurement and state preparation (corresponding to entanglement-breaking channels). These results are formulated as a resource theory of genuine quantum memories (correlated in time), mirroring the resource theory of entanglement in quantum states (correlated spatially). As the set of conditions is complete, the corresponding tests are faithful, in the sense that any non entanglement-breaking channel can be certified. Moreover, they only require the assumption of trusted inputs, known to be unavoidable for quantum channel verification. As such, the tests we propose are intrinsically different from the usual process tomography, for which the probes of both the input and the output of the channel must be trusted. An explicit construction is provided and shown to be experimentally realizable, even in the presence of arbitrarily strong losses in the memory or detectors.


title: The topological complexity of the free product

abstract: We prove the formula $TC(G\ast H)=\max\{TC(G), TC(H), cd(G\times H)\}$ for the topological complexity of the free product of discrete groups with cohomological dimension >2.


title: A unified model of supernova driven by magnetic monopoles

abstract: In this paper, we first discuss a series of important but puzzling physical mechanisms concerning the energy source, various kinds of core collapsed supernovae explosion mechanisms during central gravitational collapse in astrophysics. We also discuss the puzzle of possible association of $\gamma$-ray burst with gravitational wave perturbation, the heat source for the molten interior of the core of the earth and finally the puzzling problem of the cooling of white dwarfs. We then make use of the estimations for the space flux of magnetic monopoles (hereafter MMs) and nucleon decay induced by MMs (called the Rubakov-Callen(RC) effect) to obtain the luminosity due to the RC effect. In terms of the formula for this RC luminosity, we present a unified treatment for the heat source of the Earth's core, the energy source for the white dwarf interior, various kinds of core collapsed supernovae (Type II Supernova (SNII), Type Ib Supernova (SNIb), Type Ic Supernova (SNIc), Super luminous supernova (SLSN)), and the production mechanism for $\gamma$-ray burst. This unified model can also be used to reasonably explain the possible association of the short $\gamma$-ray burst detected by the Fermi $\gamma$-ray Burst Monitoring Satellite (GBM) with the LIGO gravitational wave event GW150914 in September 2015.


title: Chiral Topological Phases from Artificial Neural Networks

abstract: Motivated by recent progress in applying techniques from the field of artificial neural networks (ANNs) to quantum many-body physics, we investigate as to what extent the flexibility of ANNs can be used to efficiently study systems that host chiral topological phases such as fractional quantum Hall (FQH) phases. With benchmark examples, we demonstrate that training ANNs of restricted Boltzmann machine type in the framework of variational Monte Carlo can numerically solve FQH problems to good approximation. Furthermore, we show by explicit construction how n-body correlations can be kept at an exact level with ANN wave-functions exhibiting polynomial scaling with power n in system size. Using this construction, we analytically represent the paradigmatic Laughlin wave-function as an ANN state.


title: Chandra imaging of the $\sim$kpc extended outflow in 1H 0419-577

abstract: The Seyfert 1 galaxy 1H 0419-577 hosts a $\sim$kpc extended outflow that is evident in the [\ion{O}{iii}] image and that is also detected as a warm absorber in the UV/X-ray spectrum. Here, we analyze a $\sim$30 ks Chandra-ACIS X-ray image, with the aim of resolving the diffuse extranuclear X-ray emission and of investigating its relationship with the galactic outflow. Thanks to its sub-arcsecond spatial resolution, Chandra resolves the circumnuclear X-ray emission, which extends up to a projected distance of at least $\sim$16 kpc from the center. The morphology of the diffuse X-ray emission is spherically symmetrical. We could not recover a morphological resemblance between the soft X-ray emission and the ionization bicone that is traced by the [\ion{O}{iii}] outflow. we argue that the photoionized gas nebula must be distributed mostly along the polar directions, outside our line of sight. In this geometry, the X-ray/UV warm absorber must trace a different gas component, physically disconnected from the emitting gas, and located closer to the equatorial plane.


title: New results on collectivity with ALICE

abstract: An overview of recent ALICE results aimed to understand collective phenomena in Pb-Pb collisions at the LHC is presented. These include the centrality dependence of the transverse momentum ($p_{\rm T}$) distributions of charged pions, kaons, and protons as well as results on anisotropic flow using data from Pb-Pb collisions at $\sqrt{s_{{\mathrm{NN}}}} = 5.02$ TeV. In addition, comparisons between data and model predictions are presented.


title: Local and nonlocal order parameters in the Kitaev chain

abstract: We have calculated order parameters for the phases of the Kitaev chain with interaction and dimerization at a special symmetric point applying the Jordan-Wigner and other duality transformations. We use string order parameters (SOPs) defined via the correlation functions of the Majorana string operators. The SOPs are mapped onto the local order parameters of some dual Hamiltonians and easily calculated. We have shown that the phase diagram of the interacting dimerized chain comprises the phases with the conventional local order as well as the phases with nonlocal SOPs. From the results for the critical indices we infer the 2D Ising universality class of criticality at the particular symmetry point where the model is exactly solvable.


title: Generalised ballooning theory of two dimensional tokamak modes

abstract: In this work, using solutions from a local gyrokinetic flux-tube code combined with higher order ballooning theory, a new analytical approach is developed to reconstruct the global linear mode structure with associated global mode frequency. In addition to the isolated mode (IM), which usually peaks on the outboard mid-plane, the higher order ballooning theory has also captured other types of less unstable global modes: (a) the weakly asymmetric ballooning theory (WABT) predicts a mixed mode (MM) that undergoes a small poloidal shift away from the outboard mid-plane, (b) a relatively more stable general mode (GM) balloons on the top (or bottom) of the tokamak plasma. In this paper, an analytic approach is developed to combine these disconnected analytical limits into a single generalised ballooning theory (GBT). This is used to investigate how an IM behaves under the effect of sheared toroidal flow. For small values of flow an IM initially converts into a MM where the results of WABT are recaptured, and eventually, as the flow increases, the mode asymptotically becomes a GM on the top (or bottom) of the plasma. This may be an ingredient in models for understanding why in some experimental scenarios, instead of large edge localised modes (ELMs), small ELMs are observed. Finally, our theory can have other important consequences, especially for calculations involving Reynolds stress driven intrinsic rotation through the radial asymmetry in the global mode structures. Understanding the intrinsic rotation is significant because external torque in a plasma the size of ITER is expected to be relatively low.


title: "Relative-Continuity" for Non-Lipschitz Non-Smooth Convex Optimization using Stochastic (or Deterministic) Mirror Descent

abstract: The usual approach to developing and analyzing first-order methods for non-smooth (stochastic or deterministic) convex optimization assumes that the objective function is uniformly Lipschitz continuous with parameter $M_f$. However, in many settings the non-differentiable convex function $f(\cdot)$ is not uniformly Lipschitz continuous -- for example (i) the classical support vector machine (SVM) problem, (ii) the problem of minimizing the maximum of convex quadratic functions, and even (iii) the univariate setting with $f(x) := \max\{0, x\} + x^2$. Herein we develop a notion of "relative continuity" that is determined relative to a user-specified "reference function" $h(\cdot)$ (that should be computationally tractable for algorithms), and we show that many non-differentiable convex functions are relatively continuous with respect to a correspondingly fairly-simple reference function $h(\cdot)$. We also similarly develop a notion of "relative stochastic continuity" for the stochastic setting. We analysis two standard algorithms -- the (deterministic) mirror descent algorithm and the stochastic mirror descent algorithm -- for solving optimization problems in these two new settings, and we develop for the first time computational guarantees for instances where the objective function is not uniformly Lipschitz continuous. This paper is a companion paper for non-differentiable convex optimization to the recent paper by Lu, Freund, and Nesterov, which developed similar sorts of results for differentiable convex optimization.


title: The second inner variation of energy and the Morse index of limit interfaces

abstract: In this article we study the second variation of the energy functional associated to the Allen-Cahn equation on closed manifolds. Extending well known analogies between the gradient theory of phase transitions and the theory of minimal hypersurfaces, we prove the upper semicontinuity of the eigenvalues of the stability operator and consequently obtain upper bounds for the Morse index of limit interfaces which arise from solutions with bounded energy and index without assuming any multiplicity or orientability condition on these hypersurfaces. This extends some recent results of N. Le and F. Hiesmayr.


title: Asymptotic properties of the set of systoles of arithmetic Riemann surfaces

abstract: The purpose this article is to try to understand the mysterious coincidence between the asymptotic behavior of the volumes of the Moduli Space of closed hyperbolic surfaces of genus $g$ with respect to the Weil-Petersson metric and the asymptotic behavior of the number of arithmetic closed hyperbolic surfaces of genus $g$. If the set of arithmetic surfaces is well distributed then its image for any interesting function should be well distributed too. We investigate the distribution of the function systole. We give several results indicating that the systoles of arithmetic surfaces can not be concentrated, consequently the same holds for the set of arithmetic surfaces. The proofs are based in different techniques: combinatorics (obtaining regular graphs with any girth from results of B. Bollobas and constructions with cages and Ramanujan graphs), group theory (constructing finite index subgroups of surface groups from finite index subgroups of free groups using results of G. Baumslag) and geometric group theory (linking the geometry of graphs with the geometry of coverings of a surface).


title: Cox regression analysis with missing covariates via multiple imputation

abstract: We consider the situation of estimating Cox regression in which some covariates are subject to missing, and there exists additional information (including observed event time, censoring indicator and fully observed covariates) which may be predictive of the missing covariates. We propose to use two working regression models: one for predicting the missing covariates and the other for predicting the missing probabilities. For each missing covariate observation, these two working models are used to define a nearest neighbor imputing set. This set is then used to nonparametrically impute covariate values for the missing observation. Upon the completion of imputation, Cox regression is performed on the multiply imputed datasets to estimate the regression coefficients. In a simulation study, we compare the nonparametric multiple imputation approach with the augmented inverse probability weighted (AIPW) method, which directly incorporates the two working models into estimation of Cox regression. We show that all approaches can reduce bias due to non-ignorable missing mechanism. The proposed nonparametric imputation method is robust to mis-specification of either one of the two working models and robust to mis-specification of the link function of the two working models. In contrast, the AIPW method is not robust to mis-specification of the link function of the two working models and is sensitive to the selection probability. We apply the approaches to a breast cancer dataset from Surveillance, Epidemiology and End Results (SEER) Program .


title: The inverse hull of 0-left cancellative semigroups

abstract: Given a semigroup S with zero, which is left-cancellative in the sense that st=sr \neq 0 implies that t=r, we construct an inverse semigroup called the inverse hull of S, denoted H(S). When S admits least common multiples, in a precise sense defined below, we study the idempotent semilattice of H(S), with a focus on its spectrum. When S arises as the language semigroup for a subsift X on a finite alphabet, we discuss the relationship between H(S) and several C*-algebras associated to X appearing in the literature.


title: Harnessing bistability for directional propulsion of untethered, soft robots

abstract: In most macro-scale robotics systems , propulsion and controls are enabled through a physical tether or complex on-board electronics and batteries. A tether simplifies the design process but limits the range of motion of the robot, while on-board controls and power supplies are heavy and complicate the design process. Here we present a simple design principle for an untethered, entirely soft, swimming robot with the ability to achieve preprogrammed, directional propulsion without a battery or on-board electronics. Locomotion is achieved by employing actuators that harness the large displacements of bistable elements, triggered by surrounding temperature changes. Powered by shape memory polymer (SMP) muscles, the bistable elements in turn actuates the robot's fins. Our robots are fabricated entirely using a commercially available 3D printer with no post-processing. As a proof-of-concept, we demonstrate the ability to program a vessel, which can autonomously deliver a cargo and navigate back to the deployment point.


title: Training deep neural networks for the inverse design of nanophotonic structures

abstract: Data inconsistency leads to a slow training process when deep neural networks are used for the inverse design of photonic devices, an issue that arises from the fundamental property of non-uniqueness in all inverse scattering problems. Here we show that by combining forward modeling and inverse design in a tandem architecture, one can overcome this fundamental issue, allowing deep neural networks to be effectively trained by data sets that contain non-unique electromagnetic scattering instances. This paves the way for using deep neural networks to design complex photonic structures that requires large training sets.


title: Hyperparameter Importance Across Datasets

abstract: With the advent of automated machine learning, automated hyperparameter optimization methods are by now routinely used in data mining. However, this progress is not yet matched by equal progress on automatic analyses that yield information beyond performance-optimizing hyperparameter settings. In this work, we aim to answer the following two questions: Given an algorithm, what are generally its most important hyperparameters, and what are typically good values for these? We present methodology and a framework to answer these questions based on meta-learning across many datasets. We apply this methodology using the experimental meta-data available on OpenML to determine the most important hyperparameters of support vector machines, random forests and Adaboost, and to infer priors for all their hyperparameters. The results, obtained fully automatically, provide a quantitative basis to focus efforts in both manual algorithm design and in automated hyperparameter optimization. The conducted experiments confirm that the hyperparameters selected by the proposed method are indeed the most important ones and that the obtained priors also lead to statistically significant improvements in hyperparameter optimization.


title: RAWSim-O: A Simulation Framework for Robotic Mobile Fulfillment Systems

abstract: This paper deals with a new type of warehousing system, Robotic Mobile Fulfillment Systems (RMFS). In such systems, robots are sent to carry storage units, so-called "pods", from the inventory and bring them to human operators working at stations. At the stations, the items are picked according to customers' orders. There exist new decision problems in such systems, for example, the reallocation of pods after their visits at work stations or the selection of pods to fulfill orders. In order to analyze decision strategies for these decision problems and relations between them, we develop a simulation framework called "RAWSim-O" in this paper. Moreover, we show a real-world application of our simulation framework by integrating simple robot prototypes based on vacuum cleaning robots.


title: Quantified Gamow Shell Model interaction for $psd$-shell nuclei

abstract: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interactions approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow Shell Model (GSM) provides such a framework as it is capable of describing resonant and non-resonant many-body states on equal footing. To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the $psdf$ shell model space. The resulting interaction is expected to describe nuclei with $5 \leqslant A \leqslant 12$ at the $p-sd$-shell interface. The optimized one-body potential reproduces nucleon-$^4$He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the Helium, Lithium, and Beryllium isotopes up to $A=9$. A very good agreement with experiment was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the $psd$ region of the nuclear chart.


title: Entropy operates in non-linear semifields

abstract: We set out to demonstrate that the R\'enyi entropies with parameter $\alpha$ are better thought of as operating in a type of non-linear semiring called a positive semifield. We show how the R\'enyi's postulates lead to Pap's g-calculus where the functions carrying out the domain transformation are Renyi's information function and its inverse. In its turn, Pap's g-calculus under R\'enyi's information function transforms the set of positive reals into a family of semirings where "standard" product has been transformed into sum and "standard" sum into a power-deformed sum. Consequently, the transformed product has an inverse whence the structure is actually that of a positive semifield. Instances of this construction lead into idempotent analysis and tropical algebra as well as to less exotic structures. Furthermore, shifting the definition of the $\alpha$ parameter shows in full the intimate relation of the R\'enyi entropies to the weighted generalized power means. We conjecture that this is one of the reasons why tropical algebra procedures, like the Viterbi algorithm of dynamic programming, morphological processing, or neural networks are so successful in computational intelligence applications. But also, why there seem to exist so many procedures to deal with "information" at large.


title: Analysis of Dark Matter Axion Clumps with Spherical Symmetry

abstract: Recently there has been much interest in the spatial distribution of light scalar dark matter, especially axions, throughout the universe. When the local gravitational interactions between the scalar modes are sufficiently rapid, it can cause the field to re-organize into a BEC of gravitationally bound clumps. While these clumps are stable when only gravitation is included, the picture is complicated by the presence of the axion's attractive self-interactions, which can potentially cause the clumps to collapse. Here we perform a detailed stability analysis to determine under what conditions the clumps are stable. In this paper we focus on spherical configurations, leaving aspherical configurations for future work. We identify branches of clump solutions of the axion-gravity-self-interacting system and study their stability properties. We find that clumps that are (spatially) large are stable, while clumps that are (spatially) small are unstable and may collapse. Furthermore, there is a maximum number of particles that can be in a clump. We map out the full space of solutions, which includes quasi-stable axitons, and clarify how a recent claim in the literature of a new ultra-dense branch of stable solutions rests on an invalid use of the non-relativistic approximation. We also consider repulsive self-interactions that may arise from a generic scalar dark matter candidate, finding a single stable branch that extends to arbitrary particle number.


title: Bosonic topological phases of matter: bulk-boundary correspondence, SPT invariants and gauging

abstract: We analyze $2+1d$ and $3+1d$ Bosonic Symmetry Protected Topological (SPT) phases of matter protected by onsite symmetry group $G$ by using dual bulk and boundary approaches. In the bulk we study an effective field theory which upon coupling to a background flat $G$ gauge field furnishes a purely topological response theory. The response action evaluated on certain manifolds, with appropriate choice of background gauge field, defines a set of SPT topological invariants. Further, SPTs can be gauged by summing over all isomorphism classes of flat $G$ gauge fields to obtain Dijkgraaf-Witten topological $G$ gauge theories. These topological gauge theories can be ungauged by first introducing and then proliferating defects that spoils the gauge symmetry. This mechanism is related to anyon condensation in $2+1d$ and condensing bosonic gauge charges in $3+1d$. In the dual boundary approach, we study $1+1d$ and $2+1d$ quantum field theories that have $G$ 't-Hooft anomalies that can be precisely cancelled by (the response theory of) the corresponding bulk SPT. We show how to construct/compute topological invariants for the bulk SPTs directly from the boundary theories. Further we sum over boundary partition functions with different background gauge fields to construct $G$-characters that generate topological data for the bulk topological gauge theory. Finally, we study a $2+1d$ quantum field theory with a mixed $\mathbb{Z}_2^{T/R} \times U(1)$ anomaly where $\mathbb{Z}_2^{T/R}$ is time-reversal/reflection symmetry, and the $U(1)$ could be a 0-form or 1-form symmetry depending on the choice of time reversal/reflection action. We briefly discuss the bulk effective action and topological response for a theory in $3+1d$ that cancels this anomaly. This signals the existence of SPTs in $3+1d$ protected by 0,1-form $U(1)\times \mathbb{Z}_{2}^{T,R}$.


title: Planning, Fast and Slow: A Framework for Adaptive Real-Time Safe Trajectory Planning

abstract: Motion planning is an extremely well-studied problem in the robotics community, yet existing work largely falls into one of two categories: computationally efficient but with few if any safety guarantees, or able to give stronger guarantees but at high computational cost. This work builds on a recent development called FaSTrack in which a slow offline computation provides a modular safety guarantee for a faster online planner. We introduce the notion of "meta-planning" in which a refined offline computation enables safe switching between different online planners. This provides autonomous systems with the ability to adapt motion plans to a priori unknown environments in real-time as sensor measurements detect new obstacles, and the flexibility to maneuver differently in the presence of obstacles than they would in free space, all while maintaining a strict safety guarantee. We demonstrate the meta-planning algorithm both in simulation and in hardware using a small Crazyflie 2.0 quadrotor.


title: Existence of Positive Steady States for Weakly Reversible Mass-Action Systems

abstract: We prove the following. For each weakly reversible mass-action system, there exists a positive steady state in each positive stoichiometric class.


title: A poset $\Phi_n$ whose maximal chains are in bijection with the $n \times n$ alternating sign matrices

abstract: For an integer $n\geq 1$, we display a poset $\Phi_n$ whose maximal chains are in bijection with the $n\times n$ alternating sign matrices. The Hasse diagram $\widehat \Phi_n$ is obtained from the $n$-cube by adding some edges. We show that the dihedral group $D_{2n}$ acts on $\widehat \Phi_n$ as a group of automorphisms.


title: STDP Based Pruning of Connections and Weight Quantization in Spiking Neural Networks for Energy Efficient Recognition

abstract: Spiking Neural Networks (SNNs) with a large number of weights and varied weight distribution can be difficult to implement in emerging in-memory computing hardware due to the limitations on crossbar size (implementing dot product), the constrained number of conductance levels in non-CMOS devices and the power budget. We present a sparse SNN topology where non-critical connections are pruned to reduce the network size and the remaining critical synapses are weight quantized to accommodate for limited conductance levels. Pruning is based on the power law weight-dependent Spike Timing Dependent Plasticity (STDP) model; synapses between pre- and post-neuron with high spike correlation are retained, whereas synapses with low correlation or uncorrelated spiking activity are pruned. The weights of the retained connections are quantized to the available number of conductance levels. The process of pruning non-critical connections and quantizing the weights of critical synapses is performed at regular intervals during training. We evaluated our sparse and quantized network on MNIST dataset and on a subset of images from Caltech-101 dataset. The compressed topology achieved a classification accuracy of 90.1% (91.6%) on the MNIST (Caltech-101) dataset with 3.1x (2.2x) and 4x (2.6x) improvement in energy and area, respectively. The compressed topology is energy and area efficient while maintaining the same classification accuracy of a 2-layer fully connected SNN topology.


title: On the Runtime-Efficacy Trade-off of Anomaly Detection Techniques for Real-Time Streaming Data

abstract: Ever growing volume and velocity of data coupled with decreasing attention span of end users underscore the critical need for real-time analytics. In this regard, anomaly detection plays a key role as an application as well as a means to verify data fidelity. Although the subject of anomaly detection has been researched for over 100 years in a multitude of disciplines such as, but not limited to, astronomy, statistics, manufacturing, econometrics, marketing, most of the existing techniques cannot be used as is on real-time data streams. Further, the lack of characterization of performance -- both with respect to real-timeliness and accuracy -- on production data sets makes model selection very challenging. To this end, we present an in-depth analysis, geared towards real-time streaming data, of anomaly detection techniques. Given the requirements with respect to real-timeliness and accuracy, the analysis presented in this paper should serve as a guide for selection of the "best" anomaly detection technique. To the best of our knowledge, this is the first characterization of anomaly detection techniques proposed in very diverse set of fields, using production data sets corresponding to a wide set of application domains.


title: Search for Gamma-ray Emission from Galactic Novae with the Fermi-LAT

abstract: Context. A number of novae have been found to emit high-energy gamma rays (>100 MeV). However, the origin of this emission is not yet understood. We report on the search for gamma-ray emission from 75 optically-detected Galactic novae in the first 7.4 years of operation of the Fermi Large Area Telescope using the Pass 8 data set. Aims. We compile an optical nova catalog including light curves from various resources and estimate the optical peak time and optical peak magnitude in order to search for gamma-ray emission to test if all novae are gamma-ray emitters. Methods. We repeat the analysis of the six novae previously identified as gamma-ray sources and develop a unified analysis strategy which we then apply to all novae in our catalog. We search for emission in a 15-day time window in two-day steps ranging from 20 days before to 20 days after the optical peak time. We perform a population study with Monte Carlo simulations to set constraints on the properties of the gamma-ray emission of novae. Results. Two new novae candidates have been found at 2sigma global significance. Although these two novae candidates were not detected at a significant level individually, taking them together with the other non-detected novae, we found a sub-threshold nova population with a cumulative 3sigma significance. We report the measured gamma-ray flux for detected sources and flux upper limits for novae without significant detection. Our results can be reproduced by several gamma-ray emissivity models (e.g. a power-law distribution with a slope of 2), while a constant emissivity model (i.e. assuming novae are standard candles) can be rejected.


title: Gravitational Decoupling and Picard-Lefschetz

abstract: In this work, we consider tunneling between non-metastable states in gravitational theories. Such processes arise in various contexts, e.g., in inflationary scenarios where the inflaton potential involves multiple fields or multiple branches. They are also relevant for bubble wall nucleation in some cosmological settings. However, we show that the transition amplitudes computed using the Euclidean method generally do not approach the corresponding field theory limit as $M_{p}\rightarrow \infty$. This implies that in the Euclidean framework, there is no systematic expansion in powers of $G_{N}$ for such processes. Such considerations also carry over directly to no-boundary scenarios involving Hawking-Turok instantons. In this note, we illustrate this failure of decoupling in the Euclidean approach with a simple model of axion monodromy and then argue that the situation can be remedied with a Lorentzian prescription such as the Picard-Lefschetz theory. As a proof of concept, we illustrate with a simple model how tunneling transition amplitudes can be calculated using the Picard-Lefschetz approach.


title: Characterizing the Motion of Solar Magnetic Bright Points at High Resolution

abstract: Magnetic bright points in the solar photosphere, visible in both continuum and G-band images, indicate footpoints of kilogauss magnetic flux tubes extending to the corona. The power spectrum of bright-point motion is thus also the power spectrum of Alfven wave excitation, transporting energy up flux tubes into the corona. This spectrum is a key input in coronal and heliospheric models. We produce a power spectrum of bright-point motion using radiative magnetohydrodynamic simulations, exploiting spatial resolution higher than can be obtained in present-day observations, while using automated tracking to produce large data quantities. We find slightly higher amounts of power at all frequencies compared to observation-based spectra, while confirming the spectrum shape of recent observations. This also provides a prediction for observations of bright points with DKIST, which will achieve similar resolution and high sensitivity. We also find a granule size distribution in support of an observed two-population distribution, and we present results from tracking passive tracers which show a similar power spectrum to that of bright points. Finally, we introduce a simplified, laminar model of granulation, with which we explore the roles of turbulence and of the properties of the granulation pattern in determining bright-point motion.


title: The p-centre of Yangians and shifted Yangians

abstract: We study the Yangian $Y_n$ associated to the general linear Lie algebra $\mathfrak{gl}_n$ over a field of positive characteristic, as well as its shifted analog $Y_n(\sigma)$. Our main result gives a description of the centre of $Y_n(\sigma)$: it is a polynomial algebra generated by its Harish-Chandra centre (which lifts the centre in characteristic zero) together with a large $p$-centre. Moreover, $Y_n(\sigma)$ is free as a module over its center. In future work, it will be seen that every reduced enveloping algebra $U_\chi(\mathfrak{gl}_n)$ is Morita equivalent to a quotient of an appropriate choice of shifted Yangian, and so our results will have applications in classical representation theory.


title: Structured Robust Submodular Maximization: Offline and Online Algorithms

abstract: Constrained submodular function maximization has been used in subset selection problems such as selection of most informative sensor locations. While these models have been quite popular, the solutions Constrained submodular function maximization has been used in subset selection problems such as selection of most informative sensor locations. While these models have been quite popular, the solutions obtained via this approach are unstable to perturbations in data defining the submodular functions. Robust submodular maximization has been proposed as a richer model that aims to overcome this discrepancy as well as increase the modeling scope of submodular optimization. In this work, we consider robust submodular maximization with structured combinatorial constraints and give efficient algorithms with provable guarantees. Our approach is applicable to constraints defined by single or multiple matroids, knapsack as well as distributionally robust criteria. We consider both the offline setting where the data defining the problem is known in advance as well as the online setting where the input data is revealed over time. For the offline setting, we give a general (nearly) optimal bi-criteria approximation algorithm that relies on new extensions of classical algorithms for submodular maximization. For the online version of the problem, we give an algorithm that returns a bi-criteria solution with sub-linear regret.


title: Universality of the underlying event in pp collisions

abstract: In this paper we study ATLAS results on underlying event in pp collisions at $\sqrt{s}$=0.9, 7 and 13 TeV. We show that the center-of-mass energy dependences of the charged-particle production sensitive to the underlying event ("transverse" region) and to the hardest partonic interaction ("towards" and "away" regions) in pp collisions can be both understood in terms of the change of the inclusive average multiplicity. Within uncertainties, the corresponding particle production as a function of the leading charged particle shows no significant $\sqrt{s}$-dependence for the three regions once they are scaled according to the relative change in multiplicity. The scaling properties reported here are well reproduced by PYTHIA 8.212 tune Monash 2013 and suggest an universality of the underlying event in hadronic interactions at high $\sqrt{s}$. Based on the simulations, we observed that the same scaling properties are also present in the average number of multi-partonic interactions as a function of the leading charged particle. Moreover, the multiplicity distributions associated to the underlying event exhibit a KNO scaling.


title: All functions are (locally) $s$-harmonic (up to a small error) - and applications

abstract: The classical and the fractional Laplacians exhibit a number of similarities, but also some rather striking, and sometimes surprising, structural differences. A quite important example of these differences is that any function (regardless of its shape) can be locally approximated by functions with locally vanishing fractional Laplacian, as it was recently proved by Serena Dipierro, Ovidiu Savin and myself. This informal note is an exposition of this result and of some of its consequences.


title: Identifying On-time Reward Delivery Projects with Estimating Delivery Duration on Kickstarter

abstract: In Crowdfunding platforms, people turn their prototype ideas into real products by raising money from the crowd, or invest in someone else's projects. In reward-based crowdfunding platforms such as Kickstarter and Indiegogo, selecting accurate reward delivery duration becomes crucial for creators, backers, and platform providers to keep the trust between the creators and the backers, and the trust between the platform providers and users. According to Kickstarter, 35% backers did not receive rewards on time. Unfortunately, little is known about on-time and late reward delivery projects, and there is no prior work to estimate reward delivery duration. To fill the gap, in this paper, we (i) extract novel features that reveal latent difficulty levels of project rewards; (ii) build predictive models to identify whether a creator will deliver all rewards in a project on time or not; and (iii) build a regression model to estimate accurate reward delivery duration (i.e., how long it will take to produce and deliver all the rewards). Experimental results show that our models achieve good performance -- 82.5% accuracy, 78.1 RMSE, and 0.108 NRMSE at the first 5% of the longest reward delivery duration.


title: Can the early human visual system compete with Deep Neural Networks?

abstract: We study and compare the human visual system and state-of-the-art deep neural networks on classification of distorted images. Different from previous works, we limit the display time to 100ms to test only the early mechanisms of the human visual system, without allowing time for any eye movements or other higher level processes. Our findings show that the human visual system still outperforms modern deep neural networks under blurry and noisy images. These findings motivate future research into developing more robust deep networks.


title: Self-similar groups of type FP_{n}

abstract: We construct new classes of self-similar groups : S-aritmetic groups, affine groups and metabelian groups. Most of the soluble ones are finitely presented and of type FP_{n} for appropriate n.


title: Movement-efficient Sensor Deployment in Wireless Sensor Networks

abstract: We study a mobile wireless sensor network (MWSN) consisting of multiple mobile sensors or robots. Two key issues in MWSNs - energy consumption, which is dominated by sensor movement, and sensing coverage - have attracted plenty of attention, but the interaction of these issues is not well studied. To take both sensing coverage and movement energy consumption into consideration, we model the sensor deployment problem as a constrained source coding problem. %, which can be applied to different coverage tasks, such as area coverage, target coverage, and barrier coverage. Our goal is to find an optimal sensor deployment to maximize the sensing coverage with specific energy constraints. We derive necessary conditions to the optimal sensor deployment with (i) total energy constraint and (ii) network lifetime constraint. Using these necessary conditions, we design Lloyd-like algorithms to provide a trade-off between sensing coverage and energy consumption. Simulation results show that our algorithms outperform the existing relocation algorithms.


title: Solvation in Space-Time: Pre-transition Effects in Trajectory Space

abstract: We demonstrate pre-transition effects in space-time in trajectories of systems in which the dynamics displays a first-order phase transition between distinct dynamical phases. These effects are analogous to those observed for thermodynamic first-order phase transitions, most notably the hydrophobic effect in water. Considering the (infinite temperature) East model as an elementary example, we study the properties of "space-time solvation" by examining trajectories where finite space-time regions are conditioned to be inactive in an otherwise active phase. Analogous to ideal hydrophobic solutes in water, we show that solvating an inactive region of space-time within an active trajectory shows two regimes in the dynamical equivalent of solvation free energy: an "entropic" small solute regime in which uncorrelated fluctuations are sufficient to evacuate activity from the solute, and an "energetic" large solute regime which involves the formation of a solute-induced inactive domain with an associated active-inactive interface bearing a dynamical interfacial tension. We also show that as a result of this dynamical interfacial tension there is a dynamical analog of the hydrophobic collapse that drives the assembly of large hydrophobes in water. We discuss the general relevance of these results to the properties of dynamical fluctuations in systems with slow collective relaxation such as glass formers.


title: HyperENTM: Evolving Scalable Neural Turing Machines through HyperNEAT

abstract: Recent developments within memory-augmented neural networks have solved sequential problems requiring long-term memory, which are intractable for traditional neural networks. However, current approaches still struggle to scale to large memory sizes and sequence lengths. In this paper we show how access to memory can be encoded geometrically through a HyperNEAT-based Neural Turing Machine (HyperENTM). We demonstrate that using the indirect HyperNEAT encoding allows for training on small memory vectors in a bit-vector copy task and then applying the knowledge gained from such training to speed up training on larger size memory vectors. Additionally, we demonstrate that in some instances, networks trained to copy bit-vectors of size 9 can be scaled to sizes of 1,000 without further training. While the task in this paper is simple, these results could open up the problems amendable to networks with external memories to problems with larger memory vectors and theoretically unbounded memory sizes.


title: Explaining Aviation Safety Incidents Using Deep Temporal Multiple Instance Learning

abstract: Although aviation accidents are rare, safety incidents occur more frequently and require a careful analysis to detect and mitigate risks in a timely manner. Analyzing safety incidents using operational data and producing event-based explanations is invaluable to airline companies as well as to governing organizations such as the Federal Aviation Administration (FAA) in the United States. However, this task is challenging because of the complexity involved in mining multi-dimensional heterogeneous time series data, the lack of time-step-wise annotation of events in a flight, and the lack of scalable tools to perform analysis over a large number of events. In this work, we propose a precursor mining algorithm that identifies events in the multidimensional time series that are correlated with the safety incident. Precursors are valuable to systems health and safety monitoring and in explaining and forecasting safety incidents. Current methods suffer from poor scalability to high dimensional time series data and are inefficient in capturing temporal behavior. We propose an approach by combining multiple-instance learning (MIL) and deep recurrent neural networks (DRNN) to take advantage of MIL's ability to learn using weakly supervised data and DRNN's ability to model temporal behavior. We describe the algorithm, the data, the intuition behind taking a MIL approach, and a comparative analysis of the proposed algorithm with baseline models. We also discuss the application to a real-world aviation safety problem using data from a commercial airline company and discuss the model's abilities and shortcomings, with some final remarks about possible deployment directions.


title: Generalized Gaussian Multiterminal Source Coding: The Symmetric Case

abstract: Consider a generalized multiterminal source coding system, where $\ell\choose m$ encoders, each observing a distinct size-$m$ subset of $\ell$ ($\ell\geq 2$) zero-mean unit-variance symmetrically correlated Gaussian sources with correlation coefficient $\rho$, compress their observations in such a way that a joint decoder can reconstruct the sources within a prescribed mean squared error distortion based on the compressed data. The optimal rate-distortion performance of this system was previously known only for the two extreme cases $m=\ell$ (the centralized case) and $m=1$ (the distributed case), and except when $\rho=0$, the centralized system can achieve strictly lower compression rates than the distributed system under all non-trivial distortion constraints. Somewhat surprisingly, it is established in the present paper that the optimal rate-distortion performance of the afore-described generalized multiterminal source coding system with $m\geq 2$ coincides with that of the centralized system for all distortions when $\rho\leq 0$ and for distortions below an explicit positive threshold (depending on $m$) when $\rho>0$. Moreover, when $\rho>0$, the minimum achievable rate of generalized multiterminal source coding subject to an arbitrary positive distortion constraint $d$ is shown to be within a finite gap (depending on $m$ and $d$) from its centralized counterpart in the large $\ell$ limit except for possibly the critical distortion $d=1-\rho$.


title: Central Moments-based Cascaded Lattice Boltzmann Method for Thermal Convective Flows in Three-Dimensions

abstract: Fluid motion driven by thermal effects, such as that due to buoyancy in differentially heated three-dimensional (3D) enclosures, arise in several natural settings and engineering applications. It is represented by the solutions of the Navier-Stokes equations (NSE) in conjunction with the thermal energy transport equation represented as a convection-diffusion equation (CDE) for the temperature field. In this study, we develop new 3D lattice Boltzmann (LB) methods based on central moments and using multiple relaxation times for the three-dimensional, fifteen velocity (D3Q15) lattice, as well as it subset, i.e. the three-dimensional, seven velocity (D3Q7) lattice to solve the 3D CDE for the temperature field in a double distribution function framework. Their collision operators lead to a cascaded structure involving higher order terms resulting in improved stability. In this approach, the fluid motion is solved by another 3D cascaded LB model from prior work. Owing to the differences in the number of collision invariants to represent the dynamics of flow and the transport of the temperature field, the structure of the collision operator for the 3D cascaded LB formulation for the CDE is found to be markedly different from that for the NSE. The new 3D cascaded (LB) models for thermal convective flows are validated for natural convection of air driven thermally on two vertically opposite faces in a cubic cavity enclosure at different Rayleigh numbers against prior numerical benchmark solutions. Results show good quantitative agreement of the profiles of the flow and thermal fields, and the magnitudes of the peak convection velocities as well as the heat transfer rates given in terms of the Nusselt number.


title: Vertex degree sums for perfect matchings in 3-uniform hypergraphs

abstract: We determine the minimum degree sum of two adjacent vertices that ensures a perfect matching in a 3-graph without isolated vertex. More precisely, suppose that $H$ is a 3-uniform hypergraph whose order $n$ is sufficiently large and divisible by $3$. If $H$ contains no isolated vertex and $deg(u)+ deg(v) > \frac{2}{3}n^2-\frac{8}{3}n+2$ for any two vertices $u$ and $v$ that are contained in some edge of $H$, then $H$ contains a perfect matching. This bound is tight.


title: A Robust Accelerated Optimization Algorithm for Strongly Convex Functions

abstract: This work proposes an accelerated first-order algorithm we call the Robust Momentum Method for optimizing smooth strongly convex functions. The algorithm has a single scalar parameter that can be tuned to trade off robustness to gradient noise versus worst-case convergence rate. At one extreme, the algorithm is faster than Nesterov's Fast Gradient Method by a constant factor but more fragile to noise. At the other extreme, the algorithm reduces to the Gradient Method and is very robust to noise. The algorithm design technique is inspired by methods from classical control theory and the resulting algorithm has a simple analytical form. Algorithm performance is verified on a series of numerical simulations in both noise-free and relative gradient noise cases.


title: Minimizing fractional harmonic maps on the real line in the supercritical regime

abstract: This article addresses the regularity issue for minimizing fractional harmonic maps of order $s\in(0,1/2)$ from an interval into a smooth manifold. H\"older continuity away from a locally finite set is established for a general target. If the target is the standard sphere, then H\"older continuity holds everywhere.


title: A curvature flow in the plane with a nonlocal term

abstract: We study the geometric flow of a planar curve driven by its curvature and the normal derivative of its capacity potential. Under a convexity condition that is natural to our problem, we establish long term existence and large time asymptotics of this flow.


title: Spherical particle in a nematic liquid crystal under an external field: the Saturn ring regime

abstract: We consider a nematic liquid crystal occupying the exterior region in R^3 outside of a spherical particle, with radial strong anchoring. Within the context of the Landau-de Gennes theory, we study minimizers subject to a strong external field, modelled by an additional term which favors nematic alignment parallel to the field. When the external field is high enough we obtain a scaling law for the energy. The energy scale corresponds to minimizers concentrating their energy in a boundary layer around the particle, with quadrupolar symmetry. This suggests the presence of a Saturn ring defect around the particle, rather than a dipolar director field typical of a point defect.


title: Hamilton path decompositions of complete multipartite graphs

abstract: We prove that a complete multipartite graph $K$ with $n>1$ vertices and $m$ edges can be decomposed into edge-disjoint Hamilton paths if and only if $\frac m{n-1}$ is an integer and the maximum degree of $K$ is at most $\frac {2m}{n-1}$.


title: Gate tunable third-order nonlinear optical response of massless Dirac fermions in graphene

abstract: Materials with massless Dirac fermions can possess exceptionally strong and widely tunable optical nonlinearities. Experiments on graphene monolayer have indeed found very large third-order nonlinear responses, but the reported variation of the nonlinear optical coefficient by orders of magnitude is not yet understood. A large part of the difficulty is the lack of information on how doping or chemical potential affects the different nonlinear optical processes. Here we report the first experimental study, in corroboration with theory, on third harmonic generation (THG) and four-wave mixing (FWM) in graphene that has its chemical potential tuned by ion-gel gating. THG was seen to have enhanced by ~30 times when pristine graphene was heavily doped, while difference-frequency FWM appeared just the opposite. The latter was found to have a strong divergence toward degenerate FWM in undoped graphene, leading to a giant third-order nonlinearity. These truly amazing characteristics of graphene come from the possibility to gate-control the chemical potential, which selectively switches on and off one- and multi-photon resonant transitions that coherently contribute to the optical nonlinearity, and therefore can be utilized to develop graphene-based nonlinear optoelectronic devices.


title: Bayesian Hypernetworks

abstract: We study Bayesian hypernetworks: a framework for approximate Bayesian inference in neural networks. A Bayesian hypernetwork $\h$ is a neural network which learns to transform a simple noise distribution, $p(\vec\epsilon) = \N(\vec 0,\mat I)$, to a distribution $q(\pp) := q(h(\vec\epsilon))$ over the parameters $\pp$ of another neural network (the "primary network")\@. We train $q$ with variational inference, using an invertible $\h$ to enable efficient estimation of the variational lower bound on the posterior $p(\pp | \D)$ via sampling. In contrast to most methods for Bayesian deep learning, Bayesian hypernets can represent a complex multimodal approximate posterior with correlations between parameters, while enabling cheap iid sampling of~$q(\pp)$. In practice, Bayesian hypernets can provide a better defense against adversarial examples than dropout, and also exhibit competitive performance on a suite of tasks which evaluate model uncertainty, including regularization, active learning, and anomaly detection.


title: Dynamical analysis for a scalar-tensor model with Gauss-Bonnet and non-minimal couplings

abstract: We study the autonomous system for a scalar-tensor model of dark energy with Gauss-Bonnet and non-minimal couplings. The critical points describe important stable asymptotic scenarios including quintessence, phantom and de Sitter attractor solutions. Two functional forms for the coupling functions and the scalar potential were considered: power-law and exponential functions of the scalar field. For the exponential functions the existence of stable quintessence, phantom or de Sitter solutions, allows an asymptotic behavior where the effective Newtonian coupling becomes constant. The phantom solutions could be realized without appealing to ghost degrees of freedom. Transient inflationary and radiation dominated phases can also be described.


title: Chiral surface states on the step edge in a Weyl semimetal

abstract: A Weyl semimetal with a pair of Weyl nodes accommodates chiral states on its flat surface if the Weyl nodes are projected onto two different points in the corresponding surface Brillouin zone. These surface states are collectively referred to as a Fermi arc as they appear to connect the projected Weyl nodes. This statement assumes that translational symmetry is present on the surface and hence electron momentum is a conserved quantity. It is unclear how chiral surface states are modified if the translational symmetry is broken by a particular system structure. Here, focusing on a straight step edge of finite width, we numerically analyze how chiral surface states appear on it. It is shown that the chiral surface states are algebraically (i.e., weakly) localized near the step edge. It is also shown that the appearance of chiral surface states is approximately determined by a simple condition characterized by the number of unit atomic layers constituting the step edge together with the location of the Weyl nodes.


title: On propagation of higher space regularity for non-linear Vlasov equations

abstract: This work is concerned with the broad question of propagation of regularity for smooth solutions to non-linear Vlasov equations. For a class of equations (that includes Vlasov-Poisson and relativistic Vlasov-Maxwell), we prove that higher regularity in space is propagated, locally in time, into higher regularity for the moments in velocity of the solution. This in turn can be translated into some anisotropic Sobolev higher regularity for the solution itself, which can be interpreted as a kind of weak propagation of space regularity. To this end, we adapt the methods introduced in the context of the quasineutral limit of the Vlasov-Poisson system in [D. Han-Kwan and F. Rousset, Ann. Sci. \'Ecole Norm. Sup., 2016].


title: Quench detection on a superconducting radio-frequency cavity

abstract: We study quench detection in superconducting accelerator cavities cooled with He-II. A rigorous mathematical formula is derived to localize the quench position from dynamical data over a finite time interval at a second sound detector.


title: Near-infrared photothermal response of plasmonic gold-coated nanoparticles in tissues

abstract: We propose a new approach to understand the time-dependent temperature increasing process of gold-silica core-shell nanoparticles injected into chicken tissues under near-infrared laser irradiation. Gold nanoshells strongly absorb near-infrared radiations and efficiently transform absorbed energy into heat. Temperature rise given by experiments and numerical calculations based on bioheat transfer are in good agreement. Our work improves the analysis of a recent study [Richardson et al., Nano Lett. 9, 1139 (2009)] by including effects of the medium perfusion on temperature increase. The theoretical analysis can also be used to estimate the distribution of nanoparticles in experimental samples and provide a relative accuracy prediction for the temperature profile of new systems. This methodology would provide a novel and reliable tool for speeding up photothermal investigations and designing state-of-the-art photothermal devices.

categories: cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph physics.bio-ph physics.med-ph


title: Estimating a network from multiple noisy realizations

abstract: Complex interactions between entities are often represented as edges in a network. In practice, the network is often constructed from noisy measurements and inevitably contains some errors. In this paper we consider the problem of estimating a network from multiple noisy observations where edges of the original network are recorded with both false positives and false negatives. This problem is motivated by neuroimaging applications where brain networks of a group of patients with a particular brain condition could be viewed as noisy versions of an unobserved true network corresponding to the disease. The key to optimally leveraging these multiple observations is to take advantage of network structure, and here we focus on the case where the true network contains communities. Communities are common in real networks in general and in particular are believed to be presented in brain networks. Under a community structure assumption on the truth, we derive an efficient method to estimate the noise levels and the original network, with theoretical guarantees on the convergence of our estimates. We show on synthetic networks that the performance of our method is close to an oracle method using the true parameter values, and apply our method to fMRI brain data, demonstrating that it constructs stable and plausible estimates of the population network.


title: Mode-Dependent Damping in Metallic Antiferromagnets Due to Inter-Sublattice Spin Pumping

abstract: Damping in magnetization dynamics characterizes the dissipation of magnetic energy and is essential for improving the performance of spintronics-based devices. While the damping of ferromagnets has been well studied and can be artificially controlled in practice, the damping parameters of antiferromagnetic materials are nevertheless little known for their physical mechanisms or numerical values. Here we calculate the damping parameters in antiferromagnetic dynamics using the generalized scattering theory of magnetization dissipation combined with the first-principles transport computation. For the PtMn, IrMn, PdMn and FeMn metallic antiferromagnets, the damping coefficient associated with the motion of magnetization ($\alpha_m$) is one to three orders of magnitude larger than the other damping coefficient associated with the variation of the N\'eel order ($\alpha_n$), in sharp contrast to the assumptions made in the literature.


title: Higher level Zhu algebras and modules for vertex operator algebras

abstract: Motivated by the study of indecomposable, nonsimple modules for a vertex operator algebra $V$, we study the relationship between various types of $V$-modules and modules for the higher level Zhu algebras for $V$, denoted $A_n(V)$, for $n \in \mathbb{N}$, first introduced by Dong, Li, and Mason in 1998. We resolve some issues that arise in a few theorems previously presented when these algebras were first introduced, and give examples illustrating the need for certain modifications of the statements of those theorems. We establish that whether or not $A_{n-1}(V)$ is isomorphic to a direct summand of $A_n(V)$ affects the types of indecomposable $V$-modules which can be constructed by inducing from an $A_n(V)$-module, and in particular whether there are $V$-modules induced from $A_n(V)$-modules that were not already induced by $A_0(V)$. We give some characterizations of the $V$-modules that can be constructed from such inducings, in particular as regards their singular vectors. To illustrate these results, we discuss two examples of $A_1(V)$: when $V$ is the vertex operator algebra associated to either the Heisenberg algebra or the Virasoro algebra. For these two examples, we show how the structure of $A_1(V)$ in relationship to $A_0(V)$ determines what types of indecomposable $V$-modules can be induced from a module for the level zero versus level one Zhu algebras. We construct a family of indecomposable modules for the Virasoro vertex operator algebra that are logarithmic modules and are not highest weight modules.


title: Atomic optical amplifier with narrow bandwidth optical filtering

abstract: Taking advantages of ultra-narrow bandwidth and high noise rejection performance of the Faraday anomalous dispersion optical filter (FADOF), simultaneously with the coherent amplification of atomic stimulated emission, a stimulated amplified Faraday anomalous dispersion optical filter (SAFADOF) at cesium 1470 nm is realized. The SAFADOF is able to significantly amplify very weak laser signals and reject noise in order to obtain clean signals in strong background. Experiment results show that, for a weak signal of 50 pW, the gain factor can be larger than 25000 (44 dB) within a bandwidth as narrow as 13 MHz. Having this ability to amplify weak signals with low background contribution, the SAFADOF finds outstanding potential applications in weak signal detections.


title: Classification of screening systems for lattice vertex operator algebras

abstract: We study and classify systems of certain screening operators arising in a generalized vertex operator algebra, or more generally an abelian intertwining algebra with an associated vertex operator (super)algebra. Screening pairs arising from weight one primary vectors acting commutatively on a lattice vertex operator algebra (the vacuum module) are classified into four general types, one type of which has been shown to play an important role in the construction and study of certain important families of $\mathcal{W}$-vertex algebras. These types of screening pairs we go on to study in detail through the notion of a system of screeners, which are lattice elements or `screening momenta' which give rise to screening pairs. We classify screening systems for all positive definite integral lattices of rank two, and for all positive definite even lattices of arbitrary rank when these lattices are generated by a screening system.


title: Disentangling the radio emission of the supernova remnant W51C

abstract: We simulate the evolution of supernova remnant (SNR) W51C. The simulation shows the existence of a new northeast edge. We present magnetic field structure of the W51 complex (SNR W51C and two HII regions W51A/B) by employing the 11 cm survey data of Effelsberg. This new edge is identified and overlaps with W51A along the line of sight, which gives a new angular diameter of about 37' for the quasi-circular remnant. In addition, we assemble the OH spectral lines (1612/1665/1720 MHz) towards the complex by employing the newly released THOR (The HI OH Recombination line survey of Milky Way) data. We find that the known 1720 MHz OH maser in the W51B/C overlap area is located away from the detected 1612/1665MHz absorption region. The latter is sitting at the peak of the HII region G49.2-0.35 within W51B.


title: Applications of Economic and Pricing Models for Resource Management in 5G Wireless Networks: A Survey

abstract: This paper presents a comprehensive literature review on applications of economic and pricing theory for resource management in the evolving fifth generation (5G) wireless networks. The 5G wireless networks are envisioned to overcome existing limitations of cellular networks in terms of data rate, capacity, latency, energy efficiency, spectrum efficiency, coverage, reliability, and cost per information transfer. To achieve the goals, the 5G systems will adopt emerging technologies such as massive Multiple-Input Multiple-Output (MIMO), mmWave communications, and dense Heterogeneous Networks (HetNets). However, 5G involves multiple entities and stakeholders that may have different objectives, e.g., high data rate, low latency, utility maximization, and revenue/profit maximization. This poses a number of challenges to resource management designs of 5G. While the traditional solutions may neither efficient nor applicable, economic and pricing models have been recently developed and adopted as useful tools to achieve the objectives. In this paper, we review economic and pricing approaches proposed to address resource management issues in the 5G wireless networks including user association, spectrum allocation, and interference and power management. Furthermore, we present applications of economic and pricing models for wireless caching and mobile data offloading. Finally, we highlight important challenges, open issues and future research directions of applying economic and pricing models to the 5G wireless networks.


title: Edge sampling using network local information

abstract: Edge sampling is an important topic in network analysis. It provides a natural way to reduce network size while retaining desired features of the original network. Sampling methods that only use local information are common in practice as they do not require access to the entire network and can be parallelized easily. Despite promising empirical performances, most of these methods are derived from heuristic considerations and therefore still lack theoretical justification. To address this issue, we study in this paper a simple edge sampling scheme that uses network local information. We show that when local connectivity is sufficiently strong, the sampled network satisfies a strong spectral property. We quantify the strength of local connectivity by a global parameter and relate it to more common network statistics such as the clustering coefficient and network curvature. Based on this result, we also provide sufficient conditions under which random networks and hypergraphs can be sampled efficiently.


title: Residual Connections Encourage Iterative Inference

abstract: Residual networks (Resnets) have become a prominent architecture in deep learning. However, a comprehensive understanding of Resnets is still a topic of ongoing research. A recent view argues that Resnets perform iterative refinement of features. We attempt to further expose properties of this aspect. To this end, we study Resnets both analytically and empirically. We formalize the notion of iterative refinement in Resnets by showing that residual connections naturally encourage features of residual blocks to move along the negative gradient of loss as we go from one block to the next. In addition, our empirical analysis suggests that Resnets are able to perform both representation learning and iterative refinement. In general, a Resnet block tends to concentrate representation learning behavior in the first few layers while higher layers perform iterative refinement of features. Finally we observe that sharing residual layers naively leads to representation explosion and counterintuitively, overfitting, and we show that simple existing strategies can help alleviating this problem.


title: The Law of the Iterated Logarithm for a Class of SPDEs

abstract: After establishing the moderate deviation principle by the Classical Azencott method, we prove the Strassen's compact law of the iterated logarithm (LIL) for a class of stochastic partial differential equations (SPDEs). As an application, we obtain this type of LIL for two population models known as super-Brownian motion and Fleming-Viot process. In addition, the classical LIL is shown for the class of SPDEs and the two population models.


title: Short-Packet Downlink Transmission with Non-Orthogonal Multiple Access

abstract: This work introduces downlink non-orthogonal multiple access (NOMA) into short-packet communications. NOMA has great potential to improve fairness and spectral efficiency with respect to orthogonal multiple access (OMA) for low-latency downlink transmission, thus making it attractive for the emerging Internet of Things. We consider a two-user downlink NOMA system with finite blocklength constraints, in which the transmission rates and power allocation are optimized. To this end, we investigate the trade-off among the transmission rate, decoding error probability, and the transmission latency measured in blocklength. Then, a one-dimensional search algorithm is proposed to resolve the challenges mainly due to the achievable rate affected by the finite blocklength and the unguaranteed successive interference cancellation. We also analyze the performance of OMA as a benchmark to fully demonstrate the benefit of NOMA. Our simulation results show that NOMA significantly outperforms OMA in terms of achieving a higher effective throughput subject to the same finite blocklength constraint, or incurring a lower latency to achieve the same effective throughput target. Interestingly, we further find that with the finite blocklength, the advantage of NOMA relative to OMA is more prominent when the effective throughput targets at the two users become more comparable.


title: Efficient 525nm laser generation in single or double resonant cavity

abstract: This paper reports the results of a study into highly efficient sum frequency generation from 792 and 1556 nm wavelength light to 525 nm wavelength light using either a single or double resonant ring cavity based on a periodically poled potassium titanyl phosphate crystal (PPKTP). By optimizing the cavity$'$s parameters, the maximum power achieved for the resultant 525 nm laser was 263 and 373 mW for the single and double resonant cavity, respectively. The corresponding quantum conversion efficiencies were 8 and 77\% for converting 1556 nm photons to 525 nm photons with the single and double resonant cavity, respectively. The measured intra-cavity single pass conversion efficiency for both configurations was about 5\%. The performances of the sum frequency generation in these two configurations were studied and compared in detail. This work will provide guidelines for optimizing the generation of sum frequency generated laser light for a variety of configurations. The high conversion efficiency achieved in this work will help pave the way for frequency up-conversion of non-classical quantum states, such as the squeezed vacuum and single photon states. The proposed green laser source will be used in our future experiments, which includes a plan to generate two-colour entangled photon pairs and achieve the frequency down-conversion of single photons carrying orbital angular momentum.


title: Higher Order Convergence Rates in Theory of Homogenization III: viscous Hamilton-Jacobi Equations

abstract: In this paper, we establish the higher order convergence rates in periodic homogenization of viscous Hamilton-Jacobi equations, which is convex and grows quadratically in the gradient variable. We observe that although the nonlinear structure governs the first order approximation, the nonlinear effect is absorbed as an external source term of a linear equation in the second and higher order approximation. Moreover, we find that the geometric shape of the initial data has to be chosen carefully according to the effective Hamiltonian, in order to achieve the higher order convergence rates.


title: Retinal Fluid Segmentation and Detection in Optical Coherence Tomography Images using Fully Convolutional Neural Network

abstract: As a non-invasive imaging modality, optical coherence tomography (OCT) can provide micrometer-resolution 3D images of retinal structures. Therefore it is commonly used in the diagnosis of retinal diseases associated with edema in and under the retinal layers. In this paper, a new framework is proposed for the task of fluid segmentation and detection in retinal OCT images. Based on the raw images and layers segmented by a graph-cut algorithm, a fully convolutional neural network was trained to recognize and label the fluid pixels. Random forest classification was performed on the segmented fluid regions to detect and reject the falsely labeled fluid regions. The leave-one-out cross validation experiments on the RETOUCH database show that our method performs well in both segmentation (mean Dice: 0.7317) and detection (mean AUC: 0.985) tasks.


title: Bernoulli numbers and sums of powers of integers of higher order

abstract: We give an expression of polynomials for higher sums of powers of integers via the higher order Bernoulli numbers.


title: Novel Application of Density Estimation Techniques in Muon Ionization Cooling Experiment

abstract: The international Muon Ionization Cooling Experiment (MICE) aims to demonstrate muon beam ionization cooling for the first time and constitutes a key part of the R&D towards a future neutrino factory or muon collider. Beam cooling reduces the size of the phase space volume occupied by the beam. Non-parametric density estimation techniques allow very precise calculation of the muon beam phase-space density and its increase as a result of cooling. These density estimation techniques are investigated in this paper and applied in order to estimate the reduction in muon beam size in MICE under various conditions.


title: GRB 111005A at Z = 0.0133 and the Prospect of Establishing Long-short GRB/GW Association

abstract: GRB 111005A, one long duration gamma-ray burst (GRB) occurred within a metal-rich environment that lacks massive stars with $M_{\rm ZAMS}\geq 15M_\odot$, is not coincident with supernova emission down to stringent limit and thus should be classified as a "long-short" GRB (lsGRB; also known as SN-less long GRB or hybrid GRB), like GRB 060505 and GRB 060614. In this work we show that in the neutron star merger model, the non-detection of the optical/infrared emission of GRB 111005A requires a sub-relativistic neutron-rich ejecta with the mass of $\leq 0.01~M_\odot$, (significantly) less massive than that of GRB 130603B, GRB 060614 and GRB 050709. The lsGRBs are found to have a high rate density and the neutron star merger origin model can be unambiguously tested by the joint observations of the second generation gravitational wave (GW) detectors and the full-sky gamma-ray monitors such as Fermi-GBM and the proposing GECAM. If no lsGRB/GW association is observed in 2020s, alternative scenarios have to be systematically investigated. With the detailed environmental information achievable for the very-nearby events, a novel kind of merger or explosion origin may be identified.


title: Multimodal and Multiscale Deep Neural Networks for the Early Diagnosis of Alzheimer's Disease using structural MR and FDG-PET images

abstract: Alzheimer's Disease (AD) is a progressive neurodegenerative disease. Amnestic mild cognitive impairment (MCI) is a common first symptom before the conversion to clinical impairment where the individual becomes unable to perform activities of daily living independently. Although there is currently no treatment available, the earlier a conclusive diagnosis is made, the earlier the potential for interventions to delay or perhaps even prevent progression to full-blown AD. Neuroimaging scans acquired from MRI and metabolism images obtained by FDG-PET provide in-vivo view into the structure and function (glucose metabolism) of the living brain. It is hypothesized that combining different image modalities could better characterize the change of human brain and result in a more accuracy early diagnosis of AD. In this paper, we proposed a novel framework to discriminate normal control(NC) subjects from subjects with AD pathology (AD and NC, MCI subjects convert to AD in future). Our novel approach utilizing a multimodal and multiscale deep neural network was found to deliver a 85.68\% accuracy in the prediction of subjects within 3 years to conversion. Cross validation experiments proved that it has better discrimination ability compared with results in existing published literature.


title: Retinal Vasculature Segmentation Using Local Saliency Maps and Generative Adversarial Networks For Image Super Resolution

abstract: We propose an image super resolution(ISR) method using generative adversarial networks (GANs) that takes a low resolution input fundus image and generates a high resolution super resolved (SR) image upto scaling factor of $16$. This facilitates more accurate automated image analysis, especially for small or blurred landmarks and pathologies. Local saliency maps, which define each pixel's importance, are used to define a novel saliency loss in the GAN cost function. Experimental results show the resulting SR images have perceptual quality very close to the original images and perform better than competing methods that do not weigh pixels according to their importance. When used for retinal vasculature segmentation, our SR images result in accuracy levels close to those obtained when using the original images.


title: Charge exchange in galaxy clusters

abstract: Though theoretically expected, the charge exchange emission from galaxy clusters has not yet been confidently detected. Accumulating hints were reported recently, including a rather marginal detection with the Hitomi data of the Perseus cluster. As suggested in Gu et al. (2015), a detection of charge exchange line emission from galaxy clusters would not only impact the interpretation of the newly-discovered 3.5 keV line, but also open up a new research topic on the interaction between hot and cold matter in clusters. We aim to perform the most systematic search for the O VIII charge exchange line in cluster spectra using the RGS on board XMM. We introduce a sample of 21 clusters observed with the RGS. The dominating thermal plasma emission is modeled and subtracted with a two-temperature CIE component, and the residuals are stacked for the line search. The systematic uncertainties in the fits are quantified by refitting the spectra with a varying continuum and line broadening. By the residual stacking, we do find a hint of a line-like feature at 14.82 A, the characteristic wavelength expected for oxygen charge exchange. This feature has a marginal significance of 2.8 sigma, and the average equivalent width is 2.5E-4 keV. We further demonstrate that the putative feature can be hardly affected by the systematic errors from continuum modelling and instrumental effects, or the atomic uncertainties of the neighbouring thermal lines. Assuming a realistic temperature and abundance pattern, the physical model implied by the possible oxygen line agrees well with the theoretical model proposed previously to explain the reported 3.5 keV line. If the charge exchange source indeed exists, we would expect that the oxygen abundance is potentially overestimated by 8-22% in previous X-ray measurements which assumed pure thermal lines.

categories: astro-ph.HE astro-ph.CO astro-ph.GA physics.atom-ph


title: The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

abstract: We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.


title: Tunable Intrinsic Plasmons due to Band Inversion in Topological Materials

abstract: The band inversion has led to rich physical effects in both topological insulators and topological semimetals. It has been found that the inverted band structure with the Mexican-hat dispersion could enhance the interband correlation leading to a strong intrinsic plasmon excitation. Its frequency ranges from several $\mathrm{meV}$ to tens of $\mathrm{meV}$ and can be effectively tuned by the external fields. The electron-hole asymmetric term splits the peak of the plasmon excitation into double peaks. The fate and properties of this plasmon excitation can also act as a probe to characterize the topological phases even in the lightly doped systems. We numerically demonstrate the impact of the band inversion on plasmon excitations in magnetically doped thin films of three-dimensional strong topological insulators, V- or Cr-doped (Bi, Sb)$_2$Te$_3$, which support the quantum anomalous Hall states. Our work thus sheds some new light on the potential applications of topological materials in plasmonics.


title: A note on automatic continuity

abstract: We present new results regarding automatic continuity, unifying some diagonalization concepts that have been developed over the years. For example, any homomorphism from a completely metrizable topological group to Thompson's group $F$ has open kernel. A similar claim holds when $F$ is replaced with a Baumslag-Solitar group or a torsion-free word hyperbolic group.


title: A Unified Scheme to Accelerate Adaptive Cubic Regularization and Gradient Methods for Convex Optimization

abstract: In this paper we propose a unified two-phase scheme for convex optimization to accelerate: (1) the adaptive cubic regularization methods with exact/inexact Hessian matrices, and (2) the adaptive gradient method, without any knowledge of the Lipschitz constants for the gradient or the Hessian. This is achieved by tuning the parameters used in the algorithm $\textit{adaptively}$ in its process of progression, which can be viewed as a relaxation over the existing algorithms in the literature. Under the assumption that the sub-problems can be solved approximately, we establish overall iteration complexity bounds for three newly proposed algorithms to obtain an $\epsilon$-approximate solution. Specifically, we show that the adaptive cubic regularization methods with the exact/inexact Hessian matrix both achieve an iteration complexity in the order of $O\left( 1 / \epsilon^{1/3} \right)$, which matches that of the original accelerated cubic regularization method presented in [Nesterov-2008-Accelerating] assuming the availability of the exact Hessian information and the Lipschitz constants, and that the sub-problems are solved to optimality. Under the same two-phase adaptive acceleration framework, the gradient method achieves an iteration complexity in the order $O\left( 1 / \epsilon^{1/2} \right)$, which is known to be best possible (cf. Nesterov-2013-Introductory). Our numerical experiment results show a clear effect of acceleration displayed in the adaptive Newton method with cubic regularization on a set of regularized logistic regression instances.


title: Linear guided modes and Whitham-Boussinesq model for variable topogra

abstract: In this article we study two classical linear water wave problems, i) normal modes of infinite straight channels of bounded constant cross-section, and ii) trapped longitudinal modes in domains with unbounded constant cross-section. Both problems can be stated using linearized free surface potential flow theory, and our goal is to compare known analytic solutions in the literature to numerical solutions obtained using an ad-hoc but simple approximation of the non-local Dirichlet-Neumann operator for linear waves proposed in [vargas2016whitham]. To study normal modes in channels with bounded cross-section we consider special symmetric triangular cross-sections, namely symmetric triangles with sides inclined at $45^{\circ}$ and $60^{\circ}$ to the vertical, and compare modes obtained using the non-local Dirichlet-Neumann operator to known semi-exact analytic expressions by Lamb [lamb1932hydrodynamics], Macdonald [macdonald1893waves] , Greenhill [greenhill1887wave], Packham [packham1980small], and Groves [groves1994hamiltonian]. These geometries have slopping beach boundaries that should in principle limit the applicability of our approximate Dirichlet-Neumann operator. We nevertheless see that the operator gives remarkably close results for even modes, while for odd modes we have some discrepancies near the boundary. For trapped longitudinal modes in domains with an infinite cross-section we consider a piecewise constant depth profile and compare modes computed with the nonlocal operator modes to known analytic solutions of linearized shallow water theory by Miles [miles1972wave], Lin, Juang and Tsay [lin2001anomalous], see also [mei2005theory]. This is a problem of significant geophysical interest, and the proposed model is shows to give quantitatively similar results for the lowest trapped modes.


title: Dynamics of a quasiparticle in the $\alpha$-T$_3$ model: Role of pseudospin polarization and transverse magnetic field on \textbf{\textit{zitterbewegung}}

abstract: We consider the $\alpha$-$T_3$ model which provides a smooth crossover between the honeycomb lattice with pseudospin $1/2$ and the dice lattice with pseudospin $1$ through the variation of a parameter $\alpha$. We study the dynamics of a wave packet representing a quasiparticle in the $\alpha$-T$_3$ model with zero and finite transverse magnetic field. For zero field, it is shown that the wave packet undergoes a transient $zitterbewegung$ (ZB). Various features of ZB depending on the initial pseudospin polarization of the wave packet have been revealed. For an intermediate value of the parameter $\alpha$ i.e. for $0<\alpha<1$ the resulting ZB consists of two distinct frequencies when the wave packet was located initially in $rim$ site. However, the wave packet exhibits single frequency ZB for $\alpha=0$ and $\alpha=1$. It is also unveiled that the frequency of ZB corresponding to $\alpha=1$ gets exactly half of that corresponding to the $\alpha=0$ case. On the other hand, when the initial wave packet was in $hub$ site, the ZB consists of only one frequency for all values of $\alpha$. Using stationary phase approximation we find analytical expression of velocity average which can be used to extract the associated timescale over which the transient nature of ZB persists. On the contrary the wave packet undergoes permanent ZB in presence of a transverse magnetic field. Due to the presence of large number of Landau energy levels the oscillations in ZB appear to be much more complicated. The oscillation pattern depends significantly on the initial pseudospin polarization of the wave packet. Furthermore, it is revealed that the number of the frequency components involved in ZB depends on the parameter $\alpha$.


title: Spatially Heterogeneous Dynamics in a Metallic Glass Forming Liquid Imaged by Electron Correlation Microscopy

abstract: Supercooled liquids exhibit spatial heterogeneity in the dynamics of their fluctuating atomic arrangements. The length and time scales of the heterogeneous dynamics are central to the glass transition and influence nucleation and growth of crystals from the liquid. We report direct experimental visualization of the spatially heterogeneous dynamics as a function of temperature in the supercooled liquid state of a Pt-based metallic glass, using electron correlation microscopy with sub-nanometer resolution. An experimental four point space-time intensity correlation function demonstrates a growing dynamic correlation length, $\xi$, upon cooling of the liquid toward the glass transition temperature. $\xi$ as a function of the relaxation time $\tau$ data are in the good agreement with the Adam-Gibbs, inhomogeneous mode coupling theory and random first order transition theory of the glass transition. The same experiments demonstrate the existence of a nanometer thickness near surface layer with order of magnitude shorter relaxation time than inside the bulk.


title: Sparse Weighted Canonical Correlation Analysis

abstract: Given two data matrices $X$ and $Y$, sparse canonical correlation analysis (SCCA) is to seek two sparse canonical vectors $u$ and $v$ to maximize the correlation between $Xu$ and $Yv$. However, classical and sparse CCA models consider the contribution of all the samples of data matrices and thus cannot identify an underlying specific subset of samples. To this end, we propose a novel sparse weighted canonical correlation analysis (SWCCA), where weights are used for regularizing different samples. We solve the $L_0$-regularized SWCCA ($L_0$-SWCCA) using an alternating iterative algorithm. We apply $L_0$-SWCCA to synthetic data and real-world data to demonstrate its effectiveness and superiority compared to related methods. Lastly, we consider also SWCCA with different penalties like LASSO (Least absolute shrinkage and selection operator) and Group LASSO, and extend it for integrating more than three data matrices.


title: A Brief Review of Ferroelectric Control of Magnetoresistance in Organic Spin Valves

abstract: Magnetoelectric coupling has been a trending research topic in both organic and inorganic materials and hybrids. The concept of controlling magnetism using an electric field is particularly appealing in energy efficient applications. In this spirit, ferroelectricity has been introduced to organic spin valves to manipulate the magneto transport, where the spin transport through the ferromagnet/organic spacer interfaces (spinterface) are under intensive study. The ferroelectric materials in the organic spin valves provide a knob to vary the interfacial energy alignment and the interfacial crystal structures, both are critical for the spin transport. In this review, we first go over the basic concepts of spin transport in organic spin valves. Then we introduce the recent efforts of controlling magnetoresistance of organic spin valves using ferroelectricity, where the ferroelectric material is either inserted as an interfacial layer or used as a spacer material. The realization of the ferroelectric control of magneto transport in organic spin valve, advances our understanding in the spin transport through the ferromagnet/organic interface and suggests more functionality of organic spintronic devices.