title: Lower Bound and optimality for a nonlinearly damped Timoshenko system with thermoelasticity

abstract: In this paper, we consider a vibrating nonlinear Timoshenko system with thermoelasticity with second sound. We first investigate the strong stability of this system, then we devote our efforts to obtain the strong lower energy estimates using Alabau--Boussouira's energy comparison principle introduced in \cite{2} (see also \cite{alabau}). One of the main advantages of these results is that they allows us to prove the optimality of the asymptotic results (as $t\rightarrow \infty$) obtained in \cite{ali}. We also extend to our model the nice results achieved in \cite{alabau} for the case of nonlinearly damped Timoshenko system with thermoelasticity. The optimality of our results is also investigated through some explicit examples of the nonlinear damping term. The proof of our results relies on the approach in \cite{AB1, AB2}.


title: Hopf Galois structures on symmetric and alternating extensions

abstract: By using our previous results on induced Hopf Galois structures and a recent result by Koch, Kohl, Truman and Underwood on normality, we determine which types of Hopf Galois structures occur on Galois extensions with Galois group isomorphic to alternating or symmetric groups.


title: Ultrafast Exciton Dissociation and Long-Lived Charge Separation in a Photovoltaic Pentacene MoS2 van der Waals Heterojunction

abstract: Van der Waals heterojunctions between two-dimensional (2D) layered materials and nanomaterials of different dimensions present unique opportunities for gate-tunable optoelectronic devices. Mixed dimensional p-n heterojunction diodes, such as p-type pentacene (0D) and n-type monolayer MoS2 (2D), are especially interesting for photovoltaic applications where the absorption cross-section and charge transfer processes can be tailored by rational selection from the vast library of organic molecules and 2D materials. Here, we study the kinetics of excited carriers in pentacene-MoS2 p-n type-II heterojunctions by transient absorption spectroscopy. These measurements show that the dissociation of MoS2 excitons occurs by hole transfer to pentacene on the time scale of 6.7 ps. In addition, the charge-separated state lives for 5.1 ns, up to an order of magnitude longer than the recombination lifetimes from previously reported 2D material heterojunctions. By studying the fractional amplitudes of the MoS2 decay processes, the hole transfer yield from MoS2 to pentacene is found to be approximately 50 percent, with the remaining holes undergoing trapping due to surface defects. Overall, the ultrafast charge transfer and long-lived charge-separated state in pentacene-MoS2 p-n heterojunctions suggest significant promise for mixed-dimensional van der Waals heterostructures in photovoltaics, photodetectors, and related optoelectronic technologies.


title: Loginson: a transform and load system for very large scale log analysis in large IT infrastructures

abstract: Nowadays, most systems and applications produce log records that are useful for security and monitoring purposes such as debugging programming errors, checking system status, and detecting configuration problems or even attacks. To this end, a log repository becomes necessary whereby logs can be accessed and visualized in a timely manner. This paper presents Loginson, a high-performance log centralization system for large-scale log collection and processing in large IT infrastructures. Besides log collection, Loginson provides high-level analytics through a visual interface for the purpose of troubleshooting critical incidents. We note that Loginson outperforms all of the other log centralization solutions by taking full advantage of the vertical scalability, and therefore decreasing Capital Expenditure (CAPEX) and Operating Expense (OPEX) costs for deployment scenarios with a huge volume of log data.


title: Asymptotic entanglement in quantum walks from delocalized initial states

abstract: We study the entanglement between the internal (spin) and external (position) degrees of freedom of the one-dimensional discrete time quantum walk starting from local and delocalized initial states whose time evolution is driven by Hadamard and Fourier coins. We obtain the dependence of the asymptotic entanglement with the initial dispersion of the state and establish a way to connect the asymptotic entanglement between local and delocalized states. We find out that the delocalization of the state increases the number of initial spin states which achieves maximal entanglement from two states (local) to a continuous set of spin states (delocalized) given by a simple relation between the angles of the initial spin state. We also carry out numerical simulations of the average entanglement along the time to confront with our analytical results.


title: Universal dark halo scaling relation for the dwarf spheroidal satellites

abstract: Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the values of their surface densities, $\Sigma_{V_{\rm max}}$, are in good agreement with those for the observed dwarf spheroidal satellites even without employing any fitting procedures. This implies that this surface density would not be largely affected by any baryonic feedbacks and thus universal. Moreover, all subhalos on the small scales of dwarf satellites are expected to obey the relation $\Sigma_{V_{\rm max}}\propto V_{\rm max}$, irrespective of differences in their orbital evolutions, host halo properties, and observed redshifts. Therefore, we find that the universal scaling relation for dark halos on dwarf galaxy mass scales surely exists and provides us important clues to understanding fundamental properties of dark halos. We also investigate orbital and dynamical evolutions of subhalos to understand the origin of this universal dark halo relation and find that most of subhalos evolve generally along the $r_{\rm max}\propto V_{\rm max}$ sequence, even though these subhalos have undergone different histories of mass assembly and tidal stripping. This sequence, therefore, should be the key feature to understand the nature of the universality of $\Sigma_{V_{\rm max}}$.


title: Asteroseismology of the {\it Kepler} target KIC\,9204718

abstract: The high precision data obtained by the {\it Kepler} satellite allows us to detect hybrid type pulsator candidates more accurately than the data obtained by ground-based observations. In this study, we present preliminary results on the new analysis of the {\it Kepler} light curve and high resolution spectroscopic observations of pulsating Am star KIC\,9204718. Our tentative analysis therefore show that the star has hybrid pulsational characteristics.


title: Quantum diffusion of electromagnetic fields of ultrarelativistic spin-half particles

abstract: We compute electromagnetic fields created by a relativistic charged spin-half particle in empty space at distances comparable to the particle Compton wavelength. The particle is described as a wave packet evolving according to the Dirac equation. It produces the electromagnetic field that is essentially different from the Coulomb field due to the quantum diffusion effect.


title: Integral Field Spectroscopy of the Low-Mass Companion HD984B with the Gemini Planet Imager

abstract: We present new observations of the low-mass companion to HD 984 taken with the Gemini Planet Imager as a part of the Gemini Planet Imager Exoplanet Survey campaign. Images of HD 984 B were obtained in the J (1.12--1.3 micron) and H (1.50--1.80 micron) bands. Combined with archival epochs from 2012 and 2014, we fit the first orbit to the companion to find an 18 AU (70 year) orbit with a 68% confidence interval between 14 and 28 AU, an eccentricity of 0.18 with a 68% confidence interval between 0.05 and 0.47, and an inclination of 119 degrees with a 68% confidence interval between 114 degrees and 125 degrees. To address considerable spectral covariance in both spectra, we present a method of splitting the spectra into low and high frequencies to analyze the spectral structure at different spatial frequencies with the proper spectral noise correlation. Using the split spectra, we compare to known spectral types using field brown dwarf and low-mass star spectra and find a best fit match of a field gravity M6.5+/-1.5 spectral type with a corresponding temperature of 2730+120 K. Photometry of the companion yields a luminosity of log(L_bol/L_sun) = -2.88+/-0.07 dex, using DUSTY models. Mass estimates, again from DUSTY models, find an age-dependent mass of 34+/-1 to 95+/-4 M_Jup. These results are consistent with previous measurements of the object.


title: Spot dynamics in a reaction-diffusion model of plant root hair initiation

abstract: We study pattern formation in a 2-D reaction-diffusion (RD) sub-cellular model characterizing the effect of a spatial gradient of a plant hormone distribution on a family of G-proteins associated with root-hair (RH) initiation in the plant cell Arabidopsis thaliana. The activation of these G-proteins, known as the Rho of Plants (ROPs), by the plant hormone auxin, is known to promote certain protuberances on root hair cells, which are crucial for both anchorage and the uptake of nutrients from the soil. Our mathematical model for the activation of ROPs by the auxin gradient is an extension of the model of Payne and Grierson [PLoS ONE, 12(4), (2009)], and consists of a two-component Schnakenberg-type RD system with spatially heterogeneous coefficients on a 2-D domain. The nonlinear kinetics in this RD system model the nonlinear interactions between the active and inactive forms of ROPs. By using a singular perturbation analysis to study 2-D localized spatial patterns of active ROPs, it is shown that the spatial variations in the nonlinear reaction kinetics, due to the auxin gradient, lead to a slow spatial alignment of the localized regions of active ROPs along the longitudinal midline of the plant cell. Numerical bifurcation analysis, together with time-dependent numerical simulations of the RD system are used to illustrate both 2-D localized patterns in the model, and the spatial alignment of localized structures.


title: The nil Temperley--Lieb algebra of type affine C

abstract: We introduce a type affine $C$ analogue of the nil Temperley--Lieb algebra, in terms of generators and relations. We show that this algebra $T(n)$, which is a quotient of the positive part of a Kac--Moody algebra of type $D_{n+1}^{(2)}$, has an easily described faithful representation as an algebra of creation and annihilation operators on particle configurations, reminiscent of the open TASEP model in statistical physics. The centre of $T(n)$ consists of polynomials in a certain element $Q$, and $T(n)$ is a free module of finite rank over its centre. We show how to localize $T(n)$ by adjoining an inverse of $Q$, and prove that the resulting algebra is a full matrix ring over a ring of Laurent polynomials over a field. Although $T(n)$ has wild representation type, over an algebraically closed field we can classify all the finite dimensional indecomposable representations of $T(n)$ in which $Q$ acts invertibly.


title: Multi-Robot Active Information Gathering with Periodic Communication

abstract: A team of robots sharing a common goal can benefit from coordination of the activities of team members, helping the team to reach the goal more reliably or quickly. We address the problem of coordinating the actions of a team of robots with periodic communication capability executing an information gathering task. We cast the problem as a multi-agent optimal decision-making problem with an information theoretic objective function. We show that appropriate techniques for solving decentralized partially observable Markov decision processes (Dec-POMDPs) are applicable in such information gathering problems. We quantify the usefulness of coordinated information gathering through simulation studies, and demonstrate the feasibility of the method in a real-world target tracking domain.


