title: Data-driven identification of dissipative linear models for nonlinear systems

abstract: We consider the problem of identifying a dissipative linear model of an unknown nonlinear system that is known to be dissipative, from time domain input-output data. We first learn an approximate linear model of the nonlinear system using standard system identification techniques and then perturb the system matrices of the linear model to enforce dissipativity, while closely approximating the dynamical behavior of the nonlinear system. Further, we provide an analytical relationship between the size of the perturbation and the radius in which the dissipativity of the linear model guarantees local dissipativity of the unknown nonlinear system. We demonstrate the application of this identification technique to the problem of learning a dissipative model of a microgrid with high penetration of variable renewable energy sources.


title: Deuterium permeation in Er$_2$O$_3$ thin film fabricated on a type 316L stainless steel substrate

abstract: A metal-oxide film can be used as a hydrogen-isotope permeation barrier in the fuel circulation system for nuclear fusion. We fabricated Er$_2$O$_3$ thin film on a type 316L stainless-steel substrate by using a metal-organic chemical vapor deposition technique for the purpose of hydrogen-isotope permeation barrier. Electron microscopy based imaging and energy-dispersive X-ray spectroscopy measurements indicate a sound film quality together with X-ray diffraction experiments. We also measured deuterium permeation in the film at high temperatures ranging from 600 $^{\circ}$C to 800 $^{\circ}$C. The permeation reduction was most apparent at 650 $^{\circ}$C. Above 800 $^{\circ}$C, we confirmed that the film was damaged and did not work as a permeation barrier.


title: Numerical Contribution of Thermoelectric Polyaniline Sheet

abstract: This work concerns the study of thermo-physical properties and transport phenomena in thermoelectric liquid polyaniline sheet. The electro-thermal heat transfer equations coupled to Naviers-Stokes equation, continuity equation, initial and boundary conditions are solved using the finite difference and the finite element schemes. The results of both numerical techniques show good qualitative and quantitative agreements. The voltage difference, the temperature variation, the Seebeck coefficient, the figure-of-merit and the maximum efficiency of polyaniline sheet are determined numerically. Our results show that the voltage difference is equal to 0.09 mV for non doped polyaniline and 0.20 mV for doped polyaniline DPAN/HCl(I) for a sheet dimension of 192 mm3.


title: Large poroelasto-plastic deformations due to radially outward fluid injection

abstract: Flow-induced failure of granular materials is relevant to a broad range of geomechanical applications. Plasticity, which is the inherent failure mechanism of most granular materials, enables large deformations that can invalidate linearised models. Motivated by fluid injection into a borehole, we develop a steady-state model for the large deformation of a thick-walled, partially-permeable, elastic--perfectly-plastic annulus with a pressurised inner cavity. We account for pre-existing compressive stresses, as would be present in the subsurface, by subtracting a compressed initial state from our solutions to provide the additional disturbance due to fluid injection. We also introduce a simple parameter that allows for a smooth transition from an impermeable material (i.e., subject to mechanical loading at the inner wall) to a fully permeable material (i.e., subject to an internal pore-pressure gradient), which would be relevant to coated boreholes and very-low-permeability materials. We focus on the difference between poroelastic and poroelasto-plastic deformations, the role of kinematic and constitutive nonlinearity, and the transition from impermeable to fully permeable. We find that plasticity can enable much larger deformations while predicting much smaller stresses. The former makes model choice increasingly important in the plastic region, while the elastic region remains insensitive to these changes. We also find that, for a fixed total stress at the inner wall, materials experience larger deformations and generally larger stresses as they transition from impermeable to fully permeable.


title: Sensitivity Comparison of Macro- and Micro-electrochemical Biosensors for Human Chorionic Gonadotropin Biomarker Detection

abstract: Selectivity and sensitivity are important figures of merit in the design and optimization of electrochemical biosensors. The efficiency of the fabricated immunosensing surface can easily be influenced by several factors, such as detection limit, non-specific binding, and type of sensing platform. Here, we demonstrate the effect of macro- and micro-sized planner working electrodes (4 mm and 400 um diameter, respectively) on the electrochemical behavior and the performance of the developed biosensor to detect human chorionic gonadotropin (hCG). The fabricated screen-printed sensor was constructed by modifying the carbon macro- and micro-electrodes with a linker, 1-pyrenebutyric acid-N-hydroxysuccinimide ester (PANHS) and immobilization of anti-hCG antibodies to detect specifically the hCG protein. The characterization of the developed electrodes was performed by cyclic voltammetry (CV) and square wave voltammetry (SWV). Each immunesensing system has its unique electrochemical behavior which might be attributed to arrangement of particles on the surface. However, the smaller surface area of the micro-electrode is found to show higher sensitivity (1 pg/mL) compared to the macro-electrode sensor with a lower detection limit of 100 pg/mL. The proposed assay represents a promising approach that is highly effective for specific detection of an analyte and can be exploited to target biomarkers for a variety of point-of-care diagnostic applications.


title: GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

abstract: The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.


title: Camera Exposure Control for Robust Robot Vision with Noise-Aware Image Quality Assessment

abstract: In this paper, we propose a noise-aware exposure control algorithm for robust robot vision. Our method aims to capture the best-exposed image which can boost the performance of various computer vision and robotics tasks. For this purpose, we carefully design an image quality metric which captures complementary quality attributes and ensures light-weight computation. Specifically, our metric consists of a combination of image gradient, entropy, and noise metrics. The synergy of these measures allows preserving sharp edge and rich texture in the image while maintaining a low noise level. Using this novel metric, we propose a real-time and fully automatic exposure and gain control technique based on the Nelder-Mead method. To illustrate the effectiveness of our technique, a large set of experimental results demonstrates higher qualitative and quantitative performances when compared with conventional approaches.


title: Mapping road safety features from streetview imagery: A deep learning approach

abstract: Each year, around 6 million car accidents occur in the U.S. on average. Road safety features (e.g., concrete barriers, metal crash barriers, rumble strips) play an important role in preventing or mitigating vehicle crashes. Accurate maps of road safety features is an important component of safety management systems for federal or state transportation agencies, helping traffic engineers identify locations to invest on safety infrastructure. In current practice, mapping road safety features is largely done manually (e.g., observations on the road or visual interpretation of streetview imagery), which is both expensive and time consuming. In this paper, we propose a deep learning approach to automatically map road safety features from streetview imagery. Unlike existing Convolutional Neural Networks (CNNs) that classify each image individually, we propose to further add Recurrent Neural Network (Long Short Term Memory) to capture geographic context of images (spatial autocorrelation effect along linear road network paths). Evaluations on real world streetview imagery show that our proposed model outperforms several baseline methods.


title: Multi-Agent Path Finding with Capacity Constraints

abstract: In multi-agent path finding (MAPF) the task is to navigate agents from their starting positions to given individual goals. The problem takes place in an undirected graph whose vertices represent positions and edges define the topology. Agents can move to neighbor vertices across edges. In the standard MAPF, space occupation by agents is modeled by a capacity constraint that permits at most one agent per vertex. We suggest an extension of MAPF in this paper that permits more than one agent per vertex. Propositional satisfiability (SAT) models for these extensions of MAPF are studied. We focus on modeling capacity constraints in SAT-based formulations of MAPF and evaluation of performance of these models. We extend two existing SAT-based formulations with vertex capacity constraints: MDD-SAT and SMT-CBS where the former is an approach that builds the model in an eager way while the latter relies on lazy construction of the model.


title: No More Chasing Waterfalls: A Measurement Study of the Header Bidding Ad-Ecosystem

abstract: In recent years, Header Bidding (HB) has gained popularity among web publishers, challenging the status quo in the ad ecosystem. Contrary to the traditional waterfall standard, HB aims to give back to publishers control of their ad inventory, increase transparency, fairness and competition among advertisers, resulting in higher ad-slot prices. Although promising, little is known about how this ad protocol works: What are HB's possible implementations, who are the major players, and what is its network and UX overhead? To address these questions, we design and implement HBDetector: a novel methodology to detect HB auctions on a website at real time. By crawling 35,000 top Alexa websites, we collect and analyze a dataset of 800k auctions. We find that: (i) 14.28% of top websites utilize HB. (ii) Publishers prefer to collaborate with a few Demand Partners who also dominate the waterfall market. (iii) HB latency can be significantly higher (up to 3x in median case) than waterfall.


title: How to Staff When Customers Arrive in Batches

abstract: In many different settings, requests for service can arrive in near or true simultaneity with one another. This creates batches of arrivals to the underlying queueing system. In this paper, we study the staffing problem for the batch arrival queue. We show that batches place a dangerous and deceptive stress on services, requiring a high amount of resources and exhibiting a fundamentally larger tail in those demands. This uncovers a service regime in which a system with large batch arrivals may have low utilization but will still have non-trivial waiting. Methodologically, these staffing results follow from novel large batch and large batch-and-rate limits of the multi-server queueing model. In the large batch limit, we establish the first formal connection between general multi-server queues and storage processes, another family of stochastic models. By consequence, we show that the batch scaled queue length process is not asymptotically normal, and that, in fact, the fluid and diffusion-type limits coincide. Hence, the (safety) staffing of this system must be directly proportional to the batch size just to achieve a non-degenerate probability of wait. In exhibition of the existence and insights of this large batch regime, we apply our results to data on Covid-19 contact tracing in New York City. In doing so, we identify significant benefits produced by the tracing agency's decision to staff above national recommendations, and we also demonstrate that there may have been an opportunity to further improve the operation by optimizing the arrival pattern in the public health data pipeline.


title: A Physics-Constrained Data-Driven Approach Based on Locally Convex Reconstruction for Noisy Database

