title: Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

abstract: We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms.


title: Beyond NDCG: behavioral testing of recommender systems with RecList

abstract: As with most Machine Learning systems, recommender systems are typically evaluated through performance metrics computed over held-out data points. However, real-world behavior is undoubtedly nuanced: ad hoc error analysis and deployment-specific tests must be employed to ensure the desired quality in actual deployments. In this paper, we propose RecList, a behavioral-based testing methodology. RecList organizes recommender systems by use case and introduces a general plug-and-play procedure to scale up behavioral testing. We demonstrate its capabilities by analyzing known algorithms and black-box commercial systems, and we release RecList as an open source, extensible package for the community.


title: Deep IDA: A Deep Learning Method for Integrative Discriminant Analysis of Multi-View Data with Feature Ranking -- An Application to COVID-19 severity

abstract: COVID-19 severity is due to complications from SARS-Cov-2 but the clinical course of the infection varies for individuals, emphasizing the need to better understand the disease at the molecular level. We use clinical and multiple molecular data (or views) obtained from patients with and without COVID-19 who were (or not) admitted to the intensive care unit to shed light on COVID-19 severity. Methods for jointly associating the views and separating the COVID-19 groups (i.e., one-step methods) have focused on linear relationships. The relationships between the views and COVID-19 patient groups, however, are too complex to be understood solely by linear methods. Existing nonlinear one-step methods cannot be used to identify signatures to aid in our understanding of the complexity of the disease. We propose Deep IDA (Integrative Discriminant Analysis) to address analytical challenges in our problem of interest. Deep IDA learns nonlinear projections of two or more views that maximally associate the views and separate the classes in each view, and permits feature ranking for interpretable findings. Our applications demonstrate that Deep IDA has competitive classification rates compared to other state-of-the-art methods and is able to identify molecular signatures that facilitate an understanding of COVID-19 severity.


title: Null-Controllability of a Non-Local Heat Equation

abstract: We consider the null-controllability of a non-local heat equation by interior $L^2(\Omega)$ controls. We confirm a conjecture of Lissy and Zuazua by showing that it is enough to assume that the kernel $k(x,\xi)$ is symmetric and $k(x,\xi)\in L^2(\Omega\times\Omega)$ in order to obtain the result. The result is obtained by treating the non-local linear operator $K:L^2(\Omega)\rightarrow L^2(\Omega)$ as a compact and bounded perturbation of the Dirichlet-Laplacian, $A$, which leads to useful bounds on the semigroup generated by $(A+K,\mathcal{D}(A))$. Using this approach, we are also able to estimate the cost of control.


title: Freezing of Gauge Symmetries in the Heterotic String on $T^4$

abstract: We derive a map relating the gauge symmetry groups of heterotic strings on $T^4$ to other components of the moduli space with rank reduction. This generalizes the results for $T^2$ and $T^3$ which mirror the singularity freezing mechanism of K3 surfaces in F and M-theory, respectively. The novel feature in six dimensions is that the map explicitly involves the topology of the gauge groups, in particular acting only on non-simply-connected ones. This relation is equivalent to that of connected components of the moduli space of flat $G$-bundles over $T^2$ with $G$ non-simply-connected. These results are verified with a reasonably exhaustive list of gauge groups obtained with a moduli space exploration algorithm.


title: Differentiable quantum computational chemistry with PennyLane

abstract: This work describes the theoretical foundation for all quantum chemistry functionality in PennyLane, a quantum computing software library specializing in quantum differentiable programming. We provide an overview of fundamental concepts in quantum chemistry, including the basic principles of the Hartree-Fock method. A flagship feature in PennyLane is the differentiable Hartree-Fock solver, allowing users to compute exact gradients of molecular Hamiltonians with respect to nuclear coordinates and basis set parameters. PennyLane provides specialized operations for quantum chemistry, including excitation gates as Givens rotations and templates for quantum chemistry circuits. Moreover, built-in simulators exploit sparse matrix techniques for representing molecular Hamiltonians that lead to fast simulation for quantum chemistry applications. In combination with PennyLane's existing methods for constructing, optimizing, and executing circuits, these methods allow users to implement a wide range of quantum algorithms for quantum chemistry. We discuss how PennyLane can be used to implement variational algorithms for calculating ground-state energies, excited-state energies, and energy derivatives, all of which can be differentiated with respect to both circuit and Hamiltonian parameters. We provide an example workflow describing how to jointly optimize circuit parameters, nuclear coordinates, and basis set parameters for quantum chemistry algorithms. We discuss a functionality for reducing the number of qubits by using symmetries and explain how PennyLane can be used to estimate quantum resources needed to implement several quantum algorithms. By combining insights from quantum computing, computational chemistry, and machine learning, PennyLane is the first library for differentiable quantum computational chemistry.


title: Step: a tool to perform tests of smoothness on differential distributions based on expansion of polynomials

abstract: We motivate and describe a method based on fits with polynomials to test the smoothness of differential distributions. As a demonstration, we apply the method to several measurements of inclusive jet double-differential cross section in the jet transverse momentum and rapidity at the Tevatron and LHC. This method opens new possibilities to test the quality of differential distributions used for the extraction of physics quantities such as the strong coupling.


title: Simulating thermal qubits through thermofield dynamics

abstract: Quantum computing has attracted the attention of the scientific community in the past few decades. The development of quantum computers promises one path toward safer and faster ways to treat, extract, and transfer information. However, despite the significant advantages of quantum computing, the development of quantum devices operating at room temperature has been compromised by the thermal decoherence process. In addition, in most undergraduate and graduate quantum mechanics courses, the study of thermofield dynamics is usually neglected. In this scenario, this work explores a didactic approach to simulate thermal qubit systems through Thermofield Dynamics (TFD), applied in a quantum computing setup. The results show that the Bloch sphere representation for a qubit can be written in terms of the Bogoliubov transformation, which allows a practical construction for the thermal qubits in a quantum computing setup. Therefore, this work introduces thermofield dynamics through quantum computing to teachers and curious students interested in teaching and learning this important field of studying the temperature impacts on quantum protocols using the TFD technique.


title: Search for GeV-scale Dark Matter Annihilation in the Sun with IceCube DeepCore

abstract: The Sun provides an excellent target for studying spin-dependent dark matter-proton scattering due to its high matter density and abundant hydrogen content. Dark matter particles from the Galactic halo can elastically interact with Solar nuclei, resulting in their capture and thermalization in the Sun. The captured dark matter can annihilate into Standard Model particles including an observable flux of neutrinos. We present the results of a search for low-energy ($<$ 500 GeV) neutrinos correlated with the direction of the Sun using 7 years of IceCube data. This work utilizes, for the first time, new optimized cuts to extend IceCube's sensitivity to dark matter mass down to 5 GeV. We find no significant detection of neutrinos from the Sun. Our observations exclude capture by spin-dependent dark matter-proton scattering with cross-section down to a few times $10^{-41}$ cm$^2$, assuming there is equilibrium with annihilation into neutrinos/anti-neutrinos for dark matter masses between 5 GeV and 100 GeV. These are the strongest constraints at GeV energies for dark matter annihilation directly to neutrinos.


title: Learning Robust Output Control Barrier Functions from Safe Expert Demonstrations

abstract: This paper addresses learning safe output feedback control laws from partial observations of expert demonstrations. We assume that a model of the system dynamics and a state estimator are available along with corresponding error bounds, e.g., estimated from data in practice. We first propose robust output control barrier functions (ROCBFs) as a means to guarantee safety, as defined through controlled forward invariance of a safe set. We then formulate an optimization problem to learn ROCBFs from expert demonstrations that exhibit safe system behavior, e.g., data collected from a human operator or an expert controller. When the parametrization of the ROCBF is linear, then we show that, under mild assumptions, the optimization problem is convex. Along with the optimization problem, we provide verifiable conditions in terms of the density of the data, smoothness of the system model and state estimator, and the size of the error bounds that guarantee validity of the obtained ROCBF. Towards obtaining a practical control algorithm, we propose an algorithmic implementation of our theoretical framework that accounts for assumptions made in our framework in practice. We validate our algorithm in the autonomous driving simulator CARLA and demonstrate how to learn safe control laws from simulated RGB camera images.


title: COVID-19 Detection on Chest X-Ray Images: A comparison of CNN architectures and ensembles

abstract: COVID-19 quickly became a global pandemic after only four months of its first detection. It is crucial to detect this disease as soon as possible to decrease its spread. The use of chest X-ray (CXR) images became an effective screening strategy, complementary to the reverse transcription-polymerase chain reaction (RT-PCR). Convolutional neural networks (CNNs) are often used for automatic image classification and they can be very useful in CXR diagnostics. In this paper, 21 different CNN architectures are tested and compared in the task of identifying COVID-19 in CXR images. They were applied to the COVIDx8B dataset, a large COVID-19 dataset with 16,352 CXR images coming from patients of at least 51 countries. Ensembles of CNNs were also employed and they showed better efficacy than individual instances. The best individual CNN instance results were achieved by DenseNet169, with an accuracy of 98.15% and an F1 score of 98.12%. These were further increased to 99.25% and 99.24%, respectively, through an ensemble with five instances of DenseNet169. These results are higher than those obtained in recent works using the same dataset.


title: A solvable nonlinear autonomous recursion of arbitrary order

abstract: The initial-values problem of the following nonlinear autonomous recursion of order p , z (s + p) = c product of [z (s + l)]^a_l ; with p an arbitrary positive integer, z (s) the dependent variable (possibly a complex number), s the independent variable (a non negative integer), c an arbitrarily assigned, possibly complex, number, and the p exponents a_l arbitrarily assigned integers (positive, negative or vanishing, so that the right-hand side of the recursion be univalent)|is solvable by algebraic operations, involving the solution of a system of linear algebraic equations (generally explicitly solvable) and of a single polynomial equation of degree p (hence explicitly solvable for p = 1; 2; 3; 4 ).


