title: Temporal and Angular Variations of 3D Core-Collapse Supernova Emissions and their Physical Correlations

abstract: We provide the time series and angular distributions of the neutrino and gravitational-wave emissions of eleven state-of-the-art three-dimensional non-rotating core-collapse supernova models and explore correlations between these signatures and the real-time dynamics of the shock and the proto-neutron-star core. The neutrino emissions are roughly isotropic on average, with instantaneous excursions about the mean inferred luminosity of as much as $\pm$20%. The deviation from isotropy is least for the "$\nu_{\mu}$"-type neutrinos and the lowest-mass progenitors. Instantaneous temporal luminosity variations along a given direction for exploding models average $\sim$2$-$4%, but can be as high as $\sim$10%. For non-exploding models, they can achieve $\sim$25%. The temporal variations in the neutrino emissions correlate with the temporal and angular variations in the mass accretion rate. We witness the LESA phenomenon in all our models and find that the vector direction of the LESA dipole and that of the inner Y$_\mathrm{e}$ distribution are highly correlated. For our entire set of 3D models, we find strong connections between the cumulative neutrino energy losses, the radius of the proto-neutron star, and the $f$-mode frequency of the gravitational wave emissions. When physically normalized, the progenitor-to-progenitor variation in any of these quantities is no more than $\sim$10%. Moreover, the reduced $f$-mode frequency is independent of time after bounce to better than $\sim$10%. Therefore, simultaneous measurement of gravitational waves and neutrinos from a given supernova event can be used synergistically to extract real physical quantities of the supernova core.


title: Orphaned Exomoons: Tidal Detachment and Evaporation Following an Exoplanet-Star Collision

abstract: Gravitational perturbations on an exoplanet from a massive outer body, such as the Kozai-Lidov mechanism, can pump the exoplanet's eccentricity up to values that will destroy it via a collision or strong interaction with its parent star. During the final stages of this process, any exomoons orbiting the exoplanet will be detached by the star's tidal force and placed into orbit around the star. Using ensembles of three and four-body simulations, we demonstrate that while most of these detached bodies either collide with their star or are ejected from the system, a substantial fraction, ~10%, of such "orphaned" exomoons (with initial properties similar to those of the Galilean satellites in our own solar system) will outlive their parent exoplanet. The detached exomoons generally orbit inside the ice line, so that strong radiative heating will evaporate any volatile-rich layers, producing a strong outgassing of gas and dust, analogous to a comet's perihelion passage. Small dust grains ejected from the exomoon may help generate an opaque cloud surrounding the orbiting body but are quickly removed by radiation blow-out. By contrast, larger solid particles inherit the orbital properties of the parent exomoon, feeding an eccentric disk of solids that drains more gradually onto the star via Poynting-Robertson drag, and which could result in longer-timescale dimming of the star. For characteristic exomoon evaporation times of ~ 1e5-1e6 yr, attenuation of the stellar light arising from one or more out-gassing exomoons provides a promising explanation for both the dipping and secular dimming behavior observed from KIC 8462852 (Boyajian's Star).


title: Numerical binary black hole collisions in dynamical Chern-Simons gravity

abstract: We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational wave detection. Head-on collisions are interesting in their own right, however, as they cleanly probe the quasi-normal mode spectrum of the final black hole. We thus compute the leading-order dynamical Chern-Simons modifications to the complex frequencies of the post-merger gravitational radiation. We consider equal-mass systems, with equal spins oriented along the axis of collision, resulting in remnant black holes with spin. We find modifications to the complex frequencies of the quasi-normal mode spectrum that behave as a power law with the spin of the remnant, and that are not degenerate with the frequencies associated with a Kerr black hole of any mass and spin. We discuss these results in the context of testing general relativity with gravitational wave observations.


title: Multisymplectic actions of compact Lie groups on spheres

abstract: We investigate the existence of homotopy comoment maps (comoments) for high-dimensional spheres seen as multisymplectic manifolds. Especially, we solve the existence problem for compact effective group actions on spheres and provide explicit constructions for such comoments in interesting particular cases.


title: Manipulation of confined polyelectrolyte conformations through dielectric mismatch

abstract: We demonstrate that a highly charged polyelectrolyte confined in a spherical cavity undergoes reversible transformations between amorphous conformations to a four-fold symmetry morphology as a function of dielectric mismatch between the media inside and outside the cavity. Surface polarization due to dielectric mismatch exhibits an extra "confinement" effect, which are most pronounced within a certain range of the cavity radius and the electrostatic strength between the monomers and counterions and multivalent counterions. For cavities with a charged surface, surface polarization leads to an increased amount of counterions adsorbed in the outer side, further compressing the confined polyelectrolyte into a four-fold symmetry morphology. The equilibrium conformation of the chain is dependent upon several key factors including the relative permittivities of the media inside and outside the cavity, multivalent counterion concentration, cavity radius relative to the chain length, and interface charge density. Our findings offer insights into the effects of dielectric mismatch in packaging and delivery of polyelectrolytes across media with different relative permittivities. Moreover, the reversible transformation of the polyelectrolyte conformations in response to environmental permittivity allows for potential applications in biosensing and medical monitoring.


title: When Multiple Agents Learn to Schedule: A Distributed Radio Resource Management Framework

abstract: Interference among concurrent transmissions in a wireless network is a key factor limiting the system performance. One way to alleviate this problem is to manage the radio resources in order to maximize either the average or the worst-case performance. However, joint consideration of both metrics is often neglected as they are competing in nature. In this article, a mechanism for radio resource management using multi-agent deep reinforcement learning (RL) is proposed, which strikes the right trade-off between maximizing the average and the $5^{th}$ percentile user throughput. Each transmitter in the network is equipped with a deep RL agent, receiving partial observations from the network (e.g., channel quality, interference level, etc.) and deciding whether to be active or inactive at each scheduling interval for given radio resources, a process referred to as link scheduling. Based on the actions of all agents, the network emits a reward to the agents, indicating how good their joint decisions were. The proposed framework enables the agents to make decisions in a distributed manner, and the reward is designed in such a way that the agents strive to guarantee a minimum performance, leading to a fair resource allocation among all users across the network. Simulation results demonstrate the superiority of our approach compared to decentralized baselines in terms of average and $5^{th}$ percentile user throughput, while achieving performance close to that of a centralized exhaustive search approach. Moreover, the proposed framework is robust to mismatches between training and testing scenarios. In particular, it is shown that an agent trained on a network with low transmitter density maintains its performance and outperforms the baselines when deployed in a network with a higher transmitter density.


title: Limits, standard complexes and fr-codes

abstract: For a strongly connected category $\mathcal C$ with pair-wise coproducts, we introduce a cosimplicial object, which serves as a sort of resolution for computing higher derived functors of ${\sf lim} : \mathrm{Ab}^{\mathcal C}\to \mathrm{Ab}$. Applications involve K\"unneth theorem for higher limits and ${\sf lim}$-finiteness of ${\bf fr}$-codes. A dictionary for the ${\bf fr}$-codes with words of length $\leq 3$ is given.


title: Fast Neutrino Flavor Conversion: Collective Motion vs. Decoherence

abstract: In an interacting neutrino gas, flavor coherence becomes dynamical and can propagate as a collective mode. In particular, tachyonic instabilities can appear, leading to "fast flavor conversion" that is independent of neutrino masses and mixing angles. On the other hand, without neutrino-neutrino interaction, a prepared wave packet of flavor coherence simply dissipates by kinematical decoherence of infinitely many non-collective modes. We reexamine the dispersion relation for fast flavor modes and show that for any wavenumber,there exists a continuum of non-collective modes besides a few discrete collective ones. So for any initial wave packet, both decoherence and collective motion occurs, although the latter typically dominates for a sufficiently dense gas. We derive explicit eigenfunctions for both collective and non-collective modes. If the angular mode distribution of electron-lepton number crosses between positive and negative values, two non-collective modes can merge to become a tachyonic collective mode. We explicitly calculate the interaction strength for this critical point. As a corollary we find that a single crossing always leads to a tachyonic instability. For an even number of crossings, no instability needs to occur.


title: Search for gas from the disintegrating rocky exoplanet K2-22b

abstract: [Abridged] Aims. We searched for circumplanetary sodium and ionized calcium gas around the disintegrating rocky exoplanet K2-22 b to constrain its gas-loss and sublimation processes. Methods. We observed four transits of K2-22 b with X-shooter on ESO's Very Large Telescope to obtain time-series of intermediate-resolution (R $\sim$ 11400) spectra. Our analysis focused on the two sodium D lines (588.995 nm and 589.592 nm) and the Ca$^{+}$ triplet (849.802 nm, 854.209 nm and 866.214 nm). Planet-related absorption is searched for in the velocity rest frame of the planet, which changes from $\pm$66 kms$^{-1}$ during the transit. Results. Since K2-22 b exhibits highly variable transit depths, we analyzed the individual nights and their average. By injecting signals we reached 5$\sigma$ upper-limits on the individual nights that ranged from 11% - 13% and 1.7% - 2.0% for the tail's sodium and ionized calcium absorption, respectively. Night 1 was contaminated by its companion star so we considered weighted averages with and without Night 1 and quote conservative 5$\sigma$ limits without Night 1 of 9% and 1.4%, respectively. Assuming their mass fractions to be similar to those in the Earth's crust, these limits correspond to scenarios in which 0.04% and 35% of the transiting dust is sublimated and observed as absorbing gas. However, this assumes the gas to be co-moving with the planet. We show that for the high irradiation environment of K2-22 b, sodium and ionized calcium could be quickly accelerated to 100s of km s$^{-1}$ due to radiation pressure and entrainment by the stellar wind, making them much more difficult to detect. No evidence for such possibly broad and blue-shifted signals are seen in our data. Conclusions. Future observations aimed at observing circumplanetary gas should take into account the possible broad and blue-shifted velocity field of atomic and ionized species.