title: Cellulyzer - Automated analysis and interactive visualization/simulation of select cellular processes

abstract: Here we report on a set of programs developed at the ZMBH Bio-Imaging Facility for tracking real-life images of cellular processes. These programs perform 1) automated tracking; 2) quantitative and comparative track analyses of different images in different groups; 3) different interactive visualization schemes; and 4) interactive realistic simulation of different cellular processes for validation and optimal problem-specific adjustment of image acquisition parameters (tradeoff between speed, resolution, and quality with feedback from the very final results). The collection of programs is primarily developed for the common bio-image analysis software ImageJ (as a single Java Plugin). Some programs are also available in other languages (C++ and Javascript) and may be run simply with a web-browser; even on a low-end Tablet or Smartphone. The programs are available at https://github.com/nurlicht/CellulyzerDemo

categories: physics.bio-ph cs.CV physics.data-an q-bio.SC


title: Giant and accessible conductivity of charged domain walls in lithium niobate

abstract: Ferroelectric domain walls are nm-sized interfaces between sections of different allowed values of the spontaneous polarization, the so-called domains. These walls - neutral or charged - can be created, displaced, deleted, and recreated again in ferroelectric materials. [1-3] Owing to the recent progress in the studies of ferroelectrics, they are expected to be functional active elements of the future nano-electronics.[1,4,5] Metallic-like conductivity of charged domain walls (CDWs) in insulating ferroelectrics, predicted in 1970s6 and detected recently,[7-10] is especially attractive for applications. This important effect is still in its infancy. The electric currents are small, the access to the conductivity is hampered by contact barriers, and stability is low because of sophisticated domain structures and/or proximity of the Curie point. Here, we report on giant and accessible persistent CDW conductivity in lithium niobate (LN) crystals (LiNbO3) - a vital material for photonics. [12] Our results are by far superior to the data known for other materials: Increase of LN conductivity by more than 14 orders of magnitude owing to CDWs, access to the effect via Ohmic and semi-Ohmic contacts, and its high stability for temperatures of up to 70 {\deg}C are demonstrated. It is made clear why this big effect was missed earlier [11,25] in LN. Our results demonstrate that strong conductivity of CDWs is available in simple ferroelectric materials, possessing only two allowed orientations of the spontaneous polarization, far from the Curie point. Also, CDW functionalities can be combined with linear and nonlinear optical phenomena. [14,15] Our findings allow new generations of adaptive-optical elements, of electrically controlled integrated-optical chips for quantum photonics, and of advanced LN-semiconductor hybrid optoelectronic devices.


title: Iron Intercalated Covalent-Organic Frameworks: First Crystalline Porous Thermoelectric Materials

abstract: Covalent-organic frameworks (COFs) are intriguing platforms for designing functional molecular materials. Here, we present a computational study based on van der Waals dispersion-corrected hybrid density functional theory calculations to analyze the material properties of boroxine-linked and triazine-linked intercalated-COFs. The effect of Fe atoms on the electronic band structures near the Fermi energy level of the intercalated-COFs have been investigated. The density of states (DOSs) computations have been performed to analyze the material properties of these kind of intercalated-COFs. We predict that COFs's electronic properties can be fine tuned by adding Fe atoms between two organic layers in their structures. The new COFs are predicted to be thermoelectric materials. These intercalated-COFs provide a new strategy to create thermoelectric materials within a rigid porous network in a highly controlled and predictable manner.


title: Drop-on-coilable-fibre systems exhibit negative stiffness events and transitions in coiling morphology

abstract: We investigate the mechanics of elastic fibres carrying liquid droplets. In such systems, buckling may localize inside the drop cavity if the fibre is thin enough. This so-called drop-on-coilable-fibre system exhibits a surprising liquid-like response under compression, and a solid-like response under tension. Here we analyze this unconventional behavior in further details and find theoretical, numerical and experimental evidences of negative stiffness events. We find that the first and main negative stiffness regime owes its existence to the transfer of capillary-stored energy into mechanical curvature energy. The following negative stiffness events are associated with changes in the coiling morphology of the fibre. Eventually coiling becomes tightly locked into an ordered phase where liquid and solid deformations coexist.


title: Polynomial Profits in Renewable Resources Management

abstract: A system of renewal equations on a graph provides a framework to describe the exploitation of a biological resource. In this context, we formulate an optimal control problem, prove the existence of an optimal control and ensure that the target cost function is polynomial in the control. In specific situations, further information about the form of this dependence is obtained. As a consequence, in some cases the optimal control is proved to be necessarily bang--bang, in other cases the computations necessary to find the optimal control are significantly reduced.


title: Dynamics of Abell 3266 - I. An Optical View of a Complex Merging Cluster

abstract: We present spectroscopy of 880 galaxies within a 2-degree field around the massive, merging cluster Abell 3266. This sample, which includes 704 new measurements, was combined with the existing redshifts measurements to generate a sample of over 1300 spectroscopic redshifts; the largest spectroscopic sample in the vicinity of A3266 to date. We define a cluster sub-sample of 790 redshifts which lie within a velocity range of 14,000 to 22,000 kms$^{-1}$ and within 1 degree of the cluster centre. A detailed structural analysis finds A3266 to have a complex dynamical structure containing six groups and filaments to the north of the cluster as well as a cluster core which can be decomposed into two components split along a northeast-southwest axis, consistent with previous X-ray observations. The mean redshift of the cluster core is found to be $0.0594 \pm 0.0005$ and the core velocity dispersion is given as $1462^{+99}_{-99}$ kms$^{-1}$. The overall velocity dispersion and redshift of the entire cluster and related structures are $1337^{+67}_{-67}$ kms$^{-1}$ and $0.0596 \pm 0.0002$, respectively, though the high velocity dispersion does not represent virialised motions but rather is due to relative motions of the cluster components. We posit A3266 is seen following a merger along the northeast southwest axis, however, the rich substructure in the rest of the cluster suggests that the dynamical history is more complex than just a simple merger with a range of continuous dynamical interactions taking place. It is thus likely that turbulence in A3266 is very high, even for a merging cluster.


title: Can a photon be separated from its wave?

abstract: An experiment is proposed in which the overall path taken by a photon is indicated by the timing of a twin herald photon, while a particular segment of that path is determined by interference. The experiment is to be carried out in two parts. In the first, coincident pairs of photons generated by type I spontaneous- parametric-down-conversion are diffracted by a large-width grating to increase the coherence length of their associated waves, and to produce two sources of coincident photons. Upon experimental confirmation that coincidences have been maintained, one of the sources is used to furnish timing heralds, and the other to send photons to an unequal path interferometer. A photon path through the interferometer via the short arm is indicated by count synchronization with the herald. The exit output port used and final detection location are determined by the phase in the long arm. If output port usage can thus be controlled by the phase in the photon free arm, the path of the photon as particle will have been controlled by interference with an accompanying photon-empty wave.


title: Bootstrapped Graph Diffusions: Exposing the Power of Nonlinearity

abstract: Graph-based semi-supervised learning (SSL) algorithms predict labels for all nodes based on provided labels of a small set of seed nodes. Classic methods capture the graph structure through some underlying diffusion process that propagates through the graph edges. Spectral diffusion, which includes personalized page rank and label propagation, propagates through random walks. Social diffusion propagates through shortest paths. A common ground to these diffusions is their {\em linearity}, which does not distinguish between contributions of few "strong" relations and many "weak" relations. Recently, non-linear methods such as node embeddings and graph convolutional networks (GCN) demonstrated a large gain in quality for SSL tasks. These methods introduce multiple components and greatly vary on how the graph structure, seed label information, and other features are used. We aim here to study the contribution of non-linearity, as an isolated ingredient, to the performance gain. To do so, we place classic linear graph diffusions in a self-training framework. Surprisingly, we observe that SSL using the resulting {\em bootstrapped diffusions} not only significantly improves over the respective non-bootstrapped baselines but also outperform state-of-the-art non-linear SSL methods. Moreover, since the self-training wrapper retains the scalability of the base method, we obtain both higher quality and better scalability.


title: Blow-up Continuity for Type-I, Mean-Convex Mean Curvature Flow

abstract: Under mean curvature flow, a closed, embedded hypersurface $M(t)$ becomes singular in finite time. For certain classes of mean-convex mean curvature flows, we show the continuity of the first singular time $T$ and the limit set "$M(T)$", with respect to initial data. We employ an Angenent-like neck-pinching argument to force singularities in nearby flows. However, since we cannot prescribe initial data, we combine Andrews' $\alpha$-non-collapsed condition and Colding and Minicozzi's uniqueness of tangent flows to place appropriately sized spheres in the region inside the hypersurface.


title: Linguistic Knowledge as Memory for Recurrent Neural Networks

abstract: Training recurrent neural networks to model long term dependencies is difficult. Hence, we propose to use external linguistic knowledge as an explicit signal to inform the model which memories it should utilize. Specifically, external knowledge is used to augment a sequence with typed edges between arbitrarily distant elements, and the resulting graph is decomposed into directed acyclic subgraphs. We introduce a model that encodes such graphs as explicit memory in recurrent neural networks, and use it to model coreference relations in text. We apply our model to several text comprehension tasks and achieve new state-of-the-art results on all considered benchmarks, including CNN, bAbi, and LAMBADA. On the bAbi QA tasks, our model solves 15 out of the 20 tasks with only 1000 training examples per task. Analysis of the learned representations further demonstrates the ability of our model to encode fine-grained entity information across a document.


title: Thermodynamic Stabilization of Precipitates through Interface Segregation: Chemical Effects