title: Optical Confirmation of X-ray selected Galaxy clusters from the Swift AGN and Cluster survey with MDM and Pan-STARRS Data (Paper III)

abstract: To understand structure formation in the universe and impose stronger constraints on the cluster mass function and cosmological models, it is important to have large galaxy cluster catalogs. The Swift AGN and Cluster Survey is a serendipitous X-ray survey aimed at building a large statistically selected X-ray cluster catalog with 442 cluster candidates in its first release. Our initial SDSS follow-up study confirmed $50\%$ of clusters in the SDSS footprint as z $<$ 0.5 clusters. Here, we present further optical follow-up analysis of 248 (out of 442) cluster candidates from the Swift cluster catalog using multi-band imaging from the MDM $2.4m$ telescope and the Pan-STARRS survey. We report the optical confirmation of 55 clusters with $> 3\sigma$ galaxy overdensities and detectable red sequences in the color-magnitude space. The majority of these confirmed clusters have redshifts z $<$ 0.6. The remaining candidates are potentially higher redshift clusters that are excellent targets for infrared observations. We report the X-ray luminosity and the optical richness for these confirmed clusters. We also discuss the distinction between X-ray and optical observables for the detected and non-detected cluster candidates.


title: Malfustection: Obfuscated Malware Detection and Malware Classification with Data Shortage by Combining Semi-Supervised and Contrastive Learning

abstract: With the advent of new technologies, using various formats of digital gadgets is becoming widespread. In today's world, where everyday tasks are inevitable without technology, this extensive use of computers paves the way for malicious activity. As a result, it is important to provide solutions to defend against these threats. Malware is one of the well-known and widely used means utilized for doing destructive activities by malicious attackers. Producing malware from scratch is somewhat difficult, so attackers tend to obfuscate existing malware and prepare it to become an unrecognizable program. Since creating new malware from an old one using obfuscation is a creative task, there are some drawbacks to identifying obfuscated malwares. In this research, we propose a solution to overcome this problem by converting the code to an image in the first step and then using a semi-supervised approach combined with contrastive learning. In this case, an obfuscation in the malware bytecode corresponds to an augmentation in the image. Hence, by utilizing meaningful augmentations, which simulate some obfuscation changes and combine them to generate complex ambiguity procedures, our proposed solution is able to construct, learn, and detect a wide range of obfuscations. This work addresses two issues: 1) malware classification despite the data deficiency and 2) obfuscated malware detection by training on non-obfuscated malwares. According to the results, the proposed method overcomes the data shortage problem in malware classification, as its accuracy is 90.1% when just 10% of data is used for training the model. Moreover, training on basic malwares without obfuscation achieved 96.21 percent accuracy in detecting obfuscated malware.


title: M2A: Motion Aware Attention for Accurate Video Action Recognition

abstract: Advancements in attention mechanisms have led to significant performance improvements in a variety of areas in machine learning due to its ability to enable the dynamic modeling of temporal sequences. A particular area in computer vision that is likely to benefit greatly from the incorporation of attention mechanisms in video action recognition. However, much of the current research's focus on attention mechanisms have been on spatial and temporal attention, which are unable to take advantage of the inherent motion found in videos. Motivated by this, we develop a new attention mechanism called Motion Aware Attention (M2A) that explicitly incorporates motion characteristics. More specifically, M2A extracts motion information between consecutive frames and utilizes attention to focus on the motion patterns found across frames to accurately recognize actions in videos. The proposed M2A mechanism is simple to implement and can be easily incorporated into any neural network backbone architecture. We show that incorporating motion mechanisms with attention mechanisms using the proposed M2A mechanism can lead to a +15% to +26% improvement in top-1 accuracy across different backbone architectures, with only a small increase in computational complexity. We further compared the performance of M2A with other state-of-the-art motion and attention mechanisms on the Something-Something V1 video action recognition benchmark. Experimental results showed that M2A can lead to further improvements when combined with other temporal mechanisms and that it outperforms other motion-only or attention-only mechanisms by as much as +60% in top-1 accuracy for specific classes in the benchmark.


title: Eigenvalue topology of non-Hermitian band structures in two and three dimensions

abstract: In the band theory for non-Hermitian systems, the energy eigenvalues, which are complex, can exhibit non-trivial topology which is not present in Hermitian systems. In one dimension, it was recently noted theoretically and demonstrated experimentally that the eigenvalue topology is classified by the braid group. The classification of eigenvalue topology in higher dimensions, however, remained an open question. Here, we give a complete description of eigenvalue topology in two and three dimensional systems, including the gapped and gapless cases. We reduce the topological classification problem to a purely computational problem in algebraic topology. In two dimensions, the Brillouin zone torus is punctured by exceptional points, and each nontrivial loop in the punctured torus acquires a braid group invariant. These braids satisfy the constraint that the composite of the braids around the exceptional points is equal to the commutator of the braids on the fundamental cycles of the torus. In three dimensions, there are exceptional knots and links, and the classification depends on how they are embedded in the Brillouin zone three-torus. When the exceptional link is contained in a contractible ball, the classification can be expressed in terms of the knot group of the link. Our results provide a comprehensive understanding of non-Hermitian eigenvalue topology in higher dimensional systems, and should be important for the further explorations of topologically robust open quantum and classical systems.


title: Four-valued logics of truth, non-falsity, exact truth, and material equivalence

abstract: The four-valued semantics of Belnap--Dunn logic, consisting of the truth values True, False, Neither, and Both, gives rise to several non-classical logics depending on which feature of propositions we wish to preserve: truth, non-falsity, or exact truth (truth and non-falsity). Interpreting equality of truth values in this semantics as material equivalence of propositions, we can moreover see the equational consequence relation of this four-element algebra as a logic of material equivalence. In this paper we axiomatize all combinations of these four-valued logics, for example the logic of truth and exact truth or the logic of truth and material equivalence. These combined systems are consequence relations which allow us to express implications involving more than one of these features of propositions.


title: Leggett-Garg inequalities in the quantum field theory of neutrino oscillations

abstract: We investigate Leggett-Garg temporal inequalities in flavor-mixing processes. We derive an exact flavor-mass uncertainty product and we establish that it is an upper bound to the violation of the inequalities. This finding relates temporal nonclassicality to quantum uncertainty and provides a time analog of the Tsirelson upper bound to the violation of the spatial Bell inequalities. By studying the problem both in the exact field-theoretical setting and in the limiting quantum mechanical approximation, we show that Leggett-Garg inequalities are violated more strongly in quantum field theory than in quantum mechanics.


title: An imaging polarimetry survey of Type Ia supernovae: are peculiar extinction and polarization properties produced by circumstellar or interstellar matter?

abstract: Some highly reddened Type Ia supernovae (SNe Ia) display low total-to-selective extinction ratios ($R_V \lesssim 2$) in comparison to that of typical Milky Way dust ($R_V \approx 3.3$), and polarization curves that rise steeply to blue wavelengths, with peak polarization values at short wavelengths ($\lambda_{\rm max} < 0.4$ ${\mu}$m) in comparison to the typical Galactic values ($\lambda_{\rm max} \approx 0.55$ ${\mu}$m). Understanding the source of these properties could provide insight into the progenitor systems of SNe Ia. We aim to determine whether they are the result of the host galaxy's interstellar dust or circumstellar dust. This is accomplished by analysing the continuum polarization of 66 SNe Ia in dust-rich spiral galaxies and 13 SNe Ia in dust-poor elliptical galaxies as a function of normalised galactocentric distance. We find that there is a general trend of SNe Ia in spiral galaxies displaying increased polarization values when located closer to the host galaxies' centre, while SNe Ia in elliptical host galaxies display low polarization. Furthermore, all highly polarized SNe Ia in spiral host galaxies display polarization curves rising toward blue wavelengths, while no evidence of such polarization properties is shown in elliptical host galaxies. This indicates that the source of the peculiar polarization curves is likely the result of interstellar material as opposed to circumstellar material. The peculiar polarization and extinction properties observed toward some SNe Ia may be explained by the radiative torque disruption mechanism induced by the SN or the interstellar radiation field.


title: Adaptive ratchets and the evolution of molecular complexity

abstract: Biological systems have evolved to amazingly complex states, yet we do not understand in general how evolution operates to generate increasing genetic and functional complexity. Molecular recognition sites are short genome segments or peptides binding a cognate recognition target of sufficient sequence similarity. Such sites are simple, ubiquitous modules of sequence information, cellular function, and evolution. Here we show that recognition sites, if coupled to a time-dependent target, can rapidly evolve to complex states with larger code length and smaller coding density than sites recognising a static target. The underlying fitness model contains selection for recognition, which depends on the sequence similarity between site and target, and a uniform cost per unit of code length. Site sequences are shown to evolve in a specific adaptive ratchet, which produces selection of different strength for code extensions and compressions. Ratchet evolution increases the adaptive width of evolved sites, accelerating the adaptation to moving targets and facilitating refinement and innovation of recognition functions. We apply these results to the recognition of fast-evolving antigens by the human immune system. Our analysis shows how molecular complexity can evolve as a collateral to selection for function in a dynamic environment.


title: Second-Order Mirror Descent: Convergence in Games Beyond Averaging and Discounting

abstract: In this paper, we propose a second-order extension of the continuous-time game-theoretic mirror descent (MD) dynamics, referred to as MD2, which provably converges to mere (but not necessarily strict) variationally stable states (VSS) without using common auxiliary techniques such as time-averaging or discounting. We show that MD2 enjoys no-regret as well as an exponential rate of convergence towards strong VSS upon a slight modification. MD2 can also be used to derive many novel continuous-time primal-space dynamics. We then use stochastic approximation techniques to provide a convergence guarantee of discrete-time MD2 with noisy observations towards interior mere VSS. Selected simulations are provided to illustrate our results.

