title: Finding a boundary between valid and invalid regions of the input space

abstract: In the context of robustness testing, the boundary between the valid and invalid regions of the input space can be an interesting source of erroneous inputs. Knowing where a specific software under test (SUT) has a boundary is essential for validation in relation to requirements. However, finding where a SUT actually implements the boundary is a non-trivial problem that has not gotten much attention. This paper proposes a method of finding the boundary between the valid and invalid regions of the input space. The proposed method consists of two steps. First, test data generators, directed by a search algorithm to maximise distance to known, valid test cases, generate valid test cases that are closer to the boundary. Second, these valid test cases undergo mutations to try to push them over the boundary and into the invalid part of the input space. This results in a pair of test sets, one consisting of test cases on the valid side of the boundary and a matched set on the outer side, with only a small distance between the two sets. The method is evaluated on a number of examples from the standard library of a modern programming language. We propose a method of determining the boundary between valid and invalid regions of the input space and apply it on a SUT that has a non-contiguous valid region of the input space. From the small distance between the developed pairs of test sets, and the fact that one test set contains valid test cases and the other invalid test cases, we conclude that the pair of test sets described the boundary between the valid and invalid regions of that input space. Differences of behaviour can be observed between different distances and sets of mutation operators, but all show that the method is able to identify the boundary between the valid and invalid regions of the input space. This is an important step towards more automated robustness testing.


title: Optimizing Agent Behavior over Long Time Scales by Transporting Value

abstract: Humans spend a remarkable fraction of waking life engaged in acts of "mental time travel". We dwell on our actions in the past and experience satisfaction or regret. More than merely autobiographical storytelling, we use these event recollections to change how we will act in similar scenarios in the future. This process endows us with a computationally important ability to link actions and consequences across long spans of time, which figures prominently in addressing the problem of long-term temporal credit assignment; in artificial intelligence (AI) this is the question of how to evaluate the utility of the actions within a long-duration behavioral sequence leading to success or failure in a task. Existing approaches to shorter-term credit assignment in AI cannot solve tasks with long delays between actions and consequences. Here, we introduce a new paradigm for reinforcement learning where agents use recall of specific memories to credit actions from the past, allowing them to solve problems that are intractable for existing algorithms. This paradigm broadens the scope of problems that can be investigated in AI and offers a mechanistic account of behaviors that may inspire computational models in neuroscience, psychology, and behavioral economics.


title: Borel reducibility and symmetric models

abstract: We develop a correspondence between the study of Borel equivalence relations induced by closed subgroups of $S_\infty$, and the study of symmetric models and weak choice principles, and apply it to prove a conjecture of Hjorth-Kechris-Louveau (1998). For example, we show that the equivalence relation $\cong^\ast_{\omega+1,0}$ is strictly below $\cong^\ast_{\omega+1,<\omega}$ in Borel reducibility. By results of Hjorth-Kechris-Louveau, $\cong^\ast_{\omega+1,<\omega}$ provides invariants for $\Sigma^0_{\omega+1}$ equivalence relations induced by actions of $S_\infty$, while $\cong^\ast_{\omega+1,0}$ provides invariants for $\Sigma^0_{\omega+1}$ equivalence relations induced by actions of abelian closed subgroups of $S_\infty$. We further apply these techniques to study the Friedman-Stanley jumps. For example, we find an equivalence relation $F$, Borel bireducible with $=^{++}$, so that $F\restriction C$ is not Borel reducible to $=^{+}$ for any non-meager set $C$. This answers a question of Zapletal, arising from the results of Kanovei-Sabok-Zapletal (2013). For these proofs we analyze the symmetric models $M_n$, $n<\omega$, developed by Monro (1973), and extend the construction past $\omega$, through all countable ordinals. This answers a question of Karagila (2019).


title: Robust flow field reconstruction from limited measurements via sparse representation

abstract: In many applications it is important to estimate a fluid flow field from limited and possibly corrupt measurements. Current methods in flow estimation often use least squares regression to reconstruct the flow field, finding the minimum-energy solution that is consistent with the measured data. However, this approach may be prone to overfitting and sensitive to noise. To address these challenges we instead seek a sparse representation of the data in a library of examples. Sparse representation has been widely used for image recognition and reconstruction, and it is well-suited to structured data with limited, corrupt measurements. We explore sparse representation for flow reconstruction on a variety of fluid data sets with a wide range of complexity, including vortex shedding past a cylinder at low Reynolds number, a mixing layer, and two geophysical flows. In addition, we compare several measurement strategies and consider various types of noise and corruption over a range of intensities. We find that sparse representation has considerably improved estimation accuracy and robustness to noise and corruption compared with least squares methods. We also introduce a sparse estimation procedure on local spatial patches for complex multiscale flows that preclude a global sparse representation. Based on these results, sparse representation is a promising framework for extracting useful information from complex flow fields with realistic measurements.


title: Science with an ngVLA: Organics in Disk Midplanes with the ngVLA

abstract: Planets assemble in the midplanes of protoplanetary disks. The compositions of dust and gas in the disk midplane region determine the compositions of nascent planets, including their chemical hospitality to life. In this context, the distributions of volatile organic material across the planet and comet forming zones is of special interest. These are difficult to access in the disk midplane at IR and even millimeter wavelengths due to dust opacity, which can veil the midplane, low intrinsic molecular abundances due to efficient freeze-out, and, in the case of mid-sized organics, a mismatch between expected excitation temperatures and accessible line upper energy levels. At ngVLA wavelengths, the dust is optically thin, enabling observations into the planet forming disk midplane. ngVLA also has the requisite sensitivity. Using TW Hya as a case study, we show that ngVLA will be able to map out the distributions of diagnostic organics, such as CH3CN, in nearby protoplanetary disks.


title: On Improving Service Chains Survivability through Efficient Backup Provisioning

abstract: With the growing adoption of Software Defined Networking (SDN) and Network Function Virtualization (NFV), large-scale NFV infrastructure deployments are gaining momentum. Such infrastructures are home to thousands of network Service Function Chains (SFCs), each composed of a chain of virtual network functions (VNFs) that are processing incoming traffic flows. Unfortunately, in such environments, the failure of a single node may break down several VNFs and thereby breaking many service chains at the same time. In this paper, we address this particular problem and investigate possible solutions to ensure the survivability of the affected service chains by provisioning backup VNFs that can take over in case of failure. Specifically, we propose a survivability management framework to efficiently manage SFCs and the backup VNFs. We formulate the SFC survivability problem as an integer linear program that determines the minimum number of required backups to protect all the SFCs in the system and identifies their optimal placement in the infrastructure. We also propose two heuristic algorithms to cope with the large-scale instances of the problem. Through extensive simulations of different deployment scenarios, we show that these algorithms provide near-optimal solutions with minimal computation time.


title: Geometric quantization of localized surface plasmons

abstract: We consider the quasi-static problem governing the localized surface plasmon modes and permittivity eigenvalues $\epsilon$ of smooth, arbitrarily shaped, axisymmetric inclusions. We develop an asymptotic theory for the dense part of the spectrum, i.e., close to the accumulation value $\epsilon=-1$ at which a flat interface supports surface plasmons; in this regime, the field oscillates rapidly along the surface and decays exponentially away from it on a comparable scale. With $\tau=-(\epsilon+1)$ as the small parameter, we develop a surface-ray description of the eigenfunctions in a narrow boundary layer about the interface; the fast phase variation, as well as the slowly varying amplitude and geometric phase, along the rays are determined as functions of the local geometry. We focus on modes varying at most moderately in the azimuthal direction, in which case the surface rays are meridian arcs that focus at the two poles. Asymptotically matching the diverging ray solutions with expansions valid in inner regions in the vicinities of the poles yields the quantization rule $$\frac{1}{\tau} \sim \frac{\pi n }{\Theta}+\frac{1}{2}\left(\frac{\pi}{\Theta}-1\right)+o(1),$$ where $n\gg1$ is an integer and $\Theta$ a geometric parameter given by the product of the inclusion length and the reciprocal average of its cross-sectional radius along its symmetry axis. For a sphere, $\Theta=\pi$, whereby the formula returns the exact eigenvalues $\epsilon_n=-1-1/n$. We also demonstrate good agreement with exact solutions in the case of prolate spheroids.


title: Science with an ngVLA: Accretion and Jets in Local Compact Objects

abstract: Despite the prevalence of jets in accreting systems and their impact on the surrounding medium, the fundamental physics of how they are launched and collimated is not fully understood. Radio observations of local compact objects, including accreting stellar mass black holes, neutron stars and white dwarfs, probe their jet emission. Coupled with multi-wavelength observations, this allows us to test the underlying accretion-outflow connection and to establish the relationship between the accretor properties and the jet power, which is necessary to accurately model jets. Compact accretors are nearby, numerous and come in a range of accretor properties, and hence are ideal probes for the underlying jet physics. Despite this there are a number of key outstanding questions regarding accretion-driven outflows in these objects that cannot be answered with current radio observations. The vastly improved sensitivity, polarization capabilities, spatial resolution and high-frequency coverage of the ngVLA will be crucial to answering these, and subsequently determining the fundamental physics behind accretion and jets at all physical scales.


title: Forward modelling of brightness variations in Sun-like stars I. Emergence and surface transport of magnetic flux

abstract: The latitudinal distribution of starspots deviates from the solar pattern with increasing rotation rate. Numerical simulations of magnetic flux emergence and transport can help model the observed stellar activity patterns and the associated brightness variations. We set up a composite model for the processes of flux emergence and transport on Sun-like stars, to simulate stellar brightness variations for various levels of magnetic activity and rotation rates. Assuming that the distribution of magnetic flux at the base of the convection zone follows solar scaling relations, we calculate the emergence latitudes and tilt angles of bipolar regions at the surface for various rotation rates, using thin-flux-tube simulations. Taking these two quantities as input to a surface flux transport SFT model, we simulate the diffusive-advective evolution of the radial field at the stellar surface, including effects of active region nesting. As the rotation rate increases, (1) magnetic flux emerges at higher latitudes and an inactive gap opens around the equator, reaching a half-width of $20^\circ$ for $8\Omega_\odot$, (2) the tilt angles of freshly emerged bipolar regions show stronger variations with latitude. Polar spots can form at $8\Omega_\odot$ by accumulation of follower-polarity flux from decaying bipolar regions. From $4\Omega_\odot$ to $8\Omega_\odot$, the maximum spot coverage changes from 3 to 20%, respectively, compared to 0.4% for the solar model. Nesting of activity can lead to strongly non-axisymmetric spot distributions. On Sun-like stars rotating at $8\Omega_\odot$ ($P_{\rm rot}\simeq 3$ days), polar spots can form, owing to higher levels of flux emergence rate and tilt angles. Defining spots by a threshold field strength yields global spot coverages that are roughly consistent with stellar observations.


