title: Biological Regulatory Networks are Minimally Frustrated

abstract: Characterization of the differences between biological and random networks can reveal the design principles that enable the robust realization of crucial biological functions including the establishment of different cell types. Previous studies, focusing on identifying topological features that are present in biological networks but not in random networks, have, however, provided few functional insights. We use a Boolean modeling framework and ideas from spin glass literature to identify functional differences between five real biological networks and random networks with similar topological features. We show that minimal frustration is a fundamental property that allows biological networks to robustly establish cell types and regulate cell fate choice, and this property can emerge in complex networks via Darwinian evolution. The study also provides clues regarding how the regulation of cell fate choice can go awry in a disease like cancer and lead to the emergence of aberrant cell types.


title: Variable warm dust around the Herbig Ae star HD~169142: Birth of a ring?

abstract: The Herbig Ae star HD 169142 is known to have a gaseous disk with a large inner hole, and also a photometrically variable inner dust component in the sub-au region. Following up our previous analysis, we further studied the temporal evolution of inner dust around HD 169142, which may provide information on the evolution from late-stage protoplanetary disks to debris disks. We used near-infrared interferometric observations obtained with VLTI/PIONIER to constrain the dust distribution at three epochs spanning six years. We also studied the photometric variability of HD 169142 using our optical-infrared observations and archival data. Our results indicate that a dust ring at ~0.3 au formed at some time between 2013 and 2018, and then faded (but did not completely disappear) by 2019. The short-term variability resembles that observed in extreme debris disks, and is likely related to short-lived dust of secondary origin, though variable shadowing from the inner ring could be an alternative interpretation. If confirmed, this is the first direct detection of secondary dust production inside a protoplanetary disk.


title: Omega and Sharpe ratio

abstract: Omega ratio, defined as the probability-weighted ratio of gains over losses at a given level of expected return, has been advocated as a better performance indicator compared to Sharpe and Sortino ratio as it depends on the full return distribution and hence encapsulates all information about risk and return. We compute Omega ratio for the normal distribution and show that under some distribution symmetry assumptions, the Omega ratio is oversold as it does not provide any additional information compared to Sharpe ratio. Indeed, for returns that have elliptic distributions, we prove that the optimal portfolio according to Omega ratio is the same as the optimal portfolio according to Sharpe ratio. As elliptic distributions are a weak form of symmetric distributions that generalized Gaussian distributions and encompass many fat tail distributions, this reduces tremendously the potential interest for the Omega ratio.


title: On $C$-compact orthogonally additive operators

abstract: We consider $C$-compact orthogonally additive operators in vector lattices. After providing some examples of $C$-compact orthogonally additive operators on a vector lattice with values in a Banach space we show that the set of those operators is a projection band in the Dedekind complete vector lattice of all regular orthogonally additive operators. In second part of the article we introduce a new class of vector lattices, called $C$-complete, and show that any laterally-to-norm continuous $C$-compact orthogonally additive operator from a $C$-complete vector lattice to a Banach space is narrow, which generalizes a result of Pliev and Popov.


title: Local geometric properties in quasi-normed Orlicz spaces

abstract: Several local geometric properties of Orlicz space $L_\phi$ are presented for an increasing Orlicz function $\phi$ which is not necessarily convex, and thus $L_\phi$ does not need to be a Banach space. In addition to monotonicity of $\phi$ it is supposed that $\phi(u^{1/p})$ is convex for some $p>0$ which is equivalent to that its lower Matuszewska-Orlicz index $\alpha_\phi>0$. Such spaces are locally bounded and are equipped with natural quasi-norms. Therefore many local geometric properties typical for Banach spaces can also be studied in those spaces. The techniques however have to be different, since duality theory cannot be applied in this case. In this article we present complete criteria, in terms of growth conditions of $\phi$, for $L_\phi$ to have type $0<p\le2$, cotype $q\ge 2$, to be (order) $p$-convex or $q$-concave, to have an upper $p$-estimate or a lower $q$-estimate, for $0<p,q<\infty$. We provide detailed proofs of most results, avoiding appealing to general not necessary theorems.


title: On 3-dimensional homotopy quantum field theory III: comparison of two approaches

abstract: Let G be a discrete group and C be an additive spherical G-fusion category. We prove that the state sum 3-dimensional HQFT derived from C is isomorphic to the surgery 3-dimensional HQFT derived from the G-center of C.


title: Fantastic Four: Differentiable Bounds on Singular Values of Convolution Layers

abstract: In deep neural networks, the spectral norm of the Jacobian of a layer bounds the factor by which the norm of a signal changes during forward/backward propagation. Spectral norm regularizations have been shown to improve generalization, robustness and optimization of deep learning methods. Existing methods to compute the spectral norm of convolution layers either rely on heuristics that are efficient in computation but lack guarantees or are theoretically-sound but computationally expensive. In this work, we obtain the best of both worlds by deriving {\it four} provable upper bounds on the spectral norm of a standard 2D multi-channel convolution layer. These bounds are differentiable and can be computed efficiently during training with negligible overhead. One of these bounds is in fact the popular heuristic method of Miyato et al. (multiplied by a constant factor depending on filter sizes). Each of these four bounds can achieve the tightest gap depending on convolution filters. Thus, we propose to use the minimum of these four bounds as a tight, differentiable and efficient upper bound on the spectral norm of convolution layers. We show that our spectral bound is an effective regularizer and can be used to bound either the lipschitz constant or curvature values (eigenvalues of the Hessian) of neural networks. Through experiments on MNIST and CIFAR-10, we demonstrate the effectiveness of our spectral bound in improving generalization and provable robustness of deep networks.


title: Neutrino mass generation from the perspective of presymmetry

abstract: The Standard Model (SM) with one right-handed neutrino per generation is revisited with presymmetry being the global $U(1)_{B-L}$ symmetry of an electroweak theory of leptons and quarks with initially postulated symmetric fractional charges. The cancellation of gauge anomalies and the non-perturbative normalization of lepton charges proceed through the mixing of local and topological charges, the global $B-L$ measuring the induced charge associated with a unit of topological charge, and the mathematical replacement of the original fractional charges with the experimentally observed ones. The $U(1)_{B-L}$ symmetry of the SM with Dirac neutrinos is seen as a residual presymmetry. High-scale and low-scale seesaw mechanisms proposed to explain the mass of neutrinos are examined from the perspective of presymmetry, be they of Majorana or pseudo-Dirac type. We find that the tiny mass splitting in pseudo-Dirac neutrinos and the mass of heavy neutrinos ride on the opposite ends of the seesaw. We show that pseudo-Dirac neutrinos contain extra sterile neutrinos with imprints of presymmetry and for heavy ones we get constraints favoring the low-scale linear seesaw over the inverse variant.


title: Jackiw-Rebbi states and trivial states in interfaced binary waveguide arrays with cubic-quintic nonlinearity

abstract: We systematically investigate two types of localized states - one is the optical analogue of the quantum relativistic Jackiw-Rebbi states, and the other is the trivial localized state - in interfaced binary waveguide arrays in the presence of cubic-quintic nonlinearity. By using the shooting method we can exactly calculate the profiles of these nonlinear localized states. Like in the case with Kerr nonlinearity, we demonstrate that these nonlinear localized states in the case with cubic-quintic nonlinearity also have a distinguishing feature which is completely different from all other well-known nonlinear localized structures in other media. Namely, the profiles of nonlinear localized states with higher peak amplitudes in interfaced binary waveguide arrays can totally envelope those with lower peak amplitudes. We show that high values of the saturation nonlinearity parameter can help to generate and stabilize these intense localized states during propagation, especially in the case with negative coefficient for the cubic nonlinearity term.


title: Study on stability of the first kind of soft-matter quasicrystals

abstract: Based on extended free energy of soft-matter quasicrystals and the variation principle on thermodynamic stability, this study reports the results on stability of the first kind of soft-matter quasicrystals. They are dependent only upon the material constants, and quite simple and intuitive, the material constants can be measured by experiments. The results are significant in studying thermodynamics of the matter.


title: An Algorithm for Strong Stability in the Student-Project Allocation Problem with Ties

abstract: We study a variant of the Student-Project Allocation problem with lecturer preferences over Students where ties are allowed in the preference lists of students and lecturers (SPA-ST). We investigate the concept of strong stability in this context. Informally, a matching is strongly stable if there is no student and lecturer $l$ such that if they decide to form a private arrangement outside of the matching via one of $l$'s proposed projects, then neither party would be worse off and at least one of them would strictly improve. We describe the first polynomial-time algorithm to find a strongly stable matching or to report that no such matching exists, given an instance of SPA-ST. Our algorithm runs in $O(m^2)$ time, where $m$ is the total length of the students' preference lists.


title: The p-Laplacian equation in a rough thin domain with terms concentrating on the boundary

abstract: In this work we use reiterated homogenization and unfolding operator approach to study the asymptotic behavior of the solutions of the $p$-Laplacian equation with Neumann boundary conditions set in a rough thin domain with concentrated terms on the boundary. We study weak, resonant and high roughness, respectively. In the three cases, we deduce the effective equation capturing the dependence on the geometry of the thin channel and the neighborhood where the concentrations take place.


title: Orbital Tuning of Tunnel Coupling in InAs/InP Nanowire Quantum Dots

abstract: We report results on the control of barrier transparency in InAs/InP nanowire quantum dots via the electrostatic control of the device electron states. Recent works demonstrated that barrier transparency in this class of devices displays a general trend just depending on the total orbital energy of the trapped electrons. We show that a qualitatively different regime is observed at relatively low filling numbers, where tunneling rates are rather controlled by the axial configuration of the electron orbital. Transmission rates versus filling are further modified by acting on the radial configuration of the orbitals by means of electrostatic gating, and the barrier transparency for the various orbitals is found to evolve as expected from numerical simulations. The possibility to exploit this mechanism to achieve a controlled continuous tuning of the tunneling rate of an individual Coulomb blockade resonance is discussed.