abstract: Precipitation hardening, which relies on a high density of intermetallic precipitates, is a commonly utilized technique for strengthening structural alloys. Structural alloys are commonly strengthened through a high density of small size intermetallic precipitates. At high temperatures, however, the precipitates coarsen to reduce the excess energy of the interface, resulting in a significant reduction in the strengthening provided by the precipitates. In certain ternary alloys, the secondary solute segregates to the interface and results in the formation of a high density of nanosize precipitates that provide enhanced strength and are resistant to coarsening. To understand the chemical effects involved, and to identify such systems, we develop a thermodynamic model using the framework of the regular nanocrystalline solution model. For various global compositions, temperatures and thermodynamic parameters, equilibrium configuration of Mg-Sn-Zn alloy is evaluated by minimizing the Gibbs free energy function with respect to the region-specific (bulk solid-solution, interface and precipitate) concentrations and sizes. The results show that Mg$_2$Sn precipitates can be stabilized to nanoscale sizes through Zn segregation to Mg/Mg$_2$Sn interface, and the precipitates can be stabilized against coarsening at high-temperatures by providing a larger Zn concentration in the system. Together with the inclusion of elastic strain energy effects and the input of computationally informed interface thermodynamic parameters in the future, the model is expected to provide a more realistic prediction of segregation and precipitate stabilization in ternary alloys of structural importance.


title: Online Convex Optimization with Unconstrained Domains and Losses

abstract: We propose an online convex optimization algorithm (RescaledExp) that achieves optimal regret in the unconstrained setting without prior knowledge of any bounds on the loss functions. We prove a lower bound showing an exponential separation between the regret of existing algorithms that require a known bound on the loss functions and any algorithm that does not require such knowledge. RescaledExp matches this lower bound asymptotically in the number of iterations. RescaledExp is naturally hyperparameter-free and we demonstrate empirically that it matches prior optimization algorithms that require hyperparameter optimization.


title: The G-centre and gradable derived equivalences

abstract: We propose a generalisation for the notion of the centre of an algebra in the setup of algebras graded by an arbitrary abelian group G. Our generalisation, which we call the G-centre, is designed to control the endomorphism category of the grading shift functors. We show that the G-centre is preserved by gradable derived equivalences given by tilting modules. We also discuss links with existing notions in superalgebra theory and apply our results to derived equivalences of superalgebras.


title: Exploiting Strong Convexity from Data with Primal-Dual First-Order Algorithms

abstract: We consider empirical risk minimization of linear predictors with convex loss functions. Such problems can be reformulated as convex-concave saddle point problems, and thus are well suitable for primal-dual first-order algorithms. However, primal-dual algorithms often require explicit strongly convex regularization in order to obtain fast linear convergence, and the required dual proximal mapping may not admit closed-form or efficient solution. In this paper, we develop both batch and randomized primal-dual algorithms that can exploit strong convexity from data adaptively and are capable of achieving linear convergence even without regularization. We also present dual-free variants of the adaptive primal-dual algorithms that do not require computing the dual proximal mapping, which are especially suitable for logistic regression.


title: On Sampling from Massive Graph Streams

abstract: We propose Graph Priority Sampling (GPS), a new paradigm for order-based reservoir sampling from massive streams of graph edges. GPS provides a general way to weight edge sampling according to auxiliary and/or size variables so as to accomplish various estimation goals of graph properties. In the context of subgraph counting, we show how edge sampling weights can be chosen so as to minimize the estimation variance of counts of specified sets of subgraphs. In distinction with many prior graph sampling schemes, GPS separates the functions of edge sampling and subgraph estimation. We propose two estimation frameworks: (1) Post-Stream estimation, to allow GPS to construct a reference sample of edges to support retrospective graph queries, and (2) In-Stream estimation, to allow GPS to obtain lower variance estimates by incrementally updating the subgraph count estimates during stream processing. Unbiasedness of subgraph estimators is established through a new Martingale formulation of graph stream order sampling, which shows that subgraph estimators, written as a product of constituent edge estimators are unbiased, even when computed at different points in the stream. The separation of estimation and sampling enables significant resource savings relative to previous work. We illustrate our framework with applications to triangle and wedge counting. We perform a large-scale experimental study on real-world graphs from various domains and types. GPS achieves high accuracy with less than 1% error for triangle and wedge counting, while storing a small fraction of the graph with average update times of a few microseconds per edge. Notably, for a large Twitter graph with more than 260M edges, GPS accurately estimates triangle counts with less than 1% error, while storing only 40K edges.


title: Horde of Bandits using Gaussian Markov Random Fields

abstract: The gang of bandits (GOB) model \cite{cesa2013gang} is a recent contextual bandits framework that shares information between a set of bandit problems, related by a known (possibly noisy) graph. This model is useful in problems like recommender systems where the large number of users makes it vital to transfer information between users. Despite its effectiveness, the existing GOB model can only be applied to small problems due to its quadratic time-dependence on the number of nodes. Existing solutions to combat the scalability issue require an often-unrealistic clustering assumption. By exploiting a connection to Gaussian Markov random fields (GMRFs), we show that the GOB model can be made to scale to much larger graphs without additional assumptions. In addition, we propose a Thompson sampling algorithm which uses the recent GMRF sampling-by-perturbation technique, allowing it to scale to even larger problems (leading to a "horde" of bandits). We give regret bounds and experimental results for GOB with Thompson sampling and epoch-greedy algorithms, indicating that these methods are as good as or significantly better than ignoring the graph or adopting a clustering-based approach. Finally, when an existing graph is not available, we propose a heuristic for learning it on the fly and show promising results.


title: Performance Scaling Law for Multi-Cell Multi-User Massive MIMO

abstract: This work provides a comprehensive scaling law based performance analysis for multi-cell multi-user massive multiple-input-multiple-output (MIMO) downlink systems. Imperfect channel state information (CSI), pilot contamination, and channel spatial correlation are all considered. First, a sum- rate lower bound is derived by exploiting the asymptotically deterministic property of the received signal power, while keeping the random nature of other components in the signal-to-interference-plus-noise-ratio (SINR) intact. Via a general scaling model on important network parameters, including the number of users, the channel training energy and the data transmission power, with respect to the number of base station antennas, the asymptotic scaling law of the effective SINR is obtained, which reveals quantitatively the tradeoff of the network parameters. More importantly, pilot contamination and pilot contamination elimination (PCE) are considered in the analytical framework. In addition, the applicability of the derived asymptotic scaling law in practical systems with large but finite antenna numbers are discussed. Finally, sufficient conditions on the parameter scalings for the SINR to be asymptotically deterministic in the sense of mean square convergence are provided, which covers existing results on such analysis as special cases and shows the effect of PCE explicitly.


title: Global optimization of Lipschitz functions

abstract: The goal of the paper is to design sequential strategies which lead to efficient optimization of an unknown function under the only assumption that it has a finite Lipschitz constant. We first identify sufficient conditions for the consistency of generic sequential algorithms and formulate the expected minimax rate for their performance. We introduce and analyze a first algorithm called LIPO which assumes the Lipschitz constant to be known. Consistency, minimax rates for LIPO are proved, as well as fast rates under an additional H\"older like condition. An adaptive version of LIPO is also introduced for the more realistic setup where the Lipschitz constant is unknown and has to be estimated along with the optimization. Similar theoretical guarantees are shown to hold for the adaptive LIPO algorithm and a numerical assessment is provided at the end of the paper to illustrate the potential of this strategy with respect to state-of-the-art methods over typical benchmark problems for global optimization.


title: Online Learning Without Prior Information

abstract: The vast majority of optimization and online learning algorithms today require some prior information about the data (often in the form of bounds on gradients or on the optimal parameter value). When this information is not available, these algorithms require laborious manual tuning of various hyperparameters, motivating the search for algorithms that can adapt to the data with no prior information. We describe a frontier of new lower bounds on the performance of such algorithms, reflecting a tradeoff between a term that depends on the optimal parameter value and a term that depends on the gradients' rate of growth. Further, we construct a family of algorithms whose performance matches any desired point on this frontier, which no previous algorithm reaches.


title: Probing the atmosphere of a sub-Jovian planet orbiting a cool dwarf

abstract: We derive the 0.01 $\mu$m binned transmission spectrum, between 0.74 and 1.0 $\mu$m, of WASP-80b from low resolution spectra obtained with the FORS2 instrument attached to ESO's Very Large Telescope. The combination of the fact that WASP-80 is an active star, together with instrumental and telluric factors, introduces correlated noise in the observed transit light curves, which we treat quantitatively using Gaussian Processes. Comparison of our results together with those from previous studies, to theoretically calculated models reveals an equilibrium temperature in agreement with the previously measured value of 825K, and a sub-solar metallicity, as well as an atmosphere depleted of molecular species with absorption bands in the IR ($\gg 5\sigma$). Our transmission spectrum alone shows evidence for additional absorption from the potassium core and wing, whereby its presence is detected from analysis of narrow 0.003 $\mu$m bin light curves ($\gg 5\sigma$). Further observations with visible and near-UV filters will be required to expand this spectrum and provide more in-depth knowledge of the atmosphere. These detections are only made possible through an instrument-dependent baseline model and a careful analysis of systematics in the data.


title: Transfinitely valued Euclidean domains have arbitrary indecomposable order type

abstract: We prove that for every indecomposable ordinal there exists a (transfinitely valued) Euclidean domain whose minimal Euclidean norm is of that order type. Conversely, any such norm must have indecomposable type, and so we completely characterize the norm complexity of Euclidean domains. Modifying this construction, we also find a finitely valued Euclidean domain with no multiplicative integer valued norm.


title: Dissipative Effects on Inertial-Range Statistics at High Reynolds numbers

abstract: Using the unique capabilities of the Variable Density Turbulence Tunnel at the Max Planck Institute for Dynamics and Self-Organization, G\"{o}ttingen, we report experimental result on classical grid turbulence that uncover fine, yet important details of the structure functions in the inertial range. This was made possible by measuring extremely long time series of up to $10^{10}$ samples of the turbulent fluctuating velocity, which corresponds to $\mathcal{O}\left(10^5\right)$ large eddy turnover times. These classical grid measurements were conducted in a well-controlled environment at a wide range of high Reynolds numbers from $R_\lambda=110$ up to $R_\lambda=1600$, using both traditional hot-wire probes as well as NSTAP probes developed at Princeton University. We found that deviations from ideal scaling are anchored to the small scales and that dissipation influences the inertial-range statistics at scales larger than the near-dissipation range.


title: Euler's factorial series and global relations