title: Positive mass theorem for initial data sets with corners along a hypersurface

abstract: We prove positive mass theorem with angular momentum and charges for axially symmetric, simply connected, maximal, complete initial data sets with two ends, one designated asymptotically flat and the other either (Kaluza-Klein) asymptotically flat or asymptotically cylindrical, for 4-dimensional Einstein-Maxwell theory and $5$-dimensional minimal supergravity theory which metrics fail to be $C^1$ and second fundamental forms and electromagnetic fields fail to be $C^0$ across an axially symmetric hypersurface $\Sigma$. Furthermore, we remove the completeness and simple connectivity assumptions in this result and prove it for manifold with boundary such that the mean curvature of the boundary is non-positive.


title: The REASONS Survey: Resolved Millimeter Observations of a Large Debris Disk Around the Nearby F Star HD 170773

abstract: Debris disks are extrasolar analogs to our own Kuiper Belt and they are detected around at least 17% of nearby Sun-like stars. The morphology and dynamics of a disk encode information about its history, as well as that of any exoplanets within the system. We used ALMA to obtain 1.3 mm observations of the debris disk around the nearby F5V star HD 170773. We image the face-on ring and determine its fundamental parameters by forward-modeling the interferometric visibilities through a Markov Chain Monte Carlo approach. Using a symmetric Gaussian surface density profile, we find a 71 $\pm$ 4 au wide belt with a radius of 193$^{+2}_{-3}$ au, a relatively large radius compared to most other millimeter-resolved belts around late A / early F type stars. This makes HD 170773 part of a group of four disks around A and F stars with radii larger than expected from the recently reported planetesimal belt radius - stellar luminosity relation. Two of these systems are known to host directly imaged giant planets, which may point to a connection between large belts and the presence of long-period giant planets. We also set upper limits on the presence of CO and CN gas in the system, which imply that the exocomets that constitute this belt have CO and HCN ice mass fractions of <77% and <3%, respectively, consistent with Solar System comets and other exocometary belts.


title: Buffer-gas cooling, high-resolution spectroscopy and optical cycling of barium monofluoride molecules

abstract: We demonstrate buffer-gas cooling, high-resolution spectroscopy and cycling fluorescence of cold barium monofluoride (BaF) molecules. Our source produces an intense and internally cold molecular beam containing the different BaF isotopologues with a mean forward velocity of 190 m/s. For a well-collimated beam of 138BaF we observe a flux of more than 1e10 molecules/sr/pulse in the X2Sigma, N=1 state in our downstream detection region. Studying the absorption line strength of the intermediate A'Delta state we infer a lifetime of 790+\-346 ns, significantly longer than previously estimated. Finally, highly-diagonal Franck-Condon factors and magnetic remixing of dark states allow us to realize a quasi-cycling transition in 138BaF that is suitable for future laser cooling of this heavy diatomic molecule.


title: Posterior Contraction Rates for Gaussian Cox Processes with Non-identically Distributed Data

abstract: This paper considers the posterior contraction of non-parametric Bayesian inference on non-homogeneous Poisson processes. We consider the quality of inference on a rate function $\lambda$, given non-identically distributed realisations, whose rates are transformations of $\lambda$. Such data arises frequently in practice due, for instance, to the challenges of making observations with limited resources or the effects of weather on detectability of events. We derive contraction rates for the posterior estimates arising from the Sigmoidal Gaussian Cox Process and Quadratic Gaussian Cox Process models. These are popular models where $\lambda$ is modelled as a logistic and quadratic transformation of a Gaussian Process respectively. Our work extends beyond existing analyses in several regards. Firstly, we consider non-identically distributed data, previously unstudied in the Poisson process setting. Secondly, we consider the Quadratic Gaussian Cox Process model, of which there was previously little theoretical understanding. Thirdly, we provide rates on the shrinkage of both the width of balls around the true $\lambda$ in which the posterior mass is concentrated and on the shrinkage of posterior mass outside these balls - usually only the former is explicitly given. Finally, our results hold for certain finite numbers of observations, rather than only asymptotically, and we relate particular choices of hyperparameter/prior to these results.


title: Hardness and Ease of Curing the Sign Problem for Two-Local Qubit Hamiltonians

abstract: We examine the problem of determining whether a multi-qubit two-local Hamiltonian can be made stoquastic by single-qubit unitary transformations. We prove that when such a Hamiltonian contains one-local terms, then this task can be NP-hard. This is shown by constructing a class of Hamiltonians for which performing this task is equivalent to deciding $3$-SAT. In contrast, we show that when such a Hamiltonian contains no one-local terms then this task is easy, namely we present an algorithm which decides, in a number of arithmetic operations over $\mathbb{R}$ which is polynomial in the number of qubits, whether the sign problem of the Hamiltonian can be cured by single-qubit rotations.


title: Regularization, renormalization and consistency conditions in QED with x-electric potential steps

abstract: The present article is an important addition to the nonperturbative formulation of QED with x-steps presented by Gavrilov and Gitman in Phys. Rev. D. 93, 045002 (2016). Here we propose a new renormalization and volume regularization procedures which allow one to calculate and distinguish physical parts of different matrix elements of operators of the current and of the energy-momentum tensor, at the same time relating the latter quantities with characteristics of the vacuum instability. For this purpose, a modified inner product and a parameter {\tau} of the regularization are introduced. The latter parameter can be fixed using physical considerations. In the Klein zone this parameter can be interpreted as the time of the observation of the pair production effect. In the refined formulation of QED with x-steps, we succeeded to consider the backreaction problem. In the case of an uniform electric field E confined between two capacitor plates separated by a finite distance L, we see that the smallness of the backreaction implies a restriction (the consistency condition) on the product EL from above.


title: Parameter-robust Multiphysics Algorithms for Biot Model with Application in Brain Edema Simulation

abstract: In this paper, we develop two parameter-robust numerical algorithms for Biot model and applied the algorithms in brain edema simulations. By introducing an intermediate variable, we derive a multiphysics reformulation of the Biot model. Based on the reformulation, the Biot model is viewed as a generalized Stokes subproblem combining with a reaction-diffusion subproblem. Solving the two subproblems together or separately will lead to a coupled or a decoupled algorithm. We conduct extensive numerical experiments to show that the two algorithms are robust with respect to the physics parameters. The algorithms are applied to study the brain swelling caused by abnormal accumulation of cerebrospinal fluid in injured areas. The effects of key physics parameters on brain swelling are carefully investigated. It is observe that the permeability has the greatest effect on intracranial pressure (ICP) and tissue deformation; the Young's modulus and the Poisson ratio will not affect the maximum ICP too much but will affect the tissue deformation and the developing speed of brain swelling.


title: A $\phi^6$ soliton with a long-range tail

abstract: We propose an analytically solvable sextic potential model with non-trivial soliton solutions connecting the trivial vacua. The model does not respect parity symmetry, and like $\phi^4$ theory has two minima. The soliton solutions and the consequent results are obtained in terms of the Lambert W function, i.e., the inverse function of $f(W) = We^W$. They have power-law asymptotics at one spatial infinity and exponential asymptotics at the other. We compare the solution with the kink of $\phi^4$ theory, which preserves the parity symmetry and has exponential asymptotics at both spatial infinities. Moreover, we study the full spectrum (bound and continuum states) of boson and fermion fields in the presence of the proposed soliton. We consider two types of coupling for the boson-soliton interaction and Yukawa coupling for the fermion-soliton interaction. Most results are derived analytically. This property renders the model a fertile ground for further study, including parity breaking related phenomena and long-range soliton-soliton interactions.


title: Derivation of the Variational Bayes Equations

