title: A scalable spectral Stokes solver for simulation of time-periodic flows in complex geometries

abstract: Simulation of unsteady creeping flows in complex geometries has traditionally required the use of a time-stepping procedure, which is typically costly and unscalable. To reduce the cost and allow for computations at much larger scales, we propose an alternative approach that is formulated based on the unsteady Stokes equation expressed in the time-spectral domain. This transformation results in a boundary value problem with an imaginary source term proportional to the computed mode that is discretized and solved in a complex-valued finite element solver using Bubnov-Galerkin formulation. This transformed spatio-spectral formulation presents several advantages over the traditional spatio-temporal techniques. Firstly, for cases with boundary conditions varying smoothly in time, it provides a significant saving in computational cost as it can resolve time-variation of the solution using a few modes rather than thousands of time steps. Secondly, in contrast to the traditional time integration scheme with a finite order of accuracy, this method exhibits a super convergence behavior versus the number of computed modes. Thirdly, in contrast to the stabilized finite element methods for fluid, no stabilization term is employed in our formulation, producing a solution that is consistent and more accurate. Fourthly, the proposed approach is embarrassingly parallelizable owing to the independence of the solution modes, thus enabling scalable calculations at a much larger number of processors. The comparison of the proposed technique against a standard stabilized finite element solver is performed using two- and three-dimensional canonical and complex geometries. The results show that the proposed method can produce more accurate results at 1% to 11% of the cost of the standard technique for the studied cases.


title: Decomposition algorithms for solving NP-hard problems on a quantum annealer

abstract: NP-hard problems such as the maximum clique or minimum vertex cover problems, two of Karp's 21 NP-hard problems, have several applications in computational chemistry, biochemistry and computer network security. Adiabatic quantum annealers can search for the optimum value of such NP-hard optimization problems, given the problem can be embedded on their hardware. However, this is often not possible due to certain limitations of the hardware connectivity structure of the annealer. This paper studies a general framework for a decomposition algorithm for NP-hard graph problems aiming to identify an optimal set of vertices. Our generic algorithm allows us to recursively divide an instance until the generated subproblems can be embedded on the quantum annealer hardware and subsequently solved. The framework is applied to the maximum clique and minimum vertex cover problems, and we propose several pruning and reduction techniques to speed up the recursive decomposition. The performance of both algorithms is assessed in a detailed simulation study.


title: M{\o}ller scattering in $2 + 1$ of generalized quantum electrodynamics in the heisenberg picture

abstract: In this paper, we investigate from the framework of generalized electrodynamics the differential cross section of the electron-electron scattering process $e^- e^- \rightarrow e^- e^-$, i.e., M{\o}ller scattering, in $(2+1)$ dimensions in the Heisenberg picture. To this goal, one starts within the stable and unitary framework of planar generalized electrodynamics, instead of Maxwell one. We argue the Haag's theorem strongly suggests the study of the differential cross section in the Heisenberg representation. Afterward, we explore the influence of Podolsky mass cutoff and calculate the differential cross section considering data based on condensed matter systems.


title: Custom silicon technology for SPAD-arrays with red-enhanced sensitivity and low timing jitter

abstract: Single-photon detection is an invaluable tool for many applications ranging from basic research to consumer electronics. In this respect, the Single Photon Avalanche Diode (SPAD) plays a key role in enabling a broad diffusion of these techniques thanks to its remarkable performance, room-temperature operation, and scalability. In this paper we present a silicon technology that allows the fabrication of SPAD-arrays with an unprecedented combination of low timing jitter (95 ps FWHM) and high detection efficiency at red and near infrared wavelengths (peak of 70% at 650 nm, 45% at 800 nm). We discuss the device structure, the fabrication process, and we present a thorough experimental characterization of the fabricated detectors. We think that these long-awaited results can pave the way to new exciting developments in many fields, ranging from quantum optics to single molecule spectroscopy


title: On the adjoint action of the group of symplectic diffeomorphisms

abstract: We study the action of Hamiltonian diffeomorphisms of a compact symplectic manifold ($X,\omega$) on $C^\infty(X)$ and on functions $C^\infty(X)\to \mathbb R$. We describe various properties of invariant convex functions on $C^\infty(X)$. Among other things we show that continuous convex functions $C^\infty(X)\to \mathbb R$ that are invariant under the action are automatically invariant under so called strict rearrangements and they are continuous in the sup norm topology of $C^\infty(X)$; but this is not generally true if the convexity condition is dropped.


title: General structure of Thomas$-$Whitehead gravity

abstract: Thomas-Whitehead (TW) gravity is a projectively invariant model of gravity over a d-dimensional manifold that is intimately related to string theory through reparameterization invariance. Unparameterized geodesics are the ubiquitous structure that ties together string theory and higher dimensional gravitation. This is realized through the projective geometry of Tracy Thomas. The projective connection, due to Thomas and later Whitehead, admits a component that in one dimension is in one-to-one correspondence with the coadjoint elements of the Virasoro algebra. This component is called the diffeomorphism field $\mathcal{D}_{ab }$ in the literature. It also has been shown that in four dimensions, the TW\ action collapses to the Einstein-Hilbert action with cosmological constant when $\mathcal{D}_{ab}$ is proportional to the Einstein metric. These previous results have been restricted to either particular metrics, such as the Polyakov 2D\ metric, or were restricted to coordinates that were volume preserving. In this paper, we review TW gravity and derive the gauge invariant TW action that is explicitly projectively invariant and general coordinate invariant. We derive the covariant field equations for the TW action and show how fermionic fields couple to the gauge invariant theory. The independent fields are the metric tensor $g_{ab}$, the fundamental projective invariant $\Pi^{a}_{\,\,\,bc}$, and the diffeomorphism field $\mathcal D_{ab}$.


title: Measurements of anisotropic g-factors for electrons in InSb nanowire quantum dots

abstract: We have measured the Zeeman splitting of quantum levels in few-electron quantum dots (QDs) formed in narrow bandgap InSb nanowires via the Schottky barriers at the contacts under application of different spatially orientated magnetic fields. The effective g-factor tensor extracted from the measurements is strongly anisotropic and level-dependent, which can be attributed to the presence of strong spin-orbit interaction (SOI) and asymmetric quantum confinement potentials in the QDs. We have demonstrated a successful determination of the principal values and the principal axis orientations of the g-factor tensors in an InSb nanowire QD by the measurements under rotations of a magnetic field in the three orthogonal planes. We also examine the magnetic-field evolution of the excitation spectra in an InSb nanowire QD and extract a SOI strength of $\Delta_{so}\sim 180$ $\mu$eV from an avoided level crossing between a ground state and its neighboring first excited state in the QD.


title: Efficient Transformers: A Survey

abstract: Transformer model architectures have garnered immense interest lately due to their effectiveness across a range of domains like language, vision and reinforcement learning. In the field of natural language processing for example, Transformers have become an indispensable staple in the modern deep learning stack. Recently, a dizzying number of "X-former" models have been proposed - Reformer, Linformer, Performer, Longformer, to name a few - which improve upon the original Transformer architecture, many of which make improvements around computational and memory efficiency. With the aim of helping the avid researcher navigate this flurry, this paper characterizes a large and thoughtful selection of recent efficiency-flavored "X-former" models, providing an organized and comprehensive overview of existing work and models across multiple domains.


title: Limiting Spectral Resolution of a Reflection Grating Made via Electron-Beam Lithography

abstract: Gratings enable dispersive spectroscopy from the X-ray to the optical, and feature prominently in proposed flagships and SmallSats alike. The exacting performance requirements of these future missions necessitate assessing whether the present state-of-the-art in grating manufacture will limit spectrometer performance. In this work, we manufacture a 1.5 mm thick, 1000 nm period at grating using electron-beam lithography (EBL), a promising lithographic technique for patterning gratings for future astronomical observatories. We assess the limiting spectral resolution of this grating by interferometrically measuring the diffracted wavefronts produced in +/-1st order. Our measurements show this grating has a performance of at least R ~ 14,600, and that our assessment is bounded by the error of our interferometric measurement. The impact of EBL stitching error on grating performance is quantifed, and a path to measuring the period error of customized, curved gratings is presented.


title: Variable Binding for Sparse Distributed Representations: Theory and Applications

abstract: Symbolic reasoning and neural networks are often considered incompatible approaches. Connectionist models known as Vector Symbolic Architectures (VSAs) can potentially bridge this gap. However, classical VSAs and neural networks are still considered incompatible. VSAs encode symbols by dense pseudo-random vectors, where information is distributed throughout the entire neuron population. Neural networks encode features locally, often forming sparse vectors of neural activation. Following Rachkovskij (2001); Laiho et al. (2015), we explore symbolic reasoning with sparse distributed representations. The core operations in VSAs are dyadic operations between vectors to express variable binding and the representation of sets. Thus, algebraic manipulations enable VSAs to represent and process data structures in a vector space of fixed dimensionality. Using techniques from compressed sensing, we first show that variable binding between dense vectors in VSAs is mathematically equivalent to tensor product binding between sparse vectors, an operation which increases dimensionality. This result implies that dimensionality-preserving binding for general sparse vectors must include a reduction of the tensor matrix into a single sparse vector. Two options for sparsity-preserving variable binding are investigated. One binding method for general sparse vectors extends earlier proposals to reduce the tensor product into a vector, such as circular convolution. The other method is only defined for sparse block-codes, block-wise circular convolution. Our experiments reveal that variable binding for block-codes has ideal properties, whereas binding for general sparse vectors also works, but is lossy, similar to previous proposals. We demonstrate a VSA with sparse block-codes in example applications, cognitive reasoning and classification, and discuss its relevance for neuroscience and neural networks.


