title: Microscopic Models for Fusion Categories

abstract: This is a PhD Thesis on the connection between subfactors (more precisely, their corresponding fusion categories) and Conformal Field Theory (CFT). Besides being a mathematically interesting topic on its own, subfactors have also attracted the attention of physicists, since there is a conjectured correspondence between these and CFTs. Although there is quite a persuasive body of evidence for this conjecture, there are some gaps: there exists a set of exceptional subfactors with no known counterpart CFT. Hence, it is necessary to develop new techniques for building a CFT from a subfactor. Here, it is useful to study the underlying mathematical structure in more detail: The even parts of every subfactor give rise to two Unitary Fusion Categories (UFCs), and it is a promising direction to study quantum spin systems constructed from these categories to find a connection to CFTs. The simplest example that requires new techniques for building a CFT is the Haagerup subfactor, since it is the smallest subfactor with index larger than 4. In this thesis, we investigate the question whether there is a CFT corresponding to the Haagerup subfactor via lattice models in one and two dimensions. The first task here is to find the F-symbols of the fusion category since these are crucial ingredients for the construction of a physical model in all of the models we consider in this thesis. We then investigate microscopic models such as the golden chain model and the Levin-Wen model in order to find evidence for a corresponding CFT. We find that there is no evidence for a corresponding CFT from the investigation of the UFCs directly and it is necessary to expand these studies to the corresponding unitary modular tensor category, which can, for instance, be obtained via the excitations of the Levin-Wen model.

categories: math-ph cond-mat.str-el math.CT math.MP quant-ph


title: Quantum-classical nonadiabatic dynamics of Floquet driven systems

abstract: We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schroedinger equation with inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.


title: An approximate analytic solution to the coupled problems of coronal heating and solar-wind acceleration

abstract: Between the base of the solar corona and the Alfven critical point, the solar-wind density decreases by approximately five orders of magnitude, but the temperature varies by a factor of only a few. In this paper, I show that such quasi-isothermal evolution out to the Alfven critical point is a generic property of outflows powered by reflection-driven Alfven-wave (AW) turbulence, in which outward-propagating AWs partially reflect, and counter-propagating AWs interact to produce a cascade of fluctuation energy to small scales, which leads to turbulent heating. Approximating the sub-Alfvenic region as isothermal, I first present a simplified calculation of the mass outflow rate, asymptotic wind speed, and coronal temperature that neglects conductive losses and the wave pressure force. I then develop a more detailed model of the transition region, corona, and solar wind that accounts for the heat flux from the coronal base into the transition region and momentum deposition by AWs. I solve analytically for this heat flux by balancing, within the transition region, conductive heating against internal-energy losses from radiation, pdV work, and advection. The density at the coronal base is determined by locally balancing turbulent heating and radiative cooling. I solve the equations of the model analytically in two different parameter regimes. Analytic and numerical solutions to the model equations agree with a number of observations.


title: Buckling bars in nearly face-on galaxies observed with MaNGA

abstract: Over half of disk galaxies are barred, yet the mechanisms for bar formation and the life-time of bar buckling remain poorly understood. In simulations, a thin bar undergoes a rapid (<1 Gyr) event called "buckling," during which the inner part of the bar is asymmetrically bent out of the galaxy plane and eventually thickens, developing a peanut/X-shaped profile when viewed side-on. Through analyzing stellar kinematics of N-body model snapshots of a galaxy before, during, and after the buckling phase, we confirm a distinct quadrupolar pattern of out-of-plane stellar velocities in nearly face-on galaxies. This kinematic signature of buckling allows us to identify five candidates of currently buckling bars among 434 barred galaxies in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Survey, an integral field unit (IFU) spectroscopic survey that measures the composition and kinematic structure of nearby galaxies. The frequency of buckling events detected is consistent with the 0.5-1 Gyr timescale predicted by simulations. The five candidates we present more than double the total number of candidate buckling bars, and are the only ones found using the kinematic signature.


title: BERT-GT: Cross-sentence n-ary relation extraction with BERT and Graph Transformer

abstract: A biomedical relation statement is commonly expressed in multiple sentences and consists of many concepts, including gene, disease, chemical, and mutation. To automatically extract information from biomedical literature, existing biomedical text-mining approaches typically formulate the problem as a cross-sentence n-ary relation-extraction task that detects relations among n entities across multiple sentences, and use either a graph neural network (GNN) with long short-term memory (LSTM) or an attention mechanism. Recently, Transformer has been shown to outperform LSTM on many natural language processing (NLP) tasks. In this work, we propose a novel architecture that combines Bidirectional Encoder Representations from Transformers with Graph Transformer (BERT-GT), through integrating a neighbor-attention mechanism into the BERT architecture. Unlike the original Transformer architecture, which utilizes the whole sentence(s) to calculate the attention of the current token, the neighbor-attention mechanism in our method calculates its attention utilizing only its neighbor tokens. Thus, each token can pay attention to its neighbor information with little noise. We show that this is critically important when the text is very long, as in cross-sentence or abstract-level relation-extraction tasks. Our benchmarking results show improvements of 5.44% and 3.89% in accuracy and F1-measure over the state-of-the-art on n-ary and chemical-protein relation datasets, suggesting BERT-GT is a robust approach that is applicable to other biomedical relation extraction tasks or datasets.


title: Visibility of Kobayashi geodesics in convex domains and related properties

abstract: Let $D\subset \mathbb C^n$ be a bounded domain. A pair of distinct boundary points $\{p,q\}$ of $D$ has the visibility property provided there exist a compact subset $K_{p,q}\subset D$ and open neighborhoods $U_p$ of $p$ and $U_q$ of $q$, such that the real geodesics for the Kobayashi metric of $D$ which join points in $U_p$ and $U_q$ intersect $K_{p,q}$. Every Gromov hyperbolic convex domain enjoys the visibility property for any couple of boundary points. The Goldilocks domains introduced by Bharali and Zimmer and the log-type domains of Liu and Wang also enjoy the visibility property. In this paper we relate the growth of the Kobayashi distance near the boundary with visibility and provide new families of convex domains where that property holds. We use the same methods to provide refinements of localization results for the Kobayashi distance, and give a localized sufficient condition for visibility. We also exploit visibility to study the boundary behavior of biholomorphic maps.


title: Wave-kinetic approach to zonal-flow dynamics: recent advances

abstract: Basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of wave-kinetic theory. Recently, this framework has been re-examined from first principles, which has led to more accurate yet still tractable "improved" wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow "curvature" (the second radial derivative of the flow velocity) on dynamics and stability of drift waves and zonal flows. We overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.


title: Acoustic emissions of nearly steady and uniform granular flows: a proxy for flow dynamics and velocity fluctuations

abstract: The seismic waves emitted during granular flows are generated by different sources: high frequencies by inter-particle shocks and low frequencies by global motion and large scale deformation. To unravel these different mechanisms, an experimental study has been performed on the seismic waves emitted by dry quasi steady granular flows. The emitted seismic waves were recorded using shock accelerometers and the flow dynamics were captured with a fast camera. The mechanical characteristics of the particle shocks were analyzed, along with the duration between shocks and the correlations in the particle motion. The high-frequency seismic waves (1-50 kHz) were found to originate from particle shocks and waves trapped in the flowing layer. The low-frequency waves (20-60 Hz) were generated by the oscillations of the particles along their trajectories, i.e. from cycles of dilation/compression during the shear. The profiles of granular temperature (i.e. the square of particle velocity fluctuations) and average velocity were measured and related to the average properties of the flow as well as to the slope angle and flow thickness. These profiles were then used in a simple steady granular flow model to predict the radiated seismic energy and the energetic efficiency, i.e. the fraction of the flow potential energy converted to seismic energy. Scaling laws relating the seismic power, the shear strain rate and the inertial number were derived. In particular, the emitted seismic power is proportional to the granular temperature, which is also related to the mean flow velocity.


title: Algebraic Spivak's theorem and applications

abstract: We prove an analogue of Lowrey--Sch\"urg's algebraic Spivak's theorem when working over a base ring $A$ that is either a field or a nice enough discrete valuation ring, and after inverting the residual characteristic exponent $e$ in the coefficients. By this result algebraic bordism groups of quasi-projective derived $A$-schemes can be generated by classical cycles, leading to vanishing results for low degree $e$-inverted bordism classes, as well as to the classification of quasi-smooth projective $A$-schemes of low virtual dimension up to $e$-inverted cobordism. As another application, we prove that $e$-inverted bordism classes can be extended from an open subset, leading to the proof of homotopy invariance of $e$-inverted bordism groups for quasi-projective derived $A$-schemes.


title: On the Practicality of Differential Privacy in Federated Learning by Tuning Iteration Times

abstract: In spite that Federated Learning (FL) is well known for its privacy protection when training machine learning models among distributed clients collaboratively, recent studies have pointed out that the naive FL is susceptible to gradient leakage attacks. In the meanwhile, Differential Privacy (DP) emerges as a promising countermeasure to defend against gradient leakage attacks. However, the adoption of DP by clients in FL may significantly jeopardize the model accuracy. It is still an open problem to understand the practicality of DP from a theoretic perspective. In this paper, we make the first attempt to understand the practicality of DP in FL through tuning the number of conducted iterations. Based on the FedAvg algorithm, we formally derive the convergence rate with DP noises in FL. Then, we theoretically derive: 1) the conditions for the DP based FedAvg to converge as the number of global iterations (GI) approaches infinity; 2) the method to set the number of local iterations (LI) to minimize the negative influence of DP noises. By further substituting the Laplace and Gaussian mechanisms into the derived convergence rate respectively, we show that: 3) The DP based FedAvg with the Laplace mechanism cannot converge, but the divergence rate can be effectively prohibited by setting the number of LIs with our method; 4) The learning error of the DP based FedAvg with the Gaussian mechanism can converge to a constant number finally if we use a fixed number of LIs per GI. To verify our theoretical findings, we conduct extensive experiments using two real-world datasets. The results not only validate our analysis results, but also provide useful guidelines on how to optimize model accuracy when incorporating DP into FL