title: Science Opportunities from Observations of the Interstellar Neutral Gas with Adjustable Boresight Direction

abstract: The interstellar neutral (ISN) gas enters the heliosphere and is detected at a few au from the Sun, as demonstrated by Ulysses and the nterstellar Boundary Explorer (IBEX) missions. Ulysses observed ISN gas from different vantage points in a polar orbit from 1994 to 2007, while IBEX has been observing in an Earth orbit in a fixed direction relative to the Sun from 2009. McComas et al. 2018 reported about an IMAP-Lo detector on board the Interstellar Mapping and Acceleration Probe (IMAP), with an ability to track the ISN flux in the sky. We present observation geometries for ISN gas for a detector with the capability to adjust the boresight direction along the Earth orbit over a year within a multichoice ISN observation scheme. We study science opportunities from the observations as a function of time during a year and the phase of solar activity. We identify observation geometries and determine the observation seasons separately for various ISN species and populations. We find that using an adjustable viewing direction allows for ISN gas observations in the upwind hemisphere, where the signal is not distorted by gravitational focusing, in addition to the viewing of ISN species throughout the entire year. Moreover, we demonstrate that with appropriately adjusted observation geometries, primary and secondary populations can be fully separated. Additionally, we show that atoms of ISN gas on indirect trajectories are accessible for detection, and we present their impact on the study of the ionization rates for ISN species.


title: A low dimensional model for bike sharing demand forecasting

abstract: Big, transport-related datasets are nowadays publicly available, which makes data-driven mobility analysis possible. Trips with their origins, destinations and travel times are collected in publicly available big databases, which allows for a deeper and richer understanding of mobility patterns. This paper proposes a low dimensional approach to combine these data sources with weather data in order to forecast the daily demand for Bike Sharing Systems (BSS). The core of this approach lies in the proposed clustering technique, which reduces the dimension of the problem and, differently from other machine learning techniques, requires limited assumptions on the model or its parameters. The proposed clustering technique synthesizes mobility data quantitatively (number of trips) and spatially (mean trip origin and destination). This allows identifying recursive mobility patterns that - when combined with weather data - provide accurate predictions of the demand. The method is tested with real-world data from New York City. We synthesize more than four million trips into vectors of movement, which are then combined with weather data to forecast the daily demand at a city-level. Results show that, already with a one-parameters model, the proposed approach provides accurate predictions.


title: Search for lepton flavour violating decays of a neutral heavy Higgs boson to $\mu\tau$ and e$\tau$ in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search for lepton flavour violating decays of a neutral non-standard-model Higgs boson in the $\mu\tau$ and e$\tau$ decay modes is presented. The search is based on proton-proton collisions at a center of mass energy $\sqrt{s} =$ 13 TeV collected with the CMS detector in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The $\tau$ leptons are reconstructed in the leptonic and hadronic decay modes. No signal is observed in the mass range 200-900 GeV. At 95% confidence level, the observed (expected) upper limits on the production cross section multiplied by the branching fraction vary from 51.9 (57.4) fb to 1.6 (2.1) fb for the $\mu\tau$ and from 94.1 (91.6) fb to 2.3 (2.3) fb for the e$\tau$ decay modes.


title: Nonvanishing of Dirichlet L-functions, II

abstract: We show that for at least $\frac{5}{13}$ of the primitive Dirichlet characters $\chi$ of large prime modulus, the central value $L(\frac{1}{2},\chi)$ does not vanish, improving on the previous best known result of $\frac{3}{8}$.


title: Design and Autonomous Stabilization of a Ballistically Launched Multirotor

abstract: Aircraft that can launch ballistically and convert to autonomous, free flying drones have applications in many areas such as emergency response, defense, and space exploration, where they can gather critical situational data using onboard sensors. This paper presents a ballistically launched, autonomously stabilizing multirotor prototype (SQUID, Streamlined Quick Unfolding Investigation Drone) with an onboard sensor suite, autonomy pipeline, and passive aerodynamic stability. We demonstrate autonomous transition from passive to vision based, active stabilization, confirming the ability of the multirotor to autonomously stabilize after a ballistic launch in a GPS denied environment.


title: Commensurate States and Pattern Switching via Liquid Crystal Skyrmions Trapped in a Square Lattice

abstract: Using continuum based simulations we show that a rich variety of skyrmion liquid crystal states can be realized in the presence of a periodic obstacle array. As a function of the number of skyrmions per obstacle we find hexagonal, square, dimer, trimer and quadrimer ordering, where the $n$-mer structures are a realization of a molecular crystal state of skyrmions. As a function of external field and obstacle radius we show that there are transitions between the different crystalline states as well as mixed and disordered structures. We discuss how these states are related to commensurate effects seen in other systems, such as vortices in type-II superconductors and colloids interacting with two dimensional substrates.


title: The evolution of Astronomical Observatory design

abstract: This work addresses the development of the astronomical observatory all through history, from an architectural point of view, as a building in relation to the observing instruments and their functioning as a heterogeneous work center. We focused on 32 observatories (in the period 1259-2007) and carefully analyzed the architectures. Considering the impact of the construction itself or its facilities on the results of the research (thermal or structural stability, poor weather protection, turbulence, etc.), there is little attention paid to theories or studies of the architectural or construction aspects of the observatories. Therefore, this work aims to present a theoretical-critical contribution that, at least, invites the reflection of those involved in the development of astronomical observatories in the future.


title: Initial state fluctuations of QCD conserved charges in heavy-ion collisions

abstract: We initialize the Quantum Chromodynamic conserved charges of baryon number, strangeness, and electric charge arising from gluon splitting into quark-antiquark pairs for the initial conditions of relativistic heavy-ion collisions. A new Monte Carlo procedure that can sample from a generic energy density profile is presented, called Initial Conserved Charges in Nuclear Geometry (ICCING), based on quark and gluon multiplicities derived within the color glass condensate (CGC) effective theory. We find that while baryon number and electric charge have nearly identical geometries to the energy density profile, the initial strangeness distribution is considerable more eccentric and is produced primarily at the hot spots corresponding to temperatures of $T\gtrsim 400$ MeV for PbPb collisions at $\sqrt{s_{NN}}=5.02$ TeV.


title: Joint Modeling of Local and Global Temporal Dynamics for Multivariate Time Series Forecasting with Missing Values

abstract: Multivariate time series (MTS) forecasting is widely used in various domains, such as meteorology and traffic. Due to limitations on data collection, transmission, and storage, real-world MTS data usually contains missing values, making it infeasible to apply existing MTS forecasting models such as linear regression and recurrent neural networks. Though many efforts have been devoted to this problem, most of them solely rely on local dependencies for imputing missing values, which ignores global temporal dynamics. Local dependencies/patterns would become less useful when the missing ratio is high, or the data have consecutive missing values; while exploring global patterns can alleviate such problems. Thus, jointly modeling local and global temporal dynamics is very promising for MTS forecasting with missing values. However, work in this direction is rather limited. Therefore, we study a novel problem of MTS forecasting with missing values by jointly exploring local and global temporal dynamics. We propose a new framework LGnet, which leverages memory network to explore global patterns given estimations from local perspectives. We further introduce adversarial training to enhance the modeling of global temporal distribution. Experimental results on real-world datasets show the effectiveness of LGnet for MTS forecasting with missing values and its robustness under various missing ratios.


title: Titan: A Parallel Asynchronous Library for Multi-Agent and Soft-Body Robotics using NVIDIA CUDA

abstract: While most robotics simulation libraries are built for low-dimensional and intrinsically serial tasks, soft-body and multi-agent robotics have created a demand for simulation environments that can model many interacting bodies in parallel. Despite the increasing interest in these fields, no existing simulation library addresses the challenge of providing a unified, highly-parallelized, GPU-accelerated interface for simulating large robotic systems. Titan is a versatile CUDA-based C++ robotics simulation library that employs a novel asynchronous computing model for GPU-accelerated simulations of robotics primitives. The innovative GPU architecture design permits simultaneous optimization and control on the CPU while the GPU runs asynchronously, enabling rapid topology optimization and reinforcement learning iterations. Kinematics are solved with a massively parallel integration scheme that incorporates constraints and environmental forces. We report dramatically improved performance over CPU-based baselines, simulating as many as 300 million primitive updates per second, while allowing flexibility for a wide range of research applications. We present several applications of Titan to high-performance simulations of soft-body and multi-agent robots.


title: A fully-coupled computational framework for large-scale simulation of fluid-driven fracture propagation on parallel computers

abstract: The propagation of cracks driven by a pressurized fluid emerges in several areas of engineering, including structural, geotechnical, and petroleum engineering. We present a robust numerical framework to simulate fluid-driven fracture propagation that addresses the challenges emerging in the simulation of this complex coupled nonlinear hydro-mechanical response. We observe that the numerical difficulties stem from the strong nonlinearities present in the fluid equations as well as those associated with crack propagation, from the quasi-static nature of the problem, and from the a priori unknown and potentially intricate crack geometries that may arise. An additional challenge is the need for large scale simulation owing to the mesh resolution requirements and the expected 3D character of the problem in practical applications. To address these challenges we model crack propagation with a high-order hybrid discontinuous Galerkin / cohesive zone model framework, which has proven massive scalability properties, and we model the lubrication flow inside the propagating cracks using continuous finite elements, furnishing a fully-coupled discretization of the solid and fluid equations. The parallel approach for solving the linearized coupled problem consists of standard iterative solvers based on domain decomposition. The resulting computational approach provides the ability to conduct highly-resolved and quasi-static simulations of fluid-driven fracture propagation with unspecified crack path. We conduct a series of numerical tests to verify the computational framework against known analytical solutions in the toughness and viscosity dominated regimes and we demonstrate its performance in terms of robustness and parallel scalability, enabling simulations of several million degrees of freedom on hundreds of processors.


title: A Laterally Vibrating Lithium Niobate MEMS Resonator Array Operating at 500{\deg}C in Air