title: Godel Diffeomorphisms

abstract: A basic problem in smooth dynamics is determining if a system can be distinguished from its inverse, i.e., whether a smooth diffeomorphism $T$ is isomorphic to $T^{-1}$. We show that this problem is sufficiently general that asking it for particular choices of $T$ is equivalent to the validity of well-known number theoretic conjectures including the Riemann Hypothesis and Goldbach's conjecture. Further one can produce computable diffeomorphisms $T$ such that the question of whether $T$ is isomorphic to $T^{-1}$ is independent of ZFC.


title: Exploring the toolkit of Jean Bourgain

abstract: Gian-Carlo Rota once asserted that "every mathematician only has a few tricks". The sheer breadth and ingenuity in the work of Jean Bourgain may at first glance appear to be a counterexample to this maxim. However, as we hope to illustrate in this article, even Bourgain relied frequently on a core set of tools, which formed the base from which problems in many disparate mathematical fields could then be attacked. We discuss a selected number of these tools here, and then perform a case study of how an argument in one of Bourgain's papers can be interpreted as a sequential application of several of these tools.


title: Lagrangian Skeleta and Plane Curve Singularities

abstract: We construct closed arboreal Lagrangian skeleta associated to links of isolated plane curve singularities. This yields closed arboreal Lagrangian skeleta for Weinstein pairs and Weinstein 4-manifolds associated to max-tb Legendrian representatives of algebraic links. We provide computations of Legendrian and Weinstein invariants, and discuss the symplectic topological nature of the Fomin-Pylyavskyy-Shustin-Thurston cluster algebra associated to a singularity. Finally, we present a conjectural ADE-classification for Lagrangian fillings of certain Legendrian links and list some related problems.


title: Homological invariants of codimension 2 contact submanifolds

abstract: Codimension 2 contact submanifolds are the natural generalization of transverse knots to contact manifolds of arbitrary dimension. In this paper, we construct new invariants of codimension 2 contact submanifolds. Our main invariant can be viewed as a deformation of the contact homology algebra of the ambient manifold. We describe various applications of these invariants to contact topology. In particular, we exhibit examples of codimension 2 contact embeddings into overtwisted and tight contact manifolds which are formally isotopic but fail to be isotopic through contact embeddings. We also give new obstructions to certain relative symplectic and Lagrangian cobordisms.


title: Predictive Synthesis of Quantum Materials by Probabilistic Reinforcement Learning

abstract: Predictive materials synthesis is the primary bottleneck in realizing new functional and quantum materials. Strategies for synthesis of promising materials are currently identified by time-consuming trial and error approaches and there are no known predictive schemes to design synthesis parameters for new materials. We use reinforcement learning to predict optimal synthesis schedules, i.e. a time-sequence of reaction conditions like temperatures and reactant concentrations, for the synthesis of a prototypical quantum material, semiconducting monolayer MoS$_{2}$, using chemical vapor deposition. The predictive reinforcement leaning agent is coupled to a deep generative model to capture the crystallinity and phase-composition of synthesized MoS$_{2}$ during CVD synthesis as a function of time-dependent synthesis conditions. This model, trained on 10000 computational synthesis simulations, successfully learned threshold temperatures and chemical potentials for the onset of chemical reactions and predicted new synthesis schedules for producing well-sulfidized crystalline and phase-pure MoS$_{2}$, which were validated by computational synthesis simulations. The model can be extended to predict profiles for synthesis of complex structures including multi-phase heterostructures and can also predict long-time behavior of reacting systems, far beyond the domain of the MD simulations used to train the model, making these predictions directly relevant to experimental synthesis.


title: On the investigation of the closure relations for Gamma-Ray Bursts observed by Swift in the post-plateau phase and the GRB fundamental plane

abstract: Gamma-Ray Bursts (GRBs) are the most explosive phenomena in the Universe after the Big Bang. A large fraction of GRB lightcurves (LCs) shows X-ray plateaus. We perform the most comprehensive analysis of all GRBs (with known and unknown redshifts) with plateau emission observed by The Neil Gehrels Swift Observatory from its launch until August 2019. We fit 455 LCs showing a plateau and explore whether these LCs follow closure relations, relations between the temporal and spectral indices of the afterglow, corresponding to 2 distinct astrophysical environments and cooling regimes within the external forward shock (ES) model, and find that the ES model works for the majority of cases. The most favored environments are a constant density interstellar or wind medium with slow cooling. We also confirm the existence of the fundamental plane relation between the rest-frame time and luminosity at the end of the plateau emission and the peak prompt luminosity for this enlarged sample, and test this relation on groups corresponding to the astrophysical environments of our known redshift sample. The plane becomes a crucial discriminant corresponding to these environments in terms of the best fitting parameters and dispersions. Most GRBs for which the closure relations are fulfilled with respect to astrophysical environments have an intrinsic scatter sigma compatible within 1 sigma of that of the Gold GRBs, a subset of long GRBs with relatively flat plateaus. We also find that GRBs satisfying closure relations indicating a fast cooling regime have a lower sigma than ever previously found in literature.


title: Explicit construction of cofree precoalgebras and coalgebras

abstract: In this paper, we will explicitly construct cofree coalgebras, by first constructing cofree precoalgebras (namely those not necessarily coassociative or counital). Our approach does not impose any condition to the coefficient ring, which means, for example, it does not need to be a field. The main technique is to generate homomorphisms from certain type of infinite trees so that those constitute cofree precoalgebras.


title: Leveraging Domain Knowledge using Machine Learning for Image Compression in Internet-of-Things

abstract: The emergent ecosystems of intelligent edge devices in diverse Internet of Things (IoT) applications, from automatic surveillance to precision agriculture, increasingly rely on recording and processing variety of image data. Due to resource constraints, e.g., energy and communication bandwidth requirements, these applications require compressing the recorded images before transmission. For these applications, image compression commonly requires: (1) maintaining features for coarse-grain pattern recognition instead of the high-level details for human perception due to machine-to-machine communications; (2) high compression ratio that leads to improved energy and transmission efficiency; (3) large dynamic range of compression and an easy trade-off between compression factor and quality of reconstruction to accommodate a wide diversity of IoT applications as well as their time-varying energy/performance needs. To address these requirements, we propose, MAGIC, a novel machine learning (ML) guided image compression framework that judiciously sacrifices visual quality to achieve much higher compression when compared to traditional techniques, while maintaining accuracy for coarse-grained vision tasks. The central idea is to capture application-specific domain knowledge and efficiently utilize it in achieving high compression. We demonstrate that the MAGIC framework is configurable across a wide range of compression/quality and is capable of compressing beyond the standard quality factor limits of both JPEG 2000 and WebP. We perform experiments on representative IoT applications using two vision datasets and show up to 42.65x compression at similar accuracy with respect to the source. We highlight low variance in compression rate across images using our technique as compared to JPEG 2000 and WebP.


title: An introduction to the Bernoulli function

abstract: We explore a variant of the zeta function interpolating the Bernoulli numbers based on an integral representation suggested by J. Jensen. The Bernoulli function $\operatorname{B}(s, v) = - s\, \zeta(1-s, v)$ can be introduced independently of the zeta function if it is based on a formula first given by Jensen in 1895. We examine the functional equation of $\operatorname{B}(s, v)$ and its representation by the Riemann $\zeta$ and $\xi$ function, and recast classical results of Hadamard, Worpitzky, and Hasse in terms of $\operatorname{B}(s, v).$ The extended Bernoulli function defines the Bernoulli numbers for odd indices basing them on rational numbers studied by Euler in 1735 that underlie the Euler and Andr\'{e} numbers. The Euler function is introduced as the difference between values of the Hurwitz-Bernoulli function. The Andr\'{e} function and the Seki function are the unsigned versions of the extended Euler resp. Bernoulli function.


title: Functional Inequalities involving Nonlocal Operators on Complete Riemannian Manifolds and Their Applications to The Fractional Porous Medium Equation

abstract: The objective of this paper is twofold. First, we conduct a careful study of various functional inequalities involving the fractional Laplacian operators, including nonlocal Sobolev-Poincar\'e, Nash, Super Poincar\'e and logarithmic Sobolev type inequalities, on complete Riemannian manifolds satisfying some mild geometric assumptions. Second, based on the derived nonlocal functional inequalities, we analyze the asymptotic behavior of the solution to the fractional porous medium equation, $\partial_t u +(-\Delta)^\sigma (|u|^{m-1}u )=0 $ with $m>0$ and $\sigma\in (0,1)$. In addition, we establish the global well-posedness of the equation on an arbitrary complete Riemannian manifold.


title: Escape from an attractor generated by recurrent exit

abstract: Kramer's theory of activation over a potential barrier consists in computing the mean exit time from the boundary of a basin of attraction of a randomly perturbed dynamical system. Here we report that for some systems, crossing the boundary is not enough, because stochastic trajectories return inside the basin with a high probability a certain number of times before escaping far away. This situation is due to a shallow potential. We compute the mean and distribution of escape times and show how this result explains the large distribution of interburst durations in neuronal networks.

categories: cond-mat.stat-mech math.PR nlin.CD q-bio.CB


title: Orbital stability of KdV multisolitons in $H^{-1}$

abstract: We prove that multisoliton solutions of the Korteweg--de Vries equation are orbitally stable in $H^{-1}(\mathbb{R})$. We introduce a variational characterization of multisolitons that remains meaningful at such low regularity and show that all optimizing sequences converge to the manifold of multisolitons. The proximity required at the initial time is uniform across the entire manifold of multisolitons; this had not been demonstrated previously, even in $H^1$.


title: Distributed Mirror Descent with Integral Feedback: Asymptotic Convergence Analysis of Continuous-time Dynamics