title: Dynamic Ordering Learning in Multivariate Forecasting

abstract: In many fields where the main goal is to produce sequential forecasts for decision making problems, the good understanding of the contemporaneous relations among different series is crucial for the estimation of the covariance matrix. In recent years, the modified Cholesky decomposition appeared as a popular approach to covariance matrix estimation. However, its main drawback relies on the imposition of the series ordering structure. In this work, we propose a highly flexible and fast method to deal with the problem of ordering uncertainty in a dynamic fashion with the use of Dynamic Order Probabilities. We apply the proposed method in two different forecasting contexts. The first is a dynamic portfolio allocation problem, where the investor is able to learn the contemporaneous relationships among different currencies improving final decisions and economic performance. The second is a macroeconomic application, where the econometrician can adapt sequentially to new economic environments, switching the contemporaneous relations among macroeconomic variables over time.


title: Elasto-capillary circumferential buckling of soft tubes under axial loading: existence and competition with localised beading and periodic axial modes

abstract: We provide an extension to previous analysis of the localised beading instability of soft slender tubes under surface tension and axial stretching. The primary questions pondered here are: under what loading conditions, if any, can bifurcation into circumferential buckling modes occur, and do such solutions dominate localisation and periodic axial modes? Three distinct boundary conditions are considered; in case 1 the tube's curved surfaces are traction free and under surface tension, whilst in cases 2 and 3 the inner and outer surfaces (respectively) are fixed to prevent radial displacement and surface tension. A linear bifurcation analysis is conducted to determine numerically the existence of circumferential mode solutions. In case 1 we focus on the tensile stress regime given the preference of slender compressed tubes towards Euler buckling over axial wrinkling. We show that tubes under several loading paths are highly sensitive to circumferential modes; in contrast, localised and periodic axial modes are absent, suggesting that the circumferential buckling is dominant by default. In case 2, circumferential mode solutions are associated with negative surface tension values and thus are physically implausible. Circumferential buckling solutions are shown to exist in case 3 for tensile and compressive axial loads, and we demonstrate for multiple loading scenarios their dominance over localisation and periodic axial modes within specific parameter regimes.


title: Conformal Scattering and the Goursat Problem for Dirac Fields in the Interior of Charged Spherically Symmetric Black Holes

abstract: We construct a conformal scattering theory for Dirac fields in the interior of a Reissner-Nordstr\"om-like black hole, between the black hole event horizon and the Cauchy horizon. The main result is a resolution of the characteristic Cauchy problem for the Dirac equation on the horizons by solving a system of wave equations and re-interpreting its solution as the components of the required Dirac solution.


title: First-Order Problem Solving through Neural MCTS based Reinforcement Learning

abstract: The formal semantics of an interpreted first-order logic (FOL) statement can be given in Tarskian Semantics or a basically equivalent Game Semantics. The latter maps the statement and the interpretation into a two-player semantic game. Many combinatorial problems can be described using interpreted FOL statements and can be mapped into a semantic game. Therefore, learning to play a semantic game perfectly leads to the solution of a specific instance of a combinatorial problem. We adapt the AlphaZero algorithm so that it becomes better at learning to play semantic games that have different characteristics than Go and Chess. We propose a general framework, Persephone, to map the FOL description of a combinatorial problem to a semantic game so that it can be solved through a neural MCTS based reinforcement learning algorithm. Our goal for Persephone is to make it tabula-rasa, mapping a problem stated in interpreted FOL to a solution without human intervention.


title: This Week's Finds in Mathematical Physics (1-50)

abstract: These are the first 50 issues of This Week's Finds of Mathematical Physics, from January 19, 1993 to March 12, 1995. These issues focus on quantum gravity, topological quantum field theory, knot theory, and applications of $n$-categories to these subjects. However, there are also digressions into Lie algebras, elliptic curves, linear logic and other subjects. They were typeset in 2020 by Tim Hosgood. If you see typos or other problems please report them. (I already know the cover page looks weird).


title: Nonadiabatic transitions in Landau-Zener grids: integrability and semiclassical theory

abstract: We demonstrate that the general model of a linearly time-dependent crossing of two energy bands is integrable. Namely, the Hamiltonian of this model has a quadratically time-dependent commuting operator. We apply this property to four-state Landau-Zener (LZ) models that have previously been used to describe the Landau-St\"uckelberg interferometry experiments with an electron shuttling between two semiconductor quantum dots. The integrability then leads to simple but nontrivial exact relations for the transition probabilities. In addition, the integrability leads to a semiclassical theory that provides analytical approximation for the transition probabilities in these models for all parameter values. The results predict a dynamic phase transition, and show that similarly-looking models belong to different topological classes.


title: Resolution-Based Distillation for Efficient Histology Image Classification

abstract: Developing deep learning models to analyze histology images has been computationally challenging, as the massive size of the images causes excessive strain on all parts of the computing pipeline. This paper proposes a novel deep learning-based methodology for improving the computational efficiency of histology image classification. The proposed approach is robust when used with images that have reduced input resolution and can be trained effectively with limited labeled data. Pre-trained on the original high-resolution (HR) images, our method uses knowledge distillation (KD) to transfer learned knowledge from a teacher model to a student model trained on the same images at a much lower resolution. To address the lack of large-scale labeled histology image datasets, we perform KD in a self-supervised manner. We evaluate our approach on two histology image datasets associated with celiac disease (CD) and lung adenocarcinoma (LUAD). Our results show that a combination of KD and self-supervision allows the student model to approach, and in some cases, surpass the classification accuracy of the teacher, while being much more efficient. Additionally, we observe an increase in student classification performance as the size of the unlabeled dataset increases, indicating that there is potential to scale further. For the CD data, our model outperforms the HR teacher model, while needing 4 times fewer computations. For the LUAD data, our student model results at 1.25x magnification are within 3% of the teacher model at 10x magnification, with a 64 times computational cost reduction. Moreover, our CD outcomes benefit from performance scaling with the use of more unlabeled data. For 0.625x magnification, using unlabeled data improves accuracy by 4% over the baseline. Thus, our method can improve the feasibility of deep learning solutions for digital pathology with standard computational hardware.


title: Maximizing Post-selected Quantum Correlations from Classical Interference in a Multi-core Fiber Beamsplitter

abstract: Fourth-order interference is an information processing primitive for photonic quantum technologies. When used in conjunction with post-selection, it forms the basis of photonic controlled logic gates, entangling measurements, and can be used to produce quantum correlations. Here, using classical weak coherent states as inputs, we study fourth-order interference in novel $4 \times 4$ multi-port beam splitters built within multi-core optical fibers. Using two mutually incoherent weak laser pulses as inputs, we observe high-quality fourth order interference between photons from different cores, as well as self-interference of a two-photon wavepacket. In addition, we show that quantum correlations, in the form of quantum discord, can be maximized by controlling the intensity ratio between the two input weak coherent states. This should allow for the exploitation of quantum correlations in future telecommunication networks.


title: Comet 2P/Encke in apparitions of 2013 and 2017: I. Imaging photometry and long-slit spectroscopy

abstract: We present the results of imaging photometric and long-slit spectroscopic observations of comet 2P/Encke performed at the heliocentric distance 0.56 au, geocentric distance 0.65 au, and phase angle 109.2 deg on November 4, 2013 and at 1.05 au, 1.34 au, and 46.8 deg on January 23, 2017. Observations were carried out at the 6-m BTA telescope of the Special Astrophysical Observatory (Russia) with the multimode focal reducer SCORPIO-2. In 2013, the direct images of comet Encke were obtained with the broad-band V filters, whereas in 2017 the narrow-band cometary BC, RC, and NH2 filters as well as the medium-band SED500 and broad-band r-sdss filters were used for observations. About 60 emissions belonging to the CN, C2, C3, NH2, CH, and CO+ molecules were identified within the range 3750-7100 {\AA}. The ratios of the production rates C2/CN and C3/CN correspond to the typical comets, not depleted in the carbon-chain. A complex structure of the coma was detected in both observational periods. In January 2017, the dust was in general concentrated near the nucleus, the dust/gas ratio was 2.9 in the r-sdss filter, however, this ratio was larger than 1 at distances 3000-40000 km from the nucleus. We found that about 75% of the flux of the reflected light in the central pixel was due to the nucleus, whereas the nucleus's flux contributed 48% in the total intensity of the 2000 km area of the coma. We found that after correction for the dust coma contamination the nucleus magnitude is 18.8m+/-0.2m.


title: Online rating system development using blockchain-based distributed ledger technology

