title: Generalising the scattered property of subspaces

abstract: Let $V$ be an $r$-dimensional $\mathbb{F}_{q^n}$-vector space. We call an $\mathbb{F}_q$-subspace $U$ of $V$ $h$-scattered if $U$ meets the $h$-dimensional $\mathbb{F}_{q^n}$-subspaces of $V$ in $\mathbb{F}_q$-subspaces of dimension at most $h$. In 2000 Blokhuis and Lavrauw proved that $\dim_{\mathbb{F}_q} U \leq rn/2$ when $U$ is $1$-scattered. Subspaces attaining this bound have been investigated intensively because of their relations with projective two-weight codes and strongly regular graphs. MRD-codes with a maximum idealiser have also been linked to $rn/2$-dimensional $1$-scattered subspaces and to $n$-dimensional $(r-1)$-scattered subspaces. In this paper we prove the upper bound $rn/(h+1)$ for the dimension of $h$-scattered subspaces, $h>1$, and construct examples with this dimension. We study their intersection numbers with hyperplanes, introduce a duality relation among them, and study the equivalence problem of the corresponding linear sets.


title: Spatial 3D Mat\'ern priors for fast whole-brain fMRI analysis

abstract: Bayesian whole-brain functional magnetic resonance imaging (fMRI) analysis with three-dimensional spatial smoothing priors has been shown to produce state-of-the-art activity maps without pre-smoothing the data. The proposed inference algorithms are computationally demanding however, and the proposed spatial priors have several less appealing properties, such as being improper and having infinite spatial range. We propose a statistical inference framework for whole-brain fMRI analysis based on the class of Mat\'ern covariance functions. The framework uses the Gaussian Markov random field (GMRF) representation of possibly anisotropic spatial Mat\'ern fields via the stochastic partial differential equation (SPDE) approach of Lindgren et al. (2011). This allows for more flexible and interpretable spatial priors, while maintaining the sparsity required for fast inference in the high-dimensional whole-brain setting. We develop an accelerated stochastic gradient descent (SGD) optimization algorithm for empirical Bayes (EB) inference of the spatial hyperparameters. Conditionally on the inferred hyperparameters, we make a fully Bayesian treatment of the brain activity. The Mat\'ern prior is applied to both simulated and experimental task-fMRI data and clearly demonstrates that it is a more reasonable choice than the previously used priors, using comparisons of activity maps, prior simulation and cross-validation.


title: Tactile Hallucinations on Artificial Skin Induced by Homeostasis in a Deep Boltzmann Machine

abstract: Perceptual hallucinations are present in neurological and psychiatric disorders and amputees. While the hallucinations can be drug-induced, it has been described that they can even be provoked in healthy subjects. Understanding their manifestation could thus unveil how the brain processes sensory information and might evidence the generative nature of perception. In this work, we investigate the generation of tactile hallucinations on biologically inspired, artificial skin. To model tactile hallucinations, we apply homeostasis, a change in the excitability of neurons during sensory deprivation, in a Deep Boltzmann Machine (DBM). We find that homeostasis prompts hallucinations of previously learned patterns on the artificial skin in the absence of sensory input. Moreover, we show that homeostasis is capable of inducing the formation of meaningful latent representations in a DBM and that it significantly increases the quality of the reconstruction of these latent states. Through this, our work provides a possible explanation for the nature of tactile hallucinations and highlights homeostatic processes as a potential underlying mechanism.


title: Surface Polariton-Like s-Polarized Waveguide Modes in Switchable Dielectric Thin-Films on Polar Crystals

abstract: Surface phonon polaritons (SPhP) and surface plasmon polaritons (SPP), evanescent modes supported by media with negative permittivity, are a fundamental building block of nanophotonics. These modes are unmatched in terms of field enhancement and spatial confinement, and dynamical all-optical control can be achieved e.g. by employing phase-change materials (PCMs). However, the excitation of surface polaritons in planar structures is intrinsically limited to p-polarization. On the contrary, waveguide modes in high-permittivity films can couple to both p- and s-polarized light, and in thin films, their confinement can become comparable to surface polaritons. Here we demonstrate that the s-polarized waveguide mode in a thin Ge$_3$Sb$_2$Te$_6$ (GST) film features a similar dispersion, confinement, and electric field enhancement as the SPhP mode of the silicon carbide (SiC) substrate, while even expanding the allowed frequency range. Moreover, we experimentally show that switching the GST film grants non-volatile control over the SPhP and the waveguide mode dispersions. We provide an analytical model for the description of the GST/SiC waveguide mode and show that our concept is applicable to the broad variety of polar crystals throughout the infrared spectral range. As such, complementarily to the polarization-limited surface polaritons, the s-polarized PCM waveguide mode constitutes a promising additional building block for nanophotonic applications.


title: Towards the Baikal Open Laboratory in Astroparticle Physics

abstract: The open science framework defined in the German-Russian Astroparticle Data Life Cycle Initiative (GRADLCI) has triggered educational and outreach activities at the Irkutsk State University (ISU), which is actively participated in the two major astroparticle facilities in the region: TAIGA observatory and Baikal-GVD neutrino telescope. We describe the ideas grew out of this unique environment and propose a new open science laboratory based on education and outreach as well as on the development and testing new methods and techniques for the multimessenger astronomy.


title: Spectro-timing analysis of MAXI J1535-571 using AstroSat

abstract: We report the results of the analysis of an AstroSat observation of the Black Hole candidate MAXI J1535-571 during its Hard Intermediate state. We studied the evolution of the spectral and timing parameters of the source during the observation. The observation covered a period of $\sim$5 days and consisted of 66 continuous segments, corresponding to individual spacecraft orbits. Each segment was analysed independently. The source count rate increased roughly linearly by $\sim$30 %. We modelled the spectra as a combination of radiation from a thermal disk component and a power-law. The timing analysis revealed the presence of strong Quasi Periodic Oscillations with centroid frequency $\nu_{\rm{QPO}}$ fluctuating in the range 1.7-3.0 Hz. We found a tight correlation between the QPO centroid frequency $\nu_{\rm{QPO}}$ and the power-law spectral index $\Gamma$, while $\nu_{\rm{QPO}}$ appeared not to be correlated with the linearly-increasing flux itself. We discuss the implications of these results on physical models of accretion.


title: A short note on Multilevel Toeplitz Matrices

abstract: Chien, Liu, Nakazato, and Tam proved that all n by n classical Toeplitz matrices (one-level Toeplitz matrices) are unitarily similar to complex symmetric matrices via two types of unitary matrices and the type of the unitary matrices only depends on the parity of n. In this paper, we extend their result to multilevel Toeplitz matrices that any multilevel Toeplitz matrix is unitarily similar to a complex symmetric matrix. We provide a method to construct the unitary matrices that uniformly turn any multilevel Toeplitz matrix to a complex symmetric matrix by taking tensor products of these two types of unitary matrices for one-level Toeplitz matrices according to the parity of each level of the multilevel Toeplitz matrices. In addition, we introduce a class of complex symmetric matrices that are unitarily similar to some p-level Toeplitz matrices.


title: Bandgap-Assisted Quantum Control of Topological Edge States in a Cavity

abstract: Quantum matter with exotic topological order has potential applications in quantum computation. However, in present experiments, the manipulations on topological states are still challenging. We here propose an architecture for optical control of topological matter. We consider a topological superconducting qubit array with Su-Schrieffer-Heeger (SSH) Hamiltonian which couples to a microwave cavity. Based on parity properties of the topological qubit array, we propose an optical spectroscopy method to observe topological phase transition, i.e., edge-to-bulk transition. This new method can be achieved by designing cavity-qubit couplings. A main purpose of this work is to understand how topological phase transition affects light-matter interaction. We find that topological bandgap plays an essential role on this issue. In topological phase, the resonant vacuum Rabi splitting of degenerate edge states coupling to the cavity field is protected from those of bulk states by the bandgap. In dispersive regime, the cavity induced coupling between edge states is dominant over couplings between edge and bulk states, due to the topological bandgap. As a result, quantum interference between topological edge states occures and enables single-photon transport through boundaries of the topological qubit array. Our work may pave a way for topological quantum state engineering.


title: Phase tracking based on GPGPU and applications in Planetary radio Science

abstract: This paper introduces a phase tracking method for planetary radio science research with computational algorithm implemented fo r NVIDIA GPUs. In contrast to the phase-locked loop (PPL) phase counting method used in traditional Doppler data processing, this method fits the tracking data signal into the shape expressed by the Taylor polynomial with optimal phase and amplitude coefficients. The Differential Evolution (DE) algorithm is employed for polynomial fitting. In order to cope with high computational intensity of the proposed phase tracking method, the graphics processing units (GPUs) are employed. As a result, the method estimates the instantaneous phase, frequency, derivative of frequency (line-of-sight acceleration) and the total count phase of different integration scales. This data can be further used in planetary radio science research to analyze the planetary occultation and gravitational fields. The method has been tested on MEX (Mars Express, ESA) and Chang'E 4 relay satellite (China) tracking data. In a real experiment with 400K data block size and $\sim$80,000 DE solver objective function evaluations we were able to acheive the target convergence threshold in 6.5 seconds and do real-time processing on NVIDIA GTX580 and 2$\times$ NVIDIA K80 GPUs, respectively. The precision of integral Doppler (60s) is 2 mrad/s and 4 mrad/s for MEX(3-way) and Chang'E 4 relay satellite(3-way) respectively.


title: Vanishing viscosity limit to vortex sheet for the isentropic compressible circularly symmetric 2D flow

abstract: In this paper, we consider the small viscosity limit problem for the isentropic compressible Navier-Stokes equations in a 2D exterior domain with impermeable boundary conditions , and the corresponding Euler equations have vortex sheet solutions.We obtain that away from the boundary and the contact discontinuous the isentropic compressible viscous flow can be approximated by the corresponding inviscid flow, near the boundary (the contact discontinuous) there is a boundary layer (vortex layer)for the angular velocity in the leading order expansion of solution, while the radial velocity and the pressure do not have boundary layers (vortex layers) in the leading order. We rigorously justify the asymptotic behavior of solutions in the $L^{\infty}$ space for the small viscosities limit in the Lagrangian coordinates.