abstract: Physics-constrained data-driven computing is an emerging hybrid approach that integrates universal physical laws with data-driven models of experimental data for scientific computing. A new data-driven simulation approach coupled with a locally convex reconstruction, termed the local convexity data-driven (LCDD) computing, is proposed to enhance accuracy and robustness against noise and outliers in data sets in the data-driven computing. In this approach, for a given state obtained by the physical simulation, the corresponding optimum experimental solution is sought by projecting the state onto the associated local convex manifold reconstructed based on the nearest experimental data. This learning process of local data structure is less sensitive to noisy data and consequently yields better accuracy. A penalty relaxation is also introduced to recast the local learning solver in the context of non-negative least squares that can be solved effectively. The reproducing kernel approximation with stabilized nodal integration is employed for the solution of the physical manifold to allow reduced stress-strain data at the discrete points for enhanced effectiveness in the LCDD learning solver. Due to the inherent manifold learning properties, LCDD performs well for high-dimensional data sets that are relatively sparse in real-world engineering applications. Numerical tests demonstrated that LCDD enhances nearly one order of accuracy compared to the standard distance-minimization data-driven scheme when dealing with noisy database, and a linear exactness is achieved when local stress-strain relation is linear.


title: How model accuracy and explanation fidelity influence user trust

abstract: Machine learning systems have become popular in fields such as marketing, financing, or data mining. While they are highly accurate, complex machine learning systems pose challenges for engineers and users. Their inherent complexity makes it impossible to easily judge their fairness and the correctness of statistically learned relations between variables and classes. Explainable AI aims to solve this challenge by modelling explanations alongside with the classifiers, potentially improving user trust and acceptance. However, users should not be fooled by persuasive, yet untruthful explanations. We therefore conduct a user study in which we investigate the effects of model accuracy and explanation fidelity, i.e. how truthfully the explanation represents the underlying model, on user trust. Our findings show that accuracy is more important for user trust than explainability. Adding an explanation for a classification result can potentially harm trust, e.g. when adding nonsensical explanations. We also found that users cannot be tricked by high-fidelity explanations into having trust for a bad classifier. Furthermore, we found a mismatch between observed (implicit) and self-reported (explicit) trust.


title: A new and consistent well model for one-phase flow in anisotropic porous media using a distributed source model

abstract: A new well model for one-phase flow in anisotropic porous media is introduced, where the mass exchange between well and a porous medium is modeled by spatially distributed source terms over a small neighborhood region. To this end, we first present a compact derivation of the exact analytical solution for an arbitrarily oriented, infinite well cylinder in an infinite porous medium with anisotropic permeability tensor in R3 , for constant well pressure and a given injection rate, using a conformal map. The analytical solution motivates the choice of a kernel function to distribute the sources. The presented model is independent from the discretization method and the choice of computational grids. In numerical experiments, the new well model is shown to be consistent and robust with respect to rotation of the well axis, rotation of the permeability tensor, and different anisotropy ratios. Finally, a comparison with a Peaceman-type well model suggests that the new scheme leads to an increased accuracy for injection (and production) rates for arbitrarily-oriented pressure-controlled wells.


title: Selection Principles for Gaia

abstract: The Gaia hypothesis considers the life-environment coupled system as a single entity that acts to regulate and maintain habitable conditions on Earth. In this paper we discuss three mechanisms which could potentially lead to Gaia: Selection by Survival, Sequential Selection and Entropic Hierarchy. We use the Tangled Nature Model of co-evolution as a common framework for investigating all three, using an extended version of the standard model to elaborate on Gaia as an example of an entropic hierarchy. This idea, which combines sequential selection together with a reservoir of diversity that acts as a 'memory', implies a tendency towards growth and increasing resilience of the Gaian system over time. We then discuss how Gaian memory could be realised in practice via the microbial seed bank, climate refugia and lateral gene transfer and conclude by discussing testable implications of an entropic hierarchy for the study of Earth history and the search for life in the universe. This paper adds to the existing taxonomy of Gaia hypotheses to suggest an "Entropic Gaia" where we argue that increasing biomass, complexity and enhanced habitability over time is a statistically likely feature of a co-evolving system.


title: Electronics without bridging components

abstract: We propose a new paradigm of electronic devices based only on two electrodes separated by a gap, i.e. without any functional element bridging them. We use a tight-binding model to show that, depending on the type of material of the electrodes and its structure, several electronic functionalities can be achieved: ohmic behaviour, rectification, negative differential resistance, spin-filtering and magnetoresistance. In particular, we show that it is possible to deliver a given functionality by changing the coupling between the surface and bulk states and between the surface states across the gap, which dramatically changes the current-voltage characteristics. These results prove that it is possible to have functional electronic and spintronic elements on the nanoscale without having physical components bridging the electrodes.


title: Improving MPI Collective I/O Performance With Intra-node Request Aggregation

abstract: Two-phase I/O is a well-known strategy for implementing collective MPI-IO functions. It redistributes I/O requests among the calling processes into a form that minimizes the file access costs. As modern parallel computers continue to grow into the exascale era, the communication cost of such request redistribution can quickly overwhelm collective I/O performance. This effect has been observed from parallel jobs that run on multiple compute nodes with a high count of MPI processes on each node. To reduce the communication cost, we present a new design for collective I/O by adding an extra communication layer that performs request aggregation among processes within the same compute nodes. This approach can significantly reduce inter-node communication congestion when redistributing the I/O requests. We evaluate the performance and compare with the original two-phase I/O on a Cray XC40 parallel computer with Intel KNL processors. Using I/O patterns from two large-scale production applications and an I/O benchmark, we show the performance improvement of up to 29 times when running 16384 MPI processes on 256 compute nodes.


title: On certain polynomial systems involving Stirling numbers of second kind

abstract: We solve a special type of linear systems with coefficients in multivariate polynomial rings. These systems arise in the computation of parametric Bernstein-Sato polynomials associated with certain hypergeometric ideals in the Weyl algebra.


title: On the Solutions of the Lucas-Uzawa Model

abstract: In a recent paper, Naz and Chaudry provided two solutions for the model of Lucas-Uzawa, via the Partial Hamiltonian Approach. The first one of these solutions coincides exactly with that determined by Chilarescu. For the second one, they claim that this is a new solution, fundamentally different than that obtained by Chilarescu. We will prove in this paper, using the existence and uniqueness theorem of nonlinear differential equations, that this is not at all true.


title: Particle Swarm Optimisation for Evolving Deep Neural Networks for Image Classification by Evolving and Stacking Transferable Blocks

abstract: Deep Convolutional Neural Networks (CNNs) have been widely used in image classification tasks, but the process of designing CNN architectures is very complex, so Neural Architecture Search (NAS), automatically searching for optimal CNN architectures, has attracted more and more research interests. However, the computational cost of NAS is often too high to apply NAS on real-life applications. In this paper, an efficient particle swarm optimisation method named EPSOCNN is proposed to evolve CNN architectures inspired by the idea of transfer learning. EPSOCNN successfully reduces the computation cost by minimising the search space to a single block and utilising a small subset of the training set to evaluate CNNs during evolutionary process. Meanwhile, EPSOCNN also keeps very competitive classification accuracy by stacking the evolved block multiple times to fit the whole dataset. The proposed EPSOCNN algorithm is evaluated on CIFAR-10 dataset and compared with 13 peer competitors comprised of deep CNNs crafted by hand, learned by reinforcement learning methods and evolved by evolutionary computation approaches, which shows very promising results by outperforming all of the peer competitors with regard to the classification accuracy, number of parameters and the computational cost.


title: Local Spectroscopies Reveal Percolative Metal in Disordered Mott Insulators

abstract: We elucidate the mechanism by which a Mott insulator transforms into a non-Fermi liquid metal upon increasing disorder at half filling. By correlating maps of the local density of states, the local magnetization and the local bond conductivity, we find a collapse of the Mott gap toward a V-shape pseudogapped density of states that occurs concomitantly with the decrease of magnetism around the highly disordered sites but an increase of bond conductivity. These metallic regions percolate to form an emergent non-Fermi liquid phase with a conductivity that increases with temperature. Bond conductivity measured via local microwave impedance combined with charge and spin local spectroscopies are ideal tools to corroborate our predictions.


title: Compatible Cycles and CHY Integrals

abstract: The CHY construction naturally associates a vector in $\mathbb{R}^{(n-3)!}$ to every 2-regular graph with $n$ vertices. Partial amplitudes in the biadjoint scalar theory are given by the inner product of vectors associated with a pair of cycles. In this work we study the problem of extending the computation to pairs of arbitrary 2-regular graphs. This requires the construction of compatible cycles, i.e. cycles such that their union with a 2-regular graph admits a Hamiltonian decomposition. We prove that there are at least $(n-2)!/4$ such cycles for any 2-regular graph. We also find a connection to breakpoint graphs when the graph only has double edges. We end with a comparison of the lower bound on the number of randomly selected cycles needed to generate a basis of $\mathbb{R}^{(n-3)!}$, using the super Catalan numbers, and our lower bound for compatible cycles.


title: Metastability for the contact process with two types of particles and priorities

abstract: We consider a symmetric finite-range contact process on $\mathbb{Z}$ with two types of particles (or infections), which propagate according to the same supercritical rate and die (or heal) at rate $1$. Particles of type 1 can occupy any site in $(-\infty, 0]$ that is empty or occupied by a particle of type 2 and, analogously, particles of type 2 can occupy any site in $[1,+\infty)$ that is empty or occupied by a particle of type 1. We consider the model restricted to a finite interval $[-N + 1,N] \cap \mathbb{Z}$. If the initial configuration is $\mathbf{1}_ {(-N,0]}+2\mathbf{1}_{[1,N)}$, we prove that this system exhibits two metastable states: one with the two species and the other one with the family that survives the competition.