title: Robust Neural Machine Translation with Joint Textual and Phonetic Embedding

abstract: Neural machine translation (NMT) is notoriously sensitive to noises, but noises are almost inevitable in practice. One special kind of noise is the homophone noise, where words are replaced by other words with similar pronunciations. We propose to improve the robustness of NMT to homophone noises by 1) jointly embedding both textual and phonetic information of source sentences, and 2) augmenting the training dataset with homophone noises. Interestingly, to achieve better translation quality and more robustness, we found that most (though not all) weights should be put on the phonetic rather than textual information. Experiments show that our method not only significantly improves the robustness of NMT to homophone noises, but also surprisingly improves the translation quality on some clean test sets.


title: SPRT Based Transceiver for Molecular Communications

abstract: Achieving precise synchronisation between transmitters and receivers is particularly challenging in diffusive molecular communication environments. To this end, point-to-point molecular communication system design is examined wherein synchronisation errors are explicitly considered. Two transceiver design questions are considered: the development of a sequential probability ratio test-based detector which allows for additional observations in the presence of uncertainty due to mis-synchronisation at the receiver, and a modulation design which is optimised for this receiver strategy. The modulation is based on optimising an approximation for the probability of error for the detection strategy and directly exploits the structure of the probability of molecules hitting a receiver within a particular time slot. The proposed receiver and modulation designs achieve strongly improved asynchronous detection performance for the same data rate as a decision feedback based receiver by a factor of 1/2.


title: Stochastic Resonance in Neural Network, Noise Color Effects

abstract: Some systems cannot be predicted by classical theories and it is required the development of combined deterministic and stochastic theories that make used of noise for dynamical prediction. Noise is not always an interfering signal which perturbs the system. On the contrary, noise signal can enhance the performance of some non-linear systems. The advantage of noise is observed in Stochastic Resonance (SR) where noise is used for amplification and subsequent detection of small signals. To detect this phenomena it is necessary that the system have a bistable potential barrier that creates a threshold, the input of the system should be a weak periodic signal which amplitude is below threshold together with an stochastic signal. In this way, the SR is detected when there are weak periodic signals that are added to different noise colors in order to be amplified and optimised. The interactions between the two signals transform the potential of the system precisely at the frequency of the weak periodic signal that is added to the system. The behaviour of the SR is detected in a neural network and it is study under noise color variations. Here I found that Pink noise amplified the sub-threshold input signal twenty times more in comparison to white noise. This is evidence of the functionality of background noise in the brain, where neurons are naturally embedded in pink noise.


title: On the small-scale structure of turbulence and its impact on the pressure field

abstract: Understanding the small-scale structure of incompressible turbulence and its implications for the non-local pressure field is one of the fundamental challenges in fluid mechanics. Intense velocity gradient structures tend to cluster on a range of scales which affects the pressure through a Poisson equation. Here we present a quantitative investigation of the spatial distribution of these structures conditional on their intensity for Taylor-based Reynolds numbers in the range [160, 380]. We find that the correlation length, the second invariant of the velocity gradient, is proportional to the Kolmogorov scale. It also is a good indicator for the spatial localization of intense enstrophy and strain-dominated regions, as well as the separation between them. We describe and quantify the differences in the two-point statistics of these regions and the impact they have on the non-locality of the pressure field as a function of the intensity of the regions. Specifically, across the examined range of Reynolds numbers, the pressure in strong rotation-dominated regions is governed by a dissipation-scale neighbourhood. In strong strain-dominated regions, on the other hand, it is determined primarily by a larger neighbourhood reaching inertial scales.


title: The milliQan Experiment: Search for milli-charged Particles at the LHC

abstract: Recently, a search for milli-charged particles produced at the LHC has been proposed. The experiment, named milliQan, is expected to obtain sensitivity to charges of $10^{- 1} - 10^{-3}e$ for masses in the 0.1 - 100 GeV range. The detector is composed of 3 stacks of 80 cm long plastic scintillator arrays read out by PMTs. It will be installed in an existing tunnel 33 m from the CMS interaction point at the LHC, with 17 m of rock shielding to suppress beam backgrounds. In the fall of 2017 a 1% scale "demonstrator" of the proposed detector was installed at the planned site in order to study the feasibility of the experiment, focusing on understanding various background sources such as radioactivity of materials, PMT dark current, cosmic rays, and beam induced backgrounds. The data from the demonstrator provides a unique opportunity to understand the backgrounds and to optimize the design of the full detector.


title: Ridge Resonance: A new Resonance Phenomenon for Silicon Photonics Harnessing Bound States in the Continuum

abstract: We present a new resonant behavior based on bound states in the continuum in a guided wave silicon platform. The continuum has the form of a collimated beam of light which is confined vertically in a TE mode of a silicon slab. The bound state is a discrete TM mode of a ridge on the silicon slab. The coupling between the slab and ridge modes results in a single sharp resonance at the wavelength where they phase match. We experimentally demonstrate this phenomenon on a silicon photonic chip using foundry compatible parameters and interface it on-chip to standard single mode silicon nanowire waveguides. The fabricated chip exhibits a single sharp resonance near 1550 nm with a line width of a few nanometer, an extinction ratio of 25 dB and a thermal stability of 19.5 pm/C. We believe that this is the first demonstration of bound states in the continuum resonance realized using guided wave components.


title: Moduli spaces of sheaves on K3 surfaces and Galois representations

abstract: We consider two K3 surfaces defined over an arbitrary field, together with a smooth proper moduli space of stable sheaves on each. When the moduli spaces have the same dimension, we prove that if the \'etale cohomology groups (with Q_ell coefficients) of the two surfaces are isomorphic as Galois representations, then the same is true of the two moduli spaces. In particular, if the field of definition is finite and the K3 surfaces have equal zeta functions, then so do the moduli spaces, even when the moduli spaces are not birational.


title: Laser stimulated deexcitation of Rydberg antihydrogen atoms

abstract: Antihydrogen atoms are routinely formed at CERN in a broad range of Rydberg states. Ground-state anti-atoms, those useful for precision measurements, are eventually produced through spontaneous decay. However given the long lifetime of Rydberg states the number of ground-state antihydrogen atoms usable is small, in particular for experiments relying on the production of a beam of antihydrogen atoms. Therefore, it is of high interest to efficiently stimulate the decay in order to retain a higher fraction of ground-state atoms for measurements. We propose a method that optimally mixes the high angular momentum states with low ones enabling to stimulate, using a broadband frequency laser, the deexcitation toward low-lying states, which then spontaneously decay to ground-state. We evaluated the method in realistic antihydrogen experimental conditions. For instance, starting with an initial distribution of atoms within the $n=20-30$ manifolds, as formed through charge exchange mechanism, we show that more than 80\% of antihydrogen atoms will be deexcited to the ground-state within 100 ns using a laser producing 2 J at 828 nm.


title: Science with an ngVLA: How Do Cold Gas Outflows Shape Galaxies?

abstract: The ngVLA will obtain breakthrough observations of the cold phases (molecular and atomic) of galactic outflows. These observations will be key to study the driving mechanisms of cold winds, the fate of the gas, the mass-loss rates, the conditions that trigger cold outflows, the fraction of gas that escapes galaxies, the effects of winds on suppressing star formation and slowing black hole growth, the statistical properties of outflows, and the relation between large scale outflows and UFOs and BAL outflows. Finally, we highlight the ngVLA complementarity with existing and planned facilities, such as ALMA/SKA, OST/SPICA, and Lynx/ATHENA.


title: Random clique covers for graphs with local density and global sparsity

abstract: Large real-world graphs tend to be sparse, but they often contain many densely connected subgraphs and exhibit high clustering coefficients. While recent random graph models can capture this sparsity, they ignore the local density, or vice versa. We develop a Bayesian nonparametric graph model based on random edge clique covers, and show that this model can capture power law degree distribution, global sparsity and non-vanishing local clustering coefficient. This distribution can be used directly as a prior on observed graphs, or as part of a hierarchical Bayesian model for inferring latent graph structures.


title: Geometric stabilisation via p-adic integration

abstract: In this article we give a new proof of Ng\^o's Geometric Stabilisation Theorem, which implies the Fundamental Lemma. This is a statement which relates the cohomology of Hitchin fibres for a quasi-split reductive group scheme $G$ to the cohomology of Hitchin fibres for the endoscopy groups $H_{\kappa}$. Our proof avoids the Decomposition and Support Theorem, instead the argument is based on results for $p$-adic integration on coarse moduli spaces of Deligne-Mumford stacks. Along the way we establish a description of the inertia stack of the (anisotropic) moduli stack of $G$-Higgs bundles in terms of endoscopic data, and extend duality for generic Hitchin fibres of Langlands dual group schemes to the quasi-split case.


title: An Illuminating Algorithm for the Light Bulb Problem