abstract: In most websites, the online rating system provides the ratings of products and services to users. Lack of trust in data integrity and its manipulation has hindered fulfilling user satisfaction. Since existing online rating systems deal with a central server, all rating data is stored on the central server. Therefore, all rating data can be removed, modified, and manipulated by the system manager to change the ratings in favor of the service or product provider. In this paper, an online rating system using distributed ledger technologies has been presented as the proposed system to solve all the weaknesses of current systems. Distributed ledger technologies are completely decentralized and there is no centralization on them by any institution. Distributed ledger technologies have different variants. Among distributed ledger technologies, blockchain technology has been used in the proposed rating system because of its support for smart contracts. In the proposed online rating system, the Ethereum platform has been chosen from different blockchain platforms that have a public permission network. In this system, the raters cannot rate unless they submit a request to the system and be authorized to take part in the online product rating process. The important feature of the Ethereum platform is its support for smart contracts, which can be used to write the rating contract in the Solidity language. Also, using Proof of Authority consensus mechanisms, all rating transactions are approved by the surveyors. Since in the real Ethereum system, each rating transaction is sent to the network by the raters, some gas must be paid for each rating transaction. However, since this method is expensive, TestNet blockchain can be used in the rating system. Finally, the proposed rating system was used for rating the restaurants of a website and its features were tested.


title: A global method for deterministic and stochastic homogenisation in $BV$

abstract: In this paper we study the deterministic and stochastic homogenisation of free-discontinuity functionals under \emph{linear} growth and coercivity conditions. The main novelty of our deterministic result is that we work under very general assumptions on the integrands which, in particular, are not required to be periodic in the space variable. Combining this result with the pointwise Subadditive Ergodic Theorem by Akcoglu and Krengel, we prove a stochastic homogenisation result, in the case of stationary random integrands. In particular, we characterise the limit integrands in terms of asymptotic cell formulas, as in the classical case of periodic homogenisation.


title: The Science of Fundamental Catalogs

abstract: This review paper discusses the science of astrometric catalogs, their current applications and future prospects for making progress in fundamental astronomy, astrophysics and gravitational physics. We discuss the concept of fundamental catalogs, their practical realizations, and future prospects. Particular attention is paid to the astrophysical implementations of the catalogs such as the measurement of the Oort constants, the secular aberration and parallax, and asteroseismology. We also consider the use of the fundamental catalogs in gravitational physics for testing general theory of relativity and detection of ultra-long gravitational waves of cosmological origin.


title: Learning with Comparison Feedback: Online Estimation of Sample Statistics

abstract: We study an online version of the noisy binary search problem where feedback is generated by a non-stochastic adversary rather than perturbed by random noise. We reframe this as maintaining an accurate estimate for the median of an adversarial sequence of integers, $x_1, x_2, \dots$, in a model where each number $x_t$ can only be accessed through a single threshold query of the form ${1(x_t \leq q_t)}$. In this online comparison feedback model, we explore estimation of general sample statistics, providing robust algorithms for median, CDF, and mean estimation with nearly matching lower bounds. We conclude with several high-dimensional generalizations.


title: Simulations of Axion Minihalos

abstract: The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. We use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro-Frenk-White profile, and the minihalos' concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. Finally, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.


title: Action Priors for Large Action Spaces in Robotics

