title: Time series of optical spectra of Nova V659 Sct

abstract: With our robotic 1.2 m TIGRE telescope, we were able to obtain eight optical spectra with intermediate resolution (R = 20,000) of the Nova V659 Sct during different phases of its outburst. We present a list of the lines found in the Nova spectra. The most common features are H I, O I, Na I, Fe II and Ca II. Studying the spectral evolution of the strong features we found that the absorption features move to higher expansion velocities before disappearing and the emission features show (different) asymmetries. Thanks to the intermediate spectral resolution we identified and analysed the interstellar medium absorption features present in the spectra. We detected atomic absorption features of Na I and Ca II. The sodium D lines show more complex substructures with three main absorption features at around a velocity of -10, 30 and 85 km/s. We identified several DIBs in the Nova V659 Sct spectra and determined their velocities and equivalent widths.


title: Some generalizations on affine invariant points

abstract: In this note we prove a more general (and topological) version of Gr\"unbaum's conjecture about affine invariant points. As an application of our result we show that, if we consider the action of the group of similarities, Gr\"unbaum's conjecture remains valid in other families of convex sets (not necessarily convex bodies).


title: Validating psychometric survey responses

abstract: We present an approach to classify user validity in survey responses by using a machine learning techniques. The approach is based on collecting user mouse activity on web-surveys and fast predicting validity of the survey in general without analysis of specific answers. Rule based approach, LSTM and HMM models are considered. The approach might be used in web-survey applications to detect suspicious users behaviour and request from them proper answering instead of false data recording.


title: Non-Markovian decoherence of a two-level system in a Lorentzian bosonic reservoir and a stochastic environment with finite correlation time

abstract: In this paper we investigate non-Markovian evolution of a two-level system (qubit) in a bosonic bath under influence of an external classical fluctuating environment. The interaction with the bath has the Lorentzian spectral density, and the fluctuating environment (stochastic field) is represented by a set of Ornstein-Uhlenbeck processes. Each of the subenvironments of the composite environment is able to induce non-Markovian dynamics of the two-level system. By means of the numerically exact method of hierarchical equations of motion, we study dependence of the steady states of the two-level system, the reduced density matrix evolution and the equilibrium emission spectrums on frequency cutoffs and coupling strengths of the subenvironments. Additionally we investigate the impact of the rotation-wave approximation (RWA) used for the interaction with the bath on accuracy of the results.


title: A Nonlinear Navigation Observer Using IMU and Generic Position Information

abstract: This paper deals with the problem of full state estimation for vehicles navigating in a three dimensional space. We assume that the vehicle is equipped with an Inertial Measurement Unit (IMU) providing body-frame measurements of the angular velocity, the specific force, and the Earth's magnetic field. Moreover, we consider available sensors that provide partial or full information about the position of the vehicle. Examples of such sensors are those which provide full position measurements (e.g., GPS), range measurements (e.g., Ultra-Wide Band (UWB) sensors), inertial-frame bearing measurements (e.g., motion capture cameras), or altitude measurements (altimeter). We propose a generic semi-globally exponentially stable nonlinear observer that estimates the position, linear velocity, linear acceleration, and attitude of the vehicle, as well as the gyro bias. We also provide a detailed observability analysis for different types of measurements. Simulation and experimental results are provided to demonstrate the effectiveness of the proposed estimation scheme.


title: Two quantum Ising algorithms for the Shortest Vector Problem: one for now and one for later

abstract: Quantum computers are expected to break today's public key cryptography within a few decades. New cryptosystems are being designed and standardised for the post-quantum era, and a significant proportion of these rely on the hardness of problems like the Shortest Vector Problem to a quantum adversary. In this paper we describe two variants of a quantum Ising algorithm to solve this problem. One variant is spatially efficient, requiring only O(NlogN) qubits where N is the lattice dimension, while the other variant is more robust to noise. Analysis of the algorithms' performance on a quantum annealer and in numerical simulations show that the more qubit-efficient variant will outperform in the long run, while the other variant is more suitable for near-term implementation.


title: Anycast Agility: Network Playbooks to Fight DDoS

abstract: IP anycast is used for services such as DNS and Content Delivery Networks (CDN) to provide the capacity to handle Distributed Denial-of-Service (DDoS) attacks. During a DDoS attack service operators redistribute traffic between anycast sites to take advantage of sites with unused or greater capacity. Depending on site traffic and attack size, operators may instead concentrate attackers in a few sites to preserve operation in others. Operators use these actions during attacks, but how to do so has not been described systematically or publicly. This paper describes several methods to use BGP to shift traffic when under DDoS, and shows that a response playbook can provide a menu of responses that are options during an attack. To choose an appropriate response from this playbook, we also describe a new method to estimate true attack size, even though the operator's view during the attack is incomplete. Finally, operator choices are constrained by distributed routing policies, and not all are helpful. We explore how specific anycast deployment can constrain options in this playbook, and are the first to measure how generally applicable they are across multiple anycast networks.


title: Distance bounds for high dimensional consistent digital rays and 2-D partially-consistent digital rays

abstract: We consider the problem of digitalizing Euclidean segments. Specifically, we look for a constructive method to connect any two points in $\mathbb{Z}^d$. The construction must be {\em consistent} (that is, satisfy the natural extension of the Euclidean axioms) while resembling them as much as possible. Previous work has shown asymptotically tight results in two dimensions with $\Theta(\log N)$ error, where resemblance between segments is measured with the Hausdorff distance, and $N$ is the $L_1$ distance between the two points. This construction was considered tight because of a $\Omega(\log N)$ lower bound that applies to any consistent construction in $\mathbb{Z}^2$. In this paper we observe that the lower bound does not directly extend to higher dimensions. We give an alternative argument showing that any consistent construction in $d$ dimensions must have $\Omega(\log^{1/(d-1)} N)$ error. We tie the error of a consistent construction in high dimensions to the error of similar {\em weak} constructions in two dimensions (constructions for which some points need not satisfy all the axioms). This not only opens the possibility for having constructions with $o(\log N)$ error in high dimensions, but also opens up an interesting line of research in the tradeoff between the number of axiom violations and the error of the construction. In order to show our lower bound, we also consider a colored variation of the concept of discrepancy of a set of points that we find of independent interest.


title: The Complete Redshift Distribution of Dusty Star-forming Galaxies from the SPT-SZ Survey

abstract: The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870\, \mathrm{\mu m}} > 25\, \mathrm{mJy}$) sample, initially selected at $1.4$ mm. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the $3$ mm spectral window, targeting carbon monoxide line emission. By combining these measurements with ancillary data, the SPT sample is now spectroscopically complete, with redshifts spanning $1.9$$<$$z$$<$$6.9$ and a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3\, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enabling the "high-redshift tail" of extremely luminous DSFGs to be measured. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts.


title: Pareto Active Learning with Gaussian Processes and Adaptive Discretization

abstract: We consider the problem of optimizing a vector-valued objective function $\boldsymbol{f}$ sampled from a Gaussian Process (GP) whose index set is a well-behaved, compact metric space $({\cal X},d)$ of designs. We assume that $\boldsymbol{f}$ is not known beforehand and that evaluating $\boldsymbol{f}$ at design $x$ results in a noisy observation of $\boldsymbol{f}(x)$. Since identifying the Pareto optimal designs via exhaustive search is infeasible when the cardinality of ${\cal X}$ is large, we propose an algorithm, called Adaptive $\boldsymbol{\epsilon}$-PAL, that exploits the smoothness of the GP-sampled function and the structure of $({\cal X},d)$ to learn fast. In essence, Adaptive $\boldsymbol{\epsilon}$-PAL employs a tree-based adaptive discretization technique to identify an $\boldsymbol{\epsilon}$-accurate Pareto set of designs in as few evaluations as possible. We provide both information-type and metric dimension-type bounds on the sample complexity of $\boldsymbol{\epsilon}$-accurate Pareto set identification. We also experimentally show that our algorithm outperforms other Pareto set identification methods on several benchmark datasets.


title: An $\ell_p$ theory of PCA and spectral clustering

abstract: Principal Component Analysis (PCA) is a powerful tool in statistics and machine learning. While existing study of PCA focuses on the recovery of principal components and their associated eigenvalues, there are few precise characterizations of individual principal component scores that yield low-dimensional embedding of samples. That hinders the analysis of various spectral methods. In this paper, we first develop an $\ell_p$ perturbation theory for a hollowed version of PCA in Hilbert spaces which provably improves upon the vanilla PCA in the presence of heteroscedastic noises. Through a novel $\ell_p$ analysis of eigenvectors, we investigate entrywise behaviors of principal component score vectors and show that they can be approximated by linear functionals of the Gram matrix in $\ell_p$ norm, which includes $\ell_2$ and $\ell_\infty$ as special examples. For sub-Gaussian mixture models, the choice of $p$ giving optimal bounds depends on the signal-to-noise ratio, which further yields optimality guarantees for spectral clustering. For contextual community detection, the $\ell_p$ theory leads to a simple spectral algorithm that achieves the information threshold for exact recovery. These also provide optimal recovery results for Gaussian mixture and stochastic block models as special cases.


title: Practical applications of metric space magnitude and weighting vectors