title: Even-to-Odd Barium Isotope Ratio in Selected Galactic Halo Stars

abstract: We have updated the Ba II model atom by taking into account the H I impact excitation with the rate coefficients from the quantum-mechanical calculations of Belyaev and Yakovleva (2018). Using high-resolution stellar spectra and the non-local thermodynamic equilibrium (non-LTE) line formation for Ba II, we have determined the fraction of barium isotopes with an odd mass number (fodd}) in four Galactic halo giants with well-known atmospheric parameters. We use a method based on the requirement that the abundances from the Ba II 4554 A resonance and Ba II 5853, 6496 A subordinate lines be equal. An accuracy of 0.04 dex in determining the barium abundance from individual lines has been achieved. In three stars (HD 2796, HD 108317, and HD 122563) fodd >= 0.4. This suggests that >= 80 % of the barium observed in these stars was synthesized in the r-process. In HD 128279 fodd = 0.27 exceeds the fraction of odd barium isotopes in the Solar system, but only slightly. The dominance of the r-process at the formation epoch of our sample stars is confirmed by the presence of a europium overabundance relative barium, with [Eu/Ba] > 0.3. We have calculated the non-LTE abundance corrections for five Ba II lines and investigated their dependence on atmospheric parameters in the ranges of effective temperatures from 4500 to 6500 K, surface gravities log g from 0.5 to 4.5, and metallicities [Fe/H] from 0 to -3.


title: Signatures of Topological Phonons in Superisostatic Lattices

abstract: Soft topological surface phonons in idealized ball-and-spring lattices with coordination number $z=2d$ in $d$ dimensions become finite-frequency surface phonons in physically realizable superisostatic lattices with $z>2d$. We study these finite-frequency modes in model lattices with added next-nearest-neighbor springs or bending forces at nodes with an eye to signatures of the topological surface modes that are retained in the physical lattices. Our results apply to metamaterial lattices, prepared with modern printing techniques, that closely approach isostaticity.


title: Pyramid: A General Framework for Distributed Similarity Search

abstract: Similarity search is a core component in various applications such as image matching, product recommendation and low-shot classification. However, single machine solutions are usually insufficient due to the large cardinality of modern datasets and stringent latency requirement of on-line query processing. We present Pyramid, a general and efficient framework for distributed similarity search. Pyramid supports search with popular similarity functions including Euclidean distance, angular distance and inner product. Different from existing distributed solutions that are based on KD-tree or locality sensitive hashing (LSH), Pyramid is based on Hierarchical Navigable Small World graph (HNSW), which is the state of the art similarity search algorithm on a single machine. To achieve high query processing throughput, Pyramid partitions a dataset into sub-datasets containing similar items for index building and assigns a query to only some of the sub-datasets for query processing. To provide the robustness required by production deployment, Pyramid also supports failure recovery and straggler mitigation. Pyramid offers a set of concise API such that users can easily use Pyramid without knowing the details of distributed execution. Experiments on large-scale datasets show that Pyramid produces quality results for similarity search, achieves high query processing throughput and is robust under node failure and straggler.


title: Total variation vs L1 regularization: a comparison of compressive sensing optimization methods for chemical detection

abstract: One of the fundamental assumptions of compressive sensing (CS) is that a signal can be reconstructed from a small number of samples by solving an optimization problem with the appropriate regularization term. Two standard regularization terms are the L1 norm and the total variation (TV) norm. We present a comparison of CS reconstruction results based on these two approaches in the context of chemical detection, and we demonstrate that optimization based on the L1 norm outperforms optimization based on the TV norm. Our comparison is driven by CS sampling, reconstruction, and chemical detection in two real-world datasets: the Physical Sciences Inc. Fabry-P\'{e}rot interferometer sensor multispectral dataset and the Johns Hopkins Applied Physics Lab FTIR-based longwave infrared sensor hyperspectral dataset. Both datasets contain the release of a chemical simulant such as glacial acetic acid, triethyl phosphate, and sulfur hexafluoride. For chemical detection we use the adaptive coherence estimator (ACE) and bulk coherence, and we propose algorithmic ACE thresholds to define the presence or absence of a chemical of interest in both un-compressed data cubes and reconstructed data cubes. The un-compressed data cubes provide an approximate ground truth. We demonstrate that optimization based on either the L1 norm or TV norm results in successful chemical detection at a compression rate of 90%, but we show that L1 optimization is preferable. We present quantitative comparisons of chemical detection on reconstructions from the two methods, with an emphasis on the number of pixels with an ACE value above the threshold.


title: Introducing a response-based duration metric and its correlation with structural damages

abstract: This study proposes a response-based parameter for strong motion duration which is computed for structures and is the total time they are nonlinear during an earthquake. Correlation between structural response and duration for structures, subjected to a set of spectrum matched ground motions, is employed to examine the efficiency of the proposed method. The spectral matching procedure ensures that the influence of amplitude and frequency content of motions on structural response variability is significantly removed. Four concrete building type systems are studied and correlation coefficients of structural response with the proposed duration definition are examined. Comparison of the proposed method with other existing definitions, the record-based and response-based metrics, shows about 15-20% improvement in the correlation values.


title: van der Waals forces stabilize low-energy polymorphism in B2O3: Implications for the crystallization anomaly

abstract: The cohesive energies and structural properties of recently predicted, and never synthesized, B$_2$O$_3$ polymorphs are investigated from first principles using density functional theory and high-accuracy many-body methods, namely, the random phase approximation and quantum Monte Carlo. We demonstrate that the van der Waals forces play a key role in making the experimentally known polymorph (B$_2$O$_3$-I) the lowest in energy, with many competing metastable structures lying only a few kcal/mol above. Remarkably, all metastable crystals are comparable in energy and density to the glass, while having anisotropic and mechanically soft structures. Furthermore, the best metastable polymorph according to our stability criteria has a structural motif found in both the glass and a recently synthesized borosulfate compound. Our findings provide new perspectives for understanding the B$_2$O$_3$ anomalous behavior, namely, its propensity to vitrify in a glassy structure drastically different from the known crystal.


title: DALI: a large Dataset of synchronized Audio, LyrIcs and notes, automatically created using teacher-student machine learning paradigm

abstract: The goal of this paper is twofold. First, we introduce DALI, a large and rich multimodal dataset containing 5358 audio tracks with their time-aligned vocal melody notes and lyrics at four levels of granularity. The second goal is to explain our methodology where dataset creation and learning models interact using a teacher-student machine learning paradigm that benefits each other. We start with a set of manual annotations of draft time-aligned lyrics and notes made by non-expert users of Karaoke games. This set comes without audio. Therefore, we need to find the corresponding audio and adapt the annotations to it. To that end, we retrieve audio candidates from the Web. Each candidate is then turned into a singing-voice probability over time using a teacher, a deep convolutional neural network singing-voice detection system (SVD), trained on cleaned data. Comparing the time-aligned lyrics and the singing-voice probability, we detect matches and update the time-alignment lyrics accordingly. From this, we obtain new audio sets. They are then used to train new SVD students used to perform again the above comparison. The process could be repeated iteratively. We show that this allows to progressively improve the performances of our SVD and get better audio-matching and alignment.


title: Newswire versus Social Media for Disaster Response and Recovery

abstract: In a disaster situation, first responders need to quickly acquire situational awareness and prioritize response based on the need, resources available and impact. Can they do this based on digital media such as Twitter alone, or newswire alone, or some combination of the two? We examine this question in the context of the 2015 Nepal Earthquakes. Because newswire articles are longer, effective summaries can be helpful in saving time yet giving key content. We evaluate the effectiveness of several unsupervised summarization techniques in capturing key content. We propose a method to link tweets written by the public and newswire articles, so that we can compare their key characteristics: timeliness, whether tweets appear earlier than their corresponding news articles, and content. A novel idea is to view relevant tweets as a summary of the matching news article and evaluate these summaries. Whenever possible, we present both quantitative and qualitative evaluations. One of our main findings is that tweets and newswire articles provide complementary perspectives that form a holistic view of the disaster situation.


title: Probing Dark-ALP Portals at Future $e^+e^-$ Colliders

abstract: We study portal interactions connecting visible and dark sectors, and involving local interactions of a photon, a dark photon and a axion-like particle (ALP) at future $e^+e^-$ colliders. These interactions, mediated by higher-dimensional effective operators, may arise at one-loop by kinetic mixing between dark and ordinary photons, or, for massless dark photons, by direct short-distance contributions. We explore these portal interactions for a heavy ALP with masses between about 10 GeV and 230 GeV by investigating the sensitivity of the production $e^+e^- \to \gamma \gamma \bar{\gamma}$ to the effective couplings, where the dark photon $\bar\gamma$ gives rise to missing momentum in the final state. We will show how an appropriate choice of missing-energy and missing-mass cuts can optimize the signal to standard-model background ratio. Exclusion regions for the effective photon-dark-photon-ALP couplings versus the ALP mass are worked out for a few representative values of the collision energy and integrated luminosity, as presently envisaged by future $e^+e^-$ projects.


title: On operators with closed range and semi-Fredholm operators over W*-algebras

abstract: In this paper we consider A-Fredholm and semi-A-Fredholm operators on Hilbert C*-modules over a W*-algebra A defined in [3],[10]. Using the assumption that A is a W*-algebra (and not an arbitrary C*-algebra), we obtain several results such as generalization of Schechter-Lebow characterization of semi-Fredholm operators and generalization of "punctured neighbourhood" theorem, as well as some other results that generalize their classical counterparts. We consider both adjointable and non adjointable semi-Fredholm operators over W*-algebras. Moreover, we also work with general bounded, adjointable operators with closed range over C*-algebras and prove a generalization to Hilbert C*-modules of the result in [1] on Hilbert spaces.


title: Non-Collapsible Dual Complexes and Fake del Pezzo Surfaces

abstract: We propose the new construction of complex surfaces with $h^{1,0} = h^{2,0} = 0$ from smoothings of normal crossing surfaces with non-collapsible dual complexes and carry it out for the simplest case of the duncehat complex, obtaining the surface with $h^{1,1} = 9$ (presumably Barlow surface).


title: (Pseudo) Random Quantum States with Binary Phase