abstract: This work addresses distributed optimization, where a network of agents wants to minimize a global strongly convex objective function. The global function can be written as a sum of local convex functions, each of which is associated with an agent. We propose a continuous-time distributed mirror descent algorithm that uses purely local information to converge to the global optimum. Unlike previous work on distributed mirror descent, we incorporate an integral feedback in the update, allowing the algorithm to converge with a constant step-size when discretized. We establish the asymptotic convergence of the algorithm using Lyapunov stability analysis. We further illustrate numerical experiments that verify the advantage of adopting integral feedback for improving the convergence rate of distributed mirror descent.


title: Interplay between complex symmetry and Koenigs eigenfunctions

abstract: We investigate the relationship between the complex symmetry of composition operators $C_{\phi}f=f\circ \phi$ induced on the classical Hardy space $H^2(\mathbb{D})$ by an analytic self-map $\phi$ of the open unit disk $\mathbb{D}$ and its Koenigs eigenfunction. A generalization of orthogonality known as conjugate-orthogonality will play a key role in this work. We show that if $\phi$ is a Schr\"{o}der map (fixes a point $a\in \mathbb{D}$ with $0<|\phi'(a)|<1$) and $\sigma$ is its Koenigs eigenfunction, then $C_{\phi}$ is complex symmetric if and only if $(\sigma^n)_{n\in \mathbb{N}}$ is complete and conjugate-orthogonal in $H^2(\mathbb{D})$. We study the conjugate-orthogonality of Koenigs sequences with some concrete examples. We use these results to show that commutants of complex symmetric composition operators with Schr\"{o}der symbols consist entirely of complex symmetric operators.


title: A pair of Jovian Trojans at the L4 Lagrange point

abstract: Asteroid pairs, two objects that are not gravitationally bound to one another, but share a common origin, have been discovered in the Main belt and Hungaria populations. Such pairs are of major interest, as the study of their evolution under a variety of dynamical influences can indicate the time since the pair was created. To date, no asteroid pairs have been found in the Jovian Trojans, despite the presence of several binaries and collisional families in the population. The search for pairs in the Jovian Trojan population is of particular interest, given the importance of the Trojans as tracers of planetary migration during the Solar system's youth. Here we report a discovery of the first pair, (258656) 2002~ES$_{76}$ and 2013~CC$_{41}$, in the Jovian Trojans. The two objects are approximately the same size and are located very close to the L4 Lagrange point. Using numerical integrations, we find that the pair is at least $360$~Myr old, though its age could be as high as several Gyrs. The existence of the (258656) 2002~ES$_{76}$--2013~CC$_{41}$ pair implies there could be many such pairs scattered through the Trojan population. Our preferred formation mechanism for the newly discovered pair is through the dissociation of an ancient binary system, triggered by a sub-catastrophic impact, but we can not rule out rotation fission of a single object driven by YORP torques. A by-product of our work is an up-to-date catalog of Jovian Trojan proper elements, which we have made available for further studies.


title: Stop the Clock: Are Timeout Effects Real?

abstract: Timeout is a short interruption during games used to communicate a change in strategy, to give the players a rest or to stop a negative flow in the game. Whatever the reason, coaches expect an improvement in their team's performance after a timeout. But how effective are these timeouts in doing so? The simple average of the differences between the scores before and after the timeouts has been used as evidence that there is an effect and that it is substantial. We claim that these statistical averages are not proper evidence and a more sound approach is needed. We applied a formal causal framework using a large dataset of official NBA play-by-play tables and drew our assumptions about the data generation process in a causal graph. Using different matching techniques to estimate the causal effect of timeouts, we concluded that timeouts have no effect on teams' performances. Actually, since most timeouts are called when the opposing team is scoring more frequently, the moments that follow resemble an improvement in the team's performance but are just the natural game tendency to return to its average state. This is another example of what statisticians call the regression to the mean phenomenon.


title: Universal Toeplitz operators on the Hardy space over the polydisk

abstract: The Invariant Subspace Problem (ISP) for Hilbert spaces asks if every bounded linear operator has a non-trivial closed invariant subspace. Due to the existence of universal operators (in the sense of Rota), the ISP may be solved by describing the invariant subspaces of these operators alone. We characterize all anaytic Toeplitz operators $T_\phi$ on the Hardy space $H^2(\mathbb{D}^n)$ over the polydisk $\mathbb{D}^n$ for $n>1$ whose adjoints satisfy the Caradus criterion for universality, that is, when $T_\phi^*$ is surjective and has infinite dimensional kernel. In particular if $\phi$ in a non-constant inner function on $\mathbb{D}^n$, or a polynomial in the ring $\mathbb{C}[z_1,\ldots,z_n]$ that has zeros in $\mathbb{D}^n$ but is zero-free on $\mathbb{T}^n$, then $T_\phi^*$ is universal for $H^2(\mathbb{D}^n)$. The analogs of these results for $n=1$ are not true.


title: Modernizing Archimedes' Construction of $\pi$

abstract: In his famous work, "Measurement of a Circle," Archimedes described a procedure for measuring both the circumference of a circle and the area it bounds. Implicit in his work is the idea that his procedure defines these quantities. Modern approaches for defining $\pi$ eschew his method and instead use arguments that are easier to justify, but they involve ideas that are not elementary. This paper makes Archimedes' measurement procedure rigorous from a modern perspective. In so doing, it brings a rigorous and geometric treatment of the differential properties of the trigonometric functions into the purview of an introductory calculus course.


title: Effect of controlled point-like disorder induced by 2.5 MeV electron irradiation on nematic resistivity anisotropy of hole-doped (Ba,K)Fe$_2$As$_2$

abstract: In-plane anisotropy of electrical resistivity was studied in samples of the hole-doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$ in the composition range $0.21 \leq x \leq 0.26$ where anisotropy changes sign. Low-temperature ($\sim$20~K) irradiation with relativistic 2.5 MeV electrons was used to control the level of disorder and residual resistivity of the samples. Modification of the stress-detwinning technique enabled measurements of the same samples before and after irradiation, leading to conclusion of anisotropic character of predominantly inelastic scattering processes. Our main finding is that the resistivity anisotropy is of the same sign irrespective of residual resistivity, and remains the same in the orthorhombic $C_2$ phase above the re-entrant tetragonal transition. Unusual $T$-linear dependence of the anisotropy $\Delta \rho \equiv \rho_a(T)-\rho_b(T)$ is found in pristine samples with $x=$0.213 and $x=$0.219, without similar signatures in either $\rho_a(T)$ or $\rho_b(T)$. We show that this feature can be reproduced by a phenomenological model of R.~M.~Fernandes {\it et al.} Phys. Rev. Lett. {\bf 107},217002 (2011). We speculate that onset of fluctuations of nematic order on approaching the instability towards the re-entrant tetragonal phase contributes to this unusual dependence.


title: Benchmarking Machine Learning Techniques with Di-Higgs Production at the LHC

abstract: Many domains of high energy physics analysis are starting to explore machine learning techniques. Powerful methods can be used to identify and measure rare processes from previously insurmountable backgrounds. One of the most profound Standard Model signatures still to be discovered at the LHC is the pair production of Higgs bosons through the Higgs self-coupling. The small cross section of this process makes detection very difficult even for the decay channel with the largest branching fraction ($hh\rightarrow b\bar{b}b\bar{b}$). This paper benchmarks a variety of approaches (boosted decision trees, various neural network architectures, semi-supervised algorithms) against one another to catalog a few of the various techniques available to high energy physicists as the era of the HL-LHC approaches.


title: Buffon's Problem determines Gaussian Curvature in three Geometries

abstract: A version of the classical Buffon problem in the plane naturally extends to the setting of any Riemannian surface with constant Gaussian curvature. The Buffon probability determines a Buffon deficit. The relationship between Gaussian curvature and the Buffon deficit is similar to the relationship that the Bertrand-Diguet-Puiseux Theorem establishes between Gaussian curvature and both circumference and area deficits.


title: Analysis of Models for Decentralized and Collaborative AI on Blockchain

abstract: Machine learning has recently enabled large advances in artificial intelligence, but these results can be highly centralized. The large datasets required are generally proprietary; predictions are often sold on a per-query basis; and published models can quickly become out of date without effort to acquire more data and maintain them. Published proposals to provide models and data for free for certain tasks include Microsoft Research's Decentralized and Collaborative AI on Blockchain. The framework allows participants to collaboratively build a dataset and use smart contracts to share a continuously updated model on a public blockchain. The initial proposal gave an overview of the framework omitting many details of the models used and the incentive mechanisms in real world scenarios. In this work, we evaluate the use of several models and configurations in order to propose best practices when using the Self-Assessment incentive mechanism so that models can remain accurate and well-intended participants that submit correct data have the chance to profit. We have analyzed simulations for each of three models: Perceptron, Na\"ive Bayes, and a Nearest Centroid Classifier, with three different datasets: predicting a sport with user activity from Endomondo, sentiment analysis on movie reviews from IMDB, and determining if a news article is fake. We compare several factors for each dataset when models are hosted in smart contracts on a public blockchain: their accuracy over time, balances of a good and bad user, and transaction costs (or gas) for deploying, updating, collecting refunds, and collecting rewards. A free and open source implementation for the Ethereum blockchain and simulations written in Python is provided at https://github.com/microsoft/0xDeCA10B. This version has updated gas costs using newer optimizations written after the original publication.


title: Simultaneous Denoising and Motion Estimation for Low-dose Gated PET using a Siamese Adversarial Network with Gate-to-Gate Consistency Learning