abstract: Metric space magnitude, an active subject of research in algebraic topology, originally arose in the context of biology, where it was used to represent the effective number of distinct species in an environment. In a more general setting, the magnitude of a metric space is a real number that aims to quantify the effective number of distinct points in the space. The contribution of each point to a metric space's global magnitude, which is encoded by the {\em weighting vector}, captures much of the underlying geometry of the original metric space. Surprisingly, when the metric space is Euclidean, the weighting vector also serves as an effective tool for boundary detection. This allows the weighting vector to serve as the foundation of novel algorithms for classic machine learning tasks such as classification, outlier detection and active learning. We demonstrate, using experiments and comparisons on classic benchmark datasets, the promise of the proposed magnitude and weighting vector-based approaches.


title: Derivatives, Eulerian polynomials and the $g$-indexes of Young tableaux

abstract: In this paper we first present summation formulas for $k$-order Eulerian polynomials and $1/k$-Eulerian polynomials. We then present combinatorial expansions of $(c(x)D)^n$ in terms of inversion sequences as well as $k$-Young tableaux, where $c(x)$ is a differentiable function in the indeterminate $x$ and $D$ is the derivative with respect to $x$. We define the $g$-indexes of $k$-Young tableaux and Young tableaux, which have important applications in combinatorics. By establishing some relations between $k$-Young tableaux and standard Young tableaux, we express Eulerian polynomials, second-order Eulerian polynomials, Andr\'e polynomials and the generating polynomials of gamma coefficients of Eulerian polynomials in terms of standard Young tableaux, which imply a deep connection among these polynomials.


title: Optical spectra of 2D monolayers from time-dependent density functional theory

abstract: The optical spectra of two-dimensional (2D) periodic systems provide a challenge for time-dependent density-functional theory (TDDFT) because of the large excitonic effects in these materials. In this work we explore how accurately these spectra can be described within a pure Kohn-Sham time-dependent density-functional framework, i.e., a framework in which no theory beyond Kohn-Sham density-functional theory, such as $GW$, is required to correct the Kohn-Sham gap. To achieve this goal we adapted a recent approach we developed for the optical spectra of 3D systems [Cavo, Berger, Romaniello, Phys. Rev. B 101, 115109 (2020)] to those of 2D systems. Our approach relies on the link between the exchange-correlation kernel of TDDFT and the derivative discontinuity of ground-state density-functional theory, which guarantees a correct quasi-particle gap, and on a generalization of the polarization functional [Berger, Phys. Rev. Lett., 115, 137402 (2015)], which describes the excitonic effects. We applied our approach to two prototypical 2D monolayers, $h$-BN and MoS$_2$. We find that our protocol gives a qualitative good description of the optical spectrum of $h$-BN, whereas improvements are needed for MoS$_2$ to describe the intensity of the excitonic peaks.


title: Expressive curves

abstract: We initiate the study of a class of real plane algebraic curves which we call expressive. These are the curves whose defining polynomial has the smallest number of critical points allowed by the topology of the set of real points of a curve. This concept can be viewed as a global version of the notion of a real morsification of an isolated plane curve singularity. We prove that a plane curve $C$ is expressive if (a) each irreducible component of $C$ can be parametrized by real polynomials (either ordinary or trigonometric), (b) all singular points of $C$ in the affine plane are ordinary hyperbolic nodes, and (c) the set of real points of $C$ in the affine plane is connected. Conversely, an expressive curve with real irreducible components must satisfy conditions (a)-(c), unless it exhibits some exotic behaviour at infinity. We describe several constructions that produce expressive curves, and discuss a large number of examples, including: arrangements of lines, parabolas, and circles; Chebyshev and Lissajous curves; hypotrochoids and epitrochoids; and much more.


title: On measuring two-party partisan bias in unbalanced states

abstract: Assuming that partisan fairness and responsiveness are important aspects of redistricting, it is important to measure them. Many measures of partisan bias are satisfactory for states that are balanced with roughly equal proportions of voters for the two major parties. It has been less clear which metrics measure fairness robustly when the proportion of the vote is unbalanced by as little as 60% to 40%. We have addressed this by analyzing past election results for four states with Democratic preferences (CA, IL, MA, and MD), three states with Republican preferences (SC, TN, and TX) and comparing those to results for four nearly balanced states (CO, NC, OH, and PA). We used many past statewide elections in each state to build statistically precise seats for votes and rank for votes graphs to which many measures of partisan bias were applied. In addition to providing values of responsiveness, we find that five of the measures of bias provide mutually consistent values in all states, thereby providing a core of usable measures for unbalanced states. Although all five measures focus on different aspects of partisan bias, normalization of the values across the eleven states provides a suitable way to compare them, and we propose that their average provides a superior measure which we call composite bias. Regarding other measures, we find that the most seemingly plausible symmetry measure fails for unbalanced states. We also consider deviations from the proportionality ideal, but using it is difficult because the political geography of a state can entangle responsiveness with total partisan bias. We do not attempt to separate intentional partisan bias from the implicit bias that results from the interaction of the map drawing rules of a state and its political geography, on the grounds that redistricting should attempt to minimize total partisan bias whatever its provenance.


title: Rapid infrared spectroscopic nano-imaging with nano-FTIR holography

abstract: Scattering-type scanning near-field optical microscopy (s-SNOM), and its derivate, Fourier transform infrared nanospectroscopy (nano-FTIR) are emerging techniques for infrared (IR) nanoimaging and spectroscopy with applications in diverse fields ranging from nanophotonics, chemical imaging and materials science. However, spectroscopic nanoimaging is still challenged by the limited acquisition rate of current nano-FTIR technology. Here we combine s-SNOM, nano-FTIR and synthetic optical holographic (SOH) to achieve infrared spectroscopic nanoimaging at unprecedented speed (8 spectroscopically resolved images in 20 min), which we demonstrate with a polymer composite sample. Beyond being fast, our method promises to enable nanoimaging in the long IR spectral range, which is covered by IR supercontinuum and synchrotron sources, but not by current quantum cascade laser technology.


title: Bounds of the sum of edge lengths in linear arrangements of trees

abstract: A fundamental problem in network science is the normalization of the topological or physical distance between vertices, that requires understanding the range of variation of the unnormalized distances. Here we investigate the limits of the variation of the physical distance in linear arrangements of the vertices of trees. In particular, we investigate various problems on the sum of edge lengths in trees of a fixed size: the minimum and the maximum value of the sum for specific trees, the minimum and the maximum in classes of trees (bistar trees and caterpillar trees) and finally the minimum and the maximum for any tree. We establish some foundations for research on optimality scores for spatial networks in one dimension.


title: Enumeration of Standard Puzzles

abstract: We introduce a large family of combinatorial objects, called standard puzzles, defined by very simple rules. We focus on the standard puzzles for which the enumeration problems can be solved by explicit formulas or by classical numbers, such as binomial coefficients, Fibonacci numbers, tangent numbers, Catalan numbers, $\ldots$


title: Strain-engineering the topological type-II Dirac semimetal NiTe$_2$

abstract: In the present work, we investigated the electronic and elastic properties in equilibrium and under strain of the type-II Dirac semimetal NiTe$_2$ using density functional theory (DFT). Our results demonstrate the tunability of Dirac nodes' energy and momentum with strain and that it is possible to bring them closer to the Fermi level, while other metallic bands are supressed. We also derive a minimal 4-band effective model for the Dirac cones which accounts for the aforementioned strain effects by means of lattice regularization, providing an inexpensive way for further theoretical investigations and easy comparison with experiments. On an equal footing, we propose the static control of the electronic structure by intercalating alkali species into the van der Waals gap, resulting in the same effects obtained by strain-engineering and removing the requirement of in situ strain. Finally, evaluating the wavefunction's symmetry evolution as the lattice is deformed, we discuss possible consequences, such as Liftshitz transitions and the coexistence of type-I and type-II Dirac cones, thus motivating future investigations.


title: Estrellas h\'ibridas con n\'ucleos superconductores de color

abstract: The aim of this work is the study of hybrid stars interiors from the equations of state describing the matter composing them. We used modern hadronic equations of state to describe the matter in the outer core of these stars. The quark phase is modeled by an extended version of the Field Correlator Method, in which vector interactions among quarks and color superconductivity are included. We analyze the effect of these two contributions on the mass-radius relationship and tidal deformability. This is a relevant quantity for neutron stars after the GW170817 event, allowing to impose additional constraints to the equation of state of these compact objetcs.


title: Breakdown of the Einstein's Equivalence Principle for a quantum body

abstract: We review our recent theoretical results about inequivalence between passive gravitational mass and energy for a composite quantum body at a macroscopic level. In particular, we consider macroscopic ensembles of the simplest composite quantum bodies - hydrogen atoms. Our results are as follows. For the most ensembles, the Einstein's Equivalence Principle is valid. On the other hand, we discuss that for some special quantum ensembles - ensembles of the coherent superpositions of the stationary quantum states in the hydrogen atoms (which we call Gravitational demons) - the Equivalence Principle between passive gravitational mass and energy is broken. We show that, for such superpositions, the expectation values of passive gravitational masses are not related to the expectation values of energies by the famous Einstein's equation, i.e, $m_g \neq \frac{E}{c^2}$. Possible experiments at the Earth's laboratories are briefly discussed, in contrast to the numerous attempts and projects to discover the possible breakdown of the Einstein's Equivalence Principle during the space missions.


title: Design And Develop Network Storage Virtualization By Using GNS3

