title: P-type $\delta$-doping with Diborane on Si(001) for STM Based Dopant Device Fabrication

abstract: Hydrogen resist lithography using the tip of a scanning tunneling microscope (STM) is employed for patterning p-type nanostructures in silicon. For this, the carrier density and mobility of boron $\delta$-layers, fabricated by gas-phase doping, are characterized with low-temperature transport experiments. Sheet resistivities as low as $300\thinspace\Omega$ are found. Adsorption, incorporation and surface diffusion of the dopants are investigated by STM imaging and result in an upper bound of 2\,nm for the lithographic resolution which is also corroborated by fabricating a 7.5\,nm wide p-type nanowire and measuring its electrical properties. Finally, to demonstrate the feasibility of bipolar dopant device fabrication with this technique, we prepared a 100\,nm wide pn junction and show that its electrical behavior is similar to that of an Esaki diode.


title: Validating the Critical Point of Spontaneous Parametric Down-Conversion for over 600 Scanning MEMS Micro Mirrors on Wafer-Level

abstract: Sensors and actuators based on resonant micro-electro-mechanical systems (MEMS), such as scanning micro mirrors, are well-established in automotive and consumer products. As the areas of application broaden, the requirements for the MEMS are increasing. Devices outside of the performance specifications have to be rejected which is costly due to the high processing times of MEMS technologies. In particular, nonlinear system behavior is often found to cause unexpected device failure or performance issues. Thus, accurate simulation or rather system models which account for nonlinear sensor dynamics can not only increase process yield, but more importantly, lead to a comprehensive understanding of the underlying physics and to improved MEMS design. We have studied [1] the possibility of a rather drastic device failure induced by nonlinearities on the example of a resonant scanning MEMS micro mirror. On the level of a few selected chips, we have carefully measured the complex nonlinear system behavior and modeled it by a nonlinear mode-coupling phenomenon known as spontaneous parametric down-conversion (SPDC). The most intriguing feature of SPDC is the sudden change from a rather linear to a nonlinear system behavior at the critical oscillation amplitude. However, the threshold only lies within the range of the mirror's operational amplitude, if certain frequency resonance conditions regarding the mechanical modes are met. Thus, the critical amplitude strongly depends on the frequency spectrum of the MEMS design which in turn is largely influenced by fabrication imperfections. We validate the dependence of the critical amplitude on the resonance condition by measuring it for over 600 micro mirrors on wafer-level. Our work does not only validate the theory of SPDC with measurements on such a large scale, but also demonstrates modeling strategies which are essential for MEMS product design.


title: Double descent in the condition number

abstract: In solving a system of $n$ linear equations in $d$ variables $Ax=b$, the condition number of the $n,d$ matrix $A$ measures how much errors in the data $b$ affect the solution $x$. Estimates of this type are important in many inverse problems. An example is machine learning where the key task is to estimate an underlying function from a set of measurements at random points in a high dimensional space and where low sensitivity to error in the data is a requirement for good predictive performance. Here we discuss the simple observation, which is known but surprisingly little quoted (see Theorem 4.2 in \cite{Brgisser:2013:CGN:2526261}): when the columns of $A$ are random vectors, the condition number of $A$ is highest if $d=n$, that is when the inverse of $A$ exists. An overdetermined system ($n>d$) as well as an underdetermined system ($n<d$), for which the pseudoinverse must be used instead of the inverse, typically have significantly better, that is lower, condition numbers. Thus the condition number of $A$ plotted as function of $d$ shows a double descent behavior with a peak at $d=n$.


title: idris-ct: A Library to do Category Theory in Idris

abstract: We introduce idris-ct, a Idris library providing verified type definitions of categorical concepts.idris-ct strives to be a bridge between academy and industry, catering both to category theorists who want to implement and try their ideas in a practical environment and to businesses and engineers who care about formalization with category theory: It is inspired by similar libraries developed for theorem proving but remains very practical, being aimed at software production in business. Nevertheless, the use of dependent types allows for a formally correct implementation of categorical concepts, so that guarantees can be made on software properties.


title: Investigating solvent effects on the magnetic properties of molybdate ions (MoO$_{4}^{2-}$) with relativistic embedding

abstract: We investigate the ability of mechanical and electronic density functional theory (DFT)-based embedding approaches to describe the solvent effects on nuclear magnetic resonance (NMR) shielding constants of the $^{95}$Mo nucleus in the molybdate ion in aqueous solution. From the description obtained from calculations with two- and four-component relativistic Hamiltonians, we find that for such systems spin-orbit coupling effects are clearly important for absolute shielding values, but for relative quantities a scalar relativistic treatment provides a sufficient estimation of the solvent effects. We find that the electronic contributions to the solvent effects are relatively modest yet decisive to provide a more accurate magnetic response of the system, when compared to reference supermolecular calculations. We analyze the errors in the embedding calculations by statistical methods as well as through a real-space representation of NMR shielding densities, which are shown to provide a clear picture of the physical processes at play.


title: Changes to the extreme and erratic behaviour of cryptocurrencies during COVID-19

abstract: This paper introduces new methods for analysing the extreme and erratic behaviour of time series to evaluate the impact of COVID-19 on cryptocurrency market dynamics. Across 51 cryptocurrencies, we examine extreme behaviour through a study of distribution extremities, and erratic behaviour through structural breaks. First, we analyse the structure of the market as a whole and observe a reduction in self-similarity as a result of COVID-19, particularly with respect to structural breaks in variance. Second, we compare and contrast these two behaviours, and identify individual anomalous cryptocurrencies. Tether (USDT) and TrueUSD (TUSD) are consistent outliers with respect to their returns, while Holo (HOT), NEXO (NEXO), Maker (MKR) and NEM (XEM) are frequently observed as anomalous with respect to both behaviours and time. Even among a market known as consistently volatile, this identifies individual cryptocurrencies that behave most irregularly in their extreme and erratic behaviour and shows these were more affected during the COVID-19 market crisis.


title: Data Exploration and Validation on dense knowledge graphs for biomedical research

abstract: Here we present a holistic approach for data exploration on dense knowledge graphs as a novel approach with a proof-of-concept in biomedical research. Knowledge graphs are increasingly becoming a vital factor in knowledge mining and discovery as they connect data using technologies from the semantic web. In this paper we extend a basic knowledge graph extracted from biomedical literature by context data like named entities and relations obtained by text mining and other linked data sources like ontologies and databases. We will present an overview about this novel network. The aim of this work was to extend this current knowledge with approaches from graph theory. This method will build the foundation for quality control, validation of hypothesis, detection of missing data and time series analysis of biomedical knowledge in general. In this context we tried to apply multiple-valued decision diagrams to these questions. In addition this knowledge representation of linked data can be used as FAIR approach to answer semantic questions. This paper sheds new lights on dense and very large knowledge graphs and the importance of a graph-theoretic understanding of these networks.


title: A mixture theory-based finite element formulation for the study of biodegradation of poroelastic scaffolds

abstract: We derive a mixture theory-based mathematical model of the degradation of a poroelastic solid immersed in a fluid bath. The evolution of the solid's mechanical and transport properties are also modeled. The inspiration for the model is the study of the temporal evolution of biodegradable Tissue Engineered Nerve Guides (TENGs), which are surgical implants supporting the alignment and re-growth of damaged nerves. The model comprises of the degrading solid, the degradation reaction products, and the fluid in which the solid is immersed. The weak formulation of the partial differential equations (PDEs) so derived is numerically implemented using a Finite Element Method (FEM). The numerical model is studied for stability and convergence rates using the Method of Manufactured Solutions.


title: X-ray reverberation lags from the 1.5 Seyfert galaxy NGC 5273

abstract: We present the results of X-ray spectral-timing analysis of a 90ks XMM-Newton observation of the nearby, broad line, early-type AGN NGC5273. The X-ray spectrum revealed the clear presence of a reflection component at high energies, with a clear signature of a narrow iron line at 6.4 keV, consistent with distant reflection. Applying a relativistic reflection model, we found only marginal evidence for a broader relativistic line component. However, cross-spectral analysis revealed that, between 4 and 6 $\times 10^{-4}$ Hz, the 5-8 keV band lagged the 2-3 keV band, implying reflection of the iron line from material close to the black hole. From the analysis of the lag-energy spectrum, we found a broad, but skewed line with a peak of $\approx$ 1000s at 7.5 keV relative to the continuum, which we interpret as the iron line in the reverberation spectrum from an illuminated accretion disc. From the asymmetry in the shape of lag-energy spectrum, we also found that the source is consistent with having an inclination $\geq 45^\circ$.


title: CRNs Exposed: Systematic Exploration of Chemical Reaction Networks

abstract: Formal methods have enabled breakthroughs in many fields, such as in hardware verification, machine learning and biological systems. The key object of interest in systems biology, synthetic biology, and molecular programming is chemical reaction networks (CRNs) which formalizes coupled chemical reactions in a well-mixed solution. CRNs are pivotal for our understanding of biological regulatory and metabolic networks, as well as for programming engineered molecular behavior. Although it is clear that small CRNs are capable of complex dynamics and computational behavior, it remains difficult to explore the space of CRNs in search for desired functionality. We use Alloy, a tool for expressing structural constraints and behavior in software systems, to enumerate CRNs with declaratively specified properties. We show how this framework can enumerate CRNs with a variety of structural constraints including biologically motivated catalytic networks and metabolic networks, and seesaw networks motivated by DNA nanotechnology. We also use the framework to explore analog function computation in rate-independent CRNs. By computing the desired output value with stoichiometry rather than with reaction rates (in the sense that $X \to Y+Y$ computes multiplication by $2$), such CRNs are completely robust to the choice of reaction rates or rate law. We find the smallest CRNs computing the max, minmax, abs and ReLU (rectified linear unit) functions in a natural subclass of rate-independent CRNs where rate-independence follows from structural network properties.