abstract: Using Pad\'e approximations to the series $E(z)=\sum_{k=0}^\infty k!(-z)^k$, we address arithmetic and analytical questions related to its values in both $p$-adic and Archimedean valuations.


title: Tunable zero and first sounds in ultracold Fermi gases with Rabi coupling

abstract: We consider a weakly-interacting fermionic gas of alkali-metal atoms characterized by two hyper- fine states which are Rabi coupled. By using a Hartree approximation for the repulsive interaction we determine the zero-temperature equation of state of this Fermi gas in D spatial dimensions (D = 1, 2, 3). Then, adopting the Landau-Vlasov equation and hydrodynamic equations, we investi- gate the speed of first sound and zero sound. We show that the two sounds, which occur respectively in collisional and collisionless regimes, crucially depend on the interplay between interaction strength and Rabi coupling. Finally, we discuss for some experimentally relevant cases the effect of a trapping harmonic potential on the density profiles of the fermionic system.


title: A Computational Model of a Single-Photon Avalanche Diode Sensor for Transient Imaging

abstract: Single-Photon Avalanche Diodes (SPAD) are affordable photodetectors, capable to collect extremely fast low-energy events, due to their single-photon sensibility. This makes them very suitable for time-of-flight-based range imaging systems, allowing to reduce costs and power requirements, without sacrifizing much temporal resolution. In this work we describe a computational model to simulate the behaviour of SPAD sensors, aiming to provide a realistic camera model for time-resolved light transport simulation, with applications on prototyping new reconstructions techniques based on SPAD time-of-flight data. Our model accounts for the major effects of the sensor on the incoming signal. We compare our model against real-world measurements, and apply it to a variety of scenarios, including complex multiply-scattered light transport.


title: Continuous and discrete one dimensional autonomous fractional ODEs

abstract: In this paper, we study 1D autonomous fractional ODEs $D_c^{\gamma}u=f(u), 0< \gamma <1$, where $u: [0,\infty)\mapsto\mathbb{R}$ is the unknown function and $D_c^{\gamma}$ is the generalized Caputo derivative introduced by Li and Liu ( arXiv:1612.05103). Based on the existence and uniqueness theorem and regularity results in previous work, we show the monotonicity of solutions to the autonomous fractional ODEs and several versions of comparison principles. We also perform a detailed discussion of the asymptotic behavior for $f(u)=Au^p$. In particular, based on an Osgood type blow-up criteria, we find relatively sharp bounds of the blow-up time in the case $A>0, p>1$. These bounds indicate that as the memory effect becomes stronger ($\gamma\to 0$), if the initial value is big, the blow-up time tends to zero while if the initial value is small, the blow-up time tends to infinity. In the case $A<0, p>1$, we show that the solution decays to zero more slowly compared with the usual derivative. Lastly, we show several comparison principles and Gr\"onwall inequalities for discretized equations, and perform some numerical simulations to confirm our analysis.


title: Effective identifiability criteria for tensors and polynomials

abstract: A tensor $T$, in a given tensor space, is said to be $h$-identifiable if it admits a unique decomposition as a sum of $h$ rank one tensors. A criterion for $h$-identifiability is called effective if it is satisfied in a dense, open subset of the set of rank $h$ tensors. In this paper we give effective $h$-identifiability criteria for a large class of tensors. We then improve these criteria for some symmetric tensors. For instance, this allows us to give a complete set of effective identifiability criteria for ternary quintic polynomial. Finally, we implement our identifiability algorithms in Macaulay2.


title: Constellation Queries over Big Data

abstract: A geometrical pattern is a set of points with all pairwise distances (or, more generally, relative distances) specified. Finding matches to such patterns has applications to spatial data in seismic, astronomical, and transportation contexts. For example, a particularly interesting geometric pattern in astronomy is the Einstein cross, which is an astronomical phenomenon in which a single quasar is observed as four distinct sky objects (due to gravitational lensing) when captured by earth telescopes. Finding such crosses, as well as other geometric patterns, is a challenging problem as the potential number of sets of elements that compose shapes is exponentially large in the size of the dataset and the pattern. In this paper, we denote geometric patterns as constellation queries and propose algorithms to find them in large data applications. Our methods combine quadtrees, matrix multiplication, and unindexed join processing to discover sets of points that match a geometric pattern within some additive factor on the pairwise distances. Our distributed experiments show that the choice of composition algorithm (matrix multiplication or nested loops) depends on the freedom introduced in the query geometry through the distance additive factor. Three clearly identified blocks of threshold values guide the choice of the best composition algorithm. Finally, solving the problem for relative distances requires a novel continuous-to-discrete transformation. To the best of our knowledge this paper is the first to investigate constellation queries at scale.


title: A Unifying Bayesian Optimization Framework for Radio Frequency Localization

abstract: We consider the problem of estimating an RF-device's location based on observations, such as received signal strength, from a set of transmitters with known locations. We survey the literature on this problem, showing that previous authors have considered implicitly or explicitly various metrics. We present a novel Bayesian framework that unifies these works and shows how to optimize the location estimation with respect to a given metric. We demonstrate how the framework can incorporate a general class of algorithms, including both model-based methods and data-driven algorithms such fingerprinting. This is illustrated by re-deriving the most popular algorithms within this framework. When used with a data-driven approach, our framework has cognitive self-improving properties in that it provably improves with increasing data compared to traditional methods. Furthermore, we propose using the error-CDF as a unified way of comparing algorithms based on two methods: (i) stochastic dominance, and (ii) an upper bound on error-CDFs. We prove that an algorithm that optimizes any distance based cost function is not stochastically dominated by any other algorithm. This suggests that in lieu of the search for a universally best localization algorithm, the community should focus on finding the best algorithm for a given well-defined objective.


title: Realizing the Aerial Robotic Worker for Inspection Operations

abstract: This report overviews a set of recent contributions in the field of path planning that were developed to enable the realization of the autonomous aerial robotic worker for inspection operations. The specific algorithmic contributions address several fundamental challenges of robotic inspection and exploration, and specifically those of optimal coverage planning given an a priori known model of the structure to be inspected, full coverage, optimized and fast inspection path planning, as well as efficient exploration of completely unknown environments and structures. All of the developed path planners support both holonomic and nonholonomic systems, and respect the on-board sensor model and constraints. An overview of the achieved results, followed by an integrating architecture in order to enable fully autonomous and highly-efficient infrastructure inspection in both known and unknown environments.


title: Don't Fear the Bit Flips: Optimized Coding Strategies for Binary Classification

abstract: After being trained, classifiers must often operate on data that has been corrupted by noise. In this paper, we consider the impact of such noise on the features of binary classifiers. Inspired by tools for classifier robustness, we introduce the same classification probability (SCP) to measure the resulting distortion on the classifier outputs. We introduce a low-complexity estimate of the SCP based on quantization and polynomial multiplication. We also study channel coding techniques based on replication error-correcting codes. In contrast to the traditional channel coding approach, where error-correction is meant to preserve the data and is agnostic to the application, our schemes specifically aim to maximize the SCP (equivalently minimizing the distortion of the classifier output) for the same redundancy overhead.


title: CMU DeepLens: Deep Learning For Automatic Image-based Galaxy-Galaxy Strong Lens Finding

abstract: Galaxy-scale strong gravitational lensing is not only a valuable probe of the dark matter distribution of massive galaxies, but can also provide valuable cosmological constraints, either by studying the population of strong lenses or by measuring time delays in lensed quasars. Due to the rarity of galaxy-scale strongly lensed systems, fast and reliable automated lens finding methods will be essential in the era of large surveys such as LSST, Euclid, and WFIRST. To tackle this challenge, we introduce CMU DeepLens, a new fully automated galaxy-galaxy lens finding method based on Deep Learning. This supervised machine learning approach does not require any tuning after the training step which only requires realistic image simulations of strongly lensed systems. We train and validate our model on a set of 20,000 LSST-like mock observations including a range of lensed systems of various sizes and signal-to-noise ratios (S/N). We find on our simulated data set that for a rejection rate of non-lenses of 99%, a completeness of 90% can be achieved for lenses with Einstein radii larger than 1.4" and S/N larger than 20 on individual $g$-band LSST exposures. Finally, we emphasize the importance of realistically complex simulations for training such machine learning methods by demonstrating that the performance of models of significantly different complexities cannot be distinguished on simpler simulations. We make our code publicly available at https://github.com/McWilliamsCenter/CMUDeepLens .


title: Limits of the Kucera-Gacs coding method

abstract: Every real is computable from a Martin-Loef random real. This well known result in algorithmic randomness was proved by Kucera and Gacs. In this survey article we discuss various approaches to the problem of coding an arbitrary real into a Martin-Loef random real,and also describe new results concerning optimal methods of coding. We start with a simple presentation of the original methods of Kucera and Gacs and then rigorously demonstrate their limitations in terms of the size of the redundancy in the codes that they produce. Armed with a deeper understanding of these methods, we then proceed to motivate and illustrate aspects of the new coding method that was recently introduced by Barmpalias and Lewis-Pye and which achieves optimal logarithmic redundancy, an exponential improvement over the original redundancy bounds.


title: Absence of chaos in Digital Memcomputing Machines with solutions

abstract: Digital memcomputing machines (DMMs) are non-linear dynamical systems designed so that their equilibrium points are solutions of the Boolean problem they solve. In a previous work [Chaos 27, 023107 (2017)] it was argued that when DMMs support solutions of the associated Boolean problem then strange attractors cannot coexist with such equilibria. In this work, we demonstrate such conjecture. In particular, we show that both topological transitivity and the strongest property of topological mixing are inconsistent with the point dissipative property of DMMs when equilibrium points are present. This is true for both the whole phase space and the global attractor. Absence of topological transitivity is enough to imply absence of chaotic behavior. In a similar vein, we prove that if DMMs do not have equilibrium points, the only attractors present are invariant tori/periodic orbits with periods that may possibly increase with system size (quasi-attractors).


title: Cost-Optimal Learning of Causal Graphs