abstract: Virtualization is an emerging and optimistic prospect in the IT industry. Its impact has a footprint widely in digital infrastructure. Many innovativeness sectors utilized the concept of virtualization to reduce the cost of frameworks. In this paper, we have designed and developed storage virtualization for physical functional solutions. It is an auspicious type of virtualization that is accessible, secure, scalable, and manageable. In the paper, we have proposed the pool storage method used the RAID-Z file system with the ZFS model which provides the duplication of site approach, compression blueprint, adequate backup methods, expansion in error-correcting techniques, and tested procedure on the real-time network location. Therefore, this study provides useful guidelines to design and develop optimized storage virtualization.


title: Higgs modulation of emergent mass as revealed in kaon and pion parton distributions

abstract: Strangeness was discovered roughly seventy years ago, lodged in a particle now known as the kaon, $K$. Kindred to the pion, $\pi$; both states are massless in the absence of Higgs-boson couplings. Kaons and pions are Nature's most fundamental Nambu-Goldstone modes. Their properties are largely determined by the mechanisms responsible for emergent mass in the standard model, but modulations applied by the Higgs are crucial to Universe evolution. Despite their importance, little is known empirically about $K$ and $\pi$ structure. This study delivers the first parameter-free predictions for all $K$ distribution functions (DFs) and comparisons with the analogous $\pi$ distributions, i.e. the one-dimensional maps that reveal how the light-front momentum of these states is shared amongst the gluons and quarks from which they are formed. The results should stimulate improved analyses of existing data and motivate new experiments sensitive to all $K$ and $\pi$ DFs.


title: The Convex Relaxation Barrier, Revisited: Tightened Single-Neuron Relaxations for Neural Network Verification

abstract: We improve the effectiveness of propagation- and linear-optimization-based neural network verification algorithms with a new tightened convex relaxation for ReLU neurons. Unlike previous single-neuron relaxations which focus only on the univariate input space of the ReLU, our method considers the multivariate input space of the affine pre-activation function preceding the ReLU. Using results from submodularity and convex geometry, we derive an explicit description of the tightest possible convex relaxation when this multivariate input is over a box domain. We show that our convex relaxation is significantly stronger than the commonly used univariate-input relaxation which has been proposed as a natural convex relaxation barrier for verification. While our description of the relaxation may require an exponential number of inequalities, we show that they can be separated in linear time and hence can be efficiently incorporated into optimization algorithms on an as-needed basis. Based on this novel relaxation, we design two polynomial-time algorithms for neural network verification: a linear-programming-based algorithm that leverages the full power of our relaxation, and a fast propagation algorithm that generalizes existing approaches. In both cases, we show that for a modest increase in computational effort, our strengthened relaxation enables us to verify a significantly larger number of instances compared to similar algorithms.


title: Time for a Background Check! Uncovering the impact of Background Features on Deep Neural Networks

abstract: With increasing expressive power, deep neural networks have significantly improved the state-of-the-art on image classification datasets, such as ImageNet. In this paper, we investigate to what extent the increasing performance of deep neural networks is impacted by background features? In particular, we focus on background invariance, i.e., accuracy unaffected by switching background features and background influence, i.e., predictive power of background features itself when foreground is masked. We perform experiments with 32 different neural networks ranging from small-size networks to large-scale networks trained with up to one Billion images. Our investigations reveal that increasing expressive power of DNNs leads to higher influence of background features, while simultaneously, increases their ability to make the correct prediction when background features are removed or replaced with a randomly selected texture-based background.


title: Machine learning the real discriminant locus

abstract: Parameterized systems of polynomial equations arise in many applications in science and engineering with the real solutions describing, for example, equilibria of a dynamical system, linkages satisfying design constraints, and scene reconstruction in computer vision. Since different parameter values can have a different number of real solutions, the parameter space is decomposed into regions whose boundary forms the real discriminant locus. This article views locating the real discriminant locus as a supervised classification problem in machine learning where the goal is to determine classification boundaries over the parameter space, with the classes being the number of real solutions. For multidimensional parameter spaces, this article presents a novel sampling method which carefully samples the parameter space. At each sample point, homotopy continuation is used to obtain the number of real solutions to the corresponding polynomial system. Machine learning techniques including nearest neighbor and deep learning are used to efficiently approximate the real discriminant locus. One application of having learned the real discriminant locus is to develop a real homotopy method that only tracks the real solution paths unlike traditional methods which track all~complex~solution~paths. Examples show that the proposed approach can efficiently approximate complicated solution boundaries such as those arising from the equilibria of the Kuramoto model.


title: Ensuring Learning Guarantees on Concept Drift Detection with Statistical Learning Theory

abstract: Concept Drift (CD) detection intends to continuously identify changes in data stream behaviors, supporting researchers in the study and modeling of real-world phenomena. Motivated by the lack of learning guarantees in current CD algorithms, we decided to take advantage of the Statistical Learning Theory (SLT) to formalize the necessary requirements to ensure probabilistic learning bounds, so drifts would refer to actual changes in data rather than by chance. As discussed along this paper, a set of mathematical assumptions must be held in order to rely on SLT bounds, which are especially controversial in CD scenarios. Based on this issue, we propose a methodology to address those assumptions in CD scenarios and therefore ensure learning guarantees. Complementary, we assessed a set of relevant and known CD algorithms from the literature in light of our methodology. As main contribution, we expect this work to support researchers while designing and evaluating CD algorithms on different domains.


title: Accelerating MRI Reconstruction on TPUs

abstract: The advanced magnetic resonance (MR) image reconstructions such as the compressed sensing and subspace-based imaging are considered as large-scale, iterative, optimization problems. Given the large number of reconstructions required by the practical clinical usage, the computation time of these advanced reconstruction methods is often unacceptable. In this work, we propose using Google's Tensor Processing Units (TPUs) to accelerate the MR image reconstruction. TPU is an application-specific integrated circuit (ASIC) for machine learning applications, which has recently been used to solve large-scale scientific computing problems. As proof-of-concept, we implement the alternating direction method of multipliers (ADMM) in TensorFlow to reconstruct images on TPUs. The reconstruction is based on multi-channel, sparsely sampled, and radial-trajectory $k$-space data with sparsity constraints. The forward and inverse non-uniform Fourier transform operations are formulated in terms of matrix multiplications as in the discrete Fourier transform. The sparsifying transform and its adjoint operations are formulated as convolutions. The data decomposition is applied to the measured $k$-space data such that the aforementioned tensor operations are localized within individual TPU cores. The data decomposition and the inter-core communication strategy are designed in accordance with the TPU interconnect network topology in order to minimize the communication time. The accuracy and the high parallel efficiency of the proposed TPU-based image reconstruction method are demonstrated through numerical examples.


title: Local structure of glassy lithium phosphorus oxynitride thin films: a combined experimental and ab initio approach

abstract: Lithium phosphorus oxynitride (LiPON) is an amorphous solid-state lithium ion conductor displaying exemplary cyclability against lithium metal anodes. There is no definitive explanation for this stability due to the limited understanding of the structure of LiPON. We provide a structural model of RF-sputtered LiPON via experimental and computational spectroscopic methods. Information about the short-range structure results from 1D and 2D solid-state nuclear magnetic resonance experiments investigating chemical shift anisotropy and dipolar interactions. These results are compared with first principles chemical shielding calculations of Li-P-O/N crystals and ab initio molecular dynamics-generated amorphous LiPON models to unequivocally identify the glassy structure as primarily isolated phosphate monomers with N incorporated in both apical and as bridging sites in phosphate dimers. Structural results suggest LiPON's stability is a result of its glassy character. Free-standing LiPON films are produced that exhibit a high degree of flexibility highlighting the unique mechanical properties of glassy materials.


title: Discontinuous Galerkin for the wave equation: a simplified a priori error analysis

abstract: Standard discontinuous Galerkin methods, based on piecewise polynomials of degree $ \qq=0,1$, are considered for temporal semi-discretization for second order hyperbolic equations. The main goal of this paper is to present a simple and straightforward a priori error analysis of optimal order with minimal regularity requirement on the solution. Uniform norm in time error estimates are also proved. To this end, energy identities and stability estimates of bthe discrete problem are proved for a slightly more general problem. These are used to prove optimal order a priori error estimates with minimal regularity requirement on the solution. The combination with the classic continuous Galerkin finite element discretization in space variable is used, to formulate a full-discrete scheme. The a priori error analysis is presented. Numerical experiments are performed to verify the theoretical results.


title: A q deformation of true-polyanalytic Bargmann transforms when q^{-1}> 1

abstract: We combine continuous $q^{-1}$-Hermite Askey polynomials with new $2D$ orthogonal polynomials introduced by Ismail and Zhang as $q$-analogs for complex Hermite polynomials to construct a new set of coherent states depending on a nonnegative integer parameter $m$. In the analytic case corresponding to $m=0$, we recover a known result on the Ar\"{\i}k-Coon oscillator for $q'=q^{-1}>1$. Our construction leads to a new $q$-deformation of the $m$-true-polyanalytic Bargmann transform on the complex plane. The obtained result may be used to introduce a $q$-deformed Ginibre-type point process.


title: Multilabel Classification by Hierarchical Partitioning and Data-dependent Grouping

