title: Step Timing Adjustment: A Step toward Generating Robust Gaits

abstract: Step adjustment for humanoid robots has been shown to improve robustness in gaits. However, step duration adaptation is often neglected in control strategies. In this paper, we propose an approach that combines both step location and timing adjustment for generating robust gaits. In this approach, step location and step timing are decided, based on feedback from the current state of the robot. The proposed approach is comprised of two stages. In the first stage, the nominal step location and step duration for the next step or a previewed number of steps are specified. In this stage which is done at the start of each step, the main goal is to specify the best step length and step duration for a desired walking speed. The second stage deals with finding the best landing point and landing time of the swing foot at each control cycle. In this stage, stability of the gaits is preserved by specifying a desired offset between the swing foot landing point and the Divergent Component of Motion (DCM) at the end of current step. After specifying the landing point of the swing foot at a desired time, the swing foot trajectory is regenerated at each control cycle to realize desired landing properties. Simulation on different scenarios shows the robustness of the generated gaits from our proposed approach compared to the case where no timing adjustment is employed.


title: Frame completions with prescribed norms: local minimizers and applications

abstract: Let $\mathcal F_0=\{f_i\}_{i\in\mathbb{I}_{n_0}}$ be a finite sequence of vectors in $\mathbb C^d$ and let $\mathbf{a}=(a_i)_{i\in\mathbb{I}_k}$ be a finite sequence of positive numbers. We consider the completions of $\cal F_0$ of the form $\cal F=(\cal F_0,\cal G)$ obtained by appending a sequence $\cal G=\{g_i\}_{i\in\mathbb{I}_k}$ of vectors in $\mathbb C^d$ such that $\|g_i\|^2=a_i$ for $i\in\mathbb{I}_k$, and endow the set of completions with the metric $d(\cal F,\tilde {\mathcal F}) =\max\{ \,\|g_i-\tilde g_i\|: \ i\in\mathbb{I}_k\}$ where $\tilde {\cal F}=(\cal F_0,\,\tilde {\cal G})$. In this context we show that local minimizers on the set of completions of a convex potential $\text{P}_\varphi$, induced by a strictly convex function $\varphi$, are also global minimizers. In case that $\varphi(x)=x^2$ then $\text{P}_\varphi$ is the so-called frame potential introduced by Benedetto and Fickus, and our work generalizes several well known results for this potential. We show that there is an intimate connection between frame completion problems with prescribed norms and frame operator distance (FOD) problems. We use this connection and our results to settle in the affirmative a generalized version of Strawn's conjecture on the FOD.


title: Interpolating sequences on the bidisk

abstract: We give a characterization of interpolating sequences for bounded analytic functions on the bidisk.


title: BBGKY chain and kinetic equations for the level dynamics in an externally perturbed quantum system

abstract: Theoretical description and simulation of large quantum coherent systems out of equilibrium remains a daunting task. Here we are developing a new approach to it based on the Pechukas-Yukawa formalism, which is especially convenient in case of an adiabatically slow external perturbation. In this formalism the dynamics of energy levels in an externally perturbed quantum system as a function of the perturbation parameter is mapped on that of a fictitious one-dimensional classical gas of particles with cubic repulsion. Equilibrium statistical mechanics of this Pechukas gas allows to reproduce the random matrix theory of energy levels. In the present work, we develop the nonequilibrium statistical mechanics of the Pechukas gas, starting with the derivation of the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) chain of equations for the appropriate generalized distribution functions. Sets of approximate kinetic equations can be consistently obtained by breaking this chain at a particular point (i.e. approximating all higher-order distribution functions by the products of the lower-order ones). When complemented by the equations for the level occupation numbers and inter-level transition amplitudes, they allow to describe the nonequilibrium evolution of the quantum state of the system, which can describe better a large quantum coherent system than the currently used approaches. In particular, we find that corrections to the factorized approximation of the distribution function scale as 1/N, where N is the number of the "Pechukas gas particles" (i.e. energy levels in the system).


title: Geometric random graphs and Rado sets in sequence spaces

abstract: We consider a random geometric graph model, where pairs of vertices are points in a metric space and edges are formed independently with fixed probability $p$ between pairs within threshold distance $\delta $. A countable dense set in a metric space is {\sl Rado} if this random model gives, with probability 1, a graph that is unique up to isomorphism. In earlier work, the first two authors proved that in finite dimensional spaces $\mathbb{R}^n$ equipped with the $\ell_{\infty}$ norm, all countable dense set satisfying a mild non-integrality condition are Rado. In this paper, we extend this result to infinite-dimensional spaces. If the underlying metric space is a separable Banach space, then we show in some cases that we can almost surely recover the Banach space from such a geometric random graph. More precisely, we show that in the sequence spaces $c$ and $c_0$, for measures $\mu$ satisfying certain conditions, $\mu^\N$-almost all countable sets are Rado. Moreover, with probability 1, in $c$ as in $c_0$, all graphs obtained from the random geometric model with a randomly chosen dense countable vertex set are isomorphic to each other. Finally, we show that representatives of the isomorphism classes obtained in this way from $c$ and $c_0$ are non-isomorphic to each other, and also non-isomorphic to their counterparts obtained from finite dimensional spaces.


title: The three point Pick problem on the Bidisk

abstract: We analyze the three point Pick interpolation problem on the bidisk


title: A simple construction of positive loops of Legendrians

abstract: We construct positive loops of Legendrian submanifolds in several instances. In particular, we partially recover G. Liu's result stating that any loose Legendrian admits a positive loop, under some mild topological assumptions on the Legendrian. Moreover, we show contractibility of the constructed loops under an extra topological assumption.


title: Hydrodynamics of bilayer membranes with diffusing transmembrane proteins

abstract: We consider the hydrodynamics of lipid bilayers containing transmembrane proteins of arbitrary shape. This biologically-motivated problem is relevant to the cell membrane, whose fluctuating dynamics play a key role in phenomena ranging from cell migration, intercellular transport, and cell communication. Using Onsager's variational principle, we derive the equations that govern the relaxation dynamics of the membrane shape, of the mass densities of the bilayer leaflets, and of the diffusing proteins' concentration. With our generic formalism, we obtain several results on membrane dynamics. We find that proteins that span the bilayer increase the intermonolayer friction coefficient. The renormalization, which can be significant, is in inverse proportion to the protein's mobility. Second, we find that asymmetric proteins couple to the membrane curvature and to the difference in monolayer densities. For practically all accessible membrane tensions ($\sigma> 10^{-8}$ N/m) we show that the protein density is the slowest relaxing variable. Furthermore, its relaxation rate decreases at small wavelengths due to the coupling to curvature. We apply our formalism to the large-scale diffusion of a concentrated protein patch. We find that the diffusion profile is not self-similar, owing to the wavevector dependence of the effective diffusion coefficient.


title: Safe and Robust Robot Maneuvers Based on Reach Control

abstract: In this paper, we investigate the synthesis of piecewise affine feedback controllers to address the problem of safe and robust controller design in robotics based on high-level controls specifications. The methodology is based on formulating the problem as a collection of reach control problems on a polytopic state space. Reach control has so far only been developed in theory and has not been tested experimentally on a real system before. Using a quadrocopter as our experimental platform, we show that these theoretical tools can achieve fast, albeit safe and robust maneuvers. In contrast to most traditional control techniques, the reach control approach does not require a predefined open-loop reference trajectory or spacial path. Experimental results on a quadrocopter show the effectiveness and robustness of this control approach. In a proof-of-concept demonstration, the reach controller is implemented in one translational direction while the other degrees of freedom are stabilized by separate controllers.


title: Charge metastability and hysteresis in the quantum Hall regime

abstract: We report simultaneous quasi-dc magnetotransport and high frequency surface acoustic wave measurements on bilayer two-dimensional electron systems in GaAs. Near strong integer quantized Hall states a strong magnetic field sweep hysteresis in the velocity of the acoustic waves is observed at low temperatures. This hysteresis indicates the presence of a metastable state with anomalously high conductivity in the interior of the sample. This non-equilibrium state is not revealed by conventional low frequency transport measurements which are dominated by dissipationless transport at the edge of the 2D system. We find that a field-cooling technique allows the equilibrium charge configuration within the interior of the sample to be established. A simple model for this behavior is discussed.


title: Final Service Provider DevOps concept and evaluation

abstract: This report presents the results of the UNIFY Service Provider DevOps activities. First, we present the final definition and assessment of the concept. SP-DevOps is realized by a combination of various functional components facilitating integrated service verification, efficient and programmable observability, and automated troubleshooting processes. Our assessment shows that SP-DevOps can help providers to reach a medium level of DevOps maturity and allows significant reduction in OPEX. Second, we focus on the evaluation of the proposed SP-DevOps components. The set of tools proposed supports ops and devs across all stages, with a focus on the deployment, operation and debugging phases, and allows to activate automated processes for operating NFV environments. Finally, we present use-cases and our demonstrators for selected process implementions, which allowed the functional validation of SP-DevOps.


title: The role of body flexibility in stroke enhancements for finite-length undulatory swimmers in viscoelastic fluids

abstract: The role of passive body dynamics on the kinematics of swimming micro-organisms in complex fluids is investigated. Asymptotic analysis of small amplitude motions of a finite-length undulatory swimmer in a Stokes-Oldroyd-B fluid is used to predict shape changes that result as body elasticity and fluid elasticity are varied. Results from the analysis are compared with numerical simulations, and the small amplitude analysis of shape changes is quantitatively accurate at both small and large amplitudes, even for strongly elastic flows. We compute a stroke-induced swimming speed that accounts for the shape changes, but not additional effects of fluid elasticity. Elastic induced shape changes lead to larger amplitude strokes for sufficiently soft swimmers in a viscoelastic fluid, and these stroke boosts can lead to swimming speed-ups, but we find that additional effects of fluid elasticity generically slow down swimmers. High amplitude strokes in strongly elastic flows lead to a qualitatively different regime in which highly concentrated elastic stresses accumulate near swimmer bodies and where dramatic slow-downs are seen.