abstract: We consider the problem of learning a causal graph over a set of variables with interventions. We study the cost-optimal causal graph learning problem: For a given skeleton (undirected version of the causal graph), design the set of interventions with minimum total cost, that can uniquely identify any causal graph with the given skeleton. We show that this problem is solvable in polynomial time. Later, we consider the case when the number of interventions is limited. For this case, we provide polynomial time algorithms when the skeleton is a tree or a clique tree. For a general chordal skeleton, we develop an efficient greedy algorithm, which can be improved when the causal graph skeleton is an interval graph.


title: System-Theoretic Performance Metrics for Low-Inertia Stability of Power Networks

abstract: As bulk synchronous generators in the power grid are replaced by distributed generation interfaced through power electronics, inertia is removed from the system, prompting concerns over grid stability. Different metrics are available for quantifying grid stability and performance; however, no theoretical results are available comparing and contrasting these metrics. This paper presents a rigorous system-theoretic study of performance metrics for low-inertia stability. For networks with uniform parameters, we derive explicit expressions for the eigenvalue damping ratios, and for the $\mathcal{H}_{2}$ and $\mathcal{H}_{\infty}$ norms of the linearized swing dynamics, from external power disturbances to different phase/frequency performance outputs.These expressions show the dependence of system performance on inertia constants, damping constants, and on the grid topology. Surprisingly, we find that the $\mathcal{H}_2$ and $\mathcal{H}_{\infty}$ norms can display contradictory behavior as functions of the system inertia, indicating that low-inertia performance depends strongly on the chosen performance metric.


title: Streaming Weak Submodularity: Interpreting Neural Networks on the Fly

abstract: In many machine learning applications, it is important to explain the predictions of a black-box classifier. For example, why does a deep neural network assign an image to a particular class? We cast interpretability of black-box classifiers as a combinatorial maximization problem and propose an efficient streaming algorithm to solve it subject to cardinality constraints. By extending ideas from Badanidiyuru et al. [2014], we provide a constant factor approximation guarantee for our algorithm in the case of random stream order and a weakly submodular objective function. This is the first such theoretical guarantee for this general class of functions, and we also show that no such algorithm exists for a worst case stream order. Our algorithm obtains similar explanations of Inception V3 predictions $10$ times faster than the state-of-the-art LIME framework of Ribeiro et al. [2016].


title: $\epsilon$-Subgradient Algorithms for Bilevel Convex Optimization

abstract: This paper introduces and studies the convergence properties of a new class of explicit $\epsilon$-subgradient methods for the task of minimizing a convex function over the set of minimizers of another convex minimization problem. The general algorithm specializes to some important cases, such as first-order methods applied to a varying objective function, which have computationally cheap iterations. We present numerical experimentation regarding certain applications where the theoretical framework encompasses efficient algorithmic techniques, enabling the use of the resulting methods to solve very large practical problems arising in tomographic image reconstruction.


title: A Curved Magnetic field in the ring-like shell of bubble N4

abstract: We report the detection of a curved magnetic field in the ring-like shell of the bubble N4, derived from near-infrared polarization of reddened diskless stars located behind this bubble. The magnetic field in the shell is curved and parallel to the ring-like shell, and its strength is estimated to be $\sim120\,\mu$G in the plane of the sky. The magnetic field strength in the shell is significantly enhanced compared to the local field strength. We calculate the mass-to-flux ratio for the submillimeter clumps in the shell and find that they are all magnetically subcritical. Our results demonstrate that the magnetic field strengthens as the interstellar medium is compressed into a shell, and suggest that the magnetic field has the potential to hinder star formation triggered by \ion{H}{2} region expansion.


title: Joint Multichannel Deconvolution and Blind Source Separation

abstract: Blind Source Separation (BSS) is a challenging matrix factorization problem that plays a central role in multichannel imaging science. In a large number of applications, such as astrophysics, current unmixing methods are limited since real-world mixtures are generally affected by extra instrumental effects like blurring. Therefore, BSS has to be solved jointly with a deconvolution problem, which requires tackling a new inverse problem: deconvolution BSS (DBSS). In this article, we introduce an innovative DBSS approach, called DecGMCA, based on sparse signal modeling and an efficient alternative projected least square algorithm. Numerical results demonstrate that the DecGMCA algorithm performs very well on simulations. It further highlights the importance of jointly solving BSS and deconvolution instead of considering these two problems independently. Furthermore, the performance of the proposed DecGMCA algorithm is demonstrated on simulated radio-interferometric data.


title: High-fidelity single-shot readout for a spin qubit via an enhanced latching mechanism

abstract: The readout of semiconductor spin qubits based on spin blockade is fast but suffers from a small charge signal. Previous work suggested large benefits from additional charge mapping processes, however uncertainties remain about the underlying mechanisms and achievable fidelity. In this work, we study the single-shot fidelity and limiting mechanisms for two variations of an enhanced latching readout. We achieve average single-shot readout fidelities > 99.3% and > 99.86% for the conventional and enhanced readout respectively, the latter being the highest to date for spin blockade. The signal amplitude is enhanced to a full one-electron signal while preserving the readout speed. Furthermore, layout constraints are relaxed because the charge sensor signal is no longer dependent on being aligned with the conventional (2, 0) - (1, 1) charge dipole. Silicon donor-quantum-dot qubits are used for this study, for which the dipole insensitivity substantially relaxes donor placement requirements. One of the readout variations also benefits from a parametric lifetime enhancement by replacing the spin-relaxation process with a charge-metastable one. This provides opportunities to further increase the fidelity. The relaxation mechanisms in the different regimes are investigated. This work demonstrates a readout that is fast, has one-electron signal and results in higher fidelity. It further predicts that going beyond 99.9% fidelity in a few microseconds of measurement time is within reach.


title: Three-dimensional lattice Boltzmann models for solid-liquid phase change

abstract: A three-dimensional (3 D) multiple-relaxation-time (MRT) and a 3 D single-relaxation-time (SRT) lattice Boltzmann (LB) models are proposed for the solid-liquid phase change. The enthalpy conservation equation can be recovered from the present models. The reasonable relationship of the relaxation times in the MRT model is discussed. Both One-dimensional (1 D) melting and solidification with analytical solutions are respectively calculated by the SRT and MRT models for validation. Compared with the SRT model, the MRT one is more accurate to capture the phase interface. The MRT model is also verified with other published two-dimensional (2 D) numerical results. The validations suggest that the present MRT approach is qualified to simulate the 3 D solid-liquid phase change process. Furthermore, the influences of Rayleigh number and Prandtl number on the 3 D melting are investigated.


title: Does a Nash theory of gravity make dark energy superfluous?

abstract: Recently Aadne and Gr{\o}n have argued that dark energy may follow naturally from a Nash theory of gravity. In this brief note I argue why this cannot be the case.


title: The strong convexity spectra of grids

abstract: Let $D$ be a connected oriented graph. A set $S \subseteq V(D)$ is convex in $D$ if, for every pair of vertices $x, y \in S$, the vertex set of every $xy$-geodesic, ($xy$ shortest directed path) and every $yx$-geodesic in $D$ is contained in $S$. The convexity number, ${\rm con}(D)$, of a non-trivial oriented graph, $D$, is the maximum cardinality of a proper convex set of $D$. The strong convexity spectrum of the graph $G$, $S_{SC} (G)$, is the set $\{{ \rm con}(D) \colon\ D {\rm \ is \ a \ strong \ orientation \ of \ } G \}$. In this paper we prove that the problem of determining the convexity number of an oriented graph is $\mathcal{NP}$-complete, even for bipartite oriented graphs of arbitrary large girth, extending previous known results for graphs. We also determine $S_{SC} (P_n \Box P_m)$, for every pair of integers $n,m \ge 2$.


title: De la structure de jauge des \'equations de Maxwell \`a l'\'etude des th\'eories au del\`a du Mod\`ele Standard gr\^ace \`a la mesure du couplage trilin\'eaire du champ de Higgs

abstract: The concept of gauge theories progressively emerged during the 20th century while the building of the Standard Model (SM) of particle physics. This work focusses on some aspects of gauge theories. First, a modern interpretation of Maxwell's equations in terms of differential forms motivates a geometric definition of gauge theories as a theory of connections on principle bundles. After presenting some historical facts about the notion of gauge in particle physics, the geometrical construction is used to study the electroweak theory and the SM. Though the SM has not been defaulted by any experiments yet, it still remains theoretically incomplete. The Higgs boson discovered in 2012 at LHC, is a particle about which we do not know lots yet. Some interesting clues about how one should extend the SM could be obtained through the observation of significant deviations between the measured value of some parameter of the Higgs boson, and the one predicted by the SM. The trilinear Higgs coupling is one of those interesting parameters. The last part of this work is a phenomenological study of the impact of variations of this coupling on the distributions one can typically get from the CMS particle detector, at LHC.


title: Exoplanet Characterization by Multi-Observatory Transit Photometry with TESS and CHEOPS

abstract: Space-based photometric surveys have discovered large numbers of planets transiting other stars, but these observe in a single band-pass and yield only the planet radius, orbital period, and transit duration. Information on the masses, compositions, and any atmospheres of these planets requires additional observations from the ground or space. The Transiting Exoplanet Survey Satellite (TESS) will yield thousands of planets around bright stars suitable for such follow-up. In the absence of spectroscopy or spectrophotometry from space, observations through the different pass-bands of multiple space telescopes provide some spectral information useful for identifying false positive signals, distinguishing between reflected light and thermal emission from hot Jupiters, and detecting Rayleigh scattering by planetary atmospheres. We calculated the expected difference in transit depths measured by the TESS and Characterizing Exoplanets Satellites (CHEOPS) missions, which will be more sensitive to redder and bluer optical wavelengths, respectively. The difference due to companion or background stars is small (<3% for main sequence companions) and likely to be negligible and undetectable. For only a few "hot" Jupiters, can combined photometry disambiguate between the reflected and thermal signals from planets. However, Rayleigh scattering by hazy atmospheres with particles sizes near 0.04 $\mu$m and at pressure altitudes above ~1 mbar can be detected for ~100 transiting planets, assuming every planet has such an atmosphere. Hazes with this characteristic particle size do not obscure observations at longer (near-infrared) wavelengths; CHEOPS follow-up of TESS-detected planets could thus identify candidates suitable for further study with the James Webb Space Telescope.