abstract: In modern multilabel classification problems, each data instance belongs to a small number of classes from a large set of classes. In other words, these problems involve learning very sparse binary label vectors. Moreover, in large-scale problems, the labels typically have certain (unknown) hierarchy. In this paper we exploit the sparsity of label vectors and the hierarchical structure to embed them in low-dimensional space using label groupings. Consequently, we solve the classification problem in a much lower dimensional space and then obtain labels in the original space using an appropriately defined lifting. Our method builds on the work of (Ubaru & Mazumdar, 2017), where the idea of group testing was also explored for multilabel classification. We first present a novel data-dependent grouping approach, where we use a group construction based on a low-rank Nonnegative Matrix Factorization (NMF) of the label matrix of training instances. The construction also allows us, using recent results, to develop a fast prediction algorithm that has a logarithmic runtime in the number of labels. We then present a hierarchical partitioning approach that exploits the label hierarchy in large scale problems to divide up the large label space and create smaller sub-problems, which can then be solved independently via the grouping approach. Numerical results on many benchmark datasets illustrate that, compared to other popular methods, our proposed methods achieve competitive accuracy with significantly lower computational costs.


title: Topological Insights into Sparse Neural Networks

abstract: Sparse neural networks are effective approaches to reduce the resource requirements for the deployment of deep neural networks. Recently, the concept of adaptive sparse connectivity, has emerged to allow training sparse neural networks from scratch by optimizing the sparse structure during training. However, comparing different sparse topologies and determining how sparse topologies evolve during training, especially for the situation in which the sparse structure optimization is involved, remain as challenging open questions. This comparison becomes increasingly complex as the number of possible topological comparisons increases exponentially with the size of networks. In this work, we introduce an approach to understand and compare sparse neural network topologies from the perspective of graph theory. We first propose Neural Network Sparse Topology Distance (NNSTD) to measure the distance between different sparse neural networks. Further, we demonstrate that sparse neural networks can outperform over-parameterized models in terms of performance, even without any further structure optimization. To the end, we also show that adaptive sparse connectivity can always unveil a plenitude of sparse sub-networks with very different topologies which outperform the dense model, by quantifying and comparing their topological evolutionary processes. The latter findings complement the Lottery Ticket Hypothesis by showing that there is a much more efficient and robust way to find "winning tickets". Altogether, our results start enabling a better theoretical understanding of sparse neural networks, and demonstrate the utility of using graph theory to analyze them.


title: The flag manifold as a tool for analyzing and comparing data sets

abstract: The shape and orientation of data clouds reflect variability in observations that can confound pattern recognition systems. Subspace methods, utilizing Grassmann manifolds, have been a great aid in dealing with such variability. However, this usefulness begins to falter when the data cloud contains sufficiently many outliers corresponding to stray elements from another class or when the number of data points is larger than the number of features. We illustrate how nested subspace methods, utilizing flag manifolds, can help to deal with such additional confounding factors. Flag manifolds, which are parameter spaces for nested subspaces, are a natural geometric generalization of Grassmann manifolds. To make practical comparisons on a flag manifold, algorithms are proposed for determining the distances between points $[A], [B]$ on a flag manifold, where $A$ and $B$ are arbitrary orthogonal matrix representatives for $[A]$ and $[B]$, and for determining the initial direction of these minimal length geodesics. The approach is illustrated in the context of (hyper) spectral imagery showing the impact of ambient dimension, sample dimension, and flag structure.


title: Tidal tails of open star clusters as probes to early gas expulsion I: A semi-analytic model

abstract: Star clusters form out of the densest parts of infrared dark clouds. The emergence of massive stars expels the residual gas, which has not formed stars yet. Gas expulsion lowers the gravitational potential of the embedded cluster, unbinding many of the cluster stars. These stars then move on their own trajectories in the external gravitational field of the Galaxy, forming a tidal tail. We investigate the formation and evolution of the tidal tail forming due to expulsion of primordial gas under various scenarios of gas expulsion to provide predictions for tidal tails around dynamically evolved (age > 100 Myr) galactic star clusters, which can be possibly detected by the Gaia mission. We provide a semi-analytical model for the tail evolution. We find that tidal tails released during gas expulsion have different kinematic properties than the tails gradually forming due to evaporation. The gas expulsion tidal tail shows non-monotonic expansion with time, where longer epochs of expansion are interspersed with shorter epochs of contraction. The tail thickness and velocity dispersions strongly, but not exactly periodically, vary with time. The times of minima of tail thickness and velocity dispersions are given only by the properties of the galactic potential, and not by the properties of the cluster. The estimates provided by the (semi-)analytical model for the extent of the tail, the minima of tail thickness, and velocity dispersions are in a very good agreement with the nbody6 simulations. This implies that the semi-analytic model can be used for estimating the properties of the gas expulsion tidal tail for a cluster of a given age and orbital parameters without the necessity of performing numerical simulations. A study with a more extended parameter space of the initial conditions is performed in the follow up paper.


title: Cheating Robot Games: A model for insider information

abstract: Combinatorial games are two-player games of pure strategy where the players, usually called Left and Right, move alternately. In this paper, we introduce Cheating Robot games. These arise from simultaneous-play combinatorial games where one player has insider information ('cheats'). Play occurs in rounds. At the beginning of a round, both players know the moves that are available to them. Left chooses a move. Knowing Left's move, Right then chooses a move. Right's move is not constrained by Left's choice. The round is not completed until both players have made a choice. A game is finished only when one or both players do not have a move at the beginning of a round. Right choosing a move, knowing Left's, makes the games deterministic, distinguishing them from simultaneous games. Also, the ending condition distinguishes this class of games from combinatorial games, since the outcomes are now Left-win, Right-win and draw. The basic theory and properties are developed, including showing that there is an equivalence relation and partial order on the games. Whilst there are no inverses in the class of all games, we show that there is a sub-class, simple hot games, in which the integers have inverses. In this sub-class, the optimal strategies are obtained by the solutions to a minimum-weight matching problem on a graph whose number of vertices equals the number of summands in the disjunctive sum.


title: Electric But Not Eclectic: Thermal Relic Dark Matter for the XENON1T Excess

abstract: The identity of dark matter is being sought with increasingly sensitive and voluminous underground detectors. Recently the XENON1T collaboration reported excess electronic recoil events, with most of these having recoil energies around $1-30$ keV. We show that a straightforward model of inelastic dark matter produced via early universe thermal freeze-out annihilation can account for the XENON1T excess. Remarkably, this dark matter model consists of a few simple elements: sub-GeV mass Dirac fermion dark matter coupled to a lighter dark photon kinetically mixed with the Standard Model photon. A scalar field charged under the dark U(1) gauge symmetry can provide a mass for the dark photon and splits the Dirac fermion component state masses by a few keV, which survive in equal abundance and interact inelastically with electrons and nuclei.


title: Neural Architecture Design for GPU-Efficient Networks