title: Thouless energy and multifractality across the many-body localization transition

abstract: Thermal and many-body localized phases are separated by a dynamical phase transition of a new kind. We analyze the distribution of off-diagonal matrix elements of local operators across the many-body localization transition (MBLT) in a disordered spin chain, and use it to characterize the breakdown of the eigenstate thermalization hypothesis and to extract the many-body Thouless energy. We find a wide critical region around the MBLT, where Thouless energy becomes smaller than the level spacing, while matrix elements show critical dependence on the energy difference. In the same region, matrix elements, viewed as amplitudes of a fictitious wave function, exhibit strong multifractality. Our findings show that the correlation length becomes larger than the accessible system sizes in a broad range of disorder strength values, and shed light on the critical behaviour of MBL systems.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el


title: Shock Wave Interactions and the Riemann-flat Condition: The Geometry behind Metric Smoothing and the Existence of Locally Inertial Frames in General Relativity

abstract: We prove that the essential smoothness of the gravitational metric at shock waves in GR, a PDE regularity issue for weak solutions of the Einstein equations, is determined by a geometrical condition which we introduce and name the {\it Riemann-flat condition}. The Riemann-flat condition determines whether or not the essential smoothness of the gravitational metric is two full derivatives more regular than the Riemann curvature tensor. This provides a geometric framework for the open problem as to whether {\it regularity singularities} (points where the curvature is in $L^\infty$ but the essential smoothness of the gravitational metric is only Lipschitz continuous) can be created by shock wave interaction in GR, or whether metrics Lipschitz at shocks can always be smoothed one level to $C^{1,1}$ by coordinate transformation. As a corollary of the ideas we give a proof that locally inertial frames always exist in a natural sense for shock wave metrics in spherically symmetric spacetimes, independent of whether the metric itself can be smoothed to $C^{1,1}$ locally. This latter result yields an explicit procedure (analogous to Riemann Normal Coordinates in smooth spacetimes) for constructing locally inertial coordinates for Lipschitz metrics, and is a new regularity result for GR solutions constructed by the Glimm scheme.


title: Grad-CAM: Visual Explanations from Deep Networks via Gradient-based Localization

abstract: We propose a technique for producing "visual explanations" for decisions from a large class of CNN-based models, making them more transparent. Our approach - Gradient-weighted Class Activation Mapping (Grad-CAM), uses the gradients of any target concept, flowing into the final convolutional layer to produce a coarse localization map highlighting important regions in the image for predicting the concept. Grad-CAM is applicable to a wide variety of CNN model-families: (1) CNNs with fully-connected layers, (2) CNNs used for structured outputs, (3) CNNs used in tasks with multimodal inputs or reinforcement learning, without any architectural changes or re-training. We combine Grad-CAM with fine-grained visualizations to create a high-resolution class-discriminative visualization and apply it to off-the-shelf image classification, captioning, and visual question answering (VQA) models, including ResNet-based architectures. In the context of image classification models, our visualizations (a) lend insights into their failure modes, (b) are robust to adversarial images, (c) outperform previous methods on localization, (d) are more faithful to the underlying model and (e) help achieve generalization by identifying dataset bias. For captioning and VQA, we show that even non-attention based models can localize inputs. We devise a way to identify important neurons through Grad-CAM and combine it with neuron names to provide textual explanations for model decisions. Finally, we design and conduct human studies to measure if Grad-CAM helps users establish appropriate trust in predictions from models and show that Grad-CAM helps untrained users successfully discern a 'stronger' nodel from a 'weaker' one even when both make identical predictions. Our code is available at https://github.com/ramprs/grad-cam/, along with a demo at http://gradcam.cloudcv.org, and a video at youtu.be/COjUB9Izk6E.


title: An Automatic System for Acoustic Microphone Geometry Calibration based on Minimal Solvers

abstract: In this paper, robust detection, tracking and geometry estimation methods are developed and combined into a system for estimating time-difference estimates, microphone localization and sound source movement. No assumptions on the 3D locations of the microphones and sound sources are made. The system is capable of tracking continuously moving sound sources in an reverberant environment. The multi-path components are explicitly tracked and used in the geometry estimation parts. The system is based on matching between pairs of channels using GCC-PHAT. Instead of taking a single maximum at each time instant from each such pair, we select the four strongest local maxima. This produce a set of hypothesis to work with in the subsequent steps, where consistency constraints between the channels and time-continuity constraints are exploited. In the paper it demonstrated how such detections can be used to estimate microphone positions, sound source movement and room geometry. The methods are tested and verified using real data from several reverberant environments. The evaluation demonstrated accuracy in the order of few millimeters.


title: New theory of diffusive and coherent nature of optical wave via a quantum walk

abstract: We propose a new theory on a relation between diffusive and coherent nature in one dimensional wave mechanics based on a quantum walk. It is known that the quantum walk in homogeneous matrices provides the coherent property of wave mechanics. Using the recent result of a localization phenomenon in a one-dimensional quantum walk (Konno, Quantum Inf. Proc. (2010) 9, 405-418), we numerically show that the randomized localized matrices suppress the coherence and give diffusive nature.


title: Color-kinematics duality and Sudakov form factor at five loops for N=4 supersymmetric Yang-Mills theory

abstract: Using color-kinematics duality, we construct for the first time the full integrand of the five-loop Sudakov form factor in N=4 super-Yang-Mills theory, including non-planar contributions. This result also provides a first manifestation of the color-kinematics duality at five loops. The integrand is explicitly ultraviolet finite when D<26/5, coincident with the known finiteness bound for amplitudes. If the double-copy method could be applied to the form factor, this would indicate an interesting ultraviolet finiteness bound for N=8 supergravity at five loops. The result is also expected to provide an essential input for computing the five-loop non-planar cusp anomalous dimension.


title: Causal-order superposition as an enabler of free will

abstract: It is often argued that bottom-up causation under a physicalist, reductionist worldview precludes free will in the libertarian sense. On the one hand, the paradigm of classical mechanics makes determinism inescapable, while on the other, the leading models that allow a role for quantum effects are noncommittal regarding how conscious agents are supposed to translate indeterminacy into self-formed choice. Recent developments, however, not only imply that self-formed decisions are possible, but actually suggest how they might come about. The cornerstone appears to be causality superposition rather than quantum-state entanglement, as is usually assumed, and the natural arena for applying these developments is (perhaps ironically) a framework that was built without any consideration for quantum effects.


title: Readout technologies for directional WIMP Dark Matter detection

abstract: The measurement of the direction of WIMP-induced nuclear recoils is a compelling but technologically challenging strategy to provide an unambiguous signature of the detection of Galactic dark matter. Most directional detectors aim to reconstruct the dark-matter-induced nuclear recoil tracks, either in gas or solid targets. The main challenge with directional detection is the need for high spatial resolution over large volumes, which puts strong requirements on the readout technologies. In this paper we review the various detector readout technologies used by directional detectors. In particular, we summarize the challenges, advantages and drawbacks of each approach, and discuss future prospects for these technologies.

categories: physics.ins-det astro-ph.CO astro-ph.IM hep-ex


title: Second Infintesimal Neighborhoods of Projective Bundle Sections

abstract: This short paper shows that a projective bundle over a smooth variety with a section can be uniquely recovered from the embedding of the section in its second infinitesimal neighborhood within the bundle.


title: Collider Signatures of Flavorful Higgs Bosons

abstract: Motivated by our limited knowledge of the Higgs couplings to first two generation fermions, we analyze the collider phenomenology of a class of two Higgs doublet models (2HDMs) with a non-standard Yukawa sector. One Higgs doublet is mainly responsible for the masses of the weak gauge bosons and the third generation fermions, while the second Higgs doublet provides mass for the lighter fermion generations. The characteristic collider signatures of this setup differ significantly from well-studied 2HDMs with natural flavor conservation, flavor alignment, or minimal flavor violation. New production mechanisms for the heavy scalar, pseudoscalar, and charged Higgs involving second generation quarks can become dominant. The most interesting decay modes include H/A -> cc, tc, {\mu}{\mu}, {\tau}{\mu} and H+ -> cb, cs, {\mu}{\nu}. Searches for low mass di-muon resonances are currently among the best probes of the heavy Higgs bosons in this setup.


title: Statistics of Dark Matter Substructure: III. Halo-to-Halo Variance

abstract: We present a study of unprecedented statistical power regarding the halo-to-halo variance of dark matter substructure. Using a combination of N-body simulations and a semi-analytical model, we investigate the variance in subhalo mass fractions and subhalo occupation numbers, with an emphasis on how these statistics scale with halo formation time. We demonstrate that the subhalo mass fraction, f_sub, is mainly a function of halo formation time, with earlier forming haloes having less substructure. At fixed formation redshift, the average f_sub is virtually independent of halo mass, and the mass dependence of f_sub is therefore mainly a manifestation of more massive haloes assembling later. We compare observational constraints on f_sub from gravitational lensing to our model predictions and simulation results. Although the inferred f_sub are substantially higher than the median LCDM predictions, they fall within the 95th percentile due to halo-to-halo variance. We show that while the halo occupation distribution of subhaloes, P(N|M), is super-Poissonian for large <N>, a well established result, it becomes sub-Poissonian for <N> < 2. Ignoring the non-Poissonity results in systematic errors of the clustering of galaxies of a few percent, and with a complicated scale- and luminosity-dependence. Earlier-formed haloes have P(N|M) closer to a Poisson distribution, suggesting that the dynamical evolution of subhaloes drives the statistics towards Poissonian. Contrary to a recent claim, the non-Poissonity of subhalo occupation statistics does not vanish by selecting haloes with fixed mass and fixed formation redshift. Finally, we use subhalo occupation statistics to put loose constraints on the mass and formation redshift of the Milky Way halo. Using observational constraints on the V_max of the most massive satellites, we infer that 0.25<M_vir/10^12M_sun/h<1.4 and 0.1<z_f<1.4 at 90% confidence.