title: Associating vectors in $\CC^n$ with rank 2 projections in $\RR^{2n}$: with applications

abstract: We will see that vectors in $\CC^n$ have natural analogs as rank 2 projections in $\RR^{2n}$ and that this association transfers many vector properties into properties of rank two projections on $\RR^{2n}$. We believe that this association will answer many open problems in $\CC^n$ where the corresponding problem in $\RR^n$ has already been answered - and vice versa. As a application, we will see that phase retrieval (respectively, phase retrieval by projections) in $\CC^n$ transfers to a variation of phase retrieval by rank 2 projections (respectively, phase retrieval by projections) on $\RR^{2n}$. As a consequence, we will answer the open problem: Give the complex version of Edidin's Theorem \cite{E} which classifies when projections do phase retrieval in $\RR^n$. As another application we answer a longstanding open problem concerning fusion frames by showing that fusion frames in $\CC^n$ associate with fusion frames in $\RR^{2n}$ with twice the dimension. As another application, we will show that a family of mutually unbiased bases in $\CC^n$ has a natural analog as a family of mutually unbiased rank 2 projections in $\RR^{2n}$. The importance here is that there are very few real mutually unbiased bases but now there are unlimited numbers of real mutually unbiased rank 2 projections to be used in their place. As another application, we will give a variaton of Edidin's theorem which gives a surprising classification of norm retrieval. Finally, we will show that equiangular and biangular frames in $\CC^n$ have an analog as equiangular and biangular rank 2 projections in $\RR^{2n}$.


title: What Would You Do? Acting by Learning to Predict

abstract: We propose to learn tasks directly from visual demonstrations by learning to predict the outcome of human and robot actions on an environment. We enable a robot to physically perform a human demonstrated task without knowledge of the thought processes or actions of the human, only their visually observable state transitions. We evaluate our approach on two table-top, object manipulation tasks and demonstrate generalisation to previously unseen states. Our approach reduces the priors required to implement a robot task learning system compared with the existing approaches of Learning from Demonstration, Reinforcement Learning and Inverse Reinforcement Learning.


title: Crystallization and glass transition in supercooled binary Lennard-Jones liquids

abstract: The classic Kob-Andersen (KA) binary Lennard-Jones mixtures which are designed to prevent crystallization has been extensively studied in simulation of slow dynamics. Although crystallization can occur if a liquid system is cooled slowly, so far the KA model has not been crystallized. Here we report using molecular simulation the observation of crystal growth in the supercooled KA liquids. The onset of crystallization is observed occurring at temperature $T_c= 0.55$ which is higher than the glass transition temperature of $T_g=0.40$. We further examine the statistical distribution of single particle displacements in crystallization and close to glass transition. The displacement distribution for crystallization exhibits a power-law decay, whereas the distribution for glassy relaxation reflects a Gaussian center, terminated with an exponential tail (namely dynamic heterogeneity). Finally, we predict in order to crystalize KA liquids the cooling rate is approximately equal to $10^{-22}$, which is about 15 odder lower than the typical MD cooling rate.


title: Towards Generalization and Simplicity in Continuous Control

abstract: This work shows that policies with simple linear and RBF parameterizations can be trained to solve a variety of continuous control tasks, including the OpenAI gym benchmarks. The performance of these trained policies are competitive with state of the art results, obtained with more elaborate parameterizations such as fully connected neural networks. Furthermore, existing training and testing scenarios are shown to be very limited and prone to over-fitting, thus giving rise to only trajectory-centric policies. Training with a diverse initial state distribution is shown to produce more global policies with better generalization. This allows for interactive control scenarios where the system recovers from large on-line perturbations; as shown in the supplementary video.


title: A discovery of young stellar objects in older clusters of the Large Magellanic Cloud

abstract: Recent studies have shown that an extended main-sequence turn-off is a common feature among intermediate-age clusters (1--3 Gyr) in the Magellanic Clouds. Multiple-generation star formation and stellar rotation or interacting binaries have been proposed to explain the feature. However, it remains controversial in the field of stellar populations. Here we present the main results of an ongoing star formation among older star clusters in the Large Magellanic Cloud. Cross-matching the positions of star clusters and young stellar objects has yielded 15 matches with 7 located in the cluster center. We demonstrate that this is not by chance by estimating local number densities of young stellar objects for each star cluster. This method is not based on isochrone fitting, which leads to some uncertainties in age estimation and methods of background subtraction. We also find no direct correlation between atomic hydrogen and the clusters. This suggests that gas accretion for fueling the star formation must be happening in situ. These findings support for the multiple-generations scenario as a plausible explanation for the extended main-sequence turn-off.


title: Structural Data Recognition with Graph Model Boosting

abstract: This paper presents a novel method for structural data recognition using a large number of graph models. In general, prevalent methods for structural data recognition have two shortcomings: 1) Only a single model is used to capture structural variation. 2) Naive recognition methods are used, such as the nearest neighbor method. In this paper, we propose strengthening the recognition performance of these models as well as their ability to capture structural variation. The proposed method constructs a large number of graph models and trains decision trees using the models. This paper makes two main contributions. The first is a novel graph model that can quickly perform calculations, which allows us to construct several models in a feasible amount of time. The second contribution is a novel approach to structural data recognition: graph model boosting. Comprehensive structural variations can be captured with a large number of graph models constructed in a boosting framework, and a sophisticated classifier can be formed by aggregating the decision trees. Consequently, we can carry out structural data recognition with powerful recognition capability in the face of comprehensive structural variation. The experiments shows that the proposed method achieves impressive results and outperforms existing methods on datasets of IAM graph database repository.


title: A general formalism for logarithmic structures

abstract: We extend the formalism of "log spaces" of arXiv:1507.06752 to topoi equipped with a sheaf of monoids, and discuss Deligne--Faltings structures and root stacks in this context.


title: Software Defined Space-Air-Ground Integrated Vehicular Networks: Challenges and Solutions

abstract: This article proposes a software defined space-air-ground integrated network architecture for supporting diverse vehicular services in a seamless, efficient, and cost-effective manner. Firstly, the motivations and challenges for integration of space-air-ground networks are reviewed. Secondly, a software defined network architecture with a layered structure is presented. To protect the legacy services in satellite, aerial, and territorial segments, resources in each segment are sliced through network slicing to achieve service isolation. Then, available resources are put into a common and dynamic space-air-ground resource pool, which is managed by hierarchical controllers to accommodate vehicular services. Finally, a case study is carried out, followed by discussion on some open research topics.


title: Algebraic and topological properties of big mapping class groups

abstract: Let $S$ be an orientable, connected surface with infinitely-generated fundamental group. The main theorem states that if the genus of $S$ is finite and at least 4, then the isomorphism type of the pure mapping class group associated to $S$, denoted $\mathrm{PMap}(S)$, detects the homeomorphism type of $S$. As a corollary, every automorphism of $\mathrm{PMap}(S)$ is induced by a homeomorphism, which extends a theorem of Ivanov from the finite-type setting. In the process of proving these results, we show that $\mathrm{PMap}(S)$ is residually finite if and only if $S$ has finite genus, demonstrating that the algebraic structure of $\mathrm{PMap}(S)$ can distinguish finite- and infinite-genus surfaces. As an independent result, we also show that $\mathrm{Map}(S)$ fails to be residually finite for any infinite-type surface $S$. In addition, we give a topological generating set for $\mathrm{PMap}(S)$ equipped with the compact-open topology. In particular, if $S$ has at most one end accumulated by genus, then $\mathrm{PMap}(S)$ is topologically generated by Dehn twists, otherwise the Dehn twists along with handle shifts topologically generate.


title: Enhanced entanglement of two optical modes in optomechanical systems via an optical parametric amplifier

abstract: We investigate the effect of a degenerate optical parametric amplifier (OPA) placed inside an optomechanical cavity on the steady-state entanglement of two cavity modes, which jointly interact with a mechanical resonator. Two cavity modes are respectively driven at the red and blue sideband associated with the mechanical resonator, which generates entanglement between them in the limit of resolved sideband. The OPA gives rise to single-mode squeezing of the cavity fields, which results in significant improvement of the two-mode entanglement. It is found that an optimal nonlinear gain of the OPA exists, depending on the system temperatures, which yields the maximum entanglement. The improvement is particularly remarkable for the system at cryogenic temperatures.


title: Essential role of long non-coding RNAs in de novo chromatin modifications: The genomic address code hypothesis

abstract: The epigenome, i.e. the whole of chromatin modifications, is transferred from mother to daughter cells during cell differentiation. When de novo chromatin modifications (establishment or erasure of, respectively, new or pre-existing DNA methylations and/or histone modifications) are made in a daughter cell, however, it has a different epigenome than its mother cell. Although de novo chromatin modifications are an important event that comprises elementary processes of cell differentiation, its molecular mechanism remains poorly understood. We argue in this Letter that a key to solving this problem lies in understanding the role of long non-coding RNAs (lncRNAs)- a type of RNA that is becoming increasingly prominent in epigenetic studies. Many studies show that lncRNAs form ribonucleo-protein complexes in the nucleus and are involved in chromatin modifications. However, chromatin-modifying enzymes lack the information about genomic positions on which they act. It is known, on the other hand, that a single-stranded RNA in general can bind to a double-stranded DNA to form a triple helix. If each lncRNA forms a ribonucleo-protein complex with chromatin-modifying enzymes on one hand and, at the same time, a triple helix with a genomic region based on its specific nucleotide sequence on the other hand, it can induce de novo chromatin modifications at specific sites. Thus, the great variety of lncRNAs can be explained by the requirement for the diversity of "genomic address codes" specific to their cognate genomic regions where de novo chromatin modifications take place.


title: Rational Dyck Paths in the Non Relatively Prime Case