abstract: In robotics, it is often not possible to learn useful policies using pure model-free reinforcement learning without significant reward shaping or curriculum learning. As a consequence, many researchers rely on expert demonstrations to guide learning. However, acquiring expert demonstrations can be expensive. This paper proposes an alternative approach where the solutions of previously solved tasks are used to produce an action prior that can facilitate exploration in future tasks. The action prior is a probability distribution over actions that summarizes the set of policies found solving previous tasks. Our results indicate that this approach can be used to solve robotic manipulation problems that would otherwise be infeasible without expert demonstrations. Source code is available at \url{https://github.com/ondrejba/action_priors}.


title: Shocks in the Stacked Sunyaev-Zel'dovich Profiles of Clusters I: Analysis with the Three Hundred Simulations

abstract: Gas infalling into the gravitational potential wells of massive galaxy clusters is expected to experience one or more shocks on its journey to becoming part of the intracluster medium (ICM). These shocks are important for setting the thermodynamic properties of the ICM and can therefore impact cluster observables such as X-ray emission and the Sunyaev-Zel'dovich (SZ) effect. We investigate the possibility of detecting signals from cluster shocks in the averaged thermal SZ profiles of galaxy clusters. Using zoom-in hydrodynamic simulations of massive clusters from the Three Hundred Project, we show that if cluster SZ profiles are stacked as a function of $R/R_{200m}$, shock-induced features appear in the averaged SZ profile. These features are not accounted for in standard fitting formulae for the SZ profiles of galaxy clusters. We show that the shock features should be detectable with samples of clusters from ongoing and future SZ surveys. We also demonstrate that the location of these features is correlated with the cluster accretion rate, as well as the location of the cluster splashback radius. Analyses of ongoing and future surveys, such as SPT-3G, AdvACT, Simons Observatory and CMB-S4, that include gas shocks will gain a new handle on the properties and dynamics of the outskirts of massive halos, both in gas and in mass.


title: Analysis of deviation from neoclassical ion equilibrium against electron and ion temperature profiles in T-10 tokamak

abstract: The results of the analysis of the deviation of the force equilibrium for ions from the neoclassical theory prediction, calculated using the direct measurements of the radial electric field, in the view of its possible local and nonlocal correlation with the profiles of electron, Te, and ion, Ti, temperatures in the T-10 tokamak are presented. Local correlations are analyzed by means of the Pearson's correlation. Nonlocal correlations are treated with an inverse problem under the assumption of an integral equation relationship between the deviation and Te and Ti profiles. The discharges with zero, weak and strong auxiliary heating (electron cyclotron resonance heating) are analyzed. It is found that the electrons substantially (not less than ions) contribute to the deviation of the ion equilibrium from the neoclassical theory prediction both in the local and nonlocal models.


title: A note on error analysis for a nonconforming discretisation of the tri-Helmholtz equation with singular data

abstract: We apply the nonconforming discretisation of Wu and Xu (2019) to the tri-Helmholtz equation on the plane where the source term is a functional evaluating the test function on a one-dimensional mesh-aligned embedded curve. We present error analysis for the convergence of the discretisation and linear convergence as a function of mesh size is recovered almost everywhere away from the embedded curve which aligns with classic regularity theory.


title: Random projections for conic programs

abstract: We discuss the application of random projections to conic programming: notably linear, second-order and semidefinite programs. We prove general approximation results on feasibility and optimality using the framework of formally real Jordan algebras. We then discuss some computational experiments on randomly generated semidefinite programs in order to illustrate the practical applicability of our ideas


title: Detecting seasonal episodic-like spatiotemporal memory patterns using animal movement modelling

abstract: 1. Spatial memory plays a role in the way animals perceive their environments, resulting in memory-informed movement patterns that are observable to ecologists. Developing mathematical techniques to understand how animals use memory in their environments allows for an increased understanding of animal cognition. 2. Here we describe a model that accounts for the memory of seasonal or ephemeral qualities of an animal's environment. The model captures multiple behaviors at once by allowing for resource selection in the present time as well as long-distance navigations to previously visited locations within an animal's home range. 3. We performed a set of analyses on simulated data to test our model, determining that it can provide informative results from as little as one year of discrete-time location data. We also show that the accuracy of model selection and parameter estimation increases with more location data. 4. This model has potential to identify a specific mechanism in which animals use memory to optimize their foraging, by revisiting temporally and predictably variable resources at consistent time lags.


title: Asymptotics of the Number of Endpoints of a Random Walk on a Certain Class of Directed Metric Graphs

abstract: A certain class of directed metric graphs is considered. Asymptotics for a number of possible endpoints of a random walk at large times is found.


title: PEng4NN: An Accurate Performance Estimation Engine for Efficient Automated Neural Network Architecture Search

abstract: Neural network (NN) models are increasingly used in scientific simulations, AI, and other high performance computing (HPC) fields to extract knowledge from datasets. Each dataset requires tailored NN model architecture, but designing structures by hand is a time-consuming and error-prone process. Neural architecture search (NAS) automates the design of NN architectures. NAS attempts to find well-performing NN models for specialized datsets, where performance is measured by key metrics that capture the NN capabilities (e.g., accuracy of classification of samples in a dataset). Existing NAS methods are resource intensive, especially when searching for highly accurate models for larger and larger datasets. To address this problem, we propose a performance estimation strategy that reduces the resources for training NNs and increases NAS throughput without jeopardizing accuracy. We implement our strategy via an engine called PEng4NN that plugs into existing NAS methods; in doing so, PEng4NN predicts the final accuracy of NNs early in the training process, informs the NAS of NN performance, and thus enables the NAS to terminate training NNs early. We assess our engine on three diverse datasets (i.e., CIFAR-100, Fashion MNIST, and SVHN). By reducing the training epochs needed, our engine achieves substantial throughput gain; on average, our engine saves 61% to 82% of training epochs, increasing throughput by a factor of 2.5 to 5 compared to a state-of-the-art NAS method. We achieve this gain without compromising accuracy, as we demonstrate with two key outcomes. First, across all our tests, between 74% and 97% of the ground truth best models lie in our set of predicted best models. Second, the accuracy distributions of the ground truth best models and our predicted best models are comparable, with the mean accuracy values differing by at most .7 percentage points across all tests.


title: Fractional solitons in non-Euclidian elastic plates

abstract: We show that minimal-surface non-Euclidean elastic plates share the same low-energy effective theory as Haldane's dimerized quantum spin chain. As a result, such elastic plates support fractional excitations, which take the form of charge-$1/2$ solitons between degenerate states of the plates, in strong analogy to their quantum counterpart. These fractional solitons exhibit properties similar to fractional excitations in quantum fractional topological states, including deconfinement and braiding, as well as unique new features such as holographic properties and diode-like nonlinear response, demonstrating great potentials for applications as mechanical metamaterials.


title: Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

abstract: We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $\mu$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $\mu$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.


title: Downflowing umbral flashes as an evidence of standing waves in sunspot umbrae

abstract: Umbral flashes are sudden brightenings commonly visible in the core of chromospheric lines. Theoretical and numerical modeling suggest that they are produced by the propagation of shock waves. According to these models and early observations, umbral flashes are associated with upflows. However, recent studies have reported umbral flashes in downflowing atmospheres. We aim to understand the origin of downflowing umbral flashes. We explore how the existence of standing waves in the umbral chromosphere impacts the generation of flashed profiles. We performed numerical simulations of wave propagation in a sunspot umbra with the code MANCHA. The Stokes profiles of the Ca II 8542 \AA\ line were synthesized with NICOLE. For freely-propagating waves, the chromospheric temperature enhancements of the oscillations are in phase with velocity upflows. In this case, the intensity core of the Ca II 8542 \AA\ atmosphere is heated during the upflowing stage of the oscillation. If we consider a different scenario with a resonant cavity, the wave reflections at the sharp temperature gradient of the transition region lead to standing oscillations. In this situation, temperature fluctuations are shifted backward and temperature enhancements partially coincide with the downflowing stage of the oscillation. In umbral flashes produced by standing oscillations, the reversal of the emission feature is produced when the oscillation is downflowing. The chromospheric temperature keeps increasing while the atmosphere is changing from a downflow to an upflow. During the appearance of flashed Ca II 8542 \AA\ cores, the atmosphere is upflowing most of the time, and only 38\% of the flashed profiles are associated with downflows. We find a scenario that remarkably explains the recent empirical findings of downflowing umbral flashes as a natural consequence of the presence of standing oscillations above sunspot umbrae.


title: Assignment of Freight Traffic in a Large-Scale Intermodal Network under Uncertainty

abstract: This paper presents a methodology for freight traffic assignment in a large-scale road-rail intermodal network under uncertainty. Network uncertainties caused by natural disasters have dramatically increased in recent years. Several of these disasters (e.g., Hurricane Sandy, Mississippi River Flooding, and Hurricane Harvey) severely disrupted the U.S. freight transportation network, and consequently, the supply chain. To account for these network uncertainties, a stochastic freight traffic assignment model is formulated. An algorithmic framework, involving the sample average approximation and gradient projection algorithm, is proposed to solve this challenging problem. The developed methodology is tested on the U.S. intermodal network with freight flow data from the Freight Analysis Framework. The experiments consider three types of natural disasters that have different risks and impacts on transportation networks: earthquakes, hurricanes, and floods. It is found that for all disaster scenarios, freight ton-miles are higher compared to the base case without uncertainty. The increase in freight ton-miles is the highest under the flooding scenario; this is because there are more states in the flood-risk areas, and they are scattered throughout the U.S.


title: Distant Jupiter family Comet P/2011 P1 (McNaught)

abstract: The spectra and images obtained through broadband BVR filters for Jupiter famaly comet P/2011 P1 were analyzed. We observed the comet on November 24, 2011, when its heliocentric distance was 5.43 AU. Two dimensional long slit spectra and photometric images were obtained using the focal reducer SCORPIO attached to the prime focus of the 6-m telescope BTA (SAO RAS, Russia). The spectra cover the wavelength range of 4200-7000 {\AA}. No emissions of C2 and CO+, which are expected in this wavelength region, were detected above 3 sigma level. An upper limit in gas production rate of C2 is expected to be 1.1*10^24 mol*s^-1. The continuum shows a reddening effect with the normalized gradient of reflectivity along dispersion of 5.1% per 1000 {\AA}. The color indices (B-V)=0.89 and (V-R) = 0.42 for the nucleus region or (B-V)=0.68 and (V-R)=0.39 for the coma region, which are derived from the photometric data, also evidence that the color of the cometary nucleus and dust are redder with respect to the Sun. The dust coma like a spiral galaxy adge on was fitted using a Monte Carlo model.


title: Field-angle dependent vortex lattice phase diagram in MgB2

abstract: Using small-angle neutron scattering we have studied the superconducting vortex lattice (VL) phase diagram in MgB2 as the applied magnetic field is rotated away from the c axis and towards the basal plane. The field rotation gradually suppresses the intermediate VL phase which exists between end states aligned with two high symmetry directions in the hexagonal basal plane for H || c. Above a critical angle, the intermediate state disappears, and the previously continuous transition becomes discontinuous. The evolution towards the discontinuous transition can be parameterized by a vanishing twelvefold anisotropy term in the VL free energy.


title: Quantum Consensus: an overview

abstract: We review the literature about reaching agreement in quantum networks, also called quantum consensus. After a brief introduction to the key feature of quantum computing, allowing the reader with no quantum theory background to have minimal tools to understand, we report a formal definition of quantum consensus and the protocols proposed. Proposals are classified according to the quantum feature used to achieve agreement.


title: Comet 2P/Encke in apparition of 2017: II. Polarization and color

abstract: We present results of imaging polarimetry of comet 2P/Encke performed on January 23, 2017 at the heliocentric (1.052 au) and geocentric (1.336 au) distances and phase angle 46.8 deg, 46 days before perihelion. Observations were made through the medium-band SED500 ({\lambda}5019/246 {\AA}) and broadband r-sdss ({\lambda}6200/1200 {\AA}) filters with the multimode focal reducer SCORPIO-2 at the 6-m BTA telescope of the Special Astrophysical Observatory (Russia). Dust in comet 2P/Encke was mainly concentrated in the near-nucleus region of the coma: the maximum dust/gas ratios were 1.5 and 2.9 in the SED500 and the r-sdss filters near the nucleus but dropped sharply to ~0.2 and ~1 at the distance ~2500 km, respectively. Then these ratios began to increase at distances ~12000 km from the nucleus, the ratio was ~0.3 (SED500) and ~1.3 (r-sdds). There were significant variations of polarization over the coma, which correlated with the variations in the dust color and dust/gas ratio. Changes in polarization and color across the 2P/Encke coma indicate changes in physical properties of the dust particles with the distance from the nucleus. Our Sh-matrix computer simulations of light scattering by Gaussian particles allow us to suggest that the observed trends in color and polarization are mainly result from changing particle size.


title: Protecting Big Data Privacy Using Randomized Tensor Network Decomposition and Dispersed Tensor Computation

abstract: Data privacy is an important issue for organizations and enterprises to securely outsource data storage, sharing, and computation on clouds / fogs. However, data encryption is complicated in terms of the key management and distribution; existing secure computation techniques are expensive in terms of computational / communication cost and therefore do not scale to big data computation. Tensor network decomposition and distributed tensor computation have been widely used in signal processing and machine learning for dimensionality reduction and large-scale optimization. However, the potential of distributed tensor networks for big data privacy preservation have not been considered before, this motivates the current study. Our primary intuition is that tensor network representations are mathematically non-unique, unlinkable, and uninterpretable; tensor network representations naturally support a range of multilinear operations for compressed and distributed / dispersed computation. Therefore, we propose randomized algorithms to decompose big data into randomized tensor network representations and analyze the privacy leakage for 1D to 3D data tensors. The randomness mainly comes from the complex structural information commonly found in big data; randomization is based on controlled perturbation applied to the tensor blocks prior to decomposition. The distributed tensor representations are dispersed on multiple clouds / fogs or servers / devices with metadata privacy, this provides both distributed trust and management to seamlessly secure big data storage, communication, sharing, and computation. Experiments show that the proposed randomization techniques are helpful for big data anonymization and efficient for big data storage and computation.


title: The genus-zero five-vertex model

abstract: We study the free energy and limit shape problem for the five-vertex model with periodic "genus zero" weights. We derive the exact phase diagram, free energy and surface tension for this model. We show that its surface tension has trivial potential and use this to give explicit parameterizations of limit shapes.


title: Decomposition of the SU(2) gauge field in the Maximal Abelian gauge

abstract: We study decomposition of $SU(2)$ gauge field into monopole and monopoleless components. After fixing the Maximal Abelian gauge in $SU(2)$ lattice gauge theory we decompose the nonabelian gauge field into the Abelian field created by monopoles and the modified nonabelian field with monopoles removed. We then calculate respective static potentialis and show that the potential due to the modified nonabelian field is nonconfining while, as is well known, the Abelian field produces linear potential. We further find that the sum of these potentials approximates the nonabelian static potential with good precision at all distances considered. We conclude that at large distances the monopole field potential describes the classical energy of the hadronic string while the static potential due to the modified nonabelian field describes the string fluctuations energy.


title: Clustering Word Embeddings with Self-Organizing Maps. Application on LaRoSeDa -- A Large Romanian Sentiment Data Set

abstract: Romanian is one of the understudied languages in computational linguistics, with few resources available for the development of natural language processing tools. In this paper, we introduce LaRoSeDa, a Large Romanian Sentiment Data Set, which is composed of 15,000 positive and negative reviews collected from one of the largest Romanian e-commerce platforms. We employ two sentiment classification methods as baselines for our new data set, one based on low-level features (character n-grams) and one based on high-level features (bag-of-word-embeddings generated by clustering word embeddings with k-means). As an additional contribution, we replace the k-means clustering algorithm with self-organizing maps (SOMs), obtaining better results because the generated clusters of word embeddings are closer to the Zipf's law distribution, which is known to govern natural language. We also demonstrate the generalization capacity of using SOMs for the clustering of word embeddings on another recently-introduced Romanian data set, for text categorization by topic.


title: Coherent scattering of an optically-modulated electron beam by atoms

abstract: Recent technological advances allowed the coherent optical manipulation of high-energy electron wavepackets with attosecond precision. Here we theoretically investigate the collision of optically-modulated pulsed electron beams with atomic targets and reveal a quantum interference associated with different momentum components of the incident broadband electron pulse, which coherently modulates both the elastic and inelastic scattering cross sections. We show that the quantum interference has a high spatial sensitivity at the level of Angstroms, offering potential applications in high-resolution ultrafast electron microscopy. Our findings are rationalized by a simple model.


title: Where you live matters: a spatial analysis of COVID-19 mortality

abstract: The COVID-19 pandemic has caused ~ 2 million fatalities. Significant progress has been made in advancing our understanding of the disease process, one of the unanswered questions, however, is the anomaly in the case/mortality ratio with Mexico as a clear example. Herein, this anomaly is explored by spatial analysis and whether mortality varies locally according to local factors. To address this, hexagonal cartogram maps (hexbin) used to spatially map COVID-19 mortality and visualise association with patient-level data on demographics and pre-existing health conditions. This was further interrogated at local Mexico City level by choropleth mapping. Our data show that the use of hexagonal cartograms is a better approach for spatial mapping of COVID-19 data in Mexico as it addresses bias in area size and population. We report sex/age-related spatial relationship with mortality amongst the Mexican states and a trend between health conditions and mortality at the state level. Within Mexico City, there is a clear south, north divide with higher mortality in the northern municipalities. Deceased patients in these northern municipalities have the highest pre-existing health conditions. Taken together, this study provides an improved presentation of COVID-19 mapping in Mexico and demonstrates spatial divergence of the mortality in Mexico.


title: Smartajweed Automatic Recognition of Arabic Quranic Recitation Rules

abstract: Tajweed is a set of rules to read the Quran in a correct Pronunciation of the letters with all its Qualities, while Reciting the Quran. which means you have to give every letter in the Quran its due of characteristics and apply it to this particular letter in this specific situation while reading, which may differ in other times. These characteristics include melodic rules, like where to stop and for how long, when to merge two letters in pronunciation or when to stretch some, or even when to put more strength on some letters over other. Most of the papers focus mainly on the main recitation rules and the pronunciation but not (Ahkam AL Tajweed) which give different rhythm and different melody to the pronunciation with every different rule of (Tajweed). Which is also considered very important and essential in Reading the Quran as it can give different meanings to the words. In this paper we discuss in detail full system for automatic recognition of Quran Recitation Rules (Tajweed) by using support vector machine and threshold scoring system


title: A Minimal Timescale for the Continuum in 4U 1608-52 and Aql X-1

abstract: Similar to black-hole X-ray binary (BHXRB) transients, hysteresis-like state transitions are also seen in some neutron-star X-ray binaries (NSXRBs). Using a method based on wavelets and lightcurves constructed from archival RXTE observations, we extract a minimal timescale over the complete range of transitions for 4U 1608-52 during the 2002 and 2007 outbursts and the 1999 and 2000 outbursts for Aql X-1. We present evidence for a strong positive correlation between this minimal timescale and a similar timescale extracted from the corresponding power spectra of these sources.


title: Role of Crown in Tree Resistance Against High Winds

abstract: Rather than using wooden sticks to simulate the breakage of trees in high winds as in most research, we employed fresh samples with branches and leaves to certify the crucial role played by the tree crown. By using the blowdown wind tunnel with a maximum wind speed of 60 m/s, we purposely reduce the number of leaves and show that the drag force will drop by as much as two thirds when half pruned. Based on real observations, we model the leaf by an open and full cone in the presence of light and strong wind, and calculate how their corresponding cross-sectional area and drag force vary with wind speed. Different power-law relations are predicted and confirmed by experiments for these properties before and after the formation of a full cone. Compared to the empirical value of 1/3 and 3/4, our simple model gave 2/5 and 2/3 for the power-law exponent of cross-sectional area at low and high winds. Discrepancy can be accounted for by including further details, such as the reorientation of open cones and the movement of branches.


title: Outbursting Young Stellar Object PGIR 20dci in the Perseus Arm

abstract: We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance $\sim 2.8$ kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed ($1.2~\mu$m) and the Near Earth Object Widefield Infrared Survey Explorer ($3.4$ and $4.6~\mu$m) photometric time-domain surveys. Recent near-infrared imaging reveals a new, extended scattered light nebula. Recent near-infrared spectroscopy confirms the similarity of PGIR 20dci to FU Ori type sources, based on strong molecular absorption in CO, H$_2$O, and OH, weak absorption in several atomic lines, and a warm wind/outflow as indicated by a P Cygni profile in the HeI 10830 A line. This is a rare case of an FU Ori star with a well-measured long term photometric rise before a sharper outburst, and the second instance of an FU Ori star with a documented two-step brightening in the mid-infrared.


title: Freezing sandpiles and Boolean threshold networks: equivalence and complexity

abstract: The NC versus P-hard classification of the prediction problem for sandpiles on the two dimensional grid with von Neumann neighborhood is a famous open problem. In this paper we make two kinds of progresses, by studying its freezing variant. First, it enables to establish strong connections with other well known prediction problems on networks of threshold Boolean functions such as majority. Second, we can highlight some necessary and sufficient elements to the dynamical complexity of sandpiles, with a surprisingly crucial role of cells with two grains.


title: Exceptional times when the KPZ fixed point violates Johansson's conjecture on maximizer uniqueness

abstract: In 2002, Johansson conjectured that the maximum of the Airy$_2$ process minus the parabola $x^2$ is almost surely achieved at a unique location. This result was proved a decade later by Corwin and Hammond; Moreno Flores, Quastel and Remenik; and Pimentel. Up to scaling, the Airy$_2$ process minus the parabola $x^2$ arises as the fixed time spatial marginal of the KPZ fixed point when started from narrow wedge initial data. We extend this maximizer uniqueness result to the fixed time spatial marginal of the KPZ fixed point when begun from any element of a very broad class of initial data. None of these results rules out the possibility that at random times, the KPZ fixed point spatial marginal violates maximizer uniqueness. To understand this possibility, we study the probability that the KPZ fixed point has, at a given time, two or more locations where its value is close to the maximum, obtaining quantitative upper and lower bounds in terms of the degree of closeness for a very broad class of initial data. We also compute a quantity akin to the joint density of the locations of two maximizers and the maximum value. As a consequence, the set of times of maximizer non-uniqueness almost surely has Hausdorff dimension at most two-thirds. Our analysis relies on the exact formula for the distribution function of the KPZ fixed point obtained by Matetski, Quastel and Remenik, the variational formula for the KPZ fixed point involving the Airy sheet constructed by Dauvergne, Ortmann and Vir\'ag, and the Brownian Gibbs property for the Airy$_2$ process minus the parabola $x^2$ demonstrated by Corwin and Hammond.


title: TrackMPNN: A Message Passing Graph Neural Architecture for Multi-Object Tracking

abstract: This study follows many classical approaches to multi-object tracking (MOT) that model the problem using dynamic graphical data structures, and adapts this formulation to make it amenable to modern neural networks. Our main contributions in this work are the creation of a framework based on dynamic undirected graphs that represent the data association problem over multiple timesteps, and a message passing graph neural network (MPNN) that operates on these graphs to produce the desired likelihood for every association therein. We also provide solutions and propositions for the computational problems that need to be addressed to create a memory-efficient, real-time, online algorithm that can reason over multiple timesteps, correct previous mistakes, update beliefs, and handle missed/false detections. To demonstrate the efficacy of our approach, we only use the 2D box location and object category ID to construct the descriptor for each object instance. Despite this, our model performs on par with state-of-the-art approaches that make use of additional sensors, as well as multiple hand-crafted and/or learned features. This illustrates that given the right problem formulation and model design, raw bounding boxes (and their kinematics) from any off-the-shelf detector are sufficient to achieve competitive tracking results on challenging MOT benchmarks.


title: On the classification of non-big Ulrich vector bundles on surfaces and threefolds

abstract: We classify Ulrich vector bundles that are not big on smooth complex surfaces and threefolds.


title: Asymptotically exact constants in natural convergence rate estimates in the Lindeberg theorem

abstract: Following (Shevtsova, 2010) we introduce detailed classification of the asymptotically exact constants in natural estimates of the rate of convergence in the Lindeberg central limit theorem, namely in Esseen's, Rozovskii's, and Wang-Ahmad's inequalities and their structural improvements obtained in our previous works. The above inequalities involve algebraic truncated third-order moments and the classical Lindeberg fraction and assume finiteness only the second-order moments of random summands. We present lower bounds for the introduced asymptotically exact constants as well as for the universal and for the most optimistic constants which turn to be not far from the upper ones.


title: PyHealth: A Python Library for Health Predictive Models

abstract: Despite the explosion of interest in healthcare AI research, the reproducibility and benchmarking of those research works are often limited due to the lack of standard benchmark datasets and diverse evaluation metrics. To address this reproducibility challenge, we develop PyHealth, an open-source Python toolbox for developing various predictive models on healthcare data. PyHealth consists of data preprocessing module, predictive modeling module, and evaluation module. The target users of PyHealth are both computer science researchers and healthcare data scientists. With PyHealth, they can conduct complex machine learning pipelines on healthcare datasets with fewer than ten lines of code. The data preprocessing module enables the transformation of complex healthcare datasets such as longitudinal electronic health records, medical images, continuous signals (e.g., electrocardiogram), and clinical notes into machine learning friendly formats. The predictive modeling module provides more than 30 machine learning models, including established ensemble trees and deep neural network-based approaches, via a unified but extendable API designed for both researchers and practitioners. The evaluation module provides various evaluation strategies (e.g., cross-validation and train-validation-test split) and predictive model metrics. With robustness and scalability in mind, best practices such as unit testing, continuous integration, code coverage, and interactive examples are introduced in the library's development. PyHealth can be installed through the Python Package Index (PyPI) or https://github.com/yzhao062/PyHealth .


title: General Hannan and Quinn Criterion for Common Time Series

abstract: This paper aims to study data driven model selection criteria for a large class of time series, which includes ARMA or AR($\infty$) processes, as well as GARCH or ARCH($\infty$), APARCH and many others processes. We tackled the challenging issue of designing adaptive criteria which enjoys the strong consistency property. When the observations are generated from one of the aforementioned models, the new criteria, select the true model almost surely asymptotically. The proposed criteria are based on the minimization of a penalized contrast akin to the Hannan and Quinn's criterion and then involved a term which is known for most classical time series models and for more complex models, this term can be data driven calibrated. Monte-Carlo experiments and an illustrative example on the CAC 40 index are performed to highlight the obtained results.


title: Ganymede's Surface Properties from Millimeter and Infrared Thermal Emission

abstract: We present thermal observations of Ganymede from the Atacama Large Millimeter Array (ALMA) in 2016-2019 at a spatial resolution of 300-900 km (0.1-0.2'' angular resolution) and frequencies of 97.5, 233, and 343.5 GHz (wavelengths of 3, 1.3, and 0.87 mm); the observations collectively covered all Ganymede longitudes. We determine the global thermophysical properties using a thermal model that considers subsurface emission and depth- and temperature-dependent thermophysical and dielectric properties, in combination with a retrieval algorithm. The data are sensitive to emission from the upper $\sim$0.5 meter of the surface, and we find a millimeter emissivity of 0.75-0.78 and (sub)surface porosities of 10-40%, corresponding to effective thermal inertias of 400-800 J m^{-2} K^{-1} s^{-1/2}. Combined with past infrared results, as well as modeling presented here of a previously-unpublished Galileo PPR nighttime infrared observation, the multi-wavelength constraints are consistent with a compaction profile whereby the porosity drops from ~85% at the surface to 10{+30/-10}% at depth over a compaction length scale of tens of cm. We present maps of temperature residuals from the best-fit global models which indicate localized variations in thermal surface properties at some (but not all) dark terrains and at impact craters, which appear 5-8 K colder than the model. Equatorial regions are warmer than predicted by the model, in particular near the centers of the leading and trailing hemispheres, while the mid-latitudes (~30-60 degrees) are generally colder than predicted; these trends are suggestive of an exogenic origin.


title: Relativistic wave--particle duality for spinors

abstract: We propose that relativistic wave--particle duality can be embodied in a relation $\langle u^i\rangle=\bar{\psi}\gamma^i\psi/\bar{\psi}\psi$, which determines the mean four-velocity of the fermion particle associated with a Dirac wave function. We use the Einstein--Cartan theory of gravity with torsion, which incorporates the spin-orbit interaction in curved spacetime. This relation is satisfied by a spinor plane wave and it is consistent with the energy-momentum tensor for particles. We suggest that spacetime guides the evolution of a wave function, that in turn guides the mean motion of the associated particle. Consequently, spacetime guides the motion of particles. The exact motion is limited by the uncertainty principle.


title: Generalized saturation problems for cliques, paths, and stars

abstract: A graph $G$ is $F$-saturated if it does not contain any copy of $F$, but the addition of any missing edge in $G$ creates at least one copy of $F$. Inspired by work of Alon and Shikhelman regarding a similar question for $F$-free graphs, Kritschgau, Methuku, Tait, and Timmons introduced the parameter of $\text{sat}_H(n,F)$ to denote the minimum number of copies of some subgraph $H$ in an $F$-saturated graph on $n$ vertices. In this paper, we address this generalized saturation problem with special focus on $\text{sat}_{K_r}(n,S_t)$ and $\text{sat}_{S_r}(n,S_t)$ This relates to recent work by Chakraborti and Loh regarding $\text{sat}_{K_r}(n,K_t)$ and by Ergemlidze, Methuku, Tait, and Timmons regarding $\text{sat}_{S_r}(n,K_t)$. We also provide some results regarding paths and arbitrary trees.


title: On the stability of boundary equilibria in Filippov systems

abstract: The leading-order approximation to a Filippov system $f$ about a generic boundary equilibrium $x^*$ is a system $F$ that is affine one side of the boundary and constant on the other side. We prove $x^*$ is exponentially stable for $f$ if and only if it is exponentially stable for $F$ when the constant component of $F$ is not tangent to the boundary. We then show exponential stability and asymptotic stability are in fact equivalent for $F$. We also show exponential stability is preserved under small perturbations to the pieces of $F$. Such results are well known for homogeneous systems. To prove the results here additional techniques are required because the two components of $F$ have different degrees of homogeneity. The primary function of the results is to reduce the problem of the stability of $x^*$ from the general Filippov system $f$ to the simpler system $F$. Yet in general this problem remains difficult. We provide a four-dimensional example of $F$ for which orbits appear to converge to $x^*$ in a chaotic fashion. By utilising the presence of both homogeneity and sliding motion the dynamics of $F$ can in this case be reduced to the combination of a one-dimensional return map and a scalar function.


title: Multimodal Engagement Analysis from Facial Videos in the Classroom

abstract: Student engagement is a key construct for learning and teaching. While most of the literature explored the student engagement analysis on computer-based settings, this paper extends that focus to classroom instruction. To best examine student visual engagement in the classroom, we conducted a study utilizing the audiovisual recordings of classes at a secondary school over one and a half month's time, acquired continuous engagement labeling per student (N=15) in repeated sessions, and explored computer vision methods to classify engagement levels from faces in the classroom. We trained deep embeddings for attentional and emotional features, training Attention-Net for head pose estimation and Affect-Net for facial expression recognition. We additionally trained different engagement classifiers, consisting of Support Vector Machines, Random Forest, Multilayer Perceptron, and Long Short-Term Memory, for both features. The best performing engagement classifiers achieved AUCs of .620 and .720 in Grades 8 and 12, respectively. We further investigated fusion strategies and found score-level fusion either improves the engagement classifiers or is on par with the best performing modality. We also investigated the effect of personalization and found that using only 60-seconds of person-specific data selected by margin uncertainty of the base classifier yielded an average AUC improvement of .084. 4.Our main aim with this work is to provide the technical means to facilitate the manual data analysis of classroom videos in research on teaching quality and in the context of teacher training.


title: Cloud-RAN functional split for an efficient fronthaul network

abstract: The evolution of telecommunication network towards cloud-native environments enables flexible centralization of the base band processing of radio signals. There is however a trade-off between the centralization benefits and the fronthaul cost for carrying the radio data between distributed antennas and data processing centers, which host the virtual RAN functions. In this paper, we present a specific split solution for an efficient fronthaul, which enables reducing the consumed bandwidth while being compliant with advanced cooperative radio technologies (interference reduction and data rate improvements). The proposed split has been implemented on the basis of Open Air Interface code and shows important gains in the required fronthaul bandwidth as well as significant latency reduction in the processing of radio frames. \\ \textbf{Publisher:} IEEE \\ \textbf{ISBN:}978-1-7281-3130-6


title: A $\beta$-Sturm Liouville problem associated with the general quantum operator

abstract: Let $\,I\subseteq\mathbb{R}\,$ be an interval and $\,\beta:\,I\rightarrow\,I\,$ a strictly increasing and continuous function with a unique fixed point $\,s_0\in I\,$ which satisfies $\,(s_0-t)(\beta(t)-t)\geq 0\,$ for all $\,t\in I$, where the equality holds only when $\,t=s_0$. The general quantum operator defined in 2015 by Hamza et al., $\,D_{\beta}[f](t):=\displaystyle\frac{f\big(\beta(t)\big)-f(t)}{\beta(t)-t}\,$ if $\,t\neq s_0\,$ and $\,D_{\beta}[f](s_0):=f^{\prime}(s_0)\,$ if $\,t=s_0,$ generalizes the Jackson $\,q$-operator $\,D_{q}\,$ and also the Hahn $\,(q,\omega)$-operator, $\,D_{q,\omega}$. Regarding a $\beta-$Sturm Liouville eigenvalue problem associated with the above operator $\,D_{\beta}\,$, we construct the $\beta-$Lagrange's identity, show that it is self-adjoint in $\,\mathscr{L}_{\beta}^2([a,b]),$ and exhibit some properties for the corresponding eigenvalues and eigenfunctions.


title: Tauberian theorems for statistical Ces\`{a}ro and statistical logarithmic summability of sequences in intuitionistic fuzzy normed spaces

abstract: We define statistical Ces\`{a}ro and statistical logarithmic summability methods of sequences in intuitionistic fuzzy normed spaces($IFNS$) and give slowly oscillating type and Hardy type Tauberian conditions under which statistical Ces\`{a}ro summability and statistical logarithmic summability imply convergence in $IFNS$. Besides, we obtain analogous results for the higher order summability methods as corollaries. Also, two theorems concerning the convergence of statically convergent sequences in $IFNS$ are proved in the paper.


title: Interpolation of Power Mappings

abstract: Let $(M_j)_{j=1}^\infty\in\mathbb{N}$ and $(r_j)_{j=1}^\infty\in\mathbb{R}^+$ be increasing sequences satisfying some mild rate of growth conditions. We prove that there is an entire function $f: \mathbb{C} \rightarrow\mathbb{C}$ whose behavior in the large annuli $\{ z\in\mathbb{C} : r_{j}\cdot\exp(\pi/M_{j})\leq|z|\leq r_{j+1}\}$ is given by a perturbed rescaling of $z\mapsto z^{M_j}$, such that the only singular values of $f$ are rescalings of $\pm r_j^{M_j}$. We describe several applications to the dynamics of entire functions.


title: Current-driven domain wall dynamics in ferrimagnetic Ni-doped Mn4N films : very large domain wall velocities and reversal of motion direction across the magnetic compensation point

abstract: Spin-transfer torque (STT) and spin-orbit torque (SOT) are spintronic phenomena allowing magnetization manipulation using electrical currents. Beyond their fundamental interest, they allow developing new classes of magnetic memories and logic devices, in particular based on domain wall (DW) motion. In this work, we report the study of STT driven DW motion in ferrimagnetic manganese nickel nitride (Mn4-xNixN) films, in which a fine adjustment of the Ni content allows setting the magnetic compensation at room temperature. The reduced magnetization, combined with the large spin polarization of conduction electrons, strongly enhances the STT so that domain wall velocities approaching 3000 m/s can be obtained for Ni compositions close to the compensation point. In addition, a reversal of the domain wall motion direction is observed when the magnetic compensation composition is crossed. This striking feature, related to the change of direction of the spin polarization with respect to that of the net magnetization, is clarified by ab initio band structure calculations.


title: On the Blow-up criterion of Navier-Stokes equation associated with the Weinstein operator

abstract: In this paper we give Navier-Stokes system associated with the Weinstein operator $(NSW)$ (see Eq.\eqref{11}), We study the existence and uniqueness of solutions to equations (NSW) in $L_{\alpha}^{p}\left(\mathbb{R}_{+}^{d+1}\right), 2 \alpha+d+2<p \leq \infty$, and we proved some properties of the maximal solution of equation. If the maximum time $T^*$is finite, we establish that the growth of $\left\| u ( {t}) \right\|_{L ^ {p}_{\alpha}} $ is at least of the order of ${\left(T^{*}-t\right)^{{-\frac{2 p}{p-2 \alpha-d-2}}}},$ fo rall $t$ in $\left[0, T^{*}\right]$, also we give some blow-up results.


title: Best-response dynamics, playing sequences, and convergence to equilibrium in random games

abstract: We analyze the performance of the best-response dynamic across all normal-form games using a random games approach. The playing sequence -- the order in which players update their actions -- is essentially irrelevant in determining whether the dynamic converges to a Nash equilibrium in certain classes of games (e.g. in potential games) but, when evaluated across all possible games, convergence to equilibrium depends on the playing sequence in an extreme way. Our main asymptotic result shows that the best-response dynamic converges to a pure Nash equilibrium in a vanishingly small fraction of all (large) games when players take turns according to a fixed cyclic order. By contrast, when the playing sequence is random, the dynamic converges to a pure Nash equilibrium if one exists in almost all (large) games.


title: Exploiting Multiple Timescales in Hierarchical Echo State Networks

abstract: Echo state networks (ESNs) are a powerful form of reservoir computing that only require training of linear output weights whilst the internal reservoir is formed of fixed randomly connected neurons. With a correctly scaled connectivity matrix, the neurons' activity exhibits the echo-state property and responds to the input dynamics with certain timescales. Tuning the timescales of the network can be necessary for treating certain tasks, and some environments require multiple timescales for an efficient representation. Here we explore the timescales in hierarchical ESNs, where the reservoir is partitioned into two smaller linked reservoirs with distinct properties. Over three different tasks (NARMA10, a reconstruction task in a volatile environment, and psMNIST), we show that by selecting the hyper-parameters of each partition such that they focus on different timescales, we achieve a significant performance improvement over a single ESN. Through a linear analysis, and under the assumption that the timescales of the first partition are much shorter than the second's (typically corresponding to optimal operating conditions), we interpret the feedforward coupling of the partitions in terms of an effective representation of the input signal, provided by the first partition to the second, whereby the instantaneous input signal is expanded into a weighted combination of its time derivatives. Furthermore, we propose a data-driven approach to optimise the hyper-parameters through a gradient descent optimisation method that is an online approximation of backpropagation through time. We demonstrate the application of the online learning rule across all the tasks considered.


title: Challenges and approaches to time-series forecasting in data center telemetry: A Survey

abstract: Time-series forecasting has been an important research domain for so many years. Its applications include ECG predictions, sales forecasting, weather conditions, even COVID-19 spread predictions. These applications have motivated many researchers to figure out an optimal forecasting approach, but the modeling approach also changes as the application domain changes. This work has focused on reviewing different forecasting approaches for telemetry data predictions collected at data centers. Forecasting of telemetry data is a critical feature of network and data center management products. However, there are multiple options of forecasting approaches that range from a simple linear statistical model to high capacity deep learning architectures. In this paper, we attempted to summarize and evaluate the performance of well known time series forecasting techniques. We hope that this evaluation provides a comprehensive summary to innovate in forecasting approaches for telemetry data.


title: Hankel determinants of linear combinations of moments of orthogonal polynomials, II

abstract: We present a formula that expresses the Hankel determinants of a linear combination of length $d+1$ of moments of orthogonal polynomials in terms of a $d\times d$ determinant of the orthogonal polynomials. This formula exists somehow hidden in the folklore of the theory of orthogonal polynomials but deserves to be better known, and be presented correctly and with full proof. We present four fundamentally different proofs, one that uses classical formulae from the theory of orthogonal polynomials, one that uses a vanishing argument and is due to Elouafi [J. Math. Anal. Appl. 431} (2015), 1253-1274] (but given in an incomplete form there), one that is inspired by random matrix theory and is due to Br\'ezin and Hikami [Comm. Math. Phys. 214 (2000), 111-135], and one that uses (Dodgson) condensation. We give two applications of the formula. In the first application, we explain how to compute such Hankel determinants in a singular case. The second application concerns the linear recurrence of such Hankel determinants for a certain class of moments that covers numerous classical combinatorial sequences, including Catalan numbers, Motzkin numbers, central binomial coefficients, central trinomial coefficients, central Delannoy numbers, Schr\"oder numbers, Riordan numbers, and Fine numbers.


title: DBTagger: Multi-Task Learning for Keyword Mapping in NLIDBs Using Bi-Directional Recurrent Neural Networks

abstract: Translating Natural Language Queries (NLQs) to Structured Query Language (SQL) in interfaces deployed in relational databases is a challenging task, which has been widely studied in database community recently. Conventional rule based systems utilize series of solutions as a pipeline to deal with each step of this task, namely stop word filtering, tokenization, stemming/lemmatization, parsing, tagging, and translation. Recent works have mostly focused on the translation step overlooking the earlier steps by using ad-hoc solutions. In the pipeline, one of the most critical and challenging problems is keyword mapping; constructing a mapping between tokens in the query and relational database elements (tables, attributes, values, etc.). We define the keyword mapping problem as a sequence tagging problem, and propose a novel deep learning based supervised approach that utilizes POS tags of NLQs. Our proposed approach, called \textit{DBTagger} (DataBase Tagger), is an end-to-end and schema independent solution, which makes it practical for various relational databases. We evaluate our approach on eight different datasets, and report new state-of-the-art accuracy results, $92.4\%$ on the average. Our results also indicate that DBTagger is faster than its counterparts up to $10000$ times and scalable for bigger databases.


title: A deep learning modeling framework to capture mixing patterns in reactive-transport systems

abstract: Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, we will present in this paper a deep-learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework -- placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training.


title: Measurement of the Flux-Averaged Inclusive Charged-Current Electron Neutrino and Antineutrino Cross Section on Argon using the NuMI Beam and the MicroBooNE Detector

abstract: We present a measurement of the combined $\nu_e$ + $\bar{\nu}_e$ flux-averaged charged-current inclusive cross section on argon using data from the MicroBooNE liquid argon time projection chamber (LArTPC) at Fermilab. Using the off-axis flux from the NuMI beam, MicroBooNE has reconstructed 214 candidate $\nu_e$ + $\bar{\nu}_e$ interactions with an estimated exposure of 2.4$\times10^{20}$ protons on target. Given the estimated purity of 38.6\%, this implies the observation of 80 $\nu_e$ + $\bar{\nu}_e$ events in argon, the largest such sample to date. The analysis includes the first demonstration of a fully automated application of a dE/dx-based particle discrimination technique of electron and photon induced showers in a LArTPC neutrino detector. We measure the $\nu_e + \bar{\nu}_e$ flux-averaged charged-current total cross section to be ${6.84\pm\!1.51~\textrm{(stat.)}\pm\!2.33~\textrm{(sys.)}\!\times\!10^{-39}~\textrm{cm}^{2}/~\textrm{nucleon}}$, for neutrino energies above 250 MeV and an average neutrino flux energy of 905 MeV when this threshold is applied. The measurement is sensitive to neutrino events where the final state electron momentum is above 48 MeV/c, includes the entire angular phase space of the electron, and is in agreement with the theoretical predictions from \texttt{GENIE} and \texttt{NuWro}. This measurement is also the first demonstration of electron neutrino reconstruction in a surface LArTPC in the presence of cosmic ray backgrounds, which will be a crucial task for surface experiments like those that comprise the Short-Baseline Neutrino (SBN) Program at Fermilab.


title: Implicit Unlikelihood Training: Improving Neural Text Generation with Reinforcement Learning

abstract: Likelihood training and maximization-based decoding result in dull and repetitive generated texts even when using powerful language models (Holtzman et al., 2019). Adding a loss function for regularization was shown to improve text generation output by helping avoid unwanted properties, such as contradiction or repetition (Li at al., 2020). In this work, we propose fine-tuning a language model by using policy gradient reinforcement learning, directly optimizing for better generation. We apply this approach to minimizing repetition in generated text, and show that, when combined with unlikelihood training (Welleck et al., 2020), our method further reduces repetition without impacting the language model quality. We also evaluate other methods for improving generation at training and decoding time, and compare them using various metrics aimed at control for better text generation output.


title: Explaining the Black-box Smoothly- A Counterfactual Approach

abstract: We propose a BlackBox Counterfactual Explainer, designed to explain image classification models for medical applications. Classical approaches (e.g., saliency maps) that assess feature importance do not explain "how" imaging features in important anatomical regions are relevant to the classification decision. Our framework explains the decision for a target class by gradually "exaggerating" the semantic effect of the class in a query image. We adopted a Generative Adversarial Network (GAN) to generate a progressive set of perturbations to a query image, such that the classification decision changes from its original class to its negation. We used counterfactual explanations from our framework to audit a classifier trained on a chest x-ray dataset with multiple labels. We proposed clinically-relevant quantitative metrics such as cardiothoracic ratio and the score of a healthy costophrenic recess to evaluate our explanations. We conducted a human-grounded experiment with diagnostic radiology residents to compare different styles of explanations (no explanation, saliency map, cycleGAN explanation, and our counterfactual explanation) by evaluating different aspects of explanations: (1) understandability, (2) classifier's decision justification, (3) visual quality, (d) identity preservation, and (5) overall helpfulness of an explanation to the users. Our results show that our counterfactual explanation was the only explanation method that significantly improved the users' understanding of the classifier's decision compared to the no-explanation baseline. Our metrics established a benchmark for evaluating model explanation methods in medical images. Our explanations revealed that the classifier relied on clinically relevant radiographic features for its diagnostic decisions, thus making its decision-making process more transparent to the end-user.


title: Partial regularity for the optimal $p$-compliance problem with length penalization

abstract: We establish a partial $C^{1,\alpha}$ regularity result for minimizers of the optimal $p$-compliance problem with length penalization in any spatial dimension $N\geq 2$, extending some of the results obtained in [Chambolle-Lamboley-Lemenant-Stepanov 17], [Bulanyi-Lemenant 20]. The key feature is that the $C^{1,\alpha}$ regularity of minimizers for some free boundary type problem is investigated with a free boundary set of codimension $N-1$. We prove that every optimal set cannot contain closed loops, and it is $C^{1,\alpha}$ regular at $\mathcal{H}^{1}$-a.e. point for every $p\in (N-1,+\infty)$.


title: The First Dynamical Mass Determination of a Nitrogen-rich Wolf-Rayet Star using a Combined Visual and Spectroscopic Orbit

abstract: We present the first visual orbit for the nitrogen-rich Wolf-Rayet binary, WR 133 (WN5o + O9I) based on observations made with the CHARA Array and the MIRC-X combiner. This orbit represents the first visual orbit for a WN star and only the third Wolf-Rayet star with a visual orbit. The orbit has a period of 112.8 d, a moderate eccentricity of 0.36, and a separation of $a$= 0.79 mas on the sky. We combine the visual orbit with an SB2 orbit and Gaia parallax to find that the derived masses of the component stars are $M_{\rm WR}$ = $9.3\pm1.6 M_\odot$ and $M_{\rm O}$ = $22.6\pm 3.2 M_\odot$, with the large errors owing to the nearly face-on geometry of the system combined with errors in the spectroscopic parameters. We also derive an orbital parallax that is identical to the {\it Gaia}-determined distance. We present a preliminary spectral analysis and atmosphere models of the component stars, and find the mass-loss rate in agreement with polarization variability and our orbit. However, the derived masses are low compared to the spectral types and spectral model. Given the close binary nature, we suspect that WR 133 should have formed through binary interactions, and represents an ideal target for testing evolutionary models given its membership in the cluster NGC 6871.


title: Independent Policy Gradient Methods for Competitive Reinforcement Learning

abstract: We obtain global, non-asymptotic convergence guarantees for independent learning algorithms in competitive reinforcement learning settings with two agents (i.e., zero-sum stochastic games). We consider an episodic setting where in each episode, each player independently selects a policy and observes only their own actions and rewards, along with the state. We show that if both players run policy gradient methods in tandem, their policies will converge to a min-max equilibrium of the game, as long as their learning rates follow a two-timescale rule (which is necessary). To the best of our knowledge, this constitutes the first finite-sample convergence result for independent policy gradient methods in competitive RL; prior work has largely focused on centralized, coordinated procedures for equilibrium computation.


title: Empirical Model of 10-130 MeV Solar Energetic Particle Spectra at 1 AU Based on Coronal Mass Ejection Speed and Direction

abstract: We present a new empirical model to predict solar energetic particle (SEP) event-integrated and peak intensity spectra between 10 and 130 MeV at 1 AU, based on multi-point spacecraft measurements from the Solar TErrestrial RElations Observatory (STEREO), the Geostationary Operational Environmental Satellites (GOES) and the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment. The analyzed data sample includes 32 SEP events occurring between 2010 and 2014, with a statistically significant proton signal at energies in excess of a few tens of MeV, unambiguously recorded at three spacecraft locations. The spatial distributions of SEP intensities are reconstructed by assuming an energy-dependent 2D Gaussian functional form, and accounting for the correlation between the intensity and the speed of the parent coronal mass ejection (CME), and the magnetic field line connection angle. The CME measurements used are from the Space Weather Database Of Notifications, Knowledge, Information (DONKI). The model performance, including its extrapolations to lower/higher energies, is tested by comparing with the spectra of 20 SEP events not used to derive the model parameters. Despite the simplicity of the model, the observed and predicted event-integrated and peak intensities at Earth and at the STEREO spacecraft for these events show remarkable agreement, both in the spectral shapes and their absolute values.


title: Flexible Validity Conditions for the Multivariate Mat\'ern Covariance in any Spatial Dimension and for any Number of Components

abstract: This paper addresses the problem of finding parametric constraints that ensure the validity of the multivariate Mat{\'e}rn covariance for modeling the spatial correlation structure of coregionalized variables defined in an Euclidean space. To date, much attention has been given to the bivariate setting, while the multivariate setting has been explored to a limited extent only. The existing conditions often imply severe restrictions on the upper bounds for the collocated correlation coefficients, which makes the multivariate Mat{\'e}rn model appealing for the case of weak spatial cross-dependence only. We provide a collection of sufficient validity conditions for the multivariate Mat{\'e}rn covariance that allows for more flexible parameterizations than those currently available, and prove that one can attain considerably higher upper bounds for the collocated correlation coefficients in comparison with our competitors. We conclude with an illustration on a trivariate geochemical data set and show that our enlarged parametric space yields better fitting performances.


title: Improving entanglement generation rates in trapped ion quantum networks using nondestructive photon measurement and storage

abstract: Long range quantum information processing will require the integration of different technologies to form hybrid architectures combining the strengths of multiple quantum systems. In this work, we propose a hybrid networking architecture designed to improve entanglement rates in quantum networks based on trapped ions. Trapped ions are excellent candidates as network nodes but photon losses make long-distance networking difficult. To overcome some losses and extend the range of trapped-ion-based networks, we propose including neutral-atom-based non-destructive single-photon detection and single photon storage in between networking nodes, forming a hybrid network. This work builds on recently demonstrated optical frequency conversion of single photons emitted by trapped ions. We derive the average two-node entanglement rate for this proposed network architecture as a function of distance. Using reasonable experimental parameters, we show this proposed quantum network can generate remote entanglement rates up to a factor of 100 larger than that of an equivalent homogeneous network at both near-IR and C-band wavelengths for distances up to 50 km.


title: Solving Common-Payoff Games with Approximate Policy Iteration

abstract: For artificially intelligent learning systems to have widespread applicability in real-world settings, it is important that they be able to operate decentrally. Unfortunately, decentralized control is difficult -- computing even an epsilon-optimal joint policy is a NEXP complete problem. Nevertheless, a recently rediscovered insight -- that a team of agents can coordinate via common knowledge -- has given rise to algorithms capable of finding optimal joint policies in small common-payoff games. The Bayesian action decoder (BAD) leverages this insight and deep reinforcement learning to scale to games as large as two-player Hanabi. However, the approximations it uses to do so prevent it from discovering optimal joint policies even in games small enough to brute force optimal solutions. This work proposes CAPI, a novel algorithm which, like BAD, combines common knowledge with deep reinforcement learning. However, unlike BAD, CAPI prioritizes the propensity to discover optimal joint policies over scalability. While this choice precludes CAPI from scaling to games as large as Hanabi, empirical results demonstrate that, on the games to which CAPI does scale, it is capable of discovering optimal joint policies even when other modern multi-agent reinforcement learning algorithms are unable to do so. Code is available at https://github.com/ssokota/capi .


title: Categories of Nets

abstract: We present a unified framework for Petri nets and various variants, such as pre-nets and Kock's whole-grain Petri nets. Our framework is based on a less well-studied notion that we call $\Sigma$-nets, which allow finer control over whether tokens are treated using the collective or individual token philosophy. We describe three forms of execution semantics in which pre-nets generate strict monoidal categories, $\Sigma$-nets (including whole-grain Petri nets) generate symmetric strict monoidal categories, and Petri nets generate commutative monoidal categories, all by left adjoint functors. We also construct adjunctions relating these categories of nets to each other, in particular showing that all kinds of net can be embedded in the unifying category of $\Sigma$-nets, in a way that commutes coherently with their execution semantics.