title: On the Dynamical Origins of Retrograde Jupiter Trojans and their Connection to High-Inclination TNOs

abstract: Over the course of the last decade, observations of highly-inclined (orbital inclination i > 60{\deg}) Trans-Neptunian Objects (TNOs) have posed an important challenge to current models of solar system formation (Levison et al. 2008; Nesvorn\'y 2015). These remarkable minor planets necessitate the presence of a distant reservoir of strongly-out-of-plane TNOs, which itself requires some dynamical production mechanism (Gladman et al. 2009; Gomes et al. 2015; Batygin and Brown 2016). A notable recent addition to the census of high-i minor bodies in the solar system is the retrograde asteroid 514107 Ka'epaoka'awela, which currently occupies a 1:-1 mean motion resonance with Jupiter at i = 163{\deg} (Wiegert et al. 2017). In this work, we delineate a direct connection between retrograde Jupiter Trojans and high-i Centaurs. First, we back-propagate a large sample of clones of Ka'epaoka'awela for 100 Ma numerically, and demonstrate that long-term stable clones tend to decrease their inclination steadily until it concentrates between 90{\deg} and 135{\deg}, while their eccentricity and semi-major axis increase, placing many of them firmly into the trans-Neptunian domain. Importantly, the clones show significant overlap with the synthetic high-i Centaurs generated in Planet 9 studies (Batygin et al. 2019), and hint at the existence of a relatively prominent, steady-state population of minor bodies occupying polar trans-Saturnian orbits. Second, through direct numerical forward-modeling, we delineate the dynamical pathway through which conventional members of the Kuiper Belt's scattered disk population can become retrograde Jovian Trojan resonators in presence of Planet 9.


title: Extending fibre nonlinear interference power modelling to account for general dual-polarisation 4D modulation formats

abstract: In optical communications, four-dimensional (4D) modulation formats encode information onto the quadrature components of two arbitrary orthogonal states of polarisation of the optical field. These formats have recently regained attention due their potential power efficiency, nonlinearity tolerance, and ultimately to their still unexplored shaping gains. As in the fibre-optic channel the shaping gain is closely related to the nonlinearity tolerance of a given modulation format, predicting the effect of nonlinearity is key to effectively optimise the transmitted constellation. Many analytical models available in the optical communication literature allow, within a first-order perturbation framework, the computation of the average power of the nonlinear interference (NLI) accumulated in coherent fibre-optic transmission systems. However, all current models only operate under the assumption of a transmitted polarisation-multiplexed, two-dimensional (PM-2D) modulation format. PM-2D formats represent a limited subset of the possible dual-polarisation 4D formats, namely, only those where data transmitted on each polarisation channel are mutually independent and identically distributed. This document presents a step-by-step mathematical derivation of the extension of existing NLI models to the class of arbitrary dual-polarisation 4D modulation formats. In particular, the methodology adopted follows the one of the popular enhanced Gaussian noise model, albeit dropping most assumptions on the geometry and statistic of the transmitted 4D modulation format. The resulting expressions show that, whilst in the PM-2D case the NLI power depends only on different statistical high-order moments of each polarisation component, for a general 4D constellation also several others cross-polarisation correlations need to be taken into account.


title: Machine learning applied in the multi-scale 3D stress modelling

abstract: This paper proposes a methodology to estimate stress in the subsurface by a hybrid method combining finite element modeling and neural networks. This methodology exploits the idea of obtaining a multi-frequency solution in the numerical modeling of systems whose behavior involves a wide span of length scales. One low-frequency solution is obtained via inexpensive finite element modeling at a coarse scale. The second solution provides the fine-grained details introduced by the heterogeneity of the free parameters at the fine scale. This high-frequency solution is estimated via neural networks -trained with partial solutions obtained in high-resolution finite-element models. When the coarse finite element solutions are combined with the neural network estimates, the results are within a 2\% error of the results that would be computed with high-resolution finite element models. This paper discusses the benefits and drawbacks of the method and illustrates their applicability via a worked example.


title: Deep Networks and the Multiple Manifold Problem

abstract: We study the multiple manifold problem, a binary classification task modeled on applications in machine vision, in which a deep fully-connected neural network is trained to separate two low-dimensional submanifolds of the unit sphere. We provide an analysis of the one-dimensional case, proving for a simple manifold configuration that when the network depth $L$ is large relative to certain geometric and statistical properties of the data, the network width $n$ grows as a sufficiently large polynomial in $L$, and the number of i.i.d. samples from the manifolds is polynomial in $L$, randomly-initialized gradient descent rapidly learns to classify the two manifolds perfectly with high probability. Our analysis demonstrates concrete benefits of depth and width in the context of a practically-motivated model problem: the depth acts as a fitting resource, with larger depths corresponding to smoother networks that can more readily separate the class manifolds, and the width acts as a statistical resource, enabling concentration of the randomly-initialized network and its gradients. The argument centers around the neural tangent kernel and its role in the nonasymptotic analysis of training overparameterized neural networks; to this literature, we contribute essentially optimal rates of concentration for the neural tangent kernel of deep fully-connected networks, requiring width $n \gtrsim L\,\mathrm{poly}(d_0)$ to achieve uniform concentration of the initial kernel over a $d_0$-dimensional submanifold of the unit sphere $\mathbb{S}^{n_0-1}$, and a nonasymptotic framework for establishing generalization of networks trained in the NTK regime with structured data. The proof makes heavy use of martingale concentration to optimally treat statistical dependencies across layers of the initial random network. This approach should be of use in establishing similar results for other network architectures.


title: Magnetism and stability of all primitive stacking patterns in bilayer chromium trihalides

abstract: Chromium trihalides, CrX$_3$ (with X = Cl, Br, I), are a family of layered magnetic materials that can be easily exfoliated to provide ferromagnetic monolayers. When two layers are stacked together to form a bilayer the interlayer exchange coupling can be either ferromagnetic or antiferromagnetic depending on the stacking sequence. Here we combine crystallographic arguments based on the close-packing condition with first-principles simulations to enumerate all possible stacking patterns in CrX$_3$ bilayers that preserve the spatial periodicity of each layer. We recover all configurations observed in bulk crystals and disclose stacking sequences with no bulk counterpart where the two layers have opposite chirality. Stacking sequences are ranked according to their relative stability and a preferential interlayer magnetic ordering is assigned to each of them. Simulations provide a consistent picture to frame all current experimental observations on bulk and exfoliated CrX$_3$ crystals, with interesting implications for future measurements, including synthetic bilayers with non-standard stacking patterns.


title: Photonic reservoir computer based on frequency multiplexing

abstract: Reservoir computing is a brain inspired approach for information processing, well suited to analogue implementations. We report a photonic implementation of a reservoir computer that exploits frequency domain multiplexing to encode neuron states. The system processes 25 comb lines simultaneously (i.e. 25 neurons), at a rate of 20 MHz. We illustrate performances on two standard benchmark tasks: channel equalization and time series forecasting. We also demonstrate that frequency multiplexing allows output weights to be implemented in the optical domain, through optical attenuation. We discuss the perspectives for high speed high performance low footprint implementations.


title: Shifted functors of linear representations are semisimple

abstract: We prove that, for any fields $k$ and $\mathbb{F}$ of characteristic $0$ and any finite group $T$, the category of modules over the shifted Green biset functor $(kR_{\mathbb{F}})_T$ is semisimple.


title: SOAR: Simultaneous Or of And Rules for Classification of Positive & Negative Classes

abstract: Algorithmic decision making has proliferated and now impacts our daily lives in both mundane and consequential ways. Machine learning practitioners make use of a myriad of algorithms for predictive models in applications as diverse as movie recommendations, medical diagnoses, and parole recommendations without delving into the reasons driving specific predictive decisions. Machine learning algorithms in such applications are often chosen for their superior performance, however popular choices such as random forest and deep neural networks fail to provide an interpretable understanding of the predictive model. In recent years, rule-based algorithms have been used to address this issue. Wang et al. (2017) presented an or-of-and (disjunctive normal form) based classification technique that allows for classification rule mining of a single class in a binary classification; this method is also shown to perform comparably to other modern algorithms. In this work, we extend this idea to provide classification rules for both classes simultaneously. That is, we provide a distinct set of rules for both positive and negative classes. In describing this approach, we also present a novel and complete taxonomy of classifications that clearly capture and quantify the inherent ambiguity in noisy binary classifications in the real world. We show that this approach leads to a more granular formulation of the likelihood model and a simulated-annealing based optimization achieves classification performance competitive with comparable techniques. We apply our method to synthetic as well as real world data sets to compare with other related methods that demonstrate the utility of our proposal.


title: What Do We Actually Learn from Evaluations in the "Heroic Era" of Visualization?

abstract: We often point to the relative increase in the amount and sophistication of evaluations of visualization systems versus the earliest days of the field as evidence that we are maturing as a field. I am not so convinced. In particular, I feel that evaluations of visualizations, as they are ordinarily performed in the field or asked for by reviewers, fail to tell us very much that is useful or transferable about visualization systems, regardless of the statistical rigor or ecological validity of the evaluation. Through a series of thought experiments, I show how our current conceptions of visualization evaluations can be incomplete, capricious, or useless for the goal of furthering the field, more in line with the "heroic age" of medical science than the rigorous evidence-based field we might aspire to be. I conclude by suggesting that our models for designing evaluations, and our priorities as a field, should be revisited.


title: Modeling Cell Populations Measured By Flow Cytometry With Covariates Using Sparse Mixture of Regressions

abstract: The ocean is filled with microscopic microalgae called phytoplankton, which together are responsible for as much photosynthesis as all plants on land combined. Our ability to predict their response to the warming ocean relies on understanding how the dynamics of phytoplankton populations is influenced by changes in environmental conditions. One powerful technique to study the dynamics of phytoplankton is flow cytometry, which measures the optical properties of thousands of individual cells per second. Today, oceanographers are able to collect flow cytometry data in real-time onboard a moving ship, providing them with fine-scale resolution of the distribution of phytoplankton across thousands of kilometers. One of the current challenges is to understand how these small and large scale variations relate to environmental conditions, such as nutrient availability, temperature, light and ocean currents. In this paper, we propose a novel sparse mixture of multivariate regressions model to estimate the time-varying phytoplankton subpopulations while simultaneously identifying the specific environmental covariates that are predictive of the observed changes to these subpopulations. We demonstrate the usefulness and interpretability of the approach using both synthetic data and real observations collected on an oceanographic cruise conducted in the north-east Pacific in the spring of 2017.


title: On the propagation of regularity for solutions of the Zakharov-Kuznetsov equation

abstract: In this work, we study some special properties of smoothness concerning to the initial value problem associated with the Zakharov-Kuznetsov-(ZK) equation in the $n-$ dimensional setting, $n\geq 2.$ It is known that the solutions of the ZK equation in the $2d$ and $3d$ cases verify special regularity properties. More precisely, the regularity of the initial data on a family of half-spaces propagates with infinite speed. Our objective in this work is to extend this analysis to the case in that the regularity of the initial data is measured on a fractional scale. To describe this phenomenon we present new localization formulas that allow us to portray the regularity of the solution on a certain class of subsets of the euclidean space.


title: The square and add Markov chain

abstract: Squaring and adding $\pm 1$ mod p generates a curiously intractable random walk. A similar process over the finite field $\mathbf{F}_q$ (with $q=2^d$) leads to novel connections between elementary Galois theory and probability.


title: Temporal Action Localization with Variance-Aware Networks

abstract: This work addresses the problem of temporal action localization with Variance-Aware Networks (VAN), i.e., DNNs that use second-order statistics in the input and/or the output of regression tasks. We first propose a network (VANp) that when presented with the second-order statistics of the input, i.e., each sample has a mean and a variance, it propagates the mean and the variance throughout the network to deliver outputs with second order statistics. In this framework, both the input and the output could be interpreted as Gaussians. To do so, we derive differentiable analytic solutions, or reasonable approximations, to propagate across commonly used NN layers. To train the network, we define a differentiable loss based on the KL-divergence between the predicted Gaussian and a Gaussian around the ground truth action borders, and use standard back-propagation. Importantly, the variances propagation in VANp does not require any additional parameters, and during testing, does not require any additional computations either. In action localization, the means and the variances of the input are computed at pooling operations, that are typically used to bring arbitrarily long videos to a vector with fixed dimensions. Second, we propose two alternative formulations that augment the first (respectively, the last) layer of a regression network with additional parameters so as to take in the input (respectively, predict in the output) both means and variances. Results in the action localization problem show that the incorporation of second order statistics improves over the baseline network, and that VANp surpasses the accuracy of virtually all other two-stage networks without involving any additional parameters.


title: Some Remarks on H-stability of syzygy bundle on algebraic surface

abstract: Let $L$ be a globally generated line bundle over a smooth irreducible complex projective surface $X$. The syzygy bundle $M_{L}$ is the kernel of the evaluation map $H^0(L)\otimes\mathcal O_X\to L$. We prove the $L$-stability of $M_L$ for Hirzebruch surfaces, del Pezzo surfaces and Enriques surfaces. The $(-K_X)$-stability of syzygy bundles $M_L$ over del Pezzo surfaces is also obtained.


title: Differentiating a Tensor Language

abstract: How does one compile derivatives of tensor programs, such that the resulting code is purely functional (hence easier to optimize and parallelize) and provably efficient relative to the original program? We show that naively differentiating tensor code---as done in popular systems like Tensorflow and PyTorch---can cause asymptotic slowdowns in pathological cases, violating the Cheap Gradients Principle. However, all existing automatic differentiation methods that guarantee this principle (for variable size data) do so by relying on += mutation through aliases/pointers---which complicates downstream optimization. We provide the first purely functional, provably efficient, adjoint/reverse-mode derivatives of array/tensor code by explicitly accounting for sparsity. We do this by focusing on the indicator function from Iverson's APL. We also introduce a new "Tensor SSA" normal form and a new derivation of reverse-mode automatic differentiation based on the universal property of inner-products.


title: The Impact of Indirect Machine Translation on Sentiment Classification

abstract: Sentiment classification has been crucial for many natural language processing (NLP) applications, such as the analysis of movie reviews, tweets, or customer feedback. A sufficiently large amount of data is required to build a robust sentiment classification system. However, such resources are not always available for all domains or for all languages. In this work, we propose employing a machine translation (MT) system to translate customer feedback into another language to investigate in which cases translated sentences can have a positive or negative impact on an automatic sentiment classifier. Furthermore, as performing a direct translation is not always possible, we explore the performance of automatic classifiers on sentences that have been translated using a pivot MT system. We conduct several experiments using the above approaches to analyse the performance of our proposed sentiment classification system and discuss the advantages and drawbacks of classifying translated sentences.


title: Say "Sul Sul!" to SimSim, A Sims-Inspired Platform for Sandbox Game AI

abstract: This paper proposes environment design in the life simulation game The Sims as a novel platform and challenge for testing divergent search algorithms. In this domain, which includes a minimal viability criterion, the goal is to furnish a house with objects that satisfy the physical needs of a simulated agent. Importantly, the large number of objects available to the player (whether human or automated) affords a wide variety of solutions to the underlying design problem. Empirical studies in a novel open source simulator called SimSim investigate the ability of novelty-based evolutionary algorithms to effectively generate viable environment designs.


title: Bipolaron formed through electron-hole excitation

abstract: We investigate the electronic properties and electron correlations of the bipolaron formed by the electron-hole excitations in the presence of Yukawa-type coupling (between nonrelativistic fermions) in three spatial dimension. The electron-hole excitation, which is necessary to the formation of bipolaron, leads to imaginary particle-hole order parameter, and provide finite boson field mass to the single-polaron dispersion in a broken-symmetry phase. We found that the bipolaron exhibits fermi-liquid features as long as the long-range strong interaction is suppressed, and it behave differently compared to the single-polaron. The bosonic momentum determines the mass of boson field propagator and the gap function, and it also related to the self-energies and the single-particle Green's functions. The Thouless criterion is also used during the calculation of gap equation at critical temperature (which become lower in weak-coupling regime), which corresponds to the pole (instability) of the pair propagator in zero center-of-mass freamwork. The mean field term and the bosonic fluctuation-induced contribution to free energy of bipolaron are also studied.


title: Topological Terms and Diffeomorphism Anomalies in Fluid Dynamics and Sigma Models

abstract: The requirement of diffeomorphism symmetry for the target space can lead to anomalous commutators for the energy-momentum tensor for sigma models and for fluid dynamics, if certain topological terms are added to the action. We analyze several examples . A particular topological term is shown to lead to the known effective hydrodynamics of a dense collection of vortices, i.e. the vortex fluid theory in 2+1 dimensions. The possibility of a similar vortex fluid in 3+1 dimensions, as well as a fluid of knots and links, with possible extended diffeomorphism algebras is also discussed.


title: Matter-gravity coupling for fuzzy geometry and the Landau-Hall problem

abstract: We consider a set of physical degrees of freedom coupled to a finite-dimensional Hilbert space, which may be taken as modeling a fuzzy space or as the lowest Landau level of a Landau-Hall problem. These may be viewed as matter fields on a fuzzy space. Sequentially generalizing to arbitrary backgrounds, we argue that the effective action is given by the Chern-Simons form associated with the Dirac index density (with gauge and gravitational fields), with an abelian gauge field shifted by the Poincar\'e-Cartan form for matter dynamics. The result is an action for matter fields where the Lagrangian is integrated with a density which is a specific polynomial in the curvatures.


title: Student Teamwork on Programming Projects: What can GitHub logs show us?

abstract: Teamwork, often mediated by version control systems such as Git and Apache Subversion (SVN), is central to professional programming. As a consequence, many colleges are incorporating both collaboration and online development environments into their curricula even in introductory courses. In this research, we collected GitHub logs from two programming projects in two offerings of a CS2 Java programming course for computer science majors. Students worked in pairs for both projects (one optional, the other mandatory) in each year. We used the students' GitHub history to classify the student teams into three groups, collaborative, cooperative, or solo-submit, based on the division of labor. We then calculated different metrics for students' teamwork including the total number and the average number of commits in different parts of the projects and used these metrics to predict the students' teamwork style. Our findings show that we can identify the students' teamwork style automatically from their submission logs. This work helps us to better understand novices' habits while using version control systems. These habits can identify the harmful working styles among them and might lead to the development of automatic scaffolds for teamwork and peer support in the future.


title: $L^{2}$-hard Lefschetz complete symplectic manifolds

abstract: For a complete symplectic manifold $M^{2n}$, we define the $L^{2}$-hard Lefschetz property on $M^{2n}$. We also prove that the complete symplectic manifold $M^{2n}$ satisfies $L^{2}$-hard Lefschetz property if and only if every class of $L^{2}$-harmonic forms contains a $L^{2}$ symplectic harmonic form. As an application, we get if $M^{2n}$ is a closed symplectic parabolic manifold which satisfies the hard Lefschetz property, then its Euler characteristic satisfies the inequality $(-1)^{n}\chi(M^{2n})\geq0$.


title: Electromagnetic radiation and the self torque of an oscillating magnetic dipole

abstract: A uniformly-charged spherical shell of radius $R$, mass $m$, and total electrical charge $q$, having an oscillatory angular velocity $\Omega(t)$ around a fixed axis, is a model for a magnetic dipole that radiates an electromagnetic field into its surrounding free space at a fixed oscillation frequency $\omega$. An exact solution of the Maxwell-Lorentz equations of classical electrodynamics yields the self-torque of radiation resistance acting on the spherical shell as a function of $R$, $q$, and $\omega$. Invoking the Newtonian equation of motion for the shell, we relate its angular velocity $\Omega(t)$ to an externally applied torque, and proceed to examine the response of the magnetic dipole to an impulsive torque applied at a given instant of time, say, $t=0$. The impulse response of the dipole is found to be causal down to extremely small values of $R$ (i.e., as $R \to 0$) so long as the exact expression of the self-torque is used in the dynamical equation of motion of the spherical shell.


title: Delay-Efficient and Reliable Data Relaying in Ultra Dense Networks using Rateless Codes

abstract: We investigate the problem of delay-efficient and reliable data delivery in ultra-dense networks (UDNs) that constitute macro base stations (MBSs), small base stations (SBSs), and mobile users. Considering a two-hop data delivery system, we propose a partial decode-and-forward (PDF) relaying strategy together with a simple and intuitive amicable encoding scheme for rateless codes to significantly improve user experience in terms of end-to-end delay. Simulation results verify that our amicable encoding scheme is efficient in improving the intermediate performance of rateless codes. It also verifies that our proposed PDF significantly improves the performance of the decode-and-forward (DF) strategy, and that PDF is much more robust against channel degradation. Overall, the proposed strategy and encoding scheme are efficient towards delay-sensitive data delivery in the UDN scenarios.


title: The quasi-particle picture and its breakdown after local quenches: mutual information, negativity, and reflected entropy

abstract: We study the dynamics of (R\'enyi) mutual information, logarithmic negativity, and (R\'enyi) reflected entropy after exciting the ground state by a local operator. Together with recent results from Ref. [1], we are able to conjecture a close-knit structure between the three quantities that emerges in states excited above the vacuum, including both local and global quantum quenches. This structure intimately depends on the chaoticity of the theory i.e. there exist distinct sets of equivalences for integrable and chaotic theories. For rational conformal field theories (RCFT), we find all quantities to compute the quantum dimension of the primary operator inserted. In contrast, we find the correlation measures to grow (logarithmically) without bound in all $c > 1$ conformal field theories with a finite twist gap. In comparing the calculations in the two classes of theories, we are able to identify the dynamical mechanism for the breakdown of the quasi-particle picture in 2D conformal field theories. Intriguingly, we also find preliminary evidence that our general lessons apply to quantum systems considerably distinct from conformal field theories, such as integrable and chaotic spin chains, suggesting a universality of entanglement dynamics in non-equilibrium systems.


title: General theory of lifting spaces

abstract: In his classical textbook on algebraic topology Edwin Spanier developed the theory of covering spaces within a more general framework of lifting spaces (i.e., Hurewicz fibrations with unique path-lifting property). Among other, Spanier proved that for every space $X$ there exists a universal lifting space, which however need not be simply connected, unless the base space $X$ is semi-locally simply connected. The question on what exactly is the fundamental group of the universal space was left unanswered. The main source of lifting spaces are inverse limits of covering spaces over $X$, or more generally, over some inverse system of spaces converging to $X$. Every metric space $X$ can be obtained as a limit of an inverse system of polyhedra, and so inverse limits of covering spaces over the system yield lifting spaces over $X$. They are related to the geometry (in particular the fundamental group) of $X$ in a similar way as the covering spaces over polyhedra are related to the fundamental group of their base. Thus lifting spaces appear as a natural replacement for the concept of covering spaces over base spaces with bad local properties. In this paper we develop a general theory of lifting spaces and prove that they are preserved by products, inverse limits and other important constructions. We show that maps from $X$ to polyhedra give rise to coverings over $X$ and use that to prove that for a connected, locally path connected and paracompact $X$, the fundamental group of the above-mentioned Spanier's universal space is precisely the intersection of all Spanier groups associated to open covers of $X$, and that the later coincides with the shape kernel of $X$. Furthermore, we examine in more detail lifting spaces over $X$ that arise as inverse limits of coverings over some approximations of $X$.


title: Classifying minimal vanishing sums of roots of unity

abstract: A vanishing sum of roots of unity is called minimal if no proper, nonempty sub-sum of it vanishes. This paper classifies all minimal vanishing sums of roots of unity of weight at most 16 by hand and then develops and codes an algorithm to extend the classification further, up to weight no greater than 21.


title: Deep Neural Network for 3D Surface Segmentation based on Contour Tree Hierarchy

abstract: Given a 3D surface defined by an elevation function on a 2D grid as well as non-spatial features observed at each pixel, the problem of surface segmentation aims to classify pixels into contiguous classes based on both non-spatial features and surface topology. The problem has important applications in hydrology, planetary science, and biochemistry but is uniquely challenging for several reasons. First, the spatial extent of class segments follows surface contours in the topological space, regardless of their spatial shapes and directions. Second, the topological structure exists in multiple spatial scales based on different surface resolutions. Existing widely successful deep learning models for image segmentation are often not applicable due to their reliance on convolution and pooling operations to learn regular structural patterns on a grid. In contrast, we propose to represent surface topological structure by a contour tree skeleton, which is a polytree capturing the evolution of surface contours at different elevation levels. We further design a graph neural network based on the contour tree hierarchy to model surface topological structure at different spatial scales. Experimental evaluations based on real-world hydrological datasets show that our model outperforms several baseline methods in classification accuracy.


title: Fluorescence-Detected Circular Dichroism of a Chiral Molecular Monolayer with Dielectric Metasurfaces

abstract: Strong enhancement of molecular circular dichroism has the potential to enable efficient asymmetric photolysis, a method of chiral separation that has conventionally been impeded by insufficient yield and low enantiomeric excess. Here, we study experimentally how predicted enhancements in optical chirality density near resonant silicon nanodisks boost circular dichroism. We use fluorescence-detected circular dichroism spectroscopy to measure indirectly the differential absorption of circularly polarized light by a monolayer of optically active molecules functionalized to silicon nanodisk arrays. Importantly, the molecules and nanodisk antennas have spectrally-coincident resonances, and our fluorescence technique allows us to deconvolute absorption in the nanodisks from the molecules. We find that enhanced fluorescence-detected circular dichroism signals depend on nanophotonic resonances in good agreement with simulated differential absorption and optical chirality density, while no signal is detected from molecules adsorbed on featureless silicon surfaces. These results verify the potential of nanophotonic platforms to be used for asymmetric photolysis with lower energy requirements


title: Rotating Neutron Stars in F(R) Gravity with Axions

abstract: We investigate equilibrium configurations of uniformly rotating neutron stars in $R^2$ gravity with axion scalar field for GM1 equation of state (EoS) for nuclear matter. The mass-radius diagram, mass-central energy density are presented for some frequencies in comparison with static stars. We also compute equatorial and polar radii and moment of inertia for stars. For axion field $\phi$ the coupling in the form $\sim R^2\phi$ is assumed. Several interesting results follow from our consideration. Maximal possible star mass with given EoS increases due to the contribution of coupling term. We discovered the possibility to increase maximal frequency of the rotation in comparison with General Relativity. As a consequence the lower bound on mass of the fast rotating stars decreases. For frequency $f=700$ Hz neutron stars with masses $\sim M_\odot$ can exist for some choice of parameters (in General Relativity for same EoS this limit is around $1.2 M_{\odot}$). Another feature of our solutions is relatively small increase of stars radii for high frequencies in comparison with static case. Thus, eventually the new class of neutron stars in $R^2$ gravity with axions is discovered namely fast rotating compact stars with intermediate masses.


title: How to construct a upper triangular matrix that satisfy the quadratic polynomial equation with different roots

abstract: Let $R$ be an associative ring with identity $1$. We describe all matrices in $T_n(R)$ the ring of $n\times n$ upper triangular matrices over $R$ ($n\in \mathbb{N}$), and $T_{\infty}(R)$ the ring of infinite upper triangular matrices over $R$, satisfying the quadratic polynomial equation $x^2-rx+s=0$. For such propose we assume that the above polynomial have two different roots in $R$. Moreover, in the case that $R$ in finite, we compute the number of all matrices to solves the matrix equation $A^2-rA+sI=0,$ where $I$ is the identity matrix.


title: Independent Vector Analysis with Deep Neural Network Source Priors

abstract: This paper studies the density priors for independent vector analysis (IVA) with convolutive speech mixture separation as the exemplary application. Most existing source priors for IVA are too simplified to capture the fine structures of speeches. Here, we first time show that it is possible to efficiently estimate the derivative of speech density with universal approximators like deep neural networks (DNN) by optimizing certain proxy separation related performance indices. Experimental results suggest that the resultant neural network density priors consistently outperform previous ones in convergence speed for online implementation and signal-to-interference ratio (SIR) for batch implementation.


title: Structure exploiting methods for fast uncertainty quantification in multiphase flow through heterogeneous media

abstract: We present a computational framework for dimension reduction and surrogate modeling to accelerate uncertainty quantification in computationally intensive models with high-dimensional inputs and function-valued outputs. Our driving application is multiphase flow in saturated-unsaturated porous media in the context of radioactive waste storage. For fast input dimension reduction, we utilize an approximate global sensitivity measure, for function-value outputs, motivated by ideas from the active subspace methods. The proposed approach does not require expensive gradient computations. We generate an efficient surrogate model by combining a truncated Karhunen-Lo\'{e}ve (KL) expansion of the output with polynomial chaos expansions, for the output KL modes, constructed in the reduced parameter space. We demonstrate the effectiveness of the proposed surrogate modeling approach with a comprehensive set of numerical experiments, where we consider a number of function-valued (temporally or spatially distributed) QoIs.


title: Formula to Determine the Countries Equilibrium Exchange Rate With the Dollar and Proposal for a Second Bretton Woods Conference

abstract: This paper presents the way in which can be determined the exchange rates that simultaneously balance the trade balances of all countries that trade with each other within a common market. A mathematical synthesis between the theory of comparative advantages of Ricardo and Mill and the New Theory on International Trade of Paul Krugman is also presented in this paper. This mathematical synthesis shows that these theories are complementary. Also presented in this paper is a proposal to organize a common market of the American hemisphere. This economic alliance would allow to establish a political alliance for the common defense of the entire American hemisphere. The formula we developed in this paper to determine the exchange rates of the countries solves the problem that Mexico, Canada, Europe, Japan and China currently experience with the United States in relation to the deficits and surpluses in the trade balance of the countries and their consequent impediment so that stable growth of international trade can be achieved.


title: Linking Machine Learning with Multiscale Numerics: Data-Driven Discovery of Homogenized Equations

abstract: The data-driven discovery of partial differential equations (PDEs) consistent with spatiotemporal data is experiencing a rebirth in machine learning research. Training deep neural networks to learn such data-driven partial differential operators requires extensive spatiotemporal data. For learning coarse-scale PDEs from computational fine-scale simulation data, the training data collection process can be prohibitively expensive. We propose to transformatively facilitate this training data collection process by linking machine learning (here, neural networks) with modern multiscale scientific computation (here, equation-free numerics). These equation-free techniques operate over sparse collections of small, appropriately coupled, space-time subdomains ("patches"), parsimoniously producing the required macro-scale training data. Our illustrative example involves the discovery of effective homogenized equations in one and two dimensions, for problems with fine-scale material property variations. The approach holds promise towards making the discovery of accurate, macro-scale effective materials PDE models possible by efficiently summarizing the physics embodied in "the best" fine-scale simulation models available.


title: Detecting non-local effects in the electronic structure of a simple covalent system with machine learning methods

abstract: Using methods borrowed from machine learning we detect in a fully algorithmic way long range effects on local physical properties in a simple covalent system of carbon atoms. The fact that these long range effects exist for many configurations implies that atomistic simulation methods, such as force fields or modern machine learning schemes, that are based on locality assumptions, are limited in accuracy. We show that the basic driving mechanism for the long range effects is charge transfer. If the charge transfer is known, locality can be recovered for certain quantities such as the band structure energy.


title: An Adversarial Attack Defending System for Securing In-Vehicle Networks

abstract: In a modern vehicle, there are over seventy Electronics Control Units (ECUs). For an in-vehicle network, ECUs communicate with each other by following a standard communication protocol, such as Controller Area Network (CAN). However, an attacker can easily access the in-vehicle network to compromise ECUs through a WLAN or Bluetooth. Though there are various deep learning (DL) methods suggested for securing in-vehicle networks, recent studies on adversarial examples have shown that attackers can easily fool DL models. In this research, we further explore adversarial examples in an in-vehicle network. We first discover and implement two adversarial attack models that are harmful to a Long Short Term Memory (LSTM)-based detection model used in the in-vehicle network. Then, we propose an Adversarial Attack Defending System (AADS) for securing an in-vehicle network. Specifically, we focus on brake-related ECUs in an in-vehicle network. Our experimental results demonstrate that adversaries can easily attack the LSTM-based detection model with a success rate of over 98%, and the proposed AADS achieves over 99% accuracy for detecting adversarial attacks.


title: Isotropic Cellular Automata: the DDLab iso-rule paradigm

abstract: To respect physics and nature, cellular automata (CA) models of self-organisation, emergence, computation and logical universality should be isotropic, having equivalent dynamics in all directions. We present a novel paradigm, the iso-rule, a concise expression for isotropic CA by the output table for each isotropic neighborhood group, allowing an efficient method of navigating and exploring iso-rule-space. We describe new functions and tools in DDLab to generate iso-groups and iso-rules, for multi-value as well as binary, in one, two and three dimensions. These methods include filing, filtering, mutating, analysing dynamics by input-frequency and entropy, identifying the critical iso-groups for glider-gun/eater dynamics, and automatically classifying iso-rule-space. We illustrate these ideas and methods for two dimensional CA on square and hexagonal lattices.


title: Synergies and Prospects for Early Resolution of the Neutrino Mass Ordering

abstract: The measurement of neutrino Mass Ordering (MO) is a fundamental element for the understanding of leptonic flavour sector of the Standard Model of Particle Physics. Its determination relies on the precise measurement of $\Delta m^2_{31}$ and $\Delta m^2_{32}$ using either neutrino vacuum oscillations, such as the ones studied by medium baseline reactor experiments, or matter effect modified oscillations such as those manifesting in long-baseline neutrino beams (LB$\nu$B) or atmospheric neutrino experiments. Despite existing MO indication today, a fully resolved MO measurement ($\geq$5$\sigma$) is most likely to await for the next generation of neutrino experiments: JUNO, whose stand-alone sensitivity is $\sim$3$\sigma$, or LB$\nu$B experiments (DUNE and Hyper-Kamiokande). Upcoming atmospheric neutrino experiments are also expected to provide precious information. In this work, we study the possible context for the earliest full MO resolution. A firm resolution is possible even before 2028, exploiting mainly vacuum oscillation, upon the combination of JUNO and the current generation of LB$\nu$B experiments (NOvA and T2K). This opportunity is possible thanks to a powerful synergy boosting the overall sensitivity where the sub-percent precision of $\Delta m^2_{32}$ by LB$\nu$B experiments is found to be the leading order term for the MO earliest discovery. We also found that the comparison between matter and vacuum driven oscillation results enables unique discovery potential for physics beyond the Standard Model.


title: Accelerating Federated Learning in Heterogeneous Data and Computational Environments

abstract: There are situations where data relevant to a machine learning problem are distributed among multiple locations that cannot share the data due to regulatory, competitiveness, or privacy reasons. For example, data present in users' cellphones, manufacturing data of companies in a given industrial sector, or medical records located at different hospitals. Moreover, participating sites often have different data distributions and computational capabilities. Federated Learning provides an approach to learn a joint model over all the available data in these environments. In this paper, we introduce a novel distributed validation weighting scheme (DVW), which evaluates the performance of a learner in the federation against a distributed validation set. Each learner reserves a small portion (e.g., 5%) of its local training examples as a validation dataset and allows other learners models to be evaluated against it. We empirically show that DVW results in better performance compared to established methods, such as FedAvg, both under synchronous and asynchronous communication protocols in data and computationally heterogeneous environments.


title: Explainable Spatial Clustering: Leveraging Spatial Data in Radiation Oncology

abstract: Advances in data collection in radiation therapy have led to an abundance of opportunities for applying data mining and machine learning techniques to promote new data-driven insights. In light of these advances, supporting collaboration between machine learning experts and clinicians is important for facilitating better development and adoption of these models. Although many medical use-cases rely on spatial data, where understanding and visualizing the underlying structure of the data is important, little is known about the interpretability of spatial clustering results by clinical audiences. In this work, we reflect on the design of visualizations for explaining novel approaches to clustering complex anatomical data from head and neck cancer patients. These visualizations were developed, through participatory design, for clinical audiences during a multi-year collaboration with radiation oncologists and statisticians. We distill this collaboration into a set of lessons learned for creating visual and explainable spatial clustering for clinical users.


title: Cosmology Intertwined I: Perspectives for the Next Decade

abstract: The standard $\Lambda$ Cold Dark Matter cosmological model provides an amazing description of a wide range of astrophysical and astronomical data. However, there are a few big open questions, that make the standard model look like a first-order approximation to a more realistic scenario that still needs to be fully understood. In this Letter of Interest we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances present between the different cosmological probes, as the Hubble constant $H_0$ value, the $\sigma_8 - S_8$ tension, and the anomalies present in the Planck results. Finally, we will give an overview of upgraded experiments and next-generation space-missions and facilities on Earth, that will be of crucial importance to address all these questions.


title: Cosmology Intertwined II: The Hubble Constant Tension

abstract: The current cosmological probes have provided a fantastic confirmation of the standard $\Lambda$ Cold Dark Matter cosmological model, that has been constrained with unprecedented accuracy. However, with the increase of the experimental sensitivity a few statistically significant tensions between different independent cosmological datasets emerged. While these tensions can be in portion the result of systematic errors, the persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the need for new physics. In this Letter of Interest we will focus on the $4.4\sigma$ tension between the Planck estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we will list a few interesting new physics models that could solve this tension and discuss how the next decade experiments will be crucial.


title: Cosmology Intertwined III: $f \sigma_8$ and $S_8$

abstract: The standard $\Lambda$ Cold Dark Matter cosmological model provides a wonderful fit to current cosmological data, but a few tensions and anomalies became statistically significant with the latest data analyses. While these anomalies could be due to the presence of systematic errors in the experiments, they could also indicate the need for new physics beyond the standard model. In this Letter of Interest we focus on the tension of the Planck data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models for solving this tension, and we discuss the importance of trying to fit with a single model a full array of data and not just one parameter at a time.


title: Cosmology Intertwined IV: The Age of the Universe and its Curvature

abstract: A precise measurement of the curvature of the Universe is of primeval importance for cosmology since it could not only confirm the paradigm of primordial inflation but also help in discriminating between different early Universe scenarios. The recent observations, while broadly consistent with a spatially flat standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) model, are showing tensions that still allow (and, in some cases, even suggest) a few percent deviations from a flat universe. In particular, the Planck Cosmic Microwave Background power spectra, assuming the nominal likelihood, prefer a closed universe at more than 99\% confidence level. While new physics could be in action, this anomaly may be the result of an unresolved systematic error or just a statistical fluctuation. However, since a positive curvature allows a larger age of the Universe, an accurate determination of the age of the oldest objects provides a smoking gun in confirming or falsifying the current flat $\Lambda$CDM model.


title: Measuring the spectral index of turbulent gas with deep learning from projected density maps

abstract: Turbulence plays a key role in star formation in molecular clouds, affecting star cluster primordial properties. As modelling present-day objects hinges on our understanding of their initial conditions, better constraints on turbulence can result in windfalls in Galactic archaeology, star cluster dynamics and star formation. Observationally, constraining the spectral index of turbulent gas usually involves computing spectra from velocity maps. Here we suggest that information on the spectral index might be directly inferred from column density maps (possibly obtained by dust emission/absorption) through deep learning. We generate mock density maps from a large set of adaptive mesh refinement turbulent gas simulations using the hydro-simulation code RAMSES. We train a convolutional neural network (CNN) on the resulting images to predict the turbulence index, optimize hyper-parameters in validation and test on a holdout set. Our adopted CNN model achieves a mean squared error of 0.024 in its predictions on our holdout set, over underlying spectral indexes ranging from 3 to 4.5. We also perform robustness tests by applying our model to altered holdout set images, and to images obtained by running simulations at different resolutions. This preliminary result on simulated density maps encourages further developments on real data, where observational biases and other issues need to be taken into account.


title: Can a protophobic vector boson explain the ATOMKI anomaly?

abstract: In 2016, the ATOMKI collaboration announced [PRL {\bf 116}, 042501 (2016)] observing an unexpected enhancement of the $e^+-e^-$ pair production signal in one of the $^8$Be nuclear transitions induced by an incident proton beam on a $^7$Li target. Many beyond-standard-model physics explanations have subsequently been proposed. One popular theory is that the anomaly is caused by the creation of a protophobic vector boson ($X$) with a mass around 17 MeV [e.g., PRL\ {\bf 117}, 071803 (2016)] in the nuclear transition. We study this hypothesis by deriving an isospin relation between photon and $X$ couplings to nucleons. This allows us to find simple relations between protophobic $X$-production cross sections and those for measured photon production. The net result is that $X$ production is dominated by direct transitions induced by $E1^X$ and $L1^X$ (transverse and longitudinal electric dipoles) and $C1^X$ (charge dipole) without going through any nuclear resonance (i.e. Bremsstrahlung radiation) with a smooth energy dependence that occurs for all proton beam energies above threshold. This contradicts the experimental observations and invalidates the protophobic vector boson explanation.


title: Robust Character Labeling in Movie Videos: Data Resources and Self-supervised Feature Adaptation

abstract: Robust face clustering is a vital step in enabling computational understanding of visual character portrayal in media. Face clustering for long-form content is challenging because of variations in appearance and lack of supporting large-scale labeled data. Our work in this paper focuses on two key aspects of this problem: the lack of domain-specific training or benchmark datasets, and adapting face embeddings learned on web images to long-form content, specifically movies. First, we present a dataset of over 169,000 face tracks curated from 240 Hollywood movies with weak labels on whether a pair of face tracks belong to the same or a different character. We propose an offline algorithm based on nearest-neighbor search in the embedding space to mine hard-examples from these tracks. We then investigate triplet-loss and multiview correlation-based methods for adapting face embeddings to hard-examples. Our experimental results highlight the usefulness of weakly labeled data for domain-specific feature adaptation. Overall, we find that multiview correlation-based adaptation yields more discriminative and robust face embeddings. Its performance on downstream face verification and clustering tasks is comparable to that of the state-of-the-art results in this domain. We also present the SAIL-Movie Character Benchmark corpus developed to augment existing benchmarks. It consists of racially diverse actors and provides face-quality labels for subsequent error analysis. We hope that the large-scale datasets developed in this work can further advance automatic character labeling in videos. All resources are available freely at https://sail.usc.edu/~ccmi/multiface.


title: Extractive Summarizer for Scholarly Articles

abstract: We introduce an extractive method that will summarize long scientific papers. Our model uses presentation slides provided by the authors of the papers as the gold summary standard to label the sentences. The sentences are ranked based on their novelty and their importance as estimated by deep neural networks. Our window-based extractive labeling of sentences results in the improvement of at least 4 ROUGE1-Recall points.


title: On Structured-Closed-Loop versus Structured-Controller Design: the Case of Relative Measurement Feedback

abstract: We consider the optimal distributed controller design problem subject to two structural requirements: locality, i.e. available measurements and sub-controllers' interactions are governed by a graph structure, and relative feedback, i.e. only differences of measurements are available to the controller. We formalize controller locality in terms of the controller's transfer function, state-space realization, or resulting closed-loop mapping. We demonstrate that the relative feedback requirement can be written as a convex constraint on the controller and (in special cases) on the resulting closed-loop, and we characterize the allowable structures of relative feedback controllers. We prove that sparse closed-loop design is a convex relaxation of structured controller state-space design, even in the continuous time IIR setting. This formalizes and extends results of the recently developed System Level Synthesis framework. We take a first step toward quantifying the performance gap associated with this convex relaxation by constructing a class of examples (based on relative feedback requirements) for which the difference in performance, measured by an H2 norm, is infinite. The results presented are used to contrast several issues of structural constraints in distributed control design that remain as open problems.


title: Flip Paths Between Lattice Triangulations

abstract: We present a $O(n^{\frac{3}{2}})$-time algorithm for the \emph{shortest (diagonal) flip path problem} for \emph{lattice} triangulations with $n$ points, improving over previous $O(n^2)$-time algorithms. For a large, natural class of inputs, our bound is tight in the sense that our algorithm runs in time linear in the number of flips in the output flip path. Our results rely on an independently interesting structural elucidation of shortest flip paths as the linear orderings of a unique partially ordered set, called a \emph{minimum flip plan}, constructed by a novel use of Farey sequences from elementary number theory. Flip paths between general (not necessarily lattice) triangulations have been studied in the combinatorial setting for nearly a century. In the Euclidean geometric setting, finding a shortest flip path between two triangulations is NP-complete. However, for lattice triangulations, which are studied as spin systems, there are known $O\left(n^2\right)$-time algorithms to find shortest flip paths. These algorithms, as well as ours, apply to \emph{constrained} flip paths that ensure a set of \emph{constraint} edges are present in every triangulation along the path. Implications for determining simultaneously flippable edges, i.e. finding optimal simultaneous flip paths between lattice triangulations, and for counting lattice triangulations are discussed.


title: Generating (Factual?) Narrative Summaries of RCTs: Experiments with Neural Multi-Document Summarization

abstract: We consider the problem of automatically generating a narrative biomedical evidence summary from multiple trial reports. We evaluate modern neural models for abstractive summarization of relevant article abstracts from systematic reviews previously conducted by members of the Cochrane collaboration, using the authors conclusions section of the review abstract as our target. We enlist medical professionals to evaluate generated summaries, and we find that modern summarization systems yield consistently fluent and relevant synopses, but that they are not always factual. We propose new approaches that capitalize on domain-specific models to inform summarization, e.g., by explicitly demarcating snippets of inputs that convey key findings, and emphasizing the reports of large and high-quality trials. We find that these strategies modestly improve the factual accuracy of generated summaries. Finally, we propose a new method for automatically evaluating the factuality of generated narrative evidence syntheses using models that infer the directionality of reported findings.


title: Measuring the Capabilities of Quantum Computers

abstract: A quantum computer has now solved a specialized problem believed to be intractable for supercomputers, suggesting that quantum processors may soon outperform supercomputers on scientifically important problems. But flaws in each quantum processor limit its capability by causing errors in quantum programs, and it is currently difficult to predict what programs a particular processor can successfully run. We introduce techniques that can efficiently test the capabilities of any programmable quantum computer, and we apply them to twelve processors. Our experiments show that current hardware suffers complex errors that cause structured programs to fail up to an order of magnitude earlier - as measured by program size - than disordered ones. As a result, standard error metrics inferred from random disordered program behavior do not accurately predict performance of useful programs. Our methods provide efficient, reliable, and scalable benchmarks that can be targeted to predict quantum computer performance on real-world problems.


title: Concept Extraction Using Pointer-Generator Networks

abstract: Concept extraction is crucial for a number of downstream applications. However, surprisingly enough, straightforward single token/nominal chunk-concept alignment or dictionary lookup techniques such as DBpedia Spotlight still prevail. We propose a generic open-domain OOV-oriented extractive model that is based on distant supervision of a pointer-generator network leveraging bidirectional LSTMs and a copy mechanism. The model has been trained on a large annotated corpus compiled specifically for this task from 250K Wikipedia pages, and tested on regular pages, where the pointers to other pages are considered as ground truth concepts. The outcome of the experiments shows that our model significantly outperforms standard techniques and, when used on top of DBpedia Spotlight, further improves its performance. The experiments furthermore show that the model can be readily ported to other datasets on which it equally achieves a state-of-the-art performance.


title: Sequences of well-distributed vertices on graphs and spectral bounds on optimal transport

abstract: Given a graph $G=(V,E)$, suppose we are interested in selecting a sequence of vertices $(x_j)_{j=1}^n$ such that $\left\{x_1, \dots, x_k\right\}$ is `well-distributed' uniformly in $k$. We describe a greedy algorithm motivated by potential theory and corresponding developments in the continuous setting. The algorithm performs nicely on graphs and may be of use for sampling problems. We can interpret the algorithm as trying to greedily minimize a negative Sobolev norm; we explain why this is related to Wasserstein distance by establishing a purely spectral bound on the Wasserstein distance on graphs that mirrors R. Peyre's estimate in the continuous setting. We illustrate this with many examples and discuss several open problems.


title: Transductive Information Maximization For Few-Shot Learning

abstract: We introduce Transductive Infomation Maximization (TIM) for few-shot learning. Our method maximizes the mutual information between the query features and their label predictions for a given few-shot task, in conjunction with a supervision loss based on the support set. Furthermore, we propose a new alternating-direction solver for our mutual-information loss, which substantially speeds up transductive-inference convergence over gradient-based optimization, while yielding similar accuracy. TIM inference is modular: it can be used on top of any base-training feature extractor. Following standard transductive few-shot settings, our comprehensive experiments demonstrate that TIM outperforms state-of-the-art methods significantly across various datasets and networks, while used on top of a fixed feature extractor trained with simple cross-entropy on the base classes, without resorting to complex meta-learning schemes. It consistently brings between 2% and 5% improvement in accuracy over the best performing method, not only on all the well-established few-shot benchmarks but also on more challenging scenarios,with domain shifts and larger numbers of classes.


title: Rethinking Non-idealities in Memristive Crossbars for Adversarial Robustness in Neural Networks

abstract: Deep Neural Networks (DNNs) have been shown to be prone to adversarial attacks. Memristive crossbars, being able to perform Matrix-Vector-Multiplications (MVMs) efficiently, are used to realize DNNs on hardware. However, crossbar non-idealities have always been devalued since they cause errors in performing MVMs, leading to computational accuracy losses in DNNs. Several software-based defenses have been proposed to make DNNs adversarially robust. However, no previous work has demonstrated the advantage conferred by the crossbar non-idealities in unleashing adversarial robustness. We show that the intrinsic hardware non-idealities yield adversarial robustness to the mapped DNNs without any additional optimization. We evaluate the adversarial resilience of state-of-the-art DNNs (VGG8 & VGG16 networks) using benchmark datasets (CIFAR-10, CIFAR-100 & Tiny Imagenet) across various crossbar sizes. We find that crossbar non-idealities unleash significantly greater adversarial robustness (>10-20%) in crossbar-mapped DNNs than baseline software DNNs. We further assess the performance of our approach with other state-of-the-art efficiency-driven adversarial defenses and find that our approach performs significantly well in terms of reducing adversarial loss.


title: Prospects for Characterizing the Haziest Sub-Neptune Exoplanets with High Resolution Spectroscopy

abstract: Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R $\sim$ 25,000 to 100,000) enables one to probe the regions in these atmospheres above the clouds where the cores of the strongest spectral lines are formed. We present models of transmission spectra for a suite of GJ1214b-like planets with thick photochemical hazes covering 1 - 5 $\mu$m at a range of resolutions relevant to current and future ground-based spectrographs. Furthermore, we compare the utility of the cross-correlation function that is typically used with a more formal likelihood-based approach, finding that only the likelihood based method is sensitive to the presence of haze opacity. We calculate the signal-to-noise of these spectra, including telluric contamination, required to robustly detect a host of molecules such as CO, CO$_{2}$, H$_{2}$O, and CH$_{4}$, and photochemical products like HCN, as a function of wavelength range and spectral resolution. Spectra in M band require the lowest S/N$_{res}$ to detect multiple molecules simultaneously. CH$_{4}$ is only observable for the coolest models ($T_{\rm{eff}} =$ 412 K) and only in the L band. We quantitatively assess how these requirements compare to what is achievable with current and future instruments, demonstrating that characterization of small cool worlds with ground-based high resolution spectroscopy is well within reach.


title: Likelihood Landscapes: A Unifying Principle Behind Many Adversarial Defenses

abstract: Convolutional Neural Networks have been shown to be vulnerable to adversarial examples, which are known to locate in subspaces close to where normal data lies but are not naturally occurring and of low probability. In this work, we investigate the potential effect defense techniques have on the geometry of the likelihood landscape - likelihood of the input images under the trained model. We first propose a way to visualize the likelihood landscape leveraging an energy-based model interpretation of discriminative classifiers. Then we introduce a measure to quantify the flatness of the likelihood landscape. We observe that a subset of adversarial defense techniques results in a similar effect of flattening the likelihood landscape. We further explore directly regularizing towards a flat landscape for adversarial robustness.


title: Estimating conditional probabilities of historical migrations in the transatlantic slave trade using kriging and Markov decision process models

abstract: Intra-African conflicts during the collapse of the kingdom of Oyo from 1817 to 1836 resulted in the enslavement of an estimated 121,000 people who were then transported to coastal ports via complex trade networks and loaded onto slave ships destined for the Americas. Historians have a good record of where these people went across the Atlantic, but little is known about where individuals were from or enslaved \textit{within} Africa. In this work, we develop a novel two-step statistical approach to describe the enslavement of people given documented violent conflict, the transport of enslaved peoples from their location of capture to their port of departure, and---given an enslaved individual's location of departure---that person's probability of origin. We combine spatial prediction of conflict density via Kriging with a Markov decision process characterising intra-African transportation. The results of this model can be visualised using an interactive web application, plotting estimated conditional probabilities of historical migrations during the African diaspora. These results help trace the uncertain origins of people enslaved in this region of Africa during this time period: using the two-step statistical methodology developed here provides a probabilistic answer to this question.


title: Hall Lie algebras of toric monoid schemes

abstract: We associate to a projective $n$-dimensional toric variety $X_{\Delta}$ a pair of co-commutative (but generally non-commutative) Hopf algebras $H^{\alpha}_X, H^{T}_X$. These arise as Hall algebras of certain categories $\Coh^{\alpha}(X), \Coh^T(X)$ of coherent sheaves on $X_{\Delta}$ viewed as a monoid scheme - i.e. a scheme obtained by gluing together spectra of commutative monoids rather than rings. When $X_{\Delta}$ is smooth, the category $\Coh^T(X)$ has an explicit combinatorial description as sheaves whose restriction to each $\mathbb{A}^n$ corresponding to a maximal cone $\sigma \in \Delta$ is determined by an $n$-dimensional generalized skew shape. The (non-additive) categories $\Coh^{\alpha}(X), \Coh^T(X)$ are treated via the formalism of proto-exact/proto-abelian categories developed by Dyckerhoff-Kapranov. The Hall algebras $H^{\alpha}_X, H^{T}_X$ are graded and connected, and so enveloping algebras $H^{\alpha}_X \simeq U(\n^{\alpha}_X)$, $H^{T}_X \simeq U(\n^{T}_X)$, where the Lie algebras $\n^{\alpha}_X, \n^{T}_X$ are spanned by the indecomposable coherent sheaves in their respective categories. We explicitly work out several examples, and in some cases are able to relate $\n^T_X$ to known Lie algebras. In particular, when $X = \mathbb{P}^1$, $\n^T_X$ is isomorphic to a non-standard Borel in $\mathfrak{gl}_2 [t,t^{-1}]$. When $X$ is the second infinitesimal neighborhood of the origin inside $\mathbb{A}^2$, $\n^T_X$ is isomorphic to a subalgebra of $\mathfrak{gl}_2[t]$. We also consider the case $X=\mathbb{P}^2$, where we give a basis for $\n^T_X$ by describing all indecomposable sheaves in $\Coh^T(X)$.


title: Integrated Cutting and Packing Heterogeneous Precast Beams Multiperiod Production Planning Problem

abstract: We introduce a novel variant of cutting production planning problems named Integrated Cutting and Packing Heterogeneous Precast Beams Multiperiod Production Planning (ICP-HPBMPP). We propose an integer linear programming model for the ICP-HPBMPP, as well as a lower bound for its optimal objective function value, which is empirically shown to be closer to the optimal solution value than the bound obtained from the linear relaxation of the model. We also propose a genetic algorithm approach for the ICP-HPBMPP as an alternative solution method. We discuss computational experiments and propose a parameterization for the genetic algorithm using D-optimal experimental design. We observe good performance of the exact approach when solving small-sized instances, although there are difficulties in finding optimal solutions for medium and large-sized problems, or even in finding feasible solutions for large instances. On the other hand, the genetic algorithm could find good-quality solutions for large-sized instances within short computing times.


title: On quiver representations over $\mathbb{F}_1$

abstract: We study the category $\textrm{Rep}(Q,\mathbb{F}_1)$ of representations of a quiver $Q$ over "the field with one element", denoted by $\mathbb{F}_1$, and the Hall algebra of $\textrm{Rep}(Q,\mathbb{F}_1)$. Representations of $Q$ over $\mathbb{F}_1$ often reflect combinatorics of those over $\mathbb{F}_q$, but show some subtleties - for example, we prove that a connected quiver $Q$ is of finite representation type over $\mathbb{F}_1$ if and only if $Q$ is a tree. Then, to each representation $\mathbb{V}$ of $Q$ over $\mathbb{F}_1$ we associate a coefficient quiver $\Gamma_\mathbb{V}$ possessing the same information as $\mathbb{V}$. This allows us to translate representations over $\mathbb{F}_1$ purely in terms of combinatorics of associated coefficient quivers. We also explore the growth of indecomposable representations of $Q$ over $\mathbb{F}_1$ - there are also similarities to representations over a field, but with some subtle differences. Finally, we link the Hall algebra of the category of nilpotent representations of an $n$-loop quiver over $\mathbb{F}_1$ with the Hopf algebra of skew shapes introduced by Szczesny.


title: Long-term word frequency dynamics derived from Twitter are corrupted: A bespoke approach to detecting and removing pathologies in ensembles of time series

abstract: Maintaining the integrity of long-term data collection is an essential scientific practice. As a field evolves, so too will that field's measurement instruments and data storage systems, as they are invented, improved upon, and made obsolete. For data streams generated by opaque sociotechnical systems which may have episodic and unknown internal rule changes, detecting and accounting for shifts in historical datasets requires vigilance and creative analysis. Here, we show that around 10\% of day-scale word usage frequency time series for Twitter collected in real time for a set of roughly 10,000 frequently used words for over 10 years come from tweets with, in effect, corrupted language labels. We describe how we uncovered problematic signals while comparing word usage over varying time frames. We locate time points where Twitter switched on or off different kinds of language identification algorithms, and where data formats may have changed. We then show how we create a statistic for identifying and removing words with pathological time series. While our resulting process for removing `bad' time series from ensembles of time series is particular, the approach leading to its construction may be generalizeable.


title: Transverse linear subspaces to hypersurfaces over finite fields

abstract: Ballico proved that a smooth projective variety $X$ of degree $d$ over a finite field of $q$ elements admits a smooth hyperplane section if $q\geq d(d-1)^{\dim X}$. In this paper, we refine this criterion for higher codimensional linear sections on smooth hypersurfaces and for hyperplane sections on Frobenius classical hypersurfaces. We also prove a similar result for the existence of reduced hyperplane sections on reduced hypersurfaces.


title: Extended multi-scalar field theories in (1+1) dimensions

abstract: We present the explicit construction of some multi-scalar field theories in (1+1) dimensions supporting BPS (Bogomol'nyi--Prasad--Sommerfield) kink solutions. The construction is based on the ideas of the so-called extension method. In particular, several new interesting two-scalar and three-scalar field theories are explicitly constructed from non-trivial couplings between well-known one-scalar field theories. The BPS solutions of the original one-field systems will be also BPS solutions of the multi-scalar system by construction, and therefore we will analyse their linear stability properties for the constructed models.


title: Identifying Coordinated Accounts on Social Media through Hidden Influence and Group Behaviours

abstract: Disinformation campaigns on social media, involving coordinated activities from malicious accounts towards manipulating public opinion, have become increasingly prevalent. Existing approaches to detect coordinated accounts either make very strict assumptions about coordinated behaviours, or require part of the malicious accounts in the coordinated group to be revealed in order to detect the rest. To address these drawbacks, we propose a generative model, AMDN-HAGE (Attentive Mixture Density Network with Hidden Account Group Estimation) which jointly models account activities and hidden group behaviours based on Temporal Point Processes (TPP) and Gaussian Mixture Model (GMM), to capture inherent characteristics of coordination which is, accounts that coordinate must strongly influence each other's activities, and collectively appear anomalous from normal accounts. To address the challenges of optimizing the proposed model, we provide a bilevel optimization algorithm with theoretical guarantee on convergence. We verified the effectiveness of the proposed method and training algorithm on real-world social network data collected from Twitter related to coordinated campaigns from Russia's Internet Research Agency targeting the 2016 U.S. Presidential Elections, and to identify coordinated campaigns related to the COVID-19 pandemic. Leveraging the learned model, we find that the average influence between coordinated account pairs is the highest.On COVID-19, we found coordinated group spreading anti-vaccination, anti-masks conspiracies that suggest the pandemic is a hoax and political scam.


title: Ensemble-filtered vortex modeling of strongly disturbed aerodynamic flows

abstract: The task of dynamic flow estimation is to construct an approximation of an evolving flow---and particularly, its response to disturbances---using measurements from available sensors. Building from previous work by Darakananda et al.~(Phys Rev Fluids 2018), we further develop an ensemble Kalman filter (EnKF) framework for aerodynamic flows based on an ensemble of randomly-perturbed inviscid vortex models of flow about an infinitely-thin plate. In the forecast step, vortex elements in each ensemble member are advected by the flow and new elements are released from each edge of the plate; the elements are aggregated to maintain an efficient representation. The vortex elements and leading edge constraint are corrected in the analysis step by assimilating the surface pressure differences across the plate measured from the truth system. We show that the overall framework can be physically interpreted as a series of adjustments to the position and shape of an elliptical region of uncertainty associated with each vortex element. In this work, we compare the previously-used stochastic EnKF with the ensemble transform Kalman filter (ETKF), which uses a deterministic analysis step. We examine the response of the flat plate at $20^\circ$ in two perturbed flows, with truth data obtained from high-fidelity simulation at Reynolds number 500. In the first case, we apply a sequence of large-amplitude pulses near the leading edge of the plate to mimic flow actuation. In the second, we place the plate in a vortex street wake behind a cylinder. In both cases, we show that the vortex-based framework accurately estimates the pressure distribution and normal force, with no {\em a priori} knowledge of the perturbations. We show that the ETKF is consistently more robust than the stochastic EnKF. Finally, we examine the mapping from measurements to state update in the analysis step through SVD of the Kalman gain.


title: Why Shouldn't All Charts Be Scatter Plots? Beyond Precision-Driven Visualizations

abstract: A central concept in information visualization research and practice is the notion of visual variable effectiveness, or the perceptual precision at which values are decoded given visual channels of encoding. Formative work from Cleveland & McGill has shown that position along a common axis is the most effective visual variable for comparing individual values. One natural conclusion is that any chart that is not a dot plot or scatterplot is deficient and should be avoided. In this paper we refute a caricature of this "scatterplots only" argument as a way to call for new perspectives on how information visualization is researched, taught, and evaluated.


title: Symmetry and Art

abstract: We use some fundamental ideas from complex analysis to create symmetric images and animations. Using a domain coloring algorithm, we generate mappings to the entire complex plane or the hyperbolic upper half-plane. The resulting designs can have rotational, translational, or mirror symmetry according to our chosen mapping functions. An appealing feature of these designs is how they reveal important properties of Euclidean and non-Euclidean geometries. We can also generate animations of our designs. Our goal is to create designs and animations having significant artistic content.


title: Theoretical study of the $D^0 \to K^- \pi^+ \eta$ reaction

abstract: We develop a model to study the $D^0 \to K^- \pi^+ \eta$ weak decay, starting with the color favored external emission and Cabibbo favored mode at the quark level. A less favored internal emission decay mode is also studied as a source of small corrections. Some pairs of quarks are allowed to hadronize producing two pseudoscalar mesons, which posteriorly are allowed to interact to finally provide the $K^- \pi^+ \eta$ state. The chiral unitary approach is used to take into account the final state interaction of pairs of mesons, which has as a consequence the production of the $\kappa$ ($K^*_0(700)$) and the $a_0(980)$ resonances, very well visible in the invariant mass distributions. We also introduce the $\bar{K}^{*0} \eta$ production in a phenomenological way and show that the $s$-wave pseudoscalar interaction together with this vector excitation mode are sufficient to provide a fair reproduction of the experimental data. The agreement with the data, in particular the relative weight of the $a_0(980)$ to the $\kappa$ excitation, provides extra support to the picture used, in which these two resonances are a consequence of the interaction of pseudoscalar mesons and not ordinary $q \bar{q}$ mesons.


title: Soft X-Ray Observations of Quiescent Solar Active Regions using Novel Dual-zone Aperture X-ray Solar Spectrometer (DAXSS)

abstract: The Dual-zone Aperture X-ray Solar Spectrometer (DAXSS) was flown on 2018 June18 on the NASA 36.336 sounding rocket flight and obtained the highest resolution to date for solar soft X-ray (SXR) spectra over a broad energy range. This observation was during a time with quiescent (non-flaring) small active regions on the solar disk and when the 10.7 cm radio flux (F10.7) was 75 solar flux units (1 sfu = 10-22 W/m^2/Hz). The DAXSS instrument consists of a LASP-developed dual-zone aperture and a commercial X-ray spectrometer from Amptek that measures solar full-disk irradiance from 0.5-20 keV with a resolving power of 20 near 1 keV. This paper discusses the novel design of the spectrometer and the instrument characterization techniques. Additionally,the solar measurements obtained from the 2018 sounding rocket flight are analyzed using CHIANTI spectral models to fit the temperatures, emission measures, and relative elemental abundances of the solar corona plasma. The abundance of iron was found to be 35 percent higher than expected in the quiescent sun's corona suggesting either that our spectral models require additional sophistication or that the underlying atomic database may require updates. Future long-term systematic observations of this spectral range are needed. DAXSS will fly on the INSPIRESat-1 CubeSat in late-2020, and its SXR spectral data could provide further insight into the sources of coronal heating through modeling the changes of relative elemental abundances during developments of active regions and solar flaring events.


title: Generalized Passarino-Veltman reduction for relativistic quantum field theories at finite temperature and finite density

abstract: The \emph{conventional} Passarino-Veltman reduction (CPVR) is a systematic procedure based on the Lorentz symmetry, which can efficiently reduce the one-loop Feynman integrals (OLFIs) in the relativistic quantum field theories (RQFTs) at zero temperature and zero density. However, at finite temperature and finite density, the Lorentz symmetry is explicitly broken due to a rest reference frame of the many-body system in which the temperature and density are measured, rendering the CPVR not applicable any more to reduce OLFIs therein. In this paper, we report a \emph{generalized} Passarino-Veltman reduction (GPVR), which can simplify the OLFIs in the RQFTs at finite temperature and finite density. The GPVR can analyze the OLFIs in a wide range of physical systems described by the RQFTs at finite temperature and finite density, such as quark-gluon plasma in nuclear physics.

categories: hep-ph cond-mat.mes-hall cond-mat.str-el nucl-th


title: Inapproximability of Diameter in super-linear time: Beyond the 5/3 ratio

abstract: We show, assuming the Strong Exponential Time Hypothesis, that for every $\varepsilon > 0$, approximating directed Diameter on $m$-arc graphs within ratio $7/4 - \varepsilon$ requires $m^{4/3 - o(1)}$ time. Our construction uses nonnegative edge weights but even holds for sparse digraphs, i.e., for which the number of vertices $n$ and the number of arcs $m$ satisfy $m = n \log^{O(1)} n$. This is the first result that conditionally rules out a near-linear time $5/3$-approximation for Diameter.


title: DQSEGDB: A time-interval database for storing gravitational wave observatory metadata

abstract: The Data Quality Segment Database (DQSEGDB) software is a database service, backend API, frontend graphical web interface, and client package used by the Laser Interferometer Gravitational-Wave Observatory (LIGO), Virgo, GEO600 and the Kamioka Gravitational wave detector for storing and accessing metadata describing the status of their detectors. The DQSEGDB has been used in the analysis of all published detections of gravitational waves in the advanced detector era. The DQSEGDB currently stores roughly 600 million metadata entries and responds to roughly 600,000 queries per day with an average response time of 0.223 ms.


title: Empirical Understanding of Deletion Privacy: Experiences, Expectations, and Measures

abstract: Social platforms are heavily used by individuals to share their thoughts and personal information. However, due to regret over time about posting inappropriate social content, embarrassment, or even life or relationship changes, some past posts might also pose serious privacy concerns for them. To cope with these privacy concerns, social platforms offer deletion mechanisms that allow users to remove their contents. Quite naturally, these deletion mechanisms are really useful for removing past posts as and when needed. However, these same mechanisms also leave the users potentially vulnerable to attacks by adversaries who specifically seek the users' damaging content and exploit the act of deletion as a strong signal for identifying such content. Unfortunately, today user experiences and contextual expectations regarding such attacks on deletion privacy and deletion privacy in general are not well understood. To that end, in this paper, we conduct a user survey-based exploration involving 191 participants to unpack their prior deletion experiences, their expectations of deletion privacy, and how effective they find the current deletion mechanisms. We find that more than 80% of the users have deleted at least a social media post, and users self-reported that, on average, around 35% of their deletions happened after a week of posting. While the participants identified the irrelevancy (due to time passing) as the main reason for content removal, most of them believed that deletions indicate that the deleted content includes some damaging information to the owner. Importantly, the participants are significantly more concerned about their deletions being noticed by large-scale data collectors (e.g., the government) than individuals from their social circle. Finally, the participants felt that popular deletion mechanisms are not very effective in protecting the privacy of those deletions.


title: Quantum walks with quantum chaotic coins: Of the Loschmidt echo, classical limit and thermalization

abstract: Coined discrete-time quantum walks are studied using simple deterministic dynamical systems as coins whose classical limit can range from being integrable to chaotic. It is shown that a Loschmidt echo like fidelity plays a central role and when the coin is chaotic this is approximately the characteristic function of a classical random walker. Thus the classical binomial distribution arises as a limit of the quantum walk and the walker exhibits diffusive growth before eventually becoming ballistic. The coin-walker entanglement growth is shown to be logarithmic in time as in the case of many-body localization and coupled kicked rotors, and saturates to a value that depends on the relative coin and walker space dimensions. In a coin dominated scenario, the chaos can thermalize the quantum walk to typical random states such that the entanglement saturates at the Haar averaged Page value, unlike in a walker dominated case when atypical states seem to be produced.


title: TradAO: A Visual Analytics System for Trading Algorithm Optimization

abstract: With the wide applications of algorithmic trading, it has become critical for traders to build a winning trading algorithm to beat the market. However, due to the lack of efficient tools, traders mainly rely on their memory to manually compare the algorithm instances of a trading algorithm and further select the best trading algorithm instance for the real trading deployment. We work closely with industry practitioners to discover and consolidate user requirements and develop an interactive visual analytics system for trading algorithm optimization. Structured expert interviews are conducted to evaluateTradAOand a representative case study is documented for illustrating the system effectiveness. To the best of our knowledge, previous financial data visual analyses have mainly aimed to assist investment managers in investment portfolio analysis but have neglected the need of traders in developing trading algorithms for portfolio execution.TradAOis the first visual analytics system that assists users in comprehensively exploring the performances of a trading algorithm with different parameter settings.


title: CYRA: the cryogenic infrared spectrograph for the Goode Solar Telescope in Big Bear

abstract: CYRA (CrYogenic solar spectrogRAph) is a facility instrument of the 1.6-meter Goode Solar Telescope (GST) at the Big Bear Solar Observatory (BBSO). CYRA focuses on the study of the near-infrared solar spectrum between 1 and 5 microns, a under explored region which is not only a fertile ground for photospheric magnetic diagnostics, but also allows a unique window into the chromosphere lying atop the photosphere. CYRA is the first ever fully cryogenic spectrograph in any solar observatory with its two predecessors, on the McMath-Pierce and Mees Telescopes, being based on warm optics except for the detectors and order sorting filters. CYRA is used to probe magnetic fields in various solar features and the quiet photosphere. CYRA measurements will allow new and better 3D extrapolations of the solar magnetic field and will provide more accurate boundary conditions for solar activity models. Superior spectral resolution of 150,000 and better allows enhanced observations of the chromosphere in the carbon monoxide (CO) spectral bands and will yield a better understanding of energy transport in the solar atmosphere. CYRA is divided into two optical sub-systems: The Fore-Optics Module and the Spectrograph. The Spectrograph is the heart of the instrument and contains the IR detector, grating, slits, filters, and imaging optics all in a cryogenically cooled Dewar (cryostat). The detector a 2048 by 2048 pixel HAWAII 2 array produced by Teledyne Scientific & Imaging, LLC. The interior of the cryostat and the readout electronics are maintained at 90 Kelvin by helium refrigerant based cryo-coolers, while the IR array is cooled to 30 Kelvin. The Fore-Optics Module de-rotates and stabilizes the solar image, provides scanning capabilities, and transfers the GST image to the Spectrograph. CYRA has been installed and is undergoing its commissioning phase.


title: High-Performance Parallel Graph Coloring with Strong Guarantees on Work, Depth, and Quality

abstract: We develop the first parallel graph coloring heuristics with strong theoretical guarantees on work and depth and coloring quality. The key idea is to design a relaxation of the vertex degeneracy order, a well-known graph theory concept, and to color vertices in the order dictated by this relaxation. This introduces a tunable amount of parallelism into the degeneracy ordering that is otherwise hard to parallelize. This simple idea enables significant benefits in several key aspects of graph coloring. For example, one of our algorithms ensures polylogarithmic depth and a bound on the number of used colors that is superior to all other parallelizable schemes, while maintaining work-efficiency. In addition to provable guarantees, the developed algorithms have competitive run-times for several real-world graphs, while almost always providing superior coloring quality. Our degeneracy ordering relaxation is of separate interest for algorithms outside the context of coloring.


title: Simplified calculations of plasma oscillations with non-extensive statistics

abstract: We use the exponential parametrization of the nonextensive distribution to calculate the dielectric constant in an electron gas obeying the nonextensive statistics. As we show, the exponential parametrization allows us to make such calculations in a straightforward way, bypassing the use of intricate formulas obtained from integral tables and/or numerical methods. For illustrative purposes, we apply first the method to the calculation of the permittivity and the corresponding dispersion relation in the ultrarelativistic limit of the electron gas, and verify that it reproduces in a simple way the results that had been obtained previously by other authors using the standard parametrization of the nonextensive distribution. In the same spirit we revisit the calculation of the same quantities for a non-relativistic gas, in the high frequency limit, which has been previously carried out, first by Lima, Silva and Santos, and subsequently revised by Chen and Li. Our own results agree with those obtained by Chen and Li. For completeness, we also apply the method the low frequency limit in the non-relativistic case, which has been previously considered by Dai, Chen and Li in the context of the stream plasma instability. We discuss some features of the results obtained in each case and their interpretation of terms of generalized nonextensive quantities, such as the Debye length $\lambda_{D}^{(q)}$, the plasma frequency $\omega_{p}^{(q)}$ and the ultra-relativistic frequency $\Omega^{(q)}_{e,rel}$. In the limit $q \rightarrow 1$ such quantities reduce to their classical value and the classical result of the dispersion relations are reproduced.


title: Enriched infinity categories I: enriched presheaves

abstract: This is the first of a series of papers on enriched infinity categories, seeking to reduce enriched higher category theory to the higher algebra of presentable infinity categories, which is better understood and can be approached via universal properties. In this paper, we introduce enriched presheaves on an enriched infinity category. We prove analogues of most familiar properties of presheaves. For example, we compute limits and colimits of presheaves, prove that all presheaves are colimits of representable presheaves, and prove a version of the Yoneda lemma.


title: Decay widthes of $^3 P_J$ charmonium to $DD,DD^*,D^*D^*$ and corresponding mass shifts of $^3 P_J$ charmonium

abstract: In this work, we calculate the amplitudes of the processes $c\bar c({^3P_J}) \rightarrow DD,DD^*, D^*D^* \rightarrow c\bar c({^3P_J})$ in the leading order of the nonrelativistic expansion. The imaginary parts of the amplitudes are corresponding to the branch decay widthes of the charmonium $c\bar c({^3P_J}) \rightarrow DD,DD^*, D^*D^*$ and the real parts are corresponding to the mass shifts of the charmonium $c\bar c({^3P_J})$ due to these decay channels. After absorbing the polynomial contributions which are pure real and include the UV divergences, the ratios between the branch decay widthes and the corresponding mass shifts are only dependent on the center-of-mass energy. We find the decay widthes and the mass shifts of the $^3P_2$ states are exact zero in the leading order. The ratios between the branch decay widthes and the mass shifts for the $^3P_0, {^3P_1}$ states are larger than 5 when the center-of-mass energy is above the $DD,DD^*, D^*D^*$ threshold. The dependence of the mass shifts on the center-of-mass energy is nontrivial especially when the center-of-mass energy is below the threshold. The analytic results can be extended to the $b$ quark sector directly.


title: Gaussian Primes in Narrow Sectors

abstract: We generalize a Theorem of Ricci and count Gaussian primes $\mathfrak{p}$ with short interval restrictions on both the norm and the argument of $\mathfrak{p}$.


title: 8 Steps to 3.7 TFLOP/s on NVIDIA V100 GPU: Roofline Analysis and Other Tricks

abstract: Performance optimization can be a daunting task especially as the hardware architecture becomes more and more complex. This paper takes a kernel from the Materials Science code BerkeleyGW, and demonstrates a few performance analysis and optimization techniques. Despite challenges such as high register usage, low occupancy, complex data access patterns, and the existence of several long-latency instructions, we have achieved 3.7 TFLOP/s of double-precision performance on an NVIDIA V100 GPU, with 8 optimization steps. This is 55% of the theoretical peak, 6.7 TFLOP/s, at nominal frequency 1312 MHz, and 70% of the more customized peak based on our 58% FMA ratio, 5.3 TFLOP/s. An array of techniques used to analyze this OpenACC kernel and optimize its performance are shown, including the use of hierarchical Roofline performance model and the performance tool Nsight Compute. This kernel exhibits computational characteristics that are commonly seen in many high-performance computing (HPC) applications, and are expected to be very helpful to a general audience of HPC developers and computational scientists, as they pursue more performance on NVIDIA GPUs.