abstract: We study the relationship between rational slope Dyck paths and invariant subsets of $\mathbb Z,$ extending the work of the first two authors in the relatively prime case. We also find a bijection between $(dn,dm)$--Dyck paths and $d$-tuples of $(n,m)$-Dyck paths endowed with certain gluing data. These are the first steps towards understanding the relationship between rational slope Catalan combinatorics and the geometry of affine Springer fibers and knot invariants in the non relatively prime case.


title: Critical temperature and upper critical field of Li$_2$Pd$_{3-x}$Cu$_x$B (x=0.0, 0.1, 0.2)superconductors

abstract: We studied the effects of substitution of Pd by Cu on the upper critical field of the noncentrosymmetric superconductor Li$_2$Pd$_{3-x}$Cu$_x$B, with x=0.0, 0.1 and 0.2. The upper critical field as a function of temperature was determined by resistance measurements at different magnetic fields. We found that the superconducting transition temperature decreases as the Cu content increases. Moreover, the temperature dependence of the upper critical field is linear in the range of the temperature studied and, at low temperature, is enhanced compared with the prediction of the Werthamer-Helfan-Hohenberg theory. This indicates that the breaking of Cooper pairs by spin orbit scattering and Pauli paramagnetism is negligible, and that the upper critical field enhancement is mainly because the electron-phonon coupling and disorder.


title: The Spherical Bolometric Albedo of Planet Mercury

abstract: Published reflectance data covering several different wavelength intervals has been combined and analyzed in order to determine the spherical bolometric albedo of Mercury. The resulting value of 0.088 +/- 0.003 spans wavelengths from 0 to 4 {\mu}m which includes over 99% of the solar flux. This bolometric result is greater than the value determined between 0.43 and 1.01 {\mu}m by Domingue et al. (2011, Planet. Space Sci., 59, 1853-1872). The difference is due to higher reflectivity at wavelengths beyond 1.01 {\mu}m. The average effective blackbody temperature of Mercury corresponding to the newly determined albedo is 436.3 K. This temperature takes into account the eccentricity of the planet's orbit (M\'endez and Rivera-Valet\'in. 2017. ApJL, 837, L1).


title: Symmetric Assembly Puzzles are Hard, Beyond a Few Pieces

abstract: We study the complexity of symmetric assembly puzzles: given a collection of simple polygons, can we translate, rotate, and possibly flip them so that their interior-disjoint union is line symmetric? On the negative side, we show that the problem is strongly NP-complete even if the pieces are all polyominos. On the positive side, we show that the problem can be solved in polynomial time if the number of pieces is a fixed constant.


title: Intrinsic group behaviour: dependence of pedestrian dyad dynamics on principal social and personal features

abstract: Being determined by human social behaviour, pedestrian group dynamics depends on "intrinsic properties" of the group such as the purpose of the pedestrians, their personal relation, their gender, age, and body size. In this work we quantitatively study the dynamical properties of pedestrian dyads by analysing a large data set of automatically tracked pedestrian trajectories in an unconstrained "ecological" setting (a shopping mall), whose relational group properties have been coded by three different human observers. We observed that females walk slower and closer than males, that workers walk faster, at a larger distance and more abreast than leisure oriented people, and that inter group relation has a strong effect on group structure, with couples walking very close and abreast, colleagues walking at a larger distance, and friends walking more abreast than family members. Pedestrian height (obtained automatically through our tracking system) influences velocity and abreast distance, both growing functions of the average group height. Results regarding pedestrian age show as expected that elderly people walk slowly, while active age adults walk at the maximum velocity. Groups with children have a strong tendency to walk in a non abreast formation, with a large distance (despite a low abreast distance). A cross-analysis of the interplay between these intrinsic features, taking in account also the effect of extrinsic crowd density, confirms these major effects but reveals also a richer structure. An interesting and unexpected result, for example, is that the velocity of groups with children {\it increases} with density, at least in the low-medium density range found under normal conditions in shopping malls. Children also appear to behave differently according to the gender of the parent.


title: Improved Quark Coalescence for a Multi-Phase Transport Model

abstract: The string melting version of a multi-phase transport model is often applied to high-energy heavy-ion collisions since the dense matter thus formed is expected to be in parton degrees of freedom. In this work we improve its quark coalescence component, which describes the hadronization of the partonic matter to a hadronic matter. We removed the previous constraint that forced the numbers of mesons, baryons, and antibaryons in an event to be separately conserved through the quark coalescence process. A quark now could form either a meson or a baryon depending on the distance to its coalescence partner(s). We then compare results from the improved model with the experimental data on hadron $dN/dy$, $p_{_{\rm T}}$ spectra, and $v_2$ in heavy-ion collisions from $\sqrt{s_{_{\rm NN}}}=62.4$ GeV to $5.02$ TeV. We show that, besides being able to describe these observables for low-$p_{_{\rm T}}$ pions and kaons, the improved model also better describes the low-$p_{_{\rm T}}$ baryon observables in general, especially the baryon $p_{_{\rm T}}$ spectra and antibaryon-to-baryon ratios for multistrange baryons.


title: Polynomial Time Algorithms for Dual Volume Sampling

abstract: We study dual volume sampling, a method for selecting k columns from an n x m short and wide matrix (n <= k <= m) such that the probability of selection is proportional to the volume spanned by the rows of the induced submatrix. This method was proposed by Avron and Boutsidis (2013), who showed it to be a promising method for column subset selection and its multiple applications. However, its wider adoption has been hampered by the lack of polynomial time sampling algorithms. We remove this hindrance by developing an exact (randomized) polynomial time sampling algorithm as well as its derandomization. Thereafter, we study dual volume sampling via the theory of real stable polynomials and prove that its distribution satisfies the "Strong Rayleigh" property. This result has numerous consequences, including a provably fast-mixing Markov chain sampler that makes dual volume sampling much more attractive to practitioners. This sampler is closely related to classical algorithms for popular experimental design methods that are to date lacking theoretical analysis but are known to empirically work well.


title: Compositeness of hadron resonances in finite volume

abstract: We develop a theoretical framework to quantify the structure of unstable hadron resonances. With the help of the corresponding system in a finite volume, we define the compositeness of resonance states which can be interpreted as a probability. This framework is used to study the structure of the scalar mesons f_0(980) and a_0(980). In both mesons, the Kbar K component dominates about a half of the wave function. The method is also applied to the Lambda(1405) resonance. We argue that a single energy level in finite volume represents the two eigenstates in infinite volume. The KbarK N component of Lambda(1405), including contributions from both eigenstates, is found to be 58%, and the rest is composed of the pi Sigma and other channels.


title: Statistical Analysis of Supernova Remnants in the Large Magellanic Cloud

abstract: We construct the most complete sample of supernova remnants (SNRs) in any galaxy - the Large Magellanic Cloud (LMC) SNR sample. We study their various properties such as spectral index ($\alpha$), size and surface-brightness. We suggest an association between the spatial distribution, environment density of LMC SNRs and their tendency to be located around supergiant shells. We find evidence that the 16 known type Ia LMC SNRs are expanding in a lower density environment compared to the Core-Collapse (CC) type. The mean diameter of our entire population (74) is 41 pc, which is comparable to nearby galaxies. We didn't find any correlation between the type of SN explosion, ovality or age. The $N(<D)$ relationship of $a={0.96}$ implies that the randomised diameters are readily mimicking such an exponent. The rate of SNe occurring in the LMC is estimated to be $\sim$1 per 200 yr. The mean $\alpha$ of the entire LMC SNR population is $\alpha=-0.52$, which is typical of most SNRs. However, our estimates show a clear flattening of the synchrotron $\alpha$ as the remnants age. As predicted, our CC SNRs sample are significantly brighter radio emitters than the type Ia remnants. We also estimate the $\Sigma - D$ relation for the LMC to have a slope $\sim3.8$ which is comparable with other nearby galaxies. We also find the residency time of electrons in the galaxy ($4.0-14.3$ Myr), implying that SNRs should be the dominant mechanism for the production and acceleration of CRs.


title: Recent status of the understanding of neutrino-nucleus cross section

abstract: In this work we have presented current understanding of neutrino-nucleon/nucleus cross sections in the few GeV energy region relevant for a precise determination of neutrino oscillation parameters and CP violation in the leptonic sector. In this energy region various processes like quasielastic and inelastic production of single and multipion production, coherent pion production, kaon, eta, hyperon production, associated particle production as well as deep inelastic scattering processes contribute to the neutrino event rates.


title: Phase retrieval by hyperplanes

abstract: We show that a scalable frame does phase retrieval if and only if the hyperplanes of its orthogonal complements do phase retrieval. We then show this result fails in general by giving an example of a frame for $\mathbb R^3$ which does phase retrieval but its induced hyperplanes fail phase retrieval. Moreover, we show that such frames always exist in $\mathbb R^d$ for any dimension $d$. We also give an example of a frame in $\mathbb R^3$ which fails phase retrieval but its perps do phase retrieval. We will also see that a family of hyperplanes doing phase retrieval in $\mathbb R^d$ must contain at least $2d-2$ hyperplanes. Finally, we provide an example of six hyperplanes in $\mathbb R^4$ which do phase retrieval.


title: Performance Bounds for Graphical Record Linkage

abstract: Record linkage involves merging records in large, noisy databases to remove duplicate entities. It has become an important area because of its widespread occurrence in bibliometrics, public health, official statistics production, political science, and beyond. Traditional linkage methods directly linking records to one another are computationally infeasible as the number of records grows. As a result, it is increasingly common for researchers to treat record linkage as a clustering task, in which each latent entity is associated with one or more noisy database records. We critically assess performance bounds using the Kullback-Leibler (KL) divergence under a Bayesian record linkage framework, making connections to Kolchin partition models. We provide an upper bound using the KL divergence and a lower bound on the minimum probability of misclassifying a latent entity. We give insights for when our bounds hold using simulated data and provide practical user guidance.

categories: math.ST cs.IT math.IT stat.ME stat.ML stat.TH


title: A large deviation principle for empirical measures on Polish spaces: Application to singular Gibbs measures on manifolds