abstract: Many mission-critical systems are based on GPU for inference. It requires not only high recognition accuracy but also low latency in responding time. Although many studies are devoted to optimizing the structure of deep models for efficient inference, most of them do not leverage the architecture of \textbf{modern GPU} for fast inference, leading to suboptimal performance. To address this issue, we propose a general principle for designing GPU-efficient networks based on extensive empirical studies. This design principle enables us to search for GPU-efficient network structures effectively by a simple and lightweight method as opposed to most Neural Architecture Search (NAS) methods that are complicated and computationally expensive. Based on the proposed framework, we design a family of GPU-Efficient Networks, or GENets in short. We did extensive evaluations on multiple GPU platforms and inference engines. While achieving $\geq 81.3\%$ top-1 accuracy on ImageNet, GENet is up to $6.4$ times faster than EfficienNet on GPU. It also outperforms most state-of-the-art models that are more efficient than EfficientNet in high precision regimes. Our source code and pre-trained models are available from \url{https://github.com/idstcv/GPU-Efficient-Networks}.


title: Learning Reward Functions from Diverse Sources of Human Feedback: Optimally Integrating Demonstrations and Preferences

abstract: Reward functions are a common way to specify the objective of a robot. As designing reward functions can be extremely challenging, a more promising approach is to directly learn reward functions from human teachers. Importantly, data from human teachers can be collected either passively or actively in a variety of forms: passive data sources include demonstrations, (e.g., kinesthetic guidance), whereas preferences (e.g., comparative rankings) are actively elicited. Prior research has independently applied reward learning to these different data sources. However, there exist many domains where multiple sources are complementary and expressive. Motivated by this general problem, we present a framework to integrate multiple sources of information, which are either passively or actively collected from human users. In particular, we present an algorithm that first utilizes user demonstrations to initialize a belief about the reward function, and then actively probes the user with preference queries to zero-in on their true reward. This algorithm not only enables us combine multiple data sources, but it also informs the robot when it should leverage each type of information. Further, our approach accounts for the human's ability to provide data: yielding user-friendly preference queries which are also theoretically optimal. Our extensive simulated experiments and user studies on a Fetch mobile manipulator demonstrate the superiority and the usability of our integrated framework.


title: A Relaxation/Finite Difference discretization of a 2D Semilinear Heat Equation over a rectangular domain

abstract: We consider an initial and Dirichlet boundary value problem for a semilinear, two dimensional heat equation over a rectangular domain. The problem is discretized in time by a version of the Relaxation Scheme proposed by C. Besse (C. R. Acad. Sci. Paris S\'er. I, vol. 326 (1998)) for the nonlinear Schr\"odinger equation and in space by a standard second order finite difference method. The proposed method is unconditionally well-posed and its convergence is established by proving an optimal second order error estimate allowing a mild mesh condition to hold.


title: A Linear-Time Algorithm for Discrete Radius Optimally Augmenting Paths in a Metric Space

abstract: Let $P$ be a path graph of $n$ vertices embedded in a metric space. We consider the problem of adding a new edge to $P$ so that the radius of the resulting graph is minimized, where any center is constrained to be one of the vertices of $P$. Previously, the "continuous" version of the problem where a center may be a point in the interior of an edge of the graph was studied and a linear-time algorithm was known. Our "discrete" version of the problem has not been studied before. We present a linear-time algorithm for the problem.


title: Asymptotic method of moving planes for fractional parabolic equations

abstract: In this paper, we develop a systematical approach in applying an asymptotic method of moving planes to investigate qualitative properties of positive solutions for fractional parabolic equations. We first obtain a series of needed key ingredients such as narrow region principles, and various asymptotic maximum and strong maximum principles for antisymmetric functions in both bounded and unbounded domains. Then we illustrate how this new method can be employed to obtain asymptotic radial symmetry and monotonicity of positive solutions in a unit ball and on the whole space. Namely, we show that no matter what the initial data are, the solutions will eventually approach to radially symmetric functions. We firmly believe that the ideas and methods introduced here can be conveniently applied to study a variety of nonlocal parabolic problems with more general operators and more general nonlinearities.


title: Random Forests applied to High Precision Photometry Analysis with Spitzer IRAC

abstract: We present a new method employing machine learning techniques for measuring astrophysical features by correcting systematics in IRAC high precision photometry using Random Forests. The main systematic in IRAC light curve data is position changes due to unavoidable telescope motions coupled with an intrapixel response function. We aim to use the large amount of publicly available calibration data for the single pixel used for this type of work (the sweet spot pixel) to make a fast, easy to use, accurate correction to science data. This correction on calibration data has the advantage of using an independent dataset instead of using the science data on itself, which has the disadvantage of including astrophysical variations. After focusing on feature engineering and hyperparameter optimization, we show that a boosted random forest model can reduce the data such that we measure the median of ten archival eclipse observations of XO-3b to be 1459 +- 200 parts per million. This is a comparable depth to the average of those in the literature done by seven different methods, however the spread in measurements is 30-100% larger than those literature values, depending on the reduction method. We also caution others attempting similar methods to check their results with the fiducial dataset of XO-3b as we were also able to find models providing initially great scores on their internal test datasets but whose results significantly underestimated the eclipse depth of that planet.


title: Raman Scattering Cross Section of Confined Carbyne

abstract: We experimentally quantify the Raman scattering from individual carbyne chains confined in double-walled carbon nanotubes. We find that the resonant differential Raman cross section of confined carbyne is on the order of $10^{-22}$ cm$^{2}\,$sr$^{-1}$ per atom, making it the strongest Raman scatterer ever reported.


title: TV-based Reconstruction of Periodic Functions

abstract: We introduce a general framework for the reconstruction of periodic multivariate functions from finitely many and possibly noisy linear measurements. The reconstruction task is formulated as a penalized convex optimization problem, taking the form of a sum between a convex data fidelity functional and a sparsity-promoting total variation based penalty involving a suitable spline-admissible regularizing operator L. In this context, we establish a periodic representer theorem, showing that the extreme-point solutions are periodic L-splines with less knots than the number of measurements. The main results are specified for the broadest classes of measurement functionals, spline-admissible operators, and convex data fidelity functionals. We exemplify our results for various regularization operators and measurement types (e.g., spatial sampling, Fourier sampling, or square-integrable functions). We also consider the reconstruction of both univariate and multivariate periodic functions.


title: Rotational transport via spontaneous symmetry breaking in vibrated disk packings

abstract: It is shown that vibrated packings of frictional disks self-organize cooperatively onto a rotational-transport state where the long-time angular velocity $\bar\omega_i$ of each disk $i$ is nonzero. Steady rotation is mediated by the spontaneous breaking of local reflection symmetry, arising when the cages in which disks are constrained by their neighbors acquire quenched disorder at large packing densities. Experiments and numerical simulation of this unexpected phenomenon show excellent agreement with each other, revealing two rotational phases as a function of excitation intensity, respectively the low-drive (LDR) and the moderate-drive (MDR) regimes. In the LDR, interdisk contacts are persistent and rotation happens due to frictional sliding. In the MDR, disks bounce against each other, still forming a solid phase. In the LDR, simple energy-dissipation arguments are provided, that support the observed dependence of the typical rotational velocity on excitation strength.


title: AutoCP: Automated Pipelines for Accurate Prediction Intervals

abstract: Successful application of machine learning models to real-world prediction problems, e.g. financial forecasting and personalized medicine, has proved to be challenging, because such settings require limiting and quantifying the uncertainty in the model predictions, i.e. providing valid and accurate prediction intervals. Conformal Prediction is a distribution-free approach to construct valid prediction intervals in finite samples. However, the prediction intervals constructed by Conformal Prediction are often (because of over-fitting, inappropriate measures of nonconformity, or other issues) overly conservative and hence inadequate for the application(s) at hand. This paper proposes an AutoML framework called Automatic Machine Learning for Conformal Prediction (AutoCP). Unlike the familiar AutoML frameworks that attempt to select the best prediction model, AutoCP constructs prediction intervals that achieve the user-specified target coverage rate while optimizing the interval length to be accurate and less conservative. We tested AutoCP on a variety of datasets and found that it significantly outperforms benchmark algorithms.


title: A Kingman-like Theorem

abstract: We provide a Kingman-like Theorem for arbitrary finite measures and a version of Birkhoff's Theorem for bounded observable. As an application, we show that Birkhoff's limit exists for some continuous observable, in an example of Bowen.


title: Representer Theorems in Banach Spaces: Minimum Norm Interpolation, Regularized Learning and Semi-Discrete Inverse Problems

abstract: Constructing or learning a function from a finite number of sampled data points (measurements) is a fundamental problem in science and engineering. This is often formulated as a minimum norm interpolation problem, regularized learning problem or, in general, a semi-discrete inverse problem, in certain functional spaces. The choice of an appropriate space is crucial for solutions of these problems. Motivated by sparse representations of the reconstructed functions such as compressed sensing and sparse learning, much of the recent research interest has been directed to considering these problems in certain Banach spaces in order to obtain their sparse solutions, which is a feasible approach to overcome challenges coming from the big data nature of most practical applications. It is the goal of this paper to provide a systematic study of the representer theorems for these problems in Banach spaces. There are a few existing results for these problems in a Banach space, with all of them regarding implicit representer theorems. We aim at obtaining explicit representer theorems based on which convenient solution methods will then be developed. For the minimum norm interpolation, the explicit representer theorems enable us to express the infimum in terms of the norm of the linear combination of the interpolation functionals. For the purpose of developing efficient computational algorithms, we establish the fixed-point equation formulation of solutions of these problems. We reveal that unlike in a Hilbert space, in general, solutions of these problems in a Banach space may not be able to be reduced to truly finite dimensional problems (with certain infinite dimensional components hidden). We demonstrate how this obstacle can be removed, reducing the original problem to a truly finite dimensional one, in the special case when the Banach space is $\ell_1(\mathbb{N})$.


title: Using public clinical trial reports to evaluate observational study methods

abstract: Observational studies are valuable for estimating the effects of various medical interventions, but are notoriously difficult to evaluate because the methods used in observational studies require many untestable assumptions. This lack of verifiability makes it difficult both to compare different observational study methods and to trust the results of any particular observational study. In this work, we propose TrialVerify, a new approach for evaluating observational study methods based on ground truth sourced from clinical trial reports. We process trial reports into a denoised collection of known causal relationships that can then be used to estimate the precision and recall of various observational study methods. We then use TrialVerify to evaluate multiple observational study methods in terms of their ability to identify the known causal relationships from a large national insurance claims dataset. We found that inverse propensity score weighting is an effective approach for accurately reproducing known causal relationships and outperforms other observational study methods. TrialVerify is made freely available for others to evaluate observational study methods.


title: Simulation Methodology for Electron Transfer in CMOS Quantum Dots

abstract: The construction of quantum computer simulators requires advanced software which can capture the most significant characteristics of the quantum behavior and quantum states of qubits in such systems. Additionally, one needs to provide valid models for the description of the interface between classical circuitry and quantum core hardware. In this study, we model electron transport in semiconductor qubits based on an advanced CMOS technology. Starting from 3D simulations, we demonstrate an order reduction and the steps necessary to obtain ordinary differential equations on probability amplitudes in a multi-particle system. We compare numerical and semi-analytical techniques concluding this paper by examining two case studies: the electron transfer through multiple quantum dots and the construction of a Hadamard gate simulated using a numerical method to solve the time-dependent Schrodinger equation and the tight-binding formalism for a time-dependent Hamiltonian.


title: Identify Influential Nodes in Online Social Network for Brand Communication

abstract: Online social networks have become incredibly popular in recent years, which prompts an increasing number of companies to promote their brands and products through social media. This paper presents an approach for identifying influential nodes in online social network for brand communication. We first construct a weighted network model for the users and their relationships extracted from the brand-related contents. We quantitatively measure the individual value of the nodes in the community from both the network structure and brand engagement aspects. Then an algorithm for identifying the influential nodes from the virtual brand community is proposed. The algorithm evaluates the importance of the nodes by their individual values as well as the individual values of their surrounding nodes. We extract and construct a virtual brand community for a specific brand from a real-life online social network as the dataset and empirically evaluate the proposed approach. The experimental results have shown that the proposed approach was able to identify influential nodes in online social network. We can get an identification result with higher ratio of verified users and user coverage by using the approach.


title: Block-matching in FPGA

abstract: Block-matching and 3D filtering (BM3D) is an image denoising algorithm that works in two similar steps. Both of these steps need to perform grouping by block-matching. We implement the block-matching in an FPGA, leveraging its ability to perform parallel computations. Our goal is to enable other researchers to use our solution in the future for real-time video denoising in video cameras that use FPGAs (such as the AXIOM Beta).


title: C^*-algebras from k group representations

abstract: We introduce certain $C^*$-algebras and $k$-graphs associated to $k$ finite dimensional unitary representations $\rho_1,...,\rho_k$ of a compact group $G$. We define a higher rank Doplicher-Roberts algebra $\mathcal{O}_{\rho_1,...,\rho_k}$, constructed from intertwiners of tensor powers of these representations. Under certain conditions, we show that this $C^*$-algebra is isomorphic to a corner in the $C^*$-algebra of a row finite rank $k$ graph $\Lambda$ with no sources. For $G$ finite and $\rho_i$ faithful of dimension at least $2$, this graph is irreducible, it has vertices $\hat{G}$ and the edges are determined by $k$ commuting matrices obtained from the character table of the group. We illustrate with some examples when $\mathcal{O}_{\rho_1,...,\rho_k}$ is simple and purely infinite, and with some $K$-theory computations.


title: Kinematic-Structure-Preserved Representation for Unsupervised 3D Human Pose Estimation

abstract: Estimation of 3D human pose from monocular image has gained considerable attention, as a key step to several human-centric applications. However, generalizability of human pose estimation models developed using supervision on large-scale in-studio datasets remains questionable, as these models often perform unsatisfactorily on unseen in-the-wild environments. Though weakly-supervised models have been proposed to address this shortcoming, performance of such models relies on availability of paired supervision on some related tasks, such as 2D pose or multi-view image pairs. In contrast, we propose a novel kinematic-structure-preserved unsupervised 3D pose estimation framework, which is not restrained by any paired or unpaired weak supervisions. Our pose estimation framework relies on a minimal set of prior knowledge that defines the underlying kinematic 3D structure, such as skeletal joint connectivity information with bone-length ratios in a fixed canonical scale. The proposed model employs three consecutive differentiable transformations named as forward-kinematics, camera-projection and spatial-map transformation. This design not only acts as a suitable bottleneck stimulating effective pose disentanglement but also yields interpretable latent pose representations avoiding training of an explicit latent embedding to pose mapper. Furthermore, devoid of unstable adversarial setup, we re-utilize the decoder to formalize an energy-based loss, which enables us to learn from in-the-wild videos, beyond laboratory settings. Comprehensive experiments demonstrate our state-of-the-art unsupervised and weakly-supervised pose estimation performance on both Human3.6M and MPI-INF-3DHP datasets. Qualitative results on unseen environments further establish our superior generalization ability.


title: Optimal Control Concerns Regarding the COVID-19 (SARS-CoV-2) Pandemic in Bahia and Santa Catarina, Brazil

abstract: The COVID-19 pandemic is the profoundest health crisis of the 21rst century. The SARS-CoV-2 virus arrived in Brazil around March, 2020 and its social and economical backlashes are catastrophic. In this paper, it is investigated how Model Predictive Control (MPC) could be used to plan appropriate social distancing policies to mitigate the pandemic effects in Bahia and Santa Catarina, two states of different regions, culture, and population demography in Brazil. In addition, the parameters of Susceptible-Infected-Recovered-Deceased (SIRD) models for these two states are identified using an optimization procedure. The control input to the process is a social isolation guideline passed to the population. Two MPC strategies are designed: a) a centralized MPC, which coordinates a single control policy for both states; and b) a decentralized strategy, for which one optimization is solved for each state. Simulation results are shown to illustrate and compare both control strategies. The framework serves as guidelines to deals with such pandemic phenomena.