abstract: We prove a quantum information-theoretic conjecture due to Ji, Liu and Song (CRYPTO 2018) which suggested that a uniform superposition with random \emph{binary} phase is statistically indistinguishable from a Haar random state. That is, any polynomial number of copies of the aforementioned state is within exponentially small trace distance from the same number of copies of a Haar random state. As a consequence, we get a provable elementary construction of \emph{pseudorandom} quantum states from post-quantum pseudorandom functions. Generating pseduorandom quantum states is desirable for physical applications as well as for computational tasks such as quantum money. We observe that replacing the pseudorandom function with a $(2t)$-wise independent function (either in our construction or in previous work), results in an explicit construction for \emph{quantum state $t$-designs} for all $t$. In fact, we show that the circuit complexity (in terms of both circuit size and depth) of constructing $t$-designs is bounded by that of $(2t)$-wise independent functions. Explicitly, while in prior literature $t$-designs required linear depth (for $t > 2$), this observation shows that polylogarithmic depth suffices for all $t$. We note that our constructions yield pseudorandom states and state designs with only real-valued amplitudes, which was not previously known. Furthermore, generating these states require quantum circuit of restricted form: applying one layer of Hadamard gates, followed by a sequence of Toffoli gates. This structure may be useful for efficiency and simplicity of implementation.


title: A unified sparse optimization framework to learn parsimonious physics-informed models from data