title: Robust sparse image reconstruction of radio interferometric observations with purify

abstract: Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, to realise these goals significant challenges in image and data processing need to be overcome. The standard methods in radio interferometry for reconstructing images, such as CLEAN, have served the community well over the last few decades and have survived largely because they are pragmatic. However, they produce reconstructed inter\-ferometric images that are limited in quality and scalability for big data. In this work we apply and evaluate alternative interferometric reconstruction methods that make use of state-of-the-art sparse image reconstruction algorithms motivated by compressive sensing, which have been implemented in the PURIFY software package. In particular, we implement and apply the proximal alternating direction method of multipliers (P-ADMM) algorithm presented in a recent article. First, we assess the impact of the interpolation kernel used to perform gridding and degridding on sparse image reconstruction. We find that the Kaiser-Bessel interpolation kernel performs as well as prolate spheroidal wave functions, while providing a computational saving and an analytic form. Second, we apply PURIFY to real interferometric observations from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA) and find images recovered by PURIFY are higher quality than those recovered by CLEAN. Third, we discuss how PURIFY reconstructions exhibit additional advantages over those recovered by CLEAN. The latest version of PURIFY, with developments presented in this work, is made publicly available.


title: Radiative shocks create environments for dust formation in novae

abstract: Classical novae commonly show evidence of rapid dust formation within months of the outburst. However, it is unclear how molecules and grains are able to condense within the ejecta, given the potentially harsh environment created by ionizing radiation from the white dwarf. Motivated by the evidence for powerful radiative shocks within nova outflows, we propose that dust formation occurs within the cool, dense shell behind these shocks. We incorporate a simple molecular chemistry network and classical nucleation theory with a model for the thermodynamic evolution of the post-shock gas, in order to demonstrate the formation of both carbon and forsterite ($\rm Mg_2SiO_4$) grains. The high densities due to radiative shock compression ($n \sim 10^{14}$ cm$^{-3}$) result in CO saturation and rapid dust nucleation. Grains grow efficiently to large sizes $\gtrsim 0.1\mu$m, in agreement with IR observations of dust-producing novae, and with total dust masses sufficient to explain massive extinction events such as V705 Cas. As in dense stellar winds, dust formation is CO-regulated, with carbon-rich flows producing carbon-rich grains and oxygen-rich flows primarily forming silicates. CO is destroyed by non-thermal particles accelerated at the shock, allowing additional grain formation at late times, but the efficiency of this process appears to be low. Given observations showing that individual novae produce both carbonaceous and silicate grains, we concur with previous works attributing this bimodality to chemical heterogeneity of the ejecta. Nova outflows are diverse and inhomogeneous, and the observed variety of dust formation events can be reconciled by different abundances, the range of shock properties, and the observer viewing angle. The latter may govern the magnitude of extinction, with the deepest extinction events occurring for observers within the binary equatorial plane.


title: The Search for Failed Supernovae with the Large Binocular Telescope: Constraints from 7 Years of Data

abstract: We report updated results for the first 7 years of our program to monitor 27 galaxies within 10 Mpc using the Large Binocular Telescope to search for failed supernovae (SNe) -- core-collapses of massive stars that form black holes without luminous supernovae. In the new data, we identify no new compelling candidates and confirm the existing candidate. Given the 6 successful core-collapse SNe in the sample and one likely failed SN, the implied fraction of core-collapses that result in failed SNe is $f = 0.14^{+0.33}_{-0.10}$ at 90% confidence. If the current candidate is a failed SN, the fraction of failed SN naturally explains the missing high-mass red supergiant SN progenitors and the black hole mass function. If the current candidate is ultimately rejected, the data implies a 90% confidence upper limit on the failed SN fraction of $f < 0.35$.


title: Warped K\"ahler potentials and fluxes

abstract: The four-dimensional effective theory for type IIB warped flux compactifications proposed in [1] is completed by taking into account the backreaction of the K\"ahler moduli on the three-form fluxes. The only required modification consists in a flux-dependent contribution to the chiral fields parametrising the K\"ahler moduli. The resulting supersymmetric effective theory satisfies the no-scale condition and consistently combines previous partial results present in the literature. Similar results hold for M-theory warped compactifications on Calabi-Yau fourfolds, whose effective field theory and K\"ahler potential are also discussed.


title: Ultraviolet Halos Around Spiral Galaxies. I. Morphology

abstract: We examine ultraviolet halos around a sample of highly inclined galaxies within 25 Mpc to measure their morphology and luminosity. Despite contamination from galactic light scattered into the wings of the point-spread function, we find that UV halos occur around each galaxy in our sample. Around most galaxies the halos form a thick, diffuse disk-like structure, but starburst galaxies with galactic superwinds have qualitatively different halos that are more extensive and have filamentary structure. The spatial coincidence of the UV halos above star-forming regions, the lack of consistent association with outflows or extraplanar ionized gas, and the strong correlation between the halo and galaxy UV luminosity suggest that the UV light is an extragalactic reflection nebula. UV halos may thus represent 1-10 million solar masses of dust within 2-10 kpc of the disk, whose properties may change with height in starburst galaxies.


title: Advanced LIGO Constraints on Neutron Star Mergers and R-Process Sites