abstract: Gating is commonly used in PET imaging to reduce respiratory motion blurring and facilitate more sophisticated motion correction methods. In the applications of low dose PET, however, reducing injection dose causes increased noise and reduces signal-to-noise ratio (SNR), subsequently corrupting the motion estimation/correction steps, causing inferior image quality. To tackle these issues, we first propose a Siamese adversarial network (SAN) that can efficiently recover high dose gated image volume from low dose gated image volume. To ensure the appearance consistency between the recovered gated volumes, we then utilize a pre-trained motion estimation network incorporated into SAN that enables the constraint of gate-to-gate (G2G) consistency. With high-quality recovered gated volumes, gate-to-gate motion vectors can be simultaneously outputted from the motion estimation network. Comprehensive evaluations on a low dose gated PET dataset of 29 subjects demonstrate that our method can effectively recover the low dose gated PET volumes, with an average PSNR of 37.16 and SSIM of 0.97, and simultaneously generate robust motion estimation that could benefit subsequent motion corrections.


title: Accelerating Bayesian microseismic event location with deep learning

abstract: We present a series of new open source deep learning algorithms to accelerate Bayesian full waveform point source inversion of microseismic events. Inferring the joint posterior probability distribution of moment tensor components and source location is key for rigorous uncertainty quantification. However, the inference process requires forward modelling of microseismic traces for each set of parameters explored by the sampling algorithm, which makes the inference very computationally intensive. In this paper we focus on accelerating this process by training deep learning models to learn the mapping between source location and seismic traces, for a given 3D heterogeneous velocity model, and a fixed isotropic moment tensor for the sources. These trained emulators replace the expensive solution of the elastic wave equation in the inference process. We compare our results with a previous study that used emulators based on Gaussian Processes to invert microseismic events. We show that all of our models provide more accurate predictions and $\sim 100$ times faster predictions than the method based on Gaussian Processes, and a $\mathcal{O}(10^5)$ speed-up factor over a pseudo-spectral method for waveform generation. For example, a 2-s long synthetic trace can be generated in $\sim 10$ ms on a common laptop processor, instead of $\sim$ 1 hr using a pseudo-spectral method on a high-profile Graphics Processing Units card. We also show that our inference results are in excellent agreement with those obtained from traditional location methods based on travel time estimates. The speed, accuracy and scalability of our open source deep learning models pave the way for extensions of these emulators to generic source mechanisms and application to joint Bayesian inversion of moment tensor components and source location using full waveforms.


title: Transition region contribution to AIA observations in the context of coronal heating

abstract: We investigate the relative contributions from the transition region and corona of coronal loops observed by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO). Using EBTEL (Enthalpy-Based Thermal Evolution of Loops) hydrodynamic simulations, we model loops with multiple lengths and energy fluxes heated randomly by events drawn from power-law distributions with different slopes and minimum delays between events to investigate how each of these parameters influences observable loop properties. We generate AIA intensities from the corona and transition region for each realization. The variations within and between models generated with these different parameters illustrate the sensitivity of narrowband imaging to the details of coronal heating. We then analyze the transition region and coronal emission from a number of observed active regions and find broad agreement with the trends in the models. In both models and observations, the transition region brightness is significant, often greater than the coronal brightness in all six "coronal" AIA channels. We also identify an inverse relationship, consistent with heating theories, between the slope of the differential emission measure (DEM) coolward of the peak temperature and the observed ratio of coronal to transition region intensity. These results highlight the use of narrowband observations and the importance of properly considering the transition region in investigations of coronal heating.


title: A method for finding anomalous astronomical light curves and their analogs