abstract: Machine learning (ML) is redefining what is possible in data-intensive fields of science and engineering. However, applying ML to problems in the physical sciences comes with a unique set of challenges: scientists want physically interpretable models that can (i) generalize to predict previously unobserved behaviors, (ii) provide effective forecasting predictions (extrapolation), and (iii) be certifiable. Autonomous systems will necessarily interact with changing and uncertain environments, motivating the need for models that can accurately extrapolate based on physical principles (e.g. Newton's universal second law for classical mechanics, $F=ma$). Standard ML approaches have shown impressive performance for predicting dynamics in an interpolatory regime, but the resulting models often lack interpretability and fail to generalize. We introduce a unified sparse optimization framework that learns governing dynamical systems models from data, selecting relevant terms in the dynamics from a library of possible functions. The resulting models are parsimonious, have physical interpretations, and can generalize to new parameter regimes. Our framework allows the use of non-convex sparsity promoting regularization functions and can be adapted to address key challenges in scientific problems and data sets, including outliers, parametric dependencies, and physical constraints. We show that the approach discovers parsimonious dynamical models on several example systems. This flexible approach can be tailored to the unique challenges associated with a wide range of applications and data sets, providing a powerful ML-based framework for learning governing models for physical systems from data.


title: Advances in Implementation, Theoretical Motivation, and Numerical Results for the Nested Iteration with Range Decomposition Algorithm

abstract: This paper studies a low-communication algorithm for solving elliptic partial differential equations (PDE's) on high-performance machines, the nested iteration with range decomposition algorithm (NIRD). Previous work has shown that NIRD converges to a high level of accuracy within a small, fixed number of iterations (usually one or two) when applied to simple elliptic problems. This paper makes some improvements to the NIRD algorithm (including the addition of adaptivity during preprocessing, wider choice of partitioning functions, and modified error measurement) that enhance the method's accuracy and scalability, especially on more difficult problems. In addition, an updated convergence proof is presented based on heuristic assumptions that are supported by numerical evidence. Furthermore, a new performance model is developed that shows increased performance benefits for NIRD when problems are more expensive to solve using traditional methods. Finally, extensive testing on a variety of elliptic problems provides additional insight into the behavior of NIRD and additional evidence that NIRD achieves excellent convergence on a wide class of elliptic PDE's and, as such, should be a very competitive method for solving PDE's on large parallel computers.


title: On the spectral radii of the unicyclic hypergraphs with fixed matching number

abstract: We determine the unique hypergraphs with maximum spectral radius among all connected $k$-uniform ($k\geq 3$) unicyclic hypergraphs with matching number at least $z$, and among all connected $k$-uniform ($k\geq 3$) unicyclic hypergraphs with a given matching number, respectively.


title: Krivine diffusions attain the Goemans--Williamson approximation ratio

abstract: Answering a question of Abbasi-Zadeh, Bansal, Guruganesh, Nikolov, Schwartz and Singh (2018), we prove the existence of a slowed-down sticky Brownian motion whose induced rounding for MAXCUT attains the Goemans--Williamson approximation ratio. This is an especially simple particular case of the general rounding framework of Krivine diffusions that we investigate elsewhere.


title: Quantization of Poisson Hopf algebras

abstract: We describe a method for quantization of Poisson Hopf algebras in $\mathbb Q$-linear symmetric monoidal categories. It is compatible with tensor products and can also be used to produce braided Hopf algebras. The main idea comes from the fact that nerves of groups are symmetric simplicial sets. Nerves of Hopf algebras then turn out to be braided rather than symmetric and nerves of Poisson Hopf algebras to be infinitesimally braided. The problem is thus solved via the standard machinery of Drinfeld associators.


title: Galactic Tide and Local Stellar Perturbations on the Oort Cloud: Creation of Interstellar Comets

abstract: Comets in the Oort cloud evolve under the influence of internal and external perturbations, such as giant planets, stellar passages, and the galactic tidal field. We aim to study the dynamical evolution of the comets in the Oort cloud, accounting for external perturbations (passing stars and the galactic tide). We first construct an analytical model of stellar encounters. We find that individual perturbations do not modify the dynamics of the comets in the cloud unless very close (< 0.5pc) encounters occur. Using proper motions, parallaxes, and radial velocities from Gaia DR2, we construct an astrometric catalogue of 14,659 stars that are within 50pc from the Sun. For all these stars we calculate the time and the closest distance to the Sun. We find that the cumulative effect of relatively distant ($\leq1$ pc) passing stars can perturb the comets in the Oort cloud. Finally, we study the dynamical evolution of the comets in the Oort cloud under the influence of multiple stellar encounters within 2.5pc from the Sun and the galactic tidal field over $\pm10$Myr. We considered two models for the Oort cloud, compact (a $\leq$0.25 pc) and extended (a$ \leq0.5$ pc). We find that the cumulative effect of stellar encounters is the major perturber of the Oort cloud for a compact configuration while for the extended, the galactic tide is the major perturber. In both cases, the effect of passing stars and the galactic tide raises the semi-major axis of $\sim1.1$\% of the comets at the edge of the cloud up to interstellar regions ($a >0.5$pc). This leads to the creation of transitional interstellar comets, which might become interstellar objects due to external perturbations. This raises the question about the existence of a cloud of objects in the interstellar space which might overlap with our Oort cloud if we consider that other planetary systems face similar processes for the ejection of comets.


title: The Zeta Quotient $\zeta(3)/ \pi^3$ is Irrational

abstract: This note proves that the first odd zeta value does not have a closed form formula $\zeta(3)\ne r \pi^3$ for any rational number $r \in \mathbb{Q}$. Furthermore, assuming the irrationality of the second odd zeta value $\zeta(5)$, it is shown that $\zeta(5)\ne r\pi^5$ for any rational number $r \in \mathbb{Q}$.


title: Higher Segal spaces via higher excision

abstract: We show that the various higher Segal conditions of Dyckerhoff and Kapranov can all be characterized in purely categorical terms by higher excision conditions (in the spirit of Goodwillie-Weiss manifold calculus) on the simplex category $\Delta$ and the cyclic category $\Lambda$.


title: Isometry-Dual Flags of AG Codes

abstract: Consider a complete flag $\{0\} = C_0 < C_1 < \cdots < C_n = \mathbb{F}^n$ of one-point AG codes of length $n$ over the finite field $\mathbb{F}$. The codes are defined by evaluating functions with poles at a given point $Q$ in points $P_1,\dots,P_n$ distinct from $Q$. A flag has the isometry-dual property if the given flag and the corresponding dual flag are the same up to isometry. For several curves, including the projective line, Hermitian curves, Suzuki curves, Ree curves, and the Klein curve over the field of eight elements, the maximal flag, obtained by evaluation in all rational points different from the point $Q$, is self-dual. More generally, we ask whether a flag obtained by evaluation in a proper subset of rational points is isometry-dual. In [3] it is shown, for a curve of genus $g$, that a flag of one-point AG codes defined with a subset of $n > 2g+2$ rational points is isometry-dual if and only if the last code $C_n$ in the flag is defined with functions of pole order at most $n+2g-1$. Using a different approach, we extend this characterization to all subsets of size $n \geq 2g+2$. Moreover we show that this is best possible by giving examples of isometry-dual flags with $n=2g+1$ such that $C_n$ is generated by functions of pole order at most $n+2g-2$. We also prove a necessary condition, formulated in terms of maximum sparse ideals of the Weierstrass semigroup of $Q$, under which a flag of punctured one-point AG codes inherits the isometry-dual property from the original unpunctured flag.


title: Steady-state optimization of an exhaustive Levy storage process with intermittent output and random output rate

abstract: Consider a regenerative storage process with a nondecreasing L\'evy input (subordinator) such that every cycle may be split into two periods. In the first (off) the output is shut off and the workload accumulates. This continues until some stopping time. In the second (on), the process evolves like a subordinator minus a positive drift (output rate) until it hits the origin. In addition, we assume that the output rate of every on period is a random variable which is determined at the beginning of this period. For example, at each period, the output rate may depend on the workload level at the beginning of the corresponding busy period. We derive the Laplace-Stieltjes transform of the steady state distribution of the workload process and then apply this result to solve a steady-state cost minimization problem with holding, setup and output capacity costs. It is shown that the optimal output rate is a nondecreasing deterministic function of the workload level at the beginning of the corresponding on period.


title: The ampsys tool of pde2path

abstract: The computation of coefficients of amplitude systems for Turing bifurcations is a straightforward but sometimes elaborate task, in particular for 2D or 3D wave vector lattices. The Matlab tool ampsys automates such computations for two classes of problems, namely scalar equations of Swift-Hohenberg type and generalizations, and reaction-diffusion systems with an arbitrary number of components. The tool is designed to require minimal user input, and for a number of cases can also deal with symbolic computations. After a brief review of the setup of amplitude systems we explain the tool by a number of 1D, 2D and 3D examples over various wave vector lattices.


title: Emotion Recognition Using Fusion of Audio and Video Features

abstract: In this paper we propose a fusion approach to continuous emotion recognition that combines visual and auditory modalities in their representation spaces to predict the arousal and valence levels. The proposed approach employs a pre-trained convolution neural network and transfer learning to extract features from video frames that capture the emotional content. For the auditory content, a minimalistic set of parameters such as prosodic, excitation, vocal tract, and spectral descriptors are used as features. The fusion of these two modalities is carried out at a feature level, before training a single support vector regressor (SVR) or at a prediction level, after training one SVR for each modality. The proposed approach also includes preprocessing and post-processing techniques which contribute favorably to improving the concordance correlation coefficient (CCC). Experimental results for predicting spontaneous and natural emotions on the RECOLA dataset have shown that the proposed approach takes advantage of the complementary information of visual and auditory modalities and provides CCCs of 0.749 and 0.565 for arousal and valence, respectively.


title: On Capital Allocation under Information Constraints

abstract: Attempts to allocate capital across a selection of different investments are often hampered by the fact that investors' decisions are made under limited information (no historical return data) and during an extremely limited timeframe. Nevertheless, in some cases, rational investors with a certain level of experience are able to ordinally rank investment alternatives through relative assessments of the probabilities that investments will be successful. However, to apply traditional portfolio optimization models, analysts must use historical (or simulated/expected) return data as the basis for their calculations. This paper develops an alternative portfolio optimization framework that is able to handle this kind of information (given by an ordinal ranking of investment alternatives) and to calculate an optimal capital allocation based on a Cobb-Douglas function, which we call the Sorted Weighted Portfolio (SWP). Considering risk-neutral investors, we show that the results of this portfolio optimization model usually outperform the output generated by the (intuitive) Equally Weighted Portfolio (EWP) of different investment alternatives, which is the result of optimization when one is unable to incorporate additional data (the ordinal ranking of the alternatives). To further extend this work, we show that our model can also address risk-averse investors to capture correlation effects.


title: Antiforensic techniques deployed by custom developed malware in evading anti-virus detection

abstract: Both malware and antivirus detection tools advance in their capabilities. Malware aim is to evade the detection while antivirus is to detect the malware. Over time, the detection techniques evolved from simple static signature matching over antiheuristic analysis to machine learning assisted algorithms. This thesis describes several layers of anti-virus evasion deployed by the malware and conducts the analysis of the evasion success rate. The scientific contribution of this research is in the following techniques the malware used -- the new algorithm for identifying the Windows operating system functions, a new custom developed obfuscation and de-obfuscation routine and the usage of USB and sound devices enumeration in the anti-heuristic detection. The new PE mutation engine facilitates the malware static signature variation. In the next stage of the assessment, anti-virus engines then test the malware evasion capabilities. The locally installed antivirus applications and the two multi-scanner online engines inspect the submitted malware samples. The thesis examines the results and discusses the strengths and weaknesses of each evasion technique.


title: On compactifications of affine homology 3-cells into quadric fibrations

abstract: In this paper we deal with compactifications of affine homology $3$-cells into quadric fibrations such that the boundary divisors contain fibers. We show that all such affine homology $3$-cells are isomorphic to the affine $3$-space $\mathbb{A}^{3}$. Moreover, we show that all such compactifications can be connected by explicit elementary links preserving $\mathbb{A}^{3}$ to the projective $3$-space $\mathbb{P}^{3}$.


title: Dense networks that do not synchronize and sparse ones that do

abstract: For any network of identical Kuramoto oscillators with identical positive coupling, there is a critical connectivity above which the system is guaranteed to converge to the in-phase synchronous state, for almost all initial conditions. But the precise value of this critical connectivity remains unknown. In 2018, Ling, Xu, and Bandeira proved that if each oscillator is coupled to at least 79.29 percent of all the others, global synchrony is ensured. In 2019, Lu and Steinerberger improved this bound to 78.89 percent. Here, by focusing on circulant networks, we find clues that the critical connectivity may be exactly 75 percent. Our methods yield a slight improvement on the best known lower bound on the critical connectivity, from $68.18\%$ to $68.28\%$. We also consider the opposite end of the connectivity spectrum, where the networks are sparse rather than dense. In this regime, we ask how few edges one needs to add to a ring of $n$ oscillators to turn it into a globally synchronizing network. We prove a partial result: all the twisted states in a ring of size $n=2^m$ can be destabilized by adding just $\mathcal{O}(n \log_2 n)$ edges. To finish the proof, one also needs to rule out all other candidate attractors. We have done this for $n=8$ with computational algebraic geometry, but the problem remains open for larger $n$. Thus, even for systems as simple as Kuramoto oscillators, much remains to be learned about dense networks that do not globally synchronize and sparse ones that do.


title: Brezis-Nirenberg type result for Kohn Laplacian with critical Choquard Nonlinearity

abstract: In this article, we are study the following Dirichlet problem with Choquard type non linearity \[ -\Delta_{\mathbb{H}} u = a u+ \left(\int_{\Omega}\frac{|u(\eta)|^{Q^*_\lambda}}{|\eta^{-1}\xi|^{\lambda}}d\eta\right)|u|^{Q^*_\lambda-2}u \; \text{in}\; \Omega,\quad u = 0 \; \text{ on } \partial \Omega , \] where $\Omega$ is a smooth bounded subset of the Heisenberg group $\mathbb{H}^N, N\in \mathbb N$ with $C^2$ boundary and $\Delta_{\mathbb{H}}$ is the Kohn Laplacian on the Heisenberg group $\mathbb{H}^N$. Here, $Q^*_\lambda=\frac{2Q-\lambda}{Q-2},\; Q= 2N+2$ and $a$ is a positive real parameter. We derive the Brezis-Nirenberg type result for the above problem. Moreover, we also prove the regularity of solutions and nonexistence of solutions depending on the range of $a$.


title: Optical Reddening, Integrated HI Optical Depth, Total Hydrogen Column Density

abstract: Despite the vastly different angular scales on which they are measured, the integrated $\lambda$21 cm \HI\ optical depth measured interferometrically, \WHI, is a good proxy for the optical reddening derived from IR dust emission, with %\WHI\ $\equiv \int\tau({\rm \HI}) ~{\rm dv} = (13.8\pm0.7)$\EBV$^{(1.10\pm0.03)}$ \kms\ \WHI\ $\propto$ \EBV$^{1.10}$ for 0.04 mag $\la$ \EBV\ $\la$ 4 mag. For \EBV\ $\la 0.04$ mag or \WHI\ $< 0.7$ \kms, less-absorbent warm and ionized gases assert themselves and $\tau$(HI) is a less reliable tracer of \EBV. The \WHI-\EBV\ relationship can be inverted to give a broken power-law relationship between the total hydrogen column density N(H) and \WHI\ such that knowledge of \WHI\ alone predicts N(H) with an accuracy of a factor 1.5 ($\pm 0.18$ dex) across two orders of magnitude variation of \WHI. The \WHI-N(H) relation is invariant under a linear rescaling of the reddening measure used in the analysis and does not depend on knowing properties of the HI such as the spin temperature.


title: Timing Properties of Shocked Accretion Flows around Neutron Stars -- II. Viscous Disks and Boundary Layers

abstract: We use Smoothed Particle Hydrodynamics to study viscous accretion flows around a weakly magnetic neutron star. We show the formation of multiple ``boundary" layers in presence of both cooling and viscosity. We find that with the introduction of a small viscosity in a sub-Keplerian flow, much like the wind accretion in HMXBs such as Cir X-1, only a single Normal Boundary Layer (NBOL) forms to adjust the rotational velocity component. With the increase of viscosity, the region extends radially and beyond some critical value, a RAdiative KEplerian Disk/layer (RAKED) forms between the sub-Keplerian flow and the NBOL. When viscosity is increased further only NBOL and RAKED remain. In all such cases, the CENtrifugal pressure dominated BOundary Layer (CENBOL) is formed, away from the star, as in the case of black holes. This is the first self-consistent study where such a transition from sub-Keplerian flows has been reported for neutron stars. We also identify the connection between accretion and ejection of matter, following the Two-Component Advective Flow for black holes, for neutron stars. The results are crucial in the understanding of the formation of disks, boundary layers and outflows in wind dominated neutron star systems.


title: Towards global fits in EFT's and New Physics implications

abstract: I discuss recent progress on fits to dimension-six operators in the Standard Model Effective Theory (SMEFT). I focus on the top quark sector of the SMEFT, as well as the theoretical advances made in computing SMEFT effects through to next-to-leading order in QCD and the use of these calculations in global fits. I also discuss fits performed to the Higgs and electroweak sectors of the SMEFT and the possibility for performing global fits to multiple sectors simultaneously.


title: Unraveling the Mott-Peierls intrigue in Vanadium dioxide

abstract: Vanadium dioxide is one of the most studied strongly correlated materials. Nonetheless, the intertwining between electronic correlation and lattice effects has precluded a comprehensive description of the rutile metal to monoclinic insulator transition, in turn triggering a longstanding "the chicken or the egg" debate about which comes first, the Mott localisation or the Peierls distortion. Here, we suggest that this problem is in fact ill-posed: the electronic correlations and the lattice vibrations conspire to stabilise the monoclinic insulator, and so they must be both considered not to miss relevant pieces of the VO$_2$ physics. Specifically, we design a minimal model for VO$_2$ that includes all the important physical ingredients: the electronic correlations, the multi-orbital character, and the two components antiferrodistortive mode that condenses in the monoclinic insulator. We solve this model by dynamical mean-field theory within the adiabatic Born-Oppenheimer approximation. Consistently with the first-order character of the metal-insulator transition, the Born-Oppenheimer potential has a rich landscape, with minima corresponding to the undistorted phase and to the four equivalent distorted ones, and which translates into an equally rich thermodynamics that we uncover by the Monte Carlo method. Remarkably, we find that a distorted metal phase intrudes between the low-temperature distorted insulator and high-temperature undistorted metal, which sheds new light on the debated experimental evidence of a monoclinic metallic phase.


title: Cohomogeneity one Kaehler and Kaehler-Einstein manifolds with one singular orbit, II

abstract: F. Podest\`a and A. Spiro introduced a class of $G$-manifolds $M$ with a cohomogeneity one action of a compact semisimple Lie group $G$ which admit an invariant Kaehler structure $(g,J)$ (``standard $G$-manifolds") and studied invariant Kaehler and Kaehler-Einstein metrics on $M$. In the first part of this paper, we gave a combinatoric description of the standard non compact $G$-manifolds as the total space $M_{\varphi}$ of the homogeneous vector bundle $M = G\times_H V \to S_0 =G/H$ over a flag manifold $S_0$ and we gave necessary and sufficient conditions for the existence of an invariant Kaehler-Einstein metric $g$ on such manifolds $M$ in terms of the existence of an interval in the $T$-Weyl chamber of the flag manifold $F = G \times _H PV$ which satisfies some linear condition. In this paper, we consider standard cohomogeneity one manifolds of a classical simply connected Lie group $G = SU_n, Sp_n. Spin_n$ and reformulate these necessary and sufficient conditions in terms of easily checked arithmetic properties of the Koszul numbers associated with the flag manifold $S_0 = G/H$. If this conditions is fulfilled, the explicit construction of the Kaehler-Einstein metric reduces to the calculation of the inverse function to a given function of one variable.


title: Universal transient behavior in large dynamical systems on networks

abstract: We analyze how the transient dynamics of large dynamical systems in the vicinity of a stationary point, modeled by a set of randomly coupled linear differential equations, depends on the network topology. We characterize the transient response of a system through the evolution in time of the squared norm of the state vector, which is averaged over different realizations of the initial perturbation. We develop a mathematical formalism that computes this quantity for graphs that are locally tree-like. We show that for unidirectional networks the theory simplifies and general analytical results can be derived. For example, we derive analytical expressions for the average squared norm for random directed graphs with a prescribed degree distribution. These analytical results reveal that unidirectional systems exhibit a high degree of universality in the sense that the average squared norm only depends on a single parameter encoding the average interaction strength between the individual constituents. In addition, we derive analytical expressions for the average squared norm for unidirectional systems with fixed diagonal disorder and with bimodal diagonal disorder. We illustrate these results with numerical experiments on large random graphs and on real-world networks.

categories: nlin.AO cond-mat.dis-nn cond-mat.stat-mech math.DS


title: Widening Perspectives: The Intellectual and Social Benefits of Astrobiology, Big History, and the Exploration of Space

abstract: Astrobiology is the field of science devoted to searching for life elsewhere in the Universe. It is inherently interdisciplinary, integrating results from multiple fields of science, and in this respect has strong synergies with 'big history'. I argue that big history and astrobiology are both acting to widen human perspectives in intellectually and socially beneficial directions, especially by enhancing public awareness of cosmic and evolutionary worldviews. I will further argue that these perspectives have important implications for the social and political organisation of humanity, including the eventual political unification of our planet. Astrobiology and big history are also concerned with the future of humanity, and I will argue that this future will be culturally and intellectually enriched if it includes the exploration of the universe around us.


title: Measuring the frequency response of optically pumped metal-clad nanolasers

abstract: We report on our initial attempt to characterize the intrinsic frequency response of metal-clad nanolasers. The probed nanolaser is optically biased and modulated, allowing the emitted signal to be detected using a high-speed photodiode at each modulation frequency. Based on this technique, the prospect of high-speed operation of nanolasers is evaluated by measuring the D-factor, which is the ratio of the resonance frequency to the square root of its output power (fR/P_out^2). Our measurements show that for nanolasers, this factor is an order of magnitude greater than that of other state-of-the-art directly modulated semiconductor lasers. The theoretical analysis, based on the rate equation model and finite element method simulations of the cavity is in full agreement with the measurement results.


title: Coding for Crowdsourced Classification with XOR Queries

abstract: This paper models the crowdsourced labeling/classification problem as a sparsely encoded source coding problem, where each query answer, regarded as a code bit, is the XOR of a small number of labels, as source information bits. In this paper we leverage the connections between this problem and well-studied codes with sparse representations for the channel coding problem to provide querying schemes with almost optimal number of queries, each of which involving only a constant number of labels. We also extend this scenario to the case where some workers can be unresponsive. For this case, we propose querying schemes where each query involves only log n items, where n is the total number of items to be labeled. Furthermore, we consider classification of two correlated labeling systems and provide two-stage querying schemes with almost optimal number of queries each involving a constant number of labels.


title: Confluent conformal blocks and the Teukolsky master equation

abstract: Quasinormal modes of usual, four dimensional, Kerr black holes are described by certain solutions of a confluent Heun differential equation. In this work, we express these solutions in terms of the connection matrices for a Riemann-Hilbert problem, which was recently solved in terms of the Painlev\'e V transcendent. We use this formulation to generate the small-frequency expansion for the angular spheroidal harmonic eigenvalue, and derive conditions on the monodromy properties for the radial modes. Using exponentiation, we relate the accessory parameter to a semi-classical conformal description and discuss the properties of the operators involved. For the radial equation, while the operators at the horizons have Liouville momenta proportional to the entropy intake, we find that spatial infinity is described by a Whittaker operator.


title: Beyond the swap test: optimal estimation of quantum state overlap

abstract: We study the estimation of the overlap between two unknown pure quantum states of a finite dimensional system, given $M$ and $N$ copies of each type. This is a fundamental primitive in quantum information processing that is commonly accomplished from the outcomes of $N$ swap-tests, a joint measurement on one copy of each type whose outcome probability is a linear function of the squared overlap. We show that a more precise estimate can be obtained by allowing for general collective measurements on all copies. We derive the statistics of the optimal measurement and compute the optimal mean square error in the asymptotic pointwise and finite Bayesian estimation settings. Besides, we consider two strategies relying on the estimation of one or both the states, and show that, although they are suboptimal, they outperform the swap test. In particular, the swap test is extremely inefficient for small values of the overlap, which become exponentially more likely as the dimension increases. Finally, we show that the optimal measurement is less invasive than the swap test and study the robustness to depolarizing noise for qubit states.


title: SOS: Safe, Optimal and Small Strategies for Hybrid Markov Decision Processes

abstract: For hybrid Markov decision processes, UPPAAL Stratego can compute strategies that are safe for a given safety property and (in the limit) optimal for a given cost function. Unfortunately, these strategies cannot be exported easily since they are computed as a very long list. In this paper, we demonstrate methods to learn compact representations of the strategies in the form of decision trees. These decision trees are much smaller, more understandable, and can easily be exported as code that can be loaded into embedded systems. Despite the size compression and actual differences to the original strategy, we provide guarantees on both safety and optimality of the decision-tree strategy. On the top, we show how to obtain yet smaller representations, which are still guaranteed safe, but achieve a desired trade-off between size and optimality.


title: Micro Air Vehicle Link (MAVLink) in a Nutshell: A Survey

abstract: The Micro Air Vehicle Link (MAVLink in short) is a communication protocol for unmanned systems (e.g., drones, robots). It specifies a comprehensive set of messages exchanged between unmanned systems and ground stations. This protocol is used in major autopilot systems, mainly ArduPilot and PX4, and provides powerful features not only for monitoring and controlling unmanned systems missions but also for their integration into the Internet. However, there is no technical survey and/or tutorial in the literature that presents these features or explains how to make use of them. Most of the references are online tutorials and basic technical reports, and none of them presents comprehensive and systematic coverage of the protocol. In this paper, we address this gap, and we propose an overview of the MAVLink protocol, the difference between its versions, and its potential in enabling Internet connectivity to unmanned systems. We also discuss the security aspects of MAVLink. To the best of our knowledge, this is the first technical survey and tutorial on the MAVLink protocol, which represents an important reference for unmanned systems users and developers.


title: Validating Coordination Schemes between Transmission and Distribution System Operators using a Laboratory-Based Approach

abstract: The secure operation of future power systems will rely on better coordination between transmission system and distribution system operators. Increasing integration of renewables throughout the whole system is challenging the traditional operation. To tackle this problem, the SmartNet project proposes and evaluates five different coordination schemes between system operators using three benchmark scenarios from Denmark, Italy, and Spain. In the project, field tests in each of the benchmark countries are complemented with a number of laboratory validation tests, to cover scenarios that cannot be tested in field trials. This paper presents the outcome of these laboratory tests. Three tests are shown, focusing on controller validation, analysis of communication impacts, and how well price-based controls can integrate with the SmartNet coordination schemes. The results demonstrate important indications for the field tests and also show some of the limitations with the current implementations of the coordinations schemes.


title: A Review on Deep Learning in Medical Image Reconstruction

abstract: Medical imaging is crucial in modern clinics to guide the diagnosis and treatment of diseases. Medical image reconstruction is one of the most fundamental and important components of medical imaging, whose major objective is to acquire high-quality medical images for clinical usage at minimal cost and risk to the patients. Mathematical models in medical image reconstruction or, more generally, image restoration in computer vision, have been playing a prominent role. Earlier mathematical models are mostly designed by human knowledge or hypothesis on the image to be reconstructed, and we shall call these models handcrafted models. Later, handcrafted plus data-driven modeling started to emerge which still mostly relies on human designs, while part of the model is learned from the observed data. More recently, as more data and computation resources are made available, deep learning based models (or deep models) pushed data-driven modeling to the extreme where the models are mostly based on learning with minimal human designs. Both handcrafted and data-driven modeling have their own advantages and disadvantages. One of the major research trends in medical imaging is to combine handcrafted modeling with deep modeling so that we can enjoy benefits from both approaches. The major part of this article is to provide a conceptual review of some recent works on deep modeling from the unrolling dynamics viewpoint. This viewpoint stimulates new designs of neural network architectures with inspiration from optimization algorithms and numerical differential equations. Given the popularity of deep modeling, there are still vast remaining challenges in the field, as well as opportunities which we shall discuss at the end of this article.


title: Effects of Filler Configuration and Moisture on Dissipation Factor and Critical Electric Field of Epoxy Composites for HV-ICs Encapsulation

abstract: Molding compounds (MCs) have been used extensively as an encapsulation material for integrated circuits, however, MCs are susceptible to moisture and charge spreading over time. The increase in dissipation factor due to the increase of parasitic electrical conductivity ({\sigma}) and the decrease in dielectric strength (E_MC^Crit) restrict their applications. Thus, a fundamental understanding of moisture transport will suggest strategies to suppress moisture diffusion and broaden their applications. In this paper, we 1) propose a generalized effective medium and solubility (GEMS) Langmuir model to quantify water uptake as a function of filler configuration and relative humidity; 2) investigate the dominant impact of reacted-water on {\sigma} through numerical simulations, mass-uptake, and DC conductivity measurements; 3) investigate electric field distribution to explain how moisture ingress reduces E_MC^Crit; and finally 4) optimize the filler configuration to lower the dissipation factor, and enhance E_MC^Crit. The GEMS-Langmuir model can be used for any application (e.g., photovoltaics, biosensors) where moisture diffusion leads to reliability challenges.


title: Central Limit Theorems for Compound Paths on the 2-Dimensional Lattice

abstract: Zeckendorf proved that every integer can be written uniquely as a sum of non-consecutive Fibonacci numbers $\{F_n\}$, and later researchers showed that the distribution of the number of summands needed for such decompositions of integers in $[F_n, F_{n+1})$ converges to a Gaussian as $n\to\infty$. Decomposition problems have been studied extensively for a variety of different sequences and notions of a legal decompositions; for the Fibonacci numbers, a legal decomposition is one for which each summand is used at most once and no two consecutive summands may be chosen. Recently, Chen et al. [CCGJMSY] generalized earlier work to $d$-dimensional lattices of positive integers; there, a legal decomposition is a path such that every point chosen had each component strictly less than the component of the previous chosen point in the path. They were able to prove Gaussianity results despite the lack of uniqueness of the decompositions; however, their results should hold in the more general case where some components are identical. The strictly decreasing assumption was needed in that work to obtain simple, closed form combinatorial expressions, which could then be well approximated and led to the limiting behavior. In this work we remove that assumption through inclusion-exclusion arguments. These lead to more involved combinatorial sums; using generating functions and recurrence relations we obtain tractable forms in $2$ dimensions and prove Gaussianity again; a more involved analysis should work in higher dimensions.


title: Water vapor assisted sintering of silver nanoparticle inks for printed electronics

abstract: In printed electronics, conductive traces are often produced by printing inks of silver nanoparticles dispersed in solvents. A sintering process is usually needed to make the printed inks conductive by removing the organic dispersants and allowing metal-to-metal contacts among nanoparticles for atomic diffusion and neck formation. It has been challenging to sinter silver nanoparticle inks in a thermal oven at a temperature < 150 C to avoid thermal damage to the plastic substrate while achieving desired conductivity. This work presents a simple yet effective way to sinter a silver nanoparticle ink below 120 C (even at 80 C) by exposing the printed ink to water vapor in the oven. The results consistently show a significant reduction of line resistivity for the samples sintered in a moist oven compared to those sintered in a dry oven. Hence, solvent vapor-assisted sintering of metal nanoparticle inks can become an enabling approach to broaden product range of printed electronics.


title: Parameter Identification in Viscoplasticity using Transitional Markov Chain Monte Carlo Method

abstract: To evaluate the cyclic behavior under different loading conditions using the kinematic and isotropic hardening theory of steel, a Chaboche viscoplastic material model is employed. The parameters of a constitutive model are usually identified by minimization of the difference between model response and experimental data. However, measurement errors and differences in the specimens lead to deviations in the determined parameters. In this article, the Choboche model is used and a stochastic simulation technique is applied to generate artificial data which exhibit the same stochastic behavior as experimental data. Then the model parameters are identified by applying an estimation using Bayes's theorem. The Transitional Markov Chain Monte Carlo method (TMCMC) is introduced and employed to estimate the model parameters in the Bayesian setting. The uniform distributions of the parameters representing their priors are considered which literally means no knowledge of the parameters is available. Identified parameters are compared with the true parameters in the simulation, and the efficiency of the identification method is discussed. In fact, the main purpose of this study is to observe the possibility of identifying the model and hardening parameters of a viscoplastic model as a very high non-linear model with only a surface displacement measurement vector in the Bayesian setting using TMCMC and evaluate the number of measurements needed for a very acceptable estimation of the uncertain parameters of the model.


title: Global injectivity of differentiable maps via W-condition in R^2

abstract: In this paper, we study the intrinsic relation between the global injectivity of differentiable local homeomorphisms $F$ and the rate that tends to zero of $Spec(F)$ in $\mathbb{R}^2$, where $Spec(F)$ denotes the set of all (complex) eigenvalues of $DF(x)$, for all $x\in \mathbb{R}^2$. This depends on the $W$-condition deeply, which extends the $*$-condition and $B$-condition. The $W$-condition reveals the rate that tends to zero of real eigenvalues of $DF$ can not exceed $\displaystyle O\Big(x\ln x(\ln \frac{\ln x}{\ln\ln x})^2\Big)^{-1}$ by the half-Reeb component method. This improves the theorems of Guti\'{e}rrez-Nguyen \cite{GN07} and Rabanal \cite{RR10}. The $W$-condition is optimal for the half-Reeb component method in this paper setting.


title: Tunable Isolated Attosecond X-ray Pulses with Gigawatt Peak Power from a Free-Electron Laser

abstract: The quantum mechanical motion of electrons in molecules and solids occurs on the sub-femtosecond timescale. Consequently, the study of ultrafast electronic phenomena requires the generation of laser pulses shorter than 1 fs and of sufficient intensity to interact with their target with high probability. Probing these dynamics with atomic-site specificity requires the extension of sub-femtosecond pulses to the soft X-ray spectral region. Here we report the generation of isolated GW-scale soft X-ray attosecond pulses with an X-ray free-electron laser. Our source has a pulse energy that is six orders of magnitude larger than any other source of isolated attosecond pulses in the soft X-ray spectral region, with a peak power in the tens of gigawatts. This unique combination of high intensity, high photon energy and short pulse duration enables the investigation of electron dynamics with X-ray non-linear spectroscopy and single-particle imaging.


title: CHANG-ES XIX: Galaxy NGC 4013 -- a diffusion-dominated radio halo with plane-parallel disk and vertical halo magnetic fields

abstract: Using the Continuum HAloes in Nearby Galaxies - an EVLA Survey (CHANG-ES) radio continuum data from the Karl G. Jansky Very Large Array (VLA) in two frequency bands (C-band, L-band), we analyzed the radio properties, including polarization and the transport processes of the CR electrons (CREs), in the edge-on spiral galaxy NGC 4013. Supplementary LOw-Frequency ARray (LOFAR) data at 150MHz are used to study the low-frequency properties of this galaxy and X-ray (Chandra, XMM-Newton) data are used to investigate the central region. The central point source dominates the radio continuum in both CHANG-ES bands, but no clear AGN classification is possible at this time. The scale height analysis shows that Gaussian fits, with halo scale heights of 1.2 kpc in C-band, 2.0 kpc in L-band, and 3.1 kpc at 150 MHz, better represent the intensity profiles than do exponential fits. The radio continuum halo of NGC 4013 in C-band is rather small, while the low-frequency LOFAR data reveal a large halo. The polarization data reveal plane-parallel, regular magnetic fields within the entire disk and vertical halo components out to heights of about 6 kpc indicating the presence of an axisymmetric field having a radial component pointing outwards. The mean magnetic field strength of the disk of NGC 4013 of 6.6 $\mu$G (using the revised equipartition formula) is rather small. The interaction and the low star formation rate (SFR) across the disk of NGC 4013 probably influence the appearance of its radio continuum. Several observable quantities give consistent evidence that the CR transport in the halo of NGC 4013 is diffusive: the frequency dependence of the synchrotron scale height, the disk/halo flux density ratio, the vertical profile of the synchrotron spectral index, the small propagation speed measured modeled with spinnaker, and the low temperature of the X-ray emitting hot gas.


title: Interpretable Image Recognition with Hierarchical Prototypes

abstract: Vision models are interpretable when they classify objects on the basis of features that a person can directly understand. Recently, methods relying on visual feature prototypes have been developed for this purpose. However, in contrast to how humans categorize objects, these approaches have not yet made use of any taxonomical organization of class labels. With such an approach, for instance, we may see why a chimpanzee is classified as a chimpanzee, but not why it was considered to be a primate or even an animal. In this work we introduce a model that uses hierarchically organized prototypes to classify objects at every level in a predefined taxonomy. Hence, we may find distinct explanations for the prediction an image receives at each level of the taxonomy. The hierarchical prototypes enable the model to perform another important task: interpretably classifying images from previously unseen classes at the level of the taxonomy to which they correctly relate, e.g. classifying a hand gun as a weapon, when the only weapons in the training data are rifles. With a subset of ImageNet, we test our model against its counterpart black-box model on two tasks: 1) classification of data from familiar classes, and 2) classification of data from previously unseen classes at the appropriate level in the taxonomy. We find that our model performs approximately as well as its counterpart black-box model while allowing for each classification to be interpreted.


title: Monte Carlo Gradient Estimation in Machine Learning

abstract: This paper is a broad and accessible survey of the methods we have at our disposal for Monte Carlo gradient estimation in machine learning and across the statistical sciences: the problem of computing the gradient of an expectation of a function with respect to parameters defining the distribution that is integrated; the problem of sensitivity analysis. In machine learning research, this gradient problem lies at the core of many learning problems, in supervised, unsupervised and reinforcement learning. We will generally seek to rewrite such gradients in a form that allows for Monte Carlo estimation, allowing them to be easily and efficiently used and analysed. We explore three strategies--the pathwise, score function, and measure-valued gradient estimators--exploring their historical development, derivation, and underlying assumptions. We describe their use in other fields, show how they are related and can be combined, and expand on their possible generalisations. Wherever Monte Carlo gradient estimators have been derived and deployed in the past, important advances have followed. A deeper and more widely-held understanding of this problem will lead to further advances, and it is these advances that we wish to support.


title: Temperature dependence of normalized sensitivity of Love wave sensor with unidirectional carbon fiber epoxy composite/Mn-doped 0.24PIN-0.46PMN-0.30PT ternary single crystal configuration

abstract: We have derived a general formula for sensitivity optimization of gravimetric sensors and use it to design a high precision and high sensitivity gravimetric sensor using unidirectional carbon fiber epoxy composite (CFEC) guiding layer on single crystal Mn-doped yPb(In1/2Nb1/2)O3-(1-x-y)Pb(Mg1/3Nb2/3)O3-xPbTiO3 (Mn: PIN-PMN-PT) piezoelectric substrate. The normalized maximum sensitivity exhibits a decreasing tendency with temperature up to 55 degrees Celsius. For the CFEC-on-Mn: PIN-PMN-PT sensor configuration with wavelength 24 {mu}m at 25 degrees Celsius, the maximum sensitivity can reach as high as 760.88 cm2/g, which is nearly twice that of traditional SiO2/ST quartz configuration gravimetric sensor.


title: ReachNN: Reachability Analysis of Neural-Network Controlled Systems

abstract: Applying neural networks as controllers in dynamical systems has shown great promises. However, it is critical yet challenging to verify the safety of such control systems with neural-network controllers in the loop. Previous methods for verifying neural network controlled systems are limited to a few specific activation functions. In this work, we propose a new reachability analysis approach based on Bernstein polynomials that can verify neural-network controlled systems with a more general form of activation functions, i.e., as long as they ensure that the neural networks are Lipschitz continuous. Specifically, we consider abstracting feedforward neural networks with Bernstein polynomials for a small subset of inputs. To quantify the error introduced by abstraction, we provide both theoretical error bound estimation based on the theory of Bernstein polynomials and more practical sampling based error bound estimation, following a tight Lipschitz constant estimation approach based on forward reachability analysis. Compared with previous methods, our approach addresses a much broader set of neural networks, including heterogeneous neural networks that contain multiple types of activation functions. Experiment results on a variety of benchmarks show the effectiveness of our approach.


title: Complexity of Highly Parallel Non-Smooth Convex Optimization

abstract: A landmark result of non-smooth convex optimization is that gradient descent is an optimal algorithm whenever the number of computed gradients is smaller than the dimension $d$. In this paper we study the extension of this result to the parallel optimization setting. Namely we consider optimization algorithms interacting with a highly parallel gradient oracle, that is one that can answer $\mathrm{poly}(d)$ gradient queries in parallel. We show that in this case gradient descent is optimal only up to $\tilde{O}(\sqrt{d})$ rounds of interactions with the oracle. The lower bound improves upon a decades old construction by Nemirovski which proves optimality only up to $d^{1/3}$ rounds (as recently observed by Balkanski and Singer), and the suboptimality of gradient descent after $\sqrt{d}$ rounds was already observed by Duchi, Bartlett and Wainwright. In the latter regime we propose a new method with improved complexity, which we conjecture to be optimal. The analysis of this new method is based upon a generalized version of the recent results on optimal acceleration for highly smooth convex optimization.


title: A Unified Beamforming and A/D Self-Interference Cancellation Design for Full Duplex MIMO Radios

abstract: In this paper, we focus on reduced complexity full duplex Multiple-Input Multiple-Output (MIMO) systems and present a joint design of digital transmit and receive beamforming with Analog and Digital (A/D) self-interference cancellation. We capitalize on a recently proposed multi-tap analog canceller architecture, whose number of taps does not scale with the number of transceiver antennas, and consider practical transmitter impairments for the full duplex operation. Particularly, transmitter IQ imbalance and nonlinear power amplification are assumed via relevant realistic models. Aiming at suppressing the residual linear and nonlinear self-interference signal below the noise floor, we propose a novel digital self-interference cancellation technique that is jointly designed with the configuration of the analog taps and digital beamformers. Differently from the state of the art, we design pilot-assisted estimation of all involved wireless channels. Our representative Monte Carlo simulation results demonstrate that our unified full duplex MIMO design exhibits higher self-interference cancellation capability with less analog taps compared to available techniques, which results in improved achievable rate and bit error performance.


title: A Receding Horizon Framework for Autonomy in Unmanned Vehicles

abstract: In this article we present a unified framework based on receding horizon techniques that can be used to design the three tasks (guidance, navigation and path-planning) which are involved in the autonomy of unmanned vehicles. This tasks are solved using model predictive control and moving horizon estimation techniques, which allows us to include physical and dynamical constraints at the design stage, thus leading to optimal and feasible results. In order to demonstrate the capabilities of the proposed framework, we have used Gazebo simulator in order to drive a Jackal unmanned ground vehicle (UGV) along a desired path computed by the path-planning task. The results we have obtained are successful as the estimation and guidance errors are small and the Jackal UGV is able to follow the desired path satisfactorily and it is also capable to avoid the obstacles which are in its way.


title: The ideal structures of self-similar $k$-graph C*-algebras

abstract: Let $(G, \Lambda)$ be a self-similar $k$-graph with a possibly infinite vertex set $\Lambda^0$. We associate a universal C*-algebra $\mathcal{O}_{G,\Lambda}$ to $(G,\Lambda)$. The main purpose of this paper is to investigate the ideal structures of $\mathcal{O}_{G,\Lambda}$. We prove that there exists a one-to-one correspondence between the set of all $G$-hereditary and $G$-saturated subsets of $\Lambda^0$ and the set of all gauge-invariant and diagonal-invariant ideals of $\mathcal{O}_{G,\Lambda}$. Under some conditions, we characterize all primitive ideas of $\mathcal{O}_{G,\Lambda}$. Moreover, we describe the Jacobson topology of some concrete examples, which includes the C*-algebra of the product of odometers. On the way to our main results, we study self-similar $P$-graph C*-algebras in depth.


title: Dispersive sensing in hybrid InAs/Al nanowires

abstract: Dispersive charge sensing is realized in hybrid semiconductor-superconductor nanowires in gate-defined single- and double-island device geometries. Signal-to-noise ratios (SNRs) were measured both in the frequency and time domain. Frequency-domain measurements were carried out as a function of frequency and power and yield a charge sensitivity of $1 \times 10^{-3} e/\sqrt{\rm Hz}$ for an 11 MHz measurement bandwidth. Time-domain measurements yield SNR > 1 for 20 $\mu$s integration time. At zero magnetic field, photon-assisted tunneling was detected dispersively in a double-island geometry, indicating coherent hybridization of the two superconducting islands. At an axial magnetic field of 0.6 T, subgap states are detected dispersively, demonstrating the suitability of the method for sensing in the topological regime.


title: UVIT view of dwarf irregular galaxy IC 2574 : Is the star formation triggered due to expanding H$~$I shells?

abstract: Star forming dwarf galaxies in the local volume are diverse and are ideal test beds to understand details of star formation in a variety of environments. Here, we present a deep FUV imaging study of a nearby dwarf irregular galaxy IC 2574 using the Ultraviolet Imaging Telescope (UVIT). We identified 419 FUV bright regions with radii between 15 - 285 pc in the galaxy and found that 28.6\% of them to be located in H~I shells, 12.6\% inside holes and 60.1\% to be away from the holes. The H~I column density is found to be more than $10^{21} cm^{-2}$ for 82.3\% of the identified regions. 30 out of the 48 H~I holes show triggered star formation in their shells while 16 holes do not show any related FUV emission. Cross-matching with H$\alpha$ emission, we found that 23 holes have both FUV and H$\alpha$ emission in their shells signifying very recent trigger. Therefore, star formation in the galaxy has been partly triggered due to the expanding H~I holes whereas in majority of the sites it is driven by other mechanisms. Irrespective of the location, larger star forming complexes were found to have multiple sub-structures. We report two resolved components for the remnant cluster of the super giant shell and estimated their masses. The star formation rate of IC 2574 is found to be 0.57 $M_{\odot}$/yr, which is slightly higher compared to the average value of other nearby dwarf irregular galaxies.


title: Thermal Expansion Coefficient and Phonon Dynamics in Coexisting Allotropes of Monolayer WS2 Probed by Raman Scattering

abstract: We report a comprehensive temperature dependent Raman measurements on three different phases of monolayer WS2 from 4K to 330K in a wide spectral range. Our studies revels the anomalous nature of the first as well as the higher order combination modes reflected in the disappearance of the few modes and anomalous temperature evaluation of the phonon self-energy parameters attributed to the detuning of resonance condition and development of strain due to thermal expansion mismatch with the underlying substrate. Our detailed temperature dependence studies also decipher the ambiguity about assignment of the two modes in literature near ~ 297 cm-1 and 325 cm-1. Mode near 297 cm-1 is assigned as first order Raman mode, which is forbidden in the backscattering geometry and 325 cm-1 is assigned to the combination of and mode. We also estimated thermal expansion coefficient by systematically disentangling the substrate effect in the temperature range of 4K to 330K and probed its temperature dependence in 1H, 1T and 1T' phases.


title: Saving the universe with finite volume effects

abstract: Setting aside anthropic arguments, there is no reason why the universe should initially favour a net expanding phase rather than one experiencing a net contraction. However, a collapsing universe containing "normal" matter will end at a singularity in a finite time. We point out that there is a mechanism, derived from non-perturbative effects in Quantum Field Theory in a finite volume, which may provide a bias towards expansion when the spacetime volume shrinks, by dynamically violating the null energy condition, without the need for modified gravity or exotic matter. We describe a scalar field component subjected to this non-perturbative effect in a cosmological background and consider its impact on a contracting phase. We discuss how this could dynamically generate the necessary initial conditions for inflation to get started, or form part of the mechanism for a non-singular cosmological bounce.


title: Experimental demonstration of tunable graphene hyperbolic metamaterial

abstract: Tuning the macroscopic dielectric response on demand holds potential for actively tunable metaphotonics and optical devices. In recent years, graphene has been extensively investigated as a tunable element in nanophotonics. Significant theoretical work has been devoted on the tuning the hyperbolic properties of graphene/dielectric heterostructures, however, until now, such a motif has not been demonstrated experimentally. Here we focus on a graphene/polaritonic dielectric metamaterial, with strong optical resonances arising from the polar response of the dielectric, which are, in general, difficult to actively control. By controlling the doping level of graphene via external bias we experimentally demonstrate a wide range of tunability from a Fermi level of EF = 0 eV to EF = 0.5 eV, which yields an effective epsilon-near-zero crossing and tunable dielectric properties, verified through spectroscopic ellipsometry and transmission measurements.


title: Straggler Mitigation at Scale

abstract: Runtime performance variability at the servers has been a major issue, hindering the predictable and scalable performance in modern distributed systems. Executing requests or jobs redundantly over multiple servers has been shown to be effective for mitigating variability, both in theory and practice. Systems that employ redundancy has drawn significant attention, and numerous papers have analyzed the pain and gain of redundancy under various service models and assumptions on the runtime variability. This paper presents a cost (pain) vs. latency (gain) analysis of executing jobs of many tasks by employing replicated or erasure coded redundancy. Tail heaviness of service time variability is decisive on the pain and gain of redundancy and we quantify its effect by deriving expressions for the cost and latency. Specifically, we try to answer four questions: 1) How do replicated and coded redundancy compare in the cost vs. latency tradeoff? 2) Can we introduce redundancy after waiting some time and expect to reduce the cost? 3) Can relaunching the tasks that appear to be straggling after some time help to reduce cost and/or latency? 4) Is it effective to use redundancy and relaunching together? We validate the answers we found for each of the questions via simulations that use empirical distributions extracted from a Google cluster data.


title: The Aw-Rascle traffic model: Enskog-type kinetic derivation and generalisations

abstract: We study the derivation of second order macroscopic traffic models from kinetic descriptions. In particular, we recover the celebrated Aw-Rascle model as the hydrodynamic limit of an Enskog-type kinetic equation out of a precise characterisation of the microscopic binary interactions among the vehicles. Unlike other derivations available in the literature, our approach unveils the multiscale physics behind the Aw-Rascle model. This further allows us to generalise it to a new class of second order macroscopic models complying with the Aw-Rascle consistency condition, namely the fact that no wave should travel faster than the mean traffic flow.

categories: cond-mat.stat-mech cs.NA math.NA physics.soc-ph


title: Automatic Conversion from Flip-flop to 3-phase Latch-based Designs

abstract: Latch-based designs have many benefits over their flip-flop based counterparts but have limited use partially because most RTL specifications are flop-centric and automatic conversion of FF to latch-based designs is challenging. Conventional conversion algorithms target master-slave latch-based designs with two non-overlapping clocks. This paper presents a novel automated design flow that converts flip-flop to 3-phase latch-based designs. The resulting circuits have the same performance as the master-slave based designs but require significantly less latches. Our experimental results demonstrate the potential for savings in the number of latches (21.3%), area (5.8%), and power (16.3%) on a variety of ISCAS, CEP, and CPU benchmark circuits, compared to the master-slave conversions.


title: Reinforcement Learning with Competitive Ensembles of Information-Constrained Primitives

abstract: Reinforcement learning agents that operate in diverse and complex environments can benefit from the structured decomposition of their behavior. Often, this is addressed in the context of hierarchical reinforcement learning, where the aim is to decompose a policy into lower-level primitives or options, and a higher-level meta-policy that triggers the appropriate behaviors for a given situation. However, the meta-policy must still produce appropriate decisions in all states. In this work, we propose a policy design that decomposes into primitives, similarly to hierarchical reinforcement learning, but without a high-level meta-policy. Instead, each primitive can decide for themselves whether they wish to act in the current state. We use an information-theoretic mechanism for enabling this decentralized decision: each primitive chooses how much information it needs about the current state to make a decision and the primitive that requests the most information about the current state acts in the world. The primitives are regularized to use as little information as possible, which leads to natural competition and specialization. We experimentally demonstrate that this policy architecture improves over both flat and hierarchical policies in terms of generalization.


title: Discrete logarithms in quasi-polynomial time in finite fields of fixed characteristic

abstract: We prove that the discrete logarithm problem can be solved in quasi-polynomial expected time in the multiplicative group of finite fields of fixed characteristic. More generally, we prove that it can be solved in the field of cardinality $p^n$ in expected time $(pn)^{2\log_2(n) + O(1)}$.


title: Structural Design Using Laplacian Shells

abstract: We introduce a method to design lightweight shell objects that are structurally robust under the external forces they may experience during use. Given an input 3D model and a general description of the external forces, our algorithm generates a structurally-sound minimum weight shell object. Our approach works by altering the local shell thickness repeatedly based on the stresses that develop inside the object. A key issue in shell design is that large thickness values might result in self-intersections on the inner boundary creating a significant computational challenge during optimization. To address this, we propose a shape parametrization based on the solution to the Laplace's equation that guarantees smooth and intersection-free shell boundaries. Combined with our gradient-free optimization algorithm, our method provides a practical solution to the structural design of hollow objects with a single inner cavity. We demonstrate our method on a variety of problems with arbitrary 3D models under complex force configurations and validate its performance with physical experiments.


title: Improving performance of deep learning models with axiomatic attribution priors and expected gradients

abstract: Recent research has demonstrated that feature attribution methods for deep networks can themselves be incorporated into training; these attribution priors optimize for a model whose attributions have certain desirable properties -- most frequently, that particular features are important or unimportant. These attribution priors are often based on attribution methods that are not guaranteed to satisfy desirable interpretability axioms, such as completeness and implementation invariance. Here, we introduce attribution priors to optimize for higher-level properties of explanations, such as smoothness and sparsity, enabled by a fast new attribution method formulation called expected gradients that satisfies many important interpretability axioms. This improves model performance on many real-world tasks where previous attribution priors fail. Our experiments show that the gains from combining higher-level attribution priors with expected gradients attributions are consistent across image, gene expression, and health care data sets. We believe this work motivates and provides the necessary tools to support the widespread adoption of axiomatic attribution priors in many areas of applied machine learning. The implementations and our results have been made freely available to academic communities.


title: Explaining Deep Learning Models with Constrained Adversarial Examples

abstract: Machine learning algorithms generally suffer from a problem of explainability. Given a classification result from a model, it is typically hard to determine what caused the decision to be made, and to give an informative explanation. We explore a new method of generating counterfactual explanations, which instead of explaining why a particular classification was made explain how a different outcome can be achieved. This gives the recipients of the explanation a better way to understand the outcome, and provides an actionable suggestion. We show that the introduced method of Constrained Adversarial Examples (CADEX) can be used in real world applications, and yields explanations which incorporate business or domain constraints such as handling categorical attributes and range constraints.


title: Local-global principles for tori over arithmetic curves

abstract: In this paper we study local-global principles for tori over semi-global fields, which are one variable function fields over complete discretely valued fields. In particular, we show that for principal homogeneous spaces for tori over the underlying discrete valuation ring, the obstruction to a local-global principle with respect to discrete valuations can be computed using methods coming from patching. We give a sufficient condition for the vanishing of the obstruction, as well as examples were the obstruction is nontrivial or even infinite. A major tool is the notion of a flasque resolution of a torus.


title: Learning Fair and Transferable Representations

abstract: Developing learning methods which do not discriminate subgroups in the population is a central goal of algorithmic fairness. One way to reach this goal is by modifying the data representation in order to meet certain fairness constraints. In this work we measure fairness according to demographic parity. This requires the probability of the possible model decisions to be independent of the sensitive information. We argue that the goal of imposing demographic parity can be substantially facilitated within a multitask learning setting. We leverage task similarities by encouraging a shared fair representation across the tasks via low rank matrix factorization. We derive learning bounds establishing that the learned representation transfers well to novel tasks both in terms of prediction performance and fairness metrics. We present experiments on three real world datasets, showing that the proposed method outperforms state-of-the-art approaches by a significant margin.


title: Outlier eigenvalues for non-Hermitian polynomials in independent i.i.d. matrices and deterministic matrices

abstract: We consider a square random matrix of size $N$ of the form $P(Y,A)$ where $P$ is a noncommutative polynomial, $A$ is a tuple of deterministic matrices converging in $\ast$-distribution, when $N$ goes to infinity, towards a tuple $a$ in some $\mathcal{C}^*$-probability space and $Y$ is a tuple of independent matrices with i.i.d. centered entries with variance $1/N$. We investigate the eigenvalues of $P(Y,A)$ outside the spectrum of $P(c,a)$ where $c$ is a circular system which is free from $a$. We provide a sufficient condition to guarantee that these eigenvalues coincide asymptotically with those of $P(0,A)$.


title: Computational Investigations of the Lithium Superoxide Dimer Rearrangement on Noisy Quantum Devices

abstract: Currently available noisy intermediate-scale quantum (NISQ) devices are limited by the number of qubits that can be used for quantum chemistry calculations on molecules. We show herein that the number of qubits required for simulations on a quantum computer can be reduced by limiting the number of orbitals in the active space. Thus, we have utilized ans\"atze that approximate exact classical matrix eigenvalue decomposition methods (Full Configuration Interaction). Such methods are appropriate for computations with the Variational Quantum Eigensolver algorithm to perform computational investigations on the rearrangement of the lithium superoxide dimer with both quantum simulators and quantum devices. These results demonstrate that, even with a limited orbital active space, quantum simulators are capable of obtaining energy values that are similar to the exact ones. However, calculations on quantum hardware underestimate energies even after the application of readout error mitigation.


title: Spectral responsivity and photoconductive gain in thin film black phosphorus photodetectors

abstract: We have fabricated black phosphorus photodetectors and characterized their full spectral responsivity. These devices, which are effectively in the bulk thin film limit, show broadband responsivity ranging from <400 nm to the ~3.8 $\mu$m bandgap. In the visible, an intrinsic responsivity >7 A/W can be obtained due to internal gain mechanisms. By examining the full spectral response, we identify a sharp contrast between the visible and infrared behavior. In particular, the visible responsivity shows a large photoconductive gain and gate-voltge dependence, while the infrared responsivity is nearly independent of gate voltage and incident light intensity under most conditions. This is attributed to a contribution from the surface oxide. In addition, we find that the polarization anisotropy in responsivity along armchair and zigzag directions can be as large as 103 and extends from the band edge to 500 nm. The devices were fabricated in an inert atmosphere and encapsulated by Al$_2$O$_3$ providing stable operation for more than 6 months.