title: A Constant-Factor Approximation for Directed Latency in Quasi-Polynomial Time

abstract: We give the first constant-factor approximation for the Directed Latency problem in quasi-polynomial time. Here, the goal is to visit all nodes in an asymmetric metric with a single vehicle starting at a depot $r$ to minimize the average time a node waits to be visited by the vehicle. The approximation guarantee is an improvement over the polynomial-time $O(\log n)$-approximation [Friggstad, Salavatipour, Svitkina, 2013] and no better quasi-polynomial time approximation algorithm was known. To obtain this, we must extend a recent result showing the integrality gap of the Asymmetric TSP-Path LP relaxation is bounded by a constant [K\"{o}hne, Traub, and Vygen, 2019], which itself builds on the breakthrough result that the integrality gap for standard Asymmetric TSP is also a constant [Svensson, Tarnawsi, and Vegh, 2018]. We show the standard Asymmetric TSP-Path integrality gap is bounded by a constant even if the cut requirements of the LP relaxation are relaxed from $x(\delta^{in}(S)) \geq 1$ to $x(\delta^{in}(S)) \geq \rho$ for some constant $1/2 < \rho \leq 1$. We also give a better approximation guarantee in the special case of Directed Latency in regret metrics where the goal is to find a path $P$ minimize the average time a node $v$ waits in excess of $c_{rv}$, i.e. $\frac{1}{|V|} \cdot \sum_{v \in V} (c_v(P)-c_{rv})$.


title: Greenery Segmentation In Urban Images By Deep Learning

abstract: Vegetation is a relevant feature in the urban scenery and its awareness can be measured in an image by the Green View Index (GVI). Previous approaches to estimate the GVI were based upon heuristics image processing approaches and recently by deep learning networks (DLN). By leveraging some recent DLN architectures tuned to the image segmentation problem and exploiting a weighting strategy in the loss function (LF) we improved previously reported results in similar datasets.


title: On Metrics to Assess the Transferability of Machine Learning Models in Non-Intrusive Load Monitoring

abstract: To assess the performance of load disaggregation algorithms it is common practise to train a candidate algorithm on data from one or multiple households and subsequently apply cross-validation by evaluating the classification and energy estimation performance on unseen portions of the dataset derived from the same households. With an emerging discussion of transferability in Non-Intrusive Load Monitoring (NILM), there is a need for domain-specific metrics to assess the performance of NILM algorithms on new test scenarios being unseen buildings. In this paper, we discuss several metrics to assess the generalisation ability of NILM algorithms. These metrics target different aspects of performance evaluation in NILM and are meant to complement the traditional performance evaluation approach. We demonstrate how our metrics can be utilised to evaluate NILM algorithms by means of two case studies. We conduct our studies on several energy consumption datasets and take into consideration five state-of-the-art as well as four baseline NILM solutions. Finally, we formulate research challenges for future work.


title: The $l$-loop Banana Amplitude from GKZ Systems and relative Calabi-Yau Periods

abstract: We use the GKZ description of periods and certain classes of relative periods on families of Barth-Nieto Calabi-Yau $(l-1)$-folds in order to solve the $l$-loop banana amplitudes with their general mass dependence. As examples we compute the mass dependencies of the banana amplitudes up to the three-loop case and check the results against the known results for special mass values.


title: Optimal, Truthful, and Private Securities Lending

abstract: We consider a fundamental dynamic allocation problem motivated by the problem of $\textit{securities lending}$ in financial markets, the mechanism underlying the short selling of stocks. A lender would like to distribute a finite number of identical copies of some scarce resource to $n$ clients, each of whom has a private demand that is unknown to the lender. The lender would like to maximize the usage of the resource $\mbox{---}$ avoiding allocating more to a client than her true demand $\mbox{---}$ but is constrained to sell the resource at a pre-specified price per unit, and thus cannot use prices to incentivize truthful reporting. We first show that the Bayesian optimal algorithm for the one-shot problem $\mbox{---}$ which maximizes the resource's expected usage according to the posterior expectation of demand, given reports $\mbox{---}$ actually incentivizes truthful reporting as a dominant strategy. Because true demands in the securities lending problem are often sensitive information that the client would like to hide from competitors, we then consider the problem under the additional desideratum of (joint) differential privacy. We give an algorithm, based on simple dynamics for computing market equilibria, that is simultaneously private, approximately optimal, and approximately dominant-strategy truthful. Finally, we leverage this private algorithm to construct an approximately truthful, optimal mechanism for the extensive form multi-round auction where the lender does not have access to the true joint distributions between clients' requests and demands.


title: ManiGAN: Text-Guided Image Manipulation

abstract: The goal of our paper is to semantically edit parts of an image matching a given text that describes desired attributes (e.g., texture, colour, and background), while preserving other contents that are irrelevant to the text. To achieve this, we propose a novel generative adversarial network (ManiGAN), which contains two key components: text-image affine combination module (ACM) and detail correction module (DCM). The ACM selects image regions relevant to the given text and then correlates the regions with corresponding semantic words for effective manipulation. Meanwhile, it encodes original image features to help reconstruct text-irrelevant contents. The DCM rectifies mismatched attributes and completes missing contents of the synthetic image. Finally, we suggest a new metric for evaluating image manipulation results, in terms of both the generation of new attributes and the reconstruction of text-irrelevant contents. Extensive experiments on the CUB and COCO datasets demonstrate the superior performance of the proposed method. Code is available at https://github.com/mrlibw/ManiGAN.


title: On Ricci negative Lie groups

abstract: We give an overview of what is known on Lie groups admitting a left-invariant metric of negative Ricci curvature, including many natural questions and conjectures in the solvable case. We also introduce an open and convex cone C(n) of derivations attached to each nilpotent Lie algebra n, which is defined as the image of certain moment map and parametrizes a set of solvable Lie algebras with nilradical n admitting Ricci negative metrics.


title: Local theory for spatio-temporal canards and delayed bifurcations

abstract: We present a rigorous framework for the local analysis of canards and slow passages through bifurcations in a wide class of infinite-dimensional dynamical systems with time-scale separation. The framework is applicable to models where an infinite-dimensional dynamical system for the fast variables is coupled to a finite-dimensional dynamical system for slow variables. We prove the existence of centre-manifolds for generic models of this type, and study the reduced, finite-dimensional dynamics near bifurcations of (possibly) patterned steady states in the layer problem. Theoretical results are complemented with detailed examples and numerical simulations covering systems of local- and nonlocal-reaction diffusion equations, neural field models, and delay-differential equations. We provide analytical foundations for numerical observations recently reported in literature, such as spatio-temporal canards and slow-passages through Hopf bifurcations in spatially-extended systems subject to slow parameter variations. We also provide a theoretical analysis of slow passage through a Turing bifurcation in local and nonlocal models.


title: de Sitter symmetries and inflationary scalar-vector models

abstract: In this paper, we study the correspondence between a field theory in de Sitter space in $D$-dimensions and a dual conformal field theory in a euclidean space in $(D-1)$-dimensions. In particular, we investigate the form in which this correspondence is established for a system of interacting scalar and a vector fields propagating in de Sitter space. We analyze some necessary (but not sufficient) conditions for which conformal symmetry is preserved in the dual theory in $(D-1)$-dimensions, making possible the establishment of the correspondence. The discussion that we address on this paper is framed on the context of inflationary cosmology, so, the results obtained here pose some relevant possibilities of application to the calculation of the fields's correlation functions and of the primordial curvature perturbation $\zeta$, in inflationary models including coupled scalar and vector fields.


title: From deep learning to mechanistic understanding in neuroscience: the structure of retinal prediction

abstract: Recently, deep feedforward neural networks have achieved considerable success in modeling biological sensory processing, in terms of reproducing the input-output map of sensory neurons. However, such models raise profound questions about the very nature of explanation in neuroscience. Are we simply replacing one complex system (a biological circuit) with another (a deep network), without understanding either? Moreover, beyond neural representations, are the deep network's computational mechanisms for generating neural responses the same as those in the brain? Without a systematic approach to extracting and understanding computational mechanisms from deep neural network models, it can be difficult both to assess the degree of utility of deep learning approaches in neuroscience, and to extract experimentally testable hypotheses from deep networks. We develop such a systematic approach by combining dimensionality reduction and modern attribution methods for determining the relative importance of interneurons for specific visual computations. We apply this approach to deep network models of the retina, revealing a conceptual understanding of how the retina acts as a predictive feature extractor that signals deviations from expectations for diverse spatiotemporal stimuli. For each stimulus, our extracted computational mechanisms are consistent with prior scientific literature, and in one case yields a new mechanistic hypothesis. Thus overall, this work not only yields insights into the computational mechanisms underlying the striking predictive capabilities of the retina, but also places the framework of deep networks as neuroscientific models on firmer theoretical foundations, by providing a new roadmap to go beyond comparing neural representations to extracting and understand computational mechanisms.


title: Shaping representations through communication: community size effect in artificial learning systems

abstract: Motivated by theories of language and communication that explain why communities with large numbers of speakers have, on average, simpler languages with more regularity, we cast the representation learning problem in terms of learning to communicate. Our starting point sees the traditional autoencoder setup as a single encoder with a fixed decoder partner that must learn to communicate. Generalizing from there, we introduce community-based autoencoders in which multiple encoders and decoders collectively learn representations by being randomly paired up on successive training iterations. We find that increasing community sizes reduce idiosyncrasies in the learned codes, resulting in representations that better encode concept categories and correlate with human feature norms.


title: Theoretical Analysis of Metallic-Nanodimer Thermoplasmonics for Phototactic Nanoswimmers

abstract: We assess the potentiality of several geometries of metallic nanodimers (one of the simplest thermoplasmonic systems) as candidates for active particles (nanoswimmers) propelled and controlled by light (phototaxis). The studied nanodimers are formed by two spherical nanoparticles of gold, silver, or copper with radii ranging from 20 to 100 nm. Contrary to most proposals, which assume the asymmetry of the systems as a requirement for self-propulsion, our results show that nanodimers made of identical nanoparticles are excellent candidates for phototactic self-thermophoretic systems. Nonsymmetrical nanodimers, although having a tunable effective diffusion, possess much lower or zero average thermophoretic forces. We show that the effective diffusion and the net thermophoretic force in both types of systems depend strongly on the wavelength of the incident light, which makes these properties highly tunable. Our study may be useful for the design of simple-to-make but controllable self-propelled nanoparticles. This can find numerous applications ranging from autonomous drug-carrying to controlling the self-assembly of complex nanomaterials.


title: Hierarchical Bayesian Thermonuclear Rate for the $^7$Be(n,p)$^7$Li Big Bang Nucleosynthesis Reaction

abstract: Big bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than a thousand seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the Universe. This work derives new $^7$Be(n,p)$^7$Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of $^{7}$Be and $^7$Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel R-matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the $^7$Be(n,p)$^7$Li thermonuclear rates. From the resulting posteriors, we extract R-matrix parameters ($E_r$, $\gamma^2_n$, $\gamma^2_p$) and derive excitation energies, partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our $^7$Be(n,p)$^7$Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of $\leq$1 GK. We compare our rates to previous results and find that the $^7$Be(n,p)$^7$Li rates most commonly used in big bang simulations have too optimistic uncertainties.


title: Thermal transport in nanoporous holey silicon membranes investigated with optically-induced transient thermal gratings

abstract: In this study, we use the transient thermal grating optical technique \textemdash a non-contact, laser-based thermal metrology technique with intrinsically high accuracy \textemdash to investigate room-temperature phonon-mediated thermal transport in two nanoporous holey silicon membranes with limiting dimensions of 100 nm and 250 nm respectively. We compare the experimental results to ab initio calculations of phonon-mediated thermal transport according to the phonon Boltzmann transport equation (BTE) using two different computational techniques. We find that the calculations conducted within the Casimir framework, i.e. based on the BTE with the bulk phonon dispersion and diffuse scattering from surfaces, are in quantitative agreement with the experimental data, and thus conclude that this framework is adequate for describing phonon-mediated thermal transport through holey silicon membranes with feature sizes on the order of 100 nm.


title: The 1.28 GHz MeerKAT DEEP2 Image

abstract: We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.


title: Origins Space Telescope Mission Concept Study Report

abstract: The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.


title: Encoding Knowledge Graph Entity Aliases in Attentive Neural Network for Wikidata Entity Linking

abstract: The collaborative knowledge graphs such as Wikidata excessively rely on the crowd to author the information. Since the crowd is not bound to a standard protocol for assigning entity titles, the knowledge graph is populated by non-standard, noisy, long or even sometimes awkward titles. The issue of long, implicit, and nonstandard entity representations is a challenge in Entity Linking (EL) approaches for gaining high precision and recall. Underlying KG, in general, is the source of target entities for EL approaches, however, it often contains other relevant information, such as aliases of entities (e.g., Obama and Barack Hussein Obama are aliases for the entity Barack Obama). EL models usually ignore such readily available entity attributes. In this paper, we examine the role of knowledge graph context on an attentive neural network approach for entity linking on Wikidata. Our approach contributes by exploiting the sufficient context from a KG as a source of background knowledge, which is then fed into the neural network. This approach demonstrates merit to address challenges associated with entity titles (multi-word, long, implicit, case-sensitive). Our experimental study shows approx 8% improvements over the baseline approach, and significantly outperform an end to end approach for Wikidata entity linking.


title: Topological $\mathbb{Z}_2$ RVB quantum spin liquid on the ruby lattice

abstract: We construct a short-range resonating valence-bond state (RVB) on the ruby lattice, using projected entangled-pair states (PEPS) with bond dimension $D=3$. By introducing non-local moves to the dimer patterns on the torus, we distinguish four distinct sectors in the space of dimer coverings, which is a signature of the topological nature of the RVB wave function. Furthermore, by calculating the reduced density matrix of a bipartition of the RVB state on an infinite cylinder and exploring its entanglement entropy, we confirm the topological nature of the RVB wave function by obtaining non-zero topological contribution, $\gamma=-\rm{ln}\ 2$, consistent with that of a $\mathbb{Z}_2$ topological quantum spin liquid. We also calculate the ground-state energy of the spin-$\frac{1}{2}$ antiferromagnetic Heisenberg model on the ruby lattice and compare it with the RVB energy. Finally, we construct a quantum-dimer model for the ruby lattice and discuss it as a possible parent Hamiltonian for the RVB wave function.


title: Introduction to Galaxy Formation and Evolution. From Primordial Gas to Present-Day Galaxies

abstract: Present-day elliptical, spiral and irregular galaxies are large systems made of stars, gas and dark matter. Their properties result from a variety of physical processes that have occurred during the nearly fourteen billion years since the Big Bang. This comprehensive textbook, which bridges the gap between introductory and specialized texts, explains the key physical processes of galaxy formation, from the cosmological recombination of primordial gas to the evolution of the different galaxies that we observe in the Universe today. In a logical sequence, the book introduces cosmology, illustrates the properties of galaxies in the present-day Universe, then explains the physical processes behind galaxy formation in the cosmological context, taking into account the most recent developments in this field. The text ends on how to find distant galaxies with multi-wavelength observations, and how to extract the physical and evolutionary properties based on imaging and spectroscopic data.


title: Rounding Error Analysis of Mixed Precision Block Householder QR Algorithms

abstract: Although mixed precision arithmetic has recently garnered interest for training dense neural networks, many other applications could benefit from the speed-ups and lower storage cost if applied appropriately. The growing interest in employing mixed precision computations motivates the need for rounding error analysis that properly handles behavior from mixed precision arithmetic. We develop mixed precision variants of existing Householder QR algorithms and show error analyses supported by numerical experiments.


title: YOLACT++: Better Real-time Instance Segmentation

abstract: We present a simple, fully-convolutional model for real-time (>30 fps) instance segmentation that achieves competitive results on MS COCO evaluated on a single Titan Xp, which is significantly faster than any previous state-of-the-art approach. Moreover, we obtain this result after training on only one GPU. We accomplish this by breaking instance segmentation into two parallel subtasks: (1) generating a set of prototype masks and (2) predicting per-instance mask coefficients. Then we produce instance masks by linearly combining the prototypes with the mask coefficients. We find that because this process doesn't depend on repooling, this approach produces very high-quality masks and exhibits temporal stability for free. Furthermore, we analyze the emergent behavior of our prototypes and show they learn to localize instances on their own in a translation variant manner, despite being fully-convolutional. We also propose Fast NMS, a drop-in 12 ms faster replacement for standard NMS that only has a marginal performance penalty. Finally, by incorporating deformable convolutions into the backbone network, optimizing the prediction head with better anchor scales and aspect ratios, and adding a novel fast mask re-scoring branch, our YOLACT++ model can achieve 34.1 mAP on MS COCO at 33.5 fps, which is fairly close to the state-of-the-art approaches while still running at real-time.


title: The LUVOIR Mission Concept Study Final Report

abstract: The Large UV/Optical/Infrared Surveyor (LUVOIR) mission is one of four Decadal Survey Mission Concepts studied by NASA in preparation for the US National Academies' Astro2020 Decadal Survey. This observatory has the major goal of characterizing a wide range of exoplanets, including those that might be habitable -- or even inhabited. It would simultaneously enable a great leap forward in a broad range of astrophysics -- from the epoch of reionization, through galaxy formation and evolution, to star and planet formation. Powerful remote sensing observations of Solar System bodies will also be possible. This Final Report on the LUVOIR study presents the scientific motivations and goals of the mission concept, the engineering design, and technology development information. Please refer to the LUVOIR Final Report Appendices (separate document) for additional information.


title: A comparison of the real and non-archimedean Monge-Amp\`ere operator

abstract: Let $X$ be a proper algebraic variety over a non-archimedean, non-trivially valued field. We show that the non-archimedean Monge-Amp\`ere measure of a metric arising from a convex function on an open face of some skeleton of $X^{\text{an}}$ is equal to the real Monge-Amp\`ere measure of that function up to multiplication by a constant. As a consequence we obtain a regularity result for solutions of the non-archimedean Monge-Amp\`ere problem on curves.


title: An affine reconstructed algorithm for diffusion on triangular grids using the nodal discontinuous Galerkin method

abstract: This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection-diffusion equations. However, to the best of the authors' knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. The algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.


title: On the BTZ Black Hole and the Spinning Cosmic String

abstract: We reviewed the Ba\u nados-Teitelboim-Zanelli (BTZ) black hole solution in connection with the spinning string solution. We find a new exact solution, which can be related to the $(2+1)$-dimensional spinning point particle solution. There is no need for a cosmological constant, so the solution can be up-lifted to $(3+1)$ dimensions. The exact solution in a conformal invariant gravity model, where the spacetime is written as $g_{\mu\nu}=\omega^2 \tilde g_{\mu\nu}$, is horizon free and has an ergo-circle, while $\tilde g_{\mu\nu}$ is the BTZ solution. The dilaton $\omega$ determines the scale of the model. In accordance with the spinning cosmic string solution, it is conjectured that the new solution can be linked to the mass of the interior of the spinning cosmic string.


title: Arnold's potentials and quantum catastrophes

abstract: In the well known Thom's classification, every classical catastrophe is assigned a Lyapunov function. In the one-dimensional case, due to V. I. Arnold, these functions have polynomial form $V_{(k)}(x)= x^{k+1} + c_1x^{k-1} + \ldots$. A natural question is which features of the theory survive when such a function (say, with an even value of asymptotically dominant exponent $k+1$) is used as a confining potential in Schr\"{o}dinger equation. A few answers are formulated. Firstly, it is clarified that due to the tunneling, one of the possible classes of the measurable quantum catastrophes may be sought in a phenomenon of relocalization of the dominant part of the quantum particle density between different minima. For the sake of definiteness we just consider the spatially even Arnold's potentials and in the limit of the thick barriers and deep valleys we arrive at a systematic classification of the corresponding relocalization catastrophes.


title: Cable bacteria as long-range biological semiconductors

abstract: Filamentous cable bacteria exhibit unprecedented long-range biological electron transport, which takes place in a parallel fibre structure that shows an extraordinary electrical conductivity for a biological material. Still, the underlying electron transport mechanism remains undisclosed. Here we determine the intrinsic electrical properties of individual cable bacterium filaments. We retrieve an equivalent electrical circuit model, characterising cable bacteria as resistive biological wires. Temperature dependent experiments reveal that the charge transport is thermally activated, and can be described with an Arrhenius-type relation over a broad temperature range (-196{\deg}C to +50{\deg}C), thus excluding metal-like electron transport. Furthermore, when cable bacterium filaments are utilized as the channel in a field-effect transistor, they show n-type transport, indicating that electrons rather than holes are the charge carriers. Electron mobilities are in the order of 10$^{-1}$ cm$^2$/Vs, comparable to many organic semiconductors. This new type of biological centimetre-range semiconductor with low resistivity offers new perspectives for both fundamental studies and applications in (bio)electronics.


title: Forward-backward approximation of evolution equations in finite and infinite horizon

abstract: This research is concerned with evolution equations and their forward-backward discretizations. Our first contribution is an estimation for the distance between iterates of sequences generated by forward-backward schemes, useful in the convergence and robustness analysis of iterative algorithms of widespread use in variational analysis and optimization. Our second contribution is the approximation, on a bounded time frame, of the solutions of evolution equations governed by accretive (monotone) operators with an additive structure, by trajectories defined using forward-backward sequences. This provides a short, simple and self-contained proof of existence and regularity for such solutions; unifies and extends a number of classical results; and offers a guide for the development of numerical methods. Finally, our third contribution is a mathematical methodology that allows us to deduce the behavior, as the number of iterations tends to $+\infty$, of sequences generated by forward-backward algorithms, based solely on the knowledge of the behavior, as time goes to $+\infty$, of the solutions of differential inclusions, and viceversa.


title: Practical Quantum Computation of Chemical and Nuclear Energy Levels Using Quantum Imaginary Time Evolution and Lanczos Algorithms

abstract: Various methods have been developed for the quantum computation of the ground and excited states of physical and chemical systems, but many of them require either large numbers of ancilla qubits or high-dimensional optimization. The quantum imaginary-time evolution (QITE) and quantum Lanczos (QLanczos) methods proposed in [1] eschew the aforementioned issues. In this study, we demonstrate the practical application of these algorithms to nontrivial quantum computation, using the deuteron binding energy and molecular Hydrogen binding and excited state energies as examples. With the correct choice of initial and final states, we show that the number of time steps in QITE and QLanczos can be reduced significantly, which commensurately simplifies the required quantum circuit and improves compatibility with NISQ devices. We have performed these calculations on cloud-accessible IBM-Q quantum computers. With the application of readout-error mitigation and Richardson error extrapolation, we have obtained ground and excited state energies that agree well with exact results obtained from diagonalization.


title: Theme-Matters: Fashion Compatibility Learning via Theme Attention

abstract: Fashion compatibility learning is important to many fashion markets such as outfit composition and online fashion recommendation. Unlike previous work, we argue that fashion compatibility is not only a visual appearance compatible problem but also a theme-matters problem. An outfit, which consists of a set of fashion items (e.g., shirt, suit, shoes, etc.), is considered to be compatible for a "dating" event, yet maybe not for a "business" occasion. In this paper, we aim at solving the fashion compatibility problem given specific themes. To this end, we built the first real-world theme-aware fashion dataset comprising 14K around outfits labeled with 32 themes. In this dataset, there are more than 40K fashion items labeled with 152 fine-grained categories. We also propose an attention model learning fashion compatibility given a specific theme. It starts with a category-specific subspace learning, which projects compatible outfit items in certain categories to be close in the subspace. Thanks to strong connections between fashion themes and categories, we then build a theme-attention model over the category-specific embedding space. This model associates themes with the pairwise compatibility with attention, and thus compute the outfit-wise compatibility. To the best of our knowledge, this is the first attempt to estimate outfit compatibility conditional on a theme. We conduct extensive qualitative and quantitative experiments on our new dataset. Our method outperforms the state-of-the-art approaches.


title: A new renormalon in two dimensions

abstract: According to standard lore, perturbative series of super-renormalizable theories have only instanton singularities. In this paper we show that two-dimensional scalar theories with a spontaneously broken $O(N)$ symmetry at the classical level, which are super-renormalizable, have an IR renormalon singularity at large $N$. Since perturbative expansions in these theories are made around the "false vacuum" in which the global symmetry is broken, this singularity can be regarded as a manifestation of the non-perturbative absence of Goldstone bosons. We conjecture that the Borel singularity in the ground state energy of the Lieb--Liniger model is a non-relativistic manifestation of this phenomenon. We also provide {\it en passant} a detailed perturbative calculation of the Lieb--Liniger energy up to two-loops, and we check that it agrees with the prediction of the Bethe ansatz.

categories: hep-th cond-mat.quant-gas hep-lat hep-ph


title: Optimal Two-Sided Market Mechanism Design for Large-Scale Data Sharing and Trading in Massive IoT Networks

abstract: The development of the Internet of Things (IoT) generates a significant amount of data that contains valuable knowledge for system operations and business opportunities. Since the data is the property of the IoT data owners, the access to the data requires permission from the data owners, which gives rise to a potential market opportunity for the IoT data sharing and trading to create economic values and market opportunities for both data owners and buyers. In this work, we leverage optimal mechanism design theory to develop a monopolist matching platform for data trading over massive IoT networks. The proposed mechanism is composed of a pair of matching and payment rules for each side of the market. We analyze the incentive compatibility of the market and characterize the optimal mechanism with a class of cut-off matching rules for both welfare-maximization and revenue-maximization mechanisms and study three matching behaviors including complete-matched, bottom-eliminated, and top-reserved.


title: Descent for non-archimedean analytic spaces

abstract: In this paper we study two types of descent in the category of Berkovich analytic spaces: flat descent and descent with respect to an extension of the ground field. Quite surprisingly, the deepest results in this direction seem to be of the second type, including the descent of properties of being a good analytic space and being a morphism without boundary.


title: The role of Web of Science publications in China's tenure system

abstract: Tenure provides a permanent position to faculty in higher education institutions. In North America, it is granted to those who have established a record of excellence in research, teaching and services in a limited period. However, in China, research excellence represented by the number of Web of Science publications is highly weighted in the tenure assessment compared to excellence in teaching and services, but this has never been systematically investigated. By analyzing the tenure assessment documents from Chinese universities, this study reveals the role of Web of Science publications in China tenure system and presents the landscape of the tenure assessment process in Chinese higher education institutions.


title: Robust spin and charge excitations throughout high-$T_c$-cuprate phase diagram from incipient Mottness

abstract: The generic phase diagram of lightly hole-doped high-$T_c$-cuprates hosts antiferromagnetic insulating phase with well-defined spin-wave excitations. Contrary to the weak-coupling prediction, these modes persist up to the overdoped metallic regime as intense and dispersive paramagnons. Here we report on our study of the low-energy magnetic and charge excitations within the extended Hubbard model at strong-coupling, using a modified $1/N$ expansion method with a variational state serving as the saddle point solution. Despite clear separation of magnetic and Hubbard-$U$ energy scales, we find that incipient Mottness affects qualitatively dispersions and widths of magnetic modes throughout entire phase diagram. The obtained magnetic and charge dynamical structure factors agree semi-quantitatively with recent resonant $X$-ray and neutron scattering data for $\mathrm{La_{2-\mathit{x}}Sr_\mathit{x}CuO_4}$ and $\mathrm{(Bi,Pb)_2(Sr,La)_2CuO_{6+\delta}}$ at all available doping levels. The weak-coupling random-phase-approximation fails already for underdoped samples, pointing to the non-trivial intertwining of distinct energy scales in cuprate superconductors. The existence of a discrete charge mode which splits off the electron-hole continuum is also predicted.


title: Limited smoothness conditions with mixed norms for bilinear Fourier multipliers

abstract: In this paper, the $L^2 \times L^{\infty} \to L^2$ and $L^2 \times L^2 \to L^1$ boundedness of bilinear Fourier multiplier operators is discussed under weak smoothness conditions on multipliers. As an application, we prove the $L^2 \times BMO \to L^2$ and $L^2 \times L^2 \to H^1$ boundedness of bilinear operators with multipliers of limited smoothness satisfying vanishing conditions.


title: Atomic-waveguide quantum electrodynamics

abstract: Atom arrays are a new type of quantum light-matter interface. Here, we propose to employ one-dimensional ordered arrays as atomic waveguides. These arrays support optical guided modes that do not decay into free space. We show that these modes can be harnessed to mediate tunable, long-range interactions between additional "impurity qubits" coupled to the chain, without need for photonic structures. The efficient coupling between qubits and atomic waveguides enables the realization of tunable qubit-qubit interactions, which can be short- or long-range, dissipative or coherent, as well as chiral. Moreover, owing to the two-level nature of atoms, these waveguides are intrinsically quantum. In contrast to classical waveguides, where photons do not interact with each other, atomic waveguides display strong non-linearities, which create a tunable dissipative channel for qubit-qubit interactions, and opens the door to the exploration of many-body physics between guided photons. This physics is universal as it only relies on photon interference and can also be observed with other types of quantum emitters, such as those in molecular or solid-state systems.


title: A Curious Case of Circular Polarization in the 25 GHz Methanol Maser Line toward OMC-1

abstract: We report the detection of a circular polarization signature in the Stokes $V$ profile of a 25 GHz Class I CH$_3$OH maser toward the high mass star forming region OMC-1. Such a feature usually constitutes a detection of the Zeeman effect. If due to a magnetic field in OMC-1, this would represent the first detection and discovery of the Zeeman effect in the 25 GHz Class I CH$_3$OH maser. The feature in Stokes $V$ is detected in two observations with different angular resolutions taken eight years apart with the Very Large Array (VLA); for our 2009 D-configuration observations, the fitted value for $z B_\text{los}$ is $152 \pm 12$ Hz, where $z$ is the Zeeman splitting factor and $B_\text{los}$ is the line-of-sight magnetic field. For our 2017 C-configuration observations, the fitted value for $zB_\text{los} = 149 \pm 19~$Hz, likely for the same maser spot. These correspond to $B_\text{los}$ in the range 171-214 mG, depending on which hyperfine transition is responsible for the maser line. While these $B_\text{los}$ values are high, they are not implausible. If the magnetic field increases in proportion to the molecular hydrogen density in shocked regions, then our detected fields predict values for the pre-shock magnetic field that are in agreement with observations. With $B_\text{los}$ =171-214 mG, the magnetic energy in the post-shocked regions where these 25 GHz Class I CH$_3$OH masers occur would dominate over the kinetic energy density and be at least of the order of the pressure in the shock, implying that the magnetic field would exert significant influence over the dynamics of these regions.


title: Automatic Financial Feature Construction

abstract: In automatic financial feature construction task, the state-of-the-art technic leverages reverse polish expression to represent the features, then use genetic programming (GP) to conduct its evolution process. In this paper, we propose a new framework based on neural network, alpha discovery neural network (ADNN). In this work, we made several contributions. Firstly, in this task, we make full use of neural network overwhelming advantage in feature extraction to construct highly informative features. Secondly, we use domain knowledge to design the object function, batch size, and sampling rules. Thirdly, we use pre-training to replace the GP evolution process. According to neural network universal approximation theorem, pre-training can conduct a more effective and explainable evolution process. Experiment shows that ADNN can remarkably produce more diversified and higher informative features than GP. Besides, ADNN can serve as a data augmentation algorithm. It further improves the the performance of financial features constructed by GP.


title: Precise $\alpha_s$ determination from charmonium sum rules

abstract: The strong coupling, $\alpha_s$, governs perturbative Quantum Chromodynamics (QCD) and is one of the free parameters of the Standard Model. We introduce a new method that allows a precise extraction of $\alpha_s(m_Z)$ from dimensionless ratios of roots of moments of the charm-quark vector correlator. The ratios we use in our analysis have a rather weak logarithmic quark-mass dependence, starting at $\mathcal{O}(\alpha_s^2)$, and can be obtained from experimental data with good precision, since they benefit from positive correlations among the individual experimentally determined moments. We perform a careful and conservative error analysis with special emphasis on uncertainties related to the truncation of perturbation theory, treating the renormalization scales such as to ensure order-by-order convergence. Our final result, with expressions at $\mathcal{O}(\alpha_s^3)$, is $\alpha_s(m_Z)=0.1168\pm 0.0019$.


title: Estimates and Asymptotics for the stress concentration between closely spaced stiff $C^{1, \gamma}$ inclusions in linear elasticity

abstract: This paper is concerned with the stress concentration phenomenon in elastic composite materials when the inclusions are very closely spaced. We investigate the gradient blow-up estimates for the Lam\'{e} system of linear elasticity with partially infinite coefficients to show the dependence of the degree of stress enhancement on the distance between inclusions in a composite containing densely placed stiff inclusions. In this paper we assume that the inclusions to be of $C^{1, \gamma}$, weaker than the previous $C^{2, \gamma}$ assumption. To overcome this new difficulty, we make use of $W^{1, p}$ estimates for elliptic system with right hand side in divergence form, instead of a direct $W^{2, p}$ argument for $C^{2, \gamma}$ inclusion case, and combine with the Campanato's approach to establish the optimal gradient estimates, including upper and lower bounds. Moreover, an asymptotic formula of the gradient near the blow-up point is obtained for some symmetric $C^{1, \gamma}$ inclusions.


title: Dimensionality-enhanced quantum state transfer in long-range interacting spin systems

abstract: In this work we study the single-qubit quantum state transfer in uniform long-range spin XXZ systems in high-dimensional geometries. We consider prototypical long-range spin exchanges that are relevant for experiments in cold atomic platforms: Coulomb, dipolar and van der Waals-like interactions. We find that in all these cases the fidelity increases with the dimensionality of the lattice. This can be related to the emergence of a pair of bilocalized states on the sender and receiver site due to the onset of an effective weak-coupling Hamiltonian. The enhancement of the quantum state transfer fidelity is more pronounced both with the increase of the couplings interaction range and in going from a 1D to a 2D lattice. Finally, we test our predictions in the presence of temperature-induced disorder introducing a model for the thermal displacement of the lattice sites, considered as a set of local adiabatic oscillators.


title: Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory

abstract: A fast non-convex low-rank matrix decomposition method for potential field data separation is proposed. The singular value decomposition of the large size trajectory matrix, which is also a block Hankel matrix, is obtained using a fast randomized singular value decomposition algorithm in which fast block Hankel matrix-vector multiplications are implemented with minimal memory storage. This fast block Hankel matrix randomized singular value decomposition algorithm is integrated into the \texttt{Altproj} algorithm, which is a standard non-convex method for solving the robust principal component analysis optimization problem. The improved algorithm avoids the construction of the trajectory matrix. Hence, gravity and magnetic data matrices of large size can be computed. Moreover, it is more efficient than the traditional low-rank matrix decomposition method, which is based on the use of an inexact augmented Lagrange multiplier algorithm. The presented algorithm is also robust and, hence, algorithm-dependent parameters are easily determined. The improved and traditional algorithms are contrasted for the separation of synthetic gravity and magnetic data matrices of different sizes. The presented results demonstrate that the improved algorithm is not only computationally more efficient but it is also more accurate. Moreover, it is possible to solve far larger problems. As an example, for the adopted computational environment, matrices of sizes larger than $205 \times 205$ generate "out of memory" exceptions with the traditional method, but a matrix of size $2001\times 2001$ can be calculated in $1062.29$s with the new algorithm. Finally, the improved method is applied to separate real gravity and magnetic data in the Tongling area, Anhui province, China. Areas which may exhibit mineralizations are inferred based on the separated anomalies.


title: On the root count of algebraic Kuramoto equations in cycle networks with uniform coupling

abstract: The Kuramoto model is a classical model used in the study of spontaneous synchronizations in networks of coupled oscillators. In this model, frequency synchronization configurations can be formulated as complex solutions to a system of algebraic equations. Recently, upper bounds to the number of frequency synchronization configurations in cycle networks of N oscillators were calculated under the assumption of generic non-uniform coupling. In this paper, we refine these results for the special cases of uniform coupling. In particular, we show that when, and only when, N is divisible by 4, the upper bound for the number of synchronization configurations in the uniform coupling cases is significantly less than the bound in the non-uniform coupling cases. This result also establishes an explicit formula for the gap between the birationally invariant intersection index and the Bernshtein-Kushnirenko-Khovanskii bound for the underlying algebraic equations.


title: Containment strategies and statistical measures for the control of Bovine Viral Diarrhea spread in livestock trade networks

abstract: Assessing the risk of epidemic spread on networks and developing strategies for its containment is of tremendous practical importance, both due to direct effects in public health and its impact on economies. In this work we present the numerical results of a stochastic, event-driven, hierarchical agent-based model designed to reproduce the infectious dynamics of the cattle disease called Bovine Viral Diarrhea (BVD), for which the corresponding network of movements is the main route of spreading. For the farm-node dynamics, the model takes into account a vast number of breeding, infectious and animal movement mechanisms via a susceptible-infected-recovered (SIR) type of dynamics with an additional permanently infectious class. The interaction between the farms is described by a supply and demand farm manager mechanism governing the network structure and dynamics. We discuss the disease and breeding dynamics, study numerous mitigation strategies of present and past government regulations taking Germany as a case-study country and perform a sensitivity analysis on key parameters. We argue that the model, subject to calibration, has universal predictive potential, can be extended to diseases beyond BVD and demonstrate that appropriate measures can indeed lead to eradication regarding BVD. We further present the results of the time-series of the model and conduct a statistical analysis of and among the different mitigation strategies.


title: Multi-Task Offloading over Vehicular Clouds under Graph-based Representation

abstract: Vehicular cloud computing has emerged as a promising paradigm for realizing user requirements in computation-intensive tasks in modern driving environments. In this paper, a novel framework of multi-task offloading over vehicular clouds (VCs) is introduced where tasks and VCs are modeled as undirected weighted graphs. Aiming to achieve a trade-off between minimizing task completion time and data exchange costs, task components are efficiently mapped to available virtual machines in the related VCs. The problem is formulated as a non-linear integer programming problem, mainly under constraints of limited contact between vehicles as well as available resources, and addressed in low-traffic and rush-hour scenarios. In low-traffic cases, we determine optimal solutions; in rush-hour cases, a connection-restricted randommatching-based subgraph isomorphism algorithm is proposed that presents low computational complexity. Evaluations of the proposed algorithms against greedy-based baseline methods are conducted via extensive simulations.


title: Screening of Informed and Uninformed Experts

abstract: Testing the validity of claims made by self-proclaimed experts can be impossible when testing them in isolation, even with infinite observations at the disposal of the tester. However, in a multiple expert setting it is possible to design a contract that only informed experts accept and uninformed experts reject. The tester can pit competing forecasts of future events against each other and take advantage of the uncertainty experts have about the other experts' knowledge. This contract will work even when there is only a single data point to evaluate.


title: Connectivity concerning the last two subconstituents of a Q-polynomial distance-regular graph

abstract: Let $\Gamma$ be a $Q$-polynomial distance-regular graph of diameter $d\geq 3$. Fix a vertex $\gamma$ of $\Gamma$ and consider the subgraph induced on the union of the last two subconstituents of $\Gamma$ with respect to $\gamma$. We prove that this subgraph is connected.


title: Two-dimensional quantum Yang-Mills theory and the Makeenko-Migdal equations

abstract: These notes, echoing a conference given at the Strasbourg-Zurich seminar in October 2017, are written to serve as an introduction to 2-dimensional quantum Yang-Mills theory and to the results obtained in the last five to ten years about its so-called large N limit.


title: The sum theorem for maximal monotone operators in reflexive Banach spaces revisited

abstract: Our goal is to present a new shorter proof for the maximal monotonicity of the Minkowski sum of two maximal monotone multi-valued operators defined in a reflexive Banach space under the classical interiority condition involving their domains.


title: General Information Bottleneck Objectives and their Applications to Machine Learning

abstract: We view the Information Bottleneck Principle (IBP: Tishby et al., 1999; Schwartz-Ziv and Tishby, 2017) and Predictive Information Bottleneck Principle (PIBP: Still et al., 2007; Alemi, 2019) as special cases of a family of general information bottleneck objectives (IBOs). Each IBO corresponds to a particular constrained optimization problem where the constraints apply to: (a) the mutual information between the training data and the learned model parameters or extracted representation of the data, and (b) the mutual information between the learned model parameters or extracted representation of the data and the test data (if any). The heuristics behind the IBP and PIBP are shown to yield different constraints in the corresponding constrained optimization problem formulations. We show how other heuristics lead to a new IBO, different from both the IBP and PIBP, and use the techniques from (Alemi, 2019) to derive and optimize a variational upper bound on the new IBO. We then apply the theory of general IBOs to resolve the seeming contradiction between, on the one hand, the recommendations of IBP and PIBP to maximize the mutual information between the model parameters and test data, and on the other, recent information-theoretic results (see Xu and Raginsky, 2017) suggesting that this mutual information should be minimized. The key insight is that the heuristics (and thus the constraints in the constrained optimization problems) of IBP and PIBP are not applicable to the scenario analyzed by (Xu and Raginsky, 2017) because the latter makes the additional assumption that the parameters of the trained model have been selected to minimize the empirical loss function. Aided by this insight, we formulate a new IBO that accounts for this property of the parameters of the trained model, and derive and optimize a variational bound on this IBO.


title: Shock waves in a rotating non-Maxwellian magnetized dusty plasma

abstract: A theoretical model is presented to study characteristics of dust acoustic shock in a viscous, magnetized and rotating dusty plasma at both fast and slow time scales. By employing reductive perturbation technique the nonlinear Zakharov--Kuznetsov (ZK) equation has been derived for both cases when dust is inactive and dynamic (fast and slow time scales). Both electrons and ions are considered to follow kappa/Cairns distribution. It is observed that the viscosity in both cases when dust is in background and active plays as a key role in dissipation for the propagation of acoustic shock. Magnetic field and rotation are responsible for the dispersive term. Superthermality has been found to affect significantly on the formation of shock wave along with viscous nature of plasma. The present investigation may be beneficial to understanding the rotating plasma in particular experiments being carried out.


title: Reconfigurable Intelligent Surface Enabled IoT Networks in Generalized Fading Channels

abstract: This paper studies an Internet-of-Things (IoT) network employing a reconfigurable intelligent surface (RIS) over generalized fading channels. Inspired by the promising potential of RIS-based transmission, we investigate a RIS-enabled IoT network with the source node employing a RIS-based access point. The system is modelled with reference to a receiver transmitter pair and the Fisher-Snedecor F model is adopted to analyse the composite fading and shadowing channel. Closed-form expressions are derived for the system with regards to the average capacity, average bit error rate (BER) and outage probability. Monte-Carlo simulations are provided throughout to validate the results. The results investigated and reported in this study extend early results reported in the emerging literature on RIS-enabled technologies and provides a framework for the evaluation of a basic RIS-enabled IoT network over the most common multipath fading channels. The results indicate the clear benefit of employing a RIS-enabled access point, as well as the versatility of the derived expressions in analysing the effects of fading and shadowing on the network. The results further demonstrate that for a RIS-enabled IoT network, there is the need to balance between the cost and benefit of increasing the RIS cells against other parameters such as increasing transmit power, especially at low SNR and/or high to moderate fading/shadowing severity.


title: Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids

abstract: Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, AlB$_2$, and NaZn$_{13}$ in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction-size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

categories: cond-mat.soft cond-mat.stat-mech physics.chem-ph physics.comp-ph


title: A scaling-invariant algorithm for linear programming whose running time depends only on the constraint matrix

abstract: Following the breakthrough work of Tardos in the bit-complexity model, Vavasis and Ye gave the first exact algorithm for linear programming in the real model of computation with running time depending only on the constraint matrix. For solving a linear program (LP) $\max\, c^\top x,\: Ax = b,\: x \geq 0,\: A \in \mathbb{R}^{m \times n}$, Vavasis and Ye developed a primal-dual interior point method using a 'layered least squares' (LLS) step, and showed that $O(n^{3.5} \log (\bar{\chi}_A+n))$ iterations suffice to solve (LP) exactly, where $\bar{\chi}_A$ is a condition measure controlling the size of solutions to linear systems related to $A$. Monteiro and Tsuchiya, noting that the central path is invariant under rescalings of the columns of $A$ and $c$, asked whether there exists an LP algorithm depending instead on the measure $\bar{\chi}^*_A$, defined as the minimum $\bar{\chi}_{AD}$ value achievable by a column rescaling $AD$ of $A$, and gave strong evidence that this should be the case. We resolve this open question affirmatively. Our first main contribution is an $O(m^2 n^2 + n^3)$ time algorithm which works on the linear matroid of $A$ to compute a nearly optimal diagonal rescaling $D$ satisfying $\bar{\chi}_{AD} \leq n(\bar{\chi}^*)^3$. This algorithm also allows us to approximate the value of $\bar{\chi}_A$ up to a factor $n (\bar{\chi}^*)^2$. As our second main contribution, we develop a scaling invariant LLS algorithm, together with a refined potential function based analysis for LLS algorithms in general. With this analysis, we derive an improved $O(n^{2.5} \log n\log (\bar{\chi}^*_A+n))$ iteration bound for optimally solving (LP) using our algorithm. The same argument also yields a factor $n/\log n$ improvement on the iteration complexity bound of the original Vavasis-Ye algorithm.


title: Unconstrained Facial Expression Transfer using Style-based Generator

abstract: Facial expression transfer and reenactment has been an important research problem given its applications in face editing, image manipulation, and fabricated videos generation. We present a novel method for image-based facial expression transfer, leveraging the recent style-based GAN shown to be very effective for creating realistic looking images. Given two face images, our method can create plausible results that combine the appearance of one image and the expression of the other. To achieve this, we first propose an optimization procedure based on StyleGAN to infer hierarchical style vector from an image that disentangle different attributes of the face. We further introduce a linear combination scheme that fuses the style vectors of the two given images and generate a new face that combines the expression and appearance of the inputs. Our method can create high-quality synthesis with accurate facial reenactment. Unlike many existing methods, we do not rely on geometry annotations, and can be applied to unconstrained facial images of any identities without the need for retraining, making it feasible to generate large-scale expression-transferred results.


title: Dark Energy Survey's Observation Strategy, Tactics, and Exposure Scheduler

abstract: The Dark Energy Survey is a stage III dark energy experiment, performing an optical imaging survey to measure cosmological equation of state parameters using four independent methods. The scope and complexity of the survey introduced complex strategic and tactical scheduling problems that needed to be addressed. We begin with an overview of the process used to develop DES strategy and tactics, from the inception of the project, to task forces that studied and developed strategy changes over the course of the survey, to the nightly pre-observing meeting in which immediate tactical issues were addressed. We then summarize the strategic choices made for each sub-survey, including metrics, scheduling considerations, choice of time domain fields and their sequences of exposures, and wide survey footprint and pointing layout choices. We go on to describe the detailed process that determined which specific exposures were taken at which specific times. We give a chronology of the strategic and tactical peculiarities of each year of observing, including the proposal and execution of a sixth year. We give an overview of obstac, the implementation of the DES scheduler used to simulate and evaluate strategic and tactical options, and automate exposure scheduling; and describe developments in obstac for use after DES. Appendices describe further details of data quality evaluation, tau, and t_eff; airmass calculation; and modeling of the seeing and sky brightness. The significant corpus of DES data indicates that the simple scaling relations for seeing as a function of wavelength and airmass derived from the Kolmogorov turbulence model work adequately for exposure planning purposes: deviations from these relations are modest in comparison with short time-scale seeing variations.


title: Theoretically-Efficient and Practical Parallel DBSCAN

abstract: The DBSCAN method for spatial clustering has received significant attention due to its applicability in a variety of data analysis tasks. There are fast sequential algorithms for DBSCAN in Euclidean space that take $O(n\log n)$ work for two dimensions, sub-quadratic work for three or more dimensions, and can be computed approximately in linear work for any constant number of dimensions. However, existing parallel DBSCAN algorithms require quadratic work in the worst case, making them inefficient for large datasets. This paper bridges the gap between theory and practice of parallel DBSCAN by presenting new parallel algorithms for Euclidean exact DBSCAN and approximate DBSCAN that match the work bounds of their sequential counterparts, and are highly parallel (polylogarithmic depth). We present implementations of our algorithms along with optimizations that improve their practical performance. We perform a comprehensive experimental evaluation of our algorithms on a variety of datasets and parameter settings. Our experiments on a 36-core machine with hyper-threading show that we outperform existing parallel DBSCAN implementations by up to several orders of magnitude, and achieve speedups by up to 33x over the best sequential algorithms.


title: On the equivalence between quantum and random walks on finite graphs

abstract: Quantum walks on graphs are ubiquitous in quantum computing finding a myriad of applications. Likewise, random walks on graphs are a fundamental building block for a large number of algorithms with diverse applications. While the relationship between quantum and random walks has been recently discussed in specific scenarios, this work establishes a formal equivalence between the processes on arbitrary finite graphs and general conditions for shift and coin operators. It requires empowering random walks with time heterogeneity, where the transition probability of the walker is non-uniform and time dependent. The equivalence is obtained by equating the probability of measuring the quantum walk on a given node of the graph and the probability that the random walk is at that same node, for all nodes and time steps. The result is given by the construction procedure of a matrix sequence for the random walk that yields the exact same vertex probability distribution sequence of any given quantum walk, including the scenario with multiple interfering walkers. Interestingly, these matrices allows for a different simulation approach for quantum walks where node samples respect neighbor locality and convergence is guaranteed by the law of large numbers, enabling efficient (polynomial) sampling of quantum graph trajectories (paths). Furthermore, the complexity of constructing this sequence of matrices is discussed in the general case.


title: Effects of Opacity Temperature Dependence on Radiatively Accelerated Clouds

abstract: We study how different opacity-temperature scalings affect the dynamical evolution of irradiated gas clouds using time-dependent, radiation-hydrodynamics (rad-HD) simulations. When clouds are optically thick, the bright side heats up and expands, accelerating the cloud via the rocket effect. Clouds that become more optically thick as they heat accelerate $\sim 35\%$ faster than clouds that become optically thin. An enhancement of $\sim 85\%$ in the acceleration can be achieved by having a broken powerlaw opacity profile, which allows the evaporating gas driving the cloud to become optically thin and not attenuate the driving radiation flux. We find that up to $\sim 2\%$ of incident radiation is re-emitted by accelerating clouds, which we estimate as the contribution of a single accelerating cloud to an emission or absorption line. Re-emission is suppressed by "bumps" in the opacity-temperature relation since these decrease the opacity of the hot, evaporating gas, primarily responsible for the re-radiation. If clouds are optically thin, they heat nearly uniformly, expand and form shocks. This triggers the Richtmyer-Meshkov instability, leading to cloud disruption and dissipation on thermal time-scales.


title: Mcity Data Collection for Automated Vehicles Study

abstract: The main goal of this paper is to introduce the data collection effort at Mcity targeting automated vehicle development. We captured a comprehensive set of data from a set of perception sensors (Lidars, Radars, Cameras) as well as vehicle steering/brake/throttle inputs and an RTK unit. Two in-cabin cameras record the human driver's behaviors for possible future use. The naturalistic driving on selected open roads is recorded at different time of day and weather conditions. We also perform designed choreography data collection inside the Mcity test facility focusing on vehicle to vehicle, and vehicle to vulnerable road user interactions which is quite unique among existing open-source datasets. The vehicle platform, data content, tags/labels, and selected analysis results are shown in this paper.


title: A predictive path-following controller for multi-steered articulated vehicles

abstract: Stabilizing multi-steered articulated vehicles in backward motion is a complex task for any human driver. Unless the vehicle is accurately steered, its structurally unstable joint-angle kinematics during reverse maneuvers can cause the vehicle segments to fold and enter a jack-knife state. In this work, a model predictive path-following controller is proposed enabling automatic low-speed steering control of multi-steered articulated vehicles, comprising a car-like tractor and an arbitrary number of trailers with passive or active steering. The proposed path-following controller is tailored to follow nominal paths that contains full state and control-input information, and is designed to satisfy various physical constraints on the vehicle states as well as saturations and rate limitations on the tractor's curvature and the trailer steering angles. The performance of the proposed model predictive path-following controller is evaluated in a set of simulations for a multi-steered 2-trailer with a car-like tractor where the last trailer has steerable wheels.


title: Enhancement of Water Repellence by Hierarchical Surface Structures Integrating Micro-dome and Micro-pillar Arrays with Nanoporous Coatings

abstract: Superhydrophobic surfaces with multi-scale topographies offer exceptionally high apparent water contact angles and low contact angle hysteresis by virtue of the small liquid{\textendash}solid contact fractions they enable. Natural water-repellent surfaces such as lotus leaves often feature dome-shaped micro-scale protrusions, whose lack of sharp edges also facilitates smooth droplet shedding without pinning. Engineered hydrophobic surfaces, however, have not yet fully exploited the merits of protrusions with controlled curvature. In this work, thermal re-flow of photoresist patterns followed by elastomeric casting was used to fabricate arrays of micro-domes with sizes 20{\textendash}50 {\mu}m. These microstructures were coated with a nanoporous zinc oxide film and fluorosilanized to produce hierarchical surface topographies with static water contact angles up to 169.7{\pm}0.4{\deg} and contact angle hysteresis as low as 14.7{\pm}1.3{\deg}. Performance of the micro-dome arrays significantly exceeded that of arrays of sharp-edged square pillars and flat surfaces coated with the same nanoporous film. The highest performance came from the smallest micro-domes (20 {\mu}m) and closest spacings (10 {\mu}m) investigated. For larger features, contact angles reduced and hysteresis increased {\textemdash} unexpected trends not explained by contact fraction alone. This simple fabrication technique could be adapted to manufacture large surfaces for droplet shedding, including in heat transfer applications.


title: Complex poles, spectral function and reflection positivity violation of Yang-Mills theory

abstract: We discuss the analytic continuation of the gluon propagator from the Euclidean region to the complex squared-momentum plane towards the Minkowski region from a viewpoint of gluon confinement. For this purpose, we investigate the massive Yang-Mills model with one-loop quantum corrections, which is to be identified with a low-energy effective theory of the Yang-Mills theory in the sense that the confining decoupling solution for the Euclidean gluon and ghost propagators of the Yang-Mills theory in the Landau gauge obtained by numerical simulations on the lattice are reproduced with good accuracy from the massive Yang-Mills model by taking into account one-loop quantum corrections. We show that the gluon propagator in the massive Yang-Mills model has a pair of complex conjugate poles or "tachyonic" poles of multiplicity two, in accordance with the fact that the gluon field has a negative spectral function, while the ghost propagator has at most one "unphysical" pole. These results are consistent with general relationships between the number of complex poles of a propagator and the sign of the spectral function originating from the branch cut in the Minkowski region under some assumptions on the asymptotic behaviors of the propagator. Consequently, we give an analytical proof for violation of the reflection positivity as a necessary condition for gluon confinement for any choice of the parameters in the massive Yang-Mills model, including the physical point. Moreover, the complex structure of the propagator enables us to explain why the gluon propagator in the Euclidean region is well described by the Gribov-Stingl form.


title: Extracting clinical concepts from user queries

abstract: Clinical concept extraction often begins with clinical Named Entity Recognition (NER). Often trained on annotated clinical notes, clinical NER models tend to struggle with tagging clinical entities in user queries because of the structural differences between clinical notes and user queries. User queries, unlike clinical notes, are often ungrammatical and incoherent. In many cases, user queries are compounded of multiple clinical entities, without comma or conjunction words separating them. By using as dataset a mixture of annotated clinical notes and synthesized user queries, we adapt a clinical NER model based on the BiLSTM-CRF architecture for tagging clinical entities in user queries. Our contribution are the following: 1) We found that when trained on a mixture of synthesized user queries and clinical notes, the NER model performs better on both user queries and clinical notes. 2) We provide an end-to-end and easy-to-implement framework for clinical concept extraction from user queries.


title: On an Analogue Of the Gauss Circle Problem For the Heisenberg Groups

abstract: We consider the problem of estimating the error term $\mathcal{E}_{q}(x)=\big|\mathbb{Z}^{2q+1}\cap\delta_{x}\mathcal{B}\big|-\textit{vol}\big(\mathcal{B}\big)x^{2q+2}$ which occurs in the counting of lattice points in Heisenberg dilates of the Cygan-Kor{\'a}nyi ball. We prove three type of results regarding the order of magnitude of $\mathcal{E}_{q}(x)$, which are valid for any $q\geq3$. An upper bound estimate of the form $|\mathcal{E}_{q}(x)|\ll x^{2q-2/3}$ ; A sharp second moment estimate, which shows that $\mathcal{E}_{q}(x)$ has order of magnitude $x^{2q-1}$ in mean-square ; And an $\Omega$-estimate of the form $\mathcal{E}_{q}(x)=\Omega\big(x^{2q-1}\big(\log{x}\big)^{1/4}\big(\log{\log{x}}\big)^{1/8}\big)$. Consequently, we obtain the lower bound $\kappa_{q}=\sup\big\{\alpha>0:\big|\mathcal{E}_{q}(x)\big|\ll x^{2q+2-\alpha}\big\}\geq\frac{8}{3}$ for $q\geq3$, and conjecture that $\kappa_{q}=3$


title: Optimization-based motion planning for multi-steered articulated vehicles

abstract: The task of maneuvering a multi-steered articulated vehicle in confined environments is difficult even for experienced drivers. In this work, we present an optimization-based trajectory planner targeting low-speed maneuvers in unstructured environments for multi-steered N-trailer vehicles, which are comprised of a car-like tractor and an arbitrary number of interconnected trailers with fixed or steerable wheels. The proposed trajectory planning framework is divided into two steps, where a lattice-based trajectory planner is used in a first step to compute a resolution optimal solution to a discretized version of the trajectory planning problem. The output from the lattice planner is then used in a second step to initialize an optimal control problem solver, which enables the framework to compute locally optimal trajectories that start at the vehicle's initial state and reaches the goal state exactly. The performance of the proposed optimization-based trajectory planner is evaluated in a set of practically relevant scenarios for a multi-steered 3-trailer vehicle with a car-like tractor where the last trailer is steerable.


title: Towards Disentangled Representations for Human Retargeting by Multi-view Learning

abstract: We study the problem of learning disentangled representations for data across multiple domains and its applications in human retargeting. Our goal is to map an input image to an identity-invariant latent representation that captures intrinsic factors such as expressions and poses. To this end, we present a novel multi-view learning approach that leverages various data sources such as images, keypoints, and poses. Our model consists of multiple id-conditioned VAEs for different views of the data. During training, we encourage the latent embeddings to be consistent across these views. Our observation is that auxiliary data like keypoints and poses contain critical, id-agnostic semantic information, and it is easier to train a disentangling CVAE on these simpler views to separate such semantics from other id-specific attributes. We show that training multi-view CVAEs and encourage latent-consistency guides the image encoding to preserve the semantics of expressions and poses, leading to improved disentangled representations and better human retargeting results.


title: Low activity main belt comet 133P/Elst-Pizarro: New constraints on its Albedo, Temperature and Active Mechanism from a thermophysical perspective

abstract: 133P/Elst-Pizarro is the firstly recognized main-belt comet, but we still know little about its nucleus. Firstly we use mid-infrared data of Spitzer-MIPS, Spitzer-IRS and WISE to estimate its effective diameter $D_{\rm eff}=3.9^{+0.4}_{-0.3}$ km, geometric albedo $p_{\rm v}=0.074\pm0.013$ and mean Bond albedo $A_{\rm eff,B}=0.024\pm0.004$. The albedo is used to compute 133P's temperature distribution, which shows significant seasonal variation, especially polar regions, ranging from $\sim40$ to $\sim200$ K. Based on current activity observations, the maximum water gas production rate is estimated to be $\sim1.4\times10^{23}\rm~s^{-1}$, being far weaker than $\sim10^{26}\rm~s^{-1}$ of JFC 67P at similar helio-centric distance $\sim2.7$ AU, indicating a thick dust mantle on the surface to lower down the gas production rate. The diameter of the sublimation area may be $<\sim200$ m according to our model prediction. We thus propose that 133P's activity is more likely to be caused by sublimation of regional near-surface ice patch rather than homogeneous buried ice layer. Such small near-surface ice patch might be exposed by one impact event, before which 133P may be an extinct comet (or ice-rich asteroid) with ice layer buried below $\sim$40 m depth. The proposed ice patch may be located somewhere within latitude $-50\sim50^\circ$ by comparing theoretical variation of sublimation temperature to the constraints from observations. The time scale to form such a thick dust mantle is estimated to be $>100$ Myr, indicating that 133P may be more likely to be a relatively old planetesimals or a member of an old family than a recently formed fragment of some young family.


title: Infrared permittivity of the biaxial van der Waals semiconductor $\alpha$-MoO$_3$ from near- and far-field correlative studies

abstract: The biaxial van der Waals semiconductor $\alpha$-phase molybdenum trioxide ($\alpha$-MoO$_3$) has recently received significant attention due to its ability to support highly anisotropic phonon polaritons (PhPs) -infrared (IR) light coupled to lattice vibrations in polar materials-, offering an unprecedented platform for controlling the flow of energy at the nanoscale. However, to fully exploit the extraordinary IR response of this material, an accurate dielectric function is required. Here, we report the accurate IR dielectric function of $\alpha$-MoO$_3$ by modelling far-field, polarized IR reflectance spectra acquired on a single thick flake of this material. Unique to our work, the far-field model is refined by contrasting the experimental dispersion and damping of PhPs, revealed by polariton interferometry using scattering-type scanning near-field optical microscopy (s-SNOM) on thin flakes of $\alpha$-MoO$_3$, with analytical and transfer-matrix calculations, as well as full-wave simulations. Through these correlative efforts, exceptional quantitative agreement is attained to both far- and near-field properties for multiple flakes, thus providing strong verification of the accuracy of our model, while offering a novel approach to extracting dielectric functions of nanomaterials, usually too small or inhomogeneous for establishing accurate models only from standard far-field methods. In addition, by employing density functional theory (DFT), we provide insights into the various vibrational states dictating our dielectric function model and the intriguing optical properties of $\alpha$-MoO$_3$.


title: Inferring Distributions Over Depth from a Single Image

abstract: When building a geometric scene understanding system for autonomous vehicles, it is crucial to know when the system might fail. Most contemporary approaches cast the problem as depth regression, whose output is a depth value for each pixel. Such approaches cannot diagnose when failures might occur. One attractive alternative is a deep Bayesian network, which captures uncertainty in both model parameters and ambiguous sensor measurements. However, estimating uncertainties is often slow and the distributions are often limited to be uni-modal. In this paper, we recast the continuous problem of depth regression as discrete binary classification, whose output is an un-normalized distribution over possible depths for each pixel. Such output allows one to reliably and efficiently capture multi-modal depth distributions in ambiguous cases, such as depth discontinuities and reflective surfaces. Results on standard benchmarks show that our method produces accurate depth predictions and significantly better uncertainty estimations than prior art while running near real-time. Finally, by making use of uncertainties of the predicted distribution, we significantly reduce streak-like artifacts and improves accuracy as well as memory efficiency in 3D map reconstruction.


title: Learning and Optimization with Bayesian Hybrid Models

abstract: Bayesian hybrid models fuse physics-based insights with machine learning constructs to correct for systematic bias. In this paper, we compare Bayesian hybrid models against physics-based glass-box and Gaussian process black-box surrogate models. We consider ballistic firing as an illustrative case study for a Bayesian decision-making workflow. First, Bayesian calibration is performed to estimate model parameters. We then use the posterior distribution from Bayesian analysis to compute optimal firing conditions to hit a target via a single-stage stochastic program. The case study demonstrates the ability of Bayesian hybrid models to overcome systematic bias from missing physics with less data than the pure machine learning approach. Ultimately, we argue Bayesian hybrid models are an emerging paradigm for data-informed decision-making under parametric and epistemic uncertainty.


title: An Asymptotically Compatible Formulation for Local-to-Nonlocal Coupling Problems without Overlapping Regions

abstract: In this paper we design and analyze an explicit partitioned procedure for a 2D dynamic local-to-nonlocal (LtN) coupling problem, based on a new nonlocal Robin-type transmission condition. The nonlocal subproblem is modeled by the nonlocal heat equation with a finite horizon parameter $\delta$ characterizing the range of nonlocal interactions, and the local subproblem is described by the classical heat equation. We consider a heterogeneous system where the local and nonlocal subproblems present different physical properties, and employ no overlapping region between the two subdomains. We first propose a new generalization of classical local Neumann-type condition by converting the local flux to a correction term in the nonlocal model, and show that the proposed Neumann-type boundary formulation recovers the local case as $O(\delta^2)$ in the $L^{\infty}$ norm. We then extend the nonlocal Neumann-type boundary condition to a Robin-type boundary condition, and develop a local-to-nonlocal coupling formulation with Robin-Dirichlet transmission conditions. To stabilize the explicit coupling procedure and to achieve asymptotic compatibility, the choice of the coefficient in the Robin condition is obtained via amplification factor analysis for the discretized system with coarse grids. Employing a high-order meshfree discretization method in the nonlocal solver and a linear finite element method in the local solver, the selection of optimal Robin coefficients are verified with numerical experiments on heterogeneous and complicated domains. With the developed optimal coupling strategy, we numerically demonstrate the coupling framework's asymptotic convergence to the local limit with an $O(\delta)=O(h)$ rate, when there is a fixed ratio between the horizon size $\delta$ and the spatial discretization size $h$.