title: CerebroVis: Designing an Abstract yet Spatially Contextualized Cerebral Arteries Network Visualization

abstract: Blood circulation in the human brain is supplied through a network of cerebral arteries. If a clinician suspects a patient has a stroke or other cerebrovascular condition they order imaging tests. Neuroradiologists visually search the resulting scans for abnormalities. Their visual search tasks correspond to the abstract network analysis tasks of browsing and path following. To assist neuroradiologists in identifying cerebral artery abnormalities we designed CerebroVis, a novel abstract---yet spatially contextualized---cerebral artery network visualization. In this design study, we contribute a novel framing and definition of the cerebral artery system in terms of network theory and characterize neuroradiologist domain goals as abstract visualization and network analysis tasks. Through an iterative, user-centered design process we developed an abstract network layout technique which incorporates cerebral artery spatial context. The abstract visualization enables increased domain task performance over 3D geometry representations, while including spatial context helps preserve the user's mental map of the underlying geometry. We provide open source implementations of our network layout technique and prototype cerebral artery visualization tool. We demonstrate the robustness of our technique by successfully laying out 61 open source brain scans. We evaluate the effectiveness of our layout through a mixed methods study with three neuroradiologists. In a formative controlled experiment our study participants used CerebroVis and a conventional 3D visualization to examine real cerebral artery imaging data and to identify a simulated intracranial artery stenosis. Participants were more accurate at identifying stenoses using CerebroVis (absolute risk difference 13%). A free copy of this paper, the evaluation stimuli and data, and source code are available at https://osf.io/e5sxt/.


title: CUNI Systems for the Unsupervised News Translation Task in WMT 2019

abstract: In this paper we describe the CUNI translation system used for the unsupervised news shared task of the ACL 2019 Fourth Conference on Machine Translation (WMT19). We follow the strategy of Artexte et al. (2018b), creating a seed phrase-based system where the phrase table is initialized from cross-lingual embedding mappings trained on monolingual data, followed by a neural machine translation system trained on synthetic parallel data. The synthetic corpus was produced from a monolingual corpus by a tuned PBMT model refined through iterative back-translation. We further focus on the handling of named entities, i.e. the part of vocabulary where the cross-lingual embedding mapping suffers most. Our system reaches a BLEU score of 15.3 on the German-Czech WMT19 shared task.


title: Airbnb Price Prediction Using Machine Learning and Sentiment Analysis

abstract: Pricing a rental property on Airbnb is a challenging task for the owner as it determines the number of customers for the place. On the other hand, customers have to evaluate an offered price with minimal knowledge of an optimal value for the property. This paper aims to develop a reliable price prediction model using machine learning, deep learning, and natural language processing techniques to aid both the property owners and the customers with price evaluation given minimal available information about the property. Features of the rentals, owner characteristics, and the customer reviews will comprise the predictors, and a range of methods from linear regression to tree-based models, support-vector regression (SVR), K-means Clustering (KMC), and neural networks (NNs) will be used for creating the prediction model.


title: Modeling Shared Cache Performance of OpenMP Programs using Reuse Distance

abstract: Performance modeling of parallel applications on multicore computers remains a challenge in computational co-design due to the complex design of multicore processors including private and shared memory hierarchies. We present a Scalable Analytical Shared Memory Model to predict the performance of parallel applications that runs on a multicore computer and shares the same level of cache in the hierarchy. This model uses a computationally efficient, probabilistic method to predict the reuse distance profiles, where reuse distance is a hardware architecture-independent measure of the patterns of virtual memory accesses. It relies on a stochastic, static basic block-level analysis of reuse profiles measured from the memory traces of applications ran sequentially on small instances rather than using a multi-threaded trace. The results indicate that the hit-rate predictions on the shared cache are accurate.


title: Reinforced Dynamic Reasoning for Conversational Question Generation

abstract: This paper investigates a new task named Conversational Question Generation (CQG) which is to generate a question based on a passage and a conversation history (i.e., previous turns of question-answer pairs). CQG is a crucial task for developing intelligent agents that can drive question-answering style conversations or test user understanding of a given passage. Towards that end, we propose a new approach named Reinforced Dynamic Reasoning (ReDR) network, which is based on the general encoder-decoder framework but incorporates a reasoning procedure in a dynamic manner to better understand what has been asked and what to ask next about the passage. To encourage producing meaningful questions, we leverage a popular question answering (QA) model to provide feedback and fine-tune the question generator using a reinforcement learning mechanism. Empirical results on the recently released CoQA dataset demonstrate the effectiveness of our method in comparison with various baselines and model variants. Moreover, to show the applicability of our method, we also apply it to create multi-turn question-answering conversations for passages in SQuAD.


title: Perspectives on CUR Decompositions

abstract: This note discusses an interesting matrix factorization called the CUR Decomposition. We illustrate various viewpoints of this method by comparing and contrasting them in different situations. Additionally, we offer a new characterization of CUR decompositions which synergizes these viewpoints and shows that they are indeed the same in the exact decomposition case.


title: The Challenge of Imputation in Explainable Artificial Intelligence Models

abstract: Explainable models in Artificial Intelligence are often employed to ensure transparency and accountability of AI systems. The fidelity of the explanations are dependent upon the algorithms used as well as on the fidelity of the data. Many real world datasets have missing values that can greatly influence explanation fidelity. The standard way to deal with such scenarios is imputation. This can, however, lead to situations where the imputed values may correspond to a setting which refer to counterfactuals. Acting on explanations from AI models with imputed values may lead to unsafe outcomes. In this paper, we explore different settings where AI models with imputation can be problematic and describe ways to address such scenarios.


title: Minimum coprime graph labelings

abstract: A coprime labeling of a graph $G$ is a labeling of the vertices of $G$ with distinct integers from $1$ to $k$ such that adjacent vertices have coprime labels. The minimum coprime number of $G$ is the least $k$ for which such a labeling exists. In this paper, we determine the minimum coprime number for several well-studied classes of graphs, including the coronas of complete graphs with empty graphs and the joins of two paths. In particular, we resolve a conjecture of Seoud, El Sonbaty, and Mahran and two conjectures of Asplund and Fox. We also provide an asymptotic for the minimum coprime number of the Erd\H{o}s-R\'enyi random graph.


title: Periodic Ultranarrow Rods as 1D Subwavelength Optical Lattices

abstract: We report on ground state properties of a one-dimensional, weakly-interacting Bose gas constrained by an infinite multi-rods periodic structure at zero temperature. We solve the stationary Gross-Pitaevskii equation (GPE) to obtain the Bloch wave functions from which we give a semi-analytical solution for the density profile, as well as for the phase of the wave function in terms of the Jacobi elliptic functions, and the incomplete elliptic integrals of the first, second and third kind. Then, we determine numerically the energy of the ground state, the chemical potential and the compressibility of the condensate and show their dependence on the potential height, as well as on the interaction between the bosons. We show the appearance of loops in the energy band spectrum of the system for strong enough interactions, which appear at the edges of the first Brillouin zone for odd bands and at the center for even bands. We apply our model to predict the energy band structure of the Bose gas in an optical lattice with subwavelength spatial structure. To discuss the density range of the validity of the GPE predictions, we calculate the ground state energies of the free Bose gas using the GPE, which we compare with the Lieb-Liniger exact energies.


title: Correction of Channel Sounding Clock Drift and Antenna Rotation Effects for mmWave Angular Profile Measurements

abstract: Millimeter-wave (mmWave) bands will be used for the fifth generation communication systems to support high data rates. For the proper characterization of the mmWave propagation channel, it is essential to measure the power angular-delay profile (PADP) of the channel which includes angle-of-departure (AoD) and angle-of-arrival (AoA) of the multipath components (MPCs). In this paper, we first describe in detail our 28 GHz channel sounder where directional horn antennas are placed on rotating gimbals. Then, for this specific sounder class, we describe and address the following two problems in extracting the MPCs from the measurements: 1) For the channel measurements at large distances between the transmitter (TX) and the receiver (RX), it is not possible to generate the triggering signal for the TX and the RX using a single clock (SICL). This necessitates the use of separate clocks (SECLs) which introduces a random timing drift between the clocks. 2) As the positions of the antennas change during the scanning process, total distance traveled by the same MPC differs at each measurement. These two errors together cause missing some of the MPCs and detecting MPCs that do not exist in reality. We propose an algorithm to correct the clock drift and the errors in the MPC delays due to the rotation of the antennas. We compare the MPCs from the SICL measurement and the corrected SECL measurements using a Hungarian algorithm based MPC matching method. We show that the percentage of the matched MPCs increases from 28.36% to 74.13% after the correction process.


title: Dynamics of solar Coronal Mass Ejections: forces that impact their propagation