abstract: We prove a large deviation principle for a sequence of point processes defined by Gibbs probability measures on a Polish space. This is obtained as a consequence of a more general Laplace principle for the non-normalized Gibbs measures. We consider three main applications: Conditional Gibbs measures on compact spaces, Coulomb gases on compact Riemannian manifolds and the usual Gibbs measures in the Euclidean space. Finally, we study the generalization of Fekete points and prove a deterministic version of the Laplace principle known as $\Gamma$-convergence. The approach is partly inspired by the works of Dupuis and co-authors. It is remarkably natural and general compared to the usual strategies for singular Gibbs measures.


title: High-mass twin stars with a multi-polytrope EoS

abstract: We show that in the 3-polytropes model of Hebeler et al. \cite{Hebeler:2013nza} for the neutron star equation of state at supersaturation densities a third family of compact stars can be obtained which confirms the possibility of high-mass twin stars that have coincident masses $M_1=M_2\approx 2~M_\odot$ and significantly different radii $|R_1-R_2|>2-3 $ km. We show that the causality constraint puts severe limitations on the maximum mass of the third family sequence which can be relaxed when this scheme is extended to four polytropes thus mimicking a realistic high-density matter EoS.


title: Sparse Quadratic Logistic Regression in Sub-quadratic Time

abstract: We consider support recovery in the quadratic logistic regression setting - where the target depends on both p linear terms $x_i$ and up to $p^2$ quadratic terms $x_i x_j$. Quadratic terms enable prediction/modeling of higher-order effects between features and the target, but when incorporated naively may involve solving a very large regression problem. We consider the sparse case, where at most $s$ terms (linear or quadratic) are non-zero, and provide a new faster algorithm. It involves (a) identifying the weak support (i.e. all relevant variables) and (b) standard logistic regression optimization only on these chosen variables. The first step relies on a novel insight about correlation tests in the presence of non-linearity, and takes $O(pn)$ time for $n$ samples - giving potentially huge computational gains over the naive approach. Motivated by insights from the boolean case, we propose a non-linear correlation test for non-binary finite support case that involves hashing a variable and then correlating with the output variable. We also provide experimental results to demonstrate the effectiveness of our methods.


title: The Hilbert metric on Teichm\"uller space and Earthquake

abstract: Hamenst\"adt gave a parametrization of the Teichm\"uller space of punctured surfaces such that the image under this parametrization is the interior of a polytope. In this paper, we study the Hilbert metric on the Teichm\"uller space of punctured surfaces based on this parametrization. We prove that every earthquake ray is an almost geodesic under the Hilbert metric.


title: Enlargeability, foliations, and positive scalar curvature

abstract: We extend the deep and important results of Lichnerowicz, Connes, and Gromov-Lawson which relate geometry and characteristic numbers to the existence and non-existence of metrics of positive scalar curvature (PSC). In particular, we show: that a spin foliation with Hausdorff homotopy groupoid of an enlargeable manifold admits no PSC metric; that any metric of PSC on such a foliation is bounded by a multiple of the reciprocal of the foliation K-area of the ambient manifold; and that Connes' vanishing theorem for characteristic numbers of PSC foliations extends to a vanishing theorem for Haefliger cohomology classes.


title: New results on multi-agent system consensus: A graph signal processing perspective

abstract: This paper revisits the problem of multi-agent consensus from a graph signal processing perspective. By defining the graph filter from the consensus protocol, we establish the direct relation between average consensus of multi-agent systems and filtering of graph signals. This relation not only provides new insights of the average consensus, it also turns out to be a powerful tool to design effective consensus protocols for uncertain networks, which is difficult to deal with by existing time-domain methods. In this paper, we consider two cases, one is uncertain networks modeled by an estimated Laplacian matrix and a fixed eigenvalue bound, the other is connected graphs with unknown topology. The consensus protocols are designed for both cases based on the protocol filter. Several numerical examples are given to demonstrate the effectiveness of our methods.


title: Photoionization-driven Absorption Lines Variability in Balmer Absorption Line Quasar LBQS 1206+1052

abstract: In this paper we present an analysis of absorption line variability in mini-BAL quasar LBQS 1206+1052. The SDSS spectrum demonstrates that the absorption troughs can be divided into two components of blueshift velocities of $\sim$700 km s$^{-1}$ and $\sim$1400 km s$^{-1}$ relative to the quasar rest-frame. The former component shows rare Balmer absorption, which is an indicator of high density absorbing gas, thus the quasar is worth follow-up spectroscopic observations. Our follow-up optical and near-infrared spectra using MMT, YFOSC, TripleSpec and DBSP reveal that the strengths of the absorption lines vary for both of the two components, while the velocities do not change. We reproduce all of the spectral data by assuming that only the ionization state of the absorbing gas is variable and that all other physical properties are invariable. The variation of ionization is consistent with the variation of optical continuum from the V-band light-curve. Additionally, we can not interpret the data by assuming that the variability is due to a movement of the absorbing gas. Therefore, our analysis strongly indicates that the absorption line variability in LBQS 1206+1052 is photoionization-driven. As shown from photo-ionization simulations, the absorbing gas with blueshift velocity of $\sim$700 km s$^{-1}$ has a density in the range of $10^9$ to $10^{10}$ cm$^{-3}$ and a distance of $\sim$1 pc, and the gas with blueshift velocity of $\sim$1400 km s$^{-1}$ has a density of $10^3$ cm$^{-3}$ and a distance of $\sim$1 kpc.


title: Almost isometries between Teichm\"uller spaces

abstract: We prove that the Teichm\"uller space of surfaces with given boundary lengths equipped with the arc metric (resp. the Teichm\"uller metric) is almost isometric to the Teichm\"uller space of punctured surfaces equipped with the Thurston metric (resp. the Teichm\"uller metric).


title: HYDRA: HYbrid Design for Remote Attestation (Using a Formally Verified Microkernel)

abstract: Remote Attestation (RA) allows a trusted entity (verifier) to securely measure internal state of a remote untrusted hardware platform (prover). RA can be used to establish a static or dynamic root of trust in embedded and cyber-physical systems. It can also be used as a building block for other security services and primitives, such as software updates and patches, verifiable deletion and memory resetting. There are three major classes of RA designs: hardware-based, software-based, and hybrid, each with its own set of benefits and drawbacks. This paper presents the first hybrid RA design, called HYDRA, that builds upon formally verified software components that ensure memory isolation and protection, as well as enforce access control to memory and other resources. HYDRA obtains these properties by using the formally verified seL4 microkernel. (Until now, this was only attainable with purely hardware-based designs.) Using seL4 requires fewer hardware modifications to the underlying microprocessor. Building upon a formally verified software component increases confidence in security of the overall design of HYDRA and its implementation. We instantiate HYDRA on two commodity hardware platforms and assess the performance and overhead of performing RA on such platforms via experimentation; we show that HYDRA can attest 10MB of memory in less than 500msec when using a Speck-based message authentication code (MAC) to compute a cryptographic checksum over the memory to be attested.


title: Exact MAP Inference by Avoiding Fractional Vertices

abstract: Given a graphical model, one essential problem is MAP inference, that is, finding the most likely configuration of states according to the model. Although this problem is NP-hard, large instances can be solved in practice. A major open question is to explain why this is true. We give a natural condition under which we can provably perform MAP inference in polynomial time. We require that the number of fractional vertices in the LP relaxation exceeding the optimal solution is bounded by a polynomial in the problem size. This resolves an open question by Dimakis, Gohari, and Wainwright. In contrast, for general LP relaxations of integer programs, known techniques can only handle a constant number of fractional vertices whose value exceeds the optimal solution. We experimentally verify this condition and demonstrate how efficient various integer programming methods are at removing fractional solutions.


title: Leveraging Sparsity for Efficient Submodular Data Summarization

abstract: The facility location problem is widely used for summarizing large datasets and has additional applications in sensor placement, image retrieval, and clustering. One difficulty of this problem is that submodular optimization algorithms require the calculation of pairwise benefits for all items in the dataset. This is infeasible for large problems, so recent work proposed to only calculate nearest neighbor benefits. One limitation is that several strong assumptions were invoked to obtain provable approximation guarantees. In this paper we establish that these extra assumptions are not necessary---solving the sparsified problem will be almost optimal under the standard assumptions of the problem. We then analyze a different method of sparsification that is a better model for methods such as Locality Sensitive Hashing to accelerate the nearest neighbor computations and extend the use of the problem to a broader family of similarities. We validate our approach by demonstrating that it rapidly generates interpretable summaries.


title: Searches for the Anomalous FCNC Top-Higgs Couplings with Polarized Electron Beam at the LHeC

abstract: In this paper, we study the single top and Higgs associated production $\rm e^- p\rightarrow \nu_e \bar{t} \rightarrow \nu_e h \bar{q}(h\rightarrow b\bar{b})$ in the top-Higgs FCNC couplings at the LHeC with the electron beam energy of $E_{e}$ = 60 GeV and $E_{e}$ = 120 GeV, combination of a 7 TeV and 50 TeV proton beam. With the possibility of e-beam polarization ($p_{e}$ = 0, $\pm0.6$), we distinct the Cut-based method and the Multivariate Analysis (MVA) based method, and compare with the current experimental and theoretical limits. It is shown that the branching ratio Br $\rm(t\to uh)$ can be probed to 0.113 (0.093) $\%$, 0.071 (0.057) $\%$, 0.030 (0.022) $\%$ and 0.024 (0.019) $\%$ with the Cut-based (MVA-based) analysis at ($E_{p}$, $E_{e}$) = (7 TeV, 60 GeV), ($E_{p}$, $E_{e}$) = (7 TeV, 120 GeV), ($E_{p}$, $E_{e}$) = (50 TeV, 60 GeV) and ($E_{p}$, $E_{e}$) = (50 TeV, 120 GeV) beam energy and 1$\sigma$ level. With the possibility of e-beam polarization, the expected limits can be probed down to 0.090 (0.073) $\%$, 0.056 (0.045) $\%$, 0.024 (0.018) $\%$ and 0.019 (0.015) $\%$, respectively.