abstract: Our understanding of the Universe has profited from deliberate, targeted studies of known phenomena, as well as from serendipitous, unexpected discoveries, such as the discovery of a complex variability pattern in the direction of KIC 8462852 (Boyajian's star). Upcoming surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will explore the parameter space of astrophysical transients at all time scales, and offer the opportunity to discover even more extreme examples of unexpected phenomena. We investigate strategies to identify novel objects and to contextualize them within large time-series data sets in order to facilitate the discovery of new classes of objects, as well as the physical interpretation of their anomalous nature. We develop a method that combines tree-based and manifold-learning algorithms for anomaly detection in order to perform two tasks: 1) identify and rank anomalous objects in a time-domain dataset; and 2) group those anomalies according to their similarity in order to identify analogs. We achieve the latter by combining an anomaly score from a tree-based method with a dimensionality manifold-learning reduction strategy. Clustering in the reduced space allows for the successful identification of anomalies and analogs. We also assess the impact of pre-processing and feature engineering schemes and investigate the astrophysical nature of the objects that our models identify as anomalous by augmenting the Kepler data with Gaia color and luminosity information. We find that multiple models, used in combination, are a promising strategy to identify novel light curves and light curve families.


title: Introduction to a novel T2 relaxation analysis method SAME-ECOS: Spectrum Analysis for Multiple Exponentials via Experimental Condition Oriented Simulation

abstract: We propose a novel T2 relaxation data analysis method which we have named spectrum analysis for multiple exponentials via experimental condition oriented simulation (SAME-ECOS). SAME-ECOS, which was developed based on a combination of information theory and machine learning neural network algorithms, is tailored for different MR experimental conditions, decomposing multi-exponential decay data into T2 spectra, which had been considered an ill-posed problem using conventional fitting algorithms, including the commonly used non-negative least squares (NNLS) method. Our results demonstrated that, compared with NNLS, the simulation-derived SAME-ECOS model yields much more reliable T2 spectra in a dramatically shorter time, increasing the feasibility of multi-component T2 decay analysis in clinical settings.


title: New complex network building methodology for High Level Classification based on attribute-attribute interaction

abstract: High-level classification algorithms focus on the interactions between instances. These produce a new form to evaluate and classify data. In this process, the core is the complex network building methodology because it determines the metrics to be used for classification. The current methodologies use variations of kNN to produce these graphs. However, this technique ignores some hidden pattern between attributes and require normalization to be accurate. In this paper, we propose a new methodology for network building based on attribute-attribute interactions that do not require normalization and capture the hidden patterns of the attributes. The current results show us that could be used to improve some current high-level techniques.


title: Superconductivity suppression in disordered films: Interplay of two-dimensional diffusion and three-dimensional ballistics

abstract: Suppression of the critical temperature in homogeneously disordered superconducting films is a consequence of the disorder-induced enhancement of Coulomb repulsion. We demonstrate that for the majority of thin films studied now this effect cannot be completely explained in the assumption of two-dimensional diffusive nature of electrons motion. The main contribution to the $T_c$ suppression arises from the correction to the electron-electron interaction constant coming from small scales of the order of the Fermi wavelength that leads to the critical temperature shift $\delta T_c/T_{c0} \sim - 1/k_Fl$, where $k_F$ is the Fermi momentum and $l$ is the mean free path. Thus almost for all superconducting films that follow the fermionic scenario of $T_c$ suppression with decreasing the film thickness, this effect is caused by the proximity to the three-dimensional Anderson localization threshold and is controlled by the parameter $k_F l$ rather than the sheet resistance of the film.

categories: cond-mat.supr-con cond-mat.dis-nn cond-mat.mes-hall cond-mat.mtrl-sci


title: Private data sharing between decentralized users through the privGAN architecture

abstract: More data is almost always beneficial for analysis and machine learning tasks. In many realistic situations however, an enterprise cannot share its data, either to keep a competitive advantage or to protect the privacy of the data sources, the enterprise's clients for example. We propose a method for data owners to share synthetic or fake versions of their data without sharing the actual data, nor the parameters of models that have direct access to the data. The method proposed is based on the privGAN architecture where local GANs are trained on their respective data subsets with an extra penalty from a central discriminator aiming to discriminate the origin of a given fake sample. We demonstrate that this approach, when applied to subsets of various sizes, leads to better utility for the owners than the utility from their real small datasets. The only shared pieces of information are the parameter updates of the central discriminator. The privacy is demonstrated with white-box attacks on the most vulnerable elments of the architecture and the results are close to random guessing. This method would apply naturally in a federated learning setting.


title: Negativity vs. Purity and Entropy in Witnessing Entanglement

abstract: In this work, we show that while all measures of mixedness may be used to witness entanglement, all such entangled states must have a negative partial transpose (NPT). Though computing the negativity of the partial transpose scales well at high dimension, it relies on knowing the complete quantum state, which does not. To address this, we compare different measures of mixedness over uniform ensembles of joint quantum states at varying dimension to gauge their relative success in witnessing entanglement. In doing so, we find that comparing joint and marginal purities is overwhelmingly more successful at high dimension at identifying entanglement than comparing joint and marginal von Neumann entropies, in spite of requiring fewer resources. We conclude by showing how our results impact the fundamental relationship between correlation and entanglement and related witnesses.


title: Dynamic polarization and plasmons in kekul\'e-patterned graphene: Signatures of the broken valley degeneracy

abstract: The dynamic polarization for kekul\'e-patterned graphene is studied within the Random Phase Approximation (RPA). It is shown how the breaking of the valley degeneracy by the lattice modulation is manifested through the dielectric spectrum, the plasmonic dispersion, the static screening and the optical conductivity. The valley-dependent splitting of the Fermi velocities due to the kekul\'e distortion leads to a similar splitting in the dielectric spectrum of graphene, introducing new characteristic frequencies, which are given in terms of the valley-coupling amplitude. The valley coupling also splits the plasmonic dispersion, introducing a second branch in the Landau damping region. Finally, the signatures of the broken valley degeneracy in the optical conductivity are studied. The characteristic step-like spectrum of graphene is split into two half steps due to the onset of absorption in each valley occurring at different characteristic frequencies. Also, it was found an absorption phenomenon where a resonance peak related to intervalley transport emerges at a "beat frequency", determined by the difference between the characteristic frequencies of each valley. Some of these mechanisms are expected to be present in other space-modulated 2D materials and suggest how optical or electrical response measurements can be suitable to detect spatial modulation.


title: Qutrit-inspired Fully Self-supervised Shallow Quantum Learning Network for Brain Tumor Segmentation

abstract: Classical self-supervised networks suffer from convergence problems and reduced segmentation accuracy due to forceful termination. Qubits or bi-level quantum bits often describe quantum neural network models. In this article, a novel self-supervised shallow learning network model exploiting the sophisticated three-level qutrit-inspired quantum information system referred to as Quantum Fully Self-Supervised Neural Network (QFS-Net) is presented for automated segmentation of brain MR images. The QFS-Net model comprises a trinity of a layered structure of qutrits inter-connected through parametric Hadamard gates using an 8-connected second-order neighborhood-based topology. The non-linear transformation of the qutrit states allows the underlying quantum neural network model to encode the quantum states, thereby enabling a faster self-organized counter-propagation of these states between the layers without supervision. The suggested QFS-Net model is tailored and extensively validated on Cancer Imaging Archive (TCIA) data set collected from Nature repository and also compared with state of the art supervised (U-Net and URes-Net architectures) and the self-supervised QIS-Net model. Results shed promising segmented outcome in detecting tumors in terms of dice similarity and accuracy with minimum human intervention and computational resources.


title: An all-digital associated particle imaging system for the 3D determination of isotopic distributions

abstract: Associated particle imaging (API) is a non-destructive nuclear technique for the 3D determination of isotopic distributions. By detecting the alpha particle associated with the emitted neutron in the deuterium-tritium fusion reaction with a position- and time-resolving detector, the direction of the 14.1 MeV neutron and its time of emission can be determined. Employing this method, isotope characteristic gamma rays emitted in inelastic neutron scattering events can be correlated with the neutron interaction location. An API system consisting of a sealed-type neutron generator, gamma detectors, and a position-sensitive alpha detector was designed, constructed, and characterized. The system was tested with common soil elements and shown to be sensitive to 12C, 16O, 28Si, 27Al, and 56Fe. New aspects of our approach are the use of a yttrium-aluminum-perovskite (YAP) scintillator, using a sapphire window instead of a fiber-optic faceplate for light transport to the photomultiplier, and the all-digital data acquisition system. We present a description of the system with simulations and experimental results that show a position resolution on the alpha detector of 1 mm, a depth resolution using a LaBr3 detector of 6.2 cm, and an angular resolution of 4.5 degrees. Additionally, we present single-element gamma response measurements for the elements mentioned above together with a comparison to Monte Carlo simulations (MCNP6).


title: Infinite series asymptotic expansions for decaying solutions of dissipative differential equations with non-smooth nonlinearity

abstract: We study the precise asymptotic behavior of a non-trivial solution that converges to zero, as time tends to infinity, of dissipative systems of nonlinear ordinary differential equations. The nonlinear term of the equations may not possess a Taylor series expansion about the origin. This absence technically cripples previous proofs in establishing an asymptotic expansion, as an infinite series, for such a decaying solution. In the current paper, we overcome this limitation and obtain an infinite series asymptotic expansion, as time goes to infinity. This series expansion provides large time approximations for the solution with the errors decaying exponentially at any given rates. The main idea is to shift the center of the Taylor expansions for the nonlinear term to a non-zero point. Such a point turns out to come from the non-trivial asymptotic behavior of the solution, which we prove by a new and simple method. Our result applies to different classes of non-linear equations that have not been dealt with previously.


title: Joint Subgraph-to-Subgraph Transitions -- Generalizing Triadic Closure for Powerful and Interpretable Graph Modeling

abstract: We generalize triadic closure, along with previous generalizations of triadic closure, under an intuitive umbrella generalization: the Subgraph-to-Subgraph Transition (SST). We present algorithms and code to model graph evolution in terms of collections of these SSTs. We then use the SST framework to create link prediction models for both static and temporal, directed and undirected graphs which produce highly interpretable results. Quantitatively, our models match out-of-the-box performance of state of the art graph neural network models, thereby validating the correctness and meaningfulness of our interpretable results.


title: Brieskorn module and Center conditions: pull-back of differential equations in projective space

abstract: The moduli space of algebraic foliations on P2 of a fixed degree and with a center singularity has many irreducible components. We find a basis of the Brieskorn module defined for a rational function and prove that set of pull-back foliations forms an irreducible component of the moduli space. The main tools are the Picard-Lefschetz theory of a rational function in two variables, period integrals, and the Brieskorn module.


title: Carbon Diffusion in Concentrated Fe-C Glasses

abstract: By combining atomistic simulations with a detailed analysis of individual atomic hops, we show that the diffusion of carbon in a binary Fe-C glass exhibits strong (anti-)correlations and is largely determined by the local environment. Higher local carbon concentrations lead to slower atomic mobility. Our results help explain the increasing stability of Fe-C (and other similar metal-metalloid glasses) against crystallization with increasing carbon concentration.


title: Truth or Square: Aspect Ratio Biases Recall of Position Encodings

abstract: Bar charts are among the most frequently used visualizations, in part because their position encoding leads them to convey data values precisely. Yet reproductions of single bars or groups of bars within a graph can be biased. Curiously, some previous work found that this bias resulted in an overestimation of reproduced data values, while other work found an underestimation. Across three empirical studies, we offer an explanation for these conflicting findings: this discrepancy is a consequence of the differing aspect ratios of the tested bar marks. Viewers are biased to remember a bar mark as being more similar to a prototypical square, leading to an overestimation of bars with a wide aspect ratio, and an underestimation of bars with a tall aspect ratio. Experiments 1 and 2 showed that the aspect ratio of the bar marks indeed influenced the direction of this bias. Experiment 3 confirmed that this pattern of misestimation bias was present for reproductions from memory, suggesting that this bias may arise when comparing values across sequential displays or views. We describe additional visualization designs that might be prone to this bias beyond bar charts (e.g., Mekko charts and treemaps), and speculate that other visual channels might hold similar biases toward prototypical values.


title: Decoupling of dipolar and hydrophobic motions in biological membranes

abstract: Cells use homeostatic mechanisms to maintain an optimal composition of distinct types of phospholipids in cellular membranes. The hydrophilic dipolar layer at the membrane interface, composed of phospholipid headgroups, regulates the interactions between cell membranes and incoming molecules, nanoparticles, and viruses. On the other hand, the membrane hydrophobic core determines membrane thickness and forms an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of lipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy methodology with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm where phospholipid headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.


title: Controllable neural text-to-speech synthesis using intuitive prosodic features

abstract: Modern neural text-to-speech (TTS) synthesis can generate speech that is indistinguishable from natural speech. However, the prosody of generated utterances often represents the average prosodic style of the database instead of having wide prosodic variation. Moreover, the generated prosody is solely defined by the input text, which does not allow for different styles for the same sentence. In this work, we train a sequence-to-sequence neural network conditioned on acoustic speech features to learn a latent prosody space with intuitive and meaningful dimensions. Experiments show that a model conditioned on sentence-wise pitch, pitch range, phone duration, energy, and spectral tilt can effectively control each prosodic dimension and generate a wide variety of speaking styles, while maintaining similar mean opinion score (4.23) to our Tacotron baseline (4.26).


title: On the optimal certification of von Neumann measurements

abstract: In this report we study certification of quantum measurements, which can be viewed as the extension of quantum hypotheses testing. This extension involves also the study of the input state and the measurement procedure. Here, we will be interested in two-point (binary) certification scheme in which the null and alternative hypotheses are single element sets. Our goal is to minimize the probability of the type II error given some fixed statistical significance. In this report, we begin with studying the two-point certification of pure quantum states and unitary channels to later use them to prove our main result, which is the certification of von Neumann measurements in single-shot and parallel scenarios. From our main result follow the conditions when two pure states, unitary operations and von Neumann measurements cannot be distinguished perfectly but still can be certified with a given statistical significance. Moreover, we show the connection between the certification of quantum channels or von Neumann measurements and the notion of $q$-numerical range.


title: Open Cluster Chemical Homogeneity Throughout the Milky Way

abstract: The chemical homogeneity of surviving stellar clusters contains important clues about interstellar medium (ISM) mixing efficiency, star formation, and the enrichment history of the Galaxy. Existing measurements in a handful of open clusters suggest homogeneity in several elements at the 0.03 dex level. Here we present (i) a new cluster member catalog based only on APOGEE radial velocities and Gaia-DR2 proper motions, (ii) improved abundance uncertainties for APOGEE cluster members, and (iii) the dependence of cluster homogeneity on Galactic and cluster properties, using abundances of eight elements from the APOGEE survey for ten high-quality clusters. We find that cluster homogeneity is uncorrelated with Galactocentric distance, |Z|, age, and metallicity. However, velocity dispersion, which is a proxy for cluster mass, is positively correlated with intrinsic scatter at relatively high levels of significance for [Ca/Fe] and [Mg/Fe]. We also see a possible positive correlation at a low level of significance for [Ni/Fe], [Si/Fe], [Al/Fe], and [Fe/H], while [Cr/Fe] and [Mn/Fe] are uncorrelated. The elements that show a correlation with velocity dispersion are those that are predominantly produced by core-collapse supernovae (CCSNe). However, the small sample size and relatively low correlation significance highlight the need for follow-up studies. If borne out by future studies, these findings would suggest a quantitative difference between the correlation lengths of elements produced predominantly by Type~Ia SNe versus CCSNe, which would have implications for Galactic chemical evolution models and the feasibility of chemical tagging.


title: Spatio-Temporal Top-k Similarity Search for Trajectories in Graphs

abstract: We study the problem of finding the $k$ most similar trajectories to a given query trajectory. Our work is inspired by the work of Grossi et al. [6] that considers trajectories as walks in a graph. Each visited vertex is accompanied by a time-interval. Grossi et al. define a similarity function that captures temporal and spatial aspects. We improve this similarity function to derive a new spatio-temporal distance function for which we can show that a specific type of triangle inequality is satisfied. This distance function is the basis for our index structures, which can be constructed efficiently, need only linear memory, and can quickly answer queries for the top-$k$ most similar trajectories. Our evaluation on real-world and synthetic data sets shows that our algorithms outperform the baselines with respect to indexing time by several orders of magnitude while achieving similar or better query time and quality of results.


title: Factorization at Subleading Power and Endpoint Divergences in $h\to\gamma\gamma$ Decay: II. Renormalization and Scale Evolution

abstract: Building on the recent derivation of a bare factorization theorem for the $b$-quark induced contribution to the $h\to\gamma\gamma$ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where $\lambda=m_b/M_h\ll 1$ in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that "endpoint regularization" does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $\alpha{\hspace{0.3mm}}\alpha_s^2{\hspace{0.3mm}} L^k$ in the three-loop decay amplitude, where $L=\ln(-M_h^2/m_b^2)$ and $k=6,5,4,3$. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $\sim\alpha{\hspace{0.3mm}}\alpha_s^n{\hspace{0.3mm}} L^{2n+1}$.


title: Healthcare Cost Prediction: Leveraging Fine-grain Temporal Patterns

abstract: Objective: To design and assess a method to leverage individuals' temporal data for predicting their healthcare cost. To achieve this goal, we first used patients' temporal data in their fine-grain form as opposed to coarse-grain form. Second, we devised novel spike detection features to extract temporal patterns that improve the performance of cost prediction. Third, we evaluated the effectiveness of different types of temporal features based on cost information, visit information and medical information for the prediction task. Materials and methods: We used three years of medical and pharmacy claims data from 2013 to 2016 from a healthcare insurer, where the first two years were used to build the model to predict the costs in the third year. To prepare the data for modeling and prediction, the time series data of cost, visit and medical information were extracted in the form of fine-grain features (i.e., segmenting each time series into a sequence of consecutive windows and representing each window by various statistics such as sum). Then, temporal patterns of the time series were extracted and added to fine-grain features using a novel set of spike detection features (i.e., the fluctuation of data points). Gradient Boosting was applied on the final set of extracted features. Moreover, the contribution of each type of data (i.e., cost, visit and medical) was assessed. Conclusions: Leveraging fine-grain temporal patterns for healthcare cost prediction significantly improves prediction performance. Enhancing fine-grain features with extraction of temporal cost and visit patterns significantly improved the performance. However, medical features did not have a significant effect on prediction performance. Gradient Boosting outperformed all other prediction models.


title: Pilot: Winner of the Human-Agent Negotiation Challenge at IJCAI 2020

abstract: This document describes our agent Pilot, winner of the Human-Agent Negotiation Challenge at ANAC, IJCAI 2020. Pilot is a virtual human that participates in a sequence of three negotiations with a human partner. Our system is based on the Interactive Arbitration Guide Online (IAGO) negotiation framework. We leverage prior Affective Computing and Psychology research in negotiations to guide various key principles that define the behavior and personality of our agent.


title: Analysis of Random Access in NB-IoT Networks with Three Coverage Enhancement Groups: A Stochastic Geometry Approach

abstract: NarrowBand-Internet of Things (NB-IoT) is a new 3GPP radio access technology designed to provide better coverage for Low Power Wide Area (LPWA) networks. To provide reliable connections with extended coverage, a repetition transmission scheme and up to three Coverage Enhancement (CE) groups are introduced into NB-IoT during both Random Access CHannel (RACH) procedure and data transmission procedure, where each CE group is configured with different repetition values and transmission resources. To characterize the RACH performance of the NB-IoT network with three CE groups, this paper develops a novel traffic-aware spatio-temporal model to analyze the RACH success probability, where both the preamble transmission outage and the collision events of each CE group jointly determine the traffic evolution and the RACH success probability. Based on this analytical model, we derive the analytical expression for the RACH success probability of a randomly chosen IoT device in each CE group over multiple time slots with different RACH schemes, including baseline, back-off (BO), access class barring (ACB), and hybrid ACB and BO schemes (ACB&BO). Our results have shown that the RACH success probabilities of the devices in three CE groups outperform that of a single CE group network but not for all the groups, which is affected by the choice of the categorizing parameters.This mathematical model and analytical framework can be applied to evaluate the performance of multiple group users of other networks with spatial separations.


title: Learning Hidden Patterns from Patient Multivariate Time Series Data Using Convolutional Neural Networks: A Case Study of Healthcare Cost Prediction

abstract: Objective: To develop an effective and scalable individual-level patient cost prediction method by automatically learning hidden temporal patterns from multivariate time series data in patient insurance claims using a convolutional neural network (CNN) architecture. Methods: We used three years of medical and pharmacy claims data from 2013 to 2016 from a healthcare insurer, where data from the first two years were used to build the model to predict costs in the third year. The data consisted of the multivariate time series of cost, visit and medical features that were shaped as images of patients' health status (i.e., matrices with time windows on one dimension and the medical, visit and cost features on the other dimension). Patients' multivariate time series images were given to a CNN method with a proposed architecture. After hyper-parameter tuning, the proposed architecture consisted of three building blocks of convolution and pooling layers with an LReLU activation function and a customized kernel size at each layer for healthcare data. The proposed CNN learned temporal patterns became inputs to a fully connected layer. Conclusions: Feature learning through the proposed CNN configuration significantly improved individual-level healthcare cost prediction. The proposed CNN was able to outperform temporal pattern detection methods that look for a pre-defined set of pattern shapes, since it is capable of extracting a variable number of patterns with various shapes. Temporal patterns learned from medical, visit and cost data made significant contributions to the prediction performance. Hyper-parameter tuning showed that considering three-month data patterns has the highest prediction accuracy. Our results showed that patients' images extracted from multivariate time series data are different from regular images, and hence require unique designs of CNN architectures.


title: Learning Mixtures of Permutations: Groups of Pairwise Comparisons and Combinatorial Method of Moments

abstract: In applications such as rank aggregation, mixture models for permutations are frequently used when the population exhibits heterogeneity. In this work, we study the widely used Mallows mixture model. In the high-dimensional setting, we propose a polynomial-time algorithm that learns a Mallows mixture of permutations on $n$ elements with the optimal sample complexity that is proportional to $\log n$, improving upon previous results that scale polynomially with $n$. In the high-noise regime, we characterize the optimal dependency of the sample complexity on the noise parameter. Both objectives are accomplished by first studying demixing permutations under a noiseless query model using groups of pairwise comparisons, which can be viewed as moments of the mixing distribution, and then extending these results to the noisy Mallows model by simulating the noiseless oracle.


title: Existence and nonexistence of positive solutions for singular (p,q)-equations with superdiffusive perturbation

abstract: We consider a nonlinear Dirichlet problem driven by the $(p,q)$-Laplacian and with a reaction which is parametric and exhibits the combined effects of a singular term and of a superdiffusive one. We prove an existence and nonexistence result for positive solutions depending on the value of the parameter $\lambda \in \overset{\circ}{\mathbb{R}}_+=(0,+\infty)$.


title: A Finite, Feasible, Quantifier-free Foundation for Constructive Geometry

abstract: In this paper we will develop an axiomatic foundation for the geometric study of straight edge, protractor, and compass constructions, which while being related to previous foundations, will be the first to have all axioms written and all proofs conducted in quantifier-free first order logic. All constructions within the system will be justified to be feasible by basic human faculties. No statement in the system will refer to infinitely many objects and one can posit an interpretation of the system which is in accordance to our free, creative process of geometric constructions. We are also able to capture analogous results to Euclid's work on non-planar geometry in Book XI of The Elements. This paper primarily builds on Suppes' paper Quantifier-Free Axioms for Constructive Affine Plane Geometry and draws from Beeson's article A Constructive Version of Tarski's Geometry. By further developing Suppes' work on parallel line segments, we are able to develop analogs to most theorems about parallel lines without assuming an equivalent to Euclid's Fifth Postulate which we deem as introducing non-feasible constructions. In A Constructive Version of Tarski's Geometry, Beeson defines the characteristics such a geometric foundation should have to called constructive. This work satisfies these characteristics. Additionally this work would be considered constructive as Suppes defined it.


title: VRFT with ARX controller model and constrained total least squares

abstract: The virtual reference feedback tuning (VRFT) is a non-iterative data-driven (DD) method employed to tune a controller's parameters aiming to achieve a prescribed closed-loop performance. In its most common formulation, the parameters of a linearly parametrized controller are estimated by solving a least squares (LS) problem, which in the presence of noise leads to a biased estimate of the controller's parameters. To eliminate this bias, an instrumental variable (IV) variant of the method is usual, at the cost of increasing significantly the estimate's variance. In the present work, we propose to apply the constrained total least squares (CTLS) solution to the VRFT problem. We formulate explicitly the VRFT solution with CTLS for controllers described by an autoregressive exogenous (ARX) model. The effectiveness of the proposed solution is illustrated by two case studies in which it is compared to the usual VRFT solutions and to another, statistically efficient, design method.


title: The asymptotic geometry of $\rm{G}_2$-monopoles

abstract: This article investigates the asymptotics of $\rm{G}_2$-monopoles. First, we prove that when the underlying $\rm{G}_2$-manifold is nonparabolic (i.e. admits a positive Green's function), finite intermediate energy monopoles with bounded curvature have finite mass. The second main result restricts to the case when the underlying $\rm{G}_2$-manifold is asymptotically conical. In this situation, we deduce sharp decay estimates and that the connection converges, along the end, to a pseudo-Hermitian--Yang--Mills connection over the asymptotic cone. Finally, our last result exhibits a Fredholm setup describing the moduli space of finite intermediate energy monopoles on an asymptotically conical $\rm{G}_2$-manifold.


title: The Cost of Software-Based Memory Management Without Virtual Memory

abstract: Virtual memory has been a standard hardware feature for more than three decades. At the price of increased hardware complexity, it has simplified software and promised strong isolation among colocated processes. In modern computing systems, however, the costs of virtual memory have increased significantly. With large memory workloads, virtualized environments, data center computing, and chips with multiple DMA devices, virtual memory can degrade performance and increase power usage. We therefore explore the implications of building applications and operating systems without relying on hardware support for address translation. Primarily, we investigate the implications of removing the abstraction of large contiguous memory segments. Our experiments show that the overhead to remove this reliance is surprisingly small for real programs. We expect this small overhead to be worth the benefit of reducing the complexity and energy usage of address translation. In fact, in some cases, performance can even improve when address translation is avoided.


title: First-Order Methods for Wasserstein Distributionally Robust MDP

abstract: Markov decision processes (MDPs) are known to be sensitive to parameter specification. Distributionally robust MDPs alleviate this issue by allowing for \emph{ambiguity sets} which give a set of possible distributions over parameter sets. The goal is to find an optimal policy with respect to the worst-case parameter distribution. We propose a framework for solving Distributionally robust MDPs via first-order methods, and instantiate it for several types of Wasserstein ambiguity sets. By developing efficient proximal updates, our algorithms achieve a convergence rate of $O\left(NA^{2.5}S^{3.5}\log(S)\log(\epsilon^{-1})\epsilon^{-1.5} \right)$ for the number of kernels $N$ in the support of the nominal distribution, states $S$, and actions $A$; this rate varies slightly based on the Wasserstein setup. Our dependence on $N,A$ and $S$ is significantly better than existing methods, which have a complexity of $O\left(N^{3.5}A^{3.5}S^{4.5}\log^{2}(\epsilon^{-1}) \right)$. Numerical experiments show that our algorithm is significantly more scalable than state-of-the-art approaches across several domains.


title: The simple type conjecture for mod 2 Seiberg-Witten invariants

abstract: We prove that, under a simple condition on the cohomology ring, every closed 4-manifold has mod 2 Seiberg-Witten simple type. This result shows that there exists a large class of topological 4-manifolds such that all smooth structures have mod 2 simple type, and yet some have non-vanishing (mod 2) Seiberg-Witten invariants. As corollaries, we obtain adjunction inequalities and show that, under a mild topological condition, every geometrically simply connected closed 4-manifold has the vanishing mod 2 Seiberg-Witten invariant for at least one orientation.


title: On sets of subspaces with two intersection dimensions and a geometrical junta bound

abstract: In this article, constant dimension subspace codes whose codewords have subspace distance in a prescribed set of integers, are considered. The easiest example of such an object is a {\it junta}; i.e. a subspace code in which all codewords go through a common subspace. We focus on the case when only two intersection values for the codewords, are assigned. In such a case we determine an upper bound for the dimension of the vector space spanned by the elements of a non-junta code. In addition, if the two intersection values are consecutive, we prove that such a bound is tight, and classify the examples attaining the largest possible dimension as one of four infinite families.


title: Counting ternary trees according to the number of middle edges and factorizing into $(3/2)$-ary trees

abstract: The sequence A120986 in the Encyclopedia of Integer Sequences counts ternary trees according to the number of nodes and the number of middle edges. Using a certain substition, the underlying cubic equation can be factored. This leads to an extension of the concept of $(3/2)$-ary trees, introduced by Knuth in his christmas lecture from 2014.


title: Coarse quotients of metric spaces and embeddings of uniform Roe algebras

abstract: We study embeddings of uniform Roe algebras which have "large range" in their codomain and the relation of those with coarse quotients between metric spaces. Among other results, we show that if $Y$ has property A and there is an embedding $\Phi:\mathrm{C}^*_u(X)\to \mathrm{C}^*_u(Y)$ with "large range" and so that $\Phi(\ell_\infty(X))$ is a Cartan subalgebra of $\mathrm{C}^*_u(Y)$, then there is a bijective coarse quotient $X\to Y$. This shows that the large scale geometry of $Y$ is, in some sense, controlled by the one of $X$. For instance, if $X$ has finite asymptotic dimension, so does $Y$.


title: DynamicVAE: Decoupling Reconstruction Error and Disentangled Representation Learning

abstract: This paper challenges the common assumption that the weight $\beta$, in $\beta$-VAE, should be larger than $1$ in order to effectively disentangle latent factors. We demonstrate that $\beta$-VAE, with $\beta < 1$, can not only attain good disentanglement but also significantly improve reconstruction accuracy via dynamic control. The paper removes the inherent trade-off between reconstruction accuracy and disentanglement for $\beta$-VAE. Existing methods, such as $\beta$-VAE and FactorVAE, assign a large weight to the KL-divergence term in the objective function, leading to high reconstruction errors for the sake of better disentanglement. To mitigate this problem, a ControlVAE has recently been developed that dynamically tunes the KL-divergence weight in an attempt to control the trade-off to more a favorable point. However, ControlVAE fails to eliminate the conflict between the need for a large $\beta$ (for disentanglement) and the need for a small $\beta$. Instead, we propose DynamicVAE that maintains a different $\beta$ at different stages of training, thereby decoupling disentanglement and reconstruction accuracy. In order to evolve the weight, $\beta$, along a trajectory that enables such decoupling, DynamicVAE leverages a modified incremental PI (proportional-integral) controller, and employs a moving average as well as a hybrid annealing method to evolve the value of KL-divergence smoothly in a tightly controlled fashion. We theoretically prove the stability of the proposed approach. Evaluation results on three benchmark datasets demonstrate that DynamicVAE significantly improves the reconstruction accuracy while achieving disentanglement comparable to the best of existing methods. The results verify that our method can separate disentangled representation learning and reconstruction, removing the inherent tension between the two.


title: Decoding Polar Codes with Reinforcement Learning

abstract: In this paper we address the problem of selecting factor-graph permutations of polar codes under belief propagation (BP) decoding to significantly improve the error-correction performance of the code. In particular, we formalize the factor-graph permutation selection as the multi-armed bandit problem in reinforcement learning and propose a decoder that acts like an online-learning agent that learns to select the good factor-graph permutations during the course of decoding. We use state-of-the-art algorithms for the multi-armed bandit problem and show that for a 5G polar codes of length 128 with 64 information bits, the proposed decoder has an error-correction performance gain of around 0.125 dB at the target frame error rate of 10^{-4}, when compared to the approach that randomly selects the factor-graph permutations.


title: Competing AI: How does competition feedback affect machine learning?

abstract: This papers studies how competition affects machine learning (ML) predictors. As ML becomes more ubiquitous, it is often deployed by companies to compete over customers. For example, digital platforms like Yelp use ML to predict user preference and make recommendations. A service that is more often queried by users, perhaps because it more accurately anticipates user preferences, is also more likely to obtain additional user data (e.g. in the form of a Yelp review). Thus, competing predictors cause feedback loops whereby a predictor's performance impacts what training data it receives and biases its predictions over time. We introduce a flexible model of competing ML predictors that enables both rapid experimentation and theoretical tractability. We show with empirical and mathematical analysis that competition causes predictors to specialize for specific sub-populations at the cost of worse performance over the general population. We further analyze the impact of predictor specialization on the overall prediction quality experienced by users. We show that having too few or too many competing predictors in a market can hurt the overall prediction quality. Our theory is complemented by experiments on several real datasets using popular learning algorithms, such as neural networks and nearest neighbor methods.


title: On the use of local structural properties for improving the efficiency of hierarchical community detection methods

abstract: Community detection is a fundamental problem in the analysis of complex networks. It is the analogue of clustering in network data mining. Within community detection methods, hierarchical algorithms are popular. However, their iterative nature and the need to recompute the structural properties used to split the network (i.e. edge betweenness in Girvan and Newman's algorithm), make them unsuitable for large network data sets. In this paper, we study how local structural network properties can be used as proxies to improve the efficiency of hierarchical community detection while, at the same time, achieving competitive results in terms of modularity. In particular, we study the potential use of the structural properties commonly used to perform local link prediction, a supervised learning problem where community structure is relevant, as nodes are prone to establish new links with other nodes within their communities. In addition, we check the performance impact of network pruning heuristics as an ancillary tactic to make hierarchical community detection more efficient


title: The Importance of Pessimism in Fixed-Dataset Policy Optimization

abstract: We study worst-case guarantees on the expected return of fixed-dataset policy optimization algorithms. Our core contribution is a unified conceptual and mathematical framework for the study of algorithms in this regime. This analysis reveals that for naive approaches, the possibility of erroneous value overestimation leads to a difficult-to-satisfy requirement: in order to guarantee that we select a policy which is near-optimal, we may need the dataset to be informative of the value of every policy. To avoid this, algorithms can follow the pessimism principle, which states that we should choose the policy which acts optimally in the worst possible world. We show why pessimistic algorithms can achieve good performance even when the dataset is not informative of every policy, and derive families of algorithms which follow this principle. These theoretical findings are validated by experiments on a tabular gridworld, and deep learning experiments on four MinAtar environments.


title: On a conjecture of Soundararajan

abstract: Building on recent work of A. Harper (2012), and using various results of M. C. Chang (2014) and H. Iwaniec (1974) on the zero-free regions of $L$-functions $L(s,\chi)$ for characters $\chi$ with a smooth modulus $q$, we establish a conjecture of K. Soundararajan (2008) on the distribution of smooth numbers over reduced residue classes for such moduli $q$. A crucial ingredient in our argument is that, for such $q$, there is at most one "problem character" for which $L(s,\chi)$ has a smaller zero-free region. Similarly, using the "Deuring-Heilbronn" phenomenon on the repelling nature of zeros of $L$-functions close to one, we also show that Soundararajan's conjecture holds for a family of moduli having Siegel zeros.


title: Approximate Piecewise Constant Curvature Equivalent Model and Their Application to Continuum Robot Configuration Estimation

abstract: The continuum robot has attracted more attention for its flexibility. Continuum robot kinematics models are the basis for further perception, planning, and control. The design and research of continuum robots are usually based on the assumption of piecewise constant curvature (PCC). However, due to the influence of friction, etc., the actual motion of the continuum robot is approximate piecewise constant curvature (APCC). To address this, we present a kinematic equivalent model for continuum robots, i.e. APCC 2L-5R. Using classical rigid linkages to replace the original model in kinematic, the APCC 2L-5R model effectively reduces complexity and improves numerical stability. Furthermore, based on the model, the configuration self-estimation of the continuum robot is realized by monocular cameras installed at the end of each approximate constant curvature segment. The potential of APCC 2L-5R in perception, planning, and control of continuum robots remains to be explored.


title: RBE with Non-Poisson Distribution of Radiation Induced Strand Breaks

abstract: Postulating that increasing linear energy transfer (LET) causes non-random clustering of lethal lesions to deviate from the Poisson distribution, we employ a non-Poisson approach as a more flexible alternative that accounts for overdispersion of lethal lesions. Using non-homologous end-joining (NHEJ) pathway of double-strand break repair, a customized negative binomial (NB) distribution is used to describe the distribution of lethal events in a cell nucleus. The proposed model provides a novel, mechanistically based explanation for the measured values of the biological relevant quantities, such as model parameters and relative biological effectiveness (RBE) of the surviving cells, for various light ion types and LET values. The estimated quantities are compared with the predictions of several mechanism-inspired models and experimental data at medium and high LET values. The results examined are closer to the Microdosimetric-Kinetic model predictions for helium and carbon ions but progressively lower than trends predicted by the Local Effect Model and the Repair Misrepair-Fixation model in the large LET region. The results support the view that the limitation in the increase in RBE at high LET can be accounted for entirely, or, in large part, by clustering of lethal events to cause deviation from the Poisson distribution.


title: Quantum algorithm for a chemical reaction path optimization by using a variational quantum algorithm and a reaction path generation

abstract: The search for new computational tasks of quantum chemistry that can be performed on current quantum computers is important for the development of quantum computing and quantum chemistry. Although calculations of chemical reactions have a wide range of applications in quantum chemical calculations, a quantum algorithm for obtaining activation energy $E_a$, which determines the rate of chemical reactions, has not been performed. In this study, we propose a quantum algorithm for the chemical reaction path optimization to obtain $E_a$. In our algorithm, quantum circuits can be used not only for the energy evaluation by the variational quantum eigensolver (VQE) but also for chemical reaction path generation. The chemical reaction path is obtained by encoding the initial reaction path to the circuit, operating parameterized gates, and extracting the path information by measurement. The nudged elastic band method was used for optimizing a reaction path, and the ground-state calculation for each state on the path was performed by the VQE or the exact diagonalization (ED). The proposed algorithm was applied to $\mathrm{H}_2 + \mathrm{H} \rightarrow \mathrm{H}_2 + \mathrm{H}$ reaction, and we confirmed that $E_a$ was obtained accurately in the case of both the VQE and the ED. We also obtained numerical results that the entanglement between the images accelerates the path optimization. From these results, we show the feasibility of performing fast and accurate chemical reaction calculations by using quantum algorithms.


title: Quantum droplets in two-dimensional optical lattices

abstract: We study the stability of zero-vorticity and vortex lattice quantum droplets (LQDs), which are described by a two-dimensional (2D) Gross-Pitaevskii (GP) equation with a periodic potential and Lee-Huang-Yang (LHY) term. The LQDs are divided in two types: onsite-centered and offsitecentered LQDs, the centers of which are located at the minimum and the maximum of the potential, respectively. The stability areas of these two types of LQDs with different number of sites for zerovorticity and vorticity with S = 1 are given. We found that the u-N relationship of the stable LQDs with a fixed number of sites can violate the Vakhitov-Kolokolov (VK) criterion, which is a necessary stability condition for nonlinear modes with an attractive interaction. Moreover, the u-N relationship shows that two types of vortex LQDs with the same number of sites are degenerated, while the zero-vorticity LQDs are not degenerated. It is worth mentioning that the offsite-centered LQDs with zero-vorticity and vortex LQDs with S = 1 are heterogeneous.


title: Dual topological superconducting states in the layered titanium-based oxypnictide superconductor BaTi$_2$Sb$_2$O

abstract: Topological superconductors have long been predicted to host Majorana zero modes which obey non-Abelian statistics and have potential for realizing non-decoherence topological quantum computation. However, material realization of topological superconductors is still a challenge in condensed matter physics. Utilizing high-resolution angle-resolved photoemission spectroscopy and first-principles calculations, we predict and then unveil the coexistence of topological Dirac semimetal and topological insulator states in the vicinity of Fermi energy ($E_F$) in the titanium-based oxypnictide superconductor BaTi$_2$Sb$_2$O. Further spin-resolved measurements confirm its spin-helical surface states around $E_F$, which are topologically protected and give an opportunity for realization of Majorana zero modes and Majorana flat bands in one material. Hosting dual topological superconducting states, the intrinsic superconductor BaTi$_2$Sb$_2$O is expected to be a promising platform for further investigation of topological superconductivity.


title: Incentive-compatible mechanisms for online resource allocation in mobility-as-a-service systems

abstract: In the context of `Everything-as-a-Service', the transportation sector has been evolving towards user-centric business models in which customized services and mode-agnostic mobility resources are priced in a unified framework. Yet, in the vast majority of studies on Mobility as a Service (MaaS) systems, mobility resource pricing is based on segmented travel modes, e.g. private vehicle, public transit and shared mobility services. This study attempts to address this research gap by introducing innovative auction-based online MaaS mechanisms where users can bid for any amount of mode-agnostic mobility resources based on their willingness to pay and preferences. We take the perspective of a MaaS regulator which aims to maximize social welfare by allocating mobility resources to users. We propose two mechanisms which allow users to either pay for the immediate use of mobility service (pay-as-you-go), or to subscribe to mobility service packages (pay-as-a-package). We cast the proposed auction-based mechanisms as online resource allocation problems where users compete for MaaS resources and bid for travel time per trip. We propose (integer-) linear programming formulations to accommodate user bids based on available mobility resources in an online optimization approach. We show that the proposed MaaS mechanisms are incentive-compatible, develop customized online algorithms and derive performance bounds based on competitive analysis. Extensive numerical simulations are conducted on large scale instances generated from realistic mobility data, which highlight the benefits of the proposed MaaS mechanisms and the effectiveness of the proposed online optimization approaches.


title: The Radicalization Risks of GPT-3 and Advanced Neural Language Models

abstract: In this paper, we expand on our previous research of the potential for abuse of generative language models by assessing GPT-3. Experimenting with prompts representative of different types of extremist narrative, structures of social interaction, and radical ideologies, we find that GPT-3 demonstrates significant improvement over its predecessor, GPT-2, in generating extremist texts. We also show GPT-3's strength in generating text that accurately emulates interactive, informational, and influential content that could be utilized for radicalizing individuals into violent far-right extremist ideologies and behaviors. While OpenAI's preventative measures are strong, the possibility of unregulated copycat technology represents significant risk for large-scale online radicalization and recruitment; thus, in the absence of safeguards, successful and efficient weaponization that requires little experimentation is likely. AI stakeholders, the policymaking community, and governments should begin investing as soon as possible in building social norms, public policy, and educational initiatives to preempt an influx of machine-generated disinformation and propaganda. Mitigation will require effective policy and partnerships across industry, government, and civil society.


title: Optimality of short-term synaptic plasticity in modelling certain dynamic environments

abstract: Biological neurons and their in-silico emulations for neuromorphic artificial intelligence (AI) use extraordinarily energy-efficient mechanisms, such as spike-based communication and local synaptic plasticity. It remains unclear whether these neuronal mechanisms only offer efficiency or also underlie the superiority of biological intelligence. Here, we prove rigorously that, indeed, the Bayes-optimal prediction and inference of randomly but continuously transforming environments, a common natural setting, relies on short-term spike-timing-dependent plasticity, a hallmark of biological synapses. Further, this dynamic Bayesian inference through plasticity enables circuits of the cerebral cortex in simulations to recognize previously unseen, highly distorted dynamic stimuli. Strikingly, this also introduces a biologically-modelled AI, the first to overcome multiple limitations of deep learning and outperform artificial neural networks in a visual task. The cortical-like network is spiking and event-based, trained only with unsupervised and local plasticity, on a small, narrow, and static training dataset, but achieves recognition of unseen, transformed, and dynamic data better than deep neural networks with continuous activations, trained with supervised backpropagation on the transforming data. These results link short-term plasticity to high-level cortical function, suggest optimality of natural intelligence for natural environments, and repurpose neuromorphic AI from mere efficiency to computational supremacy altogether.


title: When is the rate function of a random vector strictly convex?

abstract: We give a necessary and sufficient condition for strict convexity of the rate function of a random vector in $R^d$. This condition is always satisfied when the random vector has finite Laplace transform. We also completely describe the effective domain of the rate function under a weaker condition.