abstract: We investigate the Sun-Earth dynamics of a set of 38 well-observed Coronal Mass Ejections(CMEs) using data from the STEREO, SOHO missions and WIND instrument. We seek to quantify the relative contributions of Lorentz force and aerodynamic drag on their propagation. The CMEs are 3D reconstructed using the Graduated Cylindrical Shell(GCS) model to derive observed CME parameters.Using a microphysical prescription of the drag coefficient,we find that solar wind aerodynamic drag adequately accounts for the dynamics of the fastest CMEs from as low as 3.5 Rs. For relatively slower CMEs, however, we find that when the drag-based model is initiated below the distances ranging from 12 to 50 Rs,the observed CME trajectories cannot be accounted for.This suggests that aerodynamic drag force dominates the dynamics of slower CMEs only above these heights.To investigate CME dynamics below the heights where aerodynamic drag dominates, we consider the Torus Instability model for the driving Lorentz force. We find that for fast CMEs,Lorentz forces become negligible in comparison to aerodynamic drag as early as 3.5-4 Rs.For slow CMEs, however, they become negligible only by 12-50 Rs.This justifies the success of the drag-only model for fast CMEs. In case of slow CMEs, the Lorentz force is only slightly smaller than the drag force even beyond 12-50 Rs.For these slow events,our results suggest that some of the magnetic flux carried by CMEs might be expended in expansion or heating. These dissipation effects might be important in describing the propagation of slower CMEs. To the best of our knowledge,this is the first systematic study in this regard using a diverse CME sample. A physical understanding of the forces that affect CME propagation and how they compare with each other at various heliocentric distances is an important ingredient in building tools for describing and predicting CME trajectories.


title: One-to-X analogical reasoning on word embeddings: a case for diachronic armed conflict prediction from news texts

abstract: We extend the well-known word analogy task to a one-to-X formulation, including one-to-none cases, when no correct answer exists. The task is cast as a relation discovery problem and applied to historical armed conflicts datasets, attempting to predict new relations of type `location:armed-group' based on data about past events. As the source of semantic information, we use diachronic word embedding models trained on English news texts. A simple technique to improve diachronic performance in such task is demonstrated, using a threshold based on a function of cosine distance to decrease the number of false positives; this approach is shown to be beneficial on two different corpora. Finally, we publish a ready-to-use test set for one-to-X analogy evaluation on historical armed conflicts data.


title: Neutrino counting experiments and non-unitarity from LEP and future experiments

abstract: Non-unitarity of the neutrino mixing matrix is expected in many scenarios with physics beyond the Standard Model. Motivated by the search for deviations from unitary, we study two neutrino counting observables: the neutrino-antineutrino gamma process and the invisible $Z$ boson decay into neutrinos. We report on new constraints for non-unitarity coming from the first of this observables. We study the potential constraints that future collider experiments will give from the invisible decay of the Z boson, that will be measured with improved precision.


title: Relativistic stable processes in quasi-ballistic heat conduction in thin film semiconductors

abstract: In this article, we show how relativistic alpha stable processes can be used to explain quasi-ballistic heat conduction in semiconductors. This is a method that can fit experimental results of ultrafast laser heating in alloys. It also provides a connection to a rich literature on Feynman-Kac formalism and random processes that transition from a stable L\'evy process on short time and length scales to the Brownian motion at larger scales. This transition was captured by a heuristic truncated L\'evy distribution in earlier papers. The rigorous Feynman-Kac approach is used to derive sharp bounds for the transition kernel. Future directions are briefly discussed.


title: Cal X-1: an absolute in-orbit calibrator for current and future X-ray observatories

abstract: We describe Cal X-1, a SmallSat mission concept to establish X-ray standard candles in the sky, which will enable absolute calibration of the current and future X-ray observatories. It consists of two CubeSats flying in formation, one containing an absolutely calibrated X-ray source and another, 1-2 km away, a simple X-ray telescope.


title: Analog Errors in Quantum Annealing: Doom and Hope

abstract: Quantum annealing has the potential to provide a speedup over classical algorithms in solving optimization problems. Just as for any other quantum device, suppressing Hamiltonian control errors will be necessary before quantum annealers can achieve speedups. Such analog control errors are known to lead to $J$-chaos, wherein the probability of obtaining the optimal solution, encoded as the ground state of the intended Hamiltonian, varies widely depending on the control error. Here, we show that $J$-chaos causes a catastrophic failure of quantum annealing, in that the scaling of the time-to-solution metric becomes worse than that of a deterministic (exhaustive) classical solver. We demonstrate this empirically using random Ising spin glass problems run on the two latest generations of the D-Wave quantum annealers. We then proceed to show that this doomsday scenario can be mitigated using a simple error suppression and correction scheme known as quantum annealing correction (QAC). By using QAC, the time-to-solution scaling of the same D-Wave devices is improved to below that of the classical upper bound, thus restoring hope in the speedup prospects of quantum annealing.


title: Machine Translation Evaluation with BERT Regressor

abstract: We introduce the metric using BERT (Bidirectional Encoder Representations from Transformers) (Devlin et al., 2019) for automatic machine translation evaluation. The experimental results of the WMT-2017 Metrics Shared Task dataset show that our metric achieves state-of-the-art performance in segment-level metrics task for all to-English language pairs.


title: Mergers, Starbursts, and Quenching in the Simba Simulation

abstract: We use the Simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\leq$ 4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass $M_*$ well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate sSFR$>t_{H}^{-1}$ to sSFR$<0.2t_{H}^{-1}$, where $t_{H}$ is the Hubble time. At $z\approx 0-3$, mergers show $\sim\times 2-3$ higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\leq 1\%$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the $H_2$ content, but for $M_*\geq 10^{10.5} M_{\odot}$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift $\propto (1+z)^{3.5}$, but the quenching rate evolves much more slowly, $\propto (1+z)^{0.9}$; there are insufficient mergers to explain the quenching rate at $z\leq 1.5$. Simba first quenches galaxies at $z\geq 3$, with a number density in good agreement with observations. The quenching timescales $\tau_q$ are strongly bimodal, with `slow' quenchings ($\tau_q \sim 0.1t_{H}$) dominating overall, but `fast' quenchings ($\tau_q\sim 0.01 t_H$) dominating in $M_*\sim 10^{10}-10^{10.5}M_{\odot}$ galaxies, likely induced by Simba's jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, Simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.


title: Residual-Guided In-Loop Filter Using Convolution Neural Network

abstract: The block-based coding structure in the hybrid video coding framework inevitably introduces compression artifacts such as blocking, ringing, etc. To compensate for those artifacts, extensive filtering techniques were proposed in the loop of video codecs, which are capable of boosting the subjective and objective qualities of reconstructed videos. Recently, neural network based filters were presented with the power of deep learning from a large magnitude of data. Though the coding efficiency has been improved from traditional methods in High-Efficiency Video Coding (HEVC), the rich features and information generated by the compression pipeline has not been fully utilized in the design of neural networks. Therefore, in this paper, we propose the Residual-Reconstruction-based Convolutional Neural Network (RRNet) to further improve the coding efficiency to its full extent, where the compression features induced from bitstream in form of prediction residual is fed into the network as an additional input to the reconstructed frame. In essence, the residual signal can provide valuable information about block partitions and can aid reconstruction of edge and texture regions in a picture. Thus, more adaptive parameters can be trained to handle different texture characteristics. The experimental results show that our proposed RRNet approach presents significant BD-rate savings compared to HEVC and the state-of-the-art CNN-based schemes, indicating that residual signal plays a significant role in enhancing video frame reconstruction.


title: Conditioned two-dimensional simple random walk: Green's function and harmonic measure

abstract: We study the Doob's $h$-transform of the two-dimensional simple random walk with respect to its potential kernel, which can be thought of as the two-dimensional simple random walk conditioned on never hitting the origin. We derive an explicit formula for the Green's function of this random walk, and also prove a quantitative result on the speed of convergence of the (conditional) entrance measure to the harmonic measure (for the conditioned walk) on a finite set.


title: Proof of the circulant Hadamard conjecture

abstract: In this paper the circulant Hadamard conjecture is proved.


title: Analytical percolation theory for topological color codes under qubit loss

abstract: Quantum information theory has shown strong connections with classical statistical physics. For example, quantum error correcting codes like the surface and the color code present a tolerance to qubit loss that is related to the classical percolation threshold of the lattices where the codes are defined. Here we explore such connection to study analytically the tolerance of the color code when the protocol introduced in [Phys. Rev. Lett. $\textbf{121}$, 060501 (2018)] to correct qubit losses is applied. This protocol is based on the removal of the lost qubit from the code, a neighboring qubit, and the lattice edges where these two qubits reside. We first obtain analytically the average fraction of edges $ r(p) $ that the protocol erases from the lattice to correct a fraction $ p $ of qubit losses. Then, the threshold $ p_c $ below which the logical information is protected corresponds to the value of $ p $ at which $ r(p) $ equals the bond-percolation threshold of the lattice. Moreover, we prove that the logical information is protected if and only if the set of lost qubits does not include the entire support of any logical operator. The results presented here open a route to an analytical understanding of the effects of qubit losses in topological quantum error codes.


title: A time-dependent diffusion MRI signature of axon caliber variations and beading

abstract: MRI provides a unique non-invasive window into the brain, yet is limited to millimeter resolution, orders of magnitude coarser than cell dimensions. Here we show that diffusion MRI is sensitive to the micrometer-scale variations in axon caliber or pathological beading, by identifying a signature power-law diffusion time-dependence of the along-fiber diffusion coefficient. We observe this signature in human brain white matter, and uncover its origins by Monte Carlo simulations in realistic substrates from 3d electron microscopy of mouse corpus callosum. Simulations reveal that the time-dependence originates from axon caliber variation, rather than from mitochondria or axonal undulations. We report a decreased amplitude of time-dependence in multiple sclerosis lesions, illustrating the sensitivity of our method to axonal beading in a plethora of neurodegenerative disorders. This specificity to microstructure offers an exciting possibility of bridging across scales to image cellular-level pathology with a clinically feasible MRI technique.


title: Concentration of measure, classification of submeasures, and dynamics of $L_{0}$