title: Scalable Data Classification for Security and Privacy

abstract: Content based data classification is an open challenge. Traditional Data Loss Prevention (DLP)-like systems solve this problem by fingerprinting the data in question and monitoring endpoints for the fingerprinted data. With a large number of constantly changing data assets in Facebook, this approach is both not scalable and ineffective in discovering what data is where. This paper is about an end-to-end system built to detect sensitive semantic types within Facebook at scale and enforce data retention and access controls automatically. The approach described here is our first end-to-end privacy system that attempts to solve this problem by incorporating data signals, machine learning, and traditional fingerprinting techniques to map out and classify all data within Facebook. The described system is in production achieving a 0.9+ average F2 scores across various privacy classes while handling a large number of data assets across dozens of data stores.


title: A Model of the Fed's View on Inflation

abstract: We develop a medium-size semi-structural time series model of inflation dynamics that is consistent with the view - often expressed by central banks - that three components are important: a trend anchored by long-run expectations, a Phillips curve and temporary fluctuations in energy prices. We find that a stable long-term inflation trend and a well identified steep Phillips curve are consistent with the data, but they imply potential output declining since the new millennium and energy prices affecting headline inflation not only via the Phillips curve but also via an independent expectational channel. A high-frequency energy price cycle can be related to global factors affecting the commodity market, and often overpowers the Phillips curve thereby explaining the inflation puzzles of the last ten years.


title: Heat kernel bounds for a large class of Markov process with singular jump

abstract: Let $Z=(Z^{1}, \ldots, Z^{d})$ be the $d$-dimensional L\'evy processes where $Z^{i}$'s are independent $1$-dimensional L\'evy processes with jump kernel $J^{\phi, 1}(u,w) =|u-w|^{-1}\phi(|u-w|)^{-1}$ for $u, w\in \mathbb R$. Here $\phi$ is an increasing function with weak scaling condition of order $\underline \alpha, \overline \alpha\in (0, 2)$. Let $J(x,y) \asymp J^\phi (x,y)$ be the symmetric measurable function where \begin{align*} J^\phi(x,y):=\begin{cases} J^{\phi, 1}(x^i, y^i)\qquad&\text{ if $x^i \ne y^i$ for some $i$ and $x^j = y^j$ for all $j \ne i$}\\ 0\qquad&\text{ if $x^i \ne y^i$ for more than one index $i$.} \end{cases} \end{align*} Corresponding to the jump kernel $J$, we show the existence of non-isotropic Markov processes $X:=(X^{1}, \ldots, X^{d})$ and obtain sharp two-sided heat kernel estimates for the transition density functions.


title: A Protocol to Convert Infrastructure Data from Computer-Aided Design (CAD) to Geographic Information Systems (GIS)

abstract: While many municipalities and organizations see value in converting infrastructure data from Computer-Aided Design (CAD) to Geographic Information System (GIS) format, the process can be complex, expensive, and time-consuming. Given that municipal employees often prefer to continue performing work in both CAD and GIS, depending on the type of work required, an improved conversion process would help municipalities more fully employ GIS-based analyses. Municipalities facing budget and capacity challenges would especially benefit from an improved conversion process. With advances in GIS functionality and the promise of smart and connected cities, more emphasis is placed on the quality of data, and in this case, the potential loss of data quality from CAD to GIS formats. The goals of this article are twofold. First, to understand the common practices municipalities use to convert infrastructure CAD data to GIS and the specific challenges they face. Second, based on knowledge of those practices and challenges, this article proposes a five-step process to reduce common conversion errors and reduce the time required to correct these errors. The process is illustrated through the conversion of CAD data from the University of Illinois at Chicago (UIC) campus. The findings were validated with qualitative, semi-structured interviews conducted with GIS Analysts and Managers working in eleven municipalities across the United States who directly manage at least one of the following infrastructures: water, sanitary sewer, or stormwater sewer systems. The interviews confirmed the challenges municipalities faced with the conversion and identified solutions interviewees undertook to enable data-informed decision-making.


title: Multi-marginal optimal transport and probabilistic graphical models

abstract: We study multi-marginal optimal transport problems from a probabilistic graphical model perspective. We point out an elegant connection between the two when the underlying cost for optimal transport allows a graph structure. In particular, an entropy regularized multi-marginal optimal transport is equivalent to a Bayesian marginal inference problem for probabilistic graphical models with the additional requirement that some of the marginal distributions are specified. This relation on the one hand extends the optimal transport as well as the probabilistic graphical model theories, and on the other hand leads to fast algorithms for multi-marginal optimal transport by leveraging the well-developed algorithms in Bayesian inference. Several numerical examples are provided to highlight the results.


title: Creating QED Photon Jets with Present-Day Lasers

abstract: Large-scale, relativistic particle-in-cell simulations with quantum electrodynamics (QED) models show that high energy (1$<E_\gamma\lesssim$ 75 MeV) QED photon jets with a flux of $10^{12}$ sr$^{-1}$ can be created with present-day lasers and planar, unstructured targets. This process involves a self-forming channel in the target in response to a laser pulse focused tightly ($f$ number unity) onto the target surface. We show the self-formation of a channel to be robust to experimentally motivated variations in preplasma, angle of incidence, and laser stability, and present in simulations using historical shot data from the Texas Petawatt. We estimate that a detectable photon flux in the 10s of MeV range will require about 60 J in a 150 fs pulse.


title: Direct observation of anyonic braiding statistics at the $\nu$=1/3 fractional quantum Hall state