abstract: The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron star (NS-NS) mergers that is needed in order to reproduce the observed [Eu/Fe] vs [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution (DTD) functions for NS-NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS-NS mergers and black hole - neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS-NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 Msun, the number of NS-NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS-NS merger ejecta, and low Fe yields for massive stars. Population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide important constraints on the actual local NS-NS merger rate, will provide insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the Universe.


title: Anisotropic transport of normal metal-barrier-normal metal junctions in monolayer phosphorene

abstract: We study transport properties of a phosphorene monolayer in the presence of single and multiple potential barriers of height $U_0$ and width $d$, using both continuum and microscopic lattice models, and show that the nature of electron transport along its armchair edge ($x$ direction) is qualitatively different from its counterpart in both conventional two-dimensional electron gas with Schr\"odinger-like quasiparticles and graphene or surfaces of topological insulators hosting massless Dirac quasiparticles. We show that the transport, mediated by massive Dirac electrons, allows one to achieve collimated quasiparticle motion along $x$ and thus makes monolayer phosphorene an ideal experimental platform for studying Klein paradox. We study the dependence of the tunneling conductance $G \equiv G_{xx}$ as a function of $d$ and $U_0$, and demonstrate that for a given applied voltage $V$ its behavior changes from oscillatory to decaying function of $d$ for a range of $U_0$ with finite non-zero upper and lower bounds, and provide analytical expression for these bounds within which $G$ decays with $d$. We contrast such behavior of $G$ with that of massless Dirac electrons in graphene and also with that along the zigzag edge ($y$ direction) in phosphorene where the quasiparticles obey an effective Schr\"odinger equation at low energy. We also study transport through multiple barriers along $x$ and demonstrate that these properties hold for transport through multiple barrier as well. Finally, we suggest concrete experiments which may verify our theoretical predictions.


title: On the frequency of star-forming galaxies in the vicinity of powerful AGNs: The case of SMM J04135+10277

abstract: (Abridged) In the last decade several massive molecular gas reservoirs were found <100 kpc distance from active galactic nuclei (AGNs), residing in gas-rich companion galaxies. The study of AGN-gas-rich companion systems opens the opportunity to determine whether the stellar mass of massive local galaxies was formed in their host after a merger event or outside of their host galaxy in a close starbursting companion and later incorporated via mergers. We study the quasar-companion galaxy system of SMM J04135+10277 (z=2.84) and investigate the expected frequency of quasar-starburst galaxy pairs at high redshift using a cosmological galaxy formation model. We use archive data and new APEX ArTeMiS data to construct and model the spectral energy distribution of SMM J04135. We also carry out a comprehensive analysis of the cosmological galaxy formation model GALFORM with the aim of characterising how typical the system of SMM J04135 is and whether quasar-star-forming galaxy pairs may constitute an important stage in galaxy evolution. The companion galaxy of SMM J04135 is a heavily dust-obscured starburst galaxy with a median star formation rate (SFR) of $700\,\rm{M_{\odot}\,yr^{-1}}$, median dust mass of $5.1\times 10^9\,\rm{M_{\odot}}$ and median dust luminosity of $\textrm 9.3 \times 10^{12}\,\rm{L_{\odot}}$. Our simulations, performed at z=2.8, suggest that SMM J04135 is not unique. In fact, at a distance of <100 kpc, 22% of our simulated quasar sample have at least one companion galaxy of a stellar mass $>10^8\, \rm{M_{\odot}}$, and 0.3% have at least one highly star-forming companion ($\rm{SFR}>100\,\rm{M_{\odot}\,yr^{-1}}$). Our results suggest that quasar-gas-rich companion galaxy systems are common phenomena in the early Universe and the high incidence of companions makes the study of such systems crucial to understand the growth and hierarchical build-up of galaxies and black holes.


title: A second look at transition amplitudes in (2+1)-dimensional causal dynamical triangulations

abstract: Studying transition amplitudes in (2+1)-dimensional causal dynamical triangulations, Cooperman and Miller discovered speculative evidence for Lorentzian quantum geometries emerging from its Euclidean path integral. On the basis of this evidence, Cooperman and Miller conjectured that Lorentzian de Sitter spacetime, not Euclidean de Sitter space, dominates the ground state of the quantum geometry of causal dynamical triangulations on large scales, a scenario akin to that of the Hartle-Hawking no-boundary proposal in which Lorentzian spacetimes dominate a Euclidean path integral. We argue against this conjecture: we propose a more straightforward explanation of their findings, and we proffer evidence for the Euclidean nature of these seemingly Lorentzian quantum geometries. This explanation reveals another manner in which the Euclidean path integral of causal dynamical triangulations behaves correctly in its semiclassical limit--the implementation and interaction of multiple constraints.


title: The SCUBA-2 Cosmology Legacy Survey: the nature of bright submm galaxies from 2 deg2 of 850-um imaging

abstract: We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.


title: Profinite rigidity and surface bundles over the circle

abstract: If $M$ is a compact 3-manifold whose first betti number is 1, and $N$ is a compact 3-manifold such that $\pi_1N$ and $\pi_1M$ have the same finite quotients, then $M$ fibres over the circle if and only if $N$ does. We prove that groups of the form $F_2\rtimes\mathbb{Z}$ are distinguished from one another by their profinite completions. Thus, regardless of betti number, if $M$ and $N$ are punctured torus bundles over the circle and $M$ is not homeomorphic to $N$, then there is a finite group $G$ such that one of $\pi_1M$ and $\pi_1N$ maps onto $G$ and the other does not.


title: Colours, Star formation Rates, and Environments of Star forming and Quiescent Galaxies at the Cosmic Noon

abstract: We analyse the star formation rates (SFRs), colours, and dust extinctions of galaxies in massive (10^12.5-10^13.5 Msun) haloes at z~2 in high-resolution, cosmological zoom-in simulations as part of the Feedback in Realistic Environments (FIRE) project. The simulations do not model feedback from active galactic nuclei (AGN) but reproduce well the observed relations between stellar and halo mass and between stellar mass and SFR. About half (a third) of the simulated massive galaxies (massive central galaxies) at z~2 have broad-band colours classifying them as 'quiescent', and the fraction of quiescent centrals is steeply decreasing towards higher redshift, in agreement with observations. The progenitors of z~2 quiescent central galaxies are, on average, more massive, have lower specific SFRs, and reside in more massive haloes than the progenitors of similarly massive star forming centrals. The simulations further predict a morphological mix of galaxies that includes disk-dominated, irregular, and early-type galaxies. However, our simulations do not reproduce the reddest of the quiescent galaxies observed at z~2. We also do not find evidence for a colour bimodality, but are limited by our modest sample size. In our simulations, the star formation activity of central galaxies of moderate mass (Mstar~10^10-10^11 Msun) is affected by a combination of two distinct physical processes. Outflows powered by stellar feedback result in a short-lived (<100 Myr), but almost complete, suppression of star formation activity after which many galaxies quickly recover and continue to form stars at normal rates. In addition, galaxies residing in slowly growing haloes tend to experience a moderate reduction of their SFRs ('cosmological starvation'). The relative importance of these processes and AGN feedback is uncertain and will be explored in future work.


title: A Kernel-based Lagrangian Method for Imperfectly-mixed Chemical Reactions

abstract: Current Lagrangian (particle-tracking) algorithms used to simulate diffusion-reaction equations must employ a certain number of particles to properly emulate the system dynamics---particularly for imperfectly-mixed systems. The number of particles is tied to the statistics of the initial concentration fields of the system at hand. Systems with shorter-range correlation and/or smaller concentration variance require more particles, potentially limiting the computational feasibility of the method. For the well-known problem of bimolecular reaction, we show that using kernel-based, rather than Dirac-delta, particles can significantly reduce the required number of particles. We derive the fixed width of a Gaussian kernel for a given reduced number of particles that analytically eliminates the error between kernel and Dirac solutions at any specified time. We also show how to solve for the fixed kernel size by minimizing the squared differences between solutions over any given time interval. Numerical results show that the width of the kernel should be kept below about 12% of the domain size, and that the analytic equations used to derive kernel width suffer significantly from the neglect of higher-order moments. The simulations with a kernel width given by least squares minimization perform better than those made to match at one specific time. A heuristic time-variable kernel size, based on the previous results, performs on a par with the least squares fixed kernel size.


title: Equality of Opportunity in Supervised Learning

abstract: We propose a criterion for discrimination against a specified sensitive attribute in supervised learning, where the goal is to predict some target based on available features. Assuming data about the predictor, target, and membership in the protected group are available, we show how to optimally adjust any learned predictor so as to remove discrimination according to our definition. Our framework also improves incentives by shifting the cost of poor classification from disadvantaged groups to the decision maker, who can respond by improving the classification accuracy. In line with other studies, our notion is oblivious: it depends only on the joint statistics of the predictor, the target and the protected attribute, but not on interpretation of individualfeatures. We study the inherent limits of defining and identifying biases based on such oblivious measures, outlining what can and cannot be inferred from different oblivious tests. We illustrate our notion using a case study of FICO credit scores.


title: Indoor Space Recognition using Deep Convolutional Neural Network: A Case Study at MIT Campus

abstract: In this paper, we propose a robust and parsimonious approach using Deep Convolutional Neural Network (DCNN) to recognize and interpret interior space. DCNN has achieved incredible success in object and scene recognition. In this study we design and train a DCNN to classify a pre-zoning indoor space, and from a single phone photo to recognize the learned space features, with no need of additional assistive technology. We collect more than 600,000 images inside MIT campus buildings to train our DCNN model, and achieved 97.9% accuracy in validation dataset and 81.7% accuracy in test dataset based on spatial-scale fixed model. Furthermore, the recognition accuracy and spatial resolution can be potentially improved through multiscale classification model. We identify the discriminative image regions through Class Activating Mapping (CAM) technique, to observe the model's behavior in how to recognize space and interpret it in an abstract way. By evaluating the results with misclassification matrix, we investigate the visual spatial feature of interior space by looking into its visual similarity and visual distinctiveness, giving insights into interior design and human indoor perception and wayfinding research. The contribution of this paper is threefold. First, we propose a robust and parsimonious approach for indoor navigation using DCNN. Second, we demonstrate that DCNN also has a potential capability in space feature learning and recognition, even under severe appearance changes. Third, we introduce a DCNN based approach to look into the visual similarity and visual distinctiveness of interior space.


title: Automatic chemical design using a data-driven continuous representation of molecules

abstract: We report a method to convert discrete representations of molecules to and from a multidimensional continuous representation. This model allows us to generate new molecules for efficient exploration and optimization through open-ended spaces of chemical compounds. A deep neural network was trained on hundreds of thousands of existing chemical structures to construct three coupled functions: an encoder, a decoder and a predictor. The encoder converts the discrete representation of a molecule into a real-valued continuous vector, and the decoder converts these continuous vectors back to discrete molecular representations. The predictor estimates chemical properties from the latent continuous vector representation of the molecule. Continuous representations allow us to automatically generate novel chemical structures by performing simple operations in the latent space, such as decoding random vectors, perturbing known chemical structures, or interpolating between molecules. Continuous representations also allow the use of powerful gradient-based optimization to efficiently guide the search for optimized functional compounds. We demonstrate our method in the domain of drug-like molecules and also in the set of molecules with fewer that nine heavy atoms.


title: Non-parametric Morphologies of Mergers in the Illustris Simulation

abstract: We study non-parametric morphologies of mergers events in a cosmological context, using the Illustris project. We produce mock g-band images comparable to observational surveys from the publicly available Illustris simulation idealized mock images at $z=0$. We then measure non parametric indicators: asymmetry, Gini, $M_{20}$, clumpiness and concentration for a set of galaxies with $M_* >10^{10}$ M$_\odot$. We correlate these automatic statistics with the recent merger history of galaxies and with the presence of close companions. Our main contribution is to assess in a cosmological framework, the empirically derived non-parametric demarcation line and average time-scales used to determine the merger rate observationally. We found that 98 per cent of galaxies above the demarcation line have a close companion or have experienced a recent merger event. On average, merger signatures obtained from the $G-M_{20}$ criteria anticorrelate clearly with the elapsing time to the last merger event. We also find that the asymmetry correlates with galaxy pair separation and relative velocity, exhibiting the larger enhancements for those systems with pair separations $d < 50$ h$^{-1}$ kpc and relative velocities $V < 350$ km s$^{-1}$. We find that the $G-M_{20}$ is most sensitive to recent mergers ($\sim0.14$ Gyr) and to ongoing mergers with stellar mass ratios greater than 0.1. For this indicator, we compute a merger average observability time-scale of $\sim0.2$ Gyr, in agreement with previous results and demonstrate that the morphologically derived merger rate recovers the intrinsic total merger rate of the simulation and the merger rate as a function of stellar mass.


title: A Lefschetz Hyperplane Theorem for non-Archimedean Jacobians

abstract: We establish a Lefschetz hyperplane theorem for the Berkovich analytifications of Jacobians of curves over an algebraically closed non-Archimedean field. Let $J$ be the Jacobian of a curve $X$, and let $W_d \subset J$ be the locus of effective divisor classes of degree $d$. We show that the pair $(J^{an},W_d^{an})$ is $d$-connected, and thus in particular the inclusion of the analytification of the theta divisor $\Theta^{an}$ into $J^{an}$ satisfies a Lefschetz hyperplane theorem for $\mathbb{Z}$-cohomology groups and homotopy groups. A key ingredient in our proof is a generalization, over arbitrary characteristics and allowing arbitrary singularities on the base, of a result of Brown and Foster for the homotopy type of analytic projective bundles.


title: On the Structure Theorem of Clifford Algebras

abstract: In this paper, theory and construction of spinor representations of real Clifford algebras $\cl_{p,q}$ in minimal left ideals are reviewed. Connection with a general theory of semisimple rings is shown. The actual computations can be found in, for example, [2].


title: Charged-particle multiplicity and transverse-energy distribution using the Weibull-Glauber approach in heavy-ion collisions

abstract: The charged-particle multiplicity distribution and the transverse-energy distribution measured in heavy-ion collisions at top RHIC and LHC energies are described using the two-component model approach based on a convolution of the Monte Carlo Glauber model with the Weibull model for particle production. The model successfully describes the multiplicity and transverse-energy distribution of minimum-bias collision data for a wide range of energies. The Weibull-Glauber model can be used to determine the centrality classes in heavy-ion collisions as an alternative to the conventional negative binomial distribution (NBD)-Glauber approach.


title: Lopsidependency in the Moser-Tardos framework: Beyond the Lopsided Lovasz Local Lemma

abstract: The Lopsided Lov\'{a}sz Local Lemma (LLLL) is a powerful probabilistic principle which has been used in a variety of combinatorial constructions. While originally a general statement about probability spaces, it has recently been transformed into a variety of polynomial-time algorithms. The resampling algorithm of Moser & Tardos (2010) is the most well-known example of this. A variety of criteria have been shown for the LLLL; the strongest possible criterion was shown by Shearer, and other criteria which are easier to use computationally have been shown by Bissacot et al (2011), Pegden (2014), Kolipaka & Szegedy (2011), and Kolipaka, Szegedy, Xu (2012). We show a new criterion for the Moser-Tardos algorithm to converge. This criterion is stronger than the LLLL criterion; this is possible because it does not apply in the same generality as the original LLLL; yet, it is strong enough to cover many applications of the LLLL in combinatorics. We show a variety of new bounds and algorithms. A noteworthy application is for $k$-SAT, with bounded occurrences of variables. As shown in Gebauer, Sz\'{a}bo, and Tardos (2011), a $k$-SAT instance in which every variable appears $L \leq \frac{2^{k+1}}{e (k+1)}$ times, is satisfiable. Although this bound is asymptotically tight (in $k$), we improve it to $L \leq \frac{2^{k+1} (1 - 1/k)^k}{k-1} - \frac{2}{k}$ which can be significantly stronger when $k$ is small. We introduce a new parallel algorithm for the LLLL. While Moser & Tardos described a simple parallel algorithm for the Lov\'{a}sz Local Lemma, and described a simple sequential algorithm for a form of the Lopsided Lemma, they were not able to combine the two. Our new algorithm applies in nearly all settings in which the sequential algorithm works --- this includes settings covered by our new stronger LLLL criterion.


title: Dual symmetry in a generalized Maxwell theory

abstract: We examine Podolsky's electrodynamics, which is noninvariant under the usual duality transformation. We deduce a generalization of Hodge's star duality, which leads to a dual gauge field and restores to a certain extent the dual symmetry. The model becomes fully dual symmetric asymptotically when it reduces to the Maxwell theory. We argue that this strict dual symmetry directly implies the existence of the basic invariants of the electromagnetic fields.


title: Higher Dimensional Generalizations of the SYK Model

abstract: We discuss a 1+1 dimensional generalization of the Sachdev-Ye-Kitaev model. The model contains $N$ Majorana fermions at each lattice site with a nearest-neighbour hopping term. The SYK random interaction is restricted to low momentum fermions of definite chirality within each lattice site. This gives rise to an ordinary 1+1 field theory above some energy scale and a low energy SYK-like behavior. We exhibit a class of low-pass filters which give rise to a rich variety of hyperscaling behaviour in the IR. We also discuss another set of generalizations which describes probing an SYK system with an external fermion, together with the new scaling behavior they exhibit in the IR.


title: Large conformal metrics with prescribed scalar curvature

abstract: Let $(M,g)$ be an $n-$dimensional compact Riemannian manifold. Let $h$ be a smooth function on $M$ and assume that it has a critical point $\xi\in M$ such that $h(\xi)=0$ and which satisfies a suitable flatness assumption. We are interested in finding conformal metrics $g_\lambda=u_\lambda^\frac4{n-2}g$, with $u>0$, whose scalar curvature is the prescribed function $h_\lambda=\lambda^2+h$, where $\lambda$ is a small parameter. In the positive case, i.e. when the scalar curvature $R_g$ is strictly positive, we find a family of bubbling metrics $g_\lambda$, where $u_\lambda$ blows-up at the point $\xi$ and approaches zero far from $\xi$ as $\lambda$ goes to zero. In the general case, if in addition we assume that there exists a non-degenerate conformal metric $g_0=u_0^\frac4{n-2}g$, with $u_0>0$, whose scalar curvature is equal to $h$, then there exists a bounded family of conformal metrics $g_{0,\lambda}=u_{0,\lambda}^\frac4{n-2}g$, with $u_{0,\lambda}>0$, which satisfies $u_{0,\lambda}\to u_0$ uniformly as $\lambda\to 0$. Here, we build a second family of bubbling metrics $g_\lambda$, where $u_\lambda$ blows-up at the point $\xi$ and approaches $u_0$ far from $\xi$ as $\lambda$ goes to zero. In particular, this shows that this problem admits more than one solution.


title: Diverse Beam Search: Decoding Diverse Solutions from Neural Sequence Models

abstract: Neural sequence models are widely used to model time-series data. Equally ubiquitous is the usage of beam search (BS) as an approximate inference algorithm to decode output sequences from these models. BS explores the search space in a greedy left-right fashion retaining only the top-B candidates - resulting in sequences that differ only slightly from each other. Producing lists of nearly identical sequences is not only computationally wasteful but also typically fails to capture the inherent ambiguity of complex AI tasks. To overcome this problem, we propose Diverse Beam Search (DBS), an alternative to BS that decodes a list of diverse outputs by optimizing for a diversity-augmented objective. We observe that our method finds better top-1 solutions by controlling for the exploration and exploitation of the search space - implying that DBS is a better search algorithm. Moreover, these gains are achieved with minimal computational or memory over- head as compared to beam search. To demonstrate the broad applicability of our method, we present results on image captioning, machine translation and visual question generation using both standard quantitative metrics and qualitative human studies. Further, we study the role of diversity for image-grounded language generation tasks as the complexity of the image changes. We observe that our method consistently outperforms BS and previously proposed techniques for diverse decoding from neural sequence models.


title: Quantum Cellular Automata Models for General Dirac Equation

abstract: The goal of this study is to provide an exact unitary quantum cellular automata that, under discrete time steps, converges towards the Generalized Dirac Equation (GDE) in the continuum limit. The evolutionary rules for such a single particle walk are discussed in this paper, and it is shown that this quantum celluar automata will maintain similar properties to the GDE.


title: Interplay of Site and Bond Electron-Phonon Coupling in One Dimension

abstract: The interplay of bond and charge correlations is studied in a one-dimensional model with both Holstein and Su-Schrieffer-Heeger (SSH) couplings to quantum phonons. The problem is solved exactly by quantum Monte Carlo simulations. If one of the couplings dominates, the ground state is a Peierls insulator with long-range bond or charge order. At weak coupling, the results suggest a spin-gapped and repulsive metallic phase arising from the competing order parameters and lattice fluctuations. Such a phase is absent from the pure SSH model even for quantum phonons. At strong coupling, evidence for a continuous transition between the two Peierls states is presented.


title: The Stochastic Topic Block Model for the Clustering of Vertices in Networks with Textual Edges

abstract: Due to the significant increase of communications between individuals via social media (Facebook, Twitter, Linkedin) or electronic formats (email, web, e-publication) in the past two decades, network analysis has become a unavoidable discipline. Many random graph models have been proposed to extract information from networks based on person-to-person links only, without taking into account information on the contents. This paper introduces the stochastic topic block model (STBM), a probabilistic model for networks with textual edges. We address here the problem of discovering meaningful clusters of vertices that are coherent from both the network interactions and the text contents. A classification variational expectation-maximization (C-VEM) algorithm is proposed to perform inference. Simulated data sets are considered in order to assess the proposed approach and to highlight its main features. Finally, we demonstrate the effectiveness of our methodology on two real-word data sets: a directed communication network and a undirected co-authorship network.


title: The Calabi metric and desingularization of Einstein orbifolds

abstract: Consider an Einstein orbifold $(M_0,g_0)$ of real dimension $2n$ having a singularity with orbifold group the cyclic group of order $n$ in ${\rm{SU}}(n)$ which is generated by an $n$th root of unity times the identity. Existence of a Ricci-flat K\"ahler ALE metric with this group at infinity was shown by Calabi. There is a natural "approximate" Einstein metric on the desingularization of $M_0$ obtained by replacing a small neighborhood of the singular point of the orbifold with a scaled and truncated Calabi metric. In this paper, we identify the first obstruction to perturbing this approximate Einstein metric to an actual Einstein metric. If $M_0$ is compact, we can use this to produce examples of Einstein orbifolds which do not admit any Einstein metric in a neighborhood of the natural approximate Einstein metric on the desingularization. In the case that $(M_0,g_0)$ is asymptotically hyperbolic Einstein and non-degenerate, we show that if the first obstruction vanishes, then there does in fact exist an asymptotically hyperbolic Einstein metric on the desingularization. We also obtain a non-existence result in the asymptotically hyperbolic Einstein case, provided that the obstruction does not vanish. This work extends a construction of Biquard in the case $n =2$, in which case the Calabi metric is also known as the Eguchi-Hanson metric, but there are some key points for which the higher-dimensional case differs.


title: Kinetics of Aggregation with Choice

abstract: We generalize the ordinary aggregation process to allow for choice. In ordinary aggregation, two random clusters merge and form a larger aggregate. In our implementation of choice, a target cluster and two candidate clusters are randomly selected, and the target cluster merges with the larger of the two candidate clusters. We study the long-time asymptotic behavior, and find that as in ordinary aggregation, the size density adheres to the standard scaling form. However, aggregation with choice exhibits a number of novel features. First, the density of the smallest clusters exhibits anomalous scaling. Second, both the small-size and the large-size tails of the density are overpopulated, at the expense of the density moderate-size clusters. We also study the complementary case where the smaller candidate clusters participates in the aggregation process, and find abundance of moderate clusters at the expense of small and large clusters. Additionally, we investigate aggregation processes with choice among multiple candidate clusters, and a symmetric implementation where the choice is between two pairs of clusters.


title: Local Euler obstructions of toric varieties

abstract: We use Matsui and Takeuchi's formula for toric A-discriminants to give algorithms for computing local Euler obstructions and dual degrees of toric surfaces and 3-folds. In particular, we consider weighted projective spaces. As an application we give counterexamples to a conjecture by Matsui and Takeuchi. As another application we recover the well-known fact that the only defective normal toric surfaces are cones.


title: ResearchDoom and CocoDoom: Learning Computer Vision with Games

abstract: In this short note we introduce ResearchDoom, an implementation of the Doom first-person shooter that can extract detailed metadata from the game. We also introduce the CocoDoom dataset, a collection of pre-recorded data extracted from Doom gaming sessions along with annotations in the MS Coco format. ResearchDoom and CocoDoom can be used to train and evaluate a variety of computer vision methods such as object recognition, detection and segmentation at the level of instances and categories, tracking, ego-motion estimation, monocular depth estimation and scene segmentation. The code and data are available at http://www.robots.ox.ac.uk/~vgg/research/researchdoom.


title: Model Reduction of Linear Multi-Agent Systems by Clustering and Associated $\mathcal{H}_2$- and $\mathcal{H}_\infty$-Error Bounds

abstract: In this paper, we study a model reduction technique for leader-follower networked multi-agent systems defined on weighted, undirected graphs with arbitrary linear multivariable agent dynamics. In the network graph of this network, nodes represent the agents and edges represent communication links between the agents. Only the leaders in the network receive an external input, the followers only exchange information with their neighbors. The reduced network is obtained by partitioning the set of nodes into disjoint sets, called clusters, and associating with each cluster a single, new, node in a reduced network graph. The resulting reduced network has a weighted, symmetric, directed network graph, and inherits some of the structure of the original network. We establish a priori upper bounds on the $\mathcal{H}_2$ and $\mathcal{H}_\infty$ model reduction error for the special case that the graph partition is almost equitable. These upper bounds depend on the Laplacian eigenvalues of the original and reduced network, an auxiliary system associated with the agent dynamics, and the number of nodes that belong to the same clusters as the leaders in the network. Finally, we consider the problem of obtaining a priori upper bounds if we cluster using arbitrary, possibly non almost equitable, partitions.


title: Anomalous metapopulation dynamics on scale-free networks

abstract: We model transport of individuals across a heterogeneous scale-free network where a few weakly connected nodes exhibit heavy-tailed residence times. Using the empirical law Axiom of Cumulative Inertia and fractional analysis we show that `anomalous cumulative inertia' overpowers highly connected nodes in attracting network individuals. This fundamentally challenges the classical result that individuals tend to accumulate in high-order nodes. The derived residence time distribution has a non-trivial U-shape which we encounter empirically across human residence and employment times.

categories: physics.soc-ph cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech physics.data-an


title: Algorithmic aspects of branched coverings IV/V. Expanding maps

abstract: Thurston maps are branched self-coverings of the sphere whose critical points have finite forward orbits. We give combinatorial and algebraic characterizations of Thurston maps that are isotopic to expanding maps as "Levy-free" maps and as maps with "contracting biset". We prove that every Thurston map decomposes along a unique minimal multicurve into Levy-free and finite-order pieces, and this decomposition is algorithmically computable. Each of these pieces admits a geometric structure. We apply these results to matings of post-critically finite polynomials, extending a criterion by Mary Rees and Tan Lei: they are expanding if and only if they do not admit a cycle of periodic rays.


title: Symmetry-broken states in a system of interacting bosons on a two-leg ladder with a uniform Abelian gauge field

abstract: We study the quantum phases of bosons with repulsive contact interactions on a two-leg ladder in the presence of a uniform Abelian gauge field. The model realizes many interesting states, including Meissner phases, vortex-fluids, vortex-lattices, charge-density-waves and the biased-ladder phase. Our work focuses on the subset of these states that break a discrete symmetry. We use density matrix renormalization group simulations to demonstrate the existence of three vortex-lattice states at different vortex densities and we characterize the phase transitions from these phases into neighboring states. Furthermore, we provide an intuitive explanation of the chiral-current reversal effect that is tied to some of these vortex lattices. We also study a charge-density-wave state that exists at 1/4 particle filling at large interaction strengths and flux values close to half a flux quantum. By changing the system parameters, this state can transition into a completely gapped vortex-lattice Mott-insulating state. We elucidate the stability of these phases against nearest-neighbor interactions on the rungs of the ladder relevant for experimental realizations with a synthetic lattice dimension. A charge-density-wave state at 1/3 particle filling can be stabilized for flux values close to half a flux-quantum and for very strong on-site interactions in the presence of strong repulsion on the rungs. Finally, we analytically describe the emergence of these phases in the low-density regime, and, in particular, we obtain the boundaries of the biased-ladder phase, i.e., the phase that features a density imbalance between the legs. We make contact to recent quantum-gas experiments that realized related models and discuss signatures of these quantum states in experimentally accessible observables.


title: A convex framework for high-dimensional sparse Cholesky based covariance estimation

abstract: Covariance estimation for high-dimensional datasets is a fundamental problem in modern day statistics with numerous applications. In these high dimensional datasets, the number of variables p is typically larger than the sample size n. A popular way of tackling this challenge is to induce sparsity in the covariance matrix, its inverse or a relevant transformation. In particular, methods inducing sparsity in the Cholesky pa- rameter of the inverse covariance matrix can be useful as they are guaranteed to give a positive definite estimate of the covariance matrix. Also, the estimated sparsity pattern corresponds to a Directed Acyclic Graph (DAG) model for Gaussian data. In recent years, two useful penalized likelihood methods for sparse estimation of this Cholesky parameter (with no restrictions on the sparsity pattern) have been developed. How- ever, these methods either consider a non-convex optimization problem which can lead to convergence issues and singular estimates of the covariance matrix when p > n, or achieve a convex formulation by placing a strict constraint on the conditional variance parameters. In this paper, we propose a new penalized likelihood method for sparse estimation of the inverse covariance Cholesky parameter that aims to overcome some of the shortcomings of current methods, but retains their respective strengths. We ob- tain a jointly convex formulation for our objective function, which leads to convergence guarantees, even when p > n. The approach always leads to a positive definite and symmetric estimator of the covariance matrix. We establish high-dimensional estima- tion and graph selection consistency, and also demonstrate finite sample performance on simulated/real data.


title: The late-type stellar density profile in the Galactic Center: A statistical approach

abstract: The late-type stellar population in the Galactic Center was first predicted to reside in a dynamically relaxed cusp (power law slope ranging from 3/2 to 7/4). However, observations - which rely on models to correct for projection effects - have suggested a flat distribution instead. The need for this correction is due to the lack of information regarding the line-of-sight distances. With a two decade long baseline in astrometric measurements, we are now able to measure significant projected radial accelerations, six of which are newly reported here, that directly constrain line-of-sight distances. Here we present a statistical approach to take advantage of this information and more accurately constrain the shape of the radial density profile of the late-type stellar population in the Galactic Center.


title: Perverse Sheaves and Knot Contact Homology

abstract: In this paper, which is mostly a research announcement, we give a new algebraic construction of knot contact homology in the sense of L. Ng [Ng05a]. For a link $L$ in $ {\mathbb R}^3 $, we define a differential graded (DG) $k$-category $ \tilde{\mathscr A} $ with finitely many objects, whose quasi-equivalence class is a topological invariant of $ L $. In the case when $L$ is a knot, the endomorphism algebra of a distinguished object of $ \tilde{\mathscr A} $ coincides with the fully noncommutative knot DGA as defined by Ekholm, Etnyre, Ng and Sullivan in [EENS13a]. The input of our construction is a natural action of the braid group $B_n$ on the category of perverse sheaves on a two-dimensional disk with singularities at $n$ marked points, studied by Gelfand, MacPherson and Vilonen in [GMV96]. As an application, we show that the category of finite-dimensional representations of the link $k$-category $ \tilde{A} = H_0(\tilde{\mathscr A}) $ defined as the $0$th homology of our DG category $ \tilde{\mathscr A} $ is equivalent to the category of perverse sheaves on $ {\mathbb R}^3 $ which are singular along the link $ L $. We also obtain several generalizations of the category $ \tilde{\mathscr A} $ by extending the Gelfand-MacPherson-Vilonen braid action.


title: On the Largest Integer that is not a Sum of Distinct Positive $n$th Powers

abstract: It is known that for an arbitrary positive integer \(n\) the sequence \(S(x^n)=(1^n, 2^n, \ldots)\) is complete, meaning that every sufficiently large integer is a sum of distinct \(n\)th powers of positive integers. We prove that every integer \(m\geq (b-1)2^{n-1}(r+\frac{2}{3}(b-1)(2^{2n}-1)+2(b-2))^n-2a+ab\), where \(a=n!2^{n^2}\), \(b=2^{n^3}a^{n-1}\), \(r=2^{n^2-n}a\), is a sum of distinct \(n\)th powers of positive integers.


title: Flexible Dynamics of Two Quorum Sensing Coupled Repressilators

abstract: Genetic oscillators play important roles in cell life regulation. The regulatory efficiency usually depends strongly on the emergence of stable collective dynamic modes, which requires designing the interactions between genetic networks. We investigate the dynamics of two identical synthetic genetic repressilators coupled by an additional plasmid which implements quorum sensing (QS) in each network thereby supporting global coupling. In a basic genetic ring oscillator network in which three genes inhibit each other in unidirectional manner, QS stimulates the transcriptional activity of chosen genes providing for competition between inhibitory and stimulatory activities localized in those genes. The "promoter strength", the Hill cooperativity coefficient of transcription repression, and the coupling strength, i.e., parameters controlling the basic rates of genetic reactions, were chosen for extensive bifurcation analysis. The results are presented as a map of dynamic regimes. We found that the remarkable multistability of the anti-phase limit cycle and stable homogeneous and inhomogeneous steady states exists over broad ranges of control parameters. We studied the anti-phase limit cycle stability and the evolution of irregular oscillatory regimes in the parameter areas where the anti-phase cycle loses stability. In these regions we observed developing complex oscillations, collective chaos, and multistability between regular limit cycles and complex oscillations over uncommonly large intervals of coupling strength. QS-coupling stimulates the appearance of intrachaotic periodic windows with spatially symmetric and asymmetric partial limit cycles which, in turn, change the type of chaos from a simple anti-phase character into chaos composed of pieces of the trajectories having alternating polarity. Abstract truncated...


title: Constraining the Variability and Binary Fraction of Galactic Center Young Stars

abstract: We present constraints on the variability and binarity of young stars in the central 10 arcseconds (~0.4 pc) of the Milky Way Galactic Center (GC) using Keck Adaptive Optics data over a 12 year baseline. Given our experiment's photometric uncertainties, at least 36% of our sample's known early-type stars are variable. We identified eclipsing binary systems by searching for periodic variability. In our sample of spectroscopically confirmed and likely early-type stars, we detected the two previously discovered GC eclipsing binary systems. We derived the likely binary fraction of main sequence, early-type stars at the GC via Monte Carlo simulations of eclipsing binary systems, and find that it is at least 32% with 90% confidence.


title: InfraNotes: Inconspicuous Handwritten Trajectory Tracking for Lecture Note Recording with Infrared Sensors

abstract: Lecture notes are important for students to review and understand the key points in the class. Unfortunately, the students often miss or lose part of the lecture notes. In this paper, we design and implement an infrared sensor based system, InfraNotes, to automatically record the notes on the board by sensing and analyzing hand gestures of the lecturer. Compared with existing techniques, our system does not require special accessories with lecturers such as sensor-facilitated pens, writing surfaces or the video-taping infrastructure. Instead, it only has an infrared-sensor module on the eraser holder of black/white board to capture handwritten trajectories. With a lightweight framework for handwritten trajectory processing, clear lecture notes can be generated automatically. We evaluate the quality of lecture notes by three standard character recognition techniques. The results indicate that InfraNotes is a promising solution to create clear and complete lectures to promote the education.


title: mPDF: Framework for Watermarking PDF Files using Image Watermarking Algorithms

abstract: The advancement in digital technologies have made it possible to produce perfect copies of digital content. In this environment, malicious users reproduce the digital content and share it without compensation to the content owner. Content owners are concerned about the potential loss of revenue and reputation from piracy, especially when the content is available over the Internet. Digital watermarking has emerged as a deterrent measure towards such malicious activities. Several methods have been proposed for copyright protection and fingerprinting of digital images. However, these methods are not applicable to text documents as these documents lack rich texture information which is abundantly available in digital images. In this paper, a framework (mPDF) is proposed which facilitates the usage of digital image watermarking algorithms on text documents. The proposed method divides a text document into texture and non-texture blocks using an energy-based approach. After classification, a watermark is embedded inside the texture blocks in a content adaptive manner. The proposed method is integrated with five known image watermarking methods and its performance is studied in terms of quality and robustness. Experiments are conducted on documents in 11 different languages. Experimental results clearly show that the proposed method facilitates the usage of image watermarking algorithms on text documents and is robust against attacks such as print & scan, print screen, and skew. Also, the proposed method overcomes the drawbacks of existing text watermarking methods such as manual inspection and language dependency.


title: Computational Investigation of Half-Heusler Compounds for Spintronics Applications

abstract: We present first-principles density functional calculations of the electronic structure, magnetism, and structural stability of 378 $\textit{XYZ}$ half-Heusler compounds (with $X=$ Cr, Mn, Fe, Co, Ni, Ru, Rh, $Y=$ Ti, V, Cr, Mn, Fe, Ni, $Z=$ Al, Ga, In, Si, Ge, Sn, P, As, Sb). We find that a "Slater-Pauling density of states" with a gap or pseudogap at three states per atom below the gap in at least one spin channel is a common feature in half-Heusler compounds. We find that the presence of such a gap at the Fermi energy in one or both spin channels contributes greatly to the stability of a half-Heusler compound. We calculate the formation energy of each compound and systematically investigate its stability against all other phases in the Open Quantum Materials Database (OQMD). We represent the thermodynamic phase stability of each compound as its distance from the convex hull of stable phases in the respective chemical space and show that the hull distance of a compound is a good measure of the likelihood of its experimental synthesis. We identify 26 18-electron semiconductors, 45 half-metals, and 34 near half-metals with negative formation energy, that follow the Slater-Pauling rule of three electrons per atom. Our calculations predict new thermodynamically stable semiconducting phases NiScAs, RhTiP, and RuVAs, which merit further experimental exploration. Further, two interesting zero-moment half-metals, CrMnAs and MnCrAs, are calculated to have negative formation energy. In addition, our calculations predict a number of new, hitherto unreported, semiconducting (e.g., CoVGe, FeVAs), half-metallic (e.g., RhVSb), near half-metallic (e.g., CoFeSb, CoVP) half-Heusler compounds to lie close to the respective convex hull of stable phases, and thus may be experimentally realized under suitable synthesis conditions, resulting in potential candidates for various spintronics applications.


title: Instagram Post Data Analysis

abstract: Because of the spread of the Internet, social platforms become big data pools. From there we can learn about the trends, culture and hot topics. This project focuses on analyzing the data from Instagram. It shows the relationship of Instagram filter data with location and number of likes to give users filter suggestion on achieving more likes based on their location. It also analyzes the popular hashtags in different locations to show visual culture differences between different cities.


title: Densities of 3-vertex graphs

abstract: Let d_i(G) be the density of the 3-vertex i-edge graph in a graph G, i.e., the probability that three random vertices induce a subgraph with i edges. Let S be the set of all quadruples (d_0,d_1,d_2,d_3) that are arbitrary close to 3-vertex graph densities in arbitrary large graphs. Huang, Linial, Naves, Peled and Sudakov have recently determined the projection of the set S to the (d_0,d_3) plane. We determine the projection of the set S to all the remaining planes.


title: Review: Nonstationary Spatial Modeling, with Emphasis on Process Convolution and Covariate-Driven Approaches

abstract: In many environmental applications involving spatially-referenced data, limitations on the number and locations of observations motivate the need for practical and efficient models for spatial interpolation, or kriging. A key component of models for continuously-indexed spatial data is the covariance function, which is traditionally assumed to belong to a parametric class of stationary models. While convenient, the assumption of stationarity is rarely realistic; as a result, there is a rich literature on alternative methodologies which capture and model the nonstationarity present in most environmental processes. This review document provides a rigorous and concise description of the existing literature on nonstationary methods, paying particular attention to process convolution (also called kernel smoothing or moving average) approaches. A summary is also provided of more recent methods which leverage covariate information and yield both interpretational and computational benefits. Note: the article is borrowed from Chapters 1 and 2 of the author's Ph.D. dissertation, joint with Catherine A. Calder.


title: A topological approach to the problem of emergence in complex systems

abstract: Emergent patterns in complex systems are related to many intriguing phenomena in modern science and philosophy. Several conceptions such as weak, strong and robust emergence have been proposed to emphasize different epistemological and ontological aspects of the problem. One of the most important concerns is whether emergence is an intrinsic property of the reality we observe, or it is rather a consequence of epistemological limitations. To elucidate this question, we propose a novel approximation through constructive topology, a framework that allow us to map the space of observed objects (ontology) with the knowledge subject conceptual apparatus (epistemology). Focusing in a particular type of emergent processes, namely those accessible through experiments and from which we have still no clue on the mechanistic processes yielding its formation, we analyse how a knowledge subject would build a conceptual explanatory framework. Working on these systems, we identify concept disjunction as a critical logical operation needed to identify the constraints of the system. Next, focusing on a three-bits synthetic system, we show how the number and scope of the constraints hinder the development of such scheme. Interestingly, we observe that our framework is unable to identify global constraints, clearly linking the epistemological limits of the framework with an ontological feature of the system. This allows us to propose a definition of emergence strength which we make compatible with the scientific method through the active intervention of the observer on the system. We think that this definition reconciles previous attempts to classify emergent processes, at least for the specific kind we discuss here. The paper finishes discussing the relevance of global constraints in biological systems, understood as a downward causal influence exerted by natural selection.


title: The asymptotic trace norm of random circulants and the graph energy

abstract: We compute the expected normalized trace norm (matrix/graph energy) of random symmetric band circulant matrices and graphs in the limit of large sizes, and obtain explicit bounds on the rate of convergence to the limit, and on the probabilities of large deviations. We also show that random symmetric band Toeplitz matrices have the same limit norm assuming that their band widths remain small relative to their sizes. We compare the limit norms across a range of related random matrix and graph ensembles.


title: Fe3O4 thin films: controlling and manipulating an elusive quantum material

abstract: Fe3O4 (magnetite) is one of the most elusive quantum materials and at the same time one of the most studied transition metal oxide materials for thin film applications. The theoretically expected half-metallic behavior generates high expectations that it can be used in spintronic devices. Yet, despite the tremendous amount of work devoted to preparing thin films, the enigmatic first order metal-insulator transition and the hall mark of magnetite known as the Verwey transition, is in thin films extremely broad and occurs at substantially lower temperatures as compared to that in high quality bulk single crystals. Here we have succeeded in finding and making a particular class of substrates that allows the growth of magnetite thin films with the Verwey transition as sharp as in the bulk. Moreover, we are now able to tune the transition temperature and, using tensile strain, increase it to substantially higher values than in the bulk.


title: Atom Interferometry in a Warm Vapor

abstract: We demonstrate matterwave interference in a warm vapor of rubidium atoms. Established approaches to light pulse atom interferometry rely on laser cooling to concentrate a large ensemble of atoms into a velocity class resonant with the atom optical light pulse. In our experiment, we show that clear interference signals may be obtained without laser cooling. This effect relies on the Doppler selectivity of the atom interferometer resonance. This interferometer may be configured to measure accelerations, and we demonstrate that multiple interferometers may be operated simultaneously by addressing multiple velocity classes.


title: Flagella bending affects macroscopic properties of bacterial suspensions

abstract: To survive in harsh conditions, motile bacteria swim in complex environment and respond to the surrounding flow. Here we develop a PDE model describing how the flagella bending affects macroscopic properties of bacterial suspensions. First, we show how the flagella bending contributes to the decrease of the effective viscosity observed in dilute suspension. Our results do not impose tumbling (random re-orientation) as it was done previously to explain the viscosity reduction. Second, we demonstrate a possibility of bacterium escape from the wall entrapment due to the self-induced buckling of flagella. Our results shed light on the role of flexible bacterial flagella in interactions of bacteria with shear flow and walls or obstacles.


title: A construction of certain weak colimits and an exactness property of the 2-category of categories

abstract: Given a 2-category $\mathcal{A}$, a $2$-functor $\mathcal{A} \overset {F} {\longrightarrow} \mathcal{C}at$ and a distinguished 1-subcategory $\Sigma \subset \mathcal{A}$ containing all the objects, a $\sigma$-cone for $F$ (with respect to $\Sigma$) is a lax cone such that the structural $2$-cells corresponding to the arrows of $\Sigma$ are invertible. The conical $\sigma$-limit is the universal (up to isomorphism) $\sigma$-cone. The notion of $\sigma$-limit generalises the well known notions of pseudo and lax limit. We consider the fundamental notion of $\sigma$-filtered} pair $(\mathcal{A}, \, \Sigma)$ which generalises the notion of 2-filtered 2-category. We give an explicit construction of $\sigma$-filtered $\sigma$-colimits of categories, construction which allows computations with these colimits. We then state and prove a basic exactness property of the 2-category of categories, namely, that $\sigma$-filtered $\sigma$-colimits commute with finite weighted pseudo (or bi) limits. An important corollary of this result is that a $\sigma$-filtered $\sigma$-colimit of exact category valued 2-functors is exact. This corollary is essential in the 2-dimensional theory of flat and pro-representable 2-functors, that we develop elsewhere.


title: Learning What and Where to Draw

abstract: Generative Adversarial Networks (GANs) have recently demonstrated the capability to synthesize compelling real-world images, such as room interiors, album covers, manga, faces, birds, and flowers. While existing models can synthesize images based on global constraints such as a class label or caption, they do not provide control over pose or object location. We propose a new model, the Generative Adversarial What-Where Network (GAWWN), that synthesizes images given instructions describing what content to draw in which location. We show high-quality 128 x 128 image synthesis on the Caltech-UCSD Birds dataset, conditioned on both informal text descriptions and also object location. Our system exposes control over both the bounding box around the bird and its constituent parts. By modeling the conditional distributions over part locations, our system also enables conditioning on arbitrary subsets of parts (e.g. only the beak and tail), yielding an efficient interface for picking part locations. We also show preliminary results on the more challenging domain of text- and location-controllable synthesis of images of human actions on the MPII Human Pose dataset.


title: Approximate Nearest Neighbor Search on High Dimensional Data --- Experiments, Analyses, and Improvement (v1.0)

abstract: Approximate Nearest neighbor search (ANNS) is fundamental and essential operation in applications from many domains, such as databases, machine learning, multimedia, and computer vision. Although many algorithms have been continuously proposed in the literature in the above domains each year, there is no comprehensive evaluation and analysis of their performances. In this paper, we conduct a comprehensive experimental evaluation of many state-of-the-art methods for approximate nearest neighbor search. Our study (1) is cross-disciplinary (i.e., including 16 algorithms in different domains, and from practitioners) and (2) has evaluated a diverse range of settings, including 20 datasets, several evaluation metrics, and different query workloads. The experimental results are carefully reported and analyzed to understand the performance results. Furthermore, we propose a new method that achieves both high query efficiency and high recall empirically on majority of the datasets under a wide range of settings.


title: Volatility Smile as Relativistic Effect

abstract: We give an explicit formula for the probability distribution based on a relativistic extension of Brownian motion. The distribution 1) is properly normalized and 2) obeys the tower law (semigroup property), so we can construct martingales and self-financing hedging strategies and price claims (options). This model is a 1-constant-parameter extension of the Black-Scholes-Merton model. The new parameter is the analog of the speed of light in Special Relativity. However, in the financial context there is no "speed limit" and the new parameter has the meaning of a characteristic diffusion speed at which relativistic effects become important and lead to a much softer asymptotic behavior, i.e., fat tails, giving rise to volatility smiles. We argue that a nonlocal stochastic description of such (Levy) processes is inadequate and discuss a local description from physics. The presentation is intended to be pedagogical.


title: Clustering of Local Group distances: publication bias or correlated measurements? IV. The Galactic Center

abstract: Aiming at deriving a statistically well-justified Galactic Center distance, $R_0$, and reducing any occurrence of publication bias, we compiled the most comprehensive and most complete database of Galactic Center distances available to date, containing 273 new or revised $R_0$ estimates published since records began in October 1918 until June 2016. We separate our $R_0$ compilation into direct and indirect distance measurements. The latter include a large body of estimates that rely on centroid determinations for a range of tracer populations as well as measurements based on kinematic observations of objects at the solar circle, combined with a mass and/or rotational model of the Milky Way. Careful assessment of the Galactic Center distances resulting from orbital modeling and statistical parallax measurements in the Galactic nucleus yields our final Galactic Center distance recommendation of $R_0 = 8.3 \pm 0.2 \mbox{ (statistical)} \pm 0.4 \mbox{ (systematic)}$ kpc. The centroid-based distances are in good agreement with this recommendation. Neither the direct measurements nor the post-1990 centroid-based distance determinations suggest that publication bias may be important. The kinematics-based distance estimates are affected by significantly larger uncertainties, but they can be used to constrain the Galaxy's rotation velocity at the solar Galactocentric distance, $\Theta_0$. Our results imply that the International Astronomical Union-recommended Galactic Center distance ($R_0^{\rm IAU} = 8.5$ kpc) needs a downward adjustment, while its $\Theta_0$ recommendation ($\Theta_0 = 220$ km s$^{-1}$) requires a substantial upward revision.


title: Interstellar-Medium Mapping in M82 Through Light Echoes Around Supernova 2014J

abstract: We present multiple-epoch measurements of the size and surface brightness of the light echoes from supernova (SN) 2014J in the nearby starburst galaxy M82. Hubble Space Telescope (HST) ACS/WFC images were taken ~277 and ~416 days after B-band maximum in the filters F475W, F606W, and F775W. Observations with HST WFC3/UVIS images at epochs ~216 and ~365 days (Crotts 2015) are included for a more complete analysis. The images reveal the temporal evolution of at least two major light-echo components. The first one exhibits a filled ring structure with position-angle-dependent intensity. This radially extended, diffuse echo indicates the presence of an inhomogeneous interstellar dust cloud ranging from ~100 pc to ~500 pc in the foreground of the SN. The second echo component appears as an unresolved luminous quarter-circle arc centered on the SN. The wavelength dependence of scattering measured in different dust components suggests that the dust producing the luminous arc favors smaller grain sizes, while that causing the diffuse light echo may have sizes similar to those of the Milky Way dust. Smaller grains can produce an optical depth consistent with that along the supernova-Earth line of sight measured by previous studies around maximum light. Therefore, it is possible that the dust slab, from which the luminous arc arises, is also responsible for most of the extinction towards SN 2014J. The optical depths determined from the Milky Way-like dust in the scattering matters are lower than that produced by the dust slab.


title: Anomalous Nernst Effect in Dirac Semimetal Cd3As2

abstract: Dirac and Weyl semimetals display a host of novel properties. In Cd$_3$As$_2$, the Dirac nodes lead to a protection mechanism that strongly suppresses backscattering in zero magnetic field, resulting in ultrahigh mobility ($\sim$ 10$^7$ cm$^2$ V$^{-1}$ s$^{-1}$). In applied magnetic field, an anomalous Nernst effect is predicted to arise from the Berry curvature associated with the Weyl nodes. We report observation of a large anomalous Nernst effect in Cd$_3$As$_2$. Both the anomalous Nernst signal and transport relaxation time $\tau_{tr}$ begin to increase rapidly at $\sim$ 50 K. This suggests a close relation between the protection mechanism and the anomalous Nernst effect. In a field, the quantum oscillations of bulk states display a beating effect, suggesting that the Dirac nodes split into Weyl states, allowing the Berry curvature to be observed as an anomalous Nernst effect.


title: A dimension reduction method with applications for coefficient inversion of diffusion equations

abstract: In this paper, we present a dimension reduction method to reduce the dimension of parameter space and state space and efficiently solve inverse problems. To this end, proper orthogonal decomposition (POD) and radial basis function (RBF) are combined to represent the solution of forward model with a form of variable separation. This POD-RBF method can be used to efficiently evaluate the model's output. A gradient regularization method is presented to solve the inverse problem with fast convergence. A generalized cross validation method is suggested to select the regularization parameter and differential step size for the gradient computation. Because the regularization method needs many model's evaluations. This is desirable for POD-RBF method. Thus, the POD-RBF method is integrated with the gradient regularization method to provide an efficient approach to solve inverse problems. We focus on the coefficient inversion of diffusion equations using the proposed approach. Based on different types of measurement data and different basis functions for coefficients, we present a few numerical examples for the coefficient inversion. The numerical results show that accurate reconstruction for the coefficient can be achieved efficiently.