abstract: Exhibiting a new type of measure concentration, we prove uniform concentration bounds for measurable Lipschitz functions on product spaces, where Lipschitz is taken with respect to the metric induced by a weighted covering of the index set of the product. Our proof combines the Herbst argument with an analogue of Shearer's lemma for differential entropy. We give a quantitative "geometric" classification of diffuse submeasures into elliptic, parabolic, and hyperbolic. We prove that any non-elliptic submeasure (for example, any measure, or any pathological submeasure) has a property that we call covering concentration. Our results have strong consequences for the dynamics of the corresponding topological $L_{0}$-groups.


title: Let's Make It Personal, A Challenge in Personalizing Medical Inter-Human Communication

abstract: Current AI approaches have frequently been used to help personalize many aspects of medical experiences and tailor them to a specific individuals' needs. However, while such systems consider medically-relevant information, they ignore socially-relevant information about how this diagnosis should be communicated and discussed with the patient. The lack of this capability may lead to mis-communication, resulting in serious implications, such as patients opting out of the best treatment. Consider a case in which the same treatment is proposed to two different individuals. The manner in which this treatment is mediated to each should be different, depending on the individual patient's history, knowledge, and mental state. While it is clear that this communication should be conveyed via a human medical expert and not a software-based system, humans are not always capable of considering all of the relevant aspects and traversing all available information. We pose the challenge of creating Intelligent Agents (IAs) to assist medical service providers (MSPs) and consumers in establishing a more personalized human-to-human dialogue. Personalizing conversations will enable patients and MSPs to reach a solution that is best for their particular situation, such that a relation of trust can be built and commitment to the outcome of the interaction is assured. We propose a four-part conceptual framework for personalized social interactions, expand on which techniques are available within current AI research and discuss what has yet to be achieved.


title: Isolation solution for extreme environmental vibrations for quantum-enabling cryogenic setups installed on raised frames

abstract: Cryogenic quantum sensing techniques are developing alongside the ever-increasing requirements for noiseless experimental environments. For instance, several groups have isolated internal system vibrations from cold heads in closed-cycle dilution refrigerators. However, these solutions often do not account for external vibrations, necessitating novel strategies to isolate the entire cryogenic systems from their environments in a particular set of raised cryostats. Here, we introduce a dual-stage external active vibration-isolation solution in conjunction with a closed-cycle dilution refrigerator that isolates it from the environment. This dual stage includes two sets of active attenuators and a customized steel tower for supporting experimental probes at heights of 3 m from the floor. Both stages achieve 20-40 dB of attenuation with the active systems engaged, corresponding to levels of vibration in the VC-G range (a standardized Vibration Criterion appropriate for extremely quiet research spaces) on the cryostat's room temperature baseplate and the steel tower. Our unique vibration isolation solution therefore expands the applications of modern cryogenic equipment beyond exclusively quiet specialty buildings, rendering such equipment suitable for interdisciplinary, open-floor research centers.


title: Van der Waals--Allen--Cahn--Hilliard equation with a volume constraint

abstract: We give multiplicity results for the solutions of a nonlinear elliptic equation, with an asymmetric double well potential of Van der Waals-Allen--Cahn--Hilliard type, satisfying a linear volume constraint, on a bounded Lipschitz domain $\Omega\subset\mathds R^N$. The number of solutions is estimated in terms of topological and homological invariants of the underlying domain $\Omega$.


title: Control of nonlinear, complex and black-boxed greenhouse system with reinforcement learning

abstract: Modern control theories such as systems engineering approaches try to solve nonlinear system problems by revelation of causal relationship or co-relationship among the components; most of those approaches focus on control of sophisticatedly modeled white-boxed systems. We suggest an application of actor-critic reinforcement learning approach to control a nonlinear, complex and black-boxed system. We demonstrated this approach on artificial green-house environment simulator all of whose control inputs have several side effects so human cannot figure out how to control this system easily. Our approach succeeded to maintain the circumstance at least 20 times longer than PID and Deep Q Learning.


title: Feasible bases for a polytope related to the Hamilton cycle problem

abstract: We study a certain polytope depending on a graph $G$ and a parameter $\beta\in(0,1)$ which arises from embedding the Hamiltonian cycle problem in a discounted Markov decision process. Eshragh \emph{et al.} conjectured a lower bound on the proportion of feasible bases corresponding to Hamiltonian cycles in the set of all feasible bases. We make progress towards a proof of the conjecture by proving results about the structure of feasible bases. In particular, we prove three main results: (1) the set of feasible bases is independent of the parameter $\beta$ when the parameter is close to 1, (2) the polytope can be interpreted as a generalized network flow polytope and (3) we deduce a combinatorial interpretation of the feasible bases. We also provide a full characterization for a special class of feasible bases, and we apply this to provide some computational support for the conjecture.


title: On the derivatives of the Heun functions

abstract: The Heun functions satisfy linear ordinary differential equations of second order with certain singularities in the complex plane. The first order derivatives of the Heun functions satisfy linear second order differential equations with one more singularity. In this paper we compare these equations with linear differential equations isomonodromy deformations of which are described by the Painlev\'e equations $P_{II}-P_{VI}$.


title: The effects of a compressive velocity pulse on a collapsing turbulent clump

abstract: High-resolution hydrodynamical simulations are presented to follow the gravitational collapse of a uniform turbulent clump, upon which a purely radial compressive velocity pulse is activated in the midst of the evolution of the clump, when its turbulent state has been fully developed. The shape of the velocity pulse is determined basically by two free parameters: the velocity $V_0$ and the initial radial position $r_0$. In the present paper, models are considered in which the velocity $V_0$ takes the values 2, 5, 10, 20, and 50 times the speed of sound of the clump $c_0$, while $r_0$ is fixed for all the models. The collapse of the model with $2 \, c_0$ goes faster as a consequence of the velocity pulse, while the cluster formed in the central region of the isolated clump mainly stays the same. In the models with greater velocity $V_0$, the evolution of the isolated clump is significantly changed, so that a dense shell of gas forms around $r_0$ and moves radially inward. The radial profile of the density and velocity as well as the mass contained in the dense shell of gas are calculated, and it is found that (i) the higher the velocity $V_0$, the less mass is contained in the shell; (ii) there is a critical velocity of the pulse, around $10 \, c_0$, such that for shock models with a lower velocity, there will be a well defined dense central region in the shocked clump surrounded by the shell.


title: Diffusive bounds for the critical density of activated random walks

abstract: We consider symmetric activated random walks on $\mathbb{Z}$, and show that the critical density $\zeta_c$ satisfies $c\sqrt{\lambda} \leq \zeta_c(\lambda) \leq C \sqrt{\lambda}$ where $\lambda$ denotes the sleep rate.


title: Nonlinear causality of general first-order relativistic viscous hydrodynamics

abstract: Effective theory arguments are used to derive the most general energy-momentum tensor of a relativistic viscous fluid with an arbitrary equation of state (in the absence of other conserved currents) that is first-order in the derivatives of the energy density and flow velocity and does not include extended variables such as in Mueller-Israel-Stewart-like theories. This energy-momentum tensor leads to a causal theory, provided one abandons the usual conventions for the out-of-equilibrium hydrodynamic variables put forward by Landau-Lifshitz and Eckart. In particular, causality requires nonzero out-of-equilibrium energy density corrections and heat flow. Conditions are found to ensure linear stability around equilibrium in flat space-time. We also prove local existence and uniqueness of solutions to the equations of motion. Our causality, existence, and uniqueness results hold in the full nonlinear regime, without symmetry assumptions, in four space-time dimensions, with or without coupling to Einstein's equations, and are mathematically rigorously established. Furthermore, a kinetic theory realization of this energy-momentum tensor is also provided.


title: Multiple transitions between normal and hyperballistic diffusion in quantum walks with time-dependent jumps

abstract: We extend to the gamut of functional forms of the probability distribution of the time-dependent step-length a previous model dubbed Elephant Quantum Walk, which considers a uniform distribution and yields hyperballistic dynamics where the variance grows cubicly with time, $\sigma ^2 \propto t^3$, and a Gaussian for the position of the walker. We investigate this proposal both locally and globally with the results showing that the time-dependent interplay between interference, memory and long-range hopping leads to multiple transitions between dynamical regimes, namely ballistic $\rightarrow$ diffusive $\rightarrow$ superdiffusive $\rightarrow$ ballistic $\rightarrow$ hyperballistic for non-hermitian coin whereas the first diffusive regime is quelled for implementations using the Hadamard coin. In addition, we observe a robust asymptotic approach to maximal coin-space entanglement.


title: Dual-FOFE-net Neural Models for Entity Linking with PageRank

abstract: This paper presents a simple and computationally efficient approach for entity linking (EL), compared with recurrent neural networks (RNNs) or convolutional neural networks (CNNs), by making use of feedforward neural networks (FFNNs) and the recent dual fixed-size ordinally forgetting encoding (dual-FOFE) method to fully encode the sentence fragment and its left/right contexts into a fixed-size representation. Furthermore, in this work, we propose to incorporate PageRank based distillation in our candidate generation module. Our neural linking models consist of three parts: a PageRank based candidate generation module, a dual-FOFE-net neural ranking model and a simple NIL entity clustering system. Experimental results have shown that our proposed neural linking models achieved higher EL accuracy than state-of-the-art models on the TAC2016 task dataset over the baseline system, without requiring any in-house data or complicated handcrafted features. Moreover, it achieves a competitive accuracy on the TAC2017 task dataset.


title: EVO* 2019 -- Late-Breaking Abstracts Volume

abstract: This volume contains the Late-Breaking Abstracts submitted to the EVO* 2019 Conference, that took place in Leipzig, from 24 to 26 of April. These papers where presented as short talks and also at the poster session of the conference together with other regular submissions. All of them present ongoing research and preliminary results investigating on the application of different approaches of Evolutionary Computation to different problems, most of them real world ones.