abstract: Utilizing an electronic Fabry-Perot interferometer in which Coulomb charging effects are suppressed, we report experimental observation of anyonic braiding statistics for the $\nu=1/3$ fractional quantum Hall state. Strong Aharonov-Bohm interference of the $\nu=1/3$ edge mode is punctuated by discrete phase slips consistent with an anyonic phase of $\theta_{anyon}=\frac{2\pi}{3}$. Our results are consistent with a recent theory of a Fabry-Perot interferometer operated in a regime in which device charging energy is small compared to the energy of formation of charged quasiparticles. Close correspondence between device operation and theoretical predictions substantiates our claim of observation of anyonic braiding.


title: Normalizing Text using Language Modelling based on Phonetics and String Similarity

abstract: Social media networks and chatting platforms often use an informal version of natural text. Adversarial spelling attacks also tend to alter the input text by modifying the characters in the text. Normalizing these texts is an essential step for various applications like language translation and text to speech synthesis where the models are trained over clean regular English language. We propose a new robust model to perform text normalization. Our system uses the BERT language model to predict the masked words that correspond to the unnormalized words. We propose two unique masking strategies that try to replace the unnormalized words in the text with their root form using a unique score based on phonetic and string similarity metrics.We use human-centric evaluations where volunteers were asked to rank the normalized text. Our strategies yield an accuracy of 86.7% and 83.2% which indicates the effectiveness of our system in dealing with text normalization.


title: Fast Learning of Graph Neural Networks with Guaranteed Generalizability: One-hidden-layer Case

abstract: Although graph neural networks (GNNs) have made great progress recently on learning from graph-structured data in practice, their theoretical guarantee on generalizability remains elusive in the literature. In this paper, we provide a theoretically-grounded generalizability analysis of GNNs with one hidden layer for both regression and binary classification problems. Under the assumption that there exists a ground-truth GNN model (with zero generalization error), the objective of GNN learning is to estimate the ground-truth GNN parameters from the training data. To achieve this objective, we propose a learning algorithm that is built on tensor initialization and accelerated gradient descent. We then show that the proposed learning algorithm converges to the ground-truth GNN model for the regression problem, and to a model sufficiently close to the ground-truth for the binary classification problem. Moreover, for both cases, the convergence rate of the proposed learning algorithm is proven to be linear and faster than the vanilla gradient descent algorithm. We further explore the relationship between the sample complexity of GNNs and their underlying graph properties. Lastly, we provide numerical experiments to demonstrate the validity of our analysis and the effectiveness of the proposed learning algorithm for GNNs.


title: The Max-Cut Decision Tree: Improving on the Accuracy and Running Time of Decision Trees

abstract: Decision trees are a widely used method for classification, both by themselves and as the building blocks of multiple different ensemble learning methods. The Max-Cut decision tree involves novel modifications to a standard, baseline model of classification decision tree construction, precisely CART Gini. One modification involves an alternative splitting metric, maximum cut, based on maximizing the distance between all pairs of observations belonging to separate classes and separate sides of the threshold value. The other modification is to select the decision feature from a linear combination of the input features constructed using Principal Component Analysis (PCA) locally at each node. Our experiments show that this node-based localized PCA with the novel splitting modification can dramatically improve classification, while also significantly decreasing computational time compared to the baseline decision tree. Moreover, our results are most significant when evaluated on data sets with higher dimensions, or more classes; which, for the example data set CIFAR-100, enable a 49% improvement in accuracy while reducing CPU time by 94%. These introduced modifications dramatically advance the capabilities of decision trees for difficult classification tasks.


title: Mod $\ell$ cohomology of some Deligne-Lusztig varieties for $\operatorname{GL}_n(q)$

abstract: In this article, we study the mod $\ell$ cohomology of some Deligne-Lusztig varieties for $\operatorname{GL}_n(q)$. We prove that the cohomology groups of these varieties are torsion-free under some conditions on the characteristic. Under the torsion-free assumption we can compute the cohomology groups explicitly and we prove that the cohomology complex satisfies partial-tilting condition, which is one of the necessary conditions in the geometric version of Brou\'{e}'s abelian defect group conjecture.


title: Tilt Map: Interactive Transitions Between Choropleth Map, Prism Map and Bar Chart in Immersive Environments

abstract: We introduce Tilt Map, a novel interaction technique for intuitively transitioning between 2D and 3D map visualisations in immersive environments. Our focus is visualising data associated with areal features on maps, for example, population density by state. Tilt Map transitions from 2D choropleth maps to 3D prism maps to 2D bar charts to overcome the limitations of each. Our paper includes two user studies. The first study compares subjects' task performance interpreting population density data using 2D choropleth maps and 3D prism maps in virtual reality (VR). We observed greater task accuracy with prism maps, but faster response times with choropleth maps. The complementarity of these views inspired our hybrid Tilt Map design. Our second study compares Tilt Map to: a side-by-side arrangement of the various views; and interactive toggling between views. The results indicate benefits for Tilt Map in user preference; and accuracy (versus side-by-side) and time (versus toggle).


title: Option Pricing: Channels, Target Zones and Sideways Markets

abstract: After a market downturn, especially in an uncertain economic environment such as the current state, there can be a relatively long period with a sideways market, where indexes, stocks, etc., move in channels with support and resistance levels. We discuss option pricing in such scenarios, in both cases of unattainable as well as attainable boundaries, and obtain closed-form option pricing formulas. Our results also apply to FX rates in target zones without interest rate pegging (USD/HKD, digital currencies, etc.).


title: Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational Wave Event S190521g

abstract: We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.


title: On Lyapunov Exponents for RNNs: Understanding Information Propagation Using Dynamical Systems Tools

abstract: Recurrent neural networks (RNNs) have been successfully applied to a variety of problems involving sequential data, but their optimization is sensitive to parameter initialization, architecture, and optimizer hyperparameters. Considering RNNs as dynamical systems, a natural way to capture stability, i.e., the growth and decay over long iterates, are the Lyapunov Exponents (LEs), which form the Lyapunov spectrum. The LEs have a bearing on stability of RNN training dynamics because forward propagation of information is related to the backward propagation of error gradients. LEs measure the asymptotic rates of expansion and contraction of nonlinear system trajectories, and generalize stability analysis to the time-varying attractors structuring the non-autonomous dynamics of data-driven RNNs. As a tool to understand and exploit stability of training dynamics, the Lyapunov spectrum fills an existing gap between prescriptive mathematical approaches of limited scope and computationally-expensive empirical approaches. To leverage this tool, we implement an efficient way to compute LEs for RNNs during training, discuss the aspects specific to standard RNN architectures driven by typical sequential datasets, and show that the Lyapunov spectrum can serve as a robust readout of training stability across hyperparameters. With this exposition-oriented contribution, we hope to draw attention to this understudied, but theoretically grounded tool for understanding training stability in RNNs.


title: Quantum flag manifold $\sigma$-models and Hermitian Ricci flow

abstract: We show that flag manifold $\sigma$-models (including $\mathbb{CP}^{n-1}$, Grassmannian models as special cases) and their deformed versions may be cast in the form of gauged bosonic Thirring/Gross-Neveu-type systems. Quantum mechanically the gauging is violated by chiral anomalies, which may be cancelled by adding fermions. We conjecture that such models are integrable and check on some examples that the trigonometrically deformed geometries satisfy the generalized Ricci flow equations.


title: Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

abstract: Membrane phase-separation is a mechanism that biological membranes often use to locally concentrate specific lipid species in order to organize diverse membrane processes. Phase separation has also been explored as a tool for the design of liposomes with heterogeneous and spatially organized surfaces. These "patchy" liposomes are promising platforms for delivery purposes, however their design and optimization through experimentation can be expensive and time-consuming. We developed a computationally efficient method based on the surface Cahn-Hilliard phase-field model to complement experimental investigations in the design of patchy liposomes. The method relies on thermodynamic considerations to set the initial state for numerical simulations. We show that our computational approach delivers not only qualitative pictures, but also accurate quantitative information about the dynamics of the membrane organization. In particular, the computational and experimental results are in excellent agreement in terms of raft area fraction, total raft perimeter over time and total number of rafts over time for two different membrane compositions (DOPC:DPPC with a 2:1 molar ratio with 20% Chol and DOPC:DPPC with a 3:1 molar ratio with 20% Chol). Thus, the computational phase-field model informed by experiments has a considerable potential to assist in the design of liposomes with spatially organized surfaces, thereby containing the cost and time required by the design process.

categories: physics.comp-ph cond-mat.soft cs.NA math.NA


title: Robust and Efficient Approximate Bayesian Computation: A Minimum Distance Approach

abstract: In many instances, the application of approximate Bayesian methods is hampered by two practical features: 1) the requirement to project the data down to low-dimensional summary, including the choice of this projection, which ultimately yields inefficient inference; 2) a possible lack of robustness to deviations from the underlying model structure. Motivated by these efficiency and robustness concerns, we construct a new Bayesian method that can deliver efficient estimators when the underlying model is well-specified, and which is simultaneously robust to certain forms of model misspecification. This new approach bypasses the calculation of summaries by considering a norm between empirical and simulated probability measures. For specific choices of the norm, we demonstrate that this approach can deliver point estimators that are as efficient as those obtained using exact Bayesian inference, while also simultaneously displaying robustness to deviations from the underlying model assumptions.


title: Novel dynamic critical phenomena induced by superfluidity and the chiral magnetic effect in Quantum Chromodynamics