abstract: The Light Bulb Problem is one of the most basic problems in data analysis. One is given as input $n$ vectors in $\{-1,1\}^d$, which are all independently and uniformly random, except for a planted pair of vectors with inner product at least $\rho \cdot d$ for some constant $\rho > 0$. The task is to find the planted pair. The most straightforward algorithm leads to a runtime of $\Omega(n^2)$. Algorithms based on techniques like Locality-Sensitive Hashing achieve runtimes of $n^{2 - O(\rho)}$; as $\rho$ gets small, these approach quadratic. Building on prior work, we give a new algorithm for this problem which runs in time $O(n^{1.582} + nd),$ regardless of how small $\rho$ is. This matches the best known runtime due to Karppa et al. Our algorithm combines techniques from previous work on the Light Bulb Problem with the so-called `polynomial method in algorithm design,' and has a simpler analysis than previous work. Our algorithm is also easily derandomized, leading to a deterministic algorithm for the Light Bulb Problem with the same runtime of $O(n^{1.582} + nd),$ improving previous results.


title: Achieving Covert Wireless Communications Using a Full-Duplex Receiver

abstract: Covert communications hide the transmission of a message from a watchful adversary while ensuring a certain decoding performance at the receiver. In this work, a wireless communication system under fading channels is considered where covertness is achieved by using a full-duplex (FD) receiver. More precisely, the receiver of covert information generates artificial noise with a varying power causing uncertainty at the adversary, Willie, regarding the statistics of the received signals. Given that Willie's optimal detector is a threshold test on the received power, we derive a closed-form expression for the optimal detection performance of Willie averaged over the fading channel realizations. Furthermore, we provide guidelines for the optimal choice of artificial noise power range, and the optimal transmission probability of covert information to maximize the detection errors at Willie. Our analysis shows that the transmission of artificial noise, although causes self-interference, provides the opportunity of achieving covertness but its transmit power levels need to be managed carefully. We also demonstrate that the prior transmission probability of 0.5 is not always the best choice for achieving the maximum possible covertness, when the covert transmission probability and artificial noise power can be jointly optimized.


title: Assessing and Remedying Coverage for a Given Dataset

abstract: Data analysis impacts virtually every aspect of our society today. Often, this analysis is performed on an existing dataset, possibly collected through a process that the data scientists had limited control over. The existing data analyzed may not include the complete universe, but it is expected to cover the diversity of items in the universe. Lack of adequate coverage in the dataset can result in undesirable outcomes such as biased decisions and algorithmic racism, as well as creating vulnerabilities such as opening up room for adversarial attacks. In this paper, we assess the coverage of a given dataset over multiple categorical attributes. We first provide efficient techniques for traversing the combinatorial explosion of value combinations to identify any regions of attribute space not adequately covered by the data. Then, we determine the least amount of additional data that must be obtained to resolve this lack of adequate coverage. We confirm the value of our proposal through both theoretical analyses and comprehensive experiments on real data.


title: Marrying Universal Dependencies and Universal Morphology

abstract: The Universal Dependencies (UD) and Universal Morphology (UniMorph) projects each present schemata for annotating the morphosyntactic details of language. Each project also provides corpora of annotated text in many languages - UD at the token level and UniMorph at the type level. As each corpus is built by different annotators, language-specific decisions hinder the goal of universal schemata. With compatibility of tags, each project's annotations could be used to validate the other's. Additionally, the availability of both type- and token-level resources would be a boon to tasks such as parsing and homograph disambiguation. To ease this interoperability, we present a deterministic mapping from Universal Dependencies v2 features into the UniMorph schema. We validate our approach by lookup in the UniMorph corpora and find a macro-average of 64.13% recall. We also note incompatibilities due to paucity of data on either side. Finally, we present a critical evaluation of the foundations, strengths, and weaknesses of the two annotation projects.


title: Gravitational redshift in quantum-clock interferometry

abstract: The creation of delocalized coherent superpositions of quantum systems experiencing different relativistic effects is an important milestone in future research at the interface of gravity and quantum mechanics. This could be achieved by generating a superposition of quantum clocks that follow paths with different gravitational time dilation and investigating the consequences on the interference signal when they are eventually recombined. Light-pulse atom interferometry with elements employed in optical atomic clocks is a promising candidate for that purpose, but suffers from major challenges including its insensitivity to the gravitational redshift in a uniform field. All these difficulties can be overcome with a novel scheme presented here which is based on initializing the clock when the spatially separate superposition has already been generated and performing a doubly differential measurement where the differential phase shift between the two internal states is compared for different initialization times. This can be exploited to test the universality of the gravitational redshift with delocalized coherent superpositions of quantum clocks and it is argued that its experimental implementation should be feasible with a new generation of 10-meter atomic fountains that will soon become available. Interestingly, the approach also offers significant advantages for more compact set-ups based on guided interferometry or hybrid configurations. Furthermore, in order to provide a solid foundation for the analysis of the various interferometry schemes and the effects that can be measured with them, a general formalism for a relativistic description of atom interferometry in curved spacetime is developed. It can deal with freely falling atoms, but also include the effects of external forces and guiding potentials, and can be applied to a very wide range of situations.


title: Can Euroscepticism Contribute to a European Public Sphere? The Europeanization of Media Discourses about Euroscepticism across Six Countries

abstract: This study compares the media discourses about Euroscepticism in 2014 across six countries (United Kingdom, Ireland, France, Spain, Sweden, and Denmark). We assess the extent to which the mass media's reporting of Euroscepticism indicates the Europeanization of public spheres. Using a mixed-methods approach combining LDA topic modeling and qualitative coding, we find that approximately 70 per cent of print articles mentioning "Euroscepticism" or "Eurosceptic" are framed in a non-domestic (i.e. European) context. In five of the six cases studied, articles exhibiting a European context are strikingly similar in content, with the British case as the exception. However, coverage of British Euroscepticism drives Europeanization in other Member States. Bivariate logistic regressions further reveal three macro-level structural variables that significantly correlate with a Europeanized media discourse: newspaper type (tabloid or broadsheet), presence of a strong Eurosceptic party, and relationship to the EU budget (net contributor or receiver of EU funds).


title: Using Deep Reinforcement Learning for the Continuous Control of Robotic Arms

abstract: Deep reinforcement learning enables algorithms to learn complex behavior, deal with continuous action spaces and find good strategies in environments with high dimensional state spaces. With deep reinforcement learning being an active area of research and many concurrent inventions, we decided to focus on a relatively simple robotic task to evaluate a set of ideas that might help to solve recent reinforcement learning problems. We test a newly created combination of two commonly used reinforcement learning methods, whether it is able to learn more effectively than a baseline. We also compare different ideas to preprocess information before it is fed to the reinforcement learning algorithm. The goal of this strategy is to reduce training time and eventually help the algorithm to converge. The concluding evaluation proves the general applicability of the described concepts by testing them using a simulated environment. These concepts might be reused for future experiments.


title: Which domains have two-sided supporting unit spheres at every boundary point?

abstract: We prove the quantitative equivalence of two important geometrical conditions, pertaining to the regularity of a domain $\Omega\subset\mathbb{R}^N$. These are: (i) the uniform two-sided supporting sphere condition, and (ii) the Lipschitz continuity of the outward unit normal vector. In particular, the answer to the question posed in our title is: "Those domains, whose unit normal is well defined and has Lipschitz constant one." We also offer an extension to infinitely dimensional spaces $L^p$, $p\in (1,\infty)$.


title: Assessing the Contribution of Semantic Congruency to Multisensory Integration and Conflict Resolution

abstract: The efficient integration of multisensory observations is a key property of the brain that yields the robust interaction with the environment. However, artificial multisensory perception remains an open issue especially in situations of sensory uncertainty and conflicts. In this work, we extend previous studies on audio-visual (AV) conflict resolution in complex environments. In particular, we focus on quantitatively assessing the contribution of semantic congruency during an AV spatial localization task. In addition to conflicts in the spatial domain (i.e. spatially misaligned stimuli), we consider gender-specific conflicts with male and female avatars. Our results suggest that while semantically related stimuli affect the magnitude of the visual bias (perceptually shifting the location of the sound towards a semantically congruent visual cue), humans still strongly rely on environmental statistics to solve AV conflicts. Together with previously reported results, this work contributes to a better understanding of how multisensory integration and conflict resolution can be modelled in artificial agents and robots operating in real-world environments.


title: Optimally rotated coordinate systems for adaptive least-squares regression on sparse grids

abstract: For low-dimensional data sets with a large amount of data points, standard kernel methods are usually not feasible for regression anymore. Besides simple linear models or involved heuristic deep learning models, grid-based discretizations of larger (kernel) model classes lead to algorithms, which naturally scale linearly in the amount of data points. For moderate-dimensional or high-dimensional regression tasks, these grid-based discretizations suffer from the curse of dimensionality. Here, sparse grid methods have proven to circumvent this problem to a large extent. In this context, space- and dimension-adaptive sparse grids, which can detect and exploit a given low effective dimensionality of nominally high-dimensional data, are particularly successful. They nevertheless rely on an axis-aligned structure of the solution and exhibit issues for data with predominantly skewed and rotated coordinates. In this paper we propose a preprocessing approach for these adaptive sparse grid algorithms that determines an optimized, problem-dependent coordinate system and, thus, reduces the effective dimensionality of a given data set in the ANOVA sense. We provide numerical examples on synthetic data as well as real-world data to show how an adaptive sparse grid least squares algorithm benefits from our preprocessing method.


title: Surface Algebras I: Dessins D'enfants, Surface Algebras, and Dessin Orders

abstract: In this paper, a construction of an infinite dimensional associative algebra, which will be called a \emph{Surface Algebra}, is associated in a "canonical" way to a dessin d'enfant, or more generally, a cellularly embedded graph in a Riemann surface. Once the surface algebras are constructed we will see a construction of what we call here the associated \emph{Dessin Order} or more generally the \emph{Surface Order}. This provides a way of associating to every algebraic curve $X$, with function field $k(X)$ (defined over an arbitrary field $k$) the representation theory of its Surface Algebra and the lattices over Surface Orders, which are defined as pullbacks of certain matrix algebras over commutative $k$-algebras. We will then be able to prove that the center and the (noncommutative) normalization of the surface orders are invariant under the action of the absolute Galois group $\mathcal{G}(\overline{\mathbb{Q}}/\mathbb{Q})$. We will see that the surface algebras and surface orders are closely related to the fundamental group(oid) of the Riemann surfaces and the associated monodromy group. A description of the projective resolutions of the simple modules over the surface order is given and it will be shown that one can completely recover the dessin with the projective resolutions of the simple modules alone. In particular, the projective resolutions of the simple modules encode all combinatorial and topological data of the monodromy group (or cartographic group) of a dessin. Finally, as a corollary we are able to say that classifying dessins in an orbit of $\mathcal{G}(\overline{\mathbb{Q}}/\mathbb{Q})$ is equivalent to classifying dessin orders with a given normalization. We end with some further examples of surface algebras and surface orders related to the classical and geometric version of the Langlands Program.


title: Bosonic SET to SPT Transition under Dyonic LSM Theorem

abstract: Employing large-scale quantum Monte Carlo simulatoins, we study the phase diagram of a quantum spin model which is subject to the recently developed dyonic Lieb-Shultz-Mattis (LSM) theorem. The theorem predicts there are symmetry enriched/protected topological (SET/SPT) phases in the phase diagram. Our numerical results reveal a first order quantum phase transition between SET and SPT phases, consistent with an anyon condensation mechanism that enforces SPT phase according to the theorem. Also there exists in the phase diagram a symmetry-breaking phase in the form of superfluid (SF). The transition between SET and SF is continuous and that between SPT and SF is first order. Interestingly, the SET, SPT and SF phases meet at a critical endpoint, whose presence can be universally explained via theory contains emergent gauge field coupled to vortex fields, and consequently reveals the exotic feature of our model even beyond the realm of dyonic LSM.


title: On the direct images of parabolic vector bundles and parabolic connections

abstract: Let $\varphi : Y \rightarrow X$ be a finite surjective morphism between smooth complex projective curves, where $X$ is irreducible but $Y$ need not be so. Let $V_*$ be a parabolic vector bundle on $Y$. We construct a parabolic structure on the direct image $\varphi_* V$ on $X$, where $V$ is the vector bundle underlying $V_*$. The parabolic vector bundle $\varphi_* V_*$ on $X$ obtained this way has a ramified torus sub-bundle; it is a torus bundle of $\text{Ad}(\varphi_* V)$ outside the parabolic divisor for $\varphi_* V_*$ that satisfies certain conditions at the parabolic points. Conversely, given a parabolic vector bundle $E_*$ on $X$, and a ramified torus sub-bundle $\mathcal T$ for it, we construct a ramified covering $Z$ of $X$ and a parabolic vector bundle $W_*$ on $Z$, such that the parabolic bundle $E_*$ is the direct image of $W_*$. A connection on $V_*$ produces a connection on $\varphi_* V_*$. The ramified torus sub-bundle for $\varphi_* V_*$ is preserved by the logarithmic connection on $\text{End}(\varphi_* V)$ induced by this connection on $\varphi_* V_*$. If the parabolic vector bundle $E_*$ on $X$ is equipped with a connection $D$ such that the connection on the endomorphism bundle induced by it preserves the ramified torus sub-bundle $\mathcal T$, then we prove that the corresponding parabolic vector bundle $W_*$ on $Z$ has a connection that produces the connection $D$ on the direct image $E_*$.


title: The Buchdahl Stability Bound in Eddington-inspired Born-Infeld Gravity

abstract: We give the Buchdahl stability bound in Eddington-inspired Born-Infeld (EiBI) gravity. We show that this bound depends on an energy condition controlled by the model parameter $\kappa$. From this bound, we can constrain $\kappa\lesssim 10^{8}\text{m}^2$ if a neutron star with a mass around $3M_{\odot}$ is observed in the future. In addition, to avoid the potential pathologies in EiBI, a \emph{Hagedorn-like} equation of state associated with $\kappa$ at the center of a compact star is inevitable, which is similar to the Hagedorn temperature in string theory.


title: On the Peaks of a Stochastic Heat Equation on a Sphere with a Large Radius

abstract: For every $R>0$, consider the stochastic heat equation $\partial_{t} u_{R}(t\,,x)=\tfrac12 \Delta_{S_{R}^{2}}u_{R}(t\,,x)+\sigma(u_{R}(t\,,x)) \xi_{R}(t\,,x)$ on $S_{R}^{2}$, where $\xi_{R}=\dot{W_{R}}$ are centered Gaussian noises with the covariance structure given by $E [\dot{W_{R}}(t,x)\dot{W_{R}}(s,y)]=h_{R}(x,y)\delta_{0}(t-s)$, where $h_{R}$ is symmetric and semi-positive definite and there exist some fixed constants $-2< C_{h_{up}}< 2$ and $\frac{1}{2}C_{h_{up}}-1 <C_{h_{lo}}\le C_{h_{up}}$ such that for all $R>0$ and $x\,,y \in S_{R}^{2}$, $(\log R)^{C_{h_{lo}}/2}=h_{lo}(R)\leq h_{R}(x,y) \leq h_{up}(R)=(\log R)^{C_{h_{up}}/2}$, $\Delta_{S_{R}^{2}}$ denotes the Laplace-Beltrami operator defined on $S_{R}^{2}$ and $\sigma:R \mapsto R$ is Lipschitz continuous, positive and uniformly bounded away from $0$ and $\infty$. Under the assumption that $u_{R,0}(x)=u_{R}(0\,,x)$ is a nonrandom continuous function on $x \in S_{R}^{2}$ and the initial condition that there exists a finite positive $U$ such that $\sup_{R>0}\sup_{x \in S_{R}^{2}}\vert u_{R,0}(x)\vert \le U$, we prove that for every finite positive $t$, there exist finite positive constants $C_{low}(t)$ and $C_{up}(t)$ which only depend on $t$ such that as $R \to \infty$, $\sup_{x \in S_{R}^{2}}\vert u_{R}(t\,,x)\vert$ is asymptotically bounded below by $C_{low}(t)(\log R)^{1/4+C_{h_{lo}}/4-C_{h_{up}}/8}$ and asymptotically bounded above by $C_{up}(t)(\log R)^{1/2+C_{h_{up}}/4}$ with high probability.


title: Discovering Fair Representations in the Data Domain

abstract: Interpretability and fairness are critical in computer vision and machine learning applications, in particular when dealing with human outcomes, e.g. inviting or not inviting for a job interview based on application materials that may include photographs. One promising direction to achieve fairness is by learning data representations that remove the semantics of protected characteristics, and are therefore able to mitigate unfair outcomes. All available models however learn latent embeddings which comes at the cost of being uninterpretable. We propose to cast this problem as data-to-data translation, i.e. learning a mapping from an input domain to a fair target domain, where a fairness definition is being enforced. Here the data domain can be images, or any tabular data representation. This task would be straightforward if we had fair target data available, but this is not the case. To overcome this, we learn a highly unconstrained mapping by exploiting statistics of residuals - the difference between input data and its translated version - and the protected characteristics. When applied to the CelebA dataset of face images with gender attribute as the protected characteristic, our model enforces equality of opportunity by adjusting the eyes and lips regions. Intriguingly, on the same dataset we arrive at similar conclusions when using semantic attribute representations of images for translation. On face images of the recent DiF dataset, with the same gender attribute, our method adjusts nose regions. In the Adult income dataset, also with protected gender attribute, our model achieves equality of opportunity by, among others, obfuscating the wife and husband relationship. Analyzing those systematic changes will allow us to scrutinize the interplay of fairness criterion, chosen protected characteristics, and prediction performance.


title: Graphene oxide for enhanced nonlinear optics in low and high index contrast waveguides and nanowires

abstract: We demonstrate enhanced four-wave mixing (FWM) in doped silica waveguides integrated with graphene oxide (GO) layers. Owing to strong mode overlap between the integrated waveguides and GO films that have a high Kerr nonlinearity and low loss, the FWM efficiency of the hybrid integrated waveguides is significantly improved. We perform FWM measurements for different pump powers, wavelength detuning, GO coating lengths, and number of GO layers. Our experimental results show good agreement with theory, achieving up to ~9.5-dB enhancement in the FWM conversion efficiency for a 1.5-cm-long waveguide integrated with 2 layers of GO. We show theoretically that for different waveguide geometries an enhancement in FWM efficiency of ~ 20 dB can be obtained in the doped silica waveguides, and more than 30 dB in silicon nanowires and slot waveguides. This demonstrates the effectiveness of introducing GO films into integrated photonic devices in order to enhance the performance of nonlinear optical processes.


title: An analytically simple and computationally efficient Gaussian beam mode-decomposition approach to classical diffraction theory

abstract: We present a method of Gaussian-beam-mode decomposition to calculate classical diffraction of optical beams by apertures. This method offers a entirely different approach to examine the classic problem. Although our method is based on a very straightforward setup, it is surprisingly effective. We validate our method by comparing its results with those of Kirchhoff's full diffraction formula. Not only does our method have a simple and organized analytical framework, it also offers significant computational advantage.


title: Discriminator Rejection Sampling

abstract: We propose a rejection sampling scheme using the discriminator of a GAN to approximately correct errors in the GAN generator distribution. We show that under quite strict assumptions, this will allow us to recover the data distribution exactly. We then examine where those strict assumptions break down and design a practical algorithm - called Discriminator Rejection Sampling (DRS) - that can be used on real data-sets. Finally, we demonstrate the efficacy of DRS on a mixture of Gaussians and on the SAGAN model, state-of-the-art in the image generation task at the time of developing this work. On ImageNet, we train an improved baseline that increases the Inception Score from 52.52 to 62.36 and reduces the Frechet Inception Distance from 18.65 to 14.79. We then use DRS to further improve on this baseline, improving the Inception Score to 76.08 and the FID to 13.75.


title: A Block Coordinate Descent Proximal Method for Simultaneous Filtering and Parameter Estimation

abstract: We propose and analyze a block coordinate descent proximal algorithm (BCD-prox) for simultaneous filtering and parameter estimation of ODE models. As we show on ODE systems with up to d=40 dimensions, as compared to state-of-the-art methods, BCD-prox exhibits increased robustness (to noise, parameter initialization, and hyperparameters), decreased training times, and improved accuracy of both filtered states and estimated parameters. We show how BCD-prox can be used with multistep numerical discretizations, and we establish convergence of BCD-prox under hypotheses that include real systems of interest.


title: Frozen Gaussian Approximation for 3-D Elastic Wave Equation and Seismic Tomography

abstract: The purpose of this work is to generalize the frozen Gaussian approximation (FGA) theory to solve the 3-D elastic wave equation and use it as the forward modeling tool for seismic tomography with high-frequency data. FGA has been previously developed and verified as an efficient solver for high-frequency acoustic wave propagation (P-wave). The main contribution of this paper consists of three aspects: 1. We derive the FGA formulation for the 3-D elastic wave equation. Rather than standard ray-based methods (e.g. geometric optics and Gaussian beam method), the derivation requires to do asymptotic expansion in the week sense (integral form) so that one is able to perform integration by parts. Compared to the FGA theory for acoustic wave equation, the calculations in the derivation are much more technically involved due to the existence of both P- and S-waves, and the coupling of the polarized directions for SH- and SV-waves. In particular, we obtain the diabatic coupling terms for SH- and SV-waves, with the form closely connecting to the concept of Berry phase that is intensively studied in quantum mechanics and topology (Chern number). The accuracy and parallelizability of the FGA algorithm is illustrated by comparing to the spectral element method for 3-D elastic wave equation in homogeneous media; 2. We derive the interface conditions of FGA for 3-D elastic wave equation based on an Eulerian formulation and the Snell's law. We verify these conditions by simulating high-frequency elastic wave propagation in a 1-D layered Earth model. In this example, we also show that it is natural to apply the FGA algorithm to geometries with non-Cartesian coordinates; 3. We apply the developed FGA algorithm for 3-D seismic { wave-equation-based traveltime tomography and full waveform inversion, respectively


title: Extreme enhancement of spin relaxation mediated by surface magnon polaritons

abstract: Polaritons in metals, semimetals, semiconductors, and polar insulators, with their extreme confinement of electromagnetic energy, provide many promising opportunities for enhancing typically weak light-matter interactions such as multipolar radiation, multiphoton spontaneous emission, Raman scattering, and material nonlinearities. These highly confined polaritons are quasi-electrostatic in nature, with most of their energy residing in the electric field. As a result, these "electric" polaritons are far from optimized for enhancing emission of a magnetic nature, such as spin relaxation, which is typically many orders of magnitude slower than corresponding electric decays. Here, we propose using surface magnon polaritons in negative magnetic permeability materials such as MnF$_2$ and FeF$_2$ to strongly enhance spin-relaxation in nearby emitters in the THz spectral range. We find that these magnetic polaritons in 100 nm thin-films can be confined to lengths over 10,000 times smaller than the wavelength of a photon at the same frequency, allowing for a surprising twelve orders of magnitude enhancement in magnetic dipole transitions. This takes THz spin-flip transitions, which normally occur at timescales on the order of a year, and forces them to occur at sub-ms timescales. Our results suggest an interesting platform for polaritonics at THz frequencies, and more broadly, a new way to use polaritons to control light-matter interactions.


title: Convex expansion for finite distributive lattices with applications

abstract: The concept of cutting is first explicitly introduced. By the concept, a convex expansion for finite distributive lattices is considered. Thus, a more general method for drawing the Hasse diagram is given, and the rank generating function of a finite distributive lattice is obtained. In addition, we have several enumerative properties on finite distributive lattices and verify the generalized Euler formula for polyhedrons.


title: Gelfand-Tsetlin degeneration of shift of argument subalgebras in type D

abstract: The universal enveloping algebra of any semisimple Lie algebra $\mathfrak{g}$ contains a family of maximal commutative subalgebras, called shift of argument subalgebras, parametrized by regular Cartan elements of $\mathfrak{g}$. For $\mathfrak{g}=\mathfrak{gl}_n$ the Gelfand-Tsetlin commutative subalgebra in $U(\mathfrak{g})$ arises as some limit of subalgebras from this family. In our previous work (arXiv:1807.11126) we studied the analogous limit of shift of argument subalgebras for the Lie algebras $\mathfrak{g}=\mathfrak{sp}_{2n}$ and $\mathfrak{g}=\mathfrak{o}_{2n+1}$. We described the limit subalgebras in terms of Bethe subalgebras of twisted Yangians $Y^-(2)$ and $Y^+(2)$, respectively, and parametrized the eigenbases of these limit subalgebras in the finite dimensional irreducible highest weight representations by Gelfand-Tsetlin patterns of types C and B. In this note we state and prove similar results for the last case of classical Lie algebras, $\mathfrak{g}=\mathfrak{o}_{2n}$. We describe the limit shift of argument subalgebra in terms of the Bethe subalgebra in the twisted Yangian $Y^+(2)$ and give a natural indexing of its eigenbasis in any finite dimensional irreducible highest weight $\mathfrak{g}$-module by type D Gelfand-Tsetlin patterns.


title: Simple Policy Evaluation for Data-Rich Iterative Tasks

abstract: A data-based policy for iterative control task is presented. The proposed strategy is model-free and can be applied whenever safe input and state trajectories of a system performing an iterative task are available. These trajectories, together with a user-defined cost function, are exploited to construct a piecewise affine approximation to the value function. Approximated value functions are then used to evaluate the control policy by solving a linear program. We show that for linear system subject to convex cost and constraints, the proposed strategy guarantees closed-loop constraint satisfaction and performance bounds on the closed-loop trajectory. We evaluate the proposed strategy in simulations and experiments, the latter carried out on the Berkeley Autonomous Race Car (BARC) platform. We show that the proposed strategy is able to reduce the computation time by one order of magnitude while achieving the same performance as our model-based control algorithm.


title: A survey of automatic de-identification of longitudinal clinical narratives

abstract: Use of medical data, also known as electronic health records, in research helps develop and advance medical science. However, protecting patient confidentiality and identity while using medical data for analysis is crucial. Medical data can be in the form of tabular structures (i.e. tables), free-form narratives, and images. This study focuses on medical data in the free form longitudinal text. De-identification of electronic health records provides the opportunity to use such data for research without it affecting patient privacy, and avoids the need for individual patient consent. In recent years there is increasing interest in developing an accurate, robust and adaptable automatic de-identification system for electronic health records. This is mainly due to the dilemma between the availability of an abundance of health data, and the inability to use such data in research due to legal and ethical restrictions. De-identification tracks in competitions such as the 2014 i2b2 UTHealth and the 2016 CEGS N-GRID shared tasks have provided a great platform to advance this area. The primary reasons for this include the open source nature of the dataset and the fact that raw psychiatric data were used for 2016 competitions. This study focuses on noticeable trend changes in the techniques used in the development of automatic de-identification for longitudinal clinical narratives. More specifically, the shift from using conditional random fields (CRF) based systems only or rules (regular expressions, dictionary or combinations) based systems only, to hybrid models (combining CRF and rules), and more recently to deep learning based systems. We review the literature and results that arose from the 2014 and the 2016 competitions and discuss the outcomes of these systems. We also provide a list of research questions that emerged from this survey.


title: DN-ResNet: Efficient Deep Residual Network for Image Denoising

abstract: A deep learning approach to blind denoising of images without complete knowledge of the noise statistics is considered. We propose DN-ResNet, which is a deep convolutional neural network (CNN) consisting of several residual blocks (ResBlocks). With cascade training, DN-ResNet is more accurate and more computationally efficient than the state of art denoising networks. An edge-aware loss function is further utilized in training DN-ResNet, so that the denoising results have better perceptive quality compared to conventional loss function. Next, we introduce the depthwise separable DN-ResNet (DS-DN-ResNet) utilizing the proposed Depthwise Seperable ResBlock (DS-ResBlock) instead of standard ResBlock, which has much less computational cost. DS-DN-ResNet is incrementally evolved by replacing the ResBlocks in DN-ResNet by DS-ResBlocks stage by stage. As a result, high accuracy and good computational efficiency are achieved concurrently. Whereas previous state of art deep learning methods focused on denoising either Gaussian or Poisson corrupted images, we consider denoising images having the more practical Poisson with additive Gaussian noise as well. The results show that DN-ResNets are more efficient, robust, and perform better denoising than current state of art deep learning methods, as well as the popular variants of the BM3D algorithm, in cases of blind and non-blind denoising of images corrupted with Poisson, Gaussian or Poisson-Gaussian noise. Our network also works well for other image enhancement task such as compressed image restoration.


title: Approximate Fisher Information Matrix to Characterise the Training of Deep Neural Networks

abstract: In this paper, we introduce a novel methodology for characterising the performance of deep learning networks (ResNets and DenseNet) with respect to training convergence and generalisation as a function of mini-batch size and learning rate for image classification. This methodology is based on novel measurements derived from the eigenvalues of the approximate Fisher information matrix, which can be efficiently computed even for high capacity deep models. Our proposed measurements can help practitioners to monitor and control the training process (by actively tuning the mini-batch size and learning rate) to allow for good training convergence and generalisation. Furthermore, the proposed measurements also allow us to show that it is possible to optimise the training process with a new dynamic sampling training approach that continuously and automatically change the mini-batch size and learning rate during the training process. Finally, we show that the proposed dynamic sampling training approach has a faster training time and a competitive classification accuracy compared to the current state of the art.


title: SINE: Scalable Incomplete Network Embedding

abstract: Attributed network embedding aims to learn low-dimensional vector representations for nodes in a network, where each node contains rich attributes/features describing node content. Because network topology structure and node attributes often exhibit high correlation, incorporating node attribute proximity into network embedding is beneficial for learning good vector representations. In reality, large-scale networks often have incomplete/missing node content or linkages, yet existing attributed network embedding algorithms all operate under the assumption that networks are complete. Thus, their performance is vulnerable to missing data and suffers from poor scalability. In this paper, we propose a Scalable Incomplete Network Embedding (SINE) algorithm for learning node representations from incomplete graphs. SINE formulates a probabilistic learning framework that separately models pairs of node-context and node-attribute relationships. Different from existing attributed network embedding algorithms, SINE provides greater flexibility to make the best of useful information and mitigate negative effects of missing information on representation learning. A stochastic gradient descent based online algorithm is derived to learn node representations, allowing SINE to scale up to large-scale networks with high learning efficiency. We evaluate the effectiveness and efficiency of SINE through extensive experiments on real-world networks. Experimental results confirm that SINE outperforms state-of-the-art baselines in various tasks, including node classification, node clustering, and link prediction, under settings with missing links and node attributes. SINE is also shown to be scalable and efficient on large-scale networks with millions of nodes/edges and high-dimensional node features. The source code of this paper is available at https://github.com/daokunzhang/SINE.


title: An alternative stable solution for the Kepler-419 system, obtained with the use of a genetic algorithm

abstract: The mid-transit times of an exoplanet may be non-periodic. The variations in the timing of the transits with respect to a single period, that is, the transit timing variations (TTVs), can sometimes be attributed to perturbations by other exoplanets present in the system, which may or may not transit the star. Our aim is to compute the mass and the six orbital elements of an non-transiting exoplanet, given only the central times of transit of the transiting body. We also aim to recover the mass of the star and the mass and orbital elements of the transiting exoplanet, suitably modified in order to decrease the deviation between the observed and the computed transit times by as much as possible. We have applied our method, based on a genetic algorithm, to the Kepler-419 system. We were able to compute all fourteen free parameters of the system, which, when integrated in time, give transits within the observational errors. We also studied the dynamics and the long-term orbital evolution of the Kepler-419 planetary system as defined by the orbital elements computed by us, in order to determine its stability.


title: Science with an ngVLA - The Molecular High-z Universe on Large Scales: Low-surface-brightness CO and the strength of the ngVLA Core

abstract: The Next-Generation Very Large Array (ngVLA) will revolutionize our understanding of the Early Universe by tracing the coldest phase of molecular gas -the raw ingredient for star formation- in the most distant galaxies and galaxy-clusters. The km-scale core of the ngVLA will be densely packed with antennas, making it a prime instrument for imaging low-surface-brightness emission from large-scale molecular gas in the high-z circum- and inter-galactic medium (CGM/IGM). Recent studies indicate that large amounts of cold molecular gas are hiding in the 10s-100 kpc environments of distant galaxies, but that technical limitations on existing telescope arrays have prevented us from efficiently detecting these large molecular reservoirs. This may have led to a severely biased view of the molecular Universe. We present the science case for low-surface-brightness CO observations of the Early Universe, and how the core of the ngVLA will reveal the cold molecular Universe to limits and at scales not currently detectable by radio telescopes. As such, the ngVLA core will be a powerful instrument for studying the cold baryon cycle that drives the early evolution of galaxies and clusters.


title: Antisymmetry of the stochastic order on all ordered topological spaces

abstract: In this short note, we prove that the stochastic order of Radon probability measures on any ordered topological space is antisymmetric. This has been known before in various special cases. We give a simple and elementary proof of the general result.


title: Studies of two-dimensional MoS2 on enhancing the electrical performance of ultrathin copper films

abstract: Copper nanowires are widely used as on-chip interconnects due to superior conductivity. However, with aggressive Cu interconnect scaling, the diffusive surface scattering of electrons drastically increases the electrical resistivity. In this work, we studied the electrical performance of Cu thin films on different materials. By comparing the thickness dependence of Cu films resistivity on MoS2 and SiO2, we demonstrated that two-dimensional MoS2 can be used to enhance the electrical performance of ultrathin Cu films due to a partial specular surface scattering. By fitting the experimental data with the theoretical Fuchs Sondheimer model, we obtained the specularity parameter at the Cu MoS2 interface in the temperature range 2K to 300K. Furthermore, first principle calculations based on the density functional theory indicates that there are more localized states at the Cu amorphous SiO2 interface than the Cu MoS2 interface which is responsible for the higher resistivity in the Cu SiO2 heterostructure due to more severe electron scattering. Our results suggest that Cu MoS2 hybrid is a promising candidate structure for the future generations of CMOS interconnects.


title: Evolutionary Stochastic Gradient Descent for Optimization of Deep Neural Networks

abstract: We propose a population-based Evolutionary Stochastic Gradient Descent (ESGD) framework for optimizing deep neural networks. ESGD combines SGD and gradient-free evolutionary algorithms as complementary algorithms in one framework in which the optimization alternates between the SGD step and evolution step to improve the average fitness of the population. With a back-off strategy in the SGD step and an elitist strategy in the evolution step, it guarantees that the best fitness in the population will never degrade. In addition, individuals in the population optimized with various SGD-based optimizers using distinct hyper-parameters in the SGD step are considered as competing species in a coevolution setting such that the complementarity of the optimizers is also taken into account. The effectiveness of ESGD is demonstrated across multiple applications including speech recognition, image recognition and language modeling, using networks with a variety of deep architectures.


title: Intersections of subcomplexes in non-positively curved 2-dimensional complexes

abstract: Let $X$ be a contractible $2$-complex which is a union of two contractible subcomplexes $Y$ and $Z.$ Is the intersection $Y\cap Z$ contractible as well? In this note, we prove that the inclusion-induced map $\pi _{1}(Y\cap Z)\rightarrow \pi _{1}(Z)$ is injective if $Y$ is $\pi _{1}$-injective subcomplex in a locally CAT(0) 2-complex $X$. In particular, each component in the intersection of two contractible subcomplexes in a CAT(0) 2-complex is contractible.


title: Riemannian Inexact Newton Method for Structured Inverse Eigenvalue and Singular Value Problems

abstract: Inverse eigenvalue and singular value problems have been widely discussed for decades. The well-known result is the Weyl-Horn condition, which presents the relations between the eigenvalues and singular values of an arbitrary matrix. This result by Weyl-Horn then leads to an interesting inverse problem, i.e., how to construct a matrix with desired eigenvalues and singular values. In this work, we do that and more. We propose an eclectic mix of techniques from differential geometry and the inexact Newton method for solving inverse eigenvalue and singular value problems as well as additional desired characteristics such as nonnegative entries, prescribed diagonal entries, and even predetermined entries. We show theoretically that our method converges globally and quadratically, and we provide numerical examples to demonstrate the robustness and accuracy of our proposed method. {Having theoretical interest, we provide in the appendix a necessary and sufficient condition for the existence of a $2\times 2$ real matrix, or even a nonnegative matrix, with prescribed eigenvalues, singular values, and main diagonal entries.


title: Super-resolution MRI through Deep Learning

abstract: Magnetic resonance imaging (MRI) is extensively used for diagnosis and image-guided therapeutics. Due to hardware, physical and physiological limitations, acquisition of high-resolution MRI data takes long scan time at high system cost, and could be limited to low spatial coverage and also subject to motion artifacts. Super-resolution MRI can be achieved with deep learning, which is a promising approach and has a great potential for preclinical and clinical imaging. Compared with polynomial interpolation or sparse-coding algorithms, deep learning extracts prior knowledge from big data and produces superior MRI images from a low-resolution counterpart. In this paper, we adapt two state-of-the-art neural network models for CT denoising and deblurring, transfer them for super-resolution MRI, and demonstrate encouraging super-resolution MRI results toward two-fold resolution enhancement.


title: Effects of migration on vector-borne diseases with forward and backward stage progression

abstract: Is it possible to break the host-vector chain of transmission when there is an influx of infectious hosts into a na\"{i}ve population and competent vector? To address this question, a class of vector-borne disease models with an arbitrary number of infectious stages that account for immigration of infective individuals is formulated. The proposed model accounts for forward and backward progression, capturing the mitigation and aggravation to and from any stages of the infection, respectively. The model has a rich dynamic, which depends on the patterns of infected immigrant influx into the host population and connectivity of the transfer between infectious classes. We provide conditions under which the answer of the initial question is positive.


title: Double Ore extensions versus graded skew PBW extensions

abstract: In this paper we present necessary and sufficient conditions for a graded (trimmed) double Ore extension to be a graded (quasi-commutative) skew PBW extension. Using this fact, we prove that a graded skew PBW extension $A = \sigma(R)\langle x_1,x_2 \rangle$ of an Artin-Schelter regular algebra $R$ is Artin-Schelter regular. As a consequence, every graded skew PBW extension $A = \sigma(R)\langle x_1,x_2 \rangle$ of a connected skew Calabi-Yau algebra $R$ of dimension $d$ is skew Calabi-Yau of dimension $d+2$.


title: Science with an ngVLA: A six-dimensional tomographic view of Galactic star-formation

abstract: Various sign-posts of recent star-formation activity, such as water and methanol maser emission or magnetically active low-mass young stars, can be detected with Very Long Baseline Interferometry (VLBI) radio arrays. The extremely accurate astrometry already attainable with VLBI instruments implies that the trigonometric parallax and the proper motion of these objects can be measured to better than 1% for sources within a few hundred parsecs of the Sun, and better than 10% for objects at a few kiloparsecs. An ngVLA with baselines extending to several thousand km will have a sensitivity more than one order of magnitude better than current VLBI instruments, and will enable such highly accurate astrometric measurements to be performed throughout the Milky Way. This will provide a full six-dimensional view (three spatial and three velocity coordinates) of star-formation in the Galactic disk, and have a transformative impact on our understanding of both star-formation processes and Galactic structure.


title: Subalgebras, idempotents, ideals and quasi-units of two-dimensional algebras

abstract: All subalgebras, idempotents, left(right) ideals and left quasi-units of two-dimensional algebras are described. Classification of algebras with given number of subalgebras, left(right) ideals are provided. In particular, a list of isomorphism classes of all simple two-dimensional algebras is given. In the study of ideals and subalgebras the number of them depend on roots of certain system of polynomials at structure constants of the algebra. We also give explicit forms of the polynomials.


title: On the local pairing behavior of critical points and roots of random polynomials

abstract: We study the pairing between zeros and critical points of the polynomial $p_n(z) = \prod_{j=1}^n(z-X_j)$, whose roots $X_1, \ldots, X_n$ are complex-valued random variables. Under a regularity assumption, we show that if the roots are independent and identically distributed, the Wasserstein distance between the empirical distributions of roots and critical points of $p_n$ is on the order of $1/n$, up to logarithmic corrections. The proof relies on a careful construction of disjoint random Jordan curves in the complex plane, which allow us to naturally pair roots and nearby critical points. In addition, we establish asymptotic expansions to order $1/n^2$ for the locations of the nearest critical points to several fixed roots. This allows us to describe the joint limiting fluctuations of the critical points as $n$ tends to infinity, extending a recent result of Kabluchko and Seidel. Finally, we present a local law that describes the behavior of the critical points when the roots are neither independent nor identically distributed.


title: Squeezed back-to-back correlation between bosons and antibosons with different in-medium masses in high-energy heavy-ion collisions

abstract: We derive the formulas of the squeezed back-to-back correlation (SBBC) between a boson and antiboson with different in-medium masses in high-energy heavy-ion collisions. The influence of the in-medium mass difference between a boson and antiboson on the SBBC is investigated. We calculate the SBBC functions of $D$-meson pairs for the hydrodynamic sources described by the VISH2+1 code for Au+Au collisions at $\sqrt{s_{NN}} =200$~GeV. Our results indicate that the SBBC strengths of $D^+D^-$ and $D^0{\bar D}^0$ are different if there are charge-dependent in-medium interactions.


title: Tracking Pre-inflation Era from density perturbation spectra

abstract: One of the great triumphs of the inflationary model is the prediction of the flat power spectrum of the CMB fluctuation. The prediction is based on the assumption of the de-Sitter vacuum in the past infinity. However, the true past infinity of the inflation is expected to be dominated by radiation and curvature of the space. We consider pre-inflation era as dominated by radiation and curvatures as well as inflation potential. We derive the exact solutions for the scalar fields in this era and find a exact power spectra caused by the inflaton vacuum fluctuation. We show that the power spectrum is almost flat for sub-horizon scale and deviates from flat for very high super-horizon fluctuation, which is quite sensitive to the radiation and the curvature in the pre-inflation era.


title: ProMP: Proximal Meta-Policy Search

abstract: Credit assignment in Meta-reinforcement learning (Meta-RL) is still poorly understood. Existing methods either neglect credit assignment to pre-adaptation behavior or implement it naively. This leads to poor sample-efficiency during meta-training as well as ineffective task identification strategies. This paper provides a theoretical analysis of credit assignment in gradient-based Meta-RL. Building on the gained insights we develop a novel meta-learning algorithm that overcomes both the issue of poor credit assignment and previous difficulties in estimating meta-policy gradients. By controlling the statistical distance of both pre-adaptation and adapted policies during meta-policy search, the proposed algorithm endows efficient and stable meta-learning. Our approach leads to superior pre-adaptation policy behavior and consistently outperforms previous Meta-RL algorithms in sample-efficiency, wall-clock time, and asymptotic performance.


title: Nonassociative right hoops

abstract: The class of nonassociative right hoops, or narhoops for short, is defined as a subclass of right-residuated magmas, and is shown to be a variety. These algebras generalize both right quasigroups and right hoops, and we characterize the subvarieties in which the operation $x\sqcap y=(x / y)y$ is associative and/or commutative. Narhoops with a left unit are proved to have a top element if and only if $\sqcap$ is commutative, and their congruences are determined by the equivalence class of the left unit. We also show that the four identities defining narhoops are independent.


title: Graphene and Carbon Nanotubes in Space

abstract: As the fourth most abundant element in the universe, carbon plays an important role in the physical and chemical evolution of the interstellar medium. Due to its unique property to form three different types of chemical bonds through sp^1, sp^2, and sp^3 hybridizations, carbon can be stabilized in various allotropes, including amorphous carbon, graphite, diamond, polycyclic aromatic hydrocarbon, fullerenes, graphene, and carbon nanotubes.


Supplementary Material for "Estimation of a Multiplicative Correlation Structure in the Large Dimensional Case"


title: Narrow Transient Absorptions in Late-Time Optical Spectra of Type Ia Supernovae: Evidence for Large Clumps of Iron-Rich Ejecta?

abstract: An examination of late-time, optical spectra of type Ia supernovae revealed surprisingly narrow absorption features which only become visible a few months after maximum light. These features, most clearly seen in the late-time spectra of the bright, recent type Ia supernovae ASASSN-14lp and SN 2017bzc, appear as narrow absorptions at ~4840 A, ~5000 A, and as a sharp inflection at ~4760 A on the red side of the prominent late-time 4700 A feature. A survey of on-line archival data revealed similar features present in the spectra of ten other normal and 91T-like SNe Ia, including SN 2011fe. Unlike blue spectral features which exhibit progressive red-ward shifts, these narrow absorptions remain at the same wavelength from epoch to epoch for an individual SN, but can appear at slightly different wavelengths for each object. These features are also transient, appearing and then fading in one to three months. After ruling out instrumental, data reduction, and atmospheric affects, we discuss possible explanations including progenitor mass-loss material, interaction with material from previous novae events, and absorption by large discrete clumps of high-velocity Fe-rich ejecta.


title: The Big Bang is a Coordinate Singularity for $k = -1$ Inflationary FLRW Spacetimes

abstract: We show that the big bang is a coordinate singularity for a large class of $k = -1$ inflationary FLRW spacetimes which we have dubbed 'Milne-like.' By introducing a new set of coordinates, the big bang appears as a past boundary of the universe where the metric is no longer degenerate -- a result which has already been investigated in the context of vacuum decay [12]. We generalize their results and approach the problem from a more mathematical perspective. Similar to how investigating the geometrical properties of the $r = 2m$ event horizon in Schwarzschild led to a better understanding of black holes, we believe that investigating the geometrical properties of the big bang coordinate singularity for Milne-like spacetimes could lead to a better understanding of cosmology. We show how the mathematics of these spacetimes may help illuminate certain issues associated with dark energy, dark matter, and the universe's missing antimatter.


title: Science with an ngVLA: ngVLA Studies of Classical Novae

abstract: Observations with modern radio telescopes have revealed that classical novae are far from the simple, spherically symmetric events they were once assumed to be. It is now understood that novae provide excellent laboratories to study several astrophysical properties including binary interactions, stellar outflows, and shock physics. The ngVLA will provide unprecedented opportunities to study these events. It will enable us to observe more distant and fainter novae than we can today. It will allow us to simultaneously resolve both the thermal and non-thermal components in the ejecta. Finally, monitoring novae with the ngVLA will reveal the evolution of the ejecta in better detail than is possible with any current instrument.


title: Linear relations among roots of a polynomial

abstract: Let $f(x)$ be a monic polynomial over $\mathbb{Q}$ with complex roots $\alpha_1,\dots,\alpha_n$. Linear relations among them and $1$ over $\mathbb{Q}$ play an important role when we study the distribution of roots modulo a prime. We study linear relations group-theoretically. If $\mathbb{Q}(\alpha_1)$ is an abelian extension, we can understand the structure completely.


title: Appendix to a paper [arXiv:1809.08623] of B. Williams

abstract: In this appendix to a paper [arXiv:1809.08623] by B. Williams, we give birational equivalences between the models of the Hilbert modular surfaces for ${\mathbb Q}(\sqrt{29})$ and ${\mathbb Q}(\sqrt{37})$ given there and those previously found by Elkies and Kumar.


title: Learning Two-layer Neural Networks with Symmetric Inputs

abstract: We give a new algorithm for learning a two-layer neural network under a general class of input distributions. Assuming there is a ground-truth two-layer network $$ y = A \sigma(Wx) + \xi, $$ where $A,W$ are weight matrices, $\xi$ represents noise, and the number of neurons in the hidden layer is no larger than the input or output, our algorithm is guaranteed to recover the parameters $A,W$ of the ground-truth network. The only requirement on the input $x$ is that it is symmetric, which still allows highly complicated and structured input. Our algorithm is based on the method-of-moments framework and extends several results in tensor decompositions. We use spectral algorithms to avoid the complicated non-convex optimization in learning neural networks. Experiments show that our algorithm can robustly learn the ground-truth neural network with a small number of samples for many symmetric input distributions.


title: Elimination of basal-plane stacking faults in semipolar/nonpolar GaN and light emitting diodes heteroepitaxially grown on sapphire substrates

abstract: High quality semipolar and nonpolar GaN is crucial in achieving high-performance GaN-based optoelectronic devices, yet it has been very challenging to achieve large-area wafers that are free of basal-plane stacking faults (BSFs). In this work, we report an approach to prepare large-area, stacking-fault-free (SF-free) semipolar GaN on (4-inch) sapphire substrates. A root cause of the formation of BSFs is the emergence of N-polar (000-1) facets during semipolar and non-polar heteroepitaxy. Invoking the concept of kinetic Wulff plot, we succeeded in suppressing the occurrence of N-polar GaN (000-1) facets, and consequently in eliminating the stacking faults generated in (000-1) basal-planes. The result was confirmed by transmission electron microscopy, cathodoluminescence, and low-temperature photoluminescence characterizations. Furthermore, InGaN light emitting diodes with promising characteristics have been produced on the SF-free semipolar (20-21) GaN on sapphire substrates. Our work opens up a new insight about the heteroepitaxial growth of nonpolar/semipolar GaN and provides an approach of producing SF-free nonpolar/semipolar GaN material over large-area wafers which will create new opportunities in GaN optoelectronic and microelectronic research.


title: Comment on "keV Neutrino Dark Matter in a Fast Expanding Universe" by Biswas et al

abstract: Biswas et al. found that the thermal relic density of a dark matter particle freezing out while the universe's energy density is dominated by a non-standard extra component $\phi$, whose energy density redshifts faster than radiation, can be greatly suppressed. Here we show that this result, which contradicts extensive previous literature, is incorrect: the mistake lies with the assumption that the (decoupled) extra component $\phi$ contributes to the entropic degrees of freedom relevant for dark matter freeze out. If this were the case, a completely different approach would be needed to calculate the dark matter relic abundance, with dramatically different results.


title: StructVIO : Visual-inertial Odometry with Structural Regularity of Man-made Environments

abstract: We propose a novel visual-inertial odometry approach that adopts structural regularity in man-made environments. Instead of using Manhattan world assumption, we use Atlanta world model to describe such regularity. An Atlanta world is a world that contains multiple local Manhattan worlds with different heading directions. Each local Manhattan world is detected on-the-fly, and their headings are gradually refined by the state estimator when new observations are coming. With fully exploration of structural lines that aligned with each local Manhattan worlds, our visual-inertial odometry method become more accurate and robust, as well as much more flexible to different kinds of complex man-made environments. Through extensive benchmark tests and real-world tests, the results show that the proposed approach outperforms existing visual-inertial systems in large-scale man-made environments


title: A Robust Local Binary Similarity Pattern for Foreground Object Detection

abstract: Accurate and fast extraction of the foreground object is one of the most significant issues to be solved due to its important meaning for object tracking and recognition in video surveillance. Although many foreground object detection methods have been proposed in the recent past, it is still regarded as a tough problem due to illumination variations and dynamic backgrounds challenges. In this paper, we propose a robust foreground object detection method with two aspects of contributions. First, we propose a robust texture operator named Robust Local Binary Similarity Pattern (RLBSP), which shows strong robustness to illumination variations and dynamic backgrounds. Second, a combination of color and texture features are used to characterize pixel representations, which compensate each other to make full use of their own advantages. Comprehensive experiments evaluated on the CDnet 2012 dataset demonstrate that the proposed method performs favorably against state-of-the-art methods.


title: The Jellium Edge and the Size Effect of the Chemical Potential and Surface Energy in Metal Slabs

abstract: Although free electron models have been established in order to capture the essential physics of interfacial and bulk properties in metals, some issues still remain regarding the application of free electron models to thin metal films. One of the issues relates to whether the geometric edge coincides with the potential edge in order to satisfy the charge neutrality condition when the potential profile is modeled as a rectangular potential well. We show that they coincide by rigorously taking into account the quantization effect arising from electron confinement in a thin metal slab. As a result, the overall behaviors of the chemical potential and surface energy show an increasing trend by decreasing the thickness of the slab. The chemical potential and surface energy show an oscillatory thickness dependence by further taking into account the discreteness of the total number of free electrons.


title: A systematic study of Galactic infrared bubbles along the Galactic plane with AKARI and Herschel

abstract: Galactic infrared (IR) bubbles, which have shell-like structures in the mid-IR wavelengths, are known to contain massive stars near their centers. IR bubbles in inner Galactic regions ($|$l$|\leq$ 65$^{\circ}$, $|$b$|\leq$ 1$^{\circ}$) have so far been studied well to understand the massive star formation mechanisms. In this study, we expand the research area to the whole Galactic plane (0$^{\circ}\leq$ l $<$360$^{\circ}$, $|$b$|\leq$ 5$^{\circ}$), using the AKARI all-sky survey data. We limit our study on large bubbles with angular radii of $>1'$ to reliably identify and characterize them. For the 247 IR bubbles in total, we derived the radii and the covering fractions of the shells, based on the method developed in \citet{Hattori2016}. We also created their spectral energy distributions, using the AKARI and Herschel photometric data, and decomposed them with a dust model, to obtain the total IR luminosity and the luminosity of each dust component, i.e., polycyclic aromatic hydrocarbons (PAHs), warm dust and cold dust. As a result, we find that there are systematic differences in the IR properties of the bubbles between inner and outer Galactic regions. The total IR luminosities are lower in outer Galactic regions, while there is no systematic difference in the range of the shell radii between inner and outer Galactic regions. More IR bubbles tend to be observed as broken bubbles rather than closed ones and the fractional luminosities of the PAH emission are significantly higher in outer Galactic regions. We discuss the implications of these results for the massive stars and the interstellar environments associated with the Galactic IR bubbles.


title: Percolation and jamming of random sequential adsorption samples of large linear $k$-mers on a square lattice

abstract: The behavior of the percolation threshold and the jamming coverage for isotropic random sequential adsorption samples has been studied by means of numerical simulations. A parallel algorithm that is very efficient in terms of its speed and memory usage has been developed and applied to the model involving large linear $k$-mers on a square lattice with periodic boundary conditions. We have obtained the percolation thresholds and jamming concentrations for lengths of $k$-mers up to $2^{17}$. New large $k$ regime of the percolation threshold behavior has been identified. The structure of the percolating and jamming states has been investigated. The theorem of G.~Kondrat, Z.~Koza, and P.~Brzeski [Phys. Rev. E 96, 022154 (2017)] has been generalized to the case of periodic boundary conditions. We have proved that any cluster at jamming is percolating cluster and that percolation occurs before jamming.


title: Quasi-hyperbolic momentum and Adam for deep learning

abstract: Momentum-based acceleration of stochastic gradient descent (SGD) is widely used in deep learning. We propose the quasi-hyperbolic momentum algorithm (QHM) as an extremely simple alteration of momentum SGD, averaging a plain SGD step with a momentum step. We describe numerous connections to and identities with other algorithms, and we characterize the set of two-state optimization algorithms that QHM can recover. Finally, we propose a QH variant of Adam called QHAdam, and we empirically demonstrate that our algorithms lead to significantly improved training in a variety of settings, including a new state-of-the-art result on WMT16 EN-DE. We hope that these empirical results, combined with the conceptual and practical simplicity of QHM and QHAdam, will spur interest from both practitioners and researchers. Code is immediately available.


title: Efficient Continuous-Time SLAM for 3D Lidar-Based Online Mapping

abstract: Modern 3D laser-range scanners have a high data rate, making online simultaneous localization and mapping (SLAM) computationally challenging. Recursive state estimation techniques are efficient but commit to a state estimate immediately after a new scan is made, which may lead to misalignments of measurements. We present a 3D SLAM approach that allows for refining alignments during online mapping. Our method is based on efficient local mapping and a hierarchical optimization back-end. Measurements of a 3D laser scanner are aggregated in local multiresolution maps by means of surfel-based registration. The local maps are used in a multi-level graph for allocentric mapping and localization. In order to incorporate corrections when refining the alignment, the individual 3D scans in the local map are modeled as a sub-graph and graph optimization is performed to account for drift and misalignments in the local maps. Furthermore, in each sub-graph, a continuous-time representation of the sensor trajectory allows to correct measurements between scan poses. We evaluate our approach in multiple experiments by showing qualitative results. Furthermore, we quantify the map quality by an entropy-based measure.


title: Co-manifold learning with missing data

abstract: Representation learning is typically applied to only one mode of a data matrix, either its rows or columns. Yet in many applications, there is an underlying geometry to both the rows and the columns. We propose utilizing this coupled structure to perform co-manifold learning: uncovering the underlying geometry of both the rows and the columns of a given matrix, where we focus on a missing data setting. Our unsupervised approach consists of three components. We first solve a family of optimization problems to estimate a complete matrix at multiple scales of smoothness. We then use this collection of smooth matrix estimates to compute pairwise distances on the rows and columns based on a new multi-scale metric that implicitly introduces a coupling between the rows and the columns. Finally, we construct row and column representations from these multi-scale metrics. We demonstrate that our approach outperforms competing methods in both data visualization and clustering.


title: Symmetry-Based Approach to Shape Coexistence in Nuclei

abstract: A symmetry-based approach for describing shape-coexistence, is presented in the framework of the interacting boson model of nuclei. It involves a construction of a number-conserving Hamiltonian which preserves the dynamical symmetry of selected bands associated with each shape, while breaking the symmetries in other states. The resulting structure embodies multiple partial dynamical symmetries. The procedure is applied to prolate-oblate and spherical-prolate-oblate coexistence, at and slightly away from the critical points.


title: Dirichlet problem in domains with lower dimensional boundaries

abstract: The present paper pioneers the study of the Dirichlet problem with $L^q$ boundary data for second order operators with complex coefficients in domains with lower dimensional boundaries, e.g., in $\Omega := \mathbb R^n \setminus \mathbb R^d$ with $d<n-1$. Following the first results of Guy David and the two first authors, the article introduces an appropriate degenerate elliptic operator and show that the Dirichlet problem is solvable for all $q>1$ provided that the coefficients satisfy the small Carleson norm condition. Even in the context of the classical case $d=n-1$, (the analogues of) our results are new. The conditions on the coefficients are more relaxed than the previously known ones (most notably, we do not impose any restrictions whatsoever on the first $n-1$ rows of the matrix of coefficients) and the results are more general. We establish local rather than global estimates between the square function and the non-tangential maximal function and, perhaps even more importantly, we establish new Moser-type estimates at the boundary and improve the interior ones.


title: High-order Solution Transfer between Curved Triangular Meshes

abstract: The problem of solution transfer between meshes arises frequently in computational physics, e.g. in Lagrangian methods where remeshing occurs. The interpolation process must be conservative, i.e. it must conserve physical properties, such as mass. We extend previous works --- which described the solution transfer process for straight sided unstructured meshes --- by considering high-order isoparametric meshes with curved elements. To facilitate solution transfer, we numerically integrate the product of shape functions via Green's theorem along the boundary of the intersection of two curved elements. We perform a numerical experiment and confirm the expected accuracy by transferring test fields across two families of meshes.


title: Morph: Flexible Acceleration for 3D CNN-based Video Understanding

abstract: The past several years have seen both an explosion in the use of Convolutional Neural Networks (CNNs) and the design of accelerators to make CNN inference practical. In the architecture community, the lion share of effort has targeted CNN inference for image recognition. The closely related problem of video recognition has received far less attention as an accelerator target. This is surprising, as video recognition is more computationally intensive than image recognition, and video traffic is predicted to be the majority of internet traffic in the coming years. This paper fills the gap between algorithmic and hardware advances for video recognition by providing a design space exploration and flexible architecture for accelerating 3D Convolutional Neural Networks (3D CNNs) - the core kernel in modern video understanding. When compared to (2D) CNNs used for image recognition, efficiently accelerating 3D CNNs poses a significant engineering challenge due to their large (and variable over time) memory footprint and higher dimensionality. To address these challenges, we design a novel accelerator, called Morph, that can adaptively support different spatial and temporal tiling strategies depending on the needs of each layer of each target 3D CNN. We codesign a software infrastructure alongside the Morph hardware to find good-fit parameters to control the hardware. Evaluated on state-of-the-art 3D CNNs, Morph achieves up to 3.4x (2.5x average) reduction in energy consumption and improves performance/watt by up to 5.1x (4x average) compared to a baseline 3D CNN accelerator, with an area overhead of 5%. Morph further achieves a 15.9x average energy reduction on 3D CNNs when compared to Eyeriss.