title: Towards Dynamic Simulation Guided Optimal Design of Tumbling Microrobots

abstract: Design of robots at the small scale is a trial-and-error based process, which is costly and time-consuming. There are no good dynamic simulation tools to predict the motion or performance of a microrobot as it moves against a substrate. At smaller length scales, the influence of adhesion and friction, which scales with surface area, becomes more pronounced. Thus, rigid body dynamic simulators, which implicitly assume that contact between two bodies can be modeled as point contact are not suitable. In this paper, we present techniques for simulating the motion of microrobots where there can be intermittent and non-point contact between the robot and the substrate. We use this simulator to study the motion of microrobots of different shapes and select shapes that are most promising for performing a given task.


title: A Practical Guide for the Effective Evaluation of Twitter User Geolocation

abstract: Geolocating Twitter users---the task of identifying their home locations---serves a wide range of community and business applications such as managing natural crises, journalism, and public health. Many approaches have been proposed for automatically geolocating users based on their tweets; at the same time, various evaluation metrics have been proposed to measure the effectiveness of these approaches, making it challenging to understand which of these metrics is the most suitable for this task. In this paper, we propose a guide for a standardized evaluation of Twitter user geolocation by analyzing fifteen models and two baselines in a controlled experimental setting. Models are evaluated using ten metrics over four geographic granularities. We use rank correlations to assess the effectiveness of these metrics. Our results demonstrate that the choice of effectiveness metric can have a substantial impact on the conclusions drawn from a geolocation system experiment, potentially leading experimenters to contradictory results about relative effectiveness. We show that for general evaluations, a range of performance metrics should be reported, to ensure that a complete picture of system effectiveness is conveyed. Given the global geographic coverage of this task, we specifically recommend evaluation at micro versus macro levels to measure the impact of the bias in distribution over locations. Although a lot of complex geolocation algorithms have been applied in recent years, a majority class baseline is still competitive at coarse geographic granularity. We propose a suite of statistical analysis tests, based on the employed metric, to ensure that the results are not coincidental.


title: Multi-critical Behavior in Topological Phase Transitions

abstract: Topological phase transitions can be described by the theory of critical phenomena and identified by critical exponents that define their universality classes. This is a consequence of the existence of a diverging length at the transition that has been identified as the penetration depth of the surface modes in the non-trivial topological phase. In this paper, we characterize different universality classes of topological transitions by determining their correlation length exponents directly from numerical calculations of the penetration length of the edge modes as a function of the distance to the topological transition. We consider generalizations of the topological non-trivial Su-Schrieefer-Heeger (SSH) model, for the case of next nearest neighbors hopping and in the presence of a synthetic potential. The latter allows the system to transit between two universality classes with different correlation length and dynamic critical exponents. It presents a multi-critical point in its phase diagram since the behavior of the Berry connection depends on the path it is approached. We compare our results with those obtained from a scaling approach to the Berry connection.


title: On an optimal quadrature formula for approximation of Fourier integrals in the space $L_2^{(1)}$

abstract: This paper deals with the construction of an optimal quadrature formula for the approximation of Fourier integrals in the Sobolev space $L_2^{(1)}[a,b]$ of non-periodic, complex valued functions which are square integrable with first order derivative. Here the quadrature sum consists of linear combination of the given function values in a uniform grid. The difference between the integral and the quadrature sum is estimated by the norm of the error functional. The optimal quadrature formula is obtained by minimizing the norm of the error functional with respect to coefficients. Analytic formulas for optimal coefficients can also be obtained using discrete analogue of the differential operator $d^2/d x^2$. In addition, the convergence order of the optimal quadrature formula is studied. It is proved that the obtained formula is exact for all linear polynomials. Thus, it is shown that the convergence order of the optimal quadrature formula for functions of the space $C^2[a,b]$ is $O(h^2)$. Moreover, several numerical results are presented and the obtained optimal quadrature formula is applied to reconstruct the X-ray Computed Tomography image by approximating Fourier transforms.


title: The Symmetric $2\times 2$ Hypergeometric Matrix Differential Operators

abstract: We obtain an explicit classification of all $2\times 2$ real hypergeometric Bochner pairs, ie. pairs $(W(x),\mathfrak{D})$ consisting of a $2\times 2$ real hypergeometric differential operator $\mathfrak{D}$ and a $2\times 2$ weight matrix satisfying the property that $\mathfrak{D}$ is symmetric with respect to the matrix-valued inner product defined by W(x). Furthermore, we obtain a classifying space of hypergeometric Bochner pairs by describing a bijective correspondence between the collection of pairs and an open subset of a real algebraic set whose smooth paths correspond to isospectral deformations of the weight W(x) preserving a bispectral property. We also relate the hypergeometric Bochner pairs to classical Bochner pairs via noncommutative bispectral Darboux transformations.


title: Multi-Angle Point Cloud-VAE: Unsupervised Feature Learning for 3D Point Clouds from Multiple Angles by Joint Self-Reconstruction and Half-to-Half Prediction

abstract: Unsupervised feature learning for point clouds has been vital for large-scale point cloud understanding. Recent deep learning based methods depend on learning global geometry from self-reconstruction. However, these methods are still suffering from ineffective learning of local geometry, which significantly limits the discriminability of learned features. To resolve this issue, we propose MAP-VAE to enable the learning of global and local geometry by jointly leveraging global and local self-supervision. To enable effective local self-supervision, we introduce multi-angle analysis for point clouds. In a multi-angle scenario, we first split a point cloud into a front half and a back half from each angle, and then, train MAP-VAE to learn to predict a back half sequence from the corresponding front half sequence. MAP-VAE performs this half-to-half prediction using RNN to simultaneously learn each local geometry and the spatial relationship among them. In addition, MAP-VAE also learns global geometry via self-reconstruction, where we employ a variational constraint to facilitate novel shape generation. The outperforming results in four shape analysis tasks show that MAP-VAE can learn more discriminative global or local features than the state-of-the-art methods.


title: Generalizations of the Sommerfield and Schwinger models

abstract: The Sommerfield model with a massive vector field coupled to a massless fermion in 1+1 dimensions is an exactly solvable analog of a Bank-Zaks model. The "physics" of the model comprises a massive boson and an unparticle sector that survives at low energy as a conformal field theory (Thirring model). We analyze generalizations of the Sommerfield model, and the corresponding generalizations of the Schwinger model, with more massless fermions and more vector fields.


title: Model-Free Unsupervised Learning for Optimization Problems with Constraints

abstract: In many optimization problems in wireless communications, the expressions of objective function or constraints are hard or even impossible to derive, which makes the solutions difficult to find. In this paper, we propose a model-free learning framework to solve constrained optimization problems without the supervision of the optimal solution. Neural networks are used respectively for parameterizing the function to be optimized, parameterizing the Lagrange multiplier associated with instantaneous constraints, and approximating the unknown objective function or constraints. We provide learning algorithms to train all the neural networks simultaneously, and reveal the connections of the proposed framework with reinforcement learning. Numerical and simulation results validate the proposed framework and demonstrate the efficiency of model-free learning by taking power control problem as an example.


title: Spatial Modulation for Ambient Backscatter Communications: Modeling and Analysis

abstract: Multiple-antenna backscatter is emerging as a promising approach to offer high communication performance for the data-intensive applications of ambient backscatter communications (AmBC). Although much has been understood about multiple-antenna backscatter in conventional backscatter communications (CoBC), existing analytical models cannot be directly applied to AmBC due to the structural differences in RF source and tag circuit designs. This paper takes the first step to fill the gap, by exploring the use of spatial modulation (SM) in AmBC whenever tags are equipped with multiple antennas. Specifically, we present a practical multiple-antenna backscatter design for AmBC that exempts tags from the inter-antenna synchronization and mutual coupling problems while ensuring high spectral efficiency and ultra-low power consumption. We obtain an optimal detector for the joint detection of both backscatter signal and source signal based on the maximum likelihood principle. We also design a two-step algorithm to derive bounds on the bit error rate (BER) of both signals. Simulation results validate the analysis and show that the proposed scheme can significantly improve the throughput compared with traditional systems.


title: Millimeter-Wave NOMA with User Grouping, Power Allocation and Hybrid Beamforming

abstract: This paper investigates the application of non-orthogonal multiple access in millimeter-Wave communications (mmWave-NOMA). Particularly, we consider downlink transmission with a hybrid beamforming structure. A user grouping algorithm is first proposed according to the channel correlations of the users. Whereafter, a joint hybrid beamforming and power allocation problem is formulated to maximize the achievable sum rate, subject to a minimum rate constraint for each user. To solve this non-convex problem with high-dimensional variables, we first obtain the solution of power allocation under arbitrary fixed hybrid beamforming, which is divided into intra-group power allocation and inter-group power allocation. Then, given arbitrary fixed analog beamforming, we utilize the approximate zero-forcing method to design the digital beamforming to minimize the inter-group interference. Finally, the analog beamforming problem with the constant-modulus constraint is solved with a proposed boundary-compressed particle swarm optimization algorithm. Simulation results show that the proposed joint approach, including user grouping, hybrid beamforming and power allocation, outperforms the state-of-the-art schemes and the conventional mmWave orthogonal multiple access system in terms of achievable sum rate and energy efficiency.


title: Propensity score analysis with latent covariates: Measurement error bias correction using the covariate's posterior mean, aka the inclusive factor score