abstract: In this thesis, we study the novel dynamic critical phenomena induced by superfluidity and the chiral magnetic effect in Quantum Chromodynamics (QCD), respectively. In particular, we find the new dynamic universality class of the high-density QCD critical point between the nuclear and quark superfluid phases beyond the conventional Hohnenberg and Halperin's classification. We also consider the second-order chiral phase transition under an external magnetic field and find that the inclusion of the chiral magnetic effect changes the dynamic universality class from the model E into model A within the conventional classification.


title: Identifying candidate hosts for quantum defects via data mining

abstract: Atom-like defects in solid-state hosts are promising candidates for the development of quantum information systems, but despite their importance, the host substrate/defect combinations currently under study have almost exclusively been found serendipitously. Here we systematically evaluate the suitability of host materials by applying a combined four-stage data mining and manual screening process to all entries in the Materials Project database, with literature-based experimental confirmation of band gap values. We identify 580 viable host substrates for quantum defect introduction and use in quantum information systems. While this constitutes a significant increase in the number of known and potentially viable material systems, it nonetheless represents a significant (99.54%) reduction from the total number of known inorganic phases, and the application of additional selection criteria for specific applications will reduce their number even further. The screening principles outlined may easily be applied to previously unrealized phases and other technologically important materials systems.


title: Comments on D3-Brane Holography

abstract: We revisit the idea that the quantum dynamics of open strings ending on $N$ D3-branes in the large $N$ limit can be described at large `t Hooft coupling by classical closed string theory in the background created by the D3-branes in asymptotically flat spacetime. We study the resulting thermodynamics and compute the Hagedorn temperature and other properties of the D3-brane worldvolume theory in this regime. We also consider the theory in which the D3-branes are replaced by negative branes and show that its thermodynamics is well behaved. We comment on the idea that this theory can be thought of as an irrelevant deformation of $\mathcal{N}=4$ SYM, and on its relation to $T\bar T$ deformed $CFT_2$.


title: Effect of dilute magnetism in a topological insulator

abstract: Three-dimensional topological insulators (TIs) have emerged as a unique state of quantum matter and generated enormous interests in condensed matter physics. The surfaces of a three dimensional (3D) TI are composed of a massless Dirac cone, which is characterized by the Z2 topological invariant. Introduction of magnetism on the surface of TI is essential to realize the quantum anomalous Hall effect (QAHE) and other novel magneto-electric phenomena. Here, by using a combination of first principles calculations, magneto-transport, angle-resolved photoemission spectroscopy (ARPES), and time resolved ARPES (tr-ARPES), we study the electronic properties of Gadolinium (Gd) doped Sb2Te3. Our study shows that Gd doped Sb2Te3 is a spin-orbit-induced bulk band-gap material, whose surface is characterized by a single topological surface state. We further demonstrate that introducing diluted 4f-electron magnetism into the Sb2Te3 topological insulator system by the Gd doping creates surface magnetism in this system. Our results provide a new platform to investigate the interaction between dilute magnetism and topology in doped topological materials.


title: Retiree mortality forecasting: A partial age-range or a full age-range model?

abstract: An essential input of annuity pricing is the future retiree mortality. From observed age-specific mortality data, modeling and forecasting can be taken place in two routes. On the one hand, we can first truncate the available data to retiree ages and then produce mortality forecasts based on a partial age-range model. On the other hand, with all available data, we can first apply a full age-range model to produce forecasts and then truncate the mortality forecasts to retiree ages. We investigate the difference in modeling the logarithmic transformation of the central mortality rates between a partial age-range and a full age-range model, using data from mainly developed countries in the Human Mortality Database (2020). By evaluating and comparing the short-term point and interval forecast accuracies, we recommend the first strategy by truncating all available data to retiree ages and then produce mortality forecasts. However, when considering the long-term forecasts, it is unclear which strategy is better since it is more difficult to find a model and parameters that are optimal. This is a disadvantage of using methods based on time series extrapolation for long-term forecasting. Instead, an expectation approach, in which experts set a future target, could be considered, noting that this method has also had limited success in the past.


title: Coherent Solutions and Transition to Turbulence in Two-Dimensional Rayleigh-B\'{e}nard Convection

abstract: For two-dimensional Rayleigh-B\'{e}nard convection, classes of unstable, steady solutions were previously computed using numerical continuation (Waleffe, 2015; Sondak, 2015). The `primary' steady solution bifurcates from the conduction state at $Ra \approx 1708$, and has a characteristic aspect ratio (length/height) of approximately $2$. The primary solution corresponds to one pair of counterclockwise-clockwise convection rolls with a temperature updraft in between and an adjacent downdraft on the sides. By adjusting the horizontal length of the domain, (Waleffe, 2015; Sondak, 2015) also found steady, maximal heat transport solutions, with characteristic aspect ratio less than $2$ and decreasing with increasing $Ra$. Compared to the primary solutions, optimal heat transport solutions have modifications to boundary layer thickness, the horizontal length scale of the plume, and the structure of the downdrafts. The current study establishes a direct link between these (unstable) steady solutions and transition to turbulence for $Pr = 7$ and $Pr = 100$. For transitional values of $Ra$, the primary and optimal heat transport solutions both appear prominently in appropriately-sized sub-fields of the time-evolving temperature fields. For $Ra$ beyond transitional, our data analysis shows persistence of the primary solution for $Pr = 7$, while the optimal heat transport solutions are more easily detectable for $Pr = 100$. In both cases $Pr = 7$ and $Pr = 100$, the relative prevalence of primary and optimal solutions is consistent with the $Nu$ vs. $Ra$ scalings for the numerical data and the steady solutions.


title: Rare-earth-mediated opto-mechanical system in the reversed dissipation regime

abstract: Strain-mediated interaction between phonons and telecom photons is demonstrated using excited states of erbium ions embedded in a mechanical resonator. Owing to the extremely long-lived nature of rare-earth ions, the dissipation rate of the optical resonance falls below that of the mechanical one. Thus, a reversed dissipation regime is achieved in the optical frequency region. We experimentally demonstrate an opto-mechanical coupling rate 21.7 Hz, and numerically reveal that the interaction causes stimulated excitation of erbium ions. Numerical analyses further indicate the possibility of g_0 exceeding the dissipation rates of erbium and mechanical systems, thereby leading to single-photon strong coupling. This strain-mediated interaction moreover involves the spin degree of freedom, and has a potential to be extended to highly-coherent opto-electro-mechanical hybrid systems in the reversed dissipation regime.


title: Conjecture A and $\mu$-invariant for Selmer groups of supersingular elliptic curves

abstract: Let $p$ be an odd prime and let $E$ be an elliptic curve defined over a number field $F$ with good reduction at primes above $p$. In this survey article, we give an overview of some of the important results proven for the fine Selmer group and the signed Selmer groups over cyclotomic towers as well as the signed Selmer groups over $\mathbb{Z}_p^2$-extensions of an imaginary quadratic field where $p$ splits completely. We only discuss the algebraic aspects of these objects through Iwasawa theory. We also attempt to give some of the recent results implying the vanishing of the $\mu$-invariant under the hypothesis of Conjecture A. Moreover, we draw an analogy between the classical Selmer group in the ordinary reduction case and that of the signed Selmer groups of Kobayashi in the supersingular reduction case. We highlight properties of signed Selmer groups (when $E$ has good supersingular reduction) which are completely analogous to the classical Selmer group (when $E$ has good ordinary reduction). In this survey paper, we do not present any proofs, however we have tried to give references of the discussed results for the interested reader.


title: Explainable CNN-attention Networks (C-Attention Network) for Automated Detection of Alzheimer's Disease

abstract: In this work, we propose three explainable deep learning architectures to automatically detect patients with Alzheimer`s disease based on their language abilities. The architectures use: (1) only the part-of-speech features; (2) only language embedding features and (3) both of these feature classes via a unified architecture. We use self-attention mechanisms and interpretable 1-dimensional ConvolutionalNeural Network (CNN) to generate two types of explanations of the model`s action: intra-class explanation and inter-class explanation. The inter-class explanation captures the relative importance of each of the different features in that class, while the inter-class explanation captures the relative importance between the classes. Note that although we have considered two classes of features in this paper, the architecture is easily expandable to more classes because of its modularity. Extensive experimentation and comparison with several recent models show that our method outperforms these methods with an accuracy of 92.2% and F1 score of 0.952on the DementiaBank dataset while being able to generate explanations. We show by examples, how to generate these explanations using attention values.


title: Digital quantum simulation framework for energy transport in an open quantum system

abstract: Quantum effects such as the environment assisted quantum transport (ENAQT) displayed in photosynthetic Fenna-Mathews-Olson (FMO) complex has been simulated on analog quantum simulators. Digital quantum simulations offer greater universality and flexibility over analog simulations. However, digital quantum simulations of open quantum systems face a theoretical challenge; one does not know the solutions of the continuous time master equation for developing quantum gate operators. We give a theoretical framework for digital quantum simulation of ENAQT by introducing new quantum evolution operators. We develop the dynamical equation for the operators and prove that it is an analytical solution of the master equation. As an example, using the dynamical equations, we simulate the FMO complex in the digital setting, reproducing theoretical and experimental evidence of the dynamics. The framework gives an optimal method for {quantum circuit} implementation, giving a log reduction in complexity over known methods. The generic framework can be extrapolated to study other open quantum systems.