abstract: The derivation of key equations for the variational Bayes approach is well-known in certain circles. However, translating the fundamental derivations (e.g., as found in Beal's work) to Friston's notation is somewhat delicate. Further, the notion of using variational Bayes in the context of a system with a Markov blanket requires special attention. This Technical Report presents the derivation in detail. It further illustrates how the variational Bayes method provides a framework for a new computational engine, incorporating the 2-D cluster variation method (CVM), which provides a necessary free energy equation that can be minimized across both the external and representational systems' states, respectively.


title: Finding Needles in a Moving Haystack: Prioritizing Alerts with Adversarial Reinforcement Learning

abstract: Detection of malicious behavior is a fundamental problem in security. One of the major challenges in using detection systems in practice is in dealing with an overwhelming number of alerts that are triggered by normal behavior (the so-called false positives), obscuring alerts resulting from actual malicious activity. While numerous methods for reducing the scope of this issue have been proposed, ultimately one must still decide how to prioritize which alerts to investigate, and most existing prioritization methods are heuristic, for example, based on suspiciousness or priority scores. We introduce a novel approach for computing a policy for prioritizing alerts using adversarial reinforcement learning. Our approach assumes that the attackers know the full state of the detection system and dynamically choose an optimal attack as a function of this state, as well as of the alert prioritization policy. The first step of our approach is to capture the interaction between the defender and attacker in a game theoretic model. To tackle the computational complexity of solving this game to obtain a dynamic stochastic alert prioritization policy, we propose an adversarial reinforcement learning framework. In this framework, we use neural reinforcement learning to compute best response policies for both the defender and the adversary to an arbitrary stochastic policy of the other. We then use these in a double-oracle framework to obtain an approximate equilibrium of the game, which in turn yields a robust stochastic policy for the defender. Extensive experiments using case studies in fraud and intrusion detection demonstrate that our approach is effective in creating robust alert prioritization policies.


title: The Moran forest

abstract: Starting from any graph on $\{1, \ldots, n\}$, consider the Markov chain where at each time-step a uniformly chosen vertex is disconnected from all of its neighbors and reconnected to another uniformly chosen vertex. This Markov chain has a stationary distribution whose support is the set of non-empty forests on $\{1, \ldots, n\}$. The random forest corresponding to this stationary distribution has interesting connections with the uniform rooted labeled tree and the uniform attachment tree. We fully characterize its degree distribution, the distribution of its number of trees, and the limit distribution of the size of a tree sampled uniformly. We also show that the size of the largest tree is asymptotically $\alpha \log n$, where $\alpha = (1 - \log(e - 1))^{-1} \approx 2.18$, and that the degree of the most connected vertex is asymptotically $\log n / \log\log n$.


title: Detecting Two Qubit Both-Way Positive Discord States

abstract: Quantum discord plays a pragmatic role in analyzing nonclassical feature of quantum correlations beyond entanglement. It is used in several information processing protocols which lacks sufficient amount of entanglement to be used as a resource. We have provided with an analytical method of detecting quantum discord of an arbitrary two qubit state. We have formulated a set of necessary and sufficient conditions for any two qubit state to be a both-way non-zero quantum discord state. As quantum discord is asymmetric in nature, we have framed the set of if and only if conditions for a two qubit state to be classical-quantum as well for it to be quantum-classical. Interestingly, not only correlation tensor but also local Bloch vector(corresponding to the classical party) plays a role for detecting the state to be a positive discord state.


title: Observable quantum entanglement due to gravity

abstract: No experiment to date has provided evidence for quantum features of the gravitational interaction. Recently proposed tests suggest looking for the generation of quantum entanglement between massive objects as a possible route towards the observation of such features. Motivated by advances in optical cooling of mirrors, here we provide a systematic study of entanglement between two masses that are coupled gravitationally. We first consider the masses trapped at all times in harmonic potentials (optomechanics) and then the masses released from the traps. This leads to the estimate of the experimental parameters required for the observation of gravitationally induced entanglement. The optomechanical setup demands LIGO-like mirrors and squeezing or long coherence times, but the released masses can be light and accumulate detectable entanglement in a timescale shorter than their coherence times. No macroscopic quantum superposition develops during the evolution. We discuss the implications from such thought experiments regarding the nature of the gravitational coupling.


title: A Deep Reinforcement Learning Approach for Global Routing

abstract: Global routing has been a historically challenging problem in electronic circuit design, where the challenge is to connect a large and arbitrary number of circuit components with wires without violating the design rules for the printed circuit boards or integrated circuits. Similar routing problems also exist in the design of complex hydraulic systems, pipe systems and logistic networks. Existing solutions typically consist of greedy algorithms and hard-coded heuristics. As such, existing approaches suffer from a lack of model flexibility and non-optimum solutions. As an alternative approach, this work presents a deep reinforcement learning method for solving the global routing problem in a simulated environment. At the heart of the proposed method is deep reinforcement learning that enables an agent to produce an optimal policy for routing based on the variety of problems it is presented with leveraging the conjoint optimization mechanism of deep reinforcement learning. Conjoint optimization mechanism is explained and demonstrated in details; the best network structure and the parameters of the learned model are explored. Based on the fine-tuned model, routing solutions and rewards are presented and analyzed. The results indicate that the approach can outperform the benchmark method of a sequential A* method, suggesting a promising potential for deep reinforcement learning for global routing and other routing or path planning problems in general. Another major contribution of this work is the development of a global routing problem sets generator with the ability to generate parameterized global routing problem sets with different size and constraints, enabling evaluation of different routing algorithms and the generation of training datasets for future data-driven routing approaches.


title: A New Achievable Rate-Distortion Region for Distributed Source Coding

abstract: In this work, lossy distributed compression of pairs of correlated sources is considered. Conventionally, Shannon's random coding arguments -- using randomly generated unstructured codebooks whose blocklength is taken to be asymptotically large -- are used to derive achievability results. However, it was recently observed that in various multi-terminal communications scenarios, using random codes with constant finite blocklength may lead to improved achievable regions compared to the conventional approach. In other words, in some network communication scenarios, there is a finite optimal value in the blocklength of the randomly generated code used for distributed processing of information sources. Motivated by this, a coding scheme is proposed which consists of two codebook layers: i) the primary codebook which has constant finite blocklength, and ii) the secondary codebook whose blocklength is taken to be asymptotically large. The achievable region is analyzed in two steps. In the first step, a characterization of the achievable region is derived using information measures which are functions of multi-letter probability distributions. In the next step, a computable single-letter inner-bound to the achievable region is derived. It is shown through several examples that the resulting rate-distortion region is strictly larger than the Berger Tung achievable region.


title: Multimode thermal states with multiphoton subtraction: study of the photons number distribution in the selected subsystem

abstract: Thermal states of light are widely used in quantum optics due to their correlation properties. As is well known, their correlation properties and the photon number distribution as a whole are strongly dependent on the mode number selected by the detection scheme. The same changes can be caused by photon subtraction. Therefore, we describe the general case of the multimode thermal state after a multiphoton subtraction, when the photon number statistics is registered by the detector selecting a part of the initial modes. We present an analytical form of the obtained photon number distribution and its general properties and check them in the experiment.


title: Cache-Aided NOMA Mobile Edge Computing: A Reinforcement Learning Approach

abstract: A novel non-orthogonal multiple access (NOMA) based cache-aided mobile edge computing (MEC) framework is proposed. For the purpose of efficiently allocating communication and computation resources to users' computation tasks requests, we propose a long-short-term memory (LSTM) network to predict the task popularity. Based on the predicted task popularity, a long-term reward maximization problem is formulated that involves a joint optimization of the task offloading decisions, computation resource allocation, and caching decisions. To tackle this challenging problem, a single-agent Q-learning (SAQ-learning) algorithm is invoked to learn a long-term resource allocation strategy. Furthermore, a Bayesian learning automata (BLA) based multi-agent Q-learning (MAQ-learning) algorithm is proposed for task offloading decisions. More specifically, a BLA based action select scheme is proposed for the agents in MAQ-learning to select the optimal action in every state. We prove that the BLA based action selection scheme is instantaneously self-correcting and the selected action is an optimal solution for each state. Extensive simulation results demonstrate that: 1) The prediction error of the proposed LSTMs based task popularity prediction decreases with increasing learning rate. 2) The proposed framework significantly outperforms the benchmarks like all local computing, all offloading computing, and non-cache computing. 3) The proposed BLA based MAQ-learning achieves an improved performance compared to conventional reinforcement learning algorithms.


title: High-energy cosmic rays from compact galactic star clusters: particle fluxes and anisotropy

abstract: It has been shown that supernova blast waves interacting with winds from massive stars in compact star clusters may be capable of producing cosmic-ray (CR) protons to above $10^{17}$ eV. We give a brief description of the colliding-shock-flows mechanism and look at generalizations of the diffusion of ~ 100 PeV CRs in the turbulent galactic magnetic field present in the galactic disk. We calculate the temporal evolution of the CR anisotropy from a possible distribution of young compact massive star clusters assuming the sources are intermittent on time scales of a few million years, i.e., comparable to their residence time in the Milky Way. Within the confines of our model, we determine the galactic/extra-galactic fraction of high-energy CRs resulting in anisotropies consistent with observed values. We find that galactic star clusters may contribute a substantial fraction of ~ 100 PeV CRs without producing anisotropies above observed limits.


title: Exclusion Limits on Hidden-Photon Dark Matter near 2 neV from a Fixed-Frequency Superconducting Lumped-Element Resonator

abstract: We present the design and performance of a simple fixed-frequency superconducting lumped-element resonator developed for axion and hidden photon dark matter detection. A rectangular NbTi inductor was coupled to a Nb-coated sapphire capacitor and immersed in liquid helium within a superconducting shield. The resonator was transformer-coupled to a DC SQUID for readout. We measured a quality factor of $\sim$40,000 at the resonant frequency of 492.027 kHz and set a simple exclusion limit on $\sim$2 neV hidden photons with kinetic mixing angle $\varepsilon\gtrsim1.5\times10^{-9}$ based on 5.14 hours of integrated noise. This test device informs the development of the Dark Matter Radio, a tunable superconducting lumped-element resonator which will search for axions and hidden photons over the 100 Hz to 300 MHz frequency range.


title: Electroweak pseudo-observables and Z-boson form factors at two-loop accuracy

abstract: We present Standard Model predictions for the complete set of phenomenologically relevant electroweak precision pseudo-observables related to the Z-boson: the leptonic and bottom-quark effective weak mixing angles $\sin^2\theta_{\rm eff}^\ell$, $\sin^2\theta_{\rm eff}^b$, the Z-boson partial decay widths $\Gamma_f$, where $f$ indicates any charged lepton, neutrino and quark flavor (except for the top quark), as well as the total Z decay width $\Gamma_Z$, the branching ratios $R_\ell$, $R_c$, $R_b$, and the hadronic cross section $\sigma_{\rm had}^0$. The input parameters are the masses $M_Z$, $M_H$ and $m_t$, and the couplings $\alpha_s$, $\alpha$. The scheme dependence due to the choice of $M_W$ or its alternative $G_\mu$ as a last input parameter is also discussed. Recent substantial technical progress in the calculation of Minkowskian massive higher-order Feynman integrals allows the calculation of the complete electroweak two-loop radiative corrections to all the observables mentioned. QCD contributions are included appropriately. Results are provided in terms of simple and convenient parameterization formulae whose coefficients have been determined from the full numerical multi-loop calculation. The size of the missing electroweak three-loop or QCD higher-order corrections is estimated. We briefly comment on the prospects for their calculation. Finally, direct predictions for the $Z{\bar f}f$ vector and axial-vector form-factors are given, including a discussion of separate order-by-order contributions.


title: Long time asymptotics for homoenergetic solutions of the Boltzmann equation. Hyperbolic-dominated case

abstract: In this paper we continue the formal analysis of the long-time asymptotics of the homoenergetic solutions for the Boltzmann equation that we began in [18]. They have the form $f\left( x,v,t\right) =g\left(v-L\left( t\right) x,t\right) $ where $L\left( t\right) =A\left(I+tA\right) ^{-1}$ where $A$ is a constant matrix. Homoenergetic solutions satisfy an integro-differential equation which contains, in addition to the classical Boltzmann collision operator, a linear hyperbolic term. Depending on the properties of the collision kernel the collision and the hyperbolic terms might be of the same order of magnitude as $t\to\infty$, or the collision term could be the dominant one for large times, or the hyperbolic term could be the largest. The first case has been rigorously studied in [17]. Formal asymptotic expansions in the second case have been obtained in [18]. All the solutions obtained in this case can be approximated by Maxwellian distributions with changing temperature. In this paper we focus in the case where the hyperbolic terms are much larger than the collision term for large times (hyperbolic-dominated behavior). In the hyperbolic-dominated case it does not seem to be possible to describe in a simple way all the long time asymptotics of the solutions, but we discuss several physical situations and formulate precise conjectures. We give explicit formulas for the relationship between density, temperature and entropy for these solutions. These formulas differ greatly from the ones at equilibrium.


title: Self-conjugate representation SU(3) chains

abstract: It was recently argued that SU(3) chains in the $p$-box symmetric irreducible representation (irrep) exhibit a "Haldane gap" when $p$ is a multiple of 3 and are otherwise gapless [Nucl. Phys. B $\boldsymbol{924}$, 508 (2017)]. We extend this argument to the self-conjugate irreps of SU(3) with $p$ columns of length 2 and $p$ columns of length 1 in the Young tableau ($p=1$ corresponding to the adjoint irrep), arguing that they are always gapped but have spontaneously broken parity symmetry for $p$ odd but not even.


title: Pell surfaces

abstract: In 1826 Abel started the study of the polynomial Pell equation $x^2-g(u)y^2=1$. Its solvability in polynomials $x(u), y(u)$ depends on a certain torsion point on the Jacobian of the hyperelliptic curve $v^2=g(u)$. In this paper we study the affine surfaces defined by the Pell equations in 3-space with coordinates $x, y,u$, and aim to describe all affine lines on it. These are polynomial solutions of the equation $x(t)^2-g(u(t))y(t)^2=1$. Our results are rather complete when the degree of $g$ is even but the odd degree cases are left completely open. For even degrees we also describe all curves on these Pell surfaces that have only 1 place at infinity.


title: Performance Comparison Between OpenCV Built in CPU and GPU Functions on Image Processing Operations

abstract: Image Processing is a specialized area of Digital Signal Processing which contains various mathematical and algebraic operations such as matrix inversion, transpose of matrix, derivative, convolution, Fourier Transform etc. Operations like those require higher computational capabilities than daily usage purposes of computers. At that point, with increased image sizes and more complex operations, CPUs may be unsatisfactory since they use Serial Processing by default. GPUs are the solution that come up with greater speed compared to CPUs because of their Parallel Processing/Computation nature. A parallel computing platform and programming model named CUDA was created by NVIDIA and implemented by the graphics processing units (GPUs) which were produced by them. In this paper, computing performance of some commonly used Image Processing operations will be compared on OpenCV's built in CPU and GPU functions that use CUDA.


title: Non-relativistic expansion of single-nucleon Dirac equation: Comparison between Foldy-Wouthuysen transformation and similarity renormalization group

abstract: By following the Foldy-Wouthuysen (FW) transformation of the Dirac equation, we work out the exact analytic expressions up to the $1/M^4$ order for the general cases in the covariant density functional theory. These results are further compared with the corresponding ones derived from another novel non-relativistic expansion method, the similarity renormalization group (SRG). Based on that, the origin of the difference between the results obtained by the FW transformation and SRG method is explored.


title: A Hilbert-Kunz function with a periodic term that has a given period

abstract: A result of Monsky states that the Hilbert-Kunz function of a one-dimensional local ring of prime characteristic has a term $\phi$ that is eventually periodic. For example, in the case of a power series ring in one variable over a prime-characteristic field, $\phi$ is the zero function and is therefore immediately periodic with period 1. In additional examples produced by Kunz and Monsky, $\phi$ is immediately periodic with period 2. We show that, for every positive integer $\pi$, there exists a ring for which $\phi$ is immediately periodic with period $\pi$.


title: Stochastic Approximation Monte Carlo with a Dynamic Update Factor

abstract: We present a new Monte Carlo algorithm based on the Stochastic Approximation Monte Carlo (SAMC) algorithm for directly calculating the density of states. The proposed method is Stochastic Approximation with a Dynamic update factor (SAD) which dynamically adjusts the update factor $\gamma_t$ during the course of the simulation. We test this method on the square-well fluid and the 31-atom Lennard-Jones cluster and compare the convergence behavior of several related Monte Carlo methods. We find that both the SAD and $1/t$-Wang-Landau ($1/t$-WL) methods rapidly converge to the correct density of states without the need for the user to specify an arbitrary tunable parameter $t_0$ as in the case of SAMC. SAD requires as input the temperature range of interest, in contrast to $1/t$-WL, which requires that the user identify the interesting range of energies. The convergence of the $1/t$-WL method is very sensitive to the energy range chosen for the low-temperature heat capacity of the Lennard-Jones cluster. Thus, SAD is more powerful in the common case in which the range of energies is not known in advance.


title: Cross-Subject Statistical Shift Estimation for Generalized Electroencephalography-based Mental Workload Assessment

abstract: Assessment of mental workload in real-world conditions is key to ensure the performance of workers executing tasks that demand sustained attention. Previous literature has employed electroencephalography (EEG) to this end despite having observed that EEG correlates of mental workload vary across subjects and physical strain, thus making it difficult to devise models capable of simultaneously presenting reliable performance across users. Domain adaptation consists of a set of strategies that aim at allowing for improving machine learning systems performance on unseen data at training time. Such methods, however, might rely on assumptions over the considered data distributions, which typically do not hold for applications of EEG data. Motivated by this observation, in this work we propose a strategy to estimate two types of discrepancies between multiple data distributions, namely marginal and conditional shifts, observed on data collected from different subjects. Besides shedding light on the assumptions that hold for a particular dataset, the estimates of statistical shifts obtained with the proposed approach can be used for investigating other aspects of a machine learning pipeline, such as quantitatively assessing the effectiveness of domain adaptation strategies. In particular, we consider EEG data collected from individuals performing mental tasks while running on a treadmill and pedaling on a stationary bike and explore the effects of different normalization strategies commonly used to mitigate cross-subject variability. We show the effects that different normalization schemes have on statistical shifts and their relationship with the accuracy of mental workload prediction as assessed on unseen participants at training time.


title: RefleX: X-ray absorption and reflection in active galactic nuclei for arbitrary geometries

abstract: Reprocessed X-ray radiation carries important information about the structure and physical characteristics of the material surrounding the supermassive black hole (SMBH) in active galactic nuclei (AGN). We report here on a newly developed simulation platform, RefleX, which allows to reproduce absorption and reflection by quasi-arbitrary geometries. We show here the reliability of our approach by comparing the results of our simulations with existing spectral models such as pexrav, MYTorus and BNTorus. RefleX implements both Compton scattering on free electrons and Rayleigh scattering and Compton scattering on bound electrons. We show the effect of bound-electron corrections on a torus geometry simulated like in MYTorus. We release with this paper the RefleX executable, as well as RXTorus, a model that assumes absorption and reflection from a torus with a varying ratio of the minor to major axis of the torus. To allow major flexibility RXTorus is also distributed in three components: absorbed primary emission, scattered radiation and fluorescent lines. RXTorus is provided for different values of the abundance, and with (atomic configuration) or without (free-electron configuration) taking into account Rayleigh scattering and bound electrons. We apply the RXTorus model in both configurations on the XMM-Newton and NuSTAR spectrum of the Compton-thick AGN NGC 424 and find that the models are able to reproduce very well the observations, but that the assumption on the bound or free state of the electrons has significant consequences on the fit parameters.


title: Horizon Visibility Graphs and Time Series Merge Trees are Dual

abstract: In this paper we introduce the horizon visibility graph, a simple extension to the popular horizontal visibility graph representation of a time series, and show that it possesses a rigorous mathematical foundation in computational algebraic topology. This fills a longstanding gap in the literature on the horizontal visibility approach to nonlinear time series analysis which, despite a suite of successful applications across multiple domains, lacks a formal setting in which to prove general properties and develop natural extensions. The main finding is that horizon visibility graphs are dual to merge trees arising naturally over a filtered complex associated to a time series, while horizontal visibility graphs are weak duals of these trees. Immediate consequences include availability of tree-based reconstruction theorems, connections to results on the statistics of self-similar trees, and relations between visibility graphs and the emerging field of applied persistent homology.


title: Automated crater shape retrieval using weakly-supervised deep learning

abstract: Crater ellipticity determination is a complex and time consuming task that so far has evaded successful automation. We train a state of the art computer vision algorithm to identify craters in Lunar digital elevation maps and retrieve their sizes and 2D shapes. The computational backbone of the model is MaskRCNN, an "instance segmentation" general framework that detects craters in an image while simultaneously producing a mask for each crater that traces its outer rim. Our post-processing pipeline then finds the closest fitting ellipse to these masks, allowing us to retrieve the crater ellipticities. Our model is able to correctly identify 87\% of known craters in the longitude range we hid from the network during training and validation (test set), while predicting thousands of additional craters not present in our training data. Manual validation of a subset of these "new" craters indicates that a majority of them are real, which we take as an indicator of the strength of our model in learning to identify craters, despite incomplete training data. The crater size, ellipticity, and depth distributions predicted by our model are consistent with human-generated results. The model allows us to perform a large scale search for differences in crater diameter and shape distributions between the lunar highlands and maria, and we exclude any such differences with a high statistical significance. The predicted test set catalogue and trained model are available here: https://github.com/malidib/Craters_MaskRCNN/.


title: Deformable Slice-to-Volume Registration for Motion Correction in Fetal Body MRI

abstract: In in-utero MRI, motion correction for fetal body and placenta poses a particular challenge due to the presence of local non-rigid transformations of organs caused by bending and stretching. The existing slice-to-volume registration (SVR) reconstruction methods are widely employed for motion correction of fetal brain that undergoes only rigid transformation. However, for reconstruction of fetal body and placenta, rigid registration cannot resolve the issue of misregistrations due to deformable motion, resulting in degradation of features in the reconstructed volume. We propose a Deformable SVR (DSVR), a novel approach for non-rigid motion correction of fetal MRI based on a hierarchical deformable SVR scheme to allow high resolution reconstruction of the fetal body and placenta. Additionally, a robust scheme for structure-based rejection of outliers minimises the impact of registration errors. The improved performance of DSVR in comparison to SVR and patch-to-volume registration (PVR) methods is quantitatively demonstrated in simulated experiments and 20 fetal MRI datasets from 28-31 weeks gestational age (GA) range with varying degree of motion corruption. In addition, we present qualitative evaluation of 100 fetal body cases from 20-34 weeks GA range.


title: Metrics and Motivations for Earth-Space VLBI: Time-Resolving Sgr A* with the Event Horizon Telescope

abstract: Very-long-baseline interferometry (VLBI) at frequencies above 230 GHz with Earth-diameter baselines gives spatial resolution finer than the ${\sim}50 \mu$as "shadow" of the supermassive black hole at the Galactic Center, Sagittarius A* (Sgr A*). Imaging static and dynamical structure near the "shadow" provides a test of general relativity and may allow measurement of black hole parameters. However, traditional Earth-rotation synthesis is inapplicable for sources (such as Sgr A*) with intra-day variability. Expansions of ground-based arrays to include space-VLBI stations may enable imaging capability on time scales comparable to the prograde innermost stable circular orbit (ISCO) of Sgr A*, which is predicted to be 4-30 minutes, depending on black hole spin. We examine the basic requirements for space-VLBI, and we develop tools for simulating observations with orbiting stations. We also develop a metric to quantify the imaging capabilities of an array irrespective of detailed image morphology or reconstruction method. We validate this metric on example reconstructions of simulations of Sgr A* at 230 and 345 GHz, and use these results to motivate expanding the Event Horizon Telescope (EHT) to include small dishes in Low Earth Orbit (LEO). We demonstrate that high-sensitivity sites such as the Atacama Large Millimeter/Submillimeter Array (ALMA) make it viable to add small orbiters to existing ground arrays, as space-ALMA baselines would have sensitivity comparable to ground-based non-ALMA baselines. We show that LEO-enhanced arrays sample half of the diffraction-limited Fourier plane of Sgr A* in less than 30 minutes, enabling reconstructions of near-horizon structure with normalized root-mean-square error $\lesssim0.3$ on sub-ISCO timescales.


title: Data-driven prediction of a multi-scale Lorenz 96 chaotic system using deep learning methods: Reservoir computing, ANN, and RNN-LSTM

abstract: In this paper, the performance of three deep learning methods for predicting short-term evolution and for reproducing the long-term statistics of a multi-scale spatio-temporal Lorenz 96 system is examined. The methods are: echo state network (a type of reservoir computing, RC-ESN), deep feed-forward artificial neural network (ANN), and recurrent neural network with long short-term memory (RNN-LSTM). This Lorenz 96 system has three tiers of nonlinearly interacting variables representing slow/large-scale ($X$), intermediate ($Y$), and fast/small-scale ($Z$) processes. For training or testing, only $X$ is available; $Y$ and $Z$ are never known or used. We show that RC-ESN substantially outperforms ANN and RNN-LSTM for short-term prediction, e.g., accurately forecasting the chaotic trajectories for hundreds of numerical solver's time steps, equivalent to several Lyapunov timescales. The RNN-LSTM and ANN show some prediction skills as well; RNN-LSTM bests ANN. Furthermore, even after losing the trajectory, data predicted by RC-ESN and RNN-LSTM have probability density functions (PDFs) that closely match the true PDF, even at the tails. The PDF of the data predicted using ANN, however, deviates from the true PDF. Implications, caveats, and applications to data-driven and data-assisted surrogate modeling of complex nonlinear dynamical systems such as weather/climate are discussed.


title: Supercontinuum generation in angle-etched diamond waveguides

abstract: We experimentally demonstrate on-chip supercontinuum generation in the visible region in angle etched diamond waveguides. We measure an output spectrum spanning 670 nm to 920 nm in a 5mm long waveguide using 100 fs pulses with 187 pJ of incident pulse energy. Our fabrication technique, combined with diamonds broad transparency window, offers a potential route toward broadband supercontinuum generation in the UV domain.


title: Countably and entropy expansive homeomorphisms with the shadowing property

abstract: We discuss the dynamics beyond topological hyperbolicity considering homeomorphisms satisfying the shadowing property and generalizations of expansivity. It is proved that transitive countably expansive homeomorphisms satisfying the shadowing property are expansive in the set of transitive points. This is in contrast with pseudo-Anosov diffeomorphisms of the two-dimensional sphere that are transitive, cw-expansive, satisfy the shadowing property but the dynamical ball in each transitive point contains a Cantor subset. We exhibit examples of countably expansive homeomorphisms that are not finite expansive, satisfy the shadowing property and admits an infinite number of chain-recurrent classes. We further explore the relation between countable and entropy expansivity and prove that for surface homeomorphisms $f:S\to S$ satisfying the shadowing property and $\Omega(f)=S$, both countably expansive and entropy cw-expansive are equivalent to being topologically conjugate to an Anosov diffeomorphism.


title: A Flexible Pipeline for Prediction of Tropical Cyclone Paths

abstract: Hurricanes and, more generally, tropical cyclones (TCs) are rare, complex natural phenomena of both scientific and public interest. The importance of understanding TCs in a changing climate has increased as recent TCs have had devastating impacts on human lives and communities. Moreover, good prediction and understanding about the complex nature of TCs can mitigate some of these human and property losses. Though TCs have been studied from many different angles, more work is needed from a statistical approach of providing prediction regions. The current state-of-the-art in TC prediction bands comes from the National Hurricane Center of the National Oceanographic and Atmospheric Administration (NOAA), whose proprietary model provides "cones of uncertainty" for TCs through an analysis of historical forecast errors. The contribution of this paper is twofold. We introduce a new pipeline that encourages transparent and adaptable prediction band development by streamlining cyclone track simulation and prediction band generation. We also provide updates to existing models and novel statistical methodologies in both areas of the pipeline, respectively.


title: Sample phase gradient and fringe phase shift in dual phase grating X-ray interferometry

abstract: One of the key tasks in grating based x-ray phase contrast imaging is to accurately retrieve local phase gradients of a sample from measured intensity fringe shifts. To fulfill this task in dual phase grating interferometry, one needs to know the exact mathematical relationship between the two. In this work, using intuitive analysis of the sample-generated fringe shifts based on the beat pattern formation mechanism, the authors derived the formulas relating sample's phase gradients to fringe phase shifts. These formulas provide also a design optimization tool for dual phase grating interferometry.


title: Deep Learning in the Automotive Industry: Recent Advances and Application Examples

abstract: One of the most exciting technology breakthroughs in the last few years has been the rise of deep learning. State-of-the-art deep learning models are being widely deployed in academia and industry, across a variety of areas, from image analysis to natural language processing. These models have grown from fledgling research subjects to mature techniques in real-world use. The increasing scale of data, computational power and the associated algorithmic innovations are the main drivers for the progress we see in this field. These developments also have a huge potential for the automotive industry and therefore the interest in deep learning-based technology is growing. A lot of the product innovations, such as self-driving cars, parking and lane-change assist or safety functions, such as autonomous emergency braking, are powered by deep learning algorithms. Deep learning is poised to offer gains in performance and functionality for most ADAS (Advanced Driver Assistance System) solutions. Virtual sensing for vehicle dynamics application, vehicle inspection/heath monitoring, automated driving and data-driven product development are key areas that are expected to get the most attention. This article provides an overview of the recent advances and some associated challenges in deep learning techniques in the context of automotive applications.


title: Through the Big Bang in inflationary cosmology

abstract: Singularities in General Relativity are regions where the description of spacetime in terms of a pseudo-Riemannian geometry breaks down. The theory seems unable to predict the evolution of the physical degrees of freedom around and beyond such regions. In a recent paper, the author and collaborators challenged this view by providing an example of a singularity at which Einstein's equations can be rewritten in a form that satisfies an existence and uniqueness theorem, thereby predicting that each solution can be continued uniquely through the singularity. This result was obtained under the assumption of homogeneity (but not isotropy), and requires the presence of a massless free scalar field. This paper extends the result to N scalar fields with a potential, the only requirement being that it does not grow too fast. In particular, the result is compatible with inflationary potentials, e.g. Starobinsky's. This brings us one step closer to the goal of extending the original result to realistic cosmologies.


title: An Extensible Framework for Quantifying the Coverage of Defenses Against Untrusted Foundries

abstract: The transistors used to construct Integrated Circuits (ICs) continue to shrink. While this shrinkage improves performance and density, it also reduces trust: the price to build leading-edge fabrication facilities has skyrocketed, forcing even nation states to outsource the fabrication of high-performance ICs. Outsourcing fabrication presents a security threat because the black-box nature of a fabricated IC makes comprehensive inspection infeasible. Since prior work shows the feasibility of fabrication-time attackers' evasion of existing post-fabrication defenses, IC designers must be able to protect their physical designs before handing them off to an untrusted foundry. To this end, recent work suggests methods to harden IC layouts against attack. Unfortunately, no tool exists to assess the effectiveness of the proposed defenses---meaning gaps may exist. This paper presents an extensible IC layout security analysis tool called IC Attack Surface (ICAS) that quantifies defensive coverage. For researchers, ICAS identifies gaps for future defenses to target, and enables the quantitative comparison of existing and future defenses. For practitioners, ICAS enables the exploration of the impact of design decisions on an IC's resilience to fabrication-time attack. ICAS takes a set of metrics that encode the challenge of inserting a hardware Trojan into an IC layout, a set of attacks that the defender cares about, and a completed IC layout and reports the number of ways an attacker can add each attack to the design. While the ideal score is zero, practically, our experience is that lower scores correlate with increased attacker effort.


title: Minimal resolutions of monomial ideals

abstract: An explicit combinatorial minimal free resolution of an arbitrary monomial ideal $I$ in a polynomial ring in $n$ variables over a field of characteristic $0$ is defined canonically, without any choices, using higher-dimensional generalizations of combined spanning trees for cycles and cocycles ("hedges") in the upper Koszul simplicial complexes of $I$ at lattice points in $\mathbb{Z}^n$. The differentials in these "sylvan resolutions" are expressed as matrices whose entries are sums over lattice paths of weights determined combinatorially by sequences of hedges ("hedgerows") along each lattice path. This combinatorics enters via an explicit matroidal expression for the Moore-Penrose pseudoinverses of the differentials in any CW complex as weighted averages of splittings defined by hedges. This "Hedge Formula" also yields a projection formula from CW chains to boundaries. The translation from Moore-Penrose combinatorics to free resolutions relies on Wall complexes, which construct minimal free resolutions of graded ideals from vertical splittings of Koszul bicomplexes. The algebra of Wall complexes applied to individual hedgerows yields explicit but noncanonical combinatorial minimal free resolutions of arbitrary monomial ideals in any characteristic.


title: Defect-implantation for the all-electrical detection of non-collinear spin-textures

abstract: The viability of past, current and future devices for information technology hinges on their sensitivity to the presence of impurities. The latter can lead to resistivity anomalies, the so-called Kondo effect, reshape extrinsically Hall effects or reduce the efficiency of magnetoresistance effects essential in spintronics. Here we demonstrate that atomic defects enable highly efficient all-electrical detection of spin-swirling textures, in particular magnetic skyrmions, which are promising bits candidates in future spintronics devices. Impurities amplify the bare transport signal and can alter significantly the spin-mixing magnetoresistance (XMR) depending on their chemical nature. Both effects are monitored in terms of the defect-enhanced XMR (DXMR) as shown for 3d and 4d transition metal defects implanted at the vicinity of skyrmions generated in PdFe bilayer deposited on Ir(111). The ineluctability of impurities in devices promotes the implementation of DXMR in reading architectures with immediate implications in magnetic storage technologies.


title: Autonomous Navigation of MAVs in Unknown Cluttered Environments

abstract: This paper presents an autonomous navigation framework for reaching a goal in unknown 3D cluttered environments. The framework consists of three main components. First, a computationally efficient method for mapping the environment from the disparity measurements obtained from a depth sensor. Second, a stochastic method to generate a path to a given goal, taking into account field of view constraints on the space that is assumed to be safe for navigation. Third, a fast method for the online generation of motion plans, taking into account the robot's dynamic constraints, model, and environmental uncertainty and disturbances. To highlight the contribution with respect to the available literature, we provide a qualitative and quantitative comparison with the state of the art methods for reaching a goal and for exploration in unknown environments, showing the superior performance of our approach. To illustrate the effectiveness of the proposed framework, we present experiments in multiple indoors and outdoors environments running the algorithm fully on board and in real-time, using a robotic platform based on the Intel Ready to Fly drone kit, which represents the implementation in the most frugal platform for navigation in unknown cluttered environments demonstrated to date. Open source code is available at:~\url{https://github.com/IntelLabs/autonomousmavs}. The video of the experimental results can be found at~\url{https://youtu.be/Wq0e7vF6nZM}.


title: Time-independent approximations for time-dependent optical potentials

abstract: We explore the possibility of modifying the Lewis-Riesenfeld method of invariants developed originally to find exact solutions for time-dependent quantum mechanical systems for the situation in which an exact invariant can be constructed, but the subsequently resulting time-independent eigenvalue system is not solvable exactly. We propose to carry out this step in an approximate fashion, such as employing standard time-independent perturbation theory or the WKB approximation, and subsequently feeding the resulting approximated expressions back into the time-dependent scheme. We illustrate the quality of this approach by contrasting an exactly solvable solution to one obtained with a perturbatively carried out second step for two types of explicitly time-dependent optical potentials.


title: Blowup rate control for solution of Jang's equation and its application on Penrose inequality

abstract: We prove that the blowup term of a blowup solution of Jang's equation on an initial data set (M,g,k) near an arbitrary strictly stable MOTS $ \Sigma $ is exactly $ -\frac{1}{\sqrt{\lambda}}\log \tau $, where $ \tau $ is the distance from $ \Sigma $ and $ \lambda $ is the principal eigenvalue of the MOTS stability operator of $ \Sigma $. We also prove that the gradient of the solution is of order $ \tau^{-1} $. Moreover, we apply these results to get a Penrose-like inequality under additional assumptions.


title: T-TER: Defeating A2 Trojans with Targeted Tamper-Evident Routing

abstract: Since the inception of the Integrated Circuit (IC), the size of the transistors used to construct them has continually shrunk. While this advancement significantly improves computing capability, fabrication costs have skyrocketed. As a result, most IC designers must now outsource fabrication. Outsourcing, however, presents a security threat: comprehensive post-fabrication inspection is infeasible given the size of modern ICs, so it is nearly impossible to know if the foundry has altered the original design during fabrication (i.e., inserted a hardware Trojan). Defending against a foundry-side adversary is challenging because---even with as few as two gates---hardware Trojans can completely undermine software security. Researchers have attempted to both detect and prevent foundry-side attacks, but all existing defenses are ineffective against Trojans with footprints of a few gates or less. We present Targeted Tamper-Evident Routing (T-TER), a preventive layout-level defense against untrusted foundries, capable of thwarting the insertion of even the stealthiest hardware Trojans. T-TER is directed and routing-centric: it prevents foundry-side attackers from routing Trojan wires to, or directly adjacent to, security-critical wires by shielding them with guard wires. Unlike shield wires commonly deployed for cross-talk reduction, T-TER guard wires pose an additional technical challenge: they must be tamper-evident in both the digital (deletion attacks) and analog (move and jog attacks) domains. We address this challenge by developing a class of designed-in guard wires, that are added to the design specifically to protect security-critical wires. T-TER's guard wires incur minimal overhead, scale with design complexity, and provide tamper-evidence against attacks.


title: On Statistical Properties of A Veracity Scoring Method for Spatial Data

abstract: Measuring veracity or reliability of noisy data is of utmost importance, especially in the scenarios where the information are gathered through automated systems. In a recent paper, Chakraborty et. al. (2019) have introduced a veracity scoring technique for geostatistical data. The authors have used a high-quality `reference' data to measure the veracity of the varying-quality observations and incorporated the veracity scores in their analysis of mobile-sensor generated noisy weather data to generate efficient predictions of the ambient temperature process. In this paper, we consider the scenario when no reference data is available and hence, the veracity scores (referred as VS) are defined based on `local' summaries of the observations. We develop a VS-based estimation method for parameters of a spatial regression model. Under a non-stationary noise structure and fairly general assumptions on the underlying spatial process, we show that the VS-based estimators of the regression parameters are consistent. Moreover, we establish the advantage of the VS-based estimators as compared to the ordinary least squares (OLS) estimator by analyzing their asymptotic mean squared errors. We illustrate the merits of the VS-based technique through simulations and apply the methodology to a real data set on mass percentages of ash in coal seams in Pennsylvania.


title: Service Network Design Problem with Capacity-Demand Balancing

abstract: This paper addresses developing cost-effective strategies to respond to excessive demand in the service network design problem in a multi-period setting. The common assumption states that the capacity of freight carriers' assets is capable of handling all of the forecasted demand; however, we assume that there are certain periods such as holiday season in which excessive demand is observed. The demand strictly exceeds the carrier's capacity; even though, the average demand can be still fulfilled throughout the year. In this sense, we let the carrier has three options to respond to the demand: Dispersing the demand with a penalty, leasing additional asset(s) temporarily, and outsourcing some capacity. We propose a modeling and solution approach that jointly incorporates asset management and sizing, outsourcing (3PLs), and earliness/tardiness penalties. The objective is to minimize the overall operational costs by optimally selecting and scheduling the home fleet with respect to 'demand shifting' choices, selecting services from third parties, and routing the commodities on the designed network. We propose an arc-based formulation as well as valid inequalities and present a comprehensive computational study on the randomly generated instances. The formulations with valid inequalities (VIs) outperform the regular formulation in obtaining tighter lower bounds. One set of VIs can improve the CPU time elapsed by 25% on medium-instances that can be solved optimally within the time limit. Furthermore, we develop a custom multi-phase dedicate-merge-and-mix algorithm (DMaM) to solve CSSND problem with an emphasis of obtaining solutions as high-quality as possible practically in a short time in the real world. DMaM has a promising potential to obtain solutions especially for very large instances whereas the commercial solver cannot initialize the B&B algorithm due to excessive memory usage.


title: A minimum entropy principle in the compressible multicomponent Euler equations

abstract: In this work, the space of admissible entropy functions for the compressible multicomponent Euler equations is explored, following up on [Harten, \textit{J. Comput. Phys.}, 49 (1), 1983, pp. 151-164]. This effort allows us to prove a minimum entropy principle on entropy solutions, whether smooth or discrete, in the same way it was originally demonstrated for the compressible Euler equations by [Tadmor, \textit{Appl. Numer. Math.}, 49 (3-5), 1986, pp. 211-219].


title: Octonions, Albert vectors and the group $\mathrm{E}_6(F)$

abstract: We present a uniform approach to the construction of the groups of type $\mathrm{E}_6$ over arbitrary fields without using Lie theory. This gives a simple description of the group generators and some of the subgroup structure. In the finite case our approach also permits relatively straightforward computation of the group order.


title: A Signal Subspace Rotation Method for Localization of Multiple Wideband Sound Sources

abstract: In this paper, the problem of extending narrowband multichannel sound source localization algorithms to the wideband case is addressed. The DOA estimation of narrowband algorithms is based on the estimate of inter-channel phase differences (IPD) between microphones of the sound sources. A new method for wideband sound source DOA estimation based on signal subspace rotation is present. The proposed algorithm normalizes the narrowband signal statistics by rotating the estimated signal subspace to the wideband counterpart in the eigenvector domain. Then the wideband DOA estimate can be obtained by estimating the normalized IPD from these wideband signal statistics. In addition to requiring less computational complexity compared to repeating the narrowband algorithms for all relevant frequencies of wideband signals, the proposed method also does not require any additional prior knowledge. The experimental results demonstrate the efficacy and the robustness of the proposed method.


title: A new, rotating hot corino in Serpens

abstract: We have observed 29 transitions corresponding to 12 distinct species and 7 additional isotopologs toward the deeply embedded Class 0 young stellar object Ser-emb 1 with the Atacama Large Millimeter Array at 1 mm. The detected species include CH3OH and two complex organic molecules, CH3OCH3 and CH3OCHO. The emission of CH3OH and the two COMs is compact, and the CH3OH rotational temperature is 261 (46) K, implying that Ser-emb 1 hosts a hot corino. The derived CH3OH, CH3OCH3 and CH3OCHO column densities are at least 1.2 (0.4)E17 cm-2, 9.2 (3.8)E16 cm-2, and 9.1 (3.6)E16 cm-2, respectively, comparable to the values found for other Class 0 hot corinos. In addition, we observe evidence of rotation at compact scales: two of the more strongly detected lines, corresponding to C18O and H2CO present spatially resolved red- and blue-shifted compact emission orthogonal to the direction of a jet and outflow traced by CO, SiO, and several other molecules. The spatial coincidence of the hot corino emission and a possible disk in a compact region around the central protostar suggests that these structures may be physically and/or chemically related.


title: Improved Planetary Rover Inertial Navigation and Wheel Odometry Performance through Periodic Use of Zero-Type Constraints

abstract: We present an approach to enhance wheeled planetary rover dead-reckoning localization performance by leveraging the use of zero-type constraint equations in the navigation filter. Without external aiding, inertial navigation solutions inherently exhibit cubic error growth. Furthermore, for planetary rovers that are traversing diverse types of terrain, wheel odometry is often unreliable for use in localization, due to wheel slippage. For current Mars rovers, computer vision-based approaches are generally used whenever there is a high possibility of positioning error; however, these strategies require additional computational power, energy resources, and significantly slow down the rover traverse speed. To this end, we propose a navigation approach that compensates for the high likelihood of odometry errors by providing a reliable navigation solution that leverages non-holonomic vehicle constraints as well as state-aware pseudo-measurements (e.g., zero velocity and zero angular rate) updates during periodic stops. By using this, computationally expensive visual-based corrections could be performed less often. Experimental tests that compare against GPS-based localization are used to demonstrate the accuracy of the proposed approach. The source code, post-processing scripts, and example datasets associated with the paper are published in a public repository.


title: Implicitly Adaptive Importance Sampling

abstract: Adaptive importance sampling is a class of techniques for finding good proposal distributions for importance sampling. Often the proposal distributions are standard probability distributions whose parameters are adapted based on the mismatch between the current proposal and a target distribution. In this work, we present an implicit adaptive importance sampling method that applies to complicated distributions which are not available in closed form. The method iteratively matches the moments of a set of Monte Carlo draws to weighted moments based on importance weights. We apply the method to Bayesian leave-one-out cross-validation and show that it performs better than many existing parametric adaptive importance sampling methods while being computationally inexpensive.


title: A Seft-adaptive Multicellular GEP Algorithm Based On Fuzzy Control For Function Optimization

abstract: To improve the global optimization ability of traditional GEP algorithm, a Multicellular gene expression programming algorithm based on fuzzy control (Multicellular GEP Algorithm Based On Fuzzy Control, MGEP-FC) is proposed. The MGEP-FC algorithm describes the size of cross rate, mutation rate and real number mutation rate by constructing fuzzy membership function. According to the concentration and dispersion of individual fitness values in population, the crossover rate, mutation rate and real number set mutation rate of genetic operation are dynamically adjusted. In order to make the diversity of the population continue in the iterative process, a new genetic operation scheme is designed, which combines the new individuals with the parent population to build a temporary population, and the diversity of the temporary and subpopulation are optimized. The results of 12 Benchmark optimization experiments show that the MGEP-FC algorithm has been greatly improved in stability, global convergence and optimization speed.


title: A Hybrid Precipitation Prediction Method based on Multicellular Gene Expression Programming

abstract: Prompt and accurate precipitation forecast is very important for development management of regional water resource, flood disaster prevention and people's daily activity and production plan; however, non-linear and nonstationary characteristics of precipitation data and noise seriously affect forecast accuracy. This paper combines multicellular gene expression programming with more powerful function mining ability and wavelet analysis with more powerful denoising and extracting data fine feature capability for precipitation forecast modeling, proposing to estimate meteorological precipitation with WTGEPRP algorithm. Comparative result for simulation experiment with actual precipitation data in Zhengzhou, Nanning and Melbourne in Australia indicated that: fitting and forecasting performance of WTGEPRP algorithm is better than the algorithm Multicellular Gene Expression Programming-based Hybrid Model for Precipitation Prediction Coupled with EMD, Supporting Vector Regression, BP Neural Network, Multicellular Gene Expression Programming and Gene Expression Programming, and has good application prospect.


title: Hardware-friendly Neural Network Architecture for Neuromorphic Computing

abstract: The hardware-software co-optimization of neural network architectures is becoming a major stream of research especially due to the emergence of commercial neuromorphic chips such as the IBM Truenorth and Intel Loihi. Development of specific neural network architectures in tandem with the design of the neuromorphic hardware considering the hardware constraints will make a huge impact in the complete system level application. In this paper, we study various neural network architectures and propose one that is hardware-friendly for a neuromorphic hardware with crossbar array of synapses. Considering the hardware constraints, we demonstrate how one may design the neuromorphic hardware so as to maximize classification accuracy in the trained network architecture, while concurrently, we choose a neural network architecture so as to maximize utilization in the neuromorphic cores. We also proposed a framework for mapping a neural network onto a neuromorphic chip named as the Mapping and Debugging (MaD) framework. The MaD framework is designed to be generic in the sense that it is a Python wrapper which in principle can be integrated with any simulator tool for neuromorphic chips.


title: Evolving Self-taught Neural Networks: The Baldwin Effect and the Emergence of Intelligence

abstract: The so-called Baldwin Effect generally says how learning, as a form of ontogenetic adaptation, can influence the process of phylogenetic adaptation, or evolution. This idea has also been taken into computation in which evolution and learning are used as computational metaphors, including evolving neural networks. This paper presents a technique called evolving self-taught neural networks - neural networks that can teach themselves without external supervision or reward. The self-taught neural network is intrinsically motivated. Moreover, the self-taught neural network is the product of the interplay between evolution and learning. We simulate a multi-agent system in which neural networks are used to control autonomous agents. These agents have to forage for resources and compete for their own survival. Experimental results show that the interaction between evolution and the ability to teach oneself in self-taught neural networks outperform evolution and self-teaching alone. More specifically, the emergence of an intelligent foraging strategy is also demonstrated through that interaction. Indications for future work on evolving neural networks are also presented.


title: Fuzzy adaptive teaching learning-based optimization strategy for the problem of generating mixed strength t-way test suites

abstract: The teaching learning-based optimization (TLBO) algorithm has shown competitive performance in solving numerous real-world optimization problems. Nevertheless, this algorithm requires better control for exploitation and exploration to prevent premature convergence (i.e., trapped in local optima), as well as enhance solution diversity. Thus, this paper proposes a new TLBO variant based on Mamdani fuzzy inference system, called ATLBO, to permit adaptive selection of its global and local search operations. In order to assess its performances, we adopt ATLBO for the mixed strength t-way test generation problem. Experimental results reveal that ATLBO exhibits competitive performances against the original TLBO and other meta-heuristic counterparts.


title: Learning Longer-term Dependencies via Grouped Distributor Unit

abstract: Learning long-term dependencies still remains difficult for recurrent neural networks (RNNs) despite their success in sequence modeling recently. In this paper, we propose a novel gated RNN structure, which contains only one gate. Hidden states in the proposed grouped distributor unit (GDU) are partitioned into groups. For each group, the proportion of memory to be overwritten in each state transition is limited to a constant and is adaptively distributed to each group member. In other word, every separate group has a fixed overall update rate, yet all units are allowed to have different paces. Information is therefore forced to be latched in a flexible way, which helps the model to capture long-term dependencies in data. Besides having a simpler structure, GDU is demonstrated experimentally to outperform LSTM and GRU on tasks including both pathological problems and natural data set.


title: Deep Neuroevolution of Recurrent and Discrete World Models

abstract: Neural architectures inspired by our own human cognitive system, such as the recently introduced world models, have been shown to outperform traditional deep reinforcement learning (RL) methods in a variety of different domains. Instead of the relatively simple architectures employed in most RL experiments, world models rely on multiple different neural components that are responsible for visual information processing, memory, and decision-making. However, so far the components of these models have to be trained separately and through a variety of specialized training methods. This paper demonstrates the surprising finding that models with the same precise parts can be instead efficiently trained end-to-end through a genetic algorithm (GA), reaching a comparable performance to the original world model by solving a challenging car racing task. An analysis of the evolved visual and memory system indicates that they include a similar effective representation to the system trained through gradient descent. Additionally, in contrast to gradient descent methods that struggle with discrete variables, GAs also work directly with such representations, opening up opportunities for classical planning in latent space. This paper adds additional evidence on the effectiveness of deep neuroevolution for tasks that require the intricate orchestration of multiple components in complex heterogeneous architectures.


title: One-vs-All Models for Asynchronous Training: An Empirical Analysis

abstract: Any given classification problem can be modeled using multi-class or One-vs-All (OVA) architecture. An OVA system consists of as many OVA models as the number of classes, providing the advantage of asynchrony, where each OVA model can be re-trained independent of other models. This is particularly advantageous in settings where scalable model training is a consideration (for instance in an industrial environment where multiple and frequent updates need to be made to the classification system). In this paper, we conduct empirical analysis on realizing independent updates to OVA models and its impact on the accuracy of the overall OVA system. Given that asynchronous updates lead to differences in training datasets for OVA models, we first define a metric to quantify the differences in datasets. Thereafter, using Natural Language Understanding as a task of interest, we estimate the impact of three factors: (i) number of classes, (ii) number of data points and, (iii) divergences in training datasets across OVA models; on the OVA system accuracy. Finally, we observe the accuracy impact of increased asynchrony in a Spoken Language Understanding system. We analyze the results and establish that the proposed metric correlates strongly with the model performances in both the experimental settings.


title: Closing the Accuracy Gap in an Event-Based Visual Recognition Task

abstract: Mobile and embedded applications require neural networks-based pattern recognition systems to perform well under a tight computational budget. In contrast to commonly used synchronous, frame-based vision systems and CNNs, asynchronous, spiking neural networks driven by event-based visual input respond with low latency to sparse, salient features in the input, leading to high efficiency at run-time. The discrete nature of the event-based data streams makes direct training of asynchronous neural networks challenging. This paper studies asynchronous spiking neural networks, obtained by conversion from a conventional CNN trained on frame-based data. As an example, we consider a CNN trained to steer a robot to follow a moving target. We identify possible pitfalls of the conversion and demonstrate how the proposed solutions bring the classification accuracy of the asynchronous network to only 3\% below the performance of the original synchronous CNN, while requiring 12x fewer computations. While being applied to a simple task, this work is an important step towards low-power, fast, and embedded neural networks-based vision solutions for robotic applications.


title: Techniques for Inferring Context-Free Lindenmayer Systems With Genetic Algorithm

abstract: Lindenmayer systems (L-systems) are a formal grammar system, where the most notable feature is a set of rewriting rules that are used to replace every symbol in a string in parallel; by repeating this process, a sequence of strings is produced. Some symbols in the strings may be interpreted as instructions for simulation software. Thus, the sequence can be used to model the steps of a process. Currently, creating an L-system for a specific process is done by hand by experts through much effort. The inductive inference problem attempts to infer an L-system from such a sequence of strings generated by an unknown system; this can be thought of as an intermediate step to inferring from a sequence of images. This paper evaluates and analyzes different genetic algorithm encoding schemes and mathematical properties for the L-system inductive inference problem. A new tool, the Plant Model Inference Tool for Context-Free L-systems (PMIT-D0L) is implemented based on these techniques. PMIT-D0L has been successfully evaluated on 28 known L-systems, with alphabets up to 31 symbols and a total sum of 281 symbols across the rewriting rules. PMIT-D0L can infer even the largest of these L-systems in less than a few seconds.


title: Synthesizing Images from Spatio-Temporal Representations using Spike-based Backpropagation

abstract: Spiking neural networks (SNNs) offer a promising alternative to current artificial neural networks to enable low-power event-driven neuromorphic hardware. Spike-based neuromorphic applications require processing and extracting meaningful information from spatio-temporal data, represented as series of spike trains over time. In this paper, we propose a method to synthesize images from multiple modalities in a spike-based environment. We use spiking auto-encoders to convert image and audio inputs into compact spatio-temporal representations that is then decoded for image synthesis. For this, we use a direct training algorithm that computes loss on the membrane potential of the output layer and back-propagates it by using a sigmoid approximation of the neuron's activation function to enable differentiability. The spiking autoencoders are benchmarked on MNIST and Fashion-MNIST and achieve very low reconstruction loss, comparable to ANNs. Then, spiking autoencoders are trained to learn meaningful spatio-temporal representations of the data, across the two modalities - audio and visual. We synthesize images from audio in a spike-based environment by first generating, and then utilizing such shared multi-modal spatio-temporal representations. Our audio to image synthesis model is tested on the task of converting TI-46 digits audio samples to MNIST images. We are able to synthesize images with high fidelity and the model achieves competitive performance against ANNs.


title: Neural Stored-program Memory

abstract: Neural networks powered with external memory simulate computer behaviors. These models, which use the memory to store data for a neural controller, can learn algorithms and other complex tasks. In this paper, we introduce a new memory to store weights for the controller, analogous to the stored-program memory in modern computer architectures. The proposed model, dubbed Neural Stored-program Memory, augments current memory-augmented neural networks, creating differentiable machines that can switch programs through time, adapt to variable contexts and thus resemble the Universal Turing Machine. A wide range of experiments demonstrate that the resulting machines not only excel in classical algorithmic problems, but also have potential for compositional, continual, few-shot learning and question-answering tasks.


title: Particle swarm optimization model to predict scour depth around bridge pier

abstract: Scour depth around bridge piers plays a vital role in the safety and stability of the bridges. Existing methods to predict scour depth are mainly based on regression models or black box models in which the first one lacks enough accuracy while the later one does not provide a clear mathematical expression to easily employ it for other situations or cases. Therefore, this paper aims to develop new equations using particle swarm optimization as a metaheuristic approach to predict scour depth around bridge piers. To improve the efficiency of the proposed model, individual equations are derived for laboratory and field data. Moreover, sensitivity analysis is conducted to achieve the most effective parameters in the estimation of scour depth for both experimental and filed data sets. Comparing the results of the proposed model with those of existing regression-based equations reveal the superiority of the proposed method in terms of accuracy and uncertainty. Moreover, the ratio of pier width to flow depth and ratio of d50 (mean particle diameter) to flow depth for the laboratory and field data were recognized as the most effective parameters, respectively. The derived equations can be used as a suitable proxy to estimate scour depth in both experimental and prototype scales.


title: Accurate and Energy-Efficient Classification with Spiking Random Neural Network: Corrected and Expanded Version

abstract: Artificial Neural Network (ANN) based techniques have dominated state-of-the-art results in most problems related to computer vision, audio recognition, and natural language processing in the past few years, resulting in strong industrial adoption from all leading technology companies worldwide. One of the major obstacles that have historically delayed large scale adoption of ANNs is the huge computational and power costs associated with training and testing (deploying) them. In the mean-time, Neuromorphic Computing platforms have recently achieved remarkable performance running more bio-realistic Spiking Neural Networks at high throughput and very low power consumption making them a natural alternative to ANNs. Here, we propose using the Random Neural Network (RNN), a spiking neural network with both theoretical and practical appealing properties, as a general purpose classifier that can match the classification power of ANNs on a number of tasks while enjoying all the features of a spiking neural network. This is demonstrated on a number of real-world classification datasets.


title: Evolving Self-supervised Neural Networks: Autonomous Intelligence from Evolved Self-teaching

abstract: This paper presents a technique called evolving self-supervised neural networks - neural networks that can teach themselves, intrinsically motivated, without external supervision or reward. The proposed method presents some sort-of paradigm shift, and differs greatly from both traditional gradient-based learning and evolutionary algorithms in that it combines the metaphor of evolution and learning, more specifically self-learning, together, rather than treating these phenomena alternatively. I simulate a multi-agent system in which neural networks are used to control autonomous foraging agents with little domain knowledge. Experimental results show that only evolved self-supervised agents can demonstrate some sort of intelligent behaviour, but not evolution or self-learning alone. Indications for future work on evolving intelligence are also presented.


title: Towards Efficient Neural Networks On-a-chip: Joint Hardware-Algorithm Approaches

abstract: Machine learning algorithms have made significant advances in many applications. However, their hardware implementation on the state-of-the-art platforms still faces several challenges and are limited by various factors, such as memory volume, memory bandwidth and interconnection overhead. The adoption of the crossbar architecture with emerging memory technology partially solves the problem but induces process variation and other concerns. In this paper, we will present novel solutions to two fundamental issues in crossbar implementation of Artificial Intelligence (AI) algorithms: device variation and insufficient interconnections. These solutions are inspired by the statistical properties of algorithms themselves, especially the redundancy in neural network nodes and connections. By Random Sparse Adaptation and pruning the connections following the Small-World model, we demonstrate robust and efficient performance on representative datasets such as MNIST and CIFAR-10. Moreover, we present Continuous Growth and Pruning algorithm for future learning and adaptation on hardware.


title: Self-adaptive Potential-based Stopping Criteria for Particle Swarm Optimization

abstract: We study the variant of Particle Swarm Optimization (PSO) that applies random velocities in a dimension instead of the regular velocity update equations as soon as the so-called potential of the swarm falls below a certain bound in this dimension, arbitrarily set by the user. In this case, the swarm performs a forced move. In this paper, we are interested in how, by counting the forced moves, the swarm can decide for itself to stop its movement because it is improbable to find better solution candidates as it already has found. We formally prove that when the swarm is close to a (local) optimum, it behaves like a blind-searching cloud, and that the frequency of forced moves exceeds a certain, objective function-independent value. Based on this observation, we define stopping criteria and evaluate them experimentally showing that good solution candidates can be found much faster than applying other criteria.