abstract: We address measurement error bias in propensity score (PS) analysis due to covariates that are latent variables. In the setting where latent covariate $X$ is measured via multiple error-prone items $\mathbf{W}$, PS analysis using several proxies for $X$ -- the $\mathbf{W}$ items themselves, a summary score (mean/sum of the items), or the conventional factor score (cFS , i.e., predicted value of $X$ based on the measurement model) -- often results in biased estimation of the causal effect, because balancing the proxy (between exposure conditions) does not balance $X$. We propose an improved proxy: the conditional mean of $X$ given the combination of $\mathbf{W}$, the observed covariates $Z$, and exposure $A$, denoted $X_{WZA}$. The theoretical support, which applies whether $X$ is latent or not (but is unobserved), is that balancing $X_{WZA}$ (e.g., via weighting or matching) implies balancing the mean of $X$. For a latent $X$, we estimate $X_{WZA}$ by the inclusive factor score (iFS) -- predicted value of $X$ from a structural equation model that captures the joint distribution of $(X,\mathbf{W},A)$ given $Z$. Simulation shows that PS analysis using the iFS substantially improves balance on the first five moments of $X$ and reduces bias in the estimated causal effect. Hence, within the proxy variables approach, we recommend this proxy over existing ones. We connect this proxy method to known results about weighting/matching functions (Lockwood & McCaffrey, 2016; McCaffrey, Lockwood, & Setodji, 2013). We illustrate the method in handling latent covariates when estimating the effect of out-of-school suspension on risk of later police arrests using Add Health data.


title: Depth of an initial ideal

abstract: Given an arbitrary integer $d>0$, we construct a homogeneous ideal $I$ of the polynomial ring $S = K[x_1, \ldots, x_{3d}]$ in $3d$ variables over a filed $K$ for which $S/I$ is a Cohen--Macaulay ring of dimension $d$ with the property that, for each of the integers $0 \leq r \leq d$, there exists a monomial order $<_r$ on $S$ with ${\rm depth} (S/{\rm in}_{<_r}(I)) = r$, where ${\rm in}_{<_r}(I)$ is the initial ideal of $I$ with respect to $<_r$.


title: A stochastic matching model on hypergraphs

abstract: Motivated by applications to a wide range of assemble-to-order systems, operations scheduling, healthcare systems and collaborative economy applications, we introduce a stochastic matching model on hypergraphs, extending the model in [15] to the case of hypergraphical (rather than graphical) matching structures. We address a discrete-event system under a random input of single items, simply using the system as an interface to be matched by groups of two or more. We study the stability of this stochastic system, for various hypergraph geometries.


title: On the delta function broadening in the Kubo-Greenwood equation

abstract: Understanding DC electrical conductivity is crucial for the study of materials. Macroscopic DC conductivity can be calculated from first principles using the Kubo-Greenwood equation. The procedure involves finding the thermodynamic limit of the current response to an electric field that is slowly switched on, and then taking the limit of the switching rate to zero. We introduce a nonlinear extrapolation procedure executed in systems with periodic boundary conditions, which predicts conductivity close to the thermodynamic limit even for very small systems. The scheme also overcomes a large part of the usual ambiguities of the DC conductivity definition for finite systems. We numerically compare our method to the Landauer equation-based approach and find both techniques to be consistent with each other.


title: Lecture Notes on Automata, Languages, and Grammars

abstract: These lecture notes are intended as a supplement to Moore and Mertens' The Nature of Computation or as a standalone resource, and are available to anyone who wants to use them. Comments are welcome, and please let me know if you use these notes in a course. There are 61 exercises. I emphasize that automata are elementary playgrounds where we can explore the issues of deterministic and nondeterministic computation. Unlike P vs. NP, we can prove that nondeterminism is equivalent to determinism, or strictly more powerful than determinism, in finite-state and push-down automata respectively. I also correct several historical and aesthetic injustices: in particular, the Myhill-Nerode theorem and the idea of building minimal DFAs from equivalence classes of prefixes is restored to its rightful place above the Pumping Lemma for regular languages. I also discuss the Pumping Lemma for context-free languages, and briefly discuss counter automata, queue automata, and the connection between unambiguous context-free languages and algebraic generating functions.


title: $N$-phonon bundle emission via the anti-Stokes process

abstract: We demonstrate theoretically the bundle emission of $n$ strongly correlated phonons in an acoustic cavity QED system. The mechanism relies on anti-Stokes resonances that generate super-Rabi oscillations between states with a large difference in their number of excitations, which, combined with dissipation, transfer coherently pure $n$-phonon states outside of the cavity. This process works with close to perfect purity over a wide range of parameters and is tunable optically with well-resolved operation conditions. This broadens the realm of quantum phononics, with potential applications for on-chip quantum information processing, quantum metrology, and engineering of new types of quantum devices, such as optically heralded $N$-phonon guns.


title: On Learning Hamiltonian Systems from Data

abstract: Concise, accurate descriptions of physical systems through their conserved quantities abound in the natural sciences. In data science, however, current research often focuses on regression problems, without routinely incorporating additional assumptions about the system that generated the data. Here, we propose to explore a particular type of underlying structure in the data: Hamiltonian systems, where an "energy" is conserved. Given a collection of observations of such a Hamiltonian system over time, we extract phase space coordinates and a Hamiltonian function of them that acts as the generator of the system dynamics. The approach employs an autoencoder neural network component to estimate the transformation from observations to the phase space of a Hamiltonian system. An additional neural network component is used to approximate the Hamiltonian function on this constructed space, and the two components are trained jointly. As an alternative approach, we also demonstrate the use of Gaussian processes for the estimation of such a Hamiltonian. After two illustrative examples, we extract an underlying phase space as well as the generating Hamiltonian from a collection of movies of a pendulum. The approach is fully data-driven, and does not assume a particular form of the Hamiltonian function.


title: Balanced Coherence Times of Mixed-Species Atomic Qubits in a Dual $3\times3$ Magic-Intensity Optical Dipole Trap Array

abstract: In this work, we construct a polarization-mediated magic-intensity (MI) optical dipole trap (ODT) array, in which the detrimental effects of light shifts on the mixed-species qubits are efficiently mitigated so that the coherence times of the mixed-species qubits are both substantially enhanced and balanced for the first time. This mixed-species magic trapping technique relies on the tunability of the coefficient of the third-order cross term and ground state hyperpolarizability, which are inherently dependent on the degree of circular polarization of the trap laser. Experimentally, polarization of the ODT array for $^{85}$Rb qubits is finely adjusted to a definite value so that its working magnetic field required for magic trapping amounts to the one required for magically trapping $^{87}$Rb qubits in another ODT array with fully circular polarization. Ultimately, in such a polarization-mediated MI-ODT array, the coherence times of $^{87}$Rb and $^{85}$Rb qubits are respectively enhanced up to 891$\pm$47 ms and 943$\pm$35 ms. Furthermore, a new source of dephasing effect is revealed, which arises from the noise of the elliptic polarization, and the reduction in corresponding dephasing effect on the $^{85}$Rb qubits is attainable by use of shallow magic intensity. It is anticipated that the novel mixed-species MI-ODT array is a versatile platform for building scalable quantum computers with neutral atoms.


title: Optimal Dynamic Multi-Resource Management in Earth Observation Oriented Space Information Networks

abstract: Space information network (SIN) is an innovative networking architecture to achieve near-real-time mass data observation, processing and transmission over the globe. In the SIN environment, it is essential to coordinate multi-dimensional heterogeneous resources (i.e., observation resource, computation resource and transmission resource) to improve network performance. However, the time varying property of both the observation resource and transmission resource is not fully exploited in existing studies. Dynamic resource management according to instantaneous channel conditions has a potential to enhance network performance. To this end, in this paper, we study the multi-resource dynamic management problem, considering stochastic observation and transmission channel conditions in SINs. Specifically, we develop an aggregate optimization framework for observation scheduling, compression ratio selection and transmission scheduling, and formulate a flow optimization problem based on extended time expanded graph (ETEG) to maximize the sum network utility. Then, we equivalently transform the flow optimization problem on ETEG as a queue stability-related stochastic optimization problem. An online algorithm is proposed to solve the problem in a slot-by-slot manner by exploiting the Lyapunov optimization technique. Performance analysis shows that the proposed algorithm achieves close-to-optimal network utility while guaranteeing bounded queue occupancy. Extensive simulation results further validate the efficiency of the proposed algorithm and evaluate the impacts of various network parameters on the algorithm performance.


title: ReadsClean: a new approach to error correction of sequencing reads based on alignments clustering

abstract: Motivation: Next generation methods of DNA sequencing produce relatively high rate of reading errors, which interfere with de novo genome assembly of newly sequenced organisms and particularly affect the quality of SNP detection important for diagnostics of many hereditary diseases. There exists a number of programs developed for correcting errors in NGS reads. Such programs utilize various approaches and are optimized for different specific tasks, but all of them are far from being able to correct all errors, especially in sequencing reads that crossing by repeats and DNA from di/polyploid eukaryotic genomes. Results: This paper describes a novel method of error correction based on clustering of alignments of similar reads. This method is implemented in ReadsClean program, which is designed for cleaning Illumina HiSeq sequencing reads. We compared ReadsClean to other reads cleaning programs recognized to be the best by several publications. Our sequence assembly tests using actual and simulated sequencing reads show superior results achieved by ReadsClean. Availability and implementation: ReadsClean is implemented as a standalone C code. It is incorporated in an error correction pipeline and is freely available to academic users at Softberry web server www.softberry.com.


title: Triangle singularities in ${J/\psi\rightarrow\eta\pi^0\phi}$ and ${\pi^0\pi^0\phi}$