categories: math.OC cs.GT cs.LG cs.SY eess.SY math.DS


title: Towards Measuring Fairness in Speech Recognition: Casual Conversations Dataset Transcriptions

abstract: It is well known that many machine learning systems demonstrate bias towards specific groups of individuals. This problem has been studied extensively in the Facial Recognition area, but much less so in Automatic Speech Recognition (ASR). This paper presents initial Speech Recognition results on "Casual Conversations" -- a publicly released 846 hour corpus designed to help researchers evaluate their computer vision and audio models for accuracy across a diverse set of metadata, including age, gender, and skin tone. The entire corpus has been manually transcribed, allowing for detailed ASR evaluations across these metadata. Multiple ASR models are evaluated, including models trained on LibriSpeech, 14,000 hour transcribed, and over 2 million hour untranscribed social media videos. Significant differences in word error rate across gender and skin tone are observed at times for all models. We are releasing human transcripts from the Casual Conversations dataset to encourage the community to develop a variety of techniques to reduce these statistical biases.


title: Fixed points and Langlands parameter spaces for real groups

abstract: This note records that the Langlands parameter spaces, associated by Adams- Barbasch-Vogan to a real group, may be described as homotopy fixed points (fixed point stacks) of the spaces associated to the corresponding complex group.


title: DeMFI: Deep Joint Deblurring and Multi-Frame Interpolation with Flow-Guided Attentive Correlation and Recursive Boosting

abstract: In this paper, we propose a novel joint deblurring and multi-frame interpolation (DeMFI) framework, called DeMFI-Net, which accurately converts blurry videos of lower-frame-rate to sharp videos at higher-frame-rate based on flow-guided attentive-correlation-based feature bolstering (FAC-FB) module and recursive boosting (RB), in terms of multi-frame interpolation (MFI). The DeMFI-Net jointly performs deblurring and MFI where its baseline version performs feature-flow-based warping with FAC-FB module to obtain a sharp-interpolated frame as well to deblur two center-input frames. Moreover, its extended version further improves the joint task performance based on pixel-flow-based warping with GRU-based RB. Our FAC-FB module effectively gathers the distributed blurry pixel information over blurry input frames in feature-domain to improve the overall joint performances, which is computationally efficient since its attentive correlation is only focused pointwise. As a result, our DeMFI-Net achieves state-of-the-art (SOTA) performances for diverse datasets with significant margins compared to the recent SOTA methods, for both deblurring and MFI. All source codes including pretrained DeMFI-Net are publicly available at https://github.com/JihyongOh/DeMFI.


title: Boost Distribution System Restoration with Emergency Communication Vehicles Considering Cyber-Physical Interdependence

abstract: Enhancing restoration capabilities of distribution systems is one of the main strategies for resilient power systems to cope with extreme events. However, most of the existing studies assume the communication infrastructures are intact for distribution automation, which is unrealistic. Motivated by the applications of the emergency communication vehicles (ECVs) in quickly setting up wireless communication networks after disasters, in this paper, we propose an integrated distribution system restoration (DSR) framework and optimization models, which can coordinate the repair crews, the distribution system (physical sectors), and the emergency communication (cyber sectors) to pick up unserved power loads as quickly as possible. Case studies validated the effectiveness of the proposed models and proved the benefit of considering ECVs and cyber-physical interdependencies in DSR.


title: On graphs with unique geoodesics and antipodes

abstract: In 1962, Oystein Ore asked in which graphs there is exactly one geodesic between any two vertices. He called such graphs geodetic. In this paper, we systematically study properties of geodetic graphs, and also consider antipodal graphs, in which each vertex has exactly one antipode (a farthest vertex). We find necessary and sufficient conditions for a graph to be geodetic or antipodal, obtain results related to algorithmic construction, and find interesting families of Hamiltonian geodetic graphs. By introducing and describing the maximal hereditary subclasses and the minimal hereditary superclasses of the geodetic and antipodal graphs, we get close to the goal of our research -- a constructive classification of these graphs.


title: Inner workings of fractional quantum Hall parent Hamiltonians: A matrix product state point of view

abstract: We study frustration-free Hamiltonians of fractional quantum Hall (FQH) states from the point of view of the matrix product state (MPS) representation of their ground and excited states. There is a wealth of solvable models relating to FQH physics, which, however, is mostly derived and analyzed from the vantage point of first-quantized "analytic clustering properties". In contrast, one obtains long-ranged frustration-free lattice models when these Hamiltonians are studied in an orbital basis, which is the natural basis for the MPS representation of FQH states. The connection between MPS-like states and frustration-free parent Hamiltonians is the central guiding principle in the construction of solvable lattice models, but thus far, only for short-range Hamiltonians and MPSs of finite bond dimension. The situation in the FQH context is fundamentally different. Here we expose the direct link between the infinite-bond-dimension MPS structure of Laughlin-conformal field theory (CFT) states and their parent Hamiltonians. While focusing on the Laughlin state, generalizations to other CFT-MPSs will become transparent.


title: Sequential anomaly detection with sampling constraints

abstract: The problem of sequential anomaly detection is considered, where multiple data sources are monitored in real time and the goal is to identify the "anomalous" ones among them, when it is not possible to sample all sources at all times. A detection scheme in this context requires specifying not only when to stop sampling and which sources to identify as anomalous upon stopping, but also which sources to sample at each time instance until stopping. A novel formulation for this problem is proposed, in which the number of anomalous sources is not necessarily known in advance and the number of sampled sources per time instance is not necessarily fixed. Instead, an arbitrary lower bound and an arbitrary upper bound are assumed on the number of anomalous sources, and the fraction of the expected number of samples over the expected time until stopping is required to not exceed an arbitrary, user-specified level. In addition to this sampling constraint, the probabilities of at least one false alarm and at least one missed detection are controlled below user-specified tolerance levels. A general criterion is established for a policy to achieve the minimum expected time until stopping to a first-order asymptotic approximation as both familywise error rates go to zero. This criterion is used to prove the asymptotic optimality of a family of policies that sample each source at each time instance with a probability that depends on the past observations only through the current estimate of the subset of anomalous sources. In particular, the asymptotic optimality is established of a policy that requires minimal computation under any setup of the problem.


title: Gaussian Determinantal Processes: a new model for directionality in data

abstract: Determinantal point processes (a.k.a. DPPs) have recently become popular tools for modeling the phenomenon of negative dependence, or repulsion, in data. However, our understanding of an analogue of a classical parametric statistical theory is rather limited for this class of models. In this work, we investigate a parametric family of Gaussian DPPs with a clearly interpretable effect of parametric modulation on the observed points. We show that parameter modulation impacts the observed points by introducing directionality in their repulsion structure, and the principal directions correspond to the directions of maximal (i.e. the most long ranged) dependency. This model readily yields a novel and viable alternative to Principal Component Analysis (PCA) as a dimension reduction tool that favors directions along which the data is most spread out. This methodological contribution is complemented by a statistical analysis of a spiked model similar to that employed for covariance matrices as a framework to study PCA. These theoretical investigations unveil intriguing questions for further examination in random matrix theory, stochastic geometry and related topics.

categories: stat.ML cond-mat.dis-nn cs.LG math.PR math.ST stat.TH


title: Sketch-based Creativity Support Tools using Deep Learning

abstract: Sketching is a natural and effective visual communication medium commonly used in creative processes. Recent developments in deep-learning models drastically improved machines' ability in understanding and generating visual content. An exciting area of development explores deep-learning approaches used to model human sketches, opening opportunities for creative applications. This chapter describes three fundamental steps in developing deep-learning-driven creativity support tools that consumes and generates sketches: 1) a data collection effort that generated a new paired dataset between sketches and mobile user interfaces; 2) a sketch-based user interface retrieval system adapted from state-of-the-art computer vision techniques; and, 3) a conversational sketching system that supports the novel interaction of a natural-language-based sketch/critique authoring process. In this chapter, we survey relevant prior work in both the deep-learning and human-computer-interaction communities, document the data collection process and the systems' architectures in detail, present qualitative and quantitative results, and paint the landscape of several future research directions in this exciting area.


title: The Thousand-Pulsar-Array programme on MeerKAT VII: Polarisation properties of pulsars in the Magellanic Clouds

abstract: The Magellanic Clouds are the only external galaxies known to host radio pulsars. The dispersion and rotation measures of pulsars in the Clouds can aid in understanding their structure, and studies of the pulsars themselves can point to potential differences between them and their Galactic counterparts. We use the high sensitivity of the MeerKAT telescope to observe 17 pulsars in the Small and Large Magellanic Clouds in addition to five foreground (Galactic) pulsars. We provide polarization profiles for 18 of these pulsars, improved measurements of their dispersion and rotation measures, and derive the mean parallel magnetic field along the lines of sight. The results are broadly in agreement with expectations for the structure and strength of the magnetic field in the Large and Small Magellanic Clouds. The Magellanic Cloud pulsars have profiles which are narrower than expected from the period-width relationship and we show this is due to selection effects in pulsar surveys rather than any intrinsic difference between the population of Galactic and Magellanic objects.


title: Esophageal virtual disease landscape using mechanics-informed machine learning

abstract: The pathogenesis of esophageal disorders is related to the esophageal wall mechanics. Therefore, to understand the underlying fundamental mechanisms behind various esophageal disorders, it is crucial to map the esophageal wall mechanics-based parameters onto physiological and pathophysiological conditions corresponding to altered bolus transit and supraphysiologic IBP. In this work, we present a hybrid framework that combines fluid mechanics and machine learning to identify the underlying physics of the various esophageal disorders and maps them onto a parameter space which we call the virtual disease landscape (VDL). A one-dimensional inverse model processes the output from an esophageal diagnostic device called endoscopic functional lumen imaging probe (EndoFLIP) to estimate the mechanical "health" of the esophagus by predicting a set of mechanics-based parameters such as esophageal wall stiffness, muscle contraction pattern and active relaxation of esophageal walls. The mechanics-based parameters were then used to train a neural network that consists of a variational autoencoder (VAE) that generates a latent space and a side network that predicts mechanical work metrics for estimating esophagogastric junction motility. The latent vectors along with a set of discrete mechanics-based parameters define the VDL and form clusters corresponding to the various esophageal disorders. The VDL not only distinguishes different disorders but can also be used to predict disease progression in time. Finally, we also demonstrate the clinical applicability of this framework for estimating the effectiveness of a treatment and track patient condition after a treatment.