abstract: This paper is the first report of the high-temperature characteristics of a laterally vibrating piezoelectric lithium niobate (LiNbO$_{3}$) MEMS resonator array up to 500{\deg}C in air. After a high-temperature burn-in treatment, device quality factor (Q) is enhanced to 508 and the resonance shifts to a lower frequency and remains stable up to 500{\deg}C. During subsequent in situ high-temperature testing, the resonant frequencies of two coupled shear horizontal (SH0) modes in the array are 87.36 MHz and 87.21 MHz at 25{\deg}C and 84.56 MHz and 84.39 MHz at 500{\deg}C, correspondingly, representing a -3% shift in frequency over the temperature range. Upon cooling to room temperature, the resonant frequency returns to 87.36 MHz, demonstrating recoverability of device performance. The first- and second-order temperature coefficient of frequency (TCF) are found to be -95.27 ppm/{\deg}C and 57.5 ppb/{\deg}C$^{2}$ for resonant mode A, and -95.43 ppm/{\deg}C and 55.8 ppb/{\deg}C$^{2}$ for resonant mode B, respectively. The temperature-dependent quality factor (Q) and electromechanical coupling coefficient ($k_{t}^{2}$) are extracted and reported. Device Q decreases to 334 after high-temperature exposure, while $k_{t}^{2}$ increases to 12.40%. This work supports the use of piezoelectric LiNbO$_{3}$ as a material platform for harsh environment radio-frequency (RF) resonant sensors (e.g. temperature and infrared).


title: Formation of compact systems of super-Earths via dynamical instabilities and giant impacts

abstract: NASA's Kepler mission discovered $\sim700$ planets in multi-planet systems containing 3 or more transiting bodies, many of which are super-Earths and mini-Neptunes in compact configurations. Using $N$-body simulations, we examine the in situ, final stage assembly of multi-planet systems via the collisional accretion of protoplanets. Our initial conditions are constructed using a subset of the Kepler 5-planet systems as templates. Two different prescriptions for treating planetary collisions are adopted. The simulations address numerous questions: do the results depend on the accretion prescription?; do the resulting systems resemble the Kepler systems, and do they reproduce the observed distribution of planetary multiplicities when synthetically observed?; do collisions lead to significant modification of protoplanet compositions, or to stripping of gaseous envelopes?; do the eccentricity distributions agree with those inferred for the Kepler planets? We find the accretion prescription is unimportant in determining the outcomes. The final planetary systems look broadly similar to the Kepler templates adopted, but the observed distributions of planetary multiplicities or eccentricities are not reproduced, because scattering does not excite the systems sufficiently. In addition, we find that $\sim 1$% of our final systems contain a co-orbital planet pair in horseshoe or tadpole orbits. Post-processing the collision outcomes suggests they would not significantly change the ice fractions of initially ice-rich protoplanets, but significant stripping of gaseous envelopes appears likely. Hence, it may be difficult to reconcile the observation that many low mass Kepler planets have H/He envelopes with an in situ formation scenario that involves giant impacts after dispersal of the gas disc.


title: Gravitational waves from the quasi-circular inspiral of compact binaries in Einstein-aether theory

abstract: We study gravitational waves emitted by a binary system of non-spinning bodies in a quasi-circular inspiral within the framework of Einstein-aether theory. In particular, we compute explicitly and analytically the expressions for the time-domain and frequency-domain waveforms, gravitational wave polarizations, and response functions for both ground- and space-based detectors in the post-Newtonian approximation. We find that, when going beyond leading-order in the post-Newtonian approximation, the non-Einsteinian polarization modes contain terms that depend on both the first and the second harmonics of the orbital phase. We also calculate analytically the corresponding parameterized post-Einsteinian parameters, generalizing the existing framework to allow for different propagation speeds among scalar, vector and tensor modes, without assumptions about the magnitude of its coupling parameters, and meanwhile allowing the binary system to have relative motions with respect to the aether field. Such results allow for the easy construction of Einstein-aether templates that could be used in Bayesian tests of General Relativity in the future.


title: Reaching a Consensus on Random Networks: The Power of Few

abstract: A community of $n$ individuals splits into two camps, Red and Blue. The individuals are connected by a social network, which influences their colors. Everyday, each person changes his/her color according to the majority among his/her neighbors. Red (Blue) wins if everyone in the community becomes Red (Blue) at some point. We study this process when the underlying network is the random Erdos-Renyi graph $G(n, p)$. With a balanced initial state ($n/2$ person in each camp), it is clear that each color wins with the same probability. Our study reveals that for any constants $p$ and $\varepsilon$, there is a constant $C$ such that if one camp has $n/2 +C$ individuals, then it wins with probability at least $1 - \varepsilon$. The surprising key fact here is that $C$ does not depend on $n$, the population of the community. When $p=1/2$ and $\varepsilon =.1$, one can set $C$ as small as 6. If the aim of the process is to choose a candidate, then this means it takes only $6$ "defectors" to win an election unanimously with overwhelming odd.


title: Manipulation Trajectory Optimization with Online Grasp Synthesis and Selection

abstract: In robot manipulation, planning the motion of a robot manipulator to grasp an object is a fundamental problem. A manipulation planner needs to generate a trajectory of the manipulator arm to avoid obstacles in the environment and plan an end-effector pose for grasping. While trajectory planning and grasp planning are often tackled separately, how to efficiently integrate the two planning problems remains a challenge. In this work, we present a novel method for joint motion and grasp planning. Our method integrates manipulation trajectory optimization with online grasp synthesis and selection, where we apply online learning techniques to select goal configurations for grasping, and introduce a new grasp synthesis algorithm to generate grasps online. We evaluate our planning approach and demonstrate that our method generates robust and efficient motion plans for grasping in cluttered scenes. Our video can be found at https://www.youtube.com/watch?v=LIcACf8YkGU .


title: Extracting galaxy merger timescales I: Tracking haloes with WhereWolf and spinning orbits with OrbWeaver

abstract: Hierarchical models of structure formation predict that dark matter halo assembly histories are characterised by episodic mergers and interactions with other haloes. An accurate description of this process will provide insights into the dynamical evolution of haloes and the galaxies that reside in them. Using large cosmological N-body simulations, we characterise halo orbits to study the interactions between substructure haloes and their hosts, and how different evolutionary histories map to different classes of orbits. We use two new software tools - WhereWolf, which uses halo group catalogues and merger trees to ensure that haloes are tracked accurately in dense environments, and OrbWeaver, which quantifies each halo's orbital parameters. We demonstrate how WhereWolf improves the accuracy of halo merger trees, and we use OrbWeaver to quantify orbits of haloes. We assess how well analytical prescriptions for the merger timescale from the literature compare to measured merger timescales from our simulations and find that existing prescriptions perform well, provided the ratio of substructure-to-host mass is not too small. In the limit of small substructure-to-host mass ratio, we find that the prescriptions can overestimate the merger timescales substantially, such that haloes are predicted to survive well beyond the end of the simulation. This work highlights the need for a revised analytical prescription for the merger timescale that more accurately accounts for processes such as catastrophic tidal disruption.


title: Titchmarsh-Weyl formula for the spectral density of a class of Jacobi matrices in the critical case

abstract: We consider a class of Jacobi matrices with unbounded entries in the so called critical (double root, Jordan box) case. We prove a formula for the spectral density of the matrix which relates its spectral density to the asymptotics of orthogonal polynomials associated with the matrix.


title: A Measurement of Social Capital in an Open Source Software Project

abstract: The paper provides an understanding of social capital in organizations that are open membership multi-agent systems with an emphasis in our formulation on the dynamic network of social interaction that, in part, elucidate evolving structures and impromptu topologies of networks. This paper, therefore, models an open source project as an organizational network. It provides definitions of social capital for this organizational network and formulation of the mechanism to optimize the social capital for achieving its goal that is optimized productivity. A case study of an open source Apache-Hadoop project is considered and empirically evaluated. An analysis of how social capital can be created within this type of organizations and driven to a measurement for its value is provided. Finally, a verification on whether the social capital of the organizational network is proportional towards optimizing their productivity is considered.


title: Copositivity for 3rd order symmetric tensors and applications

abstract: The strict opositivity of 4th order symmetric tensor may apply to detect vacuum stability of general scalar potential. For finding analytical expressions of (strict) opositivity of 4th order symmetric tensor, we may reduce its order to 3rd order to better deal with it. So, it is provided that several analytically sufficient conditions for the copositivity of 3th order 2 dimensional (3 dimensional) symmetric tensors. Subsequently, applying these conclusions to 4th order tensors, the analytically sufficient conditions of copositivity are proved for 4th order 2 dimensional and 3 dimensional symmetric tensors. Finally, we apply these results to present analytical vacuum stability conditions for vacuum stability for $\mathbb{Z}_3$ scalar dark matter.


title: Reliability studies of electronic components for the operation at cryogenic temperature

abstract: Cold electronics is a key technology in many areas of science and technology including space exploration programs and particle physics. A major experiment with a very large number of analog and digital electronics signal processing channels to be operated at cryogenic temperatures is the next-generation neutrino experiment, the Deep Underground Neutrino Experiment (DUNE). The DUNE detector uses liquid Argon at 87K as a target material for neutrinos, and as a medium to track charged particles resulting from interactions in the detector volume. The DUNE electronics [1] consists of custom-designed ASIC (Application Specific Integrated Circuits) chips based on low power 180 nm-CMOS technology. The main risk for this technology is that the electronics components will be immersed in liquid argon for many years (20-30 years) without access. Reliability issues of ASICs may arise from thermal stress, packaging, and manufacturing-related defects: if undetected those could lead to long-term reliability and performance problems. The scope of this paper is to explore non-destructive evaluation techniques for their potential use in a comprehensive quality control process during prototyping, testing and commissioning of the DUNE cold electronics system. Specifically, we have used the Scanning Acoustic Microscopy and X-ray tomography to study permanent structural changes in the ASIC chips associated with thermal cycling between the room and cryogenic temperatures.


title: Dynamic Control of a Fiber Manufacturing Process using Deep Reinforcement Learning

abstract: This paper presents a model-free deep reinforcement learning (DRL) approach for controlling a fiber drawing system. The custom DRL-based control system predictively regulates fiber diameter and produces a fiber with a desired, constant or non-constant, diameter trajectory, i.e. diameter variation along the fiber length. Physical models of the system are not used. The system was trained and tested on a compact fiber drawing system, which has non-linear delayed dynamics and stochastic behaviors. For a reference trajectory with random step changes, after 1 hour of training, the DRL controller showed the same root mean squared error (RMSE) as an optimized PI controller; after 3 hours of training, it achieved the performance of a quadratic dynamic matrix controller (QDMC). While the PI feedback controller showed 3.5 seconds of time lag in a step response, the DRL controller showed less than a second of time lag. Controller performance tests on trajectories not used in the training process are conducted; for a sine sweep reference trajectory, the DRL controller maintained an RMSE under 40 um up to a frequency of 45 mHz, compared to 25 mHz for QDMC.


title: Training Modern Deep Neural Networks for Memory-Fault Robustness

abstract: Because deep neural networks (DNNs) rely on a large number of parameters and computations, their implementation in energy-constrained systems is challenging. In this paper, we investigate the solution of reducing the supply voltage of the memories used in the system, which results in bit-cell faults. We explore the robustness of state-of-the-art DNN architectures towards such defects and propose a regularizer meant to mitigate their effects on accuracy. Our experiments clearly demonstrate the interest of operating the system in a faulty regime to save energy without reducing accuracy.


title: On sequences associated to the invariant theory of rank two simple Lie algebras

abstract: We study two families of sequences, listed in the On-Line Encyclopedia of Integer Sequences (OEIS), which are associated to invariant theory of Lie algebras. For the first family, we prove combinatorially that the sequences A059710 and A108307 are related by a binomial transform. Based on this, we present two independent proofs of a recurrence equation for A059710, which was conjectured by Mihailovs. Besides, we also give a direct proof of Mihailovs' conjecture by the method of algebraic residues. As a consequence, closed formulae for the generating function of sequence A059710 are obtained in terms of classical Gaussian hypergeometric functions. Moreover, we show that sequences in the second family are also related by binomial transforms.


title: Convexification for a 3D inverse scattering problem with the moving point source

abstract: For the first time, we develop in this paper the globally convergent convexification numerical method for a Coefficient Inverse Problem for the 3D Helmholtz equation for the case when the backscattering data are generated by a point source running along an interval of a straight line and the wavenumber is fixed. Thus, by varying the wavenumber, one can reconstruct the dielectric constant depending not only on spatial variables but the wavenumber (i.e. frequency) as well. Our approach relies on a new derivation of a boundary value problem for a system of coupled quasilinear elliptic partial differential equations. This is done via an application of a special truncated Fourier-like method. First, we prove the Lipschitz stability estimate for this problem via a Carleman estimate. Next, using the Carleman Weight Function generated by that estimate, we construct a globally strictly convex cost functional and prove the global convergence to the exact solution of the gradient projection method. Finally, our theoretical finding is verified via several numerical tests with computationally simulated data. These tests demonstrate that we can accurately recover all three important components of targets of interest: locations, shapes and dielectric constants. In particular, large target/background contrasts in dielectric constants (up to 10:1) can be accurately calculated.


title: Scalable sim-to-real transfer of soft robot designs

abstract: The manual design of soft robots and their controllers is notoriously challenging, but it could be augmented---or, in some cases, entirely replaced---by automated design tools. Machine learning algorithms can automatically propose, test, and refine designs in simulation, and the most promising ones can then be manufactured in reality (sim2real). However, it is currently not known how to guarantee that behavior generated in simulation can be preserved when deployed in reality. Although many previous studies have devised training protocols that facilitate sim2real transfer of control polices, little to no work has investigated the simulation-reality gap as a function of morphology. This is due in part to an overall lack of tools capable of systematically designing and rapidly manufacturing robots. Here we introduce a low cost, open source, and modular soft robot design and construction kit, and use it to simulate, fabricate, and measure the simulation-reality gap of minimally complex yet soft, locomoting machines. We prove the scalability of this approach by transferring an order of magnitude more robot designs from simulation to reality than any other method. The kit and its instructions can be found here: https://github.com/skriegman/sim2real4designs


title: Invert and Defend: Model-based Approximate Inversion of Generative Adversarial Networks for Secure Inference

abstract: Inferring the latent variable generating a given test sample is a challenging problem in Generative Adversarial Networks (GANs). In this paper, we propose InvGAN - a novel framework for solving the inference problem in GANs, which involves training an encoder network capable of inverting a pre-trained generator network without access to any training data. Under mild assumptions, we theoretically show that using InvGAN, we can approximately invert the generations of any latent code of a trained GAN model. Furthermore, we empirically demonstrate the superiority of our inference scheme by quantitative and qualitative comparisons with other methods that perform a similar task. We also show the effectiveness of our framework in the problem of adversarial defenses where InvGAN can successfully be used as a projection-based defense mechanism. Additionally, we show how InvGAN can be used to implement reparameterization white-box attacks on projection-based defense mechanisms. Experimental validation on several benchmark datasets demonstrate the efficacy of our method in achieving improved performance on several white-box and black-box attacks. Our code is available at https://github.com/yogeshbalaji/InvGAN.


title: Nonlocal Poincar\'{e} Inequalities for Integral Operators with Integrable Nonhomogeneous Kernels

abstract: The paper provides two versions of nonlocal Poincar\'e-type inequalities for integral operators with a convolution-type structure and functions satisfying a zero-Dirichlet like condition. The inequalities extend existing results to a large class of nonhomogeneous kernels with supports that can vary discontinuously and need not contain a common set throughout the domain. The measure of the supports may even vanish allowing the zero-Dirichlet condition to be imposed on only a lower-dimensional manifold, with or without boundary. The conditions may be imposed on sets with co-dimension larger than one or even at just a single point. This appears to currently be the first such results in a nonlocal setting with integrable kernels. The arguments used are direct, and examples are provided demonstrating the explicit dependence of bounds for the Poincar\'{e} constant upon structural parameters of the kernel and domain.


title: A Domain Adaptive Density Clustering Algorithm for Data with Varying Density Distribution

abstract: As one type of efficient unsupervised learning methods, clustering algorithms have been widely used in data mining and knowledge discovery with noticeable advantages. However, clustering algorithms based on density peak have limited clustering effect on data with varying density distribution (VDD), equilibrium distribution (ED), and multiple domain-density maximums (MDDM), leading to the problems of sparse cluster loss and cluster fragmentation. To address these problems, we propose a Domain-Adaptive Density Clustering (DADC) algorithm, which consists of three steps: domain-adaptive density measurement, cluster center self-identification, and cluster self-ensemble. For data with VDD features, clusters in sparse regions are often neglected by using uniform density peak thresholds, which results in the loss of sparse clusters. We define a domain-adaptive density measurement method based on K-Nearest Neighbors (KNN) to adaptively detect the density peaks of different density regions. We treat each data point and its KNN neighborhood as a subgroup to better reflect its density distribution in a domain view. In addition, for data with ED or MDDM features, a large number of density peaks with similar values can be identified, which results in cluster fragmentation. We propose a cluster center self-identification and cluster self-ensemble method to automatically extract the initial cluster centers and merge the fragmented clusters. Experimental results demonstrate that compared with other comparative algorithms, the proposed DADC algorithm can obtain more reasonable clustering results on data with VDD, ED and MDDM features. Benefitting from a few parameter requirements and non-iterative nature, DADC achieves low computational complexity and is suitable for large-scale data clustering.


title: Solution Curve for Linear Control Systems on Lie Groups

abstract: The purpose of this paper is to describe explicitly the solution for linear control systems on Lie groups. In case of linear control systems with inner derivations, the solution is given basically by the product of the exponential of the associated invariant system and the exponential of the associated invariant drift field. We present the solutions in low dimensional cases and apply the results to obtain some controllability results.


title: Fundamental limits to attractive and repulsive Casimir--Polder forces

abstract: We derive upper and lower bounds on the Casimir--Polder force between an anisotropic dipolar body and a macroscopic body separated by vacuum via algebraic properties of Maxwell's equations. These bounds require only a coarse characterization of the system---the material composition of the macroscopic object, the polarizability of the dipole, and any convenient partition between the two objects---to encompass all structuring possibilities. We find that the attractive Casimir--Polder force between a polarizable dipole and a uniform planar semi-infinite bulk medium always comes within 10% of the lower bound, implying that nanostructuring is of limited use for increasing attraction. In contrast, the possibility of repulsion is observed even for isotropic dipoles, and is routinely found to be several orders of magnitude larger than any known design, including recently predicted geometries involving conductors with sharp edges. Our results have ramifications for the design of surfaces to trap, suspend, or adsorb ultracold gases.


title: Modified black hole solution with a background Kalb-Ramond field

abstract: We study the gravitation effects on a static and spherically symmetric spacetime due to the vacuum expectation value (VEV) of a Kalb-Ramond field. The Kalb-Ramond VEV is a background tensor field which produces a local Lorentz symmetry breaking (LSB) of spacetime. Considering a non-minimal coupling between the Kalb-Ramond (VEV) and the Ricci tensor, we obtain an exact parameter-dependent power-law modified black hole. For a particular choice of the LSB parameter, the Lorentz violation produces a solution similar to the Reissner-Nordstrom, despite the absence of charge. The near-horizon geometry is modified by including a new inner horizon and shifting the Schwarzschild horizon. Asymptotically, the usual Minkowski spacetime with a background tensor field is recovered. The vacuum configurations are studied considering the energy conditions and the Lorentz violating source properties. By means of the mercury perihelion test, an upper bound to the local Lorentz violation (LV) is obtained, and its corresponding effects on the black hole temperature is investigated.


title: Financial ratios and stock returns reappraised through a topological data analysis lens

abstract: Firm financials are well established as return predictors, being the inspiration for a large set of anomalies in the asset pricing literature. Employing topological data analysis we revisit the question of association between seven of the most commonly studied financial ratios and stock returns. Specifically the TDA Ball Mapper algorithm is applied to visualise the point cloud of financial ratios as an abstract two-dimensional graph readily allowing for identification of interdependencies between factors. These relationships are seldom monotonic, opportunities for investors to profitably exploit this knowledge provided by TDA abound. Clear potential offered by the tools of TDA to shed new light on asset pricing models is demonstrated. Scope for benefit is limited only by the availability of information to the analyst.


title: CoverNet: Multimodal Behavior Prediction using Trajectory Sets

abstract: We present CoverNet, a new method for multimodal, probabilistic trajectory prediction for urban driving. Previous work has employed a variety of methods, including multimodal regression, occupancy maps, and 1-step stochastic policies. We instead frame the trajectory prediction problem as classification over a diverse set of trajectories. The size of this set remains manageable due to the limited number of distinct actions that can be taken over a reasonable prediction horizon. We structure the trajectory set to a) ensure a desired level of coverage of the state space, and b) eliminate physically impossible trajectories. By dynamically generating trajectory sets based on the agent's current state, we can further improve our method's efficiency. We demonstrate our approach on public, real-world self-driving datasets, and show that it outperforms state-of-the-art methods.


title: Predator-Prey Interaction Model with Hunting Cooperation among Predators and Allee Effect in Prey

abstract: This paper investigates a dynamical predator-prey interaction model that incorporates: (a) hunting cooperation among predators; (b) Allee effect in prey. We show all possible boundary and interior solutions. In order to analyze the stability of the solution, we make use of the Jacobian matrix and the resultant characteristic polynomial. Particularly, the sign of the eigenvalue is used to determine the stability of a solution. We then provide proof for stability of the interior solution. Finally, we verify our results numerically in MATLAB by plotting: (1) predator-prey intersection graphs; (2) prey-predator vs hunting cooperation graphs; (3) initial condition trajectory for equilibrium solution. It is interesting to notice that hunting cooperation can switch the stability of coexistence equilibrium solutions. Through numerical simulations, it was verified that increasing the hunting cooperation could lead to the extinction of both prey and predator population for alpha greater than 0.96, given our choice of parameters.


title: Non-reciprocal Cavity Polariton with Atoms Strongly Coupled to Optical Cavity

abstract: Breaking the time-reversal symmetry of light is of great importance for fundamental physics and has attracted increasing interest in the study of non-reciprocal photonic devices. Here, we experimentally demonstrate a chiral cavity QED system with multiple atoms strongly coupled to a Fabry-Perot cavity. By polarizing the internal quantum state of the atoms, the time-reversal symmetry of the atom-cavity interaction is broken. The strongly coupled atom-cavity system can be described by non-reciprocal quasiparticles, i.e., the cavity polariton. When it works in the linear regime, the inherent nonreciprocity makes the system work as a single-photon-level optical isolator. Benefiting from the collective enhancement of multiple atoms, an isolation ratio exceeding 30~dB on the single-quanta level (~0.1 photon on average) is achieved. The validity of the non-reciprocal device under zero magnetic field and the reconfigurability of the isolation direction are also experimentally demonstrated. Moreover, when the cavity polariton works in the nonlinear regime, the quantum interference between polaritons with weak anharmonicity induces non-reciprocal nonclassical statistics of cavity transmission from coherent probe light.


title: Learning a Representation with the Block-Diagonal Structure for Pattern Classification

abstract: Sparse-representation-based classification (SRC) has been widely studied and developed for various practical signal classification applications. However, the performance of a SRC-based method is degraded when both the training and test data are corrupted. To counteract this problem, we propose an approach that learns Representation with Block-Diagonal Structure (RBDS) for robust image recognition. To be more specific, we first introduce a regularization term that captures the block-diagonal structure of the target representation matrix of the training data. The resulting problem is then solved by an optimizer. Last, based on the learned representation, a simple yet effective linear classifier is used for the classification task. The experimental results obtained on several benchmarking datasets demonstrate the efficacy of the proposed RBDS method.


title: Axisymmetric diffeomorphisms and ideal fluids on Riemannian 3-manifolds

abstract: We study the Riemannian geometry of 3D axisymmetric ideal fluids. We prove that the $L^2$ exponential map on the group of volume-preserving diffeomorphisms of a $3$-manifold is Fredholm along axisymmetric flows with sufficiently small swirl. Along the way, we define the notions of axisymmetric and swirl-free diffeomorphisms of any manifold with suitable symmetries and show that such diffeomorphisms form a totally geodesic submanifold of infinite $L^2$ diameter inside the space of volume-preserving diffeomorphisms whose diameter is known to be finite. As examples we derive the axisymmetric Euler equations on $3$-manifolds equipped with each of Thurston's eight model geometries.


title: New Transceiver Designs for Interleaved Frequency Division Multiple Access

abstract: This paper puts forth a class of new transceiver designs for interleaved frequency division multiple access (IFDMA) systems. These transceivers are significantly less complex than conventional IFDMA transceiver. The simple new designs are founded on a key observation that multiplexing and demultiplexing of IFDMA data streams of different sizes are coincident with the IFFTs and FFTs of different sizes embedded within the Cooley-Tukey recursive FFT decomposition scheme. For flexible resource allocation, this paper puts forth a new IFDMA resource allocation framework called Multi-IFDMA, in which a user can be allocated multiple IFDMA streams. Our new transceivers are unified designs in that they can be used in conventional IFDMA as well as multi-IFDMA systems. Two other well-known multiple-access schemes are localized FDMA (LFDMA) and orthogonal FDMA (OFDMA). In terms of flexibility in resource allocation, Multi-IFDMA, LFDMA, and OFDMA are on an equal footing. With our new transceiver designs, however, IFDMA has the following advantages (besides other known advantages not due to our new transceiver designs): 1) IFDMA/Multi-IFDMA transceivers are significantly less complex than LFDMA transceivers; in addition, IFDMA/Multi-IFDMA has better Peak-to-Average Power Ratio (PAPR) than LFDMA; 2) IFDMA/Multi-IFDMA transceivers and OFDMA transceivers are comparable in complexity; but IFDMA/Multi-IFDMA has significantly better PAPR than OFDMA.


title: Subexponential LPs Approximate Max-Cut

abstract: We show that for every $\varepsilon > 0$, the degree-$n^\varepsilon$ Sherali-Adams linear program (with $\exp(\tilde{O}(n^\varepsilon))$ variables and constraints) approximates the maximum cut problem within a factor of $(\frac{1}{2}+\varepsilon')$, for some $\varepsilon'(\varepsilon) > 0$. Our result provides a surprising converse to known lower bounds against all linear programming relaxations of Max-Cut, and hence resolves the extension complexity of approximate Max-Cut for approximation factors close to $\frac{1}{2}$ (up to the function $\varepsilon'(\varepsilon)$). Previously, only semidefinite programs and spectral methods were known to yield approximation factors better than $\frac 12$ for Max-Cut in time $2^{o(n)}$. We also show that constant-degree Sherali-Adams linear programs (with $\text{poly}(n)$ variables and constraints) can solve Max-Cut with approximation factor close to $1$ on graphs of small threshold rank: this is the first connection of which we are aware between threshold rank and linear programming-based algorithms. Our results separate the power of Sherali-Adams versus Lov\'asz-Schrijver hierarchies for approximating Max-Cut, since it is known that $(\frac{1}{2}+\varepsilon)$ approximation of Max Cut requires $\Omega_\varepsilon (n)$ rounds in the Lov\'asz-Schrijver hierarchy. We also provide a subexponential time approximation for Khot's Unique Games problem: we show that for every $\varepsilon > 0$ the degree-$(n^\varepsilon \log q)$ Sherali-Adams linear program distinguishes instances of Unique Games of value $\geq 1-\varepsilon'$ from instances of value $\leq \varepsilon'$, for some $\varepsilon'( \varepsilon) >0$, where $q$ is the alphabet size. Such guarantees are qualitatively similar to those of previous subexponential-time algorithms for Unique Games but our algorithm does not rely on semidefinite programming or subspace enumeration techniques.


title: Dynamical System Inspired Adaptive Time Stepping Controller for Residual Network Families

abstract: The correspondence between residual networks and dynamical systems motivates researchers to unravel the physics of ResNets with well-developed tools in numeral methods of ODE systems. The Runge-Kutta-Fehlberg method is an adaptive time stepping that renders a good trade-off between the stability and efficiency. Can we also have an adaptive time stepping for ResNets to ensure both stability and performance? In this study, we analyze the effects of time stepping on the Euler method and ResNets. We establish a stability condition for ResNets with step sizes and weight parameters, and point out the effects of step sizes on the stability and performance. Inspired by our analyses, we develop an adaptive time stepping controller that is dependent on the parameters of the current step, and aware of previous steps. The controller is jointly optimized with the network training so that variable step sizes and evolution time can be adaptively adjusted. We conduct experiments on ImageNet and CIFAR to demonstrate the effectiveness. It is shown that our proposed method is able to improve both stability and accuracy without introducing additional overhead in inference phase.


title: Liouville type theorems for minimal graphs over manifolds

abstract: Let $\Sigma$ be a complete Riemannian manifold with the volume doubling property and the uniform Neumann-Poincar$\mathrm{\acute{e}}$ inequality. We show that any positive minimal graphic function on $\Sigma$ is a constant.


title: A Semi-Markov Model for the Control of Thermostatically Controlled Load

abstract: Considering the potential of thermostatically controlled loads (TCLs) to provide flexibility in demand response or load control, a semi-Markov model (SMM) for the ON/OFF controlled TCL is developed in this paper. This model makes full use of the adjustment flexibility of TCLs when the control period is long and maintains the diversity of switch state in the cluster. This model also can satisfy user comfort and protect user privacy. Then, this paper adopts the cyber-physical system (CPS) to realize the coupling of the discrete control process and the continuous physical process. Finally, the proposed model is applied to the coordination of large-scale heterogenous air-conditioners (ACs) based on the equivalent thermal parameters (ETP) model. Simulation results verify that under the proposed approach, the power of TCLs cluster can track the control signal accurately, with both user comfort and diversity of TCL cluster's operation states guaranteed.


title: Some sum-product type estimates for two-variables over prime fields

abstract: In this paper, we use a recent method given by Rudnev, Shakan, and Shkredov (2018) to improve results on sum-product type problems due to Pham and Mojarrad (2018).


title: Neural Integration of Continuous Dynamics

abstract: Neural dynamical systems are dynamical systems that are described at least in part by neural networks. The class of continuous-time neural dynamical systems must, however, be numerically integrated for simulation and learning. Here, we present a compact neural circuit for two common numerical integrators: the explicit fixed-step Runge-Kutta method of any order and the semi-implicit/predictor-corrector Adams-Bashforth-Moulton method. Modeled as constant-sized recurrent networks embedding a continuous neural differential equation, they achieve fully neural temporal output. Using the polynomial class of dynamical systems, we demonstrate the equivalence of neural and numerical integration.

categories: cs.LG cs.NA math.DS math.NA nlin.CD stat.ML


title: Multi-mode CV-QKD with Non-Gaussian Operations

abstract: Non-Gaussian operations have been studied intensively in recent years due to their ability to increase the secret key rate for certain CV-QKD protocols. However, most previous studies on such protocols are carried out in a single-mode setting, even though in reality any quantum state contains multi-mode components in frequency space. In this work we investigate the use of non-Gaussian operations in a multi-mode CV-QKD system. Our main finding is that, contrary to single-mode CV-QKD systems, in generic multi-mode CV-QKD systems it is possible to use non-Gaussian operations to increase the optimized secret key rate. More specifically, we find that at losses of order 30dB, which represents a distance of order 160km and the effective maximum distance for CV-QKD, the key rate for multi-mode non-Gaussian operations can be orders of magnitude higher than single-mode operations. Our results are important for real-world CV-QKD systems especially those dependent on quantum error correction - a process that requires non-Gaussian effects.


title: Weighted Laplacian and Its Theoretical Applications

abstract: In this paper, we develop a novel weighted Laplacian method, which is partially inspired by the theory of graph Laplacian, to study recent popular graph problems, such as multilevel graph partitioning and balanced minimum cut problem, in a more convenient manner. Since the weighted Laplacian strategy inherits the virtues of spectral methods, graph algorithms designed using weighted Laplacian will necessarily possess more robust theoretical guarantees for algorithmic performances, comparing with those existing algorithms that are heuristically proposed. In order to illustrate its powerful utility both in theory and in practice, we also present two effective applications of our weighted Laplacian method to multilevel graph partitioning and balanced minimum cut problem, respectively. By means of variational methods and theory of partial differential equations (PDEs), we have established the equivalence relations among the weighted cut problem, balanced minimum cut problem and the initial clustering problem that arises in the middle stage of graph partitioning algorithms under a multilevel structure. These equivalence relations can indeed provide solid theoretical support for algorithms based on our proposed weighted Laplacian strategy. Moreover, from the perspective of the application to the balanced minimum cut problem, weighted Laplacian can make it possible for research of numerical solutions of PDEs to be a powerful tool for the algorithmic study of graph problems. Experimental results also indicate that the algorithm embedded with our strategy indeed outperforms other existing graph algorithms, especially in terms of accuracy, thus verifying the efficacy of the proposed weighted Laplacian.


title: Design and Evaluation of a Multi-Domain Trojan Detection Method on Deep Neural Networks

abstract: This work corroborates a run-time Trojan detection method exploiting STRong Intentional Perturbation of inputs, is a multi-domain Trojan detection defence across Vision, Text and Audio domains---thus termed as STRIP-ViTA. Specifically, STRIP-ViTA is the first confirmed Trojan detection method that is demonstratively independent of both the task domain and model architectures. We have extensively evaluated the performance of STRIP-ViTA over: i) CIFAR10 and GTSRB datasets using 2D CNNs, and a public third party Trojaned model for vision tasks; ii) IMDB and consumer complaint datasets using both LSTM and 1D CNNs for text tasks; and speech command dataset using both 1D CNNs and 2D CNNs for audio tasks. Experimental results based on 28 tested Trojaned models demonstrate that STRIP-ViTA performs well across all nine architectures and five datasets. In general, STRIP-ViTA can effectively detect Trojan inputs with small false acceptance rate (FAR) with an acceptable preset false rejection rate (FRR). In particular, for vision tasks, we can always achieve a 0% FRR and FAR. By setting FRR to be 3%, average FAR of 1.1% and 3.55% are achieved for text and audio tasks, respectively. Moreover, we have evaluated and shown the effectiveness of STRIP-ViTA against a number of advanced backdoor attacks whilst other state-of-the-art methods lose effectiveness in front of one or all of these advanced backdoor attacks.


title: DTF: Deep Tensor Factorization for Predicting Anticancer Drug Synergy

abstract: Motivation: Combination therapies have been widely used to treat cancers. However, it is cost- and time-consuming to experimentally screen synergistic drug pairs due to the enormous number of possible drug combinations. Thus, computational methods have become an important way to predict and prioritize synergistic drug pairs. Results: We proposed a Deep Tensor Factorization (DTF) model, which integrated a tensor factorization method and a deep neural network (DNN), to predict drug synergy. The former extracts latent features from drug synergy information while the latter constructs a binary classifier to predict the drug synergy status. Compared to the tensor-based method, the DTF model performed better in predicting drug synergy. The area under the precision-recall curve (PR AUC) was 0.57 for DTF and 0.24 for the tensor method. We also compared the DTF model with DeepSynergy and logistic regression models and found that the DTF outperformed the logistic regression model and achieved almost the same performance as DeepSynergy using several typical metrics for the classification task. Applying the DTF model to predict missing entries in our drug-cell line tensor, we identified novel synergistic drug combinations for 10 cell lines from the 5 cancer types. A literature survey showed that some of these predicted drug synergies have been identified in vivo or in vitro. Thus, the DTF model could be valuable in silico tool for prioritizing novel synergistic drug combinations.


title: Robust Q-tunable topological induced transparency in metasurface

abstract: In this study, we demonstrate first, to the best of our knowledge, robust and dynamically polarization-controlled tunable-high-Q PIT in designed nanostructures metasurface, whose sharp resonance is guaranteed by design and protected against large geometrical imperfections. By employing the explicit analysis of near-field characteristic in the reciprocal-space based on the momentum matching, and the far-field radiation features with point-scattering approach in real-space sparked from Huygens's principles, the physics of interference resonance in the spectra for plane-wave optical transmission and reflection of the metasurface is theoretically and thoroughly investigated. The experimental results verify the theory prediction. The distinctive polarization-selective and Q-tunable PIT shows robust features to performance degradations in traditional PIT system caused by inadvertent fabrication flaws or geometry asymmetry-variations, which paves way for the development of reconfigurable and flexible metasurface and, additionally, opens new avenues in robust and multifunctional nanophotonics device design and applications.


title: Surface-plasmon-enhanced near-field radiative heat transfer between planar surfaces with a thin-film plasmonic coupler

abstract: In last decade, there have been enormous efforts to experimentally show the near-field enhancement of radiative heat transfer between planar structures. Several recent experiments also have striven to achieve further enhanced heat transfer with the excitation of coupled surface polaritons by introducing nanostructures on both emitter and the receiver; however, these symmetric structures are hardly employed in real-world applications. Here, we demonstrate substantially increased near-field radiative heat transfer between asymmetric structures (i.e., doped Si and SiO$_2$) by using a thin Ti film as a plasmonic coupler. The measured near-field enhancement at vacuum gap of 380 nm is found to be 3.5 times greater than that for the case without the coupler. The enhancement mechanism is thoroughly elucidated for both polarizations and a dimensionless parameter, which can quantify the coupling strength of the surface polaritons at vacuum, is suggested. As a thin film can be readily used in many engineering applications, this study will facilitate the development of the high-performance engineering devices exploiting the near-field thermal radiation.


title: Deep graph embedding for prioritizing synergistic anticancer drug combinations

abstract: Drug combinations are frequently used for the treatment of cancer patients in order to increase efficacy, decrease adverse side effects, or overcome drug resistance. Given the enormous number of drug combinations, it is cost- and time-consuming to screen all possible drug pairs experimentally. Currently, it has not been fully explored to integrate multiple networks to predict synergistic drug combinations using recently developed deep learning technologies. In this study, we proposed a Graph Convolutional Network (GCN) model to predict synergistic drug combinations in particular cancer cell lines. Specifically, the GCN method used a convolutional neural network model to do heterogeneous graph embedding, and thus solved a link prediction task. The graph in this study was a multimodal graph, which was constructed by integrating the drug-drug combination, drug-protein interaction, and protein-protein interaction networks. We found that the GCN model was able to correctly predict cell line-specific synergistic drug combinations from a large heterogonous network. The majority (30) of the 39 cell line-specific models show an area under the receiver operational characteristic curve (AUC) larger than 0.80, resulting in a mean AUC of 0.84. Moreover, we conducted an in-depth literature survey to investigate the top predicted drug combinations in specific cancer cell lines and found that many of them have been found to show synergistic antitumor activity against the same or other cancers in vitro or in vivo. Taken together, the results indicate that our study provides a promising way to better predict and optimize synergistic drug pairs in silico.


title: PlantDoc: A Dataset for Visual Plant Disease Detection

abstract: India loses 35% of the annual crop yield due to plant diseases. Early detection of plant diseases remains difficult due to the lack of lab infrastructure and expertise. In this paper, we explore the possibility of computer vision approaches for scalable and early plant disease detection. The lack of availability of sufficiently large-scale non-lab data set remains a major challenge for enabling vision based plant disease detection. Against this background, we present PlantDoc: a dataset for visual plant disease detection. Our dataset contains 2,598 data points in total across 13 plant species and up to 17 classes of diseases, involving approximately 300 human hours of effort in annotating internet scraped images. To show the efficacy of our dataset, we learn 3 models for the task of plant disease classification. Our results show that modelling using our dataset can increase the classification accuracy by up to 31%. We believe that our dataset can help reduce the entry barrier of computer vision techniques in plant disease detection.


title: The circumstellar environment around the embedded protostar EC 53

abstract: EC53 is an embedded protostar with quasi-periodic emission in the near-IR and sub-mm. We use ALMA high-resolution observations of continuum and molecular line emission to describe the circumstellar environment of EC 53. The continuum image reveals a disk with a flux that suggests a mass of 0.075 Msun, much less than the estimated mass in the envelope, and an in-band spectral index that indicates grain growth to centimeter sizes. Molecular lines trace the outflow cavity walls, infalling and rotating envelope, and/or the Keplerian disk. The rotation profile of the C17O 3-2 line emission cannot isolate the Keplerian motion clearly although the lower limit of the protostellar mass can be calculated as 0.3 +- 0.1 Msun if the Keplerian motion is adopted. The weak CH3OH emission, which is anti-correlated with the HCO+ 4-3 line emission, indicates that the water snow line is more extended than what expected from the current luminosity, attesting to bygone outburst events. The extended snow line may persist for longer at the disk surface because the lower density increases the freeze-out timescale of methanol and water.


title: Elementary hypergeometric functions, Heun functions, and moments of MKZ operators

abstract: We consider some hypergeometric functions and prove that they are elementary functions. Consequently, the second order moments of Meyer-Konig and Zeller type operators are elementary functions. The higher order moments of these operators are expressed in terms of elementary functions and polylogarithms. Other applications are concerned with the expansion of certain Heun functions in series or finite sums of elementary hypergeometric functions.


title: Opportunities and Outcomes for Postdocs in Canada

abstract: Currently, postdoctoral fellow (PDF) researchers in Canada face challenges due to the precarious nature of their employment and their overall low compensation and benefits coverage. This report presents three themes, written as statements of need, to support an inclusive and thriving PDF community. These themes are the need for better terms of employment and conditions, the need for access to grants by non-permanent research staff, and the need for a sustainable PDF hiring model that considers the outcomes for the PDFs. We make six recommendations: R1. PDFs should be hired and compensated as skilled experts in their areas, not as trainees. R2. Standard PDF hiring practices should be revised to be more inclusive of different life circumstances. - R2.1 Allow PDFs the option of part-time employment. - R2.2 Remove years-since-PhD time limits from PDF jobs. - R2.3 Financially support PDF hires for relocation and visa expenses. R3. CASCA should form a committee to advocate for and provide support to astronomy PDFs in Canada. R4. CASCA should encourage universities to create offices dedicated to their PDFs. R5. PDFs and other PhD-holding term researchers with a host institution should be able to compete for and win grants to self-fund their own research. R6. Astronomy in Canada should hire general-purpose continuing support scientist positions instead of term PDFs to fill project or mission-specific requirements. In short, we ask for prioritization of people over production of papers.


title: Compressing Representations for Embedded Deep Learning

abstract: Despite recent advances in architectures for mobile devices, deep learning computational requirements remains prohibitive for most embedded devices. To address that issue, we envision sharing the computational costs of inference between local devices and the cloud, taking advantage of the compression performed by the first layers of the networks to reduce communication costs. Inference in such distributed setting would allow new applications, but requires balancing a triple trade-off between computation cost, communication bandwidth, and model accuracy. We explore that trade-off by studying the compressibility of representations at different stages of MobileNetV2, showing those results agree with theoretical intuitions about deep learning, and that an optimal splitting layer for network can be found with a simple PCA-based compression scheme.


title: Meta Adaptation using Importance Weighted Demonstrations

abstract: Imitation learning has gained immense popularity because of its high sample-efficiency. However, in real-world scenarios, where the trajectory distribution of most of the tasks dynamically shifts, model fitting on continuously aggregated data alone would be futile. In some cases, the distribution shifts, so much, that it is difficult for an agent to infer the new task. We propose a novel algorithm to generalize on any related task by leveraging prior knowledge on a set of specific tasks, which involves assigning importance weights to each past demonstration. We show experiments where the robot is trained from a diversity of environmental tasks and is also able to adapt to an unseen environment, using few-shot learning. We also developed a prototype robot system to test our approach on the task of visual navigation, and experimental results obtained were able to confirm these suppositions.


title: Real-time Terahertz Wave Channeling via Multifunctional Metagratings: A Sparse Array of All-Graphene Scatterers

abstract: Acquiring full control over a large number of diffraction orders can be strongly attractive in the case of realizing multifunctional devices such as multichannel reflectors. Recently, the concept of metagrating has been introduced which enables obtaining the desired diffraction pattern through a sparse periodic array of engineered scatterers. In this letter, for the first time, a tunable all-graphene multichannel meta-reflector is proposed for operating at terahertz (THz) frequencies. In the supercell level, the designed metagrating is composed of three graphene ribbons of different controllable chemical potentials which can be regarded as a five-channel THz meta-reflector. Several illustrative examples have been presented in which by choosing proper distribution of DC voltages feeding the ribbons, our design can realize different intriguing functionalities such as anomalous reflection, retro-reflection, and three-channel power splitting within a single shared aperture and with high efficiency. This work paves the way toward designing highly-efficient and tunable THz multichannel meta-reflectors with many potential applications in photonics and optoelectronics.


title: Proximality of multidimensional $\mathscr{B}$-free systems

abstract: We characterize proximality of multidimensional $\mathscr{B}$-free systems in the case of number fields and lattices in $\mathbb{Z}^m$, $m\geq2$.


title: The Nernst effect in Corbino geometry

abstract: We study the manifestation of the Nernst effect in the Corbino disk subjected to the normal external magnetic field and to the radial temperature gradient. The Corbino geometry offers a precious opportunity for the direct measurement of the magnetization currents that are masked by kinetic contributions to the Nernst current in the conventional geometry. The magnetization currents, also referred to as the edge currents, are independent on the conductivity of the sample which is why they can be conveniently described within the thermodynamic approach. They can be related to the Landau thermodynamic potential for an infinite system. We demonstrate that the observable manifestation of this, purely thermodynamic, Nernst effect consists in the strong oscillations of the magnetic field measured in the center of the disk as a function of the external field. The oscillations depend on the temperature difference at the edges of the disk. Dirac fermions and 2D electrons with a parabolic spectrum are characterized by oscillations of different phase and frequency. We predict qualitatively different power dependencies of the magnitude of the Nernst signal on the chemical potential for normal and Dirac carriers.


title: Deep learning reconstruction of ultrashort pulses from 2D spatial intensity patterns recorded by an all-in-line system in a single-shot

abstract: We propose a simple all-in-line single-shot scheme for diagnostics of ultrashort laser pulses, consisting of a multi-mode fiber, a nonlinear crystal and a CCD camera. The system records a 2D spatial intensity pattern, from which the pulse shape (amplitude and phase) are recovered, through a fast Deep Learning algorithm. We explore this scheme in simulations and demonstrate the recovery of ultrashort pulses, robustness to noise in measurements and to inaccuracies in the parameters of the system components. Our technique mitigates the need for commonly used iterative optimization reconstruction methods, which are usually slow and hampered by the presence of noise. These features make our concept system advantageous for real time probing of ultrafast processes and noisy conditions. Moreover, this work exemplifies that using deep learning we can unlock new types of systems for pulse recovery.


title: On the Breuil-Schneider conjecture II: Potentially crystalline non-generic case

abstract: This paper improves some results of the author's previous work. We will investigate the case of non-smooth points on an automorphic components and prove Breuil-Schenider conjecture. As a consequence we will see that in case when all the components are automorphic in the potentially crystalline deformation rings, the Breuil-Schenider conjecture holds unconditionally.


title: Path integral solution for a deformed radial Rosen-Morse potential

abstract: An exact path integral treatment of a particle in a deformed radial Rosen-Morse potential is presented. For this problem with the Dirichlet boundary conditions, the Green's function is constructed in a closed form by adding to V_{q}(r) a {\delta}-function perturbation and making its strength infinitely repulsive. A transcendental equation for the energy levels E_{n_{r}} and the wave functions of the bound states can then be deduced.


title: On a compact trace embedding theorem in Musielak-Sobolev spaces

abstract: By a stronger compact boundary embedding theorem in Musielak-Orlicz-Sobolev space developed in the paper, variational method is employed to deal with the nonlinear elliptic equation with the nonlinear Neumann boundary condition in the framework of Musielak-Orlicz-Sobolev space.


title: Many-body calculations for periodic materials via quantum machine learning

abstract: A state-of-the-art method that combines a quantum computational algorithm and machine learning, so-called quantum machine learning, can be a powerful approach for solving quantum many-body problems. However, the research scope in the field was mainly limited to organic molecules and simple lattice models. Here, we propose a workflow of quantum machine learning applications for periodic systems on the basis of an effective model construction from first principles. The band structures of the Hubbard model of graphene with the mean-field approximation are calculated as a benchmark, and the calculated eigenvalues show good agreement with the exact diagonalization results within a few meV by employing the transfer learning technique in quantum machine learning. The results show that the present computational scheme has the potential to solve many-body problems quickly and correctly for periodic systems using a quantum computer.


title: Prethermal memory loss in interacting quantum systems coupled to thermal baths

abstract: We study the relaxation dynamics of an extended Fermi-Hubbard chain with a strong Wannier-Stark potential tilt coupled to a bath. When the system is subjected to dephasing noise, starting from a pure initial state the system's total von Neumann entropy is found to grow monotonously. The scenario becomes rather different when the system is coupled to a thermal bath of finite temperature. Here, for sufficiently large field gradients and initial energies, the entropy peaks in time and almost reaches its largest possible value (corresponding to the maximally mixed state), long before the system relaxes to thermal equilibrium. This entropy peak signals a prethermal memory loss and, relative to the time where it occurs, the system is found to exhibit a simple scaling behavior in space and time. By comparing the system's dynamics to that of a simplified model, the underlying mechanism is found to be related to the localization property of the Wannier-Stark system, which favors dissipative coupling between eigenstates that are close in energy.


title: Approximate path integral solution for a Dirac particle in a deformed Hulth\'en potential

abstract: The problem of a Dirac particle moving in a deformed Hulthen potential is solved in the framework of the path integral formalism. With the help of the Biedenharn transformation, the construction of a closed form for the Green's function of the second-order Dirac equation is done by using a proper approximation to the centrifugal term and the Green's function of the linear Dirac equation is calculated. The energy spectrum for the bound states is obtained from the poles of the Green's function. A Dirac particle in the standard Hulthen potential for q=1 and a Dirac hydrogen-like ion when q = 1 and a that tends to infinity are considered as particular cases.


title: The generalised buoyancy/inertial forces and available energy of axisymmetric compressible stratified vortex motions

abstract: Adiabatic and inviscid axisymmetric perturbations to a stable reference vortex in gradient wind balance are known to experience two kinds of restoring forces: one that is proportional to both the perturbation density and the reference pressure gradient, and one that is purely radial and proportional to the squared angular momentum perturbation. We show that the work required to move a fluid parcel against such forces from its equilibrium to actual position is path-independent and formally equivalent to the available energy accounting for momentum constraints previously constructed by Andrews (2006) and Codoban and Shepherd (2006). Physically, this work represents the energy of the unbalanced part of the vortex and hence a form of eddy energy. It can be partitioned into available acoustic energy, slantwise available potential energy and centrifugal potential energy. We show that the conditions required for these energies to be positive definite correspond to the classical conditions for symmetric stability, valid for both small and finite-amplitude perturbations. The new available energy framework possesses various advantageous features over previous approaches that shed new light on how the thermodynamic and mechanical sinks/sources of energy control the intensification of the balanced and unbalanced parts of a warm core cyclonic vortex. These features should prove useful in the future to re-examine and clarify the links between the different existing paradigms of tropical cyclone (TC) intensification within a single unifying theoretical framework.


title: Iteratively-Refined Interactive 3D Medical Image Segmentation with Multi-Agent Reinforcement Learning

abstract: Existing automatic 3D image segmentation methods usually fail to meet the clinic use. Many studies have explored an interactive strategy to improve the image segmentation performance by iteratively incorporating user hints. However, the dynamic process for successive interactions is largely ignored. We here propose to model the dynamic process of iterative interactive image segmentation as a Markov decision process (MDP) and solve it with reinforcement learning (RL). Unfortunately, it is intractable to use single-agent RL for voxel-wise prediction due to the large exploration space. To reduce the exploration space to a tractable size, we treat each voxel as an agent with a shared voxel-level behavior strategy so that it can be solved with multi-agent reinforcement learning. An additional advantage of this multi-agent model is to capture the dependency among voxels for segmentation task. Meanwhile, to enrich the information of previous segmentations, we reserve the prediction uncertainty in the state space of MDP and derive an adjustment action space leading to a more precise and finer segmentation. In addition, to improve the efficiency of exploration, we design a relative cross-entropy gain-based reward to update the policy in a constrained direction. Experimental results on various medical datasets have shown that our method significantly outperforms existing state-of-the-art methods, with the advantage of fewer interactions and a faster convergence.


title: Attention Deep Model with Multi-Scale Deep Supervision for Person Re-Identification

abstract: In recent years, person re-identification (PReID) has become a hot topic in computer vision duo to it is an important part in intelligent surveillance. Many state-of-the-art PReID methods are attention-based or multi-scale feature learning deep models. However, introducing attention mechanism may lead to some important feature information losing issue. Besides, most of the multi-scale models embedding the multi-scale feature learning block into the feature extraction deep network, which reduces the efficiency of inference network. To address these issue, in this study, we introduce an attention deep architecture with multi-scale deep supervision for PReID. Technically, we contribute a reverse attention block to complement the attention block, and a novel multi-scale layer with deep supervision operator for training the backbone network. The proposed block and operator are only used for training, and discard in test phase. Experiments have been performed on Market-1501, DukeMTMC-reID and CUHK03 datasets. All the experiment results show that the proposed model significantly outperforms the other competitive state-of-the-art methods.


title: Hopf-Galois structures on finite extensions with almost simple Galois group

abstract: In this paper, we study the Hopf-Galois structures on a finite Galois extension whose Galois group $G$ is an almost simple group in which the socle $A$ has prime index $p$. Each Hopf-Galois structure is associated to a group $N$ of the same order as $G$. We give necessary criteria on these $N$ in terms of their group-theoretic properties, and determine the number of Hopf-Galois structures associated to $A\times C_p$, where $C_p$ is the cyclic group of order $p$.


title: Echoing the recent Google success: Foundational Roots of Quantum Supremacy

abstract: The recent Google's claim on breakthrough in quantum computing is a gong signal for further analysis of foundational roots of (possible) superiority of some quantum algorithms over the corresponding classical algorithms. This note is a step in this direction. We start with critical analysis of rather common reference to entanglement and quantum nonlocality as the basic sources of quantum superiority. We elevate the role of the Bohr's {\it principle of complementarity}\footnote{} (PCOM) by interpreting the Bell-experiments as statistical tests of this principle. (Our analysis also includes comparison of classical vs genuine quantum entanglements.) After a brief presentation of PCOM and endowing it with the information interpretation, we analyze its computational counterpart. The main implication of PCOM is that by using the quantum representation of probability, one need not compute the joint probability distribution (jpd) for observables involved in the process of computation. Jpd's calculation is exponentially time consuming. Consequently, classical probabilistic algorithms involving calculation of jpd for $n$ random variables can be over-performed by quantum algorithms (for big values of $n).$ Quantum algorithms are based on quantum probability calculus. It is crucial that the latter modifies the classical formula of total probability (FTP). Probability inference based on the quantum version of FTP leads to constructive interference of probabilities increasing probabilities of some events. We also stress the role the basic feature of the genuine quantum superposition comparing with the classical wave superposition: generation of discrete events in measurements on superposition states. Finally, the problem of superiority of quantum computations is coupled with the quantum measurement problem and linearity of dynamics of the quantum state update.


title: Interarea Oscillations & Chimera in Power Systems

abstract: This paper proposes a novel second order mathematical model in the Kuramoto framework to simulate and study low frequency oscillations in power systems. This model facilitates better understanding of the complex dynamics of a power network. A standard four generator power system with all-to-all connectivity is considered and results obtained from the proposed model are verified. It is shown, that the model simulates various properties related to low frequency oscillations in power systems which presently are obtained through small-signal analysis. Further, we provide analogy to blackouts in a power grid, by emulating chimera behavior and thereby discuss bifurcation analysis of the proposed model.


title: Forecasting vegetation condition for drought early warning systems in pastoral communities in Kenya

abstract: Droughts are a recurring hazard in sub-Saharan Africa, that can wreak huge socioeconomic costs.Acting early based on alerts provided by early warning systems (EWS) can potentially provide substantial mitigation, reducing the financial and human cost. However, existing EWS tend only to monitor current, rather than forecast future, environmental and socioeconomic indicators of drought, and hence are not always sufficiently timely to be effective in practice. Here we present a novel method for forecasting satellite-based indicators of vegetation condition. Specifically, we focused on the 3-month Vegetation Condition Index (VCI3M) over pastoral livelihood zones in Kenya, which is the indicator used by the Kenyan National Drought Management Authority(NDMA). Using data from MODIS and Landsat, we apply linear autoregression and Gaussian process modeling methods and demonstrate high forecasting skill several weeks ahead. As a benchmark we predicted the drought alert marker used by NDMA (VCI3M<35). Both of our models were able to predict this alert marker four weeks ahead with a hit rate of around 89% and a false alarm rate of around 4%, or 81% and 6% respectively six weeks ahead. The methods developed here can thus identify a deteriorating vegetation condition well and sufficiently in advance to help disaster risk managers act early to support vulnerable communities and limit the impact of a drought hazard.


title: Oriented Diameter of Star Graphs

abstract: An {\em orientation} of an undirected graph $G$ is an assignment of exactly one direction to each edge of $G$. Converting two-way traffic networks to one-way traffic networks and bidirectional communication networks to unidirectional communication networks are practical instances of graph orientations. In these contexts minimising the diameter of the resulting oriented graph is of prime interest. The $n$-star network topology was proposed as an alternative to the hypercube network topology for multiprocessor systems by Akers and Krishnamurthy [IEEE Trans. on Computers (1989)]. The $n$-star graph $S_n$ consists of $n!$ vertices, each labelled with a distinct permutation of $[n]$. Two vertices are adjacent if their labels differ exactly in the first and one other position. $S_n$ is an $(n-1)$-regular, vertex-transitive graph with diameter $\lfloor 3(n-1)/2 \rfloor$. Orientations of $S_n$, called unidirectional star graphs and distributed routing protocols over them were studied by Day and Tripathi [Information Processing Letters (1993)] and Fujita [The First International Symposium on Computing and Networking (CANDAR 2013)]. Fujita showed that the (directed) diameter of this unidirectional star graph $\overrightarrow{S_n}$ is at most $\lceil{5n/2}\rceil + 2$. In this paper, we propose a new distributed routing algorithm for the same $\overrightarrow{S_n}$ analysed by Fujita, which routes a packet from any node $s$ to any node $t$ at an undirected distance $d$ from $s$ using at most $\min\{4d+4, 2n+4\}$ hops. This shows that the (directed) diameter of $\overrightarrow{S_n}$ is at most $2n+4$. We also show that the diameter of $\overrightarrow{S_n}$ is at least $2n$ when $n \geq 7$, thereby showing that our upper bound is tight up to an additive factor.