abstract: The BESIII Collaboration recently reported the observation of the $a_0(980)^0-f_0(980)$ mixing in the isospin breaking decay $J/\psi\to \eta\pi^0\phi$. In the Dalitz plot for that decay with the $\eta$ reconstructed from two photons, there is a band around $1.4$~GeV on the $\pi^0 \phi$ distribution. In general, this peak can be due to a resonance or a kinematic effect. In this paper, we study the effects of a set of $K^*K\bar K$ triangle diagrams, and show that due to triangle singularities such diagrams can lead to a peak around 1.4~GeV in the $\pi^0\phi$ invariant mass distribution. The Dalitz plot induced by such a mechanism has a feature consistent with the BESIII observation, namely events along the band accumulate at both ends close to the Dalitz plot boundary. The effect of the same mechanism on the $J/\psi\to \pi^0\pi^0\phi$ and $J/\psi \to \eta\pi^0 K^+K^-$ decays are also investigated. We suggest to take more data for the $J/\psi\to \eta\pi^0\phi\to \eta\pi^0 K^+ K^-$ and check whether the structure around $1.4$~GeV persists for the $K^+K^-$ invariant mass away from the $\phi$ mass region. This is crucial for understanding whether the band is due to triangle singularities or due to a resonance. Were it the latter, the band should remain while it would not if it is due to the former.


title: Exploring large scale public medical image datasets

abstract: Rationale and Objectives: Medical artificial intelligence systems are dependent on well characterised large scale datasets. Recently released public datasets have been of great interest to the field, but pose specific challenges due to the disconnect they cause between data generation and data usage, potentially limiting the utility of these datasets. Materials and Methods: We visually explore two large public datasets, to determine how accurate the provided labels are and whether other subtle problems exist. The ChestXray14 dataset contains 112,120 frontal chest films, and the MURA dataset contains 40,561 upper limb radiographs. A subset of around 700 images from both datasets was reviewed by a board-certified radiologist, and the quality of the original labels was determined. Results: The ChestXray14 labels did not accurately reflect the visual content of the images, with positive predictive values mostly between 10% and 30% lower than the values presented in the original documentation. There were other significant problems, with examples of hidden stratification and label disambiguation failure. The MURA labels were more accurate, but the original normal/abnormal labels were inaccurate for the subset of cases with degenerative joint disease, with a sensitivity of 60% and a specificity of 82%. Conclusion: Visual inspection of images is a necessary component of understanding large image datasets. We recommend that teams producing public datasets should perform this important quality control procedure and include a thorough description of their findings, along with an explanation of the data generating procedures and labelling rules, in the documentation for their datasets.


title: An Unusual Mid-Infrared Flare in a Type 2 AGN: An Obscured Turning-on AGN or Tidal Disruption Event?

abstract: We report the discovery of an exceptional MIR flare in a Type 2 AGN, SDSS J165726.81+234528.1, at $z=0.059$. This object brightened by 3 mag in the Wide-field Infrared Survey Explorer (WISE) $W1$ and $W2$ bands between 2015 and 2017 (and is fading since 2018), without significant changes ($\lesssim$ 0.2 mag) in the optical over the same period of time. Based on the WISE light curves and near-IR imaging, the flare is more significant at longer wavelengths, suggesting an origin of hot dust emission. The estimated black hole mass ($\sim 10^{6.5}\,M_\odot$) from different methods places its peak bolometric luminosity around the Eddington limit. The high luminosity of the MIR flare and its multi-year timescale suggest that it most likely originated from reprocessed dust radiation in an extended torus surrounding the AGN, instead of from stellar explosions. The MIR color variability is consistent with known changing-look AGN and tidal disruption events (TDEs), but inconsistent with normal supernovae. We suggest that it is a turning-on Type 2 AGN or TDE, where the optical variability is obscured by the dust torus during the transition. This MIR flare event reveals a population of dramatic nuclear transients that are missed in the optical.


title: Revealing Energy Dependence of Quantum Defects via Two Heteronuclear Atoms in an Optical Tweezer

abstract: As a physically motivated and computationally simple model for cold atomic and molecular collisions, the multichannel quantum defect theory (MQDT) with frame transformation (FT) formalism provides an analytical treatment of scattering resonances in an arbitrary partial wave between alkali-metal atoms, leading to the experimental observation of $p-$ and $d-$wave resonances. However, the inconsistency of quantum defects for describing scattering resonances shows up when compared with experiments. Here, with two heteronuclear atoms in the ground state of an optical tweezer, the energy dependence of quantum defects is obviously revealed by comparing the measured s-wave scattering length with the prediction of MQDT-FT. By dividing the quantum defects into energy sensitive and insensitive categories, the inconsistency is ultimately removed while retaining the analytic structure of MQDT-FT. This study represents a significant improvement in the analytical MQDT-FT and demonstrates that a clean two-particle system is valuable to the test of collisional physics.


title: Euclidean Forward-Reverse Brascamp-Lieb Inequalities: Finiteness, Structure and Extremals

abstract: A new proof is given for the fact that centered gaussian functions saturate the Euclidean forward-reverse Brascamp-Lieb inequalities, extending the Brascamp-Lieb and Barthe theorems. A duality principle for best constants is also developed, which generalizes the fact that the best constants in the Brascamp-Lieb and Barthe inequalities are equal. Finally, as the title hints, the main results concerning finiteness, structure and gaussian-extremizability for the Brascamp-Lieb inequality due to Bennett, Carbery, Christ and Tao are generalized to the setting of the forward-reverse Brascamp-Lieb inequality.


title: Classical and Quantum Algorithms for Tensor Principal Component Analysis

abstract: We present classical and quantum algorithms based on spectral methods for a problem in tensor principal component analysis. The quantum algorithm achieves a quartic speedup while using exponentially smaller space than the fastest classical spectral algorithm, and a super-polynomial speedup over classical algorithms that use only polynomial space. The classical algorithms that we present are related to, but slightly different from those presented recently in Ref. 1. In particular, we have an improved threshold for recovery and the algorithms we present work for both even and odd order tensors. These results suggest that large-scale inference problems are a promising future application for quantum computers.


title: Steady-State Simulation for Combined Transmission and Distribution Systems

abstract: The future electric grid will consist of significant penetration of renewable and distributed generation that is likely to create a homogenous transmission and distribution (T&D) system, requiring tools that can model and robustly simulate the combined T&D networks. Existing tools use disparate models and formulations for simulation of transmission versus distribution grids and solving for the steady-state solution of the combined T&D networks often lacks convergence robustness and scalability to large systems. In this paper, we show that modeling both the T&D grid elements in terms of currents and voltages using an equivalent circuit framework enables simulation of combined positive sequence networks of the transmission grids with three-phase networks of the distribution grids without loss of generality. We further demonstrate that we can ensure robust convergence for these resulting large-scale complex T&D systems when the circuit simulation methods are applied to them. Our results illustrate robust convergence of combined T&D networks using a direct Newton-Raphson solver on a single machine for smaller sized systems and using a parallel Gauss-Seidel-Newton solver on multiple machines for larger sized systems with greater than million nodes.


title: Relative phantom maps and rational homotopy

abstract: We generalize some results of Gray and McGibbon-Roitberg on relations between phantom maps and rational homotopy to relative phantom maps. Since the $\lim^1$ and the profinite completion techniques do not apply to relative phantom maps, we develop new techniques.


title: Confounder-Aware Visualization of ConvNets

abstract: With recent advances in deep learning, neuroimaging studies increasingly rely on convolutional networks (ConvNets) to predict diagnosis based on MR images. To gain a better understanding of how a disease impacts the brain, the studies visualize the salience maps of the ConvNet highlighting voxels within the brain majorly contributing to the prediction. However, these salience maps are generally confounded, i.e., some salient regions are more predictive of confounding variables (such as age) than the diagnosis. To avoid such misinterpretation, we propose in this paper an approach that aims to visualize confounder-free saliency maps that only highlight voxels predictive of the diagnosis. The approach incorporates univariate statistical tests to identify confounding effects within the intermediate features learned by ConvNet. The influence from the subset of confounded features is then removed by a novel partial back-propagation procedure. We use this two-step approach to visualize confounder-free saliency maps extracted from synthetic and two real datasets. These experiments reveal the potential of our visualization in producing unbiased model-interpretation.


title: Is There Any Recovery Guarantee with Coupled Structured Matrix Factorization for Hyperspectral Super-Resolution?

abstract: Coupled structured matrix factorization (CoSMF) for hyperspectral super-resolution (HSR) has recently drawn significant interest in hyperspectral imaging for remote sensing. Presently there is very few work that studies the theoretical recovery guarantees of CoSMF. This paper makes one such endeavor by considering the CoSMF formulation by Wei et al., which, simply speaking, is similar to coupled non-negative matrix factorization. Assuming no noise, we show sufficient conditions under which the globably optimal solution to the CoSMF problem is guaranteed to deliver certain recovery accuracies. Our analysis suggests that sparsity and the pure-pixel (or separability) condition play a hidden role in enabling CoSMF to achieve some good recovery characteristics.


title: MPGD simulation in negative-ion gas for direction-sensitive dark matter searches

abstract: Negative-ion time projection chambers (NITPCs) provide the absolute z-position from the arrival-times' difference of several species of negative-ions. This is a strong background reduction method for rare event searches, such as direct dark mater searches. The characteristics of micro pattern gaseous detectors (MPGDs) in negative-ion gas were studied by both experiments and simulations to develop an NITPC for dark matter search experiments. An MPGD simulation method for the negative-ion gas with Garfield++, which has not yet been established, was developed. We present herein the MPGD simulation results in negative-ion gas, $\mathrm{SF_{6}}$ gas, together with the experimental results, such as gain, energy resolution, and signal creation from the MPGD for future NITPC optimization.