title: Improved description of di-lepton production in $\tau^-\to\nu_\tau P^-$ decays

abstract: Recently, the Belle collaboration reported the first measurements of the $\tau^- \to\nu_{\tau}\pi^-e^+e^-$ branching fraction and the spectrum of the pion-dielectron system. In an analysis previous to Belle's results, we evaluated this branching fraction which turned out to be compatible with that reported by Belle, although with a large uncertainty. This is the motivation to seek for improvements on our previous evaluation of $\tau^- \to \nu_{\tau}\pi^-\ell^+\ell^-$ decays ($\ell=e,\,\mu$). In this paper we improve our calculation of the $WP^-\gamma^*$ vertex by including flavour symmetry breaking effects in the framework of the Resonance Chiral Theory. We impose QCD short-distance behaviour to constrain most parameters and data on the $\pi^-e^+e^-$ spectrum reported by Belle to fix the remaining free ones. As a result, improved predictions for the branching ratios and hadronic/leptonic spectra are reported, in good agreement with observations. Analogous calculations for the strangeness-changing $\tau^- \to \nu_{\tau}K^-\ell^+\ell^-$ transitions are reported for the first time. Albeit one expects the $m_{\pi\mu^+\mu^-}$ spectrum to be measured in Belle-II and the observables with $\ell=e$ can be improved, it is rather unlikely that the $K$ channels can be measured due to the suppression factor $|V_{ud}/V_{us}|^2=0.05$.


title: On Numerical Considerations for Riemannian Manifold Hamiltonian Monte Carlo

abstract: Riemannian manifold Hamiltonian Monte Carlo (RMHMC) is a sampling algorithm that seeks to adapt proposals to the local geometry of the posterior distribution. The specific form of the Hamiltonian used in RMHMC necessitates {\it implicitly-defined} numerical integrators in order to sustain reversibility and volume-preservation, two properties that are necessary to establish detailed balance of RMHMC. In practice, these implicit equations are solved to a non-zero convergence tolerance via fixed-point iteration. However, the effect of these convergence thresholds on the ergodicity and computational efficiency properties of RMHMC are not well understood. The purpose of this research is to elucidate these relationships through numerous case studies. Our analysis reveals circumstances wherein the RMHMC algorithm is sensitive, and insensitive, to these convergence tolerances. Our empirical analysis examines several aspects of the computation: (i) we examine the ergodicity of the RMHMC Markov chain by employing statistical methods for comparing probability measures based on collections of samples; (ii) we investigate the degree to which detailed balance is violated by measuring errors in reversibility and volume-preservation; (iii) we assess the efficiency of the RMHMC Markov chain in terms of time-normalized ESS. In each of these cases, we investigate the sensitivity of these metrics to the convergence threshold and further contextualize our results in terms of comparison against Euclidean HMC. We propose a method by which one may select the convergence tolerance within a Bayesian inference application using techniques of stochastic approximation and we examine Newton's method, an alternative to fixed point iterations, which can eliminate much of the sensitivity of RMHMC to the convergence threshold.


title: LOLNeRF: Learn from One Look

abstract: We present a method for learning a generative 3D model based on neural radiance fields, trained solely from data with only single views of each object. While generating realistic images is no longer a difficult task, producing the corresponding 3D structure such that they can be rendered from different views is non-trivial. We show that, unlike existing methods, one does not need multi-view data to achieve this goal. Specifically, we show that by reconstructing many images aligned to an approximate canonical pose with a single network conditioned on a shared latent space, you can learn a space of radiance fields that models shape and appearance for a class of objects. We demonstrate this by training models to reconstruct object categories using datasets that contain only one view of each subject without depth or geometry information. Our experiments show that we achieve state-of-the-art results in novel view synthesis and high-quality results for monocular depth prediction.


title: Exploring glassy dynamics with Markov state models from graph dynamical neural networks

abstract: Amorphous materials exhibit structural heterogeneities that relax only on long timescales. Using machine learning techniques, we construct a Markov state model (MSM) for model glass formers that coarse-grains the dynamics into a low-dimensional space, in which transitions occur with rates corresponding to the slowest modes of the system. The transition timescale between states is more than an order of magnitude larger than the conventional alpha-relaxation time, and reveals a fragile to strong crossover at the glass transition. The learned map of states assigned to the particles exhibits correlations of a few molecular diameters both at liquid and glassy temperatures. We show that the MSM effectively constructs a map of scaled excess Voronoi volume, and the free energy difference between the two states is given exactly by the entropy of the these distributions. These results resonate with classic free volume theories of the glass transition, singling out local packing fluctuations as the slowest relaxing features.


title: Sub-10 ps Minimum Ionizing Particle Detection with Geiger-Mode APDs

abstract: Major advances in silicon pixel detectors, with outstanding timing performance, have recently attracted significant attention in the community. In this work we present and discuss the use of state-of-the-art Geiger-mode APDs, also known as single-photon avalanche diodes (SPADs), for the detection of minimum ionizing particles (MIPs) with best-in-class timing resolution. The SPADs were implemented in standard CMOS technology and integrated with on-chip quenching and recharge circuitry. Two devices in coincidence allowed to measure the time-of-flight of 180 GeV/c momentum pions with a coincidence time resolution of 22 ps FWHM (9.5 ps Gaussian sigma). Radiation hardness measurements, also presented here, highlight the suitability of this family of devices for a wide range of high energy physics (HEP) applications.


title: TnT Attacks! Universal Naturalistic Adversarial Patches Against Deep Neural Network Systems

abstract: Deep neural networks are vulnerable to attacks from adversarial inputs and, more recently, Trojans to misguide or hijack the model's decision. We expose the existence of an intriguing class of spatially bounded, physically realizable, adversarial examples -- Universal NaTuralistic adversarial paTches -- we call TnTs, by exploring the superset of the spatially bounded adversarial example space and the natural input space within generative adversarial networks. Now, an adversary can arm themselves with a patch that is naturalistic, less malicious-looking, physically realizable, highly effective achieving high attack success rates, and universal. A TnT is universal because any input image captured with a TnT in the scene will: i) misguide a network (untargeted attack); or ii) force the network to make a malicious decision (targeted attack). Interestingly, now, an adversarial patch attacker has the potential to exert a greater level of control -- the ability to choose a location-independent, natural-looking patch as a trigger in contrast to being constrained to noisy perturbations -- an ability is thus far shown to be only possible with Trojan attack methods needing to interfere with the model building processes to embed a backdoor at the risk discovery; but, still realize a patch deployable in the physical world. Through extensive experiments on the large-scale visual classification task, ImageNet with evaluations across its entire validation set of 50,000 images, we demonstrate the realistic threat from TnTs and the robustness of the attack. We show a generalization of the attack to create patches achieving higher attack success rates than existing state-of-the-art methods. Our results show the generalizability of the attack to different visual classification tasks (CIFAR-10, GTSRB, PubFig) and multiple state-of-the-art deep neural networks such as WideResnet50, Inception-V3 and VGG-16.


title: Resonance of ellipsoidal billiard trajectories and extremal rational functions

abstract: We study resonant billiard trajectories within quadrics in the $d$-dimensional Euclidean space. We relate them to the theory of approximation, in particular the extremal rational functions on the systems of $d$ intervals on the real line. This fruitful link enables us to prove fundamental properties of the billiard dynamics and to provide a comprehensive study of a large class of non-periodic trajectories of integrable billiards. A key ingredient is a functional-polynomial relation of a generalized Pell type. Applying further these ideas and techniques to $s$-weak billiard trajectories, we come to a functional-polynomial relation of the same generalized Pell type.


title: On the error term in the explicit formula of Riemann--von Mangoldt

abstract: We provide an explicit $O(x\log x/T)$ error term for the Riemann--von Mangoldt formula by making results of Wolke (1983) and Ramar\'e (2016) explicit. We also include applications to primes between consecutive powers, the error term in the prime number theorem and an inequality of Ramanujan.


title: High-speed and single-mode FP laser based on parity-time symmetry

abstract: The ability to manipulate cavity resonant modes is of critical importance in laser physics and applications. By exploiting the parity time (PT) symmetry, we propose and experimentally realize a single-mode FP laser with improved output power and high-speed modulation have been demonstrated. The proposed PT symmetric laser consists of two coupled structurally identical FP resonators. The gain and loss in two FP resonators can be manipulated independently by changing the injection currents. In the PT symmetric FP laser, single-mode operation is accomplished by selectively breaking of PT symmetry depending solely on the relation between gain-loss and coupling. Single-mode lasing with output power of 1.7 dBm and a sidemode suppression ratio (SMSR) exceeding 24 dB is demonstrated. The 3 dB bandwidth of 7.9 GHz is achieved and clear eye-openings were obtained for 2.5 Gbps and 10Gbps NRZ operation over 10 km single-mode fibers. Furthermore, the PT symmetry breaking is experimentally confirmed with measured loss and coupling coefficient of two FP resonators. The influence of cavity length, facet reflectivity, and electrical isolation between two P-side electrodes on the side mode suppression ratio and output optical power is also been demonstrated, paving the way for further improvement of the PT symmetric FP laser.


title: Differentiable Wavetable Synthesis

abstract: Differentiable Wavetable Synthesis (DWTS) is a technique for neural audio synthesis which learns a dictionary of one-period waveforms i.e. wavetables, through end-to-end training. We achieve high-fidelity audio synthesis with as little as 10 to 20 wavetables and demonstrate how a data-driven dictionary of waveforms opens up unprecedented one-shot learning paradigms on short audio clips. Notably, we show audio manipulations, such as high quality pitch-shifting, using only a few seconds of input audio. Lastly, we investigate performance gains from using learned wavetables for realtime and interactive audio synthesis.


title: Aqueous Alteration on Asteroids Simplifies Soluble Organic Matter Mixtures

abstract: Biologically relevant abiotic extraterrestrial soluble organic matter (SOM) has been widely investigated to study the origin of life and the chemical evolution of protoplanetary disks. Synthesis of biologically relevant organics, in particular, seems to require aqueous environments in the early solar system. However, SOM in primitive meteorites includes numerous chemical species besides the biologically relevant ones, and the reaction mechanisms that comprehensively explain the complex nature of SOM are unknown. Besides, the initial reactants, which formed before asteroid accretion, were uncharacterized. We examined the mass distribution of SOM extracted from three distinct Tagish Lake meteorite fragments, which exhibit different degrees of aqueous alteration though they originated from a single asteroid. We report that mass distributions of SOM in the primordial fragments are well fit by the SchulzZimm (SZ) model for the molecular weight distribution patterns found in chain growth polymerization experiments. Also, the distribution patterns diverge further from SZ with increasing degrees of aqueous alteration. These observations imply that the complex nature of the primordial SOM (1) was established before severe alteration on the asteroid, (2) possibly existed before parent-body accretion, and (3) later became simplified on the asteroid. Therefore, aqueous reactions on asteroids are not required conditions for cultivating complex SOM. Furthermore, we found that overall H over C ratios of SOM decrease with increasing aqueous alteration, and the estimate of H loss from the SOM is 10% to 30%. Organics seem to be a significant H2 source that may have caused subsequent chemical reactions in the Tagish Lake meteorite parent body.


title: Reinforcement Learning with Adaptive Curriculum Dynamics Randomization for Fault-Tolerant Robot Control

abstract: This study is aimed at addressing the problem of fault tolerance of quadruped robots to actuator failure, which is critical for robots operating in remote or extreme environments. In particular, an adaptive curriculum reinforcement learning algorithm with dynamics randomization (ACDR) is established. The ACDR algorithm can adaptively train a quadruped robot in random actuator failure conditions and formulate a single robust policy for fault-tolerant robot control. It is noted that the hard2easy curriculum is more effective than the easy2hard curriculum for quadruped robot locomotion. The ACDR algorithm can be used to build a robot system that does not require additional modules for detecting actuator failures and switching policies. Experimental results show that the ACDR algorithm outperforms conventional algorithms in terms of the average reward and walking distance.


title: Image enhancement in acoustic-resolution photoacoustic microscopy enabled by a novel directional algorithm

abstract: Acoustic-resolution photoacoustic microscopy (AR-PAM) is a promising tool for microvascular imaging. In the focal region, resolution of AR-PAM is determined by the ultrasound transducer and ultimately limited by acoustic diffraction. In the out-of-focus region, resolution deteriorates with increasing distance from the focal plane, which restricts depth of focus (DOF). Besides, a trade-off exists between resolution and DOF. Previously, synthetic aperture focusing technique (SAFT) and/or deconvolution methods have been demonstrated to enhance AR-PAM images. However, they suffer from issues in low resolution, low signal-to-noise ratio (SNR), and/or poor image fidelity. Here, we propose a novel algorithm for AR-PAM to enhance image resolution, SNR, and fidelity. The algorithm consists of a Fourier accumulation SAFT (FA-SAFT) and a directional model-based (D-MB) deconvolution method. Inspired from Fourier denoising technique and directional SAFT, FA-SAFT mainly compensates for the defocusing effect. Besides, D-MB deconvolution enhances the resolution as well as preserves the image fidelity, especially for the objects with line patterns such as microvasculature. Full width at half maximum of 26-31 um over DOF of 1.8 mm and minimum resolvable distance of 46-49 um are experimentally achieved by imaging tungsten wire phantom. Moreover, imaging of leaf skeleton phantom and in vivo imaging of mouse blood vessels also prove that our algorithm is capable of providing high-resolution, high-SNR, and good-fidelity results for complex structures and for in vivo applications.


title: FBNetV5: Neural Architecture Search for Multiple Tasks in One Run

abstract: Neural Architecture Search (NAS) has been widely adopted to design accurate and efficient image classification models. However, applying NAS to a new computer vision task still requires a huge amount of effort. This is because 1) previous NAS research has been over-prioritized on image classification while largely ignoring other tasks; 2) many NAS works focus on optimizing task-specific components that cannot be favorably transferred to other tasks; and 3) existing NAS methods are typically designed to be "proxyless" and require significant effort to be integrated with each new task's training pipelines. To tackle these challenges, we propose FBNetV5, a NAS framework that can search for neural architectures for a variety of vision tasks with much reduced computational cost and human effort. Specifically, we design 1) a search space that is simple yet inclusive and transferable; 2) a multitask search process that is disentangled with target tasks' training pipeline; and 3) an algorithm to simultaneously search for architectures for multiple tasks with a computational cost agnostic to the number of tasks. We evaluate the proposed FBNetV5 targeting three fundamental vision tasks -- image classification, object detection, and semantic segmentation. Models searched by FBNetV5 in a single run of search have outperformed the previous stateof-the-art in all the three tasks: image classification (e.g., +1.3% ImageNet top-1 accuracy under the same FLOPs as compared to FBNetV3), semantic segmentation (e.g., +1.8% higher ADE20K val. mIoU than SegFormer with 3.6x fewer FLOPs), and object detection (e.g., +1.1% COCO val. mAP with 1.2x fewer FLOPs as compared to YOLOX).


title: The compressibility of quark matter under strong magnetic field in the NJL model

abstract: The compressibility of magnetized quark matter is investigated in the SU(2) NJL model. The increases of the chemical potential and the temperature can reduce the compressibility, and lead to the much stiffer equation of state. The variation of the compressibility with the magnetic field will depend on the phase region. Due to the anisotropic structure, the compressibility is different in the directions parallel and perpendicular to the field. The discontinuity of longitudinal compressibility with the chemical potential and the temperature captures the signature of a first-order chiral phase transition and the crossover at high temperature. Moreover, the magnetic-field-and-temperature running coupling would have an important effect on the position of the phase transition. Under the lowest landau level approximation at zero temperature, the longitudinal compressibility has a direct inverse proportional relation to the magnetic field strength and the chemical potential square as $\kappa^\parallel_{\mathrm{LLL},\chi}\propto1/(eB \mu^2)$.


title: ExoMiner: A Highly Accurate and Explainable Deep Learning Classifier that Validates 301 New Exoplanets

abstract: The kepler and TESS missions have generated over 100,000 potential transit signals that must be processed in order to create a catalog of planet candidates. During the last few years, there has been a growing interest in using machine learning to analyze these data in search of new exoplanets. Different from the existing machine learning works, ExoMiner, the proposed deep learning classifier in this work, mimics how domain experts examine diagnostic tests to vet a transit signal. ExoMiner is a highly accurate, explainable, and robust classifier that 1) allows us to validate 301 new exoplanets from the MAST Kepler Archive and 2) is general enough to be applied across missions such as the on-going TESS mission. We perform an extensive experimental study to verify that ExoMiner is more reliable and accurate than the existing transit signal classifiers in terms of different classification and ranking metrics. For example, for a fixed precision value of 99%, ExoMiner retrieves 93.6% of all exoplanets in the test set (i.e., recall=0.936) while this rate is 76.3% for the best existing classifier. Furthermore, the modular design of ExoMiner favors its explainability. We introduce a simple explainability framework that provides experts with feedback on why ExoMiner classifies a transit signal into a specific class label (e.g., planet candidate or not planet candidate).


title: UN-AVOIDS: Unsupervised and Nonparametric Approach for Visualizing Outliers and Invariant Detection Scoring

abstract: The visualization and detection of anomalies (outliers) are of crucial importance to many fields, particularly cybersecurity. Several approaches have been proposed in these fields, yet to the best of our knowledge, none of them has fulfilled both objectives, simultaneously or cooperatively, in one coherent framework. The visualization methods of these approaches were introduced for explaining the output of a detection algorithm, not for data exploration that facilitates a standalone visual detection. This is our point of departure: UN-AVOIDS, an unsupervised and nonparametric approach for both visualization (a human process) and detection (an algorithmic process) of outliers, that assigns invariant anomalous scores (normalized to $[0,1]$), rather than hard binary-decision. The main aspect of novelty of UN-AVOIDS is that it transforms data into a new space, which is introduced in this paper as neighborhood cumulative density function (NCDF), in which both visualization and detection are carried out. In this space, outliers are remarkably visually distinguishable, and therefore the anomaly scores assigned by the detection algorithm achieved a high area under the ROC curve (AUC). We assessed UN-AVOIDS on both simulated and two recently published cybersecurity datasets, and compared it to three of the most successful anomaly detection methods: LOF, IF, and FABOD. In terms of AUC, UN-AVOIDS was almost an overall winner. The article concludes by providing a preview of new theoretical and practical avenues for UN-AVOIDS. Among them is designing a visualization aided anomaly detection (VAAD), a type of software that aids analysts by providing UN-AVOIDS' detection algorithm (running in a back engine), NCDF visualization space (rendered to plots), along with other conventional methods of visualization in the original feature space, all of which are linked in one interactive environment.


title: Game-environment feedback dynamics for voluntary prisoner's dilemma games

abstract: Recently, the eco-evolutionary game theory which describes the coupled dynamics of strategies and environment have attracted great attention. At the same time, most of the current work is focused on the classic two-player two-strategy game. In this work, we study multi-strategy eco-evolutionary game theory which is an extension of the framework. For simplicity, we'll focus on the voluntary participation Prisoner's dilemma game. For the general class of payoff-dependent feedback dynamics, we show the conditions for the existence and stability of internal equilibrium by using the replicator dynamics, respectively. Where internal equilibrium points, such as, two-strategy coexistence states, three-strategy coexistence states, persistent oscillation states and interior saddle points. These states are determined by the relative feedback strength and payoff matrix, and are independent of the relative feedback speed and initial state. In particular, the three-strategy coexistence provides a new mechanism for maintaining biodiversity in biology, ecology, and sociology. Besides, we find that this three-strategy model return to the persistent oscillation state of the two-strategy model when there is no defective strategy at the initial moment.


title: Direct visualization of polaron formation in the thermoelectric SnSe

abstract: SnSe is a layered material that currently holds the record for bulk thermoelectric efficiency. The primary determinant of this high efficiency is thought to be the anomalously low thermal conductivity resulting from strong anharmonic coupling within the phonon system. Here we show that the nature of the carrier system in SnSe is also determined by strong coupling to phonons by directly visualizing polaron formation in the material. We employ ultrafast electron diffraction and diffuse scattering to track the response of phonons in both momentum and time to the photodoping of free carriers across the bandgap, observing the bimodal and anisotropic lattice distortions that drive carrier localization. Relatively large (\SI{18.7}{\angstrom}), quasi-1D polarons are formed on the \SI{300}{\femto\second} timescale with smaller (\SI{4.2}{\angstrom}) 3D polarons taking an order of magnitude longer (\SI{4}{\pico\second}) to form. This difference appears to be a consequence of the profoundly anisotropic electron-phonon coupling in SnSe, with strong Fr\"ohlich coupling only to zone center polar optical phonons. These results demonstrate that carriers in SnSe at optimal doping levels results in a high polaron density and that strong electron-phonon coupling is also critical to the thermoelectric performance of this benchmark material and potentially high-performance thermoelectrics more generally.


title: Seconds-scale coherence on nuclear spin transitions of ultracold polar molecules in 3D optical lattices

abstract: Ultracold polar molecules (UPMs) are emerging as a novel and powerful platform for fundamental applications in quantum science. Here, we report characterization of the coherence between nuclear spin levels of ultracold ground-state sodium-rubidium molecules loaded into a 3D optical lattice with a nearly photon scattering limited trapping lifetime of 9(1) seconds. After identifying and compensating the main sources of decoherence, we achieve a maximum nuclear spin coherence time of $T_2^* = 3.3(6)$~s with two-photon Ramsey spectroscopy. Furthermore, based on the understanding of the main factor limiting the coherence of the two-photon Rabi transition, we obtain a Rabi lineshape with linewidth below 0.8 Hz. The simultaneous realization of long lifetime and coherence time, and ultra-high spectroscopic resolution in our system unveils the great potentials of UPMs in quantum simulation, computation, and metrology.


title: CoCAtt: A Cognitive-Conditioned Driver Attention Dataset

abstract: The task of driver attention prediction has drawn considerable interest among researchers in robotics and the autonomous vehicle industry. Driver attention prediction can play an instrumental role in mitigating and preventing high-risk events, like collisions and casualties. However, existing driver attention prediction models neglect the distraction state and intention of the driver, which can significantly influence how they observe their surroundings. To address these issues, we present a new driver attention dataset, CoCAtt (Cognitive-Conditioned Attention). Unlike previous driver attention datasets, CoCAtt includes per-frame annotations that describe the distraction state and intention of the driver. In addition, the attention data in our dataset is captured in both manual and autopilot modes using eye-tracking devices of different resolutions. Our results demonstrate that incorporating the above two driver states into attention modeling can improve the performance of driver attention prediction. To the best of our knowledge, this work is the first to provide autopilot attention data. Furthermore, CoCAtt is currently the largest and the most diverse driver attention dataset in terms of autonomy levels, eye tracker resolutions, and driving scenarios.


title: Generalized-Hukuhara Subgradient Method for Optimization Problem with Interval-valued Functions and its Application in Lasso Problem

abstract: In this study, a \emph{$gH$-subgradient technique} is developed to obtain efficient solutions to the optimization problems with nonsmooth nonlinear convex interval-valued functions. The algorithmic implementation of the developed $gH$-subgradient technique is illustrated. As an application of the proposed \emph{$gH$-subgradient technique}, an $\ell_1$ penalized linear regression problem, known as a \emph{lasso problem}, with interval-valued features is solved.


title: On Wasserstein-1 distance in the central limit theorem for elephant random walk

abstract: Recently, the elephant random walk has attracted a lot of attentions. A wide range of literature is available for the asymptotic behavior of the process, such as the central limit theorems, functional limit theorems and the law of iterated logarithm. However, there is not result concerning Wassertein-1 distance for the normal approximations.In this paper, we show that the Wassertein-1 distance in the central limit theorem is totally different when a memory parameter $p$ belongs to one of the three cases $0< p < 1/2,$ $1/2< p<3/4$ and $p=3/4.$


title: Rethinking Query, Key, and Value Embedding in Vision Transformer under Tiny Model Constraints

abstract: A vision transformer (ViT) is the dominant model in the computer vision field. Despite numerous studies that mainly focus on dealing with inductive bias and complexity, there remains the problem of finding better transformer networks. For example, conventional transformer-based models usually use a projection layer for each query (Q), key (K), and value (V) embedding before multi-head self-attention. Insufficient consideration of semantic $Q, K$, and $V$ embedding may lead to a performance drop. In this paper, we propose three types of structures for $Q$, $K$, and $V$ embedding. The first structure utilizes two layers with ReLU, which is a non-linear embedding for $Q, K$, and $V$. The second involves sharing one of the non-linear layers to share knowledge among $Q, K$, and $V$. The third proposed structure shares all non-linear layers with code parameters. The codes are trainable, and the values determine the embedding process to be performed among $Q$, $K$, and $V$. Hence, we demonstrate the superior image classification performance of the proposed approaches in experiments compared to several state-of-the-art approaches. The proposed method achieved $71.4\%$ with a few parameters (of $3.1M$) on the ImageNet-1k dataset compared to that required by the original transformer model of XCiT-N12 ($69.9\%$). Additionally, the method achieved $93.3\%$ with only $2.9M$ parameters in transfer learning on average for the CIFAR-10, CIFAR-100, Stanford Cars datasets, and STL-10 datasets, which is better than the accuracy of $92.2\%$ obtained via the original XCiT-N12 model.


title: Superconductivity and bosonic fluid emerging from Moir\'e flat bands

abstract: Although evidence of inter-valley attraction-mediated by phonon or topological fluctuations is accumulating, the origin of superconductivity in the flat-band quantum moir\'e materials remains an open question. Here, instead of attempting to pinpoint the origin of the superconductivity, we aim at identifying universal properties of moir\'e flat bands that shall emerge in the presence of inter-valley attractions. We show that by matching the interaction strength of inter-valley attraction with intra-valley repulsion, the flat-band limit becomes exactly solvable. Away from the flat-band limit, the system can be simulated via quantum Monte Carlo (QMC) methods without sign problem for any fillings. Combining analytic solutions with large-scale numerical simulations, we show that upon increasing temperature, the superconducting phase melts into a bosonic fluid of Cooper pairs with large/diverging compressibility. In contrast to flat-band attractive Hubbard models, where similar effects arise only for on-site interactions, our study indicates this physics is a universal property of moir\'e flat bands, regardless of microscopic details such as the range of interactions and/or spin-oribt couplings. At higher temperature, the boson fluid phase gives its way to a pseudo gap phase, where some Cooper pairs are torn apart by thermal fluctuations, resulting in fermion density of states inside the gap. Unlike the superconducting transition temperature, which is very sensitive to doping and twisting angles, the gap and the temperature scale of the boson fluid phase and the pseudo gap phase are found to be nearly independent of doping level and/or flat-band bandwidth. The relevance of these phases with experimental discoveries in the flat band quantum moir\'e materials is discussed.


title: Noise spectroscopy in Rydberg atomic ensemble

abstract: The phase noise of incident light fields can be converted into amplitude noise via absorption and dispersion effects in media under electromagnetically induced transparency. The conversion process is sensitive to two-photon detuning. This effect can be used to perform a novel type of Rydberg spectroscopy. We present experimental measurements of the phase noise spectrum (PNS) with Rydberg atoms in a cesium vapor cell. Furthermore, to explore microwave field measurements, a modulation-enhanced PNS approach is proposed.


title: Inference of population-level disease transmissibility from household-structured symptom onset data

abstract: First Few X (FFX) studies collect household-stratified data in the early stages of a pandemic, in order to infer severity and transmissibility of an emerging disease. We present a Bayesian method to approximately infer population-level transmissibility for the first time from such data; previous studies have only inferred household-level transmissibility. To do this we perform the inference at two levels, assuming one transmission rate parameter for within-household infection, and another transmission rate parameter for infection between different households. We use a simplifying assumption: that between-household infections always occur in naive households; while still performing full joint inference on the within-household infection parameters. In addition, a novel technique is used to remove systematic bias when the number of new infections per day is growing or decaying, as is common in real outbreaks, without the need for contemporaneous estimates of the serial interval. The method is validated on simulated data and is shown to perform well, even when the number of infected households is relatively small.


title: Achievability and Impossibility of Exact Pairwise Ranking

abstract: We consider the problem of recovering the rank of a set of $n$ items based on noisy pairwise comparisons. We assume the SST class as the family of generative models. Our analysis gave sharp information theoretic upper and lower bound for the exact requirement, which matches exactly in the parametric limit. Our tight analysis on the algorithm induced by the moment method gave better constant in Minimax optimal rate than ~\citet{shah2017simple} and contribute to their open problem. The strategy we used in this work to obtain information theoretic bounds is based on combinatorial arguments and is of independent interest.


title: Concurrent Transmission and Multiuser Detection of LoRa Signals

abstract: This paper investigates a new model to improve the scalability of low-power long-range (LoRa) networks by allowing multiple end devices (EDs) to simultaneously communicate with multiple multi-antenna gateways on the same frequency band and using the same spreading factor. The maximum likelihood (ML) decision rule is first derived for non-coherent detection of information bits transmitted by multiple devices. To overcome the high complexity of the ML detection, we propose a sub-optimal two-stage detection algorithm to balance the computational complexity and error performance. In the first stage, we identify transmit chirps (without knowing which EDs transmit them). In the second stage, we determine the EDs that transmit the specific chirps identified from the first stage. To improve the detection performance in the second stage, we also optimize the transmit powers of EDs to minimize the similarity, measured by the Jaccard coefficient, between the received powers of any pair of EDs. As the power control optimization problem is non-convex, we use concepts from successive convex approximation to transform it to an approximate convex optimization problem that can be solved iteratively and guaranteed to reach a sub-optimal solution. Simulation results demonstrate and justify the tradeoff between transmit power penalties and network scalability of the proposed LoRa network model. In particular, by allowing concurrent transmission of 2 or 3 EDs, the uplink capacity of the proposed network can be doubled or tripled over that of a conventional LoRa network, albeit at the expense of additional 3.0 or 4.7 dB transmit power.


title: UFO: A UniFied TransfOrmer for Vision-Language Representation Learning

abstract: In this paper, we propose a single UniFied transfOrmer (UFO), which is capable of processing either unimodal inputs (e.g., image or language) or multimodal inputs (e.g., the concatenation of the image and the question), for vision-language (VL) representation learning. Existing approaches typically design an individual network for each modality and/or a specific fusion network for multimodal tasks. To simplify the network architecture, we use a single transformer network and enforce multi-task learning during VL pre-training, which includes the image-text contrastive loss, image-text matching loss, and masked language modeling loss based on the bidirectional and the seq2seq attention mask. The same transformer network is used as the image encoder, the text encoder, or the fusion network in different pre-training tasks. Empirically, we observe less conflict among different tasks and achieve new state of the arts on visual question answering, COCO image captioning (cross-entropy optimization) and nocaps (in SPICE). On other downstream tasks, e.g., image-text retrieval, we also achieve competitive performance.


title: Physical and Geometrical Parameters of the Evolved Binary System HD6009

abstract: Atmospheric modeling and dynamical analysis of the components of the visually close binary system (VCBS) HD6009 were used to estimate their individual physical and geometrical parameters. Model atmospheres were constructed using a grid of Kurucz solar metalicity blanketed models and used to compute the individual synthetic spectral energy distribution (SED) for each component separately. These SEDs were combined together to compute the entire SED for the system from the net luminosities of the components $a$ and $b$ located at a distance $d$ from the Earth. %The entire observational SED of the system was used as a reference for the comparison with the synthetic ones. We used the feedback modified parameters and iteration method to get the best fit between synthetic and observational entire SEDs. The physical and geometrical parameters of the system's components were derived as: $T_{\rm eff}^{a} =5625\pm75$\,K, $T_{\rm eff}^{b} =5575\pm75$\,K, log $g_{a}=3.75\pm0.25$, log $g_{b}=3.75\pm0.25$, $R_{a}=2.75\pm0.30 R_\odot$, $R_{b}=2.65\pm0.30 R_\odot$, $M_v^{a}= 2\fm80\pm0.30$, $M_v^{b}=2\fm93\pm0.30$, $M_a= 1.42\pm0.15 M_{\odot}$, $M_b=1.40\pm0.15 M_{\odot}$, $L_a=6.80\pm0.75 L_\odot$, $L_b=6.09\pm0.75 L_\odot$ and $\pi=14.43$mas dynamical parallax. The system is shown to be consist of G6 IV primary and G6 IV secondary components.


title: Dynamically induced conformation depending on excited normal modes of fast oscillation

abstract: We present dynamical effects on conformation in a simple bead-spring model consisting of three beads connected by two stiff springs. The conformation defined by the bending angle between the two springs is determined not only by a given potential energy function depending on the bending angle, but also fast motion of the springs which constructs the effective potential. A conformation corresponding with a local minimum of the effective potential is hence called the dynamically induced conformation. We develop a theory to derive the effective potential by using multiple-scale analysis and the averaging method. A remarkable consequence is that the effective potential depends on the excited normal modes of the springs and amount of the spring energy. Efficiency of the obtained effective potential is numerically verified.


title: IC-U-Net: A U-Net-based Denoising Autoencoder Using Mixtures of Independent Components for Automatic EEG Artifact Removal

abstract: Electroencephalography (EEG) signals are often contaminated with artifacts. It is imperative to develop a practical and reliable artifact removal method to prevent misinterpretations of neural signals and underperformance of brain-computer interfaces. This study developed a new artifact removal method, IC-U-Net, which is based on the U-Net architecture for removing pervasive EEG artifacts and reconstructing brain sources. The IC-U-Net was trained using mixtures of brain and non-brain sources decomposed by independent component analysis and employed an ensemble of loss functions to model complex signal fluctuations in EEG recordings. The effectiveness of the proposed method in recovering brain sources and removing various artifacts (e.g., eye blinks/movements, muscle activities, and line/channel noises) was demonstrated in a simulation study and three real-world EEG datasets collected at rest and while driving and walking. IC-U-Net is user-friendly and publicly available, does not require parameter tuning or artifact type designations, and has no limitations on channel numbers. Given the increasing need to image natural brain dynamics in a mobile setting, IC-U-Net offers a promising end-to-end solution for automatically removing artifacts from EEG recordings.


title: E3NE: An End-to-End Framework for Accelerating Spiking Neural Networks with Emerging Neural Encoding on FPGAs

abstract: Compiler frameworks are crucial for the widespread use of FPGA-based deep learning accelerators. They allow researchers and developers, who are not familiar with hardware engineering, to harness the performance attained by domain-specific logic. There exists a variety of frameworks for conventional artificial neural networks. However, not much research effort has been put into the creation of frameworks optimized for spiking neural networks (SNNs). This new generation of neural networks becomes increasingly interesting for the deployment of AI on edge devices, which have tight power and resource constraints. Our end-to-end framework E3NE automates the generation of efficient SNN inference logic for FPGAs. Based on a PyTorch model and user parameters, it applies various optimizations and assesses trade-offs inherent to spike-based accelerators. Multiple levels of parallelism and the use of an emerging neural encoding scheme result in an efficiency superior to previous SNN hardware implementations. For a similar model, E3NE uses less than 50% of hardware resources and 20% less power, while reducing the latency by an order of magnitude. Furthermore, scalability and generality allowed the deployment of the large-scale SNN models AlexNet and VGG.


title: Primordial black holes from the perturbations in the inflaton potential in peak theory

abstract: The primordial black hole (PBH) is an effective candidate for dark matter. In this work, the PBH abundance $f$ is calculated in peak theory, with one or two perturbations in the inflaton potential. We construct an antisymmetric perturbation that can create a perfect plateau in the inflaton potential, leading inflation to the ultra-slow-roll stage. During this stage, the power spectrum of primordial curvature perturbation is remarkably enhanced on small scales, generating abundant PBHs. The PBH abundance $f\sim 0.1$ can be achieved in one or two typical mass windows at $10^{-17}M_\odot$, $10^{-13}M_\odot$, and $30M_\odot$, without spoiling the nearly scale-invariant power spectrum on large scales. For comparison, $f$ is calculated in two approximate methods of peak theory (with different spectral moments) and also in the Press--Schechter theory. It is found that the Press--Schechter theory systematically underestimates $f$ by two or three orders of magnitude compared with peak theory.


title: Central ideals and Jacobson radicals of blocks of group algebras

abstract: We study the relationship between the central ideals and the Jacobson radicals of blocks of group algebras.In particular, we characterize the blocks with the property that the square of the radical is a non-zero central ideal. Moreover, we also consider the blocks with the property that the cube of the radical becomes a non-zero central ideal.


title: Fluid-structure instability forecasts thoracic aortic aneurysm progression

abstract: The basic mechanism driving aneurysm growth is unknown. Currently, clinical diagnosis of an aneurysm is mainly informed by retrospective tracking of its size and growth rate. However, aneurysms can rupture before reactive criteria are met or remain stable when they are exceeded. Here, we identify a fluid-structure instability that is associated with abnormal aortic dilatation. Our analysis yields a measurable dimensionless number and its analytically derived critical threshold. This threshold pinpoints the transition from stable flow to unstable aortic fluttering as a function of the physiological properties composing the dimensionless number, like blood pressure and aortic compliance. A retrospective study was then conducted with 4D-flow MRI data from 117 patients indicated for cardiac imaging and 100 healthy volunteers recruited prospectively. The difference between the dimensionless number and its critical threshold was calculated for every subject from their earliest MRI data and used as an aneurysm physiomarker to forecast future growth. As a binary predictor for abnormal growth and subsequent surgical intervention reported from follow-up imaging, the aneurysm physiomarker yielded an AUC of 0.997 in a receiving operator characteristic analysis. Though validated here for thoracic ascending aortic aneurysms, this instability mechanism may be used to understand, predict and inform patient-specific treatment of aneurysms in any location without fundamental differences.


title: Novel Real-Time EMT-TS Modeling Architecture for Feeder Blackstart Simulations

abstract: This paper presents the development and benchmarking of a novel real-time electromagnetic-transient and transient-stability (EMT-TS) modeling architecture for distribution feeder restoration studies. The work presents for the first time in the literature a real-time EMT-TS testbed in which the grid-forming unit is simulated in EMT domain, operating as the slack bus of the phasor domain while including unbalanced voltage conditions. To evaluate the performance and limitations of the proposed model, an equivalent EMT testbed is developed to serve as a benchmark. First, the co-simulation framework and the domain coupling methodology are introduced. Next, the steady-state operation of the EMT-TS and EMT models are matched. Then, tests are conducted to analyze the transient performance of the proposed EMT-TS model when compared to its EMT benchmark. Results reveal that when utilizing a battery energy storage system (BESS) as the grid-forming unit, the EMT-TS testbed can maintain high accuracy for typical load steps.


title: Meta Clustering Learning for Large-scale Unsupervised Person Re-identification

abstract: Unsupervised Person Re-identification (U-ReID) with pseudo labeling recently reaches a competitive performance compared to fully-supervised ReID methods based on modern clustering algorithms. However, such clustering-based scheme becomes computationally prohibitive for large-scale datasets. How to efficiently leverage endless unlabeled data with limited computing resources for better U-ReID is under-explored. In this paper, we make the first attempt to the large-scale U-ReID and propose a "small data for big task" paradigm dubbed Meta Clustering Learning (MCL). MCL only pseudo-labels a subset of the entire unlabeled data via clustering to save computing for the first-phase training. After that, the learned cluster centroids, termed as meta-prototypes in our MCL, are regarded as a proxy annotator to softly annotate the rest unlabeled data for further polishing the model. To alleviate the potential noisy labeling issue in the polishment phase, we enforce two well-designed loss constraints to promise intra-identity consistency and inter-identity strong correlation. For multiple widely-used U-ReID benchmarks, our method significantly saves computational cost while achieving a comparable or even better performance compared to prior works.


title: Pricing S&P 500 Index Options with L\'evy Jumps

abstract: We analyze various jumps for Heston model, non-IID model and three L\'evy jump models for S&P 500 index options. The L\'evy jump for the S&P 500 index options is inevitable from empirical studies. We estimate parameters from in-sample pricing through SSE for the BS, SV, SVJ, non-IID and L\'evy (GH, NIG, CGMY) models by the method of Bakshi et al. (1997), and utilize them for out-of-sample pricing and compare these models. The sensitivities of the call option pricing for the L\'evy models with respect to parameters are presented. Empirically, we show that the NIG model, SV and SVJ models with estimated volatilities outperform other models for both in-sample and out-of-sample periods. Using the in-sample optimized parameters, we find that the NIG model has the least SSE and outperforms the rest models on one-day prediction.


title: Limiting absorption principle and singular spectrum

abstract: In this paper I give an explicit construction of an analogue of eigenspace for points of the singular spectrum of a self-adjoint operator. This construction is based on an abstract version of homogeneous Lippmann-Schwinger equation.


title: PT and Anti-PT Symmetry Arising in Time Dependent Weak Value Measurements

abstract: Parks introduced a formulation of time dependent weak values in 2008, which is the formalism we use in this paper. In this paper we extend notions from time dependent weak values to show that Hamiltonians associated with weak value measurements can be shown to exhibit even or odd symmetric properties. They exhibit PT or anti-PT symmetry, respectively. These symmetries are manifested during the measurement process as pointer translations, which have vanishing imaginary or vanishing real parts. The consequence of this that one can characterize some of the aspects of these symmetries of the time dependent Hamiltonians that arise from time dependent weak values. This allows one to generalize some work on non-Hermitian variables in quantum mechanics (QM) due to Bender related to weak values and weak measurement. We also speculate how these symmetries might apply to distinguishing between the various two-time interpretations of QM in the Conclusions of this paper.


title: Stellar populations in local AGNs: evidence for enhanced star formation in the inner 100pc

abstract: In modern models and simulations of galactic evolution, the star formation in massive galaxies is regulated by an ad hoc active galactic nuclei (AGN) feedback process. However, the physics and the extension of such effects on the star formation history of galaxies is matter of vivid debate. In order to shed some light in the AGN effects over the star formation, we analyzed the inner 500$\times$500pc of a sample of 14 Seyfert galaxies using GMOS and MUSE integral field spectroscopy. We fitted the continuum spectra in order to derive stellar age, metallicity, velocity and velocity dispersion maps in each source. After stacking our sample and averaging their properties, we found that the contribution of young SP, as well as that of AGN featureless continuum both peak at the nucleus. The fraction of intermediate-age SPs is smaller in the nucleus if compared to outer regions, and the contribution of old SPs vary very little within our field of view (FoV). We also found no variation of velocity dispersion or metallicity within our FoV. Lastly, we detected an increase in the dust reddening towards the center of the galaxies. These results lead us to conclude that AGN phenomenon is usually related to a recent star formation episode in the circumnuclear region.


title: Explaining GNN over Evolving Graphs using Information Flow

abstract: Graphs are ubiquitous in many applications, such as social networks, knowledge graphs, smart grids, etc.. Graph neural networks (GNN) are the current state-of-the-art for these applications, and yet remain obscure to humans. Explaining the GNN predictions can add transparency. However, as many graphs are not static but continuously evolving, explaining changes in predictions between two graph snapshots is different but equally important. Prior methods only explain static predictions or generate coarse or irrelevant explanations for dynamic predictions. We define the problem of explaining evolving GNN predictions and propose an axiomatic attribution method to uniquely decompose the change in a prediction to paths on computation graphs. The attribution to many paths involving high-degree nodes is still not interpretable, while simply selecting the top important paths can be suboptimal in approximating the change. We formulate a novel convex optimization problem to optimally select the paths that explain the prediction evolution. Theoretically, we prove that the existing method based on Layer-Relevance-Propagation (LRP) is a special case of the proposed algorithm when an empty graph is compared with. Empirically, on seven graph datasets, with a novel metric designed for evaluating explanations of prediction change, we demonstrate the superiority of the proposed approach over existing methods, including LRP, DeepLIFT, and other path selection methods.


title: From word-representable graphs to altered Tverberg-type theorems

abstract: Tverberg's theorem says that a set with sufficiently many points in $\mathbb{R}^d$ can always be partitioned into $m$ parts so that the $(m-1)$-simplex is the (nerve) intersection pattern of the convex hulls of the parts. In arXiv:1808.00551v1 [math.MG] the authors investigate how other simplicial complexes arise as nerve complexes once we have a set with sufficiently many points. In this paper we relate the theory of word-representable graphs as a way of codifying $1$-skeletons of simplicial complexes to generate nerves. In particular, we show that every $2$-word-representable triangle-free graph, every circle graph, every outerplanar graph, and every bipartite graph could be induced as a nerve complex once we have a set with sufficiently many points in $\mathbb{R}^d$ for some $d$.


title: Spatio-Temporal Modeling for Flash Memory Channels Using Conditional Generative Nets

abstract: We propose a data-driven approach to modeling the spatio-temporal characteristics of NAND flash memory read voltages using conditional generative networks. The learned model reconstructs read voltages from an individual memory cell based on the program levels of the cell and its surrounding cells, as well as the specified program/erase (P/E) cycling time stamp. We evaluate the model over a range of time stamps using the cell read voltage distributions, the cell level error rates, and the relative frequency of errors for patterns most susceptible to inter-cell interference (ICI) effects. We conclude that the model accurately captures the spatial and temporal features of the flash memory channel.


title: Recursive one-way Navier Stokes equations with PSE-like cost

abstract: Spatial marching methods, in which the flow state is spatially evolved in the downstream direction, can be used to produce low-cost models of flows containing a slowly varying direction, such as mixing layers, jets, and boundary layers. The parabolized stability equations (PSE) are popular due to their extremely low cost but can only capture a single instability mode; all other modes are damped or distorted by regularization methods required to stabilize the spatial march, precluding PSE from properly capturing non-modal behavior, acoustics, and interactions between multiple instability mechanisms. The one-way Navier-Stokes (OWNS) equations properly retain all downstream-traveling modes at a cost that is a fraction of that of global methods but still one to two orders of magnitude higher than PSE. In this paper, we introduce a new variant of OWNS whose cost, both in terms of CPU time and memory requirements, approaches that of PSE while still properly capturing the contributions of all downstream-traveling modes. The method is formulated in terms of a projection operator that eliminates upstream-traveling modes. Unlike previous OWNS variants, the action of this operator on a vector can be efficiently approximated using a series of equations that can be solved recursively, i.e., successively one after the next, rather than as a coupled set. In addition to highlighting the improved cost scaling of our method, we derive explicit error expressions and elucidate the relationship with previous OWNS variants. The properties, efficiency, and accuracy of our method are demonstrated for both free-shear and wall-bounded flows.


title: Embeddings and labeling schemes for A*

abstract: A* is a classic and popular method for graphs search and path finding. It assumes the existence of a heuristic function $h(u,t)$ that estimates the shortest distance from any input node $u$ to the destination $t$. Traditionally, heuristics have been handcrafted by domain experts. However, over the last few years, there has been a growing interest in learning heuristic functions. Such learned heuristics estimate the distance between given nodes based on "features" of those nodes. In this paper we formalize and initiate the study of such feature-based heuristics. In particular, we consider heuristics induced by norm embeddings and distance labeling schemes, and provide lower bounds for the tradeoffs between the number of dimensions or bits used to represent each graph node, and the running time of the A* algorithm. We also show that, under natural assumptions, our lower bounds are almost optimal.


title: Functional equivariance and conservation laws in numerical integration

abstract: Preservation of linear and quadratic invariants by numerical integrators has been well studied. However, many systems have linear or quadratic observables that are not invariant, but which satisfy evolution equations expressing important properties of the system. For example, a time-evolution PDE may have an observable that satisfies a local conservation law, such as the multisymplectic conservation law for Hamiltonian PDEs. We introduce the concept of functional equivariance, a natural sense in which a numerical integrator may preserve the dynamics satisfied by certain classes of observables, whether or not they are invariant. After developing the general framework, we use it to obtain results on methods preserving local conservation laws in PDEs. In particular, integrators preserving quadratic invariants also preserve local conservation laws for quadratic observables, and symplectic